13-浇筑混凝土.txt

Ansys 钢筋混凝土命令流练习！ Ansys 钢筋混凝土命令流练习！

!大伙看看，有何问题，如何改进？ 
!100*80*10基岩上浇筑混凝土50*40*10，每三天浇筑1米高（第一天浇完，剩下两天窝工）， 
!共用一个月时间,分为十层（layer)浇筑，模型、条件均为假设 
!建模 
/filename,test 
/PREP7 
BLOCK,,50,,10,,40, 
BLOCK,-25,75,-10,0,-20,60 
ET,1,70 
LSEL,S,,,9,12,1 
LESIZE,ALL,,,10 
LSEL,S,,,21,24,1 
LESIZE,ALL,,,20 
LSEL,S,,,2,4,2 
LSEL,A,,,5,7,2 
LESIZE,ALL,,,10 
LSEL,S,,,14,16,2 
LSEL,A,,,17,19,2 
LESIZE,ALL,,,20 
LSEL,S,,,1,3,2 
LSEL,A,,,6,8,2 
LESIZE,ALL,,,10 
LSEL,S,,,13,15,2 
LSEL,A,,,18,20,2 
LESIZE,ALL,,,5 
allsel 
MP,DENS,1,2400.0 !混凝土材料 
MP,KXX ,1,300.89 
MP,C,1,0.94500 
MAT,1 
VMESH,1 
MP,DENS,2,2600.0 !基岩材料 
MP,KXX ,2,300.89 
MP,C,2,0.94500 
MAT,2 
VMESH,2 
nummrg,all 

!热分析 
/SOLU 
NROPT,FULL !Specifies the Newton-Raphson options in full transient analysis 
ESEL,S,MAT,,2 !MAT2 IS 基岩 
NSLE,S 
IC,all,TEMP,16.4, !加基岩初温 
ASEL,S,,,7,9,1 
ASEL,A,,,11,12,1 
NSLA,S,1 
SF,all,HFLUX,0 !加基岩绝热边界条件 
NSEL,S,LOC,Y,-0.1,0.1 
ASEL,S,,,3 
NSLA,U,1 
SF,ALL,CONV,2016,20 !加基岩第三边界条件（和空气接触面） 
!加初始条件 
ALLSEL,ALL 
ESEL,U,MAT,,2 
NSLE,S 
EKILL,ALL !KILL 混凝土 
IC,all,TEMP,35, !混凝土浇筑时温度为35度 
ANTYPE,TRANSIENT,new !分析类型 
!日子的流逝 
*DO,LAYER,1,10,1 
  !激活第N层单元 
  ALLSEL,ALL 
  NSEL,S,LOC,Y,(LAYER-1),LAYER !第一层从Y=0~1米，第二层从Y=1~2米,... 
  ESLN,S ,1 
  EALIVE,ALL 
  *DO,DAY,LAYER*3-2,LAYER*3,1 !每三天浇筑一层 
    ALLSEL,ALL 
    ASEL,S,,,1,2,1 
    ASEL,A,,,5,6,1 
    NSLA,S,1 
    NSEL,R,LOC,Y,0,LAYER 
    NSEL,A,LOC,Y,LAYER-0.1,LAYER+0.1 
    SF,ALL,CONV,2016,20 !加混凝土第三类对流边界条件 
    *DO,LOOP,1,LAYER,1 
      ALLSEL,ALL 
      NSEL,R,LOC,Y,LOOP-1,LOOP 
      ESLN,S,1 
      TDAY=DAY-(LOOP-1)*3 !已经完工的天数 
      HE00=47880*(TDAY/(0.862+TDAY)-(TDAY-1)/(0.862+TDAY-1)) 
      BFE,ALL,HGEN, ,HE00 !加水化热 
    *ENDDO 
    !计算设置 
    ALLSE 
    TRNOPT,FULL 
    LUMPM,0 
    TIME,DAY 
    AUTOTS,0 
    DELTIM,1, , ,1 
    KBC,0 
    TSRES,ERASE 
    OUTRES,ALL,ALL, 
    SOLVE 
  *ENDDO 
  ESEL,S,MAT,,1 
  NSLE,S, 
  SFEDELE,ALL,ALL,CONV !删对流边界条件 
  SFDELE,ALL,CONV 
*ENDDO 
FINISH 
   
!结构分析 
/PREP7 
ET,1,SOLID65 !变热单元为结构单元 
*DIM,EXX,ARRAY,30 !考虑变化的EX 
EXX(1)=0.325E10,0.57778E10,0.78E10,0.94545E10,0.10833E11,0.12E11,0.13E11,0.13867E11,0.14625E11,0.15294E11, 0.15889E11, 0.1642E11 ,0.169E11, 0.17333E11, 0.17727E11, 0.18087E11,0.18417E11,0.1872E11 
EXX(19)=0.190E11, 0.19259E11, 0.195E11, 0.19724E11, 0.19933E11, 0.20129E11,0.20312E11, 0.20485E11, 0.20647E11, 0.208E11, 0.20944E11, 0.21081E11 
*DO,I,1,30,1 !定义材料，混凝土浇筑完每过一天变一种材料 
    MP,DENS,I,2400 
    MP,EX,I,EXX(I) 
    MP,PRXY,I,0.167 
    MP,ALPX,I,0.9E-5 
    MP,REFT,I,15 
*ENDDO 
MP,DENS,31,2600 !定义基岩材料 
MP,EX,31,0.20800E+11 
MP,PRXY,31,0.167 
MP,ALPX,31,0.9E-5 
ESEL,S,MAT,,,2 
MPCHG,31,ALL, !!变热材料为结构材料 
ALLSEL 
/SOLU 
NROPT,FULL !Specifies the Newton-Raphson options in a full transient analysis 
NLGEOM,ON 
NSLA,S,,,7,9,1 
NSLA,A,,,11,12,1 
D,all,ALL !加基岩约束 
!KILL THE CONCRETE ELEMENT 
ESEL,U,MAT,,1 
NSLE,S 
EKILL,ALL !KILL 混凝土 
*DO,DAY,1,30,1 !随时间改变，变换材料 
  *DO,LAYER,1,10 
      *IF,DAY,GT,LAYER*3-3,THEN 
          NSEL,S,LOC,Y,LAYER-1,LAYER 
          ESLN,S,1 
          EALIVE,ALL !ACTIVE NEXT LAYER 
          MPCHG,DAY-(LAYER-1)*3,ALL, !CHANG MAT 
          ALLSEL,ALL 
      *ENDIF 
  *ENDDO 
  ALLSEL,ALL 
  ACEL,0,0,-9.8, !GRAVITY LOAD 
  LDREAD,TEMP,,,DAY, ,'TEST','rth',' ' !TEMP LOAD Reads results from the results file and applies them as loads. 
  CNVTOL,F, ,0.2,2,200, !Sets convergence values for nonlinear analyses. 
  TIME,DAY 
  AUTOTS,0 
  DELTIM,1, , ,1 
  KBC,0 
  TSRES,ERASE !Defines an array of keytimes at which the time-stepping strategy changes 
  OUTRES,ALL,ALL, 
  SOLVE 
*ENDDO 

!1、定义重力集散地是应为正值：ACEL,0,0,9.8 
!2、砼单元边长（1*4*5）相差5倍，对计算精度的影响会不会较大？ 
!3、热单元SOLID70对应的结构单元为SOLID45，改为SOLID65是否合适？ 
!4、原分析得到的结果与实际情况比较如何？ 
!5、如何体现气温的变化？