 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1855 | 回復: 7 | ||
| 【懸賞金幣】回答本帖問題，作者深邃呀將贈送您 20 個金幣 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
旋轉(zhuǎn)高斯曲面體熱源以及移動熱源及加載和后處理 已有1人參與
|
||
|
激光切割仿真，自己根據(jù)網(wǎng)上資料寫了大部分命令流，遇到問題，體熱源如何加載，bfe選取路徑附近的節(jié)點加在嘛，求解設(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)的表面效應單元 KEYOPT,2,4,1 KEYOPT,2,5,0 KEYOPT,2,8,5 !設(shè)置表面效應單元關(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 !顯示面編號 AGLUE,ALL !面粘結(jié)類似命令流LGLUE VGLUE VOFFST,1,Length VOFFST,2,-Length VOFFST,3,Length VOFFST,4,-Length !拉面成體 NUMMRG,ALL !對所有實體進行合并 NUMCMP,ALL !壓縮所定義項的編號，命令流較長時，使用 可對未使用實體的編號進行壓縮 SAVE Asel,s,loc,z,0,0 !選擇實體上表面 APLOT AATT,1,,2,0 !指定實體上表面用152單元 ESIZE,0.0001 MSHAPE,0,3D MSHKEY,1 AMESH,ALL !對上表面進行劃分 ALLSEL,ALL !選擇所有實體 VATT,1,,1,0 !指定實體單元用solid70單元 ESIZE,0.0001 !設(shè)置單元尺寸 MSHAPE,0,3D !指定劃分單元的形狀為六面體單元 MSHKEY,1 !映射劃分 VMESH,ALL /SOL antype,4 !瞬態(tài)分析 trnopt,full !瞬態(tài)分析選項（一般full） tunif,20 !對所有節(jié)點指定初始均布溫度20℃ TREF,20 !設(shè)置參考溫度 SAVE !接下來定義對流換熱邊界 esel,s,type,,2 !選擇表面效應單元設(shè)置強制對流3000 nsle,s sf,all,conv,3000 !其他表面設(shè)置空氣對流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 !切割時間 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è)想移動軌跡為沿y軸的直線 !接下來熱源的加載以及分析設(shè)置 !!NSEL,S,LOC,x,Widthglad,-Widthglad !選擇X方向位于切縫熱源中心的節(jié)點 !!NPLOT BFE,1,HGEN,,%GS2% !施加體熱源 AUTOTS,on !自動時間步長 nsubst,50,50,50 KBC,1 solve bfedele,all,hgen !*********************************** !冷卻階段 !*********************************** *DO,k,1,30,1 TIME,time+5*k !載荷步結(jié)束時間 NSUBST,25 KBC,1 SOLVE SAVE *ENDDO |
回復此樓» 猜你喜歡
26申博自薦
已經(jīng)有7人回復
-
進入個人成果庫好難，一下午都沒進去
已經(jīng)有5人回復
-
接受26屆調(diào)劑生
已經(jīng)有12人回復
-
求調(diào)劑
已經(jīng)有5人回復
-
求調(diào)劑 一志愿蘇州大學，0856化工323分 | 本科應化 | 有專利/競賽/科研助手經(jīng)歷 |
已經(jīng)有6人回復
-
求調(diào)劑
已經(jīng)有4人回復
-
308求調(diào)劑
已經(jīng)有9人回復
-
303求調(diào)劑
已經(jīng)有7人回復
-
085600材料與化工 292分 求調(diào)劑
已經(jīng)有4人回復
-
一志愿中石油，化學工程，301求調(diào)劑
已經(jīng)有4人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
【答案】應助回帖
|
將： BFE,1,HGEN,,%GS2% !施加體熱源 改為： BFE,all,HGEN,,%GS2% 是加大了加熱單元的范圍，但也有一個問題：這里的 all 是指整個結(jié)構(gòu)的全部單元嗎？好像一般情況應該只是受到熱源直接照射的表層的單元才會受熱，非表層的單元不會直接受到熱源的照射，而是通過熱傳導受熱的吧？ 因此，這里需要增加對表層單元的選擇，然后使用 BFE,all,HGEN,,%GS2% 則代表只對選中的 (表層) 施加熱載荷。施加熱載荷后，再執(zhí)行命令 allsel,all 選擇所有單元，然后進行求解。 再說一點，對于求解設(shè)置： AUTOTS,on !自動時間步長 nsubst,50,50,50 KBC,1 建議修改如下： AUTOTS,on !自動時間步長 nsubst,10,500,5 KBC,0 有利于加快收斂速度。 |
7樓2023-03-22 08:05:39
2樓2023-03-05 20:31:46
【答案】應助回帖
|
簡單看了一下你的命令流，感覺是你加熱過程有問題。好像你的加熱過程只有一個載荷步，熱源中心的位置是不變的，就是說你的熱源是固定的，沒有什么移動。 正確的做法是：將熱源的整個移動路徑劃分為若干小段，一個小段作為一個載荷步。對每個小段，先將熱源中心移動到該小段的中心處，然后對結(jié)構(gòu)施加熱源函數(shù)定義的熱載荷�？紤]到熱源中心是移動的，而函數(shù)中使用的是統(tǒng)一的坐標系編號 (不過你使用的是 *SET,_FNCCSYS,0，即總體直角坐標系，這不太合適)，因此對每個載荷步需要定義一個局部坐標系，其原點在熱源中心處。且對每個載荷步所定義的局部坐標系編號相同，如 11，然后將函數(shù)中的函數(shù)坐標系編號修改為該局部坐標系的編號，如： *SET,_FNCCSYS,11。這樣才能保證隨著熱源的移動，受熱的單元也隨之變化。 至于劃分多少載荷步比較合適，要看你的具體情況。雖然是越多越好，但是載荷步太多也會造成計算時間過長�？梢韵仁褂幂^少的載荷步，試算一下，成功后再增加載荷步，得到比較滿意的結(jié)果后可以不再增加載荷步。 |
» 本帖已獲得的紅花（最新10朵）
3樓2023-03-20 09:09:23
送紅花一朵 |
意思是 我在模型外表面建立一個工件坐標系，熱源函數(shù)的FUNCTIONS在此坐標系下生成。至于移動的話，我在函數(shù)表達式中引入了速度變量，根據(jù)仿真時間子步應該會有移動。就是不理解熱生成率怎么施加。一直比較困惑。不過感謝大佬的回答，能否加個好友，有些問題一直不好理解 發(fā)自小木蟲Android客戶端 |
4樓2023-03-20 12:25:25
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3 | ipys00zm58 2026-03-06 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 哈爾濱理工大學2026年研究生調(diào)劑，材料科學與化學工程學院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
|
[考研] 復試調(diào)劑 +7 | 呼呼？~+123456 2026-03-05 | 10/500 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +6 | cgyzqwn 2026-03-01 | 13/650 |
|
|
|
[考研] 考研282分求調(diào)劑，接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
|
[基金申請] 請問大家，研究風險與應對措施那里， 大家都怎么寫呢 ？ +3 | cauasen 2026-03-02 | 3/150 |
|
|
|
[考研] 理學，工學，農(nóng)學調(diào)劑，少走彎路，這里歡迎您！ +8 | likeihood 2026-03-02 | 11/550 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
|
[考研] 1 +7 | 黑！在干嘛 2026-02-28 | 8/400 |
|
|
|
[考研] 0854復試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 材料學調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|