请教:混凝土热分析!附命令流

tjrenxiaofeng

在如下命令流中出现如下错误:
Surface loads on unselected area *** are not transferred to elements.
急求达人指正!
也请各位提出改进之处!
万分感激!
!!!!!!!
fini
/cle
/filname,the-analysis-of-an-crossing
/title,thermal-stress-analysis
/com,thermal
*set,day,40
*set,soil_on,30                        
!前处理
/prep7   
et,1,solid70
et,2,mesh200,6   
         
mp,dens,1,2333.8
mp,kxx,1,281.5656           
mp,c,1,1.047                             
!建立模型
k,1,0,0,0
k,2,0.1,0,0
k,3,0.3,0,0
k,4,0.5,0,0
k,5,0.77,0,0
k,6,0.95,0,0
k,7,1.15,0,0
k,8,1.3,0,0
k,59,8.95,0,0
kfill,8,59,50, , ,1,
k,60,9.11,0,0
k,61,9.4,0,0
k,62,9.6,0,0

allsel
ksel,all
kgen,2,all,,,0,0.1,0,,0

allsel
ksel,s,loc,y,0.1
kgen,2,all,,,0,0.25,0,,0

allsel
ksel,s,loc,y,0.1
kgen,2,all,,,0,0.5,0,,0

k,249,0,0.85,0
k,250,0.1,0.85,0
k,251,0.3,0.85,0
k,252,0.5,0.85,0
k,253,0.77,0.85,0

ksel,s,,,6,59
kgen,2,all,,,0,0.7,0,,0

allsel
k,308,9.11,0.9,0
k,309,9.4,0.9,0
k,310,9.6,0.9,0

k,311,0,1.0,0
k,312,0.1,1.0,0
k,313,0.3,1.0,0
k,314,0.5,1.0,0
k,315,0.6,1.0,0
allsel
ksel,s,,,311,315
kgen,2,all,,,0,0.15,0,,0

k,321,9.3,1.15,0
k,322,9.4,1.15,0
k,323,9.6,1.15,0

allsel
ksel,s,loc,y,1.15
kgen,2,all,,,0,0.2,0,,0     

allsel
ksel,s,loc,y,1.35
kgen,16,all,,,0,0.25,0,,0   

allsel
ksel,s,loc,y,5.1
kgen,2,all,,,0,0.2,0,,0     

k,460,0,5.5,0
k,461,0.1,5.5,0
k,462,0.3,5.5,0
k,463,0.5,5.5,0
k,464,0.6,5.5,0
k,465,0.7,5.525,0
k,466,0.85,5.5625,0
k,467,1.0,5.6,0
k,468,1.15,5.6375,0
k,469,1.3,5.675,0
k,470,1.45,5.7125,0
k,471,1.6,5.75,0
k,472,1.8,5.8,0
k,473,1.9,5.8,0
kgen,44,473,,,0.15,0,0,,0

k,517,8.5,5.75263,0
k,518,8.65,5.70526,0
k,519,8.8,5.65790,0
k,520,8.95,5.61053,0
k,521,9.1,5.56316,0
k,522,9.3,5.5,0
k,523,9.4,5.5,0
k,524,9.6,5.5,0

k,525,0,5.9,0
k,526,0.1,5.9,0
k,527,0.3,5.9,0
k,528,0.5,5.9,0
k,529,0.6,5.9,0
k,530,0.7,5.9,0
kgen,7,530,,,0.15,0,0,,0
k,537,1.8,5.9,0
k,538,1.9,5.9,0
kgen,49,538,,,0.15,,0,,0

k,587,9.3,5.9,0
k,588,9.4,5.9,0
k,589,9.6,5.9,0

allsel
ksel,s,loc,y,5.9
kgen,2,all,,,0,0.2,0,,0      

allsel
ksel,s,loc,y,6.1
kgen,3,all,,,0,0.25,0,,0     

allsel
ksel,s,loc,y,6.6
kgen,2,all,,,0,0.1,0,,0      
/PNUM,KP,1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!   
allsel
*do,ij,1,61
a,ij,ij+1,ij+63,ij+62
*enddo
*do,ij,63,123
a,ij,ij+1,ij+63,ij+62
*enddo
*do,ij,125,185
a,ij,ij+1,ij+63,ij+62
*enddo
*do,ij,187,247
a,ij,ij+1,ij+63,ij+62
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,ij,249,252
a,ij,ij+1,ij+63,ij+62
*enddo
*do,ij,311,314
a,ij,ij+1,ij+6,ij+5
*enddo
*do,ii,0,136,8
  *do,ij,316,319
  a,ii+ij,ii+ij+1,ii+ij+9,ii+ij+8
  *enddo
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,ij,308,309
a,ij,ij+1,ij+14,ij+13
*enddo
*do,ii,0,128,8
  *do,ij,321,322
  a,ii+ij,ii+ij+1,ii+ij+9,ii+ij+8
  *enddo
*enddo
*do,ij,457,458
a,ij,ij+1,ij+66,ij+65
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,ii,0,260,65
  *do,ij,460,523
  a,ii+ij,ii+ij+1,ii+ij+66,ii+ij+65
  *enddo
*enddo

allsel,all
aglue,all
!划分混凝土网格!!
asel,all
lsla,s
lesize,all,,,1  
asel,all
type,2
mat,1
amesh,all
allsel
asel,all
extopt,esize,64,0
extopt,aclear,1,        
type,1
mat,1
vext,all,,,0,0,9.6,,,
nummrg,all
numcmp,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*dim,heatgen,array,1                                
*cfopen,gen,txt                                          
*do,t,1,day                                                
heat=400*460*0.8*exp(-0.8*t)              
*vwrite,heat
(f15.5)
*enddo
*cfclos
*vread,heatgen,gen,txt, , , , , ,
(f15.5)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*dim,air_temp,array,1
*cfopen,airtemp,txt
*do,t,1,day
temper=17+20*cos(0.01721*(t-8))           
*vwrite,temper
(f15.5)
*enddo
*cfclos
*vread,air_temp,airtemp,txt, , , , , ,
(f15.5)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*dim,under_temp,array,1
*cfopen,undertemp,txt
*do,t,1,day
temper_u=17+15*cos(0.01721*(t-8))            
*vwrite,temper_u
(f15.5)
*enddo
*cfclos
*vread,under_temp,undertemp,txt, , , , , ,
(f15.5)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/SOLU                              
ANTYPE,4                              
NROPT,FULL
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL
ESEL,S,MAT,,1                  
NSLE,S  
IC,all,TEMP,32                 
allsel
*do,t,1,day
time,t
timint,on
*vread,heatgen,gen,txt, , , , , ,t-1
(f15.5)
vsel,all
bfv,all,HGEN,,heatgen
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
allsel  
*vread,air_temp,airtemp,txt, , , , , ,t-1
(f15.5)
* if,t,ge,soil_on,then                           
asel,s,loc,y,6.7                          
SFA,all,,CONV,926.4, air_temp   
*else                                                
asel,s,loc,y,6.7
SFA,all,,CONV,600, air_temp     
*endif
allsel
asel,s,loc,y,0
SFA,all,,CONV,40, air_temp     
allsel
asel,s,loc,x,0
SFA,5,,CONV,926.4, air_temp      
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
allsel
asel,s,loc,y,1.0,5.5
asel,r,loc,x,0.6
SFA,all,,CONV,926.4, under_temp   
allsel
asel,s,loc,y,1.15,5.5
asel,r,loc,x,9.3
SFA,all,,CONV,926.4, under_temp   
allsel
asel,s,loc,x,0.95,8.95
asel,r,loc,y,0.7
SFA,all,,CONV,926.4, under_temp   
allsel
asel,s,loc,x,1.8,8.35
asel,r,loc,y,5.8
SFA,all,,CONV,926.4, under_temp   
allsel
asel,s,area,,1311
asel,a,area,,1491                                         
SFA,all,,CONV,926.4, under_temp   
allsel
asel,s,area,,1473
asel,a,area,,1743                                       
SFA,all,,CONV,926.4, under_temp
allsel
asel,s,area,,1887
asel,a,area,,1891
asel,a,area,,1895
asel,a,area,,1899
asel,a,area,,1903
asel,a,area,,1907
asel,a,area,,1911
asel,a,area,,1915                                       
sfa,all,,conv,926.4, under_temp        
allsel
asel,s,area,,2095
asel,a,area,,2099
asel,a,area,,2103
asel,a,area,,2107
asel,a,area,,2111
asel,a,area,,2115                                       
SFA,all,,conv,926.4, under_temp         
allsel
asel,s,loc,z,0
sfa,all,,conv,514.666,16                        
!!!!!!!!!!!!!!!
deltim,0.5,0.1,1.0
AUTOTS,on                                             
outres,all,last
solve
*enddo
save

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

tjrenxiaofeng

问题已搞定!

xgmSS

怎么搞定的啊？大家一起分享一下吧？

schultznj

学习先。
学习热分析命令流好资料，谢谢楼主

长风秋雁

遇到和楼主一样的问题，希望说明如何解决的哦。。。。。

lytcobo

弱弱的问一句，楼主，你的混凝土材料的传导系数，比热容和大气的对流系数在哪查的？我也做一个混凝土的模型，可是查不到相关参数