Vibrations of beams with a variable cross-section fixed on rotational rigid disks

Size: px
Start display at page:

Download "Vibrations of beams with a variable cross-section fixed on rotational rigid disks"

Transcription

1 1(13) Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks Abstract The work is focused on the probem of vibrating beams with a variabe cross-section fixed on a rotationa rigid disk. The beam is oaded by a transversa time varying force orthogona to an axis of the beam and simutaneousy parae to the disk s pane. There are many ways of usage of the technica moveabe systems composed of eements with the variabe cross-sections. The main appications are used in numerous types of turbines and pumps. The paper is a kind of introduction to the dynamic anaysis of above mentioned beam systems. The equations of motion of rotationa beams fixed on the rigid disks were derived. After introducing the Coriois forces and the centrifuga forces, the transportation effect in the mathematica mode was considered. This particuar project is the first stage research, where there were proposed certain soutions of probems connected with the inear variabe cross-sections systems. The further investigation considering the noninear systems has been proceeding. The resuts, anaysis and comparison wi be presented in the future works. Keywords variabe cross-section, beam, rotary motion, dynamic fexibiity Sawomir Zokiewski * Siesian University of Technoogy, Institute of Engineering Processes Automation and Integrated Manufacturing Systems, Division of Mechatronics and Designing of Technica Systems, Facuty of Mechanica Engineering, 18a Konarskiego Street, 44-1 Giwice, Poand * Author emai: sawomir.zokiewski@pos.p 1 INTRODUCTION The fundamenta eements of the numerous mechanisms or machines, namey the beams and rods, are most commony considered as the vibrating stationary systems. In these cases they are treated as the immovabe ones without the couping between the oca osciations and the main operationa motion. The mechanica systems are predominanty considered both in the kinematica and dynamica aspects. The main probems resuting from the anaysis are connected with the issues of controing and stabiizing mechanica systems as such. Numerous research can be found in the iterature [1, 3, 4, 7, 9, 14-16, 18-] and so far the soutions have been found by considering the transportation (considered as the main working motion) and the oca osciations separatey [, 3, 5, 6, 8, 1-1, 17]. Mesut Simsek in his paper [14] presents the noninear dynamic anaysis of a functionay graded beam with pinned-pinned supports. Simsek uses the Timoshenko beam theory with the von Karman's noninear strain-dispacement reationships. The beam was anayzed in reation to its Latin American Journa of Soids and Structures 1(13) 39 57

2 4 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks variabe thickness of a cross-section according to the power-aw formua. The system of equations of motion is derived by using the Lagrange s equations. Transverse and axia defections and rotation of the cross-sections of the beam were expressed in poynomia forms. The boundary conditions are taken into account by using the Lagrange mutipiers. The Newmark-b method in conjunction with the direct iteration method was used to sove noninear equations of motion. In the work, arge defection, veocity of the moving oad and excitation frequency on the beam dispacements were anayzed. Some detaied cases were compared with each other. Another method of anaysis of highspeed rotating beams is presented in [9] by Gunda. The author introduces a new finite eement for free vibration anaysis of rotating beams. The basis shape functions which appy a inear combination of the soution of the governing static differentia equation of a stiff-string and a cubic poynomia are used. The introduced functions depend on the anguar veocity and position of the anayzed section, what accounts for the centrifuga stiffening effect. In paper [1] the effects of rotary inertia were presented. The author shows the rotation impact on the extensiona tensie force and the eigenvaues of beams. The anayzed beam was rotated uniformy about a transverse axis taking into account ongitudina easticity. The perturbation technique and Gaerkin's method were used. In the exampe of a typica heicopter rotor bade it was indicated that the extensiona tensie force increases up to ten percent, when the rotary inertia contribution is retained in the modeing. Szefer in [15, 16] caimed that the vibrations from fexibiity of eements of the mechanica composition are much smaer than a main disocation of this composition. There is an increased range of veocities and acceerations of such a type systems. A those are caused by using more efficient drives or/and ess and ess weighting materias. In order to constraint power output of drives for motion of mechanica systems, materias with the ower mass density and stiffness are used. A those things are reasons for creating the new modes of designing technica beamike systems. The variabe crosssection beams are appied in many practica appications nowadays. Both the ineary and nonineary variabe cross-sections are used in technica systems. These systems shoud be taken under consideration with respect to a reaistic distribution of mass and fexibiity. Any variations in crosssections cannot be negected and shoud be considered in the mathematica mode assumptions. In spite of many works [e.g. 1, 9-1, 17, ], where their authors anayzes the beams in rotations and beams with variabe cross-sections, this work is devoted to genera anayses both the inear and noninear ones of vibrating transversay beams with the variabe cross-section. It is beieved that new equations and resuts are presented in this paper. MODEL OF THE ANALYZED BEAMS FIXED ON THE ROTATIONAL RIGID DISK In this section the way of modeing of beams in rotationa transportation is presented. The beam was assumed as the uniform one with the variabe cross-section. The beam was fixed on the rotationa disk that was assumed as a rigid eement (Fig. 1). The disk rotates round its axis and the motion was described in two reference frames, the goba stationary reference frame and the oca one. The vibrations occurring in the oca reference frame are transferred into the goba reference frame with taking into consideration so caed the transportation effect [3, 16, 19, ].

3 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 41 Fig. 1 Mode of the anayzed beam on the rotationa disk [3] Fig. Mode of the anayzed beam on the rotationa disk the XZ projection [3] where: ρ mass-density, A(x) cross-section, ength of the beam, x ocation of the anayzed cross-section, ω anguar veocity, Ω frequency, Q rotation matrix, S position vector, F harmonic force, E Young moduus, w vector of dispacement. There is defined a set of generaized coordinates aowing to make an assumption that generaized coordinates are the orthogona projections on the coordinate axes of the goba reference frame and can be written thus: q = r, 1 X (1) q = r, Y () and after differentiating (Eq. 1) and (Eq. ) can be obtained the generaized veocities in the form: q 1 = dq 1 dt = r X = v X, (3) q = dq dt = r Y = v Y, (4) Latin American Journa of Soids and Structures 1(13) 39 57

4 4 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks or in the equivaent vector form expressing the dispacements: r = [r X,r Y ] T (5) and the equivaent vector expressing the veocities: r = [ r X, r Y ] T. (6) The cross-section is dependent on the x coordinate and for instance making an assumption that the anayzed beam is a cone (a base with radius equa x/) it can be written as foows: = A 1 x A x = π 1 x, (7) The anguar veocity of the rotationa disk treated as the transportation veocity is defined as: ω = ω T, (8) The ocation of the anayzed cross-section is described by the position vector defined as foows: S = s T. (9) The xyz reference frame is the rotationa one and the matrix aowing expressing vectors is the rotation matrix. The xyz oca reference frame is fixed to the rotating beam and its z-axis coincides with the rotation axis of the rigid disk, thus the rotation matrix is defined as such: cosϕ sinϕ Q = sinϕ cosϕ. (1) 1 Different appications of such a type beam systems may have different boundary conditions and the cross-section can be considered as the ineary or nonineary variabe one. For exampe the beam can be considered as the fixed at one end, the camped-camped one or with other boundary conditions etc. For exampe beams with one end camped and the second free-free are considered (Fig. 1). The beam is fixed on the rotationa rigid disk. An anayzed section of the beam is oaded by a harmonic time varying distributed force with unitary ampitude on a direction perpendicuar to the centre ine of the beam. The forces on origins of the beams are assumed as equa zero and aso the dispacements are equa zero, because the system is camped in one end. Because of a these the specific boundary conditions for this beam shoud be written as foows:

5 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 43 =, w,t w (,t ) =, w (,t ) =, EI x ( w,t ) = F δ( x )cos( Ωt )dx, (11) in every time moment t. For the beams with the ineary variabe cross-section the geometric moment of inertia can be assumed as the function proportiona to the constant initia moment I for exampe as: 4 I ( x ) = I 1 x. (1) After soving the boundary vaue probem, the eigenfunction for the dispacement can be derived in the form: = x n C 1 J n k x X x C Y n ( k x ) C 3 I n ( k x ) C 4 K n ( k x ), (13) where [13]: J n the Besse function of the first kind, Y n the Besse function of the second kind, I n the hyperboic Besse function of the first kind, K n the hyperboic Besse function of the second kind, C 1,,3,4 the integrations constants, where: n is a mode of vibrations of the beam with a variabe cross-section. The k coefficient is a inear approximation of the reaistic eigenvaues of the anayzed systems. 3 EQUATIONS OF MOTION OF THE ANALYZED SYSTEMS The equations of motion of the anayzed beamike systems are obtained with the hep of the cassica method. The derived equations of motion can be presented in the matrix form as foows: Latin American Journa of Soids and Structures 1(13) 39 57

6 44 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks ρa x ρa x cosϕ sin ϕ EI x sin ϕ cosϕ 1 cosϕ sin ϕ sin ϕ cosϕ 1 cosϕ sin ϕ sin ϕ cosϕ 1 w ωs ρa x = EI x 3 w 3 cosϕ sin ϕ sin ϕ cosϕ EI x 1 cosϕ sin ϕ sin ϕ cosϕ 1 ω s cosϕ sin ϕ sin ϕ cosϕ 1 4 w 4 w (14) in each point of the range D = {( x,t ),x (, ),t > } the Eq. (3.1) coincides with boundary conditions and initias conditions. The projections of the equations of motion in the goba reference frame are obtained from the Eq. (14) and in the XYZ goba reference frame equations for individua axes can be expressed as: ω s cosϕ ωs sinϕ =, (15) ρa( x ) w EI ( x ) 4 w ω s sinϕ ωs cosϕ =, (16) EI x 4 3 w EI x 3 w =. (17) Third governing equation (17) is a we-known one and represents the equation of motion of the immovabe beam with the variabe cross-section. As we can see the boundary conditions and the oading direction does not provide any new aspects to the dynamica anaysis and the beam (Fig. ) can be treated as equivaent to the stationary one. Much more interesting case of the oading is the beam oaded by a transversa time varying force orthogona to an axis of the beam and simutaneousy parae to the disk s pane. This case (Fig. 3) where the beam with a variabe cross-section A(x) and a variabe moment of inertia I(x) is considered.

7 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 45 Fig. 3 Mode of the anayzed beam with the ineary variabe cross-section fixed on the rotationa disk and oaded by a transversa force orthogona to an axis of the beam and parae to the disk s pane [3] Anaogicay to the previous equations of motion it can be derived by means of the cassica method. Interesting case of the oading is the beam oaded by a transversa time varying force orthogona to an axis of the beam and simutaneousy parae to the disk s pane. The obtained equations of motion are couped ones with the Coriois eements. After oading the beam in the direction of the y axis of the oca reference frame it can be derived the equations of motion of the beam in the matrix form. This system of equations of motion is the fourth order partia differentia equations (PDE) and can be expressed as foows: ρa x ρa x cosϕ sinϕ sinϕ cosϕ 1 cosϕ sinϕ sinϕ cosϕ 1 cosϕ sinϕ = sinϕ cosϕ 1 w ρa x ω w EI x w ρa x in each point of the range D {( x, t), x (, ), t } cosϕ sinϕ sinϕ cosϕ 1 cosϕ sinϕ sinϕ cosϕ 1 ω s ω w ωw ωs = (18) = Œ > the Eq. (18) coincides with boundary conditions and initias conditions. The projection onto the X and Y axes, can be obtained, where the projection onto the X axis of the goba reference frame is as foows: EI x ρa( x ) w X ρa( x ) ωw Y ωs sinϕ w X = ρa x ρa ( x )ω w Y ω ( w X s cosϕ ), (19) Latin American Journa of Soids and Structures 1(13) 39 57

8 46 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks and for the Y axis of the goba reference frame, the projection can be written: ρa( x ) w Y EI x ρa( x ) ωw X ωs cosϕ w Y = ρa x ρa ( x )ω w X ω ( w Y s sinϕ ), () where using the rea coordinates the individua dispacements onto the axes of the goba reference frame are: w X = w sinϕ, w Y = w cosϕ. (1) It is obvious that the vibrating beams can be described by different boundary conditions for different oading and fixations types and this case can be treated as the genera one. 4 MATHEMATICAL MODEL OF THE SYSTEM The mathematica mode of the anayzed beam is provided by the orthogonaization of the equations of motion (19) and (). The orthogonaization is provided by the Eqs. (19) and () which are mutipied using the eigenfunction for the dispacements and computing integras from the equations of motion in the beams imits of integration from the origin of the beam (zero) to the end (a ength of the beam). After assuming a constant anguar veocity, an anguar acceeration equas zero and after introducing Eqs. (1), then the projecting equations of motion Eq. (18) onto the axes of the goba reference frame can be written: ρa( x ) w X X ( x )dx EI ( x ) 4 w X X ( x )dx 4 EI x 3 w X 3 X ( x )dx EI x w X X ( x )dx ρ A( x )ω w Y X ( x )dx ρa( x )ω w X X ( x )dx =, () and for the second equation, aso mutipied by the eigenfunction for dispacement and integrated aong the beam: ρa( x ) w Y X ( x )dx EI ( x ) 4 w Y X ( x )dx 4 EI x 3 w Y 3 X ( x )dx EI x w Y X ( x )dx ρ A( x )ω w X X ( x )dx ρa( x )ω w Y X ( x )dx =. (3)

9 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 47 The presented mathematica mode is universa for different functions describing the variabiity of the beam s cross-section. The cross-section of the anayzed beam can be now treated as both ineary and nonineary variabe. If the beam has the nonineary variabe cross-section (Fig. 4-5) the appropriate function describing this variation ought to be adopted. The different oading terms can be considered as we. For exampe the beam is oaded by a transversa force orthogona to an axis of the beam and orthogona to the disk s pane (Fig. 4) or the beam is oaded by a transversa force orthogona to an axis of the beam and in parae to the disk s pane (Fig. 5). Fig. 4 Mode of the anayzed beam with the nonineary variabe cross-section fixed on the rotationa disk and oaded by a transversa force orthogona to an axis of the beam and orthogona to the disk s pane These cases can be considered separatey or combined. According to the scientific iterature the case presented in figure 4 is we-known in the mathematica sense, ikewise the case in figure 1. Fig. 5 Mode of the anayzed beam with the nonineary variabe cross-section fixed on the rotationa disk and oaded by a transversa force orthogona to an axis of the beam and in parae to the pane of the disk It is assumed that eements of the Eqs. () and (3) are the continuousy differentiabe functions. Using the Eq. () and mutipe integrating both sides of this equation by parts, the foowing formua can be obtained: Latin American Journa of Soids and Structures 1(13) 39 57

10 48 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks E X ( x )I x 3 w X 3 X ( x )I x w X X ( x )I x E 3 X ( x )I ( x ) w 3 X E I ( x )X IV ( x )w X dx 4E I ( x ) X ( x )w X dx 6E I ( x ) X ( x )w X dx 4E I ( x ) X ( x )w X dx E I IV ( x )X ( x )w X dx ρa( x ) w X X ( x )dx E X x I ( x ) w X X x I ( x ) E I ( x ) X ( x )w X dx 6E I x w X X X x w X I ( x ) ( x )w X dx 6E I ( x ) X ( x )w X dx E I IV ( x )X ( x )w X dx E X x E I x I ( x ) w X X X x I ( x ) w X E I x ( x )w X dx E I IV ( x )X ( x )w X dx X ρ A( x )ω w Y X ( x )dx ρa( x )ω w X X ( x )dx =, w X ( x )w X dx (4) where the variabe cross-section A(x) and the moment of inertia I(x) depend on the geometry of the anayzed beam. There is a beam with the ineary variabe cross-section Eqs. (7) and (1) assumed for further cacuations. The eigenfunction associated with transversa vibrations X(x) of the beam is obtained by appying undamped conditions. Simiary to the first equation, the integration resuts for the second equation (3) can be appied as foows:

11 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 49 E X ( x )I x 3 w Y 3 X ( x )I x w Y X ( x )I x E 3 X ( x )I ( x ) w 3 Y E I ( x )X IV ( x )w Y dx 4E I ( x ) X ( x )w Y dx 6E I ( x ) X ( x )w Y dx 4E I ( x ) X ( x )w Y dx E I IV ( x )X ( x )w Y dx ρa( x ) w Y X ( x )dx E X x I ( x ) w Y X x I ( x ) w Y X x w Y I ( x ) E I ( x ) X ( x )w Y dx 6E I ( x ) X ( x )w Y dx 6E I ( x ) X ( x )w Y dx E I IV ( x )X ( x )w Y dx E X x E I x I ( x ) w Y X X x I ( x ) w Y E I x ( x )w Y dx E I IV ( x )X ( x )w Y dx X ρ A( x )ω w X X ( x )dx ρa( x )ω w Y X ( x )dx =. w Y ( x )w Y dx (5) It is possibe to simpify the obtained equations as shown in the Eq. (6). Additionay, after organizing the eements in order, the Eqs. (4) and (5) can be expressed as: Latin American Journa of Soids and Structures 1(13) 39 57

12 5 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks E X ( x )I x 3 w X 3 X ( x )I x w X X ( x )I x w X E 3 X ( x )I ( x ) w 3 X E I ( x )X IV ( x )w X dx E I ( x ) X ( x )w X dx ρa( x ) w X X ( x )dx E I x E X x E X x I ( x ) w X I ( x ) w X X x X x X I ( x ) I ( x ) ( x )w X dx w X w X X x I ( x ) ρ A( x )ω w Y X ( x )dx ρa( x )ω w X X ( x )dx =, w X (6) and the couped second equation: E X ( x )I x 3 w Y 3 X ( x )I x w Y X ( x )I x w Y E 3 X ( x )I ( x ) w 3 Y E I ( x )X IV ( x )w Y dx E I ( x ) X ( x )w Y dx E I ( x ) X ( x )w Y dx ρa( x ) w Y X ( x )dx E X x E X x I ( x ) w Y I ( x ) w Y X x I ( x ) I X x ρa( x )ω w Y X ( x )dx =. ( x ) w Y w Y X x I ( x ) ρ A( x )ω w X w Y X ( x )dx (7) There are introduced the equivaent notation for the eigenfunctions, cross-sections and the moments of inertia as: X ( x ) = X, A( x ) = A, I ( x ) = I. (8)

13 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 51 The Eqs. (6) and (7) after substituting (8) are simpified to: E XI 3 w X X I X I 3 w X w X X I X I X I E ( X I 3 X I 3 X I I X )w X E IX IV w Y dx E I X w dx E I X w dx ρa w X X X X dx E X I w X X I X I w X E X I w X ( X I X I )w X ( X I X I I X )w X ρ Aω w Y X dx ρaω w X X dx = (9) and anaogicay to the second equation (7) for the Y axis of the goba reference frame can be obtained as: E XI 3 w Y X I X I 3 w Y w Y X I X I X I E ( X I 3 X I 3 X I I X )w Y E IX IV w Y dx E I X w dx E I X w dx ρa w Y Y Y X dx E X I w Y X I X I w Y E X I w Y ( X I X I )w Y ( X I X I I X )w Y ρ Aω w X X dx ρaω w Y X dx =. (3) The Eqs. (9) and (3) can be simpified as we and the foowing two equations yied: Latin American Journa of Soids and Structures 1(13) 39 57

14 5 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks E XI 3 w X X I X I 3 w X X I w X E ( X I X I )w X E IX IV w X dx E I X w dx E I X w dx ρa w X X X X dx ρ Aω w Y X dx ρaω w X X dx =, (31) the second equation (3) additionay after simpifying the integras: E XI 3 w Y X I X I 3 w Y X I w Y E ( X I X I )w Y E I X w Y dx ρa w Y X dx ρ Aω w X X dx ρaω w Y X dx =. (3) By introducing the acting externa force, it can be noticed: ( ) ( ) EI 3 w X E I w X =, 3 EI 3 w Y E I w Y = F 3 δ( x )cos( Ωt )dx. (33) Based on the boundary conditions (11) it can be written: ( ) w X = w Y =, XI 3 w X 3 XI 3 w Y 3 XI XI w X =, w Y = (34)

15 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 53 and E X I w X X I E X I w Y X I w X w Y =, =. (35) By introducing (34-35) the first equation yieds: E IX IV w X dx E I X w dx E I X w dx ρa w X X X X dx ρ Aω w Y X dx ρaω w X X dx =. (36) Where the searched soutions are defined as dispacements in the goba reference frame and the ones are separabe in time and space domains. The system is symmetric and the motion of the anayzed cross-section is given by the mutipication of the projected harmonic time dependent eigenfunctions with the frequency Ω and the eigenfunctions for dispacements with ampitudes XY. Both the variabes coincide with the responses of two decouped harmonic osciators. In the present case the searched soutions are: w X = n=1 a X X ( x )sin( Ωt ), w Y = a Y X ( x )cos( Ωt ). n=1 (37) The norm for Eqs. (36) is assumed as: γ n = AX dx, (38) for the beam with the ineary variabe cross-section (Fig. 3) the norm is cacuated as integra. Both I 1 and I coefficients are connected with the moments of inertia and the eigenfunctions for dispacements. The I and c (veocity) coefficients are as foows: Latin American Journa of Soids and Structures 1(13) 39 57

16 54 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks I 1 = 1 IX dx, γ n I = 1 γ I X dx, n (39) c = E ρ. It can be aso noticed that: X I dx =. (4) The a coefficients (ampitudes) can be obtained as: a X = ργ n c 4 k 4 I I 1 k ωωf X ( ) ( ω Ω ) c k ( I I 1 k )( ω Ω ), (41) a Y = F X ργ n 1 I k c I 1 k 4 ω Ω I k c I 1 k 4 ω Ω 1. (4) Knowing these ampitudes the tota dispacements can be cacuated and the soution in the goba reference frame can be obtained. 5 DYNAMICAL FLEXIBILITIES According to the author the dynamica fexibiity is understood as the ampitude of generaized dispacement in the direction of i generaized coordinate changed by generaized harmonic force with ampitude equas 1 in the direction of k generaized coordinate. This definition can be written in the mathematica form as: 1 s i = Y ik s k. (43) Finay the dynamica fexibiity of the rotating beam with the variabe cross-section is presented as: 1 X ( x )X ( ) Y = ργ n I k c I 1 k 4 ( ω Ω ) X ( x )X ( ). I k c I 1 k 4 ( ω Ω ) (44) n=1

17 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 55 For comparison, the dynamica fexibiity of the rotating beam with a constant cross-section is as foows: X ( x )X ( ) c g k 4 ω Ω Y g =. ρa const γ gn c 4 g k 8 ( ω Ω ) c g k 4 ( ω Ω n=1 ) On the other side the dynamica fexibiity of the immovabe beam is presented as: (45) X ( x )X ( ) Y g =. n=1 ρa const γ gn ( c g k 4 Ω (46) ) where for the Eqs. (44) and (45): and the norm is: c g = γ gn EI ρa const (47) = X dx. (48) Both the derived dynamica fexibiities, the Eq. (44) and Eq. (46), are presented in figures (Fig. 6-7). For numerica exampes the beam of 1m ength is assumed. The beam is the camped-free one and fixed on the rotationa disk. Distinctivey from the beam with a constant cross-section, the beam with the variabe cross-section has different inconstant bifurcations for the individua modes of vibrations (Fig. 7). Fig. 6 Sampe dynamica fexibiity of the anayzed immovabe beam with the variabe cross-section Latin American Journa of Soids and Structures 1(13) 39 57

18 56 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks Fig. 7 Dynamica fexibiity of the rotating beam with the variabe cross-section. The disk is being rotated with the anguar veocity equas 1 radians per second Derivations of suitabe dynamica fexibiities are appied to the systems in the transportabe motion both vibrating ongitudinay and transversay. The soutions can be presented on charts of ampitude (dynamica fexibiity) in function of frequency. There are many numerica appications that can be used for this purpose, e.g. GENTA S DynRot or author s program the Modyfit [18]. Concusions The paper concerns a very crucia probem of beams vibrating in transportation with the variabe geometric parameters such as: tapered or coned beams, bades with various shape functions, and others. The anayzed systems coud be appied in many technica impementations, obviousy after proper adjustment of mathematica modes and considering the necessary parameters and conditions. The most popuar appications coud be connected with pumps and rotors; especiay the high speed ones, bades of heicopters and arms of robots. The obtained soutions were aso presented in the form of equations of motion. The dynamica anaysis was provided using the dynamica fexibiity method and the Gaerkin's method. The dynamica fexibiity method appies soutions presentation as the dynamica characteristics of ampitude in the function of frequency. The soutions were determined as the compact mathematica formuae. The dynamica characteristics were presented on charts. The transportation effect and shape functions describing the cross-section variation were taken into consideration. The centrifuga forces and the Coriois forces were taken into consideration in the obtained mathematica mode. The natura frequencies and the minima ampitudes can be easiy read from the dynamica fexibiity characteristics. In these characteristics the beam with the constant cross-section and the beam with the variabe cross-section are put together. The juxtaposition of the rotating beam and the immovabe one is presented. The resuts presented are very significant and enabe controing of such a type rotating systems. In future works more numerica exam-

19 Sawomir Zokiewski / Vibrations of beams with a variabe cross-section fixed on rotationa rigid disks 57 pes for beamike systems with the variabe cross-section wi be presented and the rotating beam with damping forces considered wi be anayzed. The rotating beam with some imperfections and defects or cracks can be the potentia materia for the further anaysis and research. Acknowedgments: the paper was presented during 11th Conference on Dynamica Systems - Theory and Appications. References A-Ansary M. D.: Fexura vibrations of rotating beams considering rotary inertia. Computers and Structures, 69, 1998, Aucieo N. M.: Comment on. A note on vibrating tapered beams. Journa of Sound and Vibration, 1995, Bokhonsky A. I., Zokiewski S.: Modeing and Anaysis of Eastic Systems in Motion. Monograph, Giwice, Siesian University of Technoogy Press, 11, p Chang T.-P., Chang H.-C.: Vibration and bucking anaysis of rectanguar pates with noninear eastic end restraints against rotation. Internationa Journa of Soids and Structures, 34, 18, June 1997, Cheung Y. K., Zhou D.: The free vibrations of tapered rectanguar pates using a new set of beam functions with the Rayeigh Ritz method. Journa of Sound and Vibration, 3, 5, 1999, Dyniewicz B., Bajer C. I.: New feature of the soution of a Timoshenko beam carrying the moving mass partice. Archives of Mechanics, 6, 5, 1, Genta G.: Dynamics of Rotating Systems. New York, Springer, 5. Grossi RO, Bhat RB. A note on vibrating tapered beams. Journa of Sound and Vibration 1991, Gunda J. B., Gupta R. K., Gangu R.: Hybrid stiff-string poynomia basis functions for vibration anaysis of high speed rotating beams. Computers and Structures, 87, 9, Kumar R., Kansa T.: Effect of rotation on Rayeigh waves in an isotropic generaized thermoeastic diffusive hafspace. Archives of Mechanics, 6, 5, 8, Ozturk H.: In-pane free vibration of a pre-stressed curved beam obtained from a arge defected cantiever beam. Finite Eements in Anaysis and Design, 47, 11, Rahman M. A., Kowser M. A., Hossain S. M. M.: Large defection of the cantiever stee beams of uniform strength - experiment and noninear anaysis. Internationa Journa of Theoretica and Appied Mechanics, 1, 6, Soecki R., Szymkiewicz J.: The rod and surface systems. Warszawa, Arkady Budownictwo, Sztuka, Architektura, 1964 (in Poish). Simsek M.: Non-inear vibration anaysis of a functionay graded Timoshenko beam under action of a moving harmonic oad. Composite Structures, 9, 1, Szefer G.: Dynamics of eastic bodies undergoing arge motions and uniatera contact. Journa of Technica Physics. XLI, 4,. Szefer G.: Dynamics of eastic bodies in terms of pane frictiona motion. Journa of Theoretica and Appied Mechanics,, 39, 1. Zhou D.: Vibrations of point-supported rectanguar pates with variabe thickness using a set of static tapered beam functions. Internationa Journa of Mechanica Sciences, 44,, Zokiewski S.: Numerica Appication for Dynamica Anaysis of Rod and Beam Systems in Transportation. Soid State Phenomena, 164, 1, Zokiewski S.: Dynamic fexibiity of the supported-camped beam in transportation. Journa of Vibroengineering. 13, 4, December 11, Zokiewski S.: Dynamica Fexibiity of Compex Damped Systems Vibrating Transversay in Transportation. Soid State Phenomena, 164, Zokiewski S.: Damped Vibrations Probem Of Beams Fixed On The Rotationa Disk. Internationa Journa of Bifurcation and Chaos, 1, 1, October 11, Zone-Ching L., Don-Tsun L.: Dynamic defection anaysis of a panar robot Computers and Structures, 53, 4, 17 November 1994, Latin American Journa of Soids and Structures 1(13) 39 57

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

VTU-NPTEL-NMEICT Project

VTU-NPTEL-NMEICT Project MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid

More information

Modal analysis of a multi-blade system undergoing rotational motion

Modal analysis of a multi-blade system undergoing rotational motion Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa

More information

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18 Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

2.1. Cantilever The Hooke's law

2.1. Cantilever The Hooke's law .1. Cantiever.1.1 The Hooke's aw The cantiever is the most common sensor of the force interaction in atomic force microscopy. The atomic force microscope acquires any information about a surface because

More information

LECTURE 10. The world of pendula

LECTURE 10. The world of pendula LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive

More information

On a geometrical approach in contact mechanics

On a geometrical approach in contact mechanics Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

Identification of macro and micro parameters in solidification model

Identification of macro and micro parameters in solidification model BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University

More information

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE 3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses

More information

Dynamic Stability of an Axially Moving Sandwich Composite Web

Dynamic Stability of an Axially Moving Sandwich Composite Web Mechanics and Mechanica Engineering Vo. 7 No. 1 (2004) 53-68 c Technica University of Lodz Dynamic Stabiity of an Axiay Moving Sandwich Composite Web Krzysztof MARYNOWSKI Department of Dynamics of Machines

More information

830. Nonlinear dynamic characteristics of SMA simply supported beam in axial stochastic excitation

830. Nonlinear dynamic characteristics of SMA simply supported beam in axial stochastic excitation 8. Noninear dynamic characteristics of SMA simpy supported beam in axia stochastic excitation Zhi-Wen Zhu 1, Wen-Ya Xie, Jia Xu 1 Schoo of Mechanica Engineering, Tianjin University, 9 Weijin Road, Tianjin

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

Keywords: Functionally Graded Materials, Conical shell, Rayleigh-Ritz Method, Energy Functional, Vibration.

Keywords: Functionally Graded Materials, Conical shell, Rayleigh-Ritz Method, Energy Functional, Vibration. Journa of American Science, ;8(3) Comparison of wo Kinds of Functionay Graded Conica Shes with Various Gradient Index for Vibration Anaysis Amirhossein Nezhadi *, Rosan Abdu Rahman, Amran Ayob Facuty of

More information

MA 201: Partial Differential Equations Lecture - 10

MA 201: Partial Differential Equations Lecture - 10 MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary

More information

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2 MATEC Web of Conferences 159, 01005 (018) IJCAET & ISAMPE 017 https://doi.org/10.1051/matecconf/01815901005 Dynamic Response of RC Cantiever Beam by Equivaent Singe Degree of Freedom Method on Eastic Anaysis

More information

MECHANICAL ENGINEERING

MECHANICAL ENGINEERING 1 SSC-JE SFF SELECION COMMISSION MECHNICL ENGINEERING SUDY MERIL Cassroom Posta Correspondence est-series16 Rights Reserved www.sscje.com C O N E N 1. SIMPLE SRESSES ND SRINS 3-3. PRINCIPL SRESS ND SRIN

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

3.10 Implications of Redundancy

3.10 Implications of Redundancy 118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.

More information

Technical Data for Profiles. Groove position, external dimensions and modular dimensions

Technical Data for Profiles. Groove position, external dimensions and modular dimensions Technica Data for Profies Extruded Profie Symbo A Mg Si 0.5 F 25 Materia number.206.72 Status: artificiay aged Mechanica vaues (appy ony in pressing direction) Tensie strength Rm min. 245 N/mm 2 Yied point

More information

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes... : Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................

More information

SEA Subsystem Attribute Correction Based on Stiffness Multipliers

SEA Subsystem Attribute Correction Based on Stiffness Multipliers Internationa Conference on Eectrica, Mechanica and Industria Engineering (ICEMIE 016) SEA Subsystem Attribute Correction Based on Stiffness Mutipiers Jie Gao1,* and Yong Yang 1 Coege of hysics and Information

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,

More information

PHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I

PHYSICS LOCUS / / d dt. ( vi) mass, m moment of inertia, I. ( ix) linear momentum, p Angular momentum, l p mv l I 6 n terms of moment of inertia, equation (7.8) can be written as The vector form of the above equation is...(7.9 a)...(7.9 b) The anguar acceeration produced is aong the direction of appied externa torque.

More information

XI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com

XI PHYSICS. M. Affan Khan LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com XI PHYSICS M. Affan Khan LECTURER PHYSICS, AKHSS, K affan_414@ive.com https://promotephysics.wordpress.com [TORQUE, ANGULAR MOMENTUM & EQUILIBRIUM] CHAPTER NO. 5 Okay here we are going to discuss Rotationa

More information

Experimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure

Experimental Investigation and Numerical Analysis of New Multi-Ribbed Slab Structure Experimenta Investigation and Numerica Anaysis of New Muti-Ribbed Sab Structure Jie TIAN Xi an University of Technoogy, China Wei HUANG Xi an University of Architecture & Technoogy, China Junong LU Xi

More information

Application of the Finite Fourier Sine Transform Method for the Flexural-Torsional Buckling Analysis of Thin-Walled Columns

Application of the Finite Fourier Sine Transform Method for the Flexural-Torsional Buckling Analysis of Thin-Walled Columns IOSR Journa of Mechanica and Civi Engineering (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 3-334X, Voume 14, Issue Ver. I (Mar. - Apr. 17), PP 51-6 www.iosrjournas.org Appication of the Finite Fourier Sine Transform

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Bending Analysis of Continuous Castellated Beams

Bending Analysis of Continuous Castellated Beams Bending Anaysis of Continuous Casteated Beams * Sahar Eaiwi 1), Boksun Kim ) and Long-yuan Li 3) 1), ), 3) Schoo of Engineering, Pymouth University, Drake Circus, Pymouth, UK PL4 8AA 1) sahar.eaiwi@pymouth.ac.uk

More information

Candidate Number. General Certificate of Education Advanced Level Examination January 2012

Candidate Number. General Certificate of Education Advanced Level Examination January 2012 entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE Juan Huang, Ronghui Wang and Tao Tang Coege of Traffic and Communications, South China University of Technoogy, Guangzhou, Guangdong 51641,

More information

Measurement of acceleration due to gravity (g) by a compound pendulum

Measurement of acceleration due to gravity (g) by a compound pendulum Measurement of acceeration due to gravity (g) by a compound penduum Aim: (i) To determine the acceeration due to gravity (g) by means of a compound penduum. (ii) To determine radius of gyration about an

More information

CABLE SUPPORTED STRUCTURES

CABLE SUPPORTED STRUCTURES CABLE SUPPORTED STRUCTURES STATIC AND DYNAMIC ANALYSIS OF CABLES 3/22/2005 Prof. dr Stanko Brcic 1 Cabe Supported Structures Suspension bridges Cabe-Stayed Bridges Masts Roof structures etc 3/22/2005 Prof.

More information

SECTION A. Question 1

SECTION A. Question 1 SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second

More information

Faculty of Machine Building. Technical University of Cluj Napoca

Faculty of Machine Building. Technical University of Cluj Napoca Facuty of Machine Buiding Technica University of Cuj Napoca CONTRIBUTIONS TO THE CALCULATION AND ANALYSIS OF DYNAMIC ABSORBERS, WITH APPLICATIONS ON BALANCING MECHANICAL SYSTEMS PhD THESIS 11 Scientific

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

Finite element method for structural dynamic and stability analyses

Finite element method for structural dynamic and stability analyses Finite eement method for structura dynamic and stabiity anayses Modue-9 Structura stabiity anaysis Lecture-33 Dynamic anaysis of stabiity and anaysis of time varying systems Prof C S Manohar Department

More information

STABILITY OF THERMOELASTIC LAYERED COMPOSITE IN AXIAL MOVEMENT

STABILITY OF THERMOELASTIC LAYERED COMPOSITE IN AXIAL MOVEMENT 6th European Conference on Computationa Mechanics (ECCM 6) 7th European Conference on Computationa Fuid Dynamics (ECFD 7) 11-15 June 18, Gasgow, UK STABILITY OF THERMOELASTIC LAYERED COMPOSITE IN AXIAL

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p

More information

ISSN: [Alkhatib* et al., 5(10): October, 2016] Impact Factor: 4.116

ISSN: [Alkhatib* et al., 5(10): October, 2016] Impact Factor: 4.116 [Akhatib* et a., 5(10): October, 2016] Impact Factor: 4.116 IC Vaue: 3.00 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY VIBRATION ANALYSIS, SHAPE AND DESIGN OF MOTORCYCLE MOUNTING

More information

A Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem

A Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Numerical methods for PDEs FEM - abstract formulation, the Galerkin method

Numerical methods for PDEs FEM - abstract formulation, the Galerkin method Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona

More information

Dynamic equations for curved submerged floating tunnel

Dynamic equations for curved submerged floating tunnel Appied Mathematics and Mechanics Engish Edition, 7, 8:99 38 c Editoria Committee of App. Math. Mech., ISSN 53-487 Dynamic equations for curved submerged foating tunne DONG Man-sheng, GE Fei, ZHANG Shuang-yin,

More information

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution

More information

On some problems of synthesis of spatial five-bar hinged mechanisms with two degrees of freedom

On some problems of synthesis of spatial five-bar hinged mechanisms with two degrees of freedom Internationa Journa of Mechanica Engineering and Appications 014 (6): 104-110 Puished onine ecemer 18, 014 (http://www.sciencepuishinggroup.com/j/ijmea) doi: 10.11648/j.ijmea.014006.14 ISSN: 330-03X (Print)

More information

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics

Published in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of

More information

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation

Applied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation 22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements

More information

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams Unit 48: Structura Behaviour and Detaiing for Construction 4.1 Introduction Defection of Beams This topic investigates the deformation of beams as the direct effect of that bending tendency, which affects

More information

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier

More information

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006

arxiv:nlin/ v2 [nlin.cd] 30 Jan 2006 expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis Laboratory Exercise 1: Penduum Acceeration Measurement and Prediction Laboratory Handout AME 20213: Fundamentas of Measurements and Data Anaysis Prepared by: Danie Van Ness Date exercises to be performed:

More information

221B Lecture Notes Notes on Spherical Bessel Functions

221B Lecture Notes Notes on Spherical Bessel Functions Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,

More information

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled. imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,

More information

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА Милко Стоянов Милошев 1, Константин Савков Казаков 2 Висше Строително Училище Л. Каравелов - София COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS

More information

Modelling and simulation of a feed drive using the transmission line modelling technique

Modelling and simulation of a feed drive using the transmission line modelling technique Transactions on Engineering Sciences vo 44, 2003 WT Press, www.witpress.com, SSN 1743-3533 Modeing and simuation of a feed drive using the transmission ine modeing technique V.Y. Moreno-Castaneda, C. Pisaru,

More information

Transverse Anisotropy in Softwoods

Transverse Anisotropy in Softwoods Transverse Anisotropy in Softwoods Modeing and Experiments CARL MODÉN Licenciate Thesis Stockhom, Sweden 2006 TRITA-AVE 2006:30 ISSN 1651-7660 ISBN 91-7178-385-7 KTH Engineering Sciences Department of

More information

Strain Energy in Linear Elastic Solids

Strain Energy in Linear Elastic Solids Strain Energ in Linear Eastic Soids CEE L. Uncertaint, Design, and Optimiation Department of Civi and Environmenta Engineering Duke Universit Henri P. Gavin Spring, 5 Consider a force, F i, appied gradua

More information

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7 Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Numerical simulation of javelin best throwing angle based on biomechanical model

Numerical simulation of javelin best throwing angle based on biomechanical model ISSN : 0974-7435 Voume 8 Issue 8 Numerica simuation of javein best throwing ange based on biomechanica mode Xia Zeng*, Xiongwei Zuo Department of Physica Education, Changsha Medica University, Changsha

More information

1. Measurements and error calculus

1. Measurements and error calculus EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the

More information

COUPLED FLEXURAL TORSIONAL VIBRATION AND STABILITY ANALYSIS OF PRE-LOADED BEAMS USING CONVENTIONAL AND DYNAMIC FINITE ELEMENT METHODS

COUPLED FLEXURAL TORSIONAL VIBRATION AND STABILITY ANALYSIS OF PRE-LOADED BEAMS USING CONVENTIONAL AND DYNAMIC FINITE ELEMENT METHODS COUPLED FLEXURAL TORSIONAL VIBRATION AND STABILITY ANALYSIS OF PRE-LOADED BEAMS USING CONVENTIONAL AND DYNAMIC FINITE ELEMENT METHODS by Heenkenda Jayasinghe, B. Eng Aeronautica Engineering City University

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Post-buckling behaviour of a slender beam in a circular tube, under axial load

Post-buckling behaviour of a slender beam in a circular tube, under axial load Computationa Metho and Experimenta Measurements XIII 547 Post-bucking behaviour of a sender beam in a circuar tube, under axia oad M. Gh. Munteanu & A. Barraco Transivania University of Brasov, Romania

More information

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z

Bohr s atomic model. 1 Ze 2 = mv2. n 2 Z Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated

More information

ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS. Abstract

ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS. Abstract ANALYTICAL AND EXPERIMENTAL STUDY OF FRP-STRENGTHENED RC BEAM-COLUMN JOINTS Dr. Costas P. Antonopouos, University of Patras, Greece Assoc. Prof. Thanasis C. Triantafiou, University of Patras, Greece Abstract

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

Malaysian Journal of Civil Engineering 30(2): (2018)

Malaysian Journal of Civil Engineering 30(2): (2018) Maaysian Journa of Ci Engineering 3():331-346 (18) BUBNOV-GALERKIN METHOD FOR THE ELASTIC BUCKLING OF EULER COLUMNS Ofondu I.O. 1, Ikwueze E. U. & Ike C. C. * 1 Dept. of Mechanica and Production Engineering,

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur odue 2 naysis of Staticay ndeterminate Structures by the atri Force ethod Version 2 E T, Kharagpur esson 12 The Three-oment Equations- Version 2 E T, Kharagpur nstructiona Objectives fter reading this

More information

Research Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems

Research Article Building Infinitely Many Solutions for Some Model of Sublinear Multipoint Boundary Value Problems Abstract and Appied Anaysis Voume 2015, Artice ID 732761, 4 pages http://dx.doi.org/10.1155/2015/732761 Research Artice Buiding Infinitey Many Soutions for Some Mode of Subinear Mutipoint Boundary Vaue

More information

9. EXERCISES ON THE FINITE-ELEMENT METHOD

9. EXERCISES ON THE FINITE-ELEMENT METHOD 9. EXERCISES O THE FIITE-ELEMET METHOD Exercise Thickness: t=; Pane strain proem (ν 0): Surface oad Voume oad; 4 p f ( x, ) ( x ) 0 E D 0 0 0 ( ) 4 p F( xy, ) Interna constrain: rigid rod etween D and

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached

More information

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS

More information

Physics 506 Winter 2006 Homework Assignment #6 Solutions

Physics 506 Winter 2006 Homework Assignment #6 Solutions Physics 506 Winter 006 Homework Assignment #6 Soutions Textbook probems: Ch. 10: 10., 10.3, 10.7, 10.10 10. Eectromagnetic radiation with eiptic poarization, described (in the notation of Section 7. by

More information

Parallel-Axis Theorem

Parallel-Axis Theorem Parae-Axis Theorem In the previous exampes, the axis of rotation coincided with the axis of symmetry of the object For an arbitrary axis, the paraeaxis theorem often simpifies cacuations The theorem states

More information

Formulas for Angular-Momentum Barrier Factors Version II

Formulas for Angular-Momentum Barrier Factors Version II BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

Research Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method

Research Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method Science and Technoogy of Nucear Instaations Voume 213, Artice ID 261327, 6 pages http://dx.doi.org/1.1155/213/261327 Research Artice Soution of Point Reactor Neutron Kinetics Equations with Temperature

More information

1 Equations of Motion 3: Equivalent System Method

1 Equations of Motion 3: Equivalent System Method 8 Mechanica Vibrations Equations of Motion : Equivaent System Method In systems in which masses are joined by rigid ins, evers, or gears and in some distributed systems, various springs, dampers, and masses

More information

ANALYSIS OF THERMOELASTIC WAVES ON THE BASIS OF AN EXACT THEORY

ANALYSIS OF THERMOELASTIC WAVES ON THE BASIS OF AN EXACT THEORY . K. L. VERMA ANALYSIS OF THERMOELASTIC WAVES ON THE BASIS OF AN EXACT THEORY. DEPARTMENT OF MATHEMATICS GOVERNMENT POST GRADUATE COLLEGE HAMIRPUR (H.P.) 7700, INDIA ABSTRACT: In this artice wave propagation

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

International Journal of Advance Engineering and Research Development

International Journal of Advance Engineering and Research Development Scientific Journa of Impact Factor (SJIF): 4.4 Internationa Journa of Advance Engineering and Research Deveopment Voume 3, Issue 3, March -206 e-issn (O): 2348-4470 p-issn (P): 2348-6406 Study and comparison

More information