聶 睿， 李天勻， 朱 翔， 陳 旭

(1.華中科技大學 船舶與海洋工程學院，武漢 430074； 2.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240；3.船舶與海洋水動力湖北省重點實驗室，武漢 430074)

以圓柱殼為主體的結(jié)構(gòu)廣泛存在于船舶與海洋工程、航空航天工程等領域中。其振動特性被國內(nèi)外學者進行了廣泛研究。受螺旋槳脈動力激勵而產(chǎn)生的推進軸系-艇體整體結(jié)構(gòu)振動是近年來水下航行器振動特性研究重點。推進軸系和艇體會相互傳遞振動，存在著明顯耦合效應，建立兩者耦合模型進行振動特性計算十分必要。但是在以往工作中，水下航行器推進軸系與殼體往往不考慮彼此的耦合作用，僅僅只探討了隔離后軸系或殼體的振動問題[1-5]。如欲對軸-艇耦合系統(tǒng)的振動特性進行更全面分析，對組合結(jié)構(gòu)振動這類問題開展更深入的研究勢在必行。

對于組合結(jié)構(gòu)，通常先將整個結(jié)構(gòu)劃分成多個子結(jié)構(gòu)，再利用結(jié)構(gòu)之間的連續(xù)協(xié)調(diào)條件體現(xiàn)結(jié)構(gòu)間耦合作用。周海軍等[6]利用改進傅里葉級數(shù)方法和波傳播法分別求取了梁和圓柱殼的導納，再利用阻抗綜合法，建立了耦合計算模型。Qu等[7]采用廣義變分法模擬梁、殼等結(jié)構(gòu)單元，再基于Hamilton原理，計算了梁-球殼-柱殼組合結(jié)構(gòu)振動響應。郭亮等[8]采用解析法對錐柱組合殼結(jié)構(gòu)及多跨梁結(jié)構(gòu)分別建模，并在軸承處建立橫向連續(xù)條件，形成了有效的系統(tǒng)橫向振動解析求解模型。李攀碩[9]分別利用解析法和有限元法計算子結(jié)構(gòu)導納，再通過頻響函數(shù)綜合法建立了軸與加筋殼體的耦合系統(tǒng)模型。Merz等[10]采用有限元和邊界元方法建立水下加筋殼體模型，并通過簡化的彈簧阻尼系統(tǒng)與用桿結(jié)構(gòu)模擬的推進軸系耦合，計算了系統(tǒng)振動與聲特性。

總的來說，已有文獻中軸-艇耦合系統(tǒng)振動分析方法主要有數(shù)值法和解析法兩種。有限元等數(shù)值方法可以用于分析各類耦合系統(tǒng)結(jié)構(gòu)振動問題，但它們難以從機理上對系統(tǒng)振動行為提供合理解釋。而解析法能給出較精準的解析結(jié)果，易于從機理上解釋某些研究現(xiàn)象，但其往往又存在求解范圍有限的缺點。而實際應用中存在大量復雜的邊界約束，解析法針對不同邊界需要配置不同位移函數(shù)和邊值條件，難以建立統(tǒng)一的動力學分析模型。

Rayleigh-Ritz法作為一種基于能量變分原理的(近似)解析法，由于其假定振型函數(shù)滿足結(jié)構(gòu)位移邊界條件和振動方程，由此得到的函數(shù)組合往往接近真實振型，因此有很高精度；同時，由于假定振型函數(shù)之間具有正交性，實際運算成本極低。此外，Rayleigh-Ritz法還十分容易處理附加質(zhì)量或附加剛度。只需添加附加動能項或附加應變能項，即可在不改變結(jié)構(gòu)原振型函數(shù)的前提下，求取新系統(tǒng)振動特性。但傳統(tǒng)Rayleigh-Ritz 法的限制是：對復雜邊界條件，難以選取恰當?shù)募僭O振型函數(shù)。Li[11]提出了一種改進傅里葉級數(shù)，該級數(shù)通過在傳統(tǒng)傅里葉級數(shù)增加輔助多項式函數(shù)，保證了梁結(jié)構(gòu)在邊界上的連續(xù)性，結(jié)合兩類彈簧模擬一般邊界的思想，成功建立了梁的自由振動計算模型。隨后，國內(nèi)外學者借鑒該研究思路，結(jié)合Rayleigh-Ritz法，研究了梁、板、殼等各類結(jié)構(gòu)的振動特性[12-15]。這些研究證明，引入改進傅里葉級數(shù)的Rayleigh-Ritz法，在面對復雜邊界條件時，可構(gòu)建準確統(tǒng)一的分析模型。

當對艇體中段進行截斷處理后，軸-艇耦合結(jié)構(gòu)可簡化為梁-圓柱殼耦合模型[16]。在對軸-艇系統(tǒng)進行合理簡化后，本文將以梁-圓柱殼結(jié)構(gòu)為研究對象，基于能量變分原理，采用改進傅里葉級數(shù)統(tǒng)一描述梁結(jié)構(gòu)和圓柱殼結(jié)構(gòu)的振動位移；通過人工約束彈簧組模擬各類邊界條件，耦合彈簧組來滿足梁殼之間的連續(xù)條件，建立了梁-圓柱殼耦合系統(tǒng)振動計算模型。如前文所述，本文方法是對Rayleigh-Ritz 法的改進，故所得計算模型不僅可用于圓柱殼取各類經(jīng)典邊界條件下的梁-圓柱殼耦合結(jié)構(gòu)，同時由于級數(shù)良好的收斂性亦只需極低運算成本即獲得準確結(jié)果。通過與有限元軟件的計算結(jié)果對比，驗證了本文計算模型正確性。最后，分析了系統(tǒng)在部分參數(shù)改變時，振動響應性能的變化規(guī)律。

1 梁-圓柱殼耦合系統(tǒng)振動計算模型

1.1 模型描述

當進行截斷處理后，軸-艇耦合系統(tǒng)可合理簡化為梁-圓柱殼結(jié)構(gòu)。即推進軸系可簡化為帶集中質(zhì)量點的彈性梁，艇體可簡化為兩端受彈簧約束的圓柱殼，中間軸承等結(jié)構(gòu)可簡化為連接梁結(jié)構(gòu)與殼結(jié)構(gòu)的耦合彈簧組。

梁-圓柱殼耦合模型如圖1所示，全局坐標系為直角坐標系O-xyz，坐標原點O位于梁結(jié)構(gòu)左端點處；圓柱殼的長度和截面半徑分別為L和R，殼體中面軸向位移、周向位移、徑向位移及繞y軸轉(zhuǎn)角分別為ucy、vcy、cy、βcy(βcy=?cy/?x)；梁通過兩處節(jié)點與圓柱殼連接，每段長度分別為l1和l2，梁的自由端有一附加集中質(zhì)量點，其質(zhì)量為m；kL和kR分別為左連接點與右連接點的空間點對點耦合彈簧組(由三向線位移彈簧kx、ky、kz和兩向轉(zhuǎn)角彈簧Ky、Kz構(gòu)成)；圓柱殼兩端引入兩組連續(xù)線分布的約束彈簧用于模擬各類邊界(由三種位移彈簧ku、kv、k和一種轉(zhuǎn)角彈簧Kβ構(gòu)成)；fx、fy、fz為作用在質(zhì)量點處的簡諧激勵力。

圖1 梁-圓柱殼耦合模型

1.2 理論推導

如前文所述，與直接求解結(jié)構(gòu)控制微分方程的方法不同，本文將基于能量變分原理建立梁-圓柱殼耦合模型。耦合結(jié)構(gòu)能量泛函Π的形式如下

Π=Tb-Ub+b+Tcy-Ucy-Uc

(1)

式中:Tb為梁結(jié)構(gòu)的總動能;Ub為梁結(jié)構(gòu)的總勢能;b為梁結(jié)構(gòu)上外部激勵力的外力功；Tcy為圓柱殼結(jié)構(gòu)的總動能;Ucy為圓柱殼結(jié)構(gòu)的總勢能；Uc為結(jié)構(gòu)耦合彈簧儲存的彈性勢能。

本文計算模型中梁結(jié)構(gòu)為含集中質(zhì)量點的彈性歐拉梁，其自身局部坐標系為xbybzb，坐標原點Ob位于梁左端，如圖2所示。

圖2 含集中質(zhì)量點的梁結(jié)構(gòu)

梁結(jié)構(gòu)的總動能Tb由梁結(jié)構(gòu)自身動能Tb,s和集中質(zhì)量點的動能Tb,m構(gòu)成

Tb=Tb，s+Tb，m

(2)

式中:ub，vb和b分別為梁結(jié)構(gòu)在縱向(x)，橫向(y)和垂向(z)的位移;ρb和Ab分別為梁的密度和截面面積;m為集中質(zhì)量點的質(zhì)量。集中質(zhì)量點附加在梁結(jié)構(gòu)左端面上，即梁結(jié)構(gòu)坐標系中原點處。

由于本文計算模型中梁結(jié)構(gòu)無外部約束，其勢能Ub僅由自身應變勢能構(gòu)成，形式如下

(3)

式中:Eb為梁材料的楊氏模量;Ib,y和Ib,z是梁的截面慣性矩。

集中質(zhì)量點處作用有縱向、橫向及垂向三個方向的簡諧點激勵力，由此產(chǎn)生的外力功b如下

b=fxub(xb,t)|xb=0+fyvb(xb,t)|xb=0+

fzb(xb,t)|xb=0

(4)

圓柱殼結(jié)構(gòu)采用柱坐標系(R,θ,xcy)，xcy和θ分別代表軸向和周向坐標,坐標原點Ocy位于圓柱殼左端面圓心處。圓柱殼結(jié)構(gòu)如圖3所示。

圖3 圓柱殼結(jié)構(gòu)示意圖

采用Ressiner薄殼理論計算圓柱殼的動能與應變能。動能Tcy具體形式如下

(5)

式中:ρcy為殼體密度;hcy為殼體厚度。

在圓柱殼兩端引入兩組沿端面圓周方向均勻分布的線性彈簧組用于模擬各類邊界條件。因此，圓柱殼的總勢能Ucy由自身應變能Ucy,p和存儲在邊界彈簧的彈性勢能Ucy,e組成。具體計算形式如下

Ucy=Ucy,p+Ucy,e

(6)

式中:Ecy和μcy分別為圓柱殼的楊氏模量和泊松比;ku0，kv0，k0和Kβ0為圓柱殼左端約束彈簧;ku1，kv1，k1和Kβ1為圓柱殼右端約束彈簧。連續(xù)線分布的位移彈簧和轉(zhuǎn)角彈簧的單位分別為N/m2和N/(rad·m)。

梁結(jié)構(gòu)和圓柱殼結(jié)構(gòu)將通過分布彈簧系統(tǒng)實現(xiàn)耦合。每對耦合節(jié)點之間的耦合彈簧組均包含縱(x)向、橫(y)向和垂(z)向三種線位移彈簧與y向及z向兩種轉(zhuǎn)角彈簧。圓柱殼的振動位移分量將從局部柱坐標轉(zhuǎn)換到全局笛卡爾坐標下，以便構(gòu)造耦合彈簧組的能量泛函。梁結(jié)構(gòu)和圓柱殼結(jié)構(gòu)的振動位移如圖4所示。

圖4 梁和圓柱殼的振動位移

vcy=vcycosθ-cysinθ,cy=cycosθ+vcysinθ

(7)

式中，vcy和cy為圓柱殼結(jié)構(gòu)轉(zhuǎn)化到全局坐標系下的振動線位移分量。

通常，薄殼理論在解決殼體問題時，殼體模型僅考慮三向線位移和繞y軸旋轉(zhuǎn)位移。但本文工作中為保證梁、殼結(jié)構(gòu)的耦合振動分量數(shù)目相當，需要再加入繞z軸旋轉(zhuǎn)位移αcy。引入額外旋轉(zhuǎn)位移并進行轉(zhuǎn)化后的形式如下

Rotcy,y=βcycosθ-αcysinθ

Rotcy,z=αcycosθ+βcysinθ

(8)

式中，Rotcy,y和Rotcy,z為圓柱殼結(jié)構(gòu)轉(zhuǎn)化到全局坐標系下的旋轉(zhuǎn)位移分量。

如前文所述，本文計算模型含kL、kR兩組梁殼耦合彈簧組。耦合彈簧存儲的彈性勢能Uc具體形式如下：

Uc=Uc,L+Uc,R

(9)

式中:Uc,L為左耦合彈簧組儲存的彈性勢能;kL,x，kL,y和kL,z分別代表x向，y向和z向的線位移耦合彈簧;KL,y和KL,z分別代表y向和z向的旋轉(zhuǎn)位移耦合彈簧;(xL，θL)為左耦合彈簧組上圓柱殼耦合節(jié)點在柱坐標系上的位置坐標。下標R代表右耦合彈簧組對應的變量，故不贅述。

本文將引入改進傅里葉級數(shù)作為梁結(jié)構(gòu)振動位移函數(shù)及圓柱殼結(jié)構(gòu)軸向振動位移函數(shù)。同時，由于完整的圓柱殼擁有軸對稱外形，其周向振動位移函數(shù)可由傳統(tǒng)的兩組單傅里葉級數(shù)描述。具體形式如下。

梁結(jié)構(gòu)：

(10)

式中:Mb為位移函數(shù)的截斷項數(shù);Am,u和Bl,u是梁軸向位移容許函數(shù)中的未知傅里葉系數(shù);eiωt被作為時間簡諧因子引入來描述梁在不同時刻的位移;Am,v、Bl,v、Am,和Bl,是梁橫向及垂向位移容許函數(shù)中的未知傅里葉系數(shù)。式(10)中的位移容許函數(shù)即為一種改進傅里葉級數(shù)，其中，λmb=mπx/(l1+l2)，λlb=lπx/(l1+l2)。它在傳統(tǒng)的單余弦傅里葉級數(shù)基礎上引入輔助正弦函數(shù)，使得容許函數(shù)滿足軸向振動位移一階導數(shù)或彎曲振動位移三階導數(shù)在梁上的任意一點存在且連續(xù)，保證滿足位移邊界條件，同時加速了級數(shù)自身的收斂速度。

殼結(jié)構(gòu)：

(11)

式中:Ncy、Mcy為位移函數(shù)截斷項數(shù);Amn等為未知傅里葉系數(shù)，λmcy=mπx/L，λlcy=lπx/L。

將式(1)～(11)聯(lián)立，并根據(jù)最小勢能原理，對能量泛函求取極值

(12)

式中，q為未知系數(shù)向量，q=[Am,u，Bl,u，…，Amn，al，…，Cmn，cl]。

最終得到的梁-圓柱殼耦合系統(tǒng)振動方程組矩陣化形式如下

([K]-ω2[M]){q}={F}

(13)

式中:[K]和[M]分別為耦合系統(tǒng)的剛度矩陣和質(zhì)量矩陣;ω為圓頻率;{F}為集中質(zhì)量點處激勵力的廣義力向量。廣義力向量取0時，亦可由式(13)求取系統(tǒng)的固有振動特性。

2 計算模型驗證

為驗證本文方法收斂性和計算模型準確性，下面將從傅里葉級數(shù)的截斷項數(shù)出發(fā)，分析方法收斂性。同時，將本文計算模型計算結(jié)果與有限元軟件計算結(jié)果進行比對，驗證本文計算模型正確性。有限元模型共有31 566個單元。其中，殼單元31 400個，梁單元165個，質(zhì)量點單元1個。

本文計算模型如前文所述，具體參數(shù)如表1所示。

表1 耦合模型參數(shù)

2.1 傅里葉級數(shù)截斷項數(shù)收斂性

改進傅里葉級數(shù)的截斷項數(shù)對計算結(jié)果精度有直接影響，其展開的數(shù)目越多，計算結(jié)果越趨于真實值，但過大的截斷項數(shù)會降低求解效率。一般來說，截斷項數(shù)達到一定值后，計算結(jié)果即可穩(wěn)定(收斂)。本文以圓柱殼為自由邊界條件(殼體兩端約束彈簧組取值為0)的梁-圓柱殼模型為例，探討計算結(jié)果收斂性隨Mb、Mcy的變化。而圓柱殼周向截斷數(shù)Ncy并不影響計算結(jié)果的收斂，僅僅會影響固定頻率范圍內(nèi)結(jié)果是否會缺少某些模態(tài)。

驗證模型參數(shù)見表1，Ncy=6，圓柱殼上的連接點為(xL，θL)=(0.5，π), (xR，θR)=(1.5，π)，梁的連接點為xb=l1和xb=l1+l2處。前8階固有頻率隨截斷項數(shù)Mb、Mcy的變化結(jié)果如表2所示。

表2 固有頻率與截斷項數(shù)關系

從表2可以看出，隨著截斷項數(shù)的增加，本文方法的計算結(jié)果很快就趨于收斂了，而與有限元法(finite element method,FEM)軟件計算結(jié)果的對比證明了本文方法有很高精度。綜合計算成本和精度，本文其余計算工作均取截斷項數(shù)為10。

2.2 彈簧模擬經(jīng)典邊界條件

如前所述，本文采用兩類彈簧模擬圓柱殼邊界。彈簧剛度增大，儲存的彈性勢能也越大，圓柱殼邊界處約束越強。但由剛度增大而產(chǎn)生的邊界約束效應必定會趨于完全約束(位移或轉(zhuǎn)角)。因此，取值達到一定大小的彈簧組可以用于模擬各類經(jīng)典邊界，例如固支邊界即應為所有彈簧組取足夠大值。下面將驗證該假設。

仍取前文計算模型，但假定ku0、kv0、k0、Kβ0值取0，ku1、kv1、k1、Kβ1值取Q。顯然，此時耦合模型圓柱殼左端邊界條件仍為自由。圖5將顯示耦合系統(tǒng)計算模型固有頻率隨彈簧剛度取值Q的變化關系。

圖5 邊界彈簧組剛度取值Q對左端自由系統(tǒng)固有頻率影響

從圖5可知，耦合計算模型的前10階固有頻率隨彈簧剛度增加都有著明顯的“緩慢增加—快速增加—穩(wěn)定不變”趨勢，這證明了前文彈簧剛度收斂性假設正確可靠。此外，雖然不同階固有頻率的曲線變化規(guī)律略有不同，但在剛度值取1×1012左右均可以達到穩(wěn)定。因此，本文其余工作均認為彈簧組剛度值取1×1012時，其對應方向上的殼體位移或轉(zhuǎn)動為完全約束。由于本文計算模型在y(橫)和z(垂)向是基本對稱的，因此前10階固有頻率為兩兩一組，成對出現(xiàn)。

下面假定ku0、kv0、k0取值為Q，Kβ0為0，ku1、kv1、k1、Kβ1取值為1×1012。此時，圓柱殼右端所有位移均被完全約束，即對應固支條件。圖6為此條件下耦合系統(tǒng)計算模型固有頻率隨彈簧剛度取值Q的變化關系。

圖6 邊界彈簧組剛度取值Q對右端固支系統(tǒng)固有頻率影響

從圖6可以看出，在右端為固支條件時，增加圓柱殼左端線彈簧組的剛度，耦合系統(tǒng)固有頻率也有類似的規(guī)律，彈簧剛度取值達到1×1012左右也趨于收斂。對應地，ku0、kv0、k0取1×1012且Kβ0取0，這一條件下，圓柱殼左端應視為簡支邊界條件。為進一步證明該假設的正確性，下面將這兩種邊界組合下的耦合系統(tǒng)固有頻率計算值與有限元計算結(jié)果進行對比。結(jié)果如表3所示。

表3 本文方法計算結(jié)果與有限元結(jié)果對比

從表3的結(jié)果可以看出，本文方法和FEM計算結(jié)果十分接近，使用彈簧模擬邊界的計算方法是正確可靠的。

2.3 計算模型正確性

在截斷項數(shù)收斂性和彈簧剛度收斂性的探討中，本文計算模型正確性已經(jīng)從固有頻率的比對中得到了初步驗證。下面將進一步從對應模態(tài)振型及受迫振動響應的對比角度出發(fā)，驗證本文計算模型正確性。

這里的驗證模型取圓柱殼邊界條件為簡支-簡支組合，其余參數(shù)保持不變。首先，對比梁結(jié)構(gòu)z(垂)向振動的前5階模態(tài)振型，結(jié)果如表4所示。

從表4可以看出，本文計算模型的梁結(jié)構(gòu)z向振動前2階模態(tài)，圓柱殼基本可視為剛體結(jié)構(gòu)，不參與系統(tǒng)的振動。進一步觀察得知，這2階模態(tài)中，耦合點處于梁結(jié)構(gòu)振型的不動點處，因此梁結(jié)構(gòu)和圓柱殼結(jié)構(gòu)難以通過耦合點進行振動傳遞。這也可以說明圖6中，為何圓柱殼約束增強，但耦合系統(tǒng)前4階(包括2階垂向和2階橫向)振動的固有頻率幾乎沒有變化。因為此時梁結(jié)構(gòu)與圓柱殼結(jié)構(gòu)可以視為解耦狀態(tài)，圓柱殼可等效為剛性基礎，其約束的增強并不能對系統(tǒng)產(chǎn)生明顯影響。而后面的模態(tài)中，圓柱殼結(jié)構(gòu)參與了整個系統(tǒng)的振動，因此其約束的增強會使得整個系統(tǒng)勢能增加，固有頻率也會隨之增大。

然后，對比本文方法和FEM的受迫振動響應。在集中質(zhì)量點處加z向單位幅值激勵力，掃頻范圍為1～200 Hz，提取梁結(jié)構(gòu)在質(zhì)量點、左耦合點和右耦合點處的z向位移響應及殼結(jié)構(gòu)在左耦合點和右耦合點處的z向位移響應。取參考位移A0=1×10-12m，A為提取點的z向位移響應幅值，則振動位移級響應為P=20lg(A/A0)。

從圖7可以看出，本文計算模型與FEM軟件計算結(jié)果十分吻合，這進一步驗證本文方法的正確性。在190 Hz附近兩種方法得到的位移級響應峰值位置存在微小偏差，表明兩者計算得到的系統(tǒng)固有頻率存在差別，這也與表3中內(nèi)容相符。

(a) 梁結(jié)構(gòu)

3 參數(shù)對耦合系統(tǒng)振動特性影響

3.1 耦合彈簧組剛度對系統(tǒng)振動特性影響

由前文理論推導部分可知，耦合彈簧組儲存的彈性勢能是構(gòu)成整個系統(tǒng)的重要部分。顯然，耦合彈簧剛度值將直接影響其儲存的彈性勢能大小，進一步影響耦合系統(tǒng)的振動特性。

在前文的計算中，耦合彈簧組選取了較大的耦合剛度值1×108。參照邊界彈簧收斂性討論部分的內(nèi)容，這種情況下梁結(jié)構(gòu)和圓柱殼結(jié)構(gòu)耦合點之間應該為一種“強耦合”的狀態(tài)，即梁節(jié)點與圓柱殼節(jié)點的對應位移嚴格相等。當耦合彈簧的剛度減弱時，耦合點之間可能會轉(zhuǎn)化為一種“弱耦合”的狀態(tài)。當這種狀態(tài)發(fā)生的時候，耦合系統(tǒng)振動特性應該會產(chǎn)生明顯的變化。下面將保持系統(tǒng)其余參數(shù)不變且圓柱殼取簡支-簡支邊界條件，通過改變耦合彈簧組的線彈簧剛度取值，驗證這一猜想。令kL,x=0，其余各耦合線彈簧剛度采用相同取值。

從圖8中可以看出，與邊界彈簧隨剛度值增加的變化規(guī)律類似，耦合彈簧剛度值增大到一定量之后，系統(tǒng)固有振動特性將會趨于穩(wěn)定。此時，耦合剛度的改變對系統(tǒng)振動性能影響可以忽略不計，具體表現(xiàn)為剛度取值為1×108和1×1010的響應曲線有著高度重合性。而前文中關于“強弱”耦合狀態(tài)轉(zhuǎn)變時耦合系統(tǒng)將產(chǎn)生振動特性變化的猜想也在圖中得到了明顯體現(xiàn)。系統(tǒng)振動特性的變化可以分為兩個方面。首先，從圖8和圖9(a)、9(b)可以看出，在“弱耦合”階段，梁結(jié)構(gòu)與殼結(jié)構(gòu)在對應耦合點的位移級響應曲線雖然變化趨勢基本吻合但存在明顯“落差”，說明由梁結(jié)構(gòu)質(zhì)量點處激勵力產(chǎn)生的振動在耦合點處產(chǎn)生了明顯衰減；隨著耦合剛度增加，達到“強耦合”階段后，對應耦合點處兩者位移級響應曲線基本重合了，表明此狀態(tài)下兩者振動傳遞幾乎沒有損失。其次，在耦合剛度較小的“弱耦合”階段，隨著耦合彈簧剛度增加，系統(tǒng)的勢能增大，固有頻率整體有增大的趨勢，在圖8中表現(xiàn)為受迫振動位移響應共振峰整體右移，直至達到“強耦合”階段后，共振峰不再產(chǎn)生明顯移動。

圖8 耦合彈簧剛度取值變化時梁及殼結(jié)構(gòu)左耦合點處

圖9 不同耦合彈簧剛度值下梁及殼結(jié)構(gòu)耦合點處位移級響應對比

3.2 旋轉(zhuǎn)耦合彈簧對系統(tǒng)振動特性影響

以上工作僅考慮梁-殼結(jié)構(gòu)之間的線位移耦合，但實際上兩者之間還存在耦合彎矩作用，即還應考慮本文耦合彈簧組中旋轉(zhuǎn)彈簧的影響。仍然取圓柱殼為簡支-簡支條件下的計算模型，而各耦合彈簧如存在則剛度值取1×108。表5為各耦合彈簧對系統(tǒng)固有振動特性影響。

如表5所示，如考慮梁-殼結(jié)構(gòu)之間的耦合彎矩作用，即添加旋轉(zhuǎn)耦合彈簧后，由于耦合彈性勢能的增加，系統(tǒng)的固有頻率將會增大。此外，由于本文模型在垂(z)向及橫(y)向上無耦合，因此耦合旋轉(zhuǎn)彈簧也僅會對對應方向上的振動產(chǎn)生影響。

表5 不同耦合彈簧組合下系統(tǒng)前10階固有頻率

3.3 集中質(zhì)量點質(zhì)量對系統(tǒng)振動特性影響

由本文計算模型推導部分可知，附加集中質(zhì)量的大小將直接影響梁結(jié)構(gòu)的附加動能項，進而引起系統(tǒng)振動性能變化。下面分別取耦合彈簧剛度值為1×104和1×108且圓柱殼為簡支-簡支條件下的計算模型，研究僅改變附加質(zhì)量大小時耦合系統(tǒng)振動性能變化。

從圖10(a)及10(b)中可以看出，在本文探討的變化范圍內(nèi)，集中質(zhì)量點質(zhì)量的增加對梁結(jié)構(gòu)的振動位移響應有一定的抑制作用，響應曲線有下移的趨勢；此外，當集中點質(zhì)量增大時，附加動能項能量增加，系統(tǒng)固有頻率減小，在圖中反映為共振峰的左移。實際上，本文質(zhì)量點為螺旋槳的簡化，其質(zhì)量大小的改變可視為由螺旋槳與水耦合作用產(chǎn)生的附連水質(zhì)量效應，而本文計算結(jié)果也與附連水導致含槳系統(tǒng)固有頻率降低的實際情況相符。

圖10 不同耦合剛度值下集中質(zhì)量點質(zhì)量變化對結(jié)構(gòu)左耦合點處位移響應影響

4 結(jié) 論

本文基于Rayleigh-Ritz法，引入改進傅里葉級數(shù)作為結(jié)構(gòu)位移函數(shù)，結(jié)合邊界彈簧組和耦合彈簧組，建立了梁-圓柱殼耦合動力學計算模型。通過改變邊界彈簧組剛度值即可模擬各種經(jīng)典圓柱殼邊界條件，構(gòu)建了適用范圍廣的動力學分析模型。與有限元法(FEM)軟件計算結(jié)果的對比，驗證了本文方法在計算固有特性及受迫振動響應時的正確性。隨后，進一步分析了耦合剛度值和集中點質(zhì)量等參數(shù)對梁-圓柱殼耦合系統(tǒng)振動特性的影響，得出了以下結(jié)論：

(1) 針對本文研究的梁-圓柱殼耦合系統(tǒng)，耦合系統(tǒng)的前4階模態(tài)可視為支撐于剛體上的梁結(jié)構(gòu)彎曲振動，梁結(jié)構(gòu)與圓柱殼結(jié)構(gòu)可視為解耦狀態(tài)。此時，殼體邊界約束條件變化對系統(tǒng)固有頻率幾乎沒有影響。

(2) 耦合彈簧剛度值處于較小的范圍時，由于對耦合節(jié)點處對應方向的連續(xù)性約束較弱，梁殼耦合結(jié)構(gòu)將在耦合節(jié)點處出現(xiàn)位移不連續(xù)的現(xiàn)象。此時，兩者的振動位移響應曲線雖然趨勢相同，但響應幅值出現(xiàn)落差，表明振動在結(jié)構(gòu)中的傳遞產(chǎn)生了衰減。

(3) 此外，由于集中質(zhì)量點附加質(zhì)量對結(jié)構(gòu)振動動能項的影響，梁端點處的集中質(zhì)量點質(zhì)量對耦合系統(tǒng)振動特性有一定的影響，系統(tǒng)會出現(xiàn)一定的固有頻率改變(響應曲線共振峰偏移)，對應測點處的振動位移響應結(jié)果也會產(chǎn)生一定的變化。