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.     

General analytical solution to bending of composite laminated beams with delaminations

Based on the first-order shear deformable beam theory, a refined model for composite beams containing a through-the-width delamination is presented, and the deformation at the delamination front is considered. Different from the ordinary delami- nated beam theory, each of the perfectly bonded portions of the new model is constructed as two separated beams along the interface without assuming a plane section at the de- lamination front. The governing equations of the delaminated portions and bonded ones are established, combined with continuity conditions of displacements and internal forces. Solutions of delaminated composite beams with different boundary conditions, delamina- tion locations and sizes axe shown in excellent agreement with the finite element results, showing efficiency and applicability of the present model.     

复合材料层合梁的屈曲

本文在铁摩辛柯梁理论基础上,利用迭合刚度方法及Hamilton原理建立了层合梁屈曲问题控制方程,并用此控制方程求解了在具体边界条件下层合梁的屈曲问题,得出了无论在什么边界条件下层合梁的最小屈曲载荷不会大于等效剪切刚度系数C的结论.     

Thermal buckling and postbuckling of laminated composite beams with temperature-dependent properties

The thermal buckling and postbuckling analysis of laminated composite beams with temperature-dependent material properties is presented. The governing equations are based on the first-order shear deformation beam theory (FSDT) and the geometrical nonlinearity is modeled using Green's strain tensor in conjunction with the von Karman assumptions. The differential quadrature method (DQM) as an accurate, simple and computationally efficient numerical tool is adopted to discretize the governing equations and the related boundary conditions. A direct iterative method is employed to obtain the critical temperature (bifurcation point) as well as the nonlinear equilibrium path (the postbuckling behavior) of symmetrically laminated beams. The applicability, rapid rate of convergence and high accuracy of the method are established via different examples and by comparing the results with those of existing in literature. Then, the effects of temperature dependence of the material properties, boundary conditions, length-to-thickness ratios, number of layers and ply angle on the thermal buckling and postbuckling characteristic of symmetrically laminated beams are investigated.     

Exact buckling solution for two-layer Timoshenko beams with interlayer slip

This paper deals with the buckling behavior of two-layer shear-deformable beams with partial interaction. The Timoshenko kinematic hypotheses are considered for both layers and the shear connection (no uplift is permitted) is represented by a continuous relationship between the interface shear flow and the corresponding slip. A set of differential equations is obtained from a general 3D bifurcation analysis, using the above assumptions. Original closed-form analytical solutions of the buckling load and mode of the composite beam under axial compression are derived for various boundary conditions. The new expressions of the critical loads are shown to be consistent with the ones corresponding to the Euler–Bernoulli beam theory, when transverse shear stiffnesses go to infinity. The proposed analytical formulae are validated using 2D finite element computations. Parametric analyses are performed, especially including the limiting cases of perfect bond and no bond. The effect of shear flexibility is particularly emphasized.     

A shearable and thickness stretchable finite strain beam model for soft structures

Soft materials and structures have recently attracted lots of research interests as they provide paramount potential applications in diverse fields including soft robotics, wearable devices, stretchable electronics and biomedical engineering. In a previous work, an Euler–Bernoulli finite strain beam model with thickness stretching effect was proposed for soft thin structures subject to stiff constraint in the width direction. By extending that model to account for the transverse shear effect, a Timoshenko-type finite strain beam model within the plane-strain context is developed in the present work. With some kinematic hypotheses, the finite deformation of the beam is analyzed, constitutive equations are deduced from the theory of finite elasticity, and by employing the standard variational method, the equilibrium equations and associated boundary conditions are derived. In the limit of infinitesimal strain, the new model degenerates to the classical extensible and shearable elastica model. The corresponding incremental equilibrium equations and associated boundary conditions are also obtained. Based on the new beam model, analytical solutions are given for simple deformation modes, including uniaxial tension, simple shear, pure bending, and buckling under an axial load. Furthermore, numerical solution procedures and results are presented for cantilevered beams and simply supported beams with immovable ends. The results are also compared with the previously developed finite strain Euler–Bernoulli beam model to demonstrate the significance of transverse shear effect for soft beams with a small length-to-thickness ratio. The developed beam model will contribute to the design and analysis of soft robots and soft devices.     

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.     

Free vibration analyses of axially loaded laminated composite beams based on higher-order shear deformation theory

The dynamic stiffness matrix method is introduced to solve exactly the free vibration and buckling problems of axially loaded laminated composite beams with arbitrary lay-ups. The Poisson effect, axial force, extensional deformation, shear deformation and rotary inertia are included in the mathematical formulation. The exact dynamic stiffness matrix is derived from the analytical solutions of the governing differential equations of the composite beams based on third-order shear deformation beam theory. The application of the present method is illustrated by two numerical examples, in which the effects of axial force and boundary condition on the natural frequencies, mode shapes and buckling loads are examined. Comparison of the current results to the existing solutions in the literature demonstrates the accuracy and effectiveness of the present method.     

含脱层损伤复合材料层合梁的冲击动力屈曲

针对含初始缺陷和脱层损伤的复合材料层合梁的轴向冲击动力屈曲问题进行了分析。基于Hamilton原理导出了考虑初始缺陷、轴向和横向惯性、横向剪切变形以及转动惯性影响时含脱层损伤复合材料梁的非线性动力屈曲控制方程;基于B-R准则,采用有限差分方法求解了受轴向冲击载荷作用下含脱层损伤复合材料梁的动力屈曲问题;讨论了冲击速度、初始几何缺陷、铺层角度以及脱层长度等因素对复合材料层合梁动力屈曲的影响。     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

含初缺陷裂纹损伤梁的冲击动力屈曲

由Hamilton原理导出考虑初始缺陷及横向剪切变形时裂纹梁的动力屈曲控制方程;应用断裂力学中常用的线弹簧模型将裂纹引入到屈曲控制方程中;基于B-R动力屈曲判断准则,采用数值方法求解了受轴向冲击载荷作用时裂纹梁的动力屈曲;对比讨论了不同冲击速度、初始几何缺陷大小以及分布形式等因素对梁冲击动力屈曲的影响。     

Instability of functionally graded micro-beams via micro-structure-dependent beam theory

This paper focuses on the buckling behaviors of a micro-scaled bi-directional functionally graded (FG) beam with a rectangular cross-section, which is now widely used in fabricating components of micro-nano-electro-mechanical systems (MEMS/NEMS) with a wide range of aspect ratios. Based on the modified couple stress theory and the principle of minimum potential energy, the governing equations and boundary conditions for a micro-structure-dependent beam theory are derived. The present beam theory incorporates different kinds of higher-order shear assumptions as well as the two familiar beam theories, namely, the Euler-Bernoulli and Timoshenko beam theories. A numerical solution procedure, based on a generalized differential quadrature method (GDQM), is used to calculate the results of the bi-directional FG beams. The effects of the two exponential FG indexes, the higher-order shear deformations, the length scale parameter, the geometric dimensions, and the different boundary conditions on the critical buckling loads are studied in detail, by assuming that Young’s modulus obeys an exponential distribution function in both length and thickness directions. To reach the desired critical buckling load, the appropriate exponential FG indexes and geometric shape of micro-beams can be designed according to the proposed theory.     

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 dynamic stability analysis of Euler–Bernoulli beam–columns with damping effects under thermal environment

In this study, a unified nonlinear dynamic buckling analysis for Euler–Bernoulli beam–columns subjected to constant loading rates is proposed with the incorporation of mercurial damping effects under thermal environment. Two generalized methods are developed which are competent to incorporate various beam geometries, material properties, boundary conditions, compression rates, and especially, the damping and thermal effects. The Galerkin–Force method is developed by implementing Galerkin method into force equilibrium equations. Then for solving differential equations, different buckled shape functions were introduced into force equilibrium equations in nonlinear dynamic buckling analysis. On the other hand, regarding the developed energy method, the governing partial differential equation for dynamic buckling of beams is also derived by meticulously implementing Hamilton’s principles into Lagrange’s equations. Consequently, the dynamic buckling analysis with damping effects under thermal environment can be adequately formulated as ordinary differential equations. The validity and accuracy of the results obtained by the two proposed methods are rigorously verified by the finite element method. Furthermore, comprehensive investigations on the structural dynamic buckling behavior in the presence of damping effects under thermal environment are conducted.     

热荷载作用下Timoshenko功能梯度夹层梁的静态响应

在精确考虑轴线伸长和一阶横向剪切变形的基础上建立了Timoshenko功能梯度夹层梁在热载荷作用下的几何非线性控制方程.采用打靶法数值求解所得强非线性边值问题,获得了两端固支功能梯度夹层梁在横向非均匀升温作用下的静态热过屈曲和热弯曲变形数值解.分析了功能梯度材料参数变化、不同表层厚度和升温参数对夹层梁弯曲变形、拉-弯耦...     

Geometrical nonlinear analysis of thin-walled composite beams using finite element method based on first order shear deformation theory

Based on a seven-degree-of-freedom shear deformable beam model, a geometrical nonlinear analysis of thin-walled composite beams with arbitrary lay-ups under various types of loads is presented. This model accounts for all the structural coupling coming from both material anisotropy and geometric nonlinearity. The general nonlinear governing equations are derived and solved by means of an incremental Newton–Raphson method. A displacement-based one-dimensional finite element model that accounts for the geometric nonlinearity in the von Kármán sense is developed to solve the problem. Numerical results are obtained for thin-walled composite beam under vertical load to investigate the effects of fiber orientation, geometric nonlinearity, and shear deformation on the axial–flexural–torsional response.     

CLASSICAL AND HOMOGENIZED EXPRESSIONS FOR BUCKLING SOLUTIONS OF FUNCTIONALLY GRADED MATERIAL LEVINSON BEAMS

The relationship between the critical buckling loads of functionally graded material(FGM) Levinson beams(LBs) and those of the corresponding homogeneous Euler-Bernoulli beams(HEBBs) is investigated. Properties of the beam are assumed to vary continuously in the depth direction. The governing equations of the FGM beam are derived based on the Levinson beam theory, in which a quadratic variation of the transverse shear strain through the depth is included.By eliminating the axial displacement as well as the rotational angle in the governing equations,an ordinary differential equation in terms of the deflection of the FGM LBs is derived, the form of which is the same as that of HEBBs except for the definition of the load parameter. By solving the eigenvalue problem of ordinary differential equations under different boundary conditions clamped(C), simply-supported(S), roller(R) and free(F) edges combined, a uniform analytical formulation of buckling loads of FGM LBs with S-S, C-C, C-F, C-R and S-R edges is presented for those of HEBBs with the same boundary conditions. For the C-S beam the above-mentioned equation does not hold. Instead, a transcendental equation is derived to find the critical buckling load for the FGM LB which is similar to that for HEBB with the same ends. The significance of this work lies in that the solution of the critical buckling load of a FGM LB can be reduced to that of the HEBB and calculation of three constants whose values only depend upon the throughthe-depth gradient of the material properties and the geometry of the beam. So, a homogeneous and classical expression for the buckling solution of FGM LBs is accomplished.     

基于修正偶应力和高阶剪切变形理论的变截面微梁的自由振动

基于修正偶应力和高阶剪切理论建立了仅含有一个尺度参数的Reddy变截面微梁的自由振动模型,研究了变截面微梁自由振动问题的尺度效应和横向剪切变形对自振频率计算的影响。基于哈密顿原理推导了动力学方程与边界条件,并采用微分求积法求解了各种边界条件下的自振频率。算例结果表明,基于偶应力理论预测的变截面微梁的自振频率均大于经典梁理论的预测结果,即捕捉到了尺度效应。另外,梁的几何尺寸与尺度参数越接近,尺度效应就越明显,而梁的长细比越小,横向剪切变形对自振频率的影响就越明显。     

A nonlinear planar beam formulation with stretch and shear deformations under end forces and moments

A new nonlinear planar beam formulation with stretch and shear deformations is developed in this work to study equilibria of a beam under arbitrary end forces and moments. The slope angle and stretch strain of the centroid line, and shear strain of cross-sections, are chosen as dependent variables in this formulation, and end forces and moments can be either prescribed or resultant forces and moments due to constraints. Static equations of equilibria are derived from the principle of virtual work, which consist of one second-order ordinary differential equation and two algebraic equations. These equations are discretized using the finite difference method, and equilibria of the beam can be accurately calculated. For practical, geometrically nonlinear beam problems, stretch and shear strains are usually small, and a good approximate solution of the equations can be derived from the solution of the corresponding Euler–Bernoulli beam problem. The bending deformation of the beam is the only important one in a slender beam, and stretch and shear strains can be derived from it, which give a theoretical validation of the accuracy and applicability of the nonlinear Euler–Bernoulli beam formulation. Relations between end forces and moments and relative displacements of two ends of the beam can be easily calculated. This formulation is powerful in the study of buckling of beams with various boundary conditions under compression, and can be used to calculate post-buckling equilibria of beams. Higher-order buckling modes of a long slender beam that have complex configurations are also studied using this formulation.     

Free vibration of axially loaded thin-walled composite Timoshenko beams

Based on shear-deformable beam theory, free vibration of thin-walled composite Timoshenko beams with arbitrary layups under a constant axial force is presented. This model accounts for all the structural coupling coming from material anisotropy. Governing equations for flexural-torsional-shearing coupled vibrations are derived from Hamilton’s principle. The resulting coupling is referred to as sixfold coupled vibrations. A displacement-based one-dimensional finite element model is developed to solve the problem. Numerical results are obtained for thin-walled composite beams to investigate the effects of shear deformation, axial force, fiber angle, modulus ratio on the natural frequencies, corresponding vibration mode shapes and load–frequency interaction curves.