有限元网格自动生成及绘制中的实用处理技巧

提出了超大单元，结点归一化相对距离判别法这一快速三维自动部分技术，针对大型工建筑物的特点，给出了几种常遇问题的实用处理方法，还给出了适合工程实际的超大判断立体消隐图绘制法和基于Ａｕｔｏｌｉｓｐ语言及Ａｕｔｏ－ＣＡＤ的立体图处理方法。     

数值模拟离子强度敏感水凝胶的多场特性

发展了多物理模型来研究溶液中离子强度敏感水凝胶电-力-化学多场耦合的特性. 模 型的主要控制方程包括：计算水凝胶内外离子浓度分布的Nernst-Planck化学场方程；描述 膨胀变形的力学场方程和描述电场的泊松方程. 无网格有限云团法和牛顿迭代法用来数值离 散和求解控制方程. 通过对比多场耦合的响应, 包括胶的膨胀率和胶内外离子浓度和电势的分布，探讨了影响胶体变形的主要因素. 对数值模拟结果和实验结果进行了对照.     

无网格算法在多段翼型流动计算中的应用

研究了一种求解欧拉方程的无网格算法，发展出了一套布点及点云自动生成的方法；在点云离散的基础上，采用最小二乘法求解矛盾方程的方法来求取空间导数，进而获得数值通量；采用四步龙格-库塔方法进行时间推进，并引入当地时间步长和残值光顺等加速收敛措施。通过对NA-CA0012翼型的跨音速流动和多段翼型复杂绕流的数值模拟，验证了上述无网格算法的正确性和实用性。     

改进的小损伤结构动力学有限元建模方法

在遗传算法或神经网络方法识别结构损伤位置和程度时,都是基于少量的在线测量的损伤结构振动响应数据和大量的模型仿真数据来实现的,因而建立高效和精确的损伤结构动力学有限元模型,以便仿真获得损伤结构的大量动力学响应数据是十分重要的基础前提工作。本文针对ANSYS结构分析软件在建立结构小损伤有限元动力学模型存在两个关键问题,结构损伤处直接网格划分的计算结果误差和网格节省问题,以结构损伤振动检测的实际需要为出发点,提出了建立含小损伤结构的ANSYS动力学建模技术,研究了结构局部小损伤及其位置与所在处单元刚度矩阵变化的数量关系。     

有限元计算中疏密网格过渡方法研究

工程计算中出于节省计算量的目的,往往需要在一个有限元模型中布置粗细不同的网格。为保证计算结果的准确性,必须保证网格突变情况下的位移协调问题。本文工作之一是在强天驰界面过渡单元的基础上,引入虚拟节点和子单元,在子单元中应用节理元思想,提出了基于最小势能原理的弹簧节理单元法。简化了积分运算,避免了精度要求极高的坐标转换,从而提高了方法的精度和实用性;二是提出了基于位移约束的主从自由度法,简便实用,只需简单的矩阵运算即可实现。两种方法均实现了不同尺寸网格间位移的协调性和刚度的匹配,从而使之满足有限元收敛准则,且生成的刚度阵具有对称性及带状性。算例证明两种方法精度良好,并可方便地应用于求解大规模工程问题。     

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

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

复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

影响无网格方法求解精度的因素分析

基于移动最小二乘法的无网格方法的计算精度除受到节点的分布密度和基底函数的阶次影响外，还受到其它因素的影响，其中权函数的选取、权函数影响域的大小及位移边界条件的引入对计算精度影响较大。本文分析了几种常用权函数在数值计算时的特点，包括计算精度、收敛情况、计算效率等，同时分析了影响域大小及边界条件的引入对计算精度的影响。通过分析给出了确定权函数及其影响域大小的方法。当受约束的自由度较多时，通过配点法引入位移边界条件会引起计算结果的振荡，通过施加稳定项可以消除振荡现象，通过对带孔方板的受力分析证明了其可行性。应用以上结论对J23—10曲柄压力机机身进行了受力分析，应力集中部位的计算结果得到了较高的精度。     

完全变换法在无网格伽辽金方法中的应用

由于移动最小二乘形函数一般不具有常规有限元或边界元形函数所具有的插值特征，本质边界条件的处理成为无网格伽辽金法实施中的一个难点。本文通过建立节点位移和广义位移之间的关系对移动最小二乘形函数进行修正，给出了修正的移动最小二乘形函数；以二维问题为例，对完全交换法在无网格伽辽金方法中的应用进行了研究，实现了本质边界条件在节点处的精确施加。数值计算结果表明该方法不仅简单合理，而且具有较高的精度、收敛性和稳定性。     

利用无网格方法分析斜拉破坏钢筋混凝土梁

从传统的混凝土分离裂缝模型和弥散裂缝模型出发，使用不同的模型模拟不同发展阶段的混凝土裂缝，并利用无网格方法可以灵活布置节点和边界的优点，使分离裂缝模型的使用得到很大的简化。通过在宏观裂缝表面布置基于试验结果的混凝土与混凝土界面单元，正确地模拟了裂缝表面抗剪能力的变化。算例表明，使用本文方法可以准确模拟斜拉破坏混凝土梁的破坏过程，且精度和数值稳定性要高于传统的有限元弥散裂缝方法。