| 5 | 1/1 | 返回列表 |
| 查看: 4950 | 回復: 11 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
yaojinnan金蟲 (正式寫手)
|
[交流]
【求助】ANSYS 中命令流編程(fortran) 已有11人參與
|
|||
|
大家好,我是做焊接模擬的,用的是ANSYS軟件,F(xiàn)在想做一段程序,其中想做個循環(huán)交替加載。 比如:載荷作用時間總共為10s,有兩個載荷。 其中,0-1s 給選定物體加載第一個載荷, 1-3s給選定物體加載第二個載荷 3-4s給選定物體加載第一個載荷 5-7s給選定物體加載第二個載荷 7-8s給選定物體加載第一個載荷 8-10s給選定物體加載第二個載荷 其實簡單的可以這么理解。 但是我的載荷是一個由ANSYS函數(shù)編輯器定義的一個表載荷, 并且這兩個載荷都是移動的。所以不清楚應該怎么用FORTRAN來編程。希望大家能給予幫助。定有回報。 下面就是我的命令流程序: finish /clear,start !清除數(shù)據(jù)庫以開始新的分 /filname,reverse-shuang tuo qiu,0 /title,hanhanahanhan /units,si /prep7 ET,1,SOLID70 ET,2,SOLID90 !et,2,surf152 !keyopt,2,4,0 !keyopt,2,5,1 !keyopt,2,8,3 MPTEMP,1,20,100,200,300,400,500, MPTEMP,7,600,700,800,900,1000,1500 MPDATA,KXX,1,1,152,159,164,165,168,196, MPDATA,KXX,1,7,209,162,162,162,162,162 MPDATA,DENS,1,1,2680,2662,2646,2620,2600,2580, MPDATA,DENS,1,7,2436,2384,2384,2384,2384,2384 MPDATA,C,1,1,956,963,1047,1130,1224,1308, MPDATA,C,1,7,1412,1084,1084,1084,1084,1084 MPDATA,HF,1,1,8.22,11.0,13.7,23.2,33.4,46.8, MPDATA,HF,1,7,58.0,68.5,68.5,68.5,68.5,68.5 lx=0.08 ly=0.2 lz=0.005 r=0.002 D=0.002!焊縫寬度 ah=D/2!雙橢球熱源參數(shù)定義 chf=4*ah chb=6*ah bh=0.003 rf=2*chf/(chf+chb) rb=2*chb/(chf+chb) K=0.5 功率有效系數(shù) I=130 U=30 q0=I*U q=q0*K !有效功率 pi=acos(-1) step=0.002 v=0.003 !count=ly/v count=20 !*********************************** !定義單元類型 !*********************************** !建立幾何模型 !*********************************** wpstyle,,,,,,,,0 block,0,0.006,0,ly,0,-lz block,0,0.012,0,ly,0,-lz block,0,0.024,0,ly,0,-lz block,0,lx,0,ly,0,-lz vovlap,all /view,1,1,1,1 !定義焊縫區(qū)單元尺寸0.002 lesize,7,0.001 lesize,2,0.001 lesize,12,0.001 lesize,11,0.001 lesize,6,0.001 lesize,3,0.001 lesize,8,0.001 lesize,1,0.001 lesize,4,0.001 lesize,9,0.001 lesize,5,0.001 lesize,10,0.001 !定義過渡區(qū)的單元尺寸 lesize,15,r lesize,18,r !定義遠離焊縫區(qū)的單元尺寸 lesize,34,1.25*0.001 lesize,35,1.25*0.001 lesize,47,1.25*0.001 lesize,46,1.25*0.001 lesize,39,4*r lesize,42,4*r lesize,57,4*r lesize,60,4*r lesize,27,4*r lesize,30,4*r lesize,58,4*r lesize,59,4*r !*********************************** !劃分網(wǎng)格,生生有限元模型 !*********************************** type,1 mshkey,1 mshape,0 vmesh,1 vmesh,7 type,2 mshkey,0 mshape,1 vmesh,5 vmesh,6 numcmp,area save !*********************************** !循環(huán)加載、求解 !*********************************** !施加載荷 /solu *do,i,0,count disy=i*step t=(disy+step)/v antype,4 trnopt,full tunif,293 outres,basic,last time,t nropt,full,,on autots,1 nsubst,1,,,1 kbc,1 lnsrch,1 tsres,erase !TINTP,,,,1,,,,, timint,1,therm !*********************************** !焊接過程開始 !*********************************** nsel,s,ext, nsel,u,loc,x,0 !------改 sf,all,conv,200,293 !------改 allsel,all !*********************************** !*********************************** !雙橢球熱源的加載 !*********************************** !前半球節(jié)點的選取 !*********************************** LOCAL,11,2,0,disy,0.000,,,,chf/ah,bh/ah, NSEL,S,LOC,X,0,ah NSEL,R,LOC,Y,0,90 NSEL,R,LOC,Z,-90,0 CM,FRONT,NODE LOCAL,13,0,0,disy,0.000,,,,1,1, *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNCCSYS *DEL,'TFRONT' *SET,_FNCNAME,'TFRONT' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *SET,_FNC_C1(1),q *SET,_FNC_C2(1),rf *SET,_FNC_C3(1),ah *SET,_FNC_C4(1),chf *SET,_FNC_C5(1),bh *SET,_FNCCSYS,13 *DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS% %_FNCNAME%(0,0,1)= 0.0, -999 %_FNCNAME%(2,0,1)= 0.0 %_FNCNAME%(3,0,1)= %_FNC_C1(1)% %_FNCNAME%(4,0,1)= %_FNC_C2(1)% %_FNCNAME%(5,0,1)= %_FNC_C3(1)% %_FNCNAME%(6,0,1)= %_FNC_C4(1)% %_FNCNAME%(7,0,1)= %_FNC_C5(1)% %_FNCNAME%(0,1,1)= 1.0, -1, 16, 3, 0, 0, 0 %_FNCNAME%(0,2,1)= 0.0, -2, 0, 6, 0, 0, -1 %_FNCNAME%(0,3,1)= 0, -3, 0, 1, -2, 3, -1 %_FNCNAME%(0,4,1)= 0.0, -1, 0, 1, -3, 3, 17 %_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -1, 3, 18 %_FNCNAME%(0,6,1)= 0.0, -1, 0, 0, 0, 0, 0 %_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, 0, 0, -1 %_FNCNAME%(0,8,1)= 0.0, -4, 0, 1, -1, 2, -3 %_FNCNAME%(0,9,1)= 0.0, -1, 0, 3, 0, 0, -4 %_FNCNAME%(0,10,1)= 0.0, -3, 0, 1, -4, 3, -1 %_FNCNAME%(0,11,1)= 0.0, -1, 0, 2, 0, 0, 2 %_FNCNAME%(0,12,1)= 0.0, -4, 0, 1, 2, 17, -1 %_FNCNAME%(0,13,1)= 0.0, -1, 0, 2, 0, 0, 19 %_FNCNAME%(0,14,1)= 0.0, -5, 0, 1, 19, 17, -1 %_FNCNAME%(0,15,1)= 0.0, -1, 0, 1, -4, 4, -5 %_FNCNAME%(0,16,1)= 0.0, -4, 0, 2, 0, 0, 3 %_FNCNAME%(0,17,1)= 0.0, -5, 0, 1, 3, 17, -4 %_FNCNAME%(0,18,1)= 0.0, -4, 0, 2, 0, 0, 20 %_FNCNAME%(0,19,1)= 0.0, -6, 0, 1, 20, 17, -4 %_FNCNAME%(0,20,1)= 0.0, -4, 0, 1, -5, 4, -6 %_FNCNAME%(0,21,1)= 0.0, -5, 0, 1, -1, 1, -4 %_FNCNAME%(0,22,1)= 0.0, -1, 0, 2, 0, 0, 4 %_FNCNAME%(0,23,1)= 0.0, -4, 0, 1, 4, 17, -1 %_FNCNAME%(0,24,1)= 0.0, -1, 0, 2, 0, 0, 21 %_FNCNAME%(0,25,1)= 0.0, -6, 0, 1, 21, 17, -1 %_FNCNAME%(0,26,1)= 0.0, -1, 0, 1, -4, 4, -6 %_FNCNAME%(0,27,1)= 0.0, -4, 0, 1, -5, 1, -1 %_FNCNAME%(0,28,1)= 0.0, -1, 0, 1, -3, 3, -4 %_FNCNAME%(0,29,1)= 0.0, -1, 7, 1, -1, 0, 0 %_FNCNAME%(0,30,1)= 0.0, -3, 0, 1, -2, 3, -1 %_FNCNAME%(0,31,1)= 0.0, -1, 16, 3.14159265358979310, 0, 0, 0 %_FNCNAME%(0,32,1)= 0.0, -2, 0, 3.14159265358979310, 0, 0, -1 %_FNCNAME%(0,33,1)= 0.0, -4, 0, 1, -2, 3, -1 %_FNCNAME%(0,34,1)= 0.0, -1, 0, 1, -4, 3, 20 %_FNCNAME%(0,35,1)= 0.0, -2, 0, 1, -1, 3, 19 %_FNCNAME%(0,36,1)= 0.0, -1, 0, 1, -2, 3, 21 %_FNCNAME%(0,37,1)= 0.0, -2, 0, 1, -3, 4, -1 %_FNCNAME%(0,38,1)= 0.0, 99, 0, 1, -2, 0, 0 BF,FRONT,HGEN,%TFRONT% CMDELE,FRONT ALLSEL,ALL !!!后半球節(jié)點的選取 LOCAL,12,2,0,disy,0,,,,chb/ah,bh/ah, NSEL,S,LOC,X,0,ah NSEL,R,LOC,Y,-90,0 nsel,r,loc,y,-90,0 CM,BEHIND,NODE !!! LOCAL,14,0,0,disy,0,,,,1,1, *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNCCSYS *DEL,'TBEHIND' *SET,_FNCNAME,'TBEHIND' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *SET,_FNC_C1(1),q *SET,_FNC_C2(1),rb *SET,_FNC_C3(1),ah *SET,_FNC_C4(1),chb *SET,_FNC_C5(1),bh *SET,_FNCCSYS,14 *DIM,%_FNCNAME%,TABLE,7,38,1,,,,%_FNCCSYS% %_FNCNAME%(0,0,1)= 0.0, -999 %_FNCNAME%(2,0,1)= 0.0 %_FNCNAME%(3,0,1)= %_FNC_C1(1)% %_FNCNAME%(4,0,1)= %_FNC_C2(1)% %_FNCNAME%(5,0,1)= %_FNC_C3(1)% %_FNCNAME%(6,0,1)= %_FNC_C4(1)% %_FNCNAME%(7,0,1)= %_FNC_C5(1)% %_FNCNAME%(0,1,1)= 1.0, -1, 16, 3, 0, 0, 0 %_FNCNAME%(0,2,1)= 0.0, -2, 0, 6, 0, 0, -1 %_FNCNAME%(0,3,1)= 0, -3, 0, 1, -2, 3, -1 %_FNCNAME%(0,4,1)= 0.0, -1, 0, 1, -3, 3, 17 %_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -1, 3, 18 %_FNCNAME%(0,6,1)= 0.0, -1, 0, 0, 0, 0, 0 %_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, 0, 0, -1 %_FNCNAME%(0,8,1)= 0.0, -4, 0, 1, -1, 2, -3 %_FNCNAME%(0,9,1)= 0.0, -1, 0, 3, 0, 0, -4 %_FNCNAME%(0,10,1)= 0.0, -3, 0, 1, -4, 3, -1 %_FNCNAME%(0,11,1)= 0.0, -1, 0, 2, 0, 0, 2 %_FNCNAME%(0,12,1)= 0.0, -4, 0, 1, 2, 17, -1 %_FNCNAME%(0,13,1)= 0.0, -1, 0, 2, 0, 0, 19 %_FNCNAME%(0,14,1)= 0.0, -5, 0, 1, 19, 17, -1 %_FNCNAME%(0,15,1)= 0.0, -1, 0, 1, -4, 4, -5 %_FNCNAME%(0,16,1)= 0.0, -4, 0, 2, 0, 0, 3 %_FNCNAME%(0,17,1)= 0.0, -5, 0, 1, 3, 17, -4 %_FNCNAME%(0,18,1)= 0.0, -4, 0, 2, 0, 0, 20 %_FNCNAME%(0,19,1)= 0.0, -6, 0, 1, 20, 17, -4 %_FNCNAME%(0,20,1)= 0.0, -4, 0, 1, -5, 4, -6 %_FNCNAME%(0,21,1)= 0.0, -5, 0, 1, -1, 1, -4 %_FNCNAME%(0,22,1)= 0.0, -1, 0, 2, 0, 0, 4 %_FNCNAME%(0,23,1)= 0.0, -4, 0, 1, 4, 17, -1 %_FNCNAME%(0,24,1)= 0.0, -1, 0, 2, 0, 0, 21 %_FNCNAME%(0,25,1)= 0.0, -6, 0, 1, 21, 17, -1 %_FNCNAME%(0,26,1)= 0.0, -1, 0, 1, -4, 4, -6 %_FNCNAME%(0,27,1)= 0.0, -4, 0, 1, -5, 1, -1 %_FNCNAME%(0,28,1)= 0.0, -1, 0, 1, -3, 3, -4 %_FNCNAME%(0,29,1)= 0.0, -1, 7, 1, -1, 0, 0 %_FNCNAME%(0,30,1)= 0.0, -3, 0, 1, -2, 3, -1 %_FNCNAME%(0,31,1)= 0.0, -1, 16, 3.14159265358979310, 0, 0, 0 %_FNCNAME%(0,32,1)= 0.0, -2, 0, 3.14159265358979310, 0, 0, -1 %_FNCNAME%(0,33,1)= 0.0, -4, 0, 1, -2, 3, -1 %_FNCNAME%(0,34,1)= 0.0, -1, 0, 1, -4, 3, 20 %_FNCNAME%(0,35,1)= 0.0, -2, 0, 1, -1, 3, 19 %_FNCNAME%(0,36,1)= 0.0, -1, 0, 1, -2, 3, 21 %_FNCNAME%(0,37,1)= 0.0, -2, 0, 1, -3, 4, -1 %_FNCNAME%(0,38,1)= 0.0, 99, 0, 1, -2, 0, 0 BF,BEHIND,HGEN,%TBEHIND% CMDELE,BEHIND ALLSEL,ALL SOLVE BFDELE,all,HGEN *ENDDO SAVE 其中從雙橢球熱源加載到最后就是加載的過程,這是一個載荷,還需要另外一個載荷,這里沒有加載,就像上面我所說的那樣循環(huán)交替那樣應該怎么加載,希望大家給予幫助。。應該用到怎么樣的語句呢?? |
碩博精華 |


至尊木蟲 (著名寫手)
新蟲 (初入文壇)
| 我是做焊接的,想學Ansys模擬焊接過程,有這方面的教程嗎?若有請傳我一份,郵箱:leixiaokuan_0512@1363.com |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 材料專碩326求調劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 293求調劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
[考研] 274求調劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
[考研] 08工科 320總分 求調劑 +4 | 梨花珞晚風 2026-03-17 | 4/200 |
|
|
[考研] 267一志愿南京工業(yè)大學0817化工求調劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] 環(huán)境工程調劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
[考研] 318求調劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 265求調劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 331求調劑(0703有機化學 +5 | ZY-05 2026-03-13 | 6/300 |
|
|
[考研] 材料工程調劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 333求調劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] 281求調劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 0856化學工程280分求調劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考博] 福州大學楊黃浩課題組招收2026年專業(yè)學位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
[考研] 283求調劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|