層合板頻率響應與聲功率雙目標優(yōu)化設計

吳錦武，趙飛，王縣委，李根

?

層合板頻率響應與聲功率雙目標優(yōu)化設計

吳錦武1，趙飛2，王縣委1，李根1

(1. 南昌航空大學飛行器工程學院，江西南昌 330063；2. 中航航空服務保障(天津)有限公司，天津 300301)

利用遺傳算法對復合材料層合板結構的固有頻率間隔和輻射聲功率進行雙目標優(yōu)化設計。利用分層理論結合有限元模型求解層合板的固有頻率和振速分布。通過聲輻射模態(tài)理論，計算層合板結構輻射聲功率。以鋪設角度作為設計變量，第一階與第二階固有頻率間隔和輻射聲功率作為雙目標優(yōu)化目標函數，以某4層的層合板結構為例，采用目標加權法優(yōu)化目標函數。研究了不同權重系數、不同頻率時固有頻率間隔最大化和聲功率最小化對應的優(yōu)化鋪設角度。由數值分析結果可知：不同的權重系數比下獲得的Pareto最優(yōu)解不同；在同一權重系數下，兩個優(yōu)化目標所起的作用不盡相同；隨著頻率的增加，選擇相對較大的權重系數可使Pareto最優(yōu)解較好地兼顧兩個優(yōu)化目標。

鋪設角度；固有頻率；聲功率；遺傳算法

0 引言

由于復合材料結構具有比強度大、比剛度高、材料性能可設計等特點，近年來針對復合材料層合板結構減振降噪的優(yōu)化設計，越來越引起人們的關注[1-8]。隨著纖維板聲學特性研究的深入，纖維鋪設角度是層合板結構中一個重要的結構參數，以其作為設計變量的結構優(yōu)化設計也越來越多[9]。如在減振優(yōu)化設計方面，Topal等[1]通過一階剪切理論計算層合板固有頻率，并以固有頻率最大化為目標函數，以鋪層順序為設計變量，將改進的可行方向法(Modified feasible direction method)和有限元相結合，對層合板的鋪層順序進行了優(yōu)化設計。其次，當外激勵頻率越高時，可將增大相鄰兩階固有頻率之間的間隔作為優(yōu)化目標，如Adali等[2]以最大基頻和最高頻率間隔作為設計變量在自由振動下確定了對稱層合板的最佳鋪層順序。另外Kayikci等[3]以固有頻率及固有頻率間隔作為目標函數，以鋪設角度作為設計變量進行層合板結構優(yōu)化分析。

在聲優(yōu)化設計方面，以層合板在外激勵作用下的聲輻射功率最小作為設計變量，優(yōu)化復合材料鋪層順序[4-6]。其中陳爐云等[6]利用遺傳算法以聲輻射功率最小化為目標函數，對層合板鋪設層數、厚度以及鋪層順序等進行集成優(yōu)化設計。以上研究一般基于等效單層板理論分析其動力學特性。

上述研究大多數針對單一目標進行優(yōu)化設計。針對實際需要，層合板結構多目標優(yōu)化設計的研究越來越多[7-8]。但是大多數目標函數集中在結構強度、剛度、質量、重量等。本文擬通過分層理論結合有限元模型得到層合板結構的固有頻率。基于聲輻射模態(tài)得到結構輻射聲功率，通過遺傳算法以結構固有頻率間隔與聲功率同時作為優(yōu)化目標函數，以鋪設角度作為設計變量，對層合板結構的減振與降噪進行優(yōu)化設計研究。

1 結構動力響應

基于分層理論，在沿板厚方向采用分層插值，層合板結構物理模型如圖1所示。

層合板結構的位移場可表示為[10]

(2)

式(3)中：表示層合板的第個鋪設層的厚度；為板厚的局部坐標。根據有限元理論，式(1)中，、、可表示為:

(4)

其中:

(6)

其中矩陣元素見參考文獻[11]。

根據有限元理論，層合板結構的動力學方程為

其中：為單元質量矩陣;為單元剛度矩陣;為阻尼矩陣；為單元等效節(jié)點力矩陣。由式(8)可求得結構表面法向位移、法向速度向量()。層合板的各階固有頻率可由行列式求得。

2 固有頻率間隔優(yōu)化

研究結構的固有頻率是避免結構共振的一個重要內容。由參考文獻[3]可知，為了避免結構共振，可定義固有頻率間隔作為優(yōu)化目標函數。頻率參數[3]為:

(10)

針對層合板的固有頻率，取最大相鄰階固有頻率間隔作為優(yōu)化變量。結構-頻率優(yōu)化模型[11]為

3 結構頻率間隔-聲輻射雙目標優(yōu)化模型

本文主要采用基于目標加權的遺傳算法來優(yōu)化目標函數：其核心思想是對目標中每一個目標向量分配一個權重，將多目標分量乘以各自的權重系數后再相加，構成一個新的目標函數，將其轉化成一個單目標優(yōu)化問題求解。在對每個個體取不同的權重函數的基礎上，利用遺傳算法求解得到多目標優(yōu)化問題的多個Pareto最優(yōu)解，即為多目標規(guī)劃問題的有效解。

式(12)中：1表示層合板的第一階固有頻率(基頻)，通過取合適的1、2，使1+2=1，這樣可把多目標優(yōu)化問題轉化成為單目標優(yōu)化問題。權重系數=[1,2,…,]可反映每個目標函數的重要性。通過選取不同的權重系數，可以獲得不同的Pareto最優(yōu)解。

4 數值分析

以某四層碳纖維/環(huán)氧復合材料層合板為例，層合板性能參數為：長=0.3480 m，寬=0.3048 m，厚度=0.002 m，楊氏模量，，剪切模量，，阻尼比=0.05；固定激勵頻率=200 Hz和=500 Hz，板中心激勵力大小為1，初始角度為[0o]4，利用遺傳算法對復合材料層合板鋪層順序進行單目標和雙目標優(yōu)化設計。

4.1 固有頻率間隔優(yōu)化

鋪設角度采用離散形式，相鄰兩角度之間的步長分別為90o、45o、30o、15o，初始角度為[0o]4即為[0o/0o/0o/0o]，取前兩階無量綱固有頻率間隔為優(yōu)化變量，優(yōu)化層合板鋪層順序。表1為四層層合板的優(yōu)化結果。

表1 不同角度步長時層合板的鋪設角度優(yōu)化結果

由表1可見，隨著相鄰兩角度之間的步長減小，即角度遞增量的增加，目標函數逐漸增大，從而說明角度遞增量越小，優(yōu)化效果越好；但隨著角度遞增量的減小，計算量成倍的增加，尋找全局最優(yōu)解將變得越來越困難，因此在實際應用過程中目標函數的角度遞增量要根據實際情況靈活選取。

4.2 聲功率最小化優(yōu)化

仍以上述結構為例，以結構輻射聲功率最小作為目標，采用遺傳算法優(yōu)化層合板鋪層順序，初始角度分別為[0o]4、[0o]8。表2為激勵頻率和時層合板鋪層順序優(yōu)化結果。

表2 激勵頻率wp=200 Hz和wp=500 Hz層合板鋪層順序優(yōu)化

由表2可見，通過對層合板鋪層順序的優(yōu)化，其聲輻射功率級在激勵頻率為200 Hz時降低了1.9 dB，激勵頻率為500 Hz時降低了6.2 dB。說明通過鋪設角度優(yōu)化可使層合板結構具有良好的降噪效果。這與復合材料的各向異性有關。

4.3 雙目標優(yōu)化

對于固有頻率間隔與聲功率雙目標優(yōu)化，主要的是選擇合適的權重系數1、2。本文選取如表3和表4中的12個權重系數比計算上述層合板的鋪層順序雙目標優(yōu)化的Pareto最優(yōu)解，由于實際層合板結構45o鋪設角度偏多，同時為了減少計算工作量，故層合板的鋪層角度選擇區(qū)間為[-45o，0o，45o，90o]四個離散值。

表3 不同權重系數比的雙目標優(yōu)化Pareto最優(yōu)解(200 Hz)

由表3、4可見：不同的權重系數比下獲得的Pareto最優(yōu)解不同，因此本文分別選取權重系數比1:2分別為1:1、1:10，1:20做詳細討論。

從表3和表4可知：在相同頻率下，權重系數從小于1:10開始，隨著權重系數的變小，優(yōu)化后的目標函數和()數值隨之大幅度變小。當權重系數在10:1到1:5之間，優(yōu)化后的目標函數和()結果相同，相應的優(yōu)化設計變量（即鋪設角度）也一致。當權重系數大于1:10時，優(yōu)化后的目標函數數值相差較大。

表4 不同權重系數比的雙目標優(yōu)化Pareto最優(yōu)解(500 Hz)

表5~7分別表示權重系數取1:2=1:1、1:10，1: 20時，激勵頻率=200 Hz時層合板鋪層順序多目標優(yōu)化結果。

表5 激勵頻率wp=200 Hz多目標鋪設角度優(yōu)化(l1:l2=1:1)

表6 激勵頻率wp=200 Hz多目標鋪設角度優(yōu)化(l1:l2=1:10)

表7 激勵頻率wp=200 Hz多目標鋪設角度優(yōu)化(l1:l2=1:20)

從表5~7可知：隨著1:2的增大，相同角度步長時的聲功率數值也逐漸減小。

由表5～10可知，當權重系數比1:2=1:1時，兩個優(yōu)化目標中所起的作用明顯遠遠大于聲輻射功率，以至于最終優(yōu)化結果與頻率間隔單目標優(yōu)化時的結果基本一致，導致聲輻射功率過大，激勵頻率為500 Hz時尤其明顯。當權重系數比增大到1:2=1:10甚至是1:20時，兩個優(yōu)化目標對最終結果的貢獻趨近于相等。因此在200 Hz外激勵下，選擇權重系數比1:2=1:20，而在500 Hz外激勵下，權重系數比1:2選擇1:10比較合適。

表8 激勵頻率wp=500 Hz多目標鋪設角度優(yōu)化(l1:l2=1:1)

表9 激勵頻率wp=500 Hz多目標鋪設角度優(yōu)化(l1:l2=1:10)

表10 激勵頻率wp=500 Hz多目標鋪設角度優(yōu)化(l1:l2=1:20)

層合板在200 Hz外激勵下，權重系數比1:2=1:20時的遺傳算法迭代過程如圖2所示。在500 Hz外激勵下，權重系數比1:2選擇1:10時的遺傳算法迭代過程如圖3所示。從圖2、3可知：

在角度步長一致的情況下，外激勵頻率越高，同一參數優(yōu)化計算需要的迭代次數越多。另外，角度步長越小，種群均值收斂速度越慢，其中角度步長為90°時，種群均值收斂速度最快。

(a) 90°

(b) 45°

(c) 30°

(d) 15°

圖2. 不同遞增角度的層合板多目標優(yōu)化迭代過程(=200 Hz、1:2=1:20)

Fig.2 Iteration process of double objective optimization(=200 Hz、1:2=1:20)

(a) 90°

(b) 45°

(c) 30°

(d) 15°

圖3 不同角度步長的層合板多目標遺傳算法迭代過程(=500 Hz、1:2=1:10)

Fig.3 Iteration process of double objective optimization(=500 Hz、1:2=1:10)

在200 Hz外激勵下，權重系數比1:2=1:20時的Pareto最優(yōu)解較好地兼顧了兩個優(yōu)化目標，即對應優(yōu)化角度分別為：[90o/90o/0o/0o]、[90o/0o/-45o/0o]、[90o/0o/-60o/0o]以及[90o/15o/-45o/0o]。在以上四種優(yōu)化角度下，層合板結構在全頻段內的聲功率如圖4所示。為了對比，圖4中增加了[0o]4鋪設角度下對應的聲功率。

由圖4可知：四種鋪設角度的對應的聲功率最大值基本相同。另外對照表5所示內容，由圖4也可知，[90o/15o/-45o/0o]鋪設角度下的基頻最小為82 Hz，同時頻率間隔最大。而[90o/90o/0o/0o]鋪設角度下的基頻最大為100.6 Hz，同時頻率間隔最小。

在500 Hz外激勵下，權重系數比1:2=1:10時對應的優(yōu)化角度分別為[0o]4、[90o/0o/90o/0o]以及[-30o/90o/90o/-30o]。根據上述三種優(yōu)化鋪設角度，計算可得到層合板結構前10階固有頻率，如表11所示。

表11 激勵頻率wp=500 Hz層合板前10階固有頻率(l1:l2=1:10)

另外，500 Hz三種優(yōu)化角度下結構在全頻段內的聲功率如圖5所示。

由圖5可知：三種鋪設角度的對應的聲功率最大值基本相同。[0o]4鋪設角度下的基頻最小為82 Hz，而且頻率間隔最小。[90o/90o/0o/0o]鋪設角度下的基頻最大，為100.6 Hz，頻率間隔也最大。同時當頻率超過700 Hz時，隨著頻率的增高，三種鋪設角度下的層合板聲功率基本一致。

綜上所示，選擇適當鋪設角度以及權重比，可使層合板結構在減振和降噪方面同時得到優(yōu)化。

5 結論

(1) 通過分層有限元模型可計算得到層合板結構的固有頻率。

(2) 通過遺傳算法，選擇合適的鋪設角度可進行復合材料層合板結構聲功率的最小化設計。

(3) 選擇合適的權重系數，可實現層合板結構的減振和降噪雙目標優(yōu)先，改善層合板的聲振特性，為實際復合材料的設計與應用提供參考。

[1] TOPAL U, UZMAN ü. Frequency optimization of laminated skew plates [J]. Materials and Design, 2009, 30(8): 3180-3185.

[2] ADALI S, VERIJENKO V E. Optimum stacking sequence design of symmetric hybrid laminates undergoing free vibrations[J]. Composite Structures, 2001, 54(3): 131-138.

[3] KAYIKCI R, SONMEZ F O. Design of composite laminates for optimum frequency response[J]. Journal of Sound and Vibration, 2012, 331(8): 1759-1776.

[4] DU Jianbin, OLHOFF. N. Minimization of sound radiation from vibrating bi-material structures using topology optimization. Structural and Multidisciplinary Optimization[J]. Structural and Multidisciplinary Optimization, 2007, 33(5): 305-321.

[5] Duhring M B, Jensen J S, Sigmund O. Acoustic design by topology optimization[J]. Journal of Sound and Vibration, 2008, 317(5): 557-575.

[6] 陳爐云, 張裕芳. 基于遺傳算法的復合材料結構-聲輻射優(yōu)化研究[J]. 復合材料學報, 2012, 29(3): 203-207.

CHEN Luyun, ZHANG Yufang. Composite material structural-acoustic optimization on genetic algorithm[J]. Acta Material Compositae Sinica. 2012, 29(3): 203-207.

[7] PELLETIER J L, VEL. S S, Multi-objective optimization of fiber reinforced Composite laminates for strength, stiffness and minimal mass [J]. Computers & Structures, 2006, 84(29-30): 2065-2080.

[8] RAHUL Sandeep G, CHAKRABORTY D. Multi-objective optimization of hybrid laminates subjected to transverse impact[J]. Composite Structures, 2006, 73(3): 360-369.

[9] 王佐民, 俞悟周. 鋁纖維板共振結構聲學特性的理論分析[J]. 聲學技術, 2005, 24(3): 183-185.

WANG Zuomin, YU Wuzhou. Theoretical annlysis on acoustical characteristics aluminum fiberboard[J]. Technical Acoustics, 2005, 24(3): 183-185.

[10] REDDY J N. Mechanics of laminated composite plates and shells [M]. New York: CRS Press, 2004: 96-97.

[11] 趙飛. 基于分層理論的層合板聲輻射分析與結構優(yōu)化[D]. 江西: 南昌航空大學, 2013: 11-39.

ZHAO Fei. Sound radiation analysis and structure optimization of laminated plates based on layerwise theory [D]. Jiangxi: Nanchang Hangkong University, 2013:11-39

Double objective optimization of frequency response and acoustic power for laminated composite plates

WU Jin-wu1, ZHAO Fei2, WANG Xian-wei1, LI Gen1

(1.School of Aircraft Engineering, Nanchang Hang Kong University, Nanchang330063,Jiangxi,China;2. Aviation Service Assurance Ltd, Aviation Industry Corporation of China, Tianjin300301, China)

Both the separation between two adjacent natural frequencies and the sound power of laminated composite plate are optimized by use of genetic algorithm in this paper. The layer-wise finite element model is imposed to determine the natural frequencies and the velocity distributions of laminated composite plates. Based on the acoustic radiation model, sound power is calculated. The fiber laying angles are chosen as design variables and the design objective is to minimize the sound power and to maximize the separation of two adjacent natural frequencies. A four-layer laminated plate are used as examples, the objective weighting method is used to optimize the objective function. And the optimal laying angles with sound power minimization and frequency separation maximization are discussed under different frequencies and different weight coefficients. The numerical simulations show that there are different Pareto optimal solutions under different weight coefficients. With frequency increasing, a relatively large weighting coefficient can be chosen for the Pareto optimal solution to balance the two optimization objectives.

the laying angle; natural frequency; sound power; genetic algorithm

TB332 TH113.1

A

1000-3630(2016)-02-0155-07

10.16300/j.cnki.1000-3630.2016.02.0013

2015-04-03;

2015-07-14

國家自然科學基金資助項目(51265038)

吳錦武(1976－), 男, 江西高安人, 博士, 教授, 研究方向為機械振動與聲控制。

吳錦武, E-mail: jinwu.w@hotmail.com