Gantry Type Coordinate Measuring Machines 27 frequencies of a CMM. Vermeulen [5] generated highaccuracy 3-D coordinate machines using a new conguratio
|
|
- Bernard Brett Marshall
- 5 years ago
- Views:
Transcription
1 Scientia Iranica, Vol. 13, No. 2, pp 26{216 c Sharif University of Technology, April 26 Research Note Free Vibration Analysis of Gantry Type Coordinate Measuring Machines M.T. Ahmadian, G.R. Vossoughi 1 and S. Ramezani 1 Coordinate Measuring Machines (CMMs are designed for precision inspection of complex industrial products. Even though CMMs have high accuracy in electrical systems, mechanical structures prevent very high accuracy in these machines. Mechanical accuracy of CMMs depends on both static and dynamic sources of error. In automated CMMs, one of the dynamic error sources is the vibration of probes, due to inertia forces resulting from parts deceleration. Modeling of a gantry type CMM, based on the Timoshenko beam theory, is developed and the natural frequencies of the CMM system, at dierent positions of the probe, are calculated. Findings from the analytical and nite elements method indicate high accuracy and good agreement between the results. INTRODUCTION Coordinate measuring machines are, nowadays, widely used for a large range of measuring tasks. These tasks are expected to be carried out with ever increasing accuracy, speed, exibility and the ability to operate under shop oor conditions. Research is necessary to meet these demands. CMMs are prone to many error sources. Based on functional components of a CMM, an overview has been given by Weekers [1] of the most important error sources aecting the accuracy of a CMM: Geometric errors: Limited accuracy in manufacturing, assembling and adjustment of components, like guide ways and measurement systems; Drive system: For CNC operated CMMs, the axes are equipped with drives, transmission and a servocontrol unit causing errors such as mechanical load and structural vibration; Measurement system: The actual coordinates of measuring points are derived from the values indicated by the linear scales of the CMM. The min. errors introduced by the scales are: Inaccuracy of *. Corresponding Author, Department of Mechanical Engineering, Sharif University of Technology, Tehran, I.R. Iran. 1. Department of Mechanical Engineering, Sharif University of Technology, Tehran, I.R. Iran. scale pitch, misalignment and adjustment of reading device and interpolation errors; Errors due to mechanical loads: These are errors related to static or slowly varying forces on CMM components in combination with the compliance of components and are mainly caused by the weight of moving parts; Thermally induced errors: The dierence between the temperature of the measuring scales of CMM and work piece and the temperature gradient in the machine components is a source of error; Dynamic errors: These are errors mainly caused by deceleration of moving parts before stopping. These errors depend on the CMMs structural properties, like mass distribution, component stiness and damping characteristics, as well as on the control system and disturbing forces. In this paper, the authors are interested in the dynamic errors of CMMs caused by the deceleration of moving parts. Some studies on the error analysis of CMMs have been done. Weekers and Schellekenes [2] proposed a method for compensation of the dynamic errors of CMMs using inductive position sensors for online measurement of the major dynamic errors. Barakat et al. [3] presented a kinematical and geometrical error compensation for CMMs based on experiment. Nijs et al. [4] presented a very simple model of CMM for obtaining natural
2 Gantry Type Coordinate Measuring Machines 27 frequencies of a CMM. Vermeulen [5] generated highaccuracy 3-D coordinate machines using a new conguration with fewer dynamic errors. Several researchers have applied software compensation successfully on CMMs [6,7]. Most of these researchers have considered a very simple model for their analysis while most of the studies are based on experiment. In the present study, a full CMM modeling is analyzed. All columns and guide ways are modeled as a Timoshenko beam [8,9] with exibility in all directions. All bearings are modeled as torsional springs and torsional deformation of the column and guide ways is considered. Dynamic equations of motion are derived using Hamilton's principle [1]. The derived equations are treated using a state space variables method. Natural frequencies of the CMM system are calculated at dierent positions of the probe in the system. Dynamic error analysis of CMM and optimization using a Genetic Algorithm, based on maximization of a tness function containing the eect of natural frequency, dynamic error and weight of the CMM structure, is presented by Ramezani [11]. MODELING OF CMM STRUCTURE Structural components of a gantry type CMM are shown in Figure 1. Because of the very small deformation of each component, the x-, y- and z-axes are assumed to remain in the same direction as in the undeformed state. The origin of the coordinate system is along the axis of the left column and at a height which is equal to the axis of the x-guideway. The most important problem in this modeling is that the motions of the y-guide way, x-guide way and z- pinole are relative motions. They are the absolute deformation of each beam, but, they are not the absolute motion of each beam. So, the strain energy of each beam is only a function of its deformations, however, the kinetic energy of each component is inuenced by the motion of other components. The system is assumed to be xed at each position of the probe. The modeling of each component is presented in the following sections. MODELING OF THE RIGHT COLUMN For the right column, bending, torsional and longitudinal vibration is considered. The kinetic and strain energy of the column can be written as: T c = 1 2 ( h Z ha 1 ( A_ 2 (;t+ _B 2 (;t+ _C 2 (;t + I 1z _ 2 (;t+i 1y 2 _ (;t+i 1z _ 2 1(;tid U c = 1 2 ( h Z( ; (1 E 1 I 1x 2 (;t+e 1 I 1y 2 (;t+kg 1 A + G 1 I 1z 2 1(;t+E 1 A 1 C 2 (;t d ; (2 - A(;t: Bending of column in xz-plane, - (;t: Rotation of column around y-axis, - B(;t: Bending of column in yz-plane, - (;t: Rotation of column around x-axis, - 1 (;t: Torsion of column around z-axis, - C(;t: Longitudinal vibration in z-direction, - A 1, I 1x, I 1y, I 1z, E 1, G 1, 1, h and k: Cross sectional area, actual and equivalent (for thin walled column moments of inertia, Young modulus, shear modulus, density, length of column and Timoshenko shear modication constant, respectively. Figure 1. Schematic view of a gantry type CMM. Note that the symbols ( : and ( denote derivation with respect to time and coordinates, respectively.
3 28 M.T. Ahmadian, G.R. Vossoughi and S. Ramezani MODELING OF THE RIGHT Y -GUIDE WAY For the right y-guide way, kinetic and strain energy have been considered to be as follows: T y = 2 2 ( Zys na 2 ([ _E(y; t+y _(h; t _C(h; t] 2 + _A 2 (h; t+[_d(y; t+ _B(h; t+y _ 1 (h; t] 2 + I 2z _ 2 1(y; t+i 2y [ _ 3 (y; t+ _ (h; t] 2 + I 2x _ 2 2(y; tody + ZLy ys na 2 ([ _E 1 (y; t+y _(h; t _C(h; t] 2 + _A 2 (h; t+[ _D 1 (y; t+ _B(h; t +y _ 1 (h; t] 2 +I 2x _ 2 22(y; t+i 2y [ _ 3 (y; t+ _ (h; t] 2 + I 2z _ 2 11(y; t+i 2y [ _ 3 (y; t+ _ (h; t] 2 + I 2x _ 2 2(y; tody U y = 1 2 ( Zys ; (3 ne 2 I 2z 2 1(y; t+e 2 I 2x 2 2(y; t +kg 2 A 2h 2 (y; t ZLy ys ne 2 I 2z 2 + kg 2 A 2h t + h 1 (y; + G 2 I 2y 2 3(y; tody 11(y; t+h 22 (y; 1(y; 2 22(y; t+g 2 I 2y 2 33(y; tody 1(y; t + E 2 I 2x ; (4 - y s : Location of y-carriage on the y-guide way, - D(y; t: Bending in xy-plane for <y<y s, - D 1 (y; t: Bending in xy-plane for y s < y<l y, - 1 (y; t: Rotation around z-axis <y<y s, - 11 (y; t: Rotation around z-axis y s < y<l y, - E(y; t: Bending in yz-plane <y<y s, - E 1 (y; t: Bending in yz-plane y s < y<l y, - 2 (y; t: Rotation around x-axis <y<y s, - 22 (y; t: Rotation around x-axis y s < y<l y, - 3 (y; t: Torsion around y-axis <y<y s, - 33 (y; t: Torsion around y-axis y s < y<l y, - A 2, I 2x, I 2z, I 2y, E 2, G 2, 2 and L y : cross sectional area, actual and equivalent moments of inertia, Young modulus, shear modulus and density, length of y-guide way, respectively. MODELING OF THE Y -CARRIAGE AND BEARING For the drive system eect, only the rotation of a y- carriage bearing around the y-axis is considered. The kinetic and strain energy in this member can be written as: T bc3 = 1 2 M bc3n[ _ A(h; t] 2 + [y s _ 1 (h; t+ _D(y s ; t + _B(h; t] 2 + [y s _(h; t+ _E(y s ; t _C(h; t] 2o + J bc3z _ 2 (t 2 + J bc3y [ _ 3 (y s ; t+ _ (y s ; t] 2 ; (5 U bc3 = 1 2 K bc3 2 (t 2 ; (6-2 (t: Torsion around y-axis - K bc3, M bc3, J bc3y and J bc3z : Torsional stiness of bearing and drive system, mass and mass moments of inertia of y-carriage and its bearing, respectively. MODELING OF THE X-GUIDE WAY For the x-guide way, bending and torsion kinetic and strain energy can be written as: T x = 3 2 ( Zxs na 3 ([ V _ (x; t+ A(h; _ t x _ 2 (t] 2 + [y s _ 1 (h; t+ _D(y s ; t+ _B(h; t] 2 +[ _W (x; t _C(h; t + _E(y s ; t+y s _(h; t+x _ 3 (x; t] 2 + I 3y _ 2 (x; t + I 3z _' 2 (x; t+i 3x _ 2 (x; todx + ZLx xs na 3 ([ _D(y s ; t
4 Gantry Type Coordinate Measuring Machines 29 + _B(h; t+y s _ 1 (h; t] 2 + [ _W 1 (x; t+y s _(h; t + I 3z _' 2 1(x; t+i 3x _ x(x; to 2 dx + _E(y s ; t _C(h; t+x _ 3 (x; t] 2 + I 3y _ 2 1(x; t U x = 1 2 ( Zxs ne 3 [I 3y 2 (x; t+i 3z' 2 (x; t] +G 3 I 3x 2 (x; t+kg 3 A 3h (x; t +h'(x; (x; Z ; (x; Lx odx+ ne 3 I 3y[ 2 1(x; t +I 3z ' 2 1(x; t]+g 3 I 3x 2 x(x; t+kg 3 A 3h' 1 (x; 1(x; +h 1 (x; t xs (x; 1(x; t - x s : Location of y-carriage on the y-guide way, - W (x; t: Bending in xz-plane for <x<x s, - W 1 (x; t: Bending in xz-plane for x s < x<l x, - (x; t: Rotation around y-axis <x<x s, - 1 (x; t: Rotation around y-axis x s < x<l x, - V (y; t: Bending in xy-plane <x<x s, - V 1 (y; t: Bending in xy-plane x s < x<l x, - '(x; t: Rotation around z-axis <x<x s, - ' 1 (x; t: Rotation around z-axis x s < x<l x, - (x; t: Torsion around x-axis <x<x s, - x (x; t: Torsion around x-axis x s < x<l x, (8 - A 3, I 3x, I 3y, I 3z, E 3, G 3, 3 and L x are area, actual and equivalent moments of inertia, Young modulus, shear modulus, density and length of x-guide way, respectively. MODELING OF THE BEARINGS AT Z-PINOLE The kinetic and strain energy resulting from the motion of the z-pinole and the torsion of bearings is as follows: T bcz = 1 2nM bczn[y s _ 1 (h; t+ _D(y s ; t+ _B(h; t] 2 + [y s _(h; t+ _E(y s ; t _C(h; t] 2 + [ _ A(h; t + _ V (xs ; t x s _ 2 (t] 2o+J bc4x [ _ 4 (t+ _(x s ; t] 2 +J bc5x [ _ 4 (t+ _ 5 (t+ _(x s ; t] 2 +J bc6y _ 2 6(to; (9 U bcz = 1 2 [K bc4 4 (t 2 + K bc5 5 (t 2 + K bc6 6 (t 2 ]; (1-4 (t: Torsion of x-carriage bearing around x-axis, - 5 (t: Torsion of z-axis bearing around x-axis, - 6 (t: Torsion of z-axis bearing around y-axis, - K bc4, K bc5, K bc6, M bcz, J bc4x, J bc5x and J bc6y : Torsional stiness of bearings and drive system, mass and mass moments of inertia of z-pinole assembly, respectively. MODELING OF THE Z-PINOLE Bending of the z-pinole in the xz- and yz-planes is considered. Kinetic and strain energy can be written as follows: T z = 4 2 A 4 ZLz ( [ _S(z;t+ _D(y s ; t+y s _ 1 (h; t + _B(h; t+(z + z _ 6 (t] 2 + [ _R(z;t+ _A(h; t + (z + z ( _ 4 (t+ _ 5 (t+ _(x s ; t + _ V (xs ; t x s _ 2 (t] 2 + [ _W (x s ; t+ _E(y s ; t+y s _(h; t _C(h; t+x s _ 3 (x; t+l( _ 4 (t+ _ 5 (t + _(x s ; t] 2 + I 4x _ 2 1(z;t+I 4y _ 2 2(z;t dz; (11
5 21 M.T. Ahmadian, G.R. Vossoughi and S. Ramezani U x = 1 2 ZLz ne 4 I 4x 2 1(z;t+E 4 I 4x 2 2(z;t + kg 4 A +h( 2 odz; EQUATIONS OF MOTION Using Hamilton's principle, equations of motion can be found. Zt2 t1 (T U + Wdt =; (15 where denotes the variation operator and: T = T c + T x + T y + T z + T bcz + T L ; (16 - R(z;t: Bending in yz-plane, - 1 (z;t: Rotation around x-axis, - S(z;t: Bending in xz-plane, - 2 (z;t: Rotation around y-axis - A 4, I 4x, I 4y, I 4z, E 4, G 4, 4 and L z : Cross sectional area, moments of inertia, Young modulus, shear modulus, density and length of z-axis, respectively. MODELING OF LEFT COLUMN AND Y -GUIDE WAY For the left-side support column and y-guide way, the same motions C(;t, E(y; t and 2 (y; t have been considered in a similar fashion to the right side column and y-guide way. The kinetic and strain energy can be written as follows: T L = 1 2 n 5 A 5 Z h _C 2 (;td + 6 A 6 ZLy [ _E(y; t _C(h; t] 2 dyo; (13 U L = 1 2 ne 5 A 5 Z h C 2 (;td + + kg 6 A 6h 2 (y; t Zys ne 6 I 6x 2 2(y; t odyo; where A 5, A 6, I 6x, E 5, E 6, G 5, G 6, 5 and 6 are cross sectional area, moments of inertia, Young modulus, shear modulus and density of left column and left y- guide way, respectively. Note that because of the eect of the drive system, some motions, such as the rotation of the x-carriage around the z-axis and the longitudinal vibration of the z-axis, the horizontal and vertical motion of the y-carriage have been neglected. Furthermore, the gravitational strain energy is neglected. U = U c + U x + U y + U z + U bcz + U L : (17 In the case of a xed system, neglecting gravitational eects, it can be written that: W =: (18 It is interesting to nd natural frequencies using a state space variables approach. State variables are dened as follows: A = x 1 ; A = x 2 ; = x 3 ; = x 4 ; B = x 5 ; B = x 6 ; = x 7 ; = x 8 ; C = x 9 ; C = x 1 ; 1 = x 11 ; = x 1 12 ; D = x 13 ; D = x 14 ; 1 = x 15 ; = x 1 16 ; E = x 17 ; E = x 18 2 = x 19 ; = x 2 2 ; 3 = x 21 ; = x 3 22 ; D 1 = x 23 ; D = x 1 24 ; 11 = x 25 ; = x 1 26 ; E 1 = x 27 ; E = x 1 28 ; 22 = x 29 ; = x 22 3 ; 33 = x 31 ; = x ; W = x 33 ; W = x 34 ; = x 35 ; = x 36 ; V = x 37 ; V = x 38 ; ' = x 39 ; ' = x 4 ; = x 41 ; = x 42 ; W 1 = x 43 ; W = x 1 44 ; 1 = x 45 ; = x 1 46 ; V = x 47 ; V = x 48 ; ' 1 = x 49 ; ' = x 1 5 ; x = x 51 ; x = x 52 ; R = x 53 ; R = x 54 ; 1 = x 55 ; = x 1 56 ; S = x 57 ; S = x 58 ; 2 = x 59 ; = x 2 6 : (19 Note that for free vibration, each state variable can be separated into time and spatial coordinates, so: x i (coordinates, time = X i (coordinates:e i!t : (2 Now, the equations and boundary conditions of the system are written. Note that the attractiveness of the energy method is that it gives both equations of motion and corresponding boundary conditions. EQUATIONS OF THE RIGHT COLUMN X 1( =X 2 (; X 2( =X 4 ( ( 1! 2 =kg1x 1 (; (21a (21b
6 Gantry Type Coordinate Measuring Machines 211 X3( =X 4 (; X4( = ( 1! 2 =E1X 3 (+(kg 1 A 1 =E 1 I 1x [X 3 ( X 2 (]; X5( =X 6 (; X6( =X 8 ( ( 1! 2 =kg1x 5 (; X7( =X 8 (; X8( = ( 1! 2 =E1X 7 (+(kg 1 A 1 =E 1 I 1x (21c (21d (22a (22b (22c EQUATION CORRESPONDING TO ROTATION 2! 2 [J bc3z +M bcz x 2 s +A 3 L 3 x 3=3+x 2 s A 4L z 4 ]K bc3 2! 2 x s A 4 4 [(L z + z 2 z 2 ]=2( 4 + 5! 2 [M bcz x s +A 3 L 2 x 3 =2+x s A 4 L z 4 ]X 1 (h! 2 x s [M bcz + A 4 L z 4 ]X 27 (y s A 4 4 [(L z + z 2 z 2 ]=2X 31 (x s =: (28 [X 3 ( X 2 (]; (22d X9( =X 1 (; (23a X1( =! 2 ( 1 A A 5 [E 1 A 1 + E 5 A 5 ] 1 X 9 (; (23b X11( =X 12 (; (24a X12( =! 2 1 =G 1 X 9 (; (24b EQUATIONS OF RIGHT Y -GUIDE WAY X13(y =X 14 (y; (25a X14(y =X 16 (y! 2 2 A 2 [k(g 6 A 6 + G 2 A 2 ] 1 [X 13 (y+yx 11 (h+x 5 (h]; (25b X15(y =X 16 (y; (25c X16(y =[E 6 I 6z + E 2 I 2z ] 1 f! 2 2 I 2z X 15 (y + k(g 6 A 6 + G 2 A 2 [X 15 (y X 14 (yg; (25d EQUATIONS OF X-GUIDE WAY X 23(x =X 24 (x; X 24(x =X 26 (x! 2 3 =kg 3 fx 23 (x+x 17 (y s (29a + y s X 3 (h X 9 (h]g; (29b X 25(x =X 26 (x; X 26(x =1=E 3 I 3y f! 2 3 I 3y X 25 (x+kg 3 A 3 [X 25 (x X 24 (x]g; X 27(x =X 28 (x; (29c (29d (3a X 28(x=X 3 (x! 2 3 =kg 3 [X 27 (x x 2 +X 1 (h]; (3b X 29(x =X 3 (x; X 3(x =1=E 3 I 3z f! 2 3 I 3z X 29 (x (3c + kg 3 A 3 [X 29 (x X 28 (x]g; (3d X 17(y =X 18 (y; (26a X 31(x =X 32 (x; (31a X 18(y=X 2 (y! 2 [k(g 6 A 6 +G 2 A 2 ] 1 f( 2 A 2 X 32(x =! 2 3 =G 3 X 31 (x; (31b + 6 A 6 [X 17 (y X 9 (h] + 2 A 2 yx 3 (hg; (26b X 19(y =X 2 (y; X 2(y =[E 6 I 6x + E 2 I 2x ] 1 f! 2 [ 2 I 2x + 6 I 6x ] (26c X 19 (y+k(g 6 A 6 +G 2 A 2 [X 19 (y X 18 (y]g; (26d EQUATIONS CORRESPONDING TO THE ROTATIONS 4, 5 AND 6 fk bc4! 2 [J bc4x + J bc5x + A 4 4 (L z l 2 + [(L z + z 3 z 3 ]=3]g 4! 2 fj bc5x +A 4 4 (L z l 2 + [(L z +z 3 z 3 ]=3g 6 +x s (L z +z 2 z 2 ]=2! 2 2 X 21(y =X 22 (y; X 22(y =! 2 2 =G 2 [X 21 (y+x 7 (h]: (27a (27b =! 2 [J bc4x +J bc5x +A 4 4 (L z l 2 +[(L z + z 3 z 3 ]=3]X 31 (x s +! 2 A 4 4 f[(l z + z 2
7 212 M.T. Ahmadian, G.R. Vossoughi and S. Ramezani z 2 ]=2][X 1 (h+x 27 (x s ] +[L z l 2 X 27 (x s + L z ly s X 3 (h+lx s X 21 (y s + L z lx 17 (y s + L z lx 9 (hg; (32 A 4 4 x s [(L z +z 2 z 2 ]=2]! 2 2 fj bc5x +A 4 4 (L z l 2 + [(L z + z 3 z 3 ]=3]g! 2 ( K bc5 5 =! 2 fj bc5x +A 4 4 (L z l 2 +[(L z + z 3 z 3 ]=3]g +! 2 A 4 4 f[(l z +z 2 z 2 ]=2][X 1 (h+x 27 (x s ] +L z l 2 (y s X 3 (h+x s X 21 (h X 9 (h+x 23 (x s ]g; (33 f! 2 (J bc6y + A 4 4 [(L z + z 2 z 2 ]=2] K bc6y g 6 =! 2 A 4 4 [(L z + z 2 z 2 ]=2][y s X 11 (h + X 13 (y s + X 5 (h]: (34 EQUATIONS IN Z-PINOLE X 53(z =X 54 (z; X 54(z =X 56 (z! 2 4 =kg 4 fx 53 (z+(z + z [X 41 (x s ] x s 2 + X 1 (h (35a + X 37 (h]g; (35b X 55(z =X 56 (z; (35c BOUNDARY CONDITIONS FOR LEFT COLUMN The left column is considered to be clamped at =. At = h, the boundary conditions are determined from Hamilton's principle as follows: X 1 ( = X 3 ( = X 5 ( = X 7 ( = X 9 ( = X 11 ( = ; [2A 2 L y 2 +M]X 1 (h+a 3 L 2 x 3 =2 2 kg 1 A 1 [X 3 (h X 2 (h] A 4 4 [(L z + z 2 z 2 ]=2X 31 (x s = ; [2A 2 L y 2 +M]X 1 (h+a 3 L 2 x 3=2 2 kg 1 A 1 [X 3 (h X 2 (h] A 4 4 [(L z + z 2 z 2 ]=2X 31 (x s =; [ 2 A 2 L 3 y=3+y 2 sm]x 3 (h E 1 I 1x X 4 (h[ 2 A 2 L 2 y=2 + y s M]X 9 (h+ y s MX 17 (h+y s [x s M bcz +A 3 3 L 2 x=2 + A 4 L z 4 x s ]X 21 (y s + y s [M bcz + A 4 L z 4 ]X 23 (x s +y s [x s M bcz +A 3 3 L 2 x =2+A 4L z 4 x s ]X 31 (x s = ; [ 2 A 2 L 2 y=2+y 2 sm]! 2 X 5 (h+kg 1 A 1 [X 7 (h X 6 (h] [M M bc3 ]! 2 X 13 (y s +[ 2 A 2 L 2 y=2 + y s (M M bc3 ]! 2 X 11 (h X56(z =1=E 4 I 4x f! 2 4 I 4x X 55 (z+kg 4 A 4 [X 55 (z X 54 (z]g; X57(z =X 58 (z; X 58(z =X 6 (z! 2 4 =kg 4 [X 57 (z+(z + z 6 (35d (36a = ;! 2 [ 2 I 2y L y +J bc3y ]X 7 (h+e 1 I 1y X 7 (h+j bc3y X 21 (y s = ;! 2 [2 2 A 2 L y + 6 A 6 L y +M]X 9 (E 1 A 1 +E 5 A 5 X 1 (h + y s X 11 (h+x 5 (h+x 13 (y s ]; (36b! 2 M[y s X 3 (h+x 17 (y s ] X 59(z =X 6 (z; X 6(z =1=E 4 I 4y f! 2 4 I 4y X 59 (z+kg 4 A 4 [X 59 (z X 58 (z]g: (36c (36d! 2 (M bcz +A 4 L z 4 [X 23 (x s +x s X 21 (y s ]! 2 la 4 L z 4 lx 1 (x s =; [3 2 A 2 L 2 y=2+y s M]X 5 (h+g 1 I 1z X 12 (h Now, the corresponding boundary conditions for each member are written. [2 2 A 2 L 3 y =3+y2 s M]X 11(h y s MX 11 (h =; (37
8 Gantry Type Coordinate Measuring Machines 213 where M =(M bc3 + M bcz + A 3 L x 3 + A 4 L z 4. BOUNDARY CONDITIONS OF Y -GUIDE WAY Because the motions of the y-guide way are relative to the left column, it can be assumed that it is clamped at y = h and free at y = L y. Continuity conditions must be satised at y = y s. X 13 ( = X 15 ( = X 26 (L y = X 17 ( =; X 19 ( = X 3 (L y = X 21 ( = X 32 (L y =; X 13 (y s X 23 (y s = X 24 (L y X 25 (L y =; X 15 (y s X 25 (y s = X 16 (y s X 26 (y s =; X 17 (y s X 27 (y s = X 28 (L y X 29 (L y =; X 5 (h+y s X 11 (h+x 13 (y s =X 21 (y s X 31 (y s =; M[y s X 3 (h X 9 (h]+[m bcz +A 3 3 L 2 x=2+a 4 L z 4 ] X 21 (y s +(M bcz + A 4 L z 4 X 33 (x s + la 4 L z 4 X 31 (x s =; X 19 (y s X 29 (y s = X 2 (y s X 3 (y s =; A 4 L z 4 lx s X 31 (x s + x s [M bcz + A 4 L z 4 ] [X 23 (x s + X 9 (h+x s X 17 (y s + y s X 3 (h]=: (38 BOUNDARY CONDITIONS OF THE X-GUIDE WAY The x-guide way is assumed to be a simply supported beam in the xz-plane and to be clamped-free in the xy-plane. It's torsion along the x-axis is assumed to be clamped at both ends. Continuity conditions must be satised at x = x s. X 33 (=X 43 (L x =X 35 (=X 36 (L x =X 37 (=; X 39 ( = X 5 (L x = X 41 ( = X 51 (L x =; X 33 (x s X 43 (x s = X 35 (x s X 45 (x s =; X 36 (x s X 46 (x s = X 37 (x s X 47 (x s =; la 4 L z 4 X 41 (x s +[M bcz + A 4 L z 4 ][X 17 (y s + X 33 (x s + y s X 3 (h X 9 (h]=; X 49 (L x X 48 (L x = X 39 (x s X 49 (x s =; X 4 (x s X 6 (x s = X 41 (x s X 51 (x s =; [J bc4x +J bc5x +A 4 4 (L z l 2 +[(L z +z 3 z 3 ]=3] X 31 (x s +A 4 4 [(L z +z 2 z 2 ]=2][X 1 (h + X 37 (x s ] + A 4 4 L z l[ X 9 (h+x 33 (x s + y s X 3 (h]=: (39 BOUNDARY CONDITIONS OF THE Z-PINOLE The motion of the z-pinole in the zx- and zy-planes is assumed to be clamped-free. X 53 ( = X 55 ( = X 56 (L z = X 57 ( = X 59 ( = X 6 (L z =; X 54 (L z X 55 (L z = X 58 (L z X 59 (L z =: (4 METHOD OF SOLUTION In order to obtain the natural frequencies of the system, initially, the rotations 2, 4, 5 and 6 will be found from Equations 27 and 32 to 34 and substituted in the model equations. Note that the derived equation can be listed into two categories. The rst category contains equations corresponding to the motions of the left column and may be written in the following form: fx (g =[D]fX(g; (41 where fxg = fx i ; X j ; :::;X n g t represents the vector of corresponding state variables. The sign ( denotes derivation, with respect to. Furthermore, the matrix, [D] is in the form: [D] =[K]! 2 [M]; (42 where [K] and[m] can be found readily by arranging the equations in terms of! 2. The solution of Equation 41 can be written as: fx g =[U] t [][U]fCg; (43 where [U] is the matrix of the eigenvectors of [D] and fcg is a column vector of constants. The diagonal matrix, [], is in the form:
9 214 M.T. Ahmadian, G.R. Vossoughi and S. Ramezani [] = e i ; (44 where i 's are eigenvalues of the dynamic matrix [D] depending on! 2. Taking = h, the value of X i (h in terms of C i 's can be found and will be used in other formulas. The second category of equations contains equations corresponding to the motion of members in the y- guide way, x-guide way and z-pinole. These equations are in the form: fx g =[D]fXg +[F ]; (45 where, in the above equation, the vector, fxg, may be a function of x, y or z. In order to solve the equation, initially, the matrix of the eigenvectors of [D], i.e. [U], should be found. Using the following transformation: fy g =[U]fXg; (46 it can be obtained: or: fy g =[U] t [D][U]fY g +[U] t [F ]; (47 fy g = []fy g +[U] t [F ]: (48 Because the matrix, [], is diagonal, an ordinary dierential equation for each Y i is obtained and can be solved readily. Note that for each Y i, there is a constant of integration, C i. Next, using Equation 46, each X i can be found. Taking y = y s and x = x s, the values of X i (x s or X i (y s can be found in terms of C i 's and implemented in other formulas. So far, X i 's are obtained in terms of C i 's and i 's. Now, substituting X i 's into boundary conditions and rearranging in terms of C i 's, the following equation may be obtained: [B]fCg = fg: (49 In order to have a nontrivial solution, the determinant of coecient matrix [B] must be zero, i.e.: Det[B] =: (5 In order to illustrate an example of modeling, a CMM with the following specications is modeled. i = 785 kg/m 3 ; E i = 2 Gpa; G i = 7 Gpa; h =:825 m; L y = :6 m; L x = 1:2 m; L z = :75 m; l =:16 m; z = :125 m; K bearings = 5e7 N.m/rad; A 1 = 7:6e 3m 2 ; A 2 = 4:6e 3m 2 ; A 3 = 7:6e 3m 2 ; A 4 = 1:e 3m 2 ; A 5 = 4:6e 3m 2 ; A 6 = 4:6e 3m 2 ; I 1x = 7:86e 5m 4 ; I 1y = 1:3e 5m 4 ; I 1z = 3:58e 5m 4 ; I 2x = 1:35e 5m 4 ; I 2y = 1:38e 6m 4 ; I 2z = 2:4e 5m 4 ; I 3x = 5:94e 5m 4 ; I 3y = 2:82e 5m 4 ; I 3z = 6:35e 5m 4 ; I 4x = 4:21e 7m 4 ; I 4y = 4:21e 7m 4 ; I 5x = 1:35e 5m 4 ; I 6x =1:35e 5m 4 ; I 6y =2:4e 5m 4 ; I 6z = 6:95e 5m 4 ; J bc3z = :18 kgm 4 ; J bc3z = :3 kgm 4 ; J bc4x = 2:7 kgm 4 ; J bc5x = :5 kgm 4 ; J bc6y = :5 kgm 4 ; M bc3 = 1 kg; M bcz = 35 kg: The natural frequencies of the system depend on the location of the probe. The variation of the 1st and 2nd natural frequencies versus the fxy g position of the probe is illustrated in the 3-D plots of Figures 2 and 3, when the tip of the probe is at its maximum distance from the x-guideway (z = L z. It is clear that increasing the values of x s and y s reduces the natural frequencies of the system. Note that the driving systems of the x- and y-motions are located at the right side of the CMM, consequently, the x-guide way slides on the left y-guide way, and the minimum value of the natural frequency occurs when the probe is located at the end of the x-guide way. The variation of the 1st to 4th natural frequencies of the system at dierent fxy g positions of the probe (z = L z, obtained from the This yields an equation in terms of! 2. Because the presented model is based on continuous beam modeling, the obtained equation is transcendental and has an innite number of solutions. For each value of x s and y s, it can solve for! 2. Figure 2. First natural frequency (z = L z.
10 Gantry Type Coordinate Measuring Machines 215 Figure 3. Second natural frequency (z = L z. analytical and nite elements method, are presented in Figures 4 to 7. Findings from the analytical and nite elements method indicate high accuracy and good agreement between the results. Figure 5. Natural frequencies in Y = 2 for X = to X = 12 cm. CONCLUSION The modeling of a gantry type CMM, based on Timoshenko beam theory using Hamilton's principle, was developed and the natural frequencies of the system were calculated at various positions of the probe in the system. The parameters of motion are dened so that the kinetic energy of each element is inuenced by the motion of adjacent connected elements, while the strain energy is independent of the position and velocity of these elements. Results indicate that in the rst four natural frequencies, for any y-position of the probe, the system has the lowest natural frequency when the probe moves toward the end of the x-guide way. Moreover, the values of natural frequencies again decrease when the system moves toward the end of the y-guide way. However, when the system is at the end of both x- and y-guide ways, the minimum value of the Figure 6. Natural frequencies in Y = 4 for X = to X = 12 cm. Figure 4. Natural frequencies in Y = for X =to X = 12 cm. Figure 7. Natural frequencies in Y = 5 for X = to X = 12 cm (Z = L z.
11 216 M.T. Ahmadian, G.R. Vossoughi and S. Ramezani natural frequencies occurs. This is due to the minimum stiness of the system in that position. These frequencies help us to present any resonant when the system is under operation. Any interference between the systems natural frequencies and the frequencies of the driving system should be prevented. A good agreement and high accuracy in comparison with the results found by the nite element analysis can be observed. REFERENCES 1. Weekers, W.G. \Compensation for dynamic errors of coordinate measuring machines", PhD Thesis, Eindhoven University of Technology ( Weekers, W.G. and Schellekens, P.H.J. \Compensation for dynamic errors of coordinate measuring machines", Measurement, 2, pp ( Barakat, N.A., Elbestawi, M.A. and Spence, A.D. \Kinematic and geometric error compensation of a coordinate measuring machine", Machine Tools & Manufacture, 4, pp (2. 4. De Nijs, J.F.C., Schellekens, M.G.M. and Van der Wolf, A.C.H. \Modeling of a coordinate measuring machine for analysis of its dynamic behavior", Annals of the CIRP, 37, pp ( Vermeulen, M. \High precision 3D-coordinate measuring machine, design and prototype development", PhD Thesis, Eindhoven University of Technology ( Teewwsen, J.W.M., Soons, J.A. and Schelekens, P.H. \A general method for error description of CMMs using polynomial tting procedure", Annals of the CIRP, 38, pp ( Zhang, G. \Error compensation of coordinate measuring machines", Annals of the CIRP, 34, pp ( Lee, H.P. \Transverse vibration of a Timoshenko beam acted on by an accelerating mass", Applied Acoustics, 47, pp ( Esmailzadeh, E. and Ghorashi, M. \Vibration analysis of a Timoshenko beam subjected to a traveling mass", Journal of Sound and Vibration, 199, pp ( Reddy, J.N. \Energy and variational methods in applied mechanics", John Wiley & Sons ( Ramezani, S. \Vibration analysis of gantry type CMMs and optimization using genetic algorithm", M.Sc. Thesis, Sharif University of Technology, Tehran, Iran (23.
Archive of SID. Archive of SID.
Scientia Iranica, Vol. 14, No. 3, pp 78{9 c Sharif University of Technology, June 7 Research Note Dynamic Error Analis of Gantry Type Coordinate Measuring Machines M.T. Ahmadian, G.R. Vossoughi 1 and S.
More informationDynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary Conditions
Transaction B: Mechanical Engineering Vol. 16, No. 3, pp. 273{279 c Sharif University of Technology, June 2009 Research Note Dynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationPARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationCork Institute of Technology. Summer 2007 Mechanics of Machines (Time: 3 Hours)
Cork Institute of Technology Bachelor of Engineering (Honours) in Mechanical Engineering- Award Instructions Answer FOUR questions. All questions carry equal marks. (NFQ Level 8) Summer 2007 Mechanics
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationHydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition
Fluid Structure Interaction and Moving Boundary Problems IV 63 Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition K.-H. Jeong, G.-M. Lee, T.-W. Kim & J.-I.
More informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More information34 G.R. Vossoughi and A. Karimzadeh position control system during unconstrained maneuvers (because there are no interaction forces) and accommodates/
Scientia Iranica, Vol. 14, No. 1, pp 33{45 c Sharif University of Technology, February 27 Impedance Control of a Flexible Link Robot for Constrained and Unconstrained Maneuvers Using Sliding Mode Control
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationFinite Element Analysis of Piezoelectric Cantilever
Finite Element Analysis of Piezoelectric Cantilever Nitin N More Department of Mechanical Engineering K.L.E S College of Engineering and Technology, Belgaum, Karnataka, India. Abstract- Energy (or power)
More informationEFFECT OF A CRACK ON THE DYNAMIC STABILITY OF A FREE}FREE BEAM SUBJECTED TO A FOLLOWER FORCE
Journal of Sound and
More informationApplications of Eigenvalues & Eigenvectors
Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More information19.2 Mathematical description of the problem. = f(p; _p; q; _q) G(p; q) T ; (II.19.1) g(p; q) + r(t) _p _q. f(p; v. a p ; q; v q ) + G(p; q) T ; a q
II-9-9 Slider rank 9. General Information This problem was contributed by Bernd Simeon, March 998. The slider crank shows some typical properties of simulation problems in exible multibody systems, i.e.,
More informationReview Lecture. AE1108-II: Aerospace Mechanics of Materials. Dr. Calvin Rans Dr. Sofia Teixeira De Freitas
Review Lecture AE1108-II: Aerospace Mechanics of Materials Dr. Calvin Rans Dr. Sofia Teixeira De Freitas Aerospace Structures & Materials Faculty of Aerospace Engineering Analysis of an Engineering System
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More information436 A. Barani and G.H. Rahimi assessment models have been employed to investigate the LBB of cracked pipes that are not for combined load [8]. Yun-Jae
Scientia Iranica, Vol. 4, No. 5, pp 435{44 c Sharif University of Technology, October 27 Approximate Method for Evaluation of the J-Integral for Circumferentially Semi-Elliptical-Cracked Pipes Subjected
More informationAdvanced Vibrations. Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Exact Solutions Relation between Discrete and Distributed
More informationUnit 13 Review of Simple Beam Theory
MIT - 16.0 Fall, 00 Unit 13 Review of Simple Beam Theory Readings: Review Unified Engineering notes on Beam Theory BMP 3.8, 3.9, 3.10 T & G 10-15 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More information142. Determination of reduced mass and stiffness of flexural vibrating cantilever beam
142. Determination of reduced mass and stiffness of flexural vibrating cantilever beam Tamerlan Omarov 1, Kuralay Tulegenova 2, Yerulan Bekenov 3, Gulnara Abdraimova 4, Algazy Zhauyt 5, Muslimzhan Ibadullayev
More informationMEC-E8001 FINITE ELEMENT ANALYSIS
MEC-E800 FINIE EEMEN ANAYSIS 07 - WHY FINIE EEMENS AND IS HEORY? Design of machines and structures: Solution to stress or displacement by analytical method is often impossible due to complex geometry,
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More information4. Mathematical models used in engineering structural analysis
4. Mathematical models used in engineering structural analysis In this chapter we pursue a formidable task to present the most important mathematical models in structural mechanics. In order to best situate
More informationTable 1: BEM as a solution method for a BVP dierential formulation FDM BVP integral formulation FEM boundary integral formulation BEM local view is ad
Chapter 1 Introduction to Boundary element Method - 1D Example For reference: Hong-Ki Hong and Jeng-Tzong Chen, Boundary Element Method, Chapter 1 Introduction to Boundary element Method - 1D Example,
More informationCRACK-TIP DRIVING FORCE The model evaluates the eect of inhomogeneities by nding the dierence between the J-integral on two contours - one close to th
ICF 100244OR Inhomogeneity eects on crack growth N. K. Simha 1,F.D.Fischer 2 &O.Kolednik 3 1 Department ofmechanical Engineering, University of Miami, P.O. Box 248294, Coral Gables, FL 33124-0624, USA
More informationName: Fall 2014 CLOSED BOOK
Name: Fall 2014 1. Rod AB with weight W = 40 lb is pinned at A to a vertical axle which rotates with constant angular velocity ω =15 rad/s. The rod position is maintained by a horizontal wire BC. Determine
More informationPerturbation Method in the Analysis of Manipulator Inertial Vibrations
Mechanics and Mechanical Engineering Vol. 15, No. 2 (2011) 149 160 c Technical University of Lodz Perturbation Method in the Analysis of Manipulator Inertial Vibrations Przemys law Szumiński Division of
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationwhere G is called the universal gravitational constant.
UNIT-I BASICS & STATICS OF PARTICLES 1. What are the different laws of mechanics? First law: A body does not change its state of motion unless acted upon by a force or Every object in a state of uniform
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationHong Kong Institute of Vocational Education (Tsing Yi) Higher Diploma in Civil Engineering Structural Mechanics. Chapter 2 SECTION PROPERTIES
Section Properties Centroid The centroid of an area is the point about which the area could be balanced if it was supported from that point. The word is derived from the word center, and it can be though
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationCHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS
61 CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS 4.1 INTRODUCTION The analysis of cantilever beams of small dimensions taking into the effect of fringing fields is studied and
More informationVIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS
VIBRATION AND DAMPING ANALYSIS OF FIBER REINFORCED COMPOSITE MATERIAL CONICAL SHELLS Mechanical Engineering Department, Indian Institute of Technology, New Delhi 110 016, India (Received 22 January 1992,
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationBOUNDARY CONDITION IDENTIFICATION BY SOLVING CHARACTERISTIC EQUATIONS
Journal of Sound and
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationContents. 2 Partial Derivatives. 2.1 Limits and Continuity. Calculus III (part 2): Partial Derivatives (by Evan Dummit, 2017, v. 2.
Calculus III (part 2): Partial Derivatives (by Evan Dummit, 2017, v 260) Contents 2 Partial Derivatives 1 21 Limits and Continuity 1 22 Partial Derivatives 5 23 Directional Derivatives and the Gradient
More informationModelling mechanical systems with nite elements
Modelling mechanical systems with nite elements Bastian Kanning February 26, 2010 Outline 1 Introduction 2 Virtual displacements 3 Getting started The ansatz functions 4 Mass and stiness elementwise 5
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationEE C245 ME C218 Introduction to MEMS Design
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 16: Energy
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationGroup Theory. 1. Show that Φ maps a conjugacy class of G into a conjugacy class of G.
Group Theory Jan 2012 #6 Prove that if G is a nonabelian group, then G/Z(G) is not cyclic. Aug 2011 #9 (Jan 2010 #5) Prove that any group of order p 2 is an abelian group. Jan 2012 #7 G is nonabelian nite
More informationSound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur
Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-13 WAVE PROPAGATION IN SOLIDS Longitudinal Vibrations In Thin Plates Unlike 3-D solids, thin plates have surfaces
More informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
More informationOn the use of the Cornu Method for determining Young s Modulus & Poisson s Ratio.
On the use of the Cornu Method for determining Young s Modulus & Poisson s Ratio. Jai Kai Chai Shane Duane ID: 08764522 SF Theoretical Physics November 11th, 2009 Fergal Haran Abstract The aim of this
More informationA displacement nite volume formulation for the static and dynamic analyses of shear deformable circular curved beams
Scientia Iranica A (018) 5(3), 999{1014 Sharif University of Technology Scientia Iranica Transactions A: Civil Engineering http://scientiairanica.sharif.edu A displacement nite volume formulation for the
More informationFree-edge stresses in general cross-ply laminates
Scientia Iranica B (204) 2(2), 387{402 Sharif University of Technology Scientia Iranica Transactions B: Mechanical Engineering www.scientiairanica.com Free-edge stresses in general cross-ply laminates
More informationFOR FINITE ELEMENT MODELS WITH APPLICATION TO THE RANDOM VIBRATION PROBLEM. Alexander A. Muravyov 1
DETERMINATION OF NONLINEAR STIFFNESS COEFFICIENTS FOR FINITE ELEMENT MODELS WITH APPLICATION TO THE RANDOM VIBRATION PROBLEM Alexander A. Muravyov 1 NASA Langley Research Center, Structural Acoustics Branch,
More informationComb resonator design (2)
Lecture 6: Comb resonator design () -Intro Intro. to Mechanics of Materials School of Electrical l Engineering i and Computer Science, Seoul National University Nano/Micro Systems & Controls Laboratory
More informationL bor y nnd Union One nnd Inseparable. LOW I'LL, MICHIGAN. WLDNHSDA Y. JULY ), I8T. liuwkll NATIdiNAI, liank
G k y $5 y / >/ k «««# ) /% < # «/» Y»««««?# «< >«>» y k»» «k F 5 8 Y Y F G k F >«y y
More informationDAMPING RATIO AND NATURAL FREQUENCY BIFURCATIONS IN ROTATING SYSTEMS
Journal of Sound and
More informationShape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression
15 th National Conference on Machines and Mechanisms NaCoMM011-157 Shape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression Sachindra Mahto Abstract In this work,
More informationCalculating the Applied Load
The LM Guide is capable of receiving loads and moments in all directions that are generated due to the mounting orientation, alignment, gravity center position of a traveling object, thrust position and
More informationCardan s Coupling Shaft as a Dynamic Evolutionary System
American Journal of Modern Physics and Application 2017; 4(2): 6-11 http://www.openscienceonline.com/journal/ajmpa Cardan s Coupling Shaft as a Dynamic Evolutionary System Petr Hrubý 1, Zdeněk Hlaváč 2,
More informationGravitational potential energy *
OpenStax-CNX module: m15090 1 Gravitational potential energy * Sunil Kumar Singh This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 2.0 The concept of potential
More informationSpecial Classes of Fuzzy Integer Programming Models with All-Dierent Constraints
Transaction E: Industrial Engineering Vol. 16, No. 1, pp. 1{10 c Sharif University of Technology, June 2009 Special Classes of Fuzzy Integer Programming Models with All-Dierent Constraints Abstract. K.
More informationFREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY
FREE VIBRATION OF AXIALLY LOADED FUNCTIONALLY GRADED SANDWICH BEAMS USING REFINED SHEAR DEFORMATION THEORY Thuc P. Vo 1, Adelaja Israel Osofero 1, Marco Corradi 1, Fawad Inam 1 1 Faculty of Engineering
More informationANALYSIS OF STRAINS CONCEPT OF STRAIN
ANALYSIS OF STRAINS CONCEPT OF STRAIN Concept of strain : if a bar is subjected to a direct load, and hence a stress the bar will change in length. If the bar has an original length L and changes by an
More informationMicro-meso draping modelling of non-crimp fabrics
Micro-meso draping modelling of non-crimp fabrics Oleksandr Vorobiov 1, Dr. Th. Bischoff 1, Dr. A. Tulke 1 1 FTA Forschungsgesellschaft für Textiltechnik mbh 1 Introduction Non-crimp fabrics (NCFs) are
More informationA Numerical Study of Chamber Size and Boundary Eects on CPT Tip Resistance in NC Sand
Scientia Iranica, Vol. 15, No. 5, pp 541{553 c Sharif University of Technology, October 2008 A Numerical Study of Chamber Size and Boundary Eects on CPT Tip Resistance in NC Sand M.M. Ahmadi 1; and P.K.
More informationAnalysis of Axially Loaded Non-prismatic Beams with General End Restraints Using Differential Quadrature Method
ISBN 978-93-84422-56-1 Proceedings of International Conference on Architecture, Structure and Civil Engineering (ICASCE'15 Antalya (Turkey Sept. 7-8, 2015 pp. 1-7 Analysis of Axially Loaded Non-prismatic
More informationStress Analysis Lecture 3 ME 276 Spring Dr./ Ahmed Mohamed Nagib Elmekawy
Stress Analysis Lecture 3 ME 276 Spring 2017-2018 Dr./ Ahmed Mohamed Nagib Elmekawy Axial Stress 2 Beam under the action of two tensile forces 3 Beam under the action of two tensile forces 4 Shear Stress
More informationAvailable online at ScienceDirect. Procedia IUTAM 13 (2015 ) 82 89
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 13 (215 ) 82 89 IUTAM Symposium on Dynamical Analysis of Multibody Systems with Design Uncertainties The importance of imperfections
More informationChapter 2 Governing Equations
Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is
More information1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load
1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationScientia Iranica, Vol. 13, No. 2, pp 113{123 c Sharif University of Technology, April 2006 Buckling of Discretely Stringer-Stiened Composite Cylindric
Scientia Iranica, Vol. 13, No., pp 113{13 c Sharif University of Technology, April 006 Buckling of Discretely Stringer-Stiened Composite Cylindrical Shells under Combined Axial Compression and External
More informationCO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS
COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL
More informationHomework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. Fall 2004
Homework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. 1. A beam is loaded as shown. The dimensions of the cross section appear in the insert. the figure. Draw a complete free body diagram showing an equivalent
More information(Refer Slide Time: 1: 19)
Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 4 Lecture - 46 Force Measurement So this will be lecture
More informationFree vibration and stability of axially functionally graded tapered Euler-Bernoulli beams
Shock and Vibration 18 211 683 696 683 DOI 1.3233/SAV-21-589 IOS Press Free vibration and stability of axially functionally graded tapered Euler-Bernoulli beams Ahmad Shahba a,b, Reza Attarnejad a,b, and
More informationMarch 24, Chapter 4. Deflection and Stiffness. Dr. Mohammad Suliman Abuhaiba, PE
Chapter 4 Deflection and Stiffness 1 2 Chapter Outline Spring Rates Tension, Compression, and Torsion Deflection Due to Bending Beam Deflection Methods Beam Deflections by Superposition Strain Energy Castigliano
More informationFinal EE290G Intro to MEMS 12/16/98
Final EE290G Intro to MEMS 12/16/98 Name SID 1. (20 points) Using poly1 in the MCNC/MUMPS process: Fx Fy L b L a x y cxx = c yx c xy c yy Fx F y (1) (a) Calculate c xy, the constant relating the force
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationDynamics of Nonlinear rectangular plates subjected to an orbiting mass based on shear deformation plate theory
JCAMECH Vol. XX, No. X, XXXX 2XX, pp XXX-XXX DOI: 1.2259/jcamech.217.238716.169 Dynamics of Nonlinear rectangular plates subjected to an orbiting mass based on shear deformation plate theory Reza Javidi,
More informationUnit 15 Shearing and Torsion (and Bending) of Shell Beams
Unit 15 Shearing and Torsion (and Bending) of Shell Beams Readings: Rivello Ch. 9, section 8.7 (again), section 7.6 T & G 126, 127 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics and Engineering
More information[5] Stress and Strain
[5] Stress and Strain Page 1 of 34 [5] Stress and Strain [5.1] Internal Stress of Solids [5.2] Design of Simple Connections (will not be covered in class) [5.3] Deformation and Strain [5.4] Hooke s Law
More informationNonlinear effects on the rotor driven by a motor with limited power
Applied and Computational Mechanics 1 (007) 603-61 Nonlinear effects on the rotor driven by a motor with limited power L. Pst Institute of Thermomechanics, Academy of Sciences of CR, Dolejškova 5,18 00
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationStatic and free vibration analysis of carbon nano wires based on Timoshenko beam theory using differential quadrature method
8(2011) 463 472 Static and free vibration analysis of carbon nano wires based on Timoshenko beam theory using differential quadrature method Abstract Static and free vibration analysis of carbon nano wires
More informationGeneral Relativity ASTR 2110 Sarazin. Einstein s Equation
General Relativity ASTR 2110 Sarazin Einstein s Equation Curvature of Spacetime 1. Principle of Equvalence: gravity acceleration locally 2. Acceleration curved path in spacetime In gravitational field,
More informationSOLUTION (17.3) Known: A simply supported steel shaft is connected to an electric motor with a flexible coupling.
SOLUTION (17.3) Known: A simply supported steel shaft is connected to an electric motor with a flexible coupling. Find: Determine the value of the critical speed of rotation for the shaft. Schematic and
More informationExpedient Modeling of Ball Screw Feed Drives
S. Frey a A. Dadalau a A. Verl a Expedient Modeling of Ball Screw Feed Drives Stuttgart, February 2011 a Institute for Control Engineering of Machine Tools and Manufacturing Units (ISW), University of
More informationForced Vibration Analysis of Timoshenko Beam with Discontinuities by Means of Distributions Jiri Sobotka
21 st International Conference ENGINEERING MECHANICS 2015 Svratka, Czech Republic, May 11 14, 2015 Full Text Paper #018, pp. 45 51 Forced Vibration Analysis of Timoshenko Beam with Discontinuities by Means
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More informationLoading σ Stress. Strain
hapter 2 Material Non-linearity In this chapter an overview of material non-linearity with regard to solid mechanics is presented. Initially, a general description of the constitutive relationships associated
More informationFree vibration analysis of beams by using a third-order shear deformation theory
Sādhanā Vol. 32, Part 3, June 2007, pp. 167 179. Printed in India Free vibration analysis of beams by using a third-order shear deformation theory MESUT ŞİMŞEK and TURGUT KOCTÜRK Department of Civil Engineering,
More informationHomework 6 Solutions. Solution. Note {e t, te t, t 2 e t, e 2t } is linearly independent. If β = {e t, te t, t 2 e t, e 2t }, then
Homework 6 Solutions 1 Let V be the real vector space spanned by the functions e t, te t, t 2 e t, e 2t Find a Jordan canonical basis and a Jordan canonical form of T on V dened by T (f) = f Solution Note
More information