在0号块水化热分析中定义大气温度和水化热生成怎么样定义

liudapengsnd

*set, day, 7
!计算天数
*set, time, day*24
!计算小时数
!输入外界气温表
*DIM,%qiwen %,TABLE,6,10,1
!外界气温表名为qiwen
! Begin of equation: 17.3+11.35*cos(0.01745*({TIME}-215))
% qiwen %(0,0,1)= 0.0, -999
% qiwen %(2,0,1)= 0.0
% qiwen %(3,0,1)= 0.0
% qiwen %(4,0,1)= 0.0
% qiwen %(5,0,1)= 0.0
% qiwen %(6,0,1)= 0.0
% qiwen %(0,1,1)= 1.0, -1, 0, 215, 0, 0, 1
% qiwen %(0,2,1)= 0.0, -2, 0, 1, 1, 2, -1
% qiwen %(0,3,1)= 0, -1, 0, 0.01745, 0, 0, -2
% qiwen %(0,4,1)= 0.0, -3, 0, 1, -1, 3, -2
% qiwen %(0,5,1)= 0.0, -1, 10, 1, -3, 0, 0
% qiwen %(0,6,1)= 0.0, -2, 0, 11.35, 0, 0, -1
% qiwen %(0,7,1)= 0.0, -3, 0, 1, -2, 3, -1
% qiwen %(0,8,1)= 0.0, -1, 0, 17.3, 0, 0, -3
% qiwen %(0,9,1)= 0.0, -2, 0, 1, -1, 1, -3
% qiwen %(0,10,1)= 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 17.3+11.35*cos(0.01745*({TIME}-215))
snyl=486
!单位体积水泥用量，单位kg/m3
gkdlxs=8.68

!钢模表面与空气的对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
ti=20


!混凝土浇筑时温度℃
kndlxs=6.55
!
有开口内腔内表面对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
wndlxs=4.75

!
无开口内腔内表面对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
dkdlxs=10.1
!顶面与空气接触对流系数,国际单位W/(m^2.℃)，单位kJ/(m^2.h.℃)
Q0=330
!水泥水化热常数,单位kJ/kg
!选择混凝土面加热边界条件
allsel
NSLA,S,loc,x,0
NSLA,S,loc,z,0
NSLA,S,loc,y,-6.5
SF,all,HFLUX,0
!对称面加绝热边界条件
allsel
NSLA,S,loc,y,0

SF,all,CONV, dkdlxs, %qiwen %

!
对空气接触面加第三边界条件
Allsel
NSLA,S,loc,y,0
NSLA,S,loc,z,6
NSLA,S,loc,y,-6.117,-6.5
NSLA,S,loc,x,-6.05,-3.25
NSLA,u,loc,z,0
SF,all,CONV, gkdlxs, %qiwen %
!
钢模与空气接触面加对流条件
Allsel
NSLA,S,loc,z,3,6
NSLA,u,loc,z,6
NSLA,u,loc,x,0
NSLA,u,loc,y,-6.117,-6.5
NSLA,u,loc,y,0
NSLA,u,loc,x,-3.25,-6.05
SF,all,CONV, kndlxs, %qiwen %
!
有开口内腔内表面对流条件
Allsel
NSLA,S,loc,z,0,3
NSLA,u,loc,x,0
NSLA,u,loc,z,0
NSLA,u,loc,y,0
NSLA,u,loc,x,-3.25,-6.05
NSLA,u,loc,y,-6.5


SF,all,CONV, wndlxs, %qiwen %
!
无开口内腔内表面对流条件
FINISH

!输入热生成率函数
allsel
esel,s,ENAME,,E_LAYER1

!浇筑水化热
Heoo=Q0exp{-55[0.016exp(0.039Ti)*time]-1.25}
BFE,ALL,HGEN, ,HE00
!施加加水化热
*SET,_FNCNAME,'heats'   
*DIM,%_FNCNAME%,TABLE,6,16,1
   ! Begin of equation: 2450*0.7462*227.2246*0.13*exp(-0.7426*{TIME})  
%_FNCNAME%(0,0,1)= 0.0, -999
%_FNCNAME%(2,0,1)= 0.0  
%_FNCNAME%(3,0,1)= 0.0  
%_FNCNAME%(4,0,1)= 0.0  
%_FNCNAME%(5,0,1)= 0.0  
%_FNCNAME%(6,0,1)= 0.0  
%_FNCNAME%(0,1,1)= 1.0, -1, 0, 2450, 0, 0, 0
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 0.7462, 0, 0, -1
%_FNCNAME%(0,3,1)=   1, -3, 0, 1, -1, 3, -2
%_FNCNAME%(0,4,1)= 0.0, -1, 0, 227.2246, 0, 0, -3   
%_FNCNAME%(0,5,1)= 0.0, -2, 0, 1, -3, 3, -1
%_FNCNAME%(0,6,1)= 0.0, -1, 0, 0.13, 0, 0, -2   
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,8,1)= 0.0, -1, 0, 0, 0, 0, 0   
%_FNCNAME%(0,9,1)= 0.0, -2, 0, 1, 0, 0, -1  
%_FNCNAME%(0,10,1)= 0.0, -4, 0, 1, -1, 2, -2
%_FNCNAME%(0,11,1)= 0.0, -1, 0, 0.7426, 0, 0, -4
%_FNCNAME%(0,12,1)= 0.0, -2, 0, 1, -4, 3, -1
%_FNCNAME%(0,13,1)= 0.0, -1, 0, 1, -2, 3, 1
%_FNCNAME%(0,14,1)= 0.0, -1, 7, 1, -1, 0, 0
%_FNCNAME%(0,15,1)= 0.0, -2, 0, 1, -3, 3, -1
%_FNCNAME%(0,16,1)= 0.0, 99, 0, 1, -2, 0, 0
   ! End of equation: 2450*0.7462*227.2246*0.13*exp(-0.7426*{TIME})

/SOLU                       !进入求解器
ANTYPE,4                    !选择瞬态求解                 
TRNOPT,FULL
LUMPM,0
ESEL,S,MAT,,1               !选择混凝土体的所有结点
NSLE,S  
IC,all,TEMP,20           
!混凝土的浇筑温度为20
SFA,4,1,CONV,616.32, %QIWEN%  !四周木模对流系数616.32
SFA,6,1,CONV,616.32, %QIWEN%
SFA,2,1,CONV,2054.4, %QIWEN%  !上面裸露对流系数616.32
BFV,1,HGEN, %HEATS%           !输入热生成率heats
ALLSE                        

!选择所有
TRNOPT,FULL                   !设置计算方法
*DO,time,1,7*24,1            !设置循环计算
TIME,time                      !设置荷载步时间
AUTOTS,0                      !关闭自动时间步长
KBC,0                         !设置阶梯荷载步
TSRES,ERASE                   !Defines an array of keytimes at which the time-stepping strategy changes
OUTRES,ALL,ALL                !输出所有项
SOLVE                         !开始计算
*ENDDO                        !终止循环计算

!结构应力分析
/PREP7                                !前处理
ETCHG,solid65                            !转换热分析单元为结构单元
!输入参数
!定义变化弹模EXX
*DIM,EXX,ARRAY,20
EXX(1,1,1) = 0.87e10   
EXX(2,1,1) = 1.05e10   
EXX(3,1,1) = 1.17e10   
EXX(4,1,1) = 1.26e10   
EXX(5,1,1) = 1.33e10   
EXX(6,1,1) = 1.39e10   
EXX(7,1,1) = 1.44e10   
EXX(8,1,1) = 1.48e10   
EXX(9,1,1) = 1.52e10   
EXX(10,1,1) = 1.56e10  
EXX(11,1,1) = 1.59e10  
EXX(12,1,1) = 1.61e10  
EXX(13,1,1) = 1.64e10  
EXX(14,1,1) = 1.67e10  
EXX(15,1,1) = 1.69e10  
EXX(16,1,1) = 1.71e10  
EXX(17,1,1) = 1.73e10  
EXX(18,1,1) = 1.75e10  
EXX(19,1,1) = 1.77e10  
EXX(20,1,1) = 1.79e10
!定义20种材料每8小时一种材料来反映弹模的改变
*DO,I,1,20,1                                   
    MP,DENS,I,2450        !混凝土 密度为2450
    MP,EX,I,EXX(I)        !混凝土 弹模为EXX(I)
    MP,PRXY,I,0.167       !混凝土 泊松比为0.167
    MP,ALPX,I,0.75E-5     !混凝土 热膨胀系数为0.75e-5
    MP,REFT,I,15          !定义参考温度
*ENDDO
ALLSEL                    !选择选择所有

/SOLU                     !进入求解器
NROPT,FULL                !Specifies the Newton-Raphson options in a full transient analysis
NLGEOM,ON                 !考虑大变形的影响
ASEL,S,,,1,1,1            !选择面1
NSLA,S,1                  !选择面1上的结点
D,all,ALL                  !加基础全约束
ALLSEL                    !选择选择所有
DA,3,SYMM                 !加对称约束
DA,5,SYMM                 !加对称约束
I=1
*DO,time,1,24*7,1
    ESEL, S, MAT, , 1, 1,1                    !选择混凝土材料号为1的元素
    ACEL,0,0,-9.8,                            !GRAVITY LOAD 施加重力荷载加速度
    LDREAD,TEMP,,,time, ,'themal','rth',' '    !TEMP LOAD Reads results from the results file and applies them as loads.
    TIME,time
    AUTOTS,0
    KBC,0                                     !设置阶梯荷载步
    TSRES,ERASE                               !Defines an array of keytimes at which the time-stepping strategy changes
    OUTRES,ALL,ALL,                           !输出所有的结果
    SOLVE                                     !计算
    I=I+1
    *IF,I,LE,20,THEN
        MPCHG,I,ALL
    *ENDIF
*ENDDO
请高手指点下水化热怎么加呀