 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1422 | 回復(fù): 1 | |||
| 【懸賞金幣】回答本帖問題，作者科研古惑仔將贈(zèng)送您 5 個(gè)金幣 | |||
科研古惑仔新蟲 (初入文壇)
|
[求助]
仿真冷卻時(shí)間添加 已有1人參與
|
||
|
請教大神�。。∥沂亲鲞x區(qū)激光融化熱物理仿真的，仿真過程不知道怎么在一層結(jié)束后添加冷卻，希望大佬們解答一下，不勝感激�。。。∠旅媸俏业拇a finish /clear /units,si !m-K-kg-s-J-W-Pa /nerr,100000 /UIS,MSGPOP,3 /prep7 et,1,solid70 !******************************定義熱分析物理參數(shù)************************************ MPTEMP,,,,,,,, MPTEMP,1,298.15 MPTEMP,2,673.15 MPTEMP,3,1073.15 MPTEMP,4,1473.15 MPTEMP,5,1658.15 MPTEMP,6,1723.15 MPTEMP,7,1873.15 MPTEMP,8,2973.15 MPDATA,DENS,1,,7950 MPDATA,DENS,1,,7785 MPDATA,DENS,1,,7575 MPDATA,DENS,1,,7361 MPDATA,DENS,1,,7275 MPDATA,DENS,1,,7236 MPDATA,DENS,1,,6765 MPDATA,DENS,1,,6740 MPTEMP,,,,,,,, MPTEMP,1,298.15 MPTEMP,2,673.15 MPTEMP,3,1073.15 MPTEMP,4,1473.15 MPTEMP,5,1658.15 MPTEMP,6,1723.15 MPTEMP,7,1873.15 MPTEMP,8,2973.15 MPDATA,C,1,,470 MPDATA,C,1,,560 MPDATA,C,1,,630 MPDATA,C,1,,700 MPDATA,C,1,,730 MPDATA,C,1,,830 MPDATA,C,1,,830 MPDATA,C,1,,960 MPTEMP,,,,,,,, MPTEMP,1,298.15 MPTEMP,2,673.15 MPTEMP,3,1073.15 MPTEMP,4,1473.15 MPTEMP,5,1658.15 MPTEMP,6,1723.15 MPTEMP,7,1873.15 MPTEMP,8,2973.15 MPDATA,KXX,1,,13.4 MPDATA,KXX,1,,21.8 MPDATA,KXX,1,,27.2 MPDATA,KXX,1,,30.9 MPDATA,KXX,1,,32.4 MPDATA,KXX,1,,28.5 MPDATA,KXX,1,,30.5 MPDATA,KXX,1,,32.2 MPTEMP,,,,,,,, MPTEMP,1,298.15 MPTEMP,2,673.15 MPTEMP,3,1073.15 MPTEMP,4,1473.15 MPTEMP,5,1658.15 MPTEMP,6,1723.15 MPTEMP,7,1873.15 MPTEMP,8,2973.15 MPDE,ENTH,1 MPDATA,ENTH,1,,298.15*470*7950 MPDATA,ENTH,1,,673.15*560*7785 MPDATA,ENTH,1,,1073.15*630*7575 MPDATA,ENTH,1,,1473.15*7361*700 MPDATA,ENTH,1,,1658.15*730*7275 MPDATA,ENTH,1,,1723.15*830*7236 MPDATA,ENTH,1,,1873.15*830*6765 MPDATA,ENTH,1,,2973.15*960*6740 *dim,conve,table,12,1,1,temp conve(0,1,1)=1 conve(1,0,1)=298.15 conve(2,0,1)=473.15 conve(3,0,1)=673.15 conve(4,0,1)=873.15 conve(5,0,1)=1073.15 conve(6,0,1)=1273.15 conve(7,0,1)=1773.15 conve(8,0,1)=2073.15 conve(9,0,1)=2273.15 conve(10,0,1)=2473.15 conve(11,0,1)=2673.15 conve(12,0,1)=2873.15 conve(1,1,1)=12.12 conve(2,1,1)=24.23 conve(3,1,1)=40.38 conve(4,1,1)=72.69 conve(5,1,1)=125.00 conve(6,1,1)=150.38 conve(7,1,1)=200.00 conve(8,1,1)=250.65 conve(9,1,1)=300.65 conve(10,1,1)=350.6 conve(11,1,1)=400.6 conve(12,1,1)=450.6 !******************建立模型********************** BLOCK,-0.0056,0.005,-0.0006,0.0012,0,-0.0003, BLOCK,-0.0043,0.005,-0.0006,0.0012,-0.0002,-0.0003, BLOCK,0,0.0006,0,0.0006,-0.0002,0.00003, VSBV,1,2 vovlap,all !體粘結(jié) numcmp,all,,,low !壓縮單元編號，重新編號 !******************劃分網(wǎng)格********************** allsel vsel,u,volu,,1,2 !u指不選，volu是根據(jù)實(shí)體編號選擇 wprota,,,90 !旋轉(zhuǎn)工作平面，指XOZ平面由Z至X為方向旋轉(zhuǎn)90度 wprota,,45 vsbw,all wpcsys,-1 !根據(jù)坐標(biāo)系，定義工作平面 wpoff,0.0006 !平移工作平面 wprota,,,90 wprota,,-45 vsbw,all wpcsys,-1 wprota,,,90 wpoff,,,-0.0006 vsbw,all wpoff,,,0.0012+0.0006 vsbw,all wpcsys,-1 wpoff,,,-0.0002 vsbw,all wpcsys,-1 wpoff,-0.0043 wprota,,,90 vsbw,all lsel,s,lengh,,0.0006 !選擇長度為0.006mm的直線 lesize,all,0.00002 lsel,s,lengh,,0.0018 lesize,all,,,0.0006/0.00002 lsel,s,lengh,,0.0002 lesize,all,,,4 lsel,s,lengh,,0.0001 lesize,all,,,2 vsel,s,volu,,2 ! aslv,s,1 !已選擇體中選擇面 lsla,s,1 lesize,all,0.00003 Mshkey,1 Mshape,0 Vmesh,all allsel lesize,all,0.0003 esize,0.0003 Mshkey,1 Mshape,0 Vmesh,all v=1 p=225 radius=0.06e-3 ef=0.3 pi=3.1415926 qmax=2*ef*p/(pi*radius*radius) !***激光光斑中最大的熱流密度 /solu !表示在求解器中進(jìn)行設(shè)置 antype,4 !定義分析類型，4表示分析為瞬態(tài)分析。（可在后續(xù)操作中打開或關(guān)閉） trnopt,full !完全瞬態(tài)積分法（默認(rèn)參數(shù)，可以不寫） pred,on !打開預(yù)測校正（默認(rèn)參數(shù)，可以不寫） solcontrol,on !自動(dòng)控制（默認(rèn)參數(shù)，可以不寫） nropt,full,,on !設(shè)置牛頓-拉普森方法（默認(rèn)參數(shù)，可以不寫） cutcontrol,plslimit,10, !設(shè)置在一個(gè)線性求解中時(shí)間步的減少（作用不明確，建議寫上） lnsrch,on !線性搜索（對收斂有一定影響） autots,on !自動(dòng)步長（對收斂有一定影響） tintp,,,,1,, !向后積分 tref,298.15 !參考溫度 kbc,0 !連續(xù)性載荷 neqit,50 !規(guī)定每個(gè)子步中最大迭代次數(shù),默認(rèn)為25 outres,all !將輸出結(jié)果全部顯示出來 !***********************初始化載荷步參數(shù)********************************** tinc=0.06e-3/v !*******************************************穩(wěn)態(tài)初值分析********************************************* timint,off !瞬態(tài)效應(yīng)關(guān)閉 alls !選取所有幾何 ic,all,temp,298.15 ! 所有節(jié)點(diǎn)初始溫度設(shè)定25 nsel,s,ext !選擇所有外部節(jié)點(diǎn) sf,all,conv,%conve%,298.15 !在節(jié)點(diǎn)上施加面載荷，外部所有節(jié)點(diǎn)初始298K的對流系數(shù) time,1e-6 !設(shè)置時(shí)間間隔，為載荷步設(shè)置時(shí)間 deltim,1e-6,1e-7,1e-6 !設(shè)置時(shí)間步長 kbc,1 !載荷階躍方式加載 allsel,all !選中所有項(xiàng)目 solve xmin0=0 xmax0=0.06e-3 xmin=xmin0 xmax=xmax0 xd=0.06e-3 ymin0=0 ymax0=0.06e-3 ymin=ymin0 ymax=ymax0 yd=0.06e-3 ! zmin=0 zmax=0.03e-3 zd=0.03e-3 zz=0.03e-3 !施加熱源的表面 !*******************************************瞬態(tài)分析****************************************** /solu timint,on tm=tinc lnum=0.0006/0.06e-3 !!!! *do,n,1,10,1 *do,im,1,lnum,1 ! 從左往右掃描一道 →→→→→→→ esel,s,cent,y,ymin,ymax esel,r,cent,x,xmin,xmax esel,r,cent,z,zmin,zmax eplot *get,nemax,elem,,num,max !得到所有1類單元號碼的最大值 *get,nemin,elem,,num,min a=(xmin+xmax)/2 !v*tinc*im 熱源位置(a,b,c) b=v*tinc*im c=zz !***********************************以下施加熱流密度載荷************************** *do,i,nemin,nemax,1 *if,esel(i),eq,1,then !如果該單元在上面的集合里，則 xsy=centrx(i) !讀取該單元的中心坐標(biāo) ysy=centry(i) zsy=centrz(i) rr=abs(sqrt((xsy-a)*(xsy-a)+(ysy-b)*(ysy-b)+(zsy-c)*(zsy-c))) !該單元中心距離熱源中心 *if,rr,le,radius,And,zsy,ge,0,then !該單元中心在加熱半徑范圍內(nèi) qr=qmax*exp(-2*rr*rr/(radius*radius)) !該單元中心處的熱流大小 sfe,i,5,hflux,,qr !在每個(gè)單元所指定的面上施加熱源載荷 *endif *endif *enddo !*********************************************************************************** time,tm !時(shí)間間隔 !deltim,0.01,0.01,0.2 !時(shí)間步長 !deltim,tinc/20,tinc/30,tinc/10 deltim,4.5e-6,3e-6,9e-5 !deltim,2e-4,2e-4,2e-4 alls solve *do,i,nemin,nemax,1 *if,esel(i),eq,1,then !如果該單元在上面的集合里，則 xsy=centrx(i) !讀取該單元的中心坐標(biāo) ysy=centry(i) zsy=centrz(i) rr=abs(sqrt((xsy-a)*(xsy-a)+(ysy-b)*(ysy-b))) !該單元中心距離熱源中心 *if,rr,le,radius,And,zsy,ge,0,then !該單元中心在加熱半徑范圍內(nèi) sfedele,all,all,hflux !刪除每個(gè)單元六個(gè)面上的熱載荷 *endif *endif *enddo ymin=ymin+yd ymax=ymax+yd tm=tm+tinc *enddo xmin=xmin+yd xmax=xmax+yd ymin=ymin0 ymax=ymax0 *enddo ***********************************第二部分進(jìn)行應(yīng)力場分析****************************************** !******************************定義應(yīng)力分析參數(shù)****************************** /prep7 MPTEMP,,,,,,,, MPTEMP,1,298.15 MPTEMP,2,873.15 MPTEMP,3,1273.15 MPTEMP,4,1773.15 MPTEMP,5,2273.15 MPTEMP,6,2673.15 MPTEMP,7,3473.15 MPDATA,EX,1,,2E+11 MPDATA,EX,1,,1.35E+11 MPDATA,EX,1,,2E+10 MPDATA,EX,1,,2E+8 MPDATA,EX,1,,2E+7 MPDATA,EX,1,,2E+7 MPDATA,EX,1,,2E+7 MPDATA,PRXY,1,,0.33 MPDATA,PRXY,1,,0.35 MPDATA,PRXY,1,,0.36 MPDATA,PRXY,1,,0.40 MPDATA,PRXY,1,,0.40 MPDATA,PRXY,1,,0.41 MPDATA,PRXY,1,,0.45 MPTEMP,,,,,,,, MPTEMP,1,298.15 MPTEMP,2,873.15 MPTEMP,3,1273.15 MPTEMP,4,1773.15 MPTEMP,5,2273.15 MPTEMP,6,2673.15 MPTEMP,7,3473.15 UIMP,1,REFT,,, MPDATA,ALPX,1,,12E-6 MPDATA,ALPX,1,,13E-6 MPDATA,ALPX,1,,13.4E-6 MPDATA,ALPX,1,,14E-6 MPDATA,ALPX,1,,15E-6 MPDATA,ALPX,1,,15E-6 MPDATA,ALPX,1,,15E-6 TB,BKIN,1,6,2,1 TBTEMP,298.15 TBDATA,,3.45E+08,5.89E+09,,,, TBTEMP,873.15 TBDATA,,2.12E+08,1.70E+09,,,, TBTEMP,1073.15 TBDATA,,1.99E+08,1.40E+09,,,, TBTEMP,1473.15 TBDATA,,1.00E+08,1.00E+08,,,, TBTEMP,1623.15 TBDATA,,5.00E+07,1.00E+07,,,, TBTEMP,1773.15 TBDATA,,5.00E+06,1.00E+06,,,, /com,structural /prep7 /config,nres,300000 alls lsclear,all etchg,tts nsel,s,loc,z,-0.0003 d,all,ux,0,,,,uy,uz !在結(jié)構(gòu)分析中，把節(jié)點(diǎn)載xyz三個(gè)方向上的位移全部約束掉 !***********************設(shè)置非線性求解選項(xiàng)（與溫度場大同小異）************* FINISH /solu antype,4 trnopt,full !完全瞬態(tài)積分法 lumpm,0 !使用與單元相關(guān)的質(zhì)量矩陣公式 nlgeom,on pred,on sstif,1 solcontrol,on nropt,full,,on !完全頓-拉普森方法 !cnvtol,u,,0.5, l2, -1 !設(shè)置收斂值，位移 !cnvtol,f,,0.05, l2, -1 !設(shè)置收斂值，力 cnvtol,u,,0.3 cnvtol,f,,0.03 !設(shè)置收斂值,力 cutcontrol,plslimit,10, eqslv,sparse,,2 lnsrch,on pstres,1 autots,on tintp,,,,1,, timint,on tref,298.15 kbc,0 !連續(xù)性載荷 neqit,50 OUTRES,NSOL,LAST OUTPR,NSOL,LAST !***********************初始化載荷步參數(shù)********************************** tinc=0.06e-3/v !***********************穩(wěn)態(tài)分析****************************************** alls ldread,temp,1,,,,,rth !讀第一步的溫度 /pbf,temp,,1 /replot time,1e-6 deltim,1e-6,1e-7,1e-6 solve !***********************應(yīng)力計(jì)算****************************************** timint,on ! 打開瞬態(tài)分析 !*************************生死單元********************** esel,s,cent,z,0,0.00003 ekill,all /solu j=1 time=tinc xmin0=0 xmax0=0.06e-3 xmin=xmin0 xmax=xmax0 xd=0.06e-3 ymin0=0 ymax0=0.06e-3 ymin=ymin0 ymax=ymax0 yd=0.06e-3 ! zmin=0 zmax=0.03e-3 zd=0.03e-3 zz=0.03e-3 !施加熱源的表面 *do,n,1,10,1 *do,im,1,lnum,1 ! 從左往右掃描一道 →→→→→→→ esel,s,cent,y,ymin,ymax esel,r,cent,x,xmin,xmax esel,r,cent,z,zmin,zmax ealive,all eplot ldread,temp,j+1,,,,,rth! /pbf,temp,,1 /replot autots,on time,tm deltim,0.25*tinc,0.1*tinc,0.25*tinc alls solve ymin=ymin+xd ymax=ymax+xd tm=tm+tinc j=j+1 *enddo xmin=xmin+yd xmax=xmax+yd ymin=ymin0 ymax=ymax0 *enddo |
回復(fù)此樓» 猜你喜歡
材料調(diào)劑
已經(jīng)有13人回復(fù)
-
085600材料與化工 292分 求調(diào)劑
已經(jīng)有3人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有8人回復(fù)
-
273求調(diào)劑
已經(jīng)有7人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有45人回復(fù)
-
334求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑推薦
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿211，0860總分286食品工程專業(yè)求調(diào)劑
已經(jīng)有6人回復(fù)
-
招調(diào)劑
已經(jīng)有14人回復(fù)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
你是什么時(shí)候需要進(jìn)行冷卻？是在 n 循環(huán)中，還是在 im 循環(huán)中？ 假設(shè)你是在 im 循環(huán)中要進(jìn)行冷卻，那么你可以在每個(gè) im 的加熱過程結(jié)束后，添加一個(gè)冷卻過程。此時(shí)應(yīng)該刪除所有的加熱載荷，然后添加冷卻載荷 (必須有冷卻載荷，否則溫度不變。通常是對流載荷，如果前面已經(jīng)定義過，也可以不重復(fù)定義)，如定義求解選項(xiàng) (如終止時(shí)間，等)，進(jìn)行求解。 |
2樓2020-11-19 08:40:56
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 科研古惑仔 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +12 | 帆船哥 2026-03-04 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 材料專碩290求調(diào)劑 +8 | 杰尼龜aaa 2026-03-04 | 8/400 |
|
|
[考研] 301求調(diào)劑 +7 | 朝天椒的雙馬尾 2026-03-05 | 7/350 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +13 | February_Feb 2026-03-01 | 13/650 |
|
|
|
[考研] 274求調(diào)劑 +6 | cgyzqwn 2026-03-01 | 13/650 |
|
|
|
[考研] 一志愿華中科技大學(xué)，化學(xué)專業(yè)344分，求調(diào)劑 +7 | 邢xing1 2026-03-02 | 7/350 |
|
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
 5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
|
[考研] 中國科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑，有科研和獲獎(jiǎng)經(jīng)歷 +11 | wsxw 2026-03-02 | 12/600 |
|
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
|
[考研] 化學(xué)工程求調(diào)劑 +10 | 化工人999 2026-03-04 | 10/500 |
|
|
|
[考研] 材料復(fù)試調(diào)劑 +7 | 學(xué)材料的點(diǎn) 2026-03-01 | 8/400 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +6 | 等春來， 2026-03-04 | 6/300 |
|
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[考研] 11408，學(xué)碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|