Nonlinear Flexural - Torsional Dynamic Analysis of Beams of Arbitrary Cross Section by BEM
|
|
- April Lloyd
- 5 years ago
- Views:
Transcription
1 Nonlinear Flexural - Torsional Dynamic Analysis of Beams of Arbitrary Cross Section by BEM E.J. Sapountzakis, I.C. Dikaros To cite this version: E.J. Sapountzakis, I.C. Dikaros. Nonlinear Flexural - Torsional Dynamic Analysis of Beams of Arbitrary Cross Section by BEM. International Journal of Non-Linear Mechanics, Elsevier, 011, < /j.ijnonlinmec >. <hal > HAL Id: hal Submitted on 8 Oct 01 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.
2 Author s Accepted Manuscript Nonlinear Flexural - Torsional Dynamic Analysis of Beams of Arbitrary Cross Section by BEM E.J. Sapountzakis, I.C. Dikaros PII: S (11) DOI: doi: /j.ijnonlinmec Reference: NLM 1815 To appear in: International Journal of Non- Linear Mechanics Received date: 1 September 010 Revised date: 7 December 010 Accepted date: 7 February 011 Cite this article as: E.J. Sapountzakis and I.C. Dikaros, Nonlinear Flexural - Torsional Dynamic Analysis of Beams of Arbitrary Cross Section by BEM, International Journal of Non-Linear Mechanics, doi: /j.ijnonlinmec This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting galley proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
3 NONLINEAR FLEXURAL - TORSIONAL DYNAMIC ANALYSIS OF BEAMS OF ARBITRARY CROSS SECTION BY BEM by E.J. SAPOUNTZAKIS 1 and I.C. DIKAROS cvsapoun@central.ntua.gr dikarosgiannis@gmail.com Institute of Structural Analysis and Antiseismic Research, School of Civil Engineering, National Technical University of Athens, Zografou Campus, GR Athens, Greece. Key words: Flexural-Torsional Analysis; Dynamic analysis; Wagner s coefficients; Nonlinear analysis; Shortening Effect; Boundary element method Abstract In this paper, a boundary element method is developed for the nonlinear flexuraltorsional dynamic analysis of beams of arbitrary, simply or multiply connected, constant cross section, undergoing moderately large deflections and twisting rotations under general boundary conditions, taking into account the effects of rotary and warping inertia. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading in both directions as well as to twisting 1 Associate Professor, School of Civil Engineering, National Technical University, Zografou Campus, GR , Athens, Greece. Doctoral Student, School of Civil Engineering, National Technical University, Zografou Campus, GR , Athens, Greece.
4 and/or axial loading. Four boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to the angle of twist and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique leads to a system of nonlinear coupled Differential Algebraic Equations (DAE) of motion, which is solved iteratively using the Petzold- Gear Backward Differentiation Formula (BDF), a linear multistep method for differential equations coupled to algebraic equations. The geometric, inertia, torsion and warping constants are evaluated employing the Boundary Element Method. The proposed model takes into account, both the Wagner s coefficients and the shortening effect. Numerical examples are worked out to illustrate the efficiency, wherever possible the accuracy, the range of applications of the developed method as well as the influence of the nonlinear effects to the response of the beam. 1. INTRODUCTION In engineering practice the dynamic analysis of beam-like continuous systems is frequently encountered. Such structures often undergo arbitrary external loading, leading to the formulation of the flexural-torsional vibration problem. The complexity of this problem increases significantly in the case the cross section s centroid does not coincide with its shear center (asymmetric beams). Furthermore when arbitrary torsional boundary conditions are applied either at the edges or at any other interior point of a beam due to construction requirements, the beam under the action of general twisting loading is leaded to nonuniform torsion. Moreover, since requirement of weight saving is a major aspect in the design of structures, thin-walled elements of arbitrary cross section and low flexural and/or torsional stiffness are extensively used. Treating displacements and angles of rotation of these elements as being small, leads
5 in many cases to inadequate prediction of the dynamic response; hence the occurring nonlinear effects should be taken into account. This can be achieved by retaining the nonlinear terms in the strain displacement relations (finite displacement small strain theory). When finite displacements are considered, the flexural-torsional dynamic analysis of bars becomes much more complicated, leading to the formulation of coupled and nonlinear flexural, torsional and axial equations of motion. When the displacement components of a member are small, a wide range of linear analysis tools, such as modal analysis, can be used and some analytical results are possible. As these components become larger, the induced geometric nonlinearities result in effects that are not observed in linear systems. In such situations the possibility of an analytical solution method is significantly reduced and is restricted to special cases of beam boundary conditions or loading. During the past few years, the nonlinear dynamic analysis of beams undergoing large deflections has received a good amount of attention in the literature. More specifically, Rozmarynowski and Szymczak in [1] studied the nonlinear free torsional vibrations of axially immovable thin-walled beams with doubly symmetric open cross section, employing the Finite Element Method. In this research only free vibrations are examined; the solution is provided only at points of reversal of motion (not in the time domain), no general axial, torsional or warping boundary conditions (elastic support case) are studied, while some nonlinear terms related to the finite twisting rotations as well as the axial inertia term are ignored. Crespo Da Silva in [-3] presented the nonlinear flexural-torsional-extensional vibrations of Euler-Bernoulli doubly symmetric thin-walled closed cross section beams, primarily focusing on flexural vibrations and neglecting the effect of torsional warping. Pai and Nayfeh in [4-6] studied also the nonlinear flexural-torsional-extensional vibrations of metallic
6 and composite slewing or rotating closed cross section beams, primarily focusing to flexural vibrations and neglecting again the effect of torsional warping. Simo and Vu- Quoc in [7] presented a FEM solution to a fully nonlinear (small or large strains, hyperelastic material) three dimensional rod model including the effects of transverse shear and torsion-warping deformation based on a geometrically exact description of the kinematics of deformation. Qaisi in [8] obtained nonlinear normal modes of free vibrating geometrically nonlinear beams of various edge conditions employing the harmonic balance analytical method. Moreover, Pai and Nayfeh in [9] studied a geometrically exact nonlinear curved beam model for solid composite rotor blades using the concept of local engineering stress and strain measures and taking into account the in-plane and out-of-plane warpings. Di Egidio et al. in [10-11] presented also a FEM solution to the nonlinear flexural-torsional vibrations of shear undeformable thin-walled open beams taking into account in-plane and out-of-plane warpings and neglecting warping inertia. In this paper, the torsional-extensional coupling is taken into account but the inextensionality assumption leads to the fact that the axial boundary conditions are not general. Mohri et al. in [1] proposed a FEM solution to the linear vibration problem of pre- and post- buckled thin-walled open cross section beams, neglecting warping and axial inertia, considering geometrical nonlinearity only for the static loading and presenting examples of bars subjected to free vibrations and special boundary conditions. Machado and Cortinez in [13] presented also a FEM solution to the linear free vibration analysis of bisymmetric thin-walled composite beams of open cross section, taking into account initial stresses and deformations considering geometrical nonlinearity only for the static loading and presenting examples of bars subjected to special boundary conditions. Avramov et al. [14-15] studied the free flexural-torsional vibrations of
7 beams and obtained nonlinear normal modes by expansion of the equations of motion employing the Galerkin technique and neglecting the cross section warping. Lopes and Ribeiro [16] studied also the nonlinear flexural-torsional free vibrations of beams employing a FEM solution and neglecting the longitudinal and rotary inertia as well as the cross-section warping. Duan [17] presented a FEM formulation for the nonlinear free vibration problem of thin-walled curved beams of asymmetric crosssection based on a simplified displacement field. Finally, the boundary element method has also been used for the nonlinear flexural [18-0] and torsional [1] dynamic analysis of only doubly symmetric beams. To the authors knowledge the general problem of coupled nonlinear flexural torsional free or forced vibrations of asymmetric beams has not yet been presented. In this paper, a boundary element method is developed for the nonlinear flexural-torsional dynamic analysis of beams of arbitrary, simply or multiply connected, constant cross section, undergoing moderately large deflections and twisting rotations under general boundary conditions, taking into account the effects of rotary and warping inertia. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading in both directions as well as to twisting and/or axial loading. Four boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to the angle of twist and solved using the Analog Equation Method [], a BEM based method. Application of the boundary element technique leads to a system of nonlinear coupled Differential Algebraic Equations (DAE) of motion, which is solved iteratively using the Petzold-Gear Backward Differentiation Formula (BDF) [3], a linear multistep method for differential equations coupled to algebraic equations. The geometric, inertia, torsion and warping constants are evaluated employing the Boundary Element
8 Method. The essential features and novel aspects of the present formulation compared with previous ones are summarized as follows. i. The cross section is an arbitrarily shaped thin- or thick-walled one. The formulation does not stand on the assumption of a thin-walled structure and therefore the cross section s torsional and warping rigidities are evaluated exactly in a numerical sense. ii. The beam is subjected to arbitrarily distributed or concentrated transverse loading in both directions as well as to twisting and axial loading. iii. The beam is supported by the most general boundary conditions including elastic support or restraint. iv. The effects of rotary and warping inertia are taken into account on the nonlinear flexural-torsional dynamic analysis of asymmetric beams subjected to arbitrary loading and boundary conditions. v. The transverse loading can be applied at any arbitrary point of the beam cross section. The eccentricity change of the transverse loading during the torsional beam motion, resulting in additional torsional moment is taken into account. vi. The proposed model takes into account the coupling effects of bending, axial and torsional response of the beam as well as the Wagner s coefficients and the shortening effect. vii. The proposed method employs a BEM approach (requiring boundary discretization for the cross sectional analysis) resulting in line or parabolic elements instead of area elements of the FEM solutions (requiring the whole cross section to be discretized into triangular or quadrilateral area elements), while a small number of line elements are required to achieve high accuracy.
9 Numerical examples are worked out to illustrate the efficiency, wherever possible the accuracy, the range of applications of the developed method as well as the influence of the nonlinear effects to the response of the beam.. STATEMENT OF THE PROBLEM Let us consider a prismatic beam of length l (Fig.1), of constant arbitrary cross section of area. A.. The homogeneous isotropic and linearly elastic material of the beam s cross section, with modulus of elasticity. E., shear modulus.. and Poisson s ratio.ν. occupies the two dimensional multiply connected region. Ω. of the y,z plane and is bounded by the.. Γ j ( j = 1,,...,K ).. boundary curves, which are piecewise smooth, i.e. they may have a finite number of corners. In Fig. 1 CYZ is the principal bending coordinate system through the cross section s centroid C, while y C, C z are its coordinates with respect to the Syz shear system of axes through the cross section s shear center S, with axes parallel to those of the.. system. The beam is subjected to the combined action of the arbitrarily distributed or concentrated, time dependent and conservative axial loading p = p ( X,t) along X direction, twisting moment.. along x direction and transverse loading p = p ( x,t), p = p ( x,t) X X y y z z acting along the y and z directions, applied at distances y p y, z p y and y p z, z p z, with respect to the Syz shear system of axes, respectively (Fig. 1b). Under the action of the aforementioned loading, the displacement field of an arbitrary point of the cross section can be derived with respect to those of the shear center as [1] u( x,y,z,t) = u( x,t) ( y y ) ( ) ( ) ( ) ( ) P ( ). C θ Z x,t + z z C θ Y x,t + θ x x,t ϕ S y,z. (1a)
10 ( ) = ( ) ( θ ( )) ( θ ( )) v x,y,z,t v x,t z sin x x,t y 1 cos x x,t ( ) = ( ) + ( θ ( )) ( θ ( )) w x,y,z,t w x,t ysin x x,t z 1 cos x x,t (1b) (1c) ( x,t) v ( x,t) sin( ( x,t) ) w ( x,t) cos ( ( x,t) ) θ = θ θ (1d) Y x x ( x,t) v ( x,t) cos ( ( x,t) ) w ( x,t) sin( ( x,t) ) θ = θ + θ (1e) Z x x where u, v, w are the axial and transverse beam displacement components with 1 = u x, y,z,t da denotes the A A respect to the Syz shear system of axes; u ( x,t) ( ) average axial displacement of the cross section [4] and v= v( x,t), w w( x,t) the corresponding components of the shear center S ; θ ( x,t ), ( x,t) Y Z = are θ are the angles of rotation of the cross section due to bending, with respect to its centroid; θ ( x,t) denotes the rate of change of the angle of twist ( x,t ) x θ regarded as the P torsional curvature and ϕ S is the primary warping function with respect to the shear center S [5]. Employing the strain-displacement relations of the three-dimensional elasticity for moderate displacements [6, 7] x ε γ xx xy u 1 u v w = x x x x v u u u v v w w = x y x y x y x y (a) (b) γ xz w u u u v v w w = x z x z x z x z (c) ε yy = εzz = γyz = 0 (d)
11 the following strain components can be easily obtained ε γ xx xy u 1 v w + + x x x v u v v w w x y x y x y (3a) (3b) γ xz w u v v w w x z x z x z (3c) ε yy = εzz = γyz = 0 (3d) where it has been assumed that for moderate displacements ( u ) ( u ) ( u u )( ) ( u ) u x y + x ( y) x,, x. Substitu- ( u )( u ) ( u ) + ( u ) x z x z the displacement components (1a-1e) to the strain-displacement relations (3), the strain components can be written as ting ( )( ) ( )( ) P θ x φs θ x ( zc ( v cosθx w sinθx) yc ( v sinθx w cosθx) ) ε = u + z z v sinθ w cosθ y y v cosθ + w sinθ + xx C x x C x x ( )( θ x ) 1 + v + w + y + z P φ S γ xy = εxy = z θ x y P φ S γ xz = εxz = + y θ x z (4a) (4b) (4c)
12 Considering strains to be small and employing the second Piola Kirchhoff stress tensor, the non vanishing stress components are defined in terms of the strain ones as S xx E 0 0 ε xx Sxy = 0 G 0 γ xy S 0 0 G γ xz xz (5) where * E is obtained from Hooke s stress-strain law as E = E( 1 ν ) ( 1+ ν )( 1 ν ). If the assumption of plane stress condition is made, the above expression is reduced in E E = [8], while in beam formulations, E is frequently considered instead of 1 ν * E (E E) [8, 9]. This last consideration has been followed throughout the paper, while any other reasonable expression of * E could also be used without any difficulty. Substituting eqns. (4) into eqns. (5), the non vanishing stress components are obtained as ( )( θ θ ) ( )( θ θ ) P x S x ( zc ( v cos x w sin x) yc ( v sin x w cos x) ) Sxx = E u + z zc v sin x w cos x y yc v cos x + w sin x + + θ φ θ θ + θ + θ θ + ( )( θ x ) 1 + v + w + y + z P φ S Sxy = G θ x z y P φ S Sxz = G θ x + y z (6a) (6b) (6c)
13 In order to establish the nonlinear equations of motion, the principle of virtual work δw + δw = δw (7) int mass ext where ( ) δw = S δε + S δγ + S δγ dv (8a) int V xx xx xy xy xz xz V ( ) δwmass = ρ u δu + v δv + w δw dv (8b) ( y z ) δw = p δu + p δv + p δw + m δθ dx (8c) ext l X C y p z p x x under a Total Lagrangian formulation, is employed. In the above equations, δ () denotes virtual quantities, ( ) denotes differentiation with respect to time, V is the volume of the beam, u C is the axial displacement of the centroid and v p, v y p z are the transverse displacements of the points where the loads p y, p z, respectively, are applied. It is worth here noting that the aforementioned expression of the external work (eqn. (8c)) takes into account the change of the eccentricity of the external conservative transverse loading, arising from the cross section torsional rotation, inducing additional (positive or negative) torsional moment. This effect may influence substantially the torsional response of the beam. Moreover, the stress resultants of the beam can be defined as N = S dω (9a) Ω xx
14 M = S Zd (9b) Υ Ω xx Ω M Z = S YdΩ (9c) Ω xx P P ϕ S ϕ S Mt = Sxy z Sxz y dω Ω + + y z (9d) P M w = SxxϕSdΩ Ω (9e) R ( ) M = S y + z dω Ω xx (9f) where Μ t is the primary twisting moment [5] resulting from the primary shear stress distribution S xy, S xz, M w is the warping moment due to torsional curvature and M R is a higher order stress resultant. Substituting the expressions of the stress components (6) into equations (9a-9f), the stress resultants are obtained as 1 IS N = EA u v w θ + x + + A θ θ + θ + θ θ ( z ( v cos w sin ) y ( v sin w cos )) x C x x C x x ( θ θ β θ ) Y Y x x Z x (10a) M = EI w cos v sin (10b) ( θ θ β θ ) M = EI v cos + w sin (10c) Z Z x x Y x Mt = GItθ x (10d) ( θ β θ ) M = EC + ω (10e) w S x x
15 IS MR = N EIZβY ( v cosθx + w sinθx) EIYβZ ( w cosθx v sinθx) + A 1 I S + ECSβωθ x + E IR θ x A (10f) where the area A, the polar moment of inertia I S with respect to the shear center S, the principal moments of inertia I Y, I Z with respect to the cross section s centroid, the fourth moment of inertia I R with respect to the shear center S, the torsion constant I t and the warping constant C S with respect to the shear center S, are given as A = d Ω (11a) Ω ( ) IS = y + z dω (11b) Ω IY IZ = Z dω (11c) Ω = Y dω (11d) Ω ( ) IR = y + z dω (11e) Ω C S P ( ) S = d Ω (11f) Ω ϕ P P ϕs ϕ S It = y z y z dω Ω + + (11g) z y while the Wagner s coefficients β Z, β Y and β ω are given as β Z 1 = z z y + z dω (1a) I Ω Y ( C)( )
16 β Y 1 = y y y + z dω (1b) I Ω Z ( C)( ) S ( ) 1 P βω = y + z ϕ S dω (1c) C Ω Using the expressions of strain obtained in equations (4), the definitions of the stress resultants given in equations (9) and applying the principle of virtual work (eqn. (7)), the equations of motion of the beam can be derived as N + ρ Au = p X (13a) ( NF ) v1 ( M Z cos ) x ( MY sin x ) Av F v x F v3 x ( F v4 x ) + θ + θ + ρ ρ θ + ρ θ ρ θ ( ρ θ ) ρ ( θ ) ρ ( θ ) F v5 x Fv6 xv w + Fv7 xw + v = py (13b) ( NF ) w1 ( MY cos ) x ( M Z sin x ) Aw F w x F w3 x ( F w4 x ) θ + θ + ρ ρ θ ρ θ ρ θ + ( ρ θ ) ρ ( θ ) ρ ( θ ) + F w5 x Fw6 xw + v + Fw7 w xv = pz (13c) NF + NF + NF + M F + M F M M M ( ) ( ) ( ) θ1 θ θ3 Y θ4 Z θ5 t w Rθ x ( F v F w ) F F ( v w ) F ( w v ) F ρ + + ρ θ + ρ + θ + ρ θ ρ θ θ6 θ7 θ8 x θ9 x θ10 x θ11 x ( ) z y z y ρc θ = m + p y cosθ p z cosθ p z sinθ p y sinθ (13d) S x x z p x y p x z p x y p x
17 where the expressions of the stress resultants are given from equations (10) and F vi ( i = 1,,...,7 ), F wi ( i 1,,...,7 ) = and i F θ ( i 1,,...,11) = are functions of v, w, θ x and their derivatives with respect to x, given in the Appendix A. Equations (13) are coupled and highly complicated. This set of equations can be simplified if the approximate expressions [1] x θx θ cosθx 1 = 1! 3 3 θx θx x x = x sinθ θ θ 3! 6 (14a) (14b) are employed. Thus, using the aforementioned approximations, neglecting the term ρ Au of equation (13a) denoting the axial inertia of the beam, employing the expressions of the stress resultants (eqns. (10)) and ignoring the nonlinear terms of the fourth or greater order [1], the governing partial differential equations of motion for the beam at hand can be written as IS EA u ww vv G p A θ θ = x x u1 X (15a) ( ) β β ρ( ) EI v NG + EI EI G + EI G + EI G + Av + G G = Z v1 Z Y v Z Y v3 Y Z v4 v5 v6 = py Gv7pX (15b) ( ) β β ρ( ) EI w NG + EI EI G + EI G + EI G + Aw + G G = Y w1 Z Y w Z Y w3 Y Z w4 w5 w6 = pz Gw7 px (15c)
18 3 IS ECSθx GItθx E IR θ x θ x NGθ1+ ( EIZ EIY ) Gθ + EIZβYGθ3 + A ( ) z z + EI β G + ρ G θ G v + G w + G C θ = m + p y p z + Y Z θ4 θ5 x θ6 θ7 θ8 S x x z p y p + G p + G p G p (15d) θ9 y θ10 z θ11 X where G u1, vi G ( i 1,,...,7 ) =, wi G ( i 1,,...,7 ) = and i G θ ( i 1,,...,11) = are functions of v, w, θ x and their derivatives with respect to. x., t, given in the Appendix A, while the expression of the axial stress resultant N is given as 1 IS N EA u v w x x zc w x v yc v x w A θ = θ θ + + θ ( ( ) ( )) (16) The above governing differential equations (eqns. (15)) are also subjected to the initial conditions ( x ( 0,l) ) ( ) = ( ) u( x,0) = u ( x) u x,0 u x 0 ( ) = ( ) v( x,0) = v ( x) v x,0 v x 0 ( ) = ( ) w( x,0) = w ( x) w x,0 w x θ x 0 ( x,0) = θ ( x) θ ( x,0) = θ ( x) x0 (17a,b) 0 (18a,b) 0 (19a,b) x 0 (0a,b) x0 together with the corresponding boundary conditions of the problem at hand, which are given as
19 ( ) α N ( x,t) a u x,t + = α (1) 1 3 ( ) + ( ) = ( x,t) ( x,t) β v x,t β V x,t β 1 y 3 β θ + β Μ = β (a,b) 1 Z Z 3 ( ) + ( ) = ( x,t) ( x,t) γ w x,t γ V x,t γ 1 z 3 γ θ + γ Μ = γ (3a,b) 1 Y Y 3 ( x,t) + M ( x,t) = ( x,t) M ( x,t) δ θ δ δ 1 x t 3 δ θ + δ = δ (4a,b) 1 x w 3 at the beam ends x = 0,l, where V y, V z and M Z, M Y are the reactions and bending moments with respect to y, z or to Y, Z axes, respectively, given by the following relations (ignoring again the nonlinear terms of the fourth or greater order) ( θ θ θ ) EIZ ( v θx w θ x w θx v v θxθ x βyθ xθ x ) 3 EIY ( w θx w θ x v θxθ x v θx βzθ x βzθxθ xθ x ) V = N v z y + y C x C x x ( + θ θ θ ) + EIY ( w θx v θx w v θ x w θxθ x βzθ xθ x ) 3 EIZ ( v θ x w θx v θx w θxθ x βyθxθ xθ x βyθ x ) V = N w y z z C x C x x (5a) (5b) ( θ βθ θ ) ( θ θ βθ θ) M = EI w + v v + EI w + v + (5c) Z Z x Y x x Y x x Z x x ( θ βθ θ θ ) ( βθ θ θ) M = EI w + v + EI w + w + v (5d) Y Z x Y x x x Y Z x x x
20 while M t and bar, respectively, given as M w are the torsional and warping moments at the boundaries of the IS Mt = GItθ x ECSθ x + N zcw θx + ycw ycv θx zcv + θ x A + + EI ( v w ) + EI ( w + v ) + E I A 1 IS 3 Z βy θ x θ x θx YβZ θ x θ x θx R θ x ( θ β θ ) w S x x (6a) M = EC + ω (6b) Finally, αk, βk, βk, γk, γk, δk, δ k (k = 1,,3) are time dependent functions specified at the boundaries of the bar ( x = 0,l ). The boundary conditions (1)-(4) are the most general boundary conditions for the problem at hand, including also the elastic support. It is apparent that all types of the conventional boundary conditions (clamped, simply supported, free or guided edge) can be derived from these equations by specifying appropriately these functions (e.g. for a clamped edge it is α1 = β1 = β1 = γ1 = γ1 = δ1 = δ1 = 1, α = α3 = β = β3 = γ = γ3 = δ = δ 3 = β β = γ = γ = δ = δ = ) INTEGRAL REPRESENTATIONS-NUMERICAL SOLUTION 3.1 For the axial displacement u( x,t), the transverse displacements vx,t ( ), wx,t ( ) and the angle of twist θ x ( x,t)
21 According to the precedent analysis, the nonlinear flexural-torsional vibration problem of a beam reduces in establishing the axial displacement component u( x,t ) having continuous partial derivatives up to the second order and the transverse displacement components v( x,t ), w( x,t ) and the angle of rotation ( x,t ) θ having continuous partial derivatives up to the fourth order with respect to x and up to the second order with respect to t, satisfying the nonlinear initial boundary value problem described by the coupled governing differential equations of motion (eqns. (15)) along the beam, the initial conditions (eqns. (17)-(0)) and the boundary conditions (eqns. (1)-(4)) at the beam ends x = 0,l. Eqns. (15) and (17)-(4) are solved using the Analog Equation Method [] as it is developed for hyperbolic differential equations [30]. According to this method, let u( x,t ), v( x,t ), w( x,t ) and ( x,t ) θ be the sought solutions of the aforementioned x problem. Setting as u ( x,t) = u( x,t), u ( x,t) = v( x,t), u ( x,t) w( x,t) ( ) θ ( ) 4 x 1 3 x = and u x,t = x,t and differentiating with respect to x these functions two and four times, respectively, yields u 1 = q1 ( x,t) x 4 u i 4 i x ( ) ) = q x,t, ( i =,3,4 (7a,b) Eqns. (7) are quasi-static, i.e. the time variable appears as a parameter. They indicate that the solution of eqns. (15) and (17)-(4) can be established by solving eqns. (7) under the same boundary conditions (eqns. (1)-(4)), provided that the fictitious load distributions qi ( x,t ) ( i 1,,3,4) = are first established. These distributions can be determined using BEM. Following the procedure presented in [30] and employing the
22 constant element assumption for the load distributions q i along the L internal beam elements (as the numerical implementation becomes very simple and the obtained results are of high accuracy), the integral representations of the displacement components i u ( i 1,,3,4) = and their derivatives with respect to x when applied for the beam ends (0,l ), together with the boundary conditions (1)-(4) are employed to express the unknown boundary quantities u ( ζ,t), u, ( ζ,t), u, (,t) i xxx ( ) u, ζ,t ( ζ = 0,l equations is obtained i ix ixx ζ and ) in terms of q i. Thus, the following set of 8 nonlinear algebraic E d 0 0 nl D 1 α nl D β 0 E 0 0 d + 0 = E33 0 d E44 d4 nl D3 γ nl D δ 4 (8) where E 11 F1 E1 E13 = 0 D D3 (9a)
23 E E E F3 E35 E36 E37 E38 F4 0 E46 E47 E48 = D54 D55 D56 D57 D58 D64 D65 D66 D67 D68 F7 E79 E711 E71 E713 F8 0 E811 E81 E813 = D99 D9 10 D9 11 D9 1 D9 13 D109 D1010 D1011 D101 D1013 F11 E1115 E1116 E1117 E1118 F1 0 E1 16 E1 17 E1 18 = D13 14 D13 15 D13 16 D13 17 D13 18 D14 14 D14 15 D14 16 D14 17 D14 18 (9b) (9c) (9d) T T { 3 3}, { 3 3}, { 3 3} β= β β γ = γ γ δ= δ δ (9e,f,g) T The matrices D to D 14 18, are rectangular known matrices including the values of the functions a j, β j, βj, γ j, γ j, δ j, δ j ( j = 1,) of eqns. (1)-(4); α 3, β 3, β 3, γ 3, γ 3, δ 3, δ 3 are 1 known column matrices including the boundary values of the functions a3, β3, β3, γ3, γ3, δ3, δ 3, of eqns. (1)-(4); E 1 to E 1 18 are rectangular known coefficient matrices resulting from the values of kernels at the bar ends; F1, F3, F4, F7, F8, F11, F 1 are L rectangular known matrices originating from the integration of kernels on the axis of the beam; nl ( = ) nl D 1 is a 1and D i i,3,4 are 4 1 known column matrices containing the nonlinear terms included in the expressions of the boundary conditions (eqns. (1)-(4)). Finally
24 { ˆ ˆ, } T d = q u u (30a) x T { ˆ ˆ, ˆ, ˆ, }, i,3,4) d = q u u u u ( = (30b) i i i i x i xx i xxx are generalized unknown vectors, where u ˆ = { u ( 0,t) u ( l,t) } T, ( i 1,,3,4) i i i u ( ) ( ) T i 0,t ui l,t u ˆ,, ( i 1,,3,4) i x = x x u ( ) ( ) T i 0,t ui l,t u ˆ,, ( i,3,4) = x x i xx 3 3 u ( ) ( ) T i 0,t ui l,t u ˆ,, ( i,3,4) = x x i xxx 3 3 = (31a) = (31b) = (31c) = (31d) are vectors including the two unknown time dependent boundary values of the respective boundary quantities and = { q i q i... q i } T q ( i 1,,3,4) i 1 N = are vectors including the L unknown time dependent nodal values of the fictitious loads. Discretization of the integral representations of the unknown quantities u i ( i = 1,,3,4) inside the beam ( x ( 0,l) points yields ) and application to the L collocation nodal 0 1 = ˆ 1+ 1ˆ 1, x u A q C u C u (3a) 1 1x, = ˆ 1x, u A q C u (3b)
25 u1xx, = q 1 (3c) 0 i = i i + 0ˆ i + 1ˆ i, x + ˆ i, xx + 3ˆ i, xxx u A q C u C u C u C u, ( i,3,4) 1 ix, = i i+ 0ˆ ix, + 1ˆ ixx, + ˆ ixxx, u A q C u C u C u, ( i,3,4) ixx, = i i+ 0ˆ ixx, + 1ˆ ixxx, u A q C u C u, ( i,3,4) 3 i, xxx = i + 0ˆ i, xxx u A q C u, ( i,3,4) ixxxx i = (3d) = (3e) = (3f) = (3g) u, = q, ( i=,3,4) (3h) where, i j 1 k A A ( i = 0, 1), ( j = 0, 1,, 3), ( k,3,4) C 0, C 1, C 1, C, C 3 are L known matrices and u i, are time dependent vectors including the values of u ( ) = are L L known matrices; i u ix,, i, xx u, u ixxx,, u i, xxxx x,t and their derivatives at the L nodal points. Eqns. (3a, b, d, e, f, g) can be assembled more conveniently as 0 1 = 1 1 u H d (33a) 1 1x, = 1 1 u H d (33b) 0 i = i i u H d, ( i,3,4) 1 ix, = i i u H d, ( i,3,4) i, xx= i i u H d, ( i,3,4) 3 i, xxx = i i u H d, ( i,3,4) = (33c) = (33d) = (33e) = (33f)
26 where, i j 1 k H H ( i = 0, 1), ( j = 0, 1,, 3), ( k =,3,4) are L ( L 4) L ( L+ 8) known matrices, respectively arising from C 3. i A 1, + and j A, C 0, C 1, C 1, C, Applying eqns. (15) to the L collocation points and employing eqns. (33), 4 L semidiscretized nonlinear equations of motion are formulated as d 1 d1 d d d 3 d3 d 4 d4 nl nl [ M] + [ K] + { m } + { k } = {} f (34) where M, K, f are generalized mass matrix, stiffness matrix and force vector, respectively, while nl m, nl k are nonlinear generalized mass vector and stiffness vector, respectively, containing all the nonlinear terms of the semidiscretized equations of motion. It is noted that the coefficients of the mass matrix M corresponding to the generalized vector d 1 are equal to zero as the axial inertia of the beam has been neglected. Equations (34) with equations (8) form a system of. 4 L 8 + equations with respect to the generalized unknown vectors d ( = ) i i 1,,3,4 Eqns. (33a,c) when combined with eqns. (17a)-(0a) yield the following 4 L linear equations with respect to d 1, d, d 3, d 4 for t = 0 0 Hd 1 ( 0 ) = u 0 (35a) 0 Hd ( 0 ) = v 0 (35b)
27 0 Hd 3 3 ( 0 ) = w 0 (35c) 0 Hd 4 4 ( 0 ) = θ x0 (35d) The above equations, together with eqns. (8) written for t = 0, form a set of 4 L+ 8 nonlinear algebraic equations which are solved to establish the initial conditions d 1 ( 0 ), d ( 0), d 3 ( 0), d 4 ( 0). Similarly, eqns. (33a,c) when combined with eqns. (17b)-(0b) yield the following 4 L linear equations with respect to d 1, d, d 3, d 4 for t = 0 ( ) 0 Hd = u 0 (36a) ( ) Hd 0 0 = v (36b) 0 ( ) Hd 0 0 = w (36c) ( ) Hd 0 0 = θ (36d) 4 4 x0 The above equations, together with the 8 equations resulting after differentiating eqns. (8) with respect to time and writing them for t = 0, form a set of 4 L+ 8 algebraic equations, from which the initial conditions d 1 ( 0 ), d ( 0 ), d 3 ( 0 ), d 4 ( 0 ) are established. The aforementioned initial conditions along with eqns. (8), (34) form an initial value problem of Differential-Algebraic Equations (DAE), which can be solved using any efficient solver. In this study, the Petzold-Gear Backward Differentiation Formula (BDF) [3], which is a linear multistep method for differential equations coupled to algebraic equations, is employed. For this case, the method is applied after
28 introducing new variables to reduce the order of the system [31] and after differentiating eqns. (8) with respect to time, in order to obtain an equivalent system with a value of system index ind = 1 [3]. 3. For the primary warping function P ϕ S The numerical solution for the evaluation of the displacement and rotation components assume that the warping C S and torsion I t constants given from eqns. (10e), (10d) are already established. Eqns. (11f), (11g) indicate that the evaluation of P the aforementioned constants presumes that the primary warping function ϕ S at any P interior point of the domain Ω of the cross section of the beam is known. Once ϕ S is established, C S and I t constants are evaluated by converting the domain integrals into line integrals along the boundary employing the following relations P S ϕ K+ 1 CS = B ds on Γ = Γ Γ j (37a) n P P ( ) ( ) It = yz zϕs ny + y z + yϕs n z ds Γ on Γ = Γ j (37b) j= 1 K+ 1 j= 1 and using constant boundary elements for the approximation of these line integrals. In eqns. (37a,b) n y, n z are the direction cosines, while B( y,z ) is a fictitious function defined as the solution of the following Neumann problem P B ϕ S = in Ω (38a)
29 B = 0 n K+ 1 on Γ = Γ (38b) j= 1 j P The evaluation of the primary warping function ϕ S and the fictitious function B( y,z ) is accomplished using BEM as this is presented in [5, 3, 33]. Moreover, since the torsion and warping constants of the arbitrary beam cross section are evaluated employing the boundary element method, using only boundary integration, the domain integrals for the evaluation of the area, the bending, the fourth and the polar moments of inertia and the Wagner s coefficients β Z, β Y and β ω given from expressions (1) have to be converted to boundary line integrals. This can be achieved using integration by parts, the Gauss theorem and the Green identity. Thus, the aforementioned quantities can be written as (( C) y ( C) z) 1 A = y-y n + z z n ds (39a) Γ 3 (( ) ) 1 IY = z zc nz ds (39b) Γ 3 3 (( ) ) 1 IZ = y yc ny ds (39c) Γ ( ) ( )( ) 1 1 IR = y ny + z nz + y z yny + znz ds (39d) Γ ( ) 1 IS = y ny + z nz ds (39e) Γ 3 β ( y y ) ( y y )( z z ) (39f) C C C Y = ny + nz ds + yc I Γ Z 4 3
30 β ( y y ) ( z z ) ( z z ) (39g) C C C Z = ny + nz ds + zc I Γ Y 3 4 S 3 3 P 4 4 ( y z) φs ( )( y z) βω = yn zn y z zn yn ds C Γ (39h) 4. Numerical examples On the basis of the analytical and numerical procedures presented in the previous sections, a computer program has been written and representative examples have been studied to demonstrate the validation, the efficiency, wherever possible the accuracy and the range of applications of the developed method. The numerical results have been obtained employing 1 nodal points (longitudinal discretization) and 400 boundary elements (cross section discretization). Example 1 In the first example, for comparison reasons, the forced vibration of a clamped 8 beam (Fig. a), ( E =,1 10 kn m, 7 G = 8, kn m, 3 ρ = 7,85tn m, l = 4m) of a hollow rectangular cross section, having geometric constants presented in Table 1, is examined. The beam is subjected to the suddenly applied uniformly distributed loading py ( t) = 50 kn m, ( ) p t = 500kN m at its centroid, as this is z shown in Fig. b. In Figs. 3, 4 the time histories of the midpoint displacements vl,t ( ), wl,t ( ) and angle of twist θ ( l,t ) x, respectively and in Table the maximum values of these kinematical components are presented as compared with those obtained from a BEM solution [19], noting the accuracy of the proposed
31 method. As it can be observed from these figures and table, the strong coupling between flexure and torsion may lead beams of doubly symmetric cross section undergoing biaxial transverse loading applied at their centroid to develop torsional rotation, a phenomenon that cannot be predicted from linear analysis. Example In this example, in order to investigate the response of a monosymmetric beam and the influence of the loading point upon the cross section, in nonlinear flexuraltorsional vibrations, the forced vibrations of a cantilever beam ( 8 E =, kn m, 7 G = 8, kn m, 3 ρ = 7,85tn m, l = 1m) of a thin-walled open shaped cross section (Fig. 5), under two load cases have been studied (its geometric constants are given in Table 3). More specifically, the beam is subjected to a suddenly applied concentrated force PY ( t) = 5kN either on the right (load case (i), Fig. 5b) or on the left (load case (ii), Fig. 5c) flange. In Figs. 6-9 the time histories of the axial displacement u( l, t), the transverse displacements ( t) wl, ( t) and the angle of twist ( ) vl,, θ x l,t of the cantilever beam, respectively and in Table 4 the maximum values of these kinematical components, are presented. From the aforementioned figures and table, it can easily be observed that geometrical nonlinearity affects substantially the dynamic response of the beam inducing non vanishing axial displacement and displacement with respect to z axis, while in linear analysis these kinematical components vanish. From the obtained results it can also be verified that the loading position has significant influence on the response, altering substantially the magnitude of the kinematical components. This discrepancy can be explained by the fact that in load case (ii), the change of eccentricity of the transverse
32 load during torsional rotation increases the magnitude of the twisting moment, acting adversely compared with load case (i), where the applied twisting moment is reduced during torsional rotation. Example 3 In order to demonstrate the range of applications of the developed method, in this final example the forced vibrations of a simply supported (free right end according to the axial boundary condition) steel L-shaped beam of unequal 8 (asymmetric cross section) legs (Fig. 10a), ( E =,1 10 kn m, 3 ρ = 7,85tn m, 7 G = 8, kn m, l = 1m), having the geometric constants presented in Table 5, is studied. The beam is subjected to a suddenly applied uniformly distributed twisting moment mx ( x,t) = 8kNm m (Fig. 10a,b). Due to lack of symmetry, apart from the angle of twist, the beam is expected to develop axial (axial displacement vanishes in linear analysis) and transverse displacements as well. In Figs the time histories of the axial displacement u( l,t ), the transverse displacements vl,t ( ), wl,t ( ) and the angle of twist θ ( l,t ) x, respectively, are presented, while in Table 6 the maximum values of these kinematical components, taking into account or ignoring rotary inertia effect are also shown. From these figures and table, it is observed that the geometrical nonlinearity leads to the increase of torsional stiffness decreasing the magnitude of angle of twist, while transverse displacements get significantly higher values compared with the linear ones. From Table 6 it is also noted that the influence of rotary inertia proves to be negligible on the magnitude of kinematical components.
33 Moreover, the response of a hinged-hinged beam (axially immovable ends), having the same cross section and length, under harmonic excitation is examined. More specifically, the beam is subjected to a uniformly distributed harmonic load Z ( ) = 0( ) ( ω f,lin ) p x,t p x sin t, as this is shown in Fig. 10a,c. The frequency ω f,lin is considered as ω f,lin = π f 1,lin, where f 1,lin = 118,00 Hz, is the first natural frequency of the examined beam, performing a linear analysis [34]. In Figs the time histories of the displacements v( l,t), wl,t ( ) system of axes and the angle of twist ( l,t ) x, with respect to the Syz θ are presented, noting the significant difference in response between linear and nonlinear analysis. More specifically, it is observed that only in the linear response deformations continue to increase with time, while the beating phenomenon observed in the nonlinear one is explained from the fact that large kinematical components increase the bar s fundamental natural frequency ω f (by increasing the stiffness of the bar due to the tensile axial force induced by the axially immovable ends), thereby causing a detuning of ω f with the frequency of the external loading ( ω f,lin ). After the kinematical components reach their maximum values, the amplitude of deformations decreases, leading to the reversal of the previously mentioned effects. 5. CONCLUDING REMARKS In this article a boundary element method is developed for the nonlinear flexural-torsional dynamic analysis of beams of arbitrary cross section, undergoing moderately large displacements and twisting rotations and small deformations, taking into account the effect of rotary inertia, warping inertia and change of eccentricity of transverse loads during torsional rotation. The beam is subjected to arbitrarily
34 distributed conservative transverse loads, which can be applied on any point of the cross section, and/or axial loads and twisting moments, while its edges are restrained by the most general boundary conditions. The main conclusions that can be drawn from this investigation are: a. The numerical technique presented in this investigation is well suited for computer aided analysis of beams of arbitrary simply or multiply connected cross section, supported by the most general boundary conditions and subjected to the combined action of arbitrarily distributed or concentrated time dependent loading. b. Accurate results are obtained using a relatively small number of nodal points along the beam. c. The geometrical nonlinearity leads to strong coupling between the axial, torsional and bending equilibrium equations resulting in a significantly different response of the beam compared to the one obtained by linear analysis. d. The strong coupling between flexure and torsion may lead beams of doubly symmetric cross sections, undergoing biaxial transverse loading applied on their centroid, to develop torsional rotation. e. The eccentricity change of the transverse loading during the torsional beam motion, resulting in additional torsional moment influences the beam response. f. The influence of rotary inertia, as shown in the treated examples, on the dynamic response of the beams, proves to be negligible on the magnitude of kinematical components. g. The developed procedure retains most of the advantages of a BEM solution over a FEM approach, although it requires longitudinal domain discretization. REFERENCES
35 [1]. Rozmarynowski B. and Szymczak C. (1984) Non-linear free torsional vibrations of thin-walled beams with bisymmetric cross-section, Journal of Sound and Vibration, 97, []. Crespo Da Silva M.R.M. (1988) Non-linear flexural-flexural-torsionalextensional dynamics of beams--i. Formulation, International Journal of Solids and Structures, 4, [3]. Crespo Da Silva M.R.M. (1988) Non-linear flexural-flexural-torsionalextensional dynamics of beams--ii. Response analysis, International Journal of Solids and Structures, 4, [4]. Pai P.F. and Nayfeh A.H. (1990) Three-dimensional nonlinear vibrations of composite beams--i. Equations of motion", Nonlinear Dynamics, 1, [5]. Pai P.F. and Nayfeh A.H. (1991) Three-dimensional nonlinear vibrations of composite beams--ii. flapwise excitations, Nonlinear Dynamics,, [6]. Pai P.F. and Nayfeh A.H. (1991) Three-dimensional nonlinear vibrations of composite beams--iii. Chordwise excitations, Nonlinear Dynamics,, [7]. Simo J.C. and Vu-Quoc L. (1991) A Geometrically-exact rod model incorporating shear and torsion-warping deformation, International Journal of Solids and Structures, 7, [8]. Qaisi M.I. (1993) Application of the Harmonic Balance Principle to the Nonlinear Free Vibration of Beams, Applied Acoustics, 40, [9]. Pai P.F. and Nayfeh A.H. (1994) A fully nonlinear theory of curved and twisted composite rotor blades accounting for warpings and three-dimensional stress effects, International Journal of Solids and Structures, 31,
36 [10]. Di Egidio A., Luongo A. and Vestroni F. (003) A non-linear model for the dynamics of open cross-section thin-walled beams--part I: formulation, International Journal of Non-Linear Mechanics 38 (003) [11]. Di Egidio A., Luongo A. and Vestroni F. (003) A non-linear model for the dynamics of open cross-section thin-walled beams--part II: forced motion, International Journal of Non-Linear Mechanics, 38, [1]. Mohri F., Azrar L. and Potier-Ferry M. (004) Vibration analysis of buckled thin-walled beams with open sections, Journal of Sound and Vibration, 75, [13]. Machado S.P. and Cortinez V.H. (007) Free vibration of thin-walled composite beams with static initial stresses and deformations, Engineering Structures, 9, [14]. Avramov K.V., Pierre C. and Shyriaieva N. (007) Flexural-flexural-torsional Nonlinear Vibrations of Pre-twisted Rotating Beams with Asymmetric crosssections, Journal of Vibration and Control, 13, [15]. Avramov K.V., Galas O.S. and Morachkovskii O.K. (009) Analysis of Flexural-Flexural-Torsional Nonlinear Vibrations of twisted rotating beams with cross-sectional deplanation, Strength of Materials, 41, No., [16]. Lopes Alonso R. and Ribeiro P. (008) Flexural and torsional non-linear free vibrations of beams using a p-version finite element, Computers and Structures, 86, [17]. Duan H. (008) Nonlinear free vibration analysis of asymmetric thin-walled circularly curved beams with open cross section, Thin-Walled Structures, 46,
37 [18]. Sapountzakis E.J. and Dourakopoulos J.A. (009) Nonlinear dynamic analysis of Timoshenko beams by BEM. Part I: Theory and numerical implementation, Nonlinear Dynamics, 58, [19]. Sapountzakis E.J. and Dourakopoulos J.A. (009) Nonlinear dynamic analysis of Timoshenko beams by BEM. Part II: Applications and validation, Nonlinear Dynamics, 58, [0]. Katsikadelis J.T. and Tsiatas G.C. (004) Non-linear dynamic analysis of beams with variable stiffness, Journal of Sound and Vibration, 70, [1]. Sapountzakis E.J. and Tsipiras V.J. (010) Nonlinear nonuniform torsional vibrations of bars by the boundary element method, Journal of Sound and Vibration, 39, []. Katsikadelis J.T. (00) The analog equation method. A boundary-only integral equation method for nonlinear static and dynamic problems in general bodies, Theoretical and Applied Mechanics, 7, [3]. Brenan K.E., Campbell S.L. and Petzold L.R. (1989) Numerical Solution of Initial-value Problems in Differential-Algebraic Equations, North-Holland, Amsterdam. [4]. Attard M.M. (1986) Nonlinear Theory of Non-Uniform Torsion of Thin- Walled Open Beams, Thin-Walled Structures, 4, [5]. Sapountzakis E.J. and Mokos V.G. (003) Warping shear stresses in nonuniform torsion by BEM, Computational Mechanics, 30, [6]. Ramm E. and Hofmann T.J. (1995) Stabtragwerke, Der Ingenieurbau, Ed.G. Mehlhorn, Band Baustatik/Baudynamik, Ernst & Sohn, Berlin. [7]. Rothert H. and Gensichen V. (1987) Nichtlineare Stabstatik, Springer Verlag, Berlin.
38 [8]. Vlasov V. (1963) Thin-walled elastic beams, Israel Program for Scientific Translations, Jerusalem. [9]. Armenakas A.E., (006) Advanced Mechanics of Materials and Applied Elasticity, Taylor & Francis Group, New York. [30]. Sapountzakis E.J. and Katsikadelis J.T. (000) Analysis of Plates Reinforced with Beams, Computational Mechanics, 6, [31]. Bazant Z.P. and Cedolin L. (1991) Stability of Structures: Elastic, Inelastic, Fracture and Damage Theories, Oxford University Press, New York. [3]. Sapountzakis E.J. (001) Nonuniform Torsion of Multi-Material Composite Bars by the Boundary Element Method, Computers and Structures, 79, [33]. Katsikadelis J.T. (00) Boundary Elements: Theory and Applications, Elsevier, Amsterdam-London. [34]. Sapountzakis E.J. and Dourakopoulos J.A. (009) Shear deformation effect in flexural-torsional vibrations of beams by BEM, Acta Mechanica, 03, APPENDIX A Functions F ( i 1,,...,7 ) vi =, wi F ( i 1,,...,7 ) = and i F θ ( i 1,,...,11) = of the differential equations (13), concerning the nonlinear dynamic analysis of beams of arbitrary cross section are given as: F = v y θ sinθ z θ cosθ (A.1a) v1 C x x C x x ( θ θ ) F = A y sin + z cos (A.1b) v C x C x ( θ θ ) F = A z sin y cos (A.1c) v3 C x C x ( ) ( ) Fv4 = IY IZ sinθx cosθxv IY IZ cos θx I Y w (A.1d)
39 ( ) ( ) Fv5 = IY IZ sinθx cosθxw + IY IZ cos θx I Y v (A.1e) ( ) F = I I cosθ sinθ (A.1f) v6 Y Z x x ( ) F = I I cos θ I (A.1g) v7 Y Z x Y F = w + y θ cosθ z θ sinθ (A.a) w1 C x x C x x ( θ θ ) F = A y cos + z sin (A.b) w C x C x ( θ θ ) F = A z cos + y sin (A.c) w3 C x C x ( ) ( ) Fw4 = IZ IY cosθx sinθxw + IZ IY cos θx I Z v ( ) ( ) Fw5 = IZ IY cosθx sinθxv IZ IY cos θx I Z w (A.d) (A.e) ( ) F = I I cosθ sinθ (A.f) w6 Z Y x x ( ) F = I I cos θ I (A.g) w7 Z Y x Z Fθ 1 = ycv θ x cosθx ycw θ x sinθx zcw θ x cosθx + zcv θ x sinθx (A.3a) ( ) Fθ = z v cosθ + w sinθ (A.3b) C x x ( ) Fθ = y v sinθ w cosθ (A.3c) 3 C x x ( ) Fθ = w sinθ + v cosθ (A.3d) 4 x x ( ) Fθ = w cosθ v sinθ (A.3e) 5 x x ( ) Fθ = A z cosθ + y sinθ (A.3f) 6 C x C x ( ) Fθ = A z sinθ + y cosθ (A.3g) 7 C x C x ( ) ( ) Fθ 8 = IY IZ cos θx + I Z v IZ IY vw sinθxcosθx + Is +
40 ( I I ) cos θ I w + Z Y x + Y ( ) ( ) Fθ 9 = IZ IY cos θx + I Y w + IY IZ sinθxcosθxv ( ) ( ) Fθ 10 = IZ IY cos θx I Z v + IZ IY sinθxcosθxw (A.3h) (A.3i) (A.3j) ( ) ( ) 11 Y Z x x Z Y x x Fθ = I I sinθ cosθ v + I I sinθ cosθ w + ( ) + IZ IY cos θx IZ + I Y vw (A.3k) Functions G u1, vi G ( i 1,,...,7 ) =, wi G ( i 1,,...,7 ) = and i G θ ( i 1,,...,11) = of the governing differential equations (15), concerning the nonlinear dynamic analysis of beams of arbitrary cross section are given as: ( θθ θθ θ θ θ ) G = z w + w + v + v + w + u1 C x x x x x x x ( θθ θθ θ θ θ ) + y v + v w w + v (A.4) C x x x x x x x ( ) G = z θ y θ + θ θ + v (A.5a) v1 C x C x x x G = w + w + w v 4v v v (A.5b) v θx θx θx θx θxθx θxθx θx G = (A.5c) v3 θx θx θx G = (A.5d) v4 θx θx θx θx θx θx θx v5 ( Z Y )( x x ) 1 G = I I θ v θ + v IZw A ycθx + zc zcθx θx + ( IZ IY ) ( ) θθ x xv + θxv + θ xθ x + θθ x x v + θθ x xv θ x w θxw + v θθ x x +
41 + IZ w θ x v θ xw + v θ x + v θ xθ x θ x w (A.5e) ( ) G = A y z θ θ (A.5f) v6 C C x x G = v y θ θ z θ (A.5g) v7 C x x C x ( ) G = w + y θ z θ + θ θ (A.6a) w1 C x C x x x G = v + v + v + w + 4w + w + w (A.6b) w θx θx θx θx θxθx θxθx θx G = 5 (A.6c) w3 θx θx θx θx θx θx θx G = (A.6d) w4 θx θx θx w5 ( Z Y )( x x ) 1 G = I I θ w θ w + IYv + A zcθx + yc ycθ x θx ( IZ IY ) θθ x x w + θxw + ( θ xθ x + θθ x x ) w + θθ x xw + θ xv + θxv + w θθ x x IY v θ x + w θ xv w θ x w θ xθ x θ x v (A.6e) ( ) G = A z + y θ θ (A.6f) w6 C C x x G = w + y θ z θ θ (A.6g) w7 C x C x x IS Gθ 1 = θ x + yc ( w v θx) zc ( v + w θx) (A.7a) A x x Gθ = v w v θ + w θ (A.7b) Gθ 3 = θ x w θx + w θ x + θ x v + θ x v + θ x w θx (A.7c) Gθ 4 = θ x v θx + θ x w + θ x w v θ x θ x v θx (A.7d) 5 Y Z S Gθ = v I + w I + I (A.7e)
42 1 G 6 A θ = zc + ycθ x zcθx (A.7f) 1 G 7 A θ = yc zcθx ycθx (A.7g) ( ) θ ( ) ( θ ) Υ ( θ ) G = I I ww vv + I wv + ww + I vw + vv (A.7h) θ8 Z Y x Z x x Gθ 9 = θxzp + θ y x yp θ y xypy 6 (A.7i) Gθ 10 = θxzp θ z x yp z z p θ z x 6 (A.7j) IS Gθ 11 = θ x ycv θx zcw θx zcv + ycw A (A.7k)
43 (a) (b) Fig. 1. Prismatic beam in axial - flexural - torsional loading (a) of an arbitrary cross-section occupying the two dimensional region Ω (b).
Shear Deformation Effect in Flexural-torsional Vibrations of Composite Beams by Boundary Element Method (BEM)
Shear Deformation Effect in Flexural-torsional Vibrations of Composite Beams by Boundary Element Method (BEM) E. J. SAPOUNTZAKIS J. A. DOURAKOPOULOS School of Civil Engineering, National Technical University,
More informationCalculation of Shear Areas and Torsional Constant using the Boundary Element Method with Scada Pro Software
NATIONAL TECHNICAL UNIVERSITY OF ATHENS School of Civil Engineering Institute of Structural Analysis & Antiseismic Research Dr. Sapountakis J. Evangelos Dr. Civil Engineer NTUA Professor NTUA Calculation
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 informationConsider an elastic spring as shown in the Fig.2.4. When the spring is slowly
.3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original
More informationDynamic and buckling analysis of FRP portal frames using a locking-free finite element
Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.
More informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
An integral equation formulation of the coupled vibrations of uniform Timoshenko beams Masa. Tanaka & A. N. Bercin Department of Mechanical Systems Engineering, Shinshu University 500 Wakasato, Nagano
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 informationFinite element computation of leaky modes in straight and helical elastic waveguides
Finite element computation of leaky modes in straight and helical elastic waveguides Khac-Long Nguyen, Fabien Treyssede, Christophe Hazard, Anne-Sophie Bonnet-Ben Dhia To cite this version: Khac-Long Nguyen,
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 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 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 informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More information[8] Bending and Shear Loading of Beams
[8] Bending and Shear Loading of Beams Page 1 of 28 [8] Bending and Shear Loading of Beams [8.1] Bending of Beams (will not be covered in class) [8.2] Bending Strain and Stress [8.3] Shear in Straight
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
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 information3. BEAMS: STRAIN, STRESS, DEFLECTIONS
3. BEAMS: STRAIN, STRESS, DEFLECTIONS The beam, or flexural member, is frequently encountered in structures and machines, and its elementary stress analysis constitutes one of the more interesting facets
More informationNonlocal computational methods applied to composites structures
Nonlocal computational methods applied to composites structures Norbert Germain, Frédéric Feyel, Jacques Besson To cite this version: Norbert Germain, Frédéric Feyel, Jacques Besson. Nonlocal computational
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationComb Resonator Design (2)
Lecture 6: Comb Resonator Design () -Intro. to Mechanics of Materials Sh School of felectrical ti lengineering i and dcomputer Science, Si Seoul National University Nano/Micro Systems & Controls Laboratory
More informationOn measurement of mechanical properties of sound absorbing materials
On measurement of mechanical properties of sound absorbing materials Nicolas Dauchez, Manuel Etchessahar, Sohbi Sahraoui To cite this version: Nicolas Dauchez, Manuel Etchessahar, Sohbi Sahraoui. On measurement
More informationChapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14
Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.
More informationJeff Brown Hope College, Department of Engineering, 27 Graves Pl., Holland, Michigan, USA UNESCO EOLSS
MECHANICS OF MATERIALS Jeff Brown Hope College, Department of Engineering, 27 Graves Pl., Holland, Michigan, USA Keywords: Solid mechanics, stress, strain, yield strength Contents 1. Introduction 2. Stress
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationUnfolding the Skorohod reflection of a semimartingale
Unfolding the Skorohod reflection of a semimartingale Vilmos Prokaj To cite this version: Vilmos Prokaj. Unfolding the Skorohod reflection of a semimartingale. Statistics and Probability Letters, Elsevier,
More informationFIXED BEAMS IN BENDING
FIXED BEAMS IN BENDING INTRODUCTION Fixed or built-in beams are commonly used in building construction because they possess high rigidity in comparison to simply supported beams. When a simply supported
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More information2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at A and supported at B by rod (1). What is the axial force in rod (1)?
IDE 110 S08 Test 1 Name: 1. Determine the internal axial forces in segments (1), (2) and (3). (a) N 1 = kn (b) N 2 = kn (c) N 3 = kn 2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at
More informationELASTIC WAVE PROPAGATION IN THREE-DIMENSIONAL PERIODIC COMPOSITE MATERIALS
ELASTIC WAVE PROPAGATION IN THREE-DIMENSIONAL PERIODIC COMPOSITE MATERIALS B. Auld, Y. Shui, Y. Wang To cite this version: B. Auld, Y. Shui, Y. Wang. ELASTIC WAVE PROPAGATION IN THREE-DIMENSIONAL PERI-
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationPresented By: EAS 6939 Aerospace Structural Composites
A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have
More informationThe CR Formulation: BE Plane Beam
6 The CR Formulation: BE Plane Beam 6 Chapter 6: THE CR FORMUATION: BE PANE BEAM TABE OF CONTENTS Page 6. Introduction..................... 6 4 6.2 CR Beam Kinematics................. 6 4 6.2. Coordinate
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
More informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
More informationNONLINEAR VIBRATIONS OF ROTATING 3D TAPERED BEAMS WITH ARBITRARY CROSS SECTIONS
COMPDYN 2013 4 th ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, V. Papadopoulos, V. Plevris (eds.) Kos Island, Greece, 12 14 June
More informationRHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING
RHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING Jean-François Semblat To cite this version: Jean-François Semblat. RHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING. Journal of Sound and Vibration, Elsevier,
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationUnit 18 Other Issues In Buckling/Structural Instability
Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
More informationLecture Pure Twist
Lecture 4-2003 Pure Twist pure twist around center of rotation D => neither axial (σ) nor bending forces (Mx, My) act on section; as previously, D is fixed, but (for now) arbitrary point. as before: a)
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationSome approaches to modeling of the effective properties for thermoelastic composites
Some approaches to modeling of the ective properties for thermoelastic composites Anna Nasedkina Andrey Nasedkin Vladimir Remizov To cite this version: Anna Nasedkina Andrey Nasedkin Vladimir Remizov.
More informationSub. Code:
Important Instructions to examiners: ) The answers should be examined by key words and not as word-to-word as given in the model answer scheme. ) The model answer and the answer written by candidate may
More informationMechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection
Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts
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 informationDISTORTION ANALYSIS OF TILL -WALLED BOX GIRDERS
Nigerian Journal of Technology, Vol. 25, No. 2, September 2006 Osadebe and Mbajiogu 36 DISTORTION ANALYSIS OF TILL -WALLED BOX GIRDERS N. N. OSADEBE, M. Sc., Ph. D., MNSE Department of Civil Engineering
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 informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationUnit - 7 Vibration of Continuous System
Unit - 7 Vibration of Continuous System Dr. T. Jagadish. Professor for Post Graduation, Department of Mechanical Engineering, Bangalore Institute of Technology, Bangalore Continuous systems are tore which
More informationChapter 4-b Axially Loaded Members
CIVL 222 STRENGTH OF MATERIALS Chapter 4-b Axially Loaded Members AXIAL LOADED MEMBERS Today s Objectives: Students will be able to: a) Determine the elastic deformation of axially loaded member b) Apply
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 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 informationThe sound power output of a monopole source in a cylindrical pipe containing area discontinuities
The sound power output of a monopole source in a cylindrical pipe containing area discontinuities Wenbo Duan, Ray Kirby To cite this version: Wenbo Duan, Ray Kirby. The sound power output of a monopole
More informationCHAPTER -6- BENDING Part -1-
Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and
More informationAnalytical Strip Method for Thin Isotropic Cylindrical Shells
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
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 informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
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 informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationChapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements
CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
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 informationMODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM
MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,
More informationA Boundary Integral Formulation for the Dynamic Behavior of a Timoshenko Beam
A Boundary Integral Formulation for the Dynamic Behavior of a Timoshenko Beam M. Schanz and H. Antes Institute of Applied Mechanics Technical University Braunschweig, Germany e-mail: m.schanz@tu-bs.de,
More informationWaves vibrations. Maxime Nicolas. To cite this version: HAL Id: cel Maxime Nicolas.
Waves vibrations Maxime Nicolas To cite this version: Maxime Nicolas. Waves vibrations. Engineering school. France. 2016. HAL Id: cel-01440543 https://hal.archives-ouvertes.fr/cel-01440543
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationCase report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122,
Case report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122, 244902 Juan Olives, Zoubida Hammadi, Roger Morin, Laurent Lapena To cite this version: Juan Olives,
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationVerification Examples. FEM-Design. version
FEM-Design 6.0 FEM-Design version. 06 FEM-Design 6.0 StruSoft AB Visit the StruSoft website for company and FEM-Design information at www.strusoft.com Copyright 06 by StruSoft, all rights reserved. Trademarks
More informationSimulation of Geometrical Cross-Section for Practical Purposes
Simulation of Geometrical Cross-Section for Practical Purposes Bhasker R.S. 1, Prasad R. K. 2, Kumar V. 3, Prasad P. 4 123 Department of Mechanical Engineering, R.D. Engineering College, Ghaziabad, UP,
More informationCHAPTER 4: BENDING OF BEAMS
(74) CHAPTER 4: BENDING OF BEAMS This chapter will be devoted to the analysis of prismatic members subjected to equal and opposite couples M and M' acting in the same longitudinal plane. Such members are
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More informationFlexural-Torsional Buckling of General Cold-Formed Steel Columns with Unequal Unbraced Lengths
Proceedings of the Annual Stability Conference Structural Stability Research Council San Antonio, Texas, March 21-24, 2017 Flexural-Torsional Buckling of General Cold-Formed Steel Columns with Unequal
More informationSlender Structures Load carrying principles
Slender Structures Load carrying principles Basic cases: Extension, Shear, Torsion, Cable Bending (Euler) v017-1 Hans Welleman 1 Content (preliminary schedule) Basic cases Extension, shear, torsion, cable
More informationPURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.
BENDING STRESS The effect of a bending moment applied to a cross-section of a beam is to induce a state of stress across that section. These stresses are known as bending stresses and they act normally
More informationToward a novel approach for damage identification and health monitoring of bridge structures
Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationDesign of Beams (Unit - 8)
Design of Beams (Unit - 8) Contents Introduction Beam types Lateral stability of beams Factors affecting lateral stability Behaviour of simple and built - up beams in bending (Without vertical stiffeners)
More informationChapter 5 Elastic Strain, Deflection, and Stability 1. Elastic Stress-Strain Relationship
Chapter 5 Elastic Strain, Deflection, and Stability Elastic Stress-Strain Relationship A stress in the x-direction causes a strain in the x-direction by σ x also causes a strain in the y-direction & z-direction
More informationThick-Walled Open Section Beam and Finite Nodal-Line Method
Solids and Structures (SAS) Volume 5, 206 doi: 0.4355/sas.206.05.00 www.seipub.org/sas/ hick-walled Open Section Beam and Finite Nodal-Line Method Yaoqing Gong *, Sai ao 2 School of Civil Engineering,
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
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 informationENG2000 Chapter 7 Beams. ENG2000: R.I. Hornsey Beam: 1
ENG2000 Chapter 7 Beams ENG2000: R.I. Hornsey Beam: 1 Overview In this chapter, we consider the stresses and moments present in loaded beams shear stress and bending moment diagrams We will also look at
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 informationCorrection of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams
Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Mohamed Shaat* Engineering and Manufacturing Technologies Department, DACC, New Mexico State University,
More informationLecture 7: The Beam Element Equations.
4.1 Beam Stiffness. A Beam: A long slender structural component generally subjected to transverse loading that produces significant bending effects as opposed to twisting or axial effects. MECH 40: Finite
More informationFlexure of Thick Cantilever Beam using Third Order Shear Deformation Theory
International Journal of Engineering Research and Development e-issn: 78-67X, p-issn: 78-8X, www.ijerd.com Volume 6, Issue 1 (April 13), PP. 9-14 Fleure of Thick Cantilever Beam using Third Order Shear
More information2018. Nonlinear free vibration analysis of nanobeams under magnetic field based on nonlocal elasticity theory
2018. Nonlinear free vibration analysis of nanobeams under magnetic field based on nonlocal elasticity theory Tai-Ping Chang National Kaohsiung First University of Science and Technology, Kaohsiung, Taiwan
More informationY. M. TEMIS AND V. V. KARABAN 1
J. KSIAM Vol.5, No.2, 39-51, 2001 BOUNDARY ELEMENT TECHNIQUE IN TORSION PROBLEMS OF BEAMS WITH MULTIPLY CONNECTED CROSS-SECTIONS Y. M. TEMIS AND V. V. KARABAN 1 Abstract. This paper shows how boundary
More informationMechanics of Structure
S.Y. Diploma : Sem. III [CE/CS/CR/CV] Mechanics of Structure Time: Hrs.] Prelim Question Paper Solution [Marks : 70 Q.1(a) Attempt any SIX of the following. [1] Q.1(a) Define moment of Inertia. State MI
More informationThermodynamic form of the equation of motion for perfect fluids of grade n
Thermodynamic form of the equation of motion for perfect fluids of grade n Henri Gouin To cite this version: Henri Gouin. Thermodynamic form of the equation of motion for perfect fluids of grade n. Comptes
More information1 Static Plastic Behaviour of Beams
1 Static Plastic Behaviour of Beams 1.1 Introduction Many ductile materials which are used in engineering practice have a considerable reserve capacity beyond the initial yield condition. The uniaxial
More informationChapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING )
Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) 5.1 DEFINITION A construction member is subjected to centric (axial) tension or compression if in any cross section the single distinct stress
More informationSolution to Sylvester equation associated to linear descriptor systems
Solution to Sylvester equation associated to linear descriptor systems Mohamed Darouach To cite this version: Mohamed Darouach. Solution to Sylvester equation associated to linear descriptor systems. Systems
More informationMechanics of the Fouetté turn
Mechanics of the Fouetté turn Akiko Imura, M.R. Yeadon To cite this version: Akiko Imura, M.R. Yeadon. Mechanics of the Fouetté turn. Human Movement Science, Elsevier, 2010, 29 (6), pp.947. .
More informationMechanics of Materials Primer
Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus
More information