|
计算简介:
1.施工期共分37层,每层1.5米
2.施工期及坝体建成后一个月的时间步长按天考虑,随后32个月时间步长按月计.
3.计算中多年平均气温作为岩体初始温度场,各层砼浇筑温度作为其激活时的初始温度;
4.岩体边界按绝热边界条件(第一类边界条件);大气与坝面按对流边界条件（第一类边界条件）施加;按分段线性插值函数计水化热。
====================
!!!!**********************************************************************
!!!!**********************************************************************
!!!!***************某拱坝温度施工期与运行期仿真分析数据流*****************
!!!!************************* ZH.EI ************************************
!!!!**********************************************************************
!!!**************注意事项**********************************
!!1.对流边界条件应先选节点,再将所加边界条件加在相应的单元上.
!!2.加对流条件时结点及相应单元必须都选.
!2003.03.01
!按施工单位提供的浇筑方案.浇筑温度等于月平均温度.温控措施,上游水温取河水温度，保温.
FINISH
/SOLU
LSCLEAR,ALL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
NROPT,FULL
!ESTIF,1e-0024,
!加岩石的初温
ALLSEL,ALL
LSCLEAR,ALL
VSEL,S,MAT,,3
VSEL,a,MAT,,2
ESLV,S
NSLE,S
IC,all,TEMP,17.3,
!加岩石的绝热边界条件
allsel,all
NSEL,S,LOC,Z,-39,-41
NSEL,A,LOC,X,-458,-457
NSEL,A,LOC,X,317,318
NSEL,A,LOC,y,49,51
NSEL,A,LOC,y,449,451
SF,all,HFLUX,0
!~~~~~~~~~~~~~~~~~~~~~~~~~加初始条件~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!加坝体第1层初始条件
ALLSEL,ALL
ESEL,S,MAT,,1
NSLE,S
IC,all,TEMP,10.0,
EKILL,ALL !KILL WHOLE DAM
CMSEL,S,LAYER(1,1)
EALIVE,ALL
NSLE,S
!~~~~~~~~~~~~~~~~~~~~~~~~~加边界条件并计算~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!!#######################################LAYER 1,1-3DAY##################################
*DO,T,1,3,1
!加第三类边界条件
ESEL,A,MAT,,2
ESEL,A,MAT,,3
CMSEL,S,NBIAOMIAN
CMSEL,A,JCNODE
*DO,I,1,3,1
CMSEL,A,NLAYER(1,I)
*ENDDO
!SF,ALL,CONV,200, %DAQI%
SF,ALL,CONV,2016, %DAQI%
!加水化热
CMSEL,S,LAYER(1,1)
NSLE,S
T0=0
HE01=25538*((T-T0)/(3.19+T-T0)-(T-T0-1)/(3.19+T-T0-1))
BFE,ALL,HGEN, ,HE01
ALLSEL,ALL
!计算设置
TRNOPT,FULL
LUMPM,0
TIME,T
AUTOTS,0
DELTIM,1, , ,1
KBC,0
TSRES,ERASE
OUTRES,ALL,ALL,
SOLVE
*ENDDO
!!#######################################LAYER2-37##################################
!施工弟M层计算
*DO,M,2,37,1
!删对流边界条件
/SOLU
SFEDELE,ALL,ALL,CONV
SFDELE,ALL,CONV
!激活第M层单元
CMSEL,S,LAYER(M,1)
NSLE,S
EALIVE,ALL
*DO,T,3*(M-1)+1,3*M,1
!!选出坝体活单元和岩石及纵向围堰单元,使对流边界条件可加在岩石及纵向围堰与大气相交面,
!!并使与坝体活单元相连的岩石接触节点不加上对流边界条件.
CMSEL,S,LAYER(1,1)
*DO,I,1,M,1
CMSEL,A,LAYER(I,1)
*ENDDO
ESEL,A,MAT,,2
ESEL,A,MAT,,3
!!加第三类边界条件
!选岩石与大气接触节点
CMSEL,S,NBIAOMIAN
CMSEL,A,JCNODE
!选坝体与大气接触节点
*DO,I,1,M,1
*DO,J,1,2,1
CMSEL,A,NLAYER(I,J)
*ENDDO
*ENDDO
CMSEL,A,NLAYER(I,3)
!SF,ALL,CONV,200, %DAQI%
SF,A
LL,CONV,2016, %DAQI%
!!加水化热
*DO,I,1,M,1
CMSEL,S,LAYER(I,1)
NSLE,S
T0=3*(I-1)
HE01=25538*((T-T0)/(3.19+T-T0)-(T-T0-1)/(3.19+T-T0-1))
BFE,ALL,HGEN, ,HE01
*ENDDO
ALLSEL,ALL
!计算设置
TRNOPT,FULL
LUMPM,0
TIME,T
AUTOTS,0
DELTIM,1, , ,1
KBC,0
TSRES,ERASE
OUTRES,ALL,ALL,
SOLVE
*ENDDO
*ENDDO
!!##############################LAYER1-37整体施工完毕计算一个月##########################
!!主要是为了使坝体水化热散发完毕
*DO,T,112,141,1
!!岩石及纵向围堰单元,使对流边界条件可加在岩石及纵向围堰与大气相交面.
ESEL,A,MAT,,2
ESEL,A,MAT,,3
!!加第三类边界条件
!选岩石与大气接触节点
CMSEL,S,NBIAOMIAN
!选坝体与大气接触节点
*DO,I,1,37,1
*DO,J,1,2,1
CMSEL,A,NLAYER(I,J)
*ENDDO
*ENDDO
CMSEL,A,NLAYER(37,3)
!SF,ALL,CONV,200, %DAQI%
SF,ALL,CONV,2016, %DAQI%
!!加水化热
*DO,I,1,37,1
CMSEL,S,LAYER(I,1)
NSLE,S
T0=3*(I-1)
HE01=25538*((T-T0)/(3.19+T-T0)-(T-T0-1)/(3.19+T-T0-1))
BFE,ALL,HGEN, ,HE01
*ENDDO
ALLSEL,ALL
!计算设置
TRNOPT,FULL
LUMPM,0
TIME,T
AUTOTS,0
DELTIM,1, , ,1
KBC,0
TSRES,ERASE
OUTRES,ALL,ALL,
SOLVE
*ENDDO
!!##############################LAYER1-37整体施工完毕按月计算##########################
!!主要是为了节省计算空间
*DO,M,1,32,1
T=142+(M-1)*30
!!岩石及纵向围堰单元,使对流边界条件可加在岩石及纵向围堰与大气相交面.
ESEL,A,MAT,,2
ESEL,A,MAT,,3
!!加第三类边界条件
!选岩石与大气接触节点
CMSEL,S,NBIAOMIAN
!选坝体与大气接触节点
*DO,I,1,37,1
*DO,J,1,2,1
CMSEL,A,NLAYER(I,J)
*ENDDO
*ENDDO
CMSEL,A,NLAYER(37,3)
!SF,ALL,CONV,200, %DAQI%
SF,ALL,CONV,2016, %DAQI%
ALLSEL,ALL
!计算设置
TRNOPT,FULL
LUMPM,0
TIME,T
AUTOTS,-1
NSUBST,1, , ,1
KBC,0
TSRES,ERASE
OUTRES,ALL,ALL,
SOLVE
*ENDDO
===============================
! 0 = radian, 1 = degrees
Function
17.3+11.35*cos(0.01745*({TIME}-215))
|