Nonlinear dynamic analysis of Timoshenko beams by BEM. Part I: Theory and numerical implementation

In this two-part contribution, a boundary element method is developed for the nonlinear dynamic analysis of beams of arbitrary doubly symmetric simply or multiply connected constant cross section, undergoing moderate large displacements and small deformations under general boundary conditions, taking into account the effects of shear deformation and rotary inertia. Part I is devoted to the theoretical developments and their numerical implementation and Part II discusses analytical and numerical results obtained from both analytical or numerical research efforts from the literature and the proposed method. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading and bending moments in both directions as well as to axial loading. To account for shear deformations, the concept of shear deformation coefficients is used. Five boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to two stress functions and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique yields a nonlinear coupled system of equations of motion. The solution of this system is accomplished iteratively by employing the average acceleration method in combination with the modified Newton–Raphson method. The evaluation of the shear deformation coefficients is accomplished from the aforementioned stress functions using only boundary integration. The proposed model takes into account the coupling effects of bending and shear deformations along the member, as well as the shear forces along the span induced by the applied axial loading.

     

Nonlinear dynamic analysis of Timoshenko beams by BEM. Part II: applications and validation

In this two-part contribution, a boundary element method is developed for the nonlinear dynamic analysis of beams of arbitrary doubly symmetric simply or multiply connected constant cross section, undergoing moderate large displacements and small deformations under general boundary conditions, taking into account the effects of shear deformation and rotary inertia. In Part I the governing equations of the aforementioned problem have been derived, leading to the formulation of five boundary value problems with respect to the transverse displacements, to the axial displacement and to two stress functions. These problems are numerically solved using the Analog Equation Method, a BEM based method. In this Part II, numerical examples are worked out to illustrate the efficiency, the accuracy and the range of applications of the developed method. Thus, the results obtained from the proposed method are presented as compared with those from both analytical and numerical research efforts from the literature. More specifically, the shear deformation effect in nonlinear free vibration analysis, the influence of geometric nonlinearities in forced vibration analysis, the shear deformation effect in nonlinear forced vibration analysis, the nonlinear dynamic analysis of Timoshenko beams subjected to arbitrary axial and transverse in both directions loading, the free vibration analysis of Timoshenko beams with very flexible boundary conditions and the stability under axial loading (Mathieu problem) are presented and discussed through examples of practical interest.

     

Free vibration analysis of third-order shear deformable composite beams using dynamic stiffness method

The dynamic stiffness method is introduced to investigate the free vibration of laminated composite beams based on a third-order shear deformation theory which accounts for parabolic distribution of the transverse shear strain through the thickness of the beam. The exact dynamic stiffness matrix is found directly from the analytical solutions of the basic governing differential equations of motion. The Poisson effect, shear deformation, rotary inertia, in-plane deformation are considered in the analysis. Application of the derived dynamic stiffness matrix to several particular laminated beams is discussed. The influences of Poisson effect, material anisotropy, slenderness and end condition on the natural frequencies of the beams are investigated. The numerical results are compared with the existing solutions in literature whenever possible to demonstrate and validate the present 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.     

A simple spinning laminated composite shaft model

A simple spinning composite shaft model is presented in this paper. The composite shaft contains discrete isotropic rigid disks and is supported by bearings that are modeled as springs and viscous dampers. Based on a first-order shear deformable beam theory, the strain energy of the shaft are found by adopting the three-dimensional constitutive relations of material with the help of the coordinates transformation, while the kinetic energy of the shaft system is obtained via utilizing the moving rotating coordinate systems adhered to the cross-sections of shaft. The extended Hamilton’s principle is employed to derive the governing equations. In the model the transverse shear deformation, rotary inertia and gyroscopic effects, as well as the coupling effect due to the lamination of composite layers have been incorporated. To verify the present model, the critical speeds of composite shaft systems are compared with those available in the literature. A numerical example is also given to illustrate the frequencies, mode shapes, and transient response of a particular composite shaft system.     

Vibration analysis of nonhomogeneous orthotropic cylindrical shells including combined effect of shear deformation and rotary inertia

This paper is the result of an investigation on the vibration of non-homogeneous orthotropic cylindrical shells, based on the shear deformation theory. Assume that the Young’s moduli, shear moduli and density of the orthotropic material are continuous functions of the coordinate in the thickness direction. The basic equations of non-homogeneous orthotropic cylindrical shells with the shear deformation and rotary inertia are derived in the framework of Donnell-type shell theory. The ends of a non-homogeneous orthotropic cylindrical shell are considered as simply supported. The basic equations are reduced to the sixth-order algebraic equation for the frequency using the Galerkin method. Solving this algebraic equation, the lowest values of non-dimensional frequency parameters for non-homogeneous orthotropic cylindrical shells with and without shear deformation and rotary inertia are obtained. Calculations, effects of shear stresses and rotary inertia, orthotropy, non-homogeneity and shell geometry parameters on the lowest values of non-dimensional frequency parameter are described. The results are verified by comparing the obtained values with those in the existing literature.     

Effects of rotary inertia on sub- and super-critical free vibration of an axially moving beam

The most important issue in the vibration study of an engineering system is dynamics modeling. Axially moving continua is often discussed without the inertia produced by the rotation of the continua section. The main goal of this paper is to discover the effects of rotary inertia on the free vibration characteristics of an axially moving beam in the sub-critical and super-critical regime. Specifically, an integro-partial-differential nonlinear equation is modeled for the transverse vibration of the moving beam based on the generalized Hamilton principle. Then the effects of rotary inertia on the natural frequencies, the critical speed, post-buckling vibration frequencies are presented. Two kinds of boundary conditions are also compared. In super-critical speed range, the straight configuration of the axially moving beam loses its stability. The buckling configurations are derived from the corresponding nonlinear static equilibrium equation. Then the natural frequencies of the post-buckling vibration of the super-critical moving beam are calculated by using local linearization theory. By comparing the critical speed and the vibration frequencies in the sub-critical and super-critical regime, the effects of the inertia moment due to beam section rotation are investigated. Several interesting phenomena are disclosed. For examples, without rotary inertia, the study overestimates the stability of the axially moving beam. Moreover, the relative differences between the super-critical fundamental frequencies of the two theories may increase with an increasing beam length.     

Geometrically nonlinear large deformation analysis of triangular CNT-reinforced composite plates

A first known investigation on the geometrically nonlinear large deformation behavior of triangular carbon nanotube (CNT) reinforced functionally graded composite plates under transversely distributed loads is investigated. The analysis is carried out using the element-free IMLS-Ritz method. In this study, the first-order shear deformation theory (FSDT) and von Kármán assumption are employed to account for transverse shear strains, rotary inertia and moderate rotations. A convergence study is conducted by varying the supporting size and number of nodes. The effects of transverse shear deformation, CNT distribution and CNT volume fraction on the nonlinear bending characteristics under different boundary conditions are examined.     

Response of flexure–torsion coupled composite thin-walled beams with closed cross-sections to random loads

A study of the flexure–torsion coupled random response of the composite beams with solid or thin-walled closed-sections subjected to various types of concentrated and distributed random excitations is dealt with in this paper. The effects of flexure–torsion coupling, shear deformation and rotary inertia are included in the present formulations. The random excitations are assumed to be stationary, ergodic and Gaussian. Analytical expressions for the displacement response of the composite beams are obtained by using normal mode superposition method combined with frequency response function method. The present method can produce the effective solutions for the composite Timoshenko beams with circumferentially antisymmetric (CAS) configuration and more general beam assemblages of connected beams. The influences of flexure–torsion coupling, shear deformation and rotary inertia on the random response of an appropriately chosen composite beam from the literature are demonstrated and discussed.     

Effects of shear stresses and rotary inertia on the stability and vibration of sandwich cylindrical shells with FGM core surrounded by elastic medium

The vibration and stability of axially loaded sandwich cylindrical shells with the functionally graded (FG) core with and without shear stresses and rotary inertia resting Pasternak foundation are investigated. The dynamic stability is derived based on the first order shear deformation theory (FSDT) including shear stresses. The axial load and dimensionless fundamental frequency for FG sandwich shell with shear stresses and rotary inertia and resting on the Pasternak foundation. Finally, the influences of variations of FG core, elastic foundations, shear stresses and rotary inertia on the fundamental frequencies and critical axial loads are investigated.     

Dynamic analysis of 3-D beam elements including warping and shear deformation effects

In this paper the dynamic analysis of 3-D beam elements restrained at their edges by the most general linear torsional, transverse or longitudinal boundary conditions and subjected in arbitrarily distributed dynamic twisting, bending, transverse or longitudinal loading is presented. For the solution of the problem at hand, a boundary element method is developed for the construction of the 14 × 14 stiffness matrix and the corresponding nodal load vector of a member of an arbitrarily shaped simply or multiply connected cross section, taking into account both warping and shear deformation effects, which together with the respective mass and damping matrices lead to the formulation of the equation of motion. To account for shear deformations, the concept of shear deformation coefficients is used, defining these factors using a strain energy approach. Eight boundary value problems with respect to the variable along the bar angle of twist, to the primary warping function, to a fictitious function, to the beam transverse and longitudinal displacements and to two stress functions are formulated and solved employing a pure BEM approach that is only boundary discretization is used. Both free and forced transverse, longitudinal or torsional vibrations are considered, taking also into account effects of transverse, longitudinal, rotatory, torsional and warping inertia and damping resistance. Numerical examples are presented to illustrate the method and demonstrate its efficiency and accuracy. The influence of the warping effect especially in members of open form cross section is analyzed through examples demonstrating the importance of the inclusion of the warping degrees of freedom in the dynamic analysis of a space frame. Moreover, the discrepancy in the dynamic analysis of a member of a spatial structure arising from the ignorance of the shear deformation effect necessitates the inclusion of this additional effect, especially in thick walled cross section members.     

Dynamic analysis of geometrically nonlinear robot manipulators

In this paper, a method for the dynamic analysis of geometrically nonlinear elastic robot manipulators is presented. Robot arm elasticity is introduced using a finite element method which allows for the gross arm rotations. A shape function which accounts for the combined effects of rotary inertia and shear deformation is employed to describe the arm deformation relative to a selected component reference. Geometric elastic nonlinearities are introduced into the formulation by retaining the quadratic terms in the strain-displacement relationships. This has lead to a new stiffness matrix that depends on the rotary inertia and shear deformation and which has to be iteratively updated during the dynamic simulation. Mechanical joints are introduced into the formulation using a set of nonlinear algebraic constraint equations. A set of independent coordinates is identified over each subinterval and is employed to define the system state equations. In order to exemplify the analysis, a two-armed robot manipulator is solved. In this example, the effect of introducing geometric elastic nonlinearities and inertia nonlinearities on the robot arm kinematics, deformations, joint reaction forces and end-effector trajectory are investigated.     

Dynamic stability of axially accelerating Timoshenko beam: Averaging method

This study investigates dynamic stability in transverse parametric vibrations of an axially accelerating tensioned beam of Timoshenko model on simple supports. The axial speed is assumed as a harmonic fluctuation about the constant mean speed. The Galerkin method is applied to discretize the governing equation into a finite set of ordinary differential equations. The method of averaging is applied to analyze the instability phenomena caused by subharmonic and combination resonance. Numerical examples demonstrate the effects of the mean axial speed, bending stiffness, rotary inertia and shear modulus on the instability boundaries.     

Large amplitude free flexural vibration of rings using finite element approach

Here, the large amplitude free flexural vibrations of isotropic/laminated orthotropic rings are investigated, using a shear flexible curved beam element based on field consistency principle. A laminated refined beam theory is introduced for developing the element, which satisfies the interface transverse shear stress and displacement continuity, and has a vanishing shear stress on the inner and outer surfaces of the beam. The formulation includes in-plane and rotary inertia effects, and the non-linearity due to the finite deformation of the ring. The governing equations obtained using Lagrange's equations of motion are solved through the direct integration technique. Amplitude-frequency relationships evaluated from the dynamic response history are examined. Detailed numerical results are presented considering various parameters such as radius-to-thickness ratio, circumferential wave number and ovality for isotropic and laminated orthotropic rings. The nature and degree of the participation of various modes in non-linear asymmetric vibration of oval ring brought out through the present study are useful for accurate modelling of the closed non-circular structures.     

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.     

浸没的球面各向同性球壳自由振动的中厚壳理论

采用六模态的中厚壳理论研究了可压缩流体中球面各向同性球壳的自由振动，即在分析中考虑了剪切变形，旋转惯量和横向正应变的影响，引入５个辅助变量可以理到两类振动及其频率方程，对频率方程作了简化并算例进行了相应的探讨。     

Nonlinear resonance responses of geometrically imperfect shear deformable nanobeams including surface stress effects

In this work, a thorough investigation is presented into the nonlinear resonant dynamics of geometrically imperfect shear deformable nanobeams subjected to harmonic external excitation force in the transverse direction. To this end, the Gurtin–Murdoch surface elasticity theory together with Reddy’s third-order shear deformation beam theory is utilized to take into account the size-dependent behavior of nanobeams and the effects of transverse shear deformation and rotary inertia, respectively. The kinematic nonlinearity is considered using the von Kármán kinematic hypothesis. The geometric imperfection as a slight curvature is assumed as the mode shape associated with the first vibration mode. The weak form of geometrically nonlinear governing equations of motion is derived using the variational differential quadrature (VDQ) technique and Lagrange equations. Then, a multistep numerical scheme is employed to solve the obtained governing equations in order to study the nonlinear frequency–response and force–response curves of nanobeams. Comprehensive studies into the effects of initial imperfection and boundary condition as well as geometric parameters on the nonlinear dynamic characteristics of imperfect shear deformable nanobeams are carried out through numerical results. Finally, the importance of incorporating the surface stress effects via the Gurtin–Murdoch elasticity theory, is emphasized by comparing the nonlinear dynamic responses of the nanobeams with different thicknesses.     

Non-linear dynamic response of a rotating radial Timoshenko beam with periodic pulse loading at the free-end

Consideration is given to the dynamic response of a Timoshenko beam under repeated pulse loading. Starting with the basic dynamical equations for a rotating radial cantilever Timoshenko beam clamped at the hub in a centrifugal force field, a system of equations are derived for coupled axial and lateral motions which includes the transverse shear and rotary inertia effects, as well. The hyperbolic wave equation governing the axial motion is coupled with the flexural wave equation governing the lateral motion of the beam through the velocity-dependent skew-symmetric Coriolis force terms. In the analytical formulation, Rayleigh-Ritz method with a set of sinusoidal displacement shape functions is used to determine stiffness, mass and gyroscopic matrices of the system. The tip of the rotating beam is subjected to a periodic pulse load due to local rubbing against the outer case introducing Coulomb friction in the system. Transient response of the beam with the tip deforming due to rub is discussed in terms of the frequency shift and non-linear dynamic response of the rotating beam. Numerical results are presented for this vibro-impact problem of hard rub with varying coefficients of friction and the contact-load time. The effects of beam tip rub forces transmitted through the system are considered to analyze the conditions for dynamic stability of a rotating blade with intermittent rub.     

无约束修正Timoshenko梁的冲击问题

介绍了修正后的Timoshenko梁运动方程，并比较了修正Timoshenko梁与经 典Timoshenko梁的运动方程. 推导了考虑剪切变形引起的转动惯量的修正Timoshenko 梁的正交条件，推导了集中质量对无约束修正Timoshenko梁的正碰撞对梁所引起的瞬态冲 击响应公式，并用算例进行了分析，且与集中质量对经典的无约束Timoshenko梁的正碰撞 对梁所引起的冲击响应进行了比较，另外还用算例分析了梁的刚度的变化和冲击质量比对其 冲击响应产生的影响.     

Dynamic stability analysis of shear-flexible composite beams

Dynamic stability behavior of the shear-flexible composite beams subjected to the nonconservative force is intensively investigated based on the finite element model using the Hermitian beam elements. For this, a formal engineering approach of the mechanics of the laminated composite beam is presented based on kinematic assumptions consistent with the Timoshenko beam theory, and the shear stiffness of the thin-walled composite beam is explicitly derived from the energy equivalence. An extended Hamilton’s principle is employed to evaluate the mass-, elastic stiffness-, geometric stiffness-, damping-, and load correction stiffness matrices. Evaluation procedures for the critical values of divergence and flutter loads of the nonconservative system with and without damping effects are then briefly introduced. In order to verify the validity and the accuracy of this study, the divergence and flutter loads are presented and compared with the results from other references, and the influence of various parameters on the divergence and flutter behavior of the laminated composite beams is newly addressed: (1) variation of the divergence and flutter loads with or without the effects of shear deformation and rotary inertia with respect to the nonconservativeness parameter and the fiber angle change, (2) influence of the internal and external damping on flutter loads whether to consider the shear deformation or not.