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

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

Non-linear elastic non-uniform torsion of bars of arbitrary cross-section by BEM

In this paper an elastic non-uniform torsion analysis of simply or multiply connected cylindrical bars of arbitrary cross-section taking into account the effect of geometric non-linearity is presented employing the boundary-element(BE) method. The torque-rotation relationship is computed based on the finite-displacement (finite-rotation) theory, that is the transverse displacement components are expressed so as to be valid for large rotations and the longitudinal normal strain includes the second-order geometric non-linear term often described as the “Wagner strain”. The proposed formulation does not stand on the assumption of a thin-walled structure and therefore the cross-section's torsional rigidity is evaluated exactly without using the so-called Saint-Venant's torsional constant. The torsional rigidity of the cross-section is evaluated directly employing the primary warping function of the cross-section depending on its shape. Three boundary-value problems with respect to the variable along the beam axis angle of twist, to the primary and to the secondary warping functions are formulated. The first one, employing the Analog Equation Method (a BEM-based method), yields a system of non-linear equations from which the angle of twist is computed by an iterative process. The remaining two problems are solved employing a pure BE method. Numerical results are presented to illustrate the method and demonstrate its efficiency and accuracy. The developed procedure retains most of the advantages of a BEM solution over a pure domain discretization method, although it requires domain discretization.     

Shear stresses in elastic beams: an intrinsic approach

The theory of non-uniform flexure and torsion of Saint-Venant's beam with arbitrary multiply connected cross section is revisited in a coordinate-free form to provide a computationally convenient context. Numerical implementations, by Matlab, are performed to evaluate the maximum elastic shear stresses in beams with rectangular cross sections for different Poisson's ratios. The deviations between the maximum and mean stresses are then diagrammed to adjust the results provided by Jourawski's method.     

On the torsion of a cylinder with several cracks

In this paper, based on paper [1], the analytic expression of the torsion function for a cylinder containing arbitrary oriented cracks is obtained. The problem is reduced to solve a system of singular integral equations for the unknown dislocation density functions. Using the numerical method of the singular integral equations^[2,7] the torsional rigidities and stress intensity factors are evaluated for several multicracked cylinders. Next, the creak-cutting method^[5] is firstly extended to lve the torsion problem for a rectangular prism. The numerical results show that the method presented here is successful. Projects Supported by the Science Fund of the Chinese Academy of Sciences.     

采用有限单元法计算梁任意形状截面特性

采用将梁截面离散化的方式,用数值积分计算截面的几何特性,并根据梁剪切变形和扭转理论,利用变分原理建立截面的有限元法方程,求解任意形状截面的扭转常数、剪切中心以及剪切面积修正系数等特性.本方法适用于各种形式的截面,具有计算精度高及适应性强的特点.根据上述理论编制了相应程序,按照不同的单元划分方式,分别计算出矩形截面截面特性,与理论解进行比较;又对舟山市定海长峙至岙山预应力混凝土连续箱梁截面进行了计算,并与Ansys结果进行比较,均证明采用本文的计算方法能得到满意的结果,且该方法适用于各种形状的截面形式.     

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.      

Effective properties of a composite beam subjected to torsion

The problem of finding the effective characteristics of a rectilinear beam under pure torsion is considered. The problem can be reduced to determining the torsional stress function from the solution of a boundary-value problem in a cross section of the beam for a partial differential equation with variable coefficients. Two special boundary-value problems are formulated to find the effective characteristics. It is shown that the effective coefficients are reciprocal in the case of torsion of a layer with nonuniform thickness. In the two-dimensional case, the problem is solved by a finite element method. The cases of a square beam with single and multiple inclusions are discussed. The dependence of the effective characteristics on the inclusion volume fraction is analyzed.     

Non-uniform warping including the effects of torsion and shear forces. Part II: Analytical and numerical applications

This two-part contribution presents a beam theory (BT) with a non-uniform warping (NUW) including the effects of torsion, and shear forces and valid for any homogeneous cross-section made of isotropic elastic material. In part I, the governing equations of the NUW-BT has been established and simplified-NUW-BT versions has been deduced, wherein the number of degrees of freedom is reduced. In this part II, these theories are used to analyze, for a representative set of cross-sections (CS) (solid-CS and thin-walled open/closed-CS, bi-symmetric or not), the elastic behavior of cantilever beams subjected to torsion or shear-bending. For bi-symmetrical-CS, torsion and shear-bending are analyzed separately: analytical and numerical results are given for the distributions along the beam axis of the cross-sectional displacements and stresses, for the NUW-BT and its simplified versions. Numerical results are also given for the three-dimensional stress distributions close to the embedded section: the stress predictions of the NUW-BT are compared to those obtained by three-dimensional finite elements computations. It can be drawn from all these results indications that can help to decide when the simplified theories may be applied, and hence when the warping parameters may be reduced. As specified in NUW-BT, torsion and bending are coupled for non-symmetrical-CS, even if the bending moments refer to the centroid while the torsional moment refers to the shear center. To illustrate this coupling effect, the particular example of the channel-CS presented in Kim and Kim [Kim, N.-I., Kim, M.-Y., 2005. Exact dynamic/static stiffness matrices of non-symmetric thin-walled beams considering coupled shear deformation effects. Thin-Walled Structures 43, 701–734.] is analyzed and the results are compared.     

A Universal Expression for Bounds on the Torsional Rigidity of Cylindrical Shafts Containing Fibers with Arbitrary Transverse Cross-Sections

We derive upper and lower bounds for the torsional rigidity of host shafts containing a number of cylindrical fibers. The transverse cross-sections of the host shaft and the fibers are simply connected, but could be arbitrary in shape. Utilizing the fact that the torsion solution of a homogeneous host shaft with simply connected cross-section can be known, we propose a method to construct statically and kinematically admissible fields interior to the matrix and to the fibers. Previous developments on bounding the torsional rigidity of composite shaft so far are confined to circular fibers. Here we try to simulate fibers with non-circular cross-section and incorporate the interactions of the cross-sectional shapes of the host shaft and the fibers at the same time. Proceeding from extremal principles of elasticity, together with propositions of some domain integration procedures, we provide a universal expression for bounds on the torsional rigidity of the composite shaft. The exact expressions depend on the constituent information of the fibers and the host shaft, which could offer useful information to tailor the shape and the arrangement of the constituents to achieve an optimal value.     

开口薄壁梁的扭转理论与应用

以Vlasov薄壁构件理论为基础, 推导了开口薄壁构件一阶扭转理论. 该理论考虑了翘曲剪应力对截面转角的影响, 截面的转角分为自由翘曲转角和约束剪切转角, 在约束扭转中, St.Venant扭矩仅仅与自由翘曲转角有关, 而翘曲扭矩仅与约束剪切转角有关. 利用半逆解方法求出了约束扭转中薄壁构件的St.Venant扭矩表达公式; 依据能量方法, 建立了约束剪切转角和翘曲扭矩之间的关系, 并提出了翘曲剪切系数概念, 给出了一阶扭转理论的微分方程. 为了有效求解微分方程, 给出了求解微分方程的初参数法方程和相应的影响函数矩阵; 当St.Venant扭矩可以忽略时, 得到与一阶弯曲理论(Timoshenko梁理论)相似的一阶扭转理论简化形式. 最后利用算例证明了一阶扭转理论和简化理论的有效性.      

圆截面弹性细杆的平面振动

基于Kirchhoff理论讨论圆截面弹性细杆的平面振动．以杆中心线的Frenet坐标系为参考系建立动力学方程．杆作平面运动时，其扭转振动与弯曲振动解耦．讨论任意形状杆的扭转振动和轴向受压直杆在无扭转条件下的弯曲振动，证明直杆平衡的静态Lyapunov稳定性与欧拉稳定性条件为动态稳定性的必要条件．考虑轴向力和截面转动惯性效应的影响，导出弯曲振动的固有频率．     

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

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

A quaternion-based formulation of Euler–Bernoulli beam without singularity

This paper proposes a singularity-free beam element with Euler–Bernoulli assumption, i.e., the cross section remains rigid and perpendicular to the tangent of the centerline during deformation. Each node of this two-nodal beam element has eight nodal coordinates, including three global positions and one normal strain to describe the rigid translation and flexible deformation of the centerline, respectively, four Euler parameters or quaternion to represent the attitude of cross section. Adopting quaternion instead of Eulerian angles as nodal variables avoids the traditionally encountered singularity problem. The rigid cross section assumption is automatically satisfied. To guarantee the perpendicularity of cross section to the deformed neutral axes, the position and orientation coordinates are coupled interpolated by a special method developed here. The proposed beam element allows arbitrary spatial rigid motion, and large bending, extension, and torsion deformation. The resulting governing equations include normalization constraint equations for each quaternion of the beam nodes, and can be directly solved by the available differential algebraic equation (DAE) solvers. Finally, several numerical examples are presented to verify the large deformation, natural frequencies and dynamic behavior of the proposed beam element.     

考虑剪应变的弹性薄壁杆件的振动理论

对于开口截面彈性薄壁杆件的振动問題,符拉索夫在1940年提出的理論,是基于下列两个基本假設:(1)截面具有剛硬不可变形的周綫;(2)中曲面无剪应变。其后,有些学者也曾提出过相同的理論(例如.R.海里格).由于忽略了中曲面的剪应变,我們可以估計得到,符拉索夫的理論只适用于研究低频率的固有振动和強迫振动,而不适用于研究高频瘁的振动。由于行动載荷和冲击载荷必然引起各个固有频率的振动,因此符拉索夫的理論也不适用于研究这些載荷作用下的振动問題。为了要更好地研究这些問題,必須要考慮剪应变的影响。     

Non-linear flexural–torsional dynamic analysis of beams of arbitrary cross section by BEM

In this paper, a boundary element method is developed for the non-linear flexural–torsional dynamic analysis of beams of arbitrary, simply or multiply connected, constant cross section, undergoing moderately large deflections and twisting rotations under general boundary conditions, taking into account the effects of rotary and warping inertia. The beam is subjected to the combined action of arbitrarily distributed or concentrated transverse loading in both directions as well as to twisting and/or axial loading. Four boundary value problems are formulated with respect to the transverse displacements, to the axial displacement and to the angle of twist and solved using the Analog Equation Method, a BEM based method. Application of the boundary element technique leads to a system of non-linear coupled Differential–Algebraic Equations (DAE) of motion, which is solved iteratively using the Petzold–Gear Backward Differentiation Formula (BDF), a linear multistep method for differential equations coupled to algebraic equations. The geometric, inertia, torsion and warping constants are evaluated employing the Boundary Element Method. The proposed model takes into account, both the Wagner's coefficients and the shortening effect. Numerical examples are worked out to illustrate the efficiency, wherever possible the accuracy, the range of applications of the developed method as well as the influence of the non-linear effects to the response of the beam.     

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.     

The solution of integral equations with strongly singular kernels applied to the torsion of cracked circular cylinder

In this paper,the functions of warping displacement interruption defined on the crack lines are taken for the fundamental unknown functions.The torsion problem of cracked circular cylinder is reduced to solving a system of integral equations with strongly singular kernels.Using the numerical method of these equations,the torsional rigidities and the stress intensity factors are calculated to solve the torsion problem of circular cylinder with star-type and other different types of cracks.The numerical results are satisfactory.     

用伪应力函数法求解幂硬化材料柱体的扭转问题

本文引用伪应力函数使得幂硬化材料的任意形状等截面柱体和变直径圆截面柱体扭转问题的定解方程具有弹性柱体扭转问题的相应形式,从而可用类似于求解弹性柱体扭转的方法或直接利用已知的弹性解答求解对应的幂硬化材料柱体的扭转问题,本文用这种方法求得了幂硬化材料椭圆截面柱体及含球形空腔的圆轴扭转问题的解析解。     

用伪应力函数法求解幂硬化材料柱体的扭转问题

本文引用伪应力函数使得幂硬化材料的任意形状等截面柱体和变直径圆截面柱体扭转问题的定解方程具有弹性柱体扭转问题的相应形式,从而可用类似于求解弹性柱体扭转的方法或直接利用已知的弹性解答求解对应的幂硬化材料柱体的扭转问题,本文用这种方法求得了幂硬化材料椭圆截面柱体及含球形空腔的圆轴扭转问题的解析解。