流体饱和两相多孔介质拟静态问题的有限元解法

给出基于混合物理论的流体饱和两相多孔介质模型,该模型由一可变形固体 一流体相组成。采用Ｇａｌｅｒｋｉｎ加权残值法导出求解拟静态问题的有限元公式,并编制了二维有限元程序。用程序分析了一维和二维问题,得到合理的结果。     

A new Finite Element Method for magneto-dynamical problems: two-dimensional results

A mixed Lagrange finite element technique is used to solve the Maxwell equations in the magneto-hydrodynamic (MHD) limit in an hybrid domain composed of vacuum and conducting regions. The originality of the approach is that no artificial boundary condition is enforced at the interface between the conducting and the insulating regions and the non-conducting medium is not approximated by a weakly conducting medium as is frequently done in the literature. As a first evaluation of the performance of the method, we study two-dimensional (2D) configurations, where the flow streamlines of the conducting fluid are planar, i.e., invariant in one direction, and either the magnetic field (“magnetic scalar” case) or the electric field (“electric scalar” case) is parallel to the invariant direction. Induction heating, eddy current generation, and magnetic field stretching are investigated showing the usefulness of finite element methods to solve magneto-dynamical problems with complex insulating boundaries.     

Preconditioned Iterative Methods for Linear Equations Arising from the Boundary Element Method

A precondition for the Gauss–Seidel iterative method to solve a linear system of equations arising from the boundary element method for the Laplace and convective diffusion with first-order reaction problems is presented in this paper. The present precondition is based on the elementary matrix operation. We discuss the effect of the precondition in comparison with the Gauss elimination (GE) method in some numerical experiments.     

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

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

A Hybrid Boundary/Finite Element Method for Simulating Viscous Flows and Shapes of Droplets in Electric Fields

A mixed boundary element and finite element numerical algorithm for the simultaneous prediction of the electric fields, viscous flow fields, thermal fields and surface deformation of electrically conducting droplets in an electrostatic field is described in this paper. The boundary element method is used for the computation of the electric potential distribution. This allows the boundary conditions at infinity to be directly incorporated into the boundary integral formulation, thereby obviating the need for discretization at infinity. The surface deformation is determined by solving the normal stress balance equation using the weighted residuals method. The fluid flow and thermal fields are calculated using the mixed finite element method. The computational algorithm for the simultaneous prediction of surface deformation and fluid flow involves two iterative loops, one for the electric field and surface deformation and the other for the surface tension driven viscous flows. The two loops are coupled through the droplet surface shapes for viscous fluid flow calculations and viscous stresses for updating the droplet shapes. Computing the surface deformation in a separate loop permits the freedom of applying different types of elements without complicating procedures for the internal flow and thermal calculations. Tests indicate that the quadratic, cubic spline and spectral boundary elements all give approximately the same accuracy for free surface calculations; however, the quadratic elements are preferred as they are easier to implement and also require less computing time. Linear elements, however, are less accurate. Numerical simulations are carried out for the simultaneous solution of free surface shapes and internal fluid flow and temperature distributions in droplets in electric fields under both microgravity and earthbound conditions. Results show that laser heating may induce a non-uniform temperature distribution in the droplets. This non-uniform thermal field results in a variation of surface tension along the surface of the droplet, which in turn produces a recirculating fluid flow in the droplet. The viscous stresses cause additional surface deformation by squeezing the surface areas above and below the equator plane.     

The Space-Time Finite Element Method for Parabolic Problems

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｅｑｕａｔｉｏｎｓｗｅｃｏｎｓｉｄｅｒｅｄａｒｅａｓｆｏｌｌｏｗｓｕｔ-Δｕ =ｆ(ｕ) ,　　Ω× [0 ,Ｔ] ,ｕ｜ Ω =0 ,　　　　　 Ω× [0 ,Ｔ] ,ｕ( · ,0 ) =ｕ0 ,Ω ,( 1 )ｗｈｅｒｅΩ ∈Ｒ2 ,ｔｈｅｆｕｎｃｔｉｏｎｆ(ｕ)ｓａｔｉｓｆｉｅｓ｜ｆ(ｕ)｜≤ｃ｜ｕ｜ ,　　 ｕ∈Ｃ(Ω) . ( 2 )Ａｎｄｆ(ｕ)ｉｓＬｉｐｓｃｈｉｔｚｃｏｎｔｉｎｕｏｕｓ,ｉ.ｅ.ｉｔｓａｔｉｓｆｉｅｓ｜ｆ(ｕ) -ｆ(ｖ) ｜≤Ｌ｜ｕ-ｖ｜ ,　　 ｕ ,ｖ∈Ｃ(Ω) ,( 3 )ｗｈｅｒｅＬｉｓＬｉｐｓｃｈｉｔｚｃｏｎｓｔａ…     

利用有限元法计算径向磁轴承性能

利用非线性有限元方法(FEM)计算了一种径向磁轴承(定子具有8磁极常规结构)的转子在无偏心以及偏置电流分别为0.5 A、1.0 A和1.5 A下的电磁力及其电流刚度与位置刚度,预测了径向磁轴承的承载力,给出了其磁场分布图,并与磁轴承的线性化模型参数进行了比较.结果表明,基于有限元方法得到的计算值同试验结果一致,可以为径向磁轴承的设计和分析提供依据.     

移动荷载作用下地基动力分析的有限元方法

通过对地基动力问题的基本方程进行变换,把基本方程变换到随荷载移动的运动坐标系中,通过加权残数法推导了相应的单元刚度矩阵,从而建立了移动问题的有限元格式,并发现移动荷载问题的单元刚度矩阵是对相应静力问题单元刚度矩阵的修正,在静力单元刚度矩阵的主对角元素上增加与移动速度有关的项,即可得到移动问题有限元的单元刚度矩阵,这样就将动力学问题转化为“拟静力”问题处理。文中用移动问题有限元方法计算了地基的动力响应,并与解析解进行了对比,以说明本方法具有较好的精度。     

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

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

水闸闸室结构三维样条边界元分析

用三维样边界元法分析水闸闸室结构。底板，闸墩和载水墙等为其子结构，交通桥，工作桥和胸墙等处理为内部支撑。地基和边载可以是任意的，只要能给定地表位移面力关系。在各种工况下，不论是设置平板门还是弧形门，是平底板是反拱底板，即使在稀疏剖分下也能给出高精度的位移场，应力场和地基反力场。     

A Subgrid-Scale Stabilized Finite Element Method for Multicomponent Reactive Transport through Porous Media

Standard Galerkin finite element methods (GFEM) lack stability in solving advection-dominated solute transport in porous media. They usually require prohibitively fine grids and extremely small time steps to solve for advection-dominated problems. The algebraic subgrid-scale stabilized (ASGS) finite element method has been proved to overcome such problems for single-species reactive transport. Its potential for dealing with multicomponent reactive transport has not yet been explored. Here we present a numerical formulation of ASGS for steady and transient multicomponent reactive transport. Subgrid-scale transport equations are solved first by using an ASGS approximation and their solutions are substituted back into the grid-scale equations. A sequential iteration approach (SIA) is used to solve for coupled transport and chemical equations. Coupling of ASGS and SIA, ASGS+SIA, has been implemented in a reactive transport code, CORE^2D V4, and verified for conservative solute transport. ASGS+SIA has been tested for a wide range of 1-D transient multicomponent reactive transport problems involving different types of chemical reactions such as: (1) Kinetically controlled aqueous species degradation, (2) Kinetic mineral dissolution, (3) Serial-parallel decay networks, and (4) Cation exchange and pyrite oxidation at local equilibrium. ASGS+SIA always provides accurate solutions and therefore offers an efficient option to solve for advection-dominated multicomponent reactive transport problems.     

不可压缩和近乎不可压缩粘弹性问题的有限单元法

本文首先从Herrmann泛函出发,导出了不可压缩和近乎不可压缩粘弹性问题的新的本构关系,然后根据虚功原理建立了相应的有限元公式。最后编制了平面问题的有限元分析程序VFAPINP,并计算了实例。计算结果表明所介绍的分析方法和计算程序解决了泊松比趋近或等于0.5时粘弹性应力分析问题,特别适用于固体推进剂的应力分析。     

液体三维晃动特征问题的有限元数值计算方法

本文采用有限元方法数值求解任意刚性容器内液体三维晃动的固有频率和模态。通过建立液体晃动特征问题的泛函极值原理，编制了四面体等参单元有限元程序，计算了平放圆柱腔内三维液体晃动的特征频率，并将矩形容器、球腔、带“十”字隔板球形容器内的液体三维晃动计算结果与解析解、实验结果和二维有限元数值解进行了比较，程序的正确性得到了验证。     

波纹管在子午面内整体弯曲的半解析有限元法

波纹管是一类子午线呈波纹状的旋转壳，作为弹性敏感元件和柔性连接件在航空仪表和管道工程中起着重要的作用。长期以来基于壳体理论的分析多限于轴对称变形问题。最近，虽然出现了以柔性旋转壳理论为基础的解决波纹管整体弯曲问题的解析解和数值解，但仍有必要通过别的途径加以验证和补充。为些，本文提出了波纹管在子午面内整体弯曲的半解析有限元解。把波纹管近似成有限个截顶锥壳的组合体，每个截顶锥为一个单元。将位移分量沿纬线用Fourier级数展开，沿子午线用多项式插值，截锥单元因此化为2节点的直线单元，每节点为4个自由度。显然，对于同等的RAM和CPU，线单元可比其他单元划得更小。于是，利用小单元条件，将一个单元内的壁厚和平行圆半径近似地当作不变量，给出了用显式表示的单元刚度矩阵，为直观分析结构参数对波纹管力学性能的影响提供了方便。在此基础上，计算了Ω型，C型和U型波纹管在纯弯矩作用下的变形和应力分布，所得结果和已有的解析解，数值解相符。本法不限于波纹管的计算分析。     

点式压电智能结构的模型修正与振动主动控制

提出了在点式压电智能结构中应用摄动有限元方法对结构的有限元模型进行修正 ,从而达到提高建模精度 ,改善实际结构振动主动控制效果的目的。通过对一悬臂梁在模型修正前后进行振动主动控制的不同的控制效果验证了该方法的有效性     

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

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

复合材料层合壳有限元分析的预测-修正法

对于复合材料层合壳的有限元分析，本文根据Reissner-Mindlin型的全局位移场给出了一个预测一修正法。首先按照一般的有限元分析过程(没有引入剪切修正系数)计算出全局响应(如挠度，频率和屈曲载荷等)的预测值；然后利用Lagrange插值构造横向剪应力的一般形式，使得满足层间连续和表面上为零的条件，通过最小二乘法拟合三维应力平衡方程获得横向剪应力；最后在单元上计算和引入剪切修正系数，再经过有限元分析计算出全局响应的修正值。     

A 3D direct coupling method for steady ship hydroelastic analysis

Asides from the influence of incoming waves, ships can experience steady motions, such as rigid-body sinkage and trim motions, and flexible-body vertical bending motions, due to a constant forward speed even under calm water conditions. In this paper, a novel approach to analyze steady-ship hydroelasticity, particularly for the steady-ship motions and surrounding steady-wave disturbances, is proposed using a three-dimensional (3D) direct coupling method, based on a higher-order boundary element method (HOBEM) and a higher-order shell finite element method (FEM). Within the linearized framework, a solution method is proposed based on a two-step procedure, using two types of Neumann–Kelvin (NK) linear flow models for the fluid part and a virtual work equilibrium equation for the structural part. The first step is to compute a mean position wave-resistance problem using the modified NK equation, the second step is to solve a perturbed position wave-resistance problem, by employing a classical NK model and a virtual work equation based on the first step’s solution. Detailed mathematical formulation and numerical procedures are described, and a few numerical results are illustrated. These include both rigid and flexible steady-ship motions, Von-Mises stress distributions, and wave-resistance coefficients for Froude numbers ranging from 0.15 to 0.5. Furthermore, the numerical results obtained using the present direct coupling method and a modal-based one are compared.     

A stabilized formulation for the advection–diffusion equation using the Generalized Finite Element Method

This paper presents a stable formulation for the advection–diffusion equation based on the Generalized (or eXtended) Finite Element Method, GFEM (or X‐FEM). Using enrichment functions that represent the exponential character of the exact solution, smooth numerical solutions are obtained for problems with steep gradients and high Péclet numbers in one‐ and two‐dimensions. In contrast with traditional stabilized methods that require the construction of stability parameters and stabilization terms, the present work avoids numerical instabilities by improving the classical Galerkin solution with enrichment functions (that need not be polynomials) using GFEM, which is an instance of the partition of unity framework. This work also presents a strategy for constructing enrichment functions for problems involving complex geometries by employing a global–local‐type approach. Representative numerical results are presented to illustrate the performance of the proposed method. Copyright © 2010 John Wiley & Sons, Ltd.     

A New Method for Solution of 3D Elastic-Plastic Frictional Contact Problems

ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｔａｃｔｐｒｏｂｌｅｍｉｓｏｆｐａｒｔｉｃｕｌａｒｉｎｔｅｒｅｓｔｉｎｎｕｍｅｒｏｕｓｅｎｇｉｎｅｅｒｉｎｇａｐｐｌｉｃａｔｉｏｎｓｒａｎｇｉｎｇｆｒｏｍｓｔｒｕｃｔｕｒｅ＿ｓｔｒｕｃｔｕｒｅｉｎｔｅｒａｃｔｉｏｎｔｏｍａｃｈｉｎｅｄｅｓｉｇｎａｎｄｍｅｔａｌｆｏｒｍｉｎｇｕｎｄｅｒｖａｒｉｏｕｓｌｏａｄｃａｓｅｓ.Ｃｏｎｔａｃｔｐｒｏｂｌｅｍｓａｒｅｎｏｎｌｉｎｅａｒｆｏｒｔｈｅｍｏｖｉｎｇｂｏｕｎｄａｒｙａｎｄｔｈｅｅｘｉｓｔｅｎｃｅｏｆｆｒｉｃｔｉｏｎａｌｏｎ…