我所做的ansys模拟T型管的焊接
浏览:101127
/prep7
mptemp,1,0,100,200,300,400,500 !定义材料属性温度范围
mptemp,,1200,1250,1350,1400,1450,1500
mptemp,,1505,1705,1905,2105,2305,3000
mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定义材料密度(kg/m3)
mpdata,dens,1,,7800,7800,7800,7800,7800,7800
mpdata,dens,1,,7800,7800,7800,7800,7800,7800
mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K))
mpdata,kxx,1,,30,30,30,30,30,30
mpdata,kxx,1,,30,30,30,30,30,30
mpdata,c,1,1,450,469,481,508,536,569 !考虑相变潜热的比热容(J/kg*K)
mpdata,c,1,,700,2172,5116,6589,8061,9533
mpdata,c,1,,9533,7757,5982,4206,2431,700
!******************考虑辐射和对流的总换热系数(W/m^2*℃),加载时按对流加载
*dim,conve,table,8,1,1,temp
*SET,conve(0,1,1),1
*SET,conve(1,0,1),0
*SET,conve(1,1,1),12.12
*SET,conve(2,0,1),200
*SET,conve(2,1,1),24.23
*SET,conve(3,0,1),400
*SET,conve(3,1,1),40.38
*SET,conve(4,0,1),600
*SET,conve(4,1,1),72.69
*SET,conve(5,0,1),800
*SET,conve(5,1,1),105.00
*SET,conve(6,0,1),1000
*SET,conve(6,1,1),145.38
*SET,conve(7,0,1),1200
*SET,conve(7,1,1),210.00
*SET,conve(8,0,1),1400
*SET,conve(8,1,1),278.65
!**************************
!定义单元类型并划分网格
!**************************
et,1,solid70 !温度场计算单元类型
et,2,surf152 !表面效应单元类型
keyopt,2,3,1 !设置表面效应单元选项
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,7,0
keyopt,2,8,3
keyopt,2,9,0
/prep7
/units,si
cyl4,0,0,1.413/2,90,1.439/2,270,2
local,11,0,0,-2,2, ,90
wpcsys,-1,11
cyl4,0,0,2.955/2,90,1.5,270,5
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
VDELE, 2, , ,1
cyl4,0,0,2.955/2,90,1.5,270,-5
vovlap,all
FLST,2,3,6,ORDE,3
FITEM,2,3
FITEM,2,-4
FITEM,2,7
VDELE,P51X, , ,1
VDELE, 6, , ,1
/DIST,1,0.924021086472,1
/REP,FAST
FLST,3,1,5,ORDE,1
FITEM,3,26
AGEN,2,P51X, , , , ,0.02, ,0
/DIST,1,0.924021086472,1
/REP,FAST
APLOT
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
ADELE, 1, , ,1
/DIST,1,0.924021086472,1
/REP,FAST
FLST,3,1,5,ORDE,1
FITEM,3,26
AGEN,2,P51X, , , ,0.02, , ,0
VSBA, 5, 1
FLST,5,6,4,ORDE,6
FITEM,5,2
FITEM,5,4
FITEM,5,-5
FITEM,5,7
FITEM,5,28
FITEM,5,40
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,40, , , , ,1
!*
FLST,5,4,4,ORDE,2
FITEM,5,9
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,300, , , , ,1
!*
FLST,5,5,4,ORDE,5
FITEM,5,6
FITEM,5,8
FITEM,5,14
FITEM,5,-15
FITEM,5,42
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,3, , , , ,1
!*
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,3,2,6,ORDE,2
FITEM,3,1
FITEM,3,-2
VSYMM,X,P51X, , , ,0,0
ALLSEL,ALL
ASEL,S,EXT
aatt,,,2
mshkey,1
type,2
amesh,all
esel,s,type,,2 !表面效应单元上施加对流条件
nsle,s,all
sf,all,conv,%conve%,20
FINISH
/SOL
Csys,1
Antype,4
Trnopt,full
Tunif,27
nropt,full,,on !定义完全牛顿--拉普森方法,激活自适应下降(默认)
pred,on !打开预测校正
timint,on !时间积分设置
tref,20 !参考温度20摄氏度
!----------------------------------------------
!设置载荷步参数
!----------------------------------------------
*SET,dt0,1e-6 !建立初始条件的小时间段
!****************************
vsel,s,,,2,4,2
eslv
ekill,all
/REPLOT
!稳态分析确定初始温度场
!***************************************************************
time,dt0 !初始条件,loadstep=1,time=1e-6s
timint,off !时间积分选项,关闭(稳态分析)
kbc,1 !阶越加载
nsubst,1 !子步为1步
esel,s,type,,1 !选择类型1的单元节点
nsle,s
ic,all,temp,20 !初始状态,温度20摄氏度
alls
outres,all,all
solve
!***************************
*SET,Lsize,1
*SET,V,5!速度
*SET,y0,0
*SET,tinc,lsize/v
*SET,tim,1e-6
ALLSEL
*GET,EMAX,ELEM,,NUM,MAX !获取结构上最大单元数EMAX
vsel,s,,,2,4,2
eslv
*GET,NSE,ELEM,,COUNT !获取焊缝上最大单元数NSE
*DIM,NE,,NSE ! 定义NE数组,记录所有单元编号
*DIM,NEX,,NSE ! 定义NEX数组,记录所有焊缝单元的系统坐标系的X坐标
*DIM,NEY,,NSE ! 定义NEY数组,记录所有焊缝单元的系统坐标系的Y坐标
*DIM,NEZ,,NSE ! 定义NEZ数组,记录所有焊缝单元的系统坐标系的Z坐标
*DIM,NEORDER,,NSE ! 定义单元编号数组,按几何位置排列,按层,道,角度的顺序
*SET,MINE,0 ! 定义初值,记录在循环中每次选中的单元
*SET,II,0
*DO,I,1,EMAX
*IF,ESEL(I),EQ,1,THEN
*SET,II,II+1
*SET,NE(II),I ! 焊缝上单元顺序排列
*ENDIF
*ENDDO
*DO,I,1,NSE
*GET,NEX(I),ELEM,NE(I),CENT,X
*GET,NEY(I),ELEM,NE(I),CENT,Y
*GET,NEZ(I),ELEM,NE(I),CENT,Z
*ENDDO
*do,J,1,180,1
*SET,tim,j*0.5
time,tim
antype,4 !瞬态分析
timint,on !时间积分
autots,on !自动时间步长
kbc,0 !阶越加载
*SET,YMAX,J
*SET,YMIN,(J-1)
*DO,II,1,NSE
*IF,NEY(II),Le,YMAX,THEN
*IF,NEY(II),Ge,YMIN,THEN
ealive,ne(ii)
BFE,NE(II),HGEN,,2.5E9
*ENDIF
*ENDIF
*ENDDO
ALLSEL
SOLVE
bfedele,all,hgen
*enddo
*do,J,180,1,-1
*SET,tim1,36+(181-J)*tinc
time,tim1
antype,4 !瞬态分析
timint,on !时间积分
autots,on !自动时间步长
kbc,0 !阶越加载
*SET,YMAX,-(J-1)
*SET,YMIN,-j
*DO,II,1,NSE
*IF,NEY(II),Le,YMAX,THEN
*IF,NEY(II),Ge,YMIN,THEN
ealive,ne(ii)
BFE,NE(II),HGEN,,2.5E9
*ENDIF
*ENDIF
*ENDDO
ALLSEL
SOLVE
bfedele,all,all,hgen
*enddo
以上为命令流,为何不能实现热源的移动?
mptemp,1,0,100,200,300,400,500 !定义材料属性温度范围
mptemp,,1200,1250,1350,1400,1450,1500
mptemp,,1505,1705,1905,2105,2305,3000
mpdata,dens,1,1,7820,7800,7800,7800,7800,7800 !定义材料密度(kg/m3)
mpdata,dens,1,,7800,7800,7800,7800,7800,7800
mpdata,dens,1,,7800,7800,7800,7800,7800,7800
mpdata,kxx,1,1,52,50.7,48.6,46.1,42.3,38.9 !热传导系数(W/(m*K))
mpdata,kxx,1,,30,30,30,30,30,30
mpdata,kxx,1,,30,30,30,30,30,30
mpdata,c,1,1,450,469,481,508,536,569 !考虑相变潜热的比热容(J/kg*K)
mpdata,c,1,,700,2172,5116,6589,8061,9533
mpdata,c,1,,9533,7757,5982,4206,2431,700
!******************考虑辐射和对流的总换热系数(W/m^2*℃),加载时按对流加载
*dim,conve,table,8,1,1,temp
*SET,conve(0,1,1),1
*SET,conve(1,0,1),0
*SET,conve(1,1,1),12.12
*SET,conve(2,0,1),200
*SET,conve(2,1,1),24.23
*SET,conve(3,0,1),400
*SET,conve(3,1,1),40.38
*SET,conve(4,0,1),600
*SET,conve(4,1,1),72.69
*SET,conve(5,0,1),800
*SET,conve(5,1,1),105.00
*SET,conve(6,0,1),1000
*SET,conve(6,1,1),145.38
*SET,conve(7,0,1),1200
*SET,conve(7,1,1),210.00
*SET,conve(8,0,1),1400
*SET,conve(8,1,1),278.65
!**************************
!定义单元类型并划分网格
!**************************
et,1,solid70 !温度场计算单元类型
et,2,surf152 !表面效应单元类型
keyopt,2,3,1 !设置表面效应单元选项
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,7,0
keyopt,2,8,3
keyopt,2,9,0
/prep7
/units,si
cyl4,0,0,1.413/2,90,1.439/2,270,2
local,11,0,0,-2,2, ,90
wpcsys,-1,11
cyl4,0,0,2.955/2,90,1.5,270,5
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
VDELE, 2, , ,1
cyl4,0,0,2.955/2,90,1.5,270,-5
vovlap,all
FLST,2,3,6,ORDE,3
FITEM,2,3
FITEM,2,-4
FITEM,2,7
VDELE,P51X, , ,1
VDELE, 6, , ,1
/DIST,1,0.924021086472,1
/REP,FAST
FLST,3,1,5,ORDE,1
FITEM,3,26
AGEN,2,P51X, , , , ,0.02, ,0
/DIST,1,0.924021086472,1
/REP,FAST
APLOT
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
ADELE, 1, , ,1
/DIST,1,0.924021086472,1
/REP,FAST
FLST,3,1,5,ORDE,1
FITEM,3,26
AGEN,2,P51X, , , ,0.02, , ,0
VSBA, 5, 1
FLST,5,6,4,ORDE,6
FITEM,5,2
FITEM,5,4
FITEM,5,-5
FITEM,5,7
FITEM,5,28
FITEM,5,40
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,40, , , , ,1
!*
FLST,5,4,4,ORDE,2
FITEM,5,9
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,300, , , , ,1
!*
FLST,5,5,4,ORDE,5
FITEM,5,6
FITEM,5,8
FITEM,5,14
FITEM,5,-15
FITEM,5,42
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,3, , , , ,1
!*
FLST,5,2,6,ORDE,2
FITEM,5,1
FITEM,5,-2
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VSWEEP,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,3,2,6,ORDE,2
FITEM,3,1
FITEM,3,-2
VSYMM,X,P51X, , , ,0,0
ALLSEL,ALL
ASEL,S,EXT
aatt,,,2
mshkey,1
type,2
amesh,all
esel,s,type,,2 !表面效应单元上施加对流条件
nsle,s,all
sf,all,conv,%conve%,20
FINISH
/SOL
Csys,1
Antype,4
Trnopt,full
Tunif,27
nropt,full,,on !定义完全牛顿--拉普森方法,激活自适应下降(默认)
pred,on !打开预测校正
timint,on !时间积分设置
tref,20 !参考温度20摄氏度
!----------------------------------------------
!设置载荷步参数
!----------------------------------------------
*SET,dt0,1e-6 !建立初始条件的小时间段
!****************************
vsel,s,,,2,4,2
eslv
ekill,all
/REPLOT
!稳态分析确定初始温度场
!***************************************************************
time,dt0 !初始条件,loadstep=1,time=1e-6s
timint,off !时间积分选项,关闭(稳态分析)
kbc,1 !阶越加载
nsubst,1 !子步为1步
esel,s,type,,1 !选择类型1的单元节点
nsle,s
ic,all,temp,20 !初始状态,温度20摄氏度
alls
outres,all,all
solve
!***************************
*SET,Lsize,1
*SET,V,5!速度
*SET,y0,0
*SET,tinc,lsize/v
*SET,tim,1e-6
ALLSEL
*GET,EMAX,ELEM,,NUM,MAX !获取结构上最大单元数EMAX
vsel,s,,,2,4,2
eslv
*GET,NSE,ELEM,,COUNT !获取焊缝上最大单元数NSE
*DIM,NE,,NSE ! 定义NE数组,记录所有单元编号
*DIM,NEX,,NSE ! 定义NEX数组,记录所有焊缝单元的系统坐标系的X坐标
*DIM,NEY,,NSE ! 定义NEY数组,记录所有焊缝单元的系统坐标系的Y坐标
*DIM,NEZ,,NSE ! 定义NEZ数组,记录所有焊缝单元的系统坐标系的Z坐标
*DIM,NEORDER,,NSE ! 定义单元编号数组,按几何位置排列,按层,道,角度的顺序
*SET,MINE,0 ! 定义初值,记录在循环中每次选中的单元
*SET,II,0
*DO,I,1,EMAX
*IF,ESEL(I),EQ,1,THEN
*SET,II,II+1
*SET,NE(II),I ! 焊缝上单元顺序排列
*ENDIF
*ENDDO
*DO,I,1,NSE
*GET,NEX(I),ELEM,NE(I),CENT,X
*GET,NEY(I),ELEM,NE(I),CENT,Y
*GET,NEZ(I),ELEM,NE(I),CENT,Z
*ENDDO
*do,J,1,180,1
*SET,tim,j*0.5
time,tim
antype,4 !瞬态分析
timint,on !时间积分
autots,on !自动时间步长
kbc,0 !阶越加载
*SET,YMAX,J
*SET,YMIN,(J-1)
*DO,II,1,NSE
*IF,NEY(II),Le,YMAX,THEN
*IF,NEY(II),Ge,YMIN,THEN
ealive,ne(ii)
BFE,NE(II),HGEN,,2.5E9
*ENDIF
*ENDIF
*ENDDO
ALLSEL
SOLVE
bfedele,all,hgen
*enddo
*do,J,180,1,-1
*SET,tim1,36+(181-J)*tinc
time,tim1
antype,4 !瞬态分析
timint,on !时间积分
autots,on !自动时间步长
kbc,0 !阶越加载
*SET,YMAX,-(J-1)
*SET,YMIN,-j
*DO,II,1,NSE
*IF,NEY(II),Le,YMAX,THEN
*IF,NEY(II),Ge,YMIN,THEN
ealive,ne(ii)
BFE,NE(II),HGEN,,2.5E9
*ENDIF
*ENDIF
*ENDDO
ALLSEL
SOLVE
bfedele,all,all,hgen
*enddo
以上为命令流,为何不能实现热源的移动?