本帖最后由 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 上傳

總結

以上是生活随笔為你收集整理的matlab简易编程,MATLAB简单编程的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。