| 5 | 1/1 | 返回列表 |
| 查看: 2649 | 回復: 3 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
nightrain2銅蟲 (初入文壇)
|
[求助]
焊接殘余應(yīng)力場問題 已有1人參與
|
||
|
請大家指點, 剛開始學做焊接和殘余應(yīng)力場問題。采用間接費進行 熱場與應(yīng)力場的分析,命令流執(zhí)行完成后,沒有出現(xiàn) .rth 文件,無法進行應(yīng)力場的分析,求大家指點。 命令流如下: !圓管焊接模擬 /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 ! 按照前面假設(shè),焊料的初始溫度為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 !下面是各材料在各溫度點下的屈服應(yīng)力和屈服后的彈性模量 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 !設(shè)置初始溫度 kbc,1 !設(shè)置加載方式 timint,0,struct !打開結(jié)構(gòu)瞬態(tài)分析時間積分選項 timint,1,therm !打開熱瞬態(tài)分析時間積分選項 timint,0,mag !打開磁瞬態(tài)分析時間積分選項??????? tintp,0.005,,,1,0.5,0.2 !設(shè)置瞬態(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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學調(diào)劑 ,六級已過,有科研經(jīng)歷 +14 | 曦熙兮 2026-03-15 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[考研] 生物學一志愿985,分數(shù)349求調(diào)劑 +3 | zxts12 2026-03-21 | 3/150 |
|
|
[基金申請]
學校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 材料學學碩080502 337求調(diào)劑-一志愿華中科技大學 +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 0703化學336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 東南大學364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|