| 4 | 1/1 | 返回列表 |
| 查看: 2647 | 回復: 3 | |||
[求助]
焊接殘余應力場問題 已有1人參與
|
|
請大家指點, 剛開始學做焊接和殘余應力場問題。采用間接費進行 熱場與應力場的分析,命令流執(zhí)行完成后,沒有出現(xiàn) .rth 文件,無法進行應力場的分析,求大家指點。 命令流如下: !圓管焊接模擬 /CLEAR, /config,nres,5000 /title,weld /prep7 /unit,si !單位: International system (m, kg, s, K) /PREP7 *SET,RAD,0.250 *SET,THICK,0.008 *SET,HIGH,1 !管高度,焊縫長度 *SET,enum,8 *SET,div1,2 !沿厚度方向的單元層數(shù) !*SET,E,2.06e5 !彈性模量 !*SET,velocity,1 !焊接速度 mm/s RAD2=RAD-THICK esize=THICK/div1 max_tem=1500 ! 按照前面假設,焊料的初始溫度為1500℃ dt1=1e-3 ! 用于建立初始條件的一個很小的時間段 dt=2 ! 焊接一個單元所需的時間 t=0 ! 起始時間 tcool=500 !焊縫冷卻時間 !*****單元***** ET,1,SOLID70 ET,4,SURF154 !*****材料屬性***** MPTEMP,,,,,,,, MPTEMP,1,20 !溫度數(shù)據(jù) MPTEMP,2,300 MPTEMP,3,600 MPTEMP,4,1000 MPTEMP,5,1300 MPTEMP,6,2000 MPDATA,EX,1,,2e11 !彈性模量數(shù)據(jù) MPDATA,EX,1,,1.96e11 MPDATA,EX,1,,1.45e11 MPDATA,EX,1,,2.1e10 MPDATA,EX,1,,2e7 MPDATA,EX,1,,1e6 MPDATA,PRXY,1,,0.3 MPDATA,PRXY,1,,0.3 MPDATA,PRXY,1,,0.3 MPDATA,PRXY,1,,0.3 MPDATA,PRXY,1,,0.3 MPDATA,PRXY,1,,0.3 !下面是各材料在各溫度點下的屈服應力和屈服后的彈性模量 MPTEMP,,,,,,,, TB,BKIN,1,6 TBTEMP,20,1 TBDATA,1,0.345e9,2e10 TBTEMP,500,2 TBDATA,1,0.243e9,1.45e10 TBTEMP,735,3 TBDATA,1,0.045e9,0.329e10 TBTEMP,900,4 TBDATA,1,0.017e9,0.077e10 TBTEMP,1300,5 TBDATA,1,0.017e7,0.077e8 TBTEMP,1600,6 TBDATA,1,0.017e6,0.077e7 mp,dens,1,7850 !密度 mp,alpx,1,1.4e-5 !線膨脹系數(shù) mptemp mptemp,1,20 mptemp,2,500 mptemp,3,735 mptemp,4,900 mptemp,5,1000 mptemp,6,1200 mptemp,7,1300 mptemp,8,1400 mptemp,9,1500 mptemp,10,1600 mpdata,kxx,1,,53.33 mpdata,kxx,1,,37.35 mpdata,kxx,1,,29.52 mpdata,kxx,1,,27.3 mpdata,kxx,1,,27.3 mpdata,kxx,1,,27.3 mpdata,kxx,1,,27.3 mpdata,kxx,1,,27.3 mpdata,kxx,1,,27.3 mpdata,kxx,1,,27.3 !不同溫度點的導熱系數(shù) mpdata,c,1,,439.8 mpdata,c,1,,666.5 mpdata,c,1,,1004.83 mpdata,c,1,,650.44 mpdata,c,1,,650 mpdata,c,1,,666.5 mpdata,c,1,,666.5 mpdata,c,1,,666.5 mpdata,c,1,,666.5 mpdata,c,1,,666.5 !建模 CYLIND, RAD, RAD2, 0,HIGH ,0,359 CYLIND, RAD, RAD2, 0,HIGH ,359,360 VGLUE,All, !csys , 11 !劃分網(wǎng)格 asel , s , loc , z , 0 LSLA, S !LESIZE,ALL, , , ,ENUM , , , ,1, LESIZE,ALL,esize , , , , , , ,1, AATT,1,1,4,0, MSHKEY,0 MSHAPE,0,2D amesh,all allsel VATT,1,1,1,0 ESIZE, , 30 EXTOPT,ACLEAR,1 EXTOPT,VSWE,TETS,1 vsweep,all nummrg,all,,,,low numcmp,all !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *get,nvol,volu,,NUM,max !確定焊縫實體的編號 VSEL, s, volu, , nvol eslv,s, eplot *GET, nelem, elem, 0, count !焊縫單元的數(shù)目 *dim,elem,,nelem !焊縫單元的單元號 *dim,ecordz,,nelem !焊縫單元的單元z坐標 *get,elem(1) ,elem,,num,min !焊縫單元的最小單元號 ecordz(1)=centrz(elem(1)) *do,i,2,nelem elem(i)=elnext(elem(i-1)) !獲取焊縫單元號 ecordz(i)=centrz(elem(i)) !獲取焊縫單元中心Z坐標 *enddo !將單元號按照Z坐標排序 *do,i,2,nelem elem(i)=elnext(elem(i-1)) !獲取焊縫單元號 ecordz(i)=centrz(elem(i)) !獲取焊縫單元中心Z坐標 *enddo !冒泡法 按z坐標大小排序 *do,i,1,nelem-1 *do,j,i+1,nelem *if,ecordz(j),lt,ecordz(i),then tmp=ecordz(i) ecordz(i)=ecordz(j) ecordz(j)=tmp tmp=elem(i) elem(i)=elem(j) elem(j)=tmp *endif *enddo *enddo alls *dim,eorder,,nelem/(div1*2),div1*2,, !定義單元激活數(shù)組,每次激活z坐標相同的(div1*2)個單元 *do,i,1,nelem temp1=i/(div1*2) temp2=mod(i,(div1*2)) !返回余數(shù)部分 temp3=nint(temp1+0.5) !返回整數(shù)部分!對i/(div1*2)取整,nint(x)是找與x最近的整數(shù) *if,temp2,eq,0,then eorder((temp3-1),4)=elem(i) ! temp-1 修正當i能被(div1*2)整除時temp3 會偏大 *else eorder(temp3,temp2)=elem(i) *endif *enddo !輸入初始溫度 !殺死焊縫單元 *do,i,1,nelem ekill,elem(i) *enddo esel,s,live eplot /SOLU ANTYPE,4 TRNOPT,FULL NSEL,ALL neqit,20 esel,all eplot /auto,1 /replot allsel,all outres,all,all ic,all,temp,20 !設置初始溫度 kbc,1 !設置加載方式 timint,0,struct !打開結構瞬態(tài)分析時間積分選項 timint,1,therm !打開熱瞬態(tài)分析時間積分選項 timint,0,mag !打開磁瞬態(tài)分析時間積分選項??????? tintp,0.005,,,1,0.5,0.2 !設置瞬態(tài)分析參數(shù)(除最后一個0.2外都是采用的默認參數(shù)) nsub1=2 nsub2=40 aa=nelem/(div1*2) !單元激活組的個數(shù) bb=(div1*2) !同一個z坐標下 單元的個數(shù)(以這個個數(shù)bb為一個單元激活組) esel,none t=0 *do,i,1,aa esel,none *do,j,1,bb esel,a,elem,,eorder(i,j) *enddo ealive,all ! bfe,all,hgen,,4.8e9 !施加焊縫初始溫度 esel,s,live eplot esel,all ! 下面的求解用于建立溫度的初始條件 t=t+dt1 time,t nsubst,1 *do,k,1,bb *do,j,1,4 d,nelem(eorder(i,k),j),temp,max_tem !施加溫度荷載 *enddo *enddo solve ! 下面的求解用于保證初始的升溫速度為零 ??????????????? t=t+dt1 time,t solve ! 下面的步驟用于求解溫度分布 *do,k,1,bb *do,j,1,4 ddele,nelem(eorder(i,k),j),temp !刪除溫度荷載 *enddo *enddo allsel ! ddele,all,temp t=t+dt-2*dt1 time,t nsubst,nsub1 solve *enddo *do,tt,t+2,t+tcool,2 time,tt nsubst,1 solve *enddo finish |
ANSYS命令流 |
銅蟲 (初入文壇)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗: +809 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學調劑 ,六級已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|---|---|---|---|---|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 287求調劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 328求調劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調劑 2026-03-16 | 12/600 |
|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 354求調劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 328求調劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調劑 2026-03-17 | 7/350 |
|
|
[考研] 生物學071000 329分求調劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 268求調劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 332求調劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 一志愿蘇州大學材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 085601材料工程315分求調劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 266求調劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|