求救？升温求解，怎么算出确实降温来？？瞬态温度场模拟

cfrfosh

/filename,wendu
/title,the anaslsys of xiangliang thermal
/config,nres,100000
/prep7                        !进入处理器
         
ET,1,solid70                  !定义单元类型
mp,ex,1,3.50e10               !指定混凝土材料与温度有关的属性（弹性模量）
mp,dens,1,2500                !混凝土材料属性（ 质量密度）
mp,C,1,0.880                  !混凝土材料属性（ 比热容）
mp,KXX,1,3.3                  !混凝土材料属性（ 热导率）
                     
k,1,0,0                   !建立半截面模型（关键点）
k,2,1,0
k,3,1,1
k,4,0,1
k,11,0,0,2                  !建立半截面模型（关键点）
k,12,1,0,2
k,13,1,1,2
k,14,0,1,2
v,1,2,3,4,11,12,13,14
allsel
mat,1
vsweep,1
asel,s,loc,y,0                !选择顶板顶面0
cm,mian1,area                 !生成顶板顶面元件                 

asel,s,loc,y,1                !选择顶板顶面0
cm,mian2,area
asel,s,loc,x,1                !选择顶板顶面0
cm,mian3,area
asel,s,loc,z,0                !选择顶板顶面0
cm,genbu,area               
allsel
*dim,a1,arrab,13        !!!!!!!!!定义综合温度下顶板数组为一维整数列
a1(1)=28
a1(2)=29
a1(3)=31
a1(4)=33
a1(5)=34
a1(6)=35
a1(7)=37
a1(8)=37
a1(9)=37
a1(10)=37
a1(11)=36
a1(12)=36
a1(13)=34
*dim,a2,arrab,13        
a2(1)=28
a2(2)=29
a2(3)=31
a2(4)=33
a2(5)=34
a2(6)=35
a2(7)=37
a2(8)=37
a2(9)=37
a2(10)=37
a2(11)=36
a2(12)=36
a2(13)=34
*dim,a3,arrab,13        
a3(1)=28
a3(2)=29
a3(3)=31
a3(4)=33
a3(5)=34
a3(6)=35
a3(7)=37
a3(8)=37
a3(9)=37
a3(10)=37
a3(11)=36
a3(12)=36
a3(13)=34

*dim,xwdq,arrab,13       !!!!!!!!!定义箱外大气平均温度变化
xwdq(1)=32
xwdq(2)=32
xwdq(3)=33
xwdq(4)=34
xwdq(5)=34
xwdq(6)=34
xwdq(7)=34
xwdq(8)=35
xwdq(9)=35
xwdq(10)=35
xwdq(11)=35
xwdq(12)=35
xwdq(13)=35
/solu                       !加载
antype,transient,new       !做一次新的瞬态分析
timint,off                 !先作稳态(静态)分析,确定初始条件
time,0.01                  !设定只有一个子步的时间很小的载荷步
deltim,0.01                !设定时间步长
allsel                      !选择全部
sfa,genbu,,hflux,0
                    
               !选择下顶板元件
sfa,mian1,,CONV,9.4,30
da,mian1,temp,28               
sfa,mian2,,CONV,19.6,30
da,mian2,temp,28               
sfa,mian3,,CONV,21.1,30
da,mian3,temp,28
              
                    
SOLVE                       !得到初始温度分布


sfa,genbu,,hflux,0
*do,i,1,13,1                !!!用do循环来进行瞬态求解
time,i                       !设置每一步荷载结束时的时间
timint,on                    !打开时间积分效果
                       
allsel                        !删除前一阶段边界条件
               !选择下顶板元件
sfa,mian1,,CONV,9.4,xwdq(i)     !在线上施加对流面荷载（19.6为膜系数，a1(i)为平均温度）
da,mian1,temp,a1(i)         
sfa,mian2,,CONV,19.6,xwdq(i)
da,mian2,temp,a2(i)         
sfa,mian3,,CONV,21.1,xwdq(i)
da,mian3,temp,a3(i)         
allsel
TRNOPT,FULL                  !指定完全法为瞬态分析方法
deltim,0.05,0.005,0.1       !使用自动时间跟踪设定时间步长（0.05为子步开始时间，0.005为最小时间步，0.1为最大）
autots,on                   !瞬态分析自动时间步长跟踪
KBC,0
outres,all,1                !将荷载步每隔1个子步写入到数据库中
sftran                      !转换面荷载
SOLVE                       !求解
sfadele,all,all
ddele,all,temp
*enddo
本命令流为模拟一个长方体，一面固结及hflux=0，其它边给它一个升温的边界条件，进行瞬态的分析，却计算出固结端的温度是逐渐变小的？为什么？

cfrfosh

怎么没人回复啊！！求救啊

wxpand2004

运行了，帮不了你！

tonnyw

3# wxpand2004

Why not run a simple two-dimensional model to make sure nothing is messed up?