 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1421 | 回復(fù): 1 | ||
| 【懸賞金幣】回答本帖問題，作者科研古惑仔將贈送您 5 個金幣 | ||
科研古惑仔新蟲 (初入文壇)
|
[求助]
仿真冷卻時間添加 已有1人參與
|
|
|
請教大神�。。∥沂亲鲞x區(qū)激光融化熱物理仿真的，仿真過程不知道怎么在一層結(jié)束后添加冷卻，希望大佬們解答一下，不勝感激�。。�！下面是我的代碼 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ù)實體編號選擇 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 !自動控制（默認(rèn)參數(shù)，可以不寫） nropt,full,,on !設(shè)置牛頓-拉普森方法（默認(rèn)參數(shù)，可以不寫） cutcontrol,plslimit,10, !設(shè)置在一個線性求解中時間步的減少（作用不明確，建議寫上） lnsrch,on !線性搜索（對收斂有一定影響） autots,on !自動步長（對收斂有一定影響） tintp,,,,1,, !向后積分 tref,298.15 !參考溫度 kbc,0 !連續(xù)性載荷 neqit,50 !規(guī)定每個子步中最大迭代次數(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é)點初始溫度設(shè)定25 nsel,s,ext !選擇所有外部節(jié)點 sf,all,conv,%conve%,298.15 !在節(jié)點上施加面載荷，外部所有節(jié)點初始298K的對流系數(shù) time,1e-6 !設(shè)置時間間隔，為載荷步設(shè)置時間 deltim,1e-6,1e-7,1e-6 !設(shè)置時間步長 kbc,1 !載荷階躍方式加載 allsel,all !選中所有項目 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 !在每個單元所指定的面上施加熱源載荷 *endif *endif *enddo !*********************************************************************************** time,tm !時間間隔 !deltim,0.01,0.01,0.2 !時間步長 !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 !刪除每個單元六個面上的熱載荷 *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é)點載xyz三個方向上的位移全部約束掉 !***********************設(shè)置非線性求解選項（與溫度場大同小異）************* 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)力計算****************************************** 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ù)此樓» 猜你喜歡
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有28人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有40人回復(fù)
-
070300調(diào)劑一志愿南京大學(xué)化學(xué)專業(yè)，初試273分，有論文有專利有競賽，求調(diào)劑B區(qū)211
已經(jīng)有11人回復(fù)
-
中科大材料299求調(diào)劑
已經(jīng)有10人回復(fù)
-
08工科求調(diào)劑
已經(jīng)有3人回復(fù)
-
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有12人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有12人回復(fù)
-
085602高分子方向求調(diào)劑
已經(jīng)有7人回復(fù)
-
進(jìn)入個人成果庫好難，一下午都沒進(jìn)去
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有11人回復(fù)
2樓2020-11-19 08:40:56
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 科研古惑仔 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過 a 線，總分 281 求調(diào)劑 +5 | 香香 Lu 2026-03-03 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼？~+123456 2026-03-05 | 9/450 |
|
|
[考研] 求調(diào)劑 +5 | danyyyy 2026-03-04 | 5/250 |
|
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計算機(jī)、人工智能等專業(yè)的同學(xué)報考 +10 | pin8023 2026-02-28 | 13/650 |
|
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
|
[基金申請] 沒有青基直接申請面上，感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
|
[考研] 一志愿鄭州大學(xué)，學(xué)碩，物理化學(xué)， 333求調(diào)劑 +5 | 李魔女斗篷 2026-03-04 | 5/250 |
|
|
|
[考研] 325求調(diào)劑 +5 | 學(xué)家科 2026-03-04 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研]
材料工程專碩283求調(diào)劑
 5+8
|
，��？ 2026-03-02 | 10/500 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
|
[考研] 一志愿中石油（華東）本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
|
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
|
[考研]
 尋找調(diào)劑
+4
|
LYidhsjabdj 2026-02-28 | 4/200 |
|