zzluck 发表于 2012-3-15 15:56:27

求助:施加水化热后温度升到一二百度!!!!!!!!!

finish
/config,nres,1000000
/prep7
*dim,con_time,array,32         !定义数组每个浇筑仓施工时间为7天,共32仓
*do,i,1,31,1
con_time(i)=7
*enddo
con_time(32)=63                !7天施工,56天释放期!共31*7+63=217+63=280天
con_temp=10                  !浇筑温度

/solu
Allsel,all
nsla,s,1
ddele,all,temp

nropt,full
lsclear,all                !从数据库中清除荷载及荷载步选项
Allsel,all

Allsel,all
esel,s,mat,,3             !选择基岩
nsle,s
ic,all,temp,14.3
Allsel,all



!混凝土浇筑时温度
allsel,all
esel,s,mat,,1
esel,a,mat,,2
nsle,s
ic,all,temp,con_temp
ekill,all                  !杀死混凝土单元

allsel,all
day=0                               !赋值day=0
!**************************开始计算********************************************
*do,i,1,32,1
    !实际仿真天数!!
    !激活第i层单元!!
allsel,all
cmsel,s,jzct%i%
eslv,s
ealive,all
ALLSEL,BELOW,elem
DDELE,ALL,temp

*do,ty,1,con_time(i),1    !!!con_time(i)为每仓混凝土的浇筑时间!!!!!!!!      
      day=day+1            

!删对流边界条件
allsel,all
nsel,s,,,all
sfdele,all,conv

!删绝热边界条件!!!!!!!!!
       allsel,all
       nsel,S, , ,all
       sfdele,all,HFLUX

!加基岩的绝热边界条件
Allsel,all
Cmsel,s,yjsm                            !四面基础
Cmsel,a,syjysm
Cmsel,a,xyjysm

cmsel,a,xyjym
nsla,s,1
Sf,all,hflux,0                                !(热对流密度为0)
Allsel,all

!加岩石的第三边界条件
Allsel,all
Cmsel,s,syjym
*if,i,lt,19,then
cmsel,a,jyjhm%i%
*endif
nsla,s
Sf,all,conv,2313.6,14.3+13.4*cos(6.28*(330+day-210)/365)      !按5m/s风速计算放热系数

!加混凝土的第三边界条件
allsel,all
cmsel,s,jzcm%i%
nsla,s
Sf,all,conv,2313.6,14.3+13.4*cos(6.28*(330+day-210)/365)

allsel,all
t=day
      *do,j,1,i,1         
      *if,j,eq,1,then                     !j=1
      t=t
      *else
      t=t-con_time(j-1)
      *endif

      Cmsel,s,jzct%j%
      Eslv,s
      nsle,s
      HE00=270*190*(t/(0.91+t)-(t-1)/(0.91+t-1))
      BFE,ALL,HGEN, ,HE00 !加水化热
      allsel,all
      *enddo
allsel,all

      ANTYPE,4
      trnopt,full
      LUMPM,0            
      AUTOTS,0
      DELTIM,1   !指定子步步长
      KBC,0
      NROPT,FULL      !Use full Newton-Raphson.
      time,day
      tsres,erase
      OUTPR,NSOL,last
      OUTRES,NSOL,last
      solve
    *enddo               
*enddo               


为了验证问题所在,我把水化热删了又算了一遍,最高温度30左右,怎么加水化热温度升这么高啊。
jzct指浇筑层体! 希望高手回答! 万分感谢!!!!!!!
页: [1]
查看完整版本: 求助:施加水化热后温度升到一二百度!!!!!!!!!