[原创]双曲线形式水化热模拟承台浇筑

luckprown82

采用双曲线形式水化热模拟承台浇筑，比较复合指数及指数型水化热模式，计算结果更接近实测数据。命令流提供大家，多有不当之处，望大家指正。部分命令流。
!以小时计算
!双曲线形式水化热
FINISH  
/CLEAR
/config,nproc,2                  !双核处理
/COM,  Thermal  
/filename,Cross-JHYH-17#
/PREP7
!----------------------定义热分析材料及网格划分----------------------------------
k,1,29.0/2,13.75/2,0
k,2,29.0/2,0,0
k,3,0,0,0
k,4,0,13.75/2,0
k,11,29.0/2,13.75/2,-3.5
k,12,29.0/2,0,-3.5
k,13,0,0,-3.5
k,14,0,13.75/2,-3.5
v,1,2,3,4,11,12,13,14
k,101,29.0,13.75,-3.5
k,102,29.0,0.0,-3.5
k,103,0.0,0.0,-3.5
k,104,0.0,13.75,-3.5
k,201,29.0,13.75,-3.5-5.0
k,202,29.0,0.0,-3.5-5.0
k,203,0.0,0.0,-3.5-5.0
k,204,0.0,13.75,-3.5-5.0
v,101,102,103,104,201,202,203,204
!--------------------------------建立1/4模型
snyl=303+0.25*82                 !单位体积水泥用量，单位kg/m3
XSA=0.69                         !复合指数系数a=
XSB=0.56                         !复合指数系数b=
XSN=0.862                        !双曲线系数n=
hntmd=2326.0                     !混凝土材料密度，单位kg/m3
hntdrxs=9.440                    !混凝土导热系数，国际单位W/(m.℃)，单位kJ/(m.h.℃)
hntbr=0.916                      !混凝土比热，，国际单位J/(kg.℃)单位kJ/(kg.℃)
trdrxs=2.0*3.6                   !土壤导热系数，国际单位W/(m.℃)，单位kJ/(m.h.℃)
trbr=1.80                        !土壤比热，国际单位J/(kg.℃)，单位kJ/(kg.℃)
jycswd=7.2                       !基岩施加初始温度℃
jydlxs=67.4                      !基岩表面与空气的对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
hntjzwd=14.6                     !混凝土浇筑时温度℃
cmdlxs=37.335                    !侧模对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
jzdlxs=67.4                      !浇筑顶面对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
yhdlxs=34.257                    !养护顶面对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
q0=330                           !水泥水化热常数,单位kJ/kg
hmm=0.340/24                     !常数（1/d）,浇筑温度确定5-0.295  10-0.318  15-0.34  20-0.362  25-0.384
ET,1,70
MP,DENS,1,hntmd                  !混凝土材料密度，单位kg/m3  
MP,KXX ,1,hntdrxs                !混凝土导热系数，国际单位W/(m.℃)，单位kJ/(m.h.℃)            
MP,C,1,hntbr                     !混凝土比热，国际单位J/(kg.℃)，单位kJ/(kg.℃)               
MP,DENS,2,2300.0                 !基岩材料
MP,KXX ,2,trdrxs
MP,C,2,trbr
LSEL,S,line,,5
LESIZE,ALL,,,7
LSEL,S,line,,1,
LESIZE,ALL,0.275,,
LSEL,S,line,,2
LESIZE,ALL,0.25,,
allsel
type,1
mat,1
vmesh,1,
LSEL,S,line,,17
LESIZE,ALL,1.0,,
LSEL,S,line,,13,
LESIZE,ALL,0.275,,
LSEL,S,line,,14
LESIZE,ALL,0.25,,
allsel
type,1
mat,2
vmesh,2
allsel
nummrg,all,  
numcmp,all
!***********定义环境温度及时间列表***********   
lrow=98
lcolumn=1
lplane=1
*dim,air_temp17,,lrow,lcolumn                                                              
*vread,air_temp17(1),air_temp17,dat,,ijk,lrow                                       
(f6.1)  
*dim,time_t17,,lrow,lcolumn                                                              
*vread,time_t17(1),time_t17,dat,,ijk,lrow                                       
(f6.1)  
allsel
NSEL,s,LOC,Z,-3.5,-3.0,
NSEL,a,LOC,Z,-3.1,-2.9,
ESLN,S,1,ALL           
CM,E_LAYER1,ELEM         
NSEL,S,LOC,Z,-3.0,-2.5   
NSEL,a,LOC,Z,-2.6,-2.4
ESLN,S,1,ALL              
CM,E_LAYER2,ELEM   
NSEL,S,LOC,Z,-2.5,-2.0,
NSEL,a,LOC,Z,-2.1,-1.9   
ESLN,S,1,ALL              
CM,E_LAYER3,ELEM
NSEL,S,LOC,Z,-2.0,-1.5,
NSEL,a,LOC,Z,-1.6,-1.4   
ESLN,S,1,ALL              
CM,E_LAYER4,ELEM     
NSEL,S,LOC,Z,-1.5,-1.0,
NSEL,a,LOC,Z,-1.1,-0.9   
ESLN,S,1,ALL              
CM,E_LAYER5,ELEM  
NSEL,S,LOC,Z,-1.0,-0.5,
NSEL,a,LOC,Z,-0.6,-0.4   
ESLN,S,1,ALL              
CM,E_LAYER6,ELEM  
NSEL,S,LOC,Z,-0.5,0.0,
NSEL,a,LOC,Z,-0.1,0.1   
ESLN,S,1,ALL              
CM,E_LAYER7,ELEM   
         
!--------------------------------------------热分析
/SOLU
NROPT,FULL                         !在瞬态热分析中选择FULL N-R项  Specifies the Newton-Raphson options in full transient analysis
ESEL,S,MAT,,2                      !选择基岩材料
NSLE,S                             !选择单元对应节点
IC,all,TEMP,jycswd,                !基岩施加初始温度
!
asel,s,loc,x,29.0
asel,a,loc,y,13.75
asel,a,loc,z,-3.5-5.0
NSLA,S,1
SF,all,HFLUX,0                     !加基岩绝热边界条件
ALLSEL,ALL                         !加初始条件
ESEL,U,MAT,,2
NSLE,S
EKILL,ALL                          !杀死混凝土
IC,all,TEMP,hntjzwd,               !混凝土浇筑时温度hntjzwd
ANTYPE,TRANSIENT,new               !分析类型
!!----------------------循环模拟混凝土浇筑过程----------------------------------
!****************************浇至第一层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(1)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,s,LOC,Z,-3.5,-3.0,
NSEL,a,LOC,Z,-3.1,-2.9,
ESLN,S,1,ALL
EALIVE,ALL
esel,s,live
nsle,s,all
nsel,inve
d,all,temp,hntjzwd
ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,-3.0
NSEL,a,LOC,Z,-3.1,-2.9,
SF,ALL,CONV,cmdlxs,air_temp17(1)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(1)-time_t17(1))/24)**2/24                 !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(1)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇至第二层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(2)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,S,LOC,Z,-3.0,-2.5   
NSEL,a,LOC,Z,-2.6,-2.4
ESLN,S,1,ALL            
EALIVE,ALL
ddele,all
esel,s,live
nsle,s,all
nsel,inve
d,all,temp,hntjzwd
ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,-2.5
NSEL,a,LOC,Z,-2.6,-2.4
SF,ALL,CONV,cmdlxs,air_temp17(2)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(2)-time_t17(1))/24)**2/24               !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
allsel
esel,s,ENAME,,E_LAYER2        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(2)-time_t17(2))/24)**2/24                !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热                 

ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(2)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1,
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇至第三层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(3)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,S,LOC,Z,-2.5,-2.0,
NSEL,a,LOC,Z,-2.1,-1.9
ESLN,S,1,ALL            
EALIVE,ALL
ddele,all
esel,s,live
nsle,s,all
nsel,inve
d,all,temp,hntjzwd

ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,-2.0,
NSEL,a,LOC,Z,-2.1,-1.9  
SF,ALL,CONV,cmdlxs,air_temp17(3)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(3)-time_t17(1))/24)**2/24                !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
allsel
esel,s,ENAME,,E_LAYER2        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(3)-time_t17(2))/24)**2/24              !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热                 
allsel
esel,s,ENAME,,E_LAYER3        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(3)-time_t17(3))/24)**2/24             !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     

ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(3)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1,
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇至第四层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(4)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,S,LOC,Z,-2.0,-1.5,
NSEL,a,LOC,Z,-1.6,-1.4   
ESLN,S,1,ALL            
EALIVE,ALL
ddele,all
esel,s,live
nsle,s,all
nsel,inve
d,all,temp,hntjzwd

ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,-1.5,
NSEL,a,LOC,Z,-1.6,-1.4,
SF,ALL,CONV,cmdlxs,air_temp17(4)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(4)-time_t17(1))/24)**2/24              !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
allsel
esel,s,ENAME,,E_LAYER2        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(4)-time_t17(2))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热                 
allsel
esel,s,ENAME,,E_LAYER3        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(4)-time_t17(3))/24)**2/24                 !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER4        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(4)-time_t17(4))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     

ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(4)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1,
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇至第五层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(5)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,S,LOC,Z,-1.5,-1.0,
NSEL,a,LOC,Z,-1.1,-0.9  
ESLN,S,1,ALL            
EALIVE,ALL
ddele,all
esel,s,live
nsle,s,all
nsel,inve
d,all,temp,hntjzwd

ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,-1.0,
NSEL,a,LOC,Z,-1.1,-0.9,
SF,ALL,CONV,cmdlxs,air_temp17(5)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(5)-time_t17(1))/24)**2/24               !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
allsel
esel,s,ENAME,,E_LAYER2        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(5)-time_t17(2))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热                 
allsel
esel,s,ENAME,,E_LAYER3        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(5)-time_t17(3))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER4        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(5)-time_t17(4))/24)**2/24              !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER5        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(5)-time_t17(5))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热   
ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(5)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1,
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇至第六层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(6)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,S,LOC,Z,-1.0,-0.5,
NSEL,a,LOC,Z,-0.6,-0.4   
ESLN,S,1,ALL            
EALIVE,ALL
ddele,all
esel,s,live
nsle,s,all
nsel,inve
d,all,temp,hntjzwd
ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,-0.5,
NSEL,a,LOC,Z,-0.6,-0.4,
SF,ALL,CONV,cmdlxs,air_temp17(6)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(6)-time_t17(1))/24)**2/24               !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
allsel
esel,s,ENAME,,E_LAYER2        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(6)-time_t17(2))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热                 
allsel
esel,s,ENAME,,E_LAYER3        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(6)-time_t17(3))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER4        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(6)-time_t17(4))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER5        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(6)-time_t17(5))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热  
allsel
esel,s,ENAME,,E_LAYER6        !浇筑水化热                                                
esel,s,ENAME,,E_LAYER6        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(6)-time_t17(6))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00  
  
ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(6)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1,
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇至第七层*****************************************
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(7)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,S,LOC,Z,-0.5,0.0,
NSEL,a,LOC,Z,-0.1,0.1   
ESLN,S,1,ALL            
EALIVE,ALL
ddele,all
!esel,s,live
!nsle,s,all
!nsel,inve
!d,all,temp,hntjzwd
ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,0.0,
SF,ALL,CONV,cmdlxs,air_temp17(7)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,s,LOC,Z,-0.1,0.1         !施加浇筑顶面对流系数
SF,ALL,CONV,jzdlxs,air_temp17(7)           !施加所激活的顶面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(1))/24)**2/24                 !考虑已浇筑混凝土的水化热
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热
allsel
esel,s,ENAME,,E_LAYER2        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(2))/24)**2/24            !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热                 
allsel
esel,s,ENAME,,E_LAYER3        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(3))/24)**2/24           !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER4        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(4))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热     
allsel
esel,s,ENAME,,E_LAYER5        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(5))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00                                        !施加加水化热  
allsel
esel,s,ENAME,,E_LAYER6        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(6))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00  
allsel
esel,s,ENAME,,E_LAYER7        !浇筑水化热                                                
HE00=snyl*q0*XSN/(XSN+(time_t17(7)-time_t17(7))/24)**2/24               !考虑已浇筑混凝土的水化热   
!nsle,s,all
!BF,ALL,HGEN,HE00
BFE,ALL,HGEN, ,HE00  
  
ALLSEL
TRNOPT,FULL                       !瞬态分析方法：默认
LUMPM,0                           !由单元默认的质量矩阵
TIME,time_t17(7)
AUTOTS,ON
DELTIM,2,2,12,                    !time step size =2
KBC,1                             !0=ramped  loading    1=stepped  loading
TSRES,ERASE                       !???????????????????????????????????????
OUTRES,ALL,ALL,
SOLVE
ESEL,S,MAT,,1,
NSLE,S,
SFEDELE,ALL,ALL,CONV                !删除已浇筑混凝土对流边界条件
SFDELE,ALL,CONV
allsel
NSEL,S,LOC,z,-3.4,-3.6              !删除施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV      
!****************************浇筑完成，养护阶段************************************
*DO,i,8,lrow,1
allsel
NSEL,S,LOC,z,-3.4,-3.6        !施加基岩表面对流系数        
ASEL,S,,,6
NSLA,U,1
SF,ALL,CONV,jydlxs,air_temp17(i)           !加基岩第三边界条件（和空气接触面）,地基与空气对流系数为jydlxs，空气温度为ir_temp17(i)

ALLSEL,ALL   
ASEL,S,,,2                    !施加侧面对流系数
ASEL,A,,,5
NSLA,S,1      
NSEL,R,LOC,Z,-3.5,0.0,
SF,ALL,CONV,cmdlxs,air_temp17(i)           !施加所激活的侧面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)
allsel
NSEL,s,LOC,Z,-0.1,0.1         !施加浇筑顶面对流系数
SF,ALL,CONV,yhdlxs,air_temp17(i)           !施加所激活的顶面混凝土与空气间的对流边界条件，对流系数为jydlxs，空气温度为ir_temp17(i)

allsel
esel,s,ENAME,,E_LAYER1        !浇筑水化热
HE00=snyl*q0*XSN/(XSN+(time_t17(i)-time_t17(1))/24)**2/24               !考虑已浇筑混凝土的水化热

luckprown82

本来想邀请个朋友来论坛的，结果仿真币不够，于是拿出自己的命令流想挣点仿真币。看大家只来看看，也不下载，想想大家挣币也不容易。还是取消售价，大家免费下载，觉得对你有帮助了，留个言；有不正确的地方，给小弟指正，大家共同进步。谢谢大家。

ivanjfj

学习了。。。。。。。。

sdrzlz

谢谢分享，你的参与至少能让我进步，再次表示感谢！

mybq

写得好复杂啊

739898210

谢谢楼主

moorecsu

非常感谢，新手先学习~

moorecsu

楼主，看不懂这个HE00=snyl*q0*XSN/(XSN+(time_t17(2)-time_t17(1))/24)**2/24               !考虑已浇筑混凝土的水化热。。。能不能指教一下。。还有那两组空气温度和时间的数据也看不懂，求指教~非常感谢。