Timoshenko梁的第二頻譜分析

摘? ?要：推導(dǎo)Timoshenko梁振動(dòng)微分方程的初參數(shù)解，結(jié)合邊界條件，建立簡支梁的頻率方程. 當(dāng)固有頻率小于臨界頻率時(shí)，頻率方程有雙曲正弦函數(shù)與三角正弦函數(shù)之積的因式，當(dāng)固有頻率大于臨界頻率時(shí)，此因式變成為雙三角正弦函數(shù)之積，此即Timoshenko梁產(chǎn)生第二頻譜的理論原因. 推導(dǎo)出等截面等跨徑的2～3跨連續(xù)Timoshenko梁的頻率方程，并從理論上預(yù)測存在第二頻譜現(xiàn)象的其他結(jié)構(gòu). 建立了簡支Timoshenko梁第一、二頻譜的頻率計(jì)算公式. 通過實(shí)例驗(yàn)證第二頻譜的存在. 通過微分方程求解，論證了臨界頻率是結(jié)構(gòu)固有頻率的有效組成部分，其對應(yīng)的豎向位移模態(tài)無振幅、轉(zhuǎn)角位移模態(tài)的振幅為常數(shù);指出數(shù)值分析時(shí)，由于計(jì)算機(jī)截?cái)嗾`差的影響，所預(yù)測的臨界頻率有誤差、所對應(yīng)的豎向位移模態(tài)為不規(guī)則模態(tài)等特點(diǎn).

關(guān)鍵詞：頻率分析;臨界頻率;第二頻譜;Timoshenko梁

中圖分類號：TU375? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Analysis on the Second Frequency Spectrum of Timoshenko Beam

XIA Guiyun

（School of Civil Engineering，Changsha University of Science and Technology，Changsha 410114，China）

Abstract：Incorporating the boundary conditions， initial parameter solutions of vibration differential equations for Timoshenko beam are used to derive the frequency equation of a simply-supported beam. When the natural frequency is less than the critical frequency， the frequency equation can be factorized into the hyperbolic sine function and the trigonometric sine function， while， when the natural frequency is greater than the critical frequency， the frequency equation can be factorized into double trigonometric sine functions， which is the crucial reason for the existence of the second frequency spectrum. Frequency equations for two-span and three-span continuous Timoshenko beams with uniform cross sections and equal spans are derived. Other structures with the second frequency spectrum are forecasted theoretically. The formulas for the first and second frequencies are deduced for simply-supported Timoshenko beam. The existence of the second frequency spectrum is confirmed through the examples. Through solving the differential equation of motion， the critical frequency is proven to be an efficient part of the natural frequencies for the framed structures. The corresponding mode shape of the critical frequency contributes to the displacement mode shape with zero amplitude and rotation mode shape with constant amplitude. Due to the truncation error of the computer， the critical frequency predicted by the finite element method shows error， and the mode shape of the displacement is very irregular.

Key words：frequency analysis;cutoff frequency;second frequency spectrum;Timoshenko beam

梁的彎曲振動(dòng)是土木[1]、機(jī)械、石油、化工、航空航天等領(lǐng)域的重要問題，已有多種理論. 最早、最經(jīng)典的理論是Euler梁模型，該模型考慮了梁的彎曲和截面慣性力，適應(yīng)于細(xì)長桿系結(jié)構(gòu)的分析計(jì)算，但對高跨比較大的深梁結(jié)構(gòu)，存在靜力問題計(jì)算撓度偏小[2]、動(dòng)力問題高估振動(dòng)頻率和有無限階次頻率[3]等不足. Elishakoff等[4]評述Rayleigh梁考慮了截面轉(zhuǎn)動(dòng)慣量的影響，對Euler梁進(jìn)行改進(jìn). Shear梁模型考慮結(jié)構(gòu)彎曲、剪切變形的影響和截面慣性力，但沒有考慮截面轉(zhuǎn)動(dòng)慣量影響（此模型不同于一般的不考慮結(jié)構(gòu)彎曲的簡單剪切梁模型和純剪切梁模型[5]）. 1921年Timoshenko[6-7]綜合考慮截面彎曲和剪切變形、慣性力和轉(zhuǎn)動(dòng)慣量的影響，提出經(jīng)典的Timoshenko梁模型，該模型保留梁的平截面假定，放棄直法線假定，通過引入截面剪切修正系數(shù)來彌補(bǔ)剪切本構(gòu)關(guān)系方面的不足（假定截面剪應(yīng)力不均勻、剪應(yīng)變均勻）[8-9]. 相比于Euler梁、Rayleigh梁和Shear梁，Timoshenko梁有顯著進(jìn)步，提高了計(jì)算精度，擴(kuò)大了應(yīng)用范圍，桿系結(jié)構(gòu)的靜力、動(dòng)力和穩(wěn)定問題都可基于Timoshenko梁理論進(jìn)行分析[1]. 但該理論頗具爭議，至今仍眾說紛紜的第二頻譜問題[10-12]和截面剪切修正系數(shù)定義問題[8，13-15]，還有振動(dòng)微分方程解耦后存在撓度關(guān)于時(shí)間四階導(dǎo)數(shù)項(xiàng)，其物理意義不明確的問題，因此其后出現(xiàn)眾多的修正理論. 1927年Love[16]根據(jù)梁段微元體平衡，提出Timoshenko梁的修正模型，即忽略撓度關(guān)于時(shí)間四階導(dǎo)數(shù)項(xiàng)，可稱

為Love梁. 陳镕等[17]采用雙撓度理論也推導(dǎo)出了與Love梁相同的微分方程，并認(rèn)為導(dǎo)致Timoshenko梁模型中出現(xiàn)撓度關(guān)于時(shí)間四階導(dǎo)數(shù)項(xiàng)的原因是，沒有考慮剪切轉(zhuǎn)角的轉(zhuǎn)動(dòng)慣量，如果舍棄撓度關(guān)于時(shí)間四階導(dǎo)數(shù)項(xiàng)，則可考慮截面剪切轉(zhuǎn)角的轉(zhuǎn)動(dòng)慣量影響. Elishakoff等[18-19]同樣導(dǎo)出了與Love梁相同的微分方程，并認(rèn)為此理論比Timoshenko梁理論更一致、更簡單，其命名為截?cái)郥imoshenko梁. Love梁雖然形式比Timoshenko梁簡化，微分方程求解方便，但是其沒有對應(yīng)的能量泛函，不能通過變分原理導(dǎo)出，也沒有對應(yīng)的有限元列式. Xia等[20]研究了考慮截面剪切變形和全轉(zhuǎn)動(dòng)慣量影響的Timoshenko梁振動(dòng)特性，證明此種Timoshenko梁修正理論無第二頻譜問題和結(jié)構(gòu)固有頻率有界特性.

Timoshenko梁的第二頻譜現(xiàn)象是指一種振型對應(yīng)兩個(gè)固有頻率. 兩端簡支、兩端導(dǎo)向和簡支-導(dǎo)向的單跨Timoshenko梁，多跨連續(xù)的等截面等跨徑Timoshenko梁都存在第二頻譜現(xiàn)象. Traill-Nash[21]于1953年最先發(fā)現(xiàn)和報(bào)道簡支Timoshenko梁存在第二頻譜現(xiàn)象，相繼得到Anderson[22]、Dolph[23]等學(xué)者的確認(rèn)，但也有學(xué)者[24-27]認(rèn)為第二頻譜沒有物理意義而應(yīng)舍棄. 有些學(xué)者[4，28]則認(rèn)為結(jié)構(gòu)的振型包括了豎向變形和轉(zhuǎn)角，如果將變形和轉(zhuǎn)角同等看待，則振幅不一致的振型不能認(rèn)為是同一振型，因此也就不存在第二頻譜問題. 現(xiàn)在越來越多的實(shí)驗(yàn)測試結(jié)果[29-32]證實(shí)了第二頻譜現(xiàn)象不僅存在，而且實(shí)驗(yàn)測試的結(jié)構(gòu)固有頻率與Timoshenko梁理論預(yù)測結(jié)果符合較好，因此沒有理由草率地去否定、甚至舍棄第二頻譜. 本文試圖對Timoshenko梁第二頻譜產(chǎn)生的原因進(jìn)行理論解釋，以期提高結(jié)構(gòu)模態(tài)識(shí)別精度，促進(jìn)以模態(tài)疊加法等為基礎(chǔ)的動(dòng)力分析方法的發(fā)展，謀求在Timoshenko梁第二頻譜問題上取得共識(shí).

1? ?Timoshenko梁振動(dòng)的理論解析

Timoshenko梁振動(dòng)的微分方程[1，2，20]求解采用分離變量法，設(shè)豎向位移w（x）分解為豎向位移函數(shù)W（x）和時(shí)間函數(shù)T（t），如下：

w（x，t） = W（x）·T（t）? ? ? ? （1）

式中：T（t） = a1 sin（ωt） + a2 cos（ωt）

則解耦后的振動(dòng)微分方程為：

■ + ω2■ + ■■ +

■ - ■W = 0? ? ? （2）

式中：D = EI、C = μGA，A、I和μ分別為截面的面積、抗彎慣性矩和剪切修正系數(shù);E、G和ρ分別為材料的彈性模量、剪切模量和密度;ω為結(jié)構(gòu)的圓頻率.

Timoshenko梁振動(dòng)存在臨界頻譜（或稱為移頻頻率，其定義為ωC = ■），式 （2）可根據(jù)結(jié)構(gòu)固有頻率ω與臨界頻率ωC的大小關(guān)系有3種解.

1）固有頻率小于臨界頻率時(shí)（ω < ωC）

此時(shí)結(jié)構(gòu)的豎向變形、轉(zhuǎn)角、剪力和彎矩用初參數(shù)（即x = 0時(shí)的W0、ψ0、Q0、M0）表示為：

W（x）=a1（x）W0+b1（x）ψ0+c1（x）Q0+d1（x）M0ψ（x）=a2（x）W0+b2（x）ψ0+c2（x）Q0+d2（x）M0Q（x）=a3（x）W0+b3（x）ψ0+c3（x）Q0+d3（x）M0M（x）=a4（x）W0+b4（x）ψ0+c4（x）Q0+d4（x）M0（3）

其中，波數(shù)定義為：

α=■■β=■■

相關(guān)系數(shù)定義為：

γ = α1 + ■，φ = β1 - ■? ? ?（4）

式（3）振型函數(shù)ai（x），bi（x），ci（x），di（x），i=1，2，3，4，詳見文獻(xiàn)[2，20]，為節(jié)約篇幅不再列出.

2）固有頻率大于臨界頻率時(shí)（ω > ωC）

此時(shí)結(jié)構(gòu)的豎向變形、轉(zhuǎn)角、剪力和彎矩仍可用式（3）表示，但波數(shù)、相關(guān)系數(shù)和振型函數(shù)需另行定義. 波數(shù)定義為：

α′=■■β=■■

（5）

相關(guān)系數(shù)定義為：

γ′ = α1 - ■，φ = β1 - ■? ? ?（6）

振型函數(shù)ai（x），bi（x），ci（x），di（x），i=1，2，3，4的具體表達(dá)式與文獻(xiàn)[20]一致.

3）固有頻率等于臨界頻率時(shí)（ω = ωC）

微分方程退化為：

■ + β2■ = 0? ? ? （7）

式中：波數(shù)退化為α=0，β=■.

對于常規(guī)桿系結(jié)構(gòu)，如果采用Timoshenko梁單元進(jìn)行結(jié)構(gòu)的自由振動(dòng)分析，可得到一固有頻率與臨界頻率非常接近，位移振型較特別的模態(tài)，此模態(tài)說明臨界頻率是結(jié)構(gòu)頻譜的有效組成部分，其詳細(xì)討論在本文第5節(jié)進(jìn)行.

2? ?簡支梁的第二頻譜現(xiàn)象

對于兩端簡支的Timoshenko梁，當(dāng)ω < ωC時(shí)，由邊界條件知：當(dāng)x = 0時(shí)，W0 = M0 = 0;當(dāng)x = L時(shí)，WL = ML = 0，得到頻率方程為

sinh（αL）sin（βL） = 0? ? ? ?（8）

結(jié)構(gòu)固有頻率解為sin（βL） = 0，即β = ki π/L，k = 1，2，3，…，kC . 此時(shí)結(jié)構(gòu)僅有一支固有頻率. 此支頻率的最大個(gè)數(shù)為：

kC = int■■? ? ? ?（9）

當(dāng)ω > ωC時(shí)，同理可推導(dǎo)頻率方程為：

sin（α′L）sin（βL） = 0? ? ? ?（10）

結(jié)構(gòu)固有頻率解為sin（α′L） = 0或sin（βL） = 0，對應(yīng)解為α′ = nπ/L（n = 1，2，3，…）或β = kπ/L（k = kC + 1，kC + 2，…） 結(jié)構(gòu)有2支固有頻譜.

當(dāng)k = n時(shí)，結(jié)構(gòu)振型相同（振幅歸一化處理），但頻率不同，即一種振型對應(yīng)兩種固有頻率，出現(xiàn)第二頻譜現(xiàn)象. 頻率方程式（10）出現(xiàn)兩支解是簡支Timoshenko梁存在第二頻譜的理論原因.

由于簡支Timshenko梁的豎彎振型都呈正弦波形式，含有k個(gè)半波正弦的振型所對應(yīng)的第一、第二頻率見式（11）.

利用式（11）即可快速確定振型和第一、二頻譜.

3? ?存在第二頻譜現(xiàn)象的其他結(jié)構(gòu)

當(dāng)固有頻率從小于臨界頻率變化到大于臨界頻率時(shí)，如果頻率方程有形如簡支Timoshenko梁的變化規(guī)律，則在理論上存在第二頻譜現(xiàn)象. 有此特征的結(jié)構(gòu)有兩端導(dǎo)向（豎向位移活動(dòng)、轉(zhuǎn)角固定）單跨Timoshenko梁、簡支-導(dǎo)向單跨Timoshenko梁、等截面等跨徑的多跨連續(xù)Timoshenko梁等.

3.1? ?單跨Timoshenko梁的頻率方程

兩端導(dǎo)向單跨Timoshenko梁、簡支-導(dǎo)向單跨Timoshenko梁的頻率方程如表1所示.

從表1可以看出，兩端導(dǎo)向和簡支-導(dǎo)向的單跨Timoshenko梁，其頻率方程隨固有頻率的變化都有如兩端簡支的單跨Timoshenko梁的特征，因此也存在第二頻譜現(xiàn)象.

3.2? ?等截面等跨徑的多跨連續(xù)Timoshenko梁

利用式（3）建立傳遞矩陣法，可推導(dǎo)等跨徑等截面的多跨連續(xù)Timoshenko梁的頻率方程.

1）等截面等跨徑的二跨連續(xù)Timoshenko梁

當(dāng)ω < ωC時(shí)：

sinh（αL）sin（βL）（φsinh（αL）cos（βL）-

γcosh（αL）sin（βL）） = 0? ? ?（12）

當(dāng)ω > ωC時(shí)：

sin（α′L）sin（βL）（φsin（α′L）cos（βL）-

γ′cos（αL）sin（βL）） = 0? ? ? （13）

2）等截面等跨徑的三跨連續(xù)Timoshenko梁

當(dāng)ω < ωC時(shí)：

sin（αL）sin（βL）{sinh（αL）sin（βL）[2γφ-

2γφcosh（αL）sin（βL）+（γ2-φ2）sinh（αL）sin（βL）]+

3[φsinh（αL）cos（βL）-γcosh（αL）sin（βL）]2}=0

（14）

當(dāng)ω > ωC時(shí)：

sin（α′L）sin（βL）{sin（α′L）sin（βL）[2γ′φ-

2γ′φcos（α′L）sin（βL）-（γ′2-φ2）sin（α′L）sin（βL）]+

3[φsin（α′L）cos（βL）-γ′cos（α′L）sin（βL）]2}=0

（15）

由式（12）～式（15）可知，等截面等跨徑的二跨、三跨連續(xù)Timoshenko梁中，其頻率方程中的第一個(gè)因式形如式（8）、式（10），因此理論上也存在第二頻譜現(xiàn)象，其振動(dòng)模態(tài)為單跨簡支Timoshenko梁的振動(dòng)模態(tài)在多跨連續(xù)結(jié)構(gòu)中反對稱擴(kuò)展. 對于等截面等跨徑的多于三跨的連續(xù)Timoshenko梁，頻率方程中同樣存在形如式（8）、式（10）的因式和變化規(guī)律，因此理論上也存在第二頻譜問題，但頻率方程過于復(fù)雜，此處不再列出其具體表達(dá)式. 對于等截面但跨徑不等的多跨連續(xù)Timoshenko梁，當(dāng)跨徑比滿足每跨內(nèi)有整數(shù)半波的振動(dòng)條件時(shí)，也同樣存在第二頻譜現(xiàn)象. 由于跨數(shù)、跨徑比、波長等參數(shù)變化過多，頻率方程推導(dǎo)復(fù)雜，一般只能進(jìn)行數(shù)值分析和驗(yàn)證.

4? ?算例驗(yàn)證

兩端簡支單跨Timoshenko梁計(jì)算跨徑10 m，橫截面為1.0 m（寬）×1.8 m（高）的矩形截面，其剪切修正系數(shù)為5/6，材料彈性模量為200 GPa、剪切模量為80 GPa，材料密度為7 850 kg/m3. 結(jié)構(gòu)前50階頻率的理論值（根據(jù)式（8）、式（10）求解）和Ansys數(shù)值結(jié)果如圖1所示，無量綱波數(shù)α（即αL/π）、無量綱波數(shù)β（即βL/π）值如圖2所示.

從圖1可以看出，結(jié)構(gòu)前50階頻率的理論結(jié)果與利用Ansys軟件計(jì)算[33]的數(shù)值結(jié)果（結(jié)構(gòu)劃分為200個(gè)單元、Beam3單元）符合較好，最大誤差不超過0.97%. 圖2中，無量綱波數(shù)αL/π、βL/π如果取整數(shù)k（圖2圖標(biāo)有填充時(shí)為整數(shù)，空心時(shí)為小數(shù)），其對應(yīng)于含有k個(gè)正弦半波的振型，此振型有2個(gè)固有頻率. 根據(jù)式（11），結(jié)構(gòu)第一頻譜的第1、2階振動(dòng)頻率、無量綱波數(shù)和所對應(yīng)的第二頻譜的第1、2階（按固有頻率排序，其對應(yīng)振型為第8、9階振型）振動(dòng)頻率、無量綱波數(shù)如表2所示，即結(jié)構(gòu)的第1階和第8階、第2階和第9階振動(dòng)模態(tài)為頻譜對.

將振型歸一化后，第一頻譜的第1、2階振型和所對應(yīng)的第二頻譜振型第1、2階振型（對應(yīng)振型為第8、9階振型）的位移振型、轉(zhuǎn)角振型如圖3、圖4所示.

從圖3、圖4可以看出，將振型歸一化后，第一頻譜的第1、2模態(tài)與第二頻譜的第1、2階模態(tài)完全相同，驗(yàn)證了兩端簡支的單跨Timoshenko梁存在第二頻譜現(xiàn)象.

簡支-導(dǎo)向單跨Timoshenko梁第一頻譜的第1、2階模態(tài)與第二頻譜的第1、2階模態(tài)（對應(yīng)于第7、8階模態(tài)）的位移、轉(zhuǎn)角模態(tài)如圖5、圖6所示.

從圖5、圖6可以看出，將模態(tài)歸一化后，簡支-導(dǎo)向的單跨Timoshenko梁第一頻譜的第1、2模態(tài)與第二頻譜的第1、2階模態(tài)完全相同，同樣存在第二頻譜現(xiàn)象.

5? ?臨界頻率的振型分析

臨界頻率將Timoshenko梁的振動(dòng)分析為兩區(qū)段的3種特例，其分界點(diǎn)即為臨界頻率. 臨界頻率是結(jié)構(gòu)頻譜的有效組成部分，但其對應(yīng)的模態(tài)非常特別，是Euler梁、Love梁、Shear梁所沒有的，需要特別分析. 本文以簡支Timoshenko梁為例進(jìn)行模態(tài)的理論分析.

5.1? ?臨界頻率所對應(yīng)的模態(tài)特點(diǎn)

當(dāng)結(jié)構(gòu)固有頻率等于臨界頻率時(shí)，微分方程式（7）的解[24-25]為：

W（x）=A1cos（βx）+B1 sin（βx）+C1+D1 xψ（x）=a1cos（βx）+b1 sin（βx）+c1+d1 x? （16）

根據(jù)Timoshenko梁的平衡方程要求，有：

■ = ■+ ■W? ? ? ? （17）

式（16）解中，待定系數(shù)間存在如下關(guān)系：

a1=■·B1，b1=-■·A1，D1=0，d1=■·C1? ?（18）

將待定參數(shù)代入微分方程（7）的解中，有

W（x）=A1 cos（βx）+B1 sin（βx） + C1θ（x）=A1■cos（βx）-A1■sin（βx）+c1+C1■xQ（x）=-A1■sin（βx）+B1■cos（βx）-C1■x-c1·CM（x）=B1Csin（βx）+A1Ccos（βx）-C1■

（19）

引入簡支邊界條件，即W（0）=M（0）=0，得A1=C1=0;W（L）=M（L）=0，得：

B1 sin（βL）=0，B1Csin（βL）=0? ? ? ?（20）

一般條件下sin（βL）≠0，故B1=0，有

W（x）=0，θ（x）=c1，Q（x）=c1·C，M（x）=0? ? （21）

此時(shí)結(jié)構(gòu)振動(dòng)為一種特殊模態(tài)，只有截面轉(zhuǎn)動(dòng)振動(dòng)且振幅恒定，而豎向位移振動(dòng)無振幅.

如果結(jié)構(gòu)的跨徑滿足sin（βL） = 0，即L = kπ/β（k=1，2，3，…，∞），則B1≠0，結(jié)構(gòu)則為有幅振動(dòng)，此時(shí)結(jié)構(gòu)的模態(tài)為：

W（x）=B1sin（βx），θ（x）=B1Ccos（βx）/βD+c1 （22）

5.2? ?臨界模態(tài)特征的有限元預(yù)測

利用有限元軟件進(jìn)行桿系結(jié)構(gòu)的固有頻率、振型分析時(shí)，如果采用Timoshenko梁單元，則可捕捉到臨界頻率和臨界模態(tài). 但是有限元將自由振動(dòng)分析轉(zhuǎn)為特征值問題，由于計(jì)算機(jī)存在截?cái)嗾`差，會(huì)預(yù)報(bào)與臨界頻率理論值極為接近的頻率值，所對應(yīng)的豎向位移模態(tài)未能如理論預(yù)測那樣為無振幅振動(dòng)，而是有幅振動(dòng)且不規(guī)則，但轉(zhuǎn)角位移模態(tài)振幅恒定. 如本文上述算例，采用Ansys軟件計(jì)算的臨界頻率（結(jié)構(gòu)劃分為200個(gè)單元、Beam3單元）為892.675 Hz、理論計(jì)算值為892.602 Hz. 豎向位移和轉(zhuǎn)角位移模態(tài)如圖7、圖8所示.

從圖7可以看出，豎向位移模態(tài)（自由振動(dòng)分析所對應(yīng)的特征向量）最大振幅為1.467 4×10-9，與其他模態(tài)振幅相比小6個(gè)數(shù)量級，因此可以認(rèn)為此模態(tài)的振幅理論上應(yīng)為0，但由于計(jì)算機(jī)截?cái)嗾`差的原因，出現(xiàn)一些不為0的偽數(shù)據(jù)，構(gòu)成振幅無規(guī)律的模態(tài)，理應(yīng)舍棄. 從圖8可以看出，轉(zhuǎn)角位移模態(tài)則與理論預(yù)測一致，為振幅恒定的模態(tài).

6? ?結(jié)? ?論

通過建立Timoshenko梁振動(dòng)的初參數(shù)解，利用此解對Timoshenko梁第二頻譜現(xiàn)象進(jìn)行了理論研究. 由于臨界頻率將結(jié)構(gòu)的固有頻率分為三部分，使得微分方程所對應(yīng)的特征方程根與臨界頻率有關(guān)，其性質(zhì)隨固有頻率發(fā)生變化，從而使得其頻率方程有形如式（8）、式（10）的變化規(guī)律，當(dāng)固有頻率大于臨界頻率時(shí)，頻率方程式（10）有兩支解，即第二頻譜產(chǎn)生的根本原因.

1）所有頻率方程（或者頻率方程中的因式）有形如式（8）、式（10）變化特征的結(jié)構(gòu)，都存在第二頻譜現(xiàn)象，如兩端簡支、兩端導(dǎo)向、簡支-導(dǎo)向的單跨Timoshenko梁，等截面等跨徑多跨連續(xù)Timoshenko梁、滿足每跨內(nèi)有整數(shù)半波振型的等截面不等跨徑多跨連續(xù)Timoshenko梁等.

2）理論和數(shù)值分析表明，臨界頻率是Timoshenko梁結(jié)構(gòu)固有頻率的有效組成部分. 臨界頻率所對應(yīng)的模態(tài)有無振幅豎向位移振型、有恒定振幅的轉(zhuǎn)角位移振型;在特殊條件下，如跨徑滿足L = kπ/β（k為整數(shù)）條件，臨界模態(tài)也可轉(zhuǎn)化為豎向位移有振幅的模態(tài).

3）利用Timoshenko梁單元進(jìn)行桿系結(jié)構(gòu)振動(dòng)的有限元分析時(shí)，能預(yù)測與臨界頻率極為接近的固有頻率、振幅非常小且無規(guī)則的豎向位移模態(tài). 此時(shí)應(yīng)將固有頻率視為臨界頻率、豎向位移模態(tài)視為無振幅模態(tài)，出現(xiàn)誤差的原因是計(jì)算機(jī)的截?cái)嗾`差所引起.

4）Timoshenko梁的第二頻譜現(xiàn)象目前有不同的觀點(diǎn)，但是此現(xiàn)象已被眾多實(shí)驗(yàn)所證實(shí)，也能從理論上進(jìn)行解釋和分析，不應(yīng)輕言其不合理而舍棄，而應(yīng)通過更多研究或者提出更合理的深梁結(jié)構(gòu)理論來驗(yàn)證.

參考文獻(xiàn)

[1]? ? 張玲，歐強(qiáng)，朱幸仁. 成層地基中考慮樁樁相互作用的雙排樁受力變形分析[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，47（11）：120—126.

ZHANG L，OU Q，ZHU X R. Analysis on forced deformation of double row piles considering pile-pile interaction in layered foundation[J]. Journal of Hunan University （Natural Sciences），2020，47（11）：120—126. （In Chinese）

[2]? ? 夏桂云，李傳習(xí). 考慮剪切變形影響的桿系結(jié)構(gòu)理論與應(yīng)用[M]. 北京：人民交通出版社，2008：1—5.

XIA G Y，LI C X. Frame structure theory and its applications including the shear deformation effect[M]. Beijing：China Communications Press，2008：1—5. （In Chinese）

[3]? ? BOTTEGA W J. Engineering vibrations[M]. Boca Raton：CRC Press，2006：511—675.

[4]? ? ELISHAKOFF I，KAPLUNOV J，NOLDE E. Celebrating the centenary of Timoshenkos study of effects of shear deformation and rotary inertia[J]. Applied Mechanics Reviews，2015，67（6）：060802.

[5]? ? HAN S M，BENAROYA H，WEI T. Dynamics of transversely vibrating beams using four engineering theories[J]. Journal of Sound and Vibration，1999，225（5）：935—988.

[6]? ? TIMOSHENKO S P. On the correction for shear of the differential equation for transverse vibration of prismatic bars[J]. Philosophical Magazine，1921，41（6）：744—746.

[7]? ? TIMOSHENKO S P. On the transverse vibrations of bars of uniform cross-section[J]. Philosophical Magazine，1922，43（253）：125—131.

[8]? ? COWPER G R. The shear coefficients in Timoshenko's beam theory[J]. Journal of Applied Mechanics，1966，33：335—340.

[9]? ? 胡海昌. 彈性力學(xué)的變分原理與應(yīng)用[M]. 北京：科學(xué)出版社，1983：139—147.

HU H C. Variational principles in elasticity and applications[M]. Beijing：Science Press，1981：139—147. （In Chinese）

[10]? STEPHEN N G. The second spectrum of Timoshenko beam theory-Further assessment[J]. Journal of Sound and Vibration，2006，292：372—389.

[11]? CAZZANI A，STOCHINO F，TURCO E. On the whole spectrum of Timoshenko beams. Part I：a theoretical revisitation[J]. Zeitschrift Für Angewandte Mathematik Und Physik，2016，67（2）：article 24（1—30）.

[12]? CAZZANI A，STOCHINO F，TURCO E. On the whole spectrum of Timoshenko beams. Part Ⅱ：further applications[J]. Zeitschrift Für Angewandte Mathematik Und Physik（ZAMP），2016，67（2）：article 25（1—21）.

[13]? STEPHEN N G. Timoshenkos shear coefficient from a beam subjected to gravity loading[J]. Journal of Applied Mechanics，1980，47（1）：121—127.

[14]? HUTCHINSON J R. Shear coefficients for Timoshenko beam theory[J]. Journal of Applied Mechanics，2001，68（1）：87—92.

[15]? 王樂，王亮. 一種新的計(jì)算Timoshenko梁截面剪切系數(shù)的方法[J]. 應(yīng)用數(shù)學(xué)和力學(xué)，2013，34（7）：756—763.

WANG L，WANG L. A new method of obtaining Timoshenkos shear coefficients[J]. Applied Mathematics and Mechanics，2013，34（7）：756—763.? （In Chinese）

[16]? LOVE A E. A treatise on the mathematical theory of elasticity[M]. New York：Dover，1927：314—331.

[17]? 陳镕，萬春風(fēng)，薛松濤，等. Timoshenko梁運(yùn)動(dòng)方程的修正及其影響[J]. 同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，33（6）：711—715.

CHEN R，WAN C F，XUE S T，et al. Modification of motion equation of Timoshenko beam and its effect[J]. Journal of Tongji University，2005，33（6）：711—715. （In Chinese）

[18]? ELISHAKOFF I. An equation both more consist and simpler than the Bresse-Timoshenko equation. In：Gilat R;Sills-Banks L（editors），Advances in mathematical modeling and experimental methods for materials and structures[M] Berlin：Springer，2010：249—254.

[19]? ELISHAKOFF I，HACHE F，CHALLAMEL N. Critical contrasting of three versions of vibrating Bresse-Timoshenko beam with a crack[J]. International Journal of Solids and Structures，2017，109：143—151.

[20]? XIA G Y，SHU W Y，STANCIULESCU I. Analytical and numerical studies on the slope inertia-based Timoshenko beam[J]. Journal of Sound and Vibration，2020，473：115227.

[21]? TRAILL-NASH R W，COLLAR A R. The effects of shear flexibility and rotatory inertia on the bending vibrations of beams[J]. The Quarterly Journal of Mechanics and Applied Mathematics，1953，6（2）：186—222.

[22]? ANDERSON R A. Flexural vibrations in uniform beams according to the Timoshenko theory[J]. Journal of Applied Mechanics，1953，20（4）：504—510.

[23]? DOLPH C L. On the Timoshenko theory of transverse beam vibration[J]. Quarterly of Applied Mathematics. 1954，12（2）：175—187.

[24]? BHASHYAM G R，PRATHAP G. The second frequency spectrum of Timoshenko beams[J]. Journal of Sound and Vibration，1981，76（3）：407—420.

[25]? STEPHEN N G，PUCHEGGER S. On the valid frequency range of Timoshenko beam theory[J]. Journal of Sound and Vibration， 2006，297（3-5）：1082—1087.

[26]? CHERVYAKOV A M，NESTERENKO V V. Is it possible to assign physical meaning to field theory with higher derivatives？ [J]. Physical Review D，1993，48（12）：5811.

[27]? NESTERENKO V V. A theory for transverse vibrations of the Timoshenko beam[J]. Journal of Applied Mathematics and Mechanics，1993，57（4）：669—677.

[28]? LEVINSON M，COOKE D W. On the two frequency spectra of Timoshenko beams[J]. Journal of Sound and Vibration，1982，84（3）：319—326.

[29]? D？魱AZ-DE-ANDA A，F(xiàn)LORES J，GUTI？魪RREZ L，et al. Experimental study of the Timoshenko beam theory predictions[J]. Journal of Sound and Vibration，2012，331（26）：5732—5744.

[30]? MONSIVAIS G，D？魱AZ-DE-ANDA A，F(xiàn)LORES J，et al. Experimental study of the Timoshenko beam theory predictions：Further results[J]. Journal of Sound and Vibration，2016，375：187—199.

[31]? TORRES-GUZM？魣N J C，D？魱AZ-DE-ANDA A，OTERO J A，et al. On the warping of the extreme ends of a beam under flexural oscillations[J]. Journal of Sound and Vibration，2018，435：234—245.

[32]? BR？覫NS M，THOMSEN J J. Experimental testing of Timoshenko predictions of supercritical natural frequencies and mode shapes for free-free beams[J]. Journal of Sound and Vibration，2019，459：114856.

[33]? 鄧露，李樹征，淡丹輝，等. 橋梁動(dòng)態(tài)稱重技術(shù)在中小跨徑混凝土梁橋上的適用性研究[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2020，47（3）：89—96.

DENG L，LI S Z，DAN D H，et al. Study on applicability of bridge weigh-in-motion technology in short-to medium-span concrete girder bridges[J]. Journal of Hunan University （Natural Sciences），2020，47（3）：89—96.? （In Chinese）