Study of the middle ear behavior considering the presence of fluid



==Resumo==

O objetivo deste trabalho foi estudar a influência de diferentes fluidos dentro e fora da cavidade timpânica no comportamento biomecânico do ouvido médio. O modelo computacional 3 D da cadeia tímpano‐ossicular foi construído através de imagens de tomografia axial computorizada (TAC). Após a construção do modelo da cadeia tímpano‐ossicular, para a parte interna da membrana timpânica foi criada a caixa timpânica, que inclui toda a cadeia ossicular, e para a parte externa da membrana foi criada uma passagem, simulando o canal auditivo externo.

A discretização do modelo foi feita utilizando o método dos elementos finitos, com base no software ABAQUS. As propriedades mecânicas foram extraídas de trabalhos anteriores. Foram obtidos os deslocamentos do umbo e de um ponto central da platina do estribo, para um nível de pressão sonora de 105 dB aplicada na membrana timpânica, comparando a cadeia tímpano‐ossicular e os diferentes níveis de pressão externa e interna do ouvido médio (fluido líquido e fluido gasoso). Os resultados demonstram que o comportamento biomecânico do ouvido médio é diferente para os diferentes níveis de pressões testadas neste modelo. Quando utilizada uma pressão de fluido líquido, os deslocamentos são de amplitude inferior quando comparados com os resultados obtidos com pressão de ar.

==Abstract==

The objective of this work is to study the influence of different fluids of the tympanic cavity into the biomechanical behavior of the middle ear. The 3 D computer model of the tympanic ossicular chain was built using images from computed tomography (CT). After construction of the tympanic ossicular chain model, it was adapted the simulation by placing the passage to the external auditory canal and to the tympanic cavity.

The discretization of the model was made using the finite element method based on ABAQUS software. The mechanical properties were extracted from previous works. The umbo and the stapes footplate displacements, for a sound pressure level of 105 dB applied in tympanic membrane were obtained, comparing the tympanic ossicular chain and the different levels of external and internal pressure in the middle ear (liquid and air). The results demonstrate that the biomechanical behavior of the middle ear is different for different levels of pressure tested in this model. When using liquid, displacements are of lower amplitude when compared with the results obtained with air

==Palavras‐chave==

Ouvido médio ;                             Biomecânica ;                             Pressões ;                             Método dos elementos finitos

==Keywords==

Middle ear ;                             Biomechanics ;                             Pressures ;                             Finite Element Method

==1. Introdução==

A cavidade timpânica (CT) é um espaço repleto de ar, incluída no osso temporal, que se liga ao ouvido externo pela membrana timpânica e ao ouvido interno pela janela oval. É constituída por 6 paredes (inferior, superior, anterior, posterior, lateral e medial). A CT está ligada à cavidade nasal anteriormente pela trompa de Eustáquio e posteriormente e superiormente ao antro mastoide, que são cavidades preenchidas por ar, funcionando como um reservatório de ar. Possui em seu interior a cadeia ossicular (martelo, bigorna e estribo), ligamentos e músculos e uma porção do nervo facial (corda do tímpano). Seu papel mais importante é a equalização das variações de pressão atmosférica e a do ar contido na CT.

O método dos elementos finitos é atualmente o método mais amplamente utilizado para o cálculo de complexos sistemas mecânicos, tais como o ouvido humano.

Alguns estudos do ouvido têm sido apresentados na literatura, sendo que o primeiro modelo usando elementos finitos foi feito a partir do ouvido de um gato no ano de 1978 [[#bib0110|[1]]] . Entre 1995‐1997 foram efetuados estudos experimentais nesta área para obtenção dos deslocamentos do umbo e da platina do estribo [[#bib0115|[2]]] , [[#bib0120|[3]]]  and [[#bib0125|[4]]] . Entre 1999‐2002 surgiram outros modelos computacionais que também determinavam os deslocamentos do umbo e da platina do estribo [[#bib0130|[5]]] , [[#bib0135|[6]]]  and [[#bib0140|[7]]] , seguidos de outros estudos em que comparavam os próprios resultados com os resultados anteriores [[#bib0145|[8]]] , [[#bib0150|[9]]] , [[#bib0155|[10]]] , [[#bib0160|[11]]] , [[#bib0165|[12]]]  and [[#bib0170|[13]]] .      

Estudar a influência de diferentes fluidos da CT no comportamento biomecânico do ouvido médio humano, através da utilização do método dos elementos finitos, tornou‐se de grande importância, sabendo‐se que o mau funcionamento do equilíbrio extra e intratimpânico pode causar grande desconforto aos seres humanos. Exemplo disso são as bolsas de retração, pois sabe‐se que estas são as mais importantes manifestações clínicas da diminuição da pressão no interior da CT [[#bib0175|[14]]] . A otite média é uma patologia muito comum que apresenta fluído líquido dentro da CT e que está geralmente associada a hipoacusia por acumulação dos exsudados, que podem ser serosos, mucosos, purulentos e mucopurulentos [[#bib0175|[14]]] .      

Outro exemplo de desequilíbrio em que esta pressão pode ocorrer relaciona‐se com o mergulho em profundidade. Se o mergulhador não estiver atento à compensação da pressão na CT, pode surgir como consequência hemorragia e trauma no canal auditivo externo (CAE) e membrana timpânica. Em alguns casos evoluem para uma lesão otológica provocada por barotrauma.

O objetivo principal deste trabalho foi estudar a influência de diferentes fluidos aplicados na CT e CAE no comportamento biomecânico da cadeia tímpano‐ossicular, recorrendo a um modelo de simulação numérica baseado no método dos elementos finitos. Tendo esta perceção, será mais fácil simular algumas das patologias que atingem o ouvido e comparar com o ouvido normal.

==2. Material e métodos==

O primeiro passo deste trabalho foi a construção do modelo da cadeia tímpano‐ossicular do ouvido médio (membrana timpânica e os 3 ossículos, martelo, bigorna e estribo), incluindo ligamentos (superior, lateral e anterior do martelo, superior e posterior da bigorna, ligamento anular do estribo) e 2 músculos (tensor do tímpano e estapédico) [[#bib0165|[12]]] . O modelo foi adaptado [[#bib0170|[13]]]  baseado no dimensionamento descrito na obra de Paço et al''.''[[#bib0175|[14]]] . Na [[#tbl0005|tabela 1]]  é possível verificar o número de nós e elementos utilizados na membrana timpânica e ossículos. O software utilizado foi o ABAQUS/CAE versão 6.5 [[#bib0180|[15]]] . Foi feita uma análise dinâmica em regime estacionário, utilizando o comando «*steady state dynamics» do ABAQUS.

<span id='tbl0005'></span>

Tabela 1.

Número de elementos e nós dos diferentes componentes

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"

|-

! Componentes:
! Elementos
! Nós
|-

| Membrana timpânica
| 11.165
| 15.295
|-

| Martelo
| 18.841
| 3.932
|-

| Bigorna
| 39.228
| 8.373
|-

| Estribo
| 9.218
| 2.840
|}

Baseado no trabalho de Sun et al. [[#bib0140|[7]]]  e no trabalho de Prendergast et al. [[#bib0130|[5]]] , foram estipuladas as propriedades de material para o ouvido médio. O coeficiente de Poisson foi assumido como 0,3 para todos os materiais. A membrana timpânica foi discretizada com a utilização de elementos tridimensionais hexaédricos de 8 nós, C3D8, sendo esta dividida em 3 camadas conforme a sua anatomia [[#bib0170|[13]]]  and [[#bib0175|[14]]] . A ''pars flaccida''  é considerada elástica isotrópica. Já para a ''pars tensa''  foram estipuladas propriedades diferentes para cada camada, sendo a camada externa e a camada interna considerada isotrópica e a camada central ortotrópica. Para «E», o módulo de Young, o índice ''θ''  indica direção tangencial e ''r''  a direção radial ( [[#tbl0010|tabela 2]] ). Para a discretização dos ossículos optou‐se por utilizar elementos tetraédricos, do tipo C3D4, tendo sido estabelecido um valor de 1,41 E + 10 Pa para o módulo de Young. A massa específica varia de acordo com as partes constituintes dos próprios ossículos ([[#tbl0010|tabela 2]] ).

<span id='tbl0010'></span>

{| class="wikitable" style="min-width: 60%;margin-left: auto; margin-right: auto;"
|+

Tabela 2.

Propriedades de material da membrana timpânica e ossículos

|-

! Propriedades Material
! Massa específica (Kg/m<sup>3</sup> )                                                    
! Módulo de Young (N/m<sup>2</sup> )                                                    
! 
|-

| Membrana timpânica
| 
| 1,20×103
| 
|-

| 
| Pars tensa
| 
| 
|-

| 
| Camada 1
| 
| 1,00×107
|-

| 
| Camada 2
| 
| Eθ=2,00×107; Er=3,20×107
|-

| 
| Camada 3
| 
| 1,00×107
|-

| 
| Pars flaccida
| 
| 1,00×107
|-

| Ossículos
| MarteloCabeçaPescoçoCabo
| 2,55×1034,53×1033,70×103
| 1,41×1010
|-

| 
| BigornaCorpoCurta apófiseLonga apófise
| 2,36×1032,26×1035,08×103
| 
|-

| 
| Estribo
| 2,20×103
| 
|}

O CAE e a CT foram construídos através da própria membrana timpânica. Foram criados elementos de fluido junto à membrana timpânica, tanto na parte externa como na parte interna, simulando, respetivamente, a passagem para o CAE ([[#fig0005|fig. 1]] ) e para a CT ([[#fig0010|fig. 2]] ).

<span id='fig0005'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_861468639-1-s2.0-S0213131515000334-gr1.jpg|center|357px|Simulação do CAE.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 1.

Simulação do CAE.

</span>
|}

<span id='fig0010'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_861468639-1-s2.0-S0213131515000334-gr2.jpg|center|357px|Vista da CT com a membrana timpânica à transparência, evidenciando os ossículos ...]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 2.

Vista da CT com a membrana timpânica à transparência, evidenciando os ossículos no seu interior.

</span>
|}

Os nós das extremidades livres dos ligamentos e dos músculos foram completamente fixos. Para esta simulação foram usados 7.645 elementos de fluido do tipo F3D4.

Para a passagem para o CAE e para a CT foram efetuadas 5 simulações diferentes de acordo com as propriedades estipuladas para os fluidos destas regiões:
* Só a cadeia tímpano‐ossicular.
* CAE com a presença de fluido líquido.
* CAE com a presença de ar.
* CT com presença de ar.
* CT com presença de fluido líquido.

Baseado no modelo de Yeoh [[#bib0185|[16]]] , os ligamentos foram considerados como tendo um comportamento não linear hiperelástico. As constantes ''c''<sub>''1''</sub> , ''c''<sub>''2''</sub>  e ''c''<sub>''3''</sub>  foram obtidas de referências contidas na literatura [[#bib0190|[17]]] .      

A simulação das articulações entre os ossículos, martelo/bigorna e bigorna/estribo foram efetuadas por intermédio de formulações matemáticas representativas de contacto [[#bib0185|[16]]]  and [[#bib0195|[18]]] , com um coeficiente de atrito igual a 0,9 [[#bib0200|[19]]] .      

A membrana timpânica foi fixa simulando o sulco timpânico. Anatomicamente, a ''pars flaccida''  está livre. A ''pars tensa''  foi fixa em toda a sua periferia, sendo que no quadrante póstero‐superior a membrana timpânica foi presa pela camada interna, utilizando‐se um nó para a prender. Nos outros quadrantes, a membrana timpânica foi presa pela camada interna e externa, utilizando‐se 2 nós.      

Os ossículos estão ainda suspensos por músculos [[#bib0205|[20]]]  e ligamentos [[#bib0165|[12]]]  and [[#bib0105|[21]]] ; martelo pelos ligamentos superior, lateral e anterior e pelo músculo tensor do tímpano; a bigorna pelos ligamentos superior e posterior e o estribo pelo músculo estapediano. Na periferia à volta da platina do estribo o ligamento anular foi modelado com 78 elementos de barra, do tipo T3D2, em que cada um dos nós da extremidade destes elementos foi fixo. Todos os elementos de fluido, que simulam a passagem para o CAE e para a CT, também foram fixos.      

Com o objetivo de perceber o comportamento do ouvido médio, ao longo de uma gama frequencial entre 100 Hz e10 kHz, efetuaram‐se simulações da aplicação de um nível de pressão sonora uniforme de 105 dB SPL (decibel Sound Pressure Level).

O dB SPL é uma medida logarítmica da pressão sonora (''p'' ) para um som em relação a um valor de referência (''p''<sub>0</sub>  = 20 μPa correspondente ao limiar de audibilidade), definida por:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SPL=20\times log_{10}(\frac{p}{p_0})</math>
|}
| style="text-align: right;" | ( 1)
|}

A carga aplicada na membrana timpânica foi simulada para um nível de pressão sonora uniforme de 105 dB SPL, a que corresponde uma pressão sonora de 3,56 Pa (pressão uniformemente distribuída por toda a face exterior da membrana timpânica). Esta pressão é aplicada no tempo em função sinusoidal de modo a induzir a excitação da membrana.

==3. Resultados==

A fim de validar o modelo criado para a cadeia tímpano‐ossicular, foi feito um estudo dinâmico para uma gama frequencial compreendida entre 100 Hz e 10 kHz. Para se poder comparar com os resultados obtidos do estudo experimental de Kurokawa et al. [[#bib0115|[2]]] , foi aplicada na membrana timpânica a mesma pressão sonora (3,56 Pa). As [[#fig0015|Figura 3]]  and [[#fig0020|Figura 4]] , que mostram respetivamente os deslocamentos do umbo e da platina do estribo, permitem verificar que os resultados obtidos no presente modelo estão em conformidade com os resultados obtidos por Kurokawa et al. [[#bib0115|[2]]] .

<span id='fig0015'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_861468639-1-s2.0-S0213131515000334-gr3.jpg|center|335px|Deslocamentos do umbo para 105dB SPL.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 3.

Deslocamentos do umbo para 105 dB SPL.

</span>
|}

<span id='fig0020'></span>

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; max-width: 100%;" 
|-
|


[[Image:draft_Content_861468639-1-s2.0-S0213131515000334-gr4.jpg|center|312px|Deslocamentos de um ponto central da platina do estribo para 105dB SPL.]]


|-
| <span style="text-align: center; font-size: 75%;">

Figura 4.

Deslocamentos de um ponto central da platina do estribo para 105 dB SPL.

</span>
|}

A partir das simulações efetuadas foram obtidos os deslocamentos do umbo e de um ponto central da platina do estribo. Compararam‐se os deslocamentos dessas simulações entre si, com os resultados do modelo apenas da cadeia tímpano‐ossicular (considerado modelo normal) e com os resultados obtidos por outro autor presente na literatura [[#bib0115|[2]]] .      

Verifica‐se na [[#fig0015|figura 3]]  que os deslocamentos do umbo obtidos na simulação que possui pressão de fluido líquido na CT são de amplitude inferior aos outros resultados. Em casos de otite, o interior da CT fica preenchida por líquido, situação esta que pode ser comparada com esta simulação.      

Para os outros resultados (pressão de fluido líquido, pressão de ar no CAE e pressão de ar na CT) as diferenças não são tão significativas. No entanto, nas frequências mais graves notam‐se maiores diferenças, podendo observar‐se que os deslocamentos são do menor para o maior, respetivamente, simulação com presença de fluido líquido no CAE, simulação com ar na CT e simulação com ar no CAE. Pode‐se observar nos resultados da [[#fig0020|figura 4]]  que os deslocamentos de um ponto central da platina do estribo, obtidos nestas simulações, apresentam o mesmo comportamento dos deslocamentos do umbo apresentados na [[#fig0015|figura 3]] .      

É importante salientar que sempre que há fluido líquido envolvido os deslocamentos são de amplitude inferior e quando envolve a CT são ainda menores.

Os deslocamentos obtidos com o modelo apenas com a cadeia tímpano‐ossicular são de amplitude superior aos deslocamentos obtidos com o modelo que possuiu CT e passagem para o CAE.

Quando comparamos a simulação com o modelo com a CT e o CAE, utilizando ar, vê‐se uma maior movimentação da membrana timpânica do que quando utilizado fluido líquido.

==4. Conclusões==

A presente investigação teve como objetivo o estudo de diferentes fluidos aplicados à membrana timpânica pela aplicação do método dos elementos finitos. Para o efeito, utilizou‐se um modelo computacional para a simulação biomecânica da cadeia tímpano‐ossicular do ouvido humano.

Aplicando à membrana timpânica um nível de pressão sonora de 105 dB SPL foram determinados os respetivos deslocamentos do umbo e da parte central da platina do estribo e os resultados foram comparados entre si e com os resultados de outro autor presente na literatura.

Para a simulação do CAE e da CT, consideraram‐se elementos de fluido nesta região para as diferentes partes constituintes. Simulou‐se a existência de ar e fluido líquido tanto no CAE como na CT. Foram feitas análises das diferentes pressões através da utilização de diferentes fluidos na CT. Relativamente aos deslocamentos, tanto do umbo como da platina do estribo, comparando a existência de ar e fluido líquido no CAE, os resultados obtidos com ar aproximam‐se mais da curva normal. Comparando a existência de ar e fluido líquido na CT, as diferenças são muito mais significativas, apresentando valores numa ordem de grandeza 10<sup>3</sup>  inferiores quando existe líquido dentro da CT. Pode concluir‐se que todos os resultados com a utilização de fluido são de amplitude inferior aos resultados com a utilização de pressão de ar.      

Uma visão mais aprofundada do funcionamento do ouvido humano, em busca de uma melhor intervenção futura, pode ser adquirida a partir de um melhor conhecimento da simulação do comportamento biomecânico do ouvido. Para tal, como trabalhos futuros pretende‐se a melhoria do modelo computacional, incluindo o ouvido interno, possibilitando assim uma avaliação completa do comportamento de todo o ouvido.

==Agradecimentos==

Agradece‐se o financiamento concedido ao projeto PEst‐OE/EME/LA0022/2013, da Fundação da Ciência e Tecnologia, Portugal e ao projeto «Biomechanics: contributions to the healthcare», referência NORTE‐07‐0124‐FEDER‐000035 cofinanciado pelo Programa Operacional Regional do Norte (ON.2–O Novo Norte), através do Fundo Europeu de Desenvolvimento Regional (FEDER).

==Referências==

<ol style='list-style-type: none;margin-left: 0px;'><li><span id='bib0110'></span>
[[#bib0110|[1]]] W.R.J. Funnell, C.A. Laszlo; Modeling of the cat eardrum as a thin shell using the finite element method; J Acoust Soc Am, 63 (1978), pp. 1461–1467</li>
<li><span id='bib0115'></span>
[[#bib0115|[2]]] H. Kurokawa, R. Goode; Sound pressure gain produced by the human middle ear; Otolaryngol Head Neck Surg, 113 (1995), pp. 349–355</li>
<li><span id='bib0120'></span>
[[#bib0120|[3]]] S. Nishiara, R. Goode; Measurement of tympanic membrane vibration in 99 human ears; K.B. Huttenbrink (Ed.), Middle ear mechanics in research and otosurgery, Dresden University of Technology, Dresden, Gernany (1996), pp. 91–93</li>
<li><span id='bib0125'></span>
[[#bib0125|[4]]] A. Huber, G. Ball, M. Asai, R. Goode; The vibration pattern of the tympanic membrane after placement of a total ossicular replacement prosthesis; Proceeding of the International Workshop on middle ear mechanics in research and otosurgery. Germany: Dresden (1997), pp. 219–222</li>
<li><span id='bib0130'></span>
[[#bib0130|[5]]] P. Prendergast, P. Ferris, H. Rice, A. Blayney; Vibro‐acoustic modelling of the outer and middle ear using the finite‐element method; Audiol Neurootol (1999), pp. 185–191</li>
<li><span id='bib0135'></span>
[[#bib0135|[6]]] R.Z. Gan, Q. Sun, K.J. Robert, K. Chang, K. Dormer; Three‐dimensional modeling of middle ear biomechanics and its applications; Otol Neurotol, 23 (2002), pp. 271–280</li>
<li><span id='bib0140'></span>
[[#bib0140|[7]]] Q. Sun, R. Gan, K. Chang, K. Dormer; Computer‐integrated finite element modeling of human middle ear; Biomech Model Mechanobiol, 1 (2002), pp. 109–122</li>
<li><span id='bib0145'></span>
[[#bib0145|[8]]] R.Z. Gan, B. Feng, Q. Sun; Three‐dimensional finite element modeling of human ear for sound transmission; Ann Biomed Eng, 32 (2004), pp. 847–859</li>
<li><span id='bib0150'></span>
[[#bib0150|[9]]] R.Z. Gan, Q. Sun, B. Feng, M.W. Wood; Acoustic‐structural coupled finite element analysis for sound transmission in human ear: Pressure distributions; Med Eng Phys, 28 (2006), pp. 395–404</li>
<li><span id='bib0155'></span>
[[#bib0155|[10]]] C.F. Lee, P.R. Chen, W.J. Lee, J.H. Chen, T.C. Liu; Computer aided three‐dimensional reconstruction and modeling of middle ear biomechanics by high‐resolution computed tomography and finite element analysis; Biomedical Engineering: Applications, Basis and Communications (5) (2006), pp. 214–221</li>
<li><span id='bib0160'></span>
[[#bib0160|[11]]] F. Gentil, R.M.N. Jorge, A.J.M. Ferreira, M.P.L. Parente, P.A.L.S. Martins, E. Almeida; Biomechanical simulation of middle ear using hyperelastic models; J. of Biomechanics, 39 (2006), pp. 388–389</li>
<li><span id='bib0165'></span>
[[#bib0165|[12]]] F. Gentil, M. Parente, P. Martins, C. Garbe, R.N. Jorge, J. Tavares,  ''et al.''; The influence of the mechanical behaviour of the middle ear ligaments: A finite element analysis; Part H: Journal of Engineering in Medicine, 225 (2010), pp. 68–76</li>
<li><span id='bib0170'></span>
[[#bib0170|[13]]] C. Garbe, F. Gentil, M. Parente, P. Martins, R.N. Jorge; Aplicação do método do elementos finitos no estudo da membrana timpânica; Audiologia em Revista, 2 (2009), pp. 99–106</li>
<li><span id='bib0175'></span>
[[#bib0175|[14]]] J. Paço, C. Branco, H. Estibeiro, D. Oliveira; The postero/superior quadrant of the tympanic membrane; Otolaryngol Head Neck Surg, 140 (2009), pp. 884–888</li>
<li><span id='bib0180'></span>
[[#bib0180|[15]]] D. Hibbit, B. Karlsson, P. Sorenson; ABAQUS Analysis Users Manual, version 6.5; Hibbit, Karlsson & Sorenson Inc., USA (2004)</li>
<li><span id='bib0185'></span>
[[#bib0185|[16]]] O.H. Yeoh; Characterization of elastic properties of car‐bon‐black‐filled rubber vulcanizates; Rubber Chemistry and Technology, 63 (1990), pp. 792–805</li>
<li><span id='bib0190'></span>
[[#bib0190|[17]]] P.A.L.S. Martins, R.M.N. Jorge, A.J.M. Ferreira; A comparative study of several material models for prediction of hyperelastic properties. Application to silicone‐rubber and soft tissues; Strain, 42 (2006), pp. 135–147</li>
<li><span id='bib0195'></span>
[[#bib0195|[18]]] P. Wriggers; Computational contact mechanics; John Wiley & sons Ltd (eds) (2002)</li>
<li><span id='bib0200'></span>
[[#bib0200|[19]]] F. Gentil, R.M.N. Jorge, A.J.M. Ferreira, M.P.L. Parente, M. Moreira, E. Almeida; Estudo do efeito do atrito no contacto entre os ossículos do ouvido médio; Rev. int. métodos numér. cálc. diseño ing., 23 (2) (2007), pp. 177–187</li>
<li><span id='bib0205'></span>
[[#bib0205|[20]]] F. Gentil, M. Parente, P. Martins, C. Garbe, R. Natal, A. Ferreira,  ''et al.''; The influence of muscles activation on the dynamical behaviour of the tympano‐ossicular system of the middle ear; Comput. Methods Biomech. Biomed. Engin., 16 (4) (2013), pp. 392–402</li>
<li><span id='bib0105'></span>
[[#bib0105|[21]]] F. Gentil, C. Garbe, M. Parente, P. Martins, A. Ferreira, R. Natal,  ''et al.''; Analysis of eardrum pathologies using the finite element method; J. Mech. Med. Biol., 14 (3) (2014) World Scientific Publishing Company DOI: 10.1142/S0219519414500341</li>
</ol>