双层旋转壳流固耦合振动有限元分析

本文采用旋转壳体元和流体元分析了双层旋转壳体流固耦合振动特性,基于Novozhilov壳体理论、水弹性理论以及Hamilton变分原理推导出耦合系统运动方程,计算结果与实测值比较表明本分析方法有较好的精度。     

A Bubble Element for the Compressible Euler Equations

In this paper, we present a new Galerkin finite element method with bubble function for the compressible Euler equations. This method is derived from the scaled bubble element for the advection-diffusion problems developed by Simo and his colleagues, which is based on the equivalence between the Galerkin method employing piecewise linear interpolation with bubble functions and the Streamline-Upwind/Petrov Galerkin (SUPG) finite element method using P₁ approximation in the steady advection-diffusion problem. Simo and this author have applied this approach to transient advection-diffusion problems by using a special scaled bubble function called P-scaled bubble, which is designed to work in the transient advection-diffusion problems for any Peclet number from 0 to ∞. The method presented in this paper is an application of this p-scaled bubble element to a pure hyperbolic system.     

The Finite Element Analysis of Interaction Problem Through Water,Soil, Balloon and Pile

This paper presents the finite element analysis of an interaction problem involving water, soil, balloon and pile. A building with friction piles is considered, and balloons are introduced to the top of piles. To control the vertical displacement of the building, water is injected into or removed from the balloons. The two-dimensional incompressible Navier-Stokes equation is introduced, and the ALE (Arbitrary Lagrangian Eulerian) method is applied to the water flow analysis. The FS (Fractional Step) method is also applied in the finite element formulation. The soil, which is assumed as a linear elastic body, is subjected to the deformation analysis. The balloon and pile are assumed as a linear elastic truss and a rigid frame, and the deformation analysis is also performed. All the components are discretized by the finite element method in space and are interactively solved by taking into account continuity conditions of traction and displacement.     

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

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

薄板稳定问题的边界元法

本文基于小挠度薄板弯曲问题的基本解,建立了求解薄板稳定问题的边界积分方程,并计算了若干算例,结果表明用边界元法求解薄板的稳定问题是行之有效的.     

线性船舶兴波阻力边界元新方法

采用常数边界元对船舶与流体界面进行离散，求解船舶兴波势及船舶兴波阻力。这种方法可避免在船舶与流体自由面交线上安置节点，因而避免了这些节点建立补充方程。因为满足自由面条件的Ｈａｖｅｌｏｃｋ源函数的源点和场点不能同时在自由面上，使得自由面上的节点无法用Ｈａｖｅｌｏｃｋ源函数的建立方程。如对自由面交线上的节点建立补充方程，则要对线性自由面条件中包含的未知势函数的二阶导数用差分形式表达，引入较大误差。     

一种利用静力试验数据修正有限元模型的方法

基于代数特征值逆问题理论，提出了一种利用静力试验数据修正有限元模型方法。由于修正后的模型为满阵，根据有限元模型的力学特征对修正后的模型加以约束，使其不仅是对称的，而且还保持稀疏、带状特性，其计算过程简洁，计算结果表明该方法可靠，效果令人满意。     

流体力学反命题变分有限元非线性规划法

本文在用有限元法求解变域问题时，利用变域变分公式将问题化为求泛函的驻值，但该驻值常常是一个鞍点，因此给计算带来不便，我们在这里利用最小二乘法将求驻值的问题转化为求另一泛函的极值问题，然后采用LMF方法和共轭梯度法求解极值问题。     

弹性力学轴对称问题的有限元线法

给出了解弹性力学空间轴对称问题的有限元线法的基本理论。该法包括了２－４条结线的等参数单元，沿结线方向的两点边值问题采用插值矩阵法解之。算例表明，本法具有良好的收敛性和较高的计算精度。     

A Finite Element Method for 3-D Eddy Current Problems Using an Iterative Approach

An iterative method is proposed for a finite element approximation of three-dimensional eddy current problems. The method is based on an iterative method derived from a perturbation problem of magnetostatic problems. The TEAM model and a transformer model are considered as numerical examples. In both examples, BiConjugate Gradient (BiCG) method is applicable for the complex symmetric linear system arising in each step of the iterative procedure for a rather wide range of the perturbation parameter, and the present results seem to be suitable.     

General Spatial Dynamic Problem for an Elliptic Crack under the Action of a Normal Shear Wave,with Consideration for the Contact Interaction of the Crack Faces

Consideration is given to the contact interaction of the faces of a stationary plane elliptical crack under the action of a harmonic shear wave normally incident on the crack surface. The dependence of the mode II and III stress intensity factors on the wave number is studied for different values of the friction coefficient