24小时热门版块排行榜

返回列表

国王灬咸鱼

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2.5
帖子: 5
在线: 38分钟
虫号: 3305912
注册: 2014-07-03

[交流] 我仿真一个打印过程的温度变化，但是一直出错，请各位大神帮忙指点一下已有1人参与

/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

回复此楼