本项目致力于构建一个基于MATLAB环境的高效数值计算平台，专门用于模拟和分析线接触工况下的常温弹性流体动力润滑（EHL）行为。核心算法采用全系统耦合的Newton-Raphson（牛顿-拉夫逊）迭代法，通过构建包含雷诺方程、弹性变形积分方程、载荷平衡方程以及润滑油物性方程（Roelands粘压关系和Dowson-Higginson密压关系）在内的非线性方程组，并计算其雅可比矩阵（Jacobian Matrix）进行线性化求解，从而实现了压力场、膜厚场和弹性变形场的同步更新。该方法显著优于传统的直接迭代法或逆解法，特别是在高赫兹接触压力（High Hertzian Pressure）和重载工况下具有极佳的收敛速度和数值稳定性。程序首先根据赫兹接触理论生成合理的压力和膜厚初值，利用有限差分法对控制方程进行离散化，通过多重网格积分或直接求和法计算弹性变形，并在迭代过程中精确捕捉润滑区内的特征物理现象，包括出口处的二次压力峰（Petrusevich峰）以及膜厚颈缩现象。数据后处理模块能够将计算得到的无量纲参数还原为实际物理量，自动生成高分辨率的压力分布曲线图和沿流动方向的膜厚轮廓图，并计算输出中心膜厚、最小膜厚及其位置坐标，旨在为齿轮、滚动轴承等机械摩擦副的润滑性能评估与优化设计提供可靠的理论工具。
