移动体热源ANSYS命令流Word格式文档下载.docx

1、开始，进入前处理器，定义下列属性0参数定义1定义单元类型2实常数3材料性能4几何模型5网格划分控制6其他/PREP70参数定义 模型参数几何模型尺寸me-3=mm*SET,Block_x,*SET,Block_x1, *SET,Block_x2, *SET,Block_y,*SET,Block_y1,*SET,Block_y2,*SET,Block_z,*SET,Block_z1,划分网格的大小*SET,LMeshx,3e-4 !line1 *SET,LMeshy,3e-4 !line4*SET,LMeshz,3e-4 !line9边界条件参数*SET,Tinit,30 !初始温度ET,1,S

2、OLID70ET,2,surf152 !表面效应单元r,2KEYOPT,2,4,1 !表面效应单元设置 KEYOPT,2,8,4ET,3,SOLID87NothingMPTEMP MPTEMP,1,25,400,1200,1600,2000,2500 MPDATA,DENS,1,1,3983,6908,6331 !密度MPDATA,KXX,1,1,32 !热导率 MPDATA,c,1,1,472,535,679,859,945,1034 !比热MPTEMP, MPTEMP,1,0,200,400,800,1200 MPTEMP,6,1400,1454,1800,2000 MPDATA,ENTH

3、,1,1,0,+8,+9,+9,+9 MPDATA,ENTH,1,6,+9,+9,+9,+10 !热焓MPTEMP,MPTEMP,1,20,300,600,900,1200,1500,1800,2000,3000,4000MPDATA,HF,1,1,6,50,120,200,250,378,700,850,1746,3799 !对流系数 MPTEMP, MPTEMP,1,20,200,400,600,800,900,1000,1100,1200,1420,1460 MPDATA,DENS,2,1,7966,7893,7814,7724,7630,7583,7535,7486,7436,7320

4、,7320 MPDATA,KXX,2,1,32 ! mpdata,c,2,1,470,508,550,592,634,655,676,698,719,765,765 MPTEMP,1,0,200,400,800,1200 MPTEMP,6,1400,1454,1800,2000 MPDATA,ENTH,2,1,0,+8,+9,+9,+9 MPDATA,ENTH,2,6,+9,+9,+9,+10 BLOCK,0,Block_x,0,Block_y,0,Block_z,BLOCK,Block_x1,Block_x2,Block_y1,Block_y2,0,Block_z,vsel,allVOVLA

5、P,allBLOCK,0,Block_x,0,Block_y,0,-Block_z1,vglue,allnumcmp,all划分网格lesize,8,LMeshxlesize,5,LMeshylesize,13,LMeshzvsel,s,1type,1mat,1mshkey,1vmesh,allallselesize,mshkey,0mshape,1,3dvsel,s,2type,3mat,2vsel,s,3成形件轨迹单元质心排序width=Block_y2-Block_y1wpoffs,Block_x1,Block_y1,0local,11,1,Block_x1,Block_y1,0ALLS

6、EL,ALL*GET,EMAX,ELEM,COUNT,MAXVSEL,S, , , 1ESLV,S*GET,AEMAX,ELEM,COUNT1,MAX*DIM,ANE,ARRAY,AEMAX,*DIM,ANE1,ARRAY,AEMAX,*DIM,ANE2,ARRAY,AEMAX,*DIM,ANE3,ARRAY,AEMAX,*DIM,ANE4,ARRAY,AEMAX,*DIM,ANE5,ARRAY,AEMAX,*DIM,ANEX,ARRAY,AEMAX*DIM,ANEY,ARRAY,AEMAX*DIM,ANEZ,ARRAY,AEMAX*DIM,ANEORDER,ARRAY,AEMAX,1形心排序

7、*GET,ANSEL,ELEM,COUNT1 II=0 *DO,I,1,AEMAX *IF,ESEL（I）,EQ,1,THEN II=II+1 ANE（II）=I *ENDIF *ENDDO *DO,I,1,ANSEL *GET,ANEZ（I）,ELEM,ANE（I）,CENT,Z *GET,ANEY（I）,ELEM,ANE（I）,CENT,Y *GET,ANEX（I）,ELEM,ANE（I）,CENT,Xe=0*do,i,1,AEMAX *if,ANEY（I）,ge,4*width/5+Block_y1,then e=e+1 ANE5（e）=ANE（i）*endif *enddod=0*if

8、,ANEY（I）,ge,3*width/5+Block_y1,then *if,ANEY（I）,lt,4*width/5+Block_y1,then d=d+1 ANE4（d）=ANE（i） *endifc=0*if,ANEY（I）,ge,2*width/5+Block_y1,then *if,ANEY（I）,lt,3*width/5+Block_y1,then c=c+1 ANE3（c）=ANE（i）b=0*if,ANEY（I）,ge,width/5+Block_y1,then *if,ANEY（I）,lt,2*width/5+Block_y1,then b=b+1 ANE2（b）=ANE（

9、i）*endifa=0 *if,ANEY（I）,lt,width/5+Block_y1,then a=a+1 ANE1（a）=ANE（i） *对1道焊缝区的单元按质心X坐标进行从低到高排序*esel,none *do,i,1,a *if,ANE1（i）,NE,0,then esel,a,elem,ANE1（i）MINE=0*DIM,ANEX1,ARRAY,a*DIM,ANEY1,ARRAY,a*DIM,ANEZ1,ARRAY,a*DO,I1,1, aesel,u,elem,mine *GET,ANSEL1,ELEM,COUNT1 ANE1（II）=I *DO,I,1,ANSEL1 *GET,A

10、NEZ1（I）,ELEM,ANE1（I）,CENT,Z *GET,ANEY1（I）,ELEM,ANE1（I）,CENT,Y *GET,ANEX1（I）,ELEM,ANE1（I）,CENT,X MINZ=1E20 MINY=1E20 MINX=1E20 *IF,ANEZ1（I）,LT,MINZ,THEN MINZ=ANEZ1（I） MINY=ANEY1（I） MINX=ANEX1（I） MINE=ANE1（I） *ELSE *IF,ANEZ1（I）,EQ,MINZ,THEN *IF,ANEY1（I）,LT,MINY,THEN *IF,ANEY1（I）,EQ,MINY,THEN *IF,ANEX1

11、（I）,LT,MINX,THEN ANEORDER（I1）=MINE*ENDDO*对2道焊缝区的单元按质心X坐标进行从高到低排序*MAXE=1*do,j,1,b *if,ANE2（j）,NE,0,then esel,a,elem,ANE2（j）*DIM,ANEX2,ARRAY,b*DIM,ANEY2,ARRAY,b*DIM,ANEZ2,ARRAY,b*DO,J1,1, besel,u,elem,MAXE *GET,ANSEL2,ELEM,COUNT1 JJ=0 *DO,J,1,AEMAX JJ=JJ+1 ANE2（JJ）=J *DO,J,1,ANSEL2 *GET,ANEZ2（J）,ELEM,

12、ANE2（J）,CENT,Z *GET,ANEY2（J）,ELEM,ANE2（J）,CENT,Y *GET,ANEX2（J）,ELEM,ANE2（J）,CENT,X MAXZ=1E20 MAXY=1E20 MAXX=-1E20 *IF,ANEZ2（J）,LT,MAXZ,THEN MAXZ=ANEZ2（J） MAXY=ANEY2（J） MAXX=ANEX2（J） MAXE=ANE2（J） *IF,ANEZ2（J）,EQ,MAXZ,THEN *IF,ANEY2（J）,LT,MAXY,THEN *IF,ANEY2（J）,EQ,MAXY,THEN *IF,ANEX2（J）,GT,MAXX,THEN A

13、NEORDER（J1+a）=MAXE*对3道焊缝区的单元按质心X坐标进行从低到高排序*MINE=1*do,k,1,c *if,ANE3（k）,NE,0,then esel,a,elem,ANE3（k）*DIM,ANEX3,ARRAY,c*DIM,ANEY3,ARRAY,c*DIM,ANEZ3,ARRAY,c*DO,K1,1, c*GET,ANSEL3,ELEM,COUNT1 KK=0 *DO,K,1,AEMAX *IF,ESEL（K）,EQ,1,THEN KK=KK+1 ANE3（KK）=K *DO,K,1,ANSEL3 *GET,ANEZ3（K）,ELEM,ANE3（K）,CENT,Z *G

14、ET,ANEY3（K）,ELEM,ANE3（K）,CENT,Y *GET,ANEX3（K）,ELEM,ANE3（K）,CENT,X *IF,ANEZ3（K）,LT,MINZ,THEN MINZ=ANEZ3（K） MINY=ANEY3（K） MINX=ANEX3（K） MINE=ANE3（K） *IF,ANEZ3（K）,EQ,MINZ,THEN *IF,ANEY3（K）,LT,MINY,THEN *IF,ANEY3（K）,EQ,MINY,THEN *IF,ANEX3（K）,LT,MINX,THEN ANEORDER（K1+a+b）=MINE*对4道焊缝区的单元按质心X坐标进行从高到低排序*do,

15、j,1,d *if,ANE4（j）,NE,0,then esel,a,elem,ANE4（j）*DIM,ANEX4,ARRAY,d*DIM,ANEY4,ARRAY,d*DIM,ANEZ4,ARRAY,d*DO,J1,1, d*GET,ANSEL4,ELEM,COUNT1 *DO,J,1,ANSEL4 *GET,ANEZ4（J）,ELEM,ANE4（J）,CENT,Z *GET,ANEY4（J）,ELEM,ANE4（J）,CENT,Y *GET,ANEX4（J）,ELEM,ANE4（J）,CENT,X *IF,ANEZ4（J）,LT,MAXZ,THEN MAXZ=ANEZ4（J） MAXY=AN

16、EY4（J） MAXX=ANEX4（J） MAXE=ANE4（J） *IF,ANEZ4（J）,EQ,MAXZ,THEN *IF,ANEY4（J）,LT,MAXY,THEN *IF,ANEY4（J）,EQ,MAXY,THEN *IF,ANEX4（J）,GT,MAXX,THEN ANEORDER（J1+a+b+c）=MAXE*对5道焊缝区的单元按质心X坐标进行从低到高排序*do,i,1,e *if,ANE5（i）,NE,0,then esel,a,elem,ANE5（i）*DIM,ANEX5,ARRAY,e*DIM,ANEY5,ARRAY,e*DIM,ANEZ5,ARRAY,e*DO,I1,1,

17、e*GET,ANSEL5,ELEM,COUNT1 ANE5（II）=I *DO,I,1,ANSEL5 *GET,ANEZ5（I）,ELEM,ANE5（I）,CENT,Z *GET,ANEY5（I）,ELEM,ANE5（I）,CENT,Y *GET,ANEX5（I）,ELEM,ANE5（I）,CENT,X *IF,ANEZ5（I）,LT,MINZ,THEN MINZ=ANEZ5（I） MINY=ANEY5（I） MINX=ANEX5（I） MINE=ANE5（I） *IF,ANEZ5（I）,EQ,MINZ,THEN *IF,ANEY5（I）,LT,MINY,THEN *IF,ANEY5（I）,E

18、Q,MINY,THEN *IF,ANEX5（I）,LT,MINX,THEN ANEORDER（I1+a+b+c+d）=MINE*体热源加载求解*/SOLANTYPE,4 !瞬态分析TRNOPT,FULL !指定瞬态分析选项eqslv !指定求解器solcontrol,on !求解控制（用优化求解器）ESEL,ALLEPLOT/AUTO,1/REPLOTallsel,all outres,all,allTOFFST,273tunif,24kbc,1timint,onnsub1=2 nsub2=20 asel,s,loc,z,Block_z !加对流条件asel,a,loc,x,0asel,a,l

19、oc,x,Block_xasel,a,loc,y,0asel,a,loc,y,Block_ynslasf,all,conv,-1,24dt=3e-2 ! 定义每次成形所需要的时间tim=0生死单元控制*DO,I1,1,AEMAX ESEL,S,ANEORDER（I1） EKILL,ALL ESEL,S,LIVE EPLOT*enddo j=aemax !j小于AEmax*DO,I1,1,j nsubst,1 esel,s,ANEORDER（I1） EALIVE,ANEORDER（I1-1） EPLOT *ENDIF esel,all*体载荷加载* tim=tim+dt time,timbfe,ANEORDER（I1）,hgen,1, bfe,ANEORDER（I1）,hgen,1,*! T_end=300* SOLVE bfedele,ANEORDER（I1）,hgen/soluANTYPE,restbfedele,all,hgen time,4 !step101开始完全冷却NSUBST,15,30,5solve