復(fù)雜薄壁結(jié)構(gòu)模態(tài)仿真方法

段新峰 龔純 李忠華

摘要：

為節(jié)省復(fù)雜薄壁結(jié)構(gòu)模型分析的前處理時(shí)間，探究體單元在其模態(tài)分析中的可行性，以典型的四邊固支矩形薄板為對(duì)象，詳細(xì)討論Abaqus中單元類(lèi)型、積分算法和壁厚方向的網(wǎng)格層數(shù)對(duì)模態(tài)仿真結(jié)果的影響。分析結(jié)果表明：采用體單元進(jìn)行薄壁結(jié)構(gòu)模態(tài)分析時(shí)，為保證計(jì)算精度和計(jì)算效率，應(yīng)選用二次體單元C3D10或C3D20R劃分網(wǎng)格，且模型可以采用相對(duì)較粗的網(wǎng)格進(jìn)行劃分，壁厚方向僅需1層網(wǎng)格。冰箱底板體單元模態(tài)仿真試驗(yàn)對(duì)標(biāo)良好也充分驗(yàn)證這一結(jié)論。

關(guān)鍵詞：

薄壁結(jié)構(gòu)； 模態(tài)； 殼單元； 體單元； 積分算法

中圖分類(lèi)號(hào)： TB123； TB115.1

文獻(xiàn)標(biāo)志碼： B

Modal simulation method of complex thinwall structure

DUAN Xinfeng， GONG Chun， LI Zhonghua

（Midea Corporate Research Center， Foshan 528311， Guangdong， China）

Abstract：

In order to save the preprocessing time of thinwall structure model analysis and investigate the feasibility of the body element in its modal analysis， a full clamped rectangular thin plate is selected to detailedly discuss the influence of element type， integration algorithm and number of mesh layer through thickness on its modal results. The analysis results show that the second order body element C3D10 or C3D20R should be used to analyze the thinwall structure modal to ensure the precision and efficiency. The element size could be set relatively coarse and only one layer mesh through thickness. This conclusion is fully verified by the good agreement of modal simulation with test results of a refrigerator bottom panel.

Key words：

thinwall structure； modal； shell element； body element； integration algorithm

0 引 言

家電產(chǎn)品廣泛采用各種薄壁結(jié)構(gòu)。這些薄壁結(jié)構(gòu)的模態(tài)一般都比較密集，容易被激振而產(chǎn)生振動(dòng)噪聲、疲勞等問(wèn)題。[16]因此，在產(chǎn)品設(shè)計(jì)階段需要對(duì)這些薄壁結(jié)構(gòu)方案進(jìn)行模態(tài)分析和優(yōu)化，以確保其固有頻率遠(yuǎn)離激勵(lì)頻率。[7]根據(jù)有限元分析理論，這類(lèi)薄壁件適合采用殼單元進(jìn)行模擬分析[8]，但是當(dāng)薄壁結(jié)構(gòu)形狀復(fù)雜時(shí)，比如帶有大量加強(qiáng)筋的鈑金件或者注塑件，往往需要花費(fèi)較多的時(shí)間進(jìn)行中面的抽取和修補(bǔ)工作。不等壁厚的薄壁結(jié)構(gòu)使用殼單元模擬更繁瑣。直接采用體單元模擬薄壁結(jié)構(gòu)可以省去很多前期模型處理工作量和時(shí)間，但一般認(rèn)為在薄壁件厚度方向至少需要?jiǎng)澐?～4層網(wǎng)格[910]，這使得整體單元尺寸變小、模型網(wǎng)格數(shù)量劇增、計(jì)算量大幅增加、求解時(shí)間大幅延長(zhǎng)。選擇合適的體單元類(lèi)型和網(wǎng)格密度，可以在保證計(jì)算精度的前提下提高求解效率。

Abaqus具有豐富的單元庫(kù)，如C3D8R、C3D8、C3D4和C3D10等，能夠模擬各種工程問(wèn)題。[11]本文以四邊固支矩形薄板為對(duì)象，對(duì)比分析不同單元類(lèi)型、積分算法和厚度方向網(wǎng)格層數(shù)對(duì)其模態(tài)計(jì)算結(jié)果的影響，總結(jié)薄壁結(jié)構(gòu)體單元模態(tài)計(jì)算經(jīng)驗(yàn)，并通過(guò)冰箱底板的自由模態(tài)仿真試驗(yàn)對(duì)標(biāo)驗(yàn)證這一結(jié)論。

1 矩形薄板殼單元模態(tài)分析

矩形薄板長(zhǎng)邊a=150 mm，短邊b=100 mm，厚度t=5 mm，材質(zhì)為鋼材，密度ρ=7 930 kg/m3，彈性模量E=210 GPa，泊松比μ=0.3，邊界條件為四邊固支。整個(gè)模型提取中面后，在Abaqus中采用殼單元離散，并選擇Lanczos算法進(jìn)行模態(tài)分析。不同殼單元類(lèi)型、積分算法和網(wǎng)格密度對(duì)其前3階固有頻率和求解時(shí)間的影響見(jiàn)表1。網(wǎng)絡(luò)劃分編號(hào)S13和S32的前3階模態(tài)振型見(jiàn)圖1和2。對(duì)比可知，采用二次單元S8R和較粗的網(wǎng)格可以得到較精確的模態(tài)頻率和模態(tài)振型結(jié)果，并且求解效率相對(duì)較高。這主要是因?yàn)槎螁卧m然計(jì)算量大，但是其形函數(shù)為二次多項(xiàng)式，具有更高的計(jì)算精度，所以可采用更稀疏的網(wǎng)格，降低總體計(jì)算量。

2 矩形薄板體單元模態(tài)分析

針對(duì)上述固支矩形薄板模態(tài)問(wèn)題，采用不同的體單元和網(wǎng)格密度進(jìn)行求解。在對(duì)矩形薄板進(jìn)行六面體網(wǎng)格劃分時(shí)，面內(nèi)網(wǎng)格尺寸保持2 mm不變，厚度方向分別劃分1、2、3和5層網(wǎng)格。根據(jù)有限元理論，該問(wèn)題使用殼單元模擬能夠得到較精確的模態(tài)分析結(jié)果，因此以前文中網(wǎng)格劃分編號(hào)S34的計(jì)算結(jié)果為基準(zhǔn)評(píng)估體單元模態(tài)分析精度。

不同六面體單元類(lèi)型和板厚方向網(wǎng)格層數(shù)對(duì)固支薄板前3階固有頻率的影響見(jiàn)表2。對(duì)比發(fā)現(xiàn)，板厚方向僅劃分1層網(wǎng)格時(shí)，二次單元C3D20R和C3D20具有較好的固有頻率計(jì)算精度（<2%），而一次單元C3D8R和C3D8計(jì)算的固有頻率大幅偏低，尤其是一次縮減積分單元C3D8R，板厚方向劃分3層及3層以上網(wǎng)格時(shí)，一次單元計(jì)算的固有頻率才具有一定的計(jì)算精度。比較網(wǎng)格劃分編號(hào)H14和H31的固有率可知：要使前3階固有率頻達(dá)到2%的計(jì)算精度，一次單元C3D8R需要使用5層網(wǎng)格，求解耗時(shí)10.9 s；二次單元C3D20R僅需要1層網(wǎng)格，求解耗時(shí)9.6 s，計(jì)算效率更高。實(shí)際上，采用二次單元時(shí)，同等精度下薄板面內(nèi)也可采用更粗的網(wǎng)格，求解耗時(shí)會(huì)進(jìn)一步縮短。網(wǎng)格劃分編號(hào)H31的薄板前3階模態(tài)振型見(jiàn)圖3，其與殼單元振型結(jié)果一致。

實(shí)際中的工程結(jié)構(gòu)大多形狀比較復(fù)雜，很多時(shí)候難以進(jìn)行六面體單元?jiǎng)澐?，而只能直接采用四面體單元。因此，針對(duì)該固支薄板，分別采用不同四面體單元和不同網(wǎng)格密度進(jìn)行模態(tài)分析，其固有頻率見(jiàn)表3。分析對(duì)比可以發(fā)現(xiàn)：采用一次四面體單元C3D4的固有頻率計(jì)算結(jié)果偏大，尤其是板厚方向網(wǎng)格層數(shù)較少時(shí)（板厚方向1層網(wǎng)格時(shí)，固有頻率誤差達(dá)到30%左右）；二次四面體單元C3D10具有較好的固有頻率分析精度（<1%），即使板厚方向僅有1層網(wǎng)格亦然。另外，從各工況四面體單元求解時(shí)間對(duì)比來(lái)看，同等精度下，二次四面體單元C3D10具有更高的計(jì)算效率（對(duì)比網(wǎng)格劃分編號(hào)T14和T21）。網(wǎng)格劃分編號(hào)T21的薄板前3階模態(tài)振型見(jiàn)圖4，其與殼單元振型結(jié)果一致。

3 冰箱底板體單元模態(tài)分析和試驗(yàn)驗(yàn)證

根據(jù)上文分析討論可知，對(duì)于薄壁結(jié)構(gòu)，如果因?yàn)橹忻娉槿『托扪a(bǔ)工作時(shí)間過(guò)長(zhǎng)等原因而需采用體單元進(jìn)行模態(tài)分析時(shí)，為保證計(jì)算精度和提高計(jì)算效率，應(yīng)該選用二次體單元C3D10或者C3D20R，并且模型可以采用相對(duì)較粗的網(wǎng)格進(jìn)行劃分，壁厚方向僅需1層網(wǎng)格。本節(jié)選取某款冰箱底板，按照該方法進(jìn)行體單元模態(tài)仿真，并通過(guò)模態(tài)試驗(yàn)驗(yàn)證該方法仿真結(jié)果的準(zhǔn)確性。

3.1 冰箱底板體單元模態(tài)分析

某款冰箱底板結(jié)構(gòu)示意見(jiàn)圖5，板厚1.5 mm，材質(zhì)為鍍鋅鋼板，結(jié)構(gòu)比較復(fù)雜，涉及各種加強(qiáng)筋、開(kāi)孔以及翻邊。將其3D模型導(dǎo)入到Abaqus軟件中，設(shè)定網(wǎng)格尺寸為5 mm，此時(shí)板厚方向僅1層網(wǎng)格，采用四面體單元進(jìn)行自由網(wǎng)格劃分，單元類(lèi)型選用二次單元C3D10，冰箱底板有限元模型見(jiàn)圖6。材料參數(shù)通過(guò)樣條拉伸試驗(yàn)獲取，彈性模量E=199 GPa，泊松比μ=0.3，調(diào)整材料密度使模型質(zhì)量等于實(shí)際結(jié)構(gòu)質(zhì)量，得到密度ρ=7 900 kg/m3。為便于模態(tài)試驗(yàn)驗(yàn)證，直接進(jìn)行自由模態(tài)分析，底板前3階彈性模態(tài)固有頻率見(jiàn)表4，前3階模態(tài)振型見(jiàn)圖7。

3.2 冰箱底板模態(tài)試驗(yàn)

為驗(yàn)證上述模態(tài)仿真結(jié)果的準(zhǔn)確性，對(duì)冰箱底板進(jìn)行自由模態(tài)測(cè)試。用2根橡皮繩將底板吊起，使用激振器進(jìn)行激勵(lì)，采集激振器的受力信號(hào)作為參考。使用Polytec激光測(cè)振儀PSV500對(duì)預(yù)設(shè)的掃描點(diǎn)進(jìn)行測(cè)量，獲取響應(yīng)，底板模態(tài)試驗(yàn)布置見(jiàn)圖8。模態(tài)試驗(yàn)獲得的振動(dòng)位移幅值頻率響應(yīng)曲線(xiàn)見(jiàn)圖9。結(jié)合峰值頻率對(duì)應(yīng)的振型可以發(fā)現(xiàn)，頻率響應(yīng)曲線(xiàn)中最開(kāi)始的2個(gè)峰值是懸掛系統(tǒng)導(dǎo)致的，因此該底板結(jié)構(gòu)的前3階彈性模態(tài)頻率（見(jiàn)表4）分別為26.88、122.50和162.97 Hz，與其固有頻率仿真結(jié)果誤差在4%以?xún)?nèi)。

該底板結(jié)構(gòu)前3階模態(tài)振型試驗(yàn)結(jié)果見(jiàn)圖10。圖10與圖7對(duì)比可以發(fā)現(xiàn)，兩者趨勢(shì)基本一致、吻合較好。綜合來(lái)看，該模態(tài)試驗(yàn)結(jié)果可驗(yàn)證二次體單元模態(tài)仿真結(jié)果的準(zhǔn)確性。

4 結(jié) 論

以矩形薄板為對(duì)象，詳細(xì)討論單元類(lèi)型、積分算法和網(wǎng)格密度對(duì)模態(tài)仿真結(jié)果的影響，并對(duì)冰箱底板進(jìn)行模態(tài)仿真試驗(yàn)對(duì)標(biāo)。綜合計(jì)算精度和計(jì)算效率來(lái)看，針對(duì)薄壁結(jié)構(gòu)模態(tài)分析問(wèn)題，可以給出以下幾點(diǎn)建議：

（1）對(duì)于比較容易進(jìn)行中面提取的薄壁結(jié)構(gòu)，可在提取中面后，劃分相對(duì)較粗的網(wǎng)格并選用二次殼單元S8R進(jìn)行模態(tài)分析。

（2）對(duì)于中面抽取和修補(bǔ)工作量較大的復(fù)雜薄壁結(jié)構(gòu)或者變厚度的薄壁結(jié)構(gòu)，為節(jié)省模型前處理時(shí)間，可考慮直接采用體單元進(jìn)行模態(tài)分析，此時(shí)，為保證計(jì)算精度和計(jì)算效率，應(yīng)該選用二次體單元C3D10或者C3D20R，模型采用相對(duì)較粗的網(wǎng)格進(jìn)行劃分，壁厚方向僅需1層網(wǎng)格。

參考文獻(xiàn)：

[1] 郭利， 黃道業(yè). 基于ANSYS的滾筒洗衣機(jī)振動(dòng)優(yōu)化設(shè)計(jì)[J]. 重慶理工大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2012， 26（10）： 3437. DOI： 10.3969/j.issn.16748425（z）.2012.10.009.

[2] 陳敏， 石衛(wèi)衛(wèi)， 張華偉， 等. 基于形貌優(yōu)化的滾筒洗衣機(jī)箱體的結(jié)構(gòu)設(shè)計(jì)[J]. 機(jī)械， 2016， 43（10）： 4042. DOI： 10.3969/j.issn.10060316.2016.10.009.

[3] 王勇， 田麗娟， 王珂. 基于模態(tài)仿真分析的空調(diào)系統(tǒng)的低頻噪聲優(yōu)化[J]. 家電科技， 2017（7）： 6871. DOI： 10.3969/j.issn.16720172.2017.07.019.

[4] 盧劍偉， 李曉陽(yáng)， 祖玉建. 分體空調(diào)器室外機(jī)鈑金件動(dòng)態(tài)特性分析及優(yōu)化設(shè)計(jì)[J]. 家電科技， 2015（2）： 5254. DOI： 10.3969/j.issn.16720172.2015.02.018.

[5] 白連社， 張波， 羅健. 基于模態(tài)仿真的冰箱減振降噪技術(shù)研究[J]. 電器， 2012（S1）： 127129. DOI： 10.3969/j.issn.16728823.2012.z1.026.

[6] 郭維， 劉斌， 馮濤， 等. 冰箱壓縮機(jī)機(jī)殼實(shí)驗(yàn)?zāi)B(tài)分析[J]. 噪聲與振動(dòng)控制， 2010（3）： 6770. DOI： 10.3969/j.issn.10061355.2010.03.018.

[7] 朱石堅(jiān)， 樓京俊， 何其偉， 等. 振動(dòng)理論與隔振技術(shù)[M]. 北京： 國(guó)防工業(yè)出版社， 2006： 147159.

[8] 胡于進(jìn)， 王璋奇. 有限元分析及應(yīng)用[M]. 北京： 清華大學(xué)出版社， 2009： 143164.

[9] 史冬巖， 莊重， 高山， 等. 基于Abaqus的模態(tài)分析方法對(duì)比及驗(yàn)證[J]. 計(jì)算機(jī)輔助工程， 2013， 22（S2）： 432435. DOI： 10.3969/j.issn.10060871.2013.z2.105.

[10] 曹金鳳， 石亦平. Abaqus有限元分析常見(jiàn)問(wèn)題解答[M]. 北京： 機(jī)械工業(yè)出版社， 2009： 134135.

[11] 莊茁， 張帆， 岑松， 等. Abaqus非線(xiàn)性有限元分析與實(shí)例[M]. 北京： 科學(xué)出版社， 2005： 110.

（編輯 武曉英）