九参拟协调元的直接分析

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

Specht九参数板元的分析

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

具有几何对称性的12参数矩形板元

1 引言 三角形板元中,形式最简单的是九参数元,节点参数是单元三个顶点上的函数值和两个一阶偏导数值,非协调九参三角形板元的研究取得了丰硕成果,根据不同方法已构造出多种收敛性能好的单元.相比之下,矩形板元的研究较少见报道.矩形板元中形式最简单的是12参元,节点参数是单元4个顶点上的函数值和两个一阶偏导数值,这类似于九参三角形板元.常见的12参矩形板元是ACM元,其形函数空间是完整3次多项式空间加上两个4次多项式的基函数,ACM元是C°元,但位移形函数的外法向导数平均值在单元间不连续,这类似于Zienkiewicz九参三角形板元,但由于矩形单元的特殊形状,ACM元是收敛的.龙驭球教授等在[1]中提出一种12参矩形广义协调元,其位移形函数的外法向导数平均值在     

关于拓扑度的计算及其对于非线性算子的应用

<正> 本文首先给出了关于拓扑度计算的一个一般性定理,然后利用这个定理证明了非线性算子固有元与固有值的某些性质.作为应用,本文讨论了 Hammerstein 型非线性积分方程和二阶拟线性常微分方程两点边值问题中的固有元与固有值的性质.     

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

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

W01,p(Ω)中二阶拟线性椭圆型欧拉方程的非平凡解

本文利用Mountain Pass引理,证明了二阶拟线性椭圆型方程具有非平凡的广义解,若是偶泛函时,上述的齐次问题有无穷个解。     

十二参三角形板元

九参三角形板元的研究工作已有不少,但十二参三角形板元还较少见报道。唐立民等利用他们创立的拟协调方法构造一个十二参三角形拟协调元,节点参数是单元三个顶点上的函数值和两个一阶偏导数值及三边中点上的外法向导数值,他们是用力学方法构     

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

1.引言 最近石钟慈和陈绍春使用双参数法对九参数拟协调元和九参数广义协调元进行了研究。研究结果表明,基于函数空间     

双参数十二参矩形板元的对称列式

1 引言 在位移有限元中,九参数三角形板元的研究取得了丰硕成果,根据不同方法已构造出众多收敛性能很好的单元（见[1]、[2]、[3]）。相比之下,矩形板元的研究却较少报道,ACM元及广义协调元RGC—12是其中比较成功的单元.但是ACM是C~0元。其位移形函数的外法向导数平均值在单元间不连续。广义协调元是基于势能原理建立单元协调的,其自由度(协调条件)不对称是其本身的一个弱点,陈万吉研究表明。这种不对称性会破坏单元的几何不变性。     

关于非线性狄立克雷问题的一致有效渐近解

本文研究最高阶导数项带小参数的二阶拟线性椭圆型方程的狄立克雷问题.在退化方程不存在奇点的情形下,当参数ε是充分小时,证明了解的存在性和唯一性,并在整个区域导出解的一致有效渐近近似式.     

一种用于索结构分析的悬链线单元

根据弹性悬链线的理论解析解推导出适于索结构有限元分析的悬链线单元.与常用的三节点、五节点曲线单元相比,采用该单元编制的软件具有输入数据少、计算机时省、计算精度高的特点.     

非正则条件下类Wilson元的构造及其应用

本文在非正则性条件下，研究了窄四边形上的类Wilson元。通过参考元上类Wilson元的构造，证明了由此产生的有限元对任意窄四边形剖分通过Irons分片检查，得到了二阶问题的误差估计。结果表明，该单元的收敛性质与Wilson元的类似。     

GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS

1. IntroductionIn the numerical simulation of the Navier-Stokes equations one encounters three seriousdifficulties in the case of large Reynolds numbers f the treatment of the incomPressibility con-dition divu = 0, the treatment of the noIilinear terms and the large time integration. For thetreatment of the incoInPressibility condition, one use the penalty method in the case of finiteelemellts [1--2l and for the treatmen of the noulinar terms and the large tfor integration, oneuse the nonlin…     

Rail-bridge coupling element of unequal lengths for analysing train-track-bridge interaction systems

This paper presents a rail-bridge coupling element of unequal lengths, in which the length of a bridge element is longer than that of a rail element, to investigate the dynamic problem of train-track-bridge interaction systems. The equation of motion in matrix form is given for a train-track-bridge interaction system with the proposed element. The first two numerical examples with two types of bridge models are chosen to illustrate the application of the proposed element. The results show that, for the same length of rail element, (1) the dynamic responses of train, track and bridge obtained by the proposed element are almost identical to those obtained by the rail-bridge coupling element of equal length, and (2) compared with the rail-bridge coupling element of equal length, the proposed element can help to save computer time. Furthermore, the influence of the length of rail element on the dynamic responses of rail is significant. However, the influence of the length of rail element on the dynamic responses of bridge is insignificant. Therefore, the proposed element with a shorter rail element and a longer bridge element may be adopted to study the dynamic responses of a train-track-bridge interaction system. The last numerical example is to investigate the effects of two types of track models on the dynamic responses of vehicle, rail and bridge. The results show that: (1) there are differences of the dynamic responses of vehicle, rail and bridge based on the single-layer and double-layer track models, (2) the maximum differences increase with the increase of the mass of sleeper, (3) the double-layer track model is more accurate.     

三维Poisson方程特征值的四种有限元解及比较

应用三维EQ1rot元、三维Crouzeix-Raviart元、八节点等参数元、四面体线性元计算三维Poisson方程的近似特征值.计算结果表明:三维EQ1rot元和三维Crouzeix-Raviart元特征值下逼近准确特征值,八节点等参数元、四面体线性元特征值上逼近准确特征值,三维EQr1ot元和三维Crouzeix-Raviart元外推特征值下逼近准确特征值.计算结果还表明三维Crouzeix-Raviart元是一种计算效率较高的非协调元.     

一族新的三维体元——混合刚度有限元

对于三维连续介质的有限元分析,一个通用的二次有限单元体是所谓20节点等参元.尽管这个元素已富有成效地被普遍应用,但是,它需要的节点自由度太多与其达到的二次多项式的逼近精度却十分不相称,显得计算效率很低.基于混合刚度有限元法的一     

Infinite Groups with an Anticentral Element

An element of a group is called anticentral if the conjugacy class of that element is equal to the coset of the commutator subgroup containing that element. A group is called Camina group if every element outside the commutator subgroup is anticentral. In this paper, we investigate the structure of locally finite groups with an anticentral element. Moreover, we construct some non-periodic examples of Camina groups, which are not locally solvable.     

半导体器件瞬态模拟的对称正定混合元方法

提出具有对称正定特性的混合元格式求解非稳态半导体器件瞬态模拟问题。提出一个最小二乘混合元方法、一个新的具有分裂和对称正定性质的混合元格式和一个解经典混合元方程的对称正定失窃工格式求解电场位势和电场强度方程;提出一个最小二乘混合元格式求解关于电子与空穴浓度的非稳态对流扩散方程,浓度函数和流函数被同时求解;采用标准的有限元方法求解热传导方程。建立了误差分析理论。     

Fourteen-Node Mixed Stiffness Element and Its Computational Comparisons with Twenty-Node Isoparametric Element

In this paper, a family of 3-dimensional elements different from isoparametric serendipity is developed according to the variational principle and the convergence criteria of the mixed stiffness finite element method. For the new family, which is named mixed stiffness elements, the number of nodes on the quadratic element is not 20 but 14. Theoretical analysis and various computational comparisons have found the mixed stiffness element superior over the isoparametric serendipity element, especially a substantial improvement in computational efficiency can be achieved by replacing the 20 node-isoparametric element with the 14-node mixed stiffness element.     

A postprocessed flux conserving finite element solution

We propose a local postprocessing method to get a new finite element solution whose flux is conservative element‐wise. First, we use the so‐called polynomial preserving recovery (postprocessing) technique to obtain a higher order flux which is continuous across the element boundary. Then, we use special bubble functions, which have a nonzero flux only on one face‐edge or face‐triangle of each element, to correct the finite element solution element by element, guided by the above super‐convergent flux and the element mass. The new finite element solution preserves mass element‐wise and retains the quasioptimality in approximation. The method produces a conservative flux, of high‐order accuracy, satisfying the constitutive law. Numerical tests in 2D and 3D are presented.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1859–1883, 2017