Microsoft Word - thesis.doc

Transcription

1 中華大學 碩士論文 題目 : 活體中耳聽小骨鏈動態特性分析 The vibration analysis of in-vivo human ossicular chain 系所別 : 機械與航太工程研究所 學號姓名 :M 李文賢 指導教授 : 陳 精 一 博士 余 仁 方 博士 中華民國九十六年八月

2 中文摘要 人耳主要分為外耳 中耳 內耳三個主要系統, 然而中耳的功能主要為聲音傳遞之用, 中耳腔內共有三根聽小骨, 依次為鎚骨 (malleus) 鉆骨(incus) 鐙骨(stapes) 當鼓膜(tympanic membrane) 受到刺激後, 聲音就經由聽小骨鏈將聲音擴大傳至內耳基底膜 (basilar membrane), 進而傳至大腦產生聽覺 在臨床醫學方面, 對於中耳病變患者術後聽小骨鏈及鼓膜重建, 因鼓膜 聽小骨鏈之幾何外形皆以改變, 因而可藉由有限元素模型重建加以探討 本研究採非侵入分析方式, 以有限元素模型重建探討活體中耳聽小骨鏈振動動態分析 首先利用電腦斷層掃描 (computed tomography, CT) 影像, 透過 Amira 建立其立體影像模組, 以便獲取聽小骨鏈之幾何外型尺寸, 再經由 ANSS 建立其有限元素模組 有限元素分析流程分成兩個主要步驟, 一為模態分析探討個別聽小骨鏈及其兩兩偶合後自然頻率的改變 另一則為調和力分析, 施予 90dB 聲壓值去刺激鎚骨柄部 (manubrium) 與鼓膜連接處, 檢視 25/250/500/750/k/2k/3k/4k/5k/6/7k/8k Hz 各頻率下其頻率響應, 以利探討聽小骨鏈的動態特性 關鍵詞 : 聽小骨鏈 振動分析 有限元素法 電腦斷層掃描圖 I

3 Abstract Human ears consist of three main parts including out ear, middle ear, and inner ear. Malleus, incus, stapes are involved in the middle ear cavity. When the tympanic membrane is excited by the sound, the ossicular chain will transmit the sound energy into the inner ear through the oval window. In clinic, the geometry of the ossicular chain was changed in postoperation for the patient with otosclerosis. The 3D geometric model of the ossicular chain can be reconstructed by computered tomography (CT) image. The finite element model of the ossicular chain of the patient could then be built by 3D geometric model of the ossicular chain. The finite element analysis consisted of two main phases. The first phase was to discuss the natural frequencies of malleus, incus and stapes. The second phase was to discuss the amplitude at the nodes of the footplate at the frequencies of 25, 250, 500, 750, k, 2k, 3k, 4k, 5k, 6k, 7k, 8k Hz, and sound pressure (90 db) applied at the bottom area of manubrium. Hence, the dynamic characteristic of the in-vivo human ossicular chain could be analyzed. Keywords: ossicular chain, vibration analysis, finite element method, computed tomography image II

4 誌謝 時光飛逝, 短短兩年的碩士班生活即將告一段落, 首先我最感謝我的指導教授陳精一博士, 也算是我在學業上的啟蒙老師, 兩年來的諄諄教誨, 無論在學業 待人處世方面等都深深影響著我 接著要感謝的是默默引導我從機械進入生物醫學領域的余仁方博士, 博士帶領對於生物醫學領域懵懂無知的我, 至今可以針對臨床需求, 應用所學對於發現的現象做一個完整的討論 此外, 還要感謝口試委員倪慶羽博士, 針對論文提出許多指導與建議, 使得本論文更加完善 另外, 我仍想感謝常常回實驗室關心的學長耀祥 俊諺 閔雄, 在研究生涯中給予許多經驗上的交流, 也同時給予我精神的的支持 還有, 陪我一路走來的同學 - 振忠, 每每總是在夜深人靜的夜晚一起討論, 才有至今深厚的基礎, 同時也解決不少實驗室對於 ANSS 後處理繁瑣的工作 此外, 當然不可少的是實驗室的活力泉源一宏 自豪 子翔 彥達 祥維, 沒有你們的實驗室就呈現一片暗沉, 在這一年來謝謝你們! 最後, 我最感謝能夠體諒我這兩年沒日沒夜的女朋友 - 雅惠, 謝謝你的支持讓我順利的完成我的碩士生活 III

5 目錄 中文摘要... I 英文摘要... II 誌謝...III 目錄. IV 圖目錄... V 表目錄...VI 第一章緒論... - 前言.. -2 研究動機與目的. 4-3 研究方法. 6-4 章節概要. 7 第二章文獻回顧 8 第三章有限元素分析 電腦影像前置作業 有限元素分析流程 前處理 分析方式 第四章結果與討論 IV

6 4- 模態分析結果 探討模態分析結果在臨床上的意義 調和力分析結果 調合力分析結果與文獻資料比對.. 4 第五章結論...47 V

7 圖目錄 圖 - 人耳的構造.2 圖 3- 影像處理流程...3 圖 3-2 右側顱部電腦斷層影像 (a) 鎚骨與 (b) 鉆骨.. 5 圖 3-3 重建 (a) 鎚骨 (b) 鉆骨和 (c) 鐙骨之立體模型圖...6 圖 3-4 重建 (a) 鎚骨 (b) 鉆骨與 (c) 鐙骨實體模型.6 圖 3-5 有限元素分析流程.. 7 圖 3-6 聽小骨鏈之有限元素模型.. 8 圖 3-7 鐙骨足板彈簧元素 (a) 垂直方向 (b) 水平方向. 23 圖 3-8 鐙骨足板運動模式 圖 3-9 鐙骨足板上之束制條件 圖 3-0 調合力分析負載型式 圖 4- 鎚骨 25/250/500/750/k/2k Hz 模態振形 圖 4-2 鎚骨 3k/4k/5k/6k/7k/8k Hz 模態振形 圖 4-3 鉆骨 25/250/500/750/k/2k Hz 模態振形 圖 4-4 鉆骨 3k/4k/5k/6k/7k/8k Hz 模態振形. 3 圖 4-5 鐙骨 25/250/500/750/k/2k Hz 模態振形 圖 4-6 鐙骨 3k/4k/5k/6k/7k/8k Hz 模態振形 圖 4-7 鉆骨與鎚骨 25/250/500/750/k/2k Hz 模態振形... 3 VI

8 圖 4-8 鉆骨與鎚骨 3k/4k/5k/6k/7k/8k Hz 模態振形 圖 4-9 鉆骨與鐙骨 25/250/500/750/k/2k Hz 模態振形. 33 圖 4-0 鉆骨與鐙骨 3k/4k/5k/6k/7k/8k Hz 模態振形 圖 4- 聽小骨鏈 25/250/500/750/k/2k Hz 模態振形...35 圖 4-2 聽小骨鏈 3k/4k/5k/6k/7k/8k Hz 模態振形 圖 4-3 鐙骨足板上節點及對應節點編號 圖 4-4 節點 ~3 號頻率響應圖 圖 4-5 節點 4~6 號頻率響應圖 圖 4-6 節點 7~9 號頻率響應圖 VII

9 表目錄 表 3- 聽骨材料參數.. 20 表 3-2 聽小骨各部位體積與元素數目.. 22 表 4- 六種形式聽骨組合結構與臨床測試頻率模態分析結果 表 4-2 幾何外形尺寸比較表... 4 表 4-3 材料參數比較表 (I). 42 表 4-4 材料參數比較表 (II) 表 4-5 邊界條件比較表 (I) 表 4-6 邊界條件比較表 (II) VIII

10 第一章緒論 - 前言耳朵為人體主司聽覺與身體平衡的器官, 由外而內分別為外耳 (out ear)- 從耳廓至鼓膜 ; 中耳 (middle ear)- 從鎚骨柄部延伸至鐙骨足板 (stapedial footplate); 內耳 (inner ear)- 從卵圓窗 (oval window) 至前庭系統與耳蝸迷路 空氣分子的振盪產生聲波, 經由耳道刺激鼓膜振動, 再依序由聽小骨 - 鎚骨 (malleus) 鉆骨(incus) 鐙骨(stapes) - 構成的聽小骨鏈藉由阻抗匹配效應擴大聲波振幅, 進而擾動內耳基底膜, 藉由淋巴液的流動刺激耳蝸的末梢感受器, 由耳蝸毛細胞 (hair cell) 產生電流訊號經由耳蝸神經傳送至大腦產生聽覺 中耳腔又稱為中耳空腔, 是一個充滿空氣的密閉空間, 唯一可以進出的管道是位於在中耳腔前方的耳咽管 (eustachian tube, ET), 連接中耳腔至鼻咽的側方, 如圖 - 所示 聲音自空氣傳入外耳道經中耳聽小骨傳遞至內耳之內淋巴液時, 約有 99.9% 之能量被反射而喪失, 只剩下 0.% 進入內耳, 相當於損失 30 分貝 因此, 中耳必須藉由下列三種方式加強聲音之能量 :( 一 ) 鼓膜與卵圓窗之面積比 (area ratio) 為 20:, 有效的鼓膜振動面積與卵圓窗之面積比為 4:, 可增加約 22 分貝 ( 二 ) 聽小骨之槓桿比 (lever ratio), 即鎚骨柄

11 (handle) 與鉆骨長突 (long process) 的長度比為.3:, 可增加約 2 分貝 ( 三 ) 鼓膜錐形效果 (conical effect) 可增加約 6 分貝傳輸效果 藉由以上 3 種方式, 以及外耳道之共鳴, 可補回 30 分貝由空氣 ( 外耳 ) 傳入內淋巴液 ( 內耳 ) 聲音能量之損失 圖 - 人耳的構造 中耳引起的聽力障礙主要原因為慢性中耳炎 (chronic otitis media) 漿液性中耳炎 (otitis media with effusion) 與耳硬化症 (otosclerosis) 慢性中耳炎主要發生於成人, 經由局部檢查可發現鼓膜破損 中 耳流膿 肉芽組織生成 而鼓膜破損 聽小骨因發炎而侵蝕 肉芽 組織的阻礙都會影響聲波的傳遞而導致聽力受損 除了中耳流膿可 2

12 以使用藥物治療外, 開刀清除發炎組織 修補鼓膜是治療慢性中耳炎的主要方法 但若開刀後聽力恢復不理想, 仍須靠助聽器彌補其聽力 漿液性中耳炎主要發生在小孩, 少部分在大人 大部分的病因是耳咽管阻塞, 使滲出液貯留在中耳腔中無法排除, 因而阻礙中耳聲波傳遞 小孩除了聽力減弱之外, 甚少有其他症狀, 所以極易忽略 治療早期以藥物為主, 但若長期積液不退或嚴重影響聽力, 則須手術抽除積液, 甚至植入中耳通氣管 值得注意的是, 若大人發生漿液性中耳炎, 首先要排除鼻咽癌的可能 耳硬化症主要的影響, 聲波於中耳是經由三塊小骨頭的振動而傳遞, 所以當骨頭之間的關節硬化無法振動時, 聲波就不能繼續傳遞了, 這種硬化主要發生在鐙骨和耳蝸之間, 稱之為耳硬化症 常發現於 20 歲到 40 歲之間的女性, 症狀包括漸進性聽力喪失, 有時伴隨耳鳴, 症狀可能因懷孕而加重 治療以手術為主, 若無法手術者可應用助聽器改善 臨床上傳導性聽損 (conductive hearing loss) 與聽小骨硬化症患 者, 以各種不同的外科手術補救聽力, 皆可藉由建立聽小骨鏈及其 鼓膜之有限元素模組, 作為術前替代物 (prosthesis) 植入的評估 3

13 -2 研究動機與目的在臨床上, 慢性化膿性中耳炎 (chronic suppurative otitis media) 最常見且最為惱人的症狀除了耳漏及耳膜破洞, 並且中耳的聽小骨也可能會受到波及, 造成聽小骨的破壞, 而中耳先天畸型或其他原因引起的聽小骨斷裂, 為了儘可能恢復聽力時, 聽小骨成形術 (ossiculoplasty)[] 是唯一的選擇 然而, 術後則發現聽小骨傳導力下降, 而大氣壓力所造成的正負壓亦會影響聽小骨鏈 (ossicular chain) 的振動 [2-3], 造成傳導性聽力受損, 為了解決此問題, 電路模型 (generalized circuit model)[4] 被提出用以模擬聽小骨空間振動, 或是用 malleus vibration audiometer (MVA)[5] 來評估, 但是無法有效描述聽小骨中較重的骨頭部位如 : 鎚骨頭 (malleus head) 砧骨本體(incus body), 另外使用有限元素法 [5-6] 分析聽小骨振動模態, 其中聽小骨之有限元素模型係以人體組織學圖片 (histological section images)[5] 及由解剖所取得之聽小骨的 CT 影像 [6] 建立而成的 除此之外, 中耳空腔 (middle-ear air spaces)[6-8] 也由動物實驗中發現, 會影響低頻振動, 因此, 大氣壓力所造成的正負壓以及中耳的空腔, 對活體人耳內之聽小骨鏈所造成的影響, 則是值得更進一步的探討 4

14 因此本研究運用臨床所獲得之活體聽小骨電腦斷層影像, 經過影像轉換及重建之技術, 呈現活體聽小骨之 3D 立體影像, 由此便可獲得其幾何外型, 進而建立有限元素模組, 並探討活體聽小骨鏈之模態響應與調和外力振動特性, 希望可以提供臨床醫師針對需植入聽小骨義肢患者術前植入物的評估 5

15 -3 研究方法本研究將針對患者電腦斷層掃描影像, 擷取出聽小骨鏈之幾何影像邊界, 再經由 Amira 醫學影像視覺化軟體重建聽小骨鏈立體影像模型, 再由 CAD 軟體將 STL 轉換輸出為 SAT 格式之圖檔, 最後於 ANSS 中採用自由網格方式, 建立出各聽小骨有限元素模型 有限元素分析方面, 分為模態分析與調和力分析兩部份 模態分析方面主要是探討聽骨與聽骨間偶合效應之動態特性, 了解聽小骨結構各頻段的自然頻率, 提供需進行鐙骨手術 聽小骨義肢患者術前 / 術後的評估 ; 調和力分析則探討鎚骨柄部與鼓膜連接處施予 90dB 聲壓負載下, 鐙骨足板之頻率響應 6

16 -4 章節概要本論文共分為五章, 第一章緒論, 包含前言 研究動機與目的 研究方法 ; 第二章為文獻回顧, 簡述歷年相關文獻成果 ; 第三章為有限元素分析, 分為影像前置作業 有限元素分析流程兩大部分 ; 第四章為結果與討論, 模態分析與調和力分析並與文獻中調和力分析結果相比較 ; 第五章為結論, 總結模態與調和力分析結果, 並詳述於臨床術前 / 後之意義為何 7

17 第二章文獻回顧 現今有限元素法應用於生物力學領域已經相當普遍, 也隨著醫學影像解析度與電腦計算能力大幅提升, 從人體最大的骨骼 - 股骨與最小的骨骼 - 聽小骨, 皆可透過影像擷取的方式, 如組織學切片 電腦斷層影像 核磁共振造影 等技術, 取得活體立體影像 以下, 將依照年份的演進介紹有限元素法應用於聽小骨鏈分析的文獻 992 年 Wada 與 Kobayash[9] 利用屍體 聽損患者與正常聽力三種測試人員, 提出鼓膜周圍與鐙骨足板上之環狀韌帶與為彈性邊界並提出由鎚骨前端延伸至鉆骨短突部有一假設旋轉軸 2k 年,Ferris 與 Prendergast[0] 發現聽小骨部份置換手術 (partial ossicular replacement prosthesis, PROP) 與聽小骨全置換手術 (total ossicular replacement prosthesis, TROP), 術後植入物會改變其原有的動態特性, 使得聽骨鏈有不同於原本正常結構的自然頻率 並提出植入物置換手術必須針對聽骨鏈的剛性結構設計, 使得術後聽骨鏈與原本結構自然頻率相近 2002 年 Wada 等人 [] 利用有限元素法重建中耳聽小骨鏈, 其結 構包括韌帶 (ligament) 腱 (tendons) 鉆骨與鐙骨間接點 (I-S joint) 8

18 耳蝸阻尼 外耳道 (external auditory meatus, EAM) 中耳空腔等 並假設耳蝸為一阻尼結構, 加上韌帶與腱將原本複雜的聽骨鏈運動簡化為彈簧阻尼結構 最後, 發現中耳空腔並不影響鼓膜的振動模態 然而, 聲壓經過外耳道作用於鼓膜時, 耳道入口聲壓與鼓膜表面上聲壓相差無幾, 且分布於鼓膜表面聲壓卻不受外耳道形狀 長度影響 同年,Sun 等人 [2] 利用人體顳骨組織學切片, 將聽小骨外輪廓按切片影像一一分離為各層點資料, 並使用 B-spline 方法找出每一層最擬合點資料的平滑曲線, 再經由電腦輔助設計軟體將平滑曲線拉伸為一立體體積 此法不僅僅可有效提供聽小骨鏈立體影像及量測三維幾何尺寸, 更可將其導入有限元素分析, 進一步了解其動態特性 同年,Sun 等人 [3] 將屍體的顳骨切片影像, 利用電腦輔助設計環境整合, 重建耳膜及聽小骨鏈立體影像, 導入有限元素法模擬其結構受音壓負載後的分析結果與都普勒雷射干涉儀量測鐙骨足板上的微小位移相比較, 證實電腦輔助設計整合環境便利與強大之處 2003 年,Kelly 等人 [4] 利用電腦斷層掃描影像, 重建聽小骨鏈之有限元素模型, 並使用 Kurz omed 兩種分別針對 PROP 與 TROP 組合測試, 分析發現 Kurz TROP 較 omed TROP 與正常鐙骨足板振幅相近, 但置換後之鐙骨足板, 仍較原本結構較易產生共振 9

19 現象 但已證實可以利用電腦模擬置換後聽骨鏈剛性結構, 並可預 測其振動動態特性 2004 年,Gan 等人 [5] 利用人體組織學切片影像, 利用 CAD 整合環境重建左側人耳外耳道 聽小骨 鼓膜 韌帶 肌肉 中耳空腔等結構, 最後進行有限元素分析, 發現鼓膜與聽骨鏈共振頻率約於 3kHz 至 4kHz 間 2005 年,Kokie 等人 [6] 利用有限元素法針對聽小骨鏈分析結果顯示, 鐙骨足板上環狀韌帶與鉆骨 鐙骨間接點的楊氏係數對於預測聽骨鏈動態特性有著較大的影響力, 因而欲使用有限元素分析聽骨鏈動態特性, 這兩部分將不可或缺 2006 年,Gan 等人 [7] 利用有限元素分析中耳聽骨鏈 外耳道 中耳空腔 等, 壓力分布於中耳空腔與外聽道, 發現聲壓負載無論於外耳道何處在 3.5k Hz 以下時, 皆無明顯變化 此外, 發現中耳穿孔的位置與形式會大大影響鼓膜對於頻率的敏感性 2006 年,Lee 等人 [8] 利用電腦斷層掃描影像重建聽小骨鏈之有 限元素模型, 聽小骨鏈幾何尺寸與人體組織學切片所量取的資料及 文獻皆相差無幾 鐙骨足板與鼓膜凸位移的預測皆與 Nishihara, 0

20 Huber, Gan, and Sun 的實驗資料相近 因而, 可利用電腦斷層掃描影像針對活體聽小骨鏈, 重建其有限元素模型並可預測其動態特性, 在在降低以往於聽小骨動態特性未知的情況下所耗費龐大的醫療成本 綜觀以上文獻, 針對中耳聽小骨鏈由影像資料, 經由電腦輔助設計環境整合, 再導入有限元素分析法進行聽小骨鏈之動態特性分析, 已經是一種較低成本且發展成熟的技術 因有限元素法的導入, 對於了解聽骨鏈之動態特性有極大的助益, 也因此降低龐大的醫療成本 但不難發現針對聽小骨鏈的材料特性與邊界條件, 許多文獻中仍停留在假設階段, 必須由實驗的資料加以佐證, 因而有限元素法必須與實驗數據相輔相成, 並無法完全取代實驗資料

21 第三章有限元素分析 本研究利用有限元素軟體 ANSS 針對活體人耳聽小骨, 進行聽小骨模態分析以及調合力分析 模態分析主要是針對個別聽小骨及兩兩偶合後自然頻率的改變, 可供聽小骨硬化症患者術前植入物評估, 避免術後聽小骨義肢易與外來刺激共振 調合力分析主要是探討給予鎚骨柄部聲壓 (90dB-SPL) 下, 鐙骨足板的頻率響應 其中, 調和力分析中以音壓 (sound press level, SPL) 模擬人耳聽到聲音的強度, 因人耳聽覺聲音的響應範圍相當廣大, 因而以 分貝 (decibel, db) 值 描述 然而, 分貝並非描述聲音大小的單位, 而是兩個音壓的相對強渡 其中, 以人耳平均的聽覺下限 P ref = dyne/cm 2 為基準音壓 因此, 聲音的強度以分貝表示, 如下式所示 : p SPL( db) = 20log p r ef () 本章將針對電腦影像前置作業以及有限元素分析兩部份, 並加 以說明其建立流程 2

22 3- 電腦影像前置作業 由長庚醫院取得病患顱部電腦斷層影像後, 先將 CT 影像匯入 Amira 後, 由 2-D 截面重建 3-D 立體模型, 再經由 SolidWorks 重建 實體幾何模型, 最後完成活體聽小骨影像前置作業, 其流程如圖 3- 所示 患者 CT 影像 Output *.STL Files Amira 重建立體影像 SolidWorks 重建實體幾何模型 NO Output *.SAT Files ES ANSS 建立有限元素模型 圖 3- 影像處理流程 3

23 圖 3-2 為病患右側顱部電腦斷層影像, 其中, 箭頭所指處為鎚骨和鉆骨 在影像重建方面, 先由電腦斷層影像找出每一張圖檔中鎚骨 鉆骨 及鐙骨影像邊界 利用 Amira 影像工具, 將中耳聽小骨影像一一分離, 再由 SurfaceGen 計算工具將其由所圈選之 2-D 截面推演成三角形構成之 3-D 立體影像, 如圖 3-3 所示, 箭頭所指處為聽小骨 (a) 鎚骨 (b) 鉆骨 (c) 鐙骨之空間模型 最後將各別聽小骨立體影像, 依照先前區分出的三根聽小骨, 分別輸出為 malleus.stl incus.stl stapes.stl 三個圖檔, 藉由 SolidWorks 匯入圖檔後, 可快速重建實體幾何模型 完成後即可分別輸出為 malleus.sat incus.sat stapes.sat 三種格式之圖檔, 如圖 3-4 所示, 即完成影像前置作業工作 4

24 (b) (a) 圖 3-2 右側顱部電腦斷層影像 (a) 鎚骨與 (b) 鉆骨 5

25 (b) (a) (c) 圖 3-3 重建 (a) 鎚骨 (b) 鉆骨和 (c) 鐙骨之立體模型圖 (a) (b) (c) 圖 3-4 重建 (a) 鎚骨 (b) 鉆骨與 (c) 鐙骨實體模型 6

26 3-2 有限元素分析流程 本研究採用 ANSS 有限元素分析軟體建立其有限元素模型, 模型建立分為三大步驟, 依序分別為前處理 求解以及後處理, 其 流程如圖 3-5 所示 蒐集資料 工程問題 提出改善方式 有限元素模型 幾何外型 材料性質 元素生成 邊界條件 負載條件 解題方程 前處理 求解 NO 顯示結果 後處理 ES 最佳合理解釋 圖 3-5 有限元素分析流程 7

27 3-2- 前處理 () 匯入實體模型依序匯入 *.SAT 格式檔案, 分別產生聽小骨實體體積, 並依聽小骨密度不同而劃分出不同部位, 如鎚骨 - 頭部 頸部 柄部, 鉆骨 - 軀幹部 長突部 短突部, 鐙骨, 鉆骨與鐙骨接點 (incudostapedial joint, I-S joint) 等, 如圖 3-6 所示 鉆骨 短突部 軀幹部 長突部鐙骨鐙骨鉆骨 / 鐙骨接點 鎚骨頭部頸部柄部 圖 3-6 聽小骨鏈之有限元素模型 8

28 (2) 材料參數活體人耳聽小骨有限元素模型建立, 採用 Tonndorf and Khanna [9] 當聲音 (2dB) 刺激人耳時, 耳膜最大應變小於 0-4 mm 與 Sun[2] 所提出的假設當聲音 (90dB) 從耳膜經由聽小骨傳至耳蝸, 耳膜最大應變值不超過 0-4 mm 因此, 本研究中聽小骨鏈及其韌帶 (ligament) 可視為線性微振動系統 聽小骨鏈各部位材料參數皆不盡相同, 唯整體結構之卜松比假設為 0.3, 因目前所發表文獻都非常接近此值 並且,Funnell and Laszlo[20] 提出卜松比對於分析聽小骨動態特性中, 並無顯著影響 阻尼之元素矩陣如下所示 [ c] = α[ M] + β[ K] (2) 其中 [M] [K] 分別為系統之質量矩陣與剛性矩陣,α β 則是系統 之阻尼參數 其整體結構則包括三根聽小骨 鉆骨與鐙骨接點以及 鐙骨上的彈簧元素所構成 聽小骨鏈之楊氏係數皆為 4.GPa, 並假設其為線彈性 均質 均向材料 聽小骨各部位密度不盡相同, 如表 3- 所示 9

29 鉆骨與鎚骨接點 (incudomalleolar joint, I-M joint) 採用 Vlaming[2] 所提出低頻 (<3k Hz) 時, 鎚骨與鉆骨並無相對運動 因而, 其楊氏係數假設為 4.GPa, 本研究中也將鎚骨與鉆骨相連接之節點結合, 視為同一節點 鉆骨與鐙骨接點採用 Hüttenbrink[22] 提出在於噪音與壓力遽變情況下, 此接點處並無發生剛體位移, 並且發現鉆骨與鐙骨有些許相對運動發生, 以此保護內耳不因鐙骨足板位移過大而受損 此外, 並參照 Ferris and Prendergast[0] 與 Wada 等人 [23] 假設鉆骨與鐙骨接點為均質 均向材料, 楊氏係數 0.6MPa Malleus Incus Stapes Stapes 表 3- 聽骨材料參數 Structure Density(kg/m 3 ) oung modulus(n/m 2 ) Head 2.55*0 3 Neck 4.53*0 3.4*0 Handle 3.70*0 3 Body 2.36*0 3 Short process 2.26*0 3.4*0 Long process 5.08* *0 3.4*0 0 I-S joint.2*0 3 6*0 5 Structure Damping(N*s/m) Spring constant(n/m) Stapedius annular ligament Cochlear fluid Proportional damping α:0 s - β:0.000s

30 (3) 有限元素模型本研究中因探討聽小骨動態特性為首要目標, 因而採用 SOLID85 4-Node 之四面體元素型式, 以利於聽小骨鏈不規則幾何外型建立, 並針對其結構微小應變反應及其計算時間方面是為較佳選擇, 聽小骨各部位體積與元素數目, 如表 3-2 所示 其中, 兩兩相連之聽小骨是將其相緊鄰之元素表面所共用節點壓縮為一個節點 ; 鐙骨足板上則考慮以 COMBIN4 Spring- Damper 彈簧元素使用於垂直與水平面兩方向 其垂直方向之彈簧元素為模擬聲音由中耳傳達到內耳之耳蝸流體的阻尼力, 如圖 3-7(a) 所示 ; 水平方向則模擬鐙骨足板上環狀韌帶, 考慮鐙骨足板來回擺動的運動方式, 如圖 3-7(b) 所示 (4) 設定邊界條件模態分析方面, 採用自由邊界, 分析對象僅以個別聽小骨與兩兩相連之聽小骨以及聽小骨鏈六種形式, 採用 free-free 雙自由端的邊界條件 調和力分析方面, 鐙骨足板垂直面上有 6 組 x-y 平面彈簧元素, 總彈簧常數 (K) 為 60 N/m, 總阻尼為 N*s/m; 水平面上有 9 組 3D 縱向彈簧元素則假設彈簧常數為 9 N/m 分別固定鐙骨足板上垂 2

31 直與水平兩方向之彈簧元素與聽小骨鏈相連節點的另一端 U U U 三方向位移 並且, 為了模擬聲音由聽小骨傳入內耳耳蝸流體的反應, 將其垂直方向限制為僅可縱向拉伸之 2-D 彈簧阻尼元素, 水平方向則是用以模擬鐙骨足板上實際運動, 運動方式如圖 3-8 所示, 因而假設為 3-D 彈簧阻尼元素 表 3-2 聽小骨各部位體積與元素數目 體積 (mm 3 ) 元素數目 頭部 鎚骨 頸部 柄部 短突部 鉆骨 軀幹部 長突部 鐙骨 鐙骨本體 鉆骨與鐙骨接點

32 2 (a) 2 (b) 圖 3-7 鐙骨足板彈簧元素 (a) 垂直方向 (b) 水平方向 圖 3-8 鐙骨足板運動模式 [24] 23

33 3-2-2 分析方式 本文採用模態 (modal) 與調和力 (harmonic excited) 兩種型式分析, 以下將詳細介紹兩種分析流程 : () 模態分析 : 本次分析針對個別聽小骨 聽小骨兩兩耦合以及聽小骨鏈共六種聽骨結構組合為分析對象, 主要是針對須植入聽骨義肢患者, 提供術前 / 術後評估 其中, 聽骨兩兩偶合後, 相鄰聽骨表面上將有兩個節點, 分別來自不同聽骨, 共用邊界上之雙重節點皆應用 nummrg 技巧將以合併為單一節點, 並將其視為同一結構 先以 Block Lanczos 方法, 邊界條件皆為自由邊界 (free-free), 求取其自然頻率與模態向量, 進而將系統轉換至模態座標, 求取模態解, 目的在於了解個別聽小骨及其兩兩偶合後以及聽小骨鏈結構之自然頻率以及振動模態 (2) 調合外力分析 : 本次分析對象為聽小骨鏈及鐙骨足板上環狀韌帶與阻尼彈簧元素, 並固定彈簧阻尼元素非與鐙骨足板相連接之節點三個方向的自由度, 如圖 3-9 所示 並施予鎚骨底部面積上之節點負載為 - 方向 90 db-spl, 如圖 3-0 所示 再利用模態疊加 (modal superposition, 24

34 MSUP) 理論, 轉回原座標系以求得模態疊加結果, 最後再將模態疊 加結果擴展至全結構, 目的在於了解聽骨受到聲壓負載後, 聲音由 聽小骨傳遞至內耳的反應 圖 3-9 鐙骨足板上之束制條件 25

35 圖 3-0 調合力分析負載型式 26

36 第四章結果與討論 目前已由長庚醫院取得一組頭顱電腦斷層影像, 影像前置作業已重建聽小骨 - 鎚骨 鉆骨 鐙骨 - 構成立體影像 利用 ANSS 針對聽小骨鏈進行模態分析與調和力分析, 主要是要了解其結構之自然頻率是否易與外來刺激音產生共振現象, 而臨床醫學方面則是以 25/ 250/ 500/ k/ 2k/ 4k/ 8k Hz 七個主要檢測頻率 因而, 利用 ANSS 以掃頻 (sweep-frequency) 方式求解, 主要將所重視檢測頻段加以細分為 2 個頻率作為結構物是否易與外來刺激音共振的指標 4- 模態分析結果 模態分析針對個別聽小骨與聽小骨兩兩耦合以及聽小骨鏈總共六種組合, 藉此了解個別聽小骨及其兩兩耦合後自然頻率的改變 因採用 free-free 自由邊界條件, 前六個模態頻率為剛體運動 (rigid body mode), 故頻率幾乎為零 模態分析結果求解其結構自然頻率接近於 25/ 250/ 500/ 750/ k/ 2k/ 3k/ 4k/ 5k/ 6k/ 7k/ 8k 這 2 個頻率, 聽小骨的振動模態響應及振 形 以下將六種形式聽骨結構最接近上述 2 個頻率模態振形一一呈 27

37 現 圖 4- 鎚骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形, 圖 4-2 鎚骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形, 圖 4-3 鉆骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形, 圖 4-4 鉆骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形, 圖 4-5 鐙骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形, 圖 4-6 鐙骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形, 圖 4-7 鉆骨與鎚骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形, 圖 4-8 鉆骨與鎚骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形, 圖 4-9 鉆骨與鐙骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形, 圖 4-0 鉆骨與鐙骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形, 圖 4- 聽小骨鏈 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形, 圖 4-2 聽小骨鏈 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 28

38 SUB =2 FREQ=8.959 DM =2.678 SUB =8 FREQ= DM =8.82 SUB =39 FREQ= DM =9.579 (a) (b) (c) SUB =73 FREQ= DM = SUB =20 FREQ= DM = SUB =375 FREQ=999 DM = (d) (e) (f) 圖 4- 鎚骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形 29

39 SUB =689 FREQ=3000 DM =30.25 SUB =09 FREQ=4000 DM = SUB =34 FREQ=4999 DM =7.933 (g) (h) (i) SUB =639 FREQ=6000 DM = SUB =908 FREQ=7003 DM =68.89 SUB =250 FREQ=7998 DM =88.52 (j) (k) (l) 圖 4-2 鎚骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 30

40 SUB =0 FREQ= DM =6.338 SUB =6 FREQ= DM = SUB =32 FREQ=50.94 DM =22.58 (a) (b) (c) SUB =57 FREQ= DM = SUB =92 FREQ= DM =8.93 SUB =286 FREQ=999 DM =27.5 (d) (e) (f) 圖 4-3 鉆骨 25/ 250/ 500/ 750/ k/ 2kHz 模態振形 3

41 SUB =52 FREQ=300 DM =39.47 SUB =772 FREQ=3998 DM = SUB =06 FREQ=4997 DM =57.59 (g) (h) (i) SUB =250 FREQ=6002 DM = SUB =464 FREQ=7000 DM = SUB =655 FREQ=8003 DM =8.58 (j) (k) (l) 圖 4-4 鉆骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 32

42 SUB =4 FREQ= DM = SUB =72 FREQ= DM =2258 SUB =9 FREQ= DM =4092 (a) (b) (c) SUB =99 FREQ=002 DM = SUB =95 FREQ= DM = SUB =7 FREQ=994 DM =30.42 (d) (e) (f) 圖 4-5 鐙骨 25/ 250/ 500/ 750/ k/ 2k Hz 模態振形 33

43 SUB =39 FREQ=2994 DM = SUB =64 FREQ=40 DM =43.55 SUB =9 FREQ=502 DM =36.4 (g) (h) (i) SUB =27 FREQ=5992 DM =34.22 SUB =24 FREQ=7020 DM = SUB =262 FREQ=7993 DM = (j) (k) (l) 圖 4-6 鐙骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 34

44 SUB =7 FREQ= DM =6.27 SUB =29 FREQ= DM =5.733 SUB =60 FREQ= DM =23.87 (a) (b) (c) SUB =4 FREQ= DM =4.49 SUB =90 FREQ= DM =4.954 SUB =69 FREQ=999 DM =2.89 (d) (e) (f) 圖 4-7 鉆骨與鎚骨 25/ 250/ 500/ k/ 2k Hz 模態振形 35

45 SUB =49 FREQ=300 DM = SUB =73 FREQ=400 DM = SUB =2264 FREQ=5000 DM = (g) (h) (i) SUB =2779 FREQ=600 DM =6.26 SUB =3249 FREQ=6999 DM =5.707 SUB =3672 FREQ=800 DM =52.50 (j) (k) (l) 圖 4-8 鉆骨與鎚骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 36

46 SUB =35 FREQ= DM = SUB =60 FREQ= DM =5074 SUB =90 FREQ= DM = (a) (b) (c) SUB =8 FREQ= DM = SUB =58 FREQ=00 DM = SUB =37 FREQ=200 DM = (d) (e) (f) 圖 4-9 鉆骨與鐙骨 25/ 250/ 5k/ 750/ k/ 2k Hz 模態振形 37

47 SUB =629 FREQ=3000 DM =38.24 SUB =904 FREQ=3999 DM = SUB =76 FREQ=4999 DM = (g) (h) (i) SUB =436 FREQ=6003 DM = SUB =673 FREQ=6999 DM =60.04 SUB =888 FREQ=8002 DM = (j) (k) (l) 圖 4-0 鉆骨與鐙骨 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 38

48 SUB =42 FREQ= DM = SUB =73 FREQ= DM =26.64 SUB =8 FREQ= DM =33.2 (a) (b) (c) SUB =75 FREQ= DM =4.54 SUB =255 FREQ= DM =4.582 SUB =703 FREQ=999 DM =2.94 (d) (e) (f) 圖 4- 聽小骨鏈 25/ 250/ 5k/ 750/ k/ 2kHz 模態振形 39

49 SUB =256 FREQ=300 DM = SUB =845 FREQ=400 DM = SUB =2424 FREQ=5000 DM =62.52 (g) (h) (i) SUB =2964 FREQ=6000 DM = SUB =3459 FREQ=700 DM =45.39 SUB =3904 FREQ=8000 DM = (j) (k) (l) 圖 4-2 聽小骨鏈 3k/ 4k/ 5k/ 6k/ 7k/ 8k Hz 模態振形 40

50 4-2 探討模態分析結果在臨床上的意義 中耳聽小骨鏈模態分析, 最主要是希望可以提供臨床醫師針對需植入聽小骨義肢患者術前植入物的評估 表 4- 為正常人耳右側聽小骨鏈之自然頻率與臨床實務上測試頻率是否達成共振條件的對照表 分析結果中, 在個別聽小骨部分, 鎚骨則較易與刺激音 3k Hz 4k Hz 6k Hz 產生共振現象 ; 鉆骨於 7k Hz 易產生共振現象 ; 鐙骨皆不易與測試頻率達成共振 此外, 聽骨偶合結構方面, 鉆骨與鎚骨偶合易與 250 Hz k Hz 5k Hz 頻率下達成共振 鉆骨與鐙骨偶合於 3k Hz 頻率下達成共振 三根聽小骨偶合易於 k Hz 5k Hz 6k Hz 8k Hz 等頻率下達成共振 所以, 為了避免術後植入物與原本聽骨鏈偶合後, 因聽小骨鏈易與刺激音共振而導致患者聽損無法改善, 甚至完全喪失聽力, 術前植入物的評估的確是不可或缺的首要工作 4

51 表 4- 六種形式聽骨組合結構與臨床測試頻率模態分析結果 測試頻率 組合結構 鎚骨鉆骨鐙骨鎚骨 / 鉆骨鉆骨 / 鐙骨聽骨鏈 25 Hz Hz Hz Hz k Hz k Hz k Hz k Hz k Hz k Hz k Hz k Hz 聽骨偶合結構自然頻率與測試頻率達成共振條件 (±0.5 Hz) 42

52 4-3 調和力分析結果 調和力分析則針對鎚骨柄部與鼓膜連接處, 給予 90dB (0.623Pa) 之負載作用於 3.08mm 2 鎚骨柄部底面積, 檢視鐙骨足板上特定節點, 如圖 4-3, 於 250Hz 至 8kHz 範圍內頻率與振幅的響應, 分別以 三方向檢視其結果, 如圖 4-4 至 4-6 所示 圖 4-3 鐙骨足板上節點及對應節點編號 43

53 E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No. U U U E FREQUENC (Hz) E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.2 U U U E FREQUENC (Hz) E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.3 U U U E FREQUENC (Hz) 圖 4-4 節點 ~3 號頻率響應圖 44

54 E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.4 U U U E FREQUENC (Hz) E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.5 U U U E FREQUENC (Hz) E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.6 U U U E FREQUENC (Hz) 圖 4-5 節點 4~6 號頻率響應圖 45

55 E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.7 U U U E FREQUENC (Hz) E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.8 U U U E FREQUENC (Hz) E-0 AMPLIDUTE (mm) E-03 E-05 E-07 E-09 No.9 U U U E FREQUENC (Hz) 圖 4-6 節點 7~9 號頻率響應圖 46

56 由圖 4-4~6 可以發現鐙骨足板上九個節點於 k~2khz 之間有一處明顯轉折點, 證明聽骨承受低於 k Hz 的刺激頻率時, 聽骨仍會以振動形式呈現, 其位移即刺激音傳至鐙骨足板的振幅大小 ; 然而高於 2k Hz 刺激頻率時, 聽骨幾乎不會發生位移, 而是採用聽骨直傳傳遞能量 針對, 鐙骨足板節點上頻率響應與 Sun[2] 結果出入甚大, 主要來自於負載大小與位置差異 : 本研究中以聽骨鏈為研究主體, 負載施予鎚骨柄部底面積與 Sun[2] 負載施予鼓膜上, 兩者表面積相差有 22 倍之餘, 因而無法切確得知鼓膜推動聽小骨鏈之力矩效應 但也發現針對聽小骨鏈給予負載, 並不容易推動聽骨鏈擺動, 而是需倚靠鼓膜推動聽小骨鏈進而達到擴大振幅的效果 47

57 4-4 調和力分析結果與文獻資料比對 本研究分析結果與文獻中 Wada [] Sun [2] Gan [5] 三者比 較, 針對聽小骨鏈有限元素幾何外形 ( 表 4-2) 材料參數 ( 表 4-3~4-4) 與邊界條件 ( 表 4-5~4-6) 所示 表 4-2 幾何外形尺寸比較表 Sun Gan. Wada present Malleus Length from the end of manubrium to the end of lateral process (mm) Total length (mm) *Density (EQV) (kg/m 3 ) 3.59 * * ~6.2* *0 3 Weight (mg) Incus Length along long process (mm) Length along short process (mm) *Density (EQV) (kg/m 3 ) 3.23* * ~6.2* *0 3 Weight (mg) Stapes Height (mm) Length of footplate (mm) Width of footplate (mm) Area of footplate(mm 2 ) *Density (EQV)(kg/m 3 ) 2.20 * * ~6.2* *0 3 Weight (mg) 其中, 等效密度 (equivalent density) 計算方式以鎚骨為例, 鎚骨 頭部 頸部 柄部三個部位, 分別將其體積乘上該部位密度, 並將 其加總後, 再除以鎚骨總體積, 即為鎚骨等效密度 鉆骨與鐙骨也 48

58 使用此法計算其等效密度 表 4-3 材料參數比較表 (I) Sun Gan Wada present Eardrum Density(kg/m 3 ).2*0 3.2*0 3.2*0 3 - Pars tensa (N/m 2 ) 2.0*0 7 (c) 2.0*0 7 (c) 3.2*0 7 (r) 3.5*0 7 (r) 3.34*0 7 - Pars flaccida (N/m 2 ).0*0 7.0*0 7.*0 7 - Malleus Head (kg/m 3 ) 2.55* * *0 3 Neck (kg/m 3 ) 4.53* * ~6.2* *0 3 Handle (kg/m 3 ) 3.70* * *0 3 oung s module (N/m 2 ).4*0 0.4*0 0.2*0 0.4*0 0 Incus Body (kg/m 3 ) 2.36* * *0 3 Short process (kg/m 3 ) 2.26* * ~6.2* *0 3 Long process (kg/m 3 ) 5.08* * *0 3 oung s module (N/m 2 ).4*0 0.4*0 0.2*0 0.4*0 0 Stapes Stapes(kg/m 3 ) 2.20* * ~6.2* *0 3 oung s module (N/m 2 ).4*0 0.4*0 0.2*0 0.4*0 0 其中,Wada 給予材料參數資訊只有耳膜與聽小骨兩部分,Sun & Gan 則是將其細分各自部位, 本研究也是參照 Sun 的材料參數定 義, 唯獨耳膜部份並未重建 49

59 表 4-4 材料參數比較表 (II) Sun Gan Wada present Incudomallealar joint Density (kg/m 3 ) 3.2* * oung s module (N/m 2 ).4*0 0.4* Proportional damping α=0 s - α=0 s β=5.0*0-4 s β=7.5*0-5 s Incudostapedial joint Density (kg/m 3 ).2*0 3.2* *0 3 oung s module (N/m 2 ) 6.0* * * *0 5 Proportional damping α=0 s - α=0 s - α=0 s - α=0 s - β=5.0*0-4 s β=7.5*0-5 s β=3.72*0-4 s β=5.0*0-4 s Malleus attachment on the eardrum Density (kg/m 3 ).0*0 3.0* oung s module (N/m 2 ) 4.7* * Proportional damping α=0 s - α=0 s - β=5.0*0-4 s β=7.5*0-5 s - - Poisson ratio (ALL) 其中,Sun Gan 將其上述各部位密度 楊氏係數皆相同, 唯獨 Gan 假設其比例阻尼 β 略小於 Sun 假設值,Wada 則是單獨只設定 I- S joint 此部材料參數 本研究因考慮其 I-M joint 材料參數相近於聽 小骨, 因而直接將其視為相連部位, 然而 I-S joint 材料參數與聽小 骨差異甚大, 所以考慮其影響重建於有限元素模型中 另外, 鎚骨 與鼓膜連接面則因鼓膜並未重建, 因而不加以考慮 50

60 表 4-5 邊界條件比較表 (I) Sun Gan Wada present Superior mallear ligament (C) oung s module (N/m 2 ) 4.9* * Proportional damping α=0 s β=5.0*0-4 s Lateral mallear ligament (C2) oung s module (N/m 2 ) 6.7* * Posterior mallear ligament (C3) oung s module (N/m 2 ) 6.5* *0 6 - Damping (N*s/m) Anterior mallear ligament (C4) oung s module (N/m 2 ) 2.*0 6 2.*0 7 - Damping (N*s/m) Posterior stapedial muscle (C5) oung s module (N/m 2 ) 5.2* * Tensor tympani muscle (C7) oung s module (N/m 2 ) 2.6* * *0 6 - 邊界條件部分, 表 4-5 所示,Sun Gan 部份材料參數差異甚 鉅, 如 C4 C5 C7, 差距 ~2 個數量級 然而,Wada 則捨棄 C~C5, 僅考慮 C7 的效應 本研究因無法從電腦斷層影像中擷取出 影像邊界, 因而也不考慮其邊界效應 5

61 表 4-6 邊界條件比較表 (II) Sun Gan Wada present Cochlear fluid (C6) Spring constant (N/m) Damp (N*s/m) * Stapedius annular ligament Spring constant (N/m) *0 5 (N/m 2 ) 9 Proportional damping α=0 s β=3.72*0-4 s 0 Stapedius annular muscle oung s module (N/m 2 ) *05(N/m 2 ) - Proportional damping α=0 s β=3.72*0-4 s - Tympanic annulus oung s module (N/m 2 ) 6.0* * Proportional damping α=0 s - α=0 s - β=0.000 s β=7.5*0-5 s - - Air Bulk modulus (N/m 2 ) *0 5 - Density (kg/m 3 ) - -.8*0 3 - 如表 4-6 所示, 無論是 Sun Gan Wada 以及本研究皆有考慮 C6, 模擬其內耳流體造成的阻尼力,Stapedius annular ligament 則是 模擬鐙骨足板於空間中擺動的彈簧阻尼元素 其餘, 肌肉 (muscle) 與 韌帶則是互相有無, 各家考慮不同的邊界條件 52

62 最後, 針對以上各家聽小骨鏈重建有限元素模型方法, 不難發 現以下缺失 : () 聽小骨幾何尺寸量測並無明顯分界, 臨床實務方面也並未定 義, 因而只能依照重建完畢後總體積, 與各自分野其細部材料 參數 ( 密度 ) 後, 所計算得出聽小骨重量再進行比對 (2) Sun Gan 對與聽小骨鏈邊界條件的假設是考慮每一各環節, 如聽小骨鏈接點 肌肉 韌帶 阻尼力 等 ;Wada 則是取其較為影響的鉆骨與鐙骨接點 鐙股足板上肌肉與環狀韌帶 等, 部份邊界條件作為其研究對象的主要邊界條件設定 ; 本研究也針對已重建完畢之聽小骨鏈立體影像, 在不影響重建影像幾何外形的條件下, 將聽小骨鏈影響較大的部位, 如鉆骨與鐙骨接點 鐙骨足板上的環狀韌帶與彈簧阻尼元素, 作為本研究對象之主要邊界條件的假設 53

63 第五章結論 本研究配合臨床實務上的需求, 利用人耳電腦斷層掃描圖, 針 對活體聽小骨鏈以影像重建方式還原其立體模型, 再導入有限元素 法暸解聽小骨鏈動態特性 目前研究成果有 : () 直接由電腦斷層掃描影像重建人耳活體聽小骨影像, 可快速提供臨床醫師了解患者聽骨鏈之空間分佈 (2) 立體聽小骨影像可由 CAD 軟體 SolidWorks 成功導入 CAE 軟體 ANSS, 利用有限元素法了解其動態特性, 可提供術前替代物植入的評估 (3) 個別聽小骨較不易與臨床測試頻率產生共振, 但聽骨偶合後其自然頻率將接近於臨床測試頻率, 但也不完全等於測試頻率 (4) 聽骨鏈接受低頻 (<k Hz) 刺激時, 仍以振動形式呈現所承受刺激音的振幅大小 當刺激音高於 2k Hz 時, 聽骨鏈幾乎沒有位移產生, 傳導方式改為聽骨直傳方式傳遞 54

64 研究成果之臨床實務意義 : () 需植入鉆骨 鐙骨義肢時, 須注意義肢與聽骨鏈偶合後是否易與測試頻率達到共振的條件 (2) 幫助鐙骨手術 聽骨鏈病變與聽骨老化導致聽損患者, 提供術前評估 (3) 可追蹤術後植入物仍無法挽救聽損的原因 55

65 參考文獻. Samy R.N., Pensak M.L., Revision ossiculoplasty, Otolaryngol. Clin. North Am., Vol. 39(4), pp , Petrova P., Freeman S., Sohmer H., The effects of positive and negative middle ear pressures on auditory threshold, Otol. Neurotol., Vol. 27(5), pp , Puria S., Kunda L.D., Roberson J.B. Jr., Perkins R.C., Malleus-tofootplate ossicular reconstruction prosthesis positioning: cochleovestibular pressure optimization, Otol. Neurotol., Vol. 26(3), pp , Hudde H., Weistenhofer C., Key features of the human middle ear, ORL J. Otorhinolaryngol. Relat. Spec., Vol. 68(6), pp , Maassen M.M., Rodriguez Jorge J., Herberhold S., Vonthein R., immermann R., Baumann I., Brosch S., Mauz P.S., Pfister M., alaman I.M., Lowenheim H., enner H.P., Safe and reliable sound threshold measures with direct vibration of the ossicular chain, Laryngoscope, Vol. 4(), pp , Li C.J., hu P.F., Liu.H., Wang.G., ang C., Chen H.B., Ning., hou J.H., Chen J., Comparative observation of protective effects of earplug and barrel on auditory organs of guinea pigs exposed to experimental blast underpressure, Chin. J. Traumatol., Vol. 9(4), pp , Rosowski J.J., Ravicz M.E., Songer J.E., Structures that contribute to middle-ear admittance in chinchilla, J. Comp. Physiol. A Neuroethol. Sens. Neural. Behav. Physiol., Vol. 92(2), pp ,

66 8. Willi U.B., Bronner G.N., Narins P.M., Middle ear dynamics in response to seismic stimuli in the Cape golden mole (Chrysochloris asiatica), J. Exp. Biol., Vol. 209(Pt 2), pp , Wanda H., Metoki T., Kobayashi T., Analysis of dynamic behavior of human middle ear using a finite-element method, J. Acoust. Soc. Am., Vol. 92 (6), pp , Ferris P., Prendergast P.J., Middle-ear dynamics before and after ossicular replacement, J. Biomech, Vol. 33(5), pp , Koike T., Wanda H., Kobayashi T., Modeling of human middle ear using the finite-element method, J. Acoust. Soc. Am., Vol. (3), pp , Sun Q., Gan R.., Chang K.H., Dormer K.J., Computer-integrated finite element modeling of human middle ear, Biomech Model Mechanobiol, Vol. (2), pp , Gan R.., Sun Q., Dyer R.K., Chang K.H., Dormer K.J., Threedimensional modeling of middle ear biomechanics and its applications, Otol. Neurotol., Vol. 23(3), pp , Kelly D.J., Prendergast P.J., Blayney A.W., The effect of prosthesis design on vibration of the reconstructed ossicular chain:a comparative finite element analysis of four prostheses, Otol. Neurotol., Vol. 24(), pp. -9, Gan R.., Feng B., Sun Q., Three-dimensional finite element modeling of human ear for sound transmission, Ann Biomed Eng, Vol. 32(6), pp , Koike T., Shinozaki M., Murakami S., Homma K., Kobayashi T., Wada H., Effects of Individual Differences in Size and Mobility of 57

67 the Middle Ear on Hearing, JSME International Journal, Vol. 48(4), pp , Gan R.., Sun Q., Feng B., Wood M.W., Acoustic-structural coupled finite element analysis for sound transmission in human ear -Pressure distributions, Med. Eng. Phys., Vol. 28(5), pp , Lee C.F., Chen P.R., Lee W.J., Chen J.H., Liu T.C., Computed aided three-dimensional reconstruction and modeling of middle ear biomechanics by high-resolution computed tomography and finite element analysis, Biomed. Eng., Vol. 8(5), pp. 24-2, Tonndorf J., Khanna S.M., Tympanic-membrane vibrations in human cadaver ears studied by time-averaged holography, J. Acoust. Soc. Am., Vol. 52, pp , Funnell W.R.J., Laszlo C.A., A critical review of experimental observations on eardrum structure and function, ORL. J. Otorhinolaryngol. Relat. Spec., Vol. 44, pp , Vlaming M.S.M.G., Middle ear mechanics by laser Doppler interferometry, Techniche Universiteit Delft, The Netherlands, Hüttenbrink K.B., The mechanics of the middle-ear at static air pressures: the role of the ossicular joints, the function of the middleear muscles and the behaviour of stapedial prostheses, Acta Otolaryngol Suppl, Vol. 45, pp. -35, Wada H., Koike T., Kobayashi T., Three-dimensional finite-element method (FEM) analysis of the human middle ear, In: Hüttenbrink KB (ed) Middle ear mechanics in research and otosurgery. Department of Oto-Rhino-Laryngology, Dresden University of Technology, pp ,

68 24. William A. ost, Fundamentals of Hearing: An Introduction, Academic Press, fourth edition,