二维边界元奇异积分和多域缩聚法分析

基于基本解的一种新的表达式，对二维边界元分析中奇异积分的精确求解进行了讨论，从几何方面对基本解的奇异性进行了分析，给出了超参非连续元离散位势和弹性力学问题边界积分方程时奇异积分计算的精确式，从而为判断各种近似方法的优劣和间接方法的精度提供了依据，也为精确地分析了大规模问题提供了一条有效的途径。     

Navier-Stokes方程间断Galerkin有限元方法研究

通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤.同时还给出了一种简单有效的计算方法来对物面边界进行高阶近似.为了能够对NS方程进行精度测试,采用对原始系统添加源项的方法构造精确解.二维Euler和NS系统的精度测试表明该方法达到了DG方法的理论精度.二维圆柱无黏绕流的计算结果表明关于物面边界的高阶近似方法能够保持DG方法原有的精度.卡门涡街数值模拟则进一步验证了该方法的正确性并且显示出DG方法较高的计算精度和分辨率.     

提高多体系统离散时间传递矩阵法计算精度的研究

深入探讨多体系统离散时间传递矩阵法对平面、空间刚体-光滑铰多体系统运动响应的研究。提出提高该方法计算精度和计算稳定性的方法，导出相应的多端刚体传递矩阵。设计了角坐标的迭代循环，该变量不必采用近似形式，因而提高了计算精度和计算稳定性。其它方法有：1)增加泰勒展开式的高阶项；2)合理选择将速度、加速度表示为位移的线性函数的方法(本文简称线性化方法)；3)合理确定迭代循环的初值等。对4种多体系统进行了计算机仿真研究，表明本文提出的方法有效。     

平面叶栅跨音流场气动特性的数值研究

通过对模型方程的分析，给出了一种新的隐格式构造思想。将它运用到关通量分裂格式中，可得到无近似因子分解、无矩阵运算的高效二阶精度隐式矢通量分裂差分格式，并用来直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组。数值计算标明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

近似平衡多项式加速度动力显式算法

利用已知初始时刻的信息，建立一种可以取到任意阶高精度的多项式加速度单步隐式算法。在该隐式方法中，待采解方程纽系数矩阵中质量阵的系数远远大于阻尼阵和剐度阵的系数，略去非对角阻尼阵和非对角刚度阵对方程组的影响，得到一种近似平衡多项式加速度动力显式计算方法。此方法的精度主要由加速度多项式插值的项数、步长、质量阵的每件数、质量刚度比（质量阵和刚度阵的范数之比）决定。在此基础上给出了这种算法的通式，进行了精度分析，结果表明：如果时间步长h足够短，n次加速度近似平衡动力显式算法的精度可以达到O（hn＋1）。算例采用5次加速度近似平衡显式算法，计算结果的精确性证明了本算法的可行性。     

弹性力学平面问题中一类无奇异边界积分方程

从理论上提出一种新的方法，归化出间接变量无奇异边界积分方 程. 采用Lagrange二次单元，建立一个数值求解框架系统. 此外，基于问题的计算区域的 特殊性，给出一种边界近似方法. 数值算例表明该方法所取得的数值结果与精确解相当接 近，特别是边界量的数值结果. 此外，该方法容易被推广到三维问题. 和已有的直接变量的情形相比较，具有优点：1)无需处理HFP积分. 大大降低处理问 题的复杂性，并提高了计算效率和解的精度；2)摆脱了问题的具体形式，进入纯代数操作. 这样做的好处是从理论上建立一种普遍适用的方法，不仅适用于弹性力学问题，同样可应用 于其它问题，如位势问题, Stokes问题等. 3)提供了一种计算CPV积分的方法.     

任意激励下结构动力响应的状态方程精细积分法

对只有弹性模态以及除此之外还有刚体模态的结构的瞬态响应给出了精细积分的通用公式，从而使得该方法不仅可以处理线性激励的情形，而且对激励是多项式形式或可以展开成多项式的激励也同样能够计算。对于非线性激励，只要可以用关于自变量的级数形式来近似表示，都可以用本文所给的方法进行计算，计算的精度可以通过变化级数的项数来调整。     

陷落腔的流场数值计算研究

本文对有限分析解法的近似计算方法及其对计算精度的影响进行了数值计算和分析讨论；然后采用有限分析解法，对前、后壁壁缘高度不同的三种腔内流场运动发展过程及腔口剪切层运动情况作了数值计算，得到的腔口剪切层上腔口随边相互作用情况与对应腔的流场显示实验结果进行了比较，两者一致，计算结果表明，采用本文发展的近似方法计算可靠，可以较大地提高计算效率。     

叶轮机械迴转面流场的有限元近似解方法的研究和数值计算

本文提出一种新的概念:即在应用有限元处理不同类型的流体力学、传热学问题时,应该根据各自问题的物理本质,采用不同型式的基函数(shape function),并在文[1]的基础上,应用有限元近似解方法(FEAM)对直角弯管流动进行了数值计算和精度分析,还对带长短叶片的前向离心通风机的迴转面流场进行了数值计算和分析,在迴转面流场计算中采用了一种新的简化方法,解决了在非矩形单元上求取近似解的困难,计算结果表明在稀疏网格下能保持较高精度,从而明显地提高了计算速度。     

关于几种新的析分解的计算方法

对变形樟度极分解的计算方法进行了分析，给出了极分解计算的四种新方法：（１）增量叠加法；（２）基于伸长张量不变量（近似）计算法；（３）确定主转动轴计算法；（４）坐标变换法。增加叠加极分解计算方法将为建立以伸长张量为变应度量的大变形大转动有限分析方法提供基础。本文还给出了伸长张量物质时间导数的简洁表达式。     

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.     

Exactness of penalization for exact minimax penalty function method in nonconvex programming

The exact minimax penalty function method is used to solve a noncon-vex differentiable optimization problem with both inequality and equality constraints. The conditions for exactness of the penalization for the exact minimax penalty function method are established by assuming that the functions constituting the considered con-strained optimization problem are invex with respect to the same function η (with the exception of those equality constraints for which the associated Lagrange multipliers are negative-these functions should be assumed to be incave with respect to η). Thus, a threshold of the penalty parameter is given such that, for all penalty parameters exceeding this threshold, equivalence holds between the set of optimal solutions in the considered constrained optimization problem and the set of minimizer in its associated penalized problem with an exact minimax penalty function. It is shown that coercivity is not suf-ficient to prove the results.     

无网格局部Petrov-Galerkin方法中本质边界条件的处理

鉴于无网格局部Petrov-Galerkin方法(MLPG)形函数的非插值性质,将一种新的本质边界处理方案——完全变换法与MLPG结合,通过变换矩阵修正形函数,使其满足Kronecker-δ条件,实现了本质边界的精确实施。进一步与MLPG中通常处理边界的罚方法作了比较研究,数值结果表明新方法的可靠性与精确度。     

几何精确NURBS有限元中边界条件施加方式对精度影响的三维计算分析

非均匀有理B样条(NURBS)有限元法把计算机辅助几何设计(CAGD)中的NURBS几何构形方法与有限元方法有机结合起来,有效消除了有限元离散模型的几何误差,提高了计算精度。但是由于NURBS基函数不是插值函数,直接在控制节点上施加位移边界条件会引起较大误差。本文详细讨论了NURBS基函数的插值特性,在NURBS有限元分析中采用罚函数法施加位移边界条件,提高了收敛率和计算精度。结合典型三维弹性力学问题,对两种施加位移边界条件的方法进行了对比和分析。计算结果表明,直接施加位移边界条件会导致收敛率和精度的明显降低,而基于罚函数法的NURBS有限元分析则能达到最优收敛率,并具有更高的精度。     

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.     

APPLICATION OF PENALTY FUNCTION METHOD IN ISOPARANIETRIC HYBRID FINITE ELEMENT ANALYSIS

IntroductionSince T.H.H.Pain firstly puts forward hybrid element method[1]in1964, the researchand application of hybrid element have got great development. T.H.H.Painet al.havemade important pioneer works in the domain of the research on incompatible displacementstructure hybrid elementstress pattern[2,3]in recentdecades. Reference [4] had provided theoptimizing design concept of hybrid element further and established standardization methodof incompatible displacement structure hybrid eleme…     

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.     

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.     

结构振动分析中的无网格方法

无网格法采用移动最小二乘法构造位移函数，采用罚方法满足本征边界条件，对弹性体的振动问题进行了分析。首先，对权函数中的参数进行了讨论并优化，给出了参数最优值的确定方法；在此基础上对不同边界条件下梁和板的模态进行了分析；最后计算了受突加荷载作用的简支梁以及具有初位移的筒支方板的动力响应。计算结果表明该方法在动力问题的分析中有较高的精度。     

Discontinuous penalty approach with deviation integral for global constrained minimization

In this paper, we use the discontinuous exact penalty functions to solve the constrained minimization problems with an integral approach. We examine a general form of the constrained deviation integral and its analytical properties. The optimality conditions of the penalized minimization problems are proven. To implement the al- gorithm, the cross-entropy method and the importance sampling are used based on the Monte-Carlo technique. Numerical tests show the effectiveness of the proposed algorithm.