iCAx开思网

标题: 【分享】某混凝土拱坝工程施工期及运行期温度场仿真分析数据流 [打印本页]

作者: kkkyyy 时间: 2003-11-22 06:32
标题: 【分享】某混凝土拱坝工程施工期及运行期温度场仿真分析数据流
--------------------------------------------------------------------------------


      
https://okok.org/cgi-bin/ut/topic_show.cgi?id=42320&h=1#192171

  帖一个某混凝土拱坝工程施工期及运行期温度场仿真分析数据流，供大家参考。
      附件中有两个文件：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>

欢迎光临 iCAx开思网 (https://www.icax.org/)