 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1502 | 回復: 7 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
LQgraduate銅蟲 (初入文壇)
|
[求助]
求解軸承潤滑性能時，如果承載力不收斂，應該如何調(diào)整油膜厚度？ 已有1人參與
|
||
|
問題如下（希望做過這方面工作的好心人幫忙解決一下）： 計算流體動壓推力軸承或者徑向軸承的潤滑性能時，首先需要利用迭代法使壓力收斂，從而求解雷諾方程；然后要使軸承的承載力也收斂，如果不收斂，需要修正最初假設(shè)的最小和最大油膜厚度，直到收斂為止。問題來了，就是如果承載力不收斂，應該如何調(diào)整油膜厚度？有什么方法嗎？還是純粹靠經(jīng)驗？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 機械科研經(jīng)驗專輯 |
» 猜你喜歡
機械352分求調(diào)劑，本科太原理工大學
已經(jīng)有1人回復
-
招收0805（材料）調(diào)劑
已經(jīng)有0人回復
-
工程熱物理與能源利用論文潤色/翻譯怎么收費?
已經(jīng)有292人回復
-
招收0805（材料）調(diào)劑
已經(jīng)有3人回復
-
EBSD極圖和反極圖怎么看
已經(jīng)有0人回復
-
【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】
已經(jīng)有8人回復
-
一志愿齊工大 085501機械工程 304分求調(diào)劑
已經(jīng)有2人回復
-
085501機械工程304分求調(diào)劑，一志愿齊工大
已經(jīng)有15人回復
-
AI有限融合復合材料多尺度建模
已經(jīng)有0人回復
-
304分085501機械工程 一志愿齊工大
已經(jīng)有12人回復
-
求收留
已經(jīng)有5人回復
LQgraduate
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時
- 蟲號: 1549417
- 注冊: 2011-12-23
- 性別: GG
- 專業(yè): 機械摩擦學與表面技術(shù)
5樓2013-01-10 17:06:02
Tennisgraf
銀蟲 (正式寫手)
- 應助: 5 (幼兒園)
- 金幣: 13.4
- 散金: 162
- 紅花: 2
- 帖子: 850
- 在線: 206.7小時
- 蟲號: 691801
- 注冊: 2009-01-12
- 專業(yè): 機械摩擦學與表面技術(shù)
2樓2013-01-10 14:49:31
LQgraduate
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 227
- 紅花: 1
- 帖子: 23
- 在線: 23.8小時
- 蟲號: 1549417
- 注冊: 2011-12-23
- 性別: GG
- 專業(yè): 機械摩擦學與表面技術(shù)
|
首先，謝謝你的回復。但是，問題在于，你的意思是說要靠經(jīng)驗調(diào)整了，反復嘗試對吧，知道承載力收斂為止。那么有沒有什么方法，比如某個算法，可以實現(xiàn)編程然后讓油膜厚度自動調(diào)整直到承載力收斂為止��？我在文獻上看到了一種方法，說是用二分法可以，但我試了一下，做出來啊!這樣吧，我把程序貼到上邊來，你看幫我看一下問題出在哪里吧？我加了注釋了，如果有看不懂的地方，再交流。以下是我用matlab編寫的五點差分法離散雷諾方程用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;%軸瓦長度 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è): 機械摩擦學與表面技術(shù)
4樓2013-01-10 15:31:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科求調(diào)劑286 +3 | tgs_001 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 2026年華南師范大學歡迎化學，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學ing 2026-03-28 | 3/150 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 考研化學308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 0703化學一志愿南京師范大學303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
|
[考研] 26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 上海電力大學材料防護與新材料重點實驗室招收調(diào)劑研究生（材料、化學、電化學，環(huán)境） +4 | 我愛學電池 2026-03-23 | 4/200 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|