地下电站厂房层状围岩位移的动态效应研究

1,2,3,1,1,4(11,100083;21,271018;31,100053;41,100084):,,,,,,:;;;;::Researchondynamicdisplacementoflayercompositerock2massWANGKezhong1,2,ZHANGWenlei3,CAIMeifeng1,ZHANGYanxin1,DAIRong4(11SchoolofCivilandEnvironmentalEngineering,USTB,Beijing100083;21CollegeofHydraulicandCivilEngineering,ShandongAgricultureUniversity,Taian271018;31ChineseSocietyofHydraulicEngineering,Beijing100053;41DepartmentofHydraulicEngineering,TsinghuaUniversity,Beijing100084)Abstract:Thelayercompositerock2massissoftandhardalternately,itdoesnotmeettherequirementsofhomogeneityandisotropy.Thestabilityanddeformationofrock2massaredeterminedbyusingseveralmeasurementsduringtheexcavationprocess,thenumericalmodelofundergroundhousesinstratifiedanisotropicalrocksisanalyzedbyusingthreedimensionalnonlinearelementmodel.Thechanginglawofdisplacementsandstressesoflayercompositerock2massisstudied.Keywords:rock2massmechanics;stabilityanalysis;numericalsimulation;layercompositerock2mass;mechanicalparameters:2004205212:,1965,,1,;,,,2,3km,4300MW,,,SE,NW280,NW290340NE410,14715m21175m49m,45m,13019m1614m1715m,24420058JOURNALOFHYDROELECTRICENGINEERINGVol.24No.4Aug.,2005,,,,F118,F114P4,P5,,,,,3[1,2],,,,,,,:(1),,14,nn,:Vvi,jdv=Svinjds(1):V,S,nj,vi,nj,(1)vi,j=-13V4l=1vin(l)js(l)(2):ll,(l)l(2),,9vli9t=Fli(t)ml(3):Fli(t)tli,:vli(t+t)=vli(t+tP2)+[Fli(t)Pml]t(4)(3)(3),,(3)9vli9t=Fli(t)+fli(t)ml(5)(4)[3,4]:;;11Fig.1Computationsequenceofthree2dimensionalfinitedifferenceanalysis,,:ij=-16V4l=1(vlin(l)j+vljn(l)i)S(l)(6)t,H3ij:ij=N+cijNij=H3ij(ij,ij)(7),Nij,H3ij,ij:0012005ij=-t6V4l=1(vlin(l)j+vljn(l)i)S(l)(8):cij=(ikkj-ikkj)t(9)ij=-16V4l=1(vlin(l)j-vljn(l)i)S(l),v(l)i(t+t)=v(l)i(t)+tM(l)(F(l)i+f(l)i)(10):u(l)i(t+t)=u(l)i(t)+tv(l)i(t+t)(11),,,:x(l)i(t+t)=x(l)i(t)+tv(l)i(t+t)(12)Cauchy:ij,j+bi=dvidt(13),,b,dviPdt,0ij,j+bi=0(14)ii=(vi,j+vj,i)P2,9615:,3,15,4,32@115m115m,L=8m,1600kN@415m415m,L=25m;32@115m115m,L=8m,1600kN@415m6m,L=20m:,,,,,[5,6],,4(),,x()400m,y()251m,z()350m,1299214376,,8,Mohr2Coulomb[7],,,11Table1Mechanicalandphysicalparametersofrock2massEHR(GPa)EV(GPa)c(MPa)()(MPa)Rc(MPa)(kNPm3)151001201710353101102712151001201311041310110271101006014501622201001202010,(1),2192mm,,1014:2Fig.2Cross2sectionofhouseexcavation,,,,3mm,2716mm,,,(2),F112(NE40SE78),30m,F112,2112mm,F118,P4,P4NE35NW6787,1112cm,1210MPa,NE5060,610MPa,610710MPa310410MPa411M4-1(M4-1-1M4-1-2M4-1-3M4-1-4)2(0+000),0+0402Table2Positionsofdisplacementsinboltingcave(m)M4-1-10+0021500m0125m7891402004.1.8M4-1-20+0021500m3125m7921402004.1.8M4-1-30+0021500m8125m7981402004.1.8M4-1-40+0021500m17125m8081402004.1.8,:,,33Table3Displacementsatdifferentpoints(mm)1234M4-1-10+0021500m0125m1150211021652192M4-1-20+0021500m3125m1108115411952176M4-1-30+0021500m8125m0187018801971123M4-1-40+0021500m17125m01850186019511133,20,3,,,,,,3M4-1-12192mm,2187mm0105mm20120053Fig.3Processcurveofdifferentpointsinthecavern5(1),,,(2),,:[1],,.FLAC3D[R].:,2001.[2]MortiJ,CundallPA.Mixeddiscretizationprocedureforaccuratesolutionofplasticityproblem[J].Int.J.Num.MethodinEng.,1982,(6):129139.[3],,.FLAC[J].,1996,15(3):225227[4],,.[J].,2003,22(7):10471049[5],,.[J].,2000,19():952956[6],.[M].:,1994[7]QuinteiroCarlosRoberto.Numericalmodelingoftime2dependentgroundbehavioraroundadeepmineinpreparationforlayoutdesign[J].DAI2B53P07,1993:37403014: