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

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

Using the logarithmic strain measure for solving torsion problems

The problems of free and constrained torsion of a rod of solid circular cross-section are solved numerically using a tensor linear constitutive relation written in terms of the energy compatible Cauchy stress and Hencky logarithmic strain tensors. The only function used to determine the properties of the isotropic incompressible material of the rod is a power-law function that approximates the shear diagram and corresponds to an elastoplastic material with power-law hardening. The solution obtained shows that, despite the tensorial linearity of the state law, the use of the logarithmic strain measure allows one to describe qualitatively the effect of significant elongation of the rod in free torsion (the Poynting effect) as well as the arising normal longitudinal, radial, and circumferential stresses, whose values are commensurable, at large deformations, with the maximum tangential stresses in the cross-section. Computational dependences of the torsional moment on the angle of twist in free and constrained torsion are obtained. These dependences are found to be significantly different from each other; the limitmoment and the correspondingmaximum angle of twist for free torsion are found to be considerably lower than those for constrained torsion. It follows that the shear strength, which is traditionally calculated from the maximum torsional moment, becomes indeterminate. For constrained torsion, the dependence of the longitudinal compressive force on the angle of twist is obtained.     

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

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

A new beam element for analyzing geometrical and physical nonlinearity

Based on Timoshenko＇s beam theory and Vlasov＇s thin-walled member theory, a new model of spatial thin-walled beam element is developed for analyzing geometrical and physical nonlinearity, which incorporates an interior node and independent interpolations of bending angles and warp and takes diversified factors into consideration, such as traverse shear deformation, torsional shear deformation and their coupling, coupling of flexure and torsion, and the second shear stress. The geometrical nonlinear strain is formulated in updated Lagarange （UL） and the corresponding stiffness matrix is derived. The perfectly plastic model is used to account for physical nonlinearity, and the yield rule of von Mises and incremental relationship of Prandtle-Reuss are adopted. Elastoplastic stiffness matrix is obtained by numerical integration based on the finite segment method, and a finite element program is compiled. Numerical examples manifest that the proposed model is accurate and feasible in the analysis of thin-walled structures.     

Stress bounds for bars in torsion

Upper bounds for the maximum shear stress in the St. Venant torsion problem are derived with the aid of the theory of subharmonic functions. The main result is a bound that is determined in a simple manner by the magnitude of the applied twisting moment and two parameters peculiar to the cross section: the radius of the largest circle contained in it and the minimum curvature of the curve that bounds it.

---

Zusammenfassung Mit Hilfe subharmonischer Funktionen werden obere Grenzen in dem Torsionsproblem von St. Venant erhalten. Das Hauptergebnis ist eine Grenze, die auf einfache Weise vom Drehmoment und zwei nur vom Querschnitt abhängigen Parametern bestimmt ist, und zwar von dem Radius des grössten eingeschriebenen Kreises und von der Minimalkrümmung der Begrenzungskurve des Querschnitts.

     

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.     

St. Venant torsion of composite cylinders

Using a method first applied to some fluid dynamical two-phase flow problems [1], the authors derive expressions for the warping and stress functions describing the St. Venant torsion of some compound cylinders. The construction of the cylinder must have appropriate symmetries, and composites containing both two and four different elastic phases are considered. Effective shear moduli and torsional rigidities are defined, and some tables of numerical values given for various configurations. An analogous solution for the antiplane strain of a two-phase composite is also given.     

Torsin of Functionally Graded Isotropic Linearly Elastic Bars

The purpose of this research is to investigate the effects of material inhomogeneity on the torsional response of linearly elastic isotropic bars. The work is motivated by the recent research activity on functionally graded materials (FGMs), i.e. materials with spatially varying properties tailored to satisfy particular engineering applications. The classic approach to the torsion problem for a homogenous isotropic bar of arbitrary simply-connected cross-section in terms of the Prandtl stress function is generalized to the inhomogeneous case. The special case of a circular rod with shear modulus depending on the radial coordinate only is examined. It is shown that the maximum shear stress does not, in general, occur on the boundary of the rod, in contrast to the situation for the homogeneous problem. It is shown that the material inhomogeneity may increase or decrease the torsional rigidity compared to that for the homogeneous rod. Optimal upper and lower bounds for the torsional rigidity for nonhomogeneous bars of arbitrary cross-section are established. A new formulation of the basic boundary-value problem is given. The results are illustrated using specific material models used in the literature on functionally graded elastic materials. This revised version was published online in August 2006 with corrections to the Cover Date.     

Eigenanalysis and continuum modelling of pre-twisted repetitive beam-like structures

A repetitive pin-jointed, pre-twisted structure is analysed using a state variable transfer matrix technique. Within a global coordinate system the transfer matrix is periodic, but introduction of a local coordinate system rotating with nodal cross-sections results in an autonomous transfer matrix for this Floquet system. Eigenanalysis reveals four real unity eigenvalues, indicating tension–torsion coupling, and equivalent continuum properties such as Poisson’s ratio, cross-sectional area, torsion constant and the tension–torsion coupling coefficient are determined. A variety of real and complex near diagonal Jordan decompositions are possible for the multiple (eight) complex unity eigenvalues and these are discussed in some detail. Analysis of the associated principal vectors shows that a bending moment produces curvature in the plane of the moment, together with shear deformation in the perpendicular plane, but no bending–bending coupling; the choice of a structure having an equilateral triangular cross-section is thought responsible for this unexpected behaviour, as the equivalent continuum second moments of area are equal about all cross-sectional axes. In addition, an asymmetric stiffness matrix is obtained for bending moment and shearing force coupling, and possible causes are discussed.     

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.     

波形钢腹板PC组合箱梁纯扭作用下抗扭承载力计算模型

以扭转理论和固定角软化桁架模型为基础,推导了纯扭作用下波形钢腹板PC组合箱梁的抗扭承载力计算模型.该模型针对组合箱梁极限状态时,波形钢腹板处于屈服和未屈服状态两种情况,分别建立了剪应交连续与剪力流连续方程,模型同时考虑了不同位置预应力筋对抗扭承载力的不同作用.通过四根试验梁抗扭承载力理论计算结果与试验数据的对比分析,结果表明本文提出的计算模型不仅能够有效的计算抗扭承载力,还能计算出混凝土开裂后的扭矩-扭率曲线.     

薄壁管预扭冲击拉伸实验装置的研制

材料在复合应力下的拉压或扭转行为常常与单轴拉压或者单轴扭转应力下不同,拉压以及扭转常体现出相互影响的特点。而复合应力下的塑性波可以有效地反映材料在复合应力下的本构行为特点。本文通过对霍普金森压杆进行改造,建立了一套薄壁管预扭冲击拉伸的实验装置,可以在薄壁管内产生拉扭耦合塑性波,并对率相关材料304不锈钢进行了薄壁管预扭拉伸实验研究,得到了该材料的拉扭耦合塑性波。结果表明,304不锈钢薄壁管的拉扭耦合塑性波具有明显的耦合快波和耦合慢波的双波结构,并且快波慢波之间没有恒值区间隔。同时也表明,该装置可以实现预期效果,并且可以有效避免薄壁管屈曲的产生。     

A novel class of graded cylinders neutral to host shafts of arbitrary cross-sections under torsion

A novel class of graded cylinders is proposed as neutral inclusions inside host shafts of arbitrary cross-section under Saint-Venant's torsion. The graded cylinder is made of cylindrically orthotropic materials with position varying quantities. The profiles of the two distinct shear moduli in the radial and tangential directions follow specific forms based on an arbitrarily selected function along the radial distance. We show that this type of graded cylinders can serve as universal neutral inclusions within host shafts of arbitrary cross-sections. In addition, we find that the associated warping fields can be exactly determined in terms of simple exponents of the selected function. This suggests that, by tuning the gradation parameter, one can manipulate the warping field of the inserted cylinder without disturbing the fields inside the host shaft. This finding is an original contribution to the existing solvable configurations of composite shafts under torsion.     

扇形截面杆扭转问题的差分线法

导出了扇形截面杆扭转问题偏微分方程的差分线法常微分方程组, 并解析求解了该方程组, 得到了扭转应力函数的半解析解, 计算了扭转应力及扭转刚度. 计算过程中, 用追赶法计算 常微分方程组的特解, 用公式计算三对角矩阵的特征值与特征向量, 利用实对阵矩阵的特征 向量相互正交的特性避免矩阵求逆计算, 利用复化梯形公式计算扭转刚度. 整个求解过程在 角度方向离散微分方程和用复化梯形公式进行面积积分时引入了误差, 其他求解过程是精确 的. 计算结果与已有结果进行了对比, 显示了算法的正确性. 该算法对工程中扇形截面扭 转杆的设计有一定的实用价值.     

A coupling cross-section finite element model for torsion analysis of prismatic bars

A general finite element model has been developed for the analysis of prismatic bars subject to torsional loading by modelling only a small slice of the bar. Exact analytical coupling deformation relationships between the artificial cross-sections, which are independent of the position of axis of rotation, have been formulated. Three examples from the range of analyses that have been evaluated have been selected to demonstrate the accuracy and effectiveness of the method. Analyses for an orthotropic elastic square cross-section bar, an elastic–plastic circular cross-section shaft containing a radial crack, and geometrically nonlinear deformation of a thin-walled I-section beam are presented and compared with previous results, where available.     

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

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

薄壁杆约束扭转的单肢解析化分析方法

针对薄壁杆件约柬扭转的基本受力反应分析问题,采用与符拉索夫的经典约束扭转理论截然不同的立论途径,将自由扭转刚度视为调整因素,而将剥离了该抗力效应的薄壁杆件基本体系作为主要分析对象,推导出基于翘曲理论并考虑了自由扭转刚度影响的杆件刚度方程及结问荷载的等效措施.此外,提出了基于微段薄壁杆简化单刚的有效数值化分析方法.上述研究表明,对于钢构件等具有薄壁构型截面形式的杆件,这种单肢解析化分析方法无需进行复杂的截面特性如弯、翘曲惯矩或弯、形心坐标等的分析,力与变形的表达简洁、直接,本文的抗力性能分析和表述方式可为薄壁构件的稳定及畸变等研究提供新的思路切入点.     

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

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

扇形截面杆扭转的应力分析

柱体扭转的基本方程为非齐次偏微分方程,在极坐标系下,利用分离变量法及傅立叶级数展开法,求出了扭转应力函数,进一步即可计算出扇形截面杆在外力偶作用下,扭转角和横截面上剪应力的精确解答。这种方法为精确解法,在各种机械及其他工程设备中,对受扭转作用的扇形截面杆设计,有一定实用价值。     

Prandtl’s formulation for the Saint–Venant’s torsion of homogeneous piezoelectric beams

The Saint–Venant torsional problem for homogeneous, monoclinic piezoelectric beams is formulated in terms of Prandtl’s stress function and electric displacement potential function. The analytical approach presented in this paper generalizes the known formulation of Prandtl’s solution which refers to homogeneous elastic beams. The Prandtl’s stress function and electric displacement potential function satisfy the so called coupled Dirichlet problem (CDP) in the cross-sectional domain. A direct and a variational formulation are developed. Exact analytical solutions for solid elliptical cross-section and hollow circular cross-section and an approximate solution based on a variational formulation for thin-walled closed cross-section are presented.