满足PTC的三维非协调形函数的显式处理

介绍了满足分片检验条件（ＰＴＣ）的非协调形函数生成方法，显式给出了满足ＰＴＣ的三维非协调元ＮＨ１１的内位移函数。计算表明显式处理后的单元与原单元具有相同的精度，但却在很大程度上节省了单刚形成时间和计算单元应力花费，进一步发挥了非协调元的数值潜力。     

横观各向同性轴对称问题的非协调元和杂交应力元计算

基于变分原理得出各向同性轴对称问题下的非协调元和杂交应力元方法仍然适用于分析横观各向同性轴对称问题的结论，同时对应用于各向同性问题的罚平衡优化方法进行了修改，使之能够应用于横观各向同性问题的分析。文中给出了分析算例。并对各种单元结果进行了比较，计算结果表明非协调元和杂交应力元方法不但适用于横观各向同性轴对称问题分析，而且将提高其数值解的精度，改善单元内部应力分布。     

关于轴对称问题的应力函数

本文论述按应力求解轴对称问题的协调方程和应力函数,建立了应力函数与——Neuber通解间的关系。 1.协调方程。轴对称问题的平衡方程是     

带转动自由度的内参型非协调元研究

非协调元虽然破坏了单元间位移的连续性，却能很好地反映弯曲类变形，然而在不增加单元结点自由度的情况下，非协调元的计算精度总是滞留在某一水平，无法得到较大改变。基于修正后的位移型Reissner泛函中引入独立转动场的变分原理，采用连续介质力学中的转动自由度的定义，转动场采用结点真实转角来插值，结合平面四结点单元讨论了有效附加非协调位移的合理形式，引入了适用于任何四边形单元的非协调位移函数，从而建立了一种带转动自由度的平面四结点内参型非协调元模型。本文单元能通过分片检验，并易于与带转动自由度的梁单元相容．教值算例表明具有较高的计算精度。     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法,在此基础上建立起了构造等参非协调元的新方法。作为实例,构造出两个可以给出单刚显式的四结点平面非协调新单元。     

轴对称弹性应变梯度理论公式推导及有限元实现

用张量运算推导了弹性应变梯度轴对称问题的基本公式.建立了应变梯度轴对称不协调元的弱连续条件,进一步建立了满足弱连续条件的应变梯度轴对称18-DOF三角形单元(BCIZ+ART9),其中BCIZ满足线性应变C0连续,用于计算应变ε;ART9满足常曲率C1弱连续,用于计算应变梯度η0数值结果表明该单元通过C0-1分片检验并能体现材料的尺度效应.     

含球形空腔或刚性夹杂的中厚圆板在弯曲变形时的弹性场

本文将空间轴对称问题的Папковиц-Neuber通解用复变量广义解析函数表示,推导出用复变函数法求解空间轴对称问题的基本公式,并以此为工具求得了含球形空腔或刚性夹杂的中厚圆板在轴对称弯曲变形时的完全解.     

最佳协调性四边形薄板元

本文利用最佳逼近思想，构造了一种任意四边形广义协调薄板元。文中提出了协调性度量泛函的概念。通过协调性度量泛函的极小化，使单元协调性达到最佳。该方法的特点是非常便于程式化，容易推广到其它类型非协调元的构造，数值检验表明，该类单元精度很好，且对不规则网格有较好的适应性，此外，在使用协调性度量泛函的前提下，令插值基函数在笛卡尔坐标中直接取用，可使小片检验自然满足。     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法，在此基础上建立起了构造等参非协调元的新方法，作为实例，构造出两个可以给出单刚显式的四结点平面非协调新单元。     

关于非协调位移元与广义杂交元的等价性的进一步研究

本文通过定义非协调位移元的单元特性函数,提出了一种建立非协调元与广义杂交元之间的等价关系的一般方法,进一步证明了两种有限元模型之间存在强等价性,为两种单元模型的构造提供了相互借鉴的理论依据和有效的准则。     

A class of Wilson arbitrary quadrilateral elements for an axisymmetric problem

A class of modified Wilson arbitrary quadrilateral nonconforming elements for an axisymmetric problem is proposed. Their convergence is proven by means of the strong patch test. The structure of this finite element class is investigated. Thus a general method of axisymmetric nonconforming elements with convergence properties is presented. Project supported by the Foundation of the Chinese Ministry of Machine-Building Industry     

Finite element analysis of axisymmetric elastic body problems

Linear form functions are commonly used in a long time for a toroidal volume element swept by a triangle revolved about the symmetrical axis for general axisymmetrical stress problems. It is difficult to obtain the rigidity matrix by exact integration, and instead, the method of approximate integration is used. As the locations of element close to the symmetrical axis, the accuracy of this approximation deteriorates very rapidly. The exact integration have been suggested by various authors for the calculation of rigidity matrix. However, it is shown in this paper that these exact integrations can only be used for those axisymmetric bodies with central hole. For solid axisymmetric body, it can be proved that the calculation fails due to the divergent property of rigidity matrix integration. In this paper a new form function is suggested. In this new form function, the radial displacementu vanishes as radial coordinatesr approach to zero. The calculated rigidity matrix is convergent everywhere, including these triangular toroidal element closed to the symmetrical axis. This kind of element is useful for the calculation of axisymmetric elastic solid body problems.     

THE CONSTRUCTION OF WAVELET-BASED TRUNCATED CONICAL SHELL ELEMENT USING B-SPLINE WAVELET ON THE INTERVAL

Based on B-spline wavelet on the interval (BSWI), two classes of truncated conicalshell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conicalshell element and BSWI moderately thick truncated conical shell element with independent slope-deformation interpolation. In the construction of wavelet-based element, instead of traditionalpolynomial interpolation, the scaling functions of BSWI were employed to form the shape functionsthrough the constructed elemental transformation matrix,and then construct BSWI element viathe variational principle. Unlike the process of direct wavelets adding in the wavelet Galerkinmethod, the elemental displacement field represented by the coefficients of wavelets expansionwas transformed into edges and internal modes via the constructed transformation matrix. BSWIelement combines the accuracy of B-spline function approximation and various wavelet-basedelements for structural analysis. Some static and dynamic numerical examples of conical shellswere studied to demonstrate the present element with higher efficiency and precision than thetraditional element.     

Notes on the finite element analysis of the axisymmetric elastic solid

A simple transformation of displacements considerably eases the explicit derivation of the finite element stiffness matrix for the axisymmetric elastic solid without causing a decline in the rate of convergence. The worsening of the condition of the global stiffness matrix caused by this transformation can be cured by scaling. A balanced numerical integration scheme maintaining the full rate of convergence is the one that integrates each term of the work and energy expressions to the order 2p ? 2, p being the degree of the complete polynomial in the shape functions.     

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method. The English text was polished byYunming Chen.     

轴对称热传导有限元格式

分析了目前一些有限元专著中轴对称热传导有限元方法推导中的问题,给出了轴对称热传导有限元格式的正确表达形式。     

Stress function method in shakedown analysis of axisymmetric shells

A self-equilibrated stress obtained from the stress functions of thin shells is used for the static shakedown theorem as a residual stress. In combination with the finite element method, a linear programming formulation of the shakedown analysis of axisymmetric shells is derived. The physical meaning of the stress function method is clear and its computing amount is small. Some examples of the plates and shells show that the method is reasonable and efficient.The project was supported by the National Natural Science Foundation of China     

液体自由晃动的数值分析

本文对常重条件下自旋容器内液体的对称和任意阶反对称晃动建立了边界元算法，引入了两种调和函数以确定边界元系数矩阵的主对角元，并且提出了一种对计算机程序的正确性和数值解近似程度的考核方法。     

Axisymmetric viscoplastic deformation by the boundary element method

This paper presents a boundary element formulation and numerical implementation of the problem of small axisymmetric deformation of viscoplastic bodies. While the extension from planar to axisymmetric problems can be carried out fairly simply for the finite element method (FEM), this is far from true for the boundary element method (BEM). The primary reason for this fact is that the axisymmetric kernels in the integral equations of the BEM contain elliptic functions which cannot be integrated analytically even over boundary elements and internal cells of simple shape. Thus, special methods have to be developed for the efficient and accurate numerical integration of these singular and sensitive kernels over discrete elements. The accurate determination of stress rates by differentiation of the displacement rates presents another formidable challenge.A successful numerical implementation of the boundary element method with elementwise (called the Mixed approach) or pointwise (called the pure BEM or BEM approach) determination of stress rates has been carried out. A computer program has been developed for the solution of general axisymmetric viscoplasticity problems. Comparisons of numerical results from the BEM and FEM, for several illustrative problems, are presented and discussed in the paper. It is possible to get direct solutions for the simpler class of problems for cylinders of uniform cross-section, and these solutions are also compared with the BEM and FEM results for such cases.     

基于自然单元法的轴对称结构极限上限分析

自然单元法是一种以自然邻近插值为试函数的新兴无网格数值方法,其形函数的计算不涉及矩阵求逆,也不需要任何人为参数。为了充分发挥自然单元法的优势,本文基于极限分析上限定理建立了轴对称结构极限上限分析的整套求解算法。轴对称结构的位移场由自然邻近插值构造,并且采用罚函数法处理材料的不可压条件。为了消除目标函数非光滑所引起的数值困难,采用逐步识别刚性区和塑性区,并对两者用不同方法进行处理。数值算例结果表明,本文提出的轴对称结构极限上限分析方法是行之有效的。