大体积混凝土浇筑建模的问题，急！！@@

这是我的课题，尺寸为4.5*4.5*4的大体积混凝土底板，一次性浇筑成，四周侧面是100mm塑料保温层，底下做了一个土体的底座，整体做了1/4的模型，查看温度场结果很不理想，感觉好像是混凝土、保温层和空气的边界条件有问题，但是又找不出毛病，请各位大侠帮忙！！
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

具体的命令流和相应数据文件在附件里面，操作如下：

建立一个工作目录，在这里面运行，把huanjing.txt,wensheng.txt放在这个目录下面，然后把output文件夹，和calculate.txt放在c盘下面，运行“计算命令流.txt"

计算时间大概为半个小时，查看node temperature，感觉保温层的温度是最高的，很神奇，查看混凝土温度时间变化曲线，比如143节点的，温度到后面似乎没有怎么下降，这个趋势不对，说明热量没有很好地散出去，请问我的命令流有问题么？
期待高手解决！！！

计算命令流.txt内容如下：
fini
/clea
!各参数，数组代表的意义
!Nodesum混凝土内所有结点的总数
!adiaTemp绝热温升数据数组
!enviTemp空气温度数组
!Nodetemp混凝土内各节点的温度数组
!testnode测温点的编号数组
!!!!!!!!!!!!!!!!!!!!!!!!设置参数
*set,day,60             !计算天数
*set,hour,1440     !计算小时数
*set,a,2.25          !混凝土尺寸的一半
*set,b,4           !混凝土高度
*set,c,0.1           !保温层厚度
*set,h,5           !底层土层高度
*set,t,2             !侧壁土层厚度
!!!!!!!!!!!!!!!!!!!!!!!!
/prep7               
et,1,solid70     !使用热力学单元solid70
!混凝土热力学性质
mp,dens,1,2400     !密度
mp,kxx,1,10       !导热系数
mp,c,1,0.95        !热容
!底部土的热力学性质
mp,dens,2,2000
mp,kxx,2,5
mp,c,2,1.0
!泡沫热力学性质
mp,dens,3,200
mp,kxx,3,0.05
mp,c,3,1.0
!侧部土的热力学性质
mp,dens,4,2000
mp,kxx,4,0.5
mp,c,4,1.0
!建模
block,0,a,0,a,0,b   !混凝土
block,0,a,a,a+c,b,0   !保温层1
block,a,a+c,a,a+c,b,0   !保温层2
block,a,a+c,0,a,b,0   !保温层3
block,0,a,0,a,0,-h   !底层土体1
block,0,a,a,a+c,-h,0   !底层土体2
block,a,a+c,a,a+c,-h,0   !底层土体3
block,a,a+c,0,a,-h,0   !底层土体4
block,0,a,a+c,a+c+t,0,-h   !底层土体5
block,a,a+c,a+c,a+c+t,0,-h   !底层土体6
block,a+c,a+c+t,a+c,a+c+t,0,-h   !底层土体7
block,a+c,a+c+t,a,a+c,0,-h   !底层土体8
block,a+c,a+c+t,0,a,0,-h   !底层土体9

!粘接
vglue,all
allsel
!压缩编号
numcmp,volume
numcmp,area
allsel
!划分线段长度
lsel,s,length,,a,a      !选择长度为2.25的线
lesize,all,,,10          !划分为10份
lsel,s,length,,c,c      !
lesize,all,,,1
lsel,s,length,,t,t
lesize,all,,,4
lsel,s,length,,b,b
lesize,all,,,20
lsel,s,length,,h,h
lesize,all,,,10
allsel
!划分网格
vmesh,1
vmesh,2
vmesh,3
vmesh,4
vmesh,5
vmesh,6
vmesh,7
vmesh,8
vmesh,9
vmesh,10
vmesh,11
vmesh,12
vmesh,13
allsel
!对各单元赋材料属性
!混凝土单元
VSEL,S,,,1
ESLV,S
MPCHG,1,all
allsel
!泡沫单元
VSEL,S,,,3
VSEL,a,,,4
VSEL,a,,,5
ESLV,S  
MPCHG,3,all
allsel
!土单元
vsel,s,,,2    !土单元1
eslv,s
mpcha,2,all
allsel
!!!!!!!!!!!
vsel,s,,,6
vsel,a,,,7
vsel,a,,,11
vsel,a,,,12
vsel,a,,,13
eslv,s
mpcha,2,all
allsel
!!!!!!!!!!
vsel,s,,,8
vsel,a,,,9
vsel,a,,,10
eslv,s
mpcha,4,all
allsel
!建立x对称面约束
allsel
asel,s,loc,x,0
da,all,symm
allsel
!建立y对称面约束
asel,s,loc,y,0
da,all,symm
allsel
DTRAN
!求解
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/solu           !进入求解器
antype,transient,new  !定义分析类型
nropt,full      !牛顿－拉普拉斯全解
!添加初始与边界条件
vsel,a,,,2
vsel,a,,,6
vsel,a,,,7
vsel,a,,,8
vsel,a,,,9
vsel,a,,,10
vsel,a,,,11
vsel,a,,,12
vsel,a,,,13
nslv,s,1
ic,all,temp,10   !土的初始温度为10度
allsel
vsel,s,,,1
vsel,a,,,3
vsel,a,,,4
vsel,a,,,5
nslv,s,1
ic,all,temp,16    !泡沫和混凝土的初始温度为14.7度
allsel
!添加土体边界温度约束
asel,s,,,10
asel,a,,,33
asel,a,,,11
asel,a,,,52
asel,a,,,50
asel,a,,,55
asel,s,,,30
asel,a,,,38
asel,a,,,9
asel,a,,,54
asel,a,,,46
asel,a,,,41
asel,s,,,35
asel,a,,,49
asel,a,,,58
DA,all,TEMP,10           !土体周围边界温度恒定为10度
allsel

!取得混凝土内所有节点数目
vsel,s,,,1
nslv,s,1
*get,Nodesum,NODE,0,count
allsel
!导入绝热温升数据，并写入数组adiaTemp
*DIM,adiaTemp,array,200,1,1
!*VREAD,adiaTemp,'c:\wensheng','txt'
!(f8.2)
/input,wensheng,txt
!导入气温数据，并写入数组enviTemp
*DIM,enviTemp,array,1450,1,1
!*VREAD,enviTemp,'c:\qiwenshuju','txt'
!(f8.2)
/input,huanjing,txt
!将混凝土内各节点的温度建立数组Nodetemp
*DIM,Nodetemp,array,Nodesum
*DO,Nodenum,1,Nodesum,1
   Nodetemp(Nodenum)=16
*ENDDO
*dim,te,array,10    !中转数组，用于把测温点温度写入文件

testsum=30
*dim,Testnode,array,testsum,1,1      !将测温点对应的结点号编组
   testnode(1)=143,475,480,485,1,157,741,746,741,7
   testnode(11)=122,247,252,257,2,148,570,575,580,17
   testnode(21)=202,1777,1772,1767,81,133,466,461,456,12

*cfopen,c:\output\1~10,txt,,append
    *do,testnum,1,10,1
         *get,te(testnum),node,testnode(testnum),temp
    *enddo

        *vwrite,te(1),te(2),te(3),te(4),te(5),te(6),te(7),te(8),te(9),te(10)
          (f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4)         
*cfclos
*cfopen,c:\output\11~20,txt,,append
    *do,testnum,11,20,1
         *get,te(testnum-10),node,testnode(testnum),temp
    *enddo
         *vwrite,te(1),te(2),te(3),te(4),te(5),te(6),te(7),te(8),te(9),te(10)
          (f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4)         
*cfclos
*cfopen,c:\output\21~30,txt,,append
    *do,testnum,21,30,1
         *get,te(testnum-20),node,testnode(testnum),temp
    *enddo
         *vwrite,te(1),te(2),te(3),te(4),te(5),te(6),te(7),te(8),te(9),te(10)
          (f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4,f15.4)         
*cfclos

Nstep=0    !计算次数的初始值
!前7天时间间隔为2小时
  tm_start=1
  tm_fini=179
  pace=2
  step=(tm_fini-tm_start)/pace+1
  Nstep=Nstep+step
!时间循环计算
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  *do,tm,tm_start,tm_fini,pace
     /input,'c:\calculate','txt'
*ENDDO
!7~14天时间间隔为4小时
  tm_start=181
  tm_fini=337
  pace=4
  step=(tm_fini-tm_start)/pace+1
  Nstep=Nstep+step
!时间循环计算
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,tm,tm_start,tm_fini,pace
   /input,'c:\calculate','txt'
*ENDDO
!15~60天时间间隔为8小时
  tm_start=341
  tm_fini=1445
  pace=8
  step=(tm_fini-tm_start)/pace+1
  Nstep=Nstep+step
!时间循环计算
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,tm,tm_start,tm_fini,pace
   /input,'c:\calculate','txt'
*ENDDO

*DIM,TempResult,ARRAY,NStep,Nodesum
/post26
*DO,Nodenum,1,Nodesum,1                    !Nsum为节点总数
   NSOL,2,Nodenum,TEMP                  !选择第I个节点的各计算步的温度值
   VGET,TempResult(1,Nodenum),2         !将选择的温度值存入数组TempResult的第I列
*ENDDO
PARSAV,ALL,c:\temp_RESULT          !将包括数组TempResult在内的所有参数存入外部文件

若可以进一步讨论，请加我qq：157443921