单曲拱坝

!

为查值计算坝体温度，重新定义一套数组

*DIM,TEMPUP,TABLE,LAYER_NUM !

上游面温度

*DIM,TEMPDOWN,TABLE,LAYER_NUM !

下游面温度

*DIM,T2_ARCH,TABLE,LAYER_NUM !

拱厚

*DIM,Y_RAD,TABLE,LAYER_NUM !

圆心到坝轴线距离

*DIM,RADUP,TABLE,LAYER_NUM !

外半径

TEMPUP(1)=-8.84,-7.30,-5.66,-4.94,-4.99,-5.17,-5.7 3,-6.25

TEMPUP(1,0)=0,6,14,24,29,34,44,53

TEMPDOWN(1)=-8.84,-8.21,-7. 36,-6.32,-5.86,-5.45,-4.72,-6.71

TEMPDOWN(1,0)=0,6,14,24,29,34,44,53

T2_ARCH(1)=5.00,6.68,8.27,1 0.25,11.24,12.23,14.22,16.00

T2_ARCH(1,0)=0,6,14,24,29,34,44,53

Y_RAD(1)=68.00,68.00,68.00, 68.00,68.00,68.00,68.00,68.00

Y_RAD(1,0)=0,6,14,24,29,34,44,53

RADUP(1)=68.00,68.00,68.00, 68.00,68.00,68.00,68.00,68.00

RADUP(1,0)=0,6,14,24,29,34,44,53

CSYS,0

!/PNUM,KP,1

!/PNUM,LINE,1

/PREP7

KNN=0

*DO,II,1,LAYER_NUM !

生成拱坝控制关键点

,

每层

6

个点

,

上下游个三个点

PX1=ARCH_RAD(II,1)*SIN (ARCH_ANGLE(II,1))

PY1=RAD_ CEN(II)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,1))

PX2=ARCH_RAD(II,1)*SIN(ARCH _ANGLE(II,2))

PY2=RAD_CEN(I I)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,2))

K,KNN+1,PX1, PY1,ELEVATION(II)

K,KNN+2,0,RAD_CEN(II)-ARCH_

RAD(II,1),ELEVATION(II)

K,KNN+3,PX2, PY2,ELEVATION(II)

PX3=ARCH_RAD(II ,2)*SIN(ARCH_ANGLE(II,1))

P Y3=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,1) )

PX4=ARCH_RAD(II,2)*SIN(AR CH_ANGLE(II,2))

PY4=RAD_CEN (II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,2))

K,KNN+4,PX3, PY3,ELEVATION(II)

K,KNN+5,0 ,RAD_CEN(II)-ARCH_RAD(II,2),ELEVATION(II)

K,KNN+6,PX4, PY4,ELEVATION(II)

KNN=6*II

*ENDDO

LARC,1,3,2 !

生成拱圈线

*REPEAT,LAYER_NUM*2,3,3,3

L,1,4

*REPEAT,LAYER_NUM,6,6

L,3,6

*REPEAT,LAYER_NUM,6,6

ASKIN,1,3,5,7,9,11,13,15 !

生成上游面

ASKIN,2,4,6,8,10,12,14,16 !

生成下游面

ASKIN,17,18,19,20,21,22,23,24 !

生成侧面

ASKIN,25,26,27,28,29,30,31,32

AL,1,17,2,25

AL,15,24,16,32

VA,6,1,3,2,4,5

拉伸坝体侧面、底面，并进行切割，生成坝基

!

拉伸生成坝基

ASEL,S,AREA,,4 !

选择右侧面

VEXT,ALL,,,4*H_DAM !

拉伸成基岩

ASEL,S,AREA,,3 !

选择左侧面

VEXT,ALL,,,-4*H_DAM !

拉伸成基岩

ALLSEL

ASEL,S,LO C,Z,ELEVATION(LAYER_NUM)-1,ELEVATION(LAYER_NUM)+1 !

选择底面

VEXT,ALL,,,,,-NINT(H_DAM /10)*10

ASEL,S,AREA,,10,20,5

ASEL,A,AREA,,23,27,4

VEXT,ALL,,,,4*H_DAM

ASEL,S,AREA,,8,18,5

ASEL,A,AREA,,25,29,4

VEXT,ALL,,,,-4*H_DAM

ALLSEL,ALL

SAVE

WPROTA,,90

WPOFFS,,,NINT(H_DAM*1.5/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,Y,-200,-H_DAM*1.4

VDEL,ALL,,,1

VSEL,ALL

WPOFFS,,,-NINT(H_DAM*3.5/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,Y,NINT(H_DAM*2/1 0)*10,NINT(H_DAM*2/10)*100

VDEL,ALL,,,1

VSEL,ALL

WPCSYS,,0

WPROTA,,,90

WPOFFS,,,NINT((W_DAM2+H_DAM)/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,X,NINT((W_DAM2+H _DAM)/10)*10,NINT((W_DAM2+H_DAM)/10)*50

VDEL,ALL,,,1

VSEL,ALL

WPOFFS, ,,-NINT((W_DAM2+W_DAM1+2*H_DAM)/10)*10

VSBW,ALL,,DELETE

VSEL,S,LOC,X,-NINT((W_DAM1+H_DAM)/10)*50,-NINT((W _DAM1+H_DAM)/10)*10

VDEL,ALL,,,1

WPCSYS,,0

VSEL,ALL

SAVE

用各控制高程切割整个模型，便于划分单元及计算成果和拱梁法比较

*DO,II,LAYER_NUM,3,-1

WPOFFS,,,ELEVATION(II-1)-ELEVATION(II)

VSBW,ALL,,DELETE

*ENDDO

WPCSYS,,0

ALLSEL

WPROTA,,,90

VSBW,ALL,,DELETE

WPCSYS,,0

NUMMRG,ALL

NUMCMP

,ALL

PX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))

PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))

VSEL,S,LOC,X,PX1,PX2

VSEL,R,LOC,Z,ELEVATION(LAYE R_NUM),ELEVATION(1)

CM,DAM,VOLU

VSEL,INVE

CM,BEDROCK,VOLU

ALLSEL

生成溢流堰

CMSEL,S,DAM

VSEL,R,LOC,Z,ELEVATION(1)-0 .1,ELEVATION(2)+0.1

/VIEW,1,1,1,1

VPLOT

!

生成溢流堰

,

为便于建模和施加 荷载

,

与实际不太相符

,

根据圣维南原理

,

对大坝的总体应

< br>力不会有太大影响

W_WEIR=8 !

溢流堰单孔宽

(

坝轴线处

)

NUM_WEIR=5 !

溢流堰孔数

W_WALL=2.5 !

溢流堰闸 墩宽

(

坝轴线处

)

WPOFFS,,,ELEVATION(1)-ELEVA TION(LAYER_NUM)

WPOFFS,,RAD_CEN(1)

WPROTA,,,90

*AFUN,RAD

PI=2*ACOS(0)

SIT TA1=(W_WEIR*NUM_WEIR+W_WALL*(NUM_WEIR+1))/68*180/P I !

溢流堰总中心角

SITTA2=W_W

EIR/ARCH_RAD(1,1)*180/PI !

溢流

堰单孔中心角

SITTA3=W_WALL/ARCH_RAD(1,1)*180/PI !

溢流堰单个闸墩中心角

WPROTA,,-SITTA1/2

VSBW,ALL,,DELETE

*DO,II,1,NUM_WEIR

WPROTA,,SITTA3

VSBW,ALL,,DELETE

WPROTA,,SITTA2

VSBW,ALL,,DELETE

*ENDDO

WPROTA,,SITTA3

VSBW,ALL,,DELETE

SAVE

CMSEL,S,BEDROCK

VSEL,INVE

VSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)

CSYS,11

VSEL,R,LOC,Y

,-SITTA1/2,SITTA1/2

CM,YLYY,VOLU

*DO,II,1,NUM_WEIR

CMSEL,S,YLYY

VSE L,R,LOC,Y

,-SITTA1/2+II*SITTA3+(II-1)*SI TTA2,-SITTA1/2+II*SITTA3+II*SITTA2

VDEL,ALL,,,1,,,1

*ENDDO

CSYS,0

CMSEL,S,BEDROCK

VSEL,INVE

CM,DAM,VOLU

CMSEL,S,DAM

VSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)

CSYS,11

VSEL,R,LOC,Y

,-SITTA1/2,SITTA1/2,

CM,ZHADUN,VOLU

CSYS,0

WPCSYS,,0

VSEL,ALL

NUMMRG,ALL

NUMCMP

,ALL

SAVE

划分单元

!

定义单元性质和材料

ET,1,SOLID95

MP,EX,1,1E7

MP,NUXY,1,0.25

MP,DENS,1,23 !

坝体容重

MP,ALPX,1,0.7E-5

MP,REFT,1,0

MP,EX,2,0.8E7

MP,NUXY,2,0.21

MP,ALPX,2,0.0

MP,REFT,2,0

!

坝体单元划分

CMSEL,S,DAM

CMSEL,U,ZHADUN

ASLV,S

LSLA,S

*AFUN,DEG

PX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))

PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))

PX3=ARCH_RAD(LAYER_NUM,2)*S IN(ARCH_ANGLE(LAYER_NUM,1))

PX4=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NU M,2))

LSEL,U,LOC,X,PX3,PX1

LSEL,U,LOC,X,PX4,PX2

LSEL,U,LOC,X,0

LSEL,U,LOC,Z,ELEVATION(1),ELEVATION(2)

LESIZE,ALL, , ,20,

CM,LTEMP1,LINE