68．某混凝土拱坝工程施工期及运行期温度场仿真分析
供大家参考。
附件中有两个文件：CA1*为计算数据流，DAQI.FUC为大气年变化的周期函数.
计算简介:
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,ALL,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
===============================
<EquationData>
<Comments>
</Comments>
<RegimeParameter>
<Regime Var>
</RegimeParameter>
<AngularUnits type="0">
! 0 = radian, 1 = degrees
</AngularUnits>
<FunctionType type="SINGLE">
</FunctionType>
<Regime tab="0">
<TabTitle>
Function
</TabTitle>
<Equation>
17.3+11.35*cos(0.01745*({TIME}-215))
</Equation>
<MaterialId active="0">
</MaterialId>
</Regime>
</EquationData>
<!--
!
! ANSYS input
!
:FUNC
*DIM,%_FNCNAME%,TABLE,6,10,1
!
! Begin of equation: 17.3+11.35*cos(0.01745*({TIME}-215))
%_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, 215, 0, 0, 1
%_FNCNAME%(0,2,1)= 0.0, -2, 0, 1, 1, 2, -1
%_FNCNAME%(0,3,1)= 0, -1, 0, 0.01745, 0, 0, -2
%_FNCNAME%(0,4,1)= 0.0, -3, 0, 1, -1, 3, -2
%_FNCNAME%(0,5,1)= 0.0, -1, 10, 1, -3, 0, 0
%_FNCNAME%(0,6,1)= 0.0, -2, 0, 11.35, 0, 0, -1
%_FNCNAME%(0,7,1)= 0.0, -3, 0, 1, -2, 3, -1
%_FNCNAME%(0,8,1)= 0.0, -1, 0, 17.3, 0, 0, -3
%_FNCNAME%(0,9,1)= 0.0, -2, 0, 1, -1, 1, -3
%_FNCNAME%(0,10,1)= 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 17.3+11.35*cos(0.01745*({TIME}-215))
!-->
