 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1871 | 回復(fù): 7 | |||
| 【懸賞金幣】回答本帖問題，作者深邃呀將贈(zèng)送您 20 個(gè)金幣 | |||
深邃呀新蟲 (初入文壇)
|
[求助]
旋轉(zhuǎn)高斯曲面體熱源以及移動(dòng)熱源及加載和后處理 已有1人參與
|
||
|
激光切割仿真，自己根據(jù)網(wǎng)上資料寫了大部分命令流，遇到問題，體熱源如何加載，bfe選取路徑附近的節(jié)點(diǎn)加在嘛，求解設(shè)置中，子步如何理解，附上現(xiàn)有命令流，還望大佬解答解答 /com,units:SI,m-degree c-W-kg-sec-J /PREP7 ET,1,SOLID70 !定義熱分析單元 KEYOPT,1,2,1 ET,2,SURF152 !定義結(jié)構(gòu)的表面效應(yīng)單元 KEYOPT,2,4,1 KEYOPT,2,5,0 KEYOPT,2,8,5 !設(shè)置表面效應(yīng)單元關(guān)鍵字！ !設(shè)置材料參數(shù) MPTEMP,,,,,,,, MPTEMP,1,20 MPTEMP,2,250 MPTEMP,3,500 MPTEMP,4,750 MPTEMP,5,1000 MPTEMP,6,1500 MPTEMP,7,1700 MPTEMP,8,2500 MPDATA,KXX,1,,16 MPDATA,KXX,1,,19 MPDATA,KXX,1,,22 MPDATA,KXX,1,,23 MPDATA,KXX,1,,26 MPDATA,KXX,1,,27 MPDATA,KXX,1,,28 MPDATA,KXX,1,,30 MPDATA,DENS,1,,7524 MPDATA,DENS,1,,7314 MPDATA,DENS,1,,7254 MPDATA,DENS,1,,7236 MPDATA,DENS,1,,7199 MPDATA,DENS,1,,7054 MPDATA,DENS,1,,6978 MPDATA,DENS,1,,6954 MPDATA,PRXY,1,,0.28 MPDATA,PRXY,1,,0.29 MPDATA,PRXY,1,,0.31 MPDATA,PRXY,1,,0.35 MPDATA,PRXY,1,,0.4 MPDATA,PRXY,1,,0.45 MPDATA,PRXY,1,,0.48 MPDATA,PRXY,1,,0.5 MPDATA,C,1,,530 MPDATA,C,1,,674 MPDATA,C,1,,689 MPDATA,C,1,,720 MPDATA,C,1,,786 MPDATA,C,1,,810 MPDATA,C,1,,830 MPDATA,C,1,,854 MPDATA,EX,1,,0.193e12 MPDATA,EX,1,,0.165e12 MPDATA,EX,1,,0.132e12 MPDATA,EX,1,,0.5e11 MPDATA,EX,1,,0.32e11 MPDATA,EX,1,,12 MPDATA,EX,1,,6 MPDATA,EX,1,,2 UIMP,1,REFT,,,20 MPDATA,ALPX,1,,0.172e-4 MPDATA,ALPX,1,,0.181e-4 MPDATA,ALPX,1,,0.184e-4 MPDATA,ALPX,1,,0.194e-4 MPDATA,ALPX,1,,0.201e-4 MPDATA,ALPX,1,,0.209e-4 MPDATA,ALPX,1,,0.212e-4 MPDATA,ALPX,1,,0.215e-4 *SET,Length,0.05 !工件長度 *SET,Width,0.05 !工件寬度 *SET,Height,0.03 !工件高度 *SET,Widthglad,0.005 !熱源中心區(qū)域1/2 !建立幾何模型 k,1,Widthglad,0,0 k,2,Width/2,0,0 k,3,Width/2,0,-Height k,4,Widthglad,0,-Height k,5,-Widthglad,0,0 k,6,-Width/2,0,0 k,7,-Width/2,0,-Height k,8,-Widthglad,0,-Height K,9,0,0,0 k,10,0,0,-Height a,1,2,3,4 a,5,6,7,8 A,1,4,10,9 A,5,8,10,9 L,9,10 /pnum,area,1 !顯示面編號(hào) AGLUE,ALL !面粘結(jié)類似命令流LGLUE VGLUE VOFFST,1,Length VOFFST,2,-Length VOFFST,3,Length VOFFST,4,-Length !拉面成體 NUMMRG,ALL !對(duì)所有實(shí)體進(jìn)行合并 NUMCMP,ALL !壓縮所定義項(xiàng)的編號(hào)，命令流較長時(shí)，使用 可對(duì)未使用實(shí)體的編號(hào)進(jìn)行壓縮 SAVE Asel,s,loc,z,0,0 !選擇實(shí)體上表面 APLOT AATT,1,,2,0 !指定實(shí)體上表面用152單元 ESIZE,0.0001 MSHAPE,0,3D MSHKEY,1 AMESH,ALL !對(duì)上表面進(jìn)行劃分 ALLSEL,ALL !選擇所有實(shí)體 VATT,1,,1,0 !指定實(shí)體單元用solid70單元 ESIZE,0.0001 !設(shè)置單元尺寸 MSHAPE,0,3D !指定劃分單元的形狀為六面體單元 MSHKEY,1 !映射劃分 VMESH,ALL /SOL antype,4 !瞬態(tài)分析 trnopt,full !瞬態(tài)分析選項(xiàng)（一般full） tunif,20 !對(duì)所有節(jié)點(diǎn)指定初始均布溫度20℃ TREF,20 !設(shè)置參考溫度 SAVE !接下來定義對(duì)流換熱邊界 esel,s,type,,2 !選擇表面效應(yīng)單元設(shè)置強(qiáng)制對(duì)流3000 nsle,s sf,all,conv,3000 !其他表面設(shè)置空氣對(duì)流30 SFA,8,1,CONV,30,20 SFA,16,1,CONV,30,20 SFA,20,1,CONV,30,20 SFA,13,1,CONV,30,20 SFA,10,1,CONV,30,20 SFA,19,1,CONV,30,20 SFA,15,1,CONV,30,20 SFA,5,1,CONV,30,20 SFA,1,1,CONV,30,20 SFA,3,1,CONV,30,20 SFA,4,1,CONV,30,20 SFA,2,1,CONV,30,20 SFA,12,1,CONV,30,20 SFA,7,1,CONV,30,20 V=1/60 !切割速度 time=Length/v !切割時(shí)間 n=0.6 !吸收效率 R=0.0002 !光斑半徑 Qm=20000*n !功率 H=0.0003 !熱源深度 !旋轉(zhuǎn)高斯曲面體熱源函數(shù) *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNCCSYS *SET,_FNCNAME,'GS2' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *SET,_FNC_C1(1),12000 *SET,_FNC_C2(1),0.0002 *SET,_FNC_C3(1),0.0003 *SET,_FNCCSYS,0 ! /INPUT,tst.func,,,1 *DIM,%_FNCNAME%,TABLE,6,41,1,,,,%_FNCCSYS% ! ! Begin of equation: (9*Qm/(3.1415926*R0^2*H*(1-1/2.718281828459^3)))*exp((-9* ! ({X}^2+({Y}-(1/60)*{TIME})^2))/(R0^2*ln(H/{Z}))) *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), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), 0.0 *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 9, 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.1415926, 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, -4, 3, 19 *SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 2.718281828459, 0, 0, 0 *SET,%_FNCNAME%(0,9,1), 0.0, -4, 0, 3, 0, 0, -3 *SET,%_FNCNAME%(0,10,1), 0.0, -5, 0, 1, -3, 17, -4 *SET,%_FNCNAME%(0,11,1), 0.0, -3, 0, 1, 0, 0, -5 *SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 1, -3, 4, -5 *SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, 0, 0, -4 *SET,%_FNCNAME%(0,14,1), 0.0, -5, 0, 1, -3, 2, -4 *SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 1, -1, 3, -5 *SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 1, -2, 4, -3 *SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, 0, 0, -2 *SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, -2, 2, -3 *SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 9, 0, 0, -4 *SET,%_FNCNAME%(0,21,1), 0.0, -3, 0, 1, -4, 3, -2 *SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 1, 2, 17, -2 *SET,%_FNCNAME%(0,24,1), 0.0, -2, 0, 1, 0, 0, 0 *SET,%_FNCNAME%(0,25,1), 0.0, -5, 0, 60, 0, 0, -2 *SET,%_FNCNAME%(0,26,1), 0.0, -6, 0, 1, -2, 4, -5 *SET,%_FNCNAME%(0,27,1), 0.0, -2, 0, 1, -6, 3, 1 *SET,%_FNCNAME%(0,28,1), 0.0, -5, 0, 1, 3, 2, -2 *SET,%_FNCNAME%(0,29,1), 0.0, -2, 0, 2, 0, 0, -5 *SET,%_FNCNAME%(0,30,1), 0.0, -6, 0, 1, -5, 17, -2 *SET,%_FNCNAME%(0,31,1), 0.0, -2, 0, 1, -4, 1, -6 *SET,%_FNCNAME%(0,32,1), 0.0, -4, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 2, 0, 0, 18 *SET,%_FNCNAME%(0,34,1), 0.0, -3, 0, 1, 18, 17, -2 *SET,%_FNCNAME%(0,35,1), 0.0, -2, 0, 1, 19, 4, 4 *SET,%_FNCNAME%(0,36,1), 0.0, -2, 5, 1, -2, 0, 0 *SET,%_FNCNAME%(0,37,1), 0.0, -5, 0, 1, -3, 3, -2 *SET,%_FNCNAME%(0,38,1), 0.0, -2, 0, 1, -4, 4, -5 *SET,%_FNCNAME%(0,39,1), 0.0, -2, 7, 1, -2, 0, 0 *SET,%_FNCNAME%(0,40,1), 0.0, -3, 0, 1, -1, 3, -2 *SET,%_FNCNAME%(0,41,1), 0.0, 99, 0, 1, -3, 0, 0 ! End of equation: (9*Qm/(3.1415926*R0^2*H*(1-1/2.718281828459^3)))*exp((-9* ! ({X}^2+({Y}-(1/60)*{TIME})^2))/(R0^2*ln(H/{Z}))) !--> !!存在加載問題，不清楚如何加載 設(shè)想移動(dòng)軌跡為沿y軸的直線 !接下來熱源的加載以及分析設(shè)置 !!NSEL,S,LOC,x,Widthglad,-Widthglad !選擇X方向位于切縫熱源中心的節(jié)點(diǎn) !!NPLOT BFE,1,HGEN,,%GS2% !施加體熱源 AUTOTS,on !自動(dòng)時(shí)間步長 nsubst,50,50,50 KBC,1 solve bfedele,all,hgen !*********************************** !冷卻階段 !*********************************** *DO,k,1,30,1 TIME,time+5*k !載荷步結(jié)束時(shí)間 NSUBST,25 KBC,1 SOLVE SAVE *ENDDO |
回復(fù)此樓» 猜你喜歡
085600 材料與化工 295 求調(diào)劑
已經(jīng)有11人回復(fù)
-
NSFC申報(bào)書里申請(qǐng)人簡歷中代表性論著還需要在申報(bào)書最后的附件里面再上傳一遍嗎
已經(jīng)有11人回復(fù)
-
紡織、生物、化學(xué)、材料相關(guān)專業(yè)招生了
已經(jīng)有6人回復(fù)
-
材料與化工085600調(diào)劑求老師收留
已經(jīng)有7人回復(fù)
-
面上和青基一樣限30頁不合理
已經(jīng)有5人回復(fù)
-
308 085701 四六級(jí)已過求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿河海大學(xué)085900土木水利專碩279求調(diào)劑不挑專業(yè)
已經(jīng)有5人回復(fù)
-
290求調(diào)劑
已經(jīng)有8人回復(fù)
-
伙伴們，祝我生日快樂吧
已經(jīng)有9人回復(fù)
-
281求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
深邃呀
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 80.8
- 帖子: 28
- 在線: 4.5小時(shí)
- 蟲號(hào): 31128173
- 注冊(cè): 2022-10-08
- 專業(yè): 機(jī)械工程
2樓2023-03-05 20:31:46
【答案】應(yīng)助回帖
|
簡單看了一下你的命令流，感覺是你加熱過程有問題。好像你的加熱過程只有一個(gè)載荷步，熱源中心的位置是不變的，就是說你的熱源是固定的，沒有什么移動(dòng)。 正確的做法是：將熱源的整個(gè)移動(dòng)路徑劃分為若干小段，一個(gè)小段作為一個(gè)載荷步。對(duì)每個(gè)小段，先將熱源中心移動(dòng)到該小段的中心處，然后對(duì)結(jié)構(gòu)施加熱源函數(shù)定義的熱載荷�？紤]到熱源中心是移動(dòng)的，而函數(shù)中使用的是統(tǒng)一的坐標(biāo)系編號(hào) (不過你使用的是 *SET,_FNCCSYS,0，即總體直角坐標(biāo)系，這不太合適)，因此對(duì)每個(gè)載荷步需要定義一個(gè)局部坐標(biāo)系，其原點(diǎn)在熱源中心處。且對(duì)每個(gè)載荷步所定義的局部坐標(biāo)系編號(hào)相同，如 11，然后將函數(shù)中的函數(shù)坐標(biāo)系編號(hào)修改為該局部坐標(biāo)系的編號(hào)，如： *SET,_FNCCSYS,11。這樣才能保證隨著熱源的移動(dòng)，受熱的單元也隨之變化。 至于劃分多少載荷步比較合適，要看你的具體情況。雖然是越多越好，但是載荷步太多也會(huì)造成計(jì)算時(shí)間過長�？梢韵仁褂幂^少的載荷步，試算一下，成功后再增加載荷步，得到比較滿意的結(jié)果后可以不再增加載荷步。 |
» 本帖已獲得的紅花（最新10朵）
3樓2023-03-20 09:09:23
深邃呀
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 80.8
- 帖子: 28
- 在線: 4.5小時(shí)
- 蟲號(hào): 31128173
- 注冊(cè): 2022-10-08
- 專業(yè): 機(jī)械工程
送紅花一朵 |
意思是 我在模型外表面建立一個(gè)工件坐標(biāo)系，熱源函數(shù)的FUNCTIONS在此坐標(biāo)系下生成。至于移動(dòng)的話，我在函數(shù)表達(dá)式中引入了速度變量，根據(jù)仿真時(shí)間子步應(yīng)該會(huì)有移動(dòng)。就是不理解熱生成率怎么施加。一直比較困惑。不過感謝大佬的回答，能否加個(gè)好友，有些問題一直不好理解 發(fā)自小木蟲Android客戶端 |
4樓2023-03-20 12:25:25
【答案】應(yīng)助回帖
|
你對(duì)函數(shù)的作用可能有一點(diǎn)誤解。由于你的函數(shù)中使用了 函數(shù)坐標(biāo)值 (某個(gè)局部坐標(biāo)系)，執(zhí)行函數(shù)時(shí)，使用的是最靠近函數(shù)定義時(shí)的函數(shù)坐標(biāo)系，且在計(jì)算函數(shù)時(shí)這些值不會(huì)發(fā)生變化；除非你在函數(shù)之外修改了局部坐標(biāo)系的位置，然后重新調(diào)用函數(shù)。因此，即使你的函數(shù)定義中包含了速度，相應(yīng)的速度值也只是參加函數(shù)的計(jì)算，并不會(huì)改變函數(shù)坐標(biāo)系，即不會(huì)使熱源的位置發(fā)生改變 (通常，熱源中心會(huì)設(shè)置在函數(shù)坐標(biāo)系的原點(diǎn))。所以，要想使熱源發(fā)生移動(dòng)，必須在函數(shù)定義之前，人為的改變 函數(shù)坐標(biāo)系 (局部坐標(biāo)系) 的設(shè)置，才能使熱源發(fā)生移動(dòng)。 至于施加生成熱載荷，你的命令流中有一句： BFE,1,HGEN,,%GS2% !施加體熱源 確實(shí)是可以施加熱生成載荷的，只是你的用法不對(duì)。你可以看一下幫助文件，關(guān)于 bfe 命令的格式是： BFE, ELEM, Lab, STLOC, VAL1, VAL2, VAL3, VAL4 Defines an element body force load. 其中緊跟在命令 BFE 后面的參數(shù)是 elem，即施加提載荷的單元 (編號(hào))，也可以是 all - 所有單元，或若干單元組成的 component 。而你使用的是1，則表示只對(duì)第一個(gè)單元 (編號(hào)為 1) 施加熱生成載荷。這顯然是不對(duì)的。這里，你需要先根據(jù)熱源中心位置選擇施加熱生成載荷的單元，再將這些單元定義一個(gè) component，比如名為 e-hgen，然后使用如下命令對(duì)這些單元施加熱生成載荷： BFE, e-hgen, HGEN,,%GS2% !施加體熱源 你試試看吧。 |
» 本帖已獲得的紅花（最新10朵）
5樓2023-03-21 08:46:27
深邃呀
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 80.8
- 帖子: 28
- 在線: 4.5小時(shí)
- 蟲號(hào): 31128173
- 注冊(cè): 2022-10-08
- 專業(yè): 機(jī)械工程
6樓2023-03-21 11:06:49
【答案】應(yīng)助回帖
|
將： BFE,1,HGEN,,%GS2% !施加體熱源 改為： BFE,all,HGEN,,%GS2% 是加大了加熱單元的范圍，但也有一個(gè)問題：這里的 all 是指整個(gè)結(jié)構(gòu)的全部單元嗎？好像一般情況應(yīng)該只是受到熱源直接照射的表層的單元才會(huì)受熱，非表層的單元不會(huì)直接受到熱源的照射，而是通過熱傳導(dǎo)受熱的吧？ 因此，這里需要增加對(duì)表層單元的選擇，然后使用 BFE,all,HGEN,,%GS2% 則代表只對(duì)選中的 (表層) 施加熱載荷。施加熱載荷后，再執(zhí)行命令 allsel,all 選擇所有單元，然后進(jìn)行求解。 再說一點(diǎn)，對(duì)于求解設(shè)置： AUTOTS,on !自動(dòng)時(shí)間步長 nsubst,50,50,50 KBC,1 建議修改如下： AUTOTS,on !自動(dòng)時(shí)間步長 nsubst,10,500,5 KBC,0 有利于加快收斂速度。 |
» 本帖已獲得的紅花（最新10朵）
7樓2023-03-22 08:05:39
深邃呀
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 80.8
- 帖子: 28
- 在線: 4.5小時(shí)
- 蟲號(hào): 31128173
- 注冊(cè): 2022-10-08
- 專業(yè): 機(jī)械工程
8樓2023-03-22 22:59:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 深邃呀 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +10 | dream…… 2026-03-10 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑（085602一志愿985） +9 | 化工人999 2026-03-09 | 9/450 |
|
|
|
[考研] 0856材料與化工309分求調(diào)劑 +5 | ZyZy…… 2026-03-10 | 5/250 |
|
|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
 10+5
|
穿只靴子 2026-03-07 | 22/1100 |
|
|
|
[考博] 讀博申請(qǐng) +5 | 感dd 2026-03-10 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +4 | 鶴遨予卿 2026-03-09 | 4/200 |
|
|
|
[考研] 考研材料與化工，求調(diào)劑 +7 | 戲精丹丹丹 2026-03-09 | 7/350 |
|
|
|
[考研] 420求調(diào)劑 +3 | 莫向外求11 2026-03-10 | 3/150 |
|
|
|
[考研] 材料工程專碩，一志愿中國礦業(yè)大學(xué)，總分314，求調(diào)劑 +4 | 無懈可擊的巨人 2026-03-10 | 4/200 |
|
|
|
[考研] 0856材料與化工353分求調(diào)劑 +11 | NIFFFfff 2026-03-09 | 11/550 |
|
|
|
[考研] 一志愿天津大學(xué)，英一數(shù)二305分求調(diào)劑，四六級(jí)已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
|
[考研] 求調(diào)劑！ +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+4
|
簡木ChuFront 2026-03-09 | 4/200 |
|
|
|
[考研] 337求調(diào)劑 +3 | 睡醒，。 2026-03-09 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +11 | 朝天椒的雙馬尾 2026-03-05 | 11/550 |
|
|
|
[教師之家] 交大前校長王樹國：現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室，而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|