基于解析试函数的广义协调超基膜元

利用解析试函数法构造一个带旋转自由度广义协调超基膜元。根据弹性力学平面问题的控制方程和艾雷应力函数,求出问题完备的基本解析解,然后用其作为试函数并采用广义协调条件来构造单元:ATF-SCQ4θ。数值算例表明,该类单元精度高、对网格畸变不敏感,显示出良好的性能。     

用于壳体自由振动分析的广义协调元

本文采用基于解析试函数的广义协调四边形膜单元和中厚板单元构造了平板型4节点壳体单元,并将其用于壳体振动分析。该壳体单元具有列式简单,易于编程的优点,通过数值算例表明,该单元计算精度高,非常适合工程计算     

应用基于解析试函数的广义协调四边形厚板元分析中厚板的自由振动

自从广义协调元提出以来,多种性能优异的板单元被构造且得到广泛应用,但是将广义协调元用于板的动力分析目前仅局限于薄板,该文将基于解析试函数的广义协调四边形厚板元ATF—PQ4b用于板的动力分析。该单元是将静力条件下中厚板的齐次微分方程的多项式解析解作为试函数,从而推导得到单元刚度矩阵和质量矩阵。笔者编写了相关的MATLA...     

基于解析试函数的各向异性材料厚薄通用板单元

该文采用满足Kirchhoff假设的薄板理论,推导了各向异性材料系列解析试函数,并利用该系列解析试函数构造了一个四边形应力杂交板单元。首先,该文从薄板理论的基本方程出发,推导了各向异性材料薄板中面挠度w应满足的特征微分方程。然后,从该方程出发求得w的系列特征通解,由w特征通解可进一步求得广义位移、广义应变和广义应力的解析试函数。同时,根据广义应力利用平衡条件构造了相应的横向剪力解析试函数。最后,根据已有的广义应力和横向剪力解析试函数构造了一个四边形应力杂交板单元ATF-PH4。数值算例表明:上述方法构造出的单元模型有很好的精度、收敛性,且对网格畸变不敏感,同时能较好地解决板单元的厚薄通用性问题。     

基于四边形面积坐标的广义协调轴对称单元

利用了点组合广义协调和周广义协调条件,基于四边形面积坐标方法构造了含内参的4结点四边形空间轴对称单元AQACQ6。通过进一步对内参应变矩阵进行合理修正,从而形成新单元AQACQ6M,该单元能够通过强式分片检验。两种单元的位移场都达到对整体坐标的二次完备。数值算例表明:上述轴对称单元不仅精度高,而且抗网格畸变和几乎不可压缩问题能力优于等参单元,显示了面积坐标和广义协调理论的优越性。     

广义协调六结点平面曲边单元研究

主要运用广义协调原理,针对计算平面曲边单元的有限元算法进行了研究,并且利用点、周混合协调条件构造了三种高性能六结点曲边单元。第一、二种单元在平面直角坐标内分别采用解析试函数和完全三次多项式构造,第三种单元在六结点等参单元Q6的基础上附加广义协调泡状位移而成。这三种单元均能通过强式分片试验,并且显示了良好的计算精度和抗畸变能力。     

采用面积坐标的四边形板弯曲单元

本文采用四边形面积坐标，并应用广义协调法构造出一个具有１２个自由度的四边形板弯曲单元。单元的挠度场以面积坐标多项式表示，对应于直角坐标ｘ，ｙ的完全三次式和部分四次式，因而单元是完备的广义协调的板单元。应用的１２个协调条件为挠度的四个点协调条件和四个边协调条件，以及法向转角的四个边协调条件。由于面积坐标和直角坐标之间为线性变换关系，因此单元刚度矩阵的推导相当简单。数值算例表明：本文单元具有高精度、收敛性、可靠性和对网格畸变不敏感的优点     

采用SemiLoof型约束条件的薄板矩形广义协调元

本文采用SemiLoof型约束条件,建立一个十二自由度的薄板矩形广义协调元。单元自由度只含角点位移,不含Loof结点位移、单元间的协调条件全部采用点型协调条件,不采用积分型协调条件。此单元吸取广义协调元和SemiLoof元的双重优点,消除其缺点,成为同类低阶薄板单元中的最优单元。     

弹性力学有限元分析中的平衡与协调理论

有限元分析中的单元可以遵循不同的方法构造,该文提出以单元模型的平衡性与协调性进行分类,并对弹性力学平面问题中的几种经典单元进行了分析比较,总结了协调元、非协调元和超协调元的协调性方法,以及基于解析试函数法的平衡型方法。单元的协调性理论思路包含单纯形格式协调元和非单纯形格式协调元,以及相应的非协调和超协调元格式,关注的重点是单元边界的协调。单元的平衡性理论思路包含解析试函数和权函数的高阶完备性,关注的重点研究是单元内部及边界平衡性。研究表明：针对弹性力学中平衡性和协调性要求,两类理论给出的不同单元格式各具特点,而既能保证单元内部平衡性,又能考虑单元界面协调性的单元类型给出了更精确、合理的计算结果。     

带旋转自由度的广义协调三角形膜元

应用广义协调元理论,通过增加内参位移的方法,构造了两个具有旋转自由度的三角形膜元。本文单元能通过任意三角形分片检验,符合连续介质力学关于旋转自由度的定义,并且没有引入刚性转角假设。数值算例表明其具有较高的计算精度。     

基于解析试函数的广义协调四边形膜元

本文构造三个广义协调四边形膜元。根据弹性力学平面问题的控制方程和艾雷应力函数,求出问题基本解析解,然后用其作为试函数来构造单元:ATF-Q4a、ATF-Q4b、ATF-Q4q。数值算例表明,其中两个单元ATF-Q4a和ATF-Q4q对网格畸变不敏感,显示出良好的性能。     

基于解析试函数的广义协调四边形厚板元

本文构造两个广义协调四边形厚板元ATF-PQ4a和ATF-PQ4b。根据Mindlin-Reissner厚板理论的控制方程，首先求出其基本解析解，然后用其作为试函数来构造单元。数值算例表明，这两个单元不出现剪切闭锁，显示出良好的性能。     

An incompatible and unsymmetric four-node quadrilateral plane element with high numerical performance

Recent studies show that the unsymmetric finite element method exhibits excellent performance when the discretized meshes are severely distorted. In this article, a new unsymmetric 4-noded quadrilateral plane element is presented using both incompatible test functions and trial functions. Five internal nodes, one at the elemental central and four at the middle sides, are added to ensure the quadratic completeness of the elemental displacement field. Thereafter, the total nine nodes are applied to form the shape functions of trial function, and the Lagrange interpolation functions are adopted as the incompatible test shape functions of the internal nodes. The incompatible test displacements are then revised to satisfy the patch test. Numerical tests show that the present element can provide very good numerical accuracy with badly distorted meshes. Unlike the existing unsymmetric four-node plane elements in which the analytical stress fields are employed, the present element can be extended to boundary value problems of any differential equations with no difficulties.     

Quadrilateral membrane elements with analytical element stiffness matrices formulated by the new quadrilateral area coordinate method (QACM‐II)

A novel strategy for developing low‐order membrane elements with analytical element stiffness matrices is proposed. First, some complete low‐order basic analytical solutions for plane stress problems are given in terms of the new quadrilateral area coordinates method (QACM‐II). Then, these solutions are taken as the trial functions for developing new membrane elements. Thus, the interpolation formulae for displacement fields naturally possess second‐order completeness in physical space (Cartesian coordinates). Finally, by introducing nodal conforming conditions, new 4‐node and 5‐node membrane elements with analytical element stiffness matrices are successfully constructed. The resulting models, denoted as QAC‐ATF4 and QAC‐ATF5, have high computational efficiency since the element stiffness matrices are formulated explicitly and no internal parameter is added. These two elements exhibit excellent performance in various bending problems with mesh distortion. It is demonstrated that the proposed strategy possesses advantages of both the analytical and the discrete method, and the QACM‐II is a powerful tool for constructing high‐performance quadrilateral finite element models. Copyright © 2008 John Wiley & Sons, Ltd.     

Implicit boundary method for finite element analysis using non‐conforming mesh or grid

In the finite element method (FEM), a mesh is used for representing the geometry of the analysis and for representing the test and trial functions by piece‐wise interpolation. Recently, analysis techniques that use structured grids have been developed to avoid the need for a conforming mesh. The boundaries of the analysis domain are represented using implicit equations while a structured grid is used to interpolate functions. Such a method for analysis using structured grids is presented here in which the analysis domain is constructed by Boolean combination of step functions. Implicit equations of the boundary are used in the construction of trial and test functions such that essential boundary conditions are guaranteed to be satisfied. Furthermore, these functions are constructed such that internal elements, through which no boundary passes, have the same stiffness matrix. This approach has been applied to solve linear elastostatic problems and the results are compared with analytical and finite element analysis solutions to show that the method gives solutions that are similar to the FEM in quality but is less computationally expensive for dense mesh/grid and avoids the need for a conforming mesh. Copyright © 2007 John Wiley & Sons, Ltd.