仿真3D打印过程温度变化,但一直出问题,请各位大神帮忙看一下
我这是仿真的3D打印过程温度变化的命令流。我是现将所有单元都杀死,然后一个一个的激活,激活一个单元后,给这个单元加载进行分析,分析完之后将该单元 卸载,然后给下一单元加载再进行分析,一直重复此过程,直到所有单元分析完成。但在生死单元循环中,即*DO循环中。那个RESTART命令提示必须加一 个FINISH命令才能开始RESTART,但是加了之后,仍然不对。请各位帮忙看一下我的命令流哪里有问题,谢谢啦
/BATCH
FINISH
/CLEAR,START
/CWD,G:\Exercise
/FILNAME,FDMTEMPERTURE,0
/TITLE,ANSYS OF FDM TEMPERTURE
*SET,L1,1E-2
*SET,W1,1.2E-3
*SET,H1,0.3E-3 !设置模型尺寸为10mm*1.2mm*0.3mm
!设置网格尺寸
*SET,L2,L1/25
*SET,W2,0.4E-3
*SET,H2,H1/3 !设置单元尺寸为0.4mm*0.4mm*0.1mm(喷嘴直径为0.4,分层厚度0.1mm)
ARG1=L1/L2*W1/W2*H1/H2 !计算单元总数量
/CONFIG,NRES,10000
!*************************************
!进入前处理
!*************************************
/PREP7
ET,1,SOLID70
MP,DENS,1,1150 !设置密度
MPTEMP,1,50,100,150,200 !指定四个温度
MPDATA,C,1,1,1470,1490,1710,2020 !比热容
MPDATA,KXX,1,1,3.0E-2,2.8E-2,2.9E-2,3.3E-2 !导热系数
MPDATA,EX,1,1,3.5E5,2.48E9,1.68E9,0.5E9 !弹性模量
MPDATA,NUXY,1,1,0.38,0.39,0.40,0.41 !泊松比
MPDATA,ALPX,1,1,8.51E-5,8.42E-5,8.40E-5,8.38E-5 !热膨胀系数
BLOCK,0,L1,0,W1,0,H1 !建立成型件的20mm*1mm*0.3mm的模型
!网格划分
/PNUM,AREA,1 !显示面编号
VATT,1,,0 !体单元类型solid70
/PNUM,LINE,1
MSHAPE,0,3D !单元形状六面体
MSHKEY,1 !映射网格划分
LESIZE,1,W2
LESIZE,4,L2
LESIZE,9,H2 !单元尺寸大小
VMESH,1
EPLOT
/VIEW,,-1,-2,-3
/IMAGE,SAVE,VMESH,BMP
FINISH
!**************************************
!进入载荷加载和求解
!*************************************
/SOLU
ANTYPE,TRANSIENT,NEW
AR21=0 !每一个载荷步加载时间是0.02s,单元激活时间间隔也是0.02s,扫描速度20mm/s
AR22=0.02
TIME,AR21
AUTOTS,ON !自动时间步选项(非线性分析打开加速收敛)
NROPT,FULL,,ON ! 牛顿-拉普森选项
TRNOPT,FULL !瞬态分析选项,指定最大分析模态数
LUMPM,OFF !使用一致质量矩阵
KBC,1
DELTIM,AR21/4,AR21/4,AR21/2,0 !载荷子步设定
TIMINT,ON
MAX_TEMP=220 !设定温度值220摄氏度
NSEL,S,LOC,X,0,L1,L2
SF,ALL,CONV,72,25
NSEL,S,LOC,Y,0,W1,W2
SF,ALL,CONV,72,25
NSEL,S,LOC,Z,0,H1,H2
SF,ALL,CONV,72,25 !初始化所有节点,加载对流
/REPLOT,FAST
OUTRES,ALL,ALL
OUTPR,NSOL,ALL
!生死单元循环命令
ESEL,S,ELEM,,1,ARG1,1,0
EKILL,ALL
ALLSEL,ALL
*DO,AR25,1,H1/H2,1
*DO,AR23,1,L1/L2,1
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0 !依次选择下一个单元激活
EALIVE,ALL
*DO,J,1,8,1
D,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP,MAX_TEMP !给每个单元的八个节点施加恒定温度
*ENDDO
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0
NSLE,S,ALL !选择单元中的所有节点
SF,ALL,CONV,72,25
ALLSEL,ALL
SOLVE
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP !删除每个单元的八个节点施加恒定温度
*ENDDO
ALLSEL,ALL
AR22=AR22+AR21
!*IF,AR22,GT,AR21*ARG1,*EXIT
FINISH
ANTYPE,TRANSIENT,RESTART
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+AR23,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+AR23,J),TEMP !删除每个单元的八个节点施加恒定温度
*ENDDO
TIME,AR22
*ENDDO
*DO,AR26,L1/L2,1,-1
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+L1/L2+AR26,,,0
EALIVE,ALL
*ENDDO
*DO,AR27,L1/L2,1,-1
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
EALIVE,ALL
*DO,J,1,8,1
D,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP,MAX_TEMP !给每个单元的八个节点施加恒定温度
*ENDDO
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
NSLE,S,ALL !选择单元中的所有节点
SF,ALL,CONV,72,25
ALLSEL,ALL
SOLVE
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP !删除每个单元的八个节点施加恒定温度
*ENDDO
ALLSEL,ALL
AR22=AR22+AR21
!*IF,AR22,GT,AR21*ARG1,*EXIT
FINISH
ANTYPE,TRANSIENT,RESTART
ESEL,S,ELEM,,(AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,,,0
*DO,J,1,8,1
DDELE,NELEM((AR25-1)*W1/W2*L1/L2+2*L1/L2+AR27,J),TEMP !删除每个单元的八个节点施加恒定温度
*ENDDO
TIME,AR22
*ENDDO
*ENDDO
SAVE
/POST1
PLNSOL,TEMP,,1
ANDATA,0.02,1,2,0,0,1,0,0 !生成某一范围内的结果数据的顺序梯度线动画
FINISH
!请各位大神帮帮忙呀
怎么解决了嘛 ·? 这边