用MATLAB编程，所有计算都是计入热弹性变形（热弹流润滑），计算值全部在命令行窗口输出。

1、（1）轴承参数采用表4.1的参数，载荷4kN、转速2000r/min，计算计入热弹性变形周向角度分别为30、60、90、120、150、180、210、240、170、300、330、360°，轴向中央截面的压力值（MPa）、油膜温度值（℃）；

  （2）转速2000r/min，计算载荷分别为1、2、3、4、5、6、7、8、9、10kN对应的轴颈偏心率；

2、轴承参数采用表4.2的参数，载荷恒定为4kN，轴承转速从20加速到1600、3000、4500r/min三种工况，时间取0-900s，分别计算三种工况的最大轴承温度， 各计算50个数值（均匀间隔下）。

3、轴承参数采用表4.2的参数，载荷5kN，轴颈初始转速2000r/min，升速分别为1000、2000、3000、4000、5000、6000r/min/s，时间为0-2s（可以根据实际情况选取合适的时间间隔），计算每种工况的轴颈偏心率值、最大轴承温度值、轴颈摩擦力的数值、端泄流量数值，各计算50个数值（均匀时间间隔下），计算值全部保留三位小数。

4、轴承参数采用表4.2的参数，载荷5kN，轴颈初始转速分别为1000、2000、3000、4000、5000、6000r/min，升速为2000r/min/s，时间为0-2s，计算每种工况的轴颈偏心率值、最大轴承温度值、轴颈摩擦力的数值、端泄流量数值，各计算50个数值（均匀间隔下）。

5、轴承参数采用表4.2的参数，轴颈初始转速2000r/min，时间0-1s，初始载荷5kN，升速2000r/min/s，阶跃载荷分别为0、3、6、9、12、15kN，计算每种工况的轴颈偏心率值、最大轴承温度值、轴颈摩擦力的数值、端泄流量数值，各计算50个数值（均匀间隔下）。

6、轴承参数采用表4.2的参数，恒定初始转速2000r/min，升速2000r/min/s，外载荷5kN，轴颈外径面不平衡质量分别为10、20、30、40、50、60g，计算每种工况的轴颈偏心率值、最大轴承温度值、轴颈摩擦力的数值、端泄流量数值，各计算50个数值（均匀间隔下）。