An exact element method for plane problem

In this paper, based on the step reduction method^[1] and exact analytic method^[2] anew method-exact element method for constructing finite element, is presented. Since the new method doesn ’t need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a quadrilateral noncompatible element with 8 degrees of freedom is derived for the solution of plane problem. Since Jacobi ’s transformation is not applied, the present element may degenerate into a triangle element. It is convenient to use the element in engineering. In this paper, the convergence is proved. Numerical examples are given at the endof this paper, which indicate satisfactory results of stress and displacements can be obtained and have higher numerical precision in nodes.     

一个新的八节点非协调曲边四边形平面单元

本文在文[1]和文[2]的基础上,提出构造非协调有限元的新方法。该方法不用一般的变分原理,可适用任意变系数正定和非正定偏微分方程。利用这一方法得到一个新的八节点四边形平面应力单元。与一般有限元相比,位移和应力可提高一阶收敛精度。形成单刚矩阵时,不需要进行数值积分。单元之间的协调条件容易满足,文中给出收敛性证明。文末给出数值算例,表明利用本文的方法,应力和位移均可获得满意的数值精度。     

Incompatible curved quadrilateral plate bending element with 12-degrees of freedom

This paper presents a new curved quadrilateral plate element with12-degree freedom by the exact element method.The method can be used to arbitrary non-positive and positive definite partial differential equations without variation principle.Using this method,the compatibility conditions between element can be treated very easily,if displacements and stress resultants are continuous at nodes between elements.The displacements and stress resultants obtained by the present method can converge to exact solution and have the second order convergence speed.Numerical examples are given at the end of this paper,which show the excellent precision and efficiency of the new element.     

AN EXACT ELEMENT METHOD FOR THE BENDING OF NONHOMOGENEOUS REISSNER’S PLATE

In this paper,based on the step reduction method and exact analytic method,a new method,the exact element method for constructing finite element,is presented.Since the new method doesn’t need variational principle,it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficients.By this method,a triangle noncompatible element with15 degrees of freedom is derived to solve the bending of nonhomogenous Reissner’s plate.Because the displacement parameters at the nodal point only contain deflection and rotation angle.it is convenient to deal with arbitrary boundary conditions.In this paper,the convergence of displacement and stress resultants is proved.The element obtained by the present method can be used for thin and thick plates as well,Four numerical examples are given at the end of this paper,which indicates that we can obtain satisfactory results and have higher numerical precision.     

General solution for large deflection problems of nonhomogeneous circular plates on elastic foundation

In this paper, a new method is presented based on [1]. It can be used to solve the arbitrary nonlinear system of differential equations with variable coefficients. By this method, the general solution for large deformation of nonhomogeneous circular plates resting on an elastic foundation is derived. The convergence of the solution is proved. Finally, it is only necessary to solve a set of nonlinear algebraic equations with three unknowns. The solution obtained by the present method has large convergence range and the computation is simpler and more rapid than other numerical methods.Numerical examples given at the end of this paper indicate that satisfactory results of stress resullants and displacements can be obtained by the present method. The correctness of the theory in this paper is, confirmed.     

The general solution for axial symmetrical bending of nonhomogeneous circular plates resting on an elastic foundation

In this paper,a new method,the exact analytic method,is presented on the basis of stepreduction method.By this method,the general solution for the bending of nonhomogenouscircular plates and circular plates with a circular hole at the center resting,on an elastfcfoundation is obtained under arbitrary axial symmetrical loads and boundary conditions.The uniform convergence of the solution is proved.This general solution can also be applieddirectly to the bending of circular plates without elastic foundation.Finally,it is onlynecessary to solve a set of binary linear algebraic equation.Numerical examples are givenat the end of this paper which indicate satisfactory results of stress resultants anddisplacements can be obtained by the present method.     

A new integral equation formulation of two-dimensional inclusion–crack problems

A new integral equation formulation of two-dimensional infinite isotropic medium (matrix) with various inclusions and cracks is presented in this paper. The proposed integral formulation only contains the unknown displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks. In order to solve the inclusion–crack problems, the displacement integral equation is used when the source points are acting on the inclusion–matrix interfaces, whilst the stress integral equation is adopted when the source points are being on the crack surfaces. Thus, the resulting system of equations can be formulated so that the displacements on the inclusion–matrix interfaces and the discontinuous displacements over the cracks can be obtained. Based on one point formulation, the stress intensity factors at the crack tips can be achieved. Numerical results from the present method are in excellent agreement with those from the conventional boundary element method.     

An exact analytic method applied to nonhomogeneous ring- and stringer-stiffened cylindrical shells

In this paper, the nonlinear axial symmetric deformation problem of nonhomogeneous ring- and stringer-stiffened shells is first solved by the exact analytic method. An analytic expression of displacements and stress resultants is obtained and its convergence is proved. Displacements and stress resultants converge to exact solution uniformly. Finally, it is only necessary to solve a system of linear algebraic equations with two unknowns. Four numerical examples are given at the end of the paper which indicate that satisfactory results can be obtained by the exact analytic method.     

Hahn-Tsai复合材料的非线性杂交应力有限元方法

成功建立了Hahn-Tsai复合材料模型的非线性杂交应力有限元方程,采用Newton-Raphson迭代法求解结构的非线性位移方程。在迭代过程中,为了提高计算效率可采用简单迭代法由节点位移求解单元应力场。但是,当载荷增加到一定程度以后,非线性应力场由于循环迭代而无法收敛,显然,一般的加速方法不能解决这种循环迭代的发散问题。因此,本文发展了一种确实有效的非线性应力场迭代新方法,在不增加计算工作量的情况下,不仅极大地提高了收敛速度,而且对于较大载荷也能够很好地收敛,从而解决了大载荷下非线性杂交元方法失败的关键问题。数值算例表明该方法是确实可行的。     

The general solution on nonlinear axial symmetrical deformation of nonhomogeneous cylindrical shells

In this paper, the general solution on nonlinear axial symmetrical deformation of nonhomogeneous cylindrical shells is obtained by step reduction method⁽¹⁾. The general formula of displacements and stress resultants, which is used to solve the bending problems of nonhomogeneous cylindrical shells under arbitrary axial symmetric loads, is derived. Its uniform convergence is proved. Finally, it is only necessary to solve one set of binary linear algebraic equations. A numerical example is given at the end of the paper which indicates satisfactory results of displacement and stress resultants can be obtained and converge to the exact solution.     

非均质问题中的无网格边界单元法

介绍了一种不需要内部网格计算非均匀介质问题的边界元算法.该算法是建立在一种能将任何区域积分转换成边界积分的径向积分转换法基础上,首先用对应各向同性问题的基本解来建立以正规化位移表示的非均质问题的积分方程,然后用径向积分转换法将出现在积分方程中的区域积分转换成边界积分,从而形成不需要使用内部网格来计算区域积分的纯边界元算法.与其它无网格法相比,此方法需要很少的内部点,有些问题甚至不需要内部点都能得到满意的结果,因此,可以计算大型的三维非均匀介质工程问题.由于此方法继承了边界元和无网格算法的优点,因而具有广阔的发展前景.     

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

Displacement-Based Experimental Stress Analysis of a Circularly-Perforated Asymmetrical Isotropic Structure

This paper demonstrates the ability to determine the individual full-field components of stress, strain and displacement in finite asymmetrical engineering structures from recorded values of a single displacement component while not necessitating differentiation of the measured displacement data; a process which can be ill-conditioned and adversely influenced by data noise and quality. While the present approach of processing recorded displacements with a stress function is not new, previous cases have been restricted to symmetrical scenarios. Full-field displacements are measured here using the digital image correlation (DIC) method. As well as satisfying load equilibrium, the experimental results agree with those predicted using the finite element method (FEM).     

A NEW HYBRID QUADRILATERAL FINITE ELEMENT FOR MINDLIN PLATE

ＡＮＥＷＨＹＢＲＩＤＱＵＡＤＲＩＬＡＴＥＲＡＬＦＩＮＩＴＥＥＬＥＭＥＮＴＦＯＲＭＩＮＤＬＩＮＰＬＡＴＥＣｈｉｎＹｉ（秦奕）（ＴｉａｎｊｉｎＡｒｃｈｉｔｅｃｔｕｒａｌＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅ，Ｔｉａｎｊｉｎ）ＺｈａｎｇＪｉｎｇ－ｙｕ（张敬宇）（Ｉ...     

求解接触问题的一种新的实验误差法

提出了一种带松弛因子的UZAW算法求解实验误差法中给定状态下的位移和接触力满足的等式方程，并证明了该算法是R超线性收敛的。整个区域被划分为多个子区域，不同子区域位移场的求解是独立的。还提出了一种带参数的以不完全因子分解为基础的预条件子共轭梯度法求解不同子区域位移场，该算法在块体规模较大时更加有效。     

一种提高四边形四节点平面壳单元计算精度的新方法

平面壳单元是由平面应力单元和平板弯曲单元叠加组合而成,具有简单的理论表达,但是它在计算曲面壳体结构时误差较大。为了进一步提高平面壳单元的计算精度,本文提出了一种计算平面壳单元刚度矩阵的新方法。通过该方法在高斯积分点建立多个单元局部坐标系,并保证每个局部坐标系都位于单元在高斯点处的切平面上,从而可以有效适应曲面壳体形状,达到进一步提高平面壳单元计算精度的目的。为了在这种新坐标系下计算单元刚度矩阵,给出了求解形函数对局部坐标的导数、局部到自然坐标系积分转换的雅可比、以及局部到整体坐标系的转换矩阵的新型计算方法。通过将这些新坐标系以及新计算方法运用到平面壳单元DKQ24中,可以有效提高平面壳单元尤其是在计算曲面壳体时的精度。计算结果表明,本文方法和平面壳单元相结合,不仅具有平面壳单元简单的理论表达式,还能得到满意的精度。另外,本文方法还可以应用到其他类型的平面壳单元,为提高其他类型平面壳单元的计算精度提供了一种新的途径和思路,具有广阔的应用前景。     

振荡应力奇异性及其强度系数的数值分析方法

本文以具有振荡应力奇异性的平面问题为例,提出了一种利用普通的数值分析结果（由有限地或边界元程序计算得到的应力分量或位移分量）,来确定奇异性次数以及相应的复应力强度系数的数值分析方法。为了验证该方法的有效性,应用平面应变情况下的边界元计算结果,对界面端模型进行了分析。计算结果表明,本方法可以精确地求得振动应力奇异性次数,并且与奇异性对应的复应力强度系数也可以很方便地应用外插法得到。     

弹性力学的杂交自然单元法

为了解决自然单元法无法直接求解节点应力以及应力解精度不高的问题, 将应力杂交的思想引入自然单元法中, 与弹性问题的Hellinger-Reissner变分原理结合, 提出了弹性问题的杂交自然单元法. 算例表明: 杂交自然单元法的计算结果与解析解吻合, 证明该方法是可行的; 在求解应力方面, 杂交自然单元法比自然单元法有着更高的计算精度, 而且可以直接求解出节点的应力.     

局部有限元模型边界条件的直接传递方法

提出一种将整体分析得到的节点力或节点位移直接传递到精细化局部有限元模型的方法,即部分混合单元法。沿精细化局部有限元模型周边建立一组过渡单元,该组过渡单元采用与整体模型一致的单元类型和模拟方式,其外侧边界上的节点与整体模型节点的相对坐标对应,内侧边界与精细化局部有限元模型采用基于面约束的方式连接。在外侧边界上根据节点坐标对应施加整体分析获得的节点力或节点位移,过渡单元就可直接将边界条件传递到精细化局部有限元模型。通过贵州红水河特大桥钢-混结合段的精细化有限元分析,验证了本文方法的实用性和有效性。