Structure movement mechanism for auditory ossicle chain of human middle ear based on finite element method
-
摘要:
目的 运用中耳Micro-CT扫描数据进行三维重建,采用有限元方法对听骨链运动规律及鼓膜振动特性进行生物力学研究。 方法 通过对比鼓膜凸、镫骨足板振幅实验数据验证模型,并进行频率响应分析及模态分析。 结果 在不同频率下听骨链振动与转动均不一样,三块听小骨围绕一定的转动轴进行转动和摇摆运动,接近鼓膜凸处的锤骨柄做同相位转动和平动,镫骨足板做活塞平动。在频率1000 Hz下鼓膜整体弯曲变形局部高达2.32e-006 m,出现在环韧带附近;鼓膜凸最大变形约1.0e-007 m。鼓膜环韧带在声压激励下容易发生扭曲变形。鼓膜凸附近出现应力集中,最大应力约8.33e-004 MPa。 结论 听骨链运动机理研究对人耳生命科学研究和临床手术均有一定的理论指导意义。 Abstract:Objective To reconstruct the middle ear structure in three-dimensional by Micro-CT data, and explore the motion rule of chain ossicles and tympanic membrane vibration characteristics in biomechanical theory research by using the finite element method. Methods The model was verified correctly by comparing finite element simulation data with experimental data about the amplitude of tympanic membrane umbo and stapes footplate. The conclusion was drawn by frequency response analysis and modal analysis. Results (1) The vibration and rotation of ossicular chain were different under various frequencies. Three ossicles performed rotating and rocking motion around axis of rotation. Manubrium mallei close to tympanic membrane umbo performed in same phase rotation and translation. Stapes footplate performed translation like the piston. (2) At the frequency of 1000 Hz, the total bending deformation of tympanic membrane ran up to 2.32e-006 m, which appeared in annular ligament of tympanic membrane. The maximum deformation of the membrane tympani umbo was about 1.0e-007 m. Under the excitation of sound pressure, the annular ligament of tympanic membrane was prone to distortion. The stress concentration occurred in the vicinity of membrane tympani umbo.The maximum stress was 8.33e-004 MPa. Conclusion Motion of ossicular chain will provide theoretical guidance for life science research and clinical operation of human ears. -
表 1 中耳听骨链、各韧带、各肌腱的材料参数
听小骨结构 弹性模量(MPa) 密度(kg·m-3) 软组织结构 弹性模量(MPa) 密度(kg·m-3) 锤骨 锤骨头 14100 2550 鼓膜紧张部 35 1200 锤骨颈 14100 4530 鼓膜松弛部 10 1200 锤骨柄 14100 3700 鼓膜环韧带 0.6 1200 砧骨 砧骨体 14100 2360 锤骨上悬韧带 4.9 2500 砧骨短突 14100 2260 锤骨侧韧带 6.7 2500 砧骨长突 14100 5080 锤骨前韧带 2.1 2500 镫骨 14100 2200 砧骨上悬韧带 4.9 2500 锤-砧关节 14100 3200 砧骨后韧带 6.5 2500 砧-镫关节 0.6 1200 镫骨环韧带 0.2 1200 / / / 鼓膜张肌 2.6 2500 / / / 镫骨肌 5.2 2500 -
[1] Lesser TH, Williams KR. The tympanic membrane in cross section: a finite element analysis[J]. J Laryngol Otol, 1988, 102(3): 209-14. doi: 10.1017/S0022215100104542 [2] Wada H, Metoki T, Kobayashi T. Analysis of dynamic behavior of human middle ear using a finite-element method[J]. J Acoust Soc Am, 1992, 92(6): 3157-68. doi: 10.1121/1.404211 [3] Wada H, Koike T, Kobayashi T. Middle ear mechanics in research and otosurgery[C]//Proceedings of the International Workshop on Middle Ear Mechanics, 1996: 76-81. [4] Koike T, Wada H, Kobayashi T. Modeling of the human middle ear using the finite-element method[J]. J Acoust Soc Am, 2002, 111(3): 1306-17. doi: 10.1121/1.1451073 [5] Ferris P, Prendergast PJ. Middle-ear dynamics before and after ossicular replacement[J]. J Biomech, 2000, 33(5): 581-90. doi: 10.1016/S0021-9290(99)00213-4 [6] 佚名. 标题为空[C]//Proceedings of the International Workshop on Middle Ear Mechanics, Dresden, Sept.19-22, 1996, 出版年度缺失: 页码范围缺失. [7] Gan RZ, Feng B, Sun Q. Three-dimensional finite element modeling of human ear for sound transmission[J]. Ann Biomed Eng, 2004, 32(6): 847-59. doi: 10.1023/B:ABME.0000030260.22737.53 [8] Gan RZ, Sun Q, Feng B, et al. Acoustic-structural coupled finite element analysis for sound transmission in human ear--pressure distributions[J]. Med Eng Phys, 2006, 28(5): 395-404. doi: 10.1016/j.medengphy.2005.07.018 [9] Gan RZ, Reeves BP, Wang X. Modeling of sound transmission from ear canal to cochlea[J]. Ann Biomed Eng, 2007, 35(12): 2180-95. doi: 10.1007/s10439-007-9366-y [10] Gan RZ, Wang X. Multifield coupled finite element analysis for sound transmission in otitis media with effusion[J]. J Acoust Soc Am, 2007, 122(6): 3527-38. doi: 10.1121/1.2793699 [11] Gan RZ, Cheng T, Dai C, et al. Finite element modeling of sound transmission with perforations of tympanic membrane[J]. J Acoust Soc Am, 2009, 126(1): 243-53. doi: 10.1121/1.3129129 [12] Kelly DJ, Prendergast PJ, Blayney AW. The effect of prosthesis design on vibration of the reconstructed ossicular chain: a comparative finite element analysis of four prostheses[J]. Otol Neurotol, 2003, 24(1): 11-9. doi: 10.1097/00129492-200301000-00004 [13] 刘迎曦, 李 生, 孙秀珍. 人耳鼓膜病变数值分析[J]. 医用生物力学, 2008, 23(4): 275-8. http://www.cnki.com.cn/Article/CJFDTOTAL-YISX200804006.htm [14] 刘迎曦, 李 生, 孙秀珍. 人耳传声数值模型[J]. 力学学报, 2008, 40(1): 107-13. doi: 10.6052/0459-1879-2008-1-2007-051 [15] Liu YX, Li S, Sun XZ. Numerical analysis of ossicular chain lesion of human ear[J]. Acta Mech Sin, 2009, 25(2): 241-7. doi: 10.1007/s10409-008-0206-6 [16] 姚文娟, 李 武, 付黎杰, 等. 中耳结构数值模拟及传导振动分析[J]. 系统仿真学报, 2009, 21(3): 651-4. http://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ200903008.htm [17] 姚文娟, 李晓青, 李 武, 等. 中耳病变及人工镫骨形体研究[J]. 医用生物力学, 2009, 24(2): 118-22. http://www.cnki.com.cn/Article/CJFDTOTAL-YISX200902009.htm [18] Yao W, Li B, Huang X, et al. Restoring hearing using total ossicular replacement prostheses--analysis of 3D finite element model[J]. Acta Otolaryngol, 2012, 132(2): 152-9. doi: 10.3109/00016489.2011.633229 [19] Yao W, Li B, Guo C, et al. Numerical simulation and study about partial ossicular replacement prostheses with several various materials on hearing restoration[J]. Advanc Mater Res, 2011, 238(3): 2899-904. [20] 刘后广, 塔 娜, 饶柱石. 悬浮振子对中耳声传播特性影响的数值研究[J]. 力学学报, 2010, 42(1): 109-14. doi: 10.6052/0459-1879-2010-1-2008-519 [21] 王振龙, 王学林, 胡于进, 等. 基于中耳与耳蜗集成有限元模型的耳声传递模拟[J]. 中国生物医学工程学报, 2011, 30(1): 60-6. http://www.cnki.com.cn/Article/CJFDTOTAL-ZSWY201101012.htm [22] Herrmann G, Liebowitz H. Mechanics of bone fractures[M]. new york: academic press, 1972: 772-840. [23] Kirikae I. The structure and function of the middle ear[M]. tokyo: [s.n.], 1960: 51-7.