资源描述:
《土石坝渗流(seepage of earth rock dam)》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、土石坝渗流(Seepageofearthrockdam)Seepageanalysisofearthrockdamadoptsseepageparametersofunsaturatedsoilandcalculatesseepagelineiteratively,accordingtothepreviouscalculationresults,Thecoefficientofpermeabilityandthepointofwettingescapepointareconstantlymodifieduntiltheaccuracyrequirementismet.Th
2、eearthrockdamofthisexampleisrelativelysimple.Theunsaturatedseepageiscalculated,thatis,thepermeabilitycoefficientisafunctionofthevoidpressureFirst,adatafilePPPP.TXTissetuptostorethepermeabilitycoefficientfunctionasfollows.Thefirstisthevoidpressure(waterheadM),thesecondisthepermeabilitycoef
3、ficient,andthepermeabilitycoefficientisequalto10^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-0100.0E+00Damtopwidthof4Mandthedownstreamslopeis1:2,height12M,width52M.Theupstreamwaterdepthis8M,andthed
4、ownstreamwaterisnot.FINI/TITLE,EARTHDAMSEEPAGE/FILNAME,SEEPAGE5/PLOPTS,DATE,0!Cleartimedisplay*DIM,TPRE,TABLE,11,1,1,PRESS,KKPE!Definetherelationshiparraybetweenhydraulicpressureandpermeabilitycoefficient*TREAD,TPRE,PPPP,TXT!Readthearray*DIM,NCON,ARRAY,4!Definethearrayforthefournodenumber
5、ofthestorageunit/PREP7SMRT,OFF!ClosedefaultgridpartitionANTYPE,STATIC!AstaticanalysisET,1,PLANE55!2DtemperaturefieldunitMP,KXX,1,1saturationcoefficientofsaturation,definesmaterialpropertiesMP,KXX,2,1E-4!Thecoefficientofpermeabilityunderfulldrying,assumingthatthevoidwaterpressureislessthan
6、-10MK,1,24,12K,2,24,0K,3,0,0K,4,28,12K,5,28,0K,6,52,0L,1,3L,3,2L,1,2L,4,5L,5,6L,4,6LESIZE,ALL,,24A,1,3,2A,1,2,5,4A,4,5,6MSHK,2!MAPPEDAREAMESHIFPOSSIBLEMSHA,0,2D!USINGQUADSAMESH,ALL!MESHAREASNUMMRG,NODE!MERGENODESATBOTTOMOFCAISSON!*GET,N_MAX,NODE,NUM,MAX!Getthemaximumnodenumber*GET,E_MAX,E
7、LEM,NUM,MAX!Getthemaximumunitnumber*DIM,N_TEMP,ARRAY,N_MAX!Definenodaltemperaturevariablestotalhead*DIM,N_PRE,ARRAY,N_MAX!DefinenodalpressureheadvariablesDefinethetotalheadvalueoftheupstreamfaceLSEL,S,LINE,1NSLL,S,1NSEL,R,LOC,Y,0,8D,ALL,TEMP,8!Definethetotalheadvalu