The method of solving axisymmetric problems in elastic space by complex function and some examples

This paper proves Love’s stress function of space axisymmetric problem can be represented by choosing two generalized analytic functions of complex variates reasonably^[1], and deduces the expressions of the components of stress displacements and boundary conditions in complex function. To present the feasibility of the method here and examining the truth of the formulae founded in this paper, the problem of circular shaft with globular cavity pressed on the side and pulled at the ends is solved by using power series and the result is the same as that solved by other methods. In the end, the problem of a cone sheared by uniform shear stress on the sideface is solved, and the solution of a cone acted on by gravity is given by converting constant body forces into surface forces.     

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

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

A NEW HIGH-ORDER MULTI-JOINT FINITE ELEMENTFOR THIN-WALLED BAR

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｉｎ＿ｗａｌｌｅｄｓｔｒｕｃｔｕｒｅｉｓｗｉｄｅｌｙｕｓｅｄｉｎｅｎｇｉｎｅｅｒｉｎｇ .Ｉｔｎｏｔｏｎｌｙｈａｓａｒｏｌｅｏｆｄｅｃｒｅａｓｉｎｇｑｕａｎｔｉｔｙｏｆｓｔｒｕｃｔｕｒｅ ,ｕｓｉｎｇｍａｔｅｒｉａｌｒｅａｓｏｎａｂｌｙ ,ａｎｄｉｍｐｒｏｖｉｎｇｅｃｏｎｏｍｉｃｂｅｎｅｆｉｔｓｇｒｅａｔｌｙｂｕｔａｌｓｏｈａｓａｒｏｌｅｏｆｂｅａｕｔｉｆｙｉｎｇｓｔｒｕｃｔｕｒｅ’ｓｆａｃａｄｅ.Ｎｏｗａｄａｙｓ,ｂｅｃａｕｓｅｏｆｔｈｅｓｅａｐｐａｒｅｎｔａｄｖａｎｔａｇ…     

Free torsion of thin-walled structural members of open- and closed-sections

Free torsion of thin-walled structures of open- and closed-sections is a classical elastic mechanics problem, which, in literature, is often solved by the method of membrane analogy. The method of membrane analogy, however, can be only applied to structures of a single material. If the structure consists of both open- and closed-sections, the method of membrane analogy is difficult to be applied. In this paper, a new method is presented for solving the free torsion of thin-walled structures of open- and/or closed- sections with multiple materials. By utilizing a simple statically indeterminate concept, torsional equations are derived based on the equilibrium and compatibility conditions. The method presented here not only is very simple and easy to understand but also can be applied to thin-walled structures of combined open- and closed-sections with multiple materials.     

开口薄壁杆件结构稳定分析的精确单元和两步求解算法

从控制微分方程的通解出发,构造受偏心压力作用开口薄壁杆件的精确形函数,建立用于开口薄壁杆件结构稳定性分析的精确有限元,得到了单元刚度矩阵和几何刚度矩阵的显式表达,提出了计算给定区间内各阶临界荷载以及相应失稳模态的两步计算方法。计算结果表明,与常规单元相比,采用精确单元无需进行网格细分就可以获得精确的数值结果,结合本文的两步求解算法,可以准确获得给定区间内全部临界荷载和失稳模态。     

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.     

辐射换热下瞬态热-结构分析的一种空间薄壁杆单元

发展了一种用于辐射换热条件下瞬态热-结构分析的空间薄壁杆单元,其截面形式可以是任意形状的闭口截面和单支开口截面。该单元温度场分解为平均温度和多谐摄动温度,沿杆轴方向采用两结点线性插值,沿杆截面周向用三角函数展开,每结点含多个解耦的自由度,其中结点平均温度方程同传统一维温度有限元方程为非线性,各谐摄动温度方程为线性,然后利用Wilson-θ法求解结构的瞬态温度场。本文选择了两节点Ber-noulli直梁单元得到准静态热弹性有限元方程并求解,针对非对称开口截面考虑了翘曲变形及弯扭耦合的影响。温度场引起的等效热载荷不仅包括常规的热轴力,还包括热弯矩以及热双力矩。本文针对不同截面形式的梁单元给出了瞬态温度场以及热变形的验证算例,并通过与商业程序中二维壳元计算结果的比较说明了本文所提出方法的正确性和高效性。     

ANALYSIS OF COMPOSITE LAMINATE BEAMS USING COUPLING CROSS-SECTION FINITE ELEMENT METHOD

Beams and plates manufactured from laminates of composite materials have distinct advantages in a significant number of applications. However, the anisotropy arising from these materials adds a significant degree of complexity, and thus time, to the stress and deformation analyses of such components, even using numerical approaches such as finite elements. The analysis of composite laminate beams subjected to uniform extension, bending, and/or twisting loads was performed by a novel implementation of the usual finite element method. Due to the symmetric features of the deformations, only a thin slice of the beam to be analysed needs to be modelled. Conventional three-dimensional solid finite elements were used for the structural discretization. The accurate deformation relationships were formulated and implemented through the coupling of nodal translational degrees of freedom in the numerical analysis. A sample solution for a rectangular composite laminate beam is presented to show the validity and accuracy of the proposed method.     

Computational analysis of multi-stepped beams and beams with linearly-varying heights implementing closed-form finite element formulation for multi-cracked beam elements

The model where the cracks are represented by means of internal hinges endowed with rotational springs has been shown to enable simple and effective representation of transversely-cracked slender Euler–Bernoulli beams subjected to small deflections. It, namely, provides reliable results when compared to detailed 2D and 3D models even if the basic linear moment–rotation constitutive law is adopted.This paper extends the utilisation of this model as it presents the derivation of a closed-form stiffness matrix and a load vector for slender multi-stepped beams and beams with linearly-varying heights. The principle of virtual work allows for the simple inclusion of an arbitrary number of transverse cracks. The derived at matrix and vector define an ‘exact’ finite element for the utilised simplified computational model. The presented element can be implemented for analysing multi-cracked beams by using just one finite element per structural beam member. The presented expressions for a stepped-beam are not exclusively limited to this kind of height variation, as by proper discretisation an arbitrary variation of a cross-section’s height can be adequately modelled.The accurate displacement functions presented for both types of considered beams complete the derivations. All the presented expressions can be easily utilised for achieving computationally-efficient and truthful analyses.     

Response of flexure–torsion coupled composite thin-walled beams with closed cross-sections to random loads

A study of the flexure–torsion coupled random response of the composite beams with solid or thin-walled closed-sections subjected to various types of concentrated and distributed random excitations is dealt with in this paper. The effects of flexure–torsion coupling, shear deformation and rotary inertia are included in the present formulations. The random excitations are assumed to be stationary, ergodic and Gaussian. Analytical expressions for the displacement response of the composite beams are obtained by using normal mode superposition method combined with frequency response function method. The present method can produce the effective solutions for the composite Timoshenko beams with circumferentially antisymmetric (CAS) configuration and more general beam assemblages of connected beams. The influences of flexure–torsion coupling, shear deformation and rotary inertia on the random response of an appropriately chosen composite beam from the literature are demonstrated and discussed.     

基于Hamilton体系的弹性力学问题的比例边界有限元方法

比例边界有限元方法是求解偏微分方程的一种半解析半数值解法。对于弹性力学问题,可采用基于力学相似性、基于比例坐标相似变换的加权余量法和虚功原理得到以位移为未知量的系统控制方程,属于Lagrange体系。但在求解时,又引入了表面力为未知量,控制方程属于Hamilton体系。因而,本文提出在比例边界有限元离散方法的基础上,利...     

板壳畸变单元的无网格和有限元自动耦合方法

板壳大变形时单元的严重畸变会使计算精度降低。无网格局部Petrov-Galerkin法是一种真正的无网格方法,能够消除网格畸变,但比有限元法计算效率低。根据板壳网格畸变的局部性特点,利用过渡单元法,基于板壳网格质量,建立了板壳的网格严重畸变区域由有限元分析切换为无网格分析的自动耦合算法,实现了有限元法和无网格局部彼得罗夫．迦辽金法的耦合。应用实例表明：通过自适应耦合,既能发挥有限元法计算效率高的特点,又能发挥无网格法适合大变形分析、没有网格畸变造成计算困难的特点。     

A new higher-order shear deformation theory and refined beam element of composite laminates

A new higher-order shear deformation theory based on global-local superposition technique is developed. The theory satisfies the free surface conditions and the geometric and stress continuity conditions at interfaces. The global displacement components are of the Reddy theory and local components are of the internal first to third-order terms in each layer. A two-node beam element based on this theory is proposed. The solutions are compared with 3D-elasticity solutions. Numerical results show that present beam element has higher computational efficiency and higher accuracy.The project supported by the National Natural Science Foundation of China (10172023)     

压电层合梁静动力学分析的高精度梁单元

给出了一个对复合材料压电层合梁进行数值分析的高精度压电层合梁单元。基于Shi三阶剪切变形板理论的位移场和Layer-wise理论的电势场,用力-电耦合的变分原理及Hamilton原理推导了压电层合梁单元列式。采用拟协调元方法推导了一个可显式给出单元刚度矩阵的两节点压电层合梁单元,并应用于压电层合梁的力-电耦合弯曲和自由振动分析。计算结果表明,该梁单元给出的梁挠度和固有频率与解析解吻合良好,并优于其它梁单元的计算结果,说明了本文所给压电层合梁单元的可靠性和准确性。研究结果可为力-电耦合作用下压电层合梁的力学分析提供一个简单、精确且高效的压电层合梁单元。     

基于微分求积有限元法的双层微板系统振动特性研究

基于修正的偶应力理论和两变量精化的剪切变形理论,建立了由Winkler-Pasternak连续弹性夹层连接的双层微板系统的自由振动模型,着重推导了系统异步振动的运动微分方程和势能泛函。融合Gauss-Lobatto求积准则和微分求积准则构造了具有C¹连续性的微分求积有限元。通过与已有文献进行对比,验证了数值方法的有效性。详细讨论了各种因素对系统同步和异步振动特性的影响。结果表明,系统的自由振动特性对材料尺度参数、长宽比、长厚比以及边界条件呈现出依赖性;弹性夹层刚度仅对系统异步振动产生作用;随着模态阶次的增大,材料尺度参数和弹性夹层刚度对异步振动频率和模态的影响变得显著。     

爆炸冲击波对人体胸部创伤机理的有限元方法研究

通过建立简化的人体胸部三维有限元模型,模拟爆炸冲击波与人体胸部作用,根据人体胸部各个 器官的不同,选择合理的材料模型和参数,提出LS-DYNA有限元程序局部平面波改进方法,研究爆炸冲击波 与人体胸部作用的力学过程。依据人体胸部各个器官运动的速度差,预测创伤的区域分布;给出肺模型的压、 拉应力及剪切应力的变化规律,分析肺的创伤区域的分布,与解剖实验结果基本一致。对比Bowen创伤曲 线,证明人体胸部三维有限元模型可以有效得到冲击波创伤特征。     

新型整体-局部高阶剪切变形理论及层合板精化三角形单元

基于整体-局部位移方法,建立了一种高阶剪切变形理论。整体位移部分采用的是Reddy理论的位移模式（1984）,局部位移为LIXY等（1997）建立的1,2-3理论的局部函数。这一理论使满足自由表面条件的Red@理论进一步满足层间位移、应力连续,同时有效减少了1,2—3理论的未知数个数。基于此理论深入开展了有限元法研究,建立了满足C^1连续条件的精化三节点三角形单元（每个节点参数为9个）。计算结果表明：建立的精化单元能准确计算整体位移和层间应力。     

Differential quadrature method for bending of orthotropic plates with finite deformation and transverse shear effects

Based on the Reddy ‘s theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature ( DQ ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert ( DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.     

带旋转自由度C^0类任意四边形板（壳）单元

基于Ｒｅｉｓｓｎｅｒ－Ｍｉｎｄｉｌｉｎ板弯曲理论和Ｖｏｎ－Ｋａｒｍａｎ大挠度理论,采用单元域内和边界位移插值一致性的概念,将四节点等参弯曲单元与Ａｌｌｍａｎ膜变形二次插值模式相结合,对层合板壳的大挠度分析提供了一种实用的带旋转自由度的四节点Ｃ＾０类板单元。大量算例表明：该单元对板壳结构的线性强度、稳定性和后屈曲分析都表现出良好的收敛性和足够的工程精度。