采用面积坐标和基于假设转角的薄板元

采用四边形面积坐标方法,从假设转角位移场入手构造了两个广义协调四边形4结点薄板单元AΨQ-I和AΨQ-II。通过采用边界协调条件一次项与二次项分别协调使转角场实现了三次完备。与DKQ等同类单元相比,单元的精度和抗网格畸变能力都有很大提高。     

基于解析试函数的内参型广义协调膜元

利用解析试函数法构造一个内参型四结点八自由度广义协调膜元。根据弹性力学平面问题的控制方程和艾雷应力函数,求出问题完备的基本解析解,然后用其作为试函数并采用广义协调条件来构造单元:ATF-GCQ4X。该单元采用了14个解析试函数构造了应变二次完备的内部场,同时引入6个附加边界位移模式,采用平衡力系为权函数构造相应的广义协调条件。数值算例表明,该类内参型单元能在不提高单元结点自由度的情况下提高单元精度,并显示出良好的收敛特性。     

一个有效的厚薄板通用COONS曲面矩形元

本文采用板壳有限元Coons曲面法和挠度位移与转角位移导数相关构造法构造出一个有效的厚薄板通用矩形单元,并直接给出位移场形函数;其推导简单易行,灵活多样,几何意义直观,力学概念清晰,并具有一般性。在薄板极限情况下,本文厚薄板通用单元自动退化为一个性能良好的薄板Coons曲面矩形元,不出现剪切闭锁现象。数值计算表明:本文构造的厚薄板通用单元具有较高的计算精度和良好的收敛性。     

采用Coons曲面法构造薄板广义协调元

本文应用Coons曲面法构造出一种广义协调薄板矩形单元。用此单元分析薄板弯曲和振动问题都具有良好的精度。     

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

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

具有旋转自由度的广义协调矩形膜元

本文根据广义协调的思想,在平面应力矩形单元双线性协调位移场中,引入附加广义泡状位移场,构造出一种具有平面内旋转自由度的矩形膜单元,它满足广义协调条件。数值计算结果表明,这种单元有很高的计算精度,而且计算量少,是一种能收敛于精确解的单元。     

铰接板机构运动分析的简便协调方程

一些新型空间结构的施工分析模型可以简化为铰接板机构,如Pantadome。利用三角形形状稳定性特点,将三条边长变化为零来表征三角形板单元的刚体位移。进而基于杆单元的协调方程,建立了一个简便的顶点铰接三角形板单元的机构位移协调方程,且给出了1阶和2阶协调矩阵。将四边形板单元划分为2个三角形板单元,通过引入单元四顶点共面条件,推导出平面四边形板单元的协调矩阵。理论上,利用该思路可构建任意平面多边形板单元的协调矩阵。针对该简便协调方程,进一步给出了求解铰接板机构运动路径的计算策略。对1个顶升施工的Pantadome和1个顶推施工的双坡网架的成形过程进行了数值模拟,结果表明该方法对于此类铰接板机构的运动路径分析有很高的精度。     

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

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

复合材料夹层板振动分析的精化锯齿理论和三角形板单元

基于精化锯齿理论,构造了六节点三角形协调板单元并推导了夹层板自由振动问题有限元列式。不同于已有锯齿理论,精化锯齿理论特点是面内位移不含有横向位移一阶导数,构造有限元时仅需要C0 插值函数。为验证单元性能,分析了软核夹层板自由振动问题。结果表明,该文构造的单元能准确计算软核夹层板固有频率,然而基于已有锯齿理论建立的不协调元计算结果精度较低。     

引入广义泡状型位移场影响的厚板矩形单元

本文提出了在已有的无闭锁现象的广义协调厚板矩形单元的基础之上，引入广义泡状型位移场的影响，构造厚薄板通用的优质单元的一般方法．应用此法构造出新的单元．数值算例表明，在引入适当的广义泡状位移影响后，单元的计算结果显著得到改善，精度进一步提高．     

Eight-node conforming straight-side quadrilateral element with high-order completeness (QH8-C1)

A new eight-node conforming quadrilateral element with high-order completeness, denoted as QH8-C1, is proposed in this article. First, expressions for the interpolation displacement function satisfying the requirements for high-order completeness in the global coordinate system are constructed. Second, the displacement function expression in global coordinates is transformed into isoparametric coordinates, and the relationships between the two series of coefficients for the two kinds of displacement function expressions are found. Third, the displacement function expression is modified to satisfy the requirements of nodal freedom and interelement boundary continuity. The key to the new element construction is the derivation of the linear relationship expressions among 12 coefficients of element displacement interpolation polynomials in the global and isoparametric coordinate systems. As a result, the relationship between quadratic completeness and interelement continuity is explicitly given, and a proof of the completeness and the continuity was conducted to theoretically guarantee the validity of the derivation results. Furthermore, in order to verify the correctness of the theoretical work, nine numerical examples were performed. The computation results from these examples demonstrate that QH8-C1 exhibited excellent performance, including high simulation accuracy, fast convergence, insensitivity to mesh distortion, and monotonic convergence.     

Variational principles for non‐conforming finite element methods

Variational principles with relaxed inter‐element continuity requirement for non‐conforming element methods in linear and non‐linear analyses are developed. Based on the principles, any non‐conforming element displacement can be used directly to derive the explicit expressions of non‐conforming displacement function, which can ensure the passage of the patch test C for the requirement of convergence Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

三维弹性体移动接触边界元法的一类新方案

基于对二维弹性体移动和滚动接触边界元法的前期研究,将其中的协调离散方案推广到三维问题,提出了针对给定移动方向的三维弹性体移动接触的一种边界元协调离散方案。其中在接触面上的位移和面力都能在边界元离散意义下精确满足,因此能保持边界元法在位移计算、特别是应力计算中高精度的优势。该方案将可能接触区分成大小相等的直角三角形单元,每个直角三角形单元有6个节点,其中3个角节点是固定节点,另外在每条边上有一可动节点。可动节点的位置确定于与其接触的另一面上对应单元的固定节点及其连线的位置。在每一瞬时,每个单元可由可动节点连线分成4个三角形子单元,在每个子单元上边界位移与面力线性分布。文中给出了一些算例来验证所提出算法的有效性和高精度。     

Refined non‐conforming triangular elements for analysis of shell structures

Based on the refined non‐conforming element method, simple flat triangular elements with standard nodal displacement parameters are proposed for the analysis of shell structures. For ensuring the convergence of the elements a new coupled continuity condition at the inter‐element has been established in a weaker form. A common displacement for the inter‐element, an explicit expression of refined constant strain matrix, and an adjustable constant are introduced into the formulation, in which the coupled continuity requirement at the inter‐element is satisfied in the average sense. The non‐conforming displacement function of the well‐known triangular plate element BCIZ [1] and the membrane displacement of the constant strain triangular element CST [2] are employed to derive the refined flat shell elements RTS15, and the refined flat shell elements RTS18 is derived by using the element BCIZ and the Allman's triangular plane element [3] with the drilling degrees of freedom. A simple reduced higher‐order membrane strain matrix is proposed to avoid membrane locking of the element RTS18. An alternative new reduced higher‐order strain matrix method is developed to improve the accuracy of the elements RTS15 and RTS18. Numerical examples are given to show that the present methods have improved the accuracy of the shell analysis. Copyright © 1999 John Wiley & Sons, Ltd.     

采用广义协调条件构造具有旋转自由度的四边形膜元

本文根据广义协调的概念,通过引进单元结点刚体转角,提出两种具有平面内旋转自由度的四边形膜元。单元列式简单,是能通过任意四边形分片检验的收敛单元。数值计算表明这两种单元无论是位移还是应力都有很高的计算精度。