混凝土温控问题，小女子程序求查看
各位大神，小妹做的温控程序，可以循环下去，但是总是出不来想要的结果，加了绝热边界，对流，水泥水化热，可是有时候数会很大，有时候又会出现负数，不知道哪位能给指导下，我把命令流在下方，大坝的高度是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
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 怎么没有人帮忙看看呢
小女子不甚感激
555555 呜呜~~~~(>_<)~~~~ 没有人帮助么
:'(:'(:'(:'( 555555
55555555555555 :'(:'(:'(为什么么有人帮我
找到问题所在了，符号大写了
很强，我一直在找混凝土浇筑过程的东西，我做的东西，温度应力算不出来，正好学习了。
页:
[1]