基于MatLab与C/C++混合编程的弹流润滑数值算法改进

基于等温椭圆接触弹流润滑现有的数值计算方法,选取牛顿有限差分法进行压力迭代,分析不同膜厚常数变化量ΔH对压力迭代相对精度的影响。分析表明,在ΔH初值取最小膜厚的10%时,在一定范围内,ΔH变化周期越小、减小幅度越大,相对精度整体越小;合理设置ΔH变化形式,在判断出压力收敛精度趋于稳定时,辅以相对精度极小点处进行ΔH变大处理,可中和部分迭代方法的影响,能在最低点处满足收敛条件而完成数值计算,避免陷入死循环而提高计算效率。