楊成永 馬文輝 韓薛果 程霖

摘? ?要：以矩形板的Navier解為基礎(chǔ)，采用帶補充項的傅里葉級數(shù)作為撓度函數(shù)，研究了局部均布荷載作用下四邊支承矩形薄板的彎曲問題. 推導了確定待定系數(shù)的線性代數(shù)方程組，給出了簡支邊和固支邊不同組合條件下的統(tǒng)一計算公式. 討論了帶補充項法級數(shù)解的收斂速度，并與疊加法級數(shù)解及有限元數(shù)值解分別進行了精度和計算量的對比. 結(jié)果表明，帶補充項法的級數(shù)解達到收斂的級數(shù)項數(shù)約為40項. 帶補充項法的級數(shù)解與疊加法級數(shù)解具有同樣的求解精度. 有限元解隨網(wǎng)格的細分，計算結(jié)果逐漸接近級數(shù)法解. 級數(shù)解法的計算量與有限元解法相比是微不足道的. 研究成果適于進行構(gòu)筑物頂板受局部均布荷載作用的結(jié)構(gòu)計算.

關(guān)鍵詞：矩形板;四邊支承;局部均布荷載;級數(shù)解;求解精度

中圖分類號：U411? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼：A

Internal Force Calculation of Four Edges Supported Rectangular

Plates under Local Uniformly Distributed Load

YANG Chengyong，MA Wenhui?，HAN Xueguo，CHENG Lin

（School of Civil Engineering，Beijing Jiaotong University，Beijing 100044，China）

Abstract：On the basis of Naviers solution to rectangular plates， the bending problem was studied for the four edges supported thin plates under local uniformly distributed load， where the double Fourier series with additional terms was adopted as the deflection function of the plates. Linear algebraic equations for solving the undetermined coefficients were derived. A unified solution was obtained to the rectangular plates with clamped and simply supported edges. The rate of convergence was discussed on the solution of the series method with additional terms. The proposed method was compared both with superposition series method on accuracy， and with finite element numerical method on computational cost. The results show that 40 terms should be employed for a convergence of the series. The method with additional terms shows the same accuracy of solution as superposition series method does. The solution by finite element method gradually approaches that by the series method as the mesh gets finer and finer. In comparison with finite element method， the computational time by the series method is negligible. This work is applicable for structural analysis of the top plates of underground buildings under truck wheel pressure.

Key words：rectangular plate;four supported edges;locally uniformly distributed load;series solution;solution accuracy

地鐵、熱力和燃氣等地下工程中，地下構(gòu)筑物的頂板多為四邊支承的薄板，板上常承受局部均布荷載如汽車輪壓作用. 為了確定像汽車輪壓這類荷載在板內(nèi)產(chǎn)生的最大撓度和內(nèi)力，需要進行任意位置局部均布荷載作用下?lián)隙群蛢?nèi)力的計算.

對四邊支承的矩形薄板問題，可以從四邊簡支板的Navier解出發(fā)，采用疊加方法[1-2]或加補充項的方法[3-4]解決. 如：蔡長安等[5-6]以帶附加補充項的Fourier級數(shù)作為撓度函數(shù)，求解了Winkler地基及Pasternak地基上自由邊矩形板的彎曲問題. 許琪樓等[7-8]采用一種能滿足自由角點條件的撓度表達式，解決了二鄰邊支承二鄰邊自由矩形板和二鄰邊及對角點支承矩形板的彎曲問題. 他們還采用疊加方法[9-10]，提出了四邊支承矩形板及一邊固定一角點或二角點支承的矩形板的統(tǒng)一求解方法. 岳建勇等[11-12]采用一種雙三角級數(shù)形式的撓度函數(shù)，得到了三邊固定一邊自由及兩對邊固定兩對邊自由矩形板的精確解. 鐘陽等[13]在辛幾何空間中利用分離變量法推導出了四邊固支彈性矩形薄板的精確解析表達式. 于天崇等[14]假定矩形板的抗彎剛度沿板的寬度方向按照一般冪函數(shù)形式變化，研究了四邊簡支一對邊受彎作用下面內(nèi)變剛度矩形板的彎曲問題. 肖閃閃等[15]采用載荷疊加法研究了集中載荷下四邊固支正交各向異性矩形板的線性彎曲，并討論了經(jīng)典Kirchhoff薄板假設(shè)對于集中載荷的適用性.

4.4.2? ?與文獻[16]對比

文獻[16]列出了四邊簡支板中央受局部均布荷載作用時彎矩的計算系數(shù). 為與其對比并避免查表計算中的插值，取泊松比μ = 0，輪壓x方向分布長度c = 2 m，y方向分布長度d = 1 m，其余參數(shù)采用4.2節(jié)的數(shù)據(jù)，板四邊均簡支.

根據(jù) =? = 1.4， =? = 0.4， =? = 0.2，按文獻[16] 中表4-29查得計算彎矩Mx的系數(shù)為0.148 0，My的系數(shù)為0.130 8. 然后有

計算結(jié)果列于表3.

由表3看出，兩種方法的結(jié)果，前3位有效數(shù)字相同. 由于文獻[16]表格的有效數(shù)字是4位，可以認為表3中兩種結(jié)果是一致的.

4.5? ?計算結(jié)果與有限元對比

采用4.2節(jié)的計算參數(shù)，板左邊（x = 0邊）及前邊（y = 0邊）固支，其余兩邊簡支. 按相同的參數(shù)和邊界條件采用ANSYS軟件SHELL63號單元，劃分3種粗細不同網(wǎng)格進行計算. 計算結(jié)果列于表4.

由表4可看出：

1）隨有限元網(wǎng)格的加密，計算結(jié)果逐漸趨于本文的級數(shù)解. 由此可說明，本文級數(shù)解是四邊支承板變形問題的理論解或精確解. 當有限元網(wǎng)格細到5 mm × 5 mm時，撓度及彎矩有5位有效數(shù)字與級數(shù)解相同. 可以認為這時數(shù)值解與級數(shù)解基本一致.

2）表4中的軟件運行所用時間是從數(shù)據(jù)輸入到輸出全部的計算機運行時間. 要達到較高的精度，有限元需要花費的計算機時間大大高于級數(shù)解. 就本算例來說，相差達10萬倍以上. 需要注意的是，本算例中，有限元在5 mm × 5 mm網(wǎng)格時，需要求解的方程組的階數(shù)，不少于（x方向節(jié)點數(shù)5/0.005） × （y方向節(jié)點數(shù)7/0.005） × 6個自由度 = 840萬;而級數(shù)解取40項時需要求解的方程組的階數(shù)僅為40 × 2個固支邊 = 80.

5? ?結(jié)束語

采用帶補充項的撓度函數(shù)，研究了四邊支承矩形薄板的彎曲問題. 給出了局部均布荷載作用下簡支邊和固支邊不同組合條件下的統(tǒng)一計算公式.

對比計算表明，以Navier解為基礎(chǔ)帶補充項的傅里葉級數(shù)解，達到收斂的級數(shù)項數(shù)約為40項. 該級數(shù)解與其他采用疊加法得到的傅里葉級數(shù)解，具有同樣的求解精度. 與有限元數(shù)值法相比，級數(shù)解的計算量十分微小.

值得一提的是，式（3）若換成滿布荷載、線荷載和集中力相應(yīng)的傅里葉系數(shù)，本文方法也適用.

參考文獻

[1]? ? TIMOSHENKO S. WOINOSKY-KRIEGER S. Theory of plates and shells [M]. 2nd ed. New York：McGraw-Hill Book Company，Inc，1959：180—228.

[2]? ? 張福范. 彈性薄板[M]. 2版. 北京：科學出版社，1984：58—62.

ZHANG F F. Elastic thin plates [M]. 2nd ed. Beijing：Science Press，1984：58—62. （In Chinese）

[3]? ? 嚴宗達. 結(jié)構(gòu)力學中的富里葉級數(shù)解法[M]. 天津：天津大學出版社，1989：150—197.

YAN Z D. Fourier series solutions in structural mechanics [M]. Tianjin：Tianjin University Press，1989：150—197. （In Chinese）

[4]? ? 楊端生，黃炎. 矩形板結(jié)構(gòu)的彎曲問題[J]. 湖南大學學報（自然科學版），2004，31（6）：65—69.

YANG D S，HUANG Y. The problem of bending of rectangular plate structure [J]. Journal of Hunan University （Natural Sciences），2004，31（6）：65—69. （In Chinese）

[5]? ? 蔡長安，嚴宗達. Winkler地基上自由邊矩形板橫向彎曲的Fourier級數(shù)解[J]. 貴州工學院學報，1996，25（2）：53—61.

CAI C A，YAN Z D. Fourier series solution for bending problem of the rectangular plates with free edges on Winkler foundation [J]. Journal of Guizhou Institute of Technology，1996，25（2）：53—61. （In Chinese）

[6]? ? 蔡長安. Pasternak地基上自由邊矩形板彎曲問題的Fourier級數(shù)解[J]. 貴州工業(yè)大學學報，1998，27（1）：21—31.

CAI C A. Fourier series solution for bending problem of the rectangular plates with free edges on Pasternak foundation [J]. Journal of Guizhou University of Technology，1998，27（1）：21—31. （In Chinese）

[7]? ? 許琪樓，姬同庚. 二鄰邊支承其余邊自由的矩形板在均布荷載作用下的彎曲解[J]. 土木工程學報，1995，28（3）：32—41.

XU Q L，JI T G. Bending solutions of rectangular plate with two adjacent supported edges and two free edges subjected to uniform load [J]. China Civil Engineering Journal， 1995， 28（3）： 32—41. （In Chinese）

[8]? ? 許琪樓，姜銳，唐國明. 二鄰邊支承二鄰邊自由的矩形板彎曲統(tǒng)一求解方法[J]. 東南大學學報（自然科學版），2000，30（2）：138—142.

XU Q L，JIANG R，TANG G M. Unified solution method of rectangular plate bending with two adjacent supported edges and two free edges [J]. Journal of Southeast University （Natural Science Edition），2000，30（2）：138—142. （In Chinese）

[9]? ? 許琪樓，姜銳，唐國明，等. 四邊支承矩形板彎曲統(tǒng)一求解方法——兼論納維葉解與李維解法的統(tǒng)一性[J]. 工程力學，1999，16（3）：90—99.

XU Q L，JIANG R，TANG G M，et al. Unified solution method on rectangular plate bending with four edges supported [J]. Engineering Mechanics，1999，16（3）：90—99. （In Chinese）

[10]? 許琪樓，姜銳，唐國明，等. 一邊固定一角點或二角點支承的矩形板彎曲統(tǒng)一求解方法[J]. 計算力學學報，1999，16（2）：210—215.

XU Q L，JIANG R，TANG G M，et al. Unified solution method on rectangular plate bending with one edge built-in and one or two corner point supported [J]. Chinese Journal of Computational Mechanics，1999，16（2）：210—215. （In Chinese）

[11]? 岳建勇，曲慶璋. 三邊固定一邊自由矩形板的精確解[J]. 青島建筑工程學院學報，1999，20（1）：16—21.

YUE J Y，QU Q Z. The precise solution of rectangular plate with three edges built in and the fourth edge free [J]. Journal of Qingdao Institute of Architecture and Engineering，1999，20（1）：16—21. （In Chinese）

[12]? 岳建勇，曲慶璋. 兩對邊固定兩對邊自由矩形板的精確解[J]. 青島建筑工程學院學報，2000，21（2）：12—17.

YUE J Y，QU Q Z. Exact solution of rectangular thin plate with two opposite edges clamped and the other two edges free [J]. Journal of Qingdao Institute of Architecture and Engineering，2000，21（2）：12—17. （In Chinese）

[13]? 鐘陽，李銳，劉月梅. 四邊固支矩形彈性薄板的精確解析解[J]. 力學季刊，2009，30（2）：297—303.

ZHONG Y，LI R，LIU Y M. Exact analytic solution of rectangular thin plate with four edges clamped [J]. Chinese Quarterly of Mechanics，2009，30（2）：297—303. （In Chinese）

[14]? 于天崇，聶國雋，仲政. 變剛度矩形板彎曲問題的Levy解[J]. 力學季刊，2012，33（1）：53—59.

YU T C，NIE G J，ZHONG Z. Levy-type solution for the bending of rectangular plates with variable stiffness [J]. Chinese Quarterly of Mechanics，2012，33（1）：53—59. （In Chinese）

[15]? 肖閃閃，陳普會. 集中載荷下四邊固支正交各向異性矩形板的線性彎曲問題[J]. 工程力學，2015，32（6）：28—32.

XIAO S S，CHEN P H. Analytical solutions for bending of clamped orthotropic rectangular plates under a concentrated force[J]. Engineering Mechanics，2015，32（6）：28—32. （In Chinese）

[16]? 《建筑結(jié)構(gòu)靜力計算手冊》編寫組. 建筑結(jié)構(gòu)靜力計算手冊[M]. 北京：中國建筑工業(yè)出版社，1998：227—228.

Working Group of Handbook of Building Structural Statics. Handbook of building structural statics[M]. Beijing：China Architecture & Building Press，1998：227—228. （In Chinese）