APPLICATION OF PENALTY FUNCTION METHOD IN ISOPARAMETRIC HYBRID FINITE ELEMENT ANALYSIS

By the aid of the penalty function method, the equilibrium restriction conditions were introduced to the isoparametric hybrid finite element analysis, and the concrete application course of the penalty function method in three-dimensional isoparametric hybrid finite element was discussed.The separated penalty parameters method and the optimal hybrid element model with penalty balance were also presented. The penalty balance method can effectively refrain the parasitical stress on the premise of no additional degrees of freedom. The numeric experiment shows that the presented element not only is effective in improving greatly the numeric calculation precision of distorted grids but also has the universality.     

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

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

Large deformation of incompressible rubber cylinder under plane strain

The large deformation of incompressible rubber cylinder under inner pressure is analyzed by a kind of new rubber materials strain energy function. The theory formulation for the displacement and stress is presented. The penalty finite element formulation is established in order to analyze nonlinear rubber materials, and the results of finite element method agree well with theoretical ones. A new method for controlling the calculating stability and convergence rates is put forward. The selection of the appropriate penalty factor and its influence on calculated results are discussed.     

Stream function finite element solution of Stokes flow with penalties

A finite element method for solution of the stream function formulation of Stokes flow is developed. The method involves complete cubic non-conforming (C⁰) triangular Hermite elements. This element fails the patch test. To correct the element and produce a convergent method we employ a penalty method to weakly enforce the desired continuity constraint on the normal derivative across the inter-element boundaries. Successful use of the method is demonstrated to require reduced integration of the inter-element penalty with a 1-point Gauss rule. Error estimates relate the optimal choice of penalty parameter to mesh size and are corroborated by numerical convergence studies. The need for reduced integration is interpreted using rank relations for an associated hybrid method.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

A new stable space–time formulation for two‐dimensional and three‐dimensional incompressible viscous flow

A space–time finite element method for the incompressible Navier–Stokes equations in a bounded domain in ?^d (with d=2 or 3) is presented. The method is based on the time‐discontinuous Galerkin method with the use of simplex‐type meshes together with the requirement that the space–time finite element discretization for the velocity and the pressure satisfy the inf–sup stability condition of Brezzi and Babu?ka. The finite element discretization for the pressure consists of piecewise linear functions, while piecewise linear functions enriched with a bubble function are used for the velocity. The stability proof and numerical results for some two‐dimensional problems are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Finite element analysis of creeping flows using marker particles

The penalty function formulation of the finite element method is described for the analysis of transient incompressible creeping flows. Marker particles are utilized to represent moving free surfaces and to visualize the flow patterns. For determining the movement of markers from element to element, the area coordinate system of the linear triangular element is introduced. With the method presented, a punch indentation problem and an injection problem for an L-shaped cavity are solved for Newtonian and power-law fluids.     

不连续温度场问题的间断Galerkin方法

针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允许插值函数在单元边界上存在跳变.在建立有限元方程时,通过在单元边界上引入数值通量项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出.数值算例表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高.该方法也为研究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段.     

构造满足特征值要求的杂交元应力子空间

通过引进适当参数给出了杂交元应力空间特征值与参数之间的关系,从而可以通过调节参数来构造满足特征值即模态刚度要求的杂交元应力子空间。本在位移元本征应力模式基础上引进调节参数,同时,利用矩阵H对角化方法计算杂交元应力子空间的本征应力模式,然后由此方便有效地计算特征值,从而大大提高了计算效率。本通过建立Q4杂交应力元特征值与参数之间关系说明了这一方法是确定可行的。     

Derivation of geometric stiffness and mass matrices for finite element hybrid models

An unifying approach in deriving the geometric stiffness and mass matrices for finite element hybrid models is presented. The variational formulation is based on a modified Reissner Principle. Numerical verification is illustrated through a simple beam example. The element interchangeability in the finite element method and an alternative way of deriving the loading vector for the hybrid stress model are also discussed.     

平面应变不可压缩橡胶圆柱的大变形

针对一类新的橡胶材料应变能函数，推导了受内压作用下不可压缩橡胶圆柱的大变形公式，给出了位移、应力的解析表达式（用积分形式表示）。建立了适用于分析非线性不可压缩橡胶材料的罚有限元列式，算例表明：位移与应力能很好地与理论解吻合，并提出了控制计算稳定的方法，特别是详细地讨论罚因子的选取及对计算结果的影响。     

Variable penalty method for finite element analysis of incompressible flow

A new scheme is applied for increasing the accuracy of the penalty finite element method for incompressible flow by systematically varying from element to element the sign and magnitude of the penalty parameter λ, which enters through ?.v + p/λ = 0, an approximation to the incompressibility constraint. Not only is the error in this approximation reduced beyond that achievable with a constant λ, but also digital truncation error is lowered when it is aggravated by large variations in element size, a critical problem when the discretization must resolve thin boundary layers. The magnitude of the penalty parameter can be chosen smaller than when λ is constant, which also reduces digital truncation error; hence a shorter word-length computer is more likely to succeed. Error estimates of the method are reviewed. Boundary conditions which circumvent the hazards of aphysical pressure modes are catalogued for the finite element basis set chosen here. In order to compare performance, the variable penalty method is pitted against the conventional penalty method with constant λ in several Stokes flow case studies.     

Thermal solution of cylindrical composite systems using meshless method

This paper presents the thermal solution of cylindrical composite systems using meshless element free Galerkin (EFG) method. The EFG method utilizes the moving least square approximants, which are constructed by using a weight function, a basis function and a set of non-constant coefficients to approximate the unknown function of temperature. Dirichlet (essential) boundary conditions have been enforced using Lagrange multiplier and penalty methods. Existing rational weight function has been modified and used in the present analysis. MATLAB codes have been developed to obtain the numerical solution. The EFG results have been obtained using cubicspline, quarticspline, Gaussian, quadratic, hyperbolic, exponential, rational and cosine weight functions for a model problem. The results obtained using different EFG weight functions are also compared with those obtained by finite element method. The effect of scaling and penalty parameters has also been studied in detail.     

A 3D shell-like approach using element-free Galerkin method for analysis of thin and thick plate structures

A new efficient meshless method based on the element-free Galerkin method is proposed to analyze the static deformation of thin and thick plate structures in this paper. Using the new 3D shell-like kinematics in analogy to the solid-shell concept of the finite element method, discretization is carried out by the nodes located on the upper and lower surfaces of the structures. The approximation of all unknown field variables is carried out by using the moving least squares (MLS) approximation scheme in the in-plane directions, while the linear interpolation is applied through the thickness direction. Thus, different boundary conditions are defined only using displacements and penalty method is used to enforce the essential boundary conditions. The constrained Galerkin weak form, which incorporates only displacement degrees of freedom (d.o.f.s), is derived. A modified 3D constitutive relationship is adopted in order to avoid or eliminate some self-locking effects. The numeric efficiency of the proposed meshless formulation is illustrated by the numeric examples.     

无网格RKPM法模拟平板轧制过程

将无网格再生核质点法(RKPM)用于刚塑性可压缩材料轧制过程的模拟,采用罚函数满足本质边界条件,选用矩形影响域的张量积核函数,利用有限元网格作为积分单元,对求解域内和边界上采用不同的高斯积分方案。并将计算结果与刚塑性有限元计算结果和文献中的实验数据相比较,说明RKPM方法用于刚塑性可压缩材料轧制过程的可行性和正确性。     

Velocity–pressure integrated versus penalty finite element methods for high-Reynolds-number flows

Velocity–pressure integrated and consistent penalty finite element computations of high-Reynolds-number laminar flows are presented. In both methods the pressure has been interpolated using linear shape functions for a triangular element which is contained inside the biquadratic flow element. It has been shown previously that the pressure interpolation method, when used in conjunction with the velocity-pressure integrated method, yields accurate computational results for high-Reynolds-number flows. It is shown in this paper that use of the same pressure interpolation method in the consistent penalty finite element method yields computational results which are comparable to those of the velocity–pressure integrated method for both the velocity and the pressure fields. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine grid finite difference computational results. Advantages and disadvantages of the two finite element methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow of Reynolds number 10 000, a laminar backward-facing step flow and a laminar flow through a nest of cylinders.     

Penalty finite element method for nonlinear dynamic response of viscous fluid-saturated biphase porous media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｐｏｒｏｕｓｍｅｄｉａｍｏｄｅｌｓｄｅｄｕｃｅｄｆｒｏｍｍｉｘｔｕｒｅｔｈｅｏｒｙａｒｅａｔｔｒａｃｔｉｎｇａｔｔｅｎｔｉｏｎｏｆｍｏｒｅａｎｄｍｏｒｅｒｅｓｅａｒｃｈｅｒｓｂｅｃａｕｓｅｉｔｉｓｂａｓｅｄｏｎｔｈｅｆｒａｍｅｏｆｃｏｎｔｉｎｕｕｍｍｅｃｈａｎｉｃｓ[1]ａｎｄｗｅｒｅｓｕｃｃｅｓｓｆｕｌｌｙｕｓｅｄｔｏｄｅｐｉｃｔｔｈｅｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒｓｏｆｓｏｉｌｓａｎｄｂｉｏｌｏｇｉｃａｌｓｏｆｔｔｉｓｓｕｅｓ[2 ,3].Ｗｉｔｈｔｈｉｓｔｙｐｅｏ…     

discontinuous galerkin method for discontinuous temperature field problems

针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允 许插值函数在单元边界上存在跳变. 在建立有限元方程时,通过在单元边界上引入数值通量 项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出. 数值算例 表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局 部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高. 该方法也为研 究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段.     

A review of some finite element methods to solve the stationary Navier-Stokes equations

In this paper the integrated solution approach, the penalty function approach and the solenoidal approach for the finite element solution of the stationary Navier-Stokes equations are compared. It is shown that both the penalty function approach and the solenoidal approach compare favourably to the integrated solution method. For fine meshes the solenoidal approach appears to be the cheapest method.     

橡胶楔体与刚性缺口接触的大变形分析

利用一种新的橡胶材料应变能函数，对橡胶楔体与刚性缺口接触大变形问题进行了分析。得到了接触尖点附近变形的奇异性特征，给出了奇异性指数与材料常数、橡胶楔体角度、刚性缺口角度之间的关系式。同时编制了大变形有限元程序，计算得到了与理论解一致的结论。