平面4节点广义等参单元

在文献[1]建立的广义平面矩形单元的基础上,建立了平面4节点广义等参单元,并进行了深入研究,包括分片试验、单元的变形模态以及形函数的影响等.通过典型算例,对单元的性能进行了检验,数值结果表明,常规情况下建立的广义等参单元较传统单元具有更高的计算精度.     

广义节点有限元法

应用流形方法思想,通过引入广义节点的概念,对传统有限元方法进行改进,建立了可具有任意高阶多项式托值函数的广义节点有限元方法,计算结果表明,广义节点有限元方法较之传统有限元方法有较高的精度。     

关于广义等参有限单元的讨论

研究了新近提出的广义等参单元,讨论了这一新方法与现有方法的关系,提出了一个建立广义有限单元的一般性方法。     

平面4结点等参元高斯点应力的适用性

平面４结点等参元高斯点应力的适用性黄文彬李明瑞（北京农业工程大学，北京１０００８３）在平面弹性和弹塑性有限元分析中，广泛使用３～９结点等参单元，在使用４结点等参单元时，通常用２×２高斯点积分，在应力输出时，也常输出高斯点处应力，特别在弹塑性分析中，常...     

复杂结构二十节点等参元建模法

提出用单一的二十节点等参元模拟平面、梁、板壳和板-梁结构的建模方法.通过调整板-梁接触面耦合节点的数目,分析接触面耦合节点的数目对接触面应力分布和内力的影响.结果表明,二十节点等参元完全具有模拟平面问题、梁和板壳的能力,能提供更多的相关信息;通过调整接触面耦合节点的数目和组成方式,二十节点等参元独立建模法可以模拟板-梁结构中的刚性联接,铰链联接和介于二者之间的弹簧联接的优越性能.这种建模方法不仅大为简化了复杂结构的建模过程,而且使有限元模型建立在更直观和真实的物理基础之上,避免引入传统建模方式所必须引入的联接单元,有利于有限单元法的推广和应用.     

改进四节点等参单元在一步成形模拟中的应用

在基本上不增加计算量的前提下,进一步提高一步成形逆有限元法的模拟精度,将四节点等参膜单元应用于一步成形计算,对空间四边形单元的四个节点可能不在同一平面引起的单元翘曲进行了修正,编写了计算机程序,对L形件进行了实际模拟,并将模拟结果与增量法模拟结果进行了比较,二者吻合较好.进而验证了本文方法的有效性.     

张拉结构非线性分析的五节点等参单元

本文针对张拉结构的特点，提出了一种五节点等参数单元有限元模型，采用四次多项式作为位移插值函数及单元初始形状函数，并假定索是理想柔性的且满足虎克定律，基于修正的Ｌａｇｒａｎｇｉａｎ坐标描述法，建立了非线性有限元基本方程和切线刚度矩阵，利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法进行了实例计算。结果表明：本文方法精度极高，可供大跨度索网，索穹顶，拉线塔等张拉结构分析，设计时采用。     

孔洞对平面裂纹应力强度因子影响分析的广义参数Williams单元

本文采用圆形奇异区广义参数Williams单元（W单元）建立了中心裂纹与圆孔共存的平面应力模型,奇异区外围利用ABAQUS有限元软件自动网格离散技术与FORTRAN95编程前处理相结合,克服了自主编程中网格离散的局限性．算例分析了圆孔位置和几何参数对I-II混合型裂纹尖端应力强度因子(SIFs)的影响,并与扩展有限元法(XFEM)计算结果进行比较．结果表明：靠近圆孔一侧的裂尖SIFs大于远离圆孔一侧的裂尖SIFs;控制圆孔左边缘到裂纹中心的距离,则两侧裂尖SIFs随圆孔半径的增大而增大;圆孔中心与裂纹中心水平距离越远,圆孔对裂纹扩展的影响越小．同时,基于圆形奇异区的W单元直接计算得到的裂尖SIFs与扩展有限元法得到的解吻合较好,证明了W单元对奇异区离散形状不敏感,且具有高效率和高精度．     

配点型广义节点无网格法的基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。     

带旋转自由度的精化非协调平面四边形等参元

对于C~1类不协调元文[1]提出了一种精比直接刚度法。本文进一步将其应用到C~0类问题,建立一种带旋转自由度的不协调平面四边形单元,其协调部分是用人Allman插值法建立的单元函数,不协调部分用了四个内部自由度。该单元能保证通过分片试验,保持了单变量有限元列式简单、性能可靠(无多余零能模式及坐标不变性)等长处,同时,还具备多变量有限元(杂交/拟协调元)高精度的优点。算例表明,本文提出的单元收敛、可靠、高精度且高效率。     

等参元——修改等参元——最佳等参元

本文以二维问题为模型,建立了等参元、修改等参元和最佳等参元间的相互关系.给出了最佳等参元附加节点方程的解法.其数值解与理论分析一致.显示出修改等参元的许多良好特性.     

构造几何不敏感四边形膜元的广义协调方法

许多有限元模型在规则网格下具有很高的精度，但当网格畸变程度增大时，其计算精度也随之迅速下降．如何构造出对网格畸变不敏感的单元，长期以来是人们十分关注的课题．本文应用作者早期提出的广义协调方法，构造出具有平面内旋转自由度的任意四边形膜元．该单元不仅列式简洁，而且具有对网格畸变极不敏感的优异性能，为构造对网格畸变不敏感的优质单元提供了一个通用方法     

QUADRILATERAL 4-NODE QUASI-CONFORMING PLANE ELEMENT WITH INTERNAL PARAMETERS

在拟协调框架之下,利用新的内参形函数构造了一个四边形四节点拟协调平面单元. 新的内参位移函数也可以添加到等参单元Q4 中来构造新的内参型等参单元. 新构造的拟协调单元QC6N 具有显式刚度矩阵,因而效率更高. 数值例子表明相比于四节点等参单元,新构造的单元可以提高计算精度和抗网格畸变的能力.     

结构分析中的退化壳有限元方法

本文评述了退化壳有限元的现状和发展．首先简要介绍了平板壳有限元和以壳体理论为基础的曲壳有限元,接着介绍了单变量（位移型）退化壳有限元的研究现状,最后介绍了多变量（杂交、混合、拟协调和杂交／混合型）壳有限元．      

非线性广义杂交元

选择新的三类变量建立了放松单元间连续性条件的用于非线性有限元分析的泛函，并由此建立了基于不协调模式的非线性广义杂交元方法。     

非均质材料动力分析的广义多尺度有限元法

自然界和工程中的大部分材料都具有多尺度特征,当考察尺度小到一定程度后,都将表现出非均质性.针对非均质材料的动力问题,提出了一种广义多尺度有限元方法,其基本思想是利用静态凝聚法以及罚函数法构造能够反映单元内部材料非均质特性的多尺度位移基函数.与传统扩展多尺度有限元法中的基函数构造方式不同,广义多尺度有限元法的基函数无需通过在子网格域上多次求解椭圆问题得到,而可直接通过矩阵运算获得.其主要步骤如下:利用数值基函数将一个非均质单胞等效为一个宏观单元,进而形成整个结构的等效刚度矩阵,并得到宏观网格的节点位移,最后再次利用数值基函数得到微观尺度上的位移结果.该广义多尺度有限元法是扩展多尺度有限元法的一种新的拓展,可模拟具有更加复杂几何的非均质单胞的力学行为.通过数值算例,模拟了非均质材料的静力问题、广义特征值问题以及瞬态响应问题,计算结果表明:在边界条件一样的情况下,广义多尺度有限元法的计算结果与传统有限元的计算结果保持高度一致.与传统有限元相比,该方法在保证计算精度的同时极大地提高了计算效率.研究结果表明,广义多尺度有限元法能够很好地模拟非均质单胞的力学行为,具有良好的工程应用潜力.      

基于面积坐标与B网方法的四边形样条单元

传统等参元方法中, S型等参元完备阶较低,对网格畸变敏感, L型等参元具有高阶完备性但需要使用内部节点. 另外,由于引入等参变换, 采用数值积分可能导致总刚度矩阵出现奇异性.利用三角形面积坐标与B网方法建立了一类平面四边形的样条单元函数,它们的特点是满足协调条件, 克服网格畸变敏感性.其中8节点和12节点单元分别为2次和3次样条函数,对直角坐标分别具有二阶和三阶完备性, 高于相同节点的S型等参元.通过算例测试了这些样条单元, 并与等参元和其它四边形单元比较,数值结果显示了它们的高精度和有效性.      

属性与温度相关材料的广义热弹性问题有限元分析

计及材料特性与温度的相关性,基于Lord和Shulman(L-S)广义热弹性理论,建立了此类问题的有限元控制方程. 由于材料属性的温度相关性,温度控制方程具有非线性,积分变换求解方法难以采用,因而将有限元方程直接在时间域求解. 利用所建立方法研究了材料特性与温度相关、带有孔洞的无限大体在热冲击和机械冲击作用下的广义热弹性问题. 分析表明,在时间域直接求解材料属性与温度相关的广义热弹性问题是可行的,所得结果具有很高的精度,热的波动性得到充分的展现. 同时发现,热冲击载荷作用时,材料属性与温度的相关性对结构的机械响应影响显著,对温度响应影响很小;机械载荷作用时,材料参数与温度的相关性对所有响应影响都很小. 因此,研究热冲击载荷作用的机械响应时,必须考虑材料属性的温度相关性,而研究温度响应时,无论何种冲击载荷,都可以不考虑材料属性的温度相关性.      

一维有限单元的节点位移精度

研究高次杆单元和梁单元的节点位移精度问题.首先求出一端固支均匀杆和悬臂梁在任意次多项式形式分布载荷作用下的位移精确解,然后用二次杆单元、五次欧拉梁单元和三次铁木辛柯梁单元求得了节点位移.通过比较有限元解与精确解以及利用静力凝聚方法,发现一次以上杆单元、三次以上欧拉梁单元以及三次以上铁木辛柯梁单元都可以给出精确的端点位移.     

二维任意多边形有限单元

目前常用于单晶体材料的计算细观力学分析方法是将晶体细分成许多三角形或四边形有限元，对于多晶体材料，这方法的计算量之大将难于承受，实验观测表明，多晶体材料中除角点附近局部区域外每个晶体内部的变形、滑移基本上是绶慢变化的，因此可以将每个晶体简化为一个多边形有限单元，晶界简化为界面单元，使计算量大大减少，本文导出一类任意多边形等参有限元格式，并证明这类多边形单元中的任意两个都是相互协调的，无论其边数是否相等，算例表明任意多边形单元是有效可行的．