圓拱直墻式水窖應(yīng)力集中現(xiàn)象研究

張仕華

摘要 圓拱直墻式水窖是我國(guó)西北缺水地區(qū)以及水資源分布不均衡地區(qū)的一種常用集雨設(shè)施，由于目前在水窖建設(shè)中缺乏對(duì)圓拱直墻式水窖的空間結(jié)構(gòu)分析，忽略了該水窖存在的應(yīng)力集中問(wèn)題，導(dǎo)致水窖在應(yīng)力集中部位的破壞或整個(gè)水窖采用較大厚度導(dǎo)致材料的浪費(fèi)。通過(guò)利用大型有限元分析軟件ANSYS對(duì)圓拱直墻式水窖進(jìn)行空間結(jié)構(gòu)分析，計(jì)算了該水窖的應(yīng)力分布規(guī)律，從而驗(yàn)證水窖強(qiáng)度。并通過(guò)ANSYS優(yōu)化設(shè)計(jì)理論，以水窖耗材量最小為目標(biāo)函數(shù)，以許用應(yīng)力為約束條件，得出不同容水量的水窖優(yōu)化幾何尺寸。根據(jù)該尺寸，工程施工人員可以方便地完成該型式水窖的施工，并且為新材料水窖的開發(fā)提供一定的理論基礎(chǔ)。

關(guān)鍵詞 圓拱直墻式水窖；有限元；應(yīng)力集中

中圖分類號(hào) S277 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2018）21-0185-04

Abstract The cylindrical water cellar is a common rainfall collection facility in northwest water shortage areas of China and water resources distributed unbalanced areas. The lack of analysis on water cellar spatial structure and ignoring water cellar stress concentration problems have led to damage on the water cellar stress concentration sites or a greater thickness of the waste material. In this research， we calculated the stress distribution of the water cellar to verify the water cellar strength through analyzing the water cellar structure by ANSYS software， and obtained the different volumes of water cellar geometry optimization by ANSYS optimal designing theory， so as to use the minimum water cellar materials as the objective function， and to use allowable stress of the material as the constrain condition.Then， according to the optimal sizes， the worker could build this kind of water cellar conveniently， and this research provided theoretical basis for the development of new water cellar.

Key words Circular arch and straight wall water cellar；Finite element；Stress concentration

近年來(lái)，隨著社會(huì)的發(fā)展和人口的劇增，淡水資源緊缺的程度也越來(lái)越嚴(yán)重，水窖作為收集雨水的一種典型設(shè)施得到越來(lái)越廣范的應(yīng)用，圓拱直墻式水窖作為我國(guó)常用集雨設(shè)施形式[1]，以其施工簡(jiǎn)單快捷、工程量小、窖體結(jié)構(gòu)穩(wěn)定性好等優(yōu)點(diǎn)受到廣泛建窖者與農(nóng)戶的青睞。但現(xiàn)行的水窖都是建窖者根據(jù)自己的修筑經(jīng)驗(yàn)來(lái)建造，往往在水窖使用一段時(shí)間后窖壁底部便出現(xiàn)斷裂，使水窖早早報(bào)廢，其根本原因是窖底局部應(yīng)力大于材料的許用應(yīng)力而導(dǎo)致強(qiáng)度破壞，所以，有必要對(duì)該類型水窖進(jìn)行強(qiáng)度驗(yàn)算并確定出合理的尺寸。由于水窖結(jié)構(gòu)的三維尺寸差異較小，所以在結(jié)構(gòu)計(jì)算中必須以空間結(jié)構(gòu)原理為基礎(chǔ)，計(jì)算工作量大而且復(fù)雜，必須借助于適當(dāng)?shù)挠?jì)算機(jī)軟件才能完成，ANSYS軟件的出現(xiàn)為解決這一問(wèn)題提供了有力的支持。另一方面，我國(guó)的水窖使用者和修建者是農(nóng)民，由于缺乏理論指導(dǎo)，建成的水窖或者由于局部應(yīng)力過(guò)大導(dǎo)致破壞，或者壁厚過(guò)大造成材料浪費(fèi)。截至目前，尚鮮有對(duì)水窖進(jìn)行空間結(jié)構(gòu)分析的相關(guān)研究。因此，筆者借助于大型有限元分析軟件ANSYS對(duì)圓拱直墻式水窖進(jìn)行結(jié)構(gòu)分析，對(duì)水窖底部采用直角和圓角2種情況進(jìn)行應(yīng)力比較，通過(guò)對(duì)結(jié)構(gòu)尺寸的優(yōu)化，設(shè)計(jì)出既安全又經(jīng)濟(jì)的圓拱直墻式水窖，為該水窖的建設(shè)者提供技術(shù)支撐，并為新型水窖的開發(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]。

該研究建立2種模型：模型1為圓拱直墻式水窖底部直角設(shè)計(jì)，模型2為圓拱直墻式水窖底部圓角設(shè)計(jì)。2種水窖斷面形狀如圖1所示，它們的總體結(jié)構(gòu)尺寸設(shè)計(jì)一致，由窖頸、窖身和窖底3部分組成。圖1左側(cè)為窖底與窖壁結(jié)合處為直角的水窖，右側(cè)為窖底與窖壁用半徑為R1的圓角連接。

實(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)格化有限元模型[4]，由于圓柱形水窖形狀和荷載均為對(duì)稱實(shí)體，所以，分析時(shí)只對(duì)1/4窖體進(jìn)行結(jié)構(gòu)分析。底部圓角設(shè)計(jì)水窖網(wǎng)格劃分如圖2所示。

2 水窖最不利工作情況及荷載

3 圓角和直角應(yīng)力比較

為了比較水窖底部采用直角和圓角時(shí)的應(yīng)力分布狀況，利用有限元法建立模型，然后求解，求解時(shí)需要輸入的參數(shù)如下：

水窖容積20 m3，混凝土抗拉強(qiáng)度910 kPa，混凝土抗壓強(qiáng)度7 200 kPa，混凝土極限拉應(yīng)變0.000 15 m，混凝土極限壓應(yīng)變0.003 3 m，土體容重16 kN/m3，厚度T 0.11 m，彈性模量2.2×107 kPa，泊松比0.2。

滿載情況下直角和圓角的第1主應(yīng)力[6]云圖如圖3所示。

從圖3可以看出，直角水窖滿載時(shí)，其最大應(yīng)力為779.89 kPa，集中于窖底與窖壁的結(jié)合處；而圓角水窖滿載時(shí)，其最大應(yīng)力為423.50 kPa，分散分布于圓角壁的上部，可見(jiàn)水窖圓角設(shè)計(jì)可將集中的大應(yīng)力減小分散開。空載和滿載情況下，徑向（S1）、環(huán)向（S2）和垂直向（S3）最大主應(yīng)力如表1所示，其中正值為拉應(yīng)力，負(fù)值為壓應(yīng)力。

混凝土材料可視為各向同性材料，其安全性由3種主應(yīng)力中的最大應(yīng)力控制。由表1可以看出，當(dāng)為空載工況時(shí)，水窖底部為圓角設(shè)計(jì)時(shí)受到的最大拉應(yīng)力549.88 kPa，其值小于直角設(shè)計(jì)水窖受到的最大拉應(yīng)力734.18 kPa；當(dāng)為滿載工況時(shí)，圓角設(shè)計(jì)水窖受到的最大拉應(yīng)力423.48 kPa，其值亦小于直角設(shè)計(jì)水窖受到的最大拉應(yīng)力779.88 kPa，這就體現(xiàn)了窖體底部圓角設(shè)計(jì)可以有效緩解應(yīng)力集中的情況。

從表1還可以看出，空載時(shí)，直角水窖受到的最大壓應(yīng)力2 008.00 kPa，圓角水窖受到的最大壓應(yīng)力為1 221.00 kPa；滿載時(shí)直角水窖受到的最大壓應(yīng)力為337.05 kPa，圓角水窖受到的最大壓應(yīng)力為378.21 kPa，都遠(yuǎn)遠(yuǎn)小于混凝土設(shè)計(jì)抗壓強(qiáng)度7 200 kPa。由此說(shuō)明，利用混凝土作為建窖材料，引起了壓應(yīng)力的浪費(fèi)，所以從力學(xué)角度講，混凝土并不是良好的建窖材料。實(shí)際拉應(yīng)力略小于許用拉應(yīng)力，所以，水窖將會(huì)由于拉應(yīng)力不足而破壞。

4 圓拱直墻式水窖幾何尺寸優(yōu)化

為了在滿足強(qiáng)度和剛度要求的條件下，減小材料用量，降低水窖成本，以水窖造價(jià)最低為目標(biāo)函數(shù)，以應(yīng)力和變形為約束條件，以水窖結(jié)構(gòu)尺寸決策變量，建立優(yōu)化模型。

4.1 建立優(yōu)化模型

根據(jù)《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》C15混凝土極限伸長(zhǎng)量不超過(guò)0.003 3 m，而極限拉壓縮量不超過(guò)0.000 15 m，窖壁厚度T的優(yōu)化取值范圍按以前修筑水窖窖壁厚度范圍設(shè)定為6～20 cm。圓柱半徑R的優(yōu)化取值范圍設(shè)定也按以往修筑水窖經(jīng)驗(yàn)[7]，要求為首先窖體容積固定，R在0.2H～2.0H取值，在窖容為20、40和60 m3時(shí)R取值分別為1.039 6～1.760 5、1.309 8～2.218 1和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）、用戶優(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ù)目，該實(shí)驗(yàn)將兩者設(shè)置為10和7，也可根據(jù)需要增大設(shè)計(jì)序列的數(shù)目，一般初始設(shè)計(jì)序列數(shù)目30，最大可行性設(shè)計(jì)序列10都能滿足一般優(yōu)化要求，即得到較為精確的解。

4.2 等厚度水窖不同容水量?jī)?yōu)化結(jié)果

通過(guò)ANSYS優(yōu)化得出的結(jié)果見(jiàn)表2。

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

空載時(shí)材料用量從窖容20 m3的3.38 m3增加到60 m3的12.03 m3；滿載時(shí)材料用量從20 m3的2.06 m3增加到60 m3的7.56 m3。

在空載不利工況下，水窖容積為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。在滿載不利工況下，水窖容積為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。在相同窖容下，空載情況比滿載情況下的材料用量大。所以，如果水窖埋置在地面以下使用時(shí)，2種工況均會(huì)出現(xiàn)，應(yīng)采用空載計(jì)算結(jié)果。

4.3 變厚度水窖不同容水量?jī)?yōu)化結(jié)果

以上在分析水窖時(shí)，窖體各部分厚度都一樣，而水窖的受力狀況是上部小，下部大，因此在窖體上部應(yīng)力小的部位采用小的厚度，在下部應(yīng)力大的部位采用大的厚度，將進(jìn)一步節(jié)約材料用量，變厚度窖體結(jié)構(gòu)如圖4所示。

ANSYS優(yōu)化分析結(jié)果見(jiàn)表3。從表3可以看出，水窖上部圓拱以及窖頸部分受力較小，所以這部分厚度T1較小，而窖身底部受力較大，這部分厚度T2較大，窖身部分厚度在T1到T2之間呈線性變化。從表2和表3可以看出，空載情況下，變厚度水窖與等厚度水窖相比，20、40和60 m3水窖材料用量分別減少了24%、28%和46%；在滿載情況下，材料用量變化不大。因此在滿足強(qiáng)度安全的前提下采用窖壁厚度漸變與其受力大小相對(duì)應(yīng)是一種更為理想的水窖結(jié)構(gòu)設(shè)計(jì)模式，這樣使窖體各部分承受應(yīng)力均勻，更有效地消除了水窖受力集中問(wèn)題。

5 結(jié)論與討論

（1）利用有限元分析方法對(duì)圓拱直墻式水窖的2種形式（底部采用直角和圓角）的結(jié)構(gòu)分析發(fā)現(xiàn)，當(dāng)直墻底部采用圓角時(shí)，應(yīng)力分布趨于均勻，采用直角時(shí)，底部應(yīng)力集中現(xiàn)象明顯，而且在空載和滿載2種情況下，采用圓角的最大應(yīng)力遠(yuǎn)遠(yuǎn)小于采用直角的最大應(yīng)力。

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

（3）從應(yīng)力計(jì)算結(jié)果看出，水窖許用拉應(yīng)力是控制條件，當(dāng)滿足拉應(yīng)力要求時(shí)，實(shí)際承受的壓應(yīng)力遠(yuǎn)遠(yuǎn)小于許用壓應(yīng)力。由此說(shuō)明，利用混凝土作為建窖材料，引起了壓應(yīng)力的浪費(fèi)，從力學(xué)角度考慮，混凝土并不是理想的建窖材料，所以有必要開發(fā)新材料水窖。

（4）該研究在建立ANSYS優(yōu)化模型過(guò)程中，利用APDL宏語(yǔ)言編制程序，利用GUI語(yǔ)言定制可視化界面，實(shí)際上對(duì)ANSYS進(jìn)行了二次開發(fā)，需要輸入的土壤參數(shù)、材料參數(shù)可根據(jù)需要改變。所以，該計(jì)算方法具有通用性，也可用于新材料水窖的開發(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] 闞前華，譚長(zhǎng)建，張娟，等.ANSYS高級(jí)工程應(yīng)用實(shí)例分析與二次開發(fā)[M].北京：電子工業(yè)出版社，2006：285-299.

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

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

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

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