薄壁曲梁的稳定性研究进展

 曲梁是桥梁、建筑、船舶、航空和航天工程中常见的薄壁构件，根据外载荷与主曲率平面的关系，又被称为拱或水平曲梁．随着工程材料的日益发展，如复合材料、功能梯度材料的引入，曲梁的应用范围更加广泛，进一步推进了薄壁曲梁稳定性问题的研究．本文首先对薄壁梁结构的稳定性行为进行了分类．接着简述了薄壁构件的基本假设，对比了近几十年来薄壁曲梁的基本理论，针对复合材料薄壁曲梁，总结了相应的本构关系，并对各理论间存在的分歧进行了归纳．结合最新的薄壁曲梁研究，根据平衡法、能量法和虚位移（虚功）原理推导出控制微分方程，阐述了相应的求解方法，如解析法、半解析法和数值解法．为验证薄壁曲梁理论的准确性，曲梁承载能力试验验证尤为重要，但目前国内外相关研究还很少，亟待发展．最后讨论了现阶段薄壁曲梁研究的局限性和未来发展的方向．     

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.     

Shear deformable finite beam elements for composite box beams

The shear deformable thin-walled composite beams with closed cross-sections have been developed for coupled flexural, torsional, and buckling analyses. A theoretical model applicable to the thin-walled laminated composite box beams is presented by taking into account all the structural couplings coming from the material anisotropy and the shear deformation effects. The current composite beam includes the transverse shear and the restrained warping induced shear deformation by using the first-order shear deformation beam theory. Seven governing equations are derived for the coupled axial-flexural-torsional-shearing buckling based on the principle of minimum total potential energy. Based on the present analytical model, three different types of finite composite beam elements, namely, linear, quadratic and cubic elements are developed to analyze the flexural, torsional, and buckling problems. In order to demonstrate the accuracy and superiority of the beam theory and the finite beam elements developed by this study,numerical solutions are presented and compared with the results obtained by other researchers and the detailed threedimensional analysis results using the shell elements of ABAQUS. Especially, the influences of the modulus ratio and the simplified assumptions in stress–strain relations on the deflection, twisting angle, and critical buckling loads of composite box beams are investigated.     

Refined Series Methodology for the Fully Coupled Thin-Walled Laminated Beams Considering Foundation Effects

The refined power series solutions are presented for the coupled static analysis of thin-walled laminated beams resting on elastic foundation. For this purpose, the elastic strain energy considering the material and structural coupling effects and the energy including the foundation effects are constructed. The equilibrium equations and the force-displacement relationships are derived from the extended Hamilton's principle, and the explicit expressions for displacement parameters are presented based on power series expansions of displacement components. Finally, the member stiffness matrix is determined by using the force-displacement relationships. For comparison, the finite element model based on the Hermite cubic interpolation polynomial is presented. In order to verify the accuracy and the superiority of the laminated beam element developed by this study, the numerical solutions are presented and compared with results obtained from the regular finite beam elements and the ABAQUS's shell elements. The influences of the fiber angle change and the boundary conditions on the coupled behavior of laminated beams with mono-symmetric I-sections are investigated.     

Improved torsional analysis of laminated box beams

The improved torsional analysis of the laminated box beams with single- and double-celled sections subjected to a torsional moment is performed by introducing 14 displacement parameters. For this, a thin-walled laminated box beam theory considering the effects of shear and elastic couplings is presented. The governing equations and the force-displacement relations are derived from the variation of the strain energy. The system of linear algebraic equations with non-symmetric matrix is constructed by introducing the displacement parameters and by transforming the higher order simultaneous differential equations into first order ones. This numerical technique determines eigenmodes corresponding to 12 zero and 2 non-zero eigenvalues and derives displacement functions for displacement parameters based on the undetermined parameter method. Finally, the element stiffness matrix is determined using the member force-displacement relations. The theory developed by this study is validated by comparing several torsional responses from the present approach with those from the finite element beam model using the Lagrangian interpolation polynomials and three-dimensional analysis results using the shell elements of ABAQUS for coupled laminated beams with single- and double-celled sections.     

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.     

剪切流作用下层合梁非线性振动特性研究

针对剪切流中层合梁的大变形非线性振动问题, 采用高阶剪切变形锯齿理论和冯·卡门应变描述层合梁的变形模式和几何非线性效应, 构建了大变形层合梁非线性振动有限元数值模型; 采用基于任意拉格朗日?欧拉方法的有限体积法求解不可压缩黏性流体纳维-斯托克斯方程, 结合层合梁和流体的耦合界面条件建立了剪切流作用下层合梁流固耦合非线性动力学数值模型, 采用分区并行强耦合方法对层合梁的流致非线性振动响应进行了迭代计算. 研究了不同速度分布的剪切流作用下单层梁和多层复合材料梁的振动响应特性, 并验证了本文数值建模方法的有效性. 结果表明: 剪切流作用下单层梁的振动特性与均匀流作用下的情况不同, 梁的运动轨迹受剪切流影响向下偏斜, 随着速度分布系数增加, 尾部流场中的涡结构发生改变; 刚度比对剪切流作用下层合梁的振动特性有显著影响, 随着刚度比的增加, 层合梁振动的振幅增大, 主导频率下降, 运动轨迹由‘8’字形逐渐变得不对称; 发现了不同厚度比和铺层角度情况下, 层合梁存在定点稳定模式、周期极限环振动模式和非周期振动模式三种不同的振动模式, 改变层合梁铺层角度可实现层合梁周期极限环振动模式向非周期振动模式转变.      

Stochastic response of an axially loaded composite Timoshenko beam exhibiting bending–torsion coupling

A spectral finite element method is proposed to investigate the stochastic response of an axially loaded composite Timoshenko beam with solid or thin-walled closed section exhibiting bending–torsion materially coupling under the stochastic excitations with stationary and ergodic properties. The effects of axial force, shear deformation (SD) and rotary inertia (RI) as well as bending–torsion coupling are considered in the present study. First, the damped general governing differential equations of motion of an axially loaded composite Timoshenko beam are derived. Then, the spectral finite element formulation is developed in the frequency domain using the dynamic shape functions based on the exact solutions of the governing equations in undamped free vibration, which is used to compute the mean square displacement response of axially loaded composite Timoshenko beams. Finally, the proposed method is illustrated by its application to a specific example to investigate the effects of bending–torsion coupling, axial force, SD and RI on the stochastic response of the composite beam.     

A robust hyper-elastic beam model under bi-axial normal-shear loadings

In this paper, a simple and robust constitutive model is proposed to simulate mechanical behaviors of hyper-elastic materials under bi-axial normal-shear loadings in the finite strain regime. The Mooney–Rivlin strain energy function is adopted to develop a two-dimensional (2D) normal-shear constitutive model within the framework of continuum mechanics. A motion field is first proposed for combined normal and shear deformations. The deformation gradient of the proposed field is calculated and then substituted into right Cauchy–Green deformation tensor. Constitutive equations are then derived for normal and shear deformations. They are two explicit coupled equations with high-level polynomial non-linearity. In order to examine capabilities of the developed hyper-elastic model, uniaxial tensile responses and non-linear stability behaviors of moderately thick straight and curved beams undergoing normal axial and transverse shear deformations are simulated and compared with experiments. Fused deposition modeling technique as a 3D printing technology is implemented to fabricate hyper-elastic beam structures from soft poly-lactic acid filaments. The printed specimens are tested under tensile/compressive in-plane and compressive out-of-plane forces. A finite element formulation along with the Newton–Raphson and Riks techniques is also developed to trace non-linear equilibrium path of beam structures in large defamation regimes. It is shown that the model is capable of predicting non-linear equilibrium characteristics of hyper-elastic straight and curved beams. It is found that the modeling of shear deformation and finite strain is essential toward an accurate prediction of the non-linear equilibrium responses of moderately thick hyper-elastic beams. Due to simplicity and accuracy, the model can serve in the future studies dealing with the analysis of hyper-elastic structures in which two normal and shear stress components are dominant.     

Non-linear stability analysis of imperfect thin-walled composite beams

The static non-linear behavior of thin-walled composite beams is analyzed considering the effect of initial imperfections. A simple approach is used for determining the influence of imperfection on the buckling, prebuckling and postbuckling behavior of thin-walled composite beams. The fundamental and secondary equilibrium paths of perfect and imperfect systems corresponding to a major imperfection are analyzed for the case where the perfect system has a stable symmetric bifurcation point. A geometrically non-linear theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field. An initial displacement, either in vertical or horizontal plane, is considered in presence of initial geometric imperfection. Ritz's method is applied in order to discretize the non-linear differential system and the resultant algebraic equations are solved by means of an incremental Newton-Rapshon method. The numerical results are presented for a simply supported beam subjected to axial or lateral load. It is shown in the examples that a major imperfection reduces the load-carrying capacity of thin-walled beams. The influence of this effect is analyzed for different fiber orientation angle of a symmetric balanced lamination. In addition, the postbuckling response obtained with the present beam model is compared with the results obtained with a shell finite element model (Abaqus).     

曲线梁研究进展

曲线梁在工程中的应用非常广泛, 但由于曲率的影响, 导致梁内弯矩和扭矩的耦合, 使得 曲线梁的研究相对直梁非常复杂. 本文从4个方面较为详细地评述了近年来国内外曲线梁的 研究进展情况. 从曲线梁的平衡方程、几何方程和基本微分方程出发, 概述了曲线梁静动力 学的基本理论、建模及分析方法; 面内面外振动及分析方法; 非线性问题及分析方法; 特别 评述了有限元单元法在曲线梁研究中的应用, 论述了各种曲梁单元的优缺点; 对复合曲线梁 的发展及理论研究进行介绍; 推导了空间曲线梁非线性动力学方程; 最后提出了值得进一步 研究的方向及采用的方法.     

Generalized variational principle of dynamic analysis on naturally curved and twisted box beams for anisotropic materials

An incomplete generalized variational functional for naturally curved and twisted composite box beams with complete constrained boundaries at two ends is established by means of Lagrange multiplier method. The equations of motion governing the dynamic behavior of the beams and corresponding boundary conditions are derived from the stationary condition of the functional. The non-classical influences relevant to the beams are those due to transverse shear deformations, torsion-related warping and several elastic couplings that can arise in composite beams. In order to demonstrate the correctness of the theory developed the natural frequencies and normal mode shapes of the beams under in-plane free vibration are evaluated and compared with the results using PATRAN’s beam elements.     

Nonlinear vibration of functionally graded graphene-reinforced composite laminated beams resting on elastic foundations in thermal environments

Modeling and nonlinear vibration analysis of graphene-reinforced composite (GRC) laminated beams resting on elastic foundations in thermal environments are presented. The graphene reinforcements are assumed to be aligned and are distributed either uniformly or functionally graded of piece-wise type along the thickness of the beam. The motion equations of the beams are based on a higher-order shear deformation beam theory and von Kármán strain displacement relationships. The beam–foundation interaction and thermal effects are also included. The temperature-dependent material properties of GRCs are estimated through a micromechanical model. A two-step perturbation approach is employed to determine the nonlinear-to-linear frequency ratios of GRC laminated beams. Detailed parametric studies are carried out to investigate the effects of material property gradient, temperature variation, stacking sequence as well as the foundation stiffness on the linear and nonlinear vibration characteristics of the GRC laminated 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.     

??????????——-??????????????????????

研制和开发了曲梁以及复合曲梁测量应力的创新实验装置, 通过该实验的设计、开发和应用,可以验证它的创新性和综合性,找出曲梁、复合曲梁与直 梁的诸多不同之处. 把该实验用于测试由钢制成的、钢与铜两种材料制成的,具有直角梯形 截面的简支梁在拱顶处受垂直集中力作用时的正应力和切应力,计算结果表明,理论解和实 验值吻合得很好.     

A continuum theory for viscoelastic composite beams

Summary A dynamical continuum theory is developed for laminated composite beams. Starting with an assumed displacement- and temperature field, the one-dimensional approximate theory is consistently constructed within the frame of the three-dimensional theory of linear, nonisothermal, anisotropic, coupled viscoelasticity. Each constituent of the beam may possess different constant thickness and mechanical properties. All dynamic interactions between the adjacent constituents are included. Further, the effects of transverse shear and normal strains and rotatory inertia as well as those of cross-sectional distortion are all taken into account. The resulting equations consist of the macroscopic beam equations of motion and heat conduction, the kinematical relations, the initial and boundary conditions and the constitutive equations, and they govern the extensional, flexural and torsional motions of laminated composite beams. The special cases of constituents which made of either isotropic thermoviscoelastic or anisotropic thermoelastic materials are discussed briefly.Supported by the Office of Naval Research.With 1 figure     

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.     

Nonlinear transient response of strain rate dependent composite laminated plates using multiscale simulation

The effects of strain rate dependency and inelasticity on the transient responses of composite laminated plates are investigated. A micromechanics model which accounts for the transverse shear stress effect, the effect of strain rate dependency and the effect of inelasticity is used for analyzing the mechanical responses of the fiber and matrix constituents. The accuracy of the micromechanics model under transverse shear loading is verified by comparing the results with those obtained using a general purpose finite element code. A higher order laminated plate theory is extended to capture the inelastic deformations of the composite plate and is implemented using the finite element technique. A complete micro–macro numerical procedure is developed to model the strain rate dependent behavior of inelastic composite laminates by implementing the micromechanics model into the finite element model. Parametric studies of the transient responses of composite plates are conduced. The effects of geometry, ply stacking sequence, material models, boundary conditions and loadings are investigated. The results show that the strain rate dependency and inelasticity influence the transient responses of composite plates via two significantly different mechanisms.     

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.