張仕華

摘要 針對(duì)缺水地區(qū)常用的圓柱形、球形、瓶形水窖建設(shè)中存在的因窖體受力分析不足、窖體結(jié)構(gòu)設(shè)計(jì)不合理而導(dǎo)致窖體空間結(jié)構(gòu)穩(wěn)定性差以至最后窖體破壞以及盲目增厚窖壁導(dǎo)致材料浪費(fèi)等問(wèn)題，利用基于ANSYS平臺(tái)新開(kāi)發(fā)的窖體有限元分析軟件，對(duì)圓柱形、球形、瓶形3種水窖進(jìn)行空間結(jié)構(gòu)幾何尺寸優(yōu)化。以水窖耗材量最小為目標(biāo)函數(shù)，以許用應(yīng)力為約束條件，優(yōu)化出不同容水量的水窖幾何尺寸，并將它們的最優(yōu)結(jié)果進(jìn)行對(duì)比分析。結(jié)果表明，修筑球形水窖最省材料，圓柱形水窖次之，瓶形水窖在同工況同容積下相對(duì)于另外2種水窖設(shè)計(jì)形式耗材最多，但結(jié)構(gòu)穩(wěn)定性高于圓柱形水窖，它們各有優(yōu)缺點(diǎn)，可根據(jù)當(dāng)?shù)鼐唧w情況具體用途來(lái)選擇最合適的水窖形式。

關(guān)鍵詞 水窖；有限元；結(jié)構(gòu)優(yōu)化

中圖分類(lèi)號(hào) S277 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2018）20-0186-05

Abstract In the construction of cylindrical， spherical， and bottle shaped water cellar， there are some problems in the construction of cylindrical， spherical and bottle shaped water cellar， which are due to the lack of stress analysis of the cellar body and the unreasonable design of the cellar structure， which lead to the poor stability of the space structure of the cellar， and the waste of the pit wall resulting from the blind pit wall. Based on the newly developed cellar finite element analysis software based on ANSYS platform， the geometric dimensions of 3 types of cellars， cylindrical， spherical and bottle shaped， are optimized.With the minimum consumption of water cellar as the objective function and the allowable stress as the constraint condition， the geometric sizes of water cellars with different water volume are optimized and their optimal results are compared and analyzed.The results show that the construction of the spherical water cellar is the most provincial material， the cylindrical water cellar is the second， and the bottle shaped water cellar is the most suitable for the other two kinds of water cellar design in the same volume， but the structural stability is higher than that of the cylindrical water cellar.Each has its own advantages and disadvantages. It can choose the most suitable form of water cellar according to the specific circumstances and specific uses of the site.

Key words Cellar；Finite element method；Structure optimization

近年來(lái)，隨著社會(huì)的發(fā)展和人口的劇增，淡水資源緊缺的程度也越來(lái)越嚴(yán)重，水窖作為收集雨水的一種典型設(shè)施得到越來(lái)越廣泛的應(yīng)用。圓柱形、球形和瓶形水窖是我國(guó)干旱缺水地區(qū)使用最廣泛的集雨設(shè)施形式[1]，它們以施工簡(jiǎn)單快捷、工程量小、窖體結(jié)構(gòu)穩(wěn)定性好等優(yōu)點(diǎn)受到廣大建窖者與農(nóng)戶(hù)的青睞。但現(xiàn)行的水窖都是以建窖者根據(jù)自己的修筑經(jīng)驗(yàn)來(lái)建造，沒(méi)有經(jīng)過(guò)精確的計(jì)算來(lái)確定窖體結(jié)構(gòu)幾何尺寸，往往使水窖修筑者盲目地認(rèn)為只要增加窖壁厚度就能保證窖體的使用壽命與安全，這樣修筑出來(lái)的水窖既導(dǎo)致了材料的浪費(fèi)還因窖體結(jié)構(gòu)尺寸（如：矢跨比、半徑與高度的比值、窖壁的邊坡比等）設(shè)計(jì)不合理使水窖結(jié)構(gòu)穩(wěn)定性得不到保證。所以，有必要對(duì)水窖結(jié)構(gòu)進(jìn)行強(qiáng)度驗(yàn)算并優(yōu)化出合理的尺寸。由于水窖結(jié)構(gòu)的三維尺寸差異較小，所以在結(jié)構(gòu)計(jì)算中必須以空間結(jié)構(gòu)原理為基礎(chǔ)，計(jì)算工作量大而且復(fù)雜，必須借助于適當(dāng)?shù)挠?jì)算機(jī)軟件才能完成，ANSYS軟件的出現(xiàn)為解決這一問(wèn)題提供了有力的支持。截至目前，鮮有對(duì)水窖進(jìn)行空間結(jié)構(gòu)分析的相關(guān)研究。因此，筆者借助于大型有限元分析軟件ANSYS對(duì)圓柱形、球形和瓶形3種水窖空間結(jié)構(gòu)尺寸進(jìn)行優(yōu)化分析，設(shè)計(jì)出既安全又經(jīng)濟(jì)的水窖結(jié)構(gòu)形式，并將3種水窖的優(yōu)化結(jié)果進(jìn)行對(duì)比，在滿(mǎn)足窖體結(jié)構(gòu)穩(wěn)定性的基礎(chǔ)上從節(jié)省材料出發(fā)選取一種最經(jīng)濟(jì)的水窖形式，為水窖的建設(shè)者提供技術(shù)支撐，并為新型水窖的開(kāi)發(fā)提供理論基礎(chǔ)。

1 水窖有限元法分析模型的建立

水窖有限元法分析的基本思想是將窖體結(jié)構(gòu)離散為若干個(gè)僅由節(jié)點(diǎn)相連的單元，在進(jìn)行有限元分析時(shí)，各種荷載轉(zhuǎn)移至各個(gè)單元與節(jié)點(diǎn)上，利用插值函數(shù)考慮連續(xù)邊界條件，建立有限元求解方程組，并表示為規(guī)范化的矩陣形式，求得窖體結(jié)構(gòu)各部分所受到的應(yīng)力與產(chǎn)生的位移[2-3]。

該研究首先建立圓柱形水窖、球形水窖、瓶形水窖3個(gè)實(shí)體模型，3種類(lèi)型的水窖斷面形狀如圖1所示。球形水窖半徑為R，窖壁厚度為T(mén)，窖頸高為SH，取水口直徑一般為0.7 m[4]，圈梁厚度和高度一般為15 cm；圓柱形水窖和瓶形水窖結(jié)構(gòu)尺寸見(jiàn)圖1，字母代表意義與球形水窖一致。

實(shí)體模型建立后需要?jiǎng)澐志W(wǎng)格以建立有限元模型，首先采用8節(jié)點(diǎn)四邊形單元Plane82對(duì)窖體橫截面進(jìn)行網(wǎng)格劃分，然后再采用三維20節(jié)點(diǎn)結(jié)構(gòu)實(shí)體單元Solid95通過(guò)橫截面旋轉(zhuǎn)掃略生成三維水窖實(shí)體的網(wǎng)格化有限元模型[5]。由于圓柱形、球形和瓶形3種水窖形狀和荷載均具有旋轉(zhuǎn)對(duì)稱(chēng)性，所以分析時(shí)只對(duì)1/4窖體進(jìn)行結(jié)構(gòu)分析。水窖網(wǎng)格劃分如圖2所示。

2 水窖不利工況及荷載分析

根據(jù)《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》C15混凝土極限伸長(zhǎng)量不超過(guò)0.003 3 m，而極限拉壓縮量不超過(guò)0.000 15 m，窖壁厚度T的優(yōu)化取值范圍按以前修筑水窖窖壁厚度范圍設(shè)定為1～20 cm。圓柱半徑R的優(yōu)化取值范圍設(shè)定也按以往修筑水窖經(jīng)驗(yàn)，要求為首先窖體容積固定，R在0.2H～2.0H取值，在窖容為20、40和60 m3時(shí)R取值范圍分別為1.039 6～1.760 5 m、1.309 8～2.218 1 m和1.499 3～2.539 1 m，圓角半徑R1的取值范圍為1/8R～R。

決策變量：窖壁厚度T；圓柱半徑R和窖底圓角半徑R1。

優(yōu)化方法：ANSYS優(yōu)化方法有單步運(yùn)行法（single run）、隨機(jī)法（random designs）、乘子法（factorial）、最優(yōu)梯度法（gradient）、掃描法（DV Sweeps）、子問(wèn)題法（subproblem）、一階優(yōu)化工具（firstorder）、用戶(hù)優(yōu)化算法（user optimizer）[8]。該優(yōu)化方法采用子問(wèn)題法（subproblem）對(duì)窖體結(jié)構(gòu)幾何尺寸進(jìn)行優(yōu)化，子問(wèn)題法是按照單一步長(zhǎng)在每次計(jì)算以后將設(shè)計(jì)變量在變化范圍內(nèi)加以改變，對(duì)于目標(biāo)函數(shù)和狀態(tài)變量的整體變化評(píng)估可以用本工具實(shí)現(xiàn)。子問(wèn)題法屬于ANSYS優(yōu)化工具中的零階方法，是使用所有因變量（狀態(tài)變量和目標(biāo)函數(shù)）的逼近，對(duì)結(jié)構(gòu)簡(jiǎn)單的模型能得到比較精確的解，選擇該方法后需要指定參考設(shè)計(jì)序列（即初始設(shè)計(jì)序列）和最大可行性設(shè)計(jì)序列的數(shù)目，該試驗(yàn)將兩者設(shè)置為10和7，也可根據(jù)需要增大設(shè)計(jì)序列的數(shù)目，一般初始設(shè)計(jì)序列數(shù)目30，最大可行性設(shè)計(jì)序列10都能滿(mǎn)足一般優(yōu)化要求即得到較為精確的解。

3.2 圓柱形水窖不同容水量?jī)?yōu)化結(jié)果

通過(guò)ANSYS優(yōu)化得出的最優(yōu)化結(jié)果如表1所示。

從表1優(yōu)化結(jié)果可以看出，圓柱形水窖在空載時(shí)下，容積為20、40、60 m3水窖的最大拉應(yīng)力分別為895.18、883.80、884.99 kPa；滿(mǎn)載時(shí)，容積為20、40、60 m3水窖的最大拉應(yīng)力分別為789.05、907.33、874.32 kPa。2種工況下，不同容水量窖體受到的最大拉應(yīng)力值范圍為789.05～907.33 kPa，略小于材料抗拉強(qiáng)度910 kPa；最大壓應(yīng)力值范圍296.61～1 965.6 kPa，遠(yuǎn)遠(yuǎn)小于材料抗壓強(qiáng)度7 200 kPa，同時(shí)窖體計(jì)算最大位移0.000 06～0.000 14 m都小于材料極限拉應(yīng)力值0.000 15 m，因此結(jié)構(gòu)滿(mǎn)足強(qiáng)度和剛度要求。

在空載情況下，水窖容積為20、40和60 m3時(shí)，窖體結(jié)構(gòu)最優(yōu)的幾何尺寸分別如下：厚度T=0.08、0.11、0.14 m；水窖半徑R=1.57、1.80、1.93 m；底部圓角半徑R1=0.13、0.17、0.25 m；材料用量分別為3.38、7.18、12.03 m3。

在滿(mǎn)載情況下，水窖容積為20、40和60 m3時(shí)，窖體結(jié)構(gòu)最優(yōu)的幾何尺寸分別如下：厚度T=0.06、0.08、0.10 m；水窖半徑R=1.65、2.07、2.47 m；底部圓角半徑R1=0.17、0.16、0.30 m；材料用量分別為2.06、4.53、7.56 m3。在相同窖容下，空載情況比滿(mǎn)載情況下的材料用量大。如果水窖埋置在地面以下使用時(shí)，2種工況均會(huì)出現(xiàn)，應(yīng)采用空載計(jì)算結(jié)果。

3.3 球形水窖不同容水量?jī)?yōu)化結(jié)果

通過(guò)ANSYS優(yōu)化得出的結(jié)果如表2所示。

從表2可以看出，在空載工況下容積為20、40、60 m3時(shí)，水窖受到的最大拉應(yīng)力（S1max）分別為883.56、789.23、327.98 kPa；在滿(mǎn)載工況下容積為20、40、60 m3時(shí)，水窖受到的最大拉應(yīng)力分別為740.15、906.04、888.52 kPa，都小于材料的許用拉應(yīng)力910 kPa。而在2種工況下容積為20、40、60 m3時(shí)水窖受到的最大壓應(yīng)力范圍543.30～1 700.60 kPa，遠(yuǎn)遠(yuǎn)小于材料的許用壓應(yīng)力7 200 kPa，同時(shí)窖體計(jì)算位移最大值范圍0.050～0.149 mm也略小于材料的極限應(yīng)變0.15 mm，因此滿(mǎn)足強(qiáng)度和剛度的要求。

在空載不利工況下，水窖容積為20、40和60 m3時(shí)，窖體結(jié)構(gòu)最優(yōu)的幾何尺寸分別如下：厚度T=0.037、0.048、0.073 m；水窖半徑R=1.68、2.12、2.43 m；材料用量分別為1.68、3.19、6.08 m3。

在滿(mǎn)載不利工況下，水窖容積為20、40、60 m3時(shí)，窖體結(jié)構(gòu)最優(yōu)的幾何尺寸分別如下：厚度T=0.026、0.034、0.046 m；水窖半徑R=1.68、2.12、2.43 m；材料用量分別為1.24、2.35、3.96 m3。在相同窖容下，空載情況比滿(mǎn)載情況下的材料用量大。如果水窖埋置在地面以下使用時(shí)，2種工況均會(huì)出現(xiàn)，應(yīng)采用空載計(jì)算結(jié)果。

3.4 瓶形水窖不同容水量?jī)?yōu)化結(jié)果

通過(guò)ANSYS優(yōu)化得出的結(jié)果如表3所示。

從表3可以看出，在空載工況下容積為20、40、60 m3時(shí)，水窖受到的最大拉應(yīng)力（S1max）分別為854.45、899.67、906.32 kPa；在滿(mǎn)載工況下容積為20、40、60 m3時(shí)，水窖受到的最大拉應(yīng)力分別為858.68、908.9、909.16 kPa，都接近材料的許用拉應(yīng)力910 kPa。而在2種工況下容積為20、40、60 m3時(shí)水窖受到的最大壓應(yīng)力范圍606.8～2 439.4 kPa，遠(yuǎn)遠(yuǎn)小于材料的許用壓應(yīng)力7 200 kPa，同時(shí)窖體計(jì)算位移最大值0.038～0.134 mm也小于材料的極限應(yīng)變0.15 mm，因此滿(mǎn)足強(qiáng)度和剛度的要求。

在空載不利工況下，水窖容積為20、40和60 m3時(shí)，窖體結(jié)構(gòu)最優(yōu)的幾何尺寸分別如下：厚度T=0.10、0.15、0.19 m；窖身半徑R=1.79、2.04、2.30 m；窖底半徑R1=1.69、1.83、2.05 m；窖身高度H=1.01、2.11、2.55 m；旱窖高度H1=2.50、2.92、3.38 m；材料用量分別為4.25、9.99和16.74 m3，這里由于瓶形水窖修筑材料為混凝土，所以旱窖部分也可以蓄水。

在滿(mǎn)載不利工況下，水窖容積為20、40、60 m3時(shí)，窖體結(jié)構(gòu)最優(yōu)的幾何尺寸分別如下：厚度T=0.08、0.09、0.12 m；窖身半徑R=1.45、1.55、1.79 m；窖底半徑R1=1.19、0.80、1.00 m；窖身高度H=2.65、7.49、7.91 m；旱窖高度H1=1.91、2.07、2.50 m；材料用量分別為3.35、6.99、12.21 m3。在相同窖容下，空載情況比滿(mǎn)載情況下的材料用量大。如果水窖埋置在地面以下使用時(shí)，2種工況均會(huì)出現(xiàn)，應(yīng)采用空載計(jì)算結(jié)果。

3.4 3種水窖優(yōu)化結(jié)果對(duì)比分析

從表1和表2圓柱形與球形水窖優(yōu)化結(jié)果對(duì)比可以看出，首先材料用量，不管在哪種工況下，相同容積時(shí)圓柱形水窖的材料用量都高于球形水窖。水窖處于地面下容積為20、40、60 m3時(shí)，采用球形水窖的幾何尺寸最優(yōu)設(shè)計(jì)相對(duì)圓柱形水窖來(lái)說(shuō)分別可節(jié)約材料101.19%、125.08%、97.86%；水窖放在地面上容積為20、40、60 m3時(shí)，采用球形水窖的幾何尺寸最優(yōu)設(shè)計(jì)相對(duì)圓柱形水窖來(lái)說(shuō)分別可節(jié)約材料66.13%、92.77%、90.91%。

而與圓柱形、球形水窖相比，瓶形水窖在同工況同容積時(shí)最優(yōu)設(shè)計(jì)結(jié)果材料用量最大。

從表1～3可以得出，在3種水窖設(shè)計(jì)中，在水窖滿(mǎn)足強(qiáng)度和剛度的要求下，球形水窖設(shè)計(jì)是最省材料的，但施工過(guò)程與另2種水窖相比要相對(duì)復(fù)雜點(diǎn)，如果用就地紅膠泥作為防滲材料，不考慮材料用量，那么應(yīng)優(yōu)先選擇圓柱形和瓶形水窖設(shè)計(jì)，因?yàn)樗鼈兪┕み^(guò)程簡(jiǎn)單，工程量小，由于瓶形水窖下半部分窖壁具有一定的邊坡比，所以結(jié)構(gòu)穩(wěn)定性相對(duì)于直立的圓柱形水窖來(lái)說(shuō)要強(qiáng)一些。

4 結(jié)論與討論

（1）對(duì)容水量為20、40、60 m3的圓柱形、球形、瓶形水窖在2種工況下進(jìn)行了結(jié)構(gòu)的尺寸優(yōu)化設(shè)計(jì)。以水窖造價(jià)最低為目標(biāo)函數(shù)，以許用應(yīng)力和變形量為約束條件，以水窖厚度、半徑、窖身高度等為決策變量，建立了優(yōu)化模型，并利用ANSYS軟件求解，得到不同容水量、不同工況下水窖優(yōu)化尺寸，該尺寸對(duì)水窖建設(shè)具有實(shí)用價(jià)值。

（2）通過(guò)ANSYS的優(yōu)化結(jié)果與對(duì)比分析，可以看出球形水窖設(shè)計(jì)最省材料，瓶形水窖在同工況同容積下相對(duì)于另外2種水窖設(shè)計(jì)形式需要材料最多，結(jié)構(gòu)穩(wěn)定性高于圓柱形水窖。如果不考慮建窖材料用量，比如陜北紅膠泥土水窖，用當(dāng)?shù)丶t膠泥作為水窖防滲材料，建窖應(yīng)優(yōu)先考慮瓶形水窖，因?yàn)樗慕Y(jié)構(gòu)穩(wěn)定性高于圓柱形水窖，雖然低于球形水窖，但它的施工過(guò)程比球形水窖簡(jiǎn)單，而且開(kāi)挖安全性要高于球形水窖。

（3）該研究在建立ANSYS優(yōu)化模型過(guò)程中，利用APDL宏語(yǔ)言編制程序，利用GUI語(yǔ)言定制可視化界面，實(shí)際上對(duì)Ansis進(jìn)行了二次開(kāi)發(fā)，需要輸入的土壤參數(shù)、材料參數(shù)可根據(jù)需要改變。所以，該計(jì)算方法具有通用性，也可用于新材料水窖的開(kāi)發(fā)。

參考文獻(xiàn)

[1] 許紅艷，何丙輝，李章成，等.我國(guó)黃土地區(qū)水窖的研究[J].水土保持學(xué)報(bào)，2004，18（2）：58-62.

[2] 王勖成.有限單元法[M].北京：清華大學(xué)出版社，2003：98-160.

[3] 王小民.隧道內(nèi)復(fù)合式路面有排水溝墊層應(yīng)力分析[J].山西建筑，2007，33（26）：315-316.

[4] 胡安焱，賀屹，于延鳳.渭北旱塬集雨工程類(lèi)型及可行性分析[J].中國(guó)農(nóng)村水利水電，2002（8）：17-18.

[5] 闞前華，譚長(zhǎng)建，張娟，等.ANSYS高級(jí)工程應(yīng)用實(shí)例分析與二次開(kāi)發(fā)[M].北京：電子工業(yè)出版社，2006：285-299.

[6] 范欽珊，施燮琴.工程力學(xué)[M].北京：高等教育出版社，1989：273-301.

[7] 崔云鵬，蔣定生.水土保持工程學(xué)[M].西安：陜西人民出版社，1998：29-39.

[8] 博弈創(chuàng)作室.APDL參數(shù)化有限元分析技術(shù)及其應(yīng)用實(shí)例[M].北京：中國(guó)水利水電出版社，2004：178-217.

[9] 婁宗科，黨進(jìn)謙，許永功，等.黃土地區(qū)水窯窖設(shè)計(jì)參數(shù)與施工技術(shù)研究[J].水土保持通報(bào)，2001，21（3）：14-16，19.