Triangular elements for Reissner-Mindlin plate

I-IntroductionInrecentyearsmuchresearchefforthasbeenspentonthedevel0pmentofreliableandefficientplateelementsbasedonReissner-Mindlintheory.Adifficultyisthelockingbehaviorexhibitedasl- 0(tisthethicknessoftheplate)forlowerorderc'elements.BasedonTaylorexpendi…     

三维退化的壳理论和假定应变的壳单元

本文将Reissner-Mindlin板理论推广到空间曲壳结构,可称为Reissner-Mindlin型壳理论。从这种理论出发,可直接导出C(0)连续的壳体单元,即考虑横向剪切变型的影向的壳体单元,这种单元在国外已被广泛地采用,为克服这种单元在应用中所出现的剪切和膜的锁制现象同时又防止出现任何零能模式,作者提出了一种采用假定应变的新的壳单元公式,并对这种单元进行了广泛的数值试验,结果表明这种单元具有较高的精度和良好的性能。     

Three dimensional degenerated shell theory and assumed strain shell elements

In this paper Reissner-Mindlin plate theory is extended to cater for curved shell structures. It can be considered as Reissner-Mindlin type shell theory. From this theory, the C(O) continuity formulation of shell elements of taking account the transverse shear deformation could be derived directly. These degenerated shell elements have been widely employed. To overcome the locking of shear and membrane and avoid zero energy modes the author proposed the formulation of the new elements with assumed strains. A wide range of numerical tests was conducted and the results illustrate that the assumed strain elements possess high accuracy and good performance.     

????Timoshenko????????????

将理性有限元法引入到Timoshenko梁问题中,提出了一种理性Timoshenko梁单元,克服了 剪切锁死现象. 在推导控制方程时,与传统有限元方法采用Lagrange插值不同, 理性有限元法用Timoshenko梁弯曲问题的基本解逼近单元内部场. 运用该梁单元分析 Timoshenko梁时,无需缩减积分,就能避免剪切锁死,并且极大地提高了计算精度,说明 理性有限元法具有广泛的应用前景.     

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

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

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

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

一种精确积分的四结点板弯曲单元

本文从Ｍｉｎｄｌｉｎ／Ｒｅｉｓｓｎｅｒ理论出发,采用一种新的平行四边形母单元和相应的形函数推导四结点板弯曲单元刚度矩阵的精确积分解。弯曲应变和横向剪切应变分别采用不同的插值公式构成单元刚度矩阵。理论和算例分析表明本文方法克服了“闭锁”现象并能应用于很薄的板,单元刚度矩阵计算速度比采用数值积分计算的同类单元的快四倍。     

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.     

加权残值法在钢筋混凝土拱桥非线性有限元分析中的应用

本文用圆弧梁离散拱肋;用圆柱拖带坐标、三次位移插值函数及平截面假定来描述单元位形;用加权残值配点法来消除曲梁单元的剪力与膜力闭锁。按基于连续介质力学的Ｕ．Ｌ．列式建立单元增量平衡方程,以考虑几何非线性。假定钢筋为理想弹塑性材料。按三参数各向同性强化塑性模型,建立混凝土的弹塑性本构矩阵。将拱单元分段分块,根据钢筋及砼的本构特性,建立拱单元及梁段单元的弹塑性刚度矩阵,以考虑材料非线性。用编制的程序对两座模型拱桥进行计算,计算结果与模型测试结果接近     

基于mbS模式的有限单元法

mbS模式及其有限元法是在固体和结构分析模型中引入薄膜、弯曲和剪切理论,且采用纯拉压、纯弯和纯剪单元进行分析的数值方法。在时空系中剖分物质单元和时间单元上构造以指数函数和贝塞尔函数为插入函数且按Lagrange插值条件的薄膜、弯曲和剪切等基本位移函数,由此得到更加完备和耦合的固体和结构实体单元的变形模式,根据能量泛函变分原理得到静动力有限元基本方程的一致格式。研究表明,mbS模式及其有限元法可用于梁柱和板壳等结构的静动力分析及屈曲分析。     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

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

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

Conforming radial point interpolation method for spatial shell structures on the stress-resultant shell theory

The implementation of the conforming radial point interpolation method (CRPIM) for spatial thick shell structures is presented in this paper. The formulation of the discrete system equations is derived from a stress-resultant geometrically exact theory of shear flexible shells based on the Cosserat surface. A discrete singularity-free mapping between the five degrees of freedom of the Cosserat surface and the normal formulation with six degrees of freedom is constructed by exploiting the geometry connection between the orthogonal group and the unit sphere. A radial basis function is used in both the construction of shape functions based on arbitrarily distributed nodes as well as in the surface approximation of general spatial shell geometries. The major advantage of the CRPIM is that the shape functions possess a delta function property and the interpolation function obtained passes through all the scattered points in the influence domain. Thus, essential boundary conditions can be easily imposed, as in finite element method. A range of shape parameters is studied to examine the performance of CRPIM for shells, and optimal values are proposed. The phenomena of shear locking and membrane locking are illustrated by presenting the membrane and shear energies as fractions of the total energy. Several benchmark problems for shells are analyzed to demonstrate the validity and efficiency of the present CRPIM. The convergence rate of the results using a Gaussian (EXP) radial basis is relatively high compared to those using a multi-quadric (MQ) radial basis for the shell problems.     

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.     

New method for geometric nonlinear analysis of large displacement drill strings

Introduction Drillstringsina3Dlargedisplacementwellboreareusuallyregardedas3Dcurved beamswithlargedisplacementandlargerotationbutwithsmallstrain,thus,therelationship betweenstressandstrainislinear.Thedeformationofdrillstringscanbecomposedoftwo parts:oneiscalledtheinitialdeformationcausedwhentheaxisofthedrillstringisdeformed fromitsstraightconfigurationtothepositioncoincidentwiththewellaxis,involvinglarge rotationsandwiththeimplicationoffiniterotationalproblem.Theotheristhedeformation fromthew…     

杂交应力通用壳元的模式优化

文章把杂交元的优化设计原理用于厚薄壳问题,构造了正交曲线坐标系中的Mindlin型四节点二十自由度壳元。该元通用于各种厚度、任意形状的深壳和扁壳,构造简单,精度较高且没有自锁现象。     

一个厚薄板通用的三角形板弯元

本文提出了一个新的厚薄板通用的三角形板弯元,根据Hellinger-Reissner变分原理推导出来的新单元具有独立的转角、位移、剪应变和弯曲应变的插入函数,它没有其它混合元所存在的一些缺点。一些典型例子表明这个新单元具有很好的特性,刚度矩阵简单,计算精度高,收敛速度快以及厚薄板通用性强。     

用于超弹性分析的高效杂交应变--EAS固体壳单元的研究

提出了一个用于橡胶材料分析的弱变分方程，基于一个可压缩Neo—Hookean模型，将EAS模式和杂交元法有机地结合起来，推导了一个高效、稳定的杂交应变——EAS固体壳单元，通过巧妙地选择应交插值函数和本文中提出的一个精化措施，克服了橡胶材料所表现出的大应变超弹性本构为杂交元正交化法的实施所带来的困难，保证了整个单元列式都仅采用低阶高斯积分，显著提高了计算效率，确保橡胶材料不可压缩性计算的顺利进行，并克服了固体壳元的厚度自锁问题。