Physiological Models
The physiological models are simple equivalent circuit models for parts of the human body which typically are of importance in acoustical applications, namely the skin, the full human ear, the eardrum and inner ear, and the acoustic impedance from the ear’s pinna. These models are good approximations to the active, acoustical properties of these body parts with particular relevance for development of hearing aids, mobile devices, and head phones.
Human skin model
This model has been developed as a lowest-order, reasonable approximation to experimental data in Ref. 27 (see their Figure 9). It consists of a simple serial RCL circuit applied on a transducer area of At, and has the impedance
with At = 1.5 cm2, Rs = 9.0 N s m-1, Ls = 0.53·10-3 N s2·m-1 and Cs = 5.3·10-6 m N-1.
Models related to the human ear
Four models related to the human ear have been included. These detailed and experimentally-verified circuit models describe parts of the human ear, as well as the entire human ear, see Ref. 28 to 30 for further details.
Outward Human Ear Radiation
This impedance describes the acoustic radiation losses from the pinna (also known as the auricle, this is the visible part of the ear which is exterior to the head), see Figure 2-11. For cases where you model the ear canal explicitly using pressure acoustics, this boundary condition describes the acoustic losses from the outward acoustic radiation from the ear canal and into the surrounding air.
(2-29)
where Rpar = 7.0·106 N·s·m-5, Lpar = 100 N·s2·m-5, Cpar = 1.7·10-12 m5 N-1, Q = 6, R1 = R2 = Rpar, R3 = 2Rpar, ω1 = 6000·2π Hz, ω1 = 9000·2π Hz, and ω1 = 13000·2π Hz.
Figure 2-11: Illustration of outward human ear radiation.
Human Ear Drum Impedance
This model describes the impedance of the human ear drum and the entire inner ear, that is, the acoustic impedance experienced in the ear canal when looking into the ear drum, see Figure 2-12. The model equations are given in Equation 2-30 and Equation 2-31, and the parameter values in Table 2-8.
Figure 2-12: Illustration of human eardrum impedance.
(2-30)
(2-31)
The parameter values are given in Table 2-8. Note that the value of is not reported in the papers Ref. 28 to 30, but has instead been determined during model implementation. The value ensures a continuity of the phase response.
Rtcav
2·106 N s m-5
Ctcav
Rada
1.7·106 N s m-5
Lada
880 N s2 m-5
Cant
Qmac
Cmac
ωmac
2π3500 rad s-1
Rac
4·107 N s m-5
Cac
5·10-12 N-1 m5
2.4·103 kg m-4
2π·1900 rad s-1
ωYph
2π·8000 rad s-1
sYph
A0
ωA
2π·2200 rad s-1
QA
sAph
-1.2
ωAph
2π·1500 rad s-1
Rmi
Cmi
Coss
3·10-3 m N-1
Loss
7·10-3 g
Ccpl
0.5·10-3 m N-1
Rcpl
Rfree
Lfree
12·10-3 g
Rst
18·10-3 N s m-1
Lst
3·10-3 g
Cst
1.2·10-3 m N-1
Rc
70·10-3 N s m-1/ AF2
Lc
10·10-3 g / AF2
Cc
11·10-3 m N-1·AF2
AF
Human Ear Without Pinna
This model accounts for the acoustic losses associated with the ear canal and the entire human ear, see Figure 2-13. It does not include the radiation losses associated with the pinna, the visible part of the ear which is external to the head.
(2-32)
Here, Zeardrum is the eardrum impedance defined in Equation 2-31 and Tij are the components of the transfer matrix T that describes the ear canal as a two-port. The ear canal is treated as Ntot small segments each with length Δk and radius rk so its full two port T is given by
(2-33).
where
(2-34)
In these expressions, γ is the ratio of specific heats, Γk is the propagation constant (“wave number”) of the kth ear canal segment which has the segment-specific attenuation αk,
is the segment-specific Womersley number, and
is the Prandtl number expressed in terms of the specific heat Cp, dynamic viscosity μ, and thermal conductivity k. Notice that the papers presenting the model (Ref. 28 to 30) do not exactly specify which expression for the attenuation constants αk is being used, but only refer to Ref. 31. The expression above for αk is the most general expression taken from this paper. The values for Δk and rk are listed in Table 2-9.
Figure 2-13: Illustration of the impedance of the human ear without pinna.
Table 2-9: Radii rk and lengths Δk of ear canal segments taken from Ref. 30.
rk (mm)
Δk (mm)
The pressure at the eardrum peardrum is calculated whenever this impedance boundary condition is applied. This pressure is available in postprocessing, and is calculated from the expression
,
where pt is the pressure on the boundary, Tij are the coefficients of the ear canal transfer matrix T and Zear w/o pinna is the ear impedance; T and Zear w/o pinna are defined in Equation 2-32 and Equation 2-33 above.
Human Ear, Full
This model accounts for all acoustic losses associated with the entire human ear, both the internal parts as well as the pinna, the visible, external part of the ear on the head, see Figure 2-14. The model does not include any information about the directivity of the ear (the head related transfer functions, HRTFs) which depends on the ear geometry at higher frequencies. It is valid in the low frequency limit and for normal incidence on the ear. It is given by
where Zrad and Zear w/o pinna are given by Equation 2-29 and Equation 2-32 above.
Figure 2-14: Illustration of the impedance of the full human ear including radiation losses due to the pinna.
The pressure at the eardrum peardrum is calculated whenever this impedance boundary condition model is applied. This pressure is calculated from Zear full using the expression
.