 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2648 | 回復: 3 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
焊接殘余應力場問題 已有1人參與
|
|||
|
請大家指點， 剛開始學做焊接和殘余應力場問題。采用間接費進行 熱場與應力場的分析，命令流執(zhí)行完成后，沒有出現 .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 !沿厚度方向的單元層數 !*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 !溫度數據 MPTEMP,2,300 MPTEMP,3,600 MPTEMP,4,1000 MPTEMP,5,1300 MPTEMP,6,2000 MPDATA,EX,1,,2e11 !彈性模量數據 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 !線膨脹系數 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 !不同溫度點的導熱系數 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 !劃分網格 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 !焊縫單元的數目 *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,, !定義單元激活數組，每次激活z坐標相同的（div1*2）個單元 *do,i,1,nelem temp1=i/(div1*2) temp2=mod(i,(div1*2)) !返回余數部分 temp3=nint(temp1+0.5) !返回整數部分!對i/(div1*2)取整，nint（x）是找與x最近的整數 *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)分析參數（除最后一個0.2外都是采用的默認參數） nsub1=2 nsub2=40 aa=nelem/(div1*2) !單元激活組的個數 bb=(div1*2) !同一個z坐標下 單元的個數（以這個個數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命令流 |
» 猜你喜歡
一志愿重慶大學085700資源與環(huán)境總分308求調劑
已經有6人回復
-
268求調劑
已經有7人回復
-
302求調劑
已經有11人回復
-
306求調劑
已經有4人回復
-
材料學學碩080502 337求調劑-一志愿華中科技大學
已經有4人回復
-
求調劑
已經有6人回復
-
313求調劑
已經有3人回復
-
南昌大學材料專碩311分求調劑
已經有6人回復
-
316求調劑
已經有6人回復
-
346求調劑[0856]
已經有7人回復
» 本主題相關商家推薦: (我也要在這里推廣)
shgao20
專家顧問 (職業(yè)作家)
-
專家經驗: +809 - 仿真EPI: 1
- 應助: 1311 (講師)
- 金幣: 13381.4
- 紅花: 267
- 帖子: 3531
- 在線: 473.4小時
- 蟲號: 3644138
- 注冊: 2015-01-14
- 專業(yè): 動力學與控制
- 管轄: 仿真模擬
3樓2018-10-01 09:14:37
nightrain2
銅蟲 (初入文壇)
2樓2018-09-29 07:17:01
4樓2018-10-02 11:14:45
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 南昌大學材料專碩311分求調劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 機械專碩299求調劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 265求調劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 【考研調劑】化學專業(yè) 281分，一志愿四川大學，誠心求調劑 +8 | 吃吃吃才有意義 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 354求調劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 329求調劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿武理材料工程348求調劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿華中農業(yè)071010，總分320求調劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 288求調劑，一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 復試調劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 0703化學調劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[考研] 328求調劑，英語六級551，有科研經歷 +4 | 生物工程調劑 2026-03-16 | 12/600 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研]
 11408 一志愿西電，277分求調劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 321求調劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學調劑 290分有科研經歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|