基于拟应变法的壳单元大变形分析及在冲击动力问题中的应用

为了简便有效地解决板壳结构的大变形问题，本文针对八节点相对自由度壳单元进行研究。该单元的位移场由壳的中面节点位移和上表面节点的相对位移组成，不带有转动变量。所有的研究都是基于完全的三维位移、应力、应变场。采用拟应变法，对应变场另行假设，能够改善该单元在大变形情况下的计算精度。通过引入Wilson非协调模式，构造了大变形情况下的拟应变场表达式，给出了该单元用于解决非线性动力分析问题的有限元求解方程。通过算例表明，本文针对相对自由度壳单元提出的方法及推导的公式，能够解决冲击动力问题中的大变形问题。     

适用于壳体h型自适应有限元分析的一组新单元

推导出一组适用于h型自适应分析的四边形蜕化壳元。对于大多数壳体结构,壳单元的刚度矩阵可分为薄膜、弯曲和剪切三部分。对薄膜部分本文采用杂交应力元方法进行设计,独立假设薄膜应力场以改善其精度;弯曲部分的刚度矩阵则依然由基于位移的应变来获得;而剪切部分则采用假设自然应变的方法来获得能克服薄壳下剪切自锁的新剪应变并用于计算此部...     

一种带有沙漏控制的新型EAS实体壳单元

基于 Hu-Washizu 变分原理推导了一种带有沙漏控制的实体壳单元的显式有限元列式;采用面内单点积分方案,提高了计算效率;同时引入沙漏控制力,抑制了沙漏现象.基于缩减积分方案,采用 B-bar 法消除体积自锁,并通过添加七个改善拟应变参数解决了泊松自锁和剪切自锁.采用改善拟应变法消除剪切自锁,使得表达式简洁.利用这种显式实体壳单元模型对3个非线性变形的标准算例进行了计算,并与相关参考文献和有限元软件 ABAQUS 的计算结果进行了比较.结果表明:该实体壳单元具有较高的计算精度,可有效地解决板壳非线性大变形分析问题,具有很好的工程应用前景.     

A URI 4-Node quadrilateral element by assumed strain

In this paper one-point quadrature ““““assumed strain““““ mixed element formulation based on the Hu-Washizu variational principle is presented. Special care is taken to avoid hourglass modes and volumetric locking as well as shear locking. The assumed strain fields are constructed so that those portions of the fields which lead to volumetric and shear locking phenomena are eliminated by projection, while the implementation of the proposed URI scheme is straightforward to suppress hourglass modes. In order to treat geometric nonlinearities simply and efficiently, a corotational coordinate system is used. Several numerical examples are given to demonstrate the performance of the suggested formulation, including nonlinear static/dynamic mechanical problems.     

On the Reissner theory of bending of flastic plates

The Reissner equations of elastic plates are rederived on the bases of the incomplete generalized variational principle of Complementary energy. The Stress function is naturally obtained from the variational Calculation in the form of Lagrange multiplier. The stucture of solutions of the Reissner equations is thus defined. On the bases of these discussions, a simplified theory has been put forward, in which the equations of equilibrium involving the shearing influence can be reduced into a fourth order differential equation similar to those of the Classical theory of plates.     

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

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

一种考虑横向剪切变形的三角形板弯曲单元

本文从部分协调的三角形薄板弯曲单元出发,并假设横向剪切应变在单元内线性变化,提出了一种考虑横向剪切变形的具有15个自由度的三角平板弯曲单元。该单元应用于薄板和中等厚度板分析均有较高的精度,计算效率高,可用于工程中具有复杂形状的薄板和中等厚度板结构的分析。     

一个基于膜板比拟理论的四节点二十四自由度的平板壳单元

应用膜板比拟关系 ,可以避开 c1 连续性的困难 ,为板单元的构造提供了一种新的途径 ,并已成功地构造出一系列相应的板单元。本文构造了一个四节点二十四自由度的平板壳单元 ,该单元由平面四节点理性元 RQ4(膜部分 )和由膜板比拟理论构造的一个四节点十二自由度的板单元 (弯曲部分 )构成。该单元构造简单 ,数值结果表明具有很好的收敛性和精度。     

Locking-Free Degenerated Isoparametric Shell Element

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｅｇｅｎｅｒａｔｅｄｉｓｏｐａｒａｍｅｔｒｉｃｓｈｅｌｌｅｌｅｍｅｎｔｈａｓｂｅｃｏｍｅｉｎｃｒｅａｓｉｎｇｌｙｐｏｐｕｌａｒｉｎｐｒａｃｔｉｃａｌｅｎｇｉｎｅｅｒｉｎｇｉｎａｗｉｄｅｒａｎｇｅ.ＴｈｉｓｇｅｎｅｒａｌｓｈｅｌｌｅｌｅｍｅｎｔｗａｓｏｒｉｇｉｎａｌｌｙｉｎｔｒｏｄｕｃｅｄｂｙＡｈｍａｄｅｔａｌ[1]ｆｏｒｔｈｅｌｉｎｅａｒａｎａｌｙｓｉｓｏｆｍｏｄｅｒａｔｅｌｙｔｈｉｃｋｓｈｅｌｌｓ .Ｉｔｗａｓｄｅｇｅｎｅｒａｔｅｄｆｒｏｍａｔｈｒｅｅ＿ｄｉｍｅｎｓｉ…     

New mixed plate-bending elements with shear strain interpolations

This paper is concerned with two mixed plate-bending elements with shear strain interpolations, a quadrilateral element and an 8-node serendipity-type element based on discussions on the element proposed in Ref.[1]. The shear strains and inner-forces in the natural coordinates are interpolated in an element and then transformed into Cartesian coordinates in accordance with covariant and contravariant tensor transformations, respectively. The quadrilateral element coincides with the element in Ref.[1] when it is rectangular. Numerical examples show that the two new elements are free from shear locking and spurious kinematic modes under regular and irregular meshes and have the advantages of being insensitive to element distortion and able to give fairly accurate results.The Project supported by National Natural Science Foundation of China.     

一个不闭锁和抗畸变的四边形厚板元

构造一个彻底消除剪切闭锁现象并且对网格畸变不敏感的四边形厚薄板通用单元RPAQ。在方法上有三个特点：第一,在厚板挠度和转角的试函数中,采用了合理匹配方案,从而在源头上彻底消除了剪切闭锁现象;第二,采用四边形面积坐标,以代替通常的等参坐标,从而使网格畸变时仍然保持高精度;第三,采用广义协调元做法,使协调条件的采用灵活多样,并保证单元的收敛性。进行了一系列数值例题测试,表明单元RPAQ能自动消除闭锁现象,在由薄板到厚板的不同情况下,在各种网格畸变的情况下,都能体现出良好的精度和数值稳定性。     

有限厚双边U型切口板的三维弹塑性有限元分析

针对单向拉伸载荷作用下的有限厚双边 U型切口板 ,本文对其切口根部的三维效应进行了详细的三维弹塑性大变形有限元分析 ,得到了不同于平面切口问题的新结果 :( 1 )在线弹性阶段 ,三维切口根部的应力集中因子 Kt高于二维情形 ,且 Kt在板厚 ( B)约为 8倍切口根部半径 (ρ)时达到峰值 ;进入塑性后 ,切口根部的 Kt随载荷水平的增加逐渐下降 ,并介于平面应力和平面应变情形之间 ;( 2 )切口前缘等效应变εe的最大值出现在约束最高的中面 ,且可比约束最低的自由表面处的应变值高一倍 ;( 3)薄板中 ,塑性区尺寸 xp 的最大值不是出现在约束较低的自由表面 ,而是出现在约束较高的中面 ,但随 B/ρ的增大逐渐向自由表面外推 ;( 4 )面内应力比 Tx在切口根部塑性区内的分布对板厚和载荷水平不敏感 ,可以用 Hill的滑移线解很好地近似。另外 ,对三维影响区内一些重要的三维约束参量进行了详细分析 ,总结了它们的三维分布规律。     

Geometrically non-linear analysis of laminated composite structures using a 4-node co-rotational shell element with enhanced strains

To demonstrate the solutions of linear and geometrically non-linear analysis of laminated composite plates and shells, the co-rotational non-linear formulation of the shell element is presented. The combinations of an enhanced assumed strain (EAS) in the membrane strains and assumed natural strains (ANS) in the shear strains improve the behavior of 4-node shell element. To secure computational efficiency in the incremental non-linear analysis, the present element uses the form of the resultant forces pre-integrated through the thickness. The transverse shear stiffness of the laminates is defined by an equilibrium approach instead of the shear correction factor. Numerical examples of this study show very good agreement with the references.     

弹塑性板壳结构非线性有限元分析

本文根据塑性流动理论的基本公式，由隐式积分导出了与路径无关的变量更新算法和一致切线模量。采用单元广义应力应变直接离散塑性流动定律，构造了杂交应力单元一致切线刚度矩阵的显式表达式，编制了结构有限元程序ＳＡＦＥ，数值算例表明：本文的计算方法和计算程序是正确可靠的，可用于弹塑性板壳结构的非线性分析，计算结果屈曲临界载荷和极限承载能力。     

Limit analysis of viscoplastic thick-walled cylinder and spherical shell under internal pressure using a strain gradient plasticity theory

Plastic limit load of viscoplastic thick-walled cylinder and spherical shell subjected to internal pressure is investigated analytically using a strain gradient plasticity theory. As a result, the current solutions can capture the size effect at the micron scale. Numerical results show that the smaller the inner radius of the cylinder or spherical shell, the more significant the scale effects. Results also show that the size effect is more evident with increasing strain or strain-rate sensitivity index. The classical plastic-based solutions of the same problems are shown to be a special case of the present solution.     

板壳问题的三维无网格伽辽金直接分析法

对于板壳问题,共有三种数值模拟方案:线性或非线性的板壳理论、退化连续体方案和直接三维连续体方案。无网格法近似函数可具有C1甚至更高的连续性,便于在K irchhoff-Love理论中应用。但当各种无网格法用于M ind lin-R e issner板理论时,会遇到数值锁死的困扰。对比之下,三维连续体方案是最简单,最精确但并不常用的一种方案。无网格法近似函数具有高度光滑性,在板壳的厚度方向仅布置2~5层点就可以很好地捕捉此方向场的梯度,同时还可以在一定参数范围内避免剪切和体积锁死,在处理复杂本构关系、非线性板壳等问题中更是具有很大优势。本文采用无网格伽辽金法(EFG)和三维连续体方案分析了线性板壳问题,与有限单元法做了对比,并讨论了数值锁死等问题。     

On the stability analysis of plates and shells using a quadrilateral, 16-degrees of freedom plat shell element DKQ16

The linear buckling problems of plates and shells were analysed using a recently developped quadrilateral, 16-degrees of freedom flat shell element ( called DKQ16 ). The geometrical stiffness matrix was established. Comparison of the numerical results for several typical problems shows that the DKQ16 element has a very good precision for the linear buckling problems of plates and shells.     

Novel shear deformation model for moderately thick and deep laminated composite conoidal shell

This article presents a novel mathematical model for moderately thick and deep laminated composite conoidal shell. The zero transverse shear stress at top and bottom of conoidal shell conditions is applied. Novelty in the present formulation is the inclusion of curvature effect in displacement field and cross curvature effect in strain field. This present model is suitable for deep and moderately thick conoidal shell. The peculiarity in the conoidal shell is that due to its complex geometry, its peak value of transverse deflection is not at its center like other shells. The C¹ continuity requirement associated with the present model has been suitably circumvented. A nine-node curved quadratic isoparametric element with seven nodal unknowns per node is used in finite element formulation of the proposed mathematical model. The present model results are compared with experimental, elasticity, and numerical results available in the literature. This is the first effort to solve the problem of moderately thick and deep laminated composite conoidal shell using parabolic transverse shear strain deformation across the thickness of conoidal shell. Many new numerical problems are solved for the static study of moderately thick and deep laminated composite conoidal shell considering 10 different practical boundary conditions, four types of loadings, six different hl/hh (minimum rise/maximum rise) ratios, and four different laminations.     

应用平板壳单元DKQl6分析板壳结构的几何非线性

应用新近开发的四边形十六自由度离Kirchhoff平板壳单元DKQl6，分析了板壳结构的几何非线性问题，采用Total Lagrange格式，在小应交、中等转动的假定下，建立了该单元几何刚度阵和大位移矩阵．非线性方程采用位移引导或弧长引导的牛顿-拉夫森增量迭代法求解．讨论了网格和加载步效对收敛性的影响，通过对典型算例的计算以及与其它单元的比较，说明了DKQl6单元在板壳结构几何非线性分析中也有良好的精度．     

自由边界条件轴对称薄圆柱壳体的混合有限元数值解

基于文「８」发展的混合结构问题稳定方法的框架,针对自由边界条件轴对称薄圆柱壳体问题,给出了其新型混合有限元计算模型。应用线性元,已证明该模型是强制的,其对称正定代数方程组的条件数Ｏ（ｈ＾－２）。