APDL温度场分析时出现错误，好几天都未解决，求大神帮忙！

hgdgaojinling

APDL进行温度场分析时出现Parameter NEY22 needs Subscripts。正在做课题，很着急！希望大神能给我一些解决问题的建议，不胜感激，多谢！！！

Vicly2012

建议一下，最好附上APDL程序，这样方便大神帮你发现解决问题。

小七工作室

子结构有问题吧

hgdgaojinling

FINISH
/CLEAR,START           !开始一个新的分析
/COLOR,PBAK,OFF        !关闭窗口的背景颜色
/TRIAD,OFF             !关闭整体坐标系的符号
/FILN,Temperature Field Analysis                        !指定一个文件名
/TITLE,Temperature Field Analysis of FDM RP Machine     !指定一个分析标题
/PREP7                 !进入前处理器
/UNIT,SI               !采用国际单位制
ET,1,SOLID5            !采用三维8节点六面体热分析单元SOLID5
MPTEMP,,,,,,,,
MPTEMP,1,25           !在材料特性中定义初始温度
MPDATA,EX,1,,2.1E9     !定义材料的弹性模量
MPDATA,PRXY,1,,0.3     !定义材料泊松比
MPDATA,DENS,1,,1.06    !定义材料密度
MPDATA,KXX,1,,0.16     !定义材料的热导率
MPDATA,ALPX,1,,8.8E-5  !定义材料的热膨胀系数
MPDATA,HF,1,,72        !定义材料的对流系数
MPDATA,C,1,,2          !定义材料的比热容
MPDATA,MURX,1,1
!建立几何模型
L=40                   !定义模型的长
W=8                    !定义模型的宽
H=12                   !定义模型的高
BLOCK,0,40,0,12,0,8    !建立长方体模型
!划分网格
TYPE,1                 !指定单元属性，即SOLID5单元
MAT,1                  !指定材料属性
ESIZE,1                !指定划分单元的长度
VMESH,ALL              !采用扫略的方式划分网格
/SOLU
ANTYPE,4               !定义分析类型为瞬态分析
NSEL,ALL
NSEL,S,LOC,Y,0
D,ALL,TEMP,110
!生成生死单元并添加载荷
ESEL,ALL               !选择所有单元              
*GET,NSE,ELEM,0,COUNT
*DIM,NE,,NSE           !定义NE数组，记录单元个数
L0=1                   !长度方向单元长度为1
W0=1                   !宽度方向单元长度为1
H0=1                   !高度方向单元长度为1
ESEL,ALL
*DIM,NEX21,ARRAY,40         !定义NEX21数组,记录长度方向单元个数
*DIM,NEY22,ARRAY,12        !定义NEY22数组,记录宽度方向单元个数
*DIM,NEZ23,ARRAY,8         !定义NEZ23数组,记录高度方向单元个数
MAT_TEM=220            !材料从喷嘴挤出时的温度
DT1=1E-3               !用于建立初始条件的一个很小的时间段
DT=1                   !材料堆积一个单元所需的时间
T=0                    !起始时间
ESEL,S,ELEM,,1,NSE,1,0
EKILL,ALL
ALLSEL,ALL
OUTRES,ALL,ALL         !设定输出选项
IC,ALL,TEMP,25         !施加初始温度值
KBC,1                  !设定载荷增进方式为stepped
TIMINT,0,STRUCT
TIMINT,1,THERM
!TIMINT,0,MAG
TINTP,0.005,,,1,0.5,0.2
NSUB1=2                !设定子步数1
NSUB2=40               !设定子步数2
*DO,NEY22,1,H/H0,1
  *DO,NEZ23,1,W/W0,1
    *DO,NEX21,1,L/L0,1
     ESEL,S,ELEM,,(NEY22-1)*L/L0*W/W0+(NEZ23-1)*L/L0+NEX21,,,0
     EALIVE,ALL
      *DO,J,1,8
       D,NELEM((NEY22-1)*L/L0*W/W0+(NEZ23-1)*L/L0+NEX21,J),TEMP,MAT_TEM
      *ENDDO
      ESEL,S,ELEM,,(NEY22-1)*L/L0*W/W0+(NEZ23-1)*L/L0+NEX21,,,0
      NSEL,S,1
      SF,ALL,CONV,72,25
      ALLSEL,ALL
      SOLVE
      !下面的求解用于保证初始温度的升温速度为零
      T=T+DT1
      TIME,T
      SOLVE
      ESEL,S,ELEM,,(NEY22-1)*L/L0*W/W0+(NEZ23-1)*L/L0+NEX21,,,0
      *DO,J,1,8
      DDELE,NELEM((NEY22-1)*L/L0*W/W0+(NEZ23-1)*L/L0+NEX21,J),TEMP
      *ENDDO
    *ENDDO
  *ENDDO
*ENDDO   
T=T+DT-2*DT1
TIME,T
NSUBST,NSUB1
SOLVE
T=T+5000
TIME,T
NSUBST,NSUB2
SOLVE
FINISH
/POST
SET,100                !读取第100步的计算结果
PLNSOL,TEMP,,0,1.0     !查看温度场
PLNSOL,S,EQV,2,1.0     !查看残余应力场
SET,400                !读取第400步的计算结果
PLNSOL,TEMP,,0,1.0     !查看温度场
PLNSOL,S,EQV,2,1.0     !查看残余应力场
SET,LAST               !读取最后一步的计算结果
PLNSOL,TEMP,,0,1.0     !查看温度场
PLNSOL,S,EQV,2,1.0     !查看残余应力场

hgdgaojinling

希望多多交流。。。先感谢一下你们的关注。。。