混凝土温控问题，小女子程序求查看

wb86740

各位大神，小妹做的温控程序，可以循环下去，但是总是出不来想要的结果，加了绝热边界，对流，水泥水化热，可是有时候数会很大，有时候又会出现负数，不知道哪位能给指导下，我把命令流在下方，大坝的高度是75m，一层浇筑2.5m，共30层，循环应该没问题，设置什么的 我也不太清楚，可是能循环，模拟浇筑过程，就是结果温度云图不对，求帮助，小女子不甚感激


/solu
nropt,full

allsel,all
nsel,s,loc,y,0,75       ！大坝的高度是75M
esln,s,1
ekill,all               !杀死所有的混凝土元素   

allsel,all
nsel,s,loc,y,-50,0        ！选出地基包含的节点
ic,all,temp,20        !地基的初始温度20度

allsel,all
nsel,s,loc,y,0,75         ！选出上层混凝土包含的节点
ic,all,temp,28         !混凝土的浇筑温度28度

allsel,all
asel,s,loc,y,0
asel,u,loc,x,0,63.5
sfa,all,,conv,2016,25     !地基表面与大气接触部分加对流
allsel,all
asel,s,loc,x,-50
asel,a,loc,x,130
asel,a,loc,z,0
asel,a,loc,z,20.5
asel,a,loc,y,-50
nsla,s,1                  !混凝土俩测基岩左右五面加绝热边界条件
sf,all,hflux,0

*do,layer,1,30
  allsel,all
  vsel,s,loc,y,2.5*(layer-1),2.5*layer      ！2.5m一层，选出浇筑层
  nslv,s,1
  esln,s,1
  ealive,all                       ！激活浇筑层
  allsel,all

*do,day,layer*3-2,layer*3,1

wb86740

allsel,all
    vsel,s,loc,y,(layer-1)*2.5,layer*2.5
    aslv,s,1
    asel,u,loc,z,0
    asel,u,loc,z,20.5
    asel,u,loc,x,5
    aSEL,u,LOC,Y,(LAYER-1)*2.5        ！通过布尔运算选出浇筑层与空气接触的部分
    nsla,s,1
    SF,ALL,,CONV,2016,25          !加混凝土第三类对流边界条件，气温恒定25

    *DO,LOOP,1,LAYER,1
      ALLSEL,ALL
      vSEL,s,LOC,Y,(LOOP-1)*2.5,LOOP*2.5
      nslv,s,1
      esln,s,1
      TDAY=DAY-(LOOP-1)*3         !已经完工的天数
      HE00=47880*(TDAY/(0.862+TDAY)-(TDAY-1)/(0.862+TDAY-1))
      bfe,all,hgen,,he00           !选中浇筑层，循环 加水化热
    *ENDDO
       
  trnopt,full                !瞬态分析选项
  lumpm,0                     ！块矩阵计算
  time,day                     ！时间步长
  autots,0                   ！自动步长
  deltim,1,,,1              ！确定荷载步的时间步长
  kbc,0                      ！荷载步递增方式
  tsres,erase                ！定义关键时间点
  outres,all,all            ！结果文件保存内容
  allsel,all
  solve

*enddo  
  allsel,all
  asel,s,loc,y,2.5*layer           ！选择出大坝的浇筑块上层与空气接触的地方
  nsla,s,all                       
  sfadele,all,,conv           ！浇筑新层前删除上层的对流荷载
  sfdele,all,conv
  allsel,all        

/post1  
  nsel,s,loc,y,-50,layer*2.5
  esln,s,1
  eplot
  plnsol,temp                 ！显示温度云图

/image,save,i,bmp          ！保存图片         
/solu  
*enddo
save
/post1
antime,90,0.5,,1,2,0,90          ！生成动画
finish

wb86740

怎么没有人帮忙  看看呢
  小女子不甚感激

wb86740

555555 呜呜~~~~(>_<)~~~~     没有人帮助么

wb86740

:'(:'(:'(

wb86740

555555
55555555555555

wb86740

:'(:'(为什么么有人帮我

wb86740

找到问题所在了，符号大写了

无奈和孤独的人

很强，我一直在找混凝土浇筑过程的东西，我做的东西，温度应力算不出来，正好学习了。