Buckling load of seismic isolators affected by flexibility of reinforcement

Using the thinner steel reinforcing plates in the elastomeric multilayer isolators could reduce the weight of the isolators but would have a large effect on the buckling load of an isolator, which cannot be analyzed by the Haringx theory, a traditional approach on the stability analysis of rubber bearings. The buckling load of the isolators, which includes the effect of the flexibility of the steel reinforcing plates, is analyzed by a beam theory in which shear deformation and warping of the cross-section are considered. Pressure distributions in the elastomeric layer bonded to flexible reinforcements under compression force, bending moment and warping moment are derived from an assumed displacement field, from which the warping-related parameters used in the beam theory are established. The thickness of the steel reinforcement in the isolators is determined from the buckling load of the isolators, which is solved from a cubic equation established by the beam theory.     

Non-linear buckling and postbuckling behavior of thin-walled beams considering shear deformation

The static stability of thin-walled composite beams, considering shear deformation and geometrical non-linear coupling, subjected to transverse external force has been investigated in this paper. The theory is formulated in the context of large displacements and rotations, through the adoption of a shear deformable displacement field (accounting for bending and warping shear) considering moderate bending rotations and large twist. This non-linear formulation is used for analyzing the prebuckling and postbuckling behavior of simply supported, cantilever and fixed-end beams subjected to different load condition. 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 show that the beam loses its stability through a stable symmetric bifurcation point and the postbuckling strength is in relation with the buckling load value. Classical predictions of lateral buckling are conservative when the prebuckling displacements are not negligible and the non-linear buckling analysis is required for reliable solutions. The analysis is supplemented by investigating the effects of the variation of load height parameter. In addition, the critical load values and postbuckling response obtained with the present beam model are compared with the results obtained with a shell finite element model (Abaqus).     

Lateral buckling of beams with shear deformations – A hyperelastic formulation

The equilibrium and buckling equations are derived for the lateral buckling of a prismatic straight beam. A consistent finite strain constitutive law is used, which is based on a hyperelastic model for an isotropic material. The kinematics of the cross-sectional deformations are based on a Timoshenko type beam displacement of the cross-sectional plane using Euler angles and two shear finite rotations coupled with warping taken normal to the displaced plane. Also derived are the second order approximations to the displacements, curvatures, twist and internal actions. The constitutive relationships for the internal actions reveal new coupling terms between the bending moments, torsion and bimoment, which are functions of the cross-sectional warping and shear deformations. New Wagner type nonlinear torsion terms are derived which are functions of the warping of the cross-sectional plane, and are coupled to the twisting and shear deformations of the cross-section. Solutions are determined for the lateral buckling of a prismatic monosymmetric beam under pure bending and the flexural–torsional buckling under axial compression. For the flexural–torsional buckling problem it is found that the Euler type column buckling formula is consistent with Haringx’s column buckling formula while the torsional buckling formula is different to conventional equations. The second variation of the total potential is also derived. The effects of shear deformations are explored by examining the non-dimensional lateral buckling equation for a simply supported beam.     

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

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

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.     

粘贴压电层功能梯度材料Timoshenko梁的热过屈曲分析

研究了上下表面粘贴压电层的功能梯度材料Timoshenko梁在升温及电场作用下的过屈曲行为。在精确考虑轴线伸长和一阶横向剪切变形的基础上,建立了压电功能梯度Timoshenko层合梁在热-电-机械载荷作用下的几何非线性控制方程。其中,假设功能梯度的材料性质沿厚度方向按照幂函数连续变化,压电层为各向同性均匀材料。采用打靶法数值求解所得强非线性边值问题,获得了在均匀电场和横向非均匀升温场内两端固定Timoshenko梁的静态非线性屈曲和过屈曲数值解。并给出了梁的变形随热、电载荷及材料梯度参数变化的特性曲线。结果表明,通过施加电压在压电层产生拉应力可以有效地提高梁的热屈曲临界载荷,延缓热过屈曲发生。由于材料在横向的非均匀性,即使在均匀升温和均匀电场作用下,也会产生拉-弯耦合效应。但是对于两端固定的压电-功能梯度材料梁,在横向非均匀升温下过屈曲变形仍然是分叉形的。     

Influence of shear deformation on the buckling of elastically supported beams

Summary The influence of shear deformation on the buckling behavior of a beam supported laterally by a Winkler elastic foundation is studied. A full investigation of the bifurcation points at which, under axial load, the beam becomes critical with respect to one or two simultaneous buckling modes is made. The configurations and stabilities of the equilibrium paths that bifurcate from the critical points are derived. From the results of theoretical analysis, it becomes evident that shear deformation has a considerable effect upon the equilibriums and stabilities of the post-buckling of the beam. The results for the Bernoulli-Euler beam can be obtained as a limiting case for those of the present beam by letting the shear stiffness tend to infinity.Supported by the National Natural Science Foundation of China     

基于剪切变形的矩形梁剪力滞求解方法

Timoshenko梁通过假设截面的剪切刚度和附加平均剪切转角变形的方式来近似修正初等梁中未考虑剪切变形能的问题,这与梁剪应力沿梁高变化的实际不符。本文基于材料力学剪应力计算式和相应的剪切变形理论,从剪切变形与梁的位移关系入手,导出矩形梁考虑剪切变形时的纵向位移沿梁高方向的函数关系式,证明该位移可分解为纯弯曲引起的位移和剪力引起的剪力滞翘曲位移之和。应用剪力滞广义坐标与广义力的概念,基于能量变分原理得到等截面梁剪力滞控制微分方程组及其通解形式。对均布荷载作用下矩形简支梁的算例分析表明,本文算法与弹性力学精确解对比,两者的应力和挠度剪力滞系数求解结果非常接近,本文算法有足够的精度,且比弹性力学简单。     

Non-uniform warping including the effects of torsion and shear forces. Part I: A general beam theory

This two-part contribution presents a beam theory with a non-uniform warping including the effects of torsion and shear forces, and valid for any homogeneous cross-section made of isotropic elastic material. Part I is devoted to the theoretical developments and part II discusses analytical and numerical results obtained for torsion and shear-bending of cantilever beams made of different kinds of cross-section. The theory is based on a kinematics assuming that the cross-section maintains its shape and including three independent warping parameters associated to the three warping functions corresponding to torsion and shear forces. Starting from this displacement model and using the principle of virtual work, the corresponding beam theory is derived. For this theory, closed-form results are obtained for the cross-sectional constants and the three-dimensional expressions of the normal and shear stresses. Comparison with classical beam theories is carried out and additional effects due to the non-uniformity of the warping are highlighted. In particular, the contributions of primary and secondary internal forces and the effect of the non-symmetry of the cross-section on the structural behavior of the beam are specified. Simplified versions of this theory, wherein the number of degrees of freedom is reduced, are also presented. The analytical and numerical analyzes presented in part II give responses on the quality of this non-uniform beam theory and indicate also when its simplified versions could be applied.     

Application of the flexibility influence function method in the dynamic analysis of composite beams

The flexibility influence function technique is validated as a method for calculating the displacements and the rotations of a laminated beam subjected to a dynamic load, using the first-order shear deformation laminate theory and comparing the results with those obtained by modal analysis and two finite element models (one-dimensional and three-dimensional). The movements (displacements and rotations) were calculated from a single-span beam subjected to a time-variable load with four boundary conditions: clamped-clamped, hinged-hinged, clamped-free, clamped-hinged. A carbon/epoxy cross-ply laminated beam was selected to avoid bending-torsion coupling. The maximum movements calculated by the flexibility influence function method differs very little from those calculated with the other two models accounted for by the first-order shear deformation laminate theory: modal analysis and the one-dimensional numerical model. The differences in the rotations between the three-dimensional numerical model and the flexibility influence function method are slightly bigger, and could be due to the warping of the cross-section of the beam, which is not included in the first-order shear deformation laminate theory.     

Finite strain––beam theory

An appropriate strain energy density for an isotropic hyperelastic Hookean material is proposed for finite strain from which a constitutive relationship is derived and applied to problems involving beam theory approximations. The physical Lagrangian stress normal to the surfaces of a element in the deformed state is a function of the normal component of stretch while the shear is a function of the shear component of stretch. This paper attempts to make a contribution to the controversy about who is correct, Engesser or Haringx with regard to the buckling formula for a linear elastic straight prismatic column with Timoshenko beam-type shear deformations. The derived buckling formula for a straight prismatic column including shear and axial deformations agrees with Haringx’s formula. Elastica-type equations are also derived for a three-dimensional Timoshenko beam with warping excluded. When the formulation is applied to the problem of pure torsion of a cylinder no second-order axial shortening associated with the Wagner effect is predicted which differs from conventional beam theory. When warping is included, axial shortening is predicted but the formula differs from conventional beam theory.     

基于新型广义位移的箱形梁扭转单元及应用

为了改进变截面连续箱梁桥的扭转分析理论,将截面总扭转角分解为自由翘曲扭转角和约束剪切扭转角,选取自由翘曲转角扭率作为广义位移,提出一个2节点8自由度的扭转梁段单元。从约束扭转控制微分方程出发,推导单元刚度矩阵及等效节点荷载列阵。引入应力增大系数,以反映约束扭转对初等梁应力的增大效应。数值算例验证了本文梁段单元的可靠性。最后对一个三跨变截面连续箱梁桥进行分析,结果表明,双力矩影响线与弯矩影响线较为类似,按双力矩影响线进行最不利荷载加载时最大应力值偏小;应力增大系数在集中荷载作用截面出现极值,均发生在腹板与顶板交点处;利用偏载放大系数来考虑扭转附加效应时,不宜考虑弯曲正应力较小及翘曲正应力出现极值的梁段区域。     

考虑约束扭转的薄壁梁单元刚度矩阵

推导了薄壁空间梁单元刚度矩阵 ,考虑了双向弯曲及截面约束扭转对杆件轴向变形的影响 ;计算了截面的翘曲变形 ,以及二次剪应力对翘曲变形的影响 ,可适用于任意截面 (包括开口、闭口和混合剖面 )的薄壁杆件。计算结果表明 ,考虑约束扭转的薄壁梁单元刚度矩阵有相当好的精确度 ,可以用于薄壁杆件的静动力分析。     

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.     

Solution of free vibration equations of semi-rigid connected Reddy–Bickford beams resting on elastic soil using the differential transform method

The literature regarding the free vibration analysis of Bernoulli–Euler and Timoshenko beams under various supporting conditions is plenty, but the free vibration analysis of Reddy–Bickford beams with variable cross-section on elastic soil with/without axial force effect using the Differential Transform Method (DTM) has not been investigated by any of the studies in open literature so far. In this study, the free vibration analysis of axially loaded and semi-rigid connected Reddy–Bickford beam with variable cross-section on elastic soil is carried out by using DTM. The model has six degrees of freedom at the two ends, one transverse displacement and two rotations, and the end forces are a shear force and two end moments in this study. The governing differential equations of motion of the rectangular beam in free vibration are derived using Hamilton’s principle and considering rotatory inertia. Parameters for the relative stiffness, stiffness ratio and nondimensionalized multiplication factor for the axial compressive force are incorporated into the equations of motion in order to investigate their effects on the natural frequencies. At first, the terms are found directly from the analytical solutions of the differential equations that describe the deformations of the cross-section according to the high-order theory. After the analytical solution, an efficient and easy mathematical technique called DTM is used to solve the governing differential equations of the motion. The calculated natural frequencies of semi-rigid connected Reddy–Bickford beam with variable cross-section on elastic soil using DTM are tabulated in several tables and figures and are compared with the results of the analytical solution where a very good agreement is observed.     

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.     

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.     

In-plane buckling of prismatic funicular arches with shear deformations

The in-plane buckling behavior of funicular arches is investigated numerically in this paper. A finite strain Timoshenko beam-type formulation that incorporates shear deformations is developed for generic funicular arches. The elastic constitutive relationships for the internal beam actions are based on a hyperelastic constitutive model, and the funicular arch equilibrium equations are derived. The problems of a submerged arch under hydrostatic pressure, a parabolic arch under gravity load and a catenary arch loaded by overburden are investigated. Buckling solutions are derived for the parabolic and catenary arch. Subsequent investigation addresses the effects of axial deformation prior to buckling and shear deformation during buckling. An approximate buckling solution is then obtained based on the maximum axial force in the arch. The obtained buckling solutions are compared with the numerical solutions of Dinnik (Stability of arches, 1946) [1] and the finite element package ANSYS. The effects of shear deformation are also evaluated.     

A new efficient method to evaluate exact stiffness and mass matrices of non-uniform beams resting on an elastic foundation

In this paper, a new efficient method to evaluate the exact stiffness and mass matrices of a non-uniform Bernoulli–Euler beam resting on an elastic Winkler foundation is presented. The non-uniformity may result from variable cross-section and/or from inhomogeneous linearly elastic material. It is assumed that there is no abrupt variation in the cross-section of the beam so that the Euler–Bernoulli theory is valid. The method is based on the integration of the exact shape functions which are derived from the solution of the axial deformation problem of a non-uniform bar and the bending problem of a non-uniform beam which are both formulated in terms of the two displacement components. The governing differential equations are uncoupled with variable coefficients and are solved within the framework of the analog equation concept. According to this, the two differential equations with variable coefficients are replaced by two linear ones pertaining to the axial and transverse deformation of a substitute beam with unit axial and bending stiffness, respectively, under ideal load distributions. The key point of the method is the evaluation of the two ideal loads which in this work is achieved by approximating them by two polynomials. More specifically, the axial ideal load is approximated by a linear polynomial while the transverse one by a cubic polynomial. The numerical implementation of the method is simple, and the results are compared favorably to those obtained by exact solutions available in literature.     

A new super convergent thin walled composite beam element for analysis of box beam structures

In this paper, a new composite thin wall beam element of arbitrary cross-section with open or closed contour is developed. The formulation incorporates the effect of elastic coupling, restrained warping, transverse shear deformation associated with thin walled composite structures. A first order shear deformation theory is considered with the beam deformation expressed in terms of axial, spanwise and chordwise bending, corresponding shears and twist. The formulated locking free element uses higher order interpolating polynomial obtained by solving static part of the coupled governing differential equations. The formulated element has super convergent properties as it gives the exact elemental stiffness matrix. Static and free vibration analyses are performed for various beam configuration and compared with experimental and numerical results available in current literature. Good correlation is observed in all cases with extremely small system size. The formulated element is used to study the wave propagation behavior in box beams subjected to high frequency loading such as impact. Simultaneous existence of various propagating modes are graphically captured. Here the effect of transverse shear on wave propagation characteristics in axial and transverse directions are investigated for different ply layup sequences.