移动体热源ANSYS命令流Word格式.docx
- 文档编号:16379684
- 上传时间:2022-11-23
- 格式:DOCX
- 页数:22
- 大小:18.73KB
移动体热源ANSYS命令流Word格式.docx
《移动体热源ANSYS命令流Word格式.docx》由会员分享,可在线阅读,更多相关《移动体热源ANSYS命令流Word格式.docx(22页珍藏版)》请在冰豆网上搜索。
StandardUnits
*************************************************************************************************!
开始,进入前处理器,定义下列属性-
0参数定义
1定义单元类型
2实常数
3材料性能
4几何模型
5网格划分控制
6其他
/PREP7
0参数定义
模型参数
几何模型尺寸me-3=mm
*SET,Block_x,7.8e-3
*SET,Block_x1,0.9e-3
*SET,Block_x2,6.9e-3
*SET,Block_y,3.3e-3
*SET,Block_y1,0.9e-3
*SET,Block_y2,2.4e-3
*SET,Block_z,0.3e-3
*SET,Block_z1,0.9e-3
划分网格的大小
*SET,LMeshx,3e-4!
line1
*SET,LMeshy,3e-4!
line4
*SET,LMeshz,3e-4!
line9
边界条件参数
*SET,Tinit,30!
初始温度
ET,1,SOLID70
ET,2,surf152!
表面效应单元
r,2
KEYOPT,2,4,1!
表面效应单元设置
KEYOPT,2,8,4
ET,3,SOLID87
Nothing
MPTEMP
MPTEMP,1,25,400,1200,1600,2000,2500
MPDATA,DENS,1,1,3983,4688.4,5406.8,5948.8,6908,6331!
密度
MPDATA,KXX,1,1,1.6278,2.3812,3.6396,29.8,30.6,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,1,1,0,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9
MPDATA,ENTH,1,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10!
热焓
MPTEMP,,,,,,,,
MPTEMP,1,20,300,600,900,1200,1500,1800,2000,3000,4000
MPDATA,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,7320
MPDATA,KXX,2,1,13.31,16.33,19.47,22.38,25.07,26.33,27.53,28.67,29.76,31.95,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,7.9729e+8,1.6400e+9,3.4621e+9,5.3897e+9
MPDATA,ENTH,2,6,6.3415e+9,8.2595e+9,9.8798e+9,1.0816e+10
BLOCK,0,Block_x,0,Block_y,0,Block_z,
BLOCK,Block_x1,Block_x2,Block_y1,Block_y2,0,Block_z,
vsel,all
VOVLAP,all
BLOCK,0,Block_x,0,Block_y,0,-Block_z1,
vglue,all
numcmp,all
划分网格
lesize,8,LMeshx
lesize,5,LMeshy
lesize,13,LMeshz
vsel,s,,,1
type,1
mat,1
mshkey,1
vmesh,all
allsel
esize,0.3e-3
mshkey,0
mshape,1,3d
vsel,s,,,2
type,3
mat,2
vsel,s,,,3
成形件轨迹单元质心排序
width=Block_y2-Block_y1
wpoffs,Block_x1,Block_y1,0
local,11,1,Block_x1,Block_y1,0
ALLSEL,ALL
*GET,EMAX,ELEM,,COUNT,MAX
VSEL,S,,,1
ESLV,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
形心排序
*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,X
e=0
*do,i,1,AEMAX
*if,ANEY(I),ge,4*width/5+Block_y1,then
e=e+1
ANE5(e)=ANE(i)
*endif
*enddo
d=0
*if,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)
*endif
c=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(i)
*endif
a=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,a
esel,u,elem,,mine
*GET,ANSEL1,ELEM,,COUNT1
ANE1(II)=I
*DO,I,1,ANSEL1
*GET,ANEZ1(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(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,b
esel,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,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
ANEORDER(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
*GET,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,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=ANEY4(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,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),EQ,MINY,THEN
*IF,ANEX5(I),LT,MINX,THEN
ANEORDER(I1+a+b+c+d)=MINE
*********************体热源加载求解**************************************
/SOL
ANTYPE,4!
瞬态分析
TRNOPT,FULL!
指定瞬态分析选项
eqslv!
指定求解器
solcontrol,on!
求解控制(用优化求解器)
ESEL,ALL
EPLOT
/AUTO,1
/REPLOT
allsel,all
outres,all,all
TOFFST,273
tunif,24
kbc,1
timint,on
nsub1=2
nsub2=20
asel,s,loc,z,Block_z!
加对流条件
asel,a,loc,x,0
asel,a,loc,x,Block_x
asel,a,loc,y,0
asel,a,loc,y,Block_y
nsla
sf,all,conv,-1,24
dt=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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 移动 热源 ANSYS 命令