在0号块水化热分析中定义大气温度和水化热生成怎么样定义
*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-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
请高手指点下水化热怎么加呀
页:
[1]