 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1481 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
LQgraduate銅蟲(chóng) (初入文壇)
|
[求助]
求解軸承潤(rùn)滑性能時(shí)，如果承載力不收斂，應(yīng)該如何調(diào)整油膜厚度？ 已有1人參與
|
||
|
問(wèn)題如下（希望做過(guò)這方面工作的好心人幫忙解決一下）： 計(jì)算流體動(dòng)壓推力軸承或者徑向軸承的潤(rùn)滑性能時(shí)，首先需要利用迭代法使壓力收斂，從而求解雷諾方程；然后要使軸承的承載力也收斂，如果不收斂，需要修正最初假設(shè)的最小和最大油膜厚度，直到收斂為止。問(wèn)題來(lái)了，就是如果承載力不收斂，應(yīng)該如何調(diào)整油膜厚度？有什么方法嗎？還是純粹靠經(jīng)驗(yàn)？ |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 機(jī)械科研經(jīng)驗(yàn)專(zhuān)輯 |
» 猜你喜歡
英國(guó)巴斯大學(xué) University of Bath 博士后及博士招聘
已經(jīng)有5人回復(fù)
-
求調(diào)劑收留
已經(jīng)有1人回復(fù)
-
工程熱物理與能源利用論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有154人回復(fù)
-
085500機(jī)械專(zhuān)碩 331分求調(diào)劑
已經(jīng)有1人回復(fù)
-
機(jī)械專(zhuān)碩311求調(diào)劑
已經(jīng)有0人回復(fù)
-
學(xué)碩287求調(diào)劑
已經(jīng)有0人回復(fù)
-
學(xué)碩287求調(diào)劑
已經(jīng)有5人回復(fù)
-
面上參與人員重要嗎？找不到參與人
已經(jīng)有5人回復(fù)
LQgraduate
銅蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時(shí)
- 蟲(chóng)號(hào): 1549417
- 注冊(cè): 2011-12-23
- 性別: GG
- 專(zhuān)業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
4樓2013-01-10 15:31:00
Tennisgraf
銀蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 13.4
- 散金: 162
- 紅花: 2
- 帖子: 850
- 在線: 206.7小時(shí)
- 蟲(chóng)號(hào): 691801
- 注冊(cè): 2009-01-12
- 專(zhuān)業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
2樓2013-01-10 14:49:31
LQgraduate
銅蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時(shí)
- 蟲(chóng)號(hào): 1549417
- 注冊(cè): 2011-12-23
- 性別: GG
- 專(zhuān)業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
|
首先，謝謝你的回復(fù)。但是，問(wèn)題在于，你的意思是說(shuō)要靠經(jīng)驗(yàn)調(diào)整了，反復(fù)嘗試對(duì)吧，知道承載力收斂為止。那么有沒(méi)有什么方法，比如某個(gè)算法，可以實(shí)現(xiàn)編程然后讓油膜厚度自動(dòng)調(diào)整直到承載力收斂為止��？我在文獻(xiàn)上看到了一種方法，說(shuō)是用二分法可以，但我試了一下，做出來(lái)啊!這樣吧，我把程序貼到上邊來(lái)，你看幫我看一下問(wèn)題出在哪里吧？我加了注釋了，如果有看不懂的地方，再交流。以下是我用matlab編寫(xiě)的五點(diǎn)差分法離散雷諾方程用SOR(超松弛迭代法)求解壓力分布的源程序： clear all r1=31e-3;%推力瓦內(nèi)徑 r2=71e-3;%推力瓦外徑 h1=15e-6;%假設(shè)的初始最小水膜厚度 h2=30e-6;%假設(shè)的初始最大水膜厚度 theta1=50*pi/180;%推力瓦的夾角 beita=0.001*pi/180;%推力瓦周向傾角 k=6;%推力瓦塊數(shù) F1=3000;%軸向推力 n1=1000;%軸承轉(zhuǎn)速 omiga=2*pi*n1/60;%推力軸承角速度 miu=0.001;%水的粘度 l=((r1+r2)/2)*theta1;%軸瓦長(zhǎng)度 P0=1.0064e+06;%壓力量綱系數(shù) nx=40;%周向劃分的網(wǎng)格數(shù) ny=40;%徑向劃分的網(wǎng)格數(shù) deltax=(50/nx)*(pi/180);%周向網(wǎng)格的單位長(zhǎng)度 deltay=(r2-r1)/(r1*ny);%徑向網(wǎng)格的單位長(zhǎng)度 r=zeros(nx+1,ny+1);%初始化無(wú)量綱半徑矩陣 theta=zeros(nx+1,ny+1);%初始化無(wú)量綱角度矩陣 H=zeros(nx+1,ny+1);%初始化無(wú)量綱水膜厚度矩陣 while 1 for i=1:nx+1 for j=1:ny+1 r(i,j)=1+(j-1)*deltay;%無(wú)量綱半徑矩陣 theta(i,j)=(i-1)*deltax;%無(wú)量綱角度矩陣 H(i,j)=h2/h1-((h2-h1)/h1)*(theta(i,j)/theta1);%可傾瓦水膜厚度無(wú)量綱表達(dá)式 end end h=h1*H;%有量綱潤(rùn)滑膜厚度矩陣 P=zeros(nx+1,ny+1);%初始化無(wú)量綱節(jié)點(diǎn)壓力 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);%初始化無(wú)量綱壓力矩陣P2 P3=zeros(nx+1,ny+1);%初始化無(wú)量綱壓力矩陣P3 S=0;%初始化誤差表達(dá)式的分子 T=0;%初始化誤差表達(dá)式的分母 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
銅蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時(shí)
- 蟲(chóng)號(hào): 1549417
- 注冊(cè): 2011-12-23
- 性別: GG
- 專(zhuān)業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)
5樓2013-01-10 17:06:02
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +4 | vv迷 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 | ￣^￣゜汗 2026-03-19 | 7/350 |
|
|
|
[考研] 293求調(diào)劑 +12 | zjl的號(hào) 2026-03-16 | 17/850 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|