本帖最后由 wanggh 于 2016-12-5 14:47 编辑
用数值差分、SOR迭代法求雷诺方程和用牛顿迭代法求解轴向柱塞泵滑靴副压力场的算法,
%油膜厚度场、压力场迭代 MATLAB只认弧度制,不认角度制
clear all; % (60rpm=1rad/s)
wg=1000; %1000rpm缸体公转速度; 单位: 转每分钟 rpm
wz=60; %ws为滑靴自转速度 单位: 转每分钟 rpm
RCP=45*10^(-3); %柱塞分布圆半径,单位: m
pc=24000000; %柱塞腔压力,单位:Pa(单位不能变,就是Pa) Pa
aa=0.87; %压降系数,和阻尼孔直径、长度等有关,无量纲
pr=pc/aa; %支承室压力 单位: Pa
NIANDU=0.04; %介质动压粘度, 单位: Pa*S
beita=18*pi/180; %柱塞泵斜盘倾角 单位:rad 弧度制
nx=100; %沿x方向上区间个数
ny=100;