土石坝渗流（Seepage of earth rock dam）.docVIP

土石坝渗流（Seepage of earth rock dam） Seepage analysis of earth rock dam adopts seepage parameters of unsaturated soil and calculates seepage line iteratively, according to the previous calculation results, The coefficient of permeability and the point of wetting escape point are constantly modified until the accuracy requirement is met. The earth rock dam of this example is relatively simple. The unsaturated seepage is calculated, that is, the permeability coefficient is a function of the void pressure First, a data file PPPP.TXT is set up to store the permeability coefficient function as follows. The first is the void pressure (water head M), the second is the permeability coefficient, and the permeability coefficient is equal to 10^A (M/D). -10.00 -4.0E+00 -9.00 -3.6E+00 -8.00 -3.2E+00 -7.00 -2.8E+00 -6.00 -2.4E+00 -5.00 -2.0E+00 -4.00 -1.6E+00 -3.00 -1.2E+00 -2.00 -8.0E-01 -1.00 -4.0E-01 0 0.0E+00 Dam top width of 4M and the downstream slope is 1:2, height 12M, width 52M. The upstream water depth is 8M, and the downstream water is not. FINI /TITLE, EARTHDAM SEEPAGE /FILNAME, SEEPAGE5 /PLOPTS, DATE, 0! Clear time display *DIM, TPRE, TABLE, 11,1,1, PRESS, KKPE! Define the relationship array between hydraulic pressure and permeability coefficient *TREAD, TPRE, PPPP, TXT! Read the array *DIM, NCON, ARRAY, 4! Define the array for the four node number of the storage unit /PREP7 SMRT, OFF! Close default grid partition ANTYPE, STATIC! A static analysis ET, 1, PLANE55! 2D temperature field unit MP, KXX, 1,1 saturation coefficient of saturation, defines material properties MP, KXX, 2,1E-4! The coefficient of permeability under full drying, assuming that the void water pressure is less than -10M K, 1,24,12 K, 2,24,0 K, 3,0,0 K, 4,28,12 K, 5,28,0 K, 6,52,0 L, 1,3 L, 3,2 L, 1,2 L, 4,5 L, 5,6 L, 4,6 LESIZE, ALL,, 24 A, 1,3,2 A, 1,2,5,4 A, 4,5,6 MSHK, 2! MAPPED AREA MESH IF POSSIBLE MSHA, 0,2D! USING QUADS AMESH, ALL! MESH AREAS NUMMRG, NODE! MERGE NODES AT BOTTOM OF CAISSON! *GET,