[求助]ANSYS拱坝建模命令流

johhan

请问哪位朋友有ANSYS拱坝建模命令流,谢谢!!!!!!!!!!!

AaronSpark

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

1. !坐标系原点位于拱坝轴线与拱坝参考面的交点在坝基的投影，X轴指向右岸，Y轴指向下游，Z轴铅直向上
   
2. !单位:长度-M,力-KN,加速度-M/S,密度-KN/M3
   
3. !
   
4. !
   
5. FINI
   
6. /CLEAR,START
   
7. /FILNAME,ARCHDAM,1
   
8. /PLOPTS,DATE,0
   
9. /TRIAD,LBOT
   
10. /VIEW,1,1,1,1
    
11. /VUP,1,Z
    
12. *AFUN,DEG
    
13. Z_UP=53 !上游水位
    
14. Z_DOWN=0 !下游水位
    
15. Z_SAND=38 !泥沙淤积水位
    
16. DENS_SAND=0.6 !淤积泥沙浮容重
    
17. ANG_FRI=16 !淤积泥沙内摩擦角
    
18. LAYER_NUM=8 ! 定义拱层数
    
19. *DIM,ELEVATION,ARRAY,LAYER_NUM !定义拱圈高程
    
20. *DIM,T_ARCH,ARRAY,LAYER_NUM !定义拱厚
    
21. *DIM,ARCH_RAD,ARRAY,LAYER_NUM,2 !定义圆弧半径，上下游
    
22. *DIM,RAD_CEN,ARRAY,LAYER_NUM !定义圆弧中心Y值
    
23. *DIM,ARCH_ANGLE,ARRAY,LAYER_NUM,2 !定义圆弧左右中心角
    
24. ELEVATION(1)=55.5,47,39,29,24,19,9,0 !拱圈高程
    
25. H_DAM=ELEVATION(1)-ELEVATION(LAYER_NUM) !坝高
    
26. T_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00 !拱厚
    
27. ARCH_RAD(1,1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00 !外半径
    
28. *VOPER,ARCH_RAD(1,2),ARCH_RAD(1,1),SUB,T_ARCH(1) !计算内半径
    
29. RAD_CEN(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00 !圆心位置,Y坐标
    
30. ARCH_ANGLE(1,1)=-47.72,-46.36,-45.06,-43.36,-41.88,-40.38,-37.32,-34.48 !左中心角
    
31. ARCH_ANGLE(1,2)=49.73,47.13,44.68,41.57,40.44,39.29,36.92,35.72 !右中心角
    
32. W_DAM1=-ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1)) !坝顶左弦长
    
33. W_DAM2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2)) !坝顶右弦长
    
34. LOCAL,11,1,0,RAD_CEN(1),ELEVATION(1),-90 !局部柱面坐标,原点位于顶拱圆心,X轴为拱坝中心线,指向上游,Y轴为中心角
    
35. 
    
36. !为查值计算坝体温度，重新定义一套数组
    
37. *DIM,TEMPUP,TABLE,LAYER_NUM !上游面温度
    
38. *DIM,TEMPDOWN,TABLE,LAYER_NUM !下游面温度
    
39. *DIM,T2_ARCH,TABLE,LAYER_NUM !拱厚
    
40. *DIM,Y_RAD,TABLE,LAYER_NUM !圆心到坝轴线距离
    
41. *DIM,RADUP,TABLE,LAYER_NUM !外半径
    
42. TEMPUP(1)=-8.84,-7.30,-5.66,-4.94,-4.99,-5.17,-5.73,-6.25
    
43. TEMPUP(1,0)=0,6,14,24,29,34,44,53
    
44. TEMPDOWN(1)=-8.84,-8.21,-7.36,-6.32,-5.86,-5.45,-4.72,-6.71
    
45. TEMPDOWN(1,0)=0,6,14,24,29,34,44,53
    
46. T2_ARCH(1)=5.00,6.68,8.27,10.25,11.24,12.23,14.22,16.00
    
47. T2_ARCH(1,0)=0,6,14,24,29,34,44,53
    
48. Y_RAD(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00
    
49. Y_RAD(1,0)=0,6,14,24,29,34,44,53
    
50. RADUP(1)=68.00,68.00,68.00,68.00,68.00,68.00,68.00,68.00
    
51. RADUP(1,0)=0,6,14,24,29,34,44,53
    
52. 
    
53. CSYS,0
    
54. !/PNUM,KP,1
    
55. !/PNUM,LINE,1
    
56. /PREP7
    
57. KNN=0
    
58. *DO,II,1,LAYER_NUM !生成拱坝控制关键点,每层6个点,上下游个三个点
    
59. PX1=ARCH_RAD(II,1)*SIN(ARCH_ANGLE(II,1))
    
60. PY1=RAD_CEN(II)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,1))
    
61. PX2=ARCH_RAD(II,1)*SIN(ARCH_ANGLE(II,2))
    
62. PY2=RAD_CEN(II)-ARCH_RAD(II,1)*COS(ARCH_ANGLE(II,2))
    
63. K,KNN+1,PX1, PY1,ELEVATION(II)
    
64. K,KNN+2,0,RAD_CEN(II)-ARCH_RAD(II,1),ELEVATION(II)
    
65. K,KNN+3,PX2, PY2,ELEVATION(II)
    
66. 
    
67. PX3=ARCH_RAD(II,2)*SIN(ARCH_ANGLE(II,1))
    
68. PY3=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,1))
    
69. PX4=ARCH_RAD(II,2)*SIN(ARCH_ANGLE(II,2))
    
70. PY4=RAD_CEN(II)-ARCH_RAD(II,2)*COS(ARCH_ANGLE(II,2))
    
71. K,KNN+4,PX3, PY3,ELEVATION(II)
    
72. K,KNN+5,0,RAD_CEN(II)-ARCH_RAD(II,2),ELEVATION(II)
    
73. K,KNN+6,PX4, PY4,ELEVATION(II)
    
74. KNN=6*II
    
75. *ENDDO
    
76. LARC,1,3,2 !生成拱圈线
    
77. *REPEAT,LAYER_NUM*2,3,3,3
    
78. L,1,4
    
79. *REPEAT,LAYER_NUM,6,6
    
80. L,3,6
    
81. *REPEAT,LAYER_NUM,6,6
    
82. 
    
83. ASKIN,1,3,5,7,9,11,13,15 !生成上游面
    
84. ASKIN,2,4,6,8,10,12,14,16 !生成下游面
    
85. ASKIN,17,18,19,20,21,22,23,24 !生成侧面
    
86. ASKIN,25,26,27,28,29,30,31,32
    
87. AL,1,17,2,25
    
88. AL,15,24,16,32
    
89. VA,6,1,3,2,4,5

复制代码

AaronSpark

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

1. 
   
2. !拉伸生成坝基
   
3. ASEL,S,AREA,,4 !选择右侧面
   
4. VEXT,ALL,,,4*H_DAM !拉伸成基岩
   
5. ASEL,S,AREA,,3 !选择左侧面
   
6. VEXT,ALL,,,-4*H_DAM !拉伸成基岩
   
7. ALLSEL
   
8. ASEL,S,LOC,Z,ELEVATION(LAYER_NUM)-1,ELEVATION(LAYER_NUM)+1 !选择底面
   
9. VEXT,ALL,,,,,-NINT(H_DAM /10)*10
   
10. ASEL,S,AREA,,10,20,5
    
11. ASEL,A,AREA,,23,27,4
    
12. VEXT,ALL,,,,4*H_DAM
    
13. ASEL,S,AREA,,8,18,5
    
14. ASEL,A,AREA,,25,29,4
    
15. VEXT,ALL,,,,-4*H_DAM
    
16. ALLSEL,ALL
    
17. SAVE
    
18. 
    
19. WPROTA,,90
    
20. WPOFFS,,,NINT(H_DAM*1.5/10)*10
    
21. VSBW,ALL,,DELETE
    
22. VSEL,S,LOC,Y,-200,-H_DAM*1.4
    
23. VDEL,ALL,,,1
    
24. VSEL,ALL
    
25. WPOFFS,,,-NINT(H_DAM*3.5/10)*10
    
26. VSBW,ALL,,DELETE
    
27. VSEL,S,LOC,Y,NINT(H_DAM*2/10)*10,NINT(H_DAM*2/10)*100
    
28. VDEL,ALL,,,1
    
29. VSEL,ALL
    
30. 
    
31. WPCSYS,,0
    
32. WPROTA,,,90
    
33. WPOFFS,,,NINT((W_DAM2+H_DAM)/10)*10
    
34. VSBW,ALL,,DELETE
    
35. VSEL,S,LOC,X,NINT((W_DAM2+H_DAM)/10)*10,NINT((W_DAM2+H_DAM)/10)*50
    
36. VDEL,ALL,,,1
    
37. VSEL,ALL
    
38. WPOFFS,,,-NINT((W_DAM2+W_DAM1+2*H_DAM)/10)*10
    
39. VSBW,ALL,,DELETE
    
40. VSEL,S,LOC,X,-NINT((W_DAM1+H_DAM)/10)*50,-NINT((W_DAM1+H_DAM)/10)*10
    
41. VDEL,ALL,,,1
    
42. WPCSYS,,0
    
43. VSEL,ALL
    
44. SAVE

复制代码

AaronSpark

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

1. *DO,II,LAYER_NUM,3,-1
   
2. WPOFFS,,,ELEVATION(II-1)-ELEVATION(II)
   
3. VSBW,ALL,,DELETE
   
4. *ENDDO
   
5. WPCSYS,,0
   
6. ALLSEL
   
7. WPROTA,,,90
   
8. VSBW,ALL,,DELETE
   
9. WPCSYS,,0
   
10. NUMMRG,ALL
    
11. NUMCMP,ALL
    
12. PX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))
    
13. PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))
    
14. VSEL,S,LOC,X,PX1,PX2
    
15. VSEL,R,LOC,Z,ELEVATION(LAYER_NUM),ELEVATION(1)
    
16. CM,DAM,VOLU
    
17. VSEL,INVE
    
18. CM,BEDROCK,VOLU
    
19. ALLSEL

复制代码

AaronSpark

生成溢流堰

1. CMSEL,S,DAM
   
2. VSEL,R,LOC,Z,ELEVATION(1)-0.1,ELEVATION(2)+0.1
   
3. /VIEW,1,1,1,1
   
4. VPLOT
   
5. 
   
6. !生成溢流堰,为便于建模和施加荷载,与实际不太相符,根据圣维南原理,对大坝的总体应力不会有太大影响
   
7. W_WEIR=8 !溢流堰单孔宽(坝轴线处)
   
8. NUM_WEIR=5 !溢流堰孔数
   
9. W_WALL=2.5 !溢流堰闸墩宽(坝轴线处)
   
10. WPOFFS,,,ELEVATION(1)-ELEVATION(LAYER_NUM)
    
11. WPOFFS,,RAD_CEN(1)
    
12. WPROTA,,,90
    
13. *AFUN,RAD
    
14. PI=2*ACOS(0)
    
15. SITTA1=(W_WEIR*NUM_WEIR+W_WALL*(NUM_WEIR+1))/68*180/PI !溢流堰总中心角
    
16. SITTA2=W_WEIR/ARCH_RAD(1,1)*180/PI !溢流堰单孔中心角
    
17. SITTA3=W_WALL/ARCH_RAD(1,1)*180/PI !溢流堰单个闸墩中心角
    
18. WPROTA,,-SITTA1/2
    
19. VSBW,ALL,,DELETE
    
20. *DO,II,1,NUM_WEIR
    
21. WPROTA,,SITTA3
    
22. VSBW,ALL,,DELETE
    
23. WPROTA,,SITTA2
    
24. VSBW,ALL,,DELETE
    
25. *ENDDO
    
26. WPROTA,,SITTA3
    
27. VSBW,ALL,,DELETE
    
28. SAVE
    
29. 
    
30. CMSEL,S,BEDROCK
    
31. VSEL,INVE
    
32. VSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)
    
33. CSYS,11
    
34. VSEL,R,LOC,Y,-SITTA1/2,SITTA1/2
    
35. CM,YLYY,VOLU
    
36. *DO,II,1,NUM_WEIR
    
37. CMSEL,S,YLYY
    
38. VSEL,R,LOC,Y,-SITTA1/2+II*SITTA3+(II-1)*SITTA2,-SITTA1/2+II*SITTA3+II*SITTA2
    
39. VDEL,ALL,,,1,,,1
    
40. *ENDDO
    
41. CSYS,0
    
42. 
    
43. CMSEL,S,BEDROCK
    
44. VSEL,INVE
    
45. CM,DAM,VOLU
    
46. CMSEL,S,DAM
    
47. VSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)
    
48. CSYS,11
    
49. VSEL,R,LOC,Y,-SITTA1/2,SITTA1/2,
    
50. CM,ZHADUN,VOLU
    
51. CSYS,0
    
52. WPCSYS,,0
    
53. VSEL,ALL
    
54. NUMMRG,ALL
    
55. NUMCMP,ALL
    
56. SAVE

复制代码

AaronSpark

划分单元

1. !定义单元性质和材料
   
2. ET,1,SOLID95
   
3. MP,EX,1,1E7
   
4. MP,NUXY,1,0.25
   
5. MP,DENS,1,23 !坝体容重
   
6. MP,ALPX,1,0.7E-5
   
7. MP,REFT,1,0
   
8. MP,EX,2,0.8E7
   
9. MP,NUXY,2,0.21
   
10. MP,ALPX,2,0.0
    
11. MP,REFT,2,0
    
12. 
    
13. !坝体单元划分
    
14. CMSEL,S,DAM
    
15. CMSEL,U,ZHADUN
    
16. ASLV,S
    
17. LSLA,S
    
18. *AFUN,DEG
    
19. PX1=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,1))
    
20. PX2=ARCH_RAD(1,1)*SIN(ARCH_ANGLE(1,2))
    
21. PX3=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NUM,1))
    
22. PX4=ARCH_RAD(LAYER_NUM,2)*SIN(ARCH_ANGLE(LAYER_NUM,2))
    
23. LSEL,U,LOC,X,PX3,PX1
    
24. LSEL,U,LOC,X,PX4,PX2
    
25. LSEL,U,LOC,X,0
    
26. LSEL,U,LOC,Z,ELEVATION(1),ELEVATION(2)
    
27. LESIZE,ALL, , ,20,
    
28. CM,LTEMP1,LINE
    
29. 
    
30. CSYS,0
    
31. CMSEL,S,DAM
    
32. ASLV,S
    
33. LSLA,S
    
34. CMSEL,U,LTEMP1
    
35. CM,LTEMP2,LINE
    
36. LSEL,U,LENGTH,,T_ARCH(LAYER_NUM)+0.1,1000
    
37. LSEL,R,TAN1,Z
    
38. CSYS,11
    
39. LSEL,U,LOC,X,ARCH_RAD(1,1)
    
40. LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(1)
    
41. LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(2)
    
42. CSYS,0
    
43. CM,LTEMP3,LINE
    
44. LESIZE,ALL, , ,5,
    
45. 
    
46. ALLSEL
    
47. CMSEL,S,LTEMP2
    
48. CMSEL,U,LTEMP3
    
49. LSEL,U,TAN1,Z
    
50. CM,LTEMP4,LINE
    
51. LESIZE,ALL, , ,3,
    
52. 
    
53. CMSEL,S,DAM
    
54. CMSEL,U,ZHADUN
    
55. ASLV,S
    
56. LSLA,S
    
57. LSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)
    
58. CSYS,11
    
59. LSEL,U,LOC,Y,-SITTA1/2-1,SITTA1/2+1
    
60. LSEL,U,LENGTH,,0,T_ARCH(1)+0.1
    
61. CSYS,0
    
62. CMSEL,U,LTEMP3
    
63. CMSEL,U,LTEMP4
    
64. LESIZE,ALL, , ,10,
    
65. CM,LTEMP5,LINE
    
66. 
    
67. CMSEL,S,DAM
    
68. CMSEL,U,ZHADUN
    
69. ASLV,S
    
70. LSLA,S
    
71. LSEL,R,LOC,Z,ELEVATION(1),ELEVATION(2)
    
72. CMSEL,U,LTEMP3
    
73. CMSEL,U,LTEMP4
    
74. CMSEL,U,LTEMP5
    
75. LSEL,U,LENGTH,,0,3
    
76. CM,LTEMP6,LINE
    
77. LESIZE,ALL, , ,4,
    
78. CSYS,0
    
79. 
    
80. ALLSEL
    
81. MSHAPE,0
    
82. CMSEL,S,DAM
    
83. VSWEEP,ALL
    
84. SAVE
    
85. 
    
86. !坝基单元划分
    
87. CMSEL,S,BEDROCK
    
88. ASLV,S
    
89. LSLA,S
    
90. LSEL,U,TAN1,X
    
91. LSEL,U,LOC,X,PX2,PX1
    
92. CM,LTEMP7,LINE
    
93. LESIZE,ALL, , ,5,, , , ,1
    
94. 
    
95. CMSEL,S,BEDROCK
    
96. ASLV,S
    
97. LSLA,S
    
98. LSEL,U,TAN1,Y
    
99. PY4=RAD_CEN(LAYER_NUM)-ARCH_RAD(LAYER_NUM,2)*COS(ARCH_ANGLE(LAYER_NUM,2))+20
    
100. LSEL,U,LOC,Y,-10,PY4
     
101. CM,LTEMP8,LINE
     
102. LESIZE,ALL, , ,6,4, , , ,1
     
103. 
     
104. CMSEL,S,BEDROCK
     
105. ASLV,S
     
106. LSLA,S
     
107. LSEL,U,TAN1,Z
     
108. LSEL,R,LOC,Z,ELEVATION(LAYER_NUM),ELEVATION(LAYER_NUM)-200
     
109. LESIZE,ALL, , ,5,5, , , ,1
     
110. 
     
111. CMSEL,S,BEDROCK
     
112. ASLV,S
     
113. LSLA,S
     
114. LSEL,U,TAN1,Z
     
115. LSEL,R,LOC,Z,ELEVATION(1),ELEVATION(LAYER_NUM)
     
116. LESIZE,ALL, , ,3,
     
117. 
     
118. CMSEL,S,BEDROCK
     
119. ASLV,S
     
120. LSLA,S
     
121. LSEL,R,LOC,Z,-50,-1000
     
122. LSEL,R,LOC,Y,-10,PY4
     
123. CSYS,11
     
124. LSEL,U,LOC,X,ARCH_RAD(1,1)
     
125. LSEL,U,LOC,X,ARCH_RAD(1,1)-T_ARCH(LAYER_NUM)
     
126. CSYS,0
     
127. LSEL,R,LENGTH,,T_ARCH(LAYER_NUM)-5,T_ARCH(LAYER_NUM)+3,
     
128. LESIZE,ALL, , ,5,
     
129. 
     
130. CSYS,0
     
131. CMSEL,S,BEDROCK
     
132. ASLV,S
     
133. LSLA,S
     
134. LSEL,R,LOC,Z,ELEVATION(LAYER_NUM)-50,ELEVATION(LAYER_NUM)-1000
     
135. LSEL,R,LOC,X,PX2,PX1
     
136. LESIZE,ALL, , ,20,
     
137. 
     
138. CMSEL,S,BEDROCK
     
139. MAT,2
     
140. VMESH,ALL
     
141. !单元划分完毕
     
142. 
     
143. ALLSEL
     
144. NUMMRG,ALL
     
145. NUMCMP,ALL
     
146. SAVE

复制代码

AaronSpark

定义边界条件，顺河向约束Y向位移，横河向约束X向位移，底面约束竖向Z位移

1. FINI
   
2. /SOLU
   
3. !边界位移条件
   
4. *GET,NXMAX,NODE,,MXLOC,X
   
5. *GET,NXMIN,NODE,,MNLOC,X
   
6. *GET,NYMAX,NODE,,MXLOC,Y
   
7. *GET,NYMIN,NODE,,MNLOC,Y
   
8. *GET,NZMIN,NODE,,MNLOC,Z
   
9. NSEL,S,LOC,X,NXMAX
   
10. NSEL,A,LOC,X,NXMIN
    
11. D,ALL,UX,0 !约束上下游边界X向位移
    
12. NSEL,S,LOC,Y,NYMAX
    
13. NSEL,A,LOC,Y,NYMIN
    
14. D,ALL,UY,0 !约束左右边界Y向位移
    
15. NSEL,S,LOC,Z,NZMIN
    
16. D,ALL,ALL,0 !约束底面边界Z向位移
    
17. NSEL,ALL

复制代码

AaronSpark

施加水荷载，淤沙压力，坝体自重

1. !大坝荷载
   
2. ALLSEL
   
3. CMSEL,S,DAM
   
4. CMSEL,U,ZHADUN
   
5. ASLV,S
   
6. CSYS,11
   
7. ASEL,R,LOC,X,ARCH_RAD(1,1)-0.1,ARCH_RAD(1,1)+0.1
   
8. NSLA,S,1
   
9. CM,N_DAMUP,NODE !定义大坝上游面节点
   
10. CSYS,0
    
11. 
    
12. ALLSEL
    
13. CMSEL,S,DAM
    
14. ASLV,S
    
15. ASEL,R,EXT
    
16. LSLA,S
    
17. LSEL,R,LOC,X,0
    
18. LSEL,U,LOC,Y,-0.1,0.1
    
19. LSEL,U,LOC,Z,ELEVATION(2)
    
20. LSEL,U,LOC,Z,ELEVATION(LAYER_NUM)
    
21. ASLL,S
    
22. ASEL,U,LOC,X,0
    
23. LSLA,S
    
24. NSLA,S,1
    
25. CM,N_DAMDOWN,NODE !定义大坝下游面节点
    
26. 
    
27. ALLSEL
    
28. CMSEL,S,YLYY
    
29. ASLV,S
    
30. CSYS,11
    
31. ASEL,R,LOC,X,ARCH_RAD(1,1)-0.1,ARCH_RAD(1,1)+0.1
    
32. ASEL,U,,,291,326,326-291
    
33. NSLA,S,1
    
34. CM,N_ZHADUN1,NODE !定义溢流坝中墩节点
    
35. ASEL,S,,,291,326,326-291
    
36. NSLA,S,1
    
37. CM,N_ZHADUN2,NODE !定义溢流坝边墩节点
    
38. 
    
39. CSYS,0
    
40. *IF,Z_UP,NE,ELEVATION(LAYER_NUM),THEN
    
41. CMSEL,S,N_DAMUP
    
42. NSEL,R,LOC,Z,0,Z_UP
    
43. /PSF,PRES,NORM,2,0.1
    
44. SFGRAD, PRES, 0, Z, 0,-9.81
    
45. SF,ALL,PRESS,Z_UP*9.81 !施加大坝上游面水压力
    
46. *ENDIF
    
47. 
    
48. *IF,Z_UP,GT,ELEVATION(2),THEN
    
49. CMSEL,S,N_ZHADUN1
    
50. NSEL,R,LOC,Z,ELEVATION(2),Z_UP
    
51. SFGRAD, PRES, 0, Z, 0,-(W_WEIR+W_WALL)/W_WALL*9.81
    
52. SF,ALL,PRESS,Z_UP*(W_WEIR+W_WALL)/W_WALL*9.81 !施加中墩水压力,包含闸门传递的水压力
    
53. CMSEL,S,N_ZHADUN2
    
54. NSEL,R,LOC,Z,ELEVATION(2),Z_UP
    
55. SFGRAD, PRES, 0, Z, 0,-(W_WEIR/2+W_WALL)/W_WALL*9.81
    
56. SF,ALL,PRESS,Z_UP*(W_WEIR/2+W_WALL)/W_WALL*9.81 !施加边墩水压力,包含闸门传递的水压力
    
57. *ENDIF
    
58. 
    
59. *IF,Z_SAND,GT,ELEVATION(LAYER_NUM),THEN
    
60. SFCUM,PRES,ADD !荷载为累加方式
    
61. CMSEL,S,N_DAMUP
    
62. NSEL,R,LOC,Z,0,Z_SAND
    
63. *AFUN,DEG
    
64. SAND_GRADS=DENS_SAND*(TAN(45-ANG_FRI/2))**2*9.81 !计算淤沙压力梯度
    
65. SFGRAD, PRES, 0, Z, 0,-SAND_GRADS
    
66. SF,ALL,PRESS,Z_SAND*SAND_GRADS !施加淤沙压力
    
67. SFCUM,PRES,REPL !荷载为替代方式
    
68. *ENDIF
    
69. 
    
70. ALLSEL
    
71. *IF,Z_DOWN,NE,ELEVATION(LAYER_NUM),THEN
    
72. CMSEL,S,N_DAMDOWN
    
73. NSEL,R,LOC,Z,0,Z_DOWN
    
74. SFGRAD, PRES, 0, Z, 0,-9.81
    
75. SF,ALL,PRESS,Z_DOWN*9.81 !施加大坝下游面水压力
    
76. *ENDIF
    
77. 
    
78. ALLSEL
    
79. ACEL,0,0,1 !施加大坝自重荷载
    
80. SAVE

复制代码

AaronSpark

施加温度体荷载

1. !施加大坝温度荷载
   
2. CMSEL,S,DAM
   
3. NSLV,S,1
   
4. *GET,NCOUNT_DAM, NODE,, COUNT
   
5. *GET,NMIN_DAM, NODE,,NUM,MIN
   
6. TUNIF,0,
   
7. TREF,0,
   
8. NTT=NMIN_DAM
   
9. *DO,II,1,NCOUNT_DAM
   
10. Z_NTT=NZ(NTT) !Z
    
11. Z_TEMP=ELEVATION(1)-Z_NTT !到坝基高度
    
12. TL_ARCH=T2_ARCH(Z_TEMP) !拱厚
    
13. CSYS,11
    
14. X_NTT=NX(NTT) !X:距顶拱圆心距离
    
15. CSYS,0
    
16. RAD0=Y_RAD(Z_TEMP) !圆心到坝轴线距离
    
17. RAD1=RADUP(Z_TEMP) !外半径
    
18. LXUP=RAD_CEN(1)-RAD0+RAD1-X_NTT
    
19. TTT1=TEMPUP(Z_TEMP)
    
20. TTT2=TEMPDOWN(Z_TEMP)
    
21. TTT=TTT1+(TTT2-TTT1)/TL_ARCH*LXUP
    
22. BF,NTT,TEMP,TTT
    
23. *IF,II,LT,NCOUNT_DAM,THEN
    
24. NTT=NDNEXT(NTT)
    
25. *ENDIF
    
26. *ENDDO

复制代码

AaronSpark

计算求解，采用PCG求解器
ALLSEL
SAVE
EQSLV,PCG,1E-8
SOLVE
SAVE
FINI

最大变形2cm，位于拱冠梁顶，与拱梁分载法一致。

johhan

感谢AaronSpark的大力帮助!!!

请问您有没有双曲拱坝的命令流呀??

谢谢!!!!

AaronSpark

这个我没有，不过听说抛物线双曲拱坝采用浙江大学的拱坝程序不错

wjjerry

　　

如何考虑地震动水压力?

arcok

:lol from okok.org精华贴

wzg123

谢谢啊，浙大的拱坝优化程序的确很牛啊