李鳳臣， 鄭喜亮， 張麗娜

(1.東華理工大學(xué)建筑工程學(xué)院，江西 撫州 344000;2.中冶賽迪工程技術(shù)股份有限公司，重慶 400013)

近年來，大型土工結(jié)構(gòu)的動(dòng)力學(xué)問題研究越來越受到工程界的廣泛關(guān)注，特別是建筑結(jié)構(gòu)物在動(dòng)力荷載(地震、風(fēng)、沖擊荷載等)作用下的響應(yīng)問題更是不容忽視。雖然結(jié)構(gòu)動(dòng)力響應(yīng)的計(jì)算方法很多，在實(shí)際工程中的應(yīng)用也很廣泛，但結(jié)構(gòu)波動(dòng)和振動(dòng)特性的基礎(chǔ)研究進(jìn)展卻很緩慢(錢七虎，2004;孫占學(xué)等，2004)。

求解任意動(dòng)力荷載作用下結(jié)構(gòu)的動(dòng)力響應(yīng)所采用的方法大多是直接積分法，鐘萬勰(1995)、彭建設(shè)等(1998)、于開平(2004)對(duì)此類方法進(jìn)行深入的研究和探討。但這種方法只能求出近似解，不能求出精確解，且它的計(jì)算精度與時(shí)間步長有很大關(guān)系，有比較突出的精確性、收斂性和穩(wěn)定性的問題(喻曉今，2004;楊恒山，2008)。

求解結(jié)構(gòu)動(dòng)力響應(yīng)的精確解法主要有模態(tài)疊加法(MSM)、傳遞矩陣法(MTM)和回轉(zhuǎn)射線矩陣法(MRRM)。模態(tài)疊加法是比較早的計(jì)算結(jié)構(gòu)動(dòng)力響應(yīng)的方法，最初主要是用在計(jì)算多自由度離散體系動(dòng)力學(xué)方面，在克拉夫(2006)的著作中有詳細(xì)的論述，對(duì)于分布體系，用模態(tài)疊加法求解動(dòng)力響應(yīng)較為困難。林繼德(1992)對(duì)結(jié)構(gòu)動(dòng)力計(jì)算中振型疊加法作了進(jìn)一步研究，修正了模態(tài)疊加法中計(jì)算無質(zhì)量處的位移算法。傳遞矩陣法最早由Thomson(1950)提出，用于求解非均勻梁振動(dòng)問題的矩陣方法，此外還應(yīng)用在分析聲波在多層介質(zhì)中的傳播問題(Adler，1990;Lowe，1995)。芮筱亭等(1995)將傳遞矩陣法應(yīng)用到多體系統(tǒng)的振動(dòng)問題中，提出了計(jì)算多體系統(tǒng)振動(dòng)響應(yīng)的新方法;朱勇濤等(2002)提出了變截面梁彎剪振動(dòng)的傳遞矩陣法求解方法。模態(tài)疊加法是將自振頻率的模態(tài)按照傅立葉積分定理組合而成，但由于自振頻率不完整，使得結(jié)果誤差較大，而傳遞矩陣法的矩陣中含有超越函數(shù)，又有累積誤差。

針對(duì)模態(tài)疊加法和傳遞矩陣法求解結(jié)構(gòu)動(dòng)力響應(yīng)精度方面的不足，Howard等(1998)和Pao等(1999，2003)提出了一種全新的計(jì)算方法，即回傳射線矩陣法，應(yīng)用該方法計(jì)算桁架中應(yīng)力波的傳播，并通過試驗(yàn)證明了計(jì)算的準(zhǔn)確性;之后又將該方法應(yīng)用在框架動(dòng)力響應(yīng)中和平面桁架動(dòng)力響應(yīng)計(jì)算中，均得到了很好的結(jié)果。鄭喜亮(2007)研究表明，回傳射線矩陣法在計(jì)算結(jié)構(gòu)自振頻率及動(dòng)力響應(yīng)時(shí)，不僅可以得到比傳遞矩陣法更好的結(jié)果，而且還能計(jì)算結(jié)構(gòu)的瞬態(tài)響應(yīng)，這是傳遞矩陣法很難做到的。

1 回傳射線矩陣法的基本原理

1.1 散射矩陣S

首先引入雙坐標(biāo)系統(tǒng):局部坐標(biāo)系統(tǒng)和整體坐標(biāo)系統(tǒng)，如圖1所示。對(duì)桿IJ，一個(gè)局部坐標(biāo)系統(tǒng)為x－IJ，以I點(diǎn)為原點(diǎn)，x軸的正方向?yàn)檎较?另外一個(gè)局部坐標(biāo)系統(tǒng)為x－JI，以J點(diǎn)為原點(diǎn)，x軸的負(fù)方向?yàn)檎较颉?/p>

圖1 雙坐標(biāo)系統(tǒng)Fig.1 Double coordinates system

能夠得到桿系統(tǒng)的散射關(guān)系方程如下

其中

矢量d是整體離開波的向量，代表從各個(gè)節(jié)點(diǎn)離開的波;矢量a則是整體到達(dá)波的向量，代表到達(dá)各個(gè)節(jié)點(diǎn)的波。2N×2N的矩陣S是整體散射矩陣。

1.2 相位矩陣P和回傳矩陣R

一段桿內(nèi)的同一個(gè)波，既是左端點(diǎn)的離開波，同時(shí)也是右端點(diǎn)的到達(dá)波。兩者的振幅相同，但相位不同，所以有如下相位關(guān)系

在J點(diǎn)處再引入一個(gè)新的局部矢量?dJ和一個(gè)新的整體離開波矢量?d，整體矢量?d和d的關(guān)系為

其中

故可將矢量a和?d關(guān)系寫成

其中叫做整體相位變換矩陣。

聯(lián)立方程(1)、(3)、(4)，得

此處R(ω)=SPU，稱為回傳矩陣。

2 應(yīng)用傳遞矩陣法計(jì)算Timoshenko梁的動(dòng)力響應(yīng)

下面就只以單跨梁為例，介紹傳遞矩陣法求解結(jié)構(gòu)穩(wěn)態(tài)響應(yīng)的過程。首先引入狀態(tài)矢量的概念，對(duì)于Timoshenko梁上任意點(diǎn)i，其狀態(tài)矢量可以表示為

傳遞矩陣法就是建立節(jié)點(diǎn)與節(jié)點(diǎn)之間狀態(tài)矢量的場傳遞矩陣和節(jié)點(diǎn)上左右兩端之間狀態(tài)矢量的點(diǎn)傳遞矩陣，然后將場矩陣和點(diǎn)矩陣組合起來，形成整體傳遞矩陣，最后根據(jù)邊界條件，求解出梁的自振頻率及相應(yīng)的動(dòng)力響應(yīng)。

如圖2所示的兩端鉸支的單跨Timoshenko梁，跨中施加集中簡諧荷載

圖2 簡諧集中荷載作用的單跨梁Fig.2 Single-span beam under harmonic concentrated load

梁的振動(dòng)頻率與外界荷載相同，為Ω。所以在IJ段梁上，左端I節(jié)點(diǎn)的狀態(tài)矢量為zI(0，Ω)，右端J節(jié)點(diǎn)的狀態(tài)矢量為zJ(l/2，Ω)。兩節(jié)點(diǎn)的狀態(tài)矢量可以寫成:zRI=zI(0，Ω)=MIR，zLJ=zJ(l/2，Ω)=MJR。

其中

3 應(yīng)用回傳射線矩陣法計(jì)算Timoshenko梁的穩(wěn)態(tài)響應(yīng)

3.1 單跨Timoshenko梁的穩(wěn)態(tài)響應(yīng)

Timoshenko梁在頻域的位移解為

Timoshenko梁在各點(diǎn)的轉(zhuǎn)角、剪力、彎矩等響應(yīng)的頻域解為

如圖3所示的單跨兩端固支的Timoshenko梁，跨中受集中簡諧荷載P(x，t)作用，其表達(dá)式與式(7)相同。

圖3 兩端固支單跨梁模型Fig.3 Single-span beam model with both ends fixed

整體的散射方程為:

求得參數(shù)a和d的值，在帶入到式(8)中求得梁上各點(diǎn)位移響應(yīng)的頻域解，將參數(shù)a和d的值帶入式(9)中，則可求出梁上各點(diǎn)的轉(zhuǎn)角、剪力、彎矩等響應(yīng)的頻域解，將頻域解進(jìn)行傅立葉變換，則可得到梁響應(yīng)的時(shí)域解。

3.2 兩跨Timoshenko梁的穩(wěn)態(tài)響應(yīng)

如圖4所示的兩跨連續(xù)梁，兩端為固定支撐，跨中為剛性支撐，其中一跨跨中作用簡諧荷載。

圖4 兩端固支雙跨梁模型Fig.4 Two-span beam model with both ends fixed

系統(tǒng)的整體坐標(biāo)為X-Y，將梁分成三個(gè)單元，四個(gè)節(jié)點(diǎn):1-4，局部坐標(biāo)如圖4中所示。同樣先推導(dǎo)系統(tǒng)的散射矩陣，求解響應(yīng)的方法與單跨梁穩(wěn)態(tài)響應(yīng)的求解方法相同。

4 算例分析

4.1 單跨Timoshenko梁的穩(wěn)態(tài)響應(yīng)算例

如圖3所示的兩端固定的單跨低碳鋼梁，其截面尺寸為:梁高 h=0.5 m，梁寬 b=0.04 m，梁長度L=2 m，密度 ρ=7 850 kg/m3，截面積 A=0.002 m2，彈性模量 E=206 GPa，剪切模量 G=79 GPa，剪切因子κ=0.833 3。這是一個(gè)比較典型的Timoshenko梁，表1給出了三種邊界條件下，單跨梁上跨中施加簡諧荷載，分別采用傳遞矩陣法(MTM)和回轉(zhuǎn)射線矩陣法(MRRM)計(jì)算，各點(diǎn)穩(wěn)態(tài)響應(yīng)幅值的對(duì)比。其中，跨中施加荷載的大小為P0=1 000 N，荷載頻率為Ω =100 rad/s，以梁長為長度基本單位進(jìn)行無量剛計(jì)算，得到的各點(diǎn)的響應(yīng)幅值均為無量綱結(jié)果。

表1 MTM與MRRM計(jì)算單跨梁穩(wěn)態(tài)響應(yīng)各點(diǎn)響應(yīng)幅值對(duì)比Tab.1 Comparison of the amplitude of steady-state response in each point of single-span beam computed with MTM and MRRM

4.2 兩跨Timoshenko梁的穩(wěn)態(tài)響應(yīng)算例

表2給出了兩跨連續(xù)梁在其中一跨跨中施加簡諧荷載，分別采用傳遞矩陣法(MTM)和回轉(zhuǎn)射線矩陣法(MRRM)計(jì)算，梁上各點(diǎn)穩(wěn)態(tài)響應(yīng)幅值的結(jié)果對(duì)比。其中，計(jì)算簡圖同圖4，梁的截面尺寸同算。例4.1，兩跨梁長均為1 m，左跨跨中施加荷載的大小為P=1 000 N，荷載頻率為Ω=100 rad/s，以梁長為長度基本單位進(jìn)行無量剛計(jì)算，得到的各點(diǎn)的響應(yīng)幅值均為無量綱結(jié)果。

表2 MTM與MRRM計(jì)算兩跨連續(xù)梁穩(wěn)態(tài)響應(yīng)各點(diǎn)響應(yīng)幅值對(duì)比Tab.2 Comparison of the amplitude of steady-state response in each point of two-span continuous beam computed with MTMand MRRM

5 結(jié)論

傳遞矩陣法適用于計(jì)算單跨梁、多跨連續(xù)梁的自振頻率、穩(wěn)態(tài)響應(yīng)，在計(jì)算自振頻率時(shí)，前幾階的計(jì)算精度很高，但隨后就會(huì)出現(xiàn)嚴(yán)重的數(shù)值穩(wěn)定問題，而無法繼續(xù)計(jì)算更高階的頻率?；貍魃渚€矩陣法適用于計(jì)算單跨梁、多跨連續(xù)梁的自振頻率、穩(wěn)態(tài)響應(yīng)和瞬態(tài)響應(yīng)。通過以上的對(duì)比，可以看出，回傳射線矩陣法，有適用范圍廣、計(jì)算精確、算法簡潔明了以及便于編程等特點(diǎn)。

林繼德.1992.結(jié)構(gòu)動(dòng)力計(jì)算中振型疊加法的進(jìn)一步研究[J].工程力學(xué)，9(3):131-135.

彭建設(shè)，張敬宇.1998.由Gurtin變分原理求解一維動(dòng)力響應(yīng)的半解析法[J].力學(xué)學(xué)報(bào)，26(6):708-711.

錢七虎.2004.深部巖體工程響應(yīng)的特征科學(xué)現(xiàn)象及“深部”的界定[J].東華理工學(xué)院學(xué)報(bào)，27(1):1-5.

芮筱亭，陸毓琪.1995.多體系統(tǒng)振動(dòng)的傳遞矩陣法[J].宇航學(xué)報(bào)，16(3):41-47.

孫占學(xué)，朱永剛，張文.2004.礦物-水反應(yīng)的地球化學(xué)動(dòng)力學(xué)研究進(jìn)展[J].東華理工學(xué)院學(xué)報(bào)，27(1):14-18.

克拉夫R，彭津J，王光遠(yuǎn)譯.2006.結(jié)構(gòu)動(dòng)力學(xué)[M].北京:高等教育出版社.

楊恒山.2008.向量加權(quán)平均值在橋梁控制網(wǎng)平差中的應(yīng)用[J].東華理工大學(xué)學(xué)報(bào):自然科學(xué)版，31(3):256-258.

于開平.2004.結(jié)構(gòu)動(dòng)力響應(yīng)數(shù)值分析的新的廣義-α方法的頻率域分析[J].振動(dòng)工程學(xué)報(bào)，17(4):393-398.

喻曉今.2004.求梁位移的比擬梁法[J].東華理工學(xué)院學(xué)報(bào)，27(4):398-400.

鄭喜亮.2007.基于回傳射線矩陣法的Timoshenko梁動(dòng)力響應(yīng)分析[D].杭州:浙江大學(xué)碩士論文.

鐘萬勰.1995.暫態(tài)歷程的精細(xì)計(jì)算方法[J].計(jì)算結(jié)構(gòu)力學(xué)及其應(yīng)用，12(1):1-6.

朱勇濤，易傳云.2002.變截面梁彎剪振動(dòng)的傳遞矩陣法求解[J].機(jī)械傳動(dòng)，26(2):22-23.

Adler E L.1990.Matrix methods applied to acoustic waves in multilayers[J].IEEE Transactions on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，37(6):485-490.

Howard S M，Pao Y H.1998.Analysis and experiments on stress waves in planar trusses[J].Journal of Engineering Mechanics，124(8):884-891.

Lowe M J S.1995.Matrix techniques for modeling ultrasonic waves in multilayered media[J].IEEE on Ultrasonics，F(xiàn)erroelectrics，and Frequency Control，42(4):525-542.

Pao Y H，Keh D C，Howard S M.1999.Dynamic response and wave propagation in plane trusses and frames[J].AIAA Journal，37(5):594-603.

Pao Y H，Sun G.2003.Dynamic bending strains in planar trusses with pinned of rigid joints[J].Journal of Engineering Mechanics，129(3):324-332.

Thomson W T.1950.Matrix solution for the vibration of uniform beams[J].Journal of Applied Mechanics，17:337-339.