黃 櫻，郭詠梅

基于改進傅里葉級數(shù)法的任意邊界下梁橫振特性分析

黃 櫻1，郭詠梅2

（1.河南交通職業(yè)技術(shù)學(xué)院，河南 鄭州 450005；2.長安大學(xué)，陜西 西安 710000）

【目的】船舶工程、海洋工程、汽車工程等相關(guān)領(lǐng)域涉及大量的梁結(jié)構(gòu)，其振動特性直接關(guān)系到整體結(jié)構(gòu)的振動控制?！痉椒ā繉⒘航Y(jié)構(gòu)橫向振動位移采用一種改進傅里葉級數(shù)展開，運用Rayleigh-Ritz能量法進行求解，建立任意彈性邊界條件下梁橫向振動特性統(tǒng)一求解模型；在此基礎(chǔ)上，分析不同邊界條件對梁結(jié)構(gòu)振動特性的影響規(guī)律。【結(jié)果與結(jié)論】基于改進傅里葉級數(shù)，利用Rayleigh-Ritz能量法建立的梁結(jié)構(gòu)振動的求解模型與解析解基本相同，具有良好的準(zhǔn)確性，且優(yōu)于Ansys仿真解。本研究可為船海工程等梁結(jié)構(gòu)振動機理研究與優(yōu)化提供參考。

歐拉-伯努利梁；橫振；任意邊界；改進傅里葉法

船舶工程、海洋工程、汽車工程等相關(guān)領(lǐng)域涉及大量的梁結(jié)構(gòu)，特別是船舶龍骨、肋骨、縱骨，甲板橫梁、縱桁，以及海洋平臺主體桁架結(jié)構(gòu)，其振動特性直接關(guān)系到整體結(jié)構(gòu)的動力學(xué)性能，影響整體結(jié)構(gòu)的正常使用。因此，梁結(jié)構(gòu)振動特性分析對整體結(jié)構(gòu)振動控制具有非常重要的研究意義。目前，國內(nèi)外學(xué)者對梁結(jié)構(gòu)在經(jīng)典邊界條件下的振動解析解進行大量研究，研究方法包括集中質(zhì)量法、Rayleigh-Ritz能量法、假設(shè)模態(tài)法、加權(quán)參數(shù)法（伽遼金法）、傳遞矩陣法等。鮑四元等[1]基于Rayleigh-Ritz法研究不同截面彈性多跨梁的橫向振動。許得水等[2]和周海軍等[3]采用改進傅里葉級數(shù)法分別研究任意邊界下桿扭轉(zhuǎn)振動和軸系回旋振動。周海軍等[3-4]推導(dǎo)出帶集中質(zhì)量的梁回旋方程和橫向振動方程。杜敬濤等[5]采用改進傅里葉級數(shù)研究任意邊界下桿縱向振動特性。Hozhabrossadati等[6]采用集中質(zhì)量法研究Euler-Bernoulli梁的自由振動。Kim等[7]和Li[8-9]均使用傅里葉級數(shù)法研究任意邊界梁的橫向振動。Shin等[10]應(yīng)用微分變換法分析Euler-Bernoulli梁在變軸力作用下的橫向振動。

由于梁結(jié)構(gòu)橫向振動求解時會遇到超越方程，計算其固有頻率十分困難。因此，有限差分法[11-12]、邊界元法[13]和有限元法[14]等數(shù)值方法也大量應(yīng)用于梁結(jié)構(gòu)分析。其中，有限元法是研究該類問題的主要方法，Erdelyi等[15]和錢管良[16]分別使用有限元法研究多層梁和裂紋梁的振動特性。然而，數(shù)值方法在求解過程中需要將結(jié)構(gòu)離散化，高頻精確度較差，且不適合進行結(jié)構(gòu)參數(shù)優(yōu)化等方面的研究；同時，在實際結(jié)構(gòu)中，振動結(jié)構(gòu)邊界條件并非理想的經(jīng)典邊界條件，對彈性邊界的研究更具工程價值。鑒于此，本研究提出一種基于改進傅里葉級數(shù)法的任意彈性邊界條件下梁結(jié)構(gòu)振動特性分析模型，統(tǒng)一傅里葉級數(shù)求解方法，分析任意彈性邊界下梁的振動特性，以及不同邊界條件對梁結(jié)構(gòu)振動特性的影響規(guī)律，以期為船舶、海洋平臺等骨架的梁結(jié)構(gòu)振動研究提供設(shè)計參考。

1 理論建模

1.1 彈性約束邊界條件

考慮到梁結(jié)構(gòu)厚度尺寸遠小于長度尺寸，本研究使用適用于求解細長結(jié)構(gòu)振動的Euler-Bernoulli梁理論。圖1為簡化梁結(jié)構(gòu)振動模型，兩端引入線性約束彈簧及旋轉(zhuǎn)約束彈簧以模擬彈性邊界條件，對應(yīng)的彈簧剛度分別為0、k、0和K。

當(dāng)剛度系數(shù)取某有限值時，可表示彈性約束邊界條件。將彈簧剛度系數(shù)取零或無窮大可以模擬自由或固支邊界條件（表1）。和分別為梁的長度、楊氏模量、慣性矩和密度，橫向振動位移用(,)表示。

圖1 彈性邊界約束梁橫向振動模型

表1 彈簧剛度設(shè)定值

根據(jù)梁結(jié)構(gòu)兩端的力平衡關(guān)系，邊界條件方程可寫為如下形式：

在= 0處，

在=處，

1.2 位移場改進傅里葉級數(shù)表達式

在經(jīng)典邊界條件下，振動位移函數(shù)可采用傅里葉正弦或余弦級數(shù)展開，但是在梁結(jié)構(gòu)兩端，即=0或=處，會出現(xiàn)位移導(dǎo)數(shù)始終為零的情況，無法研究彈性邊界條件下結(jié)構(gòu)的振動。為克服傳統(tǒng)傅里葉級數(shù)在彈性邊界的求導(dǎo)不連續(xù)問題，梁結(jié)構(gòu)橫向位移采用改進傅里葉級數(shù)展開為

圖2 輔助函數(shù)及其導(dǎo)數(shù)

由圖2，可知

1.3 Rayleigh-Ritz法求解

基于改進傅里葉級數(shù)，可利用Rayleigh-Ritz能量法結(jié)合Hamilton原理來確定振動位移函數(shù)的未知系數(shù)，從而建立統(tǒng)一的梁結(jié)構(gòu)振動方程求解形式。

無外部作用力時，梁結(jié)構(gòu)的拉格朗日方程為

其中，Le是拉格朗日量，是梁結(jié)構(gòu)的總勢能，是梁結(jié)構(gòu)的總動能。梁結(jié)構(gòu)的總勢能包含兩個部分，分別是結(jié)構(gòu)的勢能及邊界彈簧儲存的彈性勢能。

任意邊界條件下梁（圖1）的橫向振動總勢能表達式為

其中，梁結(jié)構(gòu)勢能beam為

模擬彈性邊界的彈簧勢能spring為

梁結(jié)構(gòu)總動能為

其中，是角頻率，是梁的橫截面積。

結(jié)合Rayleigh-Ritz法及Hamilton原理[17]，對其中每個未知系數(shù)求極值，可以得到一系列方程組

將式(5)帶入式(13)，可得矩陣表達形式

其中，

在矩陣、中，11、22和11、22分別為對應(yīng)傅里葉級數(shù)系數(shù)和輔助函數(shù)系數(shù)的對稱項，12、21和12、21分別為對應(yīng)傅里葉級數(shù)系數(shù)和輔助函數(shù)系數(shù)的耦合項。

通過以上分析，橫向振動的固有頻率和特征向量可以通過求解標(biāo)準(zhǔn)特征值問題得到。每一個特征向量包含結(jié)構(gòu)模態(tài)的所有未知系數(shù)，將系數(shù)帶入振動位移表達式中即可得到真實的模態(tài)。

2 數(shù)值算例分析

本研究選取某海洋工程平臺梁結(jié)構(gòu)進行仿真模型分析，該平臺縱梁為QSTE700TM矩形型材60 mm × 40 mm × 2.5 mm，截面的抗彎慣性矩為2.347 × 10-7m4，橫截面積為4.75 × 10-4m2，長度= 3.0 m，材料密度= 7 850 kg/m3，楊氏模量= 2.06 × 1011Pa。為方便后續(xù)推導(dǎo)及分析，引入量綱為1系數(shù)，線性約束彈簧和旋轉(zhuǎn)約束彈簧剛度系數(shù)（量綱為1）分別為

2.1 準(zhǔn)確性和收斂性驗證

本研究對梁結(jié)構(gòu)自由振動特性進行分析，并與方程解析解、有限元分析計算結(jié)果進行對比。本研究使用Ansys軟件建立梁結(jié)構(gòu)有限元模型，采用Beam 188單元，Ansys仿真結(jié)果與網(wǎng)格精細程度有關(guān)。為排除網(wǎng)格尺寸對結(jié)果的影響，比較不同網(wǎng)格尺寸下簡支梁的Ansys仿真結(jié)果（表2），可見當(dāng)單元數(shù)目達到400（單元尺寸0.75 cm）時，Ansys解已經(jīng)基本收斂，可用于本研究的分析對比。

表2 簡支邊界條件Ansys解收斂情況

將兩端旋轉(zhuǎn)約束彈簧設(shè)為0，使線性彈簧剛度逐漸增大，當(dāng)彈簧剛度足夠大時即可模擬簡支邊界條件。表3列出本研究方法模擬兩端線性約束彈簧剛度逐漸增大時前6階頻率的收斂情況（均取相同的傅里葉截斷數(shù)= 50計算）。由表2可見，當(dāng)邊界約束剛度取值達到1010時，Rayleigh-Ritz能量法結(jié)果已經(jīng)收斂，第1-3階的頻率結(jié)果與解析解相同，第4-6階頻率與解析解非常接近，精度高于Ansys解。

表3 模擬簡支邊界條件線性彈簧剛度系數(shù)收斂情況

采用改進傅里葉級數(shù)法在計算時需要對傅里葉項數(shù)進行截斷，為進一步對比不同項數(shù)下的收斂速度，表4列出兩端簡支邊界條件下，不同截斷數(shù)時梁的前6階固有頻率。由表4可見，隨著截斷數(shù)增加，固有頻率計算精度逐漸提高，當(dāng)截斷數(shù)大于20時，計算結(jié)果已基本收斂，可見改進傅里葉級數(shù)法的收斂速度非?？臁?/p>

2.2 不同邊界條件對固有頻率的影響

從表3中可看出，在彈性邊界約束剛度變化時，梁的橫向振動變化存在一個較為敏感的取值區(qū)間。

圖3給出梁從自由-自由邊界到簡支-簡支邊界，梁的前6階固有頻率變化趨勢。由圖3可見，不同階次固有頻率對邊界約束剛度變化的敏感程度不同，第1階固有頻率敏感區(qū)域在線性彈簧剛度系數(shù)（量綱為1）小于102以下，而當(dāng)剛度系數(shù)（量綱為1）取值從103到104時，第4階到第6階固有頻率變化較快。

表4 改進傅里葉級數(shù)截斷數(shù)收斂情況

圖3 固有頻率隨線性約束彈簧剛度變化

圖4給出梁從簡支-簡支邊界到固支-固支邊界，梁的前6階固有頻率變化趨勢。由圖4可見，梁的前6階固有頻率變化敏感區(qū)域均在旋轉(zhuǎn)約束彈簧剛度系數(shù)（量綱為1）取值100到102之間。

圖4 固有頻率隨旋轉(zhuǎn)約束彈簧剛度變化

圖5和圖6給出梁結(jié)構(gòu)一端固定（= 0處），一端從自由邊界到固支邊界（=處），梁結(jié)構(gòu)前2階固有頻率變化趨勢。由圖5和圖6可見，旋轉(zhuǎn)約束彈簧剛度對固有頻率的影響要明顯大于線性約束彈簧，旋轉(zhuǎn)約束彈簧的影響不能忽視。

圖5 第1階固有頻率隨約束彈簧剛度變化

圖6 第2階固有頻率隨約束彈簧剛度變化

2.3 不同邊界條件對模態(tài)振型的影響

圖7給出梁在自由-自由邊界，本研究方法、解析解和Ansys解下的第1階模態(tài)振型。由圖7可見，本研究方法和解析解、Ansys解完全一致。

圖7 自由狀態(tài)下梁第1階模態(tài)振型

同樣，為分析不同邊界約束剛度變化對梁結(jié)構(gòu)模態(tài)振型的影響，圖8和圖9給出梁一端固定，另一端從自由邊界到簡支邊界，梁結(jié)構(gòu)前2階模態(tài)振型的變化趨勢。由圖8和圖9可見，模態(tài)振型在線性彈簧剛度（量綱為1）敏感區(qū)域會出現(xiàn)劇烈變化，特別是在線性彈簧剛度（量綱為1）敏感區(qū)域101到104之間。

圖8 變化邊界條件第1階模態(tài)振型

圖9 變化邊界條件第2階模態(tài)振型

3 討論

將梁結(jié)構(gòu)橫向振動位移采用改進傅里葉級數(shù)進行展開，帶入相應(yīng)的能量方程，結(jié)合Hamilton原理確定位移函數(shù)中的未知系數(shù)，可快速建立統(tǒng)一的梁結(jié)構(gòu)振動方程求解形式，便于研究任意邊界條件下梁結(jié)構(gòu)的振動特性。本方法可以克服傳統(tǒng)傅里葉正弦或余弦級數(shù)在研究經(jīng)典邊界條件下出現(xiàn)的兩端位移導(dǎo)數(shù)不連續(xù)情況[18]，具有良好的準(zhǔn)確性和快速的收斂性。

船舶和海洋平臺結(jié)構(gòu)中各類桿件，如船舶結(jié)構(gòu)中的骨架如橫梁、肋骨、肋板、縱骨、縱桁等，均可使用本方法進行固有振動特性求解。當(dāng)然，整個船體結(jié)構(gòu)也可視為船體梁結(jié)構(gòu)[19]，使用本方法快速求得模態(tài)信息。對于不規(guī)則的梁結(jié)構(gòu)，由于無法給出其具體的能量方程，求解時需要進行一定的簡化，使用本方法時會影響其結(jié)果的精確性。

船舶結(jié)構(gòu)、海工結(jié)構(gòu)在風(fēng)浪、海流、地震等沖擊力作用下的振動特性分析也具有十分重要的價值[20]，可使用本方法進行下一步研究。對于梁結(jié)構(gòu)在外界激勵力作用下的受迫振動求解，只需將外部作用力所做的功引入方程式(8)，進行求解即可。

后續(xù)將利用改進傅里葉級數(shù)和Rayleigh-Ritz能量法建立的任意邊界下梁結(jié)構(gòu)振動求解方法拓展至海洋工程平臺板殼理論研究方面。

4 結(jié)論

本研究基于Euler-Bernoulli梁理論，利用Rayleigh-Ritz能量法建立梁結(jié)構(gòu)的振動統(tǒng)一求解模型，分析任意彈性邊界下梁的振動特性，得到以下結(jié)論：

1）利用標(biāo)準(zhǔn)傅里葉級數(shù)和輔助函數(shù)構(gòu)建梁結(jié)構(gòu)橫向振動位移函數(shù)，可消除級數(shù)求導(dǎo)過程中邊界處的不連續(xù)性，具有求解過程統(tǒng)一、收斂速度快的優(yōu)勢。

2）采用兩組旋轉(zhuǎn)約束彈簧和線性約束彈簧，可很好地模擬梁結(jié)構(gòu)的任意彈性邊界條件，不同階次固有頻率和模態(tài)振型變化敏感區(qū)域與彈簧剛度有關(guān)，前6階固有頻率和振型的剛度（量綱為1）敏感區(qū)域均出現(xiàn)在104以下。

3）基于改進傅里葉級數(shù)，利用Rayleigh-Ritz能量法建立的梁結(jié)構(gòu)振動的求解模型與解析解基本相同，具有良好的準(zhǔn)確性，且優(yōu)于Ansys仿真解。

[1] 鮑四元, 周靜.不同截面形狀下彈性支撐多跨梁振動特性分析[J].中國艦船研究, 2020, 15(1): 162-169.

[2] 許得水, 杜敬濤, 李文達, 等.任意邊界條件彈性桿結(jié)構(gòu)扭轉(zhuǎn)振動特性分析[J].振動與沖擊, 2017, 36(1): 161-166.

[3] 周海軍, 呂秉琳, 王東華, 等.一種改進傅里葉級數(shù)方法的船舶軸系回旋振動動態(tài)特性研究[J].船舶力學(xué), 2012, 16(8): 962-970.

[4] 周海軍, 呂秉琳, 杜敬濤, 等.用改進傅里葉級數(shù)的方法研究軸系橫向振動特性[J].噪聲與振動控制, 2011, 31(4): 68-72.

[5] 杜敬濤, 許得水, 呂朋, 等.任意邊界條件非局部彈性桿縱振特性分析[J].振動工程學(xué)報, 2016, 29(5): 787-794.

[6] HOZHABROSSADATI S M, SANI A A, MOFID M.Vibration of beam with elastically restrained ends and rotational spring-lumped rotary inertia system at mid-span[J].International Journal of Structural Stability and Dynamics, 2015, 15(2): 1450040.

[7] KIM H K, KIM M S.Vibration of beams with generally restrained boundary conditions using Fourier series[J].Journal of Sound and Vibration, 2001, 245(5): 771-784.

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

[9] LI W L.Comparison of Fourier sine and cosine series expansions for beams with arbitrary boundary conditions[J].Journal of Sound and Vibration, 2002, 255(1): 185-194.

[10] SHIN Y J, YUN J H.Transverse vibration of a uniform Euler-Bernoulli beam under varying axial force using differential transformation method[J].Journal of Mechanical Science and Technology, 2006, 20(2): 191-196.

[11] 許士菊, 王長華.梁振動方程的一個穩(wěn)定的有限差分近似[J].吉林化工學(xué)院學(xué)報, 2007, 24(1): 79-81.

[12] 姚振漢, 王海濤.邊界元法[M].北京: 高等教育出版社, 2010: 25-78.

[13] 韓旭炤, 黃玉美, 張永貴, 等.梁結(jié)構(gòu)振動問題的邊界元法解析[J].農(nóng)業(yè)機械學(xué)報, 2008, 39(7): 178-182.

[14] Mackerle J.Finite Element Vibration Analysis of Beams, Plates and Shells[J].Shock and Vibration, 1998, 6(2): 97-109.

[15] ERDELYI N H, HASHEMI S M.On the Finite Element Free Vibration Analysis of Delaminated Layered Beams: A New Assembly Technique[J].Shock and Vibration, 2016(4): 1-14.

[16] 錢管良.計算含裂紋圓截面梁振動的有限元法[J].機械強度, 1992, 14(1): 7-10.

[17] SOEDEL W.Vibrations of Shells and Plates(Third edition) [M].New York: Marcel Dekker Inc., 2004.

[18] ZHENG D Y, DU J T, LIU Y.Vibration characteristics analysis of an elastically restrained cylindrical shell with arbitrary thickness variation[J].Thin-Walled Structures, 2021, 165: 107930.

[19] 桑國光, 龔恢.船體梁抗彎能力的計算[J].中國造船, 1984, 25(1): 65-74.

[20] ISLAM M S, KHALIL G M.Dynamic analysis of a fixed ended beam with focus on vibration neutralization in ship hull structure[J].International Journal of Engineering Research and Technology, 2014, 3(1): 2050-2064.

Transverse Vibration Characteristics Analysis of Beam with General Boundary Conditions Based on Improved Fourier Series Method

HUANG Ying1, GUO Yong-mei2

（1.,450005,; 2.’,’710000,）

【Objective】Beam structure is widely used in various engineering occasions, such as ship engineering, ocean engineering, automobile engineering etc., and its vibration characteristics have a directs influence on the vibration control of the whole structure.【Method】The transverse vibration displacement of the beam structure is expanded by a modified Fourier series and solved by the Rayleigh-Ritz energy method to establish a unified solution model for the transverse vibration characteristics of the beam under general elastic boundary conditions.With the model developed, the influence of different boundary conditions on the vibration characteristics of beam system is studied and discussed.【Result and Conclusion】Based on the improved Fourier series, the beam vibration model established by Rayleigh-Ritz energy method is basically the same as the analytical solution, with good accuracy and surpasses the Ansys simulation solution.The study provides a very important reference for the research and vibration optimization of beam structures in ship and ocean engineering.

Euler-Bernoulli beam; transverse vibration; general boundaries; improved Fourier series method

U662.2

A

1673-9159（2022）02-0135-07

10.3969/j.issn.1673-9159.2022.02.017

2021-07-28

2019年度河南省高等職業(yè)學(xué)校青年骨干教師培養(yǎng)計劃項目（2019GZGG102）

黃櫻（1982―），女，碩士，副教授，研究方向為交通運輸規(guī)劃與管理。E-mail: topworker0@protonmail.com

黃櫻，郭詠梅.基于改進傅里葉級數(shù)法的任意邊界下梁橫振特性分析[J].廣東海洋大學(xué)學(xué)報，2022，42（2）：135-141.