 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1639 | 回復(fù): 1 | ||
hjb1109鐵蟲 (初入文壇)
|
[求助]
ANSYS用生死單元實現(xiàn)多層逐步加熱，加載中存在小問題，求解
|
|
現(xiàn)在做的工作是對一種算法的有效性演算，所以建模和加載都進行了簡化。 需要實現(xiàn)的功能是利用生死單元實現(xiàn)多層逐步加熱：就是先激活底層，在底層加載；底層加熱完成再激活上面的第二層，再進行加載；第二層加熱完成后再激活上面的第三層…… 加載用的熱源模型是錐形體熱源，就是高斯激光作用于半無限大物體表面時所產(chǎn)生的熱效應(yīng)。就是這里出了問題，加載第一層還能夠顯示出錐形加熱，但是到了第二層第三層時不再顯示出是錐形，而是圓柱形。在計算過程中ANSYS也給出了警告：在選取節(jié)點進行加載時不能識別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)備，三維表格每個維度的數(shù)量 max_x=len/lsize+1 max_y=wid/lsize+1 max_z=hth/lsize+1 !進入前處理器 /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è)置，對整個程序沒有很大影響，是計算實際需要，很長一段都可以不看 !接觸設(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 !進入求解器 /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所含單元殺死，只激活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時，加載第三層 *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 !計算結(jié)束 finish !************************************************* !下面是移動工作平面進行剖切，方便觀察內(nèi)部溫度變化過程 wpoff,len/2,wid/2,0 wprot,0,0,90 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有91人回復(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ù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有68人回復(fù)
-
華東師范大學(xué)芯片設(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ā)表 | |
|---|---|---|---|---|
|
[考研] 336求調(diào)劑 +3 | rmc8866 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
|
[考研] 354求調(diào)劑 +5 | Tyoumou 2026-03-18 | 8/400 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +3 | ￣^￣゜汗 2026-03-19 | 4/200 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|