抛物方程时间连续有限元全离散格式的超收敛性

利用张量积分解和时间方向单元正交分解,证明了线性抛物型方程的时间连续全离散有限元在单元节点和内部的特征点的超收敛性.并用连续有限元计算了非线性Schrodinger方程,验证了能量的守恒性.计算结果与理论相吻合.     

二个非协调膜元的收敛问题

为了用有限元方法求(0-2)的近似解,定义有限元空间.记P(K),?K∈?_h,为单元K上多项式组成的有限维空间,v_h∈P(K)可由K上的节点参数或其它类型的参数(例如函数的导数或函数本身在单元K上积分的数值)唯一决定.     

抛物方程初边值问题连续有限元的超收敛性

研究了一类一维抛物方程初边值问题的连续有限元方法.在空间上进行任意m次有限元半离散,在时间方向上进行二次连续有限元后,获得了一个稳定的全离散计算格式.利用单元分析法校正技术的新思想进行理论分析,连续有限元解在剖分网格节点上具有超收敛性.     

节点应力连续的四边形单元

节点应力连续的四边形单元Q4-CNS是一种基于单位分解理论的混合的有限元无网格法.Q4-CNS可以视作FE-LSPIM QUAD4的发展.Q4-CNS形函数的导数在节点处是连续的,因此可以自然的得到节点应力,而不需要使用节点应力磨平算法.数值实验表明,与传统四边形单元(QUAD4)相比,Q4-CNS具有更好的计算精度和更高的收敛速度.在扭曲网格下,Q4-CNS也能取得满意的数值精度.然而,QUAD4的数值精度则会随着网格的扭曲明显的变差.基于Kirchhoff-Love假设的非协调板单元计算中,不仅要求形函数在单元的交界面上要保持C0连续性,而且要求形函数在节点处具有C1连续性,所以在任意的四边形单元上构造满足插值条件的非协调板单元形函数较为困难.Q4-CNS形函数的导数在节点处是连续的,所以Q4-CNS在求解基于Kirchhoff-Love假设的板单元问题中具有潜在的应用价值.     

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

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

应用求解拉普拉斯方程的边值问题建立有限元网格

一、引言 应用有限元方法求解各种数学物理问题时要求对求解区域进行有限的剖分,然后在网格的节点上将连续的问题离散.如所周知,给出网格节点的坐标及其有关单元的信息是     

空间各向异性弹性问题的二十节点理性单元

常规有限元方法的插值函数通常仅仅从数学层面上考虑单元的几何性质,忽视了与物理问题相关的物性参数,因此可能降低数值分析的效果.理性有限元的构造方法与常规有限元法不同,其插值函数使用的是控制微分方程解析解的线性组合,求解过程是在物理域内直接列式,对单元的应变场和应力场同时进行插值,并在单元级别考虑分片实验的要求并直接进行修正,最终形成与问题物性参数紧密相关的单元刚度阵.该方法避免了传统方法对物理问题和数学问题的割裂,可显著提高数值分析的稳定性和精度.利用空间各向异性问题的基本解,从最小势能原理出发,构造出两种满足分片实验要求的二十节点理性块体单元.数值算例表明,所给出的理性单元不仅具有较高的求解精度,而且具有良好的数值稳定性.     

轴对称Poisson方程的Trefftz有限元解法

将径向基函数应用到一类轴对称Poisson方程的数值求解中,提出了一种Trefftz有限元计算格式.非0右端项将问题的特解引入Trefftz单元域内场,致使单元刚度方程涉及区域积分.利用径向基函数对特解近似处理,可消除区域积分,从而保持Trefftz有限元法只含边界积分的优势.为获得特解,选取求解域内所有单元的节点和形心作为基本插值点,而在求解域之外构造一个虚拟边界,在其上布置一定数目的虚拟点作为额外插值点.数值算例验证了该方法的有效性和可行性.     

有限元混合法自动网格划分技术分析贴边岔管

本文应用有限元混合法分析水电站贴边岔管于内水压力作用时的应力问题.岔管系当作为二个薄壳相交的组合壳体.由于需要计算大量的岔管,发展了自动网格划分技术用于有限元混合法.根据三个参数即二个交叉管道的轴线交角、管壁厚度及交叉口二管道直径之比,区分出八种网格,十二种单元信息,编出自动确定节点坐标及形成边界条件信息的程序,业已计算了101个岔管,计算效果优良.     

二阶常微分方程初值问题C^0有限元的超收敛

基于在一个单元上的改进的单元正交展开，得到二阶常微分方程C^0有限元Uh的误差表示式，以及节点和单元内部一些特征点上的超收敛结果。     

A discrete divergence-free basis for finite element methods

The divergence-free finite element method (DFFEM) is a method to find an approximate solution of the Navier–Stokes equations in a divergence-free space. That is, the continuity equation is satisfied a priori. DFFEM eliminates the pressure from the calculations and significantly reduces the dimension of the system to be solved at each time step. For the standard 9-node velocity and 4-node pressure DFFEM, a basis for the weakly divergence-free subspace is constructed such that each basis function has nonzero support on at most 4 contiguous elements. Given this basis, weakly divergence-free macroelements are constructed. This revised version was published online in August 2006 with corrections to the Cover Date.     

三角形REISSNER-MINDLIN板元

本文提出构造无自锁现象的Reissuer－Mindlin板元的一个一般性方法.此方法将剪切应变用它的适当的插值多项式代替,当板厚趋于零时这对应于插值点的Kirchhoff条件,因而单元无自锁现象.根据这种方法我们构造两个三角形元──一个3节点元和一个6节点元,并给出数值结果.     

Accuracy assessment to applying 20-node solid elements to pressurized composite shells

A factorial analysis is conducted to explore the impact of five factors — aspect ratio, radius-to-thickness ratio, material anisotropy, fiber angle and integration order on the finite element accuracy, with regard to the application of 20-node isoparametric solid elements to clamped composite cylindrical shells subjected to internal pressure. The maximum radial displacement at the central cross-section is employed as the indicator of accuracy. Results from finite element methods were compared with closed-form solutions. Aspect ratio, radius-to- thickness ratio, and integration order are significant main effects as predicted for isotropic materials. Material anisotropy shows no significant main effect, while the fiber angle does have strong impact on the finite element accuracy. In addition to four influential main effects, some multiple-factor interactions are shown to have significant influences on the finite element accuracy. The analysis is limited to the linear elastic range.     

采用三角形面积坐标的四边形17节点样条单元

利用二元4次样条插值基和三角形面积坐标构造17节点四边形单元．这个新单元具有4次完备阶,通过一些算例测试表明了该单元有较高精度并对网格畸变不敏感．     

局部加密的非标准混杂五点矩形元分析

In this paper we apply the nonconforming theory in Z.C.Shi to construct a hybrid 5-node finite element satisfying the F1-patch test. We show that the hybrid 5-node finitc element for two-order elliptic equation is convergent and obtain an error estimate.     

Geometrically nonlinear analysis with a 4-node membrane element formulated by the quadrilateral area coordinate method

Recently, a 4-node quadrilateral membrane element AGQ6-I, has been successfully developed for analysis of linear plane problems. Since this model is formulated by the quadrilateral area coordinate method (QACM), a new natural coordinate system for developing quadrilateral finite element models, it is much less sensitive to mesh distortion than other 4-node isoparametric elements and free of various locking problems that arise from irregular mesh geometries. In order to extend these advantages of QACM to nonlinear applications, the total Lagrangian (TL) formulations of element AGQ6-I was established in this paper, which is also the first time that a plane QACM element being applied in the implicit geometrically nonlinear analysis. Numerical examples of geometrically nonlinear analysis show that the presented formulations can prevent loss of accuracy in severely distorted meshes, and therefore, are superior to those of other 4-node isoparametric elements. The efficiency of QACM for developing simple, effective and reliable serendipity plane membrane elements in geometrically nonlinear analysis is demonstrated clearly.     

p-version hybrid/mixed finite element method

In this paper, the p-version hybrid/mixed finite element formulation using the Hu–Washizu principle is presented. The shape functions are assumed to be hierarchical for the displacement variables. The stresses and the strains are interpolated using Legendre polynomials so that the computation of the stiffness matrix will be simplified. For 2-D problems, the higher-order hierarchical shape functions start following the 8-node, quadrilateral shape functions. The geometry mapping of each element is also performed using an 8-node parametric mapping. Numerical examples for a plate with a circular hole and circular inclusion are included.     

一般壳体组合结构在均布压力下的平衡与稳定

本文采用几何非线性理论建立一般壳体组合结构的能量泛函,再根据能量极值原理提出了这类组合结构在均布压力下的平衡和稳定的非线性有限元计算方法。计算结果与实验结果符合较好。     

Linear and Nonlinear Finite Element Analysis of Active Composite Laminates

The paper presents a finite element concept for analysis of thin-walled active structures featuring fiber reinforced composite laminate as a passive structural material. The structure is rendered active by embedding piezoelectric material as a multifunctional material. A 9-node degenerated shell element based on the first order shear deformation theory is developed as a modelling tool capable of predicting the general behavior of the structure for controlling purposes. The von-Kármán type nonlinearities are considered. The solution strategy of the geometrical nonlinear analysis is based on the incremental approach using the updated Lagrangian formulation. Some numerical results are given to demonstrate the behavior of the element. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical integration over polygons using an eight-node quadrilateral spline finite element

In this paper, a cubature formula over polygons is proposed and analysed. It is based on an eight-node quadrilateral spline finite element [C.-J. Li, R.-H. Wang, A new 8-node quadrilateral spline finite element, J. Comp. Appl. Math. 195 (2006) 54–65] and is exact for quadratic polynomials on arbitrary convex quadrangulations and for cubic polynomials on rectangular partitions. The convergence of sequences of the above cubatures is proved for continuous integrand functions and error bounds are derived. Some numerical examples are given, by comparisons with other known cubatures.