| 1 | 1/1 | 返回列表 |
| 查看: 385 | 回復(fù): 0 | ||
wjj0822新蟲 (初入文壇)
|
[求助]
堆焊熱分析 ——總是提示無法收斂,請高手指教
|
|
由于項目需要,想做一個堆焊過程仿真,觀察基體溫度的變化過程,但在第二個子步及t=0.4s時提示無法收斂,不知道是什么原因,麻煩高人幫忙看一下,指出問題所在,謝謝了! /filname,weldtem /title,WELDING TEMPERATIUTE FIELD SIMULATION /units,si /prep7 et,1,solid70 keyopt,1,2,1 et,2,surf152 r,2 keyopt,2,4,1 keyopt,2,5,0 keyopt,2,8,3 mptemp,1,0,100,200,300,400,500 !定義材料1屬性溫度范圍 mptemp,,1200,1250,1350,1400,1450,1500 mptemp,,1505,1705,1905,2105,2305,2500 mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定義材料1密度(kg/m3) mpdata,dens,1,,7800,7800,7800,7800,7800,7800 mpdata,dens,1,,7800,7800,7800,7800,7800,7800 mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !熱傳導(dǎo)系數(shù)(W/(m*K)) mpdata,kxx,1,,30,30,30,30,30,30 mpdata,kxx,1,,30,30,30,30,30,30 mpdata,c,1,1,450,469,481,508.5,536,569 !考慮相變潛熱的比熱容(J/kg*K) mpdata,c,1,,700,2172,5116,6589,8061,9533 mpdata,c,1,,9533,7757,5982,4206,2431,700 !基體材料單元類型 mpdata,dens,2,1,7850,7850,7800,7800,7800,7800 !定義材料2密度(kg/m3) mpdata,dens,2,,7800,7800,7800,7800,7800,7800 mpdata,dens,2,,7800,7800,7800,7800,7800,7800 mpdata,kxx,2,1,52,50.7,48.6,46.1,42.3,38.9 !熱傳導(dǎo)系數(shù)(W/(m*K)) mpdata,kxx,2,,30,30,30,30,30,30 mpdata,kxx,2,,30,30,30,30,30,30 mpdata,c,2,1,450,469,481,508.5,536,569 !考慮相變潛熱的比熱容(J/kg*K) mpdata,c,2,,700,2172,5116,6589,8061,9533 mpdata,c,2,,9533,7757,5982,4206,2431,700 pi=3.1415926 I=160 U=23 efficient=0.78 r=0.0045 Qm=3*I*U*efficient/(pi*r**2) kk=0.85 Qmax=Qm*kk ly=0.02 lx=0.02 lz=0.01 lsize=0.001 !單元尺寸 v=0.0025 tm=lx/v !移動的總時間 tinc=lsize/v !在每個單元的移動時間 max_time=1+tm/tinc max_x=1+lx/lsize !x方向節(jié)點個數(shù) max_y=1+ly/lsize !y方向節(jié)點個數(shù) /view,1,1,1,1 block,0,lx,-ly/2,ly/2,0.0003,lz+0.0003 !建模 block,0,lx,-ly/2,ly/2,0,0.0003 vglue,all numcmp,all save esize,lsize !單元的大小 mshkey,1 !Use mapped meshing type,1 mat,1 vmesh,1 esize,lsize/4 type,2 mat,2 mshkey,1 vmesh,2 /pnum,line,off /pnum,defa eplot finish --------------- *dim,flux2,table,max_x,max_y,max_time,x,y,time *do,k,1,max_time,1 *do,i,1,max_x,1 flux2(i,0,k)=(i-1)*lsize !經(jīng)過x方向的距離 *enddo *do,j,1,max_y,1 flux2(0,j,k)=(j-1)*lsize !經(jīng)過y方向的距離 *enddo *do,i,1,max_x,1 *do,j,1,max_y,1 xcenter=v*(k-1)*tinc ycenter=0 distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2) *if,distance,le,r,then flux2(i,j,k)=Qmax*exp(-3*distance**2/(r**2)) *else flux2(i,j,k)=0 *endif *enddo *enddo flux2(0,0,k)=(k-1)*tinc *enddo /solu antype,trans tunif,15 kbc,1 timint,on autots,off esel,s,type,,2 nsle,s sf,all,conv,-1 ! -1; current flow into the element face (IN ) allsel *do,t,1e-6,tm,tinc time,t autos,0 deltim,,,,1 kbc,0 outres,all,all tsres,erase sf,all,hflux,%flux2% solve *enddo Finish /seg,off,WJ,0.2 anim,1,1 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|