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

Improvement of Numerical Algorithm for Elastohydrodynamic Lubrication Based on Mixed Programming of MatLab and C/C++

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

Based on the existing numerical calculation method for isothermal elliptical contact elastohydrodynamic lubrication,the influence of the variation of film thickness constantΔH on the relative accuracy of pressure iteration was analyzed by using the Newton finite difference method.It is found that when the initial value ofΔH is 10%of the minimum film thickness,within a certain range,the smaller the change period and the larger the reduction range ofΔH,the smaller the relative precision.After setting a reasonable variation ofΔH,when judging that the convergence accuracy of pressure tends to be stable,the partial influence of iteration method can be neutralized by enlargingΔH at the minimum of relative accuracy.It can satisfy the convergence condition at the lowest point to complete the numerical calculation,avoid falling into the dead cycle and improve the calculation efficiency.