Effects of model definitions and parameter values in finite element modeling of human middle ear mechanics

Semi-automatic algorithm to build finite element numerical models of the human hearing system from Micro-CT data

Finite Element modeling has been an extended methodology to build numerical model to simulate the behavior of the hearing system. Due to the complexity of the system and the difficulties to reduce the uncertainties of the geometric data, they result in computationally expensive models, sometimes generic, representative of average geometries. It makes it difficult to validate the model with direct experimental data from the same specimen or to establish a patient-oriented modeling strategy. In the present paper, a first attempt to automatize the process of model building is made. The source information is geometrical information obtained from CT of the different elements that compose the system. Importing that data, we have designed the complete procedure to build a model including tympanic membrane, ossicular chain and cavities. The methodology includes the proper coupling of all the elements and the generation of the corresponding finite element model. The whole automatic procedure is not complete, as we need to make some human-assisted decisions; however, the model development time is reduced from 4 weeks to approximately 3 days. The goal of the modeling algorithm is to build a Finite Element Model with a limited computational cost. Several tasks as contour identification or model decimation are designed and integrated in order to follow a semi-automated process that allows generating a patient-oriented model.

Material characterization of human middle ear using machine-learning-based surrogate models

This study aims to introduce a novel non-invasive method for rapid material characterization of middle-ear structures, taking into consideration the invaluable insights provided by the mechanical properties of ear tissues. Valuable insights into various ear pathologies can be gleaned from the mechanical properties of ear tissues, yet conventional techniques for assessing these properties often entail invasive procedures that preclude their use on living patients. In this study, in the first step, we developed machine-learning models of the middle ear to predict its responses with a significantly lower computational cost in comparison to finite-element models. Leveraging findings from prior research, we focused on the most influential model parameters: the Young's modulus and thickness of the tympanic membrane and the Young's modulus of the stapedial annular ligament. The eXtreme Gradient Boosting (XGBoost) method was implemented for creating the machine-learning models. Subsequently, we combined the created machine-learning models with Bayesian optimization (BoTorch) for fast and efficient estimation of the Young's moduli of the tympanic membrane and the stapedial annular ligament. We demonstrate that the resultant surrogate models can fairly represent the vibrational responses of the umbo, stapes footplate, and vibration patterns of the tympanic membrane at most frequencies. Also, our proposed material characterization approach successfully estimated the Young's moduli of the tympanic membrane and stapedial annular ligament (separately and simultaneously) with values of mean absolute percentage error of less than 7%. The remarkable accuracy achieved through the proposed material characterization method underscores its potential for eventual clinical applications of estimating mechanical properties of the middle-ear structures for diagnostic purposes.

Performance of personalised prosthesis under static pressure: Numerical analysis and experimental validation

This study investigates the performance of personalised middle ear prostheses under static pressure through a combined approach of numerical analysis and experimental validation. The sound transmission performances of both normal and reconstructed middle ears undergo changes under high positive or negative pressure within the middle ear cavity. This pressure fluctuation has the potential to result in prosthesis displacement/extrusion in patients. To optimise the design of middle ear prostheses, it is crucial to consider various factors, including the condition of the middle ear cavity in which the prosthesis is placed. The integration of computational modelling techniques with non-invasive imaging modalities has demonstrated significant promise and distinct prospects in middle ear surgery. In this study, we assessed the efficacy of Finite Element (FE) analysis in modelling the responses of both normal and reconstructed middle ears to elevated static pressure within the ear canal. The FE model underwent validation using experimental data derived from human cadaveric temporal bones before progressing to subsequent investigations. Afterwards, we assessed stapes and umbo displacements in the reconstructed middle ear under static pressure, with either a columella-type prosthesis or a prosthetic incus, closely resembling a healthy incus. Results indicated the superior performance of the prosthetic incus in terms of both sound transmission to the inner ear and stress distribution patterns on the TM, potentially lowering the risk of prosthesis displacement/extrusion. This study underscores the potential of computational analysis in middle ear surgery, encompassing aspects such as prosthesis design, predicting outcomes in ossicular chain reconstruction (OCR), and mitigating experimental costs.

Finite element modelling of sound transmission in the Weberian apparatus of zebrafish (Danio rerio)

Zebrafish, an essential vertebrate model, has greatly expanded our understanding of hearing. However, one area that remains unexplored is the biomechanics of the Weberian apparatus, crucial for sound conduction and perception. Using micro-computed tomography (μCT) bioimaging, we created three-dimensional finite element models of the zebrafish Weberian ossicles. These models ranged from the exact size to scaled isometric versions with constrained geometry (1 to 10 mm in ossicular chain length). Harmonic finite element analysis of all 11 models revealed that the resonance frequency of the zebrafish's Weberian ossicular chain is approximately 900 Hz, matching their optimal hearing range. Interestingly, resonance frequency negatively correlated with size, while the ratio of peak displacement and difference of resonance frequency between tripus and scaphium remained constant. This suggests the transmission efficiency of the ossicular chain and the homogeneity of resonance frequency at both ends of the chain are not size-dependent. We conclude that the Weberian apparatus's resonance frequency can explain zebrafish's best hearing frequency, and their biomechanical characteristics are not influenced by isometric ontogeny. As the first biomechanical modelling of atympanic ear and among the few non-human ear modelling, this study provides a methodological framework for further investigations into hearing mechanisms and the hearing evolution of vertebrates.

Prestrain in the eardrum investigated using laser-ablation perforation: A proof of principle study on the New Zealand white rabbit

While the presence of residual stress (also called prestress) in the tympanic membrane (TM) was hypothesized more than 150 years ago by von Helmholtz (1869), little experimental data exists to date. In this paper, a novel approach to study residual stress is presented. Using a pulsed laser, the New Zealand white rabbit TM is perforated at seven predefined locations. The subsequent retraction of the membrane around the holes is computed using digital image correlation (DIC). The amount of retraction is the so-called prestrain, which is caused by the release of prestress due to the perforation. By measuring the prestrain using DIC, we show that residual stress is clearly present over the entire rabbit TM surface. In total, fourteen TMs have been measured in this work. An automated approach allows tracking the holes' deformation during the measurement process and enables a more robust analysis than was previously possible. We find similar strains (around 5%) as reported in previous work, in which slits were created manually using flattened surgical needles. However, the new approach greatly reduces measurement time, which minimizes dehydration artifacts. To investigate the effect of perforation location on the TM, the spatial decrease of the prestrain (α) around the perforation was quantified. Perforations inferior to the umbo showed the least negative α values, i.e., the most gradual decrease around the hole, and were the most consistent. Perforations on other locations showed more negative α values, i.e., steeper decrease in strain, but were less consistent across samples. We also investigated the effect of the holes' creation sequence but did not observe a significant change in the results. Overall, the presented method allows for consistent residual stress measurements over the TM surface. The findings contribute to our fundamental knowledge of the mechanics of the rabbit TM and provide a basis for future work on human TMs.

Relative importance and interactions of parameters of finite-element models of human middle ear

In the last decades, finite-element models of the middle ear have been widely used to predict the middle-ear vibration outputs. Even with the simplest linear assumption for material properties of the structures in the middle ear, these models need tens of parameters. Due to the complexities of measurements of material properties of these structures, accurate estimations of the values of most of these parameters are not possible. In this study, we benefited from the stochastic finite-element model of the middle ear we had developed in the past, to perform global sensitivity analysis. For this aim, we implemented Sobol' sensitivity analysis which ranks the importance of all uncertain parameters and interactions among them at different frequencies. To decrease the computational costs, we found Sobol' indices from surrogate models that we created using stochastic finite-element results and the polynomial chaos expansion method. Based on the results, the Young's modulus and thickness of the tympanic membrane, Young's modulus and damping of the stapedial annular ligaments, and the Young's modulus of ossicles are among the parameters with the greatest impacts on vibrations of the umbo and stapes footplate. Furthermore, the most significant interactions happen between the Young's modulus and thickness of the tympanic membrane.

Mechanical effect of reconstructed shapes of autologous ossicles on middle ear acoustic transmission

Conductive hearing loss is caused by a variety of defects, such as chronic otitis media, osteosclerosis, and malformation of the ossicles. In such cases, the defective bones of the middle ear are often surgically reconstructed using artificial ossicles to increase the hearing ability. However, in some cases, the surgical procedure does not result in increased hearing, especially in a difficult case, for example, when only the footplate of the stapes remains and all of the other bones are destroyed. Herein, the appropriate shapes of the reconstructed autologous ossicles, which are suitable for various types of middle-ear defects, can be determined by adopting an updating calculation based on a method that combines numerical prediction of the vibroacoustic transmission and optimization. In this study, the vibroacoustic transmission characteristics were calculated for bone models of the human middle ear by using the finite element method (FEM), after which Bayesian optimization (BO) was applied. The effect of the shape of artificial autologous ossicles on the acoustic transmission characteristics of the middle ear was investigated with the combined FEM and BO method. The results suggested that the volume of the artificial autologous ossicles especially has a great influence on the numerically obtained hearing levels.

Finite-Element Modelling Based on Optical Coherence Tomography and Corresponding X-ray MicroCT Data for Three Human Middle Ears

PURPOSE

Optical coherence tomography (OCT) is an emerging imaging modality which is non-invasive, can be employed in vivo, and can record both anatomy and vibrations. The purpose here is to explore the application of finite-element (FE) modelling to OCT data.

METHODS

We recorded vibrations for three human cadaver middle ears using OCT. We also have X-ray microCT images from the same ears. Three FE models were built based on geometries obtained from the microCT images. The material properties and boundary conditions of the models were obtained from previously reported studies.

RESULTS

Tympanic-membrane (TM) vibration patterns were computed for the three models and compared with the patterns measured using OCT. Frequency responses were also computed for all three models for several locations in the middle ear and compared with the OCT displacements and with the literature. The three models were compared with each other in terms of geometry and function. Parameter sensitivity analyses were done and the results were compared among the models and with the literature. The simulated TM displacement patterns are qualitatively similar to the OCT results. The simulated displacements are closer to the OCT results for 500 Hz and 1 kHz but the differences are greater at 2 kHz.

CONCLUSION

This study provides an initial look at the combined use of OCT measurements and FE modelling based on subject-specific anatomy. The geometries and parameters of the existing FE models could be modified for individual patients in the future to help identify abnormalities in the middle ear.

Inaccuracies of deterministic finite-element models of human middle ear revealed by stochastic modelling

For over 40 years, finite-element models of the mechanics of the middle ear have been mostly deterministic in nature. Deterministic models do not take into account the effects of inter-individual variabilities on middle-ear parameters. We present a stochastic finite-element model of the human middle ear that uses variability in the model parameters to investigate the uncertainty in the model outputs (umbo, stapes, and tympanic-membrane displacements). We demonstrate: (1) uncertainties in the model parameters can be magnified by more than three times in the umbo and stapes footplate responses at frequencies above 2 kHz; (2) middle-ear models are biased and they distort the output distributions; and (3) with increased frequency, the highly-uncertain regions spatially spread out on the tympanic membrane surface. Our results assert that we should be mindful when using deterministic finite-element middle-ear models for critical tasks such as novel device developments and diagnosis.

Characterization of middle ear soft tissue damping and its role in sound transmission

Damping plays an important role in the middle ear (ME) sound transmission system. However, how to mechanically characterize the damping of ME soft tissues and the role of damping in ME sound transmission have not yet reached a consensus. In this paper, a finite element (FE) model of the partial external and ME of the human ear, considering both Rayleigh damping and viscoelastic damping for different soft tissues, is developed to quantitatively investigate the damping in soft tissues effects on the wide-frequency response of the ME sound transmission system. The model-derived results can capture the high-frequency (above 2 kHz) fluctuations and obtain the 0.9 kHz resonant frequency (RF) of the stapes velocity transfer function (SVTF) response. The results show that the damping of pars tensa (PT), stapedial annular ligament (SAL) and incudostapedial joints (ISJ) can help smooth the broadband response of the umbo and stapes footplate (SFP). It is found that, between 1 and 8 kHz, the damping of the PT increases the magnitude and phase delay of the SVTF above 2 kHz while the damping of the ISJ can avoid excessive phase delay of the SVTF, which is important in maintaining the synchronization in high-frequency vibration but has not been revealed before. Below 1 kHz, the damping of the SAL plays a more important role, and it can decrease the magnitude but increases the phase delay of the SVTF. This study has implications for a better understanding of the mechanism of ME sound transmission.

Finite element modelling of the human middle ear using synchrotron-radiation phase-contrast imaging

Finite element (FE) models of the middle ear often lack accurate geometry of soft tissue structures, such as the suspensory ligaments, as they can be difficult to discern using conventional imaging modalities, such as computed tomography. Synchrotron-radiation phase-contrast imaging (SR-PCI) is a non-destructive imaging modality that has been shown to produce excellent visualization of soft tissue structures without the need for extensive sample preparation. The objectives of the investigation were to firstly use SR-PCI to create and evaluate a biomechanical FE model of the human middle ear that includes all soft tissue structures, and secondly, to investigate how modelling assumptions and simplifications of ligament representations affect the simulated biomechanical response of the FE model. The FE model included the suspensory ligaments, ossicular chain, tympanic membrane, the incudostapedial and incudomalleal joints, and the ear canal. Frequency responses obtained from the SR-PCI-based FE model agreed well with published laser doppler vibrometer measurements on cadaveric samples. Revised models with exclusion of the superior malleal ligament (SML), simplification of the SML, and modification of the stapedial annular ligament were studied, as these revised models represented modelling assumptions that have been made in literature.

Rabbit tympanic membrane thickness distribution obtained via optical coherence tomography

Knowing the precise tympanic membrane (TM) thickness variation is crucial in understanding the functional properties of the TM and has a significant effect on the accuracy of computational models. Using optical coherence tomography, we imaged five left and five right TMs of domestic New Zealand rabbits. From these data, ten thickness distribution maps were computed. Although inter-specimen variability is present, similar features could be observed in all samples: The rabbit TM is thickest around the umbo, with values of 150 ± 32 µm. From the umbo towards the TM annulus, the thickness gradually decreases down to 38 ± 7 µm around the midway location, but increases up to 54 ± 19 µm at the TM annulus. The thickness values at the umbo are comparable to literature data for humans, but the rabbit TM is thinner at the TM annulus and in-between the umbo and annulus. Moreover, the rabbit TM thickness distribution is highly symmetrical, which is not the case for the human TM. The results improve our general understanding of TM structure in rabbits and may improve numerical models of TM dynamical behavior.

Mammalian middle ear mechanics: A review

The middle ear is part of the ear in all terrestrial vertebrates. It provides an interface between two media, air and fluid. How does it work? In mammals, the middle ear is traditionally described as increasing gain due to Helmholtz’s hydraulic analogy and the lever action of the malleus-incus complex: in effect, an impedance transformer. The conical shape of the eardrum and a frequency-dependent synovial joint function for the ossicles suggest a greater complexity of function than the traditional view. Here we review acoustico-mechanical measurements of middle ear function and the development of middle ear models based on these measurements. We observe that an impedance-matching mechanism (reducing reflection) rather than an impedance transformer (providing gain) best explains experimental findings. We conclude by considering some outstanding questions about middle ear function, recognizing that we are still learning how the middle ear works.

Exploring Mechanisms Underlying Unexplained Air-Bone Gaps Post-Myringoplasty: Temporal Bone Model and Finite Element Analysis

PURPOSE

Air-bone gap (ABG) is an essential indicator of middle ear transfer function after myringoplasty. However, there is still uncertainty about the mechanisms behind unexplained ABGs in patients post-myringoplasty. The present study investigated these mechanisms using cadaveric temporal bone (TB) measurement and finite element (FE) modeling.

METHODS

Three conditions of tympanic membrane (TM) perforation were modeled with a perforated area of 6%, 24%, and 50% of the total TM area to simulate a small, medium, or large TM perforation of TB model. A piece of paper was used to patch the TM perforation to simulate the situation post-myringoplasty. In the FE model for post-operation, the material properties at the perforation area were changed. Measurement of TM vibration at the umbo was undertaken with a laser Doppler vibrometer (LDV).

RESULTS

As the perforated area increased vibration of the TM at the umbo decreased in both the TB and FE models. But the reduction of TM vibration is more minor in the FE model than in the TB model. After the perforation was repaired, the displacement of TM at the umbo could not be recovered totally in the TB and FE models. In the FE model, the displacement of TM at the umbo decreased markedly when the cone shape of TM flattened, and the reduction was almost the same as that in the TB model in the condition of large perforation.

CONCLUSION

The material properties and the anatomical shape of the repaired TM could influence the TM's modal motion and wave motion. Except for appearance and shape current clinical instruments are unable to resolve factors that affect TM motion. Consequently the ABG seen post-myringoplasty remains unexplained.

Development of electrospun, biomimetic tympanic membrane implants with tunable mechanical and oscillatory properties for myringoplasty

Most commonly, autologous grafts are used in tympanic membrane (TM) reconstruction. However, apart from the limited availability and the increased surgical risk, they cannot replicate the full functionality of the human TM properly. Hence, biomimetic synthetic TM implants have been developed in our project to overcome these drawbacks. These innovative TM implants are made from synthetic biopolymer polycaprolactone (PCL) and silk fibroin (SF) by electrospinning technology. Static and dynamic experiments have shown that the mechanical and oscillatory behavior of the TM implants can be tuned by adjusting the solution concentration, the SF and PCL mixing ratio and the electrospinning parameters. In addition, candidates for TM implants could have comparable acousto-mechanical properties to human TMs. Finally, these candidates were further validated in in vitro experiments by performing TM reconstruction in human cadaver temporal bones. The reconstructed TM with SF-PCL blend membranes fully recovered the acoustic vibration when the perforation was smaller than 50%. Furthermore, the handling, medium adhesion and transparency of the developed TM implants were similar to those of human TMs.

Morphology of human ear canal and its effect on sound transmission

The ear canal (EC) is essential for sound transfer and crucial for hearing. Some pathological conditions may modify its morphology, leading to EC sound pressure redistribution, and stapes footplate displacement (FPD) gain alteration. However, no consensus regarding pathological EC and its impact on sound transfer has yet been achieved. To address the effect of morphology of EC on sound pressure redistribution and FPD gain. Varied pathological EC finite element (FE) models were constructed and analyzed based on FE analysis. The results indicated that canal wall down mastoidectomy decreases the second resonance frequency of the EC. The canal wall down mastoidectomy, with conchaplasty increased the first resonance frequency, but decreased the second along with the interval sound pressure gain increased, following which the FPD gain was altered. Stenosis of the EC at the internal portion decreased the second resonance frequency with minimal effect to the first part. When the stenosis moved to the outer portion of the EC, the first resonance frequency decreased, and the second one increased, along with the interval sound pressure gain decreased and FPD gain. Finally, the simplified EC model exerted a minimal effect on sound transfer. The minimal change in EC, such as simplification, straightening, canal wall down mastoidectomy, or enlargement, moderately affects the sound transfer; however, the EC stenosis deteriorates the sound transfer remarkably.

Material Identification on Thin Shells Using the Virtual Fields Method, Demonstrated on the Human Eardrum

Characterization of material parameters from experimental data remains challenging, especially on biological structures. One of such techniques allowing for the inverse determination of material parameters from measurement data is the virtual fields method (VFM). However, application of the VFM on general structures of complicated shape has not yet been extensively investigated. In this paper, we extend the framework of the VFM method to thin curved solids in three-dimensional, commonly denoted shells. Our method is then used to estimate the Young's modulus and hysteretic damping of the human eardrum. By utilizing Kirchhoff plate theory, we assume that the behavior of the shell varies linearly through the thickness. The total strain of the shell can then be separated in a bending and membrane strain. This in turn allowed for an application of the VFM based only on data of the outer surface of the shell. We validated our method on simulated and experimental data of a human eardrum made to vibrate at certain frequencies. It was shown that the identified material properties were accurately determined based only on data from the outer surface and are in agreement with literature. Additionally, we observed that neither the bending nor the membrane strain in an human eardrum can be neglected and both contribute significantly to the total strain found experimentally.

On the hearing effects of a cholesteatoma growing: A biomechanical study

Chronic otitis media enables the appearance of a benign middle ear tumor, known as a cholesteatoma, that may compromise hearing. To evaluate the influence of a cholesteatoma growth on the hearing function, a computational middle ear model based on the finite element method was used and three different size of cholesteatoma were modeled. The cholesteatoma solidification and the consequent degradation of the ossicles were also simulated as two condition that commonly occurs during cholesteatoma evolution. A sound pressure level of 80 dB SPL was applied in the tympanic membrane and a steady state analysis was performed for frequencies from 100 Hz to 10 kHz. The displacements of both the tympanic membrane and the stapes footplate were measured. The results were compared with a healthy case and it was shown that the cholesteatoma development leads to a decrease in the umbo and stapes displacements. The ossicles degradation simulation showed the higher difference comparing with the cholesteatoma in an initial stage, with lower displacements in the stapes footplate mainly for high frequencies. The observed displacement differences are directly connected to hearing loss, being possible to conclude that cholesteatoma evolution in the middle ear will lead to hearing problems, mainly in an advanced stage.

Cellulose-Based Fibrous Materials From Bacteria to Repair Tympanic Membrane Perforations

Perforation is the most common illness of the tympanic membrane (TM), which is commonly treated with surgical procedures. The success rate of the treatment could be improved by novel bioengineering approaches. In fact, a successful restoration of a damaged TM needs a supporting biomaterial or scaffold able to meet mechano-acoustic properties similar to those of the native TM, along with optimal biocompatibility. Traditionally, a large number of biological-based materials, including paper, silk, Gelfoam®, hyaluronic acid, collagen, and chitosan, have been used for TM repair. A novel biopolymer with promising features for tissue engineering applications is cellulose. It is a highly biocompatible, mechanically and chemically strong polysaccharide, abundant in the environment, with the ability to promote cellular growth and differentiation. Bacterial cellulose (BC), in particular, is produced by microorganisms as a nanofibrous three-dimensional structure of highly pure cellulose, which has thus become a popular graft material for wound healing due to a number of remarkable properties, such as water retention, elasticity, mechanical strength, thermal stability, and transparency. This review paper provides a comprehensive overview of the current experimental studies of BC, focusing on the application of BC patches in the treatment of TM perforations. In addition, computational approaches to model cellulose and TM are summarized, with the aim to synergize the available tools toward the best design and exploitation of BC patches and scaffolds for TM repair and regeneration.

Structural stiffening in the human middle ear due to static pressure: Finite-element analysis of combined static and dynamic middle-ear behavior

The vibration response of the middle ear (ME) to sound changes when static pressure gradients are applied across the tympanic membrane (TM). To date, it has not been well understood which mechanisms lead to these changes in ME vibration response. In this study, a 3D finite-element model of the human ME was developed that simulates the sound-induced ME vibration response when positive and negative static pressures of up to 4 kPa are applied to the TM. Hyperelasticity of the soft-tissue components was considered to simulate large deformations under static pressure. Some ME components were treated as viscoelastic materials to capture the difference between their static and dynamic stiffness, which was needed to replicate both static and dynamic ME behavior. The change in dynamic stiffness with static preload was simulated by linearization of the hyperelastic constitutive model around the predeformed state. For the preloaded harmonic response, we found that the statically deformed ME geometry introduced asymmetry in the vibration loss between positive and negative pressure, which was due to the TM cone shape. As opposed to previous assumptions, the prestress in the ME due to static pressure had a substantial impact on the vibration response. We also found that material nonlinearity led to a higher stiffening at the umbo but a less pronounced stiffening at the footplate compared to the linear elastic condition. The results suggest that flexibility of the incudomalleolar joint (IMJ) enhances the decoupling of static umbo and footplate displacements, and that viscosity and viscoelasticity of the IMJ could play a role in the transfer of sound-induced vibrations from the umbo to the footplate. The components of the incudostapedial joint had minimal effect on ME mechanical behavior.

[Numerical analysis of the effect of typical middle ear disease on bone conduction]

Objective:To study the effects of typical middle ear lesions on bone conduction, and to provide a theoretical basis for the optimized design of bone conduction hearing aids. Method:A finite element model including the middle ear and cochlea was established by using micro-CT scanning date from a case of adult male human-cadaver temporal bone. The reliability of the model was verified by comparison with relevant experimental data. Based on this model, the effects of several typical middle ear lesions on bone conduction were analyzed. Result:With the increase of tympanic membrane mass, the basilar membrane response at the frequency range of 0.25 kHz to 1 kHz was reduced, and the maximum decrease was 6.9 dB at the frequency of 0.5 kHz. The incus fixation reduced basilar membrane response at the frequency range of 0.2 kHz to 1 kHz, with the maximum decrease of 12.0 dB at 0.75 kHz. The ossicular chain disruption reduced the response of basilar membrane at the frequency range of 1 kHz to 4 kHz, which the decrease was 5.0 dB at 1.5 kHz in the case of incudomallear joint disruption and 11.0 dB at 2 kHz in the case of incudostapedial joint disruption. Conclusion:The results showed that the increased mass of the tympanic membrane and incus fixation mainly reduced the response of bone conduction at the low frequencies. The ossicular chain disruption mainly reduced the response of bone conduction at the middle frequencies. Moreover, the deterioration of incudostapedial disruption was more obvious.

Evaluating fidelity of CT based 3D models for Zebrafish conductive hearing system

The zebrafish Weberian apparatus is an emerging model for human conductive hearing system. Their Weberian apparatus comprises minute bones and ligamentary links, and conducts sound pressure transmission from the gas bladder to inner ear through four pairs of Weberian ossicles along the vertebral column. We herein present a methodological study using MicroCT to image the Weberian apparatus for biomechanical and morphological analysis. The aim of this work is to evaluate computational models generated from multiple MicroCT scans with different parameters, to identify the most feasible scan combination for practical (minimized scan time) yet accurate (relative to highest resolution) biomechanical simulations. We segmented and created 3D models from CT scan image stacks at 4.64 μm, 5.05 μm, 9.30 μm and 13.08 μm voxel resolutions, respectively. Then, we used geometric morphometrics analysis to quantify inter-model shape differences, as well as a series of finite element modal and harmonic analyses to simulate auditory signal vibrations. Relative to the highest resolution and most accurate model, the Model 9.30 is closest in overall geometry and biomechanical behavior of all lower resolution models. The differences in resolution and quality of the CT substantially affect the segmentation and reconstruction process of the three-dimensional model of the ossicles, and the subsequent analyses. We conclude that scan voxel resolution is a key factor influencing outcomes of biomechanical simulations of delicate and minute structures, especially when studying the harmonic response of minute ossicles connected by ligaments using finite element modeling. Furthermore, contrast variations in CT images as determined by x-ray power and scan speed, also affect fidelity in 3D models and simulation outcomes.

De novo topology optimization of total ossicular replacement prostheses

Conductive hearing loss, due to middle ear pathologies or traumas, affects more than 5% of the population worldwide. Passive prostheses to replace the ossicular chain mainly rely on piston-like titanium and/or hydroxyapatite devices, which in the long term suffer from extrusion. Although the basic shape of such devices always consists of a base for contact with the eardrum and a stem to have mechanical connection with the residual bony structures, a plethora of topologies have been proposed, mainly to help surgical positioning. In this work, we optimize the topology of a total ossicular replacement prosthesis, by maximizing the global stiffness and under the smallest possible volume constraint that ensures material continuity. This investigation optimizes the prosthesis topology in response to static displacement loads with amplitudes that normally occur during sound stimulation in a frequency range between 100 Hz and 10 kHz. Following earlier studies, we discuss how the presence and arrangement of holes on the surface of the prosthesis plate in contact with the umbo affect the overall geometry. Finally, we validate the designs through a finite-element model, in which we assess the prosthesis performance upon dynamic sound pressure loads by considering four different constitutive materials: titanium, cortical bone, silk, and collagen/hydroxyapatite. The results show that the selected prostheses present, almost independently of their constitutive material, a vibroacustic behavior close to that of the native ossicular chain, with a slight almost constant positive shift that reaches a maximum of ≈5 dB close to 1 kHz. This work represents a reference for the development of a new generation of middle ear prostheses with non-conventional topologies for fabrication via additive manufacturing technologies or ultraprecision machining in order to create patient-specific devices to recover from conductive hearing loss.

The Influence of Piezoelectric Transducer Stimulating Sites on the Performance of Implantable Middle Ear Hearing Devices: A Numerical Analysis

To overcome the inherent deficiencies of hearing aids, implantable middle ear hearing devices (IMEHDs) have emerged as a new treatment for hearing loss. However, clinical results show that the IMEHD performance varies with its transducer’s stimulating site. To numerically analyze the influence of the piezoelectric transducer’s stimulating sites on its hearing compensation performance, we constructed a human ear finite element model and confirmed its validity. Based on this finite element model, the displacement stimulation, which simulates the piezoelectric transducer’s stimulation, was applied to the umbo, the incus long process, the incus body, the stapes, and the round window membrane, respectively. Then, the stimulating site’s effect of the piezoelectric transducer was analyzed by comparing the corresponding displacements of the basilar membrane. Besides, the stimulating site’s sensitivity to the direction of excitation was also studied. The result of the finite element analysis shows that stimulating the incus body is least efficient for the piezoelectric transducer. Meanwhile, stimulating the round window membrane or the stapes generates a higher basilar membrane displacement than stimulating the eardrum or the incus long process. However, the performance of these two ideal sites’ stimulation is sensitive to the changes in the excitation’s direction. Thus, the round window membrane and the stapes is the ideal stimulating sites for the piezoelectric transducer regarding the driving efficiency. The direction of the excitation should be guaranteed for these ideal sites.

Strain distribution in rabbit eardrums under static pressure

Full-field strain maps of intact rabbit eardrums subjected to static pressures are presented. A stochastic intensity pattern was applied to 12 eardrums, and strain maps were measured at the medial site using a stereoscopic digital image correlation setup for pressures between -2 and 2 kPa. Ear canal overpressures induced circumferential orientated positive strains between manubrium and the eardrum border that increased almost linearly with pressure. Radially orientated negative strains were found at the border and manubrium. Ear canal underpressures caused negative circumferential strains between manubrium and the tympanic annulus but radially orientated positive strains at the borders. The magnitudes of these negative strains at underpressures were larger than those of positive strains at overpressures and were nonlinearly proportional to pressure. In three ears, strains were calculated with intact and removed cochlea. The effect of cochlea removal on the peak-to-peak strain was found to be no more than 3%.

How flexibility and eardrum cone shape affect sound conduction in single-ossicle ears: a dynamic model study of the chicken middle ear

It is believed that non-mammals have poor hearing at high frequencies because the sound-conduction performance of their single-ossicle middle ears declines above a certain frequency. To better understand this behavior, a dynamic three-dimensional finite-element model of the chicken middle ear was constructed. The effect of changing the flexibility of the cartilaginous extracolumella on middle-ear sound conduction was simulated from 0.125 to 8 kHz, and the influence of the outward-bulging cone shape of the eardrum was studied by altering the depth and orientation of the eardrum cone in the model. It was found that extracolumella flexibility increases the middle-ear pressure gain at low frequencies due to an enhancement of eardrum motion, but it decreases the pressure gain at high frequencies as the bony columella becomes more resistant to extracolumella movement. Similar to the inward-pointing cone shape of the mammalian eardrum, it was shown that the outward-pointing cone shape of the chicken eardrum enhances the middle-ear pressure gain compared to a flat eardrum shape. When the outward-pointing eardrum was replaced by an inward-pointing eardrum, the pressure gain decreased slightly over the entire frequency range. This decrease was assigned to an increase in bending behavior of the extracolumella and a reduction in piston-like columella motion in the model with an inward-pointing eardrum. Possibly, the single-ossicle middle ear of birds favors an outward-pointing eardrum over an inward-pointing one as it preserves a straight angle between the columella and extrastapedius and a right angle between the columella and suprastapedius, which provides the optimal transmission.

How Stapes Ankylosis and Fracture Affect Middle Ear Dynamics: A Numerical Study

Numerical models of the human middle ear have been developed throughout the last 30 years, for different purposes. While several types of pathologies have been studied, stapedial disorders were seldomly explored. This papers aims to clarify how stapes fracture and some forms of stapes ankylosis, such as stapedial tendon (ST) ossification, augmented pyramidal eminence (PE) and bony bar presence, affect the sound transmission through the middle ear. In addition, the stapes dynamics is also analyzed by means of total displacement and first principal strain. For the purpose of the study, first, a three-dimensional finite element model of the human middle ear is detailed and validated under normal (healthy) conditions. The model is then modified to represent the stapedial disorders of interest. A measure is established for evaluating how the disorders reduce sound transmission through the middle ear. Results of the reduction of sound transmission showed that the different forms of stapes ankylosis affect primarily low frequencies, while the stapes fracture mostly affects high frequency sound transmission. According to the results, an augmented PE does not restrict stapes movement unless followed by some ossification of the ST. In addition, the question whether the fracture is in the anterior or posterior crus and the distance of the fractured part from the stapes footplate have a relevant role in the reduction of the sound transmission. Finally, the analysis of total displacement and first principal strain of the stapes helped to highlight some differences among the stapedial disorders.

Influence of ossicular chain malformation on the performance of round-window stimulation: A finite element approach

As a novel application of implantable middle ear hearing device, round-window stimulation is widely used to treat hearing loss with middle ear disease, such as ossicular chain malformation. To evaluate the influence of ossicular chain malformations on the efficiency of the round-window stimulation, a human ear finite element model, which incorporates cochlear asymmetric structure, was constructed. Five groups of comparison with experimental data confirmed the model’s validity. Based on this model, we investigated the influence of three categories of ossicular chain malformations, that is, incudostapedial disconnection, incus and malleus fixation, and fixation of the stapes. These malformations’ effects were evaluated by comparing the equivalent sound pressures derived from the basilar membrane displacement. Results show that the studied ossicular chain malformations mainly affected the round-window simulation’s performance at low frequencies. In contrast to the fixation of the ossicles, which mainly deteriorates round-window simulation’s low-frequency performance, incudostapedial disconnection increases this performance, especially in the absence of incus process and stapes superstructure. Among the studied ossicular chain malformations, the stapes fixation has a much more severe impact on the round-window stimulation’s efficiency. Thus, the influence of the patients’ ossicular chain malformations should be considered in the design of the round-window stimulation’s actuator. The low-frequency output of the round-window simulation’s actuator should be enhanced, especially for treating the patients with stapes fixation.

A model and experimental approach to the middle ear transfer function related to hearing in the humpback whale (Megaptera novaeangliae)

At present, there are no direct measures of hearing for any baleen whale (Mysticeti). The most viable alternative to in vivo approaches to simulate the audiogram is through modeling outer, middle, and inner ear functions based on the anatomy and material properties of each component. This paper describes a finite element model of the middle ear for the humpback whale (Megaptera novaeangliae) to calculate the middle ear transfer function (METF) to determine acoustic energy transmission to the cochlea. The model was developed based on high resolution computed tomography imaging and direct anatomical measurements of the middle ear components for this mysticete species. Mechanical properties for the middle ear tissues were determined from experimental measurements and published values. The METF for the humpback whale predicted a better frequency range between approximately 15 Hz and 3 kHz or between 200 Hz and 9 kHz based on two potential stimulation locations. Experimental measures of the ossicular chain, tympanic membrane, and tympanic bone velocities showed frequency response characteristics consistent with the model. The predicted best sensitivity hearing ranges match well with known vocalizations of this species.

The effects of varying tympanic-membrane material properties on human middle-ear sound transmission in a three-dimensional finite-element model

An anatomically based three-dimensional finite-element human middle-ear (ME) model is used to test the sensitivity of ME sound transmission to tympanic-membrane (TM) material properties. The baseline properties produce responses comparable to published measurements of ear-canal input impedance and power reflectance, stapes velocity normalized by ear-canal pressure (P_EC), and middle-ear pressure gain (MEG), i.e., cochlear-vestibule pressure (P_V) normalized by P_EC. The mass, Young's modulus (E_TM), and shear modulus (G_TM) of the TM are varied, independently and in combination, over a wide range of values, with soft and bony TM-annulus boundary conditions. MEG is recomputed and plotted for each case, along with summaries of the magnitude and group-delay deviations from the baseline over low (below 0.75 kHz), mid (0.75-5 kHz), and high (above 5 kHz) frequencies. The MEG magnitude varies inversely with increasing TM mass at high frequencies. Increasing E_TM boosts high frequencies and attenuates low and mid frequencies, especially with a bony TM annulus and when G_TM varies in proportion to E_TM, as for an isotropic material. Increasing G_TM on its own attenuates low and mid frequencies and boosts high frequencies. The sensitivity of MEG to TM material properties has implications for model development and the interpretation of experimental observations.