基于Matlab的《液体动静压轴承》回油槽径向静压轴承图谱程序

基于matlab的根据《液体动静压轴承》编写的有回油槽径向静压轴承的可显示承载能力、压强、刚度及温升等图谱.程序已调通可直接运行。打开MATLAB就闻到机油味是怎么回事最近折腾了个有意思的玩意——基于《液体动静压轴承》教材搞的径向静压轴承仿真程序。这可不是玩具代码实测能生成正经的工程图谱连油槽回油这种细节都考虑进去了。先看核心参数设置部分这里藏着轴承的DNAR 0.05; % 轴承半径(m) L 0.1; % 轴承长度 mu 0.018; % 油液粘度(Pa·s) h0 1e-5; % 初始油膜厚度 P0 1e5; % 供油压力粘度参数mu特别容易踩坑实测发现超过0.02时温升曲线会像火箭发射。建议先用教材推荐值跑通后再慢慢调参。油膜压力计算这块用了双重循环老工程师的暴力美学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; % 杜绝负压区 endcumtrapz积分函数是MATLAB的老传统了比手工写梯形法省事十倍。不过要注意theta的采样间隔太稀疏会导致积分误差让曲线长毛刺。主程序骨架长这样theta linspace(0, 2*pi, 360); % 周向角度细分 h h0 * (1 0.5*cos(theta)); % 油膜厚度变化模型 [W, K, deltaT] deal(zeros(1,50)); % 预分配数组 for i 1:50 omega 100 200*(i-1)/49; % 转速渐变 % 调用求解器计算压力场 P oilFilmPressure(theta, h); % 承载能力积分计算 W(i) trapz(theta, P.*L.*R.*cos(theta)); % 刚度计算需要压力微分 K(i) abs(trapz(theta, gradient(P,theta).*L.*R)); % 温升模型 deltaT(i) calcTemperatureRise(P, h); end转速循环里藏了个小技巧用(i-1)/49实现0到1的线性过渡比直接写linspace省个变量。gradient函数求压力梯度时默认会用中心差分比前向差分稳定得多。基于matlab的根据《液体动静压轴承》编写的有回油槽径向静压轴承的可显示承载能力、压强、刚度及温升等图谱.程序已调通可直接运行。可视化部分用了subplot玩拼图figure(Color,w,Position,[100 100 1200 900]) subplot(2,2,1) plot(omega, W/1e3, LineWidth,2) xlabel(转速 (rad/s)); ylabel(承载力 (kN)); grid on; title(承载能力曲线) % 其他子图类似...这里把窗口预设为白色背景指定尺寸避免默认的灰底在论文配图时被导师怼。LineWidth调到2让曲线在PPT演示时不会被后排观众吐槽太细。跑起来之后能看到四个指标联动画图转速从100rad/s扫到300rad/s的过程。特别说下刚度曲线——在临界转速附近会出现拐点这个特征抓准了仿真就算成了。程序包里还塞了个animation子函数能把压力分布做成动态图给老板汇报时特别省心。想要自己魔改的话注意几个暗桩油膜厚度模型h那行可以改成更复杂的函数不过别改得太放飞自我小心积分爆炸。温升计算函数里埋了个对流换热系数不同润滑油记得改这个参数不然温升预测会飘到外太空。源码已打包成MATLAB Live Script点开就能看交互式图表。下次考虑加个GUI界面用滑块调参应该更带感。有做轴承设计的同行欢迎来交流调参心得毕竟仿真和实测的对齐是永恒的战斗啊。