平面鋼閘門三維有限元線性分析的若干問題

許中武，蔡 偉，周建方

（河海大學機電工程學院，江蘇 常州 213022）

0 引 言

閘門是用來關閉、開啟或局部開啟水工建筑物中過水孔口的活動結構，其主要作用是控制水位、調節(jié)流量。作為水工建筑物的重要組成部分，閘門的安全性和適用性在很大程度上影響著整個水工建筑物的運行效果［1］，平面鋼閘門作為閘門中使用頻率最高的類型之一，其設計計算和安全運行也更值得工程設計人員關注。對于平面鋼閘門，其設計計算方法有三維有限元方法和傳統(tǒng)的平面體系計算方法，目前各大設計院的關于平面鋼閘門設計和復核計算方法仍以平面體系計算方法為主。平面體系計算方法先將平面鋼閘門簡化分解為若干基本構件（如面板、主梁、水平次梁、豎直次梁、頂梁和底梁等），再對每一構件進行設計計算［2］。由于平面體系計算方法忽略了各構件的協(xié)同作用，且其計算模型和荷載分配均依據(jù)工程經(jīng)驗，可能會造成某些構件的材料浪費或局部偏于危險的現(xiàn)象。

隨著三維有限元方法和相關軟件的成熟和完善，許多專家學者開始將三維有限元方法應用于平面鋼閘門結構計算。同時，美國的《水工鋼結構設計》［3］、中國的SL 74-2019《水利水電鋼閘門設計規(guī)范》［4］（以下簡稱規(guī)范）和NB 35055-2015《水電工程鋼閘門設計規(guī)范》［5］等相關設計規(guī)范均指出，除了使用傳統(tǒng)的平面體系計算方法外，三維有限元方法也可供平面鋼閘門設計計算使用。但工程設計或科研人員在使用三維有限元方法對平面鋼閘門進行分析計算時并沒有相應導則可供參考。不同人員對軟件操作和三維有限元方法掌握程度的差異，可能導致對同一閘門的計算結果不同，從而降低了計算結果的可信度。為此，亟需平面鋼閘門三維有限元分析導則指導設計人員進行相關三維有限元線性分析。目前，已有國內學者對閘門三維有限元分析原則進行了研究。張雪才［3，6］等分別從閘門面板和支臂結構的合理單元類型、單元尺寸或數(shù)量等方面展開分析研究，并得出了相應的分析原則。然而，關于閘門三維有限元分析原則的研究目前仍處于初步階段，并不完善。

本文以某水電站機組進水口快速平面鋼閘門為主要研究對象，結合相關設計資料，基于ANSYS有限元分析軟件，從有限元網(wǎng)格劃分和漏水孔簡化兩個方面系統(tǒng)研究平面鋼閘門三維有限元線性分析時的若干問題，以期相關結論可供工程設計人員在對此類平面鋼閘門進行設計計算時提供參考，同時也為我國相關平面鋼閘門設計規(guī)范的修訂和三維有限元線性分析原則的制定提供依據(jù)。

1 閘門概況

本文研究的水電站機組進水口快速閘門為平面鋼閘門，孔口尺寸為8.8 m×11.0 m（寬×高），底檻高程736.00 m，設計水位827.83 m，設計水頭按92.00 m 計算。該平面鋼閘門的門葉結構布置圖如圖1所示，為方便下文闡述相關計算數(shù)據(jù)，將該平面鋼閘門布置的12 根主梁從上到下分別記為1 號到12 號主梁。其中1 號和2 號主梁為箱形結構，記為上節(jié)箱形梁；同理11 號和12 號主梁構成下節(jié)箱形梁，其余主梁為工形梁；底梁為矩形截面梁，截面尺寸為40 mm×9 960 mm×200 mm。

圖1 門葉結構布置圖Fig.1 Structural drawing of door leaf

2 有限元網(wǎng)格劃分

劃分網(wǎng)格是建立有限元模型的一個重要環(huán)節(jié)，所劃分的網(wǎng)格形式對計算精度和計算規(guī)模將產(chǎn)生直接影響［7］，因此需要考慮多方面因素。在有限元軟件中網(wǎng)格劃分流程主要包括選定單元屬性、設置網(wǎng)格尺寸和選擇網(wǎng)格劃分方式等，其中單元屬性包括單元類型、材料類型和單元坐標系等。目前，在平面鋼閘門有限元模型的單元類型選定、網(wǎng)格尺寸設置和網(wǎng)格劃分方式選擇等方面的相關依據(jù)和研究較少，科研和工程設計人員多憑經(jīng)驗劃分網(wǎng)格。本節(jié)將從上述3個方面對平面鋼閘門有限元模型的網(wǎng)格劃分進行討論。

2.1 單元類型

對于平面鋼閘門這樣典型的空間薄壁結構，為保證計算精度，在建立其有限元模型時應盡量建立完整空間薄壁結構。此時，平面鋼閘門所有構件均采用殼單元劃分網(wǎng)格。合理的殼單元類型將決定結構應力和撓度的計算精度。文獻［6］以平面鋼閘門面板區(qū)格的幾種計算模型為例，分別驗證了ANSYS軟件中分析薄壁結構常用殼單元的計算精度。然而，文獻［6］僅比較了使用不同單元類型時計算模型撓度的有限元解和理論解，實際上，在平面鋼閘門結構分析時，保證其強度符合規(guī)范要求同樣重要。另一方面，文獻［6］中的計算模型僅討論了平面鋼閘門面板區(qū)格，而實際上平面鋼閘門主梁的簡支梁模型在設計計算中也很常見。

由于ANSYS 在近幾個版本已將shell43 和shell63 單元從圖形用戶界面操作選項中淘汰，因此本文將選用當前版本中常用的shell181 和shell281 單元。shell181 和shell281 單元均適合對薄的及具有一定厚度的殼體結構進行分析，shell181 單元使用一次插值函數(shù)，由4 個結點組成，每個結點具有6 個自由度；shell281 單元則使用二次位移插值函數(shù)，由8 個節(jié)點組成，各節(jié)點自由度與shell181相同。

平面鋼閘門的平面體系計算方法一般將主梁視為簡支梁模型，將面板區(qū)格視為四邊固定矩形彈性薄板模型，下面將使用這兩種常用計算模型，比較shell181和shell281這兩種單元類型在應力和變形上的計算精度。在ANSYS 中分別建立相應的計算模型，兩種模型的材料參數(shù)均依據(jù)Q355B 確定。其中，簡支梁截面尺寸為0.3 m×0.05 m（梁高×板厚），支承長度為10 m，均布荷載為10 kN/m；四邊固定矩形彈性薄板尺寸為2 m×1 m×0.05 m（板長×板寬×板厚），均布荷載為0.5 MPa。據(jù)此建立的計算模型如圖2所示。

圖2 有限元計算模型Fig.2 Finite element calculation model

對于受均布荷載的簡支梁，其應力和最大撓度可根據(jù)材料力學公式直接求得。對于受均布荷載的四邊固定矩形彈性薄板，其應力計算可依據(jù)規(guī)范和彈性力學理論［8］；在計算其撓度時，根據(jù)彈性力學理論可知，在長邊與短邊之比為2 時，矩形薄板的最大撓度理論解w為［9］：

式中：q為均布荷載；a為短邊長度；E為彈性模量，取2.1×106MPa；h為薄板厚度。

對受均布荷載的簡支梁和四邊固定矩形彈性薄板模型，在劃分網(wǎng)格建立有限元模型時，分別選用shell181 和shell281 單元，且單元尺寸設置等均相同，計算后提取相應結果，并與理論解比較，如表1和表2所示。

表1 簡支梁計算結果Tab.1 Calculation results of simply supported beam

表2 四邊固定矩形彈性薄板計算結果Tab.2 Calculation results of rectangular elastic thin plate with four edges fixed

由表1 和表2 可知，在荷載、約束和網(wǎng)格尺寸等其他設置相同的情況下，對于簡支梁和四邊固定矩形彈性薄板模型，shell281 單元無論在應力還是變形計算上精度都更高。以簡支梁最大正應力為例，使用shell181 單元的有限元解相對于理論解誤差為-20.23%，而shell281 單元有限元解誤差僅為-0.06%。特別地，從四邊固定矩形彈性薄板的計算結果可知，選用兩種不同的單元類型對撓度計算結果影響不大，與文獻［6］的結論一致，說明了上述計算結果的正確性。但是比較應力計算結果的理論值和有限元解，shell281 單元的精度顯著提高，相差一個量級，計算時間相差不大。因此對平面鋼閘門這樣的空間薄壁結構進行有限元計算時應優(yōu)先選取shell281單元。

2.2 網(wǎng)格尺寸

有限元模型網(wǎng)格尺寸對計算精度影響較大，模型網(wǎng)格尺寸的減小在一定程度上會提高計算精度，但當網(wǎng)格尺寸減小到一定程度時，精度提高會很有限，而計算時耗卻會大大增加。此時，可以認為網(wǎng)格尺寸對計算精度的影響可以忽略不計，即有限元計算結果與網(wǎng)格尺寸之間具有無關性。在平面鋼閘門有限元計算時，獲取網(wǎng)格無關性的解是工程設計或科研人員需要關注的關鍵問題。本文將選取shell281 單元，計算分析不同網(wǎng)格尺寸下的平面鋼閘門有限元模型，探討適用于本文所研究的平面鋼閘門的合理網(wǎng)格尺寸范圍。

分別建立7 種不同網(wǎng)格尺寸的平面鋼閘門有限元模型，各模型的網(wǎng)格尺寸如表3 所示。本文在建立有限元模型時，使用ANSYS 默認的全局直角坐標系（笛卡爾全局坐標系），其中x軸沿平面鋼閘門跨度方向，y軸沿平面鋼閘門高度方向，z軸沿水流方向，坐標原點位于平面鋼閘門豎板底部中點處。

表3 模型尺寸分組Tab.3 Model size grouping

該平面鋼閘門采用Q355B 材料，計算時材料彈性模量E取2.06×105MPa，密度ρ取7 850 kg/m3，泊松比v取0.3。平面鋼閘門作用水頭按設計水頭設置。此外，1 號主梁腹板承受額外的垂直水柱壓力。根據(jù)平面鋼閘門實際工作狀況，滑塊上受水流方向約束（Uz=0），平面鋼閘門底部受鉛直方向約束（Uy=0）。為保證計算模型的幾何不變性，假定平面鋼閘門底部中點沿跨度方向無位移（Ux=0）。下文所有平面鋼閘門有限元計算的荷載、相關參數(shù)和約束設置均相同。

以網(wǎng)格尺寸為0.40 m×0.40 m 的模型1 為例，計算后得到的整體等效應力和變形云圖如圖3 所示，提取所有平面鋼閘門模型中主要構件的最大等效應力和撓度列于表4和表5。

表5 平面鋼閘門主要構件最大撓度mTab.5 Maximum deflection of main components of plane steel gate

圖3 模型1有限元計算結果Fig.3 Finite element calculation results of Model 1

由表4和表5可知：①隨著網(wǎng)格尺寸的減小，平面鋼閘門面板和各主梁的最大等效應力呈現(xiàn)增大趨勢，最大等效應力出現(xiàn)位置均相同，但網(wǎng)格尺寸對面板最大等效應力影響最大。對于面板，從模型1 到模型4，即網(wǎng)格尺寸從0.40 m×0.40 m 到0.25 m×0.25 m 時，面板最大等效應力增長平緩；但從模型4 到模型5，即網(wǎng)格尺寸從0.25 m×0.25 m 到0.20 m×0.20 m 時，其最大等效應力增長率最大，達到7.56%；從模型5 到模型7，即網(wǎng)格尺寸小于0.20 m×0.20 m 時，面板最大等效應力繼續(xù)增長平緩，趨于收斂。此時，可以認為網(wǎng)格尺寸和平面鋼閘門面板最大等效應力計算結果具有無關性。②網(wǎng)格尺寸對平面鋼閘門主要構件撓度影響很小。各模型面板和主梁最大撓度偏差較小。以面板為例，隨著網(wǎng)格尺寸的減小，各模型最大撓度增長率最大僅為0.39%，即在網(wǎng)格尺寸小于0.40 m×0.40 m 時，網(wǎng)格尺寸和平面鋼閘門主要構件最大撓度計算結果即具有網(wǎng)格無關性。綜上，在對本文所研究的平面鋼閘門進行有限元計算時，其合理網(wǎng)格尺寸應小于或等于0.20 m×0.20 m。

2.3 網(wǎng)格劃分方式

對于面，ANSYS 為用戶提供了兩種網(wǎng)格劃分方式：自由網(wǎng)格劃分和映射網(wǎng)格劃分。自由網(wǎng)格劃分通常由有限元計算軟件自動劃分，易出現(xiàn)畸變網(wǎng)格和退化網(wǎng)格［10］。而映射網(wǎng)格劃分則需滿足特定規(guī)則，具有規(guī)則形狀且明顯成排。一般映射劃分得到的網(wǎng)格質量優(yōu)于自由劃分得到的網(wǎng)格質量，映射劃分有限元模型的計算精度理論上也優(yōu)于自由劃分有限元模型的計算精度。將映射劃分應用于平面鋼閘門有限元模型的建立上，需要在前期投入大量工作，對某些復雜的面需要進行切割，對各個面也只能逐一劃分網(wǎng)格。如果映射劃分方法對平面鋼閘門在應力和撓度的計算精度上提高甚微，則會得不償失。

本文使用shell281單元，網(wǎng)格尺寸設置為0.10 m×0.10 m，分別對該平面鋼閘門進行自由網(wǎng)格劃分和映射網(wǎng)格劃分，生成的局部有限元模型如圖4 所示，其中自由網(wǎng)格劃分得到的有限元模型出現(xiàn)了大量退化網(wǎng)格。兩種網(wǎng)格劃分方式所得有限元模型的單元數(shù)和節(jié)點數(shù)如表6所示。

表6 平面鋼閘門有限元模型節(jié)點和單元數(shù)Tab.6 Node and element number of plane steel gate finite element model

圖4 平面鋼閘門局部有限元模型Fig.4 Local finite element model of plane steel gate

經(jīng)有限元計算后提取平面鋼閘門主要構件的最大等效應力和撓度列于表7。

表7 平面鋼閘門主要構件應力和撓度計算結果Tab.7 The stress and deflection calculation results of the main components of plane steel gate

兩種網(wǎng)格劃分方式下，計算得到的平面鋼閘門整體的應力和變形分布均類似，最大等效應力和最大變形值偏差較小。由表7可知：與映射劃分模型相比，平面鋼閘門各主梁和面板在自由劃分下的最大等效應力值誤差僅為-3.38%，最大撓度誤差不超過1%，即無論是應力還是變形，兩種劃分方式下的計算結果相對誤差都很小。

但實際上，在使用映射劃分方式時，平面鋼閘門有限元模型的單元數(shù)和節(jié)點數(shù)都更多，計算規(guī)模更大，并且映射劃分在前期建立模型階段還需要花費大量時間和精力。此外，雖然通過映射劃分方式得到的網(wǎng)格質量更佳，但從計算結果來看精度提高并不明顯。因此使用有限元方法對此類平面鋼閘門進行計算分析時采用自由劃分方式即可。

3 漏水孔結構簡化

由于平面鋼閘門結構復雜，構件眾多，細節(jié)特征也較多，因此在有限元軟件ANSYS中建立其模型時通常需要適當簡化，使模型與計算機性能和分析目的相匹配。簡化平面鋼閘門模型一方面可以減小有限元求解的規(guī)模，提高計算效率，另一方面也可以降低工程設計或科研人員在前期的建模工作量。

為了防止平面鋼閘門各主梁上的積水影響啟閉甚至導致其構件銹蝕，進而影響平面鋼閘門結構性能，主梁腹板上一般會開設有漏水孔。如果漏水孔對平面鋼閘門應力或變形計算影響不大，則可在建模階段忽略這些漏水孔，減小有限元分析工作量，提高效率。

本文研究的平面鋼閘門在主梁腹板上開設有若干直徑為200 mm的漏水孔，其中2號至11號各節(jié)主梁腹板分別均勻布置6 個，12 號主梁腹板均勻布置8 個。根據(jù)規(guī)范，一般主梁腹板應盡量避免開大孔，漏水孔直徑Ф與主梁截面高度（下文簡稱“梁高”）H之比（下文簡稱“徑高比”）應小于0.5，必要時還需對漏水孔采取補強措施，保證結構強度。因此漏水孔的孔徑尺寸是判斷其是否可以簡化的重要參數(shù)。

本文將針對該平面鋼閘門，在保持漏水孔位置不變的情況下，分別建立徑高比為0.1、0.15、0.2、0.25和0.3的平面鋼閘門模型，相應的孔徑和梁高參數(shù)設置如表8所示。同時，還建立了無漏水孔的平面鋼閘門模型，用以與相應的計算結果對比。開展有限元計算并提取平面鋼閘門主梁應力和撓度的計算結果，統(tǒng)計繪制出圖5和圖6。

表8 孔徑和梁高設置Tab.8 Setting of aperture and beam height

圖5 各模型主梁最大等效應力Fig.5 Maximum equivalent stress of main beam for different models

圖6 各模型主梁最大撓度Fig.6 Maximum deflection of main beam for different models

由圖5 和圖6 可知：①在徑高比為0.1 到0.15 時，各模型主梁最大等效應力值之間相差很小，出現(xiàn)的位置相同，且各組計算結果與簡化模型計算結果基本相同；在徑高比為0.2 到0.3時，各計算組主梁在兩側漏水孔附近均出現(xiàn)了較為明顯的應力集中現(xiàn)象，且奇點應力值大于簡化模型主梁最大等效應力，并隨徑高比的增大而增大，部分主梁兩側漏水孔附近奇點應力超出相應的許用應力；②不同徑高比模型的主梁最大撓度均出現(xiàn)在跨中位置，與簡化模型計算結果發(fā)生的位置相同。隨著徑高比的增大，主梁最大撓度相對偏差逐漸增大。最大偏差出現(xiàn)在徑高比為0.3模型的下節(jié)箱形梁，為10.12%，但各模型主梁撓度都滿足規(guī)范設計要求。

綜上，當徑高比小于或等于0.15 時，在建立平面鋼閘門有限元模型時可忽略主梁腹板漏水孔，簡化模型；而當徑高比大于0.15 時不可忽略漏水孔，并且需要根據(jù)有限元計算結果采取相應的補強措施。

4 結 論

為了使工程設計和科研人員在使用三維有限元方法對平面鋼閘門進行有限元線性分析時有據(jù)可依，且更加科學高效，本文從網(wǎng)格劃分和模型簡化等方面展開對平面鋼閘門有限元計算過程的若干問題進行了討論研究，主要得到了以下結論。

（1）本文以平面鋼閘門平面體系計算方法中常用的簡支梁和四邊固定矩形彈性薄板模型為算例，討論了常用殼單元的計算精度。結果表明，在同樣網(wǎng)格尺寸下，采用shell281單元的應力計算精度顯著提高，相差約一個量級，而計算時間并無明顯增加，對平面鋼閘門進行有限元計算時宜選用該單元。

（2）針對本文所研究的平面鋼閘門，在劃分網(wǎng)格建立有限元模型時，合理的網(wǎng)格尺寸設置應小于或等于0.20 m×0.20 m；映射劃分相對于自由劃分得到的網(wǎng)格質量更高，但兩種劃分方式下閘門主要構件的應力和變形計算結果相對偏差很小，鑒于映射劃分方式前期工作量較大，為提高效率，采用自由劃分即可。

（3）在對平面鋼閘門使用三維有限元方法進行計算分析時，當徑高比小于或等于0.15 時，在建立平面鋼閘門模型時可忽略漏水孔；在徑高比大于0.15 時，由于部分主梁兩側漏水孔邊緣出現(xiàn)應力集中現(xiàn)象，且奇點應力大于許用應力，此時應建立平面鋼閘門完整模型，并需要根據(jù)有限元計算結果對漏水孔采取相應的補強措施。