用Ansys－APDL建立拱坝静力线性计算模型例题详解

某砌石拱坝位于U型河谷中，坝高55.5m，为单曲等厚拱坝，顶宽5m，底宽16m，坝顶弧长115.65m，弧高比2.1。首先输入基本参数，自底向上生成坝体。

!坐标系原点位于拱坝轴线与拱坝参考面的交点在坝基的投影，X轴指向右岸，Y轴指向下游，Z轴铅直向上
!单位:长度-M,力-KN,加速度-M/S,密度-KN/M3
!
!
FINI
/CLEAR,START
/FILNAME,ARCHDAM,1
/PLOPTS,DATE,0
/TRIAD,LBOT
/VIEW,1,1,1,1
/VUP,1,Z
*AFUN,DEG
Z_UP=53 !上游水位
Z_DOWN=0 !下游水位
Z_SAND=38 !泥沙淤积水位
DENS_SAND=0.6 !淤积泥沙浮容重
ANG_FRI=16 !淤积泥沙内摩擦角
LAYER_NUM=8 ! 定义拱层数
*DIM,ELEVATION,ARRAY,LAYER_NUM !定义拱圈高程
*DIM,T_ARCH,ARRAY,LAYER_NUM !定义拱厚
*DIM,ARCH_RAD,ARRAY,LAYER_NUM,2 !定义圆弧半径，上下游
*DIM,RAD_CEN,ARRAY,LAYER_NUM !定义圆弧中心Y值
*DIM,ARCH_ANGLE,ARRAY,LAYER_NUM,2 !定义圆弧左右中心角
ELEVATION(1)=55.5,47,39,29,24,19,9,0 !拱圈高程
H_DAM=ELEVATION(1)-ELEVATION(LAYER_NUM) !坝高
T_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00 !拱厚
ARCH_RAD(1,1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00 !外半径
*VOPER,ARCH_RAD(1,2),ARCH_RAD(1,1),SUB,T_ARCH(1) !计算内半径
RAD_CEN(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00 !圆心位置,Y坐标
ARCH_ANGLE(1,1)=-47.72,-46.36,-45.06,-43.36,-41.88,-40.38,-37.32,-34.48 !左中心角
ARCH_ANGLE(1,2)=49.73,47.13,44.68,41.57,40.44,39.29,36.92,35.72 !右中心角
W_DAM1=-ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1)) !坝顶左弦长
W_DAM2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2)) !坝顶右弦长
LOCAL,11,1,0,RAD_CEN(1),ELEVATION(1),-90 !局部柱面坐标,原点位于顶拱圆心,X轴为拱坝中心线,指向上游,Y轴为中心角

!为查值计算坝体温度，重新定义一套数组
*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.73,-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,10.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(II)-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))
PY3=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,1))
PX4=ARCH_RAD(II,2)*SIN(ARCH_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