國 立 交 通 大 學
機 械 工 程 學 系
碩 士 論 文
小 功 率 平 面 揚 聲 器 之 設 計
與 製 作 研 究
Design and Fabrication of Low Power Flat-Panel Loudspeaker
研 究 生 : 黃 永 立
指 導 教 授 : 金 大 仁
小功率平面揚聲器之設計與製作研究
Design and Fabrication of Low Power Flat-Panel Loudspeaker
研 究 生:黃永立 Student : Yung-Li Huang
指導教授:金大仁 教授 Advisor : Dr. Tai-Yan Kam
國 立 交 通 大 學 機 械 工 程 學 系
碩士論文
A Thesis
Submitted to Institute of Mechanical Engineering College of Engineering
National Chiao Tung University in partial Fulfillment of the Requirements
for the Degree of Master of Science
in
Mechanical Engineering July 2007
Hsinchu, Taiwan, Republic of China
小功率平面揚聲器之設計與製作研究
研究生:黃永立 指導教授:金大仁 博士 國立交通大學機械工程學系 中文摘要 本文研究低功率平面揚聲器之設計與實作。設計方面,運用有限元素分 析軟體 ANSYS 分析揚聲板的振動行為,並設計揚聲板的結構以減少中音谷 的不良影響。實作方面,研究應用複合材料製作揚聲器的彈性支承,懸邊 的製作與揚聲器的組裝方式等。本研究製作數種揚聲器,藉測量其實驗聲 壓曲線來驗證理論聲壓分析模式的準確性,然後應用驗證後的模式研究振 動板幾何比與音圈參數對聲壓曲線的影響,並找出最佳的振動板幾何比與 音圈尺寸以提升聲壓曲線的平滑度與感度。最後本研究應用最佳設計結果 製作揚聲器並量測其聲壓曲線,結果顯示最佳設計確實能產出聲壓較高及 中音谷較小的揚聲器。Design and Fabrication of Low Power Flat Panel Loudspeaker
Student : Yung-Li Huang Advisor : Dr. Tai-Yan Kam
Institute of Mechanical Engineering National Chiao Tung University
ABSTRACT
Two parts which include design and fabrication of low power flat-panel loudspeakers are studied in the thesis. The first part is about harmonic analysis and design of radiating panels using the finite element software ANSYS. The second part is about manufacture of loudspeaker components which includes suspension and elastic support, and the assembling of the loudspeakers. In fabricating the speaker, the soft suspension and elastic support are used to design the initial resonance of the speaker so that the effective frequency range can be increased and better SPL value and loudspeaker stability can be obtained in lower frequency range. Several loudspeakers are fabricated and their experimental SPL curves are used to verify the finite element model. Then based on the finite element model, the effects of aspect ratio of panel and the parameters such as voice coil diameter on the SPL curves of the loudspeakers are studied. The optimal panel aspect ratio and the diameter of voice coil for the loudspeakers are determined. It has been shown that the optimally designed loudspeaker can have smoother SPL curve and higher sensitivity.
誌謝
在這短短的兩年碩士求學過程,充滿許多難忘的回憶,對我的人生影 響重大,要感謝我的老師、家人、女友彥汶、同學及朋友對我求學的全力 支持,讓我可以無憂無慮地完成我的學業。 感謝指導教授金大仁博士的耐心指導,讓我能在完成學業及論文外, 學習到正確的求學態度及做人處世的道理。感謝偉芬學姐、學長志明、清 榮、昌毅、于昇、崧任、加融、欣翰、慶博、建勳、哲偉和宗鴻的照顧, 提供我許多學業上及生活上的協助,並在我遇到困難時,提供寶貴的意見。 感謝同窗好友國真、魁原、聖傑和在這兩年來的互相勉勵,讓我在學 習過程中不會感到孤單,感謝學弟志傑、宏銘、政霖、建儒讓我的生活充 滿歡樂,僅以此文獻給所有關心我的人 永立 2007.7 于交大目 錄
中 文 摘 要 ……… i 英 文 摘 要 ……… ii 致 謝 ……… iii 目 錄 ……… iv 表 目 錄 ……… vi 圖 目 錄 ……… vii 第一章 緒 論………. 1 1.1 前言……….. 1 1.2 文 獻 回 顧 … … … . . 2 1.3 研究方法… … … . 3 第二章 揚聲板的振動分析………... 5 2.1 一階剪變形平板理論………. 5 2.1.1 基本假設………... 5 2.1.2 應力與應變關係………... 6 2.1.2 構成方程式(Governing equation)……….. 8 2.2 多層一階剪變形平板理論………. 9 2.2.1 位移與應變………... 10 2.2.2 應力與應變………... 12 2.3 複合材料三明治板的彈性支承……….. 15 2.4 應變能與動能………. 16 2.5 瑞雷-黎次法(Rayleigh-Ritz method)……… 17 2.6 結構系統之阻尼………. 19 2.7 聲壓值計算………. 20 2.7.1 聲壓公式推導………... 20 第三章 ANSYS 有限元素模型之建構與分析………...… 24 3.1 模型建構………. 24 3.1.1 ANSYS 模型建立步驟………... 243.2 ANSYS 模擬分析中各參數的取得……… 26 3.2.1 系統的彈性係數………... 27 3.2.2 激振力的給定………... 27 3.2.3 阻尼比的給定………... 27 3.3 ANSYS 簡諧激振分析與聲壓公式之應用探討……… 28 第四章 平面揚聲器製作與實驗程序………... 30 4.1 平面揚聲器製作………. 30 4.1.1 複合材料疊層板之製作………... 30 4.1.2 三明治振動板之製作………... 30 4.1.3 音圈與場磁鐵之製作………... 31 4.1.4 彈性懸邊之製作………... 31 4.1.5 碳纖支承之製作………... 32 4.1.6 揚聲器之組裝………... 33 4.2 頻率響應實驗………. 34 4.3 聲壓相關實驗………. 34 第五章 小型平面揚聲器之設計與討論... 36 5.1 彈性懸邊之設計………. 36 5.2 玻纖彈波與碳纖支承之設計………. 37 5.3 系統二揚聲器改變設計參數造成的影響………. 38 5.3.1 改變音圈參數造成的影響... 38 5.3.2 改變振動板長寬比造成的影響... 39 5.3.3 模擬分析討論與實驗驗證... 40 5.4 系統一揚聲器改變音圈參數造成的影響………. 43 第六章 結論與未來研究方向... 44 6.1 結論………. 44 6.2 未來研究方向………. 45 參 考 文 獻 ……… 46
表 目 錄
表 3.1 揚聲器之阻尼常數………...………... 49 表 3.2 材料常數表………...………... 49 表 4.1 音圈尺寸表………...………... 50 表 4.2 磁鐵尺寸表………...………... 50 表 4.3 系統一揚聲器參數………...………... 51 表 4.4 系統二揚聲器參數………...………... 51 表 4.5 系統三揚聲器參數………...………... 52 表 5.1 絹布材料常數………...………... 53 表 5.2 各長寬比尺寸表………...………... 53 表 5.3 長寬比2.3音谷………...………... 53 表 5.4 長寬比1.7音谷………...………... 54 表 5.5 長寬比1.28音谷………...………... 54 表 5.6 各長寬比節線位置………...………... 54 表 5.7 長寬比1.28音圈與節線比例………...……….. 55 表 5.8 長寬比1.7音圈與節線比例………. 55 表 5.9 長寬比2.3音圈與節線比例………...………… 55
圖 目 錄
圖 1.1 平面揚聲器之揚聲………... 56 圖 1.2 振動板之振形……….. 56 圖 2.1 多層一階剪變形位移場示意圖(三層)………... 57 圖 2.2 複合材料積層板之主軸與材料座標系統………... 57 圖 2.3 複合材料積層板沿厚度方向之合力與合力矩…………... 58 圖 2.4 複合材料三明治板之邊界條件………... 58 圖 3.1 系統一揚聲器 ANSYS 模型………….……… 59 圖 3.2 系統二揚聲器 ANSYS 模型……… 59 圖 3.3 系統一揚聲器實驗與模擬分析比較……… 60 圖 3.4 系統二揚聲器實驗與模擬分析比較……… 60 圖 4.1 揚聲器主要結構…………..………. 61 圖 4.2a 系統一揚聲器………..………. 61 圖 4.2b 系統二揚聲器………..………. 62 圖 4.2c 系統三揚聲器………..………. 62 圖 4.3 三明治板製作及輔助材料疊層順序……..………. 63 圖 4.4 熱壓機………... 63 圖 4.5 複合材料積層板之加熱加壓硬化成型製程圖………..……. 64 圖 4.6 複合材料三明治板完成圖…………..……… 64 圖 4.7 音圈的結構…………..……….. 65 圖 4.8 場磁鐵的結構…………..……….. 65 圖 4.9 懸邊模型…………..……….. 65 圖 4.10 矽膠彈性懸邊…………..………... 66 圖 4.11 鋁製懸邊模型…………..………... 66 圖 4.12 絹布懸邊…………..………... 67 圖 4.13 碳纖夾真空帶…………..………... 67 圖 4.14 z 字型碳纖模具………..……….... 68 圖 4.15 z 字型碳纖片…………..……….... 68圖 4.16 切割 z 字型刀具………..………... 69 圖 4.17 治具一…………..………... 69 圖 4.18 黏製彈性支承………. 70 圖 4.19 系統一半成品………. 70 圖 4.20 懸邊、振動板黏合……….. 71 圖 4.21 半成品與振動板黏合……….. 71 圖 4.22 系統一揚聲器………. 72 圖 4.23 振動板治具………. 72 圖 4.24 音圈定位………. 73 圖 4.25 懸邊定位器………. 73 圖 4.26 懸邊、振動板、框黏合……….. 74 圖 4.27 系統二揚聲器………. 74 圖 4.28 B&K Pulse 頻譜分析儀……… 75 圖 4.29 LMS 聲壓頻譜儀………. 75 圖 4.30 聲壓實驗架設………. 76 圖 5.1 傳統彈波………….……….….. 77 圖 5.2 玻纖彈波………….……….….. 77 圖 5.3 系統二、系統三實驗聲壓曲線圖………….……….. 78 圖 5.4 音圈強度對系統二揚聲器影響………….……….. 78 圖 5.5 音圈重量對系統二揚聲器影響………….……… 79 圖 5.6 長寬比 2.3 聲壓分析………….………. 79 圖 5.7 長寬比 1.7 聲壓分析………….………. 80 圖 5.8 長寬比 1.28 聲壓分析………….……… 80 圖 5.9 長寬比與節線距離關係………….………. 81 圖 5.10 長寬比 1.28 音谷落差………….……… 81 圖 5.11 長寬比 1.7 音谷落差………….……….. 82 圖 5.12 長寬比 2.3 音谷落差………….……….. 82 圖 5.13 長寬比、振動板節線與寬度關係………….……… 83
圖 5.14 長寬比 1.28 系統四揚聲器………….……… 83 圖 5.15 長寬比 1.28 分析與實驗聲壓圖………….………. 84 圖 5.16 系統一、系統四、系統五揚聲器分析曲線比較……… 84 圖 5.17 系統一、系統四、系統五揚聲器相同振形比較……… 85 圖 5.18 系統一、系統四、系統五揚聲器相同振形比較……… 85 圖 5.19 系統一、系統四、系統五揚聲器相同振形比較……… 86 圖 5.20 系統一、系統四、系統五揚聲器相同振形比較……… 87 圖 5.21 系統一、系統四、系統五揚聲器實驗聲壓曲線……… 88 圖 5.22 系統一揚聲器改變音圈重量聲壓比較圖……… 88
第一章 緒論 1.1 前言 複合材料是利用數種不同性質的材料,以協同的加工方法將材料結合 在一起,擷取各材料的特性和優點,使複合材料能發揮優異性能並滿足實 際需求。其中以纖維強化的複合材料應用較廣,此種材料在纖維方向具有 高強度、高勁度之特性,可透過不同疊層方式及不同疊層角度之設計來達 到實際的需求,更由於具有質量輕、強度高、設計之多變化性等有別於傳 統金屬材料的性質,使得目前複合材料已被廣泛的應用在航太、造船、飛 彈、汽車、休閒運動器材等受重量限制及高強度需求的結構上。 平面揚聲器主要是由激振器推動振動板,藉由振動板變形而推動空氣 來產生聲音(如圖 1.1)。激振器推動振動板,使振動板表面元素之空氣產生 速度,進而產生聲壓,由於振動板之變形有凹有凸(如圖 1.2),所以產生之 空氣速度亦有正有負,因此聲壓才有高低起伏的現象。 當激振頻率達到振動板的第一個自然頻率 時,揚聲器會因為振動板 的大變形而開始產生較高的聲壓值,而且揚聲器在此頻率之後才會有較良 好之聲壓值,所以振動板的 值可視為揚聲器有效頻寬的起點,故如何降 低 值卻又不影響揚聲器整體的穩定度與聲壓值表現,對於揚聲器而言是 相當重要的研究方向。 0 f 0 f 0 f 本文即是探討由柔軟的彈性懸邊作為振動板的邊界,降低振動板的第 一個自然頻率,以增加揚聲器之有效頻寬,讓揚聲器在低頻率也能有良好 的聲壓值與穩定度。本文也嘗試改變音圈的尺寸、重量與振動板的幾何比 例,提升聲壓曲線的平滑度與平均感度,以研製出理想之小型平面揚聲器。
1.2 文獻回顧
在分析平板的理論方面,從古典板理論(Classical plate theory,簡稱 CPT)[1]改進而成的古典積層板理論(Classical lamination theory),對於複合 材料薄板的力學分析已能得到不錯的結果,但只適用於長厚比大於 80 之平 板;對於厚板而言,由於複合材料積層板之側向剪力模數(Shear modulus) 比沿纖維方向的楊氏係數(Young’s modulus)低很多,且在厚度上較薄板高出 許多,因此容易產生側向剪變形,所以古典板理論不適合分析較厚之複合 材 料 板 。 為 此 , Mindlin 提 出 了 一 階 剪 變 形 理 論 (The First-order shear deformation theory,簡稱 FSDT)[2],首先將側向剪力的影響加以考慮,但是 因為假設側向剪力分布為常數,並不符合實際的情況,於是 Whitney[3、4] 便提出了剪力修正因子來加以修正,此種理論比較適合用在長厚比大於 15 的結構上;之後,學者又提出了各種高階剪變形的理論,雖有提高理論值 與實際狀況相比之準確性,但往往較適用於厚板結構(長厚比大於 15)之情 況,且其計算上比較複雜許多,而本文中使用之複合材料結構板並不在厚 板結構的範圍內,所以仍以一階剪變形理論為主,來分析振動板之變形行 為。
在研究三明治板的文獻方面,Reissner[5]推導 Governing equation 應用 在小變形、等向性的三明治板,文獻假設面層像薄膜,而且忽略了中心層 平行面層的應力,從此,許多文獻也在三明治板的理論慢慢的變化。後來 Liaw and Little[6]根據 Reissner 理論解出了多層三明治結構彎曲的問題; Azar[7]延伸 Liaw and Little 的結果來討論非等向性面層;O’Connor[8]提出 用有限元素來分析三明治結構,他用平面彈性元素來構建中心層,樑元素 來構建面層;Kanematsu[9]用 Ritz method 來分析矩型板的彎曲和振動, Bardell 與 Dunsdon [15]用一組新的三角函數以提供 Element p-enrichment,
如此可以讓特徵值有更快的收斂;Ding [16]利用靜態樑兩端受到彈性支承 時的形狀函數,擴充應用到板上以求出一形狀函數,再利用 Rayleigh-Rize method 求得矩形板的每個模態;Hanna 與 Leissa [17]用 Ritz method 求解 Higher order shear deformation 的薄板的自然頻率,當取多項式到 108 項與
144 項時, =0.1 h b 與 =0.5 h b 的薄板前五項頻率收斂非常快。 對於聲學與聲壓計算方面,Morse[10]中推導出了聲源在空氣中傳遞之 聲壓方程式,在 Takeo[11]中引用出有限元素之聲壓方程式,而 Tan[12]中討 論了藉由促動器(Actuator)主動控制對平板之聲場的影響。文獻[13]中陳述了 關於揚聲器量測的各參數的討論,文獻[14]中則討論了傳統半圓型或波浪型 的彈性懸邊的設計及研製過程。 1.3 研究方法 本文主要是將數種材料以不同製程參數製作的彈性懸邊作為揚聲器振 動板的邊界,目的在降低揚聲器系統之第一自然頻率,增加揚聲器在各頻 率下的穩定度,並提升低頻之聲壓值,並分析不同長寬比的振動板與不同 參數的音圈對聲壓曲線的影響。 實驗方面,製作不同材質的懸邊裝置在揚聲器上,探討各種懸邊對揚 聲器穩定度的影響,並將不同製程懸邊所製作出來的揚聲器利用 LMS 聲壓 量測系統量測聲壓及阻抗,並觀察聲壓曲線的趨勢與平均感度,作為與理 論分析的對照。 針對實務的振動系統而言,大多採用有限元素法求取振動系統的運動 方程式,但是當系統是非常複雜的結構時,就必須藉由電腦輔助以快速且 有效的獲得分析結果,ANSYS 即為此一類型的套裝工程分析軟體。本文將 以 ANSYS 有限元素分析軟體建立一個與實體幾何相近的模型,進而分析結
構之自然頻率與模態,接著再進行簡諧激振分析,計算出不同激振頻率下 振動板所有節點之振幅及相位角,輸入由 Fortran 程式所寫成的聲壓計算公 式並由 LMS 聲壓軟體繪製出模擬聲壓曲線,將模擬曲線與實驗結果相對 照,以驗證模擬的正確性,再以 ANSYS 有限元素軟體模擬改變振動板的長 寬比與音圈參數,找出有利於聲壓曲線更平滑、平均感度更高的設計參數。
第二章 揚聲板的振動分析 本章討論平板的振動分析,其中平板理論有用到一階剪變形和多層一 階剪變形兩種。由於本文振動板厚度均在 1mm 以內,因此使用一階剪變形 理論來進行分析,如有部份區域的振動板採用碳纖維加勁,則由於面層的 碳纖維強度與心層的巴沙木強度相差許多,因此使用多層一階剪變形理論 來進行分析。 2.1 一階剪變形平板理論 2.1.1 基本假設
一階剪變形理論 FSDT(The first order shear deformation theory)來作為積 層板分析的基礎,其基本假設如下: 1.板的長、寬為板的厚度的 15 倍以上。 2.板的截面變形後仍保持平面。 3.厚度仍保持不變,即εz=0。 4.板的變形量 u, v, w 很小。 其位移場的假設如下: t) y, (x, zθ t) y, (x, u u = 0 + x t) y, (x, zθ t) y, (x, v v= 0 + y t) y, w(x, w= (2.1) 其應變場如下: 0 1 0 1 x x, x 0, 1 u zθ ε zκ x u ε = + = + ∂ ∂ = 0 2 0 2 y y, y 0, 2 v zθ ε zκ y v ε = + = + ∂ ∂ =
4 w,y θy z v y w ε = + ∂ ∂ + ∂ ∂ = 5 w,x θx z u x w ε = + ∂ ∂ + ∂ ∂ = 0 6 0 6 x y, y x, x 0, y 0, 6 (u v ) z(θ θ ) ε zκ x v y u ε = + + + = + ∂ ∂ + ∂ ∂ = (2.2) 其中 、 、w 分別代表積層板在 x、y、z 方向之位移量,t 代表時間, 、 分別代表積層板中間面在 x、y 方向上的位移,而 u v u0 0 v x w θx ∂ ∂ − = 、 y w θy ∂ ∂ − = 是 剪應變所產生之旋轉量(Shear rotation)。 而中間面之曲率為: y x w 2 κ y w κ x w κ 2 0 6 2 2 0 2 2 2 0 1 ∂ ∂ ∂ − = ∂ ∂ − = ∂ ∂ − = 2.1.2 應力應變關係 複合材料積層板是由多層薄片(Lamina) 疊合而成的,而每一層的纖維 排列方向與積層板主軸方向夾一θn角,則在材料座標系 座標的構成 方程式(Constitutive equation)為: Z Y X′ ′ ′ 1 5 n 5 5 1 5 Q ε σ × = × × (2.3) 其中
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = × n 55 n 44 n 66 n 22 n 12 n 12 n 11 n 5 5 Q 0 0 0 0 0 Q 0 0 0 0 0 Q 0 0 0 0 0 Q Q 0 0 0 Q Q Q (2.4) 且 ) v v (1 E Q 21 12 n 1 n 11 − = ) v v (1 v E Q 21 12 12 n 2 n 12 − = ) v v (1 E Q 21 12 n 2 n 22 − = n 12 n 66 G Q = n 23 n 44 G Q = n 13 n 55 G Q = (2.5) 其中Q 為彈性係數,nij 為楊氏模數, 為蒲松比, 為剪力模數。 n i E vij n ij G 由座標轉換的關係,求出積層平板主軸方向XYZ座標系的應力及應變 關係 n 1 5 n 5 5 n 1 5 t 5 5 1 n n 5 5 1 n n 1 5 (T ) Q [(T ) ] ε Q ε σ × × × × − × − × = = (2.6) 其中 ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − − − = − × C S 0 0 0 S C 0 0 0 0 0 S C SC SC 0 0 2SC C S 0 0 2SC S C ) T ( 2 2 2 2 2 2 1 5 5 n (2.7) 其中 ) cos(θ C= n ) sin(θ S= n 而
⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = × n 55 n 45 n 45 n 44 n 66 n 26 n 16 n 26 n 22 n 12 n 16 n 12 n 11 n 5 5 Q Q 0 0 0 Q Q 0 0 0 0 0 Q Q Q 0 0 Q Q Q 0 0 Q Q Q Q (2.8) 其中 n 5 5 Q × 矩陣與 n 5 5 Q × 矩陣間之關係為: 4 n 22 2 2 n 66 n 12 4 n 11 n 11 Q C 2(Q 2Q )C S Q S Q = + + + ) S (C Q S )C Q 4 Q Q ( Q n 4 4 12 2 2 n 66 n 22 n 11 n 12 = + − + + n 4 22 2 2 n 12 n 12 4 n 11 n 22 Q S 2(Q 2Q )C S Q C Q = + + + n 3 66 n 12 n 11 3 n 66 n 12 n 11 n 16 (Q Q 2Q )C S (Q Q 2Q )CS Q = − − + − + Q (Q Q 2Q )CS (Q Q 2Qn )C3S 66 n 12 n 11 3 n 66 n 12 n 11 n 26 = − − + − + Q (Q Q 2Q 2Q )C S Qn (C4 S4) 66 2 2 n 66 n 12 n 22 n 11 n 66 = + − − + + n 2 55 2 n 44 n 44 Q C Q S Q = + Q Q CS Qn44CS n 55 n 45 = − n 2 55 2 n 44 n 55 Q S Q C Q = + (2.9) 2.1.3 構成方程式(Governing equation) 因複合材料積層板與一般平板理論有所差異,其差異性在於複合積層 平板必須考慮到每一層薄片的力學性質,然後再將它們累加起來,始可得
到合應力及合應力矩(Stress and moment resultants),複合積層平板之應力與
應力矩為: 0 j ij 0 j ij 0 j N 1 n N 1 n Z Z 0 j n ij Z Z ij n ij i Q ε dz Q (ε zκ )dz A ε B κ N n1 n 1 n n + = + = =
∑
∫
∑∫
= = + +∑∫
= +∑∫
= + + = + = = N 1 n Z Z N 1 n Z Z 0 j ij 0 j ij 0 j 2 0 j n ij ij n ij i 1 n n 1 n n κ D ε B )dz κ z (zε Q dz ε zQ M (2.10) 其中∑
∫
∑
∫
∑
∫
= + − = + − = + − − = = − = = − = = N 1 n 3 n 3 1 n n ij h/2 h/2 2 n ij ij N 1 n 2 n 2 1 n n ij h/2 h/2 n ij ij N 1 n n 1 n n ij h/2 h/2 n ij ij ) z (z Q 3 1 dz z Q D ) z (z Q 2 1 zdz Q B ) z (z Q dz Q A (i,j=1,2,6) (2.11) n z 為第n層薄板之表面到中間面的距離, (i,j=1,2,6)分別代表拉伸勁度矩陣(Extensional stiffness matrix)、偶合勁度矩陣(Coupling stiffness
matrix)、彎曲勁度矩陣(Bending stiffness matrix),將合力及合力矩寫成矩陣 形式如下: ij ij ij,B ,D A ⎪ ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎪ ⎨ ⎧ + + ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎪ ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎪ ⎨ ⎧ x y, y x, y y, x x, x 0, y 0, y 0, x 0, 66 26 16 66 26 16 26 22 12 26 22 12 16 12 11 16 12 11 66 26 16 66 26 16 26 22 12 26 22 12 16 12 11 16 12 11 6 2 1 3 2 1 θ θ θ θ v u v u D D D B B B D D D B B B D D D B B B B B B A A A B B B A A A B B B A A A M M M N N N (2.12) 2.2 多層一階剪變形平板理論 由於複合材料三明治板的面層與心層間材料性質差異性很大,面層是 由不同纖維方向的複合材料層板所組成,因此面層與夾心層間在同一方向 的勁度會有很大的差異,所以複合材料三明治板沿整個厚度方向的位移並 不如一階剪變形位移場平板理論所假設的整個板厚度為一斜率相同的直 線,故我們將整個複合材料三明治板的每一層當作一個一階剪變形的位移
場,並考慮每一層與層之間位移的連續性(如圖2.1)。 2.2.1 位移與應變 多層一階剪變形平板理論是將複合材料三明治層板的每一層當作一個 一階剪變形的平板來分析。假設層板任一層的位移場為 ) M , , 2 , 1 i ( ) y , x ( w ) y , x ( w w ) y , x ( z ) y , x ( v v ) y , x ( z ) y , x ( u u ) i ( ) i ( 0 ) i ( ) i ( y ) i ( ) i ( 0 ) i ( ) i ( x ) i ( ) i ( 0 ) i ( L = = = φ + = φ + = (2.13) 其中, 、 、 分別為任一層在平板參考座標中x、y、z上之位 移分量,而 、 、 分別代表任一層之層板中間面在方向 x、y、z 上 之位移量, 、 則分別代表為任一層垂直於x、y軸之截面的旋轉量。 ) i ( u (i) v (i) w ) i ( 0 u (i) 0 v (i) 0 w ) i ( x φ (i) y φ 本文中將層數設為三層,如圖2.1 所示,上下面層為碳纖維,中心層為 巴沙木,依據式(2.13),並考量每一層於交界面上的位移必須連續,可得每 一層的位移場為 ) y , x ( w ) y , x ( w w ) y , x ( z ) y , x ( v v ) y , x ( z ) y , x ( u u ) 1 ( ) 1 ( ) 1 ( y ) 1 ( ) 1 ( 0 ) 1 ( ) 1 ( x ) 1 ( ) 1 ( 0 ) 1 ( = = φ + = φ + = (2.14a) ) y , x ( w w z t 2 1 t 2 1 v z v v z t 2 1 t 2 1 u z u u ) 2 ( ) 2 ( y ) 2 ( ) 2 ( y ) 2 ( ) 1 ( y ) 1 ( ) 1 ( 0 ) 2 ( y ) 2 ( ) 2 ( 0 ) 2 ( ) 2 ( x ) 2 ( ) 2 ( x ) 2 ( ) 1 ( x ) 1 ( ) 1 ( 0 ) 2 ( x ) 2 ( ) 2 ( 0 ) 2 ( = φ + φ + φ + = φ + = φ + φ + φ + = φ + = (2.14b)
) y , x ( w w z t 2 1 t 2 1 v z v v z t 2 1 t 2 1 u z u u ) 3 ( ) 3 ( y ) 3 ( ) 3 ( y ) 3 ( ) 1 ( y ) 1 ( ) 1 ( 0 ) 3 ( y ) 3 ( ) 3 ( 0 ) 3 ( ) 3 ( x ) 3 ( ) 3 ( x ) 3 ( ) 1 ( x ) 1 ( ) 1 ( 0 ) 3 ( x ) 3 ( ) 3 ( 0 ) 3 ( = φ + φ − φ − = φ + = φ + φ − φ − = φ + = (2.14c) 其中 (i) t 為第i層厚度。 假設厚度方向應變仍保持不變εz=0,應變可表示為:
[ ]
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎧ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂∂ ∂ = ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎨ ⎧ ε ε ε ε ε = ε y w z v x w z u x v y u y v x u ) i ( ) i ( ) i ( ) i ( ) i ( ) i ( ) i ( ) i ( yz xz xy y x i=1,2,3 (2.15) 其中 x ε 、 y ε 為平板參考座標中x、y上之應變,ε 為平板參考座標中x、xy y上之剪應變,ε 為平板參考座標中y、z上之剪應變,yz 為平板參考座標 中x、z上之剪應變。 xz ε 對於層板的橫向振動(Transverse vibration)問題,可假設u(1)(x,y) 0 0 = 及 ,並將式(2.14)代入式(2.15),可得到應變張量 : 0 ) y , x ( v(1) 0 =[ ]
ε[ ]
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎧ ∂ ∂ + φ ∂ ∂ + φ ∂ φ ∂ + ∂ φ ∂ ∂ φ ∂∂ φ ∂ = ε y wx w x z y z y z x z ) 1 ( y ) 1 ( x ) 1 ( y ) 1 ( ) 1 ( x ) 1 ( ) 1 ( y ) 1 ( ) 1 ( x ) 1 ( ) 1 ( (2.16a)[ ]
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎧ ∂ ∂ + φ ∂ ∂ + φ ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∂ φ ∂ + ∂ φ ∂ + ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∂ φ ∂ + ∂ φ ∂ + ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∂ φ ∂ + ∂ φ ∂ ∂ φ ∂ + ∂ φ ∂ + ∂ φ ∂ ∂ φ ∂ + ∂ φ ∂ + ∂ φ ∂ = ε y w x w x y z x y t 2 1 x y t 2 1 y z y t 2 1 y t 2 1 x z x t 2 1 x t 2 1 ) 2 ( y ) 2 ( x ) 2 ( y ) 2 ( x ) 2 ( ) 2 ( y ) 2 ( x ) 2 ( ) 1 ( y ) 1 ( x ) 1 ( ) 2 ( y ) 2 ( ) 2 ( y ) 2 ( ) 1 ( y ) 1 ( ) 2 ( x ) 2 ( ) 2 ( x ) 2 ( ) 1 ( x ) 1 ( ) 2 ( (2.16b)[ ]
⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎧ ∂ ∂ + φ ∂ ∂ + φ ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∂ φ ∂ + ∂ φ ∂ + ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∂ φ ∂ + ∂ φ ∂ − ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ∂ φ ∂ + ∂ φ ∂ − ∂ φ ∂ + ∂ φ ∂ − ∂ φ ∂ − ∂ φ ∂ + ∂ φ ∂ − ∂ φ ∂ − = ε y wx w x y z x y t 2 1 x y t 2 1 y z y t 2 1 y t 2 1 x z x t 2 1 x t 2 1 ) 3 ( y ) 3 ( x ) 3 ( y ) 3 ( x ) 3 ( ) 3 ( y ) 3 ( x ) 3 ( ) 1 ( y ) 1 ( x ) 1 ( ) 3 ( y ) 3 ( ) 3 ( y ) 3 ( ) 1 ( y ) 1 ( ) 3 ( x ) 3 ( ) 3 ( x ) 3 ( ) 1 ( x ) 1 ( ) 3 ( (2.16c) 2.2.2 應力與應變考慮一個單層複合材料平板,在材料主軸方向的應力與應變關係可以 表示如下: ⇒ ⎪ ⎪ ⎪ ⎭ ⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎩ ⎪ ⎪ ⎪ ⎨ ⎧ γ γ γ ε ε ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎪ ⎪ ⎪ ⎭ ⎪ ⎪ ⎪ ⎬ ⎫ ⎪ ⎪ ⎪ ⎩ ⎪ ⎪ ⎪ ⎨ ⎧ τ τ τ σ σ σ 12 13 23 2 1 66 55 44 33 23 13 23 22 12 13 12 11 12 13 23 3 2 1 0 Q 0 0 0 0 0 0 Q 0 0 0 0 0 0 Q 0 0 0 0 0 0 Q Q Q 0 0 0 Q Q Q 0 0 0 Q Q Q (2.17) 下標1、2代表主軸方向,也就是各層纖維的方向,而 為彈性係數, 定義如下: ij Q Δ ν ν − = 3 2 32 23 11 E E 1 Q Δ ν ν + ν = Δ ν ν + ν = 3 1 13 32 12 3 2 23 31 21 12 E E E E Q Δ ν ν + ν = Δ ν ν + ν = 2 1 23 12 13 3 2 32 21 31 13 E E E E Q Δ ν ν − = 3 1 31 13 22 E E 1 Q Δ ν ν + ν = Δ ν ν + ν = 3 1 13 21 23 3 1 31 12 23 23 E E E E Q Δ ν ν − = 2 1 21 12 33 E E 1 Q 23 44 G Q = 、Q55 =G13、Q66 =G12 其中 3 2 1 13 32 21 23 32 31 13 21 12 E E E 2 1−ν ν −ν ν −ν ν − ν ν ν = Δ (2.18) 因複合材料是由不同角度的層板疊合而成,所以必須將主軸方向的材 料性質,轉換至層板的座標 x-y 方向,如圖(2.2)所示。則在材料主軸座 標系的構成方程式 (Constitutive equation) 為:
{ }
σ =[ ]
Q{ }
ε (2.19) 展開寫成 ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ γ γ γ ε ε ⋅ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ = ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ σ σ σ σ σ σ xy xz yz y x 66 36 26 16 55 45 45 44 36 33 23 13 26 23 22 21 16 13 12 11 xy xz yz z y x 0 Q 0 0 Q Q Q 0 Q Q 0 0 0 0 Q Q 0 0 0 Q 0 0 Q Q Q Q 0 0 Q Q Q Q 0 0 Q Q Q (2.20) 其中 4 22 2 2 66 12 4 11 11 Q C 2(Q 2Q )C S Q S Q = + + + ) S C ( Q S C ) Q 4 Q Q ( Q 12 4 4 2 2 66 22 11 12 = + − + + 2 23 2 13 13 Q C Q S Q = + 4 22 2 2 66 12 4 11 22 Q S 2(Q 2Q )C S Q C Q = + + + 2 13 2 23 23 Q C Q S Q = + 33 33 Q Q = 3 66 22 12 3 66 12 11 16 (Q Q 2Q )C S (Q Q 2Q )CS Q = − − + − + S C ) Q 2 Q Q ( CS ) Q 2 Q Q ( Q 3 66 22 12 3 66 12 11 26 = − − + − + CS ) Q Q ( Q36 = 13 − 23 ) S C ( Q S C ) Q 2 Q 2 Q Q ( Q 4 4 66 2 2 66 12 22 11 66 = + − − + + 2 55 2 44 44 Q C Q S Q = + 2 44 2 55 55 Q C Q S Q = + CS ) Q Q ( Q45 = 55 − 44 i i , S sin cos C= θ = θ (2.21) 而主軸方向與 x-y 方向夾角為θ
(圖2.2),從整個複合材料積層板而 言,必須考慮各單層板之應力沿厚度方向積分,可得複合材料積層板的合 力與合力矩如,如圖2.3所示,其結果如下:(
)
(
)
∑∫
∑∫
∑∫
∑∫
∑∫
= = = = = + + + + + κ + ε = κ + ε = σ = σ σ = κ + ε = κ + ε = σ = N 1 n z z 0 j ij 0 j ij 0 j 2 0 j n ij N 1 n z z n i i N 1 n z z n yz n xz y x N 1 n z z 0 j ij 0 j ij 0 j 0 j n ij N 1 n z z n i i 1 n n 1 n n 1 n n 1 n n 1 n n D B dz z z Q zdz M dz ) , ( ) Q , Q ( B A dz z Q dz N (2.22) 其中 ,zn則為第 n 層板之下表面至中性面的距離,A ij、B ij、Dij 分別代表拉伸勁度矩陣(Extensional stiffness matrix)、偶合勁度矩陣
(Coupling stiffness matrix)、彎曲勁度矩陣(Bending stiffness matrix),KS為剪
力修正常數,關係式如下: 6 , 2 , 1 j , i = ⎭ ⎬ ⎫ ⎩ ⎨ ⎧ θ ε ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ = ⎭ ⎬ ⎫ ⎩ ⎨ ⎧ 0 D B B A M N
[ ]
(
∑
= − − = N 1 k 1 k k k ij ij Q z z A)
(2.23)[ ]
(
∑
= − − = N 1 k 2 1 k 2 k k ij ij Q z z 2 1 B)
(2.24)[ ]
(
∑
= − = N 1 k 3 k 3 k k ij ij Q z z 3 1 D)
i,j=1,2,6 (2.25)[ ]
(
∑
= − − = N 1 k 1 k k k ij s ij k Q z z A)
i,j=4,5 (2.26) 2.3 複合材料三明治板的彈性支承 對於在平板邊界上有彈性支承的情況,如本文揚聲器系統中的懸邊元 素,可模擬成如圖(2.4)所示的連續彈簧,將懸邊元素提供的彈性力等效成 橫向(Translation)及旋轉(Rotation)彈簧加以表示,而其所引起的應變能為: dS w V 2 1 U s n T =∫
(2.27)∑ ∫
= φ = 3 1 k S ) i ( n R M dS 2 1 U (2.28) 其中,UT為橫向彈簧所引起的應變能,UR為旋轉彈簧所引起的應變能, S為彈性支承的邊界,Vn為橫向力,Mn為彎矩。 橫向力Vn及彎矩Mn可分別表示為: w K Vn = Ln (2.29)∑
= φ = 3 1 i ) i ( Rn n K M (2.30) 將式(2.29)和(2.30)代入式(2.27)及(2.28)可將應變能表示為:∫
∫
∫
∫
= = = = + + + = a 0 y b 2 4 L a 0 y 0 2 3 L b 0 x a 2 2 L b 0 x 0 2 1 L T dx w 2 K dx w 2 K dy w 2 K dy w 2 K U (2.31)( )
( )
( )
φ +( )
φ ⎥⎦⎤ + ⎢⎣ ⎡ φ + φ =∫
∫
∑
∫
∫
= = = = = a 0 y b 2 ) i ( y 4 R a 0 y 0 2 ) i ( y 3 R 3 1 i b 0 x a 2 ) i ( x 2 R b 0 x 0 2 ) i ( x 1 R R dx 2 K dx 0 2 K dy 2 K dy 2 K U (2.32) 其中, 為平板四邊的橫向剛性( Translational stiffness), 為平板四邊 的旋轉剛性(Rotational stiffness)。 Ln K KRn 2.4 應變能與動能 考慮任一個單層複合材料層板,應變能可表示為:[ ] [ ]
∫
σ ε = V ) i ( T ) i ( ) i ( p dV 2 1 U i=1,2,3 (2.33) 應用上式與(2.14)、(2.19)式,可求得 )、 及 ,複合材料層板彎 曲的應變能為: 1 ( p U (2) p U (3) p U∑
= = 3 1 k ) k ( p p U U (2.34)總應變能U為複合材料層板彎曲的應變能( )與彈性支承引起的應變 能( 、 )之總和。 p U T U UR R T p U U U U = + + (2.35) 每一單層板的動能 (i)為: T dv t w t v t u T i i v i i
∫
⎥⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ ∂ ∂ + ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ∂ ∂ + ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ∂ ∂ = 2 2 ) ( 2 ) ( ) ( ) ( 2ρ
; (2.36) 3 , 2 , 1 i= 應用上式與(2.14)式,可求得 、 及 ,複合材料層板的最大動 能為: ) 1 ( T T(2) T(3)∑
= = 3 1 i ) i ( T T (2.37) 複合材料層板的總能量泛函Π可表示為: Π = U−T (2.38) 2.5 瑞雷-黎次法(Rayleigh-Ritz method)利用Rayleigh-Ritz method ,假設位移函數(Deflection function)及橫切面
轉角函數(Cross-sectional rotation function)共有七個分別是w、 (1)、
x φ (1) Y φ 、 、 、 及 ,每一個函數以無因次化可表示為: ) 2 ( x φ (2) Y φ (3) x φ (3) Y φ
∑∑
= = η Ψ ξ Φ = η ξ 1 1 1 1 1 1 1 1 I 1 i J 1 j yj xi ) 1 ( j i ( ) ( ) C ) , ( w (2.39a)∑∑
= = η Ψ ξ Φ = η ξ φ 2 2 2 2 2 2 2 2 I 1 i J 1 j yj xi ) 2 ( j i ) 1 ( x ( , ) C ( ) ( ) (2.39b)∑∑
= = η Ψ ξ Φ = η ξ φ 3 3 3 3 3 3 3 3 I 1 i J 1 j yj xi ) 3 ( j i ) 1 ( y ( , ) C ( ) ( ) (2.39c)∑∑
= = η Ψ ξ Φ = η ξ φ 4 4 4 4 4 4 4 4 I 1 i J 1 j yj xi ) 4 ( j i ) 2 ( x ( , ) C ( ) ( ) (2.39d)∑∑
= = η Ψ ξ Φ = η ξ φ 5 5 5 5 5 5 5 5 I 1 i J 1 j yj xi ) 5 ( j i ) 2 ( y ( , ) C ( ) ( ) (2.39e)∑∑
= = η Ψ ξ Φ = η ξ φ 6 6 6 6 6 6 6 6 I 1 i J 1 j yj xi ) 6 ( j i ) 3 ( x ( , ) C ( ) ( ) (2.39f)∑∑
= = η Ψ ξ Φ = η ξ φ 7 7 7 7 7 7 7 7 I 1 i J 1 j yj xi ) 7 ( j i ) 3 ( y ( , ) C ( ) ( ) (2.39g) 其中 (1) 、 、 、 、 、 、 為未定係數,ξ、 j i11 C (2) j i22 C (3) j i33 C (4) j i44 C (5) j i55 C (6) j i66 C (7) j i77 C η為無 因次參數,且x = aξ,y = bη, x Φ 、Ψy係由Gram-Schmidt正交化法所產生的多項式函數(Polynomial functions),亦是所謂之形狀函數(Shape function)。
以Φx為例:
( ) (
1)
0 1 ξ = ξ−B Φ Φ( ) (
ξ = ξ−)
Φ − Φ( )
ξ Φk Bk k−1 Ck k−2 k≥2 (2.40) 其中∫
∫
− − − − ξ ξ Φ ξ ξ Φ ξ = 0.5 5 . 0 2 1 xk 5 . 0 5 . 0 2 1 xk k d ) ( d ) ( B∫
∫
− − − − − ξ ξ Φ ξ ξ Φ ξ Φ ξ = 0.5 5 . 0 2 2 xk 5 . 0 5 . 0 xk 1 xk 2 n d ) ( d ) ( ) ( C (2.41) 此多項式函數必符合: ij 5 . 0 5 . 0 Φxi(ξ)Φxj(ξ)dξ=δ∫
− (2.42) 其中 ij δ 是Kronecker, 令 ( ) 1及 1 x ξ = Φ Ψy1(η)=1 為求未定係數,令 0 C(1) j i11 = ∂ Π ∂ , C(2) 0 j i22 = ∂ Π ∂ , C(3) 0 j i33 = ∂ Π ∂ , C(4) 0 j i44 = ∂ Π ∂ , 0 C(5) j i55 = ∂ Π ∂ , C(6) 0 j i66 = ∂ Π ∂ , C(7) 0 j i77 = ∂ Π ∂ , (2.43) 將(2.39)代入上式可得一典型的特徵值方程:(
K−ω2M)
{ } { }
C = 0 (2.44) 其中,[M]為質量矩陣,[K]為勁度矩陣,求解上式即可得到複合材料層板振動的自然頻率ω。 2.6 結構系統之阻尼 假設次系統具有比例阻尼(Proportional damping) [C]=α[M]+β[K] (2.45) 其中 ω ξ = β ξω = α , 2 ,ξ為阻尼比。 (2.46) 構成方程式為
[M]{X&&}+[C]{X&}+[K]{X}={F}sinϖt (2.47) 其中{X}為位移向量。 因為[M]、[K]、[C]為實對稱矩陣,所以我們將其對角化為: ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = Ψ Ψ O O 0 m 0 ] ][ M [ ] [ 2 r T ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = Ψ Ψ O O 0 k 0 ] ][ K [ ] [ 2 r T (2.48) ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎣ ⎡ = Ψ Ψ O O 0 c 0 ] ][ C [ ] [ 2 r T 我們將{X}=[Ψ]{x}代入(2.37)式可改寫成為: [m ]{x} [c ]{x} [k ]{x} [ ]T{F}sin t r r r && + & + = Ψ ϖ (2.49) 所以 ⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ η + η − η + η + η − η − = 2 r 2 2 r 2 r 2 2 2 r r r ) c 2 ( ) 1 ( c 2 j ) c 2 ( ) 1 ( 1 k F x (2.50) 其中
η ω ω
= / r位移向量{X}=[Ψ]{x}=
∑
{ }
(2.51) = Φ n r r rx 1 2.7 聲壓值計算 將 2.51 式所求得的位移向量及所對應的頻率代入以下所推導的聲壓公 式,即可求得各頻率下的聲壓值。在聲壓計算或實驗的量測,因聲壓大小 與量測的距離、輸入揚聲器的功率有關,若這兩個條件不同而作聲壓比較 是沒有意義的,故本文在理論分析及實驗量測距離都以距平板中心為一公 尺,輸入功率為一瓦為基準。 2.7.1 聲壓公式推導 對平面波而言,一維波動方程式為: 2 2 2 2 t p c 1 x p ∂ ∂ = ∂ ∂ 其中,c 為聲音之速度,在室溫 240 C 時,c=343 m/s。 利用一維波動問題 D’Alembert 解,可求得: ) ct x ( f ) ct x ( f ) t , x ( p = 1 − + 2 + f1,f2為任意函數,對於諧波聲場,可直接求解: ~p 0 c dx p d 2 2 2 = ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ ω + ( (2.52) 令k =ω/c,為波數(Wave number),可求得: jkx jkx e B~ e A~ ) x ( p ~ = − + (2.53) 還原成時間域的表示法: j( t kx) j( t kx e B ) ~ e A~ ) t , x ( p = − ϖ− + ϖ− (2.54) 其中A~, 為未定複數,對於點聲源而言,3-D 波動方程式在球面座標B~(Spherical coordinate)為: 2 2 2 2 2 2 2 2 2 2 2 2 2 t p c 1 p tan r 1 r p r 2 p sin r 1 p r 1 r p ∂ ∂ = θ ∂ ∂ θ + ∂ ∂ + φ ∂ ∂ θ + θ ∂ ∂ + ∂ ∂ (2.55) 因為點聲源是一個完全對稱的聲源,所以我們可以預期它所產生的聲 場 應 該 不 會 有 方 向 性 , 所 以 可 忽 略 角 度 相 關 的 項 式 , (2.55) 式 可 簡 化成: 2 2 2 2 2 t p c 1 r p r 2 r p ∂ ∂ = ∂ ∂ + ∂ ∂ (2.56) r 為點聲源距量量測的距離 上式可改寫成 ) rp ( t c 1 r 1 ) rp ( r r 1 2 2 2 2 2 ∂ ∂ = ∂ ∂ (2.57) 故 2 2 2 2 2 t ) rp ( c 1 r ) rp ( ∂ ∂ = ∂ ∂ (2.58) 觀察(2.58)的型式和前面的一維波動方程式式一樣的,所以其解為: rp=f1(r-ct)+f2(r+ct) 所以 f (r ct) r 1 ) ct r ( f r 1 ) t , r ( p = 1 − + 2 + (2.59) 由於一個點聲源只會發出外傳波(Outgoing wave),所以 ) ct r ( f r 1 ) t , r ( p = 1 − (2.60) 對於簡諧球面波(Harmonic spherical wave)而言
( t kr j e ) r A~ ) t , r ( p = ϖ− (2.61) 對微小的振動板δs而言,距離振動板 r 的聲壓 p(r,t)可由 Rayliegh’s first integral 計算出來,積分型式如下:
( t kr j d air e r 2 s u j ) t , r ( p ⎟ ϖ− ⎠ ⎞ ⎜ ⎝ ⎛ π δ ϖρ = ) (2.62) 其中 ud 為振動板元素表面之速度,r 為量測點至振動板元素之距離, 為空氣密度,j= 3 air =1.1614kg/m ρ − ,因此對一塊振動板而言,聲壓可寫1 為: ds r e ) y , x ( u e r 2 j ) t , r ( p jkr d S t j air − ϖ
∫
π ϖρ = (2.63) 設振動板表面元素之振幅大小 j( t kr) Ae ) t , r ( w = ϖ− ,則 ) t , r ( w j Ae j ) y , x ( u i( t kr) d = ϖ = ϖ − ϖ 所以 ds r e ) t , r ( w e r 2 ) t , r ( p jkr S t j air 2 −∫
− = ϖ π ρ ϖ (2.64) 因以振動板之聲壓可改寫為 i ) kr ( j i 2 air r S e || ) t , r ( w || 2 P ⎟ i i Δ ⎠ ⎞ ⎜ ⎝ ⎛ π ϖ ρ − =∑
θ− (2.65) 其中θi為相角(Phase angle) 2 i i i 1 m k c tan ϖ − ϖ = θ − 在實務上,通常量測聲音使用的尺度主要原因是聲音的動態範圍非常 之大,同時人耳對音量的感覺也是比較接近對數尺度。所以定義聲壓位準 (Sound pressure level)) P P log( 20 SPL ref rms = ….(dB) (2.66) 其中 Prms為量測點聲壓之均方根值為: 2 / 1 2 / T 2 / T 2 rms |p(r,t)| dt T 1 P ⎥⎦ ⎤ ⎢⎣ ⎡ =
∫
− (2.67) 聲壓參考值 Pref=2x10 -5 Pa一般在量測聲壓的儀器看到的聲壓曲線即為頻率響應曲線,縱軸為 (dB),橫軸為激振頻率。
第三章 ANSYS 有限元素模型之建構與分析 根據 2.6 節的位移向量(2.51 式)可以得到振動板上各點在每個分析頻率 下的振幅、相角,本節將以有限元素軟體 ANSYS 建立揚聲器的模型,並將 模型進行簡諧激振分析求解 2.6 節的構成方程式(2.47 式),以求得振動板上 各點在每個分析頻率下的振幅、相角,再代入 2.65 公式計算聲壓,即可得 到各頻率下的聲壓值。 3.1 模型建構 文獻[7]中以 Shell 99 薄殼元素模擬複合材料薄板可獲得良好之分析結 果,此元素可設定複材板之疊層數目及各疊層之纖維角度及各疊層之材料 常數,故本文中也使用 Shell 99 薄殼元素來模擬複合材料板。根據文獻[24] 分析,彈性懸邊可以 Shell 93 薄殼元素來模擬,此元素適合用來分析柔軟之 彎曲形薄殼,且容許大變形發生,與懸邊之實際情況相當符合,故選取此 元素來模擬彈性懸邊。本文主要分析小型揚聲器,因此音圈的強度相對於 振動板的強度是必須考慮的,所以在模型建構時以 Shell 99 薄殼元素來模擬 音圈。 3.1.1 ANSYS 模型建立步驟 前處理部分:
1. Preprocessor → Element type:選擇振動板、音圈元素為 Shell 99,懸 邊元素為 Shell 93。
3. Preprocessor → Material Props → Material Models:設定振動板元素之 各材料性質。 4. Preprocessor → Modeling:由點、線、面建立揚聲器系統的模型外觀。 5. Preprocessor → MeshTool:選擇元素參數、材料性質、各元素之尺寸 大小,並分割元素。 到此前處理即算完成,接下來可以做模態分析或是聲聲分析部分。 首先說明模態分析部分:
1. Solution → Analysis Type → New Analysis:選擇分析型態,自然頻率 模態分析點選“Modal”。
2. Solution → Analysis Type →Analysis Options:No. of modes to extract 為要分析的模態個數。
3. Solution → Define Loads → Apply → Structural → Displacement → On Nodes:限制懸邊元素另一端節點的全部自由度。
4. Solution → Solve → Current Ls:求解。
5. General Postproc →Results Summary:列出所有的自然頻率。
6. General Postproc →Read Results →By Pick:選出想看 Mode Shapes 的 自然頻率。
7. General Postproc →Plot Results →Contour Plot →Nodal Solu:再選 Nodal Solution →Z-Component of displacement,即可得到自然頻率及
模態。
如果要分析聲壓,在 ANSYS 部分如下:
1. Solution → Analysis Type → New Analysis:選擇分析型態,簡諧頻率 響應分析點選“Harmonic”。
2. Solution → Define Loads → Apply → Structural → Displacement → On Nodes:限制彈簧元素另一端節點的全部自由度。
3. Solution →Define Loads → Apply → Structural → Force/Moment → On Nodes:在位於音圈位置上的節點施予 Z 方向且相角為零之推力。 4. Solution → Load Step Opts → Time/Frequenc → Damping:輸入系統阻
尼
α
、β
的值。5. Solution → Load Step Opts → Time/Frequenc → Freq and Substps:輸入 欲分析頻率響應之頻寬。
6. Solution → Solve → Current Ls:求解。
7. TimeHist Postpro → List Variables:輸出振動板模型全部節點的振幅及 相角。
由上面簡諧激振分析可得到以下資料: 1.節點編號 NodeNo(i)及位置 x(i),y(i),z(i)
2.第 i 節點在第 j 頻率的振幅 Amplitude(i,j)及相角 Phase(i,j)
ANSYS 分析模擬的各項參數和模擬音圈激振的施力皆可由實驗取 得。以下的討論將使用到參數值中的 Mms、Cms 和 BL 值,其實驗過程將 在第四章中做敘述,而系統α-damping 和β-damping 也可從阻尼量測實驗 中陳述的方法取得。 3.2.1 系統的彈性係數 模擬在振動板中央施與集中力,並將振動板材料常數調整為剛體,利 用 ANSYS 進行揚聲器模型的靜態分析,將模擬得到的位移與施力利用虎克 定律計算,求得不同懸邊之彈性係數,並將 ANSYS 靜態分析計算出來的彈 性係數與實驗參數比對,兩者誤差極小,因此彈性懸邊元素的模擬是可信 的,其中實驗所得的彈性係數可由參數表 4.3~4.5 中的 Cms 值計算出來, Cms 值是系統每牛頓力可產生的位移。 3.2.2 激振力的給定 激振力計算公式為 F=BLI,其中 I 值為流過激振器線圈的電流,B 為磁 通密度,L 為線圈長度;BL 值也可由量測參數表 4.3~表 4.5 中取得,由於 線圈阻抗 R 可由三用電錶量出,而本文量測聲壓時使用標準的量測功率為 一瓦,再由公式 W=I2R(其中 W 為功率,I 為電流,R 為阻抗)可計算出 I 值, 因此激振力即可計算出來。 3.2.3 阻尼比的給定 本 文 依 照 4.2 節 實 驗 方 法 取 了 兩 個 頻 率 來 計 算 α -damping 和 β
-damping,並藉由 Bandwidth method 來計算各激振頻率之系統阻尼比。其 中 Peak 為某一共振頻率相對應之振幅, 及f1 f2為曲線和
2
peak
值之交點所 對應的頻率,利用下式求得共振頻率之阻尼比如下式: 2 2 1 1f
f
f
f
ξ
=
−
+
(5.1) 再利用 Rayleigh damping 將頻率響應實驗計算得到的系統阻尼比,利用下式 求得系統的α-damping 和 β-damping: 2 2α
βω
ξ
ω
= + (5.2) 觀察上式可以得知阻尼比並不是一個定值,而是會隨著頻率不同以及 α-damping、β-damping 而改變。 根據文獻[25]指出,改變 α-damping 對聲壓曲線造成的最大影響是 的 部份,至於高頻部分幾乎沒什麼改變,改變β-damping 對低頻部分影響相當 的小,高頻部分才看的出其差異,因此本文取了兩個頻率來計算α-damping 和β-damping,取低頻區的第一個自然頻率且忽略β-damping 的影響來計 算α-damping,還有在高頻區 10KHz 之後找一個明顯的突起並忽略α -damping 的影響來計算β-damping,表 3.1 為本文研究之揚聲器的阻尼比。 0 f 3.3 ANSYS 簡諧激振分析與聲壓公式之應用探討 實驗方面,依照 4.1.6 節製作系統一、系統二揚聲器,分析方面,依照 上述方法依序建立揚聲器的模型如圖 3.1、圖 3.2 所示,其中振動板與加勁 用的碳纖材料常數如表 3.2 所示,由於系統一揚聲器的彈性力是由碳纖維彈 性支承以及懸邊所提供,但碳纖維彈性支承的彈性力是經由玻纖彈波傳遞 到音圈上,而音圈的直徑與振動板的寬度相當接近,所以將碳纖維彈性支承所提供的彈性力等效增加到懸邊所提供的彈性力上,因此系統一揚聲器 的模型將整體彈性力等效為懸邊元素所提供的彈性力。 將系統一、系統二揚聲器實驗所得的數據依照 3.2 節的關係式,代入模 擬所需要的參數項來進行簡諧激振分析,最後將 ANSYS 分析得到振動板的 所有點(Nodes)座標及在不同激振頻率下的振幅(Amplitude)、相位角(Phase angle)的資料存取,並套入 Fortran 程式編寫的聲壓公式計算不同頻率下之 聲壓值。將分析與實驗所得的聲壓圖比較如圖 3.3、圖 3.4 所示,可得知實 驗和分析模擬出來的兩條聲壓曲線有相同的趨勢,由此可驗證此 ANSYS 模 型是具參考價值的,往後將以此模型進行分析。
第四章 平面揚聲器製作與實驗程序 4.1 平面揚聲器製作 平面揚聲器結構可分為振動板、音圈、場磁鐵、懸邊及外框,示意圖 如圖 4.1,系統一揚聲器、系統二揚聲器、系統三揚聲器分別如圖 4.2a、圖 4.2b 與圖 4.2c 所示。 4.1.1 複合材料積層板之製作 (1) 首先由冷凍庫中取出複合材料預浸材,並在室溫下靜置 12~14 小時。 (2) 將預浸材依需要的疊層數目、角度和大小剪裁、堆疊。 (3) 在鐵板上依序放置輔助材及複合材料(如圖 4.4),並抽真空,再將鐵板 送上熱壓機(如圖 4.5),依加溫時間-溫度關係及溫度-壓力關係(如圖 4.6)來控制複材積層板成型條件。 (4) 等時間一到,關掉熱壓機的溫度,使試片在室溫下自然冷卻,待冷卻 到室溫時,再將積層板取出,即製作完成。 4.1.2 三明治振動板之製作 使用材料為碳纖預浸材與巴沙木,製作步驟如下: (1) 製作前先從冷凍庫中取出複合材料預浸材,並在室溫下回溫 12~14 小 時。 (2) 裁剪出所需預浸材之大小及角度及所需的輔助材。 (3) 依照設計好的形狀、角度、層數將碳纖預浸材置於巴沙木上。 (4) 在鐵板上依序放置輔助材及振動板半成品,並送上熱壓機(圖 4.4),依
加溫時間-溫度關係及溫度-壓力關係(圖 4.5)來控制振動板成型條件。 (5) 熱壓完成後,使振動板在室溫下自然冷卻,待溫度到達室溫後,再將 振動板取出,即製作完成。 4.1.3 音圈與場磁鐵之製作 音圈部分為設計好各尺寸高度再委託廠商製作,如圖 4.7,其中 d1為內 徑,d2為外徑,t 為厚度,tv為捲幅,h 為高度,本文所使用到的音圈尺寸 如表 4-1 所示。場磁鐵部分也是設計好各尺寸高度,委託廠商製作各零組 件,再以 AB 膠黏好成型後充磁,結構如圖 4.8,本文所使用到的激振器尺 寸如表 4-2 所示。 4.1.4 彈性懸邊之製作 分別使用不織布、矽膠、絹布三種材料,依照以下的方法壓製彈性懸 邊。 不織布熱壓形成彈性懸邊,其製作步驟如下: (1) 首先裁減下一片適當大小的不織布,於中央位置劃開 X 形狀。 (2) 將不織布放入懸邊模型如圖 4.9,以壓力 300psi、溫度 120℃熱壓 30 分鐘,放在室溫下讓其自然冷卻即完成不織布彈性懸邊。 矽膠與不織布熱壓合成彈性懸邊,其材料柔軟且富有彈性,製作步驟 如下: (2) 首先裁減下一片適當大小的鐵弗龍片,將矽膠均勻塗抹在上面。 (3) 將塗好矽膠的鐵弗龍片放上熱壓機,以 120℃的溫度預熱 10 分鐘, 將其多餘的水分揮發掉。
(4) 預熱完畢後,將其依順序蓋上適當大小的不織布與另外一片鐵弗龍 片,上下鋪好真空袋及鐵板再放上熱壓機以溫度 120℃,壓力 300psi 熱壓 20 分鐘,使矽膠穿透不織布的毛孔,均勻分佈在不織布的上下 表面。 (5) 熱壓完成之後,將矽膠與不織布合成懸邊裁剪成所需大小,放入懸邊 模型如圖 4.9 加壓成型,再放在室溫下讓其自然冷卻,完成之後如圖 4.10 所示。 絹布熱壓形成彈性懸邊,其製作步驟如下: (1) 首先裁減下一片適當大小的絹布,於中央位置劃開 X 形狀。 (2) 將絹布放入鋁制懸邊模型如圖 4.11,以壓力 300psi、溫度 120℃熱壓 30 秒鐘,放在室溫下讓其自然冷卻即絹布彈性懸邊,如圖 4.12 所示。 4.1.5 碳纖支承之製作 以碳纖維製作的長條彈性支承步驟如下: (1) 設計好支承的形狀與尺寸,於繪圖軟體繪製並列印出來。 (2) 依照 4.1.1 節的方法製作一層的碳纖維板。 (3) 將製作出的單層碳纖維板依照列印出來的長條彈性支承形狀切割成 寬度為 1.5mm 的彈性支承。 以碳纖維製作的 z 字型彈性支承步驟如下: (1) 首先由冷凍庫中取出複合材料預浸材,並在室溫下靜置 12~14 小時。 (2) 將單層碳纖預浸材上下各夾一層厚度 0.02mm 真空帶(如圖 4.13),並 放入模具壓合(如圖 4.14)。 (3) 將模具送上熱壓機(如圖 4.5),依加溫時間-溫度關係及溫度-壓力關係
(如圖 4.6)來控制 z 字型彈性支承成型條件。 (4) 等時間一到,關掉熱壓機的溫度,使模具在室溫下自然冷卻,待冷卻 到室溫時,再將 z 字型彈性支承取出(如圖 4.15)。 (5) 取出的 z 字型彈性支承以刀具(如圖 4.16)切割成 1.5mm 的寬度,再剪 出所需要的 z 字型彈性支承長度即完成。 4.1.6 揚聲器之組裝 系統一揚聲器組裝步驟如下: (1) 依序將音圈與切割好的玻纖彈波放入治具一(如圖 4.17),並將音圈與 玻纖彈波黏合。 (2) 將壓製好的寬度 1.5mm 之 z 字型碳纖維彈性支承定位在治具一的 1、 2、3、4 角落,並與玻纖彈波黏合,再將寬度 1.5mm 之碳纖維定位在 框的 5、6 位置如圖 4.18,並與 z 字型碳纖維彈性支承及框黏合完成 半成品一,如圖 4.19 所示。 (3) 將壓好的彈性懸邊裁掉多餘部份並與振動板黏合,再經由定位器精準 的與框架黏合,如圖 4.20 所示。 (4) 將音圈定位器套上激振器底座,並與半成品一結合,透過激振器底座 與框架的定位,將音圈與振動板黏合,如圖 4.21 所示。 (5) 將激振器與激振器底座黏合,並裝置在框架上黏合,即完成組裝,如 圖 4.22 所示。 系統二、系統三揚聲器組裝步驟如下: (1) 將裁好形狀的振動板放入振動板治具裡(如圖 4.23),本治具為確保振 動板的定位,使振動板邊緣與揚聲器框架保持固定間隙。 (2) 將放入振動板的治具與揚聲器框架結合,並透過音圈治具將音圈固定
在振動板上,並將音圈治具與框以熱溶膠暫時黏合,如圖 4.24 所示。 (3) 將熱壓成型的彈性懸邊裁掉多餘部份,透過治具(如圖 4.25)將懸邊分 別與振動板、框黏合(如圖 4.25),拿下音圈治具並裝上激振器,組裝 完成揚聲器(如圖 4.26)。 4.2 頻率響應實驗 頻率響應實驗以 B&K Pulse 頻譜分析儀(如圖 4.28)測定揚聲器系統之自 然頻率與阻尼係數,其原理為施予外力量測其變形,得到各頻率下振幅頻 譜,圖形中的波峰位置即為系統共振頻率。量測方式為輸入電壓激振揚聲 器,以雷射測速儀接收訊號,但此種方法所接受訊號僅為局部小面積上的 變形,因此欲知振動板整體變形,需量測不同位置,再加以平均。其實驗 程序如下: (1) 將揚聲器架設妥當,連接線路,以雷射測速儀量測對準欲量測的位 置,將訊號接受強度調整至飽和。 (2) 將 PULSE 分析軟體的環境設定完成(如:測試頻寬、解析度、激發方 式)。 (3) 啟動系統後,以掃頻的方式擷取速度振幅資料,經由頻譜分析儀計算 頻譜區域(Frequency Domain)中的頻譜。 (4) 輸出頻譜圖上峰值的振動頻率數值,計算各共振頻率阻尼係數。 (5) 將雷射測速儀對準其他量測點,並重複(2)至(4)。 (6) 平均量測頻率數值,即為此揚聲器系統的自然頻率。 4.3 聲壓相關實驗 首先使用 LMS 聲壓頻譜儀(如圖 4.29)測量揚聲器系統之系統參數,其
原理為先量測原本揚聲器的阻抗與頻率關係,將近似於揚聲器的質量的物 體另外附加於振動板中央,量測另一組阻抗與頻率,利用兩組數據比較而 得到系統的各項相關參數,表 4.3~表 4.5 為系統一~系統三揚聲器的參數。 聲壓量測實驗則是使用 LMS 聲壓測試系統(內含寬頻雜訊產生器、聲壓 頻譜分析),搭配訊號放大器(Amplifier)、個人電腦以及麥克風來進行量測, 實驗架設如圖 4.30 所示,其程序如下: (1) 將平面揚聲器固定於障板上。 (2) 麥克風放置於距離揚聲器中心一公尺之同一高度腳架上。 (3) 校正 LMS 聲壓測試系統並設定好測量相關參數。 (4) 啟動電腦發出訊號,由麥克風接收聲壓訊號,傳回電腦。即可得出頻 率響應之聲壓分貝圖。
第五章 小型平面揚聲器之設計與討論 現今的產品多講究輕、薄、短、小、不佔空間,平面揚聲器也是朝向 此方面來發展,在相同之邊界條件下,如果振動板之面積變小,則系統之 自然頻率必然會增加,導致揚聲器低頻的聲音會不明顯,所以本文藉著改 變振動板邊界條件來降低自然頻率,使揚聲器之有效頻寬往低頻增加,同 時也增加系統的穩定度,避免產生不必要的雜音。 在第一個自然頻率 之後,因為振動板的振動模態變形,造成振動板 推動的空氣有正有負互相抵消,使得聲壓會有某一段突然下降,稱之為中 音谷。傳統對振動板的加勁方式是在整面巴沙木板上下各加一片碳纖板, 但是由於加勁的面積造成增加的質量太大,進而影響整個揚聲器的效率, 本章即要探討揚聲器各參數最佳的設計,使得揚聲器有平滑的聲壓曲線與 高感度的表現。 0 f 5.1 彈性懸邊之設計 依照 4.1.4 節的方法製作出三種彈性懸邊,其中不織布熱壓而成的懸邊 效果最差,由於不織布材料的多孔性,導致懸邊無法產生氣密的效果,使 得揚聲器的音質與聲壓平均感度表現不佳,也由於不織布材料的不均勻 性,導致揚聲器的穩定度太低,容易使音圈磨到激振器,產生不必要的雜 音。 因此,改良使用矽膠與不織布熱壓合成彈性懸邊,由矽膠液體的性質, 藉以填滿不織布的孔隙,再藉由矽膠固化後的彈性增加揚聲器的穩定度與 阻尼效果,此一改良大大提升懸邊的氣密性與穩定度,揚聲器的音質與聲 壓平均感度也提升許多。 絹布是由表布與熱塑性的膠兩種材質製成,材料常數如表 5.1 所示(海
瑞電聲公司提供),由於可以更換表布以調整懸邊的剛性,又可以藉由熱塑 性的膠來固定懸邊形狀,更重要的是製程時間縮短到只需要 30 秒,所以嘗 試使用絹布來壓製彈性懸邊。表布是由互相垂直的纖維編織而成,所以均 勻性比不織布高許多,孔隙也比較密,加上熱塑性膠的塗佈使得氣密性更 足夠,實驗結果發現以絹布壓製成的懸邊更提升了揚聲器的穩定度與聲壓 平均感度。 5.2 玻纖彈波與碳纖支承之設計 傳統錐盆式喇叭單體都有彈波,用來將錐盆與音圈精準定位在框架 上,讓錐盆能夠藉音圈的控制而前後運動,並且使錐盆能在運動完畢後歸 回原位,常見的彈波截面形狀為波浪形(如圖 5.1),材質則有紙、布。傳統 截面為波浪型的彈波雖已達到降低整個振動系統彈性係數的效果,但是這 種支承在承受稍大振幅時,即呈現出明顯的非線性彈性支承力,這對揚聲 器的低頻有非常不良影響,也干涉了整體系統的穩定度。傳統波浪形彈波 由於形狀為一圓面積,佔據相當大的空間,也限制住了揚聲器的長寬比變 化,對於揚聲器朝細長形狀發展,莫過於是一大阻礙。 本文設計如圖 5.2 形狀的彈波,材質是由[紗布/玻璃纖維/紗布]熱壓而 成,系統所需要的彈性效果則由 4.1.5 節製作的碳纖支承來取代,此一設計 將傳統靠波浪形狀來提供彈性力的彈波,以碳纖支承取而代之,新的設計 如同簡支樑的效果,大大簡化傳統波浪型彈波的非線性影響;此一設計也 解決了傳統波浪形彈波對長寬比變化的限制,原本系統的彈性力是由一整 圈的波浪形彈波所提供,現在改良到只需要兩根主要的碳纖支承來提供, 不僅簡化了複雜度、減輕整體系統重量,也提升整體系統的穩定度。
5.3 系統二揚聲器改變設計參數造成的影響 由於碳纖維加勁的振動板對於聲壓沒有很明顯的幫助,而且本文所研 究的皆為小型揚聲器,因此加勁在振動板上的碳纖維重量對於聲壓曲線的 平均感度影響特別大,所以本節討論改變揚聲器系統的各項參數,以提升 聲壓曲線的平滑度與平均感度。 5.3.1 改變音圈參數造成的影響 比較系統二與系統三揚聲器的設計,此兩種設計的振動板大小、厚度 以及懸邊材料、斷面形狀均相同,唯獨兩者音圈的重量、高度、材質不同, 但是兩系統的聲壓曲線卻差異相當大;由兩系統的實驗聲壓曲線(如圖 5.3) 可以看出,系統二揚聲器在 4000Hz 的音谷聲壓值僅有 67.5dB,與 5000Hz 以後的平均聲壓值 77.5dB 有 10dB 的差距,這樣子的差距對人耳聽覺來說 是非常不舒服的,反觀系統三揚聲器聲壓曲線相當平滑,僅在 4000Hz 有約 7.5dB 的音峰,因此以下分析音圈的重量、材質對聲壓曲線的影響。 圖 5.4 為分析系統二揚聲器在改變音圈材質強度為原本的 10 倍、2 倍、 0.5 倍的聲壓曲線圖,可以發現改變音圈材質強度對於聲壓曲線並不會有幫 助的效果,因此再分析改變音圈重量對於聲壓曲線的影響,由圖 5.5 可以發 現,當逐漸增加音圈重量從 0.098g 到 0.21g 時,9000Hz 以後凸起的聲壓曲
線會逐漸趨於平緩,但 6000Hz 到 9000Hz 凸起的聲壓仍沒有改變;當增加 音圈重量到 0.28g 時, 9000Hz 以後凸起的聲壓曲線會逐漸趨於平緩,而原 本 6000Hz 到 9000Hz 凸起的聲壓往前面頻率移動,由此分析證明了改變音 圈重量的設計是可以有效平緩系統二揚聲器在高頻與中、低頻間高達 10dB 的聲壓落差的,也可以解釋改變了音圈重量的系統三揚聲器為何曲線相對 於系統二揚聲器平滑了許多。 5.3.2 改變振動板長寬比造成的影響 本節分析改變振動板的幾何比對整體聲壓的影響,為了方便模型的建 立,所以本節的振動板皆簡化以矩形來分析,與原本振動板左右兩端有半 圓面積不同。 首先在固定振動板面積為 540.89mm2下將長寬比分別設計為 2.3、1.7、 1.28,尺寸如表 5.2 所示,使用音圈半徑皆為 6.6mm,分析結果可以發現此 三種設計皆有相當深的音谷出現;根據文獻[25]指出,當音圈半徑與振動板 節線至振動板中心距離的比值越接近 1 則音谷落差越小,因此分析改變音 圈半徑以提升音谷落差,不過根據 5.3.1 節分析指出,音圈的重量變化會造 成聲壓的改變,因此先固定音圈的重量不變,再配合各長寬比設計來放大 音圈半徑,分析結果如表 5.3~表 5.5 所示,聲壓如圖 5.6~圖 5.8 所示。振 動板長寬比為 2.3 的設計,使用音圈半徑為 6.6mm 在頻率 2900Hz 的音谷有 9.747dB 的落差,當放大音圈半徑到 7.6mm 時,降低了聲壓落差到 6.502dB, 而且在 2900Hz~5833Hz 的聲壓提高了 2.69dB;振動板長寬比為 1.7 的設計, 可以發現當音圈半徑由 6.6mm 增加到 7.6mm 時,原本在頻率 3170Hz 的音
谷落差由 12.348dB 減少了 5.978dB,而原本的音谷範圍 2280Hz~7080Hz 也 因 為 在 增 加 音 圈 半 徑 後 產 生 4360Hz 的 聲 壓 凸 起 而 縮 減 為 2640Hz~4360Hz,再放大音圈半徑到 8.6mm 時,3170Hz 的音谷再減少了 5dB,而且中、低頻率的聲壓也因為音圈半徑的放大變的相當平滑;振動板 長寬比為 1.28 的設計,當音圈半徑由 6.6mm 增加到 7.6mm 時,原本 3000Hz 的音谷落差 14.906dB 減少了 4.715dB,再增加音圈半徑到 8.6mm 時,音谷 頻率往後移到 4870Hz 且落差減少 1.52dB,當音圈半徑增加到 9.2mm 時, 音谷落差再減少 3.204 dB。 5.3.3 模擬分析討論與實驗驗證 隨著振動板長寬比的改變,其橫軸彎曲模態(Bending mode)發生時所產 生的節線位置也會有所不同,本節以振動板中心為原點,將各種不同長寬 比其造成音谷之模態的節線位置整理成表 5.6 並繪成圖 5.9,可以明顯看出, 當長寬比 x 越大,節線距中心原點的距離 y 越大,以多項式的形式可找出 近似曲線: y =3.1824 x3 −16.568x2 +29.671x−6.9231 (5.1) 上述公式與振動板面積相關,因此當振動板面積 A=540.89mm2的狀況 之下才成立,將來可以用式(5.1)來預測節線的位置,進而決定該選用何種 直徑的音圈,藉以施力在節線的方式減小彎曲模態(Bending mode)對聲壓的 影響,使中音谷落差減小。 由圖 5.10~圖 5.12 可以看出來當音圈半徑與節線至中心距離的比越接 近 1 時音谷落差越小,但在面積固定之下,隨著長寬比越大,振動板的寬 度會越小,若是想使用接近節線位置的音圈,其直徑會超出板子寬度,因
此選擇適當的長寬比才能有效利用音圈的施力位置改善聲壓曲線,以下將 討論長寬比在何種範圍之內,音圈的半徑與節線位置之比例能夠接近 1,且 直徑不超過振動板寬度。 假設音圈的直徑與振動板的寬度 a 相等,y 為節線至中心的距離,則目 標為 2 2 8 . 0 a ≤ y≤ a => 0.4≤ ≤0.5 a y (5.2) 將式(5.1)無因次化之後可以得到以下式子: 3832 . 0 2483 . 1 5759 . 0 1123 . 0 3 − 2 + − = x x x a y (5.3) 將其繪成圖 5.13 可以得知,若是想達到式(5.2)的條件,則長寬比 r 必須介 於 1~1.28 之間。 根據以上分析結果,以長寬比為 1.28 並搭配音圈半徑為 9.2mm 的設計 來製作揚聲器,不過由於實驗材料音圈不容易取得的限制,因此以現有音 圈半徑為 6.6mm 為主,同步縮小振動板面積,但縮小後的長、寬與音圈半 徑彼此之間的比例仍然維持不變,以此縮小比例建模分析,實驗部份依照 4.1 節方法製作系統四揚聲器(如圖 5.14 所示),分析與實驗的聲壓曲線比較 如圖 5.15 所示,分析與實驗曲線趨勢相近,足以證明分析的正確性。 將系統一揚聲器等比例放大 1.2 倍且定義為系統五揚聲器,並將長寬比 皆為 1.28 的系統一揚聲器、系統四揚聲器、系統五揚聲器分析曲線比較如 圖 5.16 所示,可以發現三條曲線的趨勢很相近,因此可以驗證本文 5.3 式 無因次化後公式的正確性,長寬比 1.28 是較佳的設計參數,而且不同面積
大小的揚聲器只要維持本文的長寬比與音圈半徑的比例,整體的聲壓趨勢 是相似的,只是因為懸邊提供的彈性力與振動板面積大小的不同,所以導 致某些相同的振形出現在不同的頻率下。如圖 5.17 所示將三條分析曲線出 現相同振形的頻率標以 、 、 圖形分類註記,其中各圖形標記所代表的 振形如圖 5.18~圖 5.20 所示,可以發現當圖形標記 系列的振形出現時,聲 壓曲線有提升的效果,觀察振形圖可以發現板子的變形均朝同一方向,因 此產生的聲壓在此振形下不會互相抵消,所以聲壓會逐漸增加;當圖形標 記 系列的振形出現時,聲壓曲線由相對高點開始下降,觀察振形圖可以 發現板子一開始的變形均朝同一方向,但是隨著頻率的增加而開始有更多 的面積朝反向變形,因此兩種不同方向板變形所貢獻的聲壓會互相抵消, 所以聲壓就會呈現下降的趨勢;當圖形標記 系列的振形出現時,聲壓曲 線由相對低點開始上升,觀察振形圖可以發現板子的變形大部分面積均朝 同一方向,僅四邊角落朝反向變形,隨著頻率的增加而在四邊角落的反向 變形面積逐漸增加,而板中央區域也出現反向變形面積,因此聲壓曲線由 相對低點開始上升。 比較長寬比 1.28 系統四揚聲器與系統二、系統三揚聲器的實驗聲壓曲 線(如圖 5.21),長寬比 1.28 的設計在 554Hz 以後的聲壓開始緩降,直到 3590Hz 再度上升,期間的聲壓落差只有 5.6dB,因此音谷的影響遠較於其 他兩條曲線小了許多,在 554Hz~3590Hz 這段頻率間的聲壓也比其他兩條曲