Geometric nonlinear dynamic analysis of curved beams using curved beam element

Instead of using the previous straight beam element to approximate the curved beam,in this paper,a curvilinear coordinate is employed to describe the deformations,and a new curved beam element is proposed to model the curved beam.Based on exact nonlinear strain-displacement relation,virtual work principle is used to derive dynamic equations for a rotating curved beam,with the effects of axial extensibility,shear deformation and rotary inertia taken into account.The constant matrices are solved numerically utilizing the Gauss quadrature integration method.Newmark and Newton-Raphson iteration methods are adopted to solve the differential equations of the rigid-flexible coupling system.The present results are compared with those obtained by commercial programs to validate the present finite method.In order to further illustrate the convergence and efficiency characteristics of the present modeling and computation formulation,comparison of the results of the present formulation with those of the ADAMS software are made.Furthermore,the present results obtained from linear formulation are compared with those from nonlinear formulation,and the special dynamic characteristics of the curved beam are concluded by comparison with those of the straight beam.     

Free vibration and spatial stability of non-symmetric thin-walled curved beams with variable curvatures

An improved formulation for free vibration and spatial stability of non-symmetric thin-walled curved beams is presented based on the displacement field considering variable curvature effects and the second-order terms of finite-semitangential rotations. By introducing Vlasov’s assumptions and integrating over the non-symmetric cross-section, the total potential energy is consistently derived from the principle of virtual work for a continuum. In this formulation, all displacement parameters and the warping function are defined at the centroid axis and also thickness-curvature effects and Wagner effect are accurately taken into account. For F.E. analysis, a thin-walled curved beam element is developed using the third-order Hermitian polynomials. In order to illustrate the accuracy and the practical usefulness of the present method, numerical solutions by this study are presented with the results analyzed by ABAQUS’ shell elements. Particularly, the effect of arch rise to span length ratio is investigated on vibrational and buckling behaviour of non-symmetric curved beams.     

Modified shooting approach to the non-linear periodic forced response of isotropic/composite curved beams

Large amplitude periodic forced vibration of curved beams under periodic excitation is investigated using a three-noded beam element. The element is based on the higher-order shear deformation theory satisfying interlayer continuity of displacements and transverse shear stress, and top-bottom conditions on the latter. The periodic responses are obtained using shooting technique coupled with Newmark time marching and arc length continuation algorithm developed. The second order governing differential equations of motion are solved without transforming to the first order differential equations thereby resulting in a computationally more efficient algorithm. The effects of excitation amplitude, support conditions and beam curvature on the frequency versus response amplitude relation are highlighted. The typical frequency response curves for isotropic and cross-ply laminated curved beams are presented. Phenomenon of strong modal interactions is observed.     

Non-linearities in rotation and thickness deformation in a new third-order thickness deformation theory for static and dynamic analysis of isotropic and laminated doubly curved shells

A geometrically non-linear theory is developed for shells of generic shape allowing for third-order thickness and shear deformation and rotary inertia by using eight parameters; geometric imperfections are also taken into account. The geometrically non-linear strain–displacement relationships are derived retaining full non-linear terms in all the 8 parameters, i.e. in-plane and transverse displacements, rotations of the normal and thickness deformation parameters; these relationships are presented in curvilinear coordinates, ready to be implemented in computer codes. Higher order terms in the transverse coordinate are retained in the derivation so that the theory is suitable also for thick laminated shells. Three-dimensional constitutive equations are used for linear elasticity. The theory is applied to circular cylindrical shells complete around the circumference and simply supported at both ends to study initially static finite deformation. Both radially distributed forces and displacement-dependent pressure are used as load and results for different shell theories are compared. Results show that a 6 parameter non-linear shell theory is quite accurate for isotropic shells. Finally, large-amplitude forced vibrations under harmonic excitation are investigated by using the new theory and results are compared to other available theories. The new theory with non-linearity in all the 8 parameters is the only one to predict correctly the thickness deformation; it works accurately for both static and dynamics loads.     

Explicit solutions for shearing and radial stresses in curved beams

In this paper the formulae for the shearing and radial stresses in curved beams are derived analytically based on the solution for a Volterra integral equation of the second kind. These formulae satisfy both the equilibrium equations and the static boundary conditions on the surfaces of the beams. As some applications, the resulting solutions are used to calculate the shearing and radial stresses in a cantilevered curved beam subjected to a concentrated force at its free end. The numerical results are compared with other existing approximate solutions as well as the corresponding solutions based on the theory of elasticity. The calculations show a better agreement between the present solution and the one based on the theory of elasticity. The resulting formulae can be applied to more general cases of curved beams with arbitrary shapes of cross-sections.     

Approximate vibration analysis of laminated curved panel using higher-order shear deformation theory

An approximate analysis for free vibration of a laminated curved panel (shell) with four edges simply supported (SS2), is presented in this paper. The transverse shear deformation is considered by using a higher-order shear deformation theory. For solving the highly coupled partial differential governing equations and associated boundary conditions, a set of solution functions in the form of double trigonometric Fourier series, which are required to satisfy the geometry part of the considered boundary conditions, is assumed in advance. By applying the Galerkin procedure both to the governing equations and to the natural boundary conditions not satisfied by the assumed solution functions, an approximate solution, capable of providing a reliable prediction for the global response of the panel, is obtained. Numerical results of antisymmetric angle-ply as well as symmetric cross-ply and angle-ply laminated curved panels are presented and discussed.     

Out-of-plane responses of a circular curved Timoshenko beam due to a moving load

For the cases of using the finite curved beam elements and taking the effects of both the shear deformation and rotary inertias into consideration, the literature regarding either free or forced vibration analysis of the curved beams is rare. Thus, this paper tries to determine the dynamic responses of a circular curved Timoshenko beam due to a moving load using the curved beam elements. By taking account of the effect of shear deformation and that of rotary inertias due to bending and torsional vibrations, the stiffness matrix and the mass matrix of the curved beam element were obtained from the force–displacement relations and the kinetic energy equations, respectively. Since all the element property matrices for the curved beam element are derived based on the local polar coordinate system (rather than the local Cartesian one), their coefficients are invariant for any curved beam element with constant radius of curvature and subtended angle and one does not need to transform the property matrices of each curved beam element from the local coordinate system to the global one to achieve the overall property matrices for the entire curved beam structure before they are assembled. The availability of the presented approach has been verified by both the existing analytical solutions for the entire continuum curved beam and the numerical solutions for the entire discretized curved beam composed of the conventional straight beam elements based on either the consistent-mass model or the lumped-mass model. In addition to the typical circular curved beams, a hybrid curved beam composed of one curved-beam segment and two identical straight-beam segments subjected to a moving load was also studied. Influence on the dynamic responses of the curved beams of the slenderness ratio, moving-load speed, shear deformation and rotary inertias was investigated.     

Free vibration of horizontally curved thin-walled beams with rectangular hollow sections considering two compatible displacement fields

A finite element is presented for vibration analyses of horizontally curved thin-walled rectangular hollow beams. Eight cross-section deformation modes are employed to describe the mid-surface contour displacement field with the modal superposition method. Focused on the in-plane moment equilibrium condition and the displacement continuity condition, two compatible displacement fields are constructed to calculate the strain energy and the kinetic energy of the beam, respectively. With the application of Hamilton’s principle the dynamic governing equations are formulated, and then approximated for the finite element implementation. Finally, numerical examples are illustrated to verify the validity of the present theory.     

A new higher-order shear deformation theory and refined beam element of composite laminates

A new higher-order shear deformation theory based on global-local superposition technique is developed. The theory satisfies the free surface conditions and the geometric and stress continuity conditions at interfaces. The global displacement components are of the Reddy theory and local components are of the internal first to third-order terms in each layer. A two-node beam element based on this theory is proposed. The solutions are compared with 3D-elasticity solutions. Numerical results show that present beam element has higher computational efficiency and higher accuracy.The project supported by the National Natural Science Foundation of China (10172023)     

Geometrically nonlinear isogeometric analysis of laminated composite plates based on higher-order shear deformation theory

In this paper, we present an effectively numerical approach based on isogeometric analysis (IGA) and higher-order shear deformation theory (HSDT) for geometrically nonlinear analysis of laminated composite plates. The HSDT allows us to approximate displacement field that ensures by itself the realistic shear strain energy part without shear correction factors (SCFs). IGA utilizing basis functions namely B-splines or non-uniform rational B-splines (NURBS) enables to satisfy easily the stringent continuity requirement of the HSDT model without any additional variables. The nonlinearity of the plates is formed in the total Lagrange approach based on the small strain assumptions. Numerous numerical validations for the isotropic, orthotropic, cross-ply and angle-ply laminated plates are provided to demonstrate the effectiveness of the proposed method.     

Vibration analysis of a stepped laminated composite Timoshenko beam

The goal of this study is to investigate the vibration characteristics of a stepped laminated composite Timoshenko beam. Based on the first order shear deformation theory, flexural rigidity and transverse shearing rigidity of a laminated beam are determined. In order to account for the effect of shear deformation and rotary inertia of the stepped beam, Timoshenko beam theory is then used to deduce the frequency function. Graphs of the natural frequencies and mode shapes of a T300/970 laminated stepped beam are given, in order to illustrate the influence of step location parameter exerts on the dynamic behavior of the beam.     

柔性梁大变形计算的一种新方法

研究了柔性梁大变形问题。常规Lagrangian有限元格式在处理大变形问题时,由于其单元插值函数不满足位移场的协调性要求,从而需要划分较多的单元,才能得到较好的结果。本文首先推导了Lagrangian坐标描述下的位移场变量满足的协调关系式,利用此关系式给出了位移场协调的非线性单元插值函数。基于虚功原理导出了梁大变形问题的非线性控制方程,数值计算结果证明了本文方法的正确性和有效性。     

Micromorphic continuum. Part II: Finite deformation plasticity coupled with damage

It is demonstrated how a micromorphic plasticity theory may be formulated on the basis of multiplicative decompositions of the macro- and microdeformation gradient tensor, respectively. The theory exhibits non-linear isotropic and non-linear kinematic hardening. The yield function is expressed in terms of Mandel stress and double stress tensors, appropriately defined for micromorphic continua. Flow rules are derived from the postulate of Il’iushin and represent generalized normality conditions. Evolution equations for isotropic and kinematic hardening are introduced as sufficient conditions for the validity of the second law of thermodynamics in every admissible process. Finally, it is sketched how isotropic damage effects may be incorporated in the theory. This is done for the concept of effective stress combined with the hypothesis of strain equivalence.     

A new non-linear higher-order shear deformation theory for large-amplitude vibrations of laminated doubly curved shells

A consistent higher-order shear deformation non-linear theory is developed for shells of generic shape, taking geometric imperfections into account. The geometrically non-linear strain-displacement relationships are derived retaining full non-linear terms in the in-plane displacements; they are presented in curvilinear coordinates in a formulation ready to be implemented. Then, large-amplitude forced vibrations of a simply supported, laminated circular cylindrical shell are studied (i) by using the developed theory, and (ii) keeping only non-linear terms of the von Kármán type. Results show that inaccurate results are obtained by keeping only non-linear terms of the von Kármán type for vibration amplitudes of about two times the shell thickness for the studied case.     

钢管混凝土拱稳定分析的三维退化层合曲梁单元

为计算钢管混凝土拱的屈曲荷载,本文在文[1]三维退化梁单元的基础上,采用等效数值积分法,构造,出120－20结点三维退化层合曲梁单元,并考虑几何非线性影响,给出用于层合梁或拱线弹性稳定性分析的有限元列式,最后,以绍兴轻纺大桥为工程背景,计算出轻纺大枯钢管混凝土拱面内及面外屈曲的稳定系数。     

Mode III fracture behavior of laminated composite with edge crack in torsion

A three-dimensional (3-D) finite element analysis was performed on a [90,(+45/−45)_n,(−45/+45)_n,90]_s class of laminated composites under the edge crack torsion (ECT) test configuration. Finite element delamination models were established and formulas for calculating the Mode III fracture toughness from 3-D finite element models were developed. The relations between the interlaminar fracture behavior and various configuration parameters were investigated and the effects of point loads, ends, geometry, Mode II interference, and friction were evaluated. Results showed that with proper selection of ECT specimen configuration and layup, the delamination could grow in pure Mode III in the middle region of the specimen. Specimen end effect played an important role in the ECT test. A Mode II component occurred in the end regions but it did not interfere significantly with the Mode III delamination state. Specimen dimension ratio, layup, and crack length exhibited significant effect on the interlaminar fracture behavior and the calculated strain energy release rates. However, friction between crackfaces was found to have negligible effect on the interlaminar properties.     

Isogeometric analysis for the dynamic problem of curved structures including warping effects

Curved beam structures have been used in many civil, mechanical, aircraft, and aerospace constructions. The analysis is mainly based on solid and plate models due to the fact that traditional curved beam elements do not include nonuniform warping effects, especially in the dynamic analysis. In this article, independent warping parameters have been taken into account and the initial curvature effect is considered. Curved beam’s behavior becomes more complex, even for dead loading, due to the coupling between axial force, bending moments, and torque that curvature produces. In addition to these, the Isogeometric tools (b-splines or NURBS), either integrated in the Finite Element Method or in a Boundary Element–based Method called Analog Equation Method, have been employed in this contribution for the dynamic analysis of horizontally curved beams of open or closed (box-shaped) cross sections. Free vibration characteristics and responses of the stress resultants and displacements to moving loading have been studied.     

压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

Natural frequencies of composite beams with a variable fiber volume fraction including rotary inertia and shear deformation

The present study is concerned with the vibration analysis of symmetric composite beams with a variable fiber volume fraction through thickness. First-order shear deformation and rotary inertia have been included in the analysis. The solution procedure is applicable to arbitrary boundary conditions. Continuous gradation of the fiber volume fraction is modeled in the form of an m-th power polynomial of the coordinate axis in the thickness direction of the beam. By varying the fiber volume fraction within the symmetric composite beam to create a functionally graded material （FGM）, certain vibration characteristics are affected. Results are presented to demonstrate the effects of shear deformation, fiber volume fraction, and boundary conditions on the natural frequencies and mode shapes of composite beams.