基于MatLab与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.
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.
引文
[1]刘俊明,王文中,赵自强.基于弹流润滑模型的润滑油黏度确定方法[J].润滑与密封,2018,43(7):7-12.LIU J M,WANG W Z,ZHAO Z Q.Method for determining lubricant viscosity based on elastohydrodynamic lubrication theory[J].Lubrication Engineering,2018,43(7):7-12.
[2]曹蔚,王栋,王宁,等.高精密弹流润滑球轴承支承主轴系统动特性分析[J].重庆理工大学学报(自然科学),2017(2):46-51.CAO W,WANG D,WANG N,et al.Dynamic characteristic analysis of spindle system supported by high precision EHL-ball bearings[J].Journal of Chongqing University of Technology(Natural Science),2017(2):46-51.
[3]武丹丹,刘晓玲,龙涛.基于Elrod算法的线接触乏油热弹流润滑分析[J].润滑与密封,2018,43(10):51-57.WU D D,LIU X L,LONG T.Thermal analysis for starved elastohydrodynamic lubrication line contacts based on Elrod algorithm[J].Lubrication Engineering,2018,43(10):51-57.
[4]BUJURKE N M,KANTLI M H,SHETTAR B M.Wavelet preconditioned Newton-Krylov method for elastohydrodynamic lubrication of line contact problems[J].Applied Mathematical Modelling,2017,46:285-298.
[5]徐彩红,王优强,张同钢,等.载荷时变对齿轮齿条弹流润滑的影响[J].表面技术,2017,46(6):77-83.XU C H,WANG Y Q,ZHANG T G,et al.Influence of load time-varying on elastohydrodynamic lubrication of gear rack[J].Surface Technology,2017,46(6):77-83.
[6]韩兴,李昌,赵焱,等.椭圆点接触弹流润滑数值计算与影响因素分析[J].机械设计,2018,35(S1):58-64.HAN X,LI C,ZHAO Y,et al.Numerical calculation and influence factors analysis of elastohydrodynamic lubrication of elliptic point contacts[J].Journal of Machine Design,2018,35(S1):58-64.
[7]温诗铸,黄平.摩擦学原理[M].2版.北京:清华大学出版社,2002.
[8]刘维.精通MATLAB与C/C++混合程序设计[M].北京:北京航空航天大学出版社,2014.
[9]潘克家,汤井田,郑洲顺.MatLab与Fortran混合编程之DLL实现方法[J].计算机工程与应用,2011,47(31):1-3.PAN K J,TANG J T,ZHENG Z S.Mixed language programming in MatLab and Fortran via dynamic linklibrary[J].Computer Engineering and Applications,2011,47(31):1-3.
[10]陈庆远,武富春,程秀芹,等.MatLab与C混合编程在火控系统仿真中的应用[J].火力与指挥控制,2014,39(S1):156-159.CHEN Q Y,WU F C,CHENG X Q,et al.Application research on hybrid programming of MatLab and C language in the fire control system simulation[J].Fire Control&Command Control,2014,39(S1):156-159.
[11]黄平.弹性流体动压润滑数值计算方法[M].北京:清华大学出版社,2013.
[2]曹蔚,王栋,王宁,等.高精密弹流润滑球轴承支承主轴系统动特性分析[J].重庆理工大学学报(自然科学),2017(2):46-51.CAO W,WANG D,WANG N,et al.Dynamic characteristic analysis of spindle system supported by high precision EHL-ball bearings[J].Journal of Chongqing University of Technology(Natural Science),2017(2):46-51.
[3]武丹丹,刘晓玲,龙涛.基于Elrod算法的线接触乏油热弹流润滑分析[J].润滑与密封,2018,43(10):51-57.WU D D,LIU X L,LONG T.Thermal analysis for starved elastohydrodynamic lubrication line contacts based on Elrod algorithm[J].Lubrication Engineering,2018,43(10):51-57.
[4]BUJURKE N M,KANTLI M H,SHETTAR B M.Wavelet preconditioned Newton-Krylov method for elastohydrodynamic lubrication of line contact problems[J].Applied Mathematical Modelling,2017,46:285-298.
[5]徐彩红,王优强,张同钢,等.载荷时变对齿轮齿条弹流润滑的影响[J].表面技术,2017,46(6):77-83.XU C H,WANG Y Q,ZHANG T G,et al.Influence of load time-varying on elastohydrodynamic lubrication of gear rack[J].Surface Technology,2017,46(6):77-83.
[6]韩兴,李昌,赵焱,等.椭圆点接触弹流润滑数值计算与影响因素分析[J].机械设计,2018,35(S1):58-64.HAN X,LI C,ZHAO Y,et al.Numerical calculation and influence factors analysis of elastohydrodynamic lubrication of elliptic point contacts[J].Journal of Machine Design,2018,35(S1):58-64.
[7]温诗铸,黄平.摩擦学原理[M].2版.北京:清华大学出版社,2002.
[8]刘维.精通MATLAB与C/C++混合程序设计[M].北京:北京航空航天大学出版社,2014.
[9]潘克家,汤井田,郑洲顺.MatLab与Fortran混合编程之DLL实现方法[J].计算机工程与应用,2011,47(31):1-3.PAN K J,TANG J T,ZHENG Z S.Mixed language programming in MatLab and Fortran via dynamic linklibrary[J].Computer Engineering and Applications,2011,47(31):1-3.
[10]陈庆远,武富春,程秀芹,等.MatLab与C混合编程在火控系统仿真中的应用[J].火力与指挥控制,2014,39(S1):156-159.CHEN Q Y,WU F C,CHENG X Q,et al.Application research on hybrid programming of MatLab and C language in the fire control system simulation[J].Fire Control&Command Control,2014,39(S1):156-159.
[11]黄平.弹性流体动压润滑数值计算方法[M].北京:清华大学出版社,2013.