| 24小時(shí)熱門(mén)版塊排行榜 |
| 8 | 1/1 | 返回列表 |
| 查看: 2941 | 回復(fù): 7 | |||
[求助]
ansys焊接模擬過(guò)程中出現(xiàn)的錯(cuò)誤
|
|
模擬過(guò)程中出現(xiàn)了如下錯(cuò)誤: The value of TEMP at node 434 is 1.594205862E+25. It is greater than the current limit of 1.E+25. This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model is properly constrained. 查看節(jié)點(diǎn)溫度,時(shí)該節(jié)點(diǎn)的溫度超高,且該區(qū)域的溫度都很高。以下是我的命令流,求高手幫忙看看,哪里出了問(wèn)題,謝了~ /CLEAR !清除數(shù)據(jù)庫(kù),新建文件 /FILNAME, SY_FANG_02 !改變工作名 !開(kāi)始熱分析 !建立有限元模型 /NERR,,10000000 /PREP7 !進(jìn)入前處理器 ET,1,PLANE77 !平面熱單元 ET,2,SOLID5 !空間熱單元 V=0.007 !焊接線速度 LENGTH=0.035 !鋼板沿焊縫方向長(zhǎng)度 L_NUMBER=20 !沿焊縫方向單元段數(shù) TIME=0.01 TIME_INC=LENGTH/(L_NUMBER*V) !第一層焊縫時(shí)間增量 ! ! ! ! MPTEMP,1,0,300,600,1000,1300,2000 !定義材料模型,溫度表 MPDATA,DENS,1,1,7859,7770,7659,7600,6750,6750 !密度 MPDATA,DENS,2,1,7859,7770,7659,7600,6750,6750 MPDATA,C,1,1,450,514,967,1370,2800,2000 !比熱 MPDATA,C,2,1,450,514,967,1370,2800,2000 MPDATA,KXX,1,1,67,53,39,31,20,20 !傳導(dǎo)系數(shù) MPDATA,KXX,2,1,67,53,39,31,20,20 MPTEMP,1,0,300,600,1000,1300,1400,2000 !溫度表 MPDATA,EX,1,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7, 2E7 !彈性模量 MPDATA,EX,2,1,2E11,1.86E11,1.35E11,2E10,2E7, 2E7, 2E7 MPDATA,PRXY,1,1,0.33,0.33,0.35,0.36,0.4,0.4 !泊松比 MPDATA,PRXY,2,1,0.33,0.33,0.35,0.36,0.4,0.4 MPDATA,ALPX,1,1,12E-6,12.6E-6,13.0E-6,13.4E-6,13.7E-6,14E-6,13.7E-6 !線膨脹系數(shù) MPDATA,ALPX,1,1,12E-6,12.6E-6,13.0E-6,13.4E-6,13.7E-6,14E-6,13.7E-6 TB,BKIN,1,5 TBTEMP,15 TBDATA,1,280E6,2E9 !屈服極限,切向模量 TBTEMP,800 TBDATA,1,10E6,1E7 TBTEMP,2500 TBDATA,1,10E6,1E6 MPTEMP,,,,,,,, !初始化隨溫度變化的材料參數(shù)表格 MPTEMP,1,20 !定義表格的溫度數(shù)值 MPTEMP,2,250 MPTEMP,3,400 MPTEMP,4,500 MPTEMP,5,750 MPTEMP,6,800 MPTEMP,7,1000 MPTEMP,8,1200 MPTEMP,9,1500 MPTEMP,10,1700 MPTEMP,11,2500 MPDATA,HF,1,,100 !定義對(duì)應(yīng)溫度下的對(duì)流換熱系數(shù)數(shù)值 MPDATA,HF,1,,350 MPDATA,HF,1,,480 MPDATA,HF,1,,520 MPDATA,HF,1,,1000 MPDATA,HF,1,,1050 MPDATA,HF,1,,1500 MPDATA,HF,1,,2000 MPDATA,HF,1,,3000 MPDATA,HF,1,,3100 MPDATA,HF,1,,3500 MPTEMP,,,,,,,, !初始化隨溫度變化的材料參數(shù)表格 MPTEMP,1,20 !定義表格的溫度數(shù)值 MPTEMP,2,100 MPTEMP,3,200 MPTEMP,4,400 MPTEMP,5,500 MPTEMP,6,1000 MPTEMP,7,1100 MPTEMP,8,1500 MPTEMP,9,1700 MPTEMP,10,2500 MPDATA,HF,2,,120 !定義對(duì)應(yīng)溫度下的對(duì)流換熱系數(shù)值 MPDATA,HF,2,,210 MPDATA,HF,2,,300 MPDATA,HF,2,,460 MPDATA,HF,2,,540 MPDATA,HF,2,,1200 MPDATA,HF,2,,1200 MPDATA,HF,2,,2800 MPDATA,HF,2,,3000 MPDATA,HF,2,,3100 ! TBCOPY,BKIN,1,2 !復(fù)制材料BKIN數(shù)據(jù) MP,REFT,1,20 !參考溫度 MP,REFT,2,20 ! ! ! K,1,0,0 K,2,0.015,0 K,3,0.015,0.004 K,4,0,0.004 K,5,0.002,0.004 K,6,0.002,0 A,1,6,5,4 A,6,2,3,5 AGLUE, ALL LESIZE,1,,,3 LESIZE,3,,,3 LESIZE,4,,,4 LESIZE,2,,,4 LESIZE,6,,,4 LESIZE,5,,,10,4 LESIZE,7,,,10,0.25 VEXT,1,2,1,0,0,0.035 VGLUE,ALL LESIZE,14,,,30 TYPE,2 MAT,1 !2號(hào)材料為母材,1號(hào)材料為焊縫 REAL, VSWEEP,1 TYPE,2 MAT,2 REAL, VSWEEP,2 ! FINISH /SOLU !進(jìn)入求解器 DA,7,SYMM !對(duì)稱面A12施加對(duì)稱載荷!已改 DL,5,,UZ DL,1,,UZ DL,4,,UX DL,12,,UX TUNIF,20, TOFFST,273, ! ! !!!!!!!!求解設(shè)置!!!!!!!!!!!!!!!!!!!!!! ANTYPE,4 !瞬態(tài)分析 TRNOPT,FULL LUMPM,0 NLGEOM,1 PSTRES,1 !打開(kāi)應(yīng)力剛化 NSUBST,5,20,1 OUTRES,ERASE OUTRES,ALL,ALL !輸出控制 AUTOTS,1 !打開(kāi)自動(dòng)時(shí)間步長(zhǎng) PRED,ON,,ON LNSRCH,1 !打開(kāi)線性搜索 ! !開(kāi)始激活焊縫單元,進(jìn)行瞬態(tài)熱分析 TIMINT,0,MAG TIMINT,0,ELECT TIMINT,ON !打開(kāi)時(shí)間積分 ! ! !!!!!!!!!!!!!!定義數(shù)組!!!!!!!!!!!!!!!!!!!!!!! *DIM,QQ,ARRAY,8 !定義數(shù)組,存儲(chǔ)單元節(jié)點(diǎn)熱流密度 *DIM,E_NODE,ARRAY,20 !定義數(shù)組,存儲(chǔ)單元所屬節(jié)點(diǎn)編號(hào) *GET,ELEM_MAX,ELEM,0,NUM,MAX !獲取最大單元號(hào) *GET,ELEM_MIN,ELEM,0,NUM,MIN !獲取最小單元號(hào) *GET,NODE_MAX,NODE,0, NUM,MAX !獲取最大節(jié)點(diǎn)號(hào) *GET,NODE_MIN,NODE,0,NUM,MIN !獲取最小節(jié)點(diǎn)號(hào) *DIM,NODE_AVE,ARRAY,NODE_MAX !定義數(shù)組,存儲(chǔ)節(jié)點(diǎn)相鄰單元數(shù)量 ! ! ! ! ESEL,S,MAT,,1 !選擇焊縫單元!!!!!!!!!!!!!!!!! EKILL,ALL !殺死單元!!!!!!!!!!!!!!! ! ! alls!!!!!!!!!!!!!! ! ! R=0.007 !電弧有效加熱半徑 Q=140*20*0.8 !輸入熱量 QMAX=3*Q/(3.1415926*R*R) !變量 XC= 0 !第一層焊縫質(zhì)心X坐標(biāo) YC= 0.007 !第一層焊縫質(zhì)心Y坐標(biāo) ! ! ! !!!!!!!!!!!!從頭至尾激活焊縫,開(kāi)始 *DO,J,1,L_NUMBER alls SFEDELE,ALL,ALL,HFLUX !刪除以前施加的熱流密度 VSEL,S,,,1 !選擇第KK層焊縫所對(duì)應(yīng)的體!!!!!!!!!!!!!!!!!!!!!! ESLV,S,1 !選擇體上單元即第KK層焊縫單元!!!!!!!!!!!!!!!!!!!! NSEL,R,LOC,Z,(J-1)*LENGTH/L_NUMBER,(J+1)*LENGTH/L_NUMBER!!!!!!!!!!!!!!!!!!!!! !選擇第J行單元的中間節(jié)點(diǎn)!!!!!!!!!!!!!!!!!!!!! ESLN,R !選擇第KK層第J行單元!!!!!!!!!!!!!!!! EALIVE,ALL !激活單元!!!!!!!!!!!!!!!!! ALLS !!!!!!!!!!!!!!! ZC=J*LENGTH/L_NUMBER !第KK層第J行單元質(zhì)心Z坐標(biāo) *DO,I,ELEM_MIN,ELEM_MAX !計(jì)算施加單元熱流密度,循環(huán)開(kāi)始 *IF,ESEL(I),EQ,1,THEN ESEL,S,,,I !選擇第I個(gè)單元 *DO,K,1,8 !計(jì)算第I單元前8個(gè)節(jié)點(diǎn)熱流密度 *GET,E_NODE(K),ELEM,I,NODE,K !獲取第I個(gè)單元第K個(gè)節(jié)點(diǎn)的編號(hào) SS=E_NODE(K) DISP=SQRT((NX(SS)-XC)*(NX(SS)-XC)+(NY(SS)-YC)*(NY(SS)-YC)+(NZ(SS)-ZC)*(NZ(SS)-ZC)) EEE=3*DISP*DISP/(R*R) !變量 *IF,DISP,LT,R,THEN QQ(K)=QMAX*EXP(-EEE)/8 !第K個(gè)節(jié)點(diǎn)的熱流密度 *ELSE QQ(K)=0 *ENDIF *ENDDO *IF,QQ(1)+QQ(2)+QQ(3)+QQ(4),NE,0,THEN SFE,I,1,HFLUX,,QQ(1), QQ(2), QQ(3),QQ(4) !在第I個(gè)單元第一個(gè)面上加熱流密度 *ENDIF *IF,QQ(1)+QQ(2)+QQ(5)+QQ(6),NE,0,THEN SFE,I,2,HFLUX,,QQ(1), QQ(2), QQ(6),QQ(5) !在第I個(gè)單元第二個(gè)面上加熱流密度 *ENDIF *IF,QQ(6)+QQ(2)+QQ(3)+QQ(7),NE,0,THEN SFE,I,3,HFLUX,,QQ(2), QQ(3), QQ(7),QQ(6) !在第I個(gè)單元第三個(gè)面上加熱流密度 *ENDIF *IF,QQ(1)+QQ(4)+QQ(5)+QQ(7),NE,0,THEN SFE,I,5,HFLUX,,QQ(1), QQ(4), QQ(7),QQ(5) !在第I個(gè)單元第五個(gè)面上加熱流密度 *ENDIF *IF,QQ(5)+QQ(6)+QQ(7)+QQ(8),NE,0,THEN SFE,I,6,HFLUX,,QQ(5), QQ(6), QQ(7),QQ(8) !在第I個(gè)單元第六個(gè)面上加熱流密度 *ENDIF ALLS!!!!!!!!!!!!!!! *ENDIF *ENDDO !計(jì)算施加單元熱流密度,循環(huán)結(jié)束 ALLS ! ! ! TIME=TIME+TIME_INC !時(shí)間變量 ! TIME,TIME !設(shè)定分析時(shí)間 SOLVE *ENDDO SAVE !保存 FINI |
焊接有限元分析 |
新蟲(chóng) (初入文壇)
|
不知道具體哪里的問(wèn)題,僅把我認(rèn)為不妥的地方說(shuō)下。 1、solid5是8節(jié)點(diǎn)的六面體單元,是沒(méi)有中間節(jié)點(diǎn)的,你用77號(hào)單元分的面網(wǎng)格,sweep后的Z方向應(yīng)該沒(méi)有中間節(jié)點(diǎn),這個(gè)你自己試下,所以激活焊縫時(shí)候可能出錯(cuò)。 2、施加熱流密度時(shí)候?yàn)樯吨患恿?個(gè)面?有點(diǎn)不理解,第四個(gè)面呢? 3、沒(méi)有細(xì)看你的命令流,好像你沒(méi)有施加熱分析的邊界條件,這樣,熱只會(huì)進(jìn)不會(huì)出。 4、焊接這個(gè)東西,三維的話,就用順序耦合吧,看你的命令流應(yīng)該是借鑒了別人的,但又用了不同的耦合方法。 |
|
1.用的單元式solid5耦合單元,77號(hào)單元當(dāng)初是為了體掃略網(wǎng)格劃分的時(shí)候的一個(gè)過(guò)渡單元,我把77號(hào)單元換成不帶中間節(jié)點(diǎn)的面單元重新計(jì)算以后,好像還是不行。 2.熱分析邊界條件加在了材料的屬性里(具體是給材料賦予了對(duì)流換熱屬性) 3.這個(gè)命令流的確是借鑒了別人根據(jù)自己的情況改的,用順序耦合的方法的我去試試 我后來(lái)在solve之前加了 ESEL,S,LIVE NSLE,S 只對(duì)激活單元和節(jié)點(diǎn)進(jìn)行計(jì)算,上面的錯(cuò)誤倒是沒(méi)有了,但是在第二個(gè)循環(huán)以后激活的焊縫單元上面溫度總是上不去,看著就像是跟其他單元隔絕的一樣,一會(huì)兒我把圖貼出來(lái) |
銅蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)

新蟲(chóng) (初入文壇)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | 抓著星星的女孩 2026-03-10 | 3/150 |
|
|
[考研] 285 求調(diào)劑 資源與環(huán)境 一志愿北京化工大學(xué) +3 | 未名考生 2026-03-10 | 3/150 |
|
|
[考研] 求調(diào)劑(材料與化工327) +4 | 愛(ài)吃香菜啦 2026-03-11 | 4/200 |
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
[考研] 420求調(diào)劑 +4 | 莫向外求11 2026-03-10 | 6/300 |
|
|
[基金申請(qǐng)] 提交后的基金本子,已讓學(xué)校撤回了,可否換口子提交 +3 | dut_pfx 2026-03-10 | 3/150 |
|