基于Matlab的《液体动静压轴承》回油槽径向静压轴承图谱程序
打开MATLAB就闻到机油味是怎么回事？最近折腾了个有意思的玩意——基于《液体动静压轴承》教材
搞的径向静压轴承仿真程序。这可不是玩具代码，实测能生成正经的工程图谱，连油槽回油这种细节都考
虑进去了。
先看核心参数设置部分，这里藏着轴承的DNA：
```matlab
R = 0.05; % 轴承半径(m)
L = 0.1; % 轴承长度
mu = 0.018; % 油液粘度(Pa·s)
h0 = 1e-5; % 初始油膜厚度
P0 = 1e5; % 供油压力
```
粘度参数mu特别容易踩坑，实测发现超过0.02时温升曲线会像火箭发射。建议先用教材推荐值，跑
通后再慢慢调参。
油膜压力计算这块用了双重循环，老工程师的暴力美学：
```matlab
function P = oilFilmPressure(theta, h)
global R L mu omega
% 压力梯度方程求解
dPdtheta = (6*mu*omega*R)./(h.^3) .* (h - h0);
P = cumtrapz(theta, dPdtheta) + P0;
% 压力修正项处理
P(P < 0) = 0; % 杜绝负压区
end
```
cumtrapz积分函数是MATLAB的老传统了，比手工写梯形法省事十倍。不过要注意theta的采样间隔，
太稀疏会导致积分误差让曲线长毛刺。
主程序骨架长这样：
```matlab
theta = linspace(0, 2*pi, 360); % 周向角度细分
