 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1640 | 回復(fù): 1 | ||
hjb1109鐵蟲 (初入文壇)
|
[求助]
ANSYS用生死單元實(shí)現(xiàn)多層逐步加熱，加載中存在小問題，求解
|
|
現(xiàn)在做的工作是對一種算法的有效性演算，所以建模和加載都進(jìn)行了簡化。 需要實(shí)現(xiàn)的功能是利用生死單元實(shí)現(xiàn)多層逐步加熱：就是先激活底層，在底層加載；底層加熱完成再激活上面的第二層，再進(jìn)行加載；第二層加熱完成后再激活上面的第三層…… 加載用的熱源模型是錐形體熱源，就是高斯激光作用于半無限大物體表面時(shí)所產(chǎn)生的熱效應(yīng)。就是這里出了問題，加載第一層還能夠顯示出錐形加熱，但是到了第二層第三層時(shí)不再顯示出是錐形，而是圓柱形。在計(jì)算過程中ANSYS也給出了警告：在選取節(jié)點(diǎn)進(jìn)行加載時(shí)不能識別Z方向的NSEL命令。 希望能有人分析一下原因。下面是APDL程序。 /units,si /title, shengsi2204 /filn,shengsi2204 /replot !后面要用到的一些參數(shù)，如模型的長寬高，激光功率及相關(guān)參數(shù)(光斑半徑、作用深度等) len=2/1000 wid=2/1000 hth=0.5/1000 q=15 h=1.5/1000 r1=0.12/1000 pi=3.1415926 qv=2*q/(pi*h*(r1**2)) lsize=0.1/1000 !為后面建立體熱源表格做準(zhǔn)備，三維表格每個(gè)維度的數(shù)量 max_x=len/lsize+1 max_y=wid/lsize+1 max_z=hth/lsize+1 !進(jìn)入前處理器 /prep7 et,1,solid70 mp,dens,1,1000 mp,kxx,1,0.25 mp,c,1,1700 mp,emis,1,0.95 block,0,len,0,wid,0,hth block,0,len,0,wid,hth,2*hth block,0,len,0,wid,2*hth,3*hth /view,1,1,1,1 /replot allsel,all esize,0.1/1000 mat,1 vsel,s,volu,,1 vsel,a,volu,,2 vsel,a,volu,,3 vmesh,all !這里是接觸設(shè)置，對整個(gè)程序沒有很大影響，是計(jì)算實(shí)際需要，很長一段都可以不看 !接觸設(shè)置開始 /COM, contact pair creation - start allsel,all asel,s,,,2 nsla,s,1 cm,target1,node asel,s,,,7 nsla,s,1 cm,contact1,node allsel,all mp,mu,1 mat,1 r,3 real,3 et,3,170 et,4,174 keyopt,4,9,0 keyopt,4,10,2 r,3, rmore,,,,,, rmore,,2500 rmore,0 keyopt,4,1,2 nsel,s,,,target1 type,3 esln,s,0 esurf,all nsel,s,,,contact1 type,4 esln,s,0 esurf,all allsel,all asel,s,,,8 nsla,s,1 cm,target1,node asel,s,,,13 nsla,s,1 cm,contact1,node allsel,all mp,mu,1 mat,1 r,3 real,3 et,3,170 et,4,174 keyopt,4,9,0 keyopt,4,10,2 r,3, rmore,,,,,, rmore,,2500 rmore,0 keyopt,4,1,2 nsel,s,,,target1 type,3 esln,s,0 esurf,all nsel,s,,,contact1 type,4 esln,s,0 esurf,all !接觸設(shè)置結(jié)束 /COM, contact pair creation - end finish !進(jìn)入求解器 /solu !定義錐形體熱源，在水平方向上高斯分布，在深度方向上線性遞減 *dim,gene1,table,max_x,max_y,max_z,x,y,z xc=len/2 yc=wid/2 *do,k,1,max_z,1 gene1(0,0,k)=(k-1)*lsize *do,i,1,max_x,1 gene1(i,0,k)=(i-1)*lsize *enddo *do,j,1,max_y,1 gene1(0,j,k)=(j-1)*lsize *enddo *do,i,1,max_x,1 *do,j,1,max_y,1 distance= (((i-1)*lsize-xc)**2+((j-1)*lsize-yc)**2)**(1/2) hei=1+(k-max_z)*lsize/h *if,distance,le,r1,then gene1(i,j,k)=Qv*EXP(-2*(distance/r1)**2)*hei *else gene1(i,j,k)=0 *endif *enddo *enddo *enddo !定義分析類型及一些設(shè)置 antype,trans trnopt,full tunif,100 allsel *set,delt,1 timint,on autots,on !將V2、V3所含單元?dú)⑺�，只激活V1所含單元 allsel,all vsel,s,volu,,2 vsel,a,volu,,3 eslv,s,1 ekill,all esel,s,live eplot !0 !0.2<=t<0.3時(shí)，加載第三層 *do,i,1e-6,1,0.01 time,i deltim,i allsel,all *if,i,lt,0.1,then *if,yc,lt,2,then allsel,all nsel,s,x,loc,xc-r1,xc+r1 nsel,r,y,loc,yc-r1,yc+r1 nsel,r,z,loc,0,hth bf,all,hgen,%gene1% allsel,all *endif *elseif,i,ge,0.1,then allsel,all bfdele,all,hgen vsel,s,volu,,2 eslv,s,1 ealive,all eplot *if,i,lt,0.2,then *if,yc,lt,2,then allsel,all nsel,s,x,loc,xc-r1,xc+r1 nsel,r,y,loc,yc-r1,yc+r1 nsel,r,z,loc,hth,2*hth bf,all,hgen,%gene1% allsel,all *endif *elseif,i,ge,0.2,then allsel,all bfdele,all,hgen vsel,s,volu,,3 eslv,s,1 ealive,all eplot *if,i,lt,0.3,then *if,yc,lt,2,then allsel,all nsel,s,x,loc,xc-r1,xc+r1 nsel,r,y,loc,yc-r1,yc+r1 nsel,r,z,loc,2*hth,3*hth bf,all,hgen,%gene1% allsel,all *endif *elseif,i,ge,0.3,then allsel,all bfdele,all,hgen *endif *endif *endif allsel,all solve *enddo !計(jì)算結(jié)束 finish !************************************************* !下面是移動工作平面進(jìn)行剖切，方便觀察內(nèi)部溫度變化過程 wpoff,len/2,wid/2,0 wprot,0,0,90 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有269人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
2樓2013-05-07 08:42:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hjb1109 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 初試 317 +6 | 半拉月丙 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點(diǎn)0 2026-03-17 | 9/450 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|