三邊固支一邊自由矩形薄板撓度計(jì)算研究

郭翔詠

摘? 要：為求解受均布荷載作用下三邊固支一邊自由矩形薄板的撓度函數(shù)，在彈性理論的基礎(chǔ)上，建立薄板力學(xué)模型。結(jié)合三邊固支一邊自由的邊界條件，并采用Rayleigh-Ritz法計(jì)算撓度函數(shù)中的未知系數(shù)k，從而推導(dǎo)撓度函數(shù)。最后，通過(guò)ANSYS對(duì)2個(gè)算例進(jìn)行模擬，得到矩形薄板在z軸方向上的變形云圖，并對(duì)比理論計(jì)算與模擬的結(jié)果。結(jié)果的一致性為計(jì)算三邊固支一邊自由等復(fù)雜邊界條件的矩形薄板的撓度提供科學(xué)依據(jù)。

關(guān)鍵詞：三邊固支一邊自由；矩形薄板；撓度；Rayleigh-Ritz法；ANSYS

中圖分類(lèi)號(hào)：TB122? ? ? ?文獻(xiàn)標(biāo)志碼：A? ? ? ? ? 文章編號(hào)：2095-2945（2024）10-0079-04

Abstract： In order to obtain the deflection function of the rectangular thin plate with three fixed sides and one free side under the uniformly distributed load， the mechanical model for the thin plate was established based on the theory of elasticity. By considering the boundary conditions of three fixed sides and one free side， and employing the Rayleigh-Ritz method， the unknown coefficient k in the deflection function was calculated， thus deriving the deflection function. Finally， two cases were simulated using ANSYS to obtain the deformation contour plot of the rectangular thin plate in the z-axis direction， and the comparison was made between the theoretical calculations and the simulation results. The consistency of the results provides a scientific basis for calculating the deflection of the rectangular thin plate under complex boundary conditions， such as three fixed sides and one free side.

Keywords： three fixed sides and one free side; rectangular thin plate; deflection; Rayleigh-Ritz method; ANSYS

彈性薄板在工程中有廣泛的應(yīng)用，例如實(shí)體板和板梁等，而矩形薄板是其中較為常見(jiàn)的一種應(yīng)用形式[1-2]。邊界條件對(duì)薄板的穩(wěn)定性和整體行為起著關(guān)鍵作用。通過(guò)適當(dāng)選擇和施加邊界條件，可以預(yù)防薄板在受載情況下出現(xiàn)不穩(wěn)定、失穩(wěn)或振動(dòng)等問(wèn)題。目前，對(duì)于邊界條件相對(duì)簡(jiǎn)單的矩形薄板撓度函數(shù)已經(jīng)有了有效的求解方法，例如四邊固支和四邊簡(jiǎn)支矩形薄板的撓度函數(shù)。馬文強(qiáng)等[3]建立了綜放采場(chǎng)頂煤與頂板在均布載荷作用下的薄板模型，并利用Matlab編譯了薄板的撓曲面公式。馬仁香[4]基于利維解和疊加法推導(dǎo)出了均布載荷作用下的四邊固支矩形薄板的撓度表達(dá)式，并結(jié)合有限元模擬驗(yàn)證了該表達(dá)式的正確性。然而，對(duì)于求解具有復(fù)雜邊界條件的矩形薄板撓度函數(shù)而言，由于求解的困難和復(fù)雜性，此類(lèi)研究還較為有限。因此，根據(jù)已有的薄板復(fù)雜邊界條件，假設(shè)薄板的撓度函數(shù)并選擇有效的方法對(duì)其進(jìn)行求解，是工程領(lǐng)域亟待解決的問(wèn)題。

Rayleigh-Ritz法[5]用于求解薄板的撓度函數(shù)時(shí)，只需預(yù)先假設(shè)的撓度函數(shù)滿足薄板的位移邊界條件，無(wú)須滿足內(nèi)力邊界條件，這極大地簡(jiǎn)化了求解復(fù)雜邊界條件下薄板撓度函數(shù)的過(guò)程?；诖?，本文建立了薄板力學(xué)模型，使用Rayleigh-Ritz法求解三邊固支一邊自由矩形薄板的撓度函數(shù)，并利用ANSYS有限元軟件驗(yàn)證了理論計(jì)算結(jié)果的正確性。

1? 三邊固支一邊自由矩形薄板撓度函數(shù)求解

1.1? 薄板力學(xué)模型建立

一般說(shuō)來(lái)，當(dāng)矩形板的厚度與板中面最小尺寸比值在1/80和1/5之間時(shí)，稱(chēng)為薄板[6]。以o為中心建立三維直角坐標(biāo)系，薄板長(zhǎng)度為a，寬度為b，厚度為h，作用在其上表面的均布荷載為q，薄板力學(xué)模型如圖1（a）所示。三邊固支一邊自由邊界條件如圖1（b）所示。

（b）? 三邊固支一邊自由邊界條件

圖1? 薄板力學(xué)模型建立

1.2? 薄板撓度函數(shù)求解

薄板的彈性曲面微分方程為[7]

D？犖4w=q， （1）

式中：w為薄板的撓度函數(shù)；？犖4為雙調(diào)和算子，？犖4=■+2■+■；D為薄板的彎曲剛度，D=■，其中，E為薄板的彈性模量，v為薄板的泊松比。

三邊固支一邊自由薄板的邊界條件為

。（2）

取滿足位移邊界條件的撓度函數(shù)

式中：k為撓度函數(shù)的未知系數(shù)。

等厚度薄板的應(yīng)變能V？著為

依據(jù)Rayleigh-Ritz法，有

由式（5）可求得薄板撓度函數(shù)的未知系數(shù)k

將式（6）代入式（3）中可得薄板撓度函數(shù)w為

2? 模擬驗(yàn)證

2.1? 算例1

取三邊固支一邊自由矩形薄板的長(zhǎng)度a為60 m，寬度b為10 m，厚度h為1 m，作用在其上表面的均布荷載q為2 MPa，彈性模量E為20 GPa，泊松比v為0.3。由式（7）可得薄板的撓度函數(shù)圖，如圖2所示。

撓度的正負(fù)號(hào)代表方向。從圖2中可以看出，受均布荷載作用的三邊固支一邊自由矩形薄板的最大撓度出現(xiàn)在薄板自由邊的中點(diǎn)處，其值為37.194 19 mm。利用ANSYS有限元軟件建立矩形薄板彈性模型，薄板的彈性模量E為20 GPa，泊松比v為0.3。模擬計(jì)算的z軸定向變形云圖如圖3所示。

從圖3中可以看出，受均布荷載作用的三邊固支一邊自由矩形薄板的最大撓度同樣出現(xiàn)在薄板自由邊的中點(diǎn)處，與理論分析一致，其值為36.349 08 mm。本文的理論計(jì)算結(jié)果與模擬結(jié)果的相對(duì)誤差僅為2.324 9%。通常情況下，理論計(jì)算結(jié)果與模擬結(jié)果的相對(duì)誤差可以允許在5%以?xún)?nèi)[8]，因此在算例1中，理論計(jì)算得到的矩形薄板撓度函數(shù)與ANSYS的模擬結(jié)果是相符的。

2.2? 算例2

取三邊固支一邊自由矩形薄板的長(zhǎng)度a為50 m，寬度b為20 m，厚度h為2 m，作用在其上表面的均布荷載q為3 MPa，彈性模量E為20 GPa，泊松比v為0.3。由式（7）可得薄板的撓度函數(shù)圖，如圖4所示。

從圖4中可以看出，受均布荷載作用的三邊固支一邊自由矩形薄板的最大撓度出現(xiàn)在薄板自由邊的中點(diǎn)處，其值為109.022 8 mm。利用ANSYS有限元軟件建立矩形薄板彈性模型，薄板的彈性模量E為20 GPa，泊松比v為0.3。模擬計(jì)算的z軸定向變形云圖如圖5所示。

從圖5中可以看出，受均布荷載作用的三邊固支一邊自由矩形薄板的最大撓度值為109.039 2 mm。本文的理論計(jì)算結(jié)果與模擬結(jié)果的相對(duì)誤差僅為0.015 0%。對(duì)于算例2，理論計(jì)算得到的矩形薄板撓度函數(shù)同樣與ANSYS的模擬結(jié)果相吻合。這兩個(gè)算例的計(jì)算結(jié)果與模擬結(jié)果的相對(duì)誤差均在5%以?xún)?nèi)，滿足工程需求。

3? 結(jié)論

本文通過(guò)建立薄板力學(xué)模型，利用Rayleigh-Ritz法求解了三邊固支一邊自由矩形薄板的撓度函數(shù)。同時(shí)，利用ANSYS有限元軟件建立了矩形薄板的彈性模型。算例1和算例2的理論計(jì)算結(jié)果與模擬結(jié)果的相對(duì)誤差分別為2.324 9%和0.015 0%，均在5%以?xún)?nèi)，滿足工程實(shí)際要求，驗(yàn)證了理論公式的正確性。Rayleigh-Ritz法在求解具有復(fù)雜邊界條件的矩形薄板撓度函數(shù)方面具有良好的適用性。此外，該方法具有公式簡(jiǎn)單、便于工程實(shí)踐的特點(diǎn)，對(duì)工程實(shí)際具有借鑒和指導(dǎo)意義。

參考文獻(xiàn)：

[1] 孫建，胡洋.均布和靜水壓力作用下固支矩形薄板力學(xué)特性[J].應(yīng)用力學(xué)報(bào)，2015，32（6）：908-914，1096-1097.

[2] 王霂，王洪波.矩形薄板受橫向均布載荷作用彈性小撓度解[J].科技創(chuàng)新與應(yīng)用，2021（3）：82-84.

[3] 馬文強(qiáng)，王同旭，曲孔典.基于薄板模型的高韌性煤層難冒放機(jī)理分析[J].采礦與安全工程學(xué)報(bào)，2017，34（4）：644-649.

[4] 馬仁香.均布載荷下四邊固支矩形薄板的撓度[J].計(jì)算機(jī)輔助工程，2021，30（4）：22-25，31.

[5] 徐芝綸.彈性力學(xué)（上冊(cè)）[M].5版.北京：高等教育出版社，2016.

[6] 李銳.矩形板問(wèn)題的Hamilton求解方法[D].大連：大連理工大學(xué)，2012.

[7] 徐芝綸.彈性力學(xué)（下冊(cè)）[M].5版.北京：高等教育出版社，2016.

[8] 謝中敏，胡超，沈朝萍，等.基于ANSYS軟件均布載荷作用下方板撓度分析[J].數(shù)字技術(shù)與應(yīng)用，2018，36（9）：40-41.