无条件稳定数值积分方法在拟动力实验中的应用研究

用无条件稳定的隐式α－方法代替目前的条件稳定的显式方法，通过数值迭代控制方法实现了拟动力实验．在此基础上进一步应用子结构技术，解决了刚度大、自由度较多结构采用传统的条件稳定数值积分方法进行拟动力实验的问题．应用不同数值积分方法进行的三层模型结构的实验结果表明，无条件稳定的数值积分方法的实验精度优于条件稳定的数值积分方法的精度．文中还就实验误差的控制补偿问题进行了讨论，提出了采用瞬时割线刚度进行补偿的方法     

瞬态空化泡演变过程的数值模拟

采用边界积分方程方法，对无粘流体中三个空化泡以及自由面附近二个空化泡相互作用的演变过程进行了数值模拟。计算中边界用二阶有限元离散，影响系数矩阵非对角线元素用六点高斯数值积分方法计算，对第一类、第二类完全椭圆积分用高次多项式近似，对计算系数矩阵对角线元素中遇到的奇异积分进行了特殊处理。结果表明，在不同的给定参数下，空化泡的溃灭形态各异，柱状射流和环形射流都有可能发生，使空化泡演变成双泡或环形泡。     

多体系统动力学微分／代数方程组数值方法

多体系统动力学微分／代数混合方程组又称Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方程，是近十年来动力学和计算数学领域研究的热点之一．本文介绍这两个领域中引入的传统的数值积分方法与新的理论．     

权值自适应调整Unscented粒子滤波及其在组合导航中的应用

针对粒子滤波存在的重要性密度函数难以选取和粒子退化问题,提出了一种新的权值自适应调整Unscented粒子滤波算法。该算法在Unscented粒子滤波的采样过程中吸收权值自适应调整的优点,考虑最新量测影响,通过欧氏距离和反映量测噪声统计特性的精度因子来自适应的调整粒子对应权值分布,增加有用粒子的权值,降低粒子退化程度,保持粒子多样性。同时Unscented变换提高了滤波精度,使该算法能更好地适用于非线性、非高斯系统模型的计算。将提出的算法应用于GPS/DR组合导航系统进行仿真验证,结果表明,提出的权值自适应调整Unscented粒子滤波算法得到的东向定位误差控制在±5.5 m附近,北向定位误差则在±5.2 m附近,滤波性能明显优于扩展卡尔曼滤波和Unscented粒子滤波,能提高GPS/DR组合导航系统解算精度。     

数字图像相关方法中匹配及过匹配形函数的误差分析

基于图像子区的数字图像相关方法需采用合适的形函数来近似目标图像子区的真实变形。由于实际测量时目标子区的局部变形往往是未知的,实际采用的不同阶次(零阶、一阶和二阶)的形函数不可避免地产生误匹配(欠匹配和过匹配)问题,从而引入位移测量的系统或随机误差。尽管由欠匹配形函数引起的系统误差已被充分认识,由过匹配形函数引起的位移误差仍缺少理论解释。本文首先推导出采用一阶和二阶形函数的数字图像相关方法的随机误差理论公式,随后采用一系列数值实验验证了理论公式的准确性。结果显示:过匹配形函数不会引入额外的系统误差,但会增加随机误差,且二阶形函数的随机误差是一阶形函数的二倍。考虑到由欠匹配一阶形函数引入的系统误差往往远大于过匹配二阶形函数的随机误差,因此在未能确知变形的情况下,推荐使用二阶形函数。     

双颗粒沉降的动力学分析

采用格子Boltzmann-虚拟区域方法对雷诺数范围为50≤Re≤200的双颗粒自由沉降进行了直接数值模拟。首先,研究发现双颗粒最终沉降的位置分别在通道的1/4和3/4位置附近,且与颗粒的初始间距、雷诺数以及通道宽度无关。其次,重点研究了颗粒沉降过程中所受的侧向力(与沉降方向垂直),首次揭示了侧向力的振动频率与雷诺数呈二次关系,且单颗粒的结果始终小于双颗粒的结果;研究还发现侧向力的振动频率与通道宽度近似呈幂律函数关系,且幂律指数与雷诺数有关,雷诺数越大,幂律指数的绝对值越小。最后还研究了雷诺数及通道宽度对侧向力振动振幅的影响。     

刃型位错和圆弧形刚性线夹杂的干涉效应

位错和夹杂的干涉效应对于理解材料的强化和韧化机理具有十分重要的意义。文中研究了晶体材料中刃型位错和多条共圆弧刚性线夹杂的干涉作用。利用Riemann—Schwarz反照原理和复势函数的奇性主部分析技术，得到了问题的一般解答；对于只含一条刚性线夹杂的情况，给出了复势函数的封闭形式解。由Peach-Koehler公式求出了作用在刃型位错上的位错力，并讨论了圆弧形刚性线夹杂对位错力的影响规律，发现弧形刚性线对刃型位错有很强的排斥作用。本文解答不但可作为格林函数获得任意分布位错的相应解答，而且可以用于研究刚性线夹杂和任意形状裂纹的干涉效应问题。     

压电材料平面裂纹尖端场的杂交应力有限元分析

基于复势理论和杂交变分原理建立了一种适用于力电耦合分析的杂交应力有限元模 型. 给出了建立刚度矩阵的主要公式和推导过程，单元内的位移场和应力场采用满足平 衡方程的复变函数级数解，假设的复变函数级数解事先精确满足裂纹的无应力和电位移法向 分量为零的条件，单元外边界的位移场假设按抛物线变化, 单元的刚度矩阵采用Gauss积分的方法得出. 通过对力电耦合裂尖场的数值计算验证了程序 的正确性和单元的有效性，同时也用所得结果校验了理论解.     

Cr12MoV钢表面磁控溅射Ti/TiN涂层的摩擦磨损性能研究

采用非平衡磁控溅射方法在Cr12MoV钢表面制备了厚度约为3 μm的Ti/TiN涂层,测定了涂层的显微硬度,并通过划痕试验和摩擦磨损试验考察了涂层同基体的结合强度及其摩擦磨损性能.结果表明:Ti/TiN涂层能够显著提高Cr12MoV钢的表面硬度和承载能力;涂层同Cr12MoV钢基体的结合强度较高,划痕临界载荷高于60 N;与此同时,磁控溅射Ti/TiN涂层可以显著改善Cr12MoV钢的耐磨性能.这是由于磁控溅射Ti/TiN涂层硬度高且与Cr12MoV钢基体的结合强度较高所致.     

磁头/磁盘滑动接触下磁盘温度及热退磁临界条件的研究

采用二维轴对称有限元模型计算磁头/磁盘滑动接触下,铝质磁盘的稳态温度和热应力场以及热退磁临界条件.结果表明:磁盘温度在极短时间内升至摩擦稳态值,然后缓慢线性升高到最终稳态值;经过充分热传导和热交换后磁盘的温度梯度较小,此时磁层内的热应力集中分布于磁盘固定端边缘附近;磁盘的稳态温度和热应力均随速度增大而增大,且载荷越大其值增大越快;热应力小于1.2 GPa时所对应的速度和载荷为安全工况;温升大于373 K时所对应的工况将导致磁盘退磁.