基于新型降階模型的網(wǎng)格加筋圓柱殼頻率分析

李玉韋,王 博,程耿東

(大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室工程力學(xué)系,大連 116024)

基于新型降階模型的網(wǎng)格加筋圓柱殼頻率分析

李玉韋,王 博,程耿東

(大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室工程力學(xué)系,大連 116024)

針對(duì)網(wǎng)格加筋圓柱殼頻率分析計(jì)算量大的問(wèn)題，提出一種基于多項(xiàng)式及梁?jiǎn)卧魏瘮?shù)的模型降階方法，即運(yùn)用多項(xiàng)式及梁?jiǎn)卧魏瘮?shù)將復(fù)雜結(jié)構(gòu)有限元的節(jié)點(diǎn)位移轉(zhuǎn)化為主節(jié)點(diǎn)位移，實(shí)現(xiàn)模型降階。并通過(guò)控制多項(xiàng)式階數(shù)及梁?jiǎn)卧獋€(gè)數(shù)來(lái)調(diào)整降階模型的分析精度，而降階模型的計(jì)算效率較精細(xì)模型有顯著提高。以網(wǎng)格加筋圓柱殼為例，對(duì)比本文提出的模型降階方法與3-D實(shí)體模型及其他模型降階方法的頻率分析結(jié)果，結(jié)果表明，提出的模型降階方法不僅能捕捉到結(jié)構(gòu)的整體模態(tài)，還能反映出結(jié)構(gòu)部分局部模態(tài)，適用范圍廣，能夠?yàn)楣こ探Y(jié)構(gòu)設(shè)計(jì)提供簡(jiǎn)單有效的計(jì)算模型。

多項(xiàng)式函數(shù)；梁?jiǎn)卧魏瘮?shù)；模型降階；網(wǎng)格加筋圓柱殼；頻率分析

0 引言

網(wǎng)格加筋圓柱殼由于較高的比剛度和比強(qiáng)度廣泛應(yīng)用于航空航天結(jié)構(gòu)，如運(yùn)載火箭的燃料貯箱和級(jí)間段等[1-2]。運(yùn)載火箭結(jié)構(gòu)的彎曲、扭轉(zhuǎn)、縱向等整體模態(tài)及部分局部模態(tài)是載荷計(jì)算、姿態(tài)穩(wěn)定性分析、速率陀螺儀安放位置確定等工作的基礎(chǔ)。詳細(xì)描述結(jié)構(gòu)細(xì)節(jié)的有限元模型往往具有幾十萬(wàn)甚至上千萬(wàn)個(gè)自由度，幾乎不可能采用這樣的模型進(jìn)行結(jié)構(gòu)動(dòng)力特性分析，而且基于此類有限元模型得到的頻率分析結(jié)果中含有大量的局部模態(tài)，反而干擾了對(duì)結(jié)構(gòu)總體性能的認(rèn)識(shí)。因此，在結(jié)構(gòu)設(shè)計(jì)初期，選取合適的計(jì)算模型是開展相關(guān)分析的前提。

為提高網(wǎng)格加筋圓柱殼頻率分析的計(jì)算效率，邢譽(yù)峰等[3]、潘忠文等[4]研究了加筋圓柱殼彎曲剛度的面積等效方法和慣性矩等效方法，分別分析了面積等效時(shí)截面慣性矩的精度和慣性矩等效時(shí)面積的等效精度，提出了用于加筋圓柱殼頻率分析的等效厚度模型，并指出等效厚度模型可直接用于計(jì)算結(jié)構(gòu)的彎曲和縱向頻率，但扭轉(zhuǎn)頻率計(jì)算精度較差，需進(jìn)行修正。Cheng等[5]、Cai等[6-7]通過(guò)改進(jìn)漸近均勻化方法(NIAH)預(yù)測(cè)三維周期性結(jié)構(gòu)剛度特性，將復(fù)雜結(jié)構(gòu)等效為具有類似剛度特性的簡(jiǎn)單結(jié)構(gòu)。這種方法不需要復(fù)雜的數(shù)學(xué)推導(dǎo)，可以利用商業(yè)軟件中提供的各種單元和建模技術(shù)，使單胞模型保持較小的規(guī)模。鄭淑飛等[8-9]提出一種基于變形修正的動(dòng)力模型簡(jiǎn)化方法，該方法以加筋圓柱殼精細(xì)有限元模型為基礎(chǔ)，根據(jù)運(yùn)動(dòng)同步性假設(shè)將結(jié)構(gòu)分成若干同步性區(qū)域，建立結(jié)構(gòu)整體位移與局部位移模式之間的轉(zhuǎn)換關(guān)系，實(shí)現(xiàn)模型降階，但位移轉(zhuǎn)換陣中未考慮原結(jié)構(gòu)的轉(zhuǎn)角自由度，對(duì)含有梁、殼等元組件結(jié)構(gòu)的分析精度較差。在此基礎(chǔ)上，王文勝等[10，12]、Cheng等[11]將復(fù)雜結(jié)構(gòu)模型劃分為不同的梁段，通過(guò)梁平截面假設(shè)及位移插值函數(shù)將每一梁段降階為一個(gè)梁超單元，按照一定的次序組裝成超梁模型，并引入橫向剪切系數(shù)修正梁超單元的剛度陣。針對(duì)具有大開口的結(jié)構(gòu)，本文通過(guò)數(shù)值方法引入了表示截面翹曲變形的翹曲基向量，進(jìn)一步提高了扭轉(zhuǎn)頻率的計(jì)算精度，本文稱該方法為平截面梁模型方法。文獻(xiàn)[12]通過(guò)網(wǎng)格加筋結(jié)構(gòu)、含開口光筒結(jié)構(gòu)驗(yàn)證了平截面梁模型方法的有效性，但該方法僅適用于計(jì)算結(jié)構(gòu)的整體頻率，無(wú)法獲得結(jié)構(gòu)的局部模態(tài)。

本文在基于多項(xiàng)式函數(shù)的動(dòng)力模型降階方法的基礎(chǔ)上，提出一種基于多項(xiàng)式及梁?jiǎn)卧魏瘮?shù)混合的模型降階方法，該方法將結(jié)構(gòu)沿軸向劃分若干個(gè)梁?jiǎn)卧?，利用多?xiàng)式及梁?jiǎn)卧魏瘮?shù)建立所有節(jié)點(diǎn)位移與梁?jiǎn)卧?jié)點(diǎn)位移之間的轉(zhuǎn)換關(guān)系，實(shí)現(xiàn)結(jié)構(gòu)降階。建立的降階模型不僅可用于計(jì)算結(jié)構(gòu)的整體頻率，還可通過(guò)增加多項(xiàng)式項(xiàng)數(shù)及梁?jiǎn)卧獋€(gè)數(shù)獲得結(jié)構(gòu)的部分局部模態(tài)。文中首先以網(wǎng)格加筋圓柱殼為例，對(duì)比本文提出的模型降階方法與3-D實(shí)體模型及其他模型降階方法的頻率分析結(jié)果。然后，針對(duì)含開口網(wǎng)格加筋圓柱殼結(jié)構(gòu)，文中通過(guò)數(shù)值方法引入翹曲變形改善降階模型扭轉(zhuǎn)頻率的計(jì)算精度。數(shù)值結(jié)果表明，本文提出的模型降階方法計(jì)算精度較高。

1 模型降階方法

模型降階是利用原結(jié)構(gòu)有限元模型的數(shù)學(xué)或力學(xué)特征，選擇一組合適的減縮基向量，將復(fù)雜結(jié)構(gòu)動(dòng)力學(xué)問(wèn)題降階。不考慮阻尼影響，原結(jié)構(gòu)動(dòng)力分析的有限元控制方程可表示為：

(1)

引入位移轉(zhuǎn)換矩陣T(r×k)(r?k)，復(fù)雜結(jié)構(gòu)的位移可表示為：

U=TQ

(2)

其中，Q(k×1)是降階模型的位移。將式(2)帶入式(1)，并左乘TT，得到降階模型的動(dòng)力學(xué)方程：

(3)

其中，MR=TTMT,KR=TTKT是降階模型的質(zhì)量陣和剛度陣，F(xiàn)R=TTF是降階模型的載荷向量。通過(guò)求解特征方程(4)獲得降階模型的固有頻率。

(KR-λkkMR)φkk=0

(4)

其中，λkk為結(jié)構(gòu)第kk階頻率，φkk為與其對(duì)應(yīng)的特征向量。

不同模型降階方法的主要區(qū)別在于減縮基向量T的構(gòu)造，不同的減縮基構(gòu)造方法可能導(dǎo)致不同的計(jì)算精度及計(jì)算效率，文獻(xiàn)[13-16]給出了不同的模型降階方法，此處不再贅述。

2 基于多項(xiàng)式及梁?jiǎn)卧魏瘮?shù)的模型降階方法

考慮如圖1所示的復(fù)雜結(jié)構(gòu)，OXYZ為全局坐標(biāo)系，將結(jié)構(gòu)沿軸向劃分為若干個(gè)梁?jiǎn)卧∶總€(gè)梁?jiǎn)卧墓?jié)點(diǎn)為主節(jié)點(diǎn)。首先利用多項(xiàng)式建立任意一點(diǎn)與其梁軸線上的投影點(diǎn)之間的位移轉(zhuǎn)換關(guān)系，再通過(guò)梁?jiǎn)卧魏瘮?shù)建立投影點(diǎn)與主節(jié)點(diǎn)之間的位移轉(zhuǎn)換關(guān)系，通過(guò)兩次位移轉(zhuǎn)換可以建立梁?jiǎn)卧獌?nèi)任意節(jié)點(diǎn)位移與其主節(jié)點(diǎn)位移之間的轉(zhuǎn)換關(guān)系，進(jìn)而得到用于模型降階的減縮基向量。具體步驟如下：

結(jié)構(gòu)上任意一點(diǎn)與其投影點(diǎn)之間的位移轉(zhuǎn)換關(guān)系如式(5)所示：

uj=Rjqi

(5)

展開形式為：

(6)

ujx=qix1+yqix2+zqix3+y2qix4+yzqix5+z2qix6=Fτjqixτ

ujy=qiy1+yqiy2+zqiy3+y2qiy4+yzqiy5+z2qiy6=Fτjqiyτ

ujz=qiz1+yqiz2+zqiz3+y2qiz4+yzqiz5+z2qiz6=Fτjqizτ

(7)

位移轉(zhuǎn)換矩陣Rj為：

(8)

投影點(diǎn)位移與其所在梁?jiǎn)卧鞴?jié)點(diǎn)位移的轉(zhuǎn)換關(guān)系為：

(9)

表2 梁?jiǎn)卧魏瘮?shù)

由式(5)及式(9)，可知結(jié)構(gòu)上任意一點(diǎn)與其所在梁?jiǎn)卧鞴?jié)點(diǎn)之間的位移轉(zhuǎn)換陣為：

(10)

梁?jiǎn)卧猧內(nèi)所有節(jié)點(diǎn)的位移與主節(jié)點(diǎn)的位移轉(zhuǎn)換關(guān)系為：

例如：在學(xué)習(xí)乙酸乙酯的制備時(shí)，通過(guò)動(dòng)畫模擬的技術(shù)可以清晰直觀的把乙酸和乙醇反應(yīng)時(shí)的斷鍵、成鍵部位展示出來(lái)，學(xué)生會(huì)對(duì)這類多原子的化合物的反應(yīng)機(jī)理產(chǎn)生更深入的了解，提高了化學(xué)實(shí)驗(yàn)教學(xué)的教學(xué)效率

(11)

(12)

3 數(shù)值算例

通過(guò)兩個(gè)具體算例的頻率分析，比較本文提出的模型降階方法與其他模型降階方法的分析結(jié)果，證明本方法的有效性。

3.1 網(wǎng)格加筋圓柱殼算例

算例1為圖2所示的正置正交網(wǎng)格加筋圓柱

殼，其中，橫向筋條11根，軸向筋條4根。筒殼長(zhǎng)10020mm，外半徑為1000mm，蒙皮厚10mm，筋條高度為50mm，筋條厚20mm。采用鋁合金材料，彈性模量為73GPa，泊松比為0.3，密度為2.7×10-3g/mm3。采用3-D實(shí)體單元(solid45)建立有限元模型，并沿筋條高度及厚度方向局部加密，整體網(wǎng)格密度為100mm，共有17360個(gè)節(jié)點(diǎn)，52080個(gè)自由度。

為對(duì)比不同梁?jiǎn)卧獋€(gè)數(shù)對(duì)降階模型計(jì)算精度的影響，本文利用B4梁?jiǎn)卧⑼队包c(diǎn)位移與其所在梁?jiǎn)卧鞴?jié)點(diǎn)之間的位移轉(zhuǎn)換關(guān)系，采用3階多項(xiàng)式(N=3)建立截面任意一點(diǎn)位移與投影點(diǎn)位移的轉(zhuǎn)化關(guān)系。表3給出了固支-自由邊界條件下的不同梁?jiǎn)卧獋€(gè)數(shù)(Nm)對(duì)網(wǎng)格加筋圓柱殼結(jié)構(gòu)整體頻率的影響。

表3 固支-自由狀態(tài)下不同梁?jiǎn)卧獋€(gè)數(shù)對(duì)整體頻率的影響

同樣以圖2所示網(wǎng)格加筋柱殼為例，表4討論了多項(xiàng)式階數(shù)(N)對(duì)降階模型頻率計(jì)算精度的影響，假定降階模型沿軸向方向劃分20個(gè)B4梁?jiǎn)卧?。為說(shuō)明本文提出方法的有效性，采用模態(tài)置信因子對(duì)降階模型和原結(jié)構(gòu)的模態(tài)作比較。模態(tài)置信因子MAC[17]定義為：

(13)

其中，φA代表原模型計(jì)算得到的模態(tài)，φR代表降階模型的模態(tài)。一般情況下，當(dāng)MAC值大于0.7時(shí)，認(rèn)為二者相似性較好；當(dāng)MAC值小于0.5時(shí)，相似性較差。

表4 固支-自由狀態(tài)下多項(xiàng)式階數(shù)對(duì)降階模型頻率的影響(單位：Hz)

由表4可知，降階模型的計(jì)算精度能滿足工程設(shè)計(jì)的需要，得到的整體模態(tài)與原模型的模態(tài)具有極高的相似性，降階模型合理地保留了原模型的特征，計(jì)算效率有很大提高。此外，本文提出的模型降階方法不僅能捕捉到原模型的整體模態(tài)，還能獲得原結(jié)構(gòu)某些局部模態(tài)，如當(dāng)多項(xiàng)式階數(shù)N=7時(shí)，降階模型得到的局部頻率計(jì)算誤差為7.05%，計(jì)算時(shí)間為原精細(xì)模型的43%。進(jìn)一步提高多項(xiàng)式階數(shù)至N=8時(shí)，局部頻率的計(jì)算誤差為5.21%，計(jì)算時(shí)間為精細(xì)模型的63%。圖3～圖6給出部分模態(tài)對(duì)比圖。因此，本文提出的模型降階方法進(jìn)一步拓寬了降階模型的適用范圍，為設(shè)計(jì)人員提供了一個(gè)簡(jiǎn)單有效的計(jì)算模型。

3.2 含開口網(wǎng)格加筋圓柱殼算例

為滿足適用性和維修性要求，在運(yùn)載火箭結(jié)構(gòu)上要設(shè)計(jì)各種開口，含開口網(wǎng)格加筋圓柱殼發(fā)生扭轉(zhuǎn)變形時(shí)，截面上各點(diǎn)沿軸線方向有不同的位移，產(chǎn)生了翹曲變形，Vlasov[18]提出了薄壁桿件約束扭轉(zhuǎn)時(shí)截面變形規(guī)律，認(rèn)為橫截面上各點(diǎn)的翹曲軸向變形沿桿件橫截面按各點(diǎn)的扇形坐標(biāo)分布，而對(duì)于復(fù)雜結(jié)構(gòu)，計(jì)算扇形坐標(biāo)是非常困難的。本文通過(guò)數(shù)值方法[12]實(shí)現(xiàn)截面變形的精確描述，該方法通過(guò)對(duì)模型施加單位扭矩，利用靜力分析計(jì)算得到各個(gè)截面上的節(jié)點(diǎn)位移，進(jìn)而通過(guò)Schmidt正交化引入翹曲基向量。具體步驟如下：

梁?jiǎn)卧獌?nèi)所有節(jié)點(diǎn)的變形可表示為：

(14)

(15)

梁?jiǎn)卧獌?nèi)各主節(jié)點(diǎn)的位移分量可表示為：

(16)

(17)

(18)

將式(18)代入式(12)得到考慮翹曲變形的位移轉(zhuǎn)換陣T。

算例2為圖7所示的含開口網(wǎng)格加筋圓柱殼，結(jié)構(gòu)尺寸及筋條分布情況和算例1相同，矩形開口位于軸向1/2位置處，開口長(zhǎng)度L=1980mm，寬度H=692.8mm，采用鋁合金材料，彈性模量為73GPa，泊松比為0.3，密度為2.7×10-3g/mm3。有限元模型共有節(jié)點(diǎn)16912個(gè)，自由度為50736個(gè)。

表5給出了引入翹曲基向量前后的降階模型頻率計(jì)算結(jié)果與ANSYS計(jì)算結(jié)果的對(duì)比，其中第1、2、4、5階為彎曲模態(tài)，第3階為扭轉(zhuǎn)模態(tài)，降階模型軸向仍采用20個(gè)B4梁?jiǎn)卧?/p>

表5 固支-自由狀態(tài)下含開口網(wǎng)格加筋圓柱殼降階模型頻率計(jì)算結(jié)果(單位：Hz)

由表5可看出，僅通過(guò)增加多項(xiàng)式階數(shù)可提高降階模型的計(jì)算精度，當(dāng)N=7時(shí)，彎曲頻率的計(jì)算誤差最高為3.07%，但扭轉(zhuǎn)頻率的計(jì)算誤差仍高達(dá)18.21%，這是由于降階模型不能精確描述截面的扭轉(zhuǎn)變形，導(dǎo)致扭轉(zhuǎn)剛度增大。引入截面翹曲基向量后，扭轉(zhuǎn)頻率的計(jì)算誤差降為0.56%，降階模型的計(jì)算精度大幅提高，計(jì)算時(shí)間為精細(xì)模型的34.1%。圖8給出了考慮翹曲變形前后降階模型得到的扭轉(zhuǎn)模態(tài)，對(duì)比圖中紅色區(qū)域的變形可看出，引入翹曲基向量的降階模型能更好地描述具有大開口結(jié)構(gòu)的實(shí)際變形。

4 結(jié)論

本文針對(duì)網(wǎng)格加筋圓柱殼結(jié)構(gòu)提出一種新的模型降階方法，該方法將復(fù)雜結(jié)構(gòu)劃分為若干個(gè)梁?jiǎn)卧瑢⒚總€(gè)梁?jiǎn)卧墓?jié)點(diǎn)作為主節(jié)點(diǎn)，利用多項(xiàng)式函數(shù)建立結(jié)構(gòu)任意節(jié)點(diǎn)與其所在軸線上投影點(diǎn)之間的位移轉(zhuǎn)換關(guān)系，再利用梁?jiǎn)卧魏瘮?shù)建立投影點(diǎn)與主節(jié)點(diǎn)之間的位移轉(zhuǎn)換關(guān)系，進(jìn)而實(shí)現(xiàn)模型降階。針對(duì)含開口網(wǎng)格加筋結(jié)構(gòu)，文中通過(guò)數(shù)值方法引入翹曲變形，提高降階模型扭轉(zhuǎn)頻率的計(jì)算精度。

本文提出的模型降階方法不需要大量的矩陣運(yùn)算，計(jì)算效率高，可利用商業(yè)軟件作為一個(gè)黑箱來(lái)實(shí)現(xiàn)，并通過(guò)調(diào)整多項(xiàng)式階數(shù)及梁?jiǎn)卧膫€(gè)數(shù)控制降階模型的計(jì)算精度。文中首先以網(wǎng)格加筋柱殼為例驗(yàn)證該降階方法的有效性，計(jì)算結(jié)果表明該降階模型合理地保留了原結(jié)構(gòu)的特征，不僅能用于計(jì)算結(jié)構(gòu)的整體頻率，還可捕捉到結(jié)構(gòu)的某些局部模態(tài)。其次，文中將該方法應(yīng)用于含開口網(wǎng)格加筋柱殼的頻率分析中，通過(guò)數(shù)值方法引入翹曲變形提高降階模型扭轉(zhuǎn)頻率的計(jì)算精度。因此，本文提出的模型降階方法擴(kuò)展了降階模型的應(yīng)用范圍，為設(shè)計(jì)人員提供了一個(gè)簡(jiǎn)單有效的計(jì)算模型。

本文提出的模型降階方法可以很方便建立不同約束條件下降階模型，在此基礎(chǔ)上，研究多工況、多邊界條件下結(jié)構(gòu)的動(dòng)力特性分析及動(dòng)力優(yōu)化設(shè)計(jì)是下一步開展的工作。

[1] Hao P, Wang B, Tian K, et al. Efficient optimization of cylindrical stiffened shells with reinforced cutouts by curvilinear stiffeners[J]. AIAA Journal, 2016, 54 (4):1-14.

[2] Wang B, Tian K, Hao P, et al. Numerical-based smeared stiffener method for global buckling analysis of grid-stiffened composite cylindrical shells[J]. Composite Structures, 2016, 15(2): 807-815.

[3] 邢譽(yù)峰, 潘忠文,楊陽(yáng).蒙皮加筋圓柱殼彎曲頻率的三種計(jì)算模型[J].北京航空航天大學(xué)學(xué)報(bào), 2012, 38 (4):438-443.

[4] 潘忠文,王小軍,馬興瑞，等. 基于梁模型的蒙皮加筋結(jié)構(gòu)縱橫扭一體化建模研究[J].中國(guó)科學(xué): 技術(shù)科學(xué), 2014, 44 (5) :517-524.

[5] Cheng G D ,Cai Y W, Xu L. Novel implementation of homogenization method to predict effective properties of periodic materials[J]. Acta Mechanica Sinica, 2013, 29 (4):550-556.

[6] Cai Y W, Xu L, Cheng G D. Novel numerical implementation of asymptotic homogenization method for periodic plate structures[J]. International Journal of Solids and Structures, 2014, 51 (1) :284-292.

[7] 蔡園武. 周期性板結(jié)構(gòu)的漸近均勻化方法及微結(jié)構(gòu)優(yōu)化 [D]. 大連：大連理工大學(xué). 2014.

[8] Zheng S F, Ding H. A model reduction method for dynamic analysis based on deformation modification and local rigid body mode[J]. Mechanics in Engineering, 2008, 30 (9):31-35.

[9] 鄭淑飛,丁樺. 基于變形修正的局部剛體化動(dòng)力模型簡(jiǎn)化方法[J].力學(xué)與實(shí)踐, 2008, 30 (6) :31-34.

[10] 王文勝,程耿東,郝鵬. 基于超梁降階模型的蒙皮加筋圓柱殼頻率分析[J].固體火箭技術(shù), 2015, 38 (3) :401-406.

[11] Cheng G D, Wang W S. Fast dynamic analysis of complicated beam-type structure based on reduced super beam model[J]. AIAA Journal, 2014, 52 (5) :952-963.

[12] Wang W S, Cheng G D, Li Q H. Fast dynamic performance optimization of complicated beam-type structures based on two new reduced physical models[J]. Engineering Optimization, 2013, 45 (7):835-850.

[13] Guyan R J. Reduction of stiffness and mass matrices[J]. AIAA Journal, 1965, 3 (2) :380.

[14] Irons B. Structural eigenvalue problems-elimination of unwanted variables[J]. AIAA Journal, 1965, 3 (5):961-962.

[15] Carrera E, Miglioretti F, Petrolo M. Computations and evaluations of higher-order theories for free vibration analysis of beams[J]. Journal of Sound and Vibration, 2012, 331 (19) :4269-4284.

[16] Carrera E, Zappino E, Filippi M. Free vibration analysis of thin-walled cylinders reinforced with longitudinal and transversal stiffeners[J]. Journal of Vibration and Acoustics, 2013, 135(1):011019.

[17] Allemang R J, Brown, D L. A correlation coefficient for modal vector analysis[M]. Proceedings of the Proceeding of the 1st international modal analysis conference, Orlando, FL. Union College, Schenectady, NY, 1982:110-116.

[18] Vlasov V Z. Thin-Walled elastic beams(2nd edn)[M].Moscow, 1959:57-162.

Frequency Analysis of Stiffened Cylinder Based on NovelReduced Order Model

LI Yu-wei, WANG Bo, CHENG Geng-dong

(State Key Laboratory of Structural Analysis for Industrial Equipment, Department of Engineering Mechanics,Dalian University of Technology, Dalian 116024, China)

A novel reduced order model is proposed to determine the natural frequencies of the stiffened cylinder in this study, which is established by using a reduction basis along with the polynomial interpolation function and shape function of beams. The basic idea is to convert the displacements of FEM nodes to the master nodes with a few generalized DOFs. The proposed ROM can gain a significant reduction of computational cost without sacrifice of accuracy. Typical stiffened cylinder are studied and the outcomes are validated by benchmark studies. The results show that the reduced order model has the ability to identify shell lobe-type modes, which can be served as an efficient alternative model during initial design.

Polynomial interpolation function; Shape function of beams; Model reduction; Stiffened cylinder; Frequency analysis

2017-03-27；

2017-04-27基金項(xiàng)目：國(guó)家“973”計(jì)劃(2014CB049000)；國(guó)家自然科學(xué)基金(11372062, 11402049)

李玉韋(1989-)，男，在讀博士，主要研究方向：結(jié)構(gòu)動(dòng)力模型降階及優(yōu)化。E-mail:liyuwei@mail.dlut.edu.cn通訊作者：王博(1978-)，男，教授，主要從事結(jié)構(gòu)多學(xué)科優(yōu)化、航天先進(jìn)材料與結(jié)構(gòu)設(shè)計(jì)和輕質(zhì)結(jié)構(gòu)抗沖擊分析與設(shè)計(jì)。E-mail:wangbo@dlut.edu.cn

V415

A

2096-4080(2017)01-0041-08