石壩現(xiàn)狀的分析管理論文

1、石壩現(xiàn)狀的分析管理論文 現(xiàn)代碾壓面板堆石壩具有工程量小工期短造價(jià)低等優(yōu)點(diǎn)，已成為壩型選擇的主要類型之一，在我國(guó)已成大規(guī)模推廣的趨勢(shì)。國(guó)內(nèi)已建成的有關(guān)門山、柯柯亞、成屏、西北口、天生橋面板堆石壩等十余座，正在規(guī)劃設(shè)計(jì)的有10余座，廣泛分布于全國(guó)各地。20世紀(jì)50年代以前設(shè)計(jì)的堆石壩一般是不計(jì)算壩體應(yīng)力應(yīng)變的，但隨著計(jì)算技術(shù)的發(fā)展和壩體的日益增高，同時(shí)在一些壩體中發(fā)現(xiàn)裂縫，壩體應(yīng)力應(yīng)變計(jì)算才逐步為人們所重視。 從“七五”科技攻關(guān)開始，國(guó)內(nèi)專家用不同的計(jì)算模型，考慮面板受力特點(diǎn)，采用相應(yīng)的參數(shù)，對(duì)高面板壩進(jìn)行二維、三維有限元分析計(jì)算。但是到目前為止，國(guó)內(nèi)外對(duì)這種壩主要憑經(jīng)驗(yàn)設(shè)計(jì)。由于堆石體材料特性十

2、分復(fù)雜，有限元的計(jì)算結(jié)果與實(shí)際觀測(cè)值之間還存在一定的差距。可見，在這方面尚有許多研究工作有待開展。 堆石體的本構(gòu)關(guān)系表達(dá)了堆石體應(yīng)力應(yīng)變之間的關(guān)系，它無疑對(duì)應(yīng)力應(yīng)變的計(jì)算結(jié)果起決定性作用。目前建立堆石體本構(gòu)關(guān)系往往基于已有模型，再針對(duì)堆石的力學(xué)特性確定甚至調(diào)整本構(gòu)關(guān)系中各種材料參數(shù)。面板壩有限元計(jì)算時(shí)常采用非線性彈性模型和彈塑性模型。 2.1非線性彈性模型 非線性彈性模型包括鄧肯模型、內(nèi)勒模型、修正鄧肯模型和鄧肯模型等。這些模型最初是針對(duì)土、砂等一類材料在常規(guī)三軸試驗(yàn)基礎(chǔ)上提出的，僅適用于二維分析計(jì)算，其中計(jì)算結(jié)果和已建壩實(shí)測(cè)結(jié)果較為符合的是鄧肯E-B模型和內(nèi)勒模型7。 2.1.1鄧肯模型18

3、 康納(Kondner)根據(jù)常規(guī)三軸試驗(yàn)結(jié)果發(fā)現(xiàn)，關(guān)系可用雙曲線擬合。隨后，鄧肯(Duncan)引入摩爾-庫(kù)侖(MohrCoulomb)破壞準(zhǔn)則，導(dǎo)出切線彈模為 式中,為凝聚力；為初始彈模；為破壞比,是破壞時(shí)主應(yīng)力差和應(yīng)力應(yīng)變雙曲線的漸近線的縱坐標(biāo)的比值,即/；為內(nèi)摩擦角,對(duì)于粗粒土假定是的函數(shù)，即,式中和為試驗(yàn)參數(shù)，大氣壓(下同)。 詹布(Janbu)據(jù)試驗(yàn)指出初始彈模和側(cè)限壓力有如下關(guān)系：式中、為兩個(gè)由試驗(yàn)確定的參數(shù)。 相應(yīng)于卸荷再加荷情況的切線彈模為，按下式計(jì)算式中、是由試驗(yàn)確定的兩個(gè)參數(shù)，一般較為大。 1980年，鄧肯(Duncan)提出切線體積模量的計(jì)算式式中、為試驗(yàn)常數(shù)。 體模相當(dāng)

4、于修正鄧肯模型中的,考慮到只能在之間變化，故有限元計(jì)算中應(yīng)限定在范圍內(nèi)，否則得不到合理計(jì)算結(jié)果。 平面應(yīng)變條件下鄧肯模型的增量型本構(gòu)關(guān)系為式中、為法向應(yīng)力增量；為剪應(yīng)力增量；、為法向應(yīng)變?cè)隽?；為剪?yīng)變?cè)隽俊?鄧肯模型的不足之處在于它不能反映壓縮與剪切的交叉影響；不能反映各向異性；不能反映加荷卸荷對(duì)的變化；不能反映中主應(yīng)力對(duì)強(qiáng)度指標(biāo)的影響等。但總的來說，該模型反映了堆石體變形的主要規(guī)律。它反映了非線性；用于增量計(jì)算時(shí)，能反映應(yīng)力路徑對(duì)變形的影響；通過與的差別部分體現(xiàn)了加荷卸荷對(duì)變形的影響；通過假定是的函數(shù)，在一定程度上考慮了高固結(jié)壓力的影響。再加上該模型簡(jiǎn)便直觀，概念明確，并積累了相當(dāng)多的應(yīng)用經(jīng)

5、驗(yàn)，因此目前被工程設(shè)計(jì)人員廣泛使用。 2.1.2內(nèi)勒模型7 內(nèi)勒非線性彈性模型是模型中較為簡(jiǎn)單的一種，模型參數(shù)容易由常規(guī)三軸試驗(yàn)確定。內(nèi)勒(Naylor)認(rèn)為土的切線體積變形模量隨側(cè)限壓力的增加而增加，土的切線剪切模量隨剪應(yīng)力的增加而減小，在破壞時(shí)變?yōu)榱恪Ｒ虼饲芯€體積變形模量為切線剪切模量為內(nèi)勒非線性彈性模型在面板壩有限元分析中較少采用。 2.2彈塑性模型4 彈塑性模型把總的變形分成彈性變形和塑性變形兩部分，前者用虎克定律計(jì)算，后者用塑性理論求解。南京水科院1987年提出了用于堆石體的雙屈服面彈塑性模型，一般稱為南水模型。該模型可以反映應(yīng)力引起的各向異性和堆石的剪縮特性，在理論方面有其優(yōu)越性。

6、但限于試驗(yàn)設(shè)備，該模型中某些參數(shù)如、較難從試驗(yàn)中得到,因此尚未進(jìn)入工程實(shí)用階段。 3.非線性有限元求解方法16 一般來說，求解方法有迭代法和增量法兩種，較為常用的是增量法。用增量法計(jì)算時(shí)，荷載逐級(jí)遞增，可以模擬施工過程，計(jì)算結(jié)果也符合觀測(cè)結(jié)果。以鄧肯模型為例，假設(shè)現(xiàn)進(jìn)行第級(jí)加荷計(jì)算，按中點(diǎn)增量法其計(jì)算步驟如下： 1根據(jù)前一級(jí)應(yīng)力全量確定彈性常數(shù)和，并形成剛度矩陣； 2加本級(jí)荷載增量的一半于結(jié)構(gòu)，用下式求位移增量： 并計(jì)算應(yīng)力和應(yīng)變?cè)隽?，進(jìn)而累計(jì)求得應(yīng)力全量，據(jù)此計(jì)算； 3施加全荷，求位移增量： 相應(yīng)地可求出應(yīng)力應(yīng)變?cè)隽?，累?jì)則可得位移和應(yīng)力應(yīng)變?nèi)浚?重復(fù)上述步驟，可得各級(jí)荷載增量下的解答。

7、還有一個(gè)問題，即新填筑層各單元初始應(yīng)力狀態(tài)是，如果以此計(jì)算，則切線彈模，就無法進(jìn)行計(jì)算?？死?Clough)等人將新填土層視為重液體處理，即(為單元形心距土表的距離，為填土容重)。目前較為常用的處理方法式計(jì)算是引入側(cè)壓力系數(shù)，即。事實(shí)上，新填土層經(jīng)過反復(fù)碾壓，可視為超固結(jié)土，可取前期預(yù)固壓力(視碾壓輕重可取MPa)。 4算例 天生橋一級(jí)面板堆石壩位于紅水河上游南盤江干流，壩高178m，上游壩坡 11.4,下游壩坡11.29。按平面應(yīng)變問題計(jì)算。 4.1單元剖分和計(jì)算參數(shù) 面板頂厚0.3，底厚為0.9m，沿面板厚度方向剖分一層單元，沿壩高分20層，共剖分500個(gè)單元，1486個(gè)節(jié)點(diǎn)。其中堆石單

8、元459個(gè)，面板單元20個(gè)，接觸面單元20個(gè)，趾板單元1個(gè)。單元形態(tài)為等參元計(jì)算采用中點(diǎn)增量法，按20級(jí)加荷方式進(jìn)行，每級(jí)荷載又分三次施加。這樣相當(dāng)于按60級(jí)加荷計(jì)算，但又比剖分60層的情況節(jié)省計(jì)算時(shí)間。 混凝土面板和趾板采用線彈性模型，=18000MPa，=0.167。 堆石各分區(qū)對(duì)應(yīng)于模型的計(jì)算參數(shù)見表1。 4.2計(jì)算成果 從面板壩有限元的計(jì)算結(jié)果來看，依據(jù)不同模型所得的應(yīng)力結(jié)果是比較一致的，不同之處在于位移的大小和分布。而實(shí)際工程中，設(shè)計(jì)人員也更關(guān)心位移的大小和分布。 圖2和圖3分別為竣工時(shí)鉛直和水平位移等值線圖。最大鉛直位移為2.35m，約位于1/2壩高處；最大水平位移：上游為0.58

9、m，下游為1.03m，參見圖2。 5.計(jì)算結(jié)果和實(shí)測(cè)值的比較分析 除本文算例的計(jì)算結(jié)果外，還將引用文獻(xiàn)3和文獻(xiàn)4對(duì)天生橋一級(jí)面板壩的計(jì)算結(jié)果，來主要說明鄧肯模型和彈塑性模型在計(jì)算結(jié)果上的差異。 天生橋一級(jí)面板壩壩體變形觀測(cè)，沉降觀測(cè)結(jié)果規(guī)律性強(qiáng)，十分可信；水平位移的觀測(cè)結(jié)果，國(guó)內(nèi)外專家尚有不同看法。因此，下文主要針對(duì)鉛直位移進(jìn)行討論。 表2天生橋一級(jí)面板壩有限元位移計(jì)算結(jié)果 竣工時(shí)最大鉛直位移的觀測(cè)值為3.06m，約為壩高的1.72%，是同類型壩中沉降量較大的。從數(shù)值上看，鄧肯模型計(jì)算的最大鉛值位移在2.3m左右，和南水模型的計(jì)算結(jié)果1.34m相比，和實(shí)測(cè)值較為接近，但仍相差0.76m，究其原

10、因，一是計(jì)算不能考慮實(shí)際施工狀況對(duì)沉降的影響，二是按鄧肯模型計(jì)算時(shí)沒有考慮堆石體流變引起的沉降。 南水模型由于可以考慮堆石體的剪縮特性，計(jì)算的鉛值和水平位移都較小。 雖然無可靠的實(shí)測(cè)水平位移以資比較，一般認(rèn)為依據(jù)南水模型計(jì)算的水平位移較為可靠,而鄧肯模型計(jì)算的水平位移一般偏大。另外，依據(jù)南水模型計(jì)算的位移的分布規(guī)律較模型更符合實(shí)際觀測(cè)34。 雖然設(shè)計(jì)面板壩時(shí)更關(guān)心面板、周邊縫、伸縮縫，止水的變位與位移，但基本變位是堆石體的變位，它對(duì)其他構(gòu)件的變位有重大的影響。因此以堆石體變位來評(píng)價(jià)模型應(yīng)該是較為合適的。通過以上的分析可知，兩個(gè)模型都存在不足之處，面板壩有限元分析還有待深入研究。 6面板壩有限元

11、分析的展望2 a.今后一段時(shí)間內(nèi)，簡(jiǎn)單實(shí)用、概念明確的模型仍占主導(dǎo)地位。而計(jì)算結(jié)果在多大程度上和實(shí)際相符，不僅取決于模型本身，也取決于模型參數(shù)測(cè)定的準(zhǔn)確程度。以模型為例，該模型涉及七個(gè)參數(shù)，任何參數(shù)偏離都會(huì)造成計(jì)算結(jié)果的變動(dòng)。鑒于此，有必要在參數(shù)測(cè)定方面作深入研究，研制新的大型三軸試驗(yàn)儀，發(fā)展原位試驗(yàn)技術(shù)等。 b.應(yīng)該承認(rèn)，現(xiàn)存的本構(gòu)關(guān)系在描述堆石體真實(shí)特性的準(zhǔn)確性和完整性方面是遠(yuǎn)遠(yuǎn)不夠。無論是非線性模型還是彈塑性模型，都是在連續(xù)性假設(shè)基礎(chǔ)上從宏觀的唯象的角度描述堆石材料的應(yīng)力應(yīng)變關(guān)系。實(shí)際上，堆石體可被視為多種材料組成的多孔介質(zhì)，不連續(xù)性為其主要特征。因此，要深刻揭示堆石體的應(yīng)力應(yīng)變規(guī)律，就

12、必然要開展堆石材料的微觀研究，并力圖在微觀研究的基礎(chǔ)上，結(jié)合新理論，提出新模型。 c.面板壩從勘測(cè)、設(shè)計(jì)、試驗(yàn)研究到施工的各個(gè)階段無不存在許多不確定性因素，因此，計(jì)算結(jié)果和觀測(cè)值完全相符幾乎是不可能的。但是如果在計(jì)算中考慮到不確定因素的影響，則會(huì)有效改善計(jì)算結(jié)果?？紤]不確定因素基本上有兩種方法：一是采用隨機(jī)有限元，一時(shí)采用模糊有限元。對(duì)于壩工而言，要準(zhǔn)確得到不確定量的隨機(jī)分布是十分困難的。因此，建議將不確定量按模糊量處理，采用模糊有限元作分析計(jì)算，計(jì)算結(jié)果仍為模糊量，這也符合一般的工程經(jīng)驗(yàn)。 d.采用新的數(shù)值計(jì)算方法如離散元法、流形元法。 1錢家歡、殷宗澤，土工原理與計(jì)算，北京：中國(guó)水利水電出版社，1996 2何廣訥，土工的若干新理論研究與應(yīng)用，北京：水利電力出版社，1994 3混凝土面板堆石壩會(huì)議論文集，河海大學(xué)出版社，1990 4章為民、沈