薄板弯曲问题的P型杂交解析有限元方法

本文采用两套变量构造有限元试函数空间，在单元内部要求试函数精确满足平衡微分方程，在单元边界上对位移和转角分别用Ｐｅａｎｏ升阶函数插值，然后利用广义变分原理建立了一种薄板弯曲问题的Ｐ型杂交解析有限方法，与常规有限元法相比，该方法不心进行过细的网格剖分，通过增加单元插值多项式的阶数Ｐ来提高精度，此外，该方法还具有积分计算只需在单元边界上进行、单元钢度矩阵和载荷向量具有嵌入结构、协调程度可以自动控制等优     

一维区间B样条小波单元的构造研究

基于区间B样条小波及小波有限元理论,提出了一种区间B样条小波有限元方法。传统有限元多项式插值被一维区间B样条小波尺度函数取代,进而构造形状函数和单元。与小波Galer-kin方法不同,本文构造的区间B样条小波单元通过转换矩阵将无明确物理意义的小波插值系数转换到物理空间。转换矩阵在小波单元构造过程中起到关键作用,为了保证求解的稳定性,转换矩阵必须非奇异。构造了以区间B样条尺度函数为插值函数的一系列一维区间B样条小波单元。数值算例表明,本文构造的区间B样条小波单元与传统有限元方法相比,在求解变截面,变载荷等问题时具有收敛快和精度高等优势;有效地丰富了小波有限元法单元库。     

固体中短波传播单位分解有限元法的解析积分

针对固体中短波传播数值模拟的单位分解有限元法中单元矩阵积分的被积函数的强烈振荡特性,应用直角坐标系下标准有限元形函数和单元内的波动方向知识提出了一种单元矩阵的解析积分方案。它对于平面三,六,四,八和九节点的直边单位分解有限单元是完全解析的,对于与这些单元相应的曲边单元则是半解析的。数值结果显示所提出的积分方案在计算效率上比高斯-勒让德积分有大幅度提高。     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

BUCKLING ANALYSIS UNDER COMBINED LOADING OF THIN-WALLED PLATE ASSEMBLIES USING BUBBLE FUNCTIONS

Bubble functions are finite element modes that are zero on the boundary of the element but nonzero at the other point. The present paper adds bubble functions to the ordinary Complex Finite Strip Method(CFSM) to calculate the elastic local buckling stress of plates and plate assemblies. The results indicate that the use of bubble functions greatly improves the convergence of the Finite Strip Method(FSM) in terms of strip subdivision, and leads to much smaller storage required for the structure stiffness and stability matrices. Numerical examples are given, including plates and plate structures subjected to a combination of longitudinal and transverse compression, bending and shear. This study illustrates the power of bubble functions in solving stability problems of plates and plate structures. The work was partly supported by the Natural Science Foundation of Jiangxi Province of China and by the Basic Theory Research Foundation of Nanchang University.     

Application of the quadrilateral area coordinate method： a new element for laminated composite plate bending problems

Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate (QAC) method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF (5-DOF per node) quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-processing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models. The project is supported by the National Natural Science Foundation of China (10502028), the Special Foundation for the Authors of the Nationwide (China) Excellent Doctoral Dissertation (200242), and the Science Research Foundation of China Agricultural University (2004016).     

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.      

Shear stress analysis in engineering beams using deplanation field of special 2-D finite elements

This paper deals with the 2-D finite element shear stress analysis in beams, loaded by bending with shear and St. Venant’s torsion. The properties of these finite elements, like stiffness matrices as well as load vectors, are derived on the basis of their axial nodal displacements, e.g. by warping field. Proposed finite elements enable stress analysis independently of both cross-sectional member shape and material properties. Stiffness matrices and load vectors are derived for several finite element types. Material is assumed to be isotropic and linear elastic. For justification of the proposed stress analysis procedure, some examples are presented.     

多边形有限元研究进展

有限元法是数值求解偏微分方程边值问题的重要方法,采用不规则多边形单元网格, 可以方便有效地模拟材料的力学性能, 又使得区域网格剖分变得灵活方便. 特别是对于复杂的几何形状, 多边形单元网格具有更大的优势. 本文对国内外有关多边形有限元法的最新进展作了初步的总结和评述, 主要以基于位移法的多边形有限元为主.论述了多边形有限元的发展历史, 给出了多边形单元上的Wachspress插值、Laplace插值和重心坐标的一些最新研究成果. 与经典有限元法形函数为多项式形式不同, 多边形单元的形函数为有理函数或者无理函数形式. 多边形单元插值形函数满足线性完备性, 可以再现线性位移场, 像经典有限元法一样直接施加本质边界条件; 插值函数在多边形的边界上是线性的,确保不同单元间的自动协调. 不同单元的插值形函数表达公式形式统一, 方便混合单元网格计算的程序编写. 提出了多边形有限元法今后需要研究的问题.      

Acoustic simulation using a novel approach for reducing dispersion error

It is well‐known that the traditional finite element method (FEM) fails to provide accurate results to the Helmholtz equation with the increase of wave number because of the ‘pollution error’ caused by numerical dispersion. In order to overcome this deficiency, a gradient‐weighted finite element method (GW‐FEM) that combines Shepard interpolation and linear shape functions is proposed in this work. Three‐node triangular and four‐node tetrahedral elements that can be generated automatically are first used to discretize the problem domain in 2D and 3D spaces, respectively. For each independent element, a compacted support domain is then formed based on the element itself and its adjacent elements sharing common edges (or faces). With the aid of Shepard interpolation, a weighted acoustic gradient field is then formulated, which will be further used to construct the discretized system equations through the generalized Galerkin weak form. Numerical examples demonstrate that the present algorithm can significantly reduces the dispersion error in computational acoustics. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of piezoelectric ceramic multilayer actuators based on an electro-mechanical coupled meshless method

This paper presents an efficient meshless method for analyzing cracked piezoelectric structures subjected to mechanical and electrical loading. In this method, an element free Galerkin (EFG) formulation, an enriched basic function and some special shape functions that contain discontinuous derivatives are employed. Based on the moving least squares (MLS) interpolation approach, the EFG method is one of the promising methods for dealing with problems involving progressive crack growth. Since the method is meshless and no element connectivity data are needed, the burdensome remeshing procedure required in the conventional finite element method (FEM) is avoided. The numerical results show that the proposed method can yield an accurate near-tip stress field in an infinite piezoelectric plate containing an interior hole. In another example studying a ceramic multilayer actuator, the proposed model was found to be accurate in the simulation of stress and electric field concentrations arround the abrupt end of an internal electrode. The project supported by the National Natural Science Foundation of China (10025209, 10132010, and 90208002), and the Research Grants Council of the Hong Kong Special Administrative Region, China (HKU 7203/03E). The English text was polished by Yunming Chen.     

A mixed finite element scheme for viscoelastic flows with XPP model

A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.     

A meshless method of a thin plate

A meshless approach to analysis of arbitrary Kirchhoff plates by the local boundary integral equation(LBIE) method is presented. The method combines the advantageous features of, all the three methods: the Galerkin finite element method (GFEM), the boundary element method (BEM) and the element-free Galerkin method (EFGM). It is a truly meshless method, which means that the discretization is independent of geometric subdivision into elements or cells, but is only based on a set of nodes (ordered or scattered) over a domain in question. It involves only boundary integration, however, over a local boundary centered at the node in question; It poses no difficulties in satisfying the essential boundary conditions while leading to banded and sparse system matrices using the moving least square (MLS) approximations. It is shown that high accuracy can be achieved for arbitrary geometries for clamped and simply-supported edge conditions. The method is found to be simple, efficient, and attractive. Project supported by the National Science Foundation of China (No. 19972019).     

Vertex solution theorem for the upper and lower bounds on the dynamic response of structures with uncertain-but-bounded parameters

The aim of this paper is to evaluate the effects of uncertain-but-bounded parameters on the dynamic response of structures. By combining the interval mathematics and the finite element analysis, the mass matrix, damping matrix, stiffness matrix and the external loads are represented as interval matrices and vector. With the help of the optimization theory, we present the vertex solution theorem for determining both the exact upper bounds or maximum values and the exact lower bounds or minimum values of the dynamic response of structures, in which these parameters reach their extreme values on the boundary of the interval mass, damping, stiffness matrices and the interval extemal loads vector. Three examples are used to illustrate the computational aspects of the presented vertex solution theorem.     

A triangular plate element 2343 using second-order absolute-nodal-coordinate slopes: numerical computation of shape functions

In this paper, the process by which geometrical and structural matrices of plate finite elements employing absolute nodal coordinate formulation (ANCF) are constructed is studied. The kinematic and topological properties of an arbitrary plate finite element are described using universal digital code dncm that provides systematic enumeration of finite elements. This code is formed using the element’s dimension d, the number of nodes it possesses n, the number of scalar coordinates per node c, and a multiplier describing the process of transforming a conventional finite element to an ANCF element m. The detailed generation of a new type of triangular plate finite element 2343 using numerical computation of shape functions is also discussed in the paper. The new triangular element employs position vectors and slope vectors up to second-order mixed-derivative slope vector. A detailed derivation of the equations of motion of the element is also provided and examples of its numerical simulation and validation presented.     

应用分形有限元方法于外域声场计算

 应用二级分形有限元方法计算了外域声场. 用一人工边界把外域声场分为两 部分，人工边界以内使用常规有限元方法，人工边界以外的无限大区域使用分形有限元方法. 使用分形有限元方法的优点是：一方面形成几何自相似网格使得相邻层之间的单元刚度矩阵 和质量矩阵具有非常简单的关系；另一方面引用自动满足无限远辐射条件的全域插值函数把 节点自由度变换为一组广义坐标，因而计算量可以大大减少. 数值算例表明：该方法对于计 算无限大外域声场是有效的.     

非网格重剖分模拟宏观裂纹体的扩展有限单元法(Ⅰ:基础理论)

基于有限元计算网格,扩展有限单元法通过建立特殊的广义节点插值形式来描述含裂缝体的不连续位移场,避免了有限元法模拟裂缝时需要的网格重划分。进而,本文从虚功原理出发,在有限元法框架内完整地推导了能模拟宏观裂纹力学场的扩展有限元法实现公式,在理论上更全面地考虑了内部裂纹面上分布外载荷及缝内粘连材料刚度的影响,并提出了构建统一的扩展有限单元刚度阵形成模式,保证了与传统有限单元方式的协调一致。文中对方法的实现过程也做了详细阐述,给出了通用的计算公式,确保了算法的可行性。     

Transient thermal stresses in an orthotropic hollow cylinder for axisymmetric problems

For the thermoelastic dynamic axisymmetric problem of a finite orthotropic hollow cylinder, one comes closer to reality to involve the effect of axial strain than to consider the plane strain case only. However, additional mathematical difficulties should be encountered and a different solution procedure should be developed. By the separation of variables, the thermoelastic axisymmetric dynamic problem of an orthotropic hollow cylinder taking account of the axial strain is transformed to a Volterra integral equation of the second kind for a function of time, which can be solved efficiently and quickly by the interpolation method. The solutions of displacements and stresses are obtained. It is noted that the present method is suitable for an orthotropic hollow cylinder with an arbitrary thickness subjected to arbitrary axisymmetric thermal loads. Numerical comparison is made to show the effect of the axial strain on the displacements and stresses. The project supported by the National Natural Science Foundation of China (10172075) and China Postdoctoral Science Foundation (20040350712)     

A variational solution of 2-D sound-structure interaction problems

Based on the fluid-structure coupling theory, a localized variational principle for analyzing the sound radiation from elastic structure submerged in water due to harmonic excitations is presented. It will be a powerful tool to formulate various numerical methods for steady response of structural-acoustic systems. By means of this variational principle a hybrid element method, in which an analytical solution valid in most of the surrounding water is incorporated with finite elements distributed in the structure and its neighboring water, is devised. Computational examples are then given to demonstrate its high accuracy and time saving. The project supported by the National Natural Science Foundation of China (10172038) and the Doctoral Program Foundation of Institution of Higher Education of China (20040487013). The English text was polished by Yunming Chen.     

Stochastic boundary element method in elasticity

The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and randomly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displacements are derived, respectively. It is found that the randomness of material parameters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods. The project supported by the National Natural Science Foundation of China and the State Education Commission Foundation of China