重载滑动轴承非牛顿流体的三维弹性流体动力润滑分析的加权平均法

本文采用加权平均法成功地求解了由广义Ｒｅｙｎｏｌｄｓ方程（幂次型）和三维弹性方程组成的非牛顿体弹流模型，该方法较好的解决了大偏心率下的弹流迭代收敛困难的问题。该方法具有收敛快，对初值要求松且计算精度高等优点，是一种解决弹流问题较好的迭代模式。     

基于半光滑牛顿法的润滑液膜有限元空化算法

针对润滑液膜中空化问题,引入Fischer-Burmeister函数,提出一种求解满足质量守恒雷诺方程的半光滑牛顿迭代算法.该算法将空化问题的非线性互补关系转化为等式约束方程,避免了迭代计算中的不等式约束识别问题.算法可将空化约束方程与雷诺方程、力平衡方程、变形方程等同时纳入牛顿迭代方程组,有效解决了传统松弛迭代算法需要多重嵌套循环带来的效率低下问题及压力与膜厚的强耦合性带来的收敛困难问题.计算实例表明,该算法计算效率高、收敛性好,且易应用于弹流润滑分析中,在滑动轴承和机械端面密封等多种物理模型下均有良好的适用性.     

Hahn-Tsai复合材料的非线性杂交应力有限元方法

成功建立了Hahn-Tsai复合材料模型的非线性杂交应力有限元方程,采用Newton-Raphson迭代法求解结构的非线性位移方程。在迭代过程中,为了提高计算效率可采用简单迭代法由节点位移求解单元应力场。但是,当载荷增加到一定程度以后,非线性应力场由于循环迭代而无法收敛,显然,一般的加速方法不能解决这种循环迭代的发散问题。因此,本文发展了一种确实有效的非线性应力场迭代新方法,在不增加计算工作量的情况下,不仅极大地提高了收敛速度,而且对于较大载荷也能够很好地收敛,从而解决了大载荷下非线性杂交元方法失败的关键问题。数值算例表明该方法是确实可行的。     

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

粘弹流体轴对称流动的格子Boltzmann方法

基于BGK碰撞模型，通过在迁移方程中引入作用力项，建立了粘弹流体的轴对称格子Boltzmann模型．通过Chapman-Enskog展开，获得了准确的柱坐标下轴对称宏观流动方程．采用双分布函数对运动方程和本构方程进行迭代求解，模拟分析了粘弹流体管道流动，获得了流场中的速度和构型张量的分布，通过与解析解进行比较，验证了模型的准确性．研究了作为粘弹流体流动基准问题的收敛流动，对涡旋位置进行了定量分析，将回转长度的计算结果与有限体积法进行了比较，两种数值结果十分吻合．研究结果表明，模型能够准确表征粘弹流体的轴对称流动，具有较广阔的应用前景．     

一种快速求解非稳态弹流问题的直接迭代算法

本文提出一种求解弹流润滑问题的直接迭代方法。实践证明,这种算法具有较好的收敛性,极快的运算速度和相当宽的参数适应范围,特别适用于非稳态问题的研究。新的直接迭代格式还具有程序简单、占用内存少等普通直接迭代法的优点。     

回转体局部空泡绕流的非线性分析

基于面元法, 通过在回转体和空泡壁面放置源汇,对回转体定长局部空泡的绕流问题进行了分析和讨论,并提出了求解回转体局部空泡绕流“正问题”的方法。计算结果表明：所给出的方法具有快速收敛的特征,第1次迭代和最终收敛时空泡壁面切向速度的误差不超过5％;随着回转体面元总数N的增加,局部空泡的空泡数σ趋于稳定;通过比较可知,该方法得到的理论估算值与实测值的一致性较好。     

采空区流场非达西渗流一种新的迭代算法

对Bachmat非线性渗流方程求解问题，结合采空区问题的力学特点，提出基于变渗透率达西 (Darcy)渗流求解非线性渗流的迭代算法. 建立了冒落采空区非线性漏风流态的有限元数 值模型，通过对复杂边界采空区的漏风渗流的计算，得到与实际流态更接近的流动规律(风 压等值线和流函数线)和采空区漏风强度分布(速度场). 结果表明，迭代算法是振荡性收 敛的，收敛速度快；与Darcy渗流相比，采空区非线性渗流速度场趋向于平缓，在工作面附 近30\,m范围内最大风速降低到原来的0.61倍，其余大部分区域速度略有增加. 迭代方法满 足工程要求.     

面向对象的土石坝参数随机反演程序设计

将储液容器流固耦合系统中的液体和容器分别视为理想可压缩流体和线弹 性固体，采用流体压力单元和固体壳单元对流固耦合系统进行有限元离散，得到一个非对称 的大型流固耦合有限元方程. 采用Arnoldi方法求解上面这个大型有限元方程的非对称特征 值问题，以得到储液容器的动力特性. 通过移频技术避免了处理零频问题，并构造了迭代格 式计算Arnoldi向量. 数值算例表明所用解法对于流固耦合系统都是非常有效的.     

半空间弹性波双参数强非线性逆散射

半空间弹性波强散射源逆问题是一个实用而又难度很大的强非线性问题，而经典的Born迭代方法只适用于弱非线性问题．基于DBI方法（Distorted-Born Iterative方法），提出了一种求解弹性波强非线性逆散射问题的迭代方法．在数值模拟运算时利用矩量法进行离散处理，并采用正则化原则避免求解病态矩阵方程。计算结果表明该方法具有较快的收敛速度及较高的精度。     

基于格子Boltzmann方法的粘弹性流体圆柱绕流分析

将光滑界面法引入到格子Boltzmann方法中分析粘弹性流体绕流问题，分别采用单松弛模型和对流扩散模型求解运动方程和Oldroyd-B本构方程，针对圆形和椭圆内部边界条件，给出连续界面插值函数，在此基础上，运用光滑界面法将内部边界转换为作用力项施加到演化方程中。首先分析圆柱绕流问题，给出不同材料参数情况下的流场分布和阻力系数计算结果，比较发现与宏观数值模拟结果相吻合。将模型拓展到绕椭圆流动中，分析椭圆形状和材料参数对粘弹性流体绕柱流的影响，发现随着椭圆长轴与短轴比值的增加和维森伯格数的增加，阻力系数逐渐下降，并且长短轴比对迭代收敛有较大影响。     

垂直上升管内油—气—水三相弹状流压力降的研究

给出一种垂直上升油－气－水三相弹状流压力降的计算模型。该模型考虑弹状流中Ｔａｙｌｏｒ气泡周围下降液膜的变化历程。通过油－气－水弹状流的实验研究发现，该模型的数值模拟结果与低压工况下的实验值符合得较好。本模型是计算垂直油－气－水三相弹状流中液相的连续相为水相时的压力降的有效方法。     

φ205超高压釜体螺纹部分热接触有限元分析

用有限元混合法对射孔弹模拟装置中φ２０５超高压釜体螺纹部分的热弹性接触应力进行了综合分析，所用的数学模型把确定接触力的大量迭代运算凝缩到少数螺牙接触点上进行，避免了多次形成和求解原方程，大大提高了运算速度；引入伪弹簧解决了总刚奇异性问题。     

φ205超高压釜体螺纹部分热接触有限元分析

用有限元混合法对射孔弹模拟装置中φ２０５超高压釜体螺纹部分的热弹性接触应力进行了综合分析，所用的数学模型把确定接触力的大量迭代运算凝缩到少数螺牙接触点上进行，避免了多次形成和求解原方程，大大提高了运算速度，引入伪弹簧解决了总刚奇异性问题。     

基于设计空间调整的结构拓扑优化方法

提出了一种基于设计空间调整的结构拓扑优化方法,以求解有限元网格规模大的问题,且获得0/1拓扑解. 首先, 借鉴有理分式材料模型, 建立了材料刚度性质与拓扑变量的关系. 为了解决分析计算量大和需要解释得到的材料分布等问题,给出了一种不影响数学规划求解算法收敛特性的设计空间调整手段. 其次,当优化迭代求解接近结构最佳拓扑邻域后,采用了加速收敛求解的策略,并给出了一种加速收敛的启发式算法. 然后, 结合基于倒设计变量的位移函数的非完整二阶近似式,建立了一种基于设计空间调整的结构拓扑优化算法. 该方法能获得较好0--1分布特征的优化拓扑,能较好地处理多载荷和多约束的结构拓扑优化问题. 给出的算例表明通过结构分析模型规模的减小和传统的位移迭代求解法的采用,方法效率明显提高. 算例验证了该方法的正确性和有效性.      

结构可靠指标求解的一种新的迭代方法

由国际安全度联合委员会推荐的一次二阶距可靠度方法是国际工程界公认的一种较好的可靠度分析方法（简称ＪＣ法）。这种方法具有迭代格式简单，收敛快的特点，但许多实际计算表明，这种方法可适用于功能函数非线性程度不变的情况，否则不收敛，本文提出一种新的迭代方法，这种方法迭代格式与ＪＣ法同样简单，一般情况下均能保证收敛。     

变厚度载流旋转壳体的热磁弹性分析

针对磁场环境中受机械载荷作用变厚度载流旋转壳体热磁弹性问题，基于薄壳的几何方程、物理方程、运动方程、电动力学方程，建立其热磁弹性基本方程．应用线性化方法，得到了关于热磁弹性问题的线性迭代方程，转化为包括8个基本未知量的标准型方程组．关于载流旋转壳体，给出Lorentz力表达式，导出了温度场及温度场积分特征值．分析了变厚度锥形旋转壳体应力、温度及变形随外加电磁参量的变化规律，并通过数值模拟对理论分析进行了验证．研究结果可为载流旋转壳体热磁弹性问题研究提供理论参考．     

一种求算结构可靠度指标的新方法

为解决结构功能函数在结构设计点附近非线性程度较高时,一次可靠度计算方法(如JC法)不收敛的问题,提出一种新的可靠度计算方法.该算法根据结构可靠指标的几何意义,先将迭代点靠近极限状态面,再通过在极限状态面上搜索下一迭代点,逐渐减小极限状态面上的点到标准正态空间中坐标原点的距离,从而达到收敛的目的.与其他方法相比,该方法不用选取算法参数,在计算效率和鲁棒性方面都有较好的优势,尤其在极限状态函数非线性强的情况下优势更为明显.     

有限厚度、大压入的球压试验材料力学响应

针对磁场环境中受机械载荷作用变厚度载流旋转壳体热磁弹性问题，基于薄壳的几何方程、物理方程、运动方程、电动力学方程，建立其热磁弹性基本方程．应用线性化方法，得到了关于热磁弹性问题的线性迭代方程，转化为包括8个基本未知量的标准型方程组．关于载流旋转壳体，给出Lorentz力表达式，导出了温度场及温度场积分特征值．分析了变厚度锥形旋转壳体应力、温度及变形随外加电磁参量的变化规律，并通过数值模拟对理论分析进行了验证．研究结果可为载流旋转壳体热磁弹性问题研究提供理论参考．     

部分输入未知条件下结构动力复合反演的分解算法

在充分利用部分输入已确知而部分输入未知的激励特性的基础上，提出了结构动力复合反演的分解算法，该算法从源头上消除了迭代过程中参数识别与荷载反演的相互影响，降低了问题的计算规模。对于线性参数系统，该算法不经过任何迭代计算即可一次性完成结构参数识别及荷载反演。将其与松弛法结合，可解决非线性参数系统的识别问题，与文献[4]的方法比较，其收敛速度有显著提高。