王悅帆,張海峰,王彥宗

(西安建筑科技大學(xué)理學(xué)院,陜西西安710055)

薄板的屈曲問題是工程中比較常見的問題,若利用普通解法求解四邊固定正交各向異性矩形板的屈曲問題比較困難,計算復(fù)雜,不易得到正確結(jié)果。但是如果利用傅立葉級數(shù)求解,求解起來就變的比較容易。

此法首先須將撓度假設(shè)為傅立葉級數(shù)形式,然后求出撓度各階導(dǎo)數(shù)的傅立葉級數(shù)形式,再將其代入邊界條件和屈曲的微分方程進行求解。為簡化計算取1/4板為研究對象。對于四邊固定的邊界條件,可以取兩臨邊簡支,其中一邊受有彎矩來求得一個撓度表達式,再將兩個邊的邊界條件對調(diào),求得另外一個撓度表達式,此時將兩個表達式相加求得疊加解,從而求得臨界力。

1 基本方程

正交各向異性板屈曲的微分方程

當不受縱向剪力,即Nxy=0時

2 疊加解

在邊界條件對稱時(如圖1所示),其解可分為:關(guān)于x=a/2與y=b/2雙軸對稱;關(guān)于x=a/2與y=b/2雙軸反對稱;關(guān)于x=a/2對稱而關(guān)于y=b/2反對稱;關(guān)于x=a/2反對稱而關(guān)于y=b/2對稱等4個問題。矩形板屬于關(guān)于x=a/2與y=b/2雙軸對稱的問題。此問題給出兩個疊加解。所有疊加解的傅立葉級數(shù)都取如下形式:

2.1 疊加解一

圖1 邊界條件對稱板

圖2 疊加解一受力圖

如圖2所示,矩形板兩鄰邊簡支,且在x=0的邊上受有彎矩。其邊界條件為:

故x=0邊上的邊界條件可化為:

且y=0邊上的邊界條件為:

假設(shè)其三階導(dǎo)數(shù)為以下形式:

由傅里葉級數(shù)的相關(guān)性質(zhì)計算求得:

二階導(dǎo)數(shù)可直接微分得到:

將式(5)代入式(4),得:

將式(6)式代入式(3)式得三階導(dǎo)數(shù):

進而求得各階偏導(dǎo)數(shù):

將以上各階偏導(dǎo)數(shù)代入式(1)得:

由于sinαmxsinβny≠0,故

解得:

將其代入式(2),即可求得撓度表達式:

即得疊加解一。

2.2 疊加解二

如圖3所示,矩形板兩鄰邊簡支,且在y=0的邊上受有彎矩。其邊界條件為:

且x=0的邊界條件為:

圖3 疊加解二受力圖

利用與疊加解一同樣的方法可以得到:

求得疊加解二的撓度表達式為:

3 邊界條件的疊加

將疊加解一與疊加解二相加,即式(7)與式(8)相加,求得兩鄰邊簡支且兩邊都受有彎矩時的撓度表達式,即

即可求出滿足全部邊界條件的解。

對式(9)求偏導(dǎo)數(shù),得:

由上兩式聯(lián)立,只需求得方程組中的cn和dm,然后將求解得到的cn和dm回代到式(9)中便可最終求得撓度w表達式。

4 數(shù)值計算

如圖4所示的圖形,四邊固定的正交各向異性矩形薄板,它的兩對邊受有均布壓力,在板邊的每單位長度上為Px=1 kN/m,其中E1/E2=40,G12/E2=0.5,V12=0.25。

用matlab計算結(jié)果如表1

?

[1]徐芝綸.彈性力學(xué)[M].北京:高等教育出版社,1990

[2]王克林,劉俊卿,趙冬.平板的彎曲、振動和屈曲[M].北京:冶金工業(yè)出版社,2006

[3]王克林.帶有自由邊的矩形板的重三角級數(shù)解[J].西安建筑科技大學(xué)學(xué)報,1988,20(3):47-57

[4]黃義,王克林.反對稱角鋪層矩形板振動和屈曲分析[J].振動工程學(xué)報,1992,5(4):403-409.

[5]Kennedy JB,ASCE M,Prabhakaara M K.Combined-load buckling of orthotropic skew plates[J].Journal of engineering mechanics division,1979,105(1):71-79

[6]Wang C M,Liem K M,Xiang Y.Buckling of rectangular Mindlin Plates with internal line supports[J].JSolidsand Stuctures,1993,30(1):1-17