应用变域变分有限元法解平面叶栅反命题

反命题作为一种可变(未知)边界问题近年来得到了广泛的研究。本文给出了亚声速平面叶栅反命题计算的势函数变域变分有限元解法。变域变分通过把可变边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而使可变边界求解和流场分析可以完全耦合进行。本文针对亚声速平面叶栅的反命题，根据泛函的驻值必要条件，介绍了变域变分有限元方法的求解过程，最后给出了两个数值算例。这两个算例均采用四节点矩形单元的插值基函数，第一个算例用于检验程序的可靠性，第二个算例设计了一个给定叶面马赫数分布的叶型，并与试验结果进行对照。     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.     

线性区间有限元静力控制方程的组合解法

区间有限元的静力控制方程常被归结为区间方程组来求解。但实际上两者并不等价。本文根据不确定结构有限元分析的力学背景，直接从问题的基本参量的不确定性出发，将基本区间参量的边界组合与求解区间方程组的有关解法相结合，提出了线性区间有限元静力控制方程的两种组合解法－参量边界全组合法和组合迭代法。可以以较小的计算量获得或逼近位移和应力区间的准确界限。且不受基本参量变化范围的限制。算例分析表明文中方法是实用和可行的。     

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

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

用面力边界积分方程求解断裂力学J积分

证明面力边界积分方程被积函数的散度等于零，应用Stokes公式，对平面线弹性问题，将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果，推导出J积分亦可表示为边界点的积分势函数的点值计算，无需进行数值积分，实例计算说明该方法的有效性。     

用间接边界元方法求解不可压缩孤立翼型绕流的一种库塔条件处理方法

本文对不可压二维翼型势流绕流的边界元法求解作了分析,在对基本方程与边界条件进行数值离散化时,将库塔条件代入基本方程。按本方法编制的计算机程序对若干算例进行了验算。结果表明,本文提出的方法是可靠的,该法计算简捷、方便,占用计算机内存少,具有实用意义。     

非定常可压缩紊流边界层流场的计算

本文通过引入非定常流函数,在Levy—Lees坐标下导出了非定常可压边界层的动量方程,并以Box格式将它离散并求解。而非定常能量方程则在物理坐标网格点上直接进行离散,其解法较为简便。采用本文方法首先计算了一个实际的低速非定常紊流边界层;并分别与实验结果及前人的计算进行了对比。本文的计算预示了较大的速度振幅过冲量。在可压流范围内。本文计算了零攻角平板叶栅、零攻角和5°正攻角双凸叶栅三种情况在进口来流简谐扰动条件下的非定常紊流边界层。计算结果表明,递压梯度下的平均速度剖面中出现涡量向壁面的迁移现象。     

TVD格式在超音速喷管三维粘性流动求解中的应用

详细给出了任意三维曲线坐标系中Ｎｏｖｉｅｒ－Ｓｔｏｋｅｓ方程的对流项ＴＶＤ格式的构造过程，建立了数值求解三维粘性流动的计算方法，应用该方法对三维超音速喷管中有激波及无激波情况下的两种工况的层流流场进行了数值求解，并与实验做了对比。结果表明本文建立的计算方法具有较高的精度，同时也证明ＴＶＤ格式具有分辩率高，稳定收敛等优点，为进一步开展叶栅流场及紊流的研究打下了基础。     

透平叶栅跨音速流动计算中的新型有限元法

将Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法和多级有限元的思想结合起来，构成了在收敛速度和稳定性两方面均较好的新型有限元算法：多级广义有限元。利用这一方法，分别基于Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和Ｅｕｌｅｒ方程，研究了透平跨音速叶栅无粘流动和粘性流动，并将计算结果与实验结果作了比较。计算结果表明，本方法是透平机械内部跨音速流动计算的强有力的手段。     

跨声速势函数的构成与数值求解

现有的非等熵势函数方法是基于热力学第二定律建立的密度与熵之间的关系式,由动量方程间接或直接计算出熵的变化,再与势函数方程迭代求解。与这种复杂的迭代计算不同,本文建议直接从散度型动量方程计算包含了熵的影响的密度,其与势函数方程迭代便可方便地求解整个流场,采用人工密度和AF2方法,对任意迴转面叶栅内跨声速流动计算的结果表明,本文方法得到的流场与等熵计算相比,激波位置前移,强度减弱,和试验结果较为接近;与计算熵增的非等熵势函数相比,结果是相同的,而计算却比较简单。     

成层半空间出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性水平成层半空间中自由波场时域计算的一维化有 限元方法. 在进行有限元网格划分时，竖向单元取满足有限元模拟精度的任意尺寸，水平向 网格尺寸由时间离散步长和水平视波速确定，并自动进行虚拟网格划分. 基底设置人工边界， 并将波动输入转化为等效荷载施加在边界节点上. 然后将集中质量有限元法和中心差分法相 结合建立节点运动方程，并将水平方向相邻节点的运动用该节点相邻时刻的运动表示，从而 将求解节点运动的二维方程组转化为一维方程组. 求解此方程组，即得到自由场中竖向一列 节点的运动. 最后根据行波传播的特点，可方便地确定全部自由波场. 理论分析和数值算例 表明，该方法具有较高的精度和良好的稳定性.     

含液多孔介质力学问题的边界元方法

提出了一种含液多孔介质力学问题的边界元求解方法.首先将问题分解为一系列含单孔流体夹杂的子问题,然后针对每个子问题建立了流体孔体积变化率与流体压力之问的函数关系,进一步采用边界元方法建立了以各流体孔压力为基本未知量的线性代数方程组,最后根据所求出的各流体孔的压力计算含液多孔介质内各点的位移、变形和应力.为了说明方法的有效...     

On kinematic coupling equations within the scaled boundary finite-element method

The scaled boundary finite element method is a semi-analytical analysis technique, which combines the advantages of the finite element method and the boundary-element method. Assuming that the geometry of the governing structure can be represented by mapping its boundary with respect to the so-called scaling coordinate, the problem can be handled in a closed-form analytical manner in the scaling direction and by a finite-element approximation in the other directions. Thus, a discretization of the boundary is sufficient and the nodal degrees of freedom are functions of the scaling coordinate. In some situations, such as the analysis of the free-edge effect in laminated plates, it is useful to introduce kinematic coupling equations, which are valid not only on the boundary, but also within the domain. The implementation of linear kinematic coupling equations within the method is presented for the case of a three-dimensional structure with scaling in a fixed Cartesian direction. Rigid-body modes are handled by using the concept of generalized inverse matrices. In some benchmark examples the efficiency of the approach is demonstrated and comparison with the results of the finite-element method shows good accordance.     

Boundary element analysis of anisotropic Kirchhoff plates

In this paper, the radial integration method is used to obtain a boundary element formulation without any domain integral for general anisotropic plate bending problems. Two integral equations are used and the unknown variables are assumed to be constant along each boundary element. The domain integral which arises from a transversely applied load is exactly transformed into a boundary integral by a radial integration technique. Uniformly and linearly distributed loads are considered. Several computational examples concerning orthotropic and general anisotropic plate bending problems are presented. The results show good agreement with analytical and finite element results available in the literature.     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

Transient heat conduction analysis using the NURBS-enhanced scaled boundary finite element method and modified precise integration method

The Non-uniform rational B-spline(NURBS)enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper.The scaled boundary finite element method is a semi-analytical technique,which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction.In this method,only the boundary is discretized in the finite element sense leading to a reduction of the spatial dimension by one with no fundamental solution required.Nevertheless,in case of the complex geometry,a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often unavoidable in the conventional finite element approach,which leads to huge computational efforts and loss of accuracy.NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape.In the proposed methodology,the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions,while the straight part of the boundary is discretized by the conventional Lagrange shape functions.Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analysis and the solution is obtained using the modified precise integration method.The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion.Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method.The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.     

Potential flow around obstacles using the scaled boundary finite‐element method

The scaled boundary finite‐element method is a novel semi‐analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. The method works by weakening the governing differential equations in one co‐ordinate direction through the introduction of shape functions, then solving the weakened equations analytically in the other (radial) co‐ordinate direction. These co‐ordinate directions are defined by the geometry of the domain and a scaling centre. The method can be employed for both bounded and unbounded domains. This paper applies the method to problems of potential flow around streamlined and bluff obstacles in an infinite domain. The method is derived using a weighted residual approach and extended to include the necessary velocity boundary conditions at infinity. The ability of the method to model unbounded problems is demonstrated, together with its ability to model singular points in the near field in the case of bluff obstacles. Flow fields around circular and square cylinders are computed, graphically illustrating the accuracy of the technique, and two further practical examples are also presented. Comparisons are made with boundary element and finite difference solutions. Copyright © 2003 John Wiley & Sons, Ltd.     

梯度材料中矩形裂纹的对偶边界元方法分析

采用对偶边界元方法分析了梯度材料中的矩形裂纹. 该方法基于层状材料基本解,以非裂纹边界的位移和面力以及裂纹面的间断位移作为未知量. 位移边界积分方程的源点配置在非裂纹边界上,面力边界积分方程的源点配置在裂纹面上. 发展了边界积分方程中不同类型奇异积分的数值方法. 借助层状材料基本解,采用分层方法逼近梯度材料夹层沿厚度方向力学参数的变化. 与均匀介质中矩形裂纹的数值解对比,建议方法可以获得高精度的计算结果. 最后,分析了梯度材料中均匀张应力作用下矩形裂纹的应力强度因子,讨论了梯度材料非均匀参数、夹层厚度和裂纹与夹层之间相对位置对应力强度因子的影响.      

自由单元法及其在结构分析中的应用

通过吸收有限元与无网格法的优点,提出了一种新的数值方法------自由单元法.此方法在离散方面,采用有限元法中的等参单元,表征几何形状和进行物理量的插值;在算法方面,采用单元配点技术,逐点产生系统方程.主要特点是,在每个配置点只需要一个和周围自由选择的节点而形成的一个独立的等参单元,因而不需要考虑物理量在单元之间的相互连接关系与导数连续性问题. 本文介绍强形式与弱形式两种自由单元法,前者直接由控制方程和边界条件直接产生系统方程,后者通过在自由单元上建立控制方程的加权余量式产生弱形式积分式,并通过像传统有限元法中的积分过程建立系统方程组.本文提出的方法是一种单元配点法,对于域内点为了获得较高的导数精度,需要采用至少具有一个内部点的等参单元,为此除了可使用各阶次的拉格朗日四边形单元外, 还 给出了七节点三角形等参单元,用于模拟较为复杂的几何形状问题.