變形加固理論及工程應(yīng)用

1、巖土與結(jié)構(gòu)工程數(shù)值方法課程專題巖體工程穩(wěn)定分析和變形加固理論及工程應(yīng)用清華大學(xué)水利系巖石力學(xué)與工程研究中心Mar. 6, 2006TSINGHUA UNIVERSITY1 概述2 巖體工程的穩(wěn)定分析3 變形加固理論 4 大規(guī)模數(shù)值計(jì)算技術(shù)5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究TSINGHUA UNIVERSITY內(nèi)容提綱TSINGHUA UNIVERSITY1 概述我國水電資源集中在西部地區(qū)，該地區(qū)具有全球最強(qiáng)烈的現(xiàn)代地殼活動和高地應(yīng)力場，其特殊復(fù)雜的地質(zhì)、地形條件，以及超大規(guī)模長期安全、穩(wěn)定運(yùn)行的要求，使得西部地區(qū)巨型水電樞紐大量采用高拱壩、地下廠房和深埋隧洞引水發(fā)電的布置形式： 溪洛渡

2、(278m)：地下廠房 小灣(292m)：壩基地應(yīng)力卸荷、壩肩基礎(chǔ)處理 錦屏一級(305m)：高邊坡、壩肩基礎(chǔ)處理、地下廠房 錦屏二級：深埋長引水隧洞，17km，2500m埋深 .高壩地基系統(tǒng)的開裂破壞、壩基的優(yōu)化選擇、巖石高邊坡的穩(wěn)定和加固、深埋地下洞室等是這些工程中經(jīng)常遇到的問題。 TSINGHUA UNIVERSITY1 概述研究背景高拱壩壩體體型優(yōu)化建基面優(yōu)化選擇壩基加固處理高拱壩設(shè)計(jì)核心提高拱壩抗裂及整體穩(wěn)定性 工程實(shí)踐表明，大壩的失穩(wěn)或破壞大多是由于壩的基巖或壩體建基面附近存在缺陷，從而引起壩上的一定部位，尤其是上游壩踵和壩與基巖的相接處，產(chǎn)生拉壓應(yīng)力的集中，導(dǎo)致拱壩局部開裂，進(jìn)而

3、引起大壩開裂，并造成壩體的最終失穩(wěn)破壞。 TSINGHUA UNIVERSITY1 概述奧地利Kolnbrein壩在1977年建成后第一次蓄水時即出現(xiàn)上游壩踵處開裂而漏水的問題 TSINGHUA UNIVERSITY1 概述我國的二灘拱壩下游壩面出現(xiàn)的裂縫 TSINGHUA UNIVERSITY1 概述地基失穩(wěn)是引起拱壩破壞的另一重要原因1959年法國Malpasset拱壩因降雨量大和地基中深層裂隙的發(fā)展，使壩體連同地基發(fā)生深層滑動而潰壩，死亡400余人。 TSINGHUA UNIVERSITY1 概述研究背景高邊坡工程天然邊坡主要指水庫庫岸邊坡等，該類邊坡在水庫蓄水后，受水的影響，容易發(fā)生庫

4、岸滑坡，造成涌浪，影響大壩和附近居民的安全。 意大利的Vajoint大壩由于邊坡的滑動導(dǎo)致大壩的不能使用 TSINGHUA UNIVERSITY1 概述人工邊坡則包括大壩開挖形成的壩肩邊坡、船閘邊坡、電站廠房邊坡等。 錦屏一級水電站高邊坡 TSINGHUA UNIVERSITY1 概述研究背景深埋地下洞室典型問題： 高地應(yīng)力 高滲透壓力 高地溫 圍巖穩(wěn)定分析 支護(hù)設(shè)計(jì) 地應(yīng)力模擬 穩(wěn)定性評價 支護(hù)時機(jī)選擇 反饋分析 .數(shù)值分析和地質(zhì)力學(xué)模型試驗(yàn)存在的問題TSINGHUA UNIVERSITY1 概述水工巖體工程的需求常規(guī)的數(shù)值計(jì)算方法和常規(guī)的巖體本構(gòu)關(guān)系都不能滿足要求。需要建立新的數(shù)值模擬手段

5、和新的本構(gòu)關(guān)系，通過新的評價體系來解決這類問題。因此，對這些巖體工程的穩(wěn)定性和加固措施進(jìn)行深入研究，為設(shè)計(jì)和施工提供可靠的依據(jù)，對巖石力學(xué)理論在水工巖體工程中的應(yīng)用的發(fā)展都具有重大的意義。 水工巖體工程特點(diǎn)：規(guī)模大 地質(zhì)條件復(fù)雜 不確定因素多要求：使用期長 安全度高TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述高拱壩混凝土和壩基巖體的開裂分析 拱壩失效破壞可以分為兩個階段：開裂和失穩(wěn)拱壩開裂破壞的原因很復(fù)雜：干縮、凍融、溫度變化、超常荷載、水力劈裂、材料質(zhì)量惡化、地基異常變形等，都是可能引起破壞的因素。 目前，拱壩設(shè)計(jì)中沒有明確的抗裂設(shè)計(jì)方法和原則，抗裂要求主要通過控制拉應(yīng)力區(qū)

6、、損傷塑性區(qū)和優(yōu)化體型來體現(xiàn)。 斷裂力學(xué)引入巖石力學(xué)，使得分析巖體中的裂紋擴(kuò)展問題成為可能。目前主要集中于二維問題或者簡單三維試件的研究，但如何對三維裂紋擴(kuò)展進(jìn)行有效的數(shù)值模擬還遠(yuǎn)遠(yuǎn)沒有得到解決，距實(shí)現(xiàn)高拱壩壩開裂過程的仿真分析尚有一定的距離傳統(tǒng)有限單元等方法的裂縫追蹤能力不強(qiáng)，商業(yè)軟件如ANSYS，有判斷開裂區(qū)的功能，但不能對裂縫進(jìn)行定位與追蹤；有些新方法如無網(wǎng)格法雖然在開裂分析方面具有一定的優(yōu)勢，但是還不能真正有效處理三維問題，距離實(shí)際工程應(yīng)用還有很多工作要做。 TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述高拱壩混凝土和壩基巖體的開裂分析 現(xiàn)有的K1，K2，K3的高拱壩開

7、裂和整體穩(wěn)定評價體系張楚漢，工程類比法，以二灘拱壩作為不開裂的基準(zhǔn)工程來評價其他拱壩的開裂和穩(wěn)定?；谝郧敖?jīng)驗(yàn)和已有工程評價拱壩開裂和穩(wěn)定 基于網(wǎng)格加密和大規(guī)模計(jì)算技術(shù)的開裂分析方法。從思路上是可行的，但是對高拱壩這樣空間尺度很大的三維復(fù)雜結(jié)構(gòu)而言，實(shí)際應(yīng)用上卻存在計(jì)算量非常巨大、計(jì)算可能不收斂的問題。 裂紋追蹤和模擬的可行方法：結(jié)合使用網(wǎng)格自適應(yīng)方法、非局部化損傷方法，采用大規(guī)模計(jì)算技術(shù)來追蹤裂紋擴(kuò)展是可行的。 TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述 從拱壩裂縫形成的機(jī)理上來看，拱壩開裂多集中和發(fā)源于拱端附近的應(yīng)力集中區(qū)，這些應(yīng)力集中區(qū)進(jìn)入壩肩巖體內(nèi)衰減很快，形成高梯

8、度應(yīng)力區(qū)。拱肩槽開挖、隧洞開挖導(dǎo)致的巖體卸荷也會形成高梯度應(yīng)力。高梯度應(yīng)力區(qū)而非高應(yīng)力區(qū)是使混凝土和巖體發(fā)生開裂的根源。高拱壩和加固措施的設(shè)計(jì)中要盡量避免出現(xiàn)突變，不論是幾何形狀突變，還是材料裂隙性質(zhì)的突變，都是非常不利的。實(shí)際上，這種突變是導(dǎo)致高梯度應(yīng)力區(qū)的重要根源。要反應(yīng)高梯度應(yīng)力的影響，需要在巖石本構(gòu)關(guān)系中引入非線性、各向異性、多尺度、多相耦合、局部化分叉、非局部化處理及內(nèi)稟長度等概念，這在工程應(yīng)用上具有一定的困難。在目前的研究水平上，通過深化已有方法來引入這些概念，可能更具有工程使用價值。目前，關(guān)于高梯度應(yīng)力方面的研究文獻(xiàn)還相對比較少。 高拱壩混凝土和壩基巖體的開裂分析 TSINGHU

9、A UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn) 巖體工程的穩(wěn)定評價 抗滑穩(wěn)定變形穩(wěn)定 傳統(tǒng)方法：抗滑穩(wěn)定剛體極限平衡地質(zhì)力學(xué)模型試驗(yàn)和三維非線性有限元計(jì)算可以更真實(shí)的反映結(jié)構(gòu)和地基的聯(lián)合工作情況以及其中的非線性應(yīng)力調(diào)整，但是只能給出位移、應(yīng)力和點(diǎn)安全度的分布，其中點(diǎn)安全度反映的是局部點(diǎn)的安全度，而不是整體穩(wěn)定情況。 變形穩(wěn)定目前尚無明確的定義和控制指標(biāo)，而從地質(zhì)力學(xué)模型超載試驗(yàn)經(jīng)常可以觀察到大壩破壞是由于壩肩變形過大導(dǎo)致大壩喪失承載力造成的。 拱壩建基面優(yōu)化分析TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn) 巖體

10、工程的穩(wěn)定評價 壩肩和邊坡穩(wěn)定：剛體極限平衡法和基于有限元的強(qiáng)度折減法基于多重網(wǎng)格法，基于非線性有限元的邊坡穩(wěn)定分析方法，更能反映巖體的非線性變形，結(jié)合了抗滑穩(wěn)定和變形穩(wěn)定整體穩(wěn)定：能量方法的使用更為合理破壞后階段的數(shù)值模擬問題：解決方法是采用精細(xì)的荷載步、采用反映破壞后階段的幾何非線性和損傷局部化特征的數(shù)值方法進(jìn)行。大規(guī)模數(shù)值計(jì)算技術(shù)TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn) 拱壩建基面的應(yīng)力控制標(biāo)準(zhǔn) 有限元數(shù)值計(jì)算結(jié)果的網(wǎng)格依賴性要求有限元計(jì)算成果的規(guī)范化應(yīng)力的線性等效應(yīng)力的非線性等效自適應(yīng)方法解決方案其中，要在工程中得到實(shí)際應(yīng)用，非線性

11、等效和自適應(yīng)方法均需采用大規(guī)模計(jì)算技術(shù)TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn) 加固設(shè)計(jì)和加固效果的評價 常規(guī)數(shù)值方法：模擬效果甚微 錨索、墊腳加錨加固措施的效應(yīng)和加固力對結(jié)構(gòu)變形特征極為敏感 變形加固理論：變形、穩(wěn)定、加固 基于抗滑穩(wěn)定的加固設(shè)計(jì)比較明確，但滑動失穩(wěn)是一種處在超載或降強(qiáng)下的極端狀態(tài)，偏離正常工作狀態(tài)很遠(yuǎn)。著眼于正常工作狀態(tài)的基于變形穩(wěn)定的加固設(shè)計(jì)可能更具有工程實(shí)際意義，但到目前為止，還沒有一套明確的設(shè)計(jì)方法和準(zhǔn)則。 TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn) 加固設(shè)計(jì)和加

12、固效果的評價 變形（不平衡力）加固理論（楊強(qiáng)，2003）根據(jù)不平衡力的大小、方向、分布，即可有效地確定加固力 解決收斂性問題加固力大小、方向均能確定同時考慮了抗滑穩(wěn)定和變形穩(wěn)定存在問題： 理論、數(shù)值模擬結(jié)果與試驗(yàn)分析的驗(yàn)證？ 與規(guī)范的配套使用？TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn) 深埋地下洞室的穩(wěn)定和支護(hù)問題 塑性松動圈：確定和數(shù)值模擬 巖爆：判別標(biāo)準(zhǔn)，極限拉應(yīng)變、最大主壓應(yīng)力/最小主壓應(yīng)力 流變問題：工程實(shí)用上的問題 支護(hù)時機(jī)選擇：新奧法，和變形的聯(lián)系 噴錨支護(hù)的數(shù)值模擬TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述巖

13、體工程穩(wěn)定和加固分析的評價標(biāo)準(zhǔn)降強(qiáng)和超載 超載分析突出了拱端推力作為主要滑動力的地位，這和常規(guī)剛體極限平衡分析的精神是一致的。 降強(qiáng)分析則強(qiáng)調(diào)了巖體、結(jié)構(gòu)面力學(xué)參數(shù)不確定性的一面。 綜合強(qiáng)度指標(biāo)將兩者結(jié)合起來：先降強(qiáng)至一定程度，再超載。材料參數(shù)降低后還能承受多大的荷載。比較符合實(shí)際情況。 巖石的強(qiáng)度應(yīng)力比 強(qiáng)度的尺寸效應(yīng)、數(shù)值模擬中的應(yīng)力取值問題？ 對巖體工程這樣復(fù)雜的結(jié)構(gòu)體系追求單一的安全指標(biāo)是不適宜的，無法全面反映客觀存在的可能多種失穩(wěn)模式，應(yīng)該建立反映破壞多尺度特點(diǎn)的安全指標(biāo)體系，如點(diǎn)、面、塊、整體安全度。這也是所謂的“工程穩(wěn)定廣角度準(zhǔn)則”的精神（卓家壽） TSINGHUA UNIVER

14、SITY1 概述關(guān)鍵技術(shù)問題綜述地應(yīng)力場的數(shù)值模擬 參數(shù)反分析 直接的地應(yīng)力場多元回歸法應(yīng)力函數(shù)法由測點(diǎn)應(yīng)力值直接推算應(yīng)力場 計(jì)算量大收斂速度慢結(jié)果不夠精確 逆反分析方法直接反分析方法：解析方法、隨機(jī)方法和枚舉法 大規(guī)模數(shù)值計(jì)算技術(shù)TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述滲流場模擬及和應(yīng)力場的耦合分析 拱壩壩基開裂 巖石高邊坡 地下洞室的穩(wěn)定 滲流場和應(yīng)力場的耦合方法 全耦合 迭代耦合：實(shí)用方法 顯式耦合 偽耦合 關(guān)鍵問題：應(yīng)力場對滲流場的影響，對于雙重孔隙介質(zhì)模型，則主要是應(yīng)力場對裂隙巖體的滲透系數(shù)的影響滲流場的數(shù)值模擬 雙重孔隙介質(zhì)模型（DP） 多重孔隙介質(zhì)模型 滲透

15、系數(shù)張量理論 離散裂隙網(wǎng)絡(luò)模型（DFN）DP和DFN的結(jié)合TSINGHUA UNIVERSITY1 概述關(guān)鍵技術(shù)問題綜述大規(guī)模數(shù)值計(jì)算和三維數(shù)字建模技術(shù) 三維有限元的并行計(jì)算在水利工程中的應(yīng)用還不普及，大空間尺度水工結(jié)構(gòu)的有限元分析精度較低，同時，破壞后階段的精細(xì)模擬（采用很精細(xì)的荷載步長）不能實(shí)現(xiàn)。 三維數(shù)字建模技術(shù)、計(jì)算機(jī)圖形理論和虛擬現(xiàn)實(shí)技術(shù)相對比較成熟，但是在水利工程中的應(yīng)用還處于一個起步階段，造成的現(xiàn)狀是有限元的前處理工作巨大、后處理不能實(shí)現(xiàn)真正的人機(jī)交互和圖形圖像的實(shí)時處理。 由于計(jì)算能力的限制，三維有限元損傷局部化理論中的關(guān)鍵參數(shù)（內(nèi)稟長度因子、損傷演化方程中的參數(shù)等）的確定和取

16、值規(guī)律研究都難以進(jìn)行，成為限制其在工程中應(yīng)用的關(guān)鍵因素。 對并行計(jì)算的認(rèn)識誤區(qū)：認(rèn)為這是計(jì)算機(jī)專業(yè)搞的東西，其實(shí)這是非常錯誤的認(rèn)識。并行計(jì)算是和具體求解問題關(guān)系非常密切的研究方向，只有水工專業(yè)的研究人員參與，才能更有效的解決實(shí)際問題。否則，將導(dǎo)致先進(jìn)的計(jì)算機(jī)技術(shù)在水利工程中的應(yīng)用越來越落后。TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析穩(wěn)定分析內(nèi)容： 拱壩壩肩穩(wěn)定分析 邊坡的穩(wěn)定分析 地下洞室的穩(wěn)定分析 .目的：為加固分析做準(zhǔn)備TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析整體穩(wěn)定分析一般通過非線性計(jì)算是否收斂和屈服區(qū)進(jìn)行綜合判斷。但是，當(dāng)材料進(jìn)入非線性后，一般很難

17、保證絕對收斂。巖體工程的特點(diǎn)及基礎(chǔ)的不均質(zhì)特性容易造成三維畸形網(wǎng)格、高度不均勻的材料分區(qū)、高量級荷載，這都對收斂性有很不利的影響。難以從收斂性上判斷穩(wěn)定。能量方法可能是發(fā)展趨勢。評價標(biāo)準(zhǔn)比較客觀。 ？主要問題：數(shù)值計(jì)算結(jié)果的整理和使用。TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析塊體的穩(wěn)定分析極限分析兩大類方法： 剛體極限平衡法（規(guī)范中的方法）通過分析巖體滑塊在破壞那一刻的平衡來求得滑塊的滑動安全系數(shù)，來判斷邊坡是否穩(wěn)定；該方法概念清晰、計(jì)算簡單，但是不能考慮結(jié)構(gòu)的變形效應(yīng)，而且假定滑動面同時達(dá)到極限破壞狀態(tài)，不能反應(yīng)滑動面的真實(shí)狀態(tài)。 有限元法（發(fā)展趨勢） 有限元方法能給出邊

18、坡中的應(yīng)力場和位移場，但是無法給出明確的安全系數(shù)。 強(qiáng)度折減法綜合判斷方法 在一定程度上是對有限元成果的規(guī)范化和應(yīng)用多重網(wǎng)格法TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析整體穩(wěn)定分析余能范數(shù)塑性應(yīng)力增量：余能范數(shù)：用來衡量整體穩(wěn)定TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析整體穩(wěn)定分析余能范數(shù)TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析多重網(wǎng)格法的思路剛體極限平衡：通過分析巖體滑塊在破壞那一刻的平衡來求得滑塊的滑動安全系數(shù)，來判斷邊坡是否穩(wěn)定；該方法概念清晰、計(jì)算簡單，但是不能考慮結(jié)構(gòu)的變形效應(yīng)，而且假定滑動面同時達(dá)到極限破壞狀態(tài)，不能反應(yīng)滑動面

19、的真實(shí)狀態(tài)。 強(qiáng)度折減系數(shù)法：根據(jù)某種破壞判據(jù)來判定系統(tǒng)是否進(jìn)入極限平衡狀態(tài)，例如限定迭代次數(shù)、限定節(jié)點(diǎn)不平衡力與外荷載的比值大小等。而高拱壩的工作特點(diǎn)及壩基不均質(zhì)的特性容易造成三維畸形網(wǎng)格、高度不均勻的材料分區(qū)、高量級荷載，這都對收斂性有很不利的影響。難以從收斂性上判斷穩(wěn)定。 將極限分析和有限元結(jié)合起來，既可以利用有限元的結(jié)果，考慮非線性應(yīng)力調(diào)整和變形的影響，又能得到滑塊的抗滑穩(wěn)定安全系數(shù)。也就是說，能同時考慮到變形穩(wěn)定和抗滑穩(wěn)定 TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析3個結(jié)構(gòu)面（上游側(cè)滑面F1、下游側(cè)滑面F2和上游拉裂面F3）與臨空面構(gòu)成的塊體。失

20、穩(wěn)巖體主要承受的荷載包括巖塊自重W，作用在結(jié)構(gòu)面上的滲透壓力U1、U2、U3和巖體抗力。考慮滑動塊體沿交線om的滑動模式，此時可假設(shè)拉裂面nol面拉開，不存在法向反力R3和切向力S3，但存在滲透壓力U3，并假設(shè)上游側(cè)滑面和下游側(cè)滑面同時達(dá)到極限狀態(tài)且切向力S1、S2平行于交線。這些假定意味著只有3個未知數(shù)，R1，R2和S1S2，由三個力的平衡條件可以確定它們。TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析穩(wěn)定安全系數(shù)Kc 用有限元法或者有限差分法計(jì)算塊體受力時，已經(jīng)考慮了滲透壓力的作用，故而上式中的U1、U2不予計(jì)入。另外，計(jì)算出的上游側(cè)滑面及下游側(cè)滑面切向力

21、S1、S2的方向一般都不平行于交線om。因此，計(jì)算滑塊安全度時，要根據(jù)二者夾角，投影到om上再進(jìn)行計(jì)算。 TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析考慮材料分區(qū)TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析滑面的安全系數(shù)用于分析滑動面上的安全情況TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析多重網(wǎng)格法核心思想：結(jié)構(gòu)網(wǎng)格和滑面網(wǎng)格相互獨(dú)立，滑面上各個結(jié)點(diǎn)的應(yīng)力由結(jié)構(gòu)網(wǎng)格上高斯點(diǎn)的應(yīng)力插值得到 插值方法：n的取值：n2即可滿足要求。TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基

22、于有限元的極限分析多重網(wǎng)格法要求滑面P的合力，需先將每一個四邊形單元分為兩個三角形，單元力矢量為兩個三角形上力矢量的和。單元中任意一個三角形上的力矢量按下式計(jì)算： TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析多重網(wǎng)格法將滑面上各單元法向矢量加和求平均，即可得到滑面的平均法向矢量，然后將各單元的力矢量加和得到的合矢量投影到滑面平均法向上，即可得到滑面的法向合力，另一分量即為滑面的切向合力。由此可確定各滑面上的面安全度及塊體的整體安全度。 TSINGHUA UNIVERSITY2 巖體工程的穩(wěn)定分析基于有限元的極限分析精度的影響因素 計(jì)算網(wǎng)格的密度網(wǎng)格越密，精度

23、越高。地形和地質(zhì)條件的精細(xì)模擬對大規(guī)模計(jì)算技術(shù)提出很高的要求。 結(jié)構(gòu)網(wǎng)格的密度與計(jì)算網(wǎng)格尺寸相當(dāng)時精度最高。TSINGHUA UNIVERSITY3 變形加固理論 變形加固理論綜述穩(wěn)定安全分析的主要目的之一是指導(dǎo)加固設(shè)計(jì)。對巖體工程之類的高次超靜定結(jié)構(gòu)，加固措施的效應(yīng)和加固力對結(jié)構(gòu)變形特征極為敏感。如擋土墻的土壓力可視為擋土墻對填土的加固力，擋土墻土壓力的發(fā)揮在很大程度上取決于擋土墻和填土的相對位移（主動、靜止、被動土壓力）。新奧法施工適時支護(hù)原則也說明這一點(diǎn)。不同的變形狀態(tài)對應(yīng)著不同的加固力，這和剛體極限平衡法的基本假設(shè)也是不相容的。有效的加固措施必須著眼于結(jié)構(gòu)自承能力的充分發(fā)揮。結(jié)構(gòu)自承能

24、力的發(fā)揮需要通過變形來實(shí)現(xiàn)，自承能力的發(fā)揮又減少了對加固力的依賴，這也是新奧法施工的原理，擋土墻土壓力理論也說明了這一點(diǎn)。變形、穩(wěn)定、加固力三者是緊密聯(lián)系在一起的，孤立地談一方面或兩方面是毫無意義的。常規(guī)剛體極限平衡法分析不能反映巖土結(jié)構(gòu)和基礎(chǔ)的非線性相互作用，且不能考慮變形機(jī)制，而且針對的是一特定的破壞狀態(tài)，以此指導(dǎo)加固設(shè)計(jì)是不全面，也肯定不是最有效的方法。 TSINGHUA UNIVERSITY3 變形加固理論 基于有限元的極限分析變形加固理論三維非線性有限元和有限差分分析無疑是更全面的一種分析方法。但如何以變形體彈塑性分析指導(dǎo)加固設(shè)計(jì)，目前尚沒有成熟的規(guī)則可循，至少以應(yīng)力、位移、塑性區(qū)的

25、變化來指導(dǎo)加固設(shè)計(jì)現(xiàn)在看起來不是很成功（相當(dāng)大的程度上是依靠經(jīng)驗(yàn)）。這和變形體的穩(wěn)定特點(diǎn)有關(guān)。 在連續(xù)介質(zhì)力學(xué)里，平衡和穩(wěn)定條件都被連續(xù)化為成場的概念，分別對應(yīng)于微元體的平衡微分方程和屈服條件，兩者都要求在結(jié)構(gòu)內(nèi)逐點(diǎn)滿足。 有限元法則通過變分原理把連續(xù)的平衡和穩(wěn)定條件凝聚到結(jié)點(diǎn)上，要求逐節(jié)點(diǎn)滿足，所以不平衡力可以視為屈服條件在結(jié)點(diǎn)上的集中體現(xiàn)。故在有限元分析的框架內(nèi)，連續(xù)結(jié)構(gòu)體的穩(wěn)定條件就是要求所有結(jié)點(diǎn)不平衡力均為零。如果經(jīng)過反復(fù)迭代計(jì)算，仍有殘余的不平衡力，就意味著必須施加一個額外的力（它和殘余的不平衡力大小相等、方向相反）才能維持平衡，這就是加固力。 TSINGHUA UNIVERSITY

26、3 變形加固理論 基于有限元的極限分析變形加固理論TSINGHUA UNIVERSITY3 變形加固理論 基于有限元的極限分析變形加固理論上面是在某一迭代步（即某一變形狀態(tài)）下靜力容許應(yīng)力場、外荷載和加固力之間的關(guān)系。對于一給定的結(jié)構(gòu)無法承擔(dān)的外荷載，在連續(xù)的迭代過程中，變形、應(yīng)力場，加固力不斷發(fā)生變化，每一個迭代步所對應(yīng)的不平衡力都是一個可行的加固方案（相應(yīng)于一定的變形）?？梢愿鶕?jù)變形、穩(wěn)定的要求來確定加固力：一般應(yīng)選取最小的不平衡力作為加固力。加固力不是惟一的，彈塑性計(jì)算所得的殘余的不平衡力所對應(yīng)的加固力是最優(yōu)的一組加固力（加固力最?。冃慰赡茌^大；如要限制變位（如拱端變位）或滿足其它

27、約束條件，也可得到一組加固力，但肯定比最優(yōu)加固力為大。 TSINGHUA UNIVERSITY4 大規(guī)模計(jì)算技術(shù)必要性 靜力問題 巖體中裂隙和軟弱斷層等各類地質(zhì)構(gòu)造、大型地下洞室的穩(wěn)定、地下廠房的構(gòu)造、尾水管內(nèi)墩子以及地應(yīng)力的計(jì)算，都要求加密網(wǎng)格，以提高精度。動力問題 動力的時程分析計(jì)算量巨大，例如，要計(jì)算20秒鐘的時程反映，一般要進(jìn)行幾千個時步的計(jì)算。如果考慮各點(diǎn)之間的相位差以及材料和幾何非線性，計(jì)算量更是大得驚人。數(shù)值模擬問題已經(jīng)成為理論進(jìn)一步發(fā)展的障礙。結(jié)構(gòu)的破壞和局部化分析 考慮巖土材料的非均勻性和非線性，進(jìn)行大位移、開裂計(jì)算以及局部化分析，都需要進(jìn)行大規(guī)模的計(jì)算。有限元并行計(jì)算的分類

28、有限元計(jì)算中的并行性分析 數(shù)據(jù)準(zhǔn)備單元分析集成整剛求解方程組應(yīng)力求解并行性數(shù)據(jù)準(zhǔn)備可完全并行分塊法和波前法 直接法和迭代法可完全并行4 大規(guī)模計(jì)算技術(shù)有限元并行計(jì)算的分類分類 求解并行只對方程組的求解采用并行方法，而對其他步驟均采用傳統(tǒng)的串行方法區(qū)域分解 本質(zhì)上是子結(jié)構(gòu)方法。EBE（Element-by-element）方法 將一個總體矩陣的向量積轉(zhuǎn)化到一組單個矩陣的向量積的計(jì)算。不再集成整剛和整體荷載向量，所有計(jì)算均在單元一級上進(jìn)行。 優(yōu)缺點(diǎn) 求解并行并行實(shí)現(xiàn)容易，但當(dāng)規(guī)模增大時，整剛集成的并行實(shí)現(xiàn)困難區(qū)域分解 在剖分較好時并行計(jì)算效率高。但三維問題的區(qū)域分解算法還不成熟EBE方法 可以解決

29、內(nèi)存和并行實(shí)現(xiàn)的問題，使用方便，不過效率可能比區(qū)域分解方法稍低4 大規(guī)模計(jì)算技術(shù)PFEM程序區(qū)域分解區(qū)域分解算法：SHANG-HSIEN HSIEH. EVALUATION OF AUTOMATIC DOMAIN PARTITIONING ALGORITHMS FOR PARALLEL FINITE ELEMENT ANALYSIS. International Journal for Numerical Methods In Engineering, 1997,40:10251051FETI（FiniteElement Tearing and Interconnection） 方法：K. C

30、. Park, Manoel R. Juitino, JR., Carlos A. Felippa. An Algebraically Partitioned FETI Method For Parallel Structural Analysis:Algorithm Description. International Journal for Numerical Methods In Engineering, 1997,40:27172737Manoel R. Juitino., JR, K. C. Park, Carlos A. Felippa. An Algebraically Part

31、itioned FETI Method For Parallel Structural Analysis:Performance Evaluation. International Journal for Numerical Methods In Engineering, 1997,40:27392758Farhat C., Lesoinne M, LeTallec P, Pierson K, Rixen D. FETI-DP: a dual-primal unified FETI method-part I: a faster alternative to the two-level FET

32、I method. International Journal for Numerical Methods In Engineering, 2001,50(7):15231544Daniel J. Rixen. Extended preconditioner for the FETI method applied to constrained problems. International Journal for Numerical Methods In Engineering, 2002,54:1-264 大規(guī)模計(jì)算技術(shù)并行計(jì)算的研究現(xiàn)狀Cray X-MPCray Y-MPFujitsu V

33、PP500YH-2 IBM SP2 網(wǎng)絡(luò)機(jī)群高性能計(jì)算機(jī)的通訊系統(tǒng)：100M1000M以太網(wǎng)、Myrinet、SCI、Quadrics、Infiniband 高性能所機(jī)群系統(tǒng) 網(wǎng)絡(luò)機(jī)群4 大規(guī)模計(jì)算技術(shù)有限元并行計(jì)算的研究現(xiàn)狀國外的發(fā)展: BLACS（Basic Linear Algebra Communication Subprograms）：BLAS的一個并行版本 LINPACK：LAPACK的并行版本 PIM（Parallel Iterative Methods）：Fortran 77的一個并行函數(shù)庫 PETSC（Portable Extensible Toolkit for Scien

34、tific Computation）最近兩年國內(nèi)并行軟件的發(fā)展也取得很大的發(fā)展:飛箭公司：武漢巖土所：Rsm VPC集成化并行計(jì)算環(huán)境（基于Windows和PVM）石油科學(xué)院：做了很多工作，具體不詳清華大學(xué)：TFINE.Pfem程序，巖體工程結(jié)構(gòu)分析并行計(jì)算的軟件方面4 大規(guī)模計(jì)算技術(shù)有限元并行計(jì)算的研究現(xiàn)狀通用有限元軟件方面：MARC：區(qū)域分解方法，提供圖形界面，用戶可以進(jìn)行手工剖分。操作不方便，而且很難達(dá)到一個計(jì)算任務(wù)的均衡。ANSYS：多重網(wǎng)格（Multigrid）求解器，主要用于流體動力學(xué)方程。NASTRAN和ABUQUS：具體算法不詳細(xì)。4 大規(guī)模計(jì)算技術(shù)三維有限元并行計(jì)算面臨的困難

35、 使用上的考慮 使用上應(yīng)該方便 計(jì)算任務(wù)的分配 如何將計(jì)算任務(wù)分配到每個計(jì)算節(jié)點(diǎn)上，并保持計(jì)算任務(wù)的平衡？ 分布內(nèi)存并行機(jī)上的并行實(shí)現(xiàn) 如何解決內(nèi)存的使用問題？當(dāng)規(guī)模增大時，怎樣將剛度矩陣保存在每個計(jì)算節(jié)點(diǎn)上？關(guān)鍵問題：存儲和計(jì)算任務(wù)在分布存儲并行機(jī)上的實(shí)現(xiàn)。4 大規(guī)模計(jì)算技術(shù)TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序簡介TFINE.Pfem （Parallel Finite Element Method）程序是采用EBE策略的有限元程序，使用C/C+語言開發(fā)。 TSINGHUA UNIVERSITY4 大規(guī)模計(jì)算技術(shù)TFINE.Pfem的計(jì)算流程 TFINE.Pfem（Para

36、llel Finite Element Method）是基于TFINE三維有限元計(jì)算程序，采用EBE策略的有限元程序，使用C/C+語言開發(fā)而成 TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序簡介運(yùn)行環(huán)境TFINE.Pfem程序及其驗(yàn)證程序驗(yàn)證矩形截面懸臂梁 TFINE.Pfem程序及其驗(yàn)證程序驗(yàn)證厚壁圓筒 內(nèi)壓應(yīng)力分布 規(guī)律TFINE.Pfem程序及其驗(yàn)證程序驗(yàn)證厚壁圓筒 內(nèi)壓位移分布圖 （a）理論值和計(jì)算值沿徑向分布曲線，（b） 整體位移等值線圖 （a）（b）TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率普通方法EBE方法EBE方法和常規(guī)方法的迭代次數(shù)對比TFI

37、NE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率EBE方法和常規(guī)方法的計(jì)算時間對比與常規(guī)方法相比，EBE方法的迭代次數(shù)減少，同時計(jì)算量也減少TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率EBE方法和常規(guī)方法的內(nèi)存需求對比均采用8節(jié)點(diǎn)6面體單元：EBE方法：12520300（單元剛度矩陣只存儲上三角 ）二維等帶寬存儲：1475332200EBE方法只占二維等帶寬存儲的4左右單元總數(shù)：12520節(jié)點(diǎn)總數(shù)：14730 700m700m1000m294m255m305m330mTFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率并行效率分析矩形截面懸臂梁T

38、FINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率并行效率分析矩形截面懸臂梁（）TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率并行效率分析二灘拱壩地基系統(tǒng)的有限元數(shù)值模擬有限元網(wǎng)格：20879個節(jié)點(diǎn)17980個單元單元類型：8節(jié)點(diǎn)六面體單元6節(jié)點(diǎn)五面體單元壩體和地基共包含31種材料計(jì)算模擬了施工的過程：（1）使用地基自重來模擬地應(yīng)力（2）計(jì)算壩體自重情況下的應(yīng)力和位移（3）計(jì)算水載作用下的應(yīng)力和位移二灘的有限元計(jì)算網(wǎng)格（其中不同的顏色表示不同的材料） TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率并行效率分析二灘拱壩地基系統(tǒng)的有限元數(shù)值模擬

39、單機(jī)計(jì)算效率與FLAC-3D的比較（單位：s）TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率并行效率分析二灘拱壩地基系統(tǒng)的有限元數(shù)值模擬并行計(jì)算效率采用了116個CPU進(jìn)行計(jì)算 TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的計(jì)算效率并行效率分析二灘拱壩地基系統(tǒng)的有限元數(shù)值模擬EBE方法的計(jì)算任務(wù)分配的本質(zhì)問題區(qū)域分解TFINE.Pfem程序及其驗(yàn)證TFINE.Pfem程序的開裂計(jì)算數(shù)值算例四點(diǎn)彎曲梁 錦屏邊坡的穩(wěn)定和加固分析 （自由度150萬，使用16個CPU）錦屏拱壩的大規(guī)模數(shù)值計(jì)算 （自由度107萬，使用16個CPU） 水布埡地下洞室的并行數(shù)值計(jì)算 （自由度

40、70萬，使用16個CPU） 二灘拱壩的開裂分析 （自由度47萬，使用8個CPU）工程算例TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究工程概況 錦屏一級水電站位于四川省涼山彝族自治州鹽源縣和木里縣境內(nèi)，是雅礱江干流上的重要梯級電站，壩址位于鹽源縣洼里鄉(xiāng)下游普斯羅溝與手爬溝間1.5km長的河段上。攔河大壩為混凝土雙曲拱壩，壩高305m，水庫庫容77.6億m3，電站裝機(jī)容量3600MW，年發(fā)電量166.2億KW.h。 錦屏一級水電站工程規(guī)模巨大，開挖邊坡最高達(dá)200300m。壩址區(qū)斷層、層間擠壓錯動帶、節(jié)理裂隙等構(gòu)造結(jié)構(gòu)面發(fā)育，特別是壩區(qū)左岸坡體內(nèi)存在的深部裂縫，

41、走向與岸坡近于平行，陡傾坡外，松弛張開顯著，對邊坡穩(wěn)定不利。因此，詳細(xì)查明工程邊坡開挖區(qū)的工程地質(zhì)條件，研究邊坡巖體可能出現(xiàn)的破壞形式及其在地震、霧化等工況下的穩(wěn)定性，對合理選取開挖坡比和支護(hù)方案，確保施工期和運(yùn)行期工程安全具有十分重要的意義。 TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究錦屏山TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究錦屏一級水電站高邊坡穩(wěn)定和加

42、固措施研究樞紐區(qū)左岸高邊坡全貌1960m2402m2038m2146m2284m2230m2105m2056m2100m1970mf8f53(4)厚層砂巖夾板巖1961m3(6)厚層砂巖厚60-80m1930m3(5)板巖厚45-50m X2120m2050m2080m2090m2020m3(4)1970m開口線處巖性纜機(jī)平臺1960m壩頂高程1885m拱肩槽開口線PD38PD48PD42PD64PD18PD54PD66PD40PD12PD30PD36PD10順坡裂隙延伸長、連通程度較高。順坡裂隙延伸短、連通程度較低。梁下游側(cè)強(qiáng)卸荷帶表層順坡裂隙的垂直分帶性NNE組順坡裂隙的發(fā)育特征： 在強(qiáng)卸

43、荷表層順坡裂隙延伸較長、貫通程度高，使坡體結(jié)構(gòu)呈現(xiàn)順向坡；以里延伸短小、貫通程度低； 在 梁，順坡裂傾角一般較陡(65-75 )，少數(shù)較緩(45-50 )；在 梁還隨高程表現(xiàn)出一定的垂直分帶性(見右圖)。NNE組順坡裂隙3(6)厚層砂巖需清理的危石3(5)粉砂質(zhì)板巖危巖體N60W向(卸荷)陡裂隙層面裂隙PD33主洞口9#公路順傾邊坡中的楔形體滑動TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究危巖體小的危巖體需要對地形地貌進(jìn)行精細(xì)的模擬，然后采用數(shù)值方法進(jìn)行計(jì)算。TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究錦屏一級水電站高邊坡

44、穩(wěn)定和加固分析錦屏一級水電站的開挖邊坡最高達(dá)200300m。壩址區(qū)斷層、層間擠壓錯動帶、節(jié)理裂隙等構(gòu)造結(jié)構(gòu)面發(fā)育，特別是左岸坡體內(nèi)存在的深部裂縫，走向與岸坡近于平行，陡傾坡外，對邊坡穩(wěn)定不利。 139536個單元149585個節(jié)點(diǎn)計(jì)算采用Drucker-Prager屈服準(zhǔn)則，模擬了整個開挖過程，從纜機(jī)平臺開挖一直模擬到拱肩槽完全開挖完成，共包含17個開挖步。采用擬靜力法考慮橫河向地震荷載，峰值加速度為0.1g。 左岸開挖步驟16個開挖步Gravity（地應(yīng)力）2080206020402020200019801960（纜機(jī)平臺）1885（壩頂高程）187018301790175017101670

45、163016001580 計(jì)算模型和計(jì)算工況 變形拉裂體與煌斑巖脈和f5斷層 計(jì)算模型和計(jì)算工況 f42-9上的抗剪洞 變形拉裂體 TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究屈服準(zhǔn)則巖體強(qiáng)度采用Drucker-Prager屈服準(zhǔn)則（D-P準(zhǔn)則）和抗拉屈服準(zhǔn)則（混凝土和巖體為低抗拉材料，故應(yīng)有抗拉準(zhǔn)則） TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究屈服準(zhǔn)則與剛體極限平衡法不同，在三維非線性有限元穩(wěn)定分析中，穩(wěn)定性表現(xiàn)為彈塑性分析解的收斂性。高拱壩及壩肩的特點(diǎn)宜于造成三維畸形網(wǎng)格、高度不均勻的材料分區(qū)、高量級荷載，這都對收斂

46、性有很不利的影響。TFINE程序采用了基于D-P準(zhǔn)則的迭代算法和積分策略，提高了數(shù)值穩(wěn)定性。該方法可有效改善計(jì)算的收斂性，并確保計(jì)算收斂到正確解上。并在此基礎(chǔ)上發(fā)展出一種基于非線性規(guī)劃的有限元極限分析方法，可大幅度減少優(yōu)化變量，克服“維數(shù)障礙”。該算法已被集成在TFINE及TFINE.Pfem程序內(nèi)，并被用于研究拱端力在超載過程的非線性分布規(guī)律。 TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究地應(yīng)力回歸 三維初始地應(yīng)力場的多元回歸分析方法：線性疊加 回歸結(jié)果滿足屈服條件的非線性校正 對各因素影響按線彈性計(jì)算，經(jīng)線性疊加得到的回歸分析結(jié)果，并不能保證回歸所得到的

47、初始地應(yīng)力場內(nèi)各點(diǎn)的應(yīng)力狀態(tài)都能滿足塑性屈服條件。這時就需要對回歸結(jié)果進(jìn)行校正，讓初始地應(yīng)力場內(nèi)各點(diǎn)的應(yīng)力狀態(tài)都滿足塑性屈服條件，從而使初始地應(yīng)力場的分析結(jié)果更加合理。 回歸結(jié)果滿足平衡條件的非線性校正 有限元分析計(jì)算求得的應(yīng)力場不但滿足整體平衡，而且在內(nèi)部節(jié)點(diǎn)也滿足平衡條件。 材料參數(shù)的確定參數(shù)的確定完全依據(jù)成勘院可行性研究報告的工程地質(zhì)報告。根據(jù)錦屏一級壩基巖體質(zhì)量分級表和壩區(qū)結(jié)構(gòu)面按性狀分類表進(jìn)行確定（1）各類巖體參數(shù)的確定各類巖體參數(shù)的確定采用加權(quán)平均的方法，主要根據(jù)巖體質(zhì)量分級表中提供數(shù)據(jù)確定。 巖體變模的確定 將平行于結(jié)構(gòu)面和垂直于結(jié)構(gòu)面變模的數(shù)值加權(quán)后取平均值。 f和c的確定 采

48、用壩基巖體質(zhì)量分級表提供的數(shù)據(jù)。 容重的確定 主要參考壩區(qū)新鮮巖塊物理力學(xué)性質(zhì)匯總表。（2）各類結(jié)構(gòu)面參數(shù)的確定各類結(jié)構(gòu)面參數(shù)依據(jù)壩區(qū)結(jié)構(gòu)面按性狀分類表中代表性結(jié)構(gòu)面項(xiàng)和壩基巖體質(zhì)量分級表中的巖級分布位置項(xiàng)的描述，從而判斷出結(jié)構(gòu)面所屬的巖級種類，最終確定其力學(xué)參數(shù)。 材料參數(shù)的確定計(jì)算工況工況1：在自重應(yīng)力場情況下，模擬開挖過程；工況2：地基材料參數(shù)f和C降低1.3倍，進(jìn)行開挖過程的模擬；工況3：在回歸地應(yīng)力場情況下，模擬開挖過程；工況4：完全開挖后，按擬靜力法施加橫河向的動力荷載；工況5：完全開挖后，模擬滲流場施加滲壓荷載；工況6：考慮錨索加固（錨索長度為70m）情況下的開挖過程 模擬；工況

49、7：考慮錨固洞加固影響的開挖過程模擬；工況8：變形體f5外側(cè)鎖固段巖體完全挖除的敏感性分析。計(jì)算荷載（1）自重：按照地質(zhì)報告提供的參數(shù)進(jìn)行計(jì)算；（2）滲壓：上游水位為正常蓄水位1880.0m，下游水位取1640.0m，中間壩體部位考慮防滲帷幕的影響，水頭采用線性過渡；（3）動力荷載：設(shè)計(jì)地震烈度7度，只考慮橫河向加速度0.1g，按照擬靜力法施加動力荷載。 TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究錦屏一級水電站高邊坡穩(wěn)定和加固分析左岸邊坡穩(wěn)定的關(guān)鍵塊體有兩個： 滑移塊1 左岸的NNW向裂密帶煌斑巖脈f42-9斷層 滑移塊2 左岸的NNW向裂密帶卸荷裂隙f42-9斷層 TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究錦屏一級水電站高邊坡穩(wěn)定和加固分析網(wǎng)格密度對結(jié)果的精度影響分析 不加密網(wǎng)格（單元數(shù)目為17442）加密網(wǎng)格（單元數(shù)目為139536） TSINGHUA UNIVERSITY5 錦屏一級水電站高邊坡穩(wěn)定和加固措施研究錦屏