关于九参数二阶拟协调元

韩厚德最近在讨论拟协调元时,引进了一个九参数二阶拟协调元.在[1]中,对一个完全三次多项式(十个自由度)附加一个特殊的约束条件并使它满足所谓的二阶拟协调条件,即形函数及其两个一阶偏导数在单元之间的内部边界上保持积分意义下的连续性.我们证明,这样得到的九参二阶拟协调元实际上就是熟知的 de Veubeke元.其证明如下.     

关于构造单元刚度矩阵的“自由列式”方法

[1]中提出的单片检验IET(Individual Element Test)用于构造有限元的单元刚度矩阵.随后,此方法逐步改进,又有了构造单元刚度矩阵的“自由列式”方法.此公式将单元刚度矩阵对应常应变部分与高阶模态分别构造,没有交叉项,常应变部分是恒不变的.     

Wilson非协调元的V循环多重网格法 *

非协调有限元V循环多重网格法的收敛性至今仍是一个没有很好解决的问题 .给出了Wilson非协调有限元的两类V循环多重网格法的收敛性证明 .     

一类非协调元的收敛性分析

１．引言最近,［１］在求解Ｓｔｏｋｅｓ方程时提出了一类非协调元．这类非协调元定义在矩形网格上,形式简单,自由度少,比熟知的多线性元还少．例如,对于三维问题,它们只有六个自由度,是型函数在各个表面中点处的函数值或者表面上的平均值．与一些熟知的非协调元不同之处是,由于这类非协调元简单的形式及其非常少的自由度,它们不包含任何协调的部分,但是,已有的结果表明它们具有良好的数值表现［１］,［２］中已有一些理论分析,用它们来求解晶体的微结构方程,也取得较好的结果．本文将详细分析这类非协调元的收敛性,所用的工…     

Specht九参数板元的分析

Zienkiewicz 九参数的不协调三角形板元以单元顶点上的函数值和两个一阶偏导数值作为节点参数,形式简单,总体自由度少,但它只对特殊剖分才收敛,不便于实际应用.后来相继出现一些对 Zienkiewicz元的改进形式,如拟协调元,TRUNC元,广义协调元等.这是一类非常规的有限元,均为三结点九参数三角形元.从常规有限元     

关于九参数双参数元与广义协调元的对称列式

关于九参数双参数元与广义协调元的对称列式石钟慈（中国科学院计算数学与科学工程计算研究所）石东洋（西安交通大学）ＯＮＴＨＥＳＹＭＭＥＴＲＩＣＡＬＦＯＲＭＯＦＴＨＥ９－ＰＡＲＡＭＥＴＥＲＤＯＵＢＬＥＳＥＴＰＡＲＭＥＴＥＲＥＬＥＭＥＮＴＡＮＤＴＨＥＧＥＮＥ...     

一般组合弹性结构的有限元分析

提出一个求解一般组合弹性结构问题的有限元方法:对体件的位移、 板件的纵向位移和杆件的纵向位移用线性协调元离散, 对杆件的纵向转角用二次协调元离散, 对板件的横向位移和杆件的横向位移分别用Morley元和三次Hermite元离散. 在相应的非协调元空间上建立了广义Korn不等式, 进而证得有限元方法的唯一可解性. 最后导出方法在能量模意义下的最优误差估计.     

Reissner-Mindlin板问题带约束非协调旋转Q_1有限元方法

本文利用带约束非协调旋转Q_1元逼近Reissner-Mindlin板问题中旋度的两个分量.并分别选择Wilson元、双线性元和带约束非协调旋转Q_1元逼近挠度,相应地选取不连续的矢量值分片线性函数空间、最低阶旋转Raviart-Thomas元空间和矢量值分片常数函数空间为离散的剪应力空间,在矩形网格上构造了三个板元.通过证明一个离散的Korn不等式,并借助MITC4元的解构造了旋度、挠度和剪应力一个具有某种特殊且关键的可交换性的插值.再利用Helmholtz分解分析相容性误差.我们证明了这三个矩形元在能量范数意义下与板厚无关的一致最优收敛性.数值算例验证了我们的理论结果.     

两种不协调板元的统一型函数形式

一、引言 在有限元法中,选取坐标合适,会给理论分析和实际计算带来不少方便.对于矩形单元,一般采用直角坐标(总体的或局部的);对于三角形元,则采用三角形的面积坐标.     

九参拟协调元的直接分析

§1.引言 几年前,唐立民等提出一种构造弹性力学方程离散格式的非常规有限元方法,称之为拟协调元法.用这种方法构造单元刚度阵简单灵活,并有良好的数值精度.张鸿庆在[3]中首先对九参拟协调板元进行了理论分析,证明这个非常规板元实际上等价于一