Rong Z. Gan

Three-dimensional finite element modeling of human ear for sound transmission

An accurate, comprehensive finite element model of the human ear can provide better understanding of sound transmission, and can be used for assessing the influence of diseases on hearing and the treatment of hearing loss. In this study, we proposed a three-dimensional finite element model of the human ear that included the external ear canal, tympanic membrane (eardrum), ossicular bones, middle ear suspensory ligaments/muscles, and middle ear cavity. This model was constructed based on a complete set of histological section images of a left ear temporal bone. The finite element (FE) model of the human ear was validated by comparing model-predicted ossicular movements at the stapes footplate and tympanic membrane with published experimental measurements on human temporal bones. The FE model was employed to predict the effects of eardrum thickness and stiffness, incudostapedial joint material, and cochlear load on acoustic-mechanical transmission through the human ossicular chain. The acoustic-structural coupled FE analysis between the ear canal air column and middle ear ossicles was also conducted and the results revealed that the peak responses of both tympanic membrane and stapes footplate occurred between 3000 and 4000 Hz.

Modeling of Sound Transmission from Ear Canal to Cochlea

A 3-D finite element (FE) model of the human ear consisting of the external ear canal, middle ear, and cochlea is reported in this paper. The acoustic-structure-fluid coupled FE analysis was conducted on the model which included the air in the ear canal and middle ear cavity, the fluid in the cochlea, and the middle ear and cochlea structures (i.e., bones and soft tissues). The middle ear transfer function such as the movements of tympanic membrane, stapes footplate, and round window, the sound pressure gain across the middle ear, and the cochlear input impedance in response to sound stimulus applied in the ear canal were derived and compared with the published experimental measurements in human temporal bones. The frequency sensitivity of the basilar membrane motion and intracochlear pressure induced by sound pressure in the ear canal was predicted along the length of the basilar membrane from the basal turn to the apex. The satisfactory agreements between the model and experimental data in the literature indicate that the middle ear function was well simulated by the model and the simplified cochlea was able to correlate sound stimulus in the ear canal with vibration of the basilar membrane and pressure variation of the cochlear fluid. This study is the first step toward the development of a comprehensive FE model of the entire human ear for acoustic-mechanical analysis.

Human middle ear transfer function measured by double laser interferometry system.

Hypothesis: Simultaneous measurements of vibrations on the stapes footplate, incudostapedial (IS) joint, and tympanic membrane (TM) can be made in both normal and drained cochleae, and the stapes displacement transfer function (S-DTF) and TM displacement transfer function (TM-DTF) are derived. Background: A single laser Doppler interferometer previously has been used for measuring movement of the stapes or TM in temporal bones. However, there may be a limitation to optimally describing acoustic–mechanical transmission when the interferometer and temporal bone are moved frequently during experimental recordings. Simultaneous measurements of vibrations of the TM and stapes footplate, or TM and IS joint may reveal different acoustic–mechanical characteristics of the middle ear. Methods: Dual laser interferometers simultaneously measured vibrations of the TM, IS joint, and stapes in 10 temporal bones with both intact and drained cochleae. From these measurements, the middle ear transfer function was expressed as the S-DTF, TM-DTF, and displacement transmission ratio (DTR). Results: Simultaneous displacements of the TM, IS joint, and stapes footplate induced by sound pressure in the ear canal were recorded in both amplitude and phase. The middle ear transfer functions in terms of displacement ratio confirmed published single interferometer data but provided new information from drained cochlea. Conclusion: Stapes and TM displacement transfer functions were determined using dual interferometry, provided accurate amplitude and phase relationships from stapes footplate, IS joint, and TM, with new data from drained and normal cochlea.

Mass loading on the ossicles and middle ear function.

The middle ear as a levered vibrating system for sound transmission from the external to the inner ear is affected by changes in ossicular chain mass. Mass loading of the ossicles may impair ossicular dynamics and sound transmission to the inner ear. It is incumbent on otologic surgeons and researchers of middle ear mechanics to consider the mass loading effect on middle ear function in clinical and physiological applications. The residual hearing and frequency response can change after surgery or implantation of middle ear prostheses. We conducted experiments on mass loading effects on the middle ear transfer functions by using laser Doppler interferometry and a human temporal bone model. Two implant mass loading conditions were tested on 17 fresh or fresh-frozen temporal bones and compared with the unloaded condition for the frequencies 250 to 8,000 Hz. The results show that the linearity of the middle ear function did not change, although displacement of the stapes footplate decreased after the increased masses were placed on the incudostapedial joint. The greater the mass of the implant, the less displacement was measured at the stapes footplate. We conclude that there is a quantitative limit to increased mass on the ossicular chain above which the mass will remarkably impair hearing thresholds.

An advanced computer-aided geometric modeling and fabrication method for human middle ear

This paper presents a practical and systematic method for reconstructing accurate computer and physical models of the entire human middle ear. The proposed method starts with the histological section preparation of human temporal bone. Through tracing outlines of the middle ear components on the sections, a set of discrete points is obtained and employed to construct B-spline curves that represent the exterior contours of the components using a curve-fitting technique. The surface-skinning technique is then employed to quilt the B-spline curves for smooth boundary surfaces of the middle ear components using B-spline surfaces. The solid models of the middle ear components are constructed using these surfaces and then assembled to create the entire middle ear in a computer-aided design environment. This method not only provides an effective way to visualize and measure the three-dimensional structure of the middle ear, but also provides a detailed knowledge of middle ear geometry that is required for finite element analysis or multibody dynamic analysis of the human middle ear. In addition, the geometric model constructed using the proposed method is smooth and can be fabricated in various scales using solid freeform fabrication technology. The physical model of the human middle ear is extremely effective in realizing the middle ear anatomy and enhancing discussion and collaboration among researchers and physicians.

Characterization of the linearly viscoelastic behavior of human tympanic membrane by nanoindentation

Human tympanic membrane (or eardrum) is composed of three membrane layers with collagen fibers oriented in the radial and circumferential directions, and exhibits viscoelastic behavior with membrane (or in-plane) properties different from through-thickness (or out-of-plane) properties. Due to the interaction of bundled fibers and ground substance, which is inhomogeneous, these properties could change with locations. In this paper, we use nanoindentation techniques to measure the viscoelastic functions of four quadrants of tympanic membrane (TM). For measurement of in-plane Youngs relaxation modulus we fixed a sectioned quadrant of the TM on a circular hole and used a spherical nanoindenter tip to apply force at the center of the suspended circular portion of the specimen. An inverse problem solving methodology was employed using finite element method to determine the average in-plane Youngs relaxation modulus of the TM quadrant. Results indicate that the in-plane steady-state Youngs relaxation modulus for four quadrants of the TM does not vary significantly. However, a variation of the modulus from 25.73 MPa to 37.8 MPa was observed with measurements from different individuals. For measurement of Youngs relaxation modulus in the through-thickness direction a spherical indenter tip was used to indent into different locations on the surface of the TM specimen supported by a substrate. Viscoelastic contact mechanics analysis of the load-displacement curve, representative primarily of the through-thickness stiffness of the TM, was conducted to extract the Youngs relaxation modulus in the out-of-plane direction. Results indicate a wide variation in steady-state Youngs relaxation modulus, from 2 MPa to 15 MPa, in the through-thickness direction over the TM.

Finite-element analysis of middle-ear pressure effects on static and dynamic behavior of human ear

A finite-element analysis for static behavior of middle ear under variation of the middle-ear pressure was conducted in a 3D model of human ear by combining the hyperelastic Mooney-Rivlin material model and geometry nonlinearity. An empirical formula was then developed to calculate material parameters of the middle-ear soft tissues as the stress-dependent elastic modulus relative to the middle-ear pressure. Dynamic behavior of the middle ear in response to sound pressure in the ear canal was predicted under various positive and negative middle-ear pressures. The results from static analysis indicate that a positive middle ear pressure produces the static displacements of the tympanic membrane (TM) and footplate more than a negative pressure. The dynamic analysis shows that the reductions of the TM and footplate vibration magnitudes under positive middle-ear pressure are mainly determined by stress dependence of elastic modulus. The reduction of the TM and footplate vibrations under negative pressure was caused by both the geometry changes of middle-ear structures and the stress dependence of elastic modulus.

Finite element modeling of sound transmission with perforations of tympanic membrane

A three-dimensional finite element (FE) model of human ear with structures of the external ear canal, middle ear, and cochlea has been developed recently. In this paper, the FE model was used to predict the effect of tympanic membrane (TM) perforations on sound transmission through the middle ear. Two perforations were made in the posterior-inferior quadrant and inferior site of the TM in the model with areas of 1.33 and 0.82 mm(2), respectively. These perforations were also created in human temporal bones with the same size and location. The vibrations of the TM (umbo) and stapes footplate were calculated from the model and measured from the temporal bones using laser Doppler vibrometers. The sound pressure in the middle ear cavity was derived from the model and measured from the bones. The results demonstrate that the TM perforations can be simulated in the FE model with geometrical visualization. The FE model provides reasonable predictions on effects of perforation size and location on middle ear transfer function. The middle ear structure-function relationship can be revealed with multi-field coupled FE analysis.

Fixation and detachment of superior and anterior malleolar ligaments in human middle ear: Experiment and modeling

The aim of this study is to investigate the function of the superior malleolar ligament (SML) and the anterior malleolar ligament (AML) in human middle ear for sound transmission through simulations of fixation and detachment of these ligaments in human temporal bones and a finite element (FE) ear model. Two laser vibrometers were used to measure the vibrations of the tympanic membrane (TM) and stapes footplate. A 3-D FE ear model was used to predict the transfer function of the middle ear with ligament fixation and detachment. The results demonstrate that fixations and detachments of the SML and AML had different effects on TM and stapes footplate movements. Fixation of the SML resulted in a reduction of displacement of the TM (umbo) and the footplate at low frequencies (f<1000 Hz), but also caused a shift of displacement peak to higher frequencies. Fixation of both SML and AML caused a reduction of 15 dB at umbo or stapes at low frequencies. Detachment of the SML had almost no effect on TM and footplate mobility, but AML detachment had a minor effect on TM and footplate movement. The FE model was able to predict the effects of SML and AML fixation and detachment.

A Comprehensive Model of Human Ear for Analysis of Implantable Hearing Devices

A finite element (FE) model of the human ear including the ear canal, middle ear, and spiral cochlea was constructed from histological sections of human temporal bone. Multiphysics analysis of the acoustics, structure, and fluid coupling in the ear was conducted in the model. The viscoelastic material behavior was applied to the middle ear soft tissues based on dynamic measurements of tissues in our laboratory. The FE model was first validated using the experimental data obtained in human cadaver ears, and then used to investigate the efficiency of the forward and reverse mechanical driving with middle ear implant, and the passive vibration of basilar membrane (BM) with cochlear implant placed in the cochlear scala tympani. The middle ear transfer function and the cochlear function of the BM vibration were derived from the model. This comprehensive ear model provides a novel computational tool to visualize and compute the implantable hearing devices and surgical procedures.