| 查看: 3004 | 回復(fù): 11 | ||
| 【懸賞金幣】回答本帖問題,作者ywq741將贈送您 9 個金幣 | ||
ywq741新蟲 (初入文壇)
|
[求助]
ANSYS焊接生死單元問題,為什么死單元會有溫度,該怎么解決問題? 已有2人參與
|
|
|
最近正在做ANSYS焊接溫度場模擬,運用生死單元,模擬每層焊接運動。試過用d還進(jìn)行約束,激活后解除約束ddele,但是發(fā)現(xiàn)最后溫度會有負(fù)的溫度。請問有什么方法可以實現(xiàn)啊。!謝謝! /PREP7 ET,1,SOLID70 MP,DENS,1,7930 MP,C,1,502 mptemp,1,24,40,60,80,110,150 mpdata,kxx,1,1,9.33,10.42,11.58,12.7,14,18 /VIEW,1,1,1,1 !*************定義基板尺寸*********** lx=0.080 lz=0.060 ly=0.005 !*************定義焊接參數(shù)*********** weld=0.012 !每道焊道寬度 wheight=0.003 !焊道的高度 LSIZE=0.002 V=0.004 !速度 tinc=LSIZE/V !每小段的時間間隔 pi=3.1415926 U=25 I=180 Q=U*I !電源有效功率 effect=0.8 f1=0.6 f2=1.4 ah=0.003 chf=2*ah chb=8*ah bh=0.003 Qmax11=6*effect*3**(1/2)*f1*q/(pi*pi**(1/2)*ah*bh*chf) Qmax12=6*effect*3**(1/2)*f2*q/(pi*pi**(1/2)*ah*bh*chb) !*****************建立模型,生成網(wǎng)格************** block,0,lx,0,ly,0,lz wpoff,0,ly,lz/2 block,0,lx,0,wheight*5,-weld/2,weld/2 wpoff,0,0,-weld/2 vsbw,all,,delete wpoff,0,0,weld vsbw,all,,delete wpoff,0,-ly,-weld/2 block,-0.02,0,0,ly,-lz/2,lz/2 block,lx,lx+0.02,0,ly,-lz/2,lz/2 wpoff,0,ly,weld/2 vglue,all FLST,5,4,4,ORDE,4 FITEM,5,14 FITEM,5,16 FITEM,5,19 FITEM,5,39 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,40, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,23 FITEM,5,-24 FITEM,5,69 FITEM,5,71 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,13 FITEM,5,15 FITEM,5,70 FITEM,5,72 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,10, , , , ,1 !* MSHAPE,0,3D MSHKEY,1 !* CM,_Y,VOLU VSEL, , , , 8 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VMESH,_Y1 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !* ESIZE,2e-3,0, /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST /DIST,1,1.08222638492,1 /REP,FAST FLST,5,3,6,ORDE,3 FITEM,5,1 FITEM,5,7 FITEM,5,9 CM,_Y,VOLU VSEL, , , ,P51X CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VMESH,_Y1 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !* FLST,5,2,6,ORDE,2 FITEM,5,10 FITEM,5,-11 CM,_Y,VOLU VSEL, , , ,P51X CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VSWEEP,_Y1 !* CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !* /UI,MESH,OFF /solu antype,trans trnopt,full nropt,full,,on AUTOTS,on lnsrch,on tunif,25 !**************殺死焊縫單元******************** nsel,s,loc,y,ly,ly+wheight*5 esln,s,1 CM,weld_elem,ELEM ekill,all esel,s,elem,,weld_elem nsle,s cm,weld_node,node ALLSEL,ALL esel,s,live eplot !**********施加對流載荷************* ! 熔積層暴露空氣,輻射和對流 ALLSEL,ALL NSEL,S,LOC,y,0,0.020 NSEL,R,EXT SF,ALL,CONV,rad,25 !空氣 !**********約束殺死單元的溫度約束************* nsel,s,node,,weld_node D,ALL,temp,298 !************定義數(shù)組維數(shù)********** *do,j,1,5,1 MAX_X=1+lx/LSIZE !*************定義table數(shù)組****************************************** *do,i,1,MAX_X,1 XX=(I-1)*LSIZE CSYS,0 nsel,s,loc,y,ly+wheight*(j-1),ly+wheight*j nsel,r,loc,z,lz/2-weld/2,lz/2+weld/2 esln,s,1 DDELE,ALL allsel,all !!!!!!!生死單元的處理 csys,0 nsel,s,loc,y,ly+wheight*(j-1),ly+wheight*j nsel,r,loc,z,lz/2-weld/2,lz/2+weld/2 nsel,r,loc,x,XX-chf/2,XX+chf/2 esln,s,1 ealive,all eplot !!!!!!雙橢球體熱源的加載 ALLSEL,ALL !!!雙橢球熱源的加載 !!!前半球節(jié)點的選取 local,11,2,XX,ly+wheight*j,lz/2,,,,bh/chf,ah/chf, NSEL,S,LOC,X,0,chf NSEL,R,LOC,Y,270,360 NSEL,R,LOC,Z,0,180 CM,FRONT,NODE ALLSEL,ALL !!!!!!!!!!!!!!!!!!!!!!!!!!! LOCAL,13,0,XX,ly+wheight*j,lz/2, , , ,1,1, *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNCCSYS *DEL,'front' *SET,_FNCNAME,'front' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *SET,_FNC_C1(1),qmax11 *SET,_FNC_C2(1),chf *SET,_FNC_C3(1),bh *SET,_FNC_C4(1),ah *SET,_FNCCSYS,13 ! /INPUT,front.func,,,1 *DIM,%_FNCNAME%,TABLE,6,20,1,,,,%_FNCCSYS% ! ! Begin of equation: Qmax11*exp(-3(({x}/chf)^2+({y}/bh)^2+({z}/ah)^2)) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 4, 18 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, 3, 4, 19 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,9,1), 0.0, -5, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, 4, 4, 20 *SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 2, 0, 0, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -2, 17, -4 *SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -1, 1, -5 *SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 3, 0, 0, -2 *SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 1, -3, 3, -4 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 17, 3, -1 *SET,%_FNCNAME%(0,20,1), 0.0, 99, 0, 1, -2, 0, 0 ! End of equation: Qmax11*exp(-3(({x}/chf)^2+({y}/bh)^2+({z}/ah)^2)) !--> BF,FRONT,HGEN,%front% CMDELE,FRONT ALLSEL,ALL !!!后半球節(jié)點的選取 local,12,2,XX,ly+wheight*j,lz/2,,,,bh/chb,ah/chb, NSEL,S,LOC,X,0,chb NSEL,R,LOC,Y,180,270 NSEL,R,LOC,Z,0,180 CM,BEHIND,NODE ALLSEL,ALL !!!!!!!!!!!!!!!!!!!!!!!!!!! LOCAL,14,0,XX,ly+wheight*j,lz/2, , , ,1,1, *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNCCSYS *DEL,'back' *SET,_FNCNAME,'back' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *SET,_FNC_C1(1),qmax12 *SET,_FNC_C2(1),chb *SET,_FNC_C3(1),bh *SET,_FNC_C4(1),ah *SET,_FNCCSYS,0 ! /INPUT,back.func,,,1 *DIM,%_FNCNAME%,TABLE,6,20,1,,,,%_FNCCSYS% ! ! Begin of equation: Qmax12*exp(-3(({x}/chb)^2+({y}/bh)^2+({z}/ah)^2)) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 2, 4, 18 *SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, 3, 4, 19 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 2, 0, 0, -1 *SET,%_FNCNAME%(0,9,1), 0.0, -5, 0, 1, -1, 17, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 1, 4, 4, 20 *SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 2, 0, 0, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -2, 17, -4 *SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 1, -1, 1, -5 *SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 3, 0, 0, -2 *SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 1, -3, 3, -4 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, 17, 3, -1 *SET,%_FNCNAME%(0,20,1), 0.0, 99, 0, 1, -2, 0, 0 ! End of equation: Qmax12*exp(-3(({x}/chb)^2+({y}/bh)^2+({z}/ah)^2)) !--> BF,BEHIND,HGEN,%back% CMDELE,BEHIND ALLSEL,ALL time,i*tinc+50*(j-1) nsubst,10 SOLVE BFDELE,all,HGEN esel,s,live eplot OUTRES,ALL,ALL, *ENDDO *enddo *DO,ii,1,120 csys,0 nsel,s,ext sf,all,conv,30,25 ALLSEL,ALL TIME,5*ii+250 solve *ENDDO |
新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
鐵蟲 (初入文壇)
|
這樣做是可以的,但是會導(dǎo)致活單元與死單元相鄰部位的溫度結(jié)果不正確:因為你設(shè)置死單元的溫度為 0,實際上是設(shè)置死單元的所有節(jié)點的溫度為 0,而死單元的邊界節(jié)點同時又屬于相鄰的活單元,那么就相當(dāng)于給相鄰活單元的節(jié)點施加了一個固定的 0 溫度條件,從而導(dǎo)致活單元的計算結(jié)果不正確。 計算時需要死單元也參加計算,但是顯示結(jié)果是可以不顯示死單元的結(jié)果,方法很簡單:顯示結(jié)果前,先選擇所有活單元以及屬于活單元的所有節(jié)點,然后在顯示結(jié)果時就只有活單元的結(jié)果了。如果這種情況仍然有負(fù)的溫度,那就是模型中有問題了?梢韵燃せ钊繂卧,算一次看看。 |
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎211本科同學(xué),過A區(qū)國家線,A區(qū)非偏遠(yuǎn)一本,交叉學(xué)科課題組 +8 | lisimayy 2026-03-04 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)工程求調(diào)劑 +10 | 化工人999 2026-03-04 | 10/500 |
|
|
[考研] 085600材料與化工 298 調(diào)劑 +8 | 小西笑嘻嘻 2026-03-03 | 8/400 |
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計算機(jī)、人工智能等專業(yè)的同學(xué)報考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
[基金申請] 請問大家,研究風(fēng)險與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 理學(xué),工學(xué),農(nóng)學(xué)調(diào)劑,少走彎路,這里歡迎您! +8 | likeihood 2026-03-02 | 11/550 |
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
[基金申請] 成果系統(tǒng)訪問量大,請一小時后再嘗試。---NSFC啥時候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|