急死了，看看这个帖子，请高手帮忙！非常感谢

scrunize

出现一个问题，由大体积混凝土温度场转换为应力分析时，怎么混凝土表面出现的值也是负的，命令流如下：（网格划分使用智能网格划分的）
!热分析
/filename,themal
/title,themal_analysis  !标题
/PREP7                  !前处理

!定义单元类型
ET,1,SOLID70

!设置材料属性
   !材料1为混凝土  
MP,KXX,1,280  
MP,C,1,0.976            
MP,DENS,1,2500      
MAT,1                  
   !材料2为基础土  
MP,KXX,2,300        
MP,C,2,1.45            
MP,DENS,2,1670         
MAT,2   

!输入参数
   
   !输入外界气温表
*DIM,qiwen,TABLE,57,1,1,TIME, , ,  !外界气温表名为qiwen
QIWEN(1,0,1) = 0
QIWEN(1,1,1) = 29
QIWEN(2,0,1) = 0.5
QIWEN(2,1,1) =30
QIWEN(3,0,1) = 1
QIWEN(3,1,1) = 24  
QIWEN(4,0,1) =1.5
QIWEN(4,1,1) = 30
QIWEN(5,0,1) =2
QIWEN(5,1,1) = 29
QIWEN(6,0,1) = 2.5
QIWEN(6,1,1) = 28.6
QIWEN(7,0,1) = 3
QIWEN(7,1,1) = 30  
QIWEN(8,0,1) =3.5
QIWEN(8,1,1) = 29   
QIWEN(9,0,1) = 4
QIWEN(9,1,1) = 28
QIWEN(10,0,1) = 4.5
QIWEN(10,1,1) = 30  
QIWEN(11,0,1) =5
QIWEN(11,1,1) = 30
QIWEN(12,0,1) =5.5
QIWEN(12,1,1) = 31  
QIWEN(13,0,1) = 6  
QIWEN(13,1,1) = 29
QIWEN(14,0,1) = 6.5
QIWEN(14,1,1) = 30
QIWEN(15,0,1) = 7  
QIWEN(15,1,1) = 30
QIWEN(16,0,1) = 7.5
QIWEN(16,1,1) = 31
QIWEN(17,0,1) =8
QIWEN(17,1,1) = 29
QIWEN(18,0,1) = 8.5
QIWEN(18,1,1) = 26
QIWEN(19,0,1) =9   
QIWEN(19,1,1) = 24
QIWEN(20,0,1) = 9.5
QIWEN(20,1,1) = 26  
QIWEN(21,0,1) = 10  
QIWEN(21,1,1) = 29  
QIWEN(22,0,1) = 10.5  
QIWEN(22,1,1) = 29
QIWEN(23,0,1) = 11  
QIWEN(23,1,1) = 29
QIWEN(24,0,1) = 11.5  
QIWEN(24,1,1) = 29
QIWEN(25,0,1) = 12  
QIWEN(25,1,1) = 29
QIWEN(26,0,1) = 12.5  
QIWEN(26,1,1) = 29
QIWEN(27,0,1) = 13  
QIWEN(27,1,1) = 29
QIWEN(28,0,1) = 13.5  
QIWEN(28,1,1) = 29
QIWEN(29,0,1) = 14
QIWEN(29,1,1) = 29
QIWEN(30,0,1) = 14.5
QIWEN(30,1,1) = 29
QIWEN(31,0,1) = 15
QIWEN(31,1,1) = 29
QIWEN(32,0,1) = 15.5
QIWEN(32,1,1) = 29
QIWEN(33,0,1) = 16
QIWEN(33,1,1) = 29
QIWEN(34,0,1) = 16.5
QIWEN(34,1,1) = 29
QIWEN(35,0,1) = 17
QIWEN(35,1,1) = 29
QIWEN(36,0,1) = 17.5
QIWEN(36,1,1) = 29
QIWEN(37,0,1) = 18
QIWEN(37,1,1) = 29
QIWEN(38,0,1) = 18.5
QIWEN(38,1,1) = 29
QIWEN(39,0,1) = 19
QIWEN(39,1,1) = 29
QIWEN(40,0,1) = 19.5
QIWEN(40,1,1) = 29
QIWEN(41,0,1) = 20
QIWEN(41,1,1) = 29
QIWEN(42,0,1) = 20.5
QIWEN(42,1,1) = 29
QIWEN(43,0,1) = 21
QIWEN(43,1,1) = 29
QIWEN(44,0,1) = 21.5
QIWEN(44,1,1) = 29
QIWEN(45,0,1) = 22
QIWEN(45,1,1) = 29
QIWEN(46,0,1) = 22.5
QIWEN(46,1,1) = 29
QIWEN(47,0,1) = 23
QIWEN(47,1,1) = 29
QIWEN(48,0,1) = 23.5
QIWEN(48,1,1) = 29
QIWEN(49,0,1) = 24
QIWEN(49,1,1) = 29
QIWEN(50,0,1) = 24.5
QIWEN(50,1,1) = 29
QIWEN(51,0,1) = 25
QIWEN(51,1,1) = 29
QIWEN(52,0,1) = 25.5
QIWEN(52,1,1) = 29
QIWEN(53,0,1) = 26
QIWEN(53,1,1) = 29
QIWEN(54,0,1) = 26.5
QIWEN(54,1,1) = 29
QIWEN(55,0,1) = 27
QIWEN(55,1,1) = 29
QIWEN(56,0,1) = 27.5
QIWEN(56,1,1) = 29
QIWEN(57,0,1) = 28
QIWEN(57,1,1) = 29

  !输入热生成率函数
*SET,_FNCNAME,'heats'   
*DIM,%_FNCNAME%,TABLE,6,16,1
   ! Begin of equation: 2450*0.82*280*0.13*exp(-0.82*{TIME})  
%_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0  
%_FNCNAME%(3,0,1)= 0.0  
%_FNCNAME%(4,0,1)= 0.0  
%_FNCNAME%(5,0,1)= 0.0  
%_FNCNAME%(6,0,1)= 0.0  
%_FNCNAME%(0,1,1)= 1.0, -1, 0, 2500, 0, 0, 0
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 0.8, 0, 0, -1
%_FNCNAME%(0,3,1)=   1, -3, 0, 1, -1, 3, -2
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 280, 0, 0, -3   
%_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -3, 3, -1
%_FNCNAME%(0,6,1)= 0.0, -1, 0, 0.205, 0, 0, -2   
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,8,1)= 0.0, -1, 0, 0, 0, 0, 0   
%_FNCNAME%(0,9,1)= 0.0, -2, 0, 1, 0, 0, -1  
%_FNCNAME%(0,10,1)= 0.0, -4, 0, 1, -1, 2, -2
%_FNCNAME%(0,11,1)= 0.0, -1, 0, 0.80, 0, 0, -4
%_FNCNAME%(0,12,1)= 0.0, -2, 0, 1, -4, 3, -1
%_FNCNAME%(0,13,1)= 0.0, -1, 0, 1, -2, 3, 1
%_FNCNAME%(0,14,1)= 0.0, -1, 7, 1, -1, 0, 0
%_FNCNAME%(0,15,1)= 0.0, -2, 0, 1, -3, 3, -1
%_FNCNAME%(0,16,1)= 0.0, 99, 0, 1, -2, 0, 0
   ! End of equation: 2450*0.7462*227.2246*0.13*exp(-0.7426*{TIME})

!创建几何体
BLC4,-14.0,-16,14.0,16,2.6          !创建混凝土块
BLOCK,-14.0,-5.5,-16,-11,-1.8,0         
vadd,1,2
BLOCK,-14.0,10,-11,10,-5,0           !创建基础土体
BLOCK,-14.0,-5.5,-16,-11,-5,-1.8
BLOCK,-5.5,10,-16,-11,-5,0
vadd,2,4,1
vglue,all
NUMMRG,ALL                  !合并所有重复元素
NUMCMP, ALL                 !压缩所有元素的编号
CSYS,0                      !返回全局坐标系

      !进行网格划分
finish

/SOLU                       !进入求解器
ANTYPE,4                    !选择瞬态求解                 
TRNOPT,FULL
LUMPM,0
ESEL,S,MAT,,1               !选择混凝土体的所有结点
NSLE,S  
IC,all,TEMP,32        !混凝土的浇筑温度为23.6
ESEL,S,MAT,,2               !选择基础土体的所有结点
NSLE,S                     
IC,all,TEMP,20             !基础的初始温度为16
ALLSEL
SFA,4,1,CONV,280, %QIWEN%  !四周木模对流系数616.32
SFA,5,1,CONV,280, %QIWEN%
SFA,9,1,CONV,280, %QIWEN%
SFA,11,1,CONV,280, %QIWEN%
SFA,2,1,CONV,600, %QIWEN%  !上面裸露对流系数616.32
BFV,1,HGEN, %HEATS%           !输入热生成率heats
ALLSE                         !选择所有
TRNOPT,FULL                   !设置计算方法
*DO,DAY,0.5,28,0.5            !设置循环计算
TIME,DAY                      !设置荷载步时间
AUTOTS,0                      !关闭自动时间步长
KBC,0                         !设置阶梯荷载步
TSRES,ERASE                   !Defines an array of keytimes at which the time-stepping strategy changes
OUTRES,ALL,ALL                !输出所有项
SOLVE                         !开始计算
*ENDDO                        !终止循环计算
/post1
plnsol,temp
set,20
/replot

!结构应力分析
/PREP7                                !前处理
ETCHG,TTS                             !转换热分析单元为结构单元
  
  !输入参数
!定义变化弹模EXX
*DIM,EXX,ARRAY,30
EXX(1,1,1) = 0.304e10   
EXX(2,1,1) = 0.608e10   
EXX(3,1,1) = 1.216e10   
EXX(4,1,1) = 1.824e10   
EXX(5,1,1) = 2.432e10   
EXX(6,1,1) = 3.04e10   
EXX(7,1,1) = 3.144e10   
EXX(8,1,1) = 3.248e10   
EXX(9,1,1) = 3.352e10   
EXX(10,1,1) = 3.456e10  
EXX(11,1,1) = 3.506e10  
EXX(12,1,1) = 3.556e10  
EXX(13,1,1) = 3.606e10  
EXX(14,1,1) = 3.656e10  
EXX(15,1,1) = 3.663e10  
EXX(16,1,1) = 3.67e10  
EXX(17,1,1) = 3.677e10  
EXX(18,1,1) = 3.684e10  
EXX(19,1,1) = 3.691e10  
EXX(20,1,1) = 3.702e10
EXX(21,1,1) = 3.709e10   
EXX(22,1,1) = 3.716e10   
EXX(23,1,1) = 3.723e10   
EXX(24,1,1) = 3.73e10   
EXX(25,1,1) = 3.737e10   
EXX(26,1,1) = 3.743e10   
EXX(27,1,1) = 3.75e10   
EXX(28,1,1) = 3.757e10   
EXX(29,1,1) = 3.759e10   
EXX(30,1,1) = 3.760e10
!定义30种材料每0.5d一种材料来反映弹模的改变
*DO,I,1,30,1                                   
    MP,DENS,I,2500       !混凝土 密度为2500
    MP,EX,I,EXX(I)        !混凝土 弹模为EXX(I)
    MP,PRXY,I,0.167       !混凝土 泊松比为0.167
    MP,ALPX,I,0.75E-5     !混凝土 热膨胀系数为0.8e-5
    MP,REFT,I,15          !定义参考温度
*ENDDO
VCLEAR, 2, 2, 1           !清除基础土体网格         
VDELE, 2, 2, 1, 1         !删除基础土体
ALLSEL                    !选择选择所有

/SOLU                     !进入求解器
NROPT,FULL                !Specifies the Newton-Raphson options in a full transient analysis
NLGEOM,ON                 !考虑大变形的影响

ASEL,S,,,1,1,1            !选择面1
NSLA,S,1                  !选择面1上的结点
D,all,ALL                 !加基础全约束
ALLSEL                    !选择选择所有
DA,3,SYMM                 !加对称约束
DA,6,SYMM                 !加对称约束
DA,4,SYMM                 !加对称约束
DA,5,SYMM
DA,9,SYMM                 !加对称约束
DA,11,SYMM
DA,12,SYMM                 !加对称约束
DA,13,SYMM
DA,14,SYMM                 !加对称约束
DA,15,SYMM
ALLSEL
I=1
*DO,DAY,0.5,10,0.5
    ESEL, S, MAT, , I, I,1                    !选择混凝土材料号为I的元素
    ACEL,0,0,-9.8,                            !GRAVITY LOAD 施加重力荷载加速度
    LDREAD,TEMP,,,DAY, ,'themal','rth',' '    !TEMP LOAD Reads results from the results file and applies them as loads.
    TIME,DAY
    AUTOTS,0
    KBC,0                                     !设置阶梯荷载步
    TSRES,ERASE                               !Defines an array of keytimes at which the time-stepping strategy changes
    OUTRES,ALL,ALL,                           !输出所有的结果
    SOLVE                                     !计算
    I=I+1
    *IF,I,LE,20,THEN
        MPCHG,I,ALL
    *ENDIF
*ENDDO

caoer

感觉你定义的30种材料都没用上么