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.     

用高阶剪切变形理论研究双模量梁的弯曲

利用高阶剪切变形理论研究了双模量梁的弯曲变形问题，推导出了双模量梁的挠曲线方程及弯曲正应力公式．讨论分析了翘曲函数的指数n对挠度、正应力的影响．研究结果表明：拉压弹性模量的差异对梁的弯曲应力有较大影响．把高阶剪切变形理论的计算结果与弹性理论计算结果进行比较，可知该方法计算精度非常高．     

开口厚壁截面短构件的约束扭转1）

为了分析开口厚壁截面短构件的约束扭转问题,采用统一分析梁模型与有限节线法,对T形和L形厚壁截面短构件约束扭转时横截面的翘曲和应力分布情况等问题进行了分析研究.算例计算结果表明：开口厚壁截面短构件存在与其横截面形心位置不一致的扭转（弯曲）中心,构件在不过扭转中心的外力作用下会产生弯扭耦合变形,其横截面将产生不均匀翘曲,横截面上的翘曲正应力和扭转剪应力均呈非线性分布.     

Justification of the Asymptotic Expansion Method for Homogeneous Isotropic Beams by Comparison with De Saint-Venant’s Solutions

The formal asymptotic expansion method is an attractive mean to derive simplified models for problems exhibiting a small parameter, such as the elastic analysis of beam-like structures. Usually this method is rigorously justified using convergence theorems Yu and Hodges, 2004. In this paper it is illustrated how the Saint-Venant’s solution naturally arises from the lowest order terms of an asymptotic expansion of the elastic state for the case of homogeneous isotropic beams. It is also highlighted that the Saint-Venant solutions corresponding to pure traction, bending and torsion involve the solution of the first-order microscopic problems, while for the simple bending problem, the solution of the second-order microscopic problems is needed. The second-order problems provide therefore a way to characterize the transverse shear behavior and the cross-sectional warping of the beam.     

Shear stress analysis in engineering beams using deplanation field of special 2-D finite elements

This paper deals with the 2-D finite element shear stress analysis in beams, loaded by bending with shear and St. Venant’s torsion. The properties of these finite elements, like stiffness matrices as well as load vectors, are derived on the basis of their axial nodal displacements, e.g. by warping field. Proposed finite elements enable stress analysis independently of both cross-sectional member shape and material properties. Stiffness matrices and load vectors are derived for several finite element types. Material is assumed to be isotropic and linear elastic. For justification of the proposed stress analysis procedure, some examples are presented.     

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

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

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.     

Bending and free vibrational analysis of bi-directional functionally graded beams with circular cross-section

The bending and free vibrational behaviors of functionally graded (FG) cylindrical beams with radially and axially varying material inhomogeneities are investigated. Based on a high-order cylindrical beam model, where the shear deformation and rotary inertia are both considered, the two coupled governing differential motion equations for the deflection and rotation are established. The analytical bending solutions for various boundary conditions are derived. In the vibrational analysis of FG cylindrical beams, the two governing equations are firstly changed to a single equation by means of an auxiliary function, and then the vibration mode is expanded into shifted Chebyshev polynomials. Numerical examples are given to investigate the effects of the material gradient indices on the deflections, the stress distributions, and the eigenfrequencies of the cylindrical beams, respectively. By comparing the obtained numerical results with those obtained by the three-dimensional (3D) elasticity theory and the Timoshenko beam theory, the effectiveness of the present approach is verified.     

复合材料叠层梁和金属梁的固有振动特性

对根据三种梁理论得到的金属梁和复合材料叠层梁的固有振动特性进行了对比性的研究对常用的三种梁理论在弹性碰撞分析中的应用进行了分析和比较     

A nonlinear composite beam theory

Presented here is a general theory for the three-dimensional nonlinear dynamics of elastic anisotropic initially straight beams undergoing moderate displacements and rotations. The theory fully accounts for geometric nonlinearities (large rotations and displacements) by using local stress and strain measures and an exact coordinate transformation, which result in nonlinear curvature and strain-displacement expressions that contain the von Karman strains as a special case. Extensionality is included in the formulation, and transverse shear deformations are accounted for by using a third-order theory. Six third-order nonlinear partial-differential equations are derived for describing one extension, two bending, one torsion, and two shearing vibrations of composite beams. They show that laminated beams display linear elastic and nonlinear geometric couplings among all motions. The theory contains, as special cases, the Euler-Bernoulli theory, Timoshenko's beam theory, the third-order shear theory, and the von Karman type nonlinear theory.     

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.     

A new finite element of spatial thin-walled beams

Based on the theories of Timoshenko＇s beams and Vlasov＇s thin-walled members, a new spatial thin-walled beam element with an interior node is developed. By independently interpolating bending angles and warp, factors such as transverse shear deformation, torsional shear deformation and their Coupling, coupling of flexure and torsion, and second shear stress are considered. According to the generalized variational theory of Hellinger-Reissner, the element stiffness matrix is derived. Examples show that the developed model is accurate and can be applied in the finite element analysis of thinwalled structures.     

Analytical test functions for free vibration analysis of rotating non-homogeneous Timoshenko beams

In this paper, the governing equations for free vibration of a non-homogeneous rotating Timoshenko beam, having uniform cross-section, is studied using an inverse problem approach, for both cantilever and pinned-free boundary conditions. The bending displacement and the rotation due to bending are assumed to be simple polynomials which satisfy all four boundary conditions. It is found that for certain polynomial variations of the material mass density, elastic modulus and shear modulus, along the length of the beam, the assumed polynomials serve as simple closed form solutions to the coupled second order governing differential equations with variable coefficients. It is found that there are an infinite number of analytical polynomial functions possible for material mass density, shear modulus and elastic modulus distributions, which share the same frequency and mode shape for a particular mode. The derived results are intended to serve as benchmark solutions for testing approximate or numerical methods used for the vibration analysis of rotating non-homogeneous Timoshenko beams.     

Homogenized and classical expressions for static bending solutions for functionally graded material Levinson beams

The exact relationship between the bending solutions of functionally graded material (FGM) beams based on the Levinson beam theory and those of the corresponding homogenous beams based on the classical beam theory is presented for the material properties of the FGM beams changing continuously in the thickness direction. The deflection, the rotational angle, the bending moment, and the shear force of FGM Levinson beams (FGMLBs) are given analytically in terms of the deflection of the reference homogenous Euler-Bernoulli beams (HEBBs) with the same loading, geometry, and end supports. Consequently, the solution of the bending of non-homogenous Levinson beams can be simplified to the calculation of transition coefficients, which can be easily determined by variation of the gradient of material properties and the geometry of beams. This is because the classical beam theory solutions of homogenous beams can be easily determined or are available in the textbook of material strength under a variety of boundary conditions. As examples, for different end constraints, particular solutions are given for the FGMLBs under specified loadings to illustrate validity of this approach. These analytical solutions can be used as benchmarks to check numerical results in the investigation of static bending of FGM beams based on higher-order shear deformation theories.     

Coupled bending and torsional vibration of nonsymmetrical axially loaded thin-walled Bernoulli–Euler beams

The dynamic transfer matrix is formulated for a straight uniform and axially loaded thin-walled Bernoulli–Euler beam element whose elastic and inertia axes are not coincident by directly solving the governing differential equations of motion of the beam element. Bernoulli–Euler beam theory is used, and the cross section of the beam does not have any symmetrical axes. The bending vibrations in two perpendicular directions are coupled with torsional vibration and the effect of warping stiffness is included. The dynamic transfer matrix method is used for calculation of exact natural frequencies and mode shapes of the nonsymmetrical thin-walled beams. Numerical results are given for a specific example of thin-walled beam under a variety of end conditions, and exact numerical solutions are tabulated for natural frequencies and solutions calculated by the other method are also tabulated for comparison. The effects of axial force and warping stiffness are also discussed.     

用材料力学公式计算任意截面短梁正应力时的修正项研究

对材料力学中梁的弯曲应力公式增加一修正项,以反映短梁弯剪翘曲变形对应力分布的影响。提出一种根据短梁横截面边界形状及艾瑞应力函数求解应力修正项的方法,应用弹性力学空间问题的一般理论,通过应力平衡方程、应变相容方程及应力边界条件,建立了关于任意截面短梁的应力修正项及剪应力的基本方程。在所建立的基本方程基础上,导出了矩形截面和圆形截面短梁修正应力的具体计算公式,该修正应力与均布荷载大小及弹性模量与剪切模量之比均成正比,但与截面惯性矩成反比。数值算例表明,本文方法计算的应力与通用有限元软件ANSYS计算的结果吻合良好,从而验证了本文方法及其基本公式的正确性。     

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.     

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.     

Saint-Venant’s Problem for Compound Piezoelectric Beams

The solution of the Saint-Venant’s Problem for a slender compound piezoelectric beam presented in this paper generalizes the recent solution by the authors and E. Harash (J. Appl. Mech. 11:1–10, 2007) for a homogeneous piezoelectric beam and the solution for a compound elastic beam developed by O. Rand and the first author (Analytical Methods in Anisotropic Elasticity with Symbolic Computational Tools, Birkhauser, Boston, 2005). Justification for this approximation emerges from the St. Venant’s Principle. The stress, strain and (electrical) displacement components (“solution hypothesis”) are presented as a set of initially assumed expressions involving twelve tip loading parameters, six unknown weight coefficients, and three pairs of torsion/bending functions of two variables. Each pair of functions satisfies the so-called coupled non-homogeneous Neumann problem (CNNP) in the cross-sectional domain. The work develops concepts of the torsion/bending functions, the torsional rigidity and piezoelectric shear center, the tip coupling matrix, for a compound piezoelectric beam. Examples of exact and approximate solutions for rectangular laminated beams made of transtropic materials are presented.