正交各向異性矩形板的自由振動特性分析

曾軍才， 王久法， 姚　望， 于　濤

(1.中國船舶重工集團公司第七一○研究所，湖北 宜昌　443003; 2.上海帝西恩精密工具有限公司，上?！?00137)

?

正交各向異性矩形板的自由振動特性分析

曾軍才1， 王久法1， 姚望2， 于濤2

(1.中國船舶重工集團公司第七一○研究所，湖北 宜昌443003; 2.上海帝西恩精密工具有限公司，上海200137)

正交各向異性矩形板越來越廣泛應(yīng)用于航空、航天、船舶和車輛等各種現(xiàn)代工程中，由于其在各個方向上具有不同力學(xué)性質(zhì)，原有的各向同性板理論已不適用，因此對其振動特性的研究具有重要的工程應(yīng)用價值和理論價值。近年來，許多學(xué)者進行了大量研究，取得了一系列成果。正交各向異性板的求解方法主要為數(shù)值方法，如Reyleigh-Ritz法[1]，微分求積法[2]和有限元法[3]等方法；也有一些學(xué)者使用解析解進行求解，如Gorman[4]用疊加法，Kshirsagar等[5]使用無截斷無窮級數(shù)的疊加原理，Hurlebaus等[6]基于Galerkin原理，張承宗等[7]用復(fù)級數(shù)展開法，Huang等[8]使用Green函數(shù)法，黃炎等[9]和Xing等[10]采用分離變量法，Bercin等[11]用Kantorovich法，Sakata等[12]和Rahbar Ranji等[13]使用延拓Kantorovich法分析了正交各向異性板的自由振動特性。不過這些解析法研究的矩形板的邊界條件都是經(jīng)典的邊界條件，而對于更符合工程實際的彈性邊界條件卻很少涉及。

近年來，Li等[14-15]提出了一種改進的Fourier級數(shù)方法進行了板梁結(jié)構(gòu)的自由振動分析，通過將板梁結(jié)構(gòu)的位移函數(shù)表示為傅里葉余弦級數(shù)和輔助多項式或者輔助級數(shù)的線性組合，使得彈性約束邊界條件能夠得到精確滿足。

本文采用改進Fourier級數(shù)的方法，將正交各向異性矩形薄板的位移函數(shù)表示為標(biāo)準(zhǔn)的二維Fourier余弦級數(shù)和輔助Fourier級數(shù)之和的形式。結(jié)合Kirchhoff理論，建立了矩形板在任意邊界條件下的自由振動模型，推導(dǎo)出與控制方程等價的矩陣表達式，板結(jié)構(gòu)的振動模態(tài)可以通過求解矩陣特征值而得到。最后通過數(shù)值仿真，驗證了本文方法的準(zhǔn)確性和快速收斂性。

1控制微分方程求解

本文所研究的正交各向異性矩形板的模型如圖1所示，板結(jié)構(gòu)的四個邊界處分別均勻地布置橫向位移彈簧和旋轉(zhuǎn)約束彈簧，通過改變剛度值，實現(xiàn)對任意彈性邊界條件的模擬。所有的經(jīng)典邊界條件都能夠通過將彈簧系數(shù)設(shè)置為無窮大或零來簡單地獲得。例如將四邊的橫向位移約束彈簧剛度值和旋轉(zhuǎn)約束彈簧剛度值同時設(shè)置為無窮大，就相當(dāng)于模擬了四邊固支的邊界條件。

圖1　任意邊界條件下板結(jié)構(gòu)示意圖Fig.1 A plate with general elastic boundary support

根據(jù)薄板振動理論，各向異性矩形板自由振動的控制方程為：

(1)

式中：D1=E1h3/[12(1-v1v2)]，D2=E2h3/[12(1-v1v2)]，D3=D12+2D66，D12=v1D2=v2D1，D66=G12h3/12，w為撓度，ρ為密度，h為厚度，E1、E2、v1、v2和G12為板的彈性常數(shù)，ω為角頻率。

彎曲振動的位移場可以通過沿x和y軸方向的兩個分量來描述，本文中，位移采用二維改進傅里葉級數(shù)展開來表示：

(2)

式中：λm=mπ/a，λn=nπ/b，Amn、clm、dln分別為用來描述板結(jié)構(gòu)彎曲振動未知的Fourier系數(shù)和輔助級數(shù)的系數(shù)。與x相關(guān)的輔助函數(shù)分別表示為：

(3)

(4)

(5)

(6)

與y相關(guān)的輔助函數(shù)可以將式(3)-(6)中的a和x分別用b和y進行替換得到。通過輔助級數(shù)的引入，解決了振動位移導(dǎo)數(shù)在邊界不連續(xù)的問題。從而此位移函數(shù)可以同時滿足位移邊界條件和力的邊界條件。這種改進的傅里葉級數(shù)解能適用于任意的彈性邊界條件，同時也能改善級數(shù)的收斂性。

將式(2)代入式(1)中有：

(7)

將所有的輔助級數(shù)及其導(dǎo)數(shù)均展開為Fourier余弦級數(shù)，并利用方程左右兩端余弦項系數(shù)相等有：

(8)

式中：m=0,1,2,…,n=0,1,2,…，輔助級數(shù)及其導(dǎo)數(shù)的Fourier展開為：

(9)

(10)

當(dāng)所有級數(shù)展開在數(shù)值計算過程中均截斷于m=M和n=N，方程(8)可以寫為如下矩陣形式：

(BA+CP)-ρhω2(EA+FP)=0

(11)

式中：A=[A00,A01,…,A0N,A10,A11,…,A1N,…,AM0,AM1,…,AMN]，P=[c10,c11,…,c1M,c20,c21,…,c2M,…c40,c41,…,c4M,d10,d11,…,d1N, …,d40,d41,…,d4N] 。

2彈性約束邊界條件

根據(jù)薄板理論，彈性約束邊界條件寫為：

在y=0上,

(12)

(13)

在y=b上,

(14)

(15)

在x=0上,

(16)

(17)

在x=a上,

(18)

(19)

式中：kx0和Kx0(kxa和Kxa)分別表示x=0(x=a)處橫向位移和旋轉(zhuǎn)約束彈簧剛度，ky0和Ky0(kyb和Kyb)分別表示y=0(y=b)處橫向位移和旋轉(zhuǎn)約束彈簧剛度。

將式(2)代入式(12)~(19)中，并將輔助級數(shù)展開Fourier級數(shù)，由方程兩端余弦項系數(shù)相等有：

(20)

(21)

(22)

(23)

上面四個式子中，m=0,1,2,…。

(24)

(25)

(26)

(27)

式(24)~(27)中，n=0,1,2,…。

當(dāng)級數(shù)的截斷數(shù)取為m=M和n=N時，式(20)-(27)可寫成矩陣表達示為：

ΗP=QA

(28)

將式(28)代入式(11)中，可得到最終的系統(tǒng)方程為

[K-(ρhω2/D1)M]A=0

(29)

式中：K=(B+CH-1Q)/D1，M=(E+FH-1Q)。

通過求解這一標(biāo)準(zhǔn)和矩陣特征值即可得到正交各向異性矩形板的固有頻率和特征向量，每階特征向量為所對應(yīng)結(jié)構(gòu)模態(tài)形狀分布的Fourier系數(shù)。

3數(shù)值結(jié)果與討論

文中用C表示固支邊界條件，F(xiàn)表示自由邊界條件，S表示簡支邊界條件。為了驗證本文方法的準(zhǔn)確性，表1給出了不同長寬比和邊界條件下矩形薄板的前六階無量綱頻率參數(shù)Ω=ωa2(ρh/D1)1/2，矩形板的彈性參數(shù)為D2=0.5D1,D3=0.5D1。表中SSSS表示沿邊界x=0，y=0，x=a和y=b的邊界條件都為簡支，簡支邊界條件可通過設(shè)置橫向位移約束彈簧剛度值為無窮大，旋轉(zhuǎn)約束彈簧剛度值為零，在本文中，無窮大取為107×D1。同時表中也給出了文獻[12]中采用延拓Kantorovich法得到的結(jié)果；表2給出了彈性參數(shù)為D2=D1,D3=0.5D1的矩形薄板在不同長寬比和邊界條件下的無量綱頻率參數(shù)，同時也給出了文獻[10]通過分離變量法得到的計算結(jié)果。通過比較可以發(fā)現(xiàn)本文方法得到的結(jié)果與其它方法得到的結(jié)果吻合良好，在本文計算過程中，兩個方向的位移展開采用相同的截斷數(shù)，取值為M=N=14。同時，圖2給出了采用本方法得到的D2=D1,D3=0.5D1的矩形薄板在a/b=1時的前四階固有振型圖。

表1　 D2=0.5D1, D3=0.5D1時矩形板頻率的頻率參數(shù)

表2　 D2=D1, D3=0.5D1時矩形板頻率的頻率參數(shù)

圖2　SSCC矩形板前四階的振型Fig.2 The first four mode shapes for the SSCC rectangular plate

為了檢驗本文方法的收斂性，表3給出了各向異性矩形薄板在不同截斷數(shù)時的計算結(jié)果，板的長寬比a/b=1，彈性參數(shù)為D2=D1,D3=0.5D1，邊界條件為CCCC。從表中可以看出，M=N=5和M=N=20時得到的前八階無量綱固有頻率的最大偏差為0.22%，即本方法具有較好的收斂性，當(dāng)截斷數(shù)取較小的值時就能得到比較精確的結(jié)果。

本文方法不僅可以計算經(jīng)典邊界條件下矩形薄板的振動頻率，而且適用于彈性支撐的矩形板。表4給出了長寬比a/b=1的正交各向異性矩形薄板在不同彈性支撐剛度下的前8階無量綱振動頻率曲線。x=0和x=a上旋轉(zhuǎn)約束彈簧剛度值都為無窮大，y=0和y=b上旋轉(zhuǎn)約束彈簧剛度值都為零，而四邊上的橫向位移約束彈簧剛度值都為10k×D1。從表中可以看出，隨著彈簧剛度值的增加，矩形板的頻率也隨之增大，當(dāng)k≥5時，數(shù)值幾乎不在變化，頻率值變?yōu)檫吔鐥l件為CSCS下的結(jié)果。從這里也能看出，本文看中剛度值無窮大取為107×D1是合適的。

表3　正交各向異性矩形板頻率的收斂性

表4　 彈性邊界條件下板的振動頻率

4結(jié)論

(1) 本文采用橫向位移約束彈簧和旋轉(zhuǎn)約束彈簧來模擬任意的彈性邊界條件，所有經(jīng)典的邊界條件都可以通過設(shè)置四邊上的彈簧剛度值為無窮大或者零來模擬。

(2) 文中采用改進Fourier級數(shù)法建立了正交各向異性矩形薄板的自由振動模型，得到了板結(jié)構(gòu)固有頻率的解析解。振動位移函數(shù)表示為一個標(biāo)準(zhǔn)的二維Fourier余弦級數(shù)和輔助級數(shù)的線性疊加。通過輔助級數(shù)的引入，解決了振動位移的導(dǎo)數(shù)在邊界潛在的不連續(xù)性問題，從而使改進Fourier級數(shù)法適用于任意彈性邊界條件下的矩形板。

(3) 本文方法中，所有的頻率參數(shù)可以通過求解一個標(biāo)準(zhǔn)的特征值問題而得到。對于不同邊界條件下的模態(tài)參數(shù)，只需設(shè)置相應(yīng)的邊界彈簧剛度值，然后就可以通過統(tǒng)一的程序進行求解。

(4) 最后進行了數(shù)值分析，并與其它文獻的結(jié)果進行了比較，驗證了本方法的快速收斂性和準(zhǔn)確性。

參 考 文 獻

[1] Bhat R B.Natural frequencies of rectangular plates using characteristic orthogonal polynomials in the Rayleigh-Ritz method [J]. Journal of Sound and Vibration, 1985, 102(4):493-499.

[2] Jafari A A,Eftekhari S A. An efficient mixed methodology for free vibration and buckling analysis of orthotropic rectangular plates [J]. Applied Mathematics and Computation, 2011, 218:2670-2692.

[3] Bardell N S, Dunsdon J M, Langley D S.Free vibration analysis of thin coplanar rectangular plate assemblies-Part I: theory, and initial results for specially orthotropic plates[J]. Composite Structures, 1996, 34:129-143.

[4] Gorman D J. Accurate free vibration analysis of clamped orthotropic plate by the method of superposition[J]. Journal of Sound and Vibration, 1990, 140:391-411.

[5] Kshirsagar S, Bhaskar K.Accurate and elegant free vibration and buckling studies of orthotropic rectangular plates using untruncated infinite series [J]. Journal of Sound and Vibration, 2008, 314:837-850.

[6] Hurlebaus S, Gaul L, Wang J T S.An exact series solution for calculating the eigenfrequencies of orthotropic plates with completely free boundary [J]. Journal of Sound and Vibration, 2001, 244(5):747-859.

[7] 張承宗，楊光松. 各向異性板結(jié)構(gòu)橫向彎曲一般解析解[J]. 力學(xué)學(xué)報, 1996, 28(4):429-440.

ZHANG Cheng-zong, YANG Guang-song.General analytic solutions for transverse bending problem of anisotropic plate structures [J]. Acta Mechanica Sinica, 1996, 28(4):429-440.

[8] Huang M, Ma X Q, Sakiyama T, el al.Free vibration analysis of orthotropic rectangular plates with variable thickness and general boundary conditions [J], Journal of Sound and Vibration, 2005, 288:931-955.

[9] 黃炎，廖瑛，謝燕. 雙參數(shù)彈性地基上受壓的正交異性板的自由振動[J]. 工程力學(xué), 2006, 23(3):46-49.

HUANG Yan, LIAO Ying, XIE Yan. Free vibration of compressed orthotropic plate on two parameter elastic foundation [J]. Engineering Mechanics, 2006, 23(3):46-49.

[10] Xing Y F, Liu B. New exact solutions for free vibrations of thin orthotropic rectangular plates [J]. Composite Structures, 2009, 89:567-574.

[11] Bercin A N.Free vibration solution for clamped orthotropic plates using the Kantorovich method[J]. Journal of Sound and Vibration, 1996, 196(2):243-247.

[12] Sakata T, Takahashi K, Bhat R B.Natural frequencies of orthotropic rectangular plates obtained by iterative reduction of the partial differential equation [J]. Journal of Sound and Vibration, 1996, 189(1):89-101.

[13] Rahbar Ranji A, Rostami Hoseynabadi H.A semi-analytical solution for forced vibrations response of rectangular orthotropic plates with various boundary conditions[J]. Journal of Mechanical Science and Technology, 2010, 24:357-364.

[14] Li W L. Free vibrations of beams with general boundary conditions [J]. Journal of Sound and Vibration, 2000, 237(4):709-725.

[15] Li W L, Zhang X F, Du J T. An exact series solution for the transverse vibration of rectangular plates with general elastic boundary supports [J]. Journal of Sound and Vibration, 2009,321(1/2):254-269.

[16] 史冬巖，石先杰，王青山,等. T型耦合板結(jié)構(gòu)振動特性研究[J]. 振動與沖擊, 2014, 33(4): 185-189.

SHI Dong-yan,SHI Xian-jie,WANG Qing-shan,et al. Vibration analysis of a T-coupled plate structure [J].Journal of Vibration and Shock,2014,33(4):185-189.

第一作者 曾軍才 男，高級工程師，1980年10月生

摘要：采用改進Fourier級數(shù)方法，建立了正交各向異性矩形薄板的彎曲振動模型，推導(dǎo)出與振動控制方程等價的矩陣方程，得到控制方程在任意邊界條件下的解析解。彎曲振動的位移函數(shù)表示為標(biāo)準(zhǔn)的二維Fourier余弦級數(shù)和輔助Fourier級數(shù)之和，通過輔助級數(shù)的引入，解決了振動位移函數(shù)的偏導(dǎo)數(shù)在各邊界處潛在不連續(xù)的問題。矩形板的振動模態(tài)信息能夠通過求解一個標(biāo)準(zhǔn)的矩陣特征值而得到。最后進行數(shù)值計算并與現(xiàn)有的文獻結(jié)果進行比較，驗證了該方法的快速收斂性和計算精確性。

關(guān)鍵詞：正交各向異性板；改進Fourier級數(shù)；任意彈性邊界條件；解析解

Free vibration characteristics of orthotropic rectangular plates

ZENGJun-cai1,WANGJiu-fa1,YAOWang2,YUTao2(1. No. 710 R&D Institute, CSIC, Yichang 433003, China;2. TCM China Tool Consulting & Management Ltd. Shanghai 200137, China)

Abstract:An improved Fourier series method was proposed to develop the transverse vibration model of orthotropic rectangular plates and derive the matrix equation which is equivalent to governing differential equations. An analytical solution for vibration of plates with general elastic boundary conditions was provided. The vibration displacement was solved as the linear combination of a double Fourier cosine series and an auxiliary series. The use of these supplementary series is to solve the discontinuity problem encountered in the partial differentials of displacement function along the edges. The vibration mode characteristics were obtained by solving the eigen values of the matrix. Several numerical examples were given and the comparison of the results with those of the available literature validates the convergence and correctness of the method.

Key words:orthotropic plates; improved Fourier series; general elastic boundary support; analytic solution

基金項目：國家自然科學(xué)基金(51075215,51475246);江蘇省自然科學(xué)基金(BK20131402)

中圖分類號:TP533

文獻標(biāo)志碼：A DOI:10.13465/j.cnki.jvs.2015.24.021

收稿日期：2014-10-13修改稿收到日期：2014-12-12