- 积分
- 0
- 注册时间
- 2012-2-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
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指浇筑层体! 希望高手回答! 万分感谢!!!!!!! |
|