 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1347 | 回復(fù): 3 | |||
[求助]
求助高斯光源加載不上去，初學(xué)者。 已有2人參與
|
|
我做的是選擇性激光燒結(jié)溫度場(chǎng)的仿真，材料為聚苯乙烯，改了一下別人的程序后死活熱源顯示加載了，但是求解中就沒有了。 /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àng)**************************************** /solu antype,4 trnopt,full !完全瞬態(tài)積分法 pred,on !打開預(yù)測(cè)校正 solcontrol,on !自動(dòng)控制 nropt,full,,on !完全牛頓-拉普森方法 !cnvtol,u,,0.5 !設(shè)置收斂值，位移 !cnvtol,f,,0.05 !設(shè)置收斂值，力 cutcontrol,plslimit,10,!設(shè)置在一個(gè)非線性求解中時(shí)間步的減少 !eqslv,sparse,,2 !求解器類型 lnsrch,on !線性搜索 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ù)此樓» 猜你喜歡
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有32人回復(fù)
-
0856材料與化工353分求調(diào)劑
已經(jīng)有4人回復(fù)
-
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有13人回復(fù)
-
315食品工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境調(diào)劑
已經(jīng)有3人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有20人回復(fù)
-
考研材料與化工，求調(diào)劑
已經(jīng)有3人回復(fù)
-
歡迎211本科同學(xué)，過A區(qū)國(guó)家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有35人回復(fù)
-
304求調(diào)劑（085602一志愿985）
已經(jīng)有4人回復(fù)
-
294 英二數(shù)二物化 求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
klzg77541812
專家顧問 (正式寫手)
NVH-FEA工程師
- 應(yīng)助: 117 (高中生)
- 金幣: 2141.8
- 紅花: 28
- 帖子: 621
- 在線: 606.2小時(shí)
- 蟲號(hào): 1479383
- 注冊(cè): 2011-11-06
- 性別: GG
- 專業(yè): 固體力學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
!1定義局部坐標(biāo)系 2 加載命令放在循環(huán)語(yǔ)句外面，因?yàn)榧虞d的是函數(shù)，沒有必要重復(fù)加 3 熱源處網(wǎng)格太少了，加密網(wǎng)格 4 計(jì)算總時(shí)間太短，可能你還沒調(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àng)**************************************** /solu antype,4 trnopt,full !完全瞬態(tài)積分法 pred,on !打開預(yù)測(cè)校正 solcontrol,on !自動(dòng)控制 nropt,full,,on !完全牛頓-拉普森方法 !cnvtol,u,,0.5 !設(shè)置收斂值，位移 !cnvtol,f,,0.05 !設(shè)置收斂值，力 cutcontrol,plslimit,10,!設(shè)置在一個(gè)非線性求解中時(shí)間步的減少 !eqslv,sparse,,2 !求解器類型 lnsrch,on !線性搜索 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
專家顧問 (正式寫手)
NVH-FEA工程師
- 應(yīng)助: 117 (高中生)
- 金幣: 2141.8
- 紅花: 28
- 帖子: 621
- 在線: 606.2小時(shí)
- 蟲號(hào): 1479383
- 注冊(cè): 2011-11-06
- 性別: GG
- 專業(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ā)表 | |
|---|---|---|---|---|
|
[考研] 標(biāo)題:撿漏預(yù)警|08工科/09農(nóng)學(xué)調(diào)劑!英語(yǔ)要求低，過線即有機(jī)會(huì)! +4 | 馬超放煙花 2026-03-07 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿天大化工（085600）調(diào)劑總分338 +4 | 蔡大美女 2026-03-09 | 4/200 |
|
|
|
[考博] 2026博士申請(qǐng) +5 | 起泡酒 2026-03-08 | 5/250 |
|
|
|
[考研] 310 070300化學(xué)求調(diào)劑 +4 | 撲風(fēng)鈴的貓 2026-03-08 | 5/250 |
|
|
|
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
|
[考研] 招調(diào)劑 +20 | 帆船哥 2026-03-04 | 23/1150 |
|
|
|
[考研] 一志愿武漢理工085601，初試301分，請(qǐng)問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
|
[考研] 化學(xué)專業(yè)調(diào)劑 +5 | 好好好1233 2026-03-04 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
|
[考研] 一志愿清華深研院材料專碩294分，專業(yè)課111分，本科中南大學(xué)材料，有六級(jí)，有工作經(jīng)驗(yàn) +3 | H14528 2026-03-04 | 3/150 |
|
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|