亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于數(shù)值流形法的重力壩抗滑穩(wěn)定性分析

        2016-10-10 06:59:34陳遠(yuǎn)強
        長江科學(xué)院院報 2016年9期
        關(guān)鍵詞:重力壩流形勢能

        陳遠(yuǎn)強,鄭 宏,陳 濤

        (1.三峽大學(xué) 土木與建筑學(xué)院,湖北 宜昌 443002;2.北京工業(yè)大學(xué) 建筑工程學(xué)院,北京 100124)

        ?

        基于數(shù)值流形法的重力壩抗滑穩(wěn)定性分析

        陳遠(yuǎn)強1,鄭宏2,陳濤1

        (1.三峽大學(xué) 土木與建筑學(xué)院,湖北 宜昌443002;2.北京工業(yè)大學(xué) 建筑工程學(xué)院,北京100124)

        作為一種新的數(shù)值方法,數(shù)值流形法在重力壩抗滑穩(wěn)定性方面的研究較少。首先給出了水壓及揚壓力的荷載矩陣、安全系數(shù)的求解方法,然后采用數(shù)值流形法分析了重力壩沿建基面及深層雙斜面的抗滑穩(wěn)定性,得出了安全系數(shù),并與有限元接觸分析的結(jié)果進行了對比。分析結(jié)果表明:采用數(shù)值流形法和有限元接觸分析法得到的安全系數(shù)基本一致,從而驗證了數(shù)值流形法在重力壩抗滑穩(wěn)定分析中的可行性。

        數(shù)值流形法;重力壩;抗滑穩(wěn)定分析;安全系數(shù)

        1 研究背景

        混凝土重力壩是我國水利水電工程建設(shè)中廣泛采用的重要壩型之一,而抗滑穩(wěn)定分析又是重力壩設(shè)計中的一項重要內(nèi)容,因此采用合適的方法評價重力壩的穩(wěn)定性具有十分重要的意義。

        對于重力壩抗滑穩(wěn)定的分析方法,我國現(xiàn)行重力壩設(shè)計規(guī)范規(guī)定:以剛體極限平衡法計算為主,必要時可輔以有限元法等其它方法[1]。剛體極限平衡法操作簡單,具有長期實踐經(jīng)驗,計算得出的整體抗滑穩(wěn)定安全系數(shù)相對是可靠的,但不能模擬復(fù)雜地基對壩體穩(wěn)定性的影響,且無法得到壩體及地基內(nèi)應(yīng)力和位移情況[2];有限元方法能夠求得壩體及壩基內(nèi)的應(yīng)力分布、軟弱結(jié)構(gòu)面受力,以及破壞情況,但存在計算量大和計算難以收斂的問題,且有限元構(gòu)造的位移函數(shù)是協(xié)調(diào)的,無法很好地模擬不連續(xù)變形[3]。

        數(shù)值流形法(NumericalManifoldMethod,NMM)是石根華博士在塊體理論和非連續(xù)變形分析的基礎(chǔ)上創(chuàng)立的,該方法以拓?fù)淞餍魏臀⒎至餍螢榛A(chǔ),采用數(shù)學(xué)覆蓋和物理覆蓋,能夠統(tǒng)一解決工程中的連續(xù)與非連續(xù)變形問題[4]。該方法在處理各種復(fù)雜邊界和不連續(xù)界面等方面具有巨大優(yōu)勢,因此,在巖土工程中得到廣泛應(yīng)用[5-7]。

        本文推導(dǎo)了水壓及揚壓力的荷載矩陣,并基于數(shù)值流形法,對重力壩沿建基面及深層雙斜面的抗滑穩(wěn)定性進行了分析和評價。

        2 數(shù)值流形法的基本原理

        2.1有限覆蓋及覆蓋函數(shù)

        數(shù)值流形法的物理覆蓋系統(tǒng)由數(shù)學(xué)覆蓋和物理網(wǎng)格兩者組成,數(shù)學(xué)覆蓋由用戶選擇,只定義近似解的精度;物理網(wǎng)格是由材料體的邊界、裂縫、塊體和不同材料區(qū)域的交接面所構(gòu)成。物理覆蓋重疊的公共部分形成流形單元[8]。

        數(shù)值流形法在各個物理片上獨立定義局部位移函數(shù),即

        Ui=[ui(x,y)vi(x,y)]T=PiDi。

        (1)

        式中:Pi為控制局部位移函數(shù)的基函數(shù)矩陣;Di為廣義自由度。這些局部位移函數(shù)可以是常量、線性的、高階多項式或級數(shù)的有限項。

        局部位移函數(shù)通過權(quán)函數(shù)聯(lián)結(jié)在一起,構(gòu)成總體位移函數(shù),即

        (2)

        式中:we(i)(x,y)為權(quán)函數(shù);Ti為位移矩陣;q為物理片的個數(shù)。

        2.2平衡方程

        由給定的位移公式(2),可以求出系統(tǒng)的總勢能∏(包括應(yīng)變勢能、初應(yīng)力勢能、點荷載勢能、線荷載勢能、體荷載勢能、慣性力勢能、接觸彈簧的應(yīng)變勢能和摩擦力勢能),再根據(jù)最小勢能原理可建立平衡方程[9],即

        KD=F。

        (3)

        其中:

        (4)

        B=[LT1,LT2,…,LTq]。

        (5)

        式中:K為剛度矩陣;F為荷載矩陣,文獻[8]給出了常見荷載類型的計算方法;Ω為求解區(qū)域;E為彈性矩陣;L為應(yīng)變微分算子。

        3 水壓力及揚壓力作用時的荷載矩陣

        圖1 作用有水壓力的單元Fig.1 Element with water pressure

        原有的數(shù)值流形法中并沒有考慮水壓及揚壓力荷載的作用,對大壩進行抗滑穩(wěn)定分析需要考慮這些荷載的作用。在此,利用能量原理進行推導(dǎo)。

        設(shè)某流形單元的一條邊上作用有分布荷載,如圖1。

        對于分布荷載段上某點s(x,y),其參數(shù)方程可表示為

        (6)

        式中:(xk,yk),(xk+1,yk+1)分別為k,k+1點的坐標(biāo)。

        點s處的水壓力為

        (7)

        式中:qx,qy為對應(yīng)點處水壓力作用分量;ls,k為點s到點k的距離;lk,k+1為點k到點k+1的距離。

        水壓力產(chǎn)生的分布勢能可表示為

        (8)

        (9)

        當(dāng)單元上作用有揚壓力時,式(9)可同樣求出由揚壓力引起的單元節(jié)點荷載,揚壓力的分布如圖2所示。圖中H1,H2分別為上、下游水深;γ為水的重度;α為折減系數(shù)。

        (a)建基面揚壓力分布

        (b)深層雙斜面揚壓力分布圖2 揚壓力分布Fig.2 Distribution of uplift pressure

        4 安全系數(shù)的求解

        重力壩的整體安全狀況可通過抗滑穩(wěn)定整體安全系數(shù)來評價。計算安全系數(shù)時,參照傳統(tǒng)有限元方法,先運用數(shù)值流形法對重力壩進行應(yīng)力分析,根據(jù)得出的應(yīng)力計算出可能滑動面上的正應(yīng)力和剪應(yīng)力,從而可得出滑動面上的抗滑力和滑動力,各自進行疊加,再仿照傳統(tǒng)的剛體極限平衡法,可計算得到安全系數(shù)[10],即

        (10)

        式中:σi,τi為滑動面上第i個單元的正應(yīng)力和剪應(yīng)力;fi,ci為滑動面上第i個單元的抗剪斷摩擦系數(shù)和黏聚力;li為滑動面上第i個單元的長度;n為滑動面上單元的個數(shù)。

        5 工程算例

        5.1算例1:重力壩沿建基面的抗滑穩(wěn)定分析

        某混凝土重力壩壩高103m,壩頂寬14.8m,壩底寬度為70.24m,壩體斷面具體尺寸如圖3(a)所示。上游水位100m,下游無水。模型計算范圍為:壩基自壩踵向上游取1.5倍壩高,自壩址向下游取1.5倍壩高,壩基深度取1倍壩高,具體計算模型及數(shù)學(xué)覆蓋如圖3(b)所示。約束條件為:上下游基巖施加法向約束,底部施加全約束。計算荷載考慮大壩及壩基自重、水壓力以及揚壓力。重力壩壩體彈性模量30.0GPa,泊松比0.20,密度2 630kg/m3;壩基巖體彈性模量19.5GPa,泊松比0.25,密度2 720kg/m3;壩體與壩基接觸面的抗剪斷摩擦系數(shù)與黏聚力分別為0.9和1MPa。

        (a)壩體尺寸

        (b)計算模型及數(shù)學(xué)覆蓋圖3 算例1壩體尺寸與數(shù)學(xué)模型覆蓋Fig.3 Size of the dam and mathematical cover forexample 1

        數(shù)值流形法中采用加減彈簧來模擬不連續(xù)界面的接觸狀態(tài),故而本文用有限元接觸分析來與之進行結(jié)果比對,水平向與豎直向位移對比見圖4和圖5,具體結(jié)果對比見表1。由結(jié)果分析可知:2種算法得到的位移大小基本一致,最大位移出現(xiàn)的位置相同:最大水平位移均出現(xiàn)在壩頂處,最大豎向位移均出現(xiàn)在下游面處。同時,二者計算得出的安全系數(shù)差別很小,說明對重力壩采用數(shù)值流形法進行分析是可行的。

        (a)數(shù)值流形法

        (b)有限元接觸分析法圖4 不同算法水平位移對比Fig.4 Comparison of calculated horizontal displacementbetween NMM and FEM contact analysis

        (a)數(shù)值流形法

        (b)有限元接觸分析法圖5 不同算法豎直位移對比Fig.5 Comparison of calculated vertical displacementbetween NMM and FEM contact analysis表1 不同算法結(jié)果比較Table 1 Comparison of calculated results between NMMand FEM contact analysis

        算法最大水平位移/mm最大豎向位移/mm安全系數(shù)數(shù)值流形法22.8810.883.08有限元接觸分析法18.4111.053.14

        5.2算例2:重力壩深層雙斜面抗滑穩(wěn)定分析

        某重力壩非溢流壩段壩頂高程964.5m,壩高28m,壩頂寬8m,上游水位962.0m,下游水位941.0m?;鶐r內(nèi)存在發(fā)育斷層F5和小斷層f3,形成壩基深層抗滑穩(wěn)定的控制面,其中:斷層F5與水平面夾角為22°,小斷層f3與水平面夾角為28°,滑動示意圖如圖6(a)所示。

        模型計算范圍為:壩基長度自壩踵向上游取2倍壩底寬度,自壩址向下游取2倍壩底寬度,壩基深度取1倍壩高,具體計算模型及數(shù)學(xué)覆蓋見圖6(b)。約束條件為:上下游基巖施加法向約束,基巖底部施加全約束。計算荷載包括大壩及壩基自重、上下游水壓力及揚壓力,排水幕處的滲透壓力折減系數(shù)為0.35。

        壩體混凝土彈性模量為22GPa,泊松比為0.18,重度為24.0kN/m3;壩基巖體巖性單一,彈性模量為12.3GPa,泊松比為0.20,重度為25.4kN/m3;發(fā)育斷層F5和小斷層f3的抗剪斷摩擦系數(shù)分別為0.25和0.70,黏聚力分別為0.05MPa和0.40MPa。

        數(shù)值流形法得出的應(yīng)力云圖見圖7,與有限元接觸分析的具體結(jié)果對比見表2。經(jīng)比較分析:兩種算法的位移分布同樣是一致的,最大水平位移出現(xiàn)在斷層F5上部,最大豎向位移出現(xiàn)在壩頂處;安全系數(shù)相當(dāng)接近,進一步驗證了數(shù)值流形法對重力壩穩(wěn)定分析的適用性。

        (a)雙斜面滑動示意圖

        (b)計算模型及數(shù)學(xué)覆蓋圖6 算例2雙斜面滑動與數(shù)學(xué)模型覆蓋Fig.6 Sketch of double-inclined planes and mathematicalcover for example 2

        (a)水平向應(yīng)力

        (b)豎向應(yīng)力

        (c)切應(yīng)力圖7 數(shù)值流行法的應(yīng)力云圖Fig.7 Contours of stress calculated by NMM表2 不同算法結(jié)果比較Table 2 Comparison of calculated results between NMMand FEM contact analysis

        算法最大水平位移數(shù)值/mm位置最大豎向位移數(shù)值/mm位置安全系數(shù)數(shù)值流形法0.94F5上端1.75壩頂5.73有限元接觸分析法1.07F5上端1.84壩頂6.07

        6 結(jié) 論

        本文將數(shù)值流形法應(yīng)用到重力壩沿建基面及深層雙斜面抗滑穩(wěn)定分析的工程算例中,通過與有限元接觸分析的結(jié)果比較,得出以下結(jié)論:

        (1)數(shù)值流形法計算得出的位移和應(yīng)力分布與有限元接觸分析結(jié)果基本一致,體現(xiàn)了數(shù)值流形法對工程中連續(xù)與非連續(xù)問題的可靠性。

        (2)在數(shù)值流形法求解應(yīng)力和引入剛體極限平衡法求安全系數(shù)的基礎(chǔ)上,計算了重力壩的抗滑穩(wěn)定安全系數(shù),結(jié)果驗證了數(shù)值流形法在重力壩抗滑穩(wěn)定分析中的可行性。

        [1]SL319—2005,混凝土重力壩設(shè)計規(guī)范[S].北京:中國水利水電出版社,2005.

        [2]張國新,金峰.用DDA進行重力壩穩(wěn)定安全分析[J].水力發(fā)電學(xué)報,2004,(1):10-14.

        [3]王義鋒,章青.基于界面元法的向家壩重力壩深層抗滑穩(wěn)定分析[J].巖石力學(xué),2009,30(9):2691-2696.

        [4]SHIGen-hua.ManifoldMethodofMaterialAnalysis[C]∥Transactionsofthe9thArmyConferenceonAppliedMathematicsandComputing.ArmyResearchOffice,Minnesota,June18-21,1991:57-76.

        [5]王水林,葛修潤.流形元方法在模擬裂紋擴展中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報,1997,16(5):405-420.

        [6]王書法,朱維申,李術(shù)才,等.加錨巖體變形分析的數(shù)值流形方法[J].巖石力學(xué)與工程學(xué)報,2002,21(8):1120-1123.

        [7]張國新,金峰,王光綸.用基于流形元的子域奇異邊界元法模擬重力壩的地震破壞[J].工程力學(xué),2001,18(4):18-27.

        [8]石根華.數(shù)值流形方法與非連續(xù)變形分析[M].裴覺民,譯.北京:清華大學(xué)出版社,1997.

        [9]王芝銀,李云鵬.數(shù)值流形法中的幾點改進[J].巖土工程學(xué)報,1998,20(6):33-36.

        [10]張國新,金峰.重力壩抗滑穩(wěn)定分析中DDA與有限元方法的比較[J].水利發(fā)電學(xué)報,2004,23(1):11-14.

        (編輯:劉運飛)

        Analysis of Anti-sliding Stability of Gravity Dam Using Numerical Manifold Method

        CHEN Yuan-qiang1, ZHENG Hong2, CHEN Tao1

        (1.CollegeofCivilEngineering&Architecture,ChinaThreeGorgesUniversity,Yichang443002,China;2.CollegeofArchitectureandCivilEngineering,BeijingUniversityofTechnology,Beijing100124,China)

        TheNumericalManifoldMethod(NMM),asanewnumericalmethod,hasseldombeenutilizedintheanti-slidingstabilityanalysisofgravitydam.Inthisstudy,theloadmatrixofwaterpressureandupliftpressure,aswellasthesolutionofsafetyfactorweregiven.Furthermore,NMMwasadoptedtoanalyzetheanti-slidingstabilityofthefoundationplanesofgravitydamandthedeeply-laiddouble-inclinedplanes,thesafetyfactorwasthuscalculated,andwasfurthercomparedwiththecorrespondingresultsobtainedbythecontactanalysisinfiniteelementmethod(FEM).Theresultsdemonstratedthatthesafetyfactorscalculatedbytheabovetwomethods,namelytheNMMandFEM,werefundamentallyconsistent,whichverifiedthefeasibilityofapplyingNMMtotheanti-slidingstabilityanalysisofgravitydam.

        numericalmanifoldmethod;gravitydam;stabilityanalysisagainstsliding;safetyfactor

        2015-08-05;

        2015-09-09

        陳遠(yuǎn)強(1990-),男,湖北隨州人,碩士研究生,主要從事數(shù)值流行法方面的研究,(電話)18602707189(電子信箱)691205842@qq.com。

        鄭宏(1964-),男,湖北南漳人,研究員,博士生導(dǎo)師,主要從事計算巖土力學(xué)研究,(電話)13986011345(電子信箱)hzheng@whrsm.ac.cn。

        10.11988/ckyyb.20150655

        2016,33(09):133-137

        TV145

        A

        1001-5485(2016)09-0133-05

        猜你喜歡
        重力壩流形勢能
        “動能和勢能”知識鞏固
        作 品:景觀設(shè)計
        ——《勢能》
        文化縱橫(2022年3期)2022-09-07 11:43:18
        “動能和勢能”知識鞏固
        “動能和勢能”隨堂練
        緊流形上的Schr?dinger算子的譜間隙估計
        考慮各向異性滲流的重力壩深層抗滑穩(wěn)定分析
        迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
        Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
        豐滿混凝土重力壩防滲降壓灌漿處理工藝探討
        潰壩涌浪及其對重力壩影響的數(shù)值模擬
        蜜桃91精品一区二区三区| 亚洲av乱码专区国产乱码| 69搡老女人老妇女老熟妇| 亚洲六月丁香色婷婷综合久久| 免费的日本一区二区三区视频| 人妻无码一区二区视频| 亚洲午夜精品a区| 精品人妻一区二区三区蜜臀在线| 喷水白浆视频在线观看| 国产av无码专区亚洲av蜜芽| 色窝窝在线无码中文| 免费高清视频在线观看视频| 一区二区三区国产内射| 18禁无遮拦无码国产在线播放| 福利一区在线观看| 国产大陆av一区二区三区| 日韩精品视频免费在线观看网站| 国内精品久久久久久久97牛牛| 亚洲av无码乱码国产精品fc2| 校花高潮一区日韩| 三级日本理论在线观看| 亚洲乱色伦图片区小说| 伊人色综合视频一区二区三区| 久久久国产精品粉嫩av| 久久精品不卡一区二区三区| 欧美日韩精品久久久免费观看| 国产精品原创巨作AV女教师 | 蜜桃一区二区三区| 国语对白做受xxxxx在线中国 | 国产女主播在线免费看| 亚洲午夜无码毛片av久久| 午夜三级a三级三点| 亚洲不卡无码高清视频| 久久亚洲综合亚洲综合| а√天堂8资源中文在线| 99久久精品国产成人综合| 欧美日韩国产乱了伦| 人妻少妇艳情视频中文字幕| av综合网男人的天堂| 亚洲伊人久久一次| 国产精品成人久久a级片|