 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1341 | 回復(fù): 3 | ||
[求助]
求助高斯光源加載不上去，初學(xué)者。 已有2人參與
|
|
我做的是選擇性激光燒結(jié)溫度場(chǎng)的仿真，材料為聚苯乙烯，改了一下別人的程序后死活熱源顯示加載了，但是求解中就沒(méi)有了。 /PREP7 !* !======================================================== =======================================!指定單元 ET,1,SOLID70 !* !* !======================================================== =======================================!材料屬性 MPTEMP,,,,,,,, MPTEMP,1,0 MPTEMP,2,25 MPTEMP,3,100 MPTEMP,4,165 MPTEMP,5,227 MPTEMP,6,327 MPDATA,KXX,1,,0.041 MPDATA,KXX,1,,0.047 MPDATA,KXX,1,,0.048 MPDATA,KXX,1,,0.049 MPDATA,KXX,1,,0.05 MPDATA,KXX,1,,0.05 MPTEMP,, MPTEMP,1,0 MPTEMP,2,25 MPTEMP,3,100 MPTEMP,4,165 MPTEMP,5,227 MPTEMP,6,327 MPDATA,C,1,,2810 MPDATA,C,1,,2904 MPDATA,C,1,,3046 MPDATA,C,1,,3612 MPDATA,C,1,,3821 MPDATA,C,1,,3912 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,472 *dim,conve,table,8,1,1,temp !總換熱系數(shù)(W/m^2*C) conve(0,1,1)=1 conve(1,0,1)=0 conve(1,1,1)=12.12e-3 conve(2,0,1)=200 conve(2,1,1)=24.23e-3 conve(3,0,1)=400 conve(3,1,1)=40.38e-3 conve(4,0,1)=600 conve(4,1,1)=72.69e-3 conve(5,0,1)=800 conve(5,1,1)=125e-3 conve(6,0,1)=1000 conve(6,1,1)=150e-3 conve(7,0,1)=1500 conve(7,1,1)=200e-3 conve(8,0,1)=2000 conve(8,1,1)=250e-3 !======================================================== =======================================!定義常量 Width_Base=0.025 !基底寬度 Height_Base=0.002 !基底高度 Length_Base=0.05 !基底長(zhǎng)度 Length_Clad=0.05 !熔覆層長(zhǎng)度 Height_Clad=0.0025 !熔覆層高度 Width_Clad=0.025 !覆層總寬度 !======================================================== =======================================!定義常量 BLOCK,0,Length_Clad,-Width_Clad/2,Width_Clad/2,0,Height_Clad, !覆層 BLOCK,0,Length_Base,-Width_Base/2,Width_Base/2,-Height_Base,0, !基底 esize,0.001 vmesh,all SFA,3,1,CONV,%conve%,90 SFA,4,1,CONV,%conve%,90 SFA,5,1,CONV,%conve%,90 SFA,6,1,CONV,%conve%,90 SFA,7,1,CONV,%conve%,90 SFA,8,1,CONV,%conve%,90 SFA,9,1,CONV,%conve%,90 SFA,10,1,CONV,%conve%,90 SFA,11,1,CONV,%conve%,90 V=2 !掃描速度 P=20 !激光功率 r=0.00015 !激光光斑半徑 Area=3.14159265*(r**2) !激光光斑面積 Qm=2*P/Area StepDis=0.001 !每個(gè)載荷步位移 StepTime=StepDis/V !每個(gè)載荷步時(shí)間 StepNum=(Length_Clad)/StepDis !載荷步數(shù) !************************************非線(xiàn)性求解選項(xiàng)**************************************** /solu antype,4 trnopt,full !完全瞬態(tài)積分法 pred,on !打開(kāi)預(yù)測(cè)校正 solcontrol,on !自動(dòng)控制 nropt,full,,on !完全牛頓-拉普森方法 !cnvtol,u,,0.5 !設(shè)置收斂值，位移 !cnvtol,f,,0.05 !設(shè)置收斂值，力 cutcontrol,plslimit,10,!設(shè)置在一個(gè)非線(xiàn)性求解中時(shí)間步的減少 !eqslv,sparse,,2 !求解器類(lèi)型 lnsrch,on !線(xiàn)性搜索 autots,on !自動(dòng)步長(zhǎng) tintp,,,,1,, !向后積分 timint,on !瞬態(tài)效果 tref,90 !參考溫度 kbc,0 !連續(xù)性載荷 neqit,50 !規(guī)定每個(gè)子步中最大迭代次數(shù),默認(rèn)為25 !======================================================== =======================================!高斯單元 !**** 高斯光束能量公式定義（GUI方式自動(dòng)生成）**** !****GAUSS=2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)**** *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNCCSYS *SET,_FNCNAME,'GAOSI' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *SET,_FNC_C1(1),P ! *SET,_FNC_C2(1),R ! *SET,_FNCCSYS,11 !將方程坐標(biāo)系設(shè)定為自己定義的局部坐標(biāo)系 ! /INPUT,GAUSS.func,,,1 *DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS% ! ! Begin of equation: 2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), 0.0 *SET,%_FNCNAME%(6,0,1), 0.0 *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,3,1), 0, -1, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, 18, 17, -1 *SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 3, -3 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 4, -4 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, -4 *SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -4, 3, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 2, 17, -2 *SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, 3, 17, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 18, 17, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -3 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 7, 1, -2, 0, 0 *SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation:2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2) !--> *DO, k, 1, StepNum, 1 SFA,2,1,HFLUX, %GAOSI% !激活單元的上表面,加熱流密度 OUTRES,ALL,ALL, TIME,StepTime*k !載荷步結(jié)束時(shí)間 NSUBST, 1 !指定載荷子步數(shù)目 DELTIM,0.1,0.1,0.1 !載荷子步 solve SFADELE,all,HFLUX !刪除熱流密度載荷 |
回復(fù)此樓» 猜你喜歡
288求調(diào)劑
已經(jīng)有5人回復(fù)
-
0703求調(diào)劑
已經(jīng)有9人回復(fù)
-
083000環(huán)境科學(xué)與工程調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
085701環(huán)境工程專(zhuān)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料與化工290求調(diào)劑
已經(jīng)有5人回復(fù)
-
求材料讀博院校
已經(jīng)有3人回復(fù)
-
347求調(diào)劑
已經(jīng)有4人回復(fù)
-
2026考研求調(diào)劑-材料類(lèi)-本科211一志愿985-初試301分
已經(jīng)有5人回復(fù)
-
求調(diào)劑，一志愿華中科大0702，數(shù)一英一，293
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
klzg77541812
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
NVH-FEA工程師
- 應(yīng)助: 117 (高中生)
- 金幣: 2141.8
- 紅花: 28
- 帖子: 621
- 在線(xiàn): 606.2小時(shí)
- 蟲(chóng)號(hào): 1479383
- 注冊(cè): 2011-11-06
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
!1定義局部坐標(biāo)系 2 加載命令放在循環(huán)語(yǔ)句外面，因?yàn)榧虞d的是函數(shù)，沒(méi)有必要重復(fù)加 3 熱源處網(wǎng)格太少了，加密網(wǎng)格 4 計(jì)算總時(shí)間太短，可能你還沒(méi)調(diào)整 5 熱流密度的大小不知道對(duì)不對(duì)，顯示加載數(shù)值很小0.12 /PREP7 !* !======================================================== =======================================!指定單元 ET,1,SOLID70 !* !* !======================================================== =======================================!材料屬性 MPTEMP,,,,,,,, MPTEMP,1,0 MPTEMP,2,25 MPTEMP,3,100 MPTEMP,4,165 MPTEMP,5,227 MPTEMP,6,327 MPDATA,KXX,1,,0.041 MPDATA,KXX,1,,0.047 MPDATA,KXX,1,,0.048 MPDATA,KXX,1,,0.049 MPDATA,KXX,1,,0.05 MPDATA,KXX,1,,0.05 MPTEMP,, MPTEMP,1,0 MPTEMP,2,25 MPTEMP,3,100 MPTEMP,4,165 MPTEMP,5,227 MPTEMP,6,327 MPDATA,C,1,,2810 MPDATA,C,1,,2904 MPDATA,C,1,,3046 MPDATA,C,1,,3612 MPDATA,C,1,,3821 MPDATA,C,1,,3912 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,472 *dim,conve,table,8,1,1,temp !總換熱系數(shù)(W/m^2*C) conve(0,1,1)=1 conve(1,0,1)=0 conve(1,1,1)=12.12e-3 conve(2,0,1)=200 conve(2,1,1)=24.23e-3 conve(3,0,1)=400 conve(3,1,1)=40.38e-3 conve(4,0,1)=600 conve(4,1,1)=72.69e-3 conve(5,0,1)=800 conve(5,1,1)=125e-3 conve(6,0,1)=1000 conve(6,1,1)=150e-3 conve(7,0,1)=1500 conve(7,1,1)=200e-3 conve(8,0,1)=2000 conve(8,1,1)=250e-3 !======================================================== =======================================!定義常量 Width_Base=0.025 !基底寬度 Height_Base=0.002 !基底高度 Length_Base=0.05 !基底長(zhǎng)度 Length_Clad=0.05 !熔覆層長(zhǎng)度 Height_Clad=0.0025 !熔覆層高度 Width_Clad=0.025 !覆層總寬度 !======================================================== =======================================!定義常量 BLOCK,0,Length_Clad,-Width_Clad/2,Width_Clad/2,0,Height_Clad, !覆層 BLOCK,0,Length_Base,-Width_Base/2,Width_Base/2,-Height_Base,0, !基底 esize,0.001 vmesh,all SFA,3,1,CONV,%conve%,90 SFA,4,1,CONV,%conve%,90 SFA,5,1,CONV,%conve%,90 SFA,6,1,CONV,%conve%,90 SFA,7,1,CONV,%conve%,90 SFA,8,1,CONV,%conve%,90 SFA,9,1,CONV,%conve%,90 SFA,10,1,CONV,%conve%,90 SFA,11,1,CONV,%conve%,90 V=2 !掃描速度 P=20 !激光功率 r=0.00015 !激光光斑半徑 Area=3.14159265*(r**2) !激光光斑面積 Qm=2*P/Area StepDis=0.001 !每個(gè)載荷步位移 StepTime=StepDis/V !每個(gè)載荷步時(shí)間 StepNum=(Length_Clad)/StepDis !載荷步數(shù) !************************************非線(xiàn)性求解選項(xiàng)**************************************** /solu antype,4 trnopt,full !完全瞬態(tài)積分法 pred,on !打開(kāi)預(yù)測(cè)校正 solcontrol,on !自動(dòng)控制 nropt,full,,on !完全牛頓-拉普森方法 !cnvtol,u,,0.5 !設(shè)置收斂值，位移 !cnvtol,f,,0.05 !設(shè)置收斂值，力 cutcontrol,plslimit,10,!設(shè)置在一個(gè)非線(xiàn)性求解中時(shí)間步的減少 !eqslv,sparse,,2 !求解器類(lèi)型 lnsrch,on !線(xiàn)性搜索 autots,on !自動(dòng)步長(zhǎng) tintp,,,,1,, !向后積分 timint,on !瞬態(tài)效果 tref,90 !參考溫度 kbc,0 !連續(xù)性載荷 neqit,50 !規(guī)定每個(gè)子步中最大迭代次數(shù),默認(rèn)為25 CSWPLA,11,0,1,1, !定義局部坐標(biāo)系11 !======================================================== =======================================!高斯單元 !**** 高斯光束能量公式定義（GUI方式自動(dòng)生成）**** !****GAUSS=2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2)**** *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNCCSYS *SET,_FNCNAME,'GAOSI' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *SET,_FNC_C1(1),P ! *SET,_FNC_C2(1),R ! *SET,_FNCCSYS,11 !將方程坐標(biāo)系設(shè)定為自己定義的局部坐標(biāo)系 ! /INPUT,GAUSS.func,,,1 *DIM,%_FNCNAME%,TABLE,6,24,1,,,,%_FNCCSYS% ! ! Begin of equation: 2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), 0.0 *SET,%_FNCNAME%(6,0,1), 0.0 *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,3,1), 0, -1, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, 18, 17, -1 *SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3 *SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 3, -3 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 4, -4 *SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, -4 *SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, -4, 3, -2 *SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, 2, 17, -2 *SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 2, 0, 0, 3 *SET,%_FNCNAME%(0,16,1), 0.0, -5, 0, 1, 3, 17, -2 *SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -4, 1, -5 *SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, 18, 17, -2 *SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -3 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 7, 1, -2, 0, 0 *SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,24,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation:2*P/(PI*R^2)*exp(-2*({X}-V*{TIME})^2)+{Y}^2)/R^2) !--> SFA,2,1,HFLUX, %GAOSI% !激活單元的上表面,加熱流密度 *DO, k, 1, StepNum, 1 OUTRES,ALL,ALL, TIME,StepTime*k !載荷步結(jié)束時(shí)間 NSUBST, 1 !指定載荷子步數(shù)目 DELTIM,0.1,0.1,0.1 !載荷子步 solve |
2樓2019-07-03 11:51:17
klzg77541812
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
NVH-FEA工程師
- 應(yīng)助: 117 (高中生)
- 金幣: 2141.8
- 紅花: 28
- 帖子: 621
- 在線(xiàn): 606.2小時(shí)
- 蟲(chóng)號(hào): 1479383
- 注冊(cè): 2011-11-06
- 性別: GG
- 專(zhuān)業(yè): 固體力學(xué)
- 管轄: 計(jì)算模擬區(qū)
3樓2019-07-03 13:01:52
4樓2019-07-03 22:37:03
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 桂~~~ 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑，一志愿華中科大0702，數(shù)一英一，293 +4 | 小羅露一二 2026-03-07 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工求調(diào)劑 +3 | 與冬清寧 2026-03-07 | 5/250 |
|
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 6/300 |
|
|
[考研] 考研求調(diào)劑，工科，總分304 +3 | 矮子不想長(zhǎng)高了 2026-03-08 | 4/200 |
|
|
|
[考研] 083000，總分284，求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
|
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷 +4 | 熊境喆 2026-03-05 | 4/200 |
|
|
|
[考研] 08工科 +4 | li李樂(lè)成 2026-03-06 | 4/200 |
|
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +9 | llx0610 2026-03-02 | 10/500 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +13 | February_Feb 2026-03-01 | 13/650 |
|
|
|
[考研] 不限學(xué)校專(zhuān)業(yè)的調(diào)劑同學(xué)看過(guò)來(lái) +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
|
[考博] 2026申博自薦 六級(jí)440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
|
[考研] 080500材料科學(xué)與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|
|
[考研] 考研282分求調(diào)劑，接受跨專(zhuān)業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
|
[考研] 282求調(diào)劑 +7 | 夕～日 2026-03-05 | 8/400 |
|
|
|
[考研] 一志愿西交化工專(zhuān)碩288專(zhuān)業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
|
[基金申請(qǐng)] 沒(méi)有青基直接申請(qǐng)面上，感覺(jué)自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
|
[考研] 295求調(diào)劑 +6 | 等春來(lái)， 2026-03-04 | 6/300 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|