非惯性系下考虑剪切变形的柔性梁的动力学建模

研究非惯性坐标系下考虑剪切变形的柔性梁的动力学建模. 首先借鉴Euler-Bernoulli梁的几何非线性变形模式,考虑了Timoshenko梁弯曲以及剪切变形产生的几何非线性效应对纵向、横向变形位移的影响,在考虑两个方向的变形耦合项后,利用有限元法对柔性梁进行了离散,采用Lagrange方程建立了柔性梁的动力学模型,首次建立了包含变形二次耦合量的Timoshenko梁的动力学方程. 关键词： 非惯性坐标系 剪切变形 柔性梁 动力学建模     

Convergence of Galerkin truncation for dynamic response of finite beams on nonlinear foundations under a moving load

The present paper investigates the convergence of the Galerkin method for the dynamic response of an elastic beam resting on a nonlinear foundation with viscous damping subjected to a moving concentrated load. It also studies the effect of different boundary conditions and span length on the convergence and dynamic response. A train–track or vehicle–pavement system is modeled as a force moving along a finite length Euler–Bernoulli beam on a nonlinear foundation. Nonlinear foundation is assumed to be cubic. The Galerkin method is utilized in order to discretize the nonlinear partial differential governing equation of the forced vibration. The dynamic response of the beam is obtained via the fourth-order Runge–Kutta method. Three types of the conventional boundary conditions are investigated. The railway tracks on stiff soil foundation running the train and the asphalt pavement on soft soil foundation moving the vehicle are treated as examples. The dependence of the convergence of the Galerkin method on boundary conditions, span length and other system parameters are studied.     

Analysis of nonlinear dynamic response for delaminated fiber–metal laminated beam under unsteady temperature field

The nonlinear dynamic response problems of fiber–metal laminated beams with delamination are studied in this paper. Basing on the Timoshenko beam theory, and considering geometric nonlinearity, transverse shear deformation, temperature effect and contact effect, the nonlinear governing equations of motion for fiber–metal laminated beams under unsteady temperature field are established, which are solved by the differential quadrature method, Nermark-β method and iterative method. In numerical examples, the effects of delamination length, delamination depth, temperature field, geometric nonlinearity and transverse shear deformation on the nonlinear dynamic response of the glass reinforced aluminum laminated beam with delamination are discussed in details.     

Higher order representation of the beam cross section deformation in large displacement finite element analysis

Most existing beam formulations assume that the cross section of the beam remains rigid regardless of the amplitude of the displacement. The absolute nodal coordinate formulation (ANCF); however, allows for the deformation of the cross section and leads to a more general beam models that capture the coupling between different modes of displacement. This paper examines the effect of the order of interpolation on the modes of deformation of the beam cross section using ANCF finite elements. To this end, a new two-dimensional shear deformable ANCF beam element is developed. The new finite element employs a higher order of interpolation, and allows for new cross section deformation modes that cannot be captured using previously developed shear deformable ANCF beam elements. The element developed in this study relaxes the assumption of planar cross section; thereby allowing for including the effect of warping as well as for different stretch values at different points on the element cross section. The displacement field of the new element is assumed to be cubic in the axial direction and quadratic in the transverse direction. Using this displacement field, more expressions for the element extension, shear and the cross section stretch can be systematically defined. The change in the cross section area is measured using Nanson’s formula. Measures of the shear angle, extension, and cross section stretch can also be systematically defined using coordinate systems defined at the element material points. Using these local coordinate systems, expressions for a nominal shear angle are obtained. The differences between the cross section deformation modes obtained using the new higher order element and those obtained using the previously developed lower order elements are highlighted. Numerical examples are presented in order to compare the results obtained using the new finite element and the results obtained using previously developed ANCF finite elements.     

Exact closed-form solutions for the natural frequencies and stability of elastically connected multiple beam system using Timoshenko and high-order shear deformation theory

A general procedure for the determination of the natural frequencies and buckling load for a set of beam system under compressive axial loading is investigated using Timoshenko and high-order shear deformation theory. It is assumed that the set beams of the system are simply supported and continuously joined by a Winkler elastic layer. The model of beam includes the effects of axial loading, shear deformation and rotary inertia. In the special case of identical beams, explicit expressions for the natural frequencies and the critical buckling load are determined using a trigonometric method. The influences of the compressive axial loading and the number of beams in the system on the natural frequencies and the critical buckling load are discussed. These results are of considerable practical interest and have wide application in engineering practice of frameworks.     

Nonlinear nonuniform torsional vibrations of bars by the boundary element method

In this paper a boundary element method is developed for the nonuniform torsional vibration problem of bars of arbitrary doubly symmetric constant cross-section taking into account the effect of geometrical nonlinearity. The bar is subjected to arbitrarily distributed or concentrated conservative dynamic twisting and warping moments along its length, while its edges are supported by the most general torsional boundary conditions. The transverse displacement components are expressed so as to be valid for large twisting rotations (finite displacement-small strain theory), thus the arising governing differential equations and boundary conditions are in general nonlinear. The resulting coupling effect between twisting and axial displacement components is considered and torsional vibration analysis is performed in both the torsional pre- or post-buckled state. A distributed mass model system is employed, taking into account the warping, rotatory and axial inertia, leading to the formulation of a coupled nonlinear initial boundary value problem with respect to the variable along the bar angle of twist and to an “average” axial displacement of the cross-section of the bar. The numerical solution of the aforementioned initial boundary value problem is performed using the analog equation method, a BEM based method, leading to a system of nonlinear differential-algebraic equations (DAE), which is solved using an efficient time discretization scheme. Additionally, for the free vibrations case, a nonlinear generalized eigenvalue problem is formulated with respect to the fundamental mode shape at the points of reversal of motion after ignoring the axial inertia to verify the accuracy of the proposed method. The problem is solved using the direct iteration technique (DIT), with a geometrically linear fundamental mode shape as a starting vector. The validity of negligible axial inertia assumption is examined for the problem at hand.     

新型低电压大变形微驱动器数值求解及仿真

针对横向加载单向变形MEMS静电微驱动器位移过小或驱动电压过大问题,提出一种基于纵横弯曲变形原理的硅基大位移低驱动电压静电驱动器模型;基于拉格朗日-麦克斯韦机电动分析力学,建立轴向横向同时加载的微驱动器动力方程;分析温度应力、静电调节力和轴向挤压力对轴向载荷的影响;基于龙格-库塔算法和有限差分法分别将横向分布载荷和轴向载荷等效转化为横向集中载荷;仿真得到变形同驱动电压、调节电压、轴向挤压量和温差的关系;结果显示当驱动电压仅为16 V时,位移高达10.861μm,远大于传统横向加载单向变形微驱动器的位移量.实验验证了仿真结果.     

The refined theory of deep rectangular beams for symmetrical deformation

Based on elasticity theory, various one-dimensional equations for symmetrical deformation have been deduced systematically and directly from the two-dimensional theory of deep rectangular beams by using the Papkovich-Neuber solution and the Lur’e method without ad hoc assumptions, and they construct the refined theory of beams for symmetrical deformation. It is shown that the displacements and stresses of the beam can be represented by the transverse normal strain and displacement of the mid-plane. In the case of homogeneous boundary conditions, the exact solutions for the beam are derived, and the exact equations consist of two governing differential equations: the second-order equation and the transcendental equation. In the case of non-homogeneous boundary conditions, the approximate governing differential equations and solutions for the beam under normal loadings only and shear loadings only are derived directly from the refined beam theory, respectively, and the correctness of the stress assumptions in classic extension or compression problems is revised. Meanwhile, as an example, explicit expressions of analytical solutions are obtained for beams subjected to an exponentially distributed load along the length of beams. Supported by the National Natural Science Foundation of China (Grant Nos. 10702077, 10672001, and 10602001), the Beijing Natural Science Foundation (Grant No. 1083012), and the Alexander von Humboldt Foundation in Germany     

An efficient coupled layerwise theory for dynamic analysis of piezoelectric composite beams

An efficient new coupled one-dimensional model is developed for the dynamics of piezoelectric composite beams. The model combines third order zigzag approximation for the displacement with layerwise approximation of the electric field as piecewise linear for sublayers. By enforcing the conditions of zero transverse shear stress at the top and bottom and its continuity at layer interfaces, the displacement field is expressed in terms of three primary displacement variables and potentials. The governing coupled equations of stress and charge equilibrium and boundary conditions are derived from Hamilton's principle. Analytical solutions are obtained, for free vibrations and forced response under harmonic load, for simply supported hybrid beams and the results are compared with the exact three-dimensional solution and uncoupled first order shear deformation theory solution. The present results show significant improvement over the first order solution and agree very well with the exact solution for both thin and thick hybrid beams. The results demonstrate the capability of the developed theory to adequately model open and closed circuit electric boundary conditions to accurately predict their influence on the response.     

轴向加速运动黏弹性梁受迫振动中的混沌动力学

研究了黏弹性轴向运动梁在外部激励和参数激励共同作用下横向振动的混沌非线性动力学行为. 引入有限支撑刚度, 并考虑黏弹性本构关系取物质导数, 同时计入由梁轴向加速度引起的沿径向变化的轴力, 建立轴向运动黏弹性梁横向非线性振动的偏微分-积分模型. 通过Galerkin截断方法研究了外部激励的频率和因速度简谐脉动引起的参数激励的频率在不可通约关系时轴向运动连续体的非线性动力学行为, 并对不同截断阶数的数值预测进行了对比. 基于对控制方程的Galerkin截断, 得到离散化的常微分方程组, 使用四阶Runge-Kutta方法求解. 基于此数值解, 运用非线性动力学时间序列分析方法, 通过Poincaré 映射, 观察到轴向运动梁随扰动速度幅值的倍周期分岔现象, 并比较了有无外部激励对倍周期分岔的影响. 分别在低速以及近临界高速运动状态下, 从相平面图、Poincaré 映射以及频谱分析的角度识别了系统中存在的准周期运动形态. 关键词： 轴向运动梁 非线性 混沌 分岔     

Frequency analysis of finite beams on nonlinear Kelvin-Voight foundation under moving loads

The vibration of an Euler-Bernoulli beam, resting on a nonlinear Kelvin-Voight viscoelastic foundation, traversed by a moving load is studied in the frequency domain. The objective is to obtain the frequency responses of the beam and the effects of different parameters on the system response. The parameters include the magnitude and speed of the moving load and the foundation nonlinearity and its damping coefficient. The solution is obtained by using the Galerkin method in conjunction with the multiple scales method (MSM). The governing nonlinear partial differential equations of motion are discretized into sets of nonlinear ordinary differential equations. Subsequently, the solution is calculated for different harmonics by using the MSM as one of the powerful perturbation techniques. The steady-state responses of the main harmonic as well as its two super-harmonics are then obtained. As a case study, a conventional railway track is dynamically simulated and the jump phenomenon in the response is observed for three harmonics. Moreover, a thorough stability analysis of the system is carried out.     

Fundamental frequency evolution in slender beams subjected to imposed axial displacements

The influence of applied axial loads on the fundamental vibration frequency is strictly connected with the stability analysis of elastic slender beams. For this reason, the correct evaluation of the fundamental frequency is of primary importance in designing new structures and components, as well as in monitoring existing ones. At the same time, if an internal axial load arises in a slender element as the consequence of an imposed (static) axial end displacement, then a different dynamic structural response is encountered respect to the case in which a beam end is free to slide, during transverse vibration, and a (constant) axial load is applied externally. This difference is due to the change in the axial boundary condition. Moreover, the presence of an initial curvature of the beam axis may significantly affect the aforesaid response. The experimental study proposed in the present paper investigates the dependence of the fundamental frequency on the axial load in slender beams subjected to imposed axial end displacements. The considered specimens presented different geometrical imperfections (initial curvatures), and were tested in two different constraint conditions (hinged–hinged and hinged–clamped). In addition, the behaviors observed during the experiments were reproduced by numerical simulations offering a valid confirmation for test results and contributing to understand the evolution of the fundamental frequency in the analyzed slender elements subjected to imposed axial end displacements.     

Analysis of energy dissipation in an elastic moving string with a viscous damper at one end

In this paper transverse vibration of an axially moving viscoelastic string with a viscous damper at one end is investigated analytically. The string is assumed to be travelling with constant velocity and the length of string is constant or time varying. The linear and nonlinear mathematical models are derived using the Lagrangian function and implemented using a finite element method. The method considers a time varying state space function applied to the linear model, the Newmark-Beta method is used to solve the response for the nonlinear problem numerically. The case of energy dissipated by a viscoelastic damper at one end of the string for different axial string velocities is considered. When a disturbance arrives at the boundary an exact value for the damper which provides maximum energy dissipation is investigated. Finally, numerical simulations are presented to establish the feasibility of the method.     

Vibration control of beams on elastic foundation under a moving vehicle and random lateral excitations

The formulation of three-dimensional dynamic behavior of a Beam On Elastic Foundation (BOEF) under moving loads and a moving mass is considered. The weight of the vehicle is modeled as a moving point load, however the effect of the lateral excitation is considered by modeling: (case 1) a lateral moving load with random intensity for wind excitation and (case 2) a moving mass just in lateral direction of the beam for earthquake excitation. A Dirac-delta function is used to describe the position of the moving load and the moving mass along the beam. The beam foundations are considered as elastic Winkler-type in two perpendicular transverse directions. This model is proposed to investigate the bending response of the rails under the effect of traveling vehicle weight while a random excitation such as earthquake or wind takes place. The results showed the importance of considering the effect of earthquake/wind actions as in bending stress of the beam on elastic foundations. The effect of different regions (different support stiffness) and different velocities of the vehicle on the response of the beam are investigated in mentioned directions. At the end, a linear optimal control algorithm with displacement–velocity feedback is proposed as a solution to suppress the response of BOEFs. By the method of modal analyses and taking into account enough number of vibration modes, state-space equation is obtained, then sufficient number of actuators was chosen for each direction. Stochastic analyses were performed in lateral direction in order to illustrate a comprehensive view for the response of the beam under the random moving load in both controlled and uncontrolled systems. Furthermore, the efficiency of control algorithm on critical velocities is verified by parametric analyses in the vertical direction with the constant moving load for different regions.     

DYNAMIC RESPONSE OF A ROTATING FLEXIBLE ARM CARRYING A MOVING MASS

A clamped-free flexible arm rotating in a horizontal plane and carrying a moving mass is studied in this paper. The arm is modelled by the Euler-Bernoulli beam theory in which rotatory inertia and shear deformation effects are ignored. The assumed mode method in conjunction with Hamilton's principle is used to derive the equation of motion of the system which takes into account the effect of centrifugal stiffening due to the rotation of the beam. The eigenfunctions of a cantilever beam which satisfy the prescribed geometric boundary conditions are used as basis functions in the assumed mode method. The equation of motion is expressed in non-dimensional matrix form. Pre-designed transformed cosine profiles are used as trajectory inputs for the hub angle and the moving mass. The equation of motion is solved numerically using the fourth order Runge-Kutta method. Graphical results are presented to show the influence of centrifugal stiffening effect, moving mass values, mass travelling time, hub angle and mass trajectory profile on the deflection of the beam.     

Galerkin method for steady-state response of nonlinear forced vibration of axially moving beams at supercritical speeds

The present paper investigates the steady-state periodic response of an axially moving viscoelastic beam in the supercritical speed range. The straight equilibrium configuration bifurcates in multiple equilibrium positions in the supercritical regime. It is assumed that the excitation of the forced vibration is spatially uniform and temporally harmonic. Under the quasi-static stretch assumption, a nonlinear integro-partial-differential equation governs the transverse motion of the axially moving beam. The equation is cast in the standard form of continuous gyroscopic systems via introducing a coordinate transform for non-trivial equilibrium configuration. For a beam constituted by the Kelvin model, the primary resonance is analyzed via the Galerkin method under the simply supported boundary conditions. Based on the Galerkin truncation, the finite difference schemes are developed to verify the results via the method of multiple scales. Numerical simulations demonstrate that the steady-state periodic responses exist in the transverse vibration and a resonance with a softening-type behavior occurs if the external load frequency approaches the linear natural frequency in the supercritical regime. The effects of the viscoelastic damping, external excitation amplitude, and nonlinearity on the steady-state response amplitude for the first mode are illustrated.     

Dynamic instability of composite plates subjected to non-uniform in-plane loads

In this paper, the dynamic instability of a shear deformable composite plate subjected to periodic non-uniform in-plane loading is studied for four sets of boundary conditions. The static component and the dynamic component of the applied periodic in-plane loading are assumed to vary according to either parabolic or linear distributions. Initially, the plate membrane problem is solved using the Ritz method to evaluate the plate in-plane stress distributions within the prebuckling range due to the applied non-uniform in-plane edge loading. Subsequently using the evaluated stress distribution within the plate, the equations governing the plate instability boundaries are formulated via Hamilton's variational principle. Employing Galerkin's method, these partial differential equations are reduced into a set of ordinary differential equations (Mathieu type of equations) describing the plate dynamic instability behaviour. Following Bolotin's method, the instability regions are determined from the boundaries of instability, which represents the periodic solution of the differential equations with period T and 2T to the Mathieu equations. The instability regions are determined for uniform, linear and parabolic dynamic in-plane loads using first-order and second-order approximations. Numerical results are also presented to bring out the effects of span to thickness ratio, shear deformation, aspect ratio, boundary conditions and static load factor on the instability regions.     

A conserved quantity and the stability of axially moving nonlinear beams

Free nonlinear transverse vibration is investigated for an axially moving beam modeled by an integro-partial-differential equation. Based on the equation, a conserved quantity is defined and confirmed for axially moving beams with pinned or clamped ends. The conserved quantity is applied to demonstrate the Lyapunov stability of the straight equilibrium configuration in transverse nonlinear of beam with a low axial speed.     

DYNAMIC STABILITY OF CURVED PANELS WITH CUTOUTS

The parametric instability behaviour of curved panels with cutouts subjected to in-plane static and periodic compressive edge loadings are studied using finite element analysis. The first order shear deformation theory is used to model the curved panels, considering the effects of transverse shear deformation and rotary inertia. The theory used is the extension of dynamic, shear deformable theory according to Sanders' first approximation for doubly curved shells, which can be reduced to Love's and Donnell's theories by means of tracers. The effects of static and dynamic load factors, geometry, boundary conditions and the cutout parameters on the principal instability regions of curved panels with cutouts are studied in detail using Bolotin's method. Quantitative results are presented to show the effects of shell geometry and load parameters on the stability boundaries. Results for plates are also presented as special cases and are compared with those available in the literature.     

On the nonlinear axisymmetric dynamic buckling behavior of clamped functionally graded spherical caps

Here, the dynamic thermal buckling behavior of functionally graded spherical caps is studied considering geometric nonlinearity based on von Karman's assumptions. The formulation is based on first-order shear deformation theory and it includes the in-plane and rotary inertia effects. The material properties are graded in the thickness direction according to the power-law distribution in terms of volume fractions of the material constituents. The effective material properties are evaluated using homogenization method. The governing equations obtained using finite element approach are solved employing the Newmark's integration technique coupled with a modified Newton–Raphson iteration scheme. The pressure load corresponding to a sudden jump in the maximum average displacement in the time history of the shell structure is taken as the dynamic buckling load. The present model is validated against the available isotropic case. A detailed numerical study is carried out to highlight the influences of shell geometries, power law index of functional graded material and boundary conditions on the dynamic buckling load of shallow spherical shells.