地下洞室群面向?qū)ο笥邢拊椒ǖ奶接?/h1>
2016-01-13 06:45:04周永強(qiáng)付曉東劉曉敏
長江科學(xué)院院報(bào) 2015年7期
周永強(qiáng)，盛 謙，付曉東，劉曉敏
(中國科學(xué)院 a.武漢巖土力學(xué)研究所;b.巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室， 武漢 430071)
地下洞室群面向?qū)ο笥邢拊椒ǖ奶接?/p>
周永強(qiáng)a，b，盛謙a，b，付曉東a，b，劉曉敏a，b
(中國科學(xué)院 a.武漢巖土力學(xué)研究所;b.巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室， 武漢430071)
摘要：以單元為代表類簡單介紹了地下洞室群面向?qū)ο笥邢拊绦?。利用自編的面向?qū)ο笥邢拊绦驅(qū)崿F(xiàn)了開挖模擬，并分析某地下廠房的分步開挖變形問題，其結(jié)果與Phase2軟件計(jì)算結(jié)果基本一致；實(shí)現(xiàn)了面向?qū)ο笥邢拊幚韯?dòng)力響應(yīng)的人工邊界輸入問題，驗(yàn)證了其正確性，并對該地下廠房的動(dòng)力響應(yīng)規(guī)律進(jìn)行了探討。結(jié)果表明：監(jiān)測點(diǎn)響應(yīng)存在一定的滯后性，且峰值的放大系數(shù)比無洞室作用下的理論解要大；在此基礎(chǔ)上，分析了不同方法下的阻尼系數(shù)對地下洞室群動(dòng)力響應(yīng)的影響；驗(yàn)證了面向?qū)ο笥邢拊椒☉?yīng)用在地下洞室群中的可行性和適用性。研究成果為進(jìn)一步研究動(dòng)態(tài)彈塑性情況奠定了基礎(chǔ)。
關(guān)鍵詞：地下洞室群；面向?qū)ο螅挥邢拊绦?；人工邊界；阻?/p>
1研究背景
傳統(tǒng)的有限元程序一般采用結(jié)構(gòu)化設(shè)計(jì)，即采用例如Fortran編程語言的面向過程進(jìn)行編寫，然而隨著工程問題復(fù)雜性的增多和計(jì)算規(guī)模的增大，面向過程的編程語言由于其自身的特點(diǎn)，如代碼重新利用率低、調(diào)試復(fù)雜、開發(fā)周期長等，潛在的問題日益突出，因而已經(jīng)難以滿足開發(fā)和維護(hù)上的要求。
采用面向?qū)ο缶幊陶Z言的有限元程序可以彌補(bǔ)上述不足，可改善有限元程序的性能，提高有限元程序的開發(fā)效率，并能為功能模塊的集成化、擴(kuò)充提供有利條件。面向?qū)ο蟮某绦蛟O(shè)計(jì)方法與有限元方法相結(jié)合，并沒有從根本上改變有限元方法的核心思想，而是程序設(shè)計(jì)和開發(fā)軟件的一種新嘗試。
Forde等[2]在1990年首次提出了面向?qū)ο笥邢拊绦蛟O(shè)計(jì)，之后國內(nèi)外對該設(shè)計(jì)方法做了許多研究。Zimmermann等[3]提出了面向?qū)ο笥邢拊绦蛟O(shè)計(jì)控制原理，Miller等[4]提出了基于面向?qū)ο笥邢拊椒ǖ姆蔷€性動(dòng)力學(xué)分析的軟件框架，李于峰等[5]實(shí)現(xiàn)了面向?qū)ο笥邢拊⑿杏?jì)算框架PANDA，Werner等[6]研究了利用面向?qū)ο蟮哪Ｐ蛠碓O(shè)計(jì)和分析隧道工程問題，平揚(yáng)等[7]實(shí)現(xiàn)了面向?qū)ο笥邢拊谏罨庸こ讨械膽?yīng)用，楊志勇等利用面向?qū)ο笥邢拊椒ǚ治隽私Y(jié)構(gòu)抗震問題，項(xiàng)陽等[8]則提出了面向?qū)ο笥邢拊椒ㄔ趲r土工程中的應(yīng)用實(shí)例，劉小生等[9]應(yīng)用面向?qū)ο笥邢拊椒ń鉀Q大壩滲流問題，簡榮林等[10]利用C++語言開發(fā)了基于Windows平臺的巷道圍巖變形有限元分析軟件。
從上述可以看出，面向?qū)ο笥邢拊讶〉昧讼喈?dāng)大的進(jìn)展，然而該方法應(yīng)用到地下洞室群領(lǐng)域的研究少有報(bào)道。本文首先以單元類作為代表介紹地下洞室群面向?qū)ο笥邢拊绦颍黄浯卫米跃幊绦騺矸治瞿车叵聫S房的開挖變形問題；最后對該地下廠房的動(dòng)力響應(yīng)規(guī)律進(jìn)行了探討；以此為基礎(chǔ)，分析了不同方法下的阻尼系數(shù)對地下洞室群動(dòng)力響應(yīng)的影響。
2面向?qū)ο笥邢拊绦蛟O(shè)計(jì)
有限元的基本思想是先將整個(gè)計(jì)算區(qū)域離散成許多單元，然后在各個(gè)單元上計(jì)算，最后整合在一起，因此單元反映了整個(gè)有限元離散體系的拓?fù)鋵傩?。本?jié)將重點(diǎn)介紹單元類，以此為基礎(chǔ)簡單概括地下洞室群面向?qū)ο笥邢拊绦?。單元類中既包含了材料屬性、單元?jié)點(diǎn)、單元類型、單元編號等基本數(shù)據(jù)屬性，同時(shí)又具備一些計(jì)算功能，如外荷載、單元的應(yīng)力和應(yīng)變、質(zhì)量矩陣和阻尼矩陣等。該程序是用VC++語言來編寫的，單元類的程序如圖1所示。
圖1 單元類程序Fig.1 Class element
其中ComputeNewmarkLhsAt(TimeStep*)和ComputeNewmarkRhsAt(TimeStep*)分別為Newmark法計(jì)算動(dòng)力方程左右兩邊的函數(shù)；giveDampMatrix(int, FloatMatrix*, TimeStep*)是計(jì)算阻尼矩陣的函數(shù)；giveStep()則為計(jì)算開挖步的函數(shù)；而FloatMatrix，F(xiàn)loatArray和IntArray分別是矩陣類、向量類和整數(shù)向量類；GaussPoint為高斯節(jié)點(diǎn)類；TimeStep則為時(shí)間步類，一般用于動(dòng)力計(jì)算。從圖1可以看出，地下洞室群面向?qū)ο笥邢拊绦蚣扔幸话阌邢拊绦虻墓残?，又有其自身的特點(diǎn)，如開挖、分析動(dòng)力響應(yīng)以及考慮阻尼等。
已知t和t+Δt時(shí)刻的動(dòng)力方程分別為
(1)
(2)
假設(shè)Δut=ut+Δt-ut，則根據(jù)Newmark法有
(3)
結(jié)合式(3),由式(2)-式(1)可得
(4)
(5)
式中:K,M和C分別為剛度矩陣、質(zhì)量矩陣和瑞利阻尼矩陣;α和δ是Newmark法參數(shù)；γ和β是不依賴頻率的參數(shù)；ρ為材料密度。對于解的穩(wěn)定性，Newmark法所選用的Δt的值是任意的，但考慮到精度要求，Δt則需小于體系固有振型最小周期的(1/10～1/20)。根據(jù)式(4)和式(5)可知，單元類的組裝過程如圖2所示。
圖2 單元類的組裝過程Fig.2 Assembly process of class element
3地下洞室群的開挖模擬
3.1基本思路
模擬地下洞室的開挖過程，重點(diǎn)是把作用在需要開挖部分的重力去除(不考慮構(gòu)造應(yīng)力)。地下洞室群在開挖之前，重力荷載是利用質(zhì)量矩陣轉(zhuǎn)化成等效節(jié)點(diǎn)作用力來實(shí)現(xiàn)的，因此為模擬地下洞室群開挖過程中重力作用的變化，只需根據(jù)時(shí)間步將質(zhì)量矩陣進(jìn)行修改，而材料密度是影響質(zhì)量矩陣的唯一因素。為了防止開挖的單元不起作用而使總體剛度矩陣產(chǎn)生奇異，開挖單元的彈性模量不能取得太大，也不能取得太小，一般比周圍巖土體彈性模量小5至8個(gè)數(shù)量級[11]，具體思路如圖3所示。
圖3 地下洞室群開挖過程的思路Fig.3 Idea of the excavation process ofunderground caverns
圖4 半無限空間體簡化圖Fig.4 Simplified diagramof half infinite space body
為驗(yàn)證該程序的開挖思想的正確性，設(shè)有半無限空間體，如圖4所示，圖中數(shù)字表示單元節(jié)點(diǎn)。材料的彈性模量為1.2×104MPa，泊松比為0.2，密度為2.7×103kg/m3。每個(gè)單元的尺寸大小均為6.25 m×25 m。初始時(shí)該體系不開挖，之后開挖最上面的單元，則開挖后單元節(jié)點(diǎn)的位移為[11]
(6)
式中:u0為初始位移，Δu為變形回彈位移?？煞謩e按下式計(jì)算：
(7)
(8)
式中：Δh為開挖的豎向高度;υ為材料的泊松比。運(yùn)用自編程序?qū)w系進(jìn)行開挖模擬，其結(jié)果如表1所示。從表1中看出，數(shù)值解與解析解基本一致。
表1 體系開挖步1后節(jié)點(diǎn)位移的數(shù)值解和解析解
3.2實(shí)例驗(yàn)證
某水電站地下廠房某斷面的計(jì)算剖面如圖5所示，其主廠房、主變室和尾調(diào)室的尺寸大小分別為29 m×57 m，20 m×37 m，26 m×87 m，整個(gè)模型尺寸為1 500 m×600 m。材料的彈性模量為1.2×104MPa，泊松比為0.2，密度為2.7×103kg/m3。
圖5 地下洞室群某計(jì)算斷面Fig.5 Cross-section of underground cavernsfor computation
圖6 地下洞室群開挖工序及監(jiān)測點(diǎn)布置Fig.6 Excavation process of underground cavernsand layout of monitoring points
圖7 地下洞室群主變室各測點(diǎn)開挖變形Fig.7 Deformation of monitoring points in maintransformer chamber of underground caverns
該地下廠房共分為10步開挖，開挖工序及廠房的監(jiān)測點(diǎn)的布置如圖6所示。利用自編程序來模擬該地下廠房的開挖過程。利用Phase2軟件建立相同的計(jì)算模型，材料參數(shù)保持不變，得到相應(yīng)監(jiān)測點(diǎn)的結(jié)果。圖7為地下洞室群主變室各監(jiān)測點(diǎn)的位移隨開挖步數(shù)的變化圖，限于篇幅，其他洞室就不一一給出。從圖7中可知，自編程序與Phase2軟件計(jì)算的結(jié)果基本一致，說明該程序的開挖基本思路是正確的，也說明了該程序能正確分析地下洞室群的開挖問題。
4地下洞室群的波動(dòng)影響
4.1人工邊界及地震輸入
要知道地下洞室群的波動(dòng)影響，首先要解決其地震輸入問題，即動(dòng)力荷載在邊界上的輸入。動(dòng)力荷載的輸入方法有2種：位移輸入和應(yīng)力輸入，本程序采用應(yīng)力輸入，即在輸入邊界上設(shè)置人工邊界，并同時(shí)把速度時(shí)程通過式(9)和式(10)轉(zhuǎn)換成應(yīng)力時(shí)程。
(9)
(10)
式中：vn和vs分別為輸入波在模型邊界上質(zhì)點(diǎn)的法向和切向的速度分量；ρ為材料密度；CP和CS分別為P波和S波在介質(zhì)中的傳播波速。
其中人工邊界的作用是吸收邊界上的入射波，減少因邊界上存在反射波而對動(dòng)力分析的結(jié)果產(chǎn)生影響。人工邊界大致可以分為局部人工邊界和全局人工邊界[12]，局部人工邊界由于其計(jì)算時(shí)間少、效率高等特點(diǎn)得到了廣泛的應(yīng)用。局部人工邊界一般包括黏性邊界、黏彈性邊界、旁軸邊界和透射邊界等4種邊界[13]，因黏性邊界簡單，易于編程實(shí)現(xiàn)，且效果較好，綜合考慮本程序采用了黏性邊界。為體現(xiàn)面向?qū)ο蟪绦虻亩鄳B(tài)性和繼承性，將含有黏性邊界的單元發(fā)展成黏性邊界單元，即其剛度矩陣和質(zhì)量矩陣保持不變，阻尼矩陣則為含有邊界特性阻尼矩陣，即根據(jù)虛功原理，采用類似質(zhì)量矩陣的集中形式(四邊形單元)，即
(11)
在分析動(dòng)力情況之前，往往先考慮體系的靜力平衡問題。劉晶波等[14-15]雖已將二維黏彈性邊界發(fā)展成黏彈性邊界單元，又推導(dǎo)了三維等效黏彈性邊界單元，并在此基礎(chǔ)上建立了三維黏彈性靜動(dòng)力統(tǒng)一人工邊界[16]，然而用此邊界來分析靜力問題時(shí)，會(huì)有一定的誤差[17]，且參數(shù)的選取有一定的主觀性。另外計(jì)算靜力問題(或動(dòng)力問題)，選取相對應(yīng)的靜力(或動(dòng)力)分析方法，然而靜力平衡之后分析動(dòng)力問題，若不去掉約束條件，同時(shí)又不施加約束反力，則整個(gè)體系是在不平衡狀態(tài)下計(jì)算動(dòng)力荷載，因此得出的結(jié)果會(huì)出現(xiàn)很大的波動(dòng)，與實(shí)際結(jié)果不相符。因此分析靜動(dòng)力問題時(shí)，應(yīng)選擇的思路為：以靜力分析結(jié)果的應(yīng)力場作為初始條件進(jìn)行動(dòng)力分析時(shí)，為保證體系在動(dòng)力計(jì)算初始時(shí)刻處于靜力平衡狀態(tài)，可在輸入靜應(yīng)力場的同時(shí)輸入引起此應(yīng)力場的荷載和約束，即在轉(zhuǎn)換邊界時(shí)，需去掉約束條件，同時(shí)施加相應(yīng)的約束反力。
為驗(yàn)證黏性邊界和地震波輸入方法在該程序中的可行性，建立了圖8所示的 200 m×2 000 m的計(jì)算模型，模型底部設(shè)置了固定約束，兩側(cè)設(shè)置了法向約束。材料參數(shù)：彈性模量為1.2×104MPa，泊松比為0.2，密度為2.7×103kg/m3。不考慮阻尼的影響，在模型左邊輸入一個(gè)周期的正弦速度波，頻率為1 Hz，峰值為1 m/s。在模型左右兩邊設(shè)置監(jiān)測點(diǎn)A和B，分別考慮在模型左邊設(shè)置黏性邊界和不設(shè)置黏性邊界2種情況，觀察監(jiān)測點(diǎn)的速度響應(yīng)時(shí)程，其結(jié)果如圖9所示。波從模型左端傳到右端，由于自由邊界反射后左行波疊加作用，使得監(jiān)測點(diǎn)B的峰值速度是入射波峰值速度的2倍，即放大了2倍，這一現(xiàn)象與理論解一致，且與文獻(xiàn)[18]得出的結(jié)果一致，同時(shí)黏性邊界較好地模擬了左行波穿過人工邊界向無限域輻射的過程。
圖9 各監(jiān)測點(diǎn)的速度時(shí)程曲線Fig.9 Curves of velocity vs. time at monitoring points
圖8動(dòng)力計(jì)算模型圖
Fig.8Model for dynamic computation
4.2波在地下洞室群的傳播規(guī)律
圖11 主變室各監(jiān)測點(diǎn)加速度時(shí)程峰值的放大倍數(shù)Fig.11 Peak amplification factor of accelerationprocess at monitoring points of maintransformer chamber
基于上述黏性邊界理論和地震波輸入方法的正確性，以上述水電站地下廠房某斷面(地下洞室群開挖完成之后)作為實(shí)例，各種材料參數(shù)和尺寸都保持不變，不考慮阻尼的影響，在模型底端輸入一個(gè)周期的正弦加速度波(頻率為1 Hz，峰值為1 m/s2)，并施加黏性邊界。為研究地下洞室群的動(dòng)力響應(yīng)，監(jiān)測布置點(diǎn)如圖5和圖6所示。同樣利用該自編程序計(jì)算，得到地下洞室群各監(jiān)測點(diǎn)的加速度響應(yīng)見圖10和圖11。由圖10可知，各監(jiān)測點(diǎn)動(dòng)力響應(yīng)存在滯后性，離模型底端越遠(yuǎn)，滯后的時(shí)間越長，如監(jiān)測點(diǎn)D滯后了約0.05 s，而監(jiān)測點(diǎn)B則大概滯后了0.22 s；由于洞室的影響，監(jiān)測點(diǎn)B的最大放大系數(shù)為2.08，大于2.0，監(jiān)測點(diǎn)C的最大放大倍數(shù)為1.21，也超過了1.0；0.57 s是監(jiān)測點(diǎn)A的波形時(shí)變化的轉(zhuǎn)折點(diǎn)，即0.57 s之前基本與輸入波形相同，之后則完全不一致，這是輸入波在洞室中產(chǎn)生反射波與輸入波疊加作用的結(jié)果。綜上分析，洞室的存在會(huì)增加地震波對圍巖變形的影響，因而需更加關(guān)注洞室頂拱和頂板的變形。
表2 瑞利阻尼系數(shù)的確定方法
5阻尼對地下洞室群動(dòng)力響應(yīng)的影響
阻尼的產(chǎn)生主要源于材料內(nèi)部摩擦以及可能存在的接觸表面的滑動(dòng)。在實(shí)際分析中，精確確定阻尼矩陣是相當(dāng)困難的，通常是將阻尼矩陣表示為質(zhì)量矩陣和剛度矩陣的線性組合(瑞利阻尼)，即
(12)
其中α和β為瑞利阻尼系數(shù)，它們之間滿足下式：
(13)
式中ξ和ω分別為阻尼比和圓頻率。對于瑞利阻尼系數(shù)的確定，一般有如表2所示的5種方法[19]。
采用第3節(jié)中相同的計(jì)算模型和材料參數(shù)，外荷載也一致，假定該模型ξ=0.05，ω=1 Hz，則計(jì)算出各種方法的阻尼系數(shù)見表2。同樣利用該程序?qū)υ撃Ｐ瓦M(jìn)行動(dòng)力時(shí)程分析，分析不同方法下的阻尼系數(shù)對地下洞室群響應(yīng)的影響。圖12給出了監(jiān)測點(diǎn)C的加速度時(shí)程響應(yīng)，圖13表示阻尼比與圓頻率之間的關(guān)系。表3展示了該監(jiān)測點(diǎn)加速度峰值與阻尼系數(shù)的關(guān)系。從圖12可以看出，無阻尼得出的加速度時(shí)程曲線與所有方法得出的時(shí)程曲線在峰值附近區(qū)域內(nèi)有所不同，即無阻尼的有2個(gè)極大值點(diǎn)和2個(gè)極小值點(diǎn)，而所有有阻尼的極大值點(diǎn)和極小值點(diǎn)則各僅有1個(gè)，同時(shí)無阻尼的振動(dòng)周期比有阻尼的要??；阻尼對峰值有一定的折減作用，包括最大值和最小值；方法1最大值的折減量略小于其他方案，但最小值的折減量卻明顯大于其他方案，達(dá)到17.89%，可能是由于該方案高估了所有頻率范圍內(nèi)的阻尼；方法4和方法5得出的加速度峰值最接近，時(shí)程響應(yīng)也比較一致。綜合考慮，方法1最簡單，但效果最不好，其余的幾種方法大同小異，沒有一種有特別的優(yōu)勢，因此要精確知道阻尼對洞室群動(dòng)力響應(yīng)的影響還需更進(jìn)一步的研究。
圖12 主變室底板監(jiān)測點(diǎn)C的加速度時(shí)程曲線Fig.12 The acceleration process at monitoring point Cof main transformer chamber floor
圖13 阻尼比與圓頻率的關(guān)系Fig.13 Relationship between damping ratio andcircular frequency
方法峰值/(m·s-2)與無阻尼相比/%最大值最小值最大值折減最小值折減方法11.148-0.9824.8117.89方法21.138-1.0925.648.70方法31.139-1.0685.5410.65方法41.127-1.1146.556.81方法51.124-1.0996.778.08無阻尼1.206-1.19600
6結(jié)論
本文將面向?qū)ο笥邢拊椒☉?yīng)用到地下洞室群中，得出了以下結(jié)論：
(1) 從單元類方面介紹了地下洞室群面向?qū)ο笥邢拊绦?，說明了該程序既有一般有限元的共性，也有其自身的特點(diǎn)。
(2) 應(yīng)用該自編的面向?qū)ο笥邢拊绦蚍治瞿乘娬镜叵聫S房的開挖問題，得出的結(jié)果與Phase2軟件計(jì)算的結(jié)果吻合很好，說明了該程序可以很好地模擬地下洞室開挖問題。
(3) 在該面向?qū)ο笥邢拊绦驅(qū)崿F(xiàn)了人工邊界和地震波輸入，在此基礎(chǔ)上分析了波在地下廠房中的傳播規(guī)律，結(jié)果表明監(jiān)測點(diǎn)響應(yīng)峰值的放大系數(shù)比無洞室作用下的理論解要大。
(4) 在該自編程序中設(shè)置了瑞利阻尼，并分析了不同方法下的阻尼系數(shù)對地下廠房主變室底板監(jiān)測點(diǎn)動(dòng)力響應(yīng)情況的影響，結(jié)果顯示有阻尼的響應(yīng)時(shí)程與無阻尼相比存在一定區(qū)域內(nèi)的不同，且峰值有一定的折減，雖然方法1最簡單，但效果最不好，其余的幾種方法則大同小異，沒有一種有特別的優(yōu)勢，因此要精確知道阻尼對洞室群動(dòng)力響應(yīng)的影響還需更進(jìn)一步的研究。
(5) 該程序雖很好地分析了地下洞室群的靜力和動(dòng)力問題，但僅局限于靜態(tài)本構(gòu)模型問題，因此考慮動(dòng)態(tài)本構(gòu)模型將是該程序以后的重點(diǎn)方向。此外本文輸入的地震荷載都是正弦波，并不能代表實(shí)際的地震波。
參考文獻(xiàn)：
[1]楊志勇,何若全. 面向?qū)ο笥邢拊椒ㄔ诮Y(jié)構(gòu)抗震分析中的應(yīng)用[J]. 哈爾濱建筑大學(xué)學(xué)報(bào),2002,35(4):29-33. (YANG Zhi-yong, HE Ruo-quan. Use of Object-oriented Finite Element Method for Earthquake Analysis of Structure[J]. Journal of Harbin University of Civil Engineering & Architecture, 2002, 35(4): 29-33. (in Chinese))
[2]FORDE B W R, FOSCHI R O, STIEMER S F. Object-oriented Finite Element Analysis[J]. Computers & Structures,1990,34(3):355-374.
[3]ZIMMERMANN T,DUBOIS-PELERIN Y,BOMME P.Object-oriented Finite Element Programming: (I) Governing Principles[J]. Computer Methods in Applied Mechanics and Engineering, 1992,98(2): 291-303.
[4]MILLER G R.Coordinate-free Isoparametric Elements[J]. Computers & Structures,1993,49(6):1027-1035.
[5]李于峰,張亞林. 面向?qū)ο笥邢拊⑿杏?jì)算框架PANDA的并行機(jī)制[J].計(jì)算機(jī)輔助工程,2011,20(1):25-28,41. (LI Yu-feng, ZHANG Ya-lin. Parallel Strategy of Object-oriented Finite Element Parallel Computation Framework PANDA[J]. Computer Aided Engineering, 2011,20(1):25-28,41. (in Chinese))
[6]WERNER H, MACKERT M, STARK M. Object Oriented Models and Tools in Tunnel Design and Analysis∥Proceedings of the 6th International Conference on Computing in Civil and Building Engineering, Rotterdam: Balkema. Berlin, Germany, July 12-15, 1995: 107-112.
[7]平揚(yáng),劉崇權(quán),白世偉 .面向?qū)ο笥邢拊治鲈谏罨庸こ讨械膽?yīng)用[J].電力勘測,1996,(2):8-12. (PING Yang, LIU Chong-quan, BAI Shi-wei. OOP-FEM Analysis and Its Application in Deep Excavation[J]. Electric Power Survey, 1996,(2):8-12. (in Chinese))
[8]項(xiàng)陽,平揚(yáng),葛修潤. 面向?qū)ο笥邢拊椒ㄔ趲r土工程中的應(yīng)用[J]. 巖土力學(xué),2000,21(4):346-349. (XIANG Yang, PING Yang, GE Xiu-run. The Application of Object Oriented FEM(OOFEM) in Geotechnical Engineering[J]. Rock and Soil Mechanics,2000,21(4):346-349. (in Chinese))
[9]劉小生,吳祥松,袁海平. 面向?qū)ο笥邢拊椒ㄔ诖髩螡B流計(jì)算中的應(yīng)用[J]. 有色金屬(礦山部分),2003,55(2):44-46. (LIU Xiao-sheng, WU Xiang-song, YUAN Hai-ping. The Application of Object Oriented FEM(OOFEM) in Dam Seepage[J]. Nonferrous Metals (Mining Section), 2003,55(2):44-46. (in Chinese))
[10]簡榮林,李化敏. 面向?qū)ο笥邢拊椒ㄔ诓傻V工程中的應(yīng)用[J]. 有色金屬(礦山部分),2003,55(1):24-27. (JIAN Rong-lin, LI Hua-min. The Application of Object Oriented FEM(OOFEM) in Mining Engineering[J]. Nonferrous Metals (Mining Section), 2003,55(1):24-27. (in Chinese))
[11]肖專文. 巖土開挖與填筑三維有限元程序的開發(fā)[J]. 中國礦業(yè),1998,7(6):51-55. (XIAO Zhuan-wen. Development of 3-D FEM Program for Rock and Soil Excavation and Filling[J]. China Mining Magazine,1998,7(6):51-55. (in Chinese))
[12]廖振鵬. 工程波動(dòng)理論導(dǎo)論[M]. 北京：科學(xué)出版社，1996. (LIAO Zhen-peng. Introduction to the Theory of Engineering Fluctuations[M]. Beijing: Science Press, 1996. (in Chinese))
[13]李宏恩,李同春,田景元，等. 黏-彈性人工邊界在雙江口土石壩動(dòng)力分析中的應(yīng)用[J]. 巖土力學(xué),2008,29(增刊):189-192. (LI Hong-en, LI Tong-chun, TIAN Jing-yuan,etal. Application of Viscoelastic Artificial Boundary to Dynamic Response Analysis of Shuangjiangkou Rock-fill Dam[J]. Rock and Soil Mechanics, 2008,29(Sup.):189-192. (in Chinese))
[14]劉晶波,谷音,杜義欣. 一致黏彈性人工邊界及黏彈性邊界單元[J]. 巖土工程學(xué)報(bào),2006,28(9):1070-1075. (LIU Jing-bo, GU Yin, DU Yi-xin. Consistent Viscous-spring Artificial Boundaries and Viscous-spring Boundary Elements[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(9): 1070-1075. (in Chinese))
[15]谷音,劉晶波,杜義欣. 三維一致黏彈性人工邊界及等效黏彈性邊界單元[J]. 工程力學(xué),2007,24(12):31-37. (GU Yin, LIU Jing-bo, DU Yi-xin. 3D Consistent Viscous-spring Artificial Boundary and Viscous-spring Boundary Element[J]. Engineering Mechanics, 2007, 24(12): 31-37. (in Chinese))
[16]劉晶波,李彬. 三維黏彈性靜-動(dòng)力統(tǒng)一人工邊界[J]. 中國科學(xué)E輯：工程科學(xué) 材料科學(xué)，2005,35(9):966-980. (LIU Jing-bo, LI Bin. 3D Consistent Viscous-spring Artificial Boundary[J]. Science in China (Series E), 2005,35(9):966-980. (in Chinese))
[17]高峰,趙馮兵. 地下結(jié)構(gòu)靜-動(dòng)力分析中的人工邊界轉(zhuǎn)換方法研究[J]. 振動(dòng)與沖擊,2011,30(11):165-170. (GAO Feng, ZHAO Feng-bing. Study on Transformation Method for Artificial Boundaries in Static-Dynamic Analysis of Underground Structure[J]. Journal of Vibration and Shock, 2011,30(11): 165-170. (in Chinese))
[18]付曉東,盛謙,張勇慧,等. DDA方法在波動(dòng)問題計(jì)算中的初步研究[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版),2012,44(2):70-75. (FU Xiao-dong, SHENG Qian, ZHANG Yong-hui,etal. Preliminary Research on Wave Problem Calculation Using DDA Method[J]. Journal of Sichuan University (Engineering Science Edition), 2012, 44(2): 70-75. (in Chinese))
[19]鄒德高,徐斌,孔憲京. 瑞利阻尼系數(shù)確定方法對高土石壩地震反應(yīng)的影響研究[J]. 巖土力學(xué),2011,32(3):797-803. (ZOU De-gao, XU Bin, KONG Xian-jing. Study of Influence of Different Methods for Calculating Rayleigh Damping Coefficient on High Earth-rock Dam Seismic Response[J]. Rock and Soil Mechanics,2011,32(3):797-803. (in Chinese))
(編輯：劉運(yùn)飛)
Discussion on the Application of Object-Oriented FEMto Underground Caverns
ZHOU Yong-qiang1,2, SHENG Qian1,2, FU Xiao-dong1,2, LIU Xiao-min1,2
(1. Institute of Rock and Soil Mechanics,Chinese Academy of Sciences, Wuhan430071,China;
2.State Key Laboratory of Geomechanics and Geotechnical Engineering,Chinese Academy
of Sciences, Wuhan430071, China)
Abstract:Object-oriented FEM program of underground caverns is briefly introduced with class element as representative. A program code written in VC + + is presented to simulate the stepped excavation of an underground powerhouse and to analyse the deformation. The result is consistent with that of Phase2 software. The seismic input in dealing with dynamic response is implemented through artificial boundary by using object-oriented FEM. Furthermore, the law of dynamic response of the underground powerhouse is analysed and the result shows that the response of monitoring points lags behind and the amplification factor of peak is larger than theoretical solution in the absence of cavern under the same action. On this basis, the influence of damping coefficient on the dynamic response of underground caverns by different methods is researched. The feasibility and applicability of object-oriented FEM in underground caverns are validated. The research result lays a foundation for further study on elastic-plastic behavior under dynamic load.
Key words:underground caverns; object-oriented; FEM; artificial boundary; damping
2015,32(07):89-93
DOI:10.3969/j.issn.1001-5485.2015.07.016
作者簡介：蘭滔(1989-)，男，重慶綦江人，碩士研究生，主要從事港口航道與近海工程方面的研究，(電話)15111973922(電子信箱)lantao944266569@163.com。
基金項(xiàng)目：國家自然科學(xué)基金項(xiàng)目(51209242)
收稿日期：2014-01-08；修回日期：2014-01-25
中圖分類號：TV554
文獻(xiàn)標(biāo)志碼：A
文章編號：1001-5485(2015)07-0082-07
猜你喜歡
N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減數(shù)學(xué)物理學(xué)報(bào)(2021年3期)2021-07-19 06:02:32 關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程數(shù)學(xué)年刊A輯(中文版)(2021年2期)2021-07-17 08:38:02 具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子數(shù)學(xué)物理學(xué)報(bào)(2021年2期)2021-06-09 08:54:38 面向?qū)ο蟮挠?jì)算機(jī)網(wǎng)絡(luò)設(shè)計(jì)軟件系統(tǒng)的開發(fā)電子測試(2018年15期)2018-09-26 06:01:34 面向?qū)ο蟮臄?shù)據(jù)交換協(xié)議研究與應(yīng)用電測與儀表(2017年24期)2017-12-19 05:15:24 烏干達(dá)Karuma水電站地下廠房洞室群圍巖穩(wěn)定分析水利水電工程設(shè)計(jì)(2017年1期)2017-05-17 05:20:03 基于改進(jìn)的非連續(xù)變形方法的洞室圍巖穩(wěn)定性分析中南大學(xué)學(xué)報(bào)（自然科學(xué)版）(2016年2期)2017-01-19 07:37:13 具阻尼項(xiàng)的Boussinesq型方程的長時(shí)間行為數(shù)學(xué)物理學(xué)報(bào)(2016年6期)2016-04-16 04:41:06 面向?qū)ο骔eb開發(fā)編程語言的的評估方法現(xiàn)代計(jì)算機(jī)(2016年12期)2016-02-28 18:35:25 水電站地下洞室群施工通風(fēng)研究水利規(guī)劃與設(shè)計(jì)(2016年7期)2016-02-28 15:06:34
周永強(qiáng)，盛 謙，付曉東，劉曉敏
(中國科學(xué)院 a.武漢巖土力學(xué)研究所;b.巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室， 武漢 430071)
地下洞室群面向?qū)ο笥邢拊椒ǖ奶接?/p>
周永強(qiáng)a，b，盛謙a，b，付曉東a，b，劉曉敏a，b
(中國科學(xué)院 a.武漢巖土力學(xué)研究所;b.巖土力學(xué)與工程國家重點(diǎn)實(shí)驗(yàn)室， 武漢430071)
摘要：以單元為代表類簡單介紹了地下洞室群面向?qū)ο笥邢拊绦?。利用自編的面向?qū)ο笥邢拊绦驅(qū)崿F(xiàn)了開挖模擬，并分析某地下廠房的分步開挖變形問題，其結(jié)果與Phase2軟件計(jì)算結(jié)果基本一致；實(shí)現(xiàn)了面向?qū)ο笥邢拊幚韯?dòng)力響應(yīng)的人工邊界輸入問題，驗(yàn)證了其正確性，并對該地下廠房的動(dòng)力響應(yīng)規(guī)律進(jìn)行了探討。結(jié)果表明：監(jiān)測點(diǎn)響應(yīng)存在一定的滯后性，且峰值的放大系數(shù)比無洞室作用下的理論解要大；在此基礎(chǔ)上，分析了不同方法下的阻尼系數(shù)對地下洞室群動(dòng)力響應(yīng)的影響；驗(yàn)證了面向?qū)ο笥邢拊椒☉?yīng)用在地下洞室群中的可行性和適用性。研究成果為進(jìn)一步研究動(dòng)態(tài)彈塑性情況奠定了基礎(chǔ)。
關(guān)鍵詞：地下洞室群；面向?qū)ο螅挥邢拊绦?；人工邊界；阻?/p>
1研究背景
傳統(tǒng)的有限元程序一般采用結(jié)構(gòu)化設(shè)計(jì)，即采用例如Fortran編程語言的面向過程進(jìn)行編寫，然而隨著工程問題復(fù)雜性的增多和計(jì)算規(guī)模的增大，面向過程的編程語言由于其自身的特點(diǎn)，如代碼重新利用率低、調(diào)試復(fù)雜、開發(fā)周期長等，潛在的問題日益突出，因而已經(jīng)難以滿足開發(fā)和維護(hù)上的要求。
采用面向?qū)ο缶幊陶Z言的有限元程序可以彌補(bǔ)上述不足，可改善有限元程序的性能，提高有限元程序的開發(fā)效率，并能為功能模塊的集成化、擴(kuò)充提供有利條件。面向?qū)ο蟮某绦蛟O(shè)計(jì)方法與有限元方法相結(jié)合，并沒有從根本上改變有限元方法的核心思想，而是程序設(shè)計(jì)和開發(fā)軟件的一種新嘗試。
Forde等[2]在1990年首次提出了面向?qū)ο笥邢拊绦蛟O(shè)計(jì)，之后國內(nèi)外對該設(shè)計(jì)方法做了許多研究。Zimmermann等[3]提出了面向?qū)ο笥邢拊绦蛟O(shè)計(jì)控制原理，Miller等[4]提出了基于面向?qū)ο笥邢拊椒ǖ姆蔷€性動(dòng)力學(xué)分析的軟件框架，李于峰等[5]實(shí)現(xiàn)了面向?qū)ο笥邢拊⑿杏?jì)算框架PANDA，Werner等[6]研究了利用面向?qū)ο蟮哪Ｐ蛠碓O(shè)計(jì)和分析隧道工程問題，平揚(yáng)等[7]實(shí)現(xiàn)了面向?qū)ο笥邢拊谏罨庸こ讨械膽?yīng)用，楊志勇等利用面向?qū)ο笥邢拊椒ǚ治隽私Y(jié)構(gòu)抗震問題，項(xiàng)陽等[8]則提出了面向?qū)ο笥邢拊椒ㄔ趲r土工程中的應(yīng)用實(shí)例，劉小生等[9]應(yīng)用面向?qū)ο笥邢拊椒ń鉀Q大壩滲流問題，簡榮林等[10]利用C++語言開發(fā)了基于Windows平臺的巷道圍巖變形有限元分析軟件。
從上述可以看出，面向?qū)ο笥邢拊讶〉昧讼喈?dāng)大的進(jìn)展，然而該方法應(yīng)用到地下洞室群領(lǐng)域的研究少有報(bào)道。本文首先以單元類作為代表介紹地下洞室群面向?qū)ο笥邢拊绦颍黄浯卫米跃幊绦騺矸治瞿车叵聫S房的開挖變形問題；最后對該地下廠房的動(dòng)力響應(yīng)規(guī)律進(jìn)行了探討；以此為基礎(chǔ)，分析了不同方法下的阻尼系數(shù)對地下洞室群動(dòng)力響應(yīng)的影響。
2面向?qū)ο笥邢拊绦蛟O(shè)計(jì)
有限元的基本思想是先將整個(gè)計(jì)算區(qū)域離散成許多單元，然后在各個(gè)單元上計(jì)算，最后整合在一起，因此單元反映了整個(gè)有限元離散體系的拓?fù)鋵傩?。本?jié)將重點(diǎn)介紹單元類，以此為基礎(chǔ)簡單概括地下洞室群面向?qū)ο笥邢拊绦?。單元類中既包含了材料屬性、單元?jié)點(diǎn)、單元類型、單元編號等基本數(shù)據(jù)屬性，同時(shí)又具備一些計(jì)算功能，如外荷載、單元的應(yīng)力和應(yīng)變、質(zhì)量矩陣和阻尼矩陣等。該程序是用VC++語言來編寫的，單元類的程序如圖1所示。
圖1 單元類程序Fig.1 Class element
其中ComputeNewmarkLhsAt(TimeStep*)和ComputeNewmarkRhsAt(TimeStep*)分別為Newmark法計(jì)算動(dòng)力方程左右兩邊的函數(shù)；giveDampMatrix(int, FloatMatrix*, TimeStep*)是計(jì)算阻尼矩陣的函數(shù)；giveStep()則為計(jì)算開挖步的函數(shù)；而FloatMatrix，F(xiàn)loatArray和IntArray分別是矩陣類、向量類和整數(shù)向量類；GaussPoint為高斯節(jié)點(diǎn)類；TimeStep則為時(shí)間步類，一般用于動(dòng)力計(jì)算。從圖1可以看出，地下洞室群面向?qū)ο笥邢拊绦蚣扔幸话阌邢拊绦虻墓残?，又有其自身的特點(diǎn)，如開挖、分析動(dòng)力響應(yīng)以及考慮阻尼等。
已知t和t+Δt時(shí)刻的動(dòng)力方程分別為
(1)
(2)
假設(shè)Δut=ut+Δt-ut，則根據(jù)Newmark法有
(3)
結(jié)合式(3),由式(2)-式(1)可得
(4)
(5)
式中:K,M和C分別為剛度矩陣、質(zhì)量矩陣和瑞利阻尼矩陣;α和δ是Newmark法參數(shù)；γ和β是不依賴頻率的參數(shù)；ρ為材料密度。對于解的穩(wěn)定性，Newmark法所選用的Δt的值是任意的，但考慮到精度要求，Δt則需小于體系固有振型最小周期的(1/10～1/20)。根據(jù)式(4)和式(5)可知，單元類的組裝過程如圖2所示。
圖2 單元類的組裝過程Fig.2 Assembly process of class element
3地下洞室群的開挖模擬
3.1基本思路
模擬地下洞室的開挖過程，重點(diǎn)是把作用在需要開挖部分的重力去除(不考慮構(gòu)造應(yīng)力)。地下洞室群在開挖之前，重力荷載是利用質(zhì)量矩陣轉(zhuǎn)化成等效節(jié)點(diǎn)作用力來實(shí)現(xiàn)的，因此為模擬地下洞室群開挖過程中重力作用的變化，只需根據(jù)時(shí)間步將質(zhì)量矩陣進(jìn)行修改，而材料密度是影響質(zhì)量矩陣的唯一因素。為了防止開挖的單元不起作用而使總體剛度矩陣產(chǎn)生奇異，開挖單元的彈性模量不能取得太大，也不能取得太小，一般比周圍巖土體彈性模量小5至8個(gè)數(shù)量級[11]，具體思路如圖3所示。
圖3 地下洞室群開挖過程的思路Fig.3 Idea of the excavation process ofunderground caverns
圖4 半無限空間體簡化圖Fig.4 Simplified diagramof half infinite space body
為驗(yàn)證該程序的開挖思想的正確性，設(shè)有半無限空間體，如圖4所示，圖中數(shù)字表示單元節(jié)點(diǎn)。材料的彈性模量為1.2×104MPa，泊松比為0.2，密度為2.7×103kg/m3。每個(gè)單元的尺寸大小均為6.25 m×25 m。初始時(shí)該體系不開挖，之后開挖最上面的單元，則開挖后單元節(jié)點(diǎn)的位移為[11]
(6)
式中:u0為初始位移，Δu為變形回彈位移?？煞謩e按下式計(jì)算：
(7)
(8)
式中：Δh為開挖的豎向高度;υ為材料的泊松比。運(yùn)用自編程序?qū)w系進(jìn)行開挖模擬，其結(jié)果如表1所示。從表1中看出，數(shù)值解與解析解基本一致。
表1 體系開挖步1后節(jié)點(diǎn)位移的數(shù)值解和解析解
3.2實(shí)例驗(yàn)證
某水電站地下廠房某斷面的計(jì)算剖面如圖5所示，其主廠房、主變室和尾調(diào)室的尺寸大小分別為29 m×57 m，20 m×37 m，26 m×87 m，整個(gè)模型尺寸為1 500 m×600 m。材料的彈性模量為1.2×104MPa，泊松比為0.2，密度為2.7×103kg/m3。
圖5 地下洞室群某計(jì)算斷面Fig.5 Cross-section of underground cavernsfor computation
圖6 地下洞室群開挖工序及監(jiān)測點(diǎn)布置Fig.6 Excavation process of underground cavernsand layout of monitoring points
圖7 地下洞室群主變室各測點(diǎn)開挖變形Fig.7 Deformation of monitoring points in maintransformer chamber of underground caverns
該地下廠房共分為10步開挖，開挖工序及廠房的監(jiān)測點(diǎn)的布置如圖6所示。利用自編程序來模擬該地下廠房的開挖過程。利用Phase2軟件建立相同的計(jì)算模型，材料參數(shù)保持不變，得到相應(yīng)監(jiān)測點(diǎn)的結(jié)果。圖7為地下洞室群主變室各監(jiān)測點(diǎn)的位移隨開挖步數(shù)的變化圖，限于篇幅，其他洞室就不一一給出。從圖7中可知，自編程序與Phase2軟件計(jì)算的結(jié)果基本一致，說明該程序的開挖基本思路是正確的，也說明了該程序能正確分析地下洞室群的開挖問題。
4地下洞室群的波動(dòng)影響
4.1人工邊界及地震輸入
要知道地下洞室群的波動(dòng)影響，首先要解決其地震輸入問題，即動(dòng)力荷載在邊界上的輸入。動(dòng)力荷載的輸入方法有2種：位移輸入和應(yīng)力輸入，本程序采用應(yīng)力輸入，即在輸入邊界上設(shè)置人工邊界，并同時(shí)把速度時(shí)程通過式(9)和式(10)轉(zhuǎn)換成應(yīng)力時(shí)程。
(9)
(10)
式中：vn和vs分別為輸入波在模型邊界上質(zhì)點(diǎn)的法向和切向的速度分量；ρ為材料密度；CP和CS分別為P波和S波在介質(zhì)中的傳播波速。
其中人工邊界的作用是吸收邊界上的入射波，減少因邊界上存在反射波而對動(dòng)力分析的結(jié)果產(chǎn)生影響。人工邊界大致可以分為局部人工邊界和全局人工邊界[12]，局部人工邊界由于其計(jì)算時(shí)間少、效率高等特點(diǎn)得到了廣泛的應(yīng)用。局部人工邊界一般包括黏性邊界、黏彈性邊界、旁軸邊界和透射邊界等4種邊界[13]，因黏性邊界簡單，易于編程實(shí)現(xiàn)，且效果較好，綜合考慮本程序采用了黏性邊界。為體現(xiàn)面向?qū)ο蟪绦虻亩鄳B(tài)性和繼承性，將含有黏性邊界的單元發(fā)展成黏性邊界單元，即其剛度矩陣和質(zhì)量矩陣保持不變，阻尼矩陣則為含有邊界特性阻尼矩陣，即根據(jù)虛功原理，采用類似質(zhì)量矩陣的集中形式(四邊形單元)，即
(11)
在分析動(dòng)力情況之前，往往先考慮體系的靜力平衡問題。劉晶波等[14-15]雖已將二維黏彈性邊界發(fā)展成黏彈性邊界單元，又推導(dǎo)了三維等效黏彈性邊界單元，并在此基礎(chǔ)上建立了三維黏彈性靜動(dòng)力統(tǒng)一人工邊界[16]，然而用此邊界來分析靜力問題時(shí)，會(huì)有一定的誤差[17]，且參數(shù)的選取有一定的主觀性。另外計(jì)算靜力問題(或動(dòng)力問題)，選取相對應(yīng)的靜力(或動(dòng)力)分析方法，然而靜力平衡之后分析動(dòng)力問題，若不去掉約束條件，同時(shí)又不施加約束反力，則整個(gè)體系是在不平衡狀態(tài)下計(jì)算動(dòng)力荷載，因此得出的結(jié)果會(huì)出現(xiàn)很大的波動(dòng)，與實(shí)際結(jié)果不相符。因此分析靜動(dòng)力問題時(shí)，應(yīng)選擇的思路為：以靜力分析結(jié)果的應(yīng)力場作為初始條件進(jìn)行動(dòng)力分析時(shí)，為保證體系在動(dòng)力計(jì)算初始時(shí)刻處于靜力平衡狀態(tài)，可在輸入靜應(yīng)力場的同時(shí)輸入引起此應(yīng)力場的荷載和約束，即在轉(zhuǎn)換邊界時(shí)，需去掉約束條件，同時(shí)施加相應(yīng)的約束反力。
為驗(yàn)證黏性邊界和地震波輸入方法在該程序中的可行性，建立了圖8所示的 200 m×2 000 m的計(jì)算模型，模型底部設(shè)置了固定約束，兩側(cè)設(shè)置了法向約束。材料參數(shù)：彈性模量為1.2×104MPa，泊松比為0.2，密度為2.7×103kg/m3。不考慮阻尼的影響，在模型左邊輸入一個(gè)周期的正弦速度波，頻率為1 Hz，峰值為1 m/s。在模型左右兩邊設(shè)置監(jiān)測點(diǎn)A和B，分別考慮在模型左邊設(shè)置黏性邊界和不設(shè)置黏性邊界2種情況，觀察監(jiān)測點(diǎn)的速度響應(yīng)時(shí)程，其結(jié)果如圖9所示。波從模型左端傳到右端，由于自由邊界反射后左行波疊加作用，使得監(jiān)測點(diǎn)B的峰值速度是入射波峰值速度的2倍，即放大了2倍，這一現(xiàn)象與理論解一致，且與文獻(xiàn)[18]得出的結(jié)果一致，同時(shí)黏性邊界較好地模擬了左行波穿過人工邊界向無限域輻射的過程。
圖9 各監(jiān)測點(diǎn)的速度時(shí)程曲線Fig.9 Curves of velocity vs. time at monitoring points
圖8動(dòng)力計(jì)算模型圖
Fig.8Model for dynamic computation
4.2波在地下洞室群的傳播規(guī)律
圖11 主變室各監(jiān)測點(diǎn)加速度時(shí)程峰值的放大倍數(shù)Fig.11 Peak amplification factor of accelerationprocess at monitoring points of maintransformer chamber
基于上述黏性邊界理論和地震波輸入方法的正確性，以上述水電站地下廠房某斷面(地下洞室群開挖完成之后)作為實(shí)例，各種材料參數(shù)和尺寸都保持不變，不考慮阻尼的影響，在模型底端輸入一個(gè)周期的正弦加速度波(頻率為1 Hz，峰值為1 m/s2)，并施加黏性邊界。為研究地下洞室群的動(dòng)力響應(yīng)，監(jiān)測布置點(diǎn)如圖5和圖6所示。同樣利用該自編程序計(jì)算，得到地下洞室群各監(jiān)測點(diǎn)的加速度響應(yīng)見圖10和圖11。由圖10可知，各監(jiān)測點(diǎn)動(dòng)力響應(yīng)存在滯后性，離模型底端越遠(yuǎn)，滯后的時(shí)間越長，如監(jiān)測點(diǎn)D滯后了約0.05 s，而監(jiān)測點(diǎn)B則大概滯后了0.22 s；由于洞室的影響，監(jiān)測點(diǎn)B的最大放大系數(shù)為2.08，大于2.0，監(jiān)測點(diǎn)C的最大放大倍數(shù)為1.21，也超過了1.0；0.57 s是監(jiān)測點(diǎn)A的波形時(shí)變化的轉(zhuǎn)折點(diǎn)，即0.57 s之前基本與輸入波形相同，之后則完全不一致，這是輸入波在洞室中產(chǎn)生反射波與輸入波疊加作用的結(jié)果。綜上分析，洞室的存在會(huì)增加地震波對圍巖變形的影響，因而需更加關(guān)注洞室頂拱和頂板的變形。
表2 瑞利阻尼系數(shù)的確定方法
5阻尼對地下洞室群動(dòng)力響應(yīng)的影響
阻尼的產(chǎn)生主要源于材料內(nèi)部摩擦以及可能存在的接觸表面的滑動(dòng)。在實(shí)際分析中，精確確定阻尼矩陣是相當(dāng)困難的，通常是將阻尼矩陣表示為質(zhì)量矩陣和剛度矩陣的線性組合(瑞利阻尼)，即
(12)
其中α和β為瑞利阻尼系數(shù)，它們之間滿足下式：
(13)
式中ξ和ω分別為阻尼比和圓頻率。對于瑞利阻尼系數(shù)的確定，一般有如表2所示的5種方法[19]。
采用第3節(jié)中相同的計(jì)算模型和材料參數(shù)，外荷載也一致，假定該模型ξ=0.05，ω=1 Hz，則計(jì)算出各種方法的阻尼系數(shù)見表2。同樣利用該程序?qū)υ撃Ｐ瓦M(jìn)行動(dòng)力時(shí)程分析，分析不同方法下的阻尼系數(shù)對地下洞室群響應(yīng)的影響。圖12給出了監(jiān)測點(diǎn)C的加速度時(shí)程響應(yīng)，圖13表示阻尼比與圓頻率之間的關(guān)系。表3展示了該監(jiān)測點(diǎn)加速度峰值與阻尼系數(shù)的關(guān)系。從圖12可以看出，無阻尼得出的加速度時(shí)程曲線與所有方法得出的時(shí)程曲線在峰值附近區(qū)域內(nèi)有所不同，即無阻尼的有2個(gè)極大值點(diǎn)和2個(gè)極小值點(diǎn)，而所有有阻尼的極大值點(diǎn)和極小值點(diǎn)則各僅有1個(gè)，同時(shí)無阻尼的振動(dòng)周期比有阻尼的要??；阻尼對峰值有一定的折減作用，包括最大值和最小值；方法1最大值的折減量略小于其他方案，但最小值的折減量卻明顯大于其他方案，達(dá)到17.89%，可能是由于該方案高估了所有頻率范圍內(nèi)的阻尼；方法4和方法5得出的加速度峰值最接近，時(shí)程響應(yīng)也比較一致。綜合考慮，方法1最簡單，但效果最不好，其余的幾種方法大同小異，沒有一種有特別的優(yōu)勢，因此要精確知道阻尼對洞室群動(dòng)力響應(yīng)的影響還需更進(jìn)一步的研究。
圖12 主變室底板監(jiān)測點(diǎn)C的加速度時(shí)程曲線Fig.12 The acceleration process at monitoring point Cof main transformer chamber floor
圖13 阻尼比與圓頻率的關(guān)系Fig.13 Relationship between damping ratio andcircular frequency
方法峰值/(m·s-2)與無阻尼相比/%最大值最小值最大值折減最小值折減方法11.148-0.9824.8117.89方法21.138-1.0925.648.70方法31.139-1.0685.5410.65方法41.127-1.1146.556.81方法51.124-1.0996.778.08無阻尼1.206-1.19600
6結(jié)論
本文將面向?qū)ο笥邢拊椒☉?yīng)用到地下洞室群中，得出了以下結(jié)論：
(1) 從單元類方面介紹了地下洞室群面向?qū)ο笥邢拊绦?，說明了該程序既有一般有限元的共性，也有其自身的特點(diǎn)。
(2) 應(yīng)用該自編的面向?qū)ο笥邢拊绦蚍治瞿乘娬镜叵聫S房的開挖問題，得出的結(jié)果與Phase2軟件計(jì)算的結(jié)果吻合很好，說明了該程序可以很好地模擬地下洞室開挖問題。
(3) 在該面向?qū)ο笥邢拊绦驅(qū)崿F(xiàn)了人工邊界和地震波輸入，在此基礎(chǔ)上分析了波在地下廠房中的傳播規(guī)律，結(jié)果表明監(jiān)測點(diǎn)響應(yīng)峰值的放大系數(shù)比無洞室作用下的理論解要大。
(4) 在該自編程序中設(shè)置了瑞利阻尼，并分析了不同方法下的阻尼系數(shù)對地下廠房主變室底板監(jiān)測點(diǎn)動(dòng)力響應(yīng)情況的影響，結(jié)果顯示有阻尼的響應(yīng)時(shí)程與無阻尼相比存在一定區(qū)域內(nèi)的不同，且峰值有一定的折減，雖然方法1最簡單，但效果最不好，其余的幾種方法則大同小異，沒有一種有特別的優(yōu)勢，因此要精確知道阻尼對洞室群動(dòng)力響應(yīng)的影響還需更進(jìn)一步的研究。
(5) 該程序雖很好地分析了地下洞室群的靜力和動(dòng)力問題，但僅局限于靜態(tài)本構(gòu)模型問題，因此考慮動(dòng)態(tài)本構(gòu)模型將是該程序以后的重點(diǎn)方向。此外本文輸入的地震荷載都是正弦波，并不能代表實(shí)際的地震波。
參考文獻(xiàn)：
[1]楊志勇,何若全. 面向?qū)ο笥邢拊椒ㄔ诮Y(jié)構(gòu)抗震分析中的應(yīng)用[J]. 哈爾濱建筑大學(xué)學(xué)報(bào),2002,35(4):29-33. (YANG Zhi-yong, HE Ruo-quan. Use of Object-oriented Finite Element Method for Earthquake Analysis of Structure[J]. Journal of Harbin University of Civil Engineering & Architecture, 2002, 35(4): 29-33. (in Chinese))
[2]FORDE B W R, FOSCHI R O, STIEMER S F. Object-oriented Finite Element Analysis[J]. Computers & Structures,1990,34(3):355-374.
[3]ZIMMERMANN T,DUBOIS-PELERIN Y,BOMME P.Object-oriented Finite Element Programming: (I) Governing Principles[J]. Computer Methods in Applied Mechanics and Engineering, 1992,98(2): 291-303.
[4]MILLER G R.Coordinate-free Isoparametric Elements[J]. Computers & Structures,1993,49(6):1027-1035.
[5]李于峰,張亞林. 面向?qū)ο笥邢拊⑿杏?jì)算框架PANDA的并行機(jī)制[J].計(jì)算機(jī)輔助工程,2011,20(1):25-28,41. (LI Yu-feng, ZHANG Ya-lin. Parallel Strategy of Object-oriented Finite Element Parallel Computation Framework PANDA[J]. Computer Aided Engineering, 2011,20(1):25-28,41. (in Chinese))
[6]WERNER H, MACKERT M, STARK M. Object Oriented Models and Tools in Tunnel Design and Analysis∥Proceedings of the 6th International Conference on Computing in Civil and Building Engineering, Rotterdam: Balkema. Berlin, Germany, July 12-15, 1995: 107-112.
[7]平揚(yáng),劉崇權(quán),白世偉 .面向?qū)ο笥邢拊治鲈谏罨庸こ讨械膽?yīng)用[J].電力勘測,1996,(2):8-12. (PING Yang, LIU Chong-quan, BAI Shi-wei. OOP-FEM Analysis and Its Application in Deep Excavation[J]. Electric Power Survey, 1996,(2):8-12. (in Chinese))
[8]項(xiàng)陽,平揚(yáng),葛修潤. 面向?qū)ο笥邢拊椒ㄔ趲r土工程中的應(yīng)用[J]. 巖土力學(xué),2000,21(4):346-349. (XIANG Yang, PING Yang, GE Xiu-run. The Application of Object Oriented FEM(OOFEM) in Geotechnical Engineering[J]. Rock and Soil Mechanics,2000,21(4):346-349. (in Chinese))
[9]劉小生,吳祥松,袁海平. 面向?qū)ο笥邢拊椒ㄔ诖髩螡B流計(jì)算中的應(yīng)用[J]. 有色金屬(礦山部分),2003,55(2):44-46. (LIU Xiao-sheng, WU Xiang-song, YUAN Hai-ping. The Application of Object Oriented FEM(OOFEM) in Dam Seepage[J]. Nonferrous Metals (Mining Section), 2003,55(2):44-46. (in Chinese))
[10]簡榮林,李化敏. 面向?qū)ο笥邢拊椒ㄔ诓傻V工程中的應(yīng)用[J]. 有色金屬(礦山部分),2003,55(1):24-27. (JIAN Rong-lin, LI Hua-min. The Application of Object Oriented FEM(OOFEM) in Mining Engineering[J]. Nonferrous Metals (Mining Section), 2003,55(1):24-27. (in Chinese))
[11]肖專文. 巖土開挖與填筑三維有限元程序的開發(fā)[J]. 中國礦業(yè),1998,7(6):51-55. (XIAO Zhuan-wen. Development of 3-D FEM Program for Rock and Soil Excavation and Filling[J]. China Mining Magazine,1998,7(6):51-55. (in Chinese))
[12]廖振鵬. 工程波動(dòng)理論導(dǎo)論[M]. 北京：科學(xué)出版社，1996. (LIAO Zhen-peng. Introduction to the Theory of Engineering Fluctuations[M]. Beijing: Science Press, 1996. (in Chinese))
[13]李宏恩,李同春,田景元，等. 黏-彈性人工邊界在雙江口土石壩動(dòng)力分析中的應(yīng)用[J]. 巖土力學(xué),2008,29(增刊):189-192. (LI Hong-en, LI Tong-chun, TIAN Jing-yuan,etal. Application of Viscoelastic Artificial Boundary to Dynamic Response Analysis of Shuangjiangkou Rock-fill Dam[J]. Rock and Soil Mechanics, 2008,29(Sup.):189-192. (in Chinese))
[14]劉晶波,谷音,杜義欣. 一致黏彈性人工邊界及黏彈性邊界單元[J]. 巖土工程學(xué)報(bào),2006,28(9):1070-1075. (LIU Jing-bo, GU Yin, DU Yi-xin. Consistent Viscous-spring Artificial Boundaries and Viscous-spring Boundary Elements[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(9): 1070-1075. (in Chinese))
[15]谷音,劉晶波,杜義欣. 三維一致黏彈性人工邊界及等效黏彈性邊界單元[J]. 工程力學(xué),2007,24(12):31-37. (GU Yin, LIU Jing-bo, DU Yi-xin. 3D Consistent Viscous-spring Artificial Boundary and Viscous-spring Boundary Element[J]. Engineering Mechanics, 2007, 24(12): 31-37. (in Chinese))
[16]劉晶波,李彬. 三維黏彈性靜-動(dòng)力統(tǒng)一人工邊界[J]. 中國科學(xué)E輯：工程科學(xué) 材料科學(xué)，2005,35(9):966-980. (LIU Jing-bo, LI Bin. 3D Consistent Viscous-spring Artificial Boundary[J]. Science in China (Series E), 2005,35(9):966-980. (in Chinese))
[17]高峰,趙馮兵. 地下結(jié)構(gòu)靜-動(dòng)力分析中的人工邊界轉(zhuǎn)換方法研究[J]. 振動(dòng)與沖擊,2011,30(11):165-170. (GAO Feng, ZHAO Feng-bing. Study on Transformation Method for Artificial Boundaries in Static-Dynamic Analysis of Underground Structure[J]. Journal of Vibration and Shock, 2011,30(11): 165-170. (in Chinese))
[18]付曉東,盛謙,張勇慧,等. DDA方法在波動(dòng)問題計(jì)算中的初步研究[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版),2012,44(2):70-75. (FU Xiao-dong, SHENG Qian, ZHANG Yong-hui,etal. Preliminary Research on Wave Problem Calculation Using DDA Method[J]. Journal of Sichuan University (Engineering Science Edition), 2012, 44(2): 70-75. (in Chinese))
[19]鄒德高,徐斌,孔憲京. 瑞利阻尼系數(shù)確定方法對高土石壩地震反應(yīng)的影響研究[J]. 巖土力學(xué),2011,32(3):797-803. (ZOU De-gao, XU Bin, KONG Xian-jing. Study of Influence of Different Methods for Calculating Rayleigh Damping Coefficient on High Earth-rock Dam Seismic Response[J]. Rock and Soil Mechanics,2011,32(3):797-803. (in Chinese))
(編輯：劉運(yùn)飛)
Discussion on the Application of Object-Oriented FEMto Underground Caverns
ZHOU Yong-qiang1,2, SHENG Qian1,2, FU Xiao-dong1,2, LIU Xiao-min1,2
(1. Institute of Rock and Soil Mechanics,Chinese Academy of Sciences, Wuhan430071,China;
2.State Key Laboratory of Geomechanics and Geotechnical Engineering,Chinese Academy
of Sciences, Wuhan430071, China)
Abstract:Object-oriented FEM program of underground caverns is briefly introduced with class element as representative. A program code written in VC + + is presented to simulate the stepped excavation of an underground powerhouse and to analyse the deformation. The result is consistent with that of Phase2 software. The seismic input in dealing with dynamic response is implemented through artificial boundary by using object-oriented FEM. Furthermore, the law of dynamic response of the underground powerhouse is analysed and the result shows that the response of monitoring points lags behind and the amplification factor of peak is larger than theoretical solution in the absence of cavern under the same action. On this basis, the influence of damping coefficient on the dynamic response of underground caverns by different methods is researched. The feasibility and applicability of object-oriented FEM in underground caverns are validated. The research result lays a foundation for further study on elastic-plastic behavior under dynamic load.
Key words:underground caverns; object-oriented; FEM; artificial boundary; damping
2015,32(07):89-93
DOI:10.3969/j.issn.1001-5485.2015.07.016
作者簡介：蘭滔(1989-)，男，重慶綦江人，碩士研究生，主要從事港口航道與近海工程方面的研究，(電話)15111973922(電子信箱)lantao944266569@163.com。
基金項(xiàng)目：國家自然科學(xué)基金項(xiàng)目(51209242)
收稿日期：2014-01-08；修回日期：2014-01-25
中圖分類號：TV554
文獻(xiàn)標(biāo)志碼：A
文章編號：1001-5485(2015)07-0082-07
猜你喜歡
數(shù)學(xué)物理學(xué)報(bào)(2021年3期)2021-07-19 06:02:32
數(shù)學(xué)年刊A輯(中文版)(2021年2期)2021-07-17 08:38:02
數(shù)學(xué)物理學(xué)報(bào)(2021年2期)2021-06-09 08:54:38
電子測試(2018年15期)2018-09-26 06:01:34
電測與儀表(2017年24期)2017-12-19 05:15:24
水利水電工程設(shè)計(jì)(2017年1期)2017-05-17 05:20:03
中南大學(xué)學(xué)報(bào)（自然科學(xué)版）(2016年2期)2017-01-19 07:37:13
數(shù)學(xué)物理學(xué)報(bào)(2016年6期)2016-04-16 04:41:06
現(xiàn)代計(jì)算機(jī)(2016年12期)2016-02-28 18:35:25
水利規(guī)劃與設(shè)計(jì)(2016年7期)2016-02-28 15:06:34
