 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1469 | 回復: 7 | |||
LQgraduate銅蟲 (初入文壇)
|
[求助]
求解軸承潤滑性能時，如果承載力不收斂，應該如何調(diào)整油膜厚度？ 已有1人參與
|
|
問題如下（希望做過這方面工作的好心人幫忙解決一下）： 計算流體動壓推力軸承或者徑向軸承的潤滑性能時，首先需要利用迭代法使壓力收斂，從而求解雷諾方程；然后要使軸承的承載力也收斂，如果不收斂，需要修正最初假設的最小和最大油膜厚度，直到收斂為止。問題來了，就是如果承載力不收斂，應該如何調(diào)整油膜厚度？有什么方法嗎？還是純粹靠經(jīng)驗？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 機械科研經(jīng)驗專輯 |
» 猜你喜歡
0858動力工程290求調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑，認真踏實，吃苦耐勞，有科研潛力
已經(jīng)有0人回復
-
機械工程論文潤色/翻譯怎么收費?
已經(jīng)有105人回復
-
0858能源動力求B區(qū)調(diào)劑，動力工程也可以！
已經(jīng)有0人回復
-
關于Chinese Journal of Aeronautics期刊收不收版面費的投稿求助貼
已經(jīng)有0人回復
-
招收0805（材料）調(diào)劑
已經(jīng)有0人回復
-
調(diào)劑 動力工程322 一志愿中石油華東
已經(jīng)有2人回復
-
招收0805（材料）調(diào)劑
已經(jīng)有3人回復
-
求調(diào)劑，總分391
已經(jīng)有3人回復
-
EBSD極圖和反極圖怎么看
已經(jīng)有0人回復
Tennisgraf
銀蟲 (正式寫手)
- 應助: 5 (幼兒園)
- 金幣: 13.4
- 散金: 162
- 紅花: 2
- 帖子: 850
- 在線: 206.7小時
- 蟲號: 691801
- 注冊: 2009-01-12
- 專業(yè): 機械摩擦學與表面技術
2樓2013-01-10 14:49:31
LQgraduate
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時
- 蟲號: 1549417
- 注冊: 2011-12-23
- 性別: GG
- 專業(yè): 機械摩擦學與表面技術
|
首先，謝謝你的回復。但是，問題在于，你的意思是說要靠經(jīng)驗調(diào)整了，反復嘗試對吧，知道承載力收斂為止。那么有沒有什么方法，比如某個算法，可以實現(xiàn)編程然后讓油膜厚度自動調(diào)整直到承載力收斂為止��？我在文獻上看到了一種方法，說是用二分法可以，但我試了一下，做出來啊!這樣吧，我把程序貼到上邊來，你看幫我看一下問題出在哪里吧？我加了注釋了，如果有看不懂的地方，再交流。以下是我用matlab編寫的五點差分法離散雷諾方程用SOR(超松弛迭代法)求解壓力分布的源程序： clear all r1=31e-3;%推力瓦內(nèi)徑 r2=71e-3;%推力瓦外徑 h1=15e-6;%假設的初始最小水膜厚度 h2=30e-6;%假設的初始最大水膜厚度 theta1=50*pi/180;%推力瓦的夾角 beita=0.001*pi/180;%推力瓦周向傾角 k=6;%推力瓦塊數(shù) F1=3000;%軸向推力 n1=1000;%軸承轉速 omiga=2*pi*n1/60;%推力軸承角速度 miu=0.001;%水的粘度 l=((r1+r2)/2)*theta1;%軸瓦長度 P0=1.0064e+06;%壓力量綱系數(shù) nx=40;%周向劃分的網(wǎng)格數(shù) ny=40;%徑向劃分的網(wǎng)格數(shù) deltax=(50/nx)*(pi/180);%周向網(wǎng)格的單位長度 deltay=(r2-r1)/(r1*ny);%徑向網(wǎng)格的單位長度 r=zeros(nx+1,ny+1);%初始化無量綱半徑矩陣 theta=zeros(nx+1,ny+1);%初始化無量綱角度矩陣 H=zeros(nx+1,ny+1);%初始化無量綱水膜厚度矩陣 while 1 for i=1:nx+1 for j=1:ny+1 r(i,j)=1+(j-1)*deltay;%無量綱半徑矩陣 theta(i,j)=(i-1)*deltax;%無量綱角度矩陣 H(i,j)=h2/h1-((h2-h1)/h1)*(theta(i,j)/theta1);%可傾瓦水膜厚度無量綱表達式 end end h=h1*H;%有量綱潤滑膜厚度矩陣 P=zeros(nx+1,ny+1);%初始化無量綱節(jié)點壓力 A=zeros(nx+1,ny+1);%初始化系數(shù)矩陣A B=zeros(nx+1,ny+1);%初始化系數(shù)矩陣B C=zeros(nx+1,ny+1);%初始化系數(shù)矩陣C D=zeros(nx+1,ny+1);%初始化系數(shù)矩陣D E=zeros(nx+1,ny+1);%初始化系數(shù)矩陣E F=zeros(nx+1,ny+1);%初始化系數(shù)矩陣F P2=zeros(nx+1,ny+1);%初始化無量綱壓力矩陣P2 P3=zeros(nx+1,ny+1);%初始化無量綱壓力矩陣P3 S=0;%初始化誤差表達式的分子 T=0;%初始化誤差表達式的分母 wucha=1;%初始化誤差 error=1e-3;%壓力迭代終止精度 error2=0.01;%承載力迭代終止精度 count=0;%初始化迭代次數(shù) w=1.75;%松弛因子 while(wucha>=error) PP=P; for i=2:nx for j=2:ny P2(i,j)=P(i,j); A(i,j)=((H(i+1,j)+H(i,j))/2)/r(i,j); B(i,j)=((H(i-1,j)+H(i,j))/2)/r(i,j); C(i,j)=((deltax/deltay)^2)*((r(i,j+1)+r(i,j))/2)*((H(i,j+1)+H(i,j))/2)^3; D(i,j)=((deltax/deltay)^2)*((r(i,j-1)+r(i,j))/2)*((H(i,j-1)+H(i,j))/2)^3; E(i,j)=A(i,j)+B(i,j)+C(i,j)+D(i,j); F(i,j)=6*deltax*r(i,j)*((H(i+1,j)-H(i-1,j))/2); P3(i,j)=(A(i,j)*P(i+1,j)+B(i,j)*P(i-1,j)+C(i,j)*P(i,j+1)+D(i,j)*P(i,j-1)-F(i,j))/E(i,j); P(i,j)=(1-w)*P2(i,j)+w*P3(i,j); if P(i,j)<0 P(i,j)=0; end end end for i=2:nx for j=2:ny S=S+abs(P(i,j)-PP(i,j)); T=T+abs(P(i,j)); end end wucha=S/T; count=count+1; end load=0;%承載力初始值 for i=1:nx for j=1:ny ds=pi*(r(i,j+1)^2-r(i,j)^2)*(deltax/(2*pi)); load=load+k*P0*P(i,j)*ds*(r1^2); end end tt=abs(load-F1)/F1; if tt hmax=max(max(h)); break; elseif load>F1 h1=(h1+h2)/2; h2=h1+l*tan(beita); else h2=(h1+h2)/2; h1=h2-l*tan(beita); end P0=omiga*miu*r1^2/(h1^2); end figure(1); [x,y]=meshgrid((0:deltax/(pi/180):50),(0:deltay*r1:r2-r1)); mesh(x,y,h); xlabel('周向/°'); ylabel('徑向r/m'); zlabel('h/m'); figure(2); [x,y]=meshgrid((0:deltay*r1:r2-r1),(0:deltax/(pi/180):50)); mesh(x,y,P); xlabel('徑向r/m'); ylabel('周向/°'); zlabel('P'); |
3樓2013-01-10 15:23:23
LQgraduate
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時
- 蟲號: 1549417
- 注冊: 2011-12-23
- 性別: GG
- 專業(yè): 機械摩擦學與表面技術
4樓2013-01-10 15:31:00
LQgraduate
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時
- 蟲號: 1549417
- 注冊: 2011-12-23
- 性別: GG
- 專業(yè): 機械摩擦學與表面技術
5樓2013-01-10 17:06:02
wshuagfnui
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 117.5
- 帖子: 1
- 在線: 37.4小時
- 蟲號: 1200847
- 注冊: 2011-02-10
- 專業(yè): 機械摩擦學與表面技術
6樓2013-01-26 12:27:27
LQgraduate
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時
- 蟲號: 1549417
- 注冊: 2011-12-23
- 性別: GG
- 專業(yè): 機械摩擦學與表面技術
7樓2013-01-28 14:53:55
yangqijiang
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 11.8
- 紅花: 1
- 帖子: 64
- 在線: 77小時
- 蟲號: 1519679
- 注冊: 2011-12-01
- 專業(yè): 生物無機化學
8樓2014-03-15 21:54:35
相關版塊跳轉
我要訂閱樓主 LQgraduate 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +9 | zjl的號 2026-03-16 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
|
[考研] 一志愿211化學學碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|