自适应多重网格有限元网格生成器研制

基于协调三角形剖分算法、分子表数据结构和Ｚｉｅｎｋｉｅｗｉｃｚ－Ｚｈｕ误估计方法，本文研制出适用于自适应多重网格有限元的网格生成器，该网格生成器可对复杂的区域进行自适应加密。当荷载作用边界随时间变及在动力荷载作用下，网格生成器可退化与再加密网格。     

一种新的并行自适应网格有限元算法

给出了一种新的适用于流体力学问题的并行自适应有限元算法。首先，基于初始稀网格上获得的事后误差估算值，应用反复谱对剖分方法对初网格进行划分，使各子域上总体误差近似相等，从而解决并行自适应计算中的负载平衡问题。然后在各处理器上独立地求解整体问题，并进行指定子域上的网格自适应处理。最后将各子域上的自适应网格组合成一个整体网格，应用基于粘接元技术的区域分裂法在该网格上获得最终解。文末给出了数值实验结果。     

自适应多重网格有限元网格生成器研制

基于协调三角形剖分算法、分子表数据结构和Ｚｉｅｎｋｉｅｗｉｃｚ－Ｚｈｕ误差估计方法，本文研制出适用于自适应多重网格有限元的网格生成器。该网格生成器可对复杂的区域进行自适应加密。当荷载作用边界随时间变化及在动力荷载作用下，网格生成器可退化与再加密网格。     

基于SIMPLER算法的多重网格方法研究

以二维方腔顶盖驱动流为模型,将多重网格方法和SIMPLER算法进行耦合,对不同雷诺数下多重网格加速SIMPLER算法和SIMPLER算法的计算效率进行了对比,数值计算表明:多重网格加速SIMPLER算法不仅能够解决SIMPLER算法不能准确模拟较高雷诺数流场的问题,而且其计算效率远远高于SIMPLER算法.本文也对松弛因子的选取、多重网格实现形式以及网格层数对多重网格加速SIMPLER算法的影响进行了研究,从而为多重网格加速SIMPLER算法的实施提供了计算技术.     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

一种人—机结合的二维有限元网格生成方法

本文提出了一种人一机结合的二维有限元网格生成方法,可将在CAD环境下绘制的矢量型网格图形自动地转换为有限元分析时所用的结点坐标和单元组信息。此方法能完全按人的意图生成网格形状,但比人工划分大大节省时间,而且对计算机硬件要求甚低,在任何PC机上都可应用、     

一种部分连续复合网格生成方法

本文描述了一种构造部分连续复合网络的偏微分方程的数值生成方法，该方法以泊松方程为变换控制方程，构造特定的边界条件，对复杂的蜗尾船型的横剖面进行了变换，变换出来的复合网格与单一区域网格相比，生成的网格质量有明显提高，并计算了另外一个例子，证明了该方法的通用性。     

一种基于有限元网格的车身CAD曲面数据修复方法

针对车身CAD中非同构数据格式间的转换使自由形体的近似值发生变化而导致曲面片之间产生缝隙或网格叠加的情况,提出了一种基于有限元网格、应用于车身CAD曲面网格数据自动修复的方法。该方法还能够在一定程度上局部修复原来划分有缺陷的网格,从而降低对车身CAD模型有限元网格划分的质量要求,并通过工业复杂成形件的实际处理结果验证了该方法的有效性。     

三维有限元六面体网格几何自适应再生成方法

基于几何自适应的六面体网再生成方法,提出了一种以栅格法为基本方法,对有限元网格再划分过程中的网格再划分标准、几何形状的获取以及新旧网格之间物理量的传递等关键问题进行了研究.重点介绍了基于几何自适应的六面体网格再生成算法,首先对旧网格实体模型进行识别,根据实体模型几何特征建立加密源点信息场;然后采用栅格法生成核心网格并对...     

一种特殊设计约束面的有限元网格优化问题

网格设计质量的优劣是有限元分析的关键，对一些特殊的设计约束面，网格的优化设计是建模时需要首先研究解决的问题，本文对两种材料在几何尺寸差异很大情况下的有限元网格进行了优化设计并对其进行了讨论。     

FINITE ELEMENT METHOD ON NUMERICAL SIMULATION OF STRATUM CORNEUM＇S PENETRATION PROPERTY

How the outer substance could penetrate through the skin lies in the stratum corneum, because it is the main barrier in the multi-layers of the skin. Supposing the keratin cell with a special geometry as tetrakaidecahedron, the penetration property of stratum corneum was the key problem which was numerically simulated with finite element method. At first the discretization of the stratum corneum region was given in two steps： first, the discretization of the keratin cell; second, the discretization of fattiness that surrounds the keratin. Then there was the work of numerical simulation. In this procedure, the finite element method and the multi-grid method were used. The former was to obtain the discretization of basic elements; the latter was to decrease the high frequency error. At last the visualization of the numerical simulation was shown.     

微机械杠杆力放大效果的计算和提高

硅微谐振式加速度计是一种高精度的惯性传感器,它通过谐振梁刚度随惯性力的变化检测输入加速度的大小。为了得到较高的标度因数,国内外的研究机构普遍使用微机械杠杆机构来放大惯性力,但实验中微机械杠杆往往达不到理想的放大倍数,尤其是多级杠杆对力的放大作用非常有限。结合国内加工条件,微机械杠杆力放大机构的模型被合理简化,使用有限元分析运动模态的方法,分别计算微机械杠杆、检测质量支撑梁以及谐振梁的刚度,并按此仿真计算结果估算微机械杠杆的力放大倍数。通过理论推导、ANSYS仿真和对微机械杠杆实验的结果验证了这种方法的正确性。为分析微机械杠杆实际的力放大效果提供了一种实用方法,并基于此方法提出通过提高输出刚度来提高微机械杠杆力放大倍数的方法。通过对杠杆结构的优化,可将原有3.9倍的杠杆放大倍数提高到7.0倍。     

Electro-mechanical coupling analysis of MEMS structures by boundary element method

In this paper, we present the applications of Boundary Element Method (BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems (MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics. Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples. The project supported by the 973 Program (G199033108) and the national Natural Science Foundation of China (10125211)     

三维h自适应有限元分析的资料转移

研究了两类用于三维h自适应网络管理的资料转移方法。文中将降低解的精度以节省计算时间的简单方法（简法）与在资料转移中确保解的精度的严格方法（严法）进行了比较分析，并通过数值算例进行论证，结果表明：对于简单问题，两类方法都能满足精度，完成资料转移。对于复杂应力状态问题，严法更可靠，特别是在网络稀疏时，严法明显优越。     

高温MEMS剪应力传感器残余应力分析

采用理论分析与数值模拟相结合的方法,研究了残余应力对静电反馈式MEMS高温剪应力传感器的结构设计的影响.研究表明:在电极宽度和间距不变的情况下,增加电极数和减小电极长度可有效地减小由静电力引起的应力和变形;残余应力对结构应力和变形的影响远大于静电力.     

An adaptive domain decomposition method for simulation of transport in porous media

A Galerkin finite element method is used along with a self-adaptive strategy of domain discretisation to model dispersion in an axisymmetric cylindrical porous medium. A solution strategy is proposed based on the use of a Gear scheme for the time stepping and partial vectorisation of the code. The domain is highly discretised in the area of the sharp transient front, while the remainder is coarsely discretised. The area covered by the fine mesh is determined by the value of the local concentration gradients. Numerical results are presented for the one and two dimensional cases.Nomenclature a _L Longitudinal dispersivity - a _T Transverse dispersivity - [A] Global matrix - C Concentration - D Dispersion coefficient - D _o Diffusion coefficient - {F} Solicitation vector - {G} Solicitation vector - h Finite element length - k Constant - [K] Global matrix of coefficients - [M] Mass matrix - n Outward pointing normal vector - P Pressure - Pe Peclet number - t Time - v Interstitial velocity Greek Symbols Boundary - Test function - Domain     

预处理方法在低速弹箭流场计算中的应用

基于Weiss-Smith预处理矩阵和全局截断预处理参数,采用有限体积方法对雷诺平均Navier-Stokes方程进行离散。对流项离散采用二阶线性重构和AUSM ＋-up格式,时间推进方法采用多重网格下的LU-SGS方法。结合M PI消息传递方法,建立了一套计算低速流动的并行数值方法。计算了低速椭球体的流场和气动力,压力系数和切应力系数计算结果与文献实验结果对比吻合度较好。生成了末敏弹的流场计算网格,对绕末敏弹流场进行了数值模拟。对多重网格下多进程的加速比和并行效率进行了测试,显示了程序良好的并行效率。计算的气动力结果与实验结果吻合。综合结果表明：本文的数值方法能够用于低速弹箭流场和气动力计算,为新型弹箭的设计和定型提供保证。     

岔管结构的图形有限元分析

本文介绍微机辅助水电站岔管结构分析的理论体系及相应的软件系统的结构与特点。该系统由前处理、有限元分析和后处理三部分组成。前处理利用无限变换法或直接法自动生成分级网格和有限元分析属性数据,亦可形成有限元分析的输入文件和有关后处理文件。有限元分析选用精密度高的正交曲线坐标系下的等参曲壳单元,以罚单元解决相贯线处的位移协调;利用只承受压缩变形的压簧元可计及围岩的弹性抗力作用;且具有基于局部重设计(重分析)的再起动功能。后处理以直观的图形方式表示有限元分析的结果,从而实现了有关数据的图形检查和图形分析。     

FINITE ELEMENT METHOD ON NUMERICAL SIMULATION OF STRATUM CORNEUM''''S PENETRATION PROPERTY

How the outer substance could penetrate through the skin lies in the stratum corneum, because it is the main barrier in the multi-layers of the skin. Supposing the keratin cell with a special geometry as tetrakaidecahedron, the penetration property of stratum corneum was the key problem which was numerically simulated with finite element method. At first the discretization of the stratum corneum region was given in two steps: first, the discretization of the keratin cell; second, the discretization of fattiness that surrounds the keratin. Then there was the work of numerical simulation. In this procedure, the finite element method and the multi-grid method were used. The former was to obtain the discretization of basic elements; the latter was to decrease the high frequency error. At last the visualization of the numerical simulation was shown.     

Some improvements on free surface simulation by the particle finite element method

The Lagrangian method has become increasingly popular in numerical simulation of free surface problems. In this paper, after a brief review of a recent Lagrangian method, namely the particle finite element method, some issues are discussed and some improvements are made. The least‐square finite element method is adopted to simplify the solving of the Navier–Stokes equations. An adaptive time method is derived to obtain suitable time steps. A mass correction procedure is imported to improve the mass conservation in long time calculations and time discretization scheme is adopted to decrease the pressure oscillations during the calculations. Finally, the method is used to simulate a series of examples and the results are compared with the commercial FLOW3D code. Copyright © 2008 John Wiley & Sons, Ltd.