本帖最后由 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; %沿y方向上区间个数
fai=60*pi/180; %滑靴转动角度，(如果可以，fai为向量)单位： rad
R0=15.4*10^(-3); %滑靴外径 单位： m
r0=10.1*10^(-3); %支承室半径 单位 ： m
TT=100; %时间
dt=TT/(2*nx); %时间间隔
deltax=2*pi/nx; %x方向上区间长度 (极角)，单位： rad
deltay=(R0-r0)/ny; %y方向上区间长度，半径 单位： m
ms=50*10^(-3); %滑靴质量， 单位： kg
mh=56*10^(-3); %柱塞质量 单位： kg
lg=8*10^(-3); %滑靴质心与球头中心距离, 单位： m
ls=10*10^(-3); %滑靴底面到球头的距离， 单位： m
RP=12.4*10^(-3); %柱塞半径 单位： m
RFF=10.5; %柱塞孔的等效液阻,单位：Pa*s/mm^3(这个数值没用上额)
rou=RCP*sqrt(1+(tan(beita))^2*(cos(fai))^2); %公式(数组) m
vs=wg*rou*2*pi/60; %滑靴速度(数组) m/s
f0=0;
MOX=0;
MOY=0; %三个函数初始化
P=0.75*pr*ones(nx+1,ny+1); %压力矩阵，初始值 Pa
P2=ones(nx+1,ny+1); %压力矩阵，初始值
P3=ones(nx+1,ny+1); %压力矩阵，初始值
H=ones(nx+1,ny+1); %油膜迭代初始矩阵 油膜厚度
HL=ones(nx+1,ny+1); %油膜厚度对时间求导矩阵
W=ones(nx+1,ny+1); %系数初始矩阵，西
E=ones(nx+1,ny+1); %系数初始矩阵，东
S=ones(nx+1,ny+1); %系数初始矩阵，南
N=ones(nx+1,ny+1); %系数初始矩阵，北
G=ones(nx+1,ny+1); %系数初始矩阵，东西南北四个系数之和
B=ones(nx+1,ny+1); %原雷诺方程右边常系数，，矩阵初始化
FO1=pi*RP^2*pc; %柱塞腔作用力 N=m^2*Pa
FO2=(ms+mh)*(wg/60)^2*RCP*cos(fai)*tan(beita); %惯性力 (wg/60)为角速度 N
MX=ms*(wg/60)^2*lg*RCP*cos(beita); %滑靴油膜的离心力矩 N*m
con=0; %第二次迭代次数初值
count=0; %压力迭代次数初始值(第三层)
coo=0; %第一层迭代次数初始值
h0=10^(-3)*[0.008 0.01 0.009]' ; %油膜初值 单位： m
hl=10^(-3)*[0.001 0.002 0.003]'; %油膜变化率初值
2016-11-23 21:56 上传
2016-11-23 21:56 上传
