Effects of poroelastic coefficients on normal vibration modes in vocal-fold tissues

The vocal-fold tissue is treated as a transversally isotropic fluid-saturated porous material. Effects of poroelastic coefficients on eigenfrequencies and eigenmodes of the vocal-fold vibration are investigated using the Ritz method. The study demonstrates that the often-used elastic model is only a particular case of the poroelastic model with an infinite fluid-solid mass coupling parameter. The elastic model may be considered appropriate for the vocal-fold tissue when the absolute value of the fluid-solid mass coupling parameter is larger than 10(5) kg/m(3). Otherwise, the poroelastic model may be more accurate. The degree of compressibility of the vocal tissue can also been described by the poroelastic coefficients. Finally, it is revealed that the liquid and solid components in a poroelastic model could have different modal shapes when the coupling between them is weak. The mode decoupling could cause desynchronization and irregular vibration of the folds.     

A mechanical model of vocal-fold collision with high spatial and temporal resolution

The tissue mechanics governing vocal-fold closure and collision during phonation are modeled in order to evaluate the role of elastic forces in glottal closure and in the development of stresses that may be a risk factor for pathology development. The model is a nonlinear dynamic contact problem that incorporates a three-dimensional, linear elastic, finite-element representation of a single vocal fold, a rigid midline surface, and quasistatic air pressure boundary conditions. Qualitative behavior of the model agrees with observations of glottal closure during normal voice production. The predicted relationship between subglottal pressure and peak collision force agrees with published experimental measurements. Accurate predictions of tissue dynamics during collision suggest that elastic forces play an important role during glottal closure and are an important determinant of aerodynamic variables that are associated with voice quality. Model predictions of contact force between the vocal folds are directly proportional to compressive stress (r2 = 0.79), vertical shear stress (r2 = 0.69), and Von Mises stress (r2 = 0.83) in the tissue. These results guide the interpretation of experimental measurements by relating them to a quantity that is important in tissue damage.     

A nonlinear finite-element model of the newborn ear canal

A three-dimensional nonlinear finite-element model of a 22-day-old newborn ear canal is presented. The geometry is based on a clinical x-ray CT scan. A nonlinear hyperelastic constitutive law is applied to model large deformations. The Young's modulus of the soft tissue is found to have a significant effect on the ear-canal volume change, which ranges from approximately 27% to 75% over the static-pressure range of +/-3kPa. The effects of Poisson's ratio and of the ratio C10: C01 in the hyperelastic model are found to be small. The volume changes do not reach a plateau at high pressures, which implies that the newborn ear-canal wall would not be rigid in tympanometric measurements. The displacements and volume changes calculated from the model are compared with available experimental data.     

A nonlinear finite-element model of the newborn middle ear

A three-dimensional static nonlinear finite-element model of a 22-day-old newborn middle ear is presented. The model includes the tympanic membrane (TM), malleus, incus, and two ligaments. The effects of the middle-ear cavity are taken into account indirectly. The geometry is based on a computed-tomography scan and on the published literature, supplemented by histology. A nonlinear hyperelastic constitutive law is applied to model large deformations. The middle-ear cavity and the Young's modulus of the TM have significant effects on TM volume displacements. The TM volume displacement and its nonlinearity and asymmetry increase as the middle-ear cavity volume increases. The effects of the Young's moduli of the ligaments and ossicles are found to be small. The simulated TM volume changes do not reach a plateau when the pressure is varied to either -3 kPa or +3 kPa, which is consistent with the nonflat tails often found in tympanograms in newborns. The simulated TM volume displacements, by themselves and also together with previous ear-canal model results, are compared with equivalent-volume differences derived from tympanometric measurements in newborns. The results suggest that the canal-wall volume displacement makes a major contribution to the total canal volume change, and may be larger than the TM volume displacement.     

Direct-numerical simulation of the glottal jet and vocal-fold dynamics in a three-dimensional laryngeal model

An immersed-boundary method based flow solver coupled with a finite-element solid dynamics solver is employed in order to conduct direct-numerical simulations of phonatory dynamics in a three-dimensional model of the human larynx. The computed features of the glottal flow including mean and peak flow rates, and the open and skewness quotients are found to be within the normal physiological range. The flow-induced vibration pattern shows the classical "convergent-divergent" glottal shape, and the vibration amplitude is also found to be typical for human phonation. The vocal fold motion is analyzed through the method of empirical eigenfunctions and this analysis indicates a 1:1 modal entrainment between the "adduction-abduction" mode and the "mucosal wave" mode. The glottal jet is found to exhibit noticeable cycle-to-cycle asymmetric deflections and the mechanism underlying this phenomenon is examined.     

A model of audio-frequency vibration of buildings

A model is proposed wherein audio-frequency vibrations in buildings are transmitted by column motion in the fixed-fixed mode. Transmission and attenuation depend on the relative tuning of neighboring storey columns.     

A model for vibration transmission of light-heavy structures

The vibration transmission of light-heavy structures is investigated in this paper. The light-heavy structure consists of a thin beam and a mass block. Based on numerical simulations with the finite element method and experiments, the block's effect on the thin beam is defined. A theoretical model for this beam-block structure is successfully developed, which is validated and agrees very well with the numerical and experimental models. Two kinds of transfer functions of velocities between any two points on the beam-block structure are studied experimentally and theoretically. The theoretical transfer functions agree well with the experimental results. There are peaks and valleys in the transfer functions, where the peaks occur at the anti-resonant frequencies of the second point and the valleys at the anti-resonant frequencies of the first point. Away from these peaks and valleys the magnitude of the transfer functions are about 0 dB for two points on the beam, and about 20 dB in our experiments for a point on the beam and another point on the block (close to the theoretical prediction of 18 dB determined by the mass ratio of the beam and the block). With these transfer functions, new techniques might be developed for indirect measurement of the vibration of the thin beam by measuring the vibration of the block.     

Irregular vocal-fold vibration--high-speed observation and modeling

Direct observations of nonstationary asymmetric vocal-fold oscillations are reported. Complex time series of the left and the right vocal-fold vibrations are extracted from digital high-speed image sequences separately. The dynamics of the corresponding high-speed glottograms reveals transitions between low-dimensional attractors such as subharmonic and quasiperiodic oscillations. The spectral components of either oscillation are given by positive linear combinations of two fundamental frequencies. Their ratio is determined from the high-speed sequences and is used as a parameter of laryngeal asymmetry in model calculations. The parameters of a simplified asymmetric two-mass model of the larynx are preset by using experimental data. Its bifurcation structure is explored in order to fit simulations to the observed time series. Appropriate parameter settings allow the reproduction of time series and differentiated amplitude contours with quantitative agreement. In particular, several phase-locked episodes ranging from 4:5 to 2:3 rhythms are generated realistically with the model.     

具有直线结构的三原子分子振动模型

给出了直线结构的三原子分子振动的一个模型,该模型通过角簧实现弯曲振动模式.建立了该模型的振动微分方程组.求解方程组得到了3个本征振动频率和对应的振动模式.就二氧化碳、羰基硫和氰化氢3个三原子分子,使用文献数据验证了所给出的模型.结果表明按该模型的3个振动频率与文献结果吻合良好.     

A microstructure model for finite-element simulation of 3D rectangular braided composite under ballistic penetration

The non-delamination feature of 3D braided composites under transverse impact leads to their potential application in the field of ballistic impact protection. One of the effective ways to investigate the ballistic impact damage of the 3D braided composite is to simulate the penetration process by numerical method, such as finite element method. However the numerical simulations of ballistic impact damage are seldom conducted based on the microstructure level. This paper presents a microstructure model for simulating ballistic impact damage of 4-step 3D braided rectangular composite penetrated by a rigid steel projectile. The microstructure model is based on the same yarn spatial configuration with that of the braided composite and also on the assumptions of the braided yarns appear straight inside the braided preform, bending and then change to other directions only at the surface. The ballistic perforation of the braided composite specimen by a cylindrical-conically steel projectile has been simulated with finite element method. The comparisons between FEA and experimental results show the validity of the microstructure model, especially for the penetration resistance and impact damage of the composite. Compared with the other continuum models of the braided composite, the microstructure model can simulate impact damage more precisely. The velocity history and acceleration history of projectile, and impact damage development of the composite in FEM simulation indicate the different damage and energy absorption mechanisms of the braided composite compared with those of laminated composite.     

A general model of the axle vibration in piezoelectric motors

In the present work a general model of the vibrational behavior of the axle of a piezoelectric motor is proposed. In this motor, a cylinder-shaped permanent magnet, which act as a rotor, is pressed in contact with an end of a steel axle by means of the magnetic forces. The other end of the axle is fitted at the center of a rotating traveling wave generator. A piezoelectric membrane, vibrating in a flexural anti-symmetrical mode, or a thick disk, vibrating in a radial anti-symmetrical mode, can be exploited as traveling wave generators. In the first case a bending moment, in the second case a transverse force is applied to the axle. In both cases, if the driving frequency coincides with a resonance frequency of the axle, the axle acts as a resonant displacement amplifier; a continuous slipping takes place between the axle and the rotor, and a torque is transmitted to the rotor. The proposed model is able to describe the axle vibrational behavior when it is excited by a bending moment, by a transverse force, and also when these two excitations are simultaneously applied. The axle is modeled as a four-port system and all its transfer functions, as well as the transversal displacement along the axle at each frequency can be easily computed. Computed results have been compared with experimental measurements carried out on two motor prototypes that exploit as traveling wave generators a membrane and a disk, respectively. A good agreement was obtained by properly taking into account the loading effect of the generator on the axle.     

Simulation of vocal fold impact pressures with a self-oscillating finite-element model

Vocal fold impact pressures were studied using a self-oscillating finite-element model capable of simulating vocal fold vibration and airflow. The calculated airflow pressure is applied on the vocal fold as the driving force. The airflow region is then adjusted according to the calculated vocal fold displacement. The interaction between airflow and the vocal folds produces a self-oscillating solution. Lung pressures between 0.2 and 2.5 kPa were used to drive this self-oscillating model. The spatial distribution of the impact pressure was studied. Studies revealed that the tissue collision during phonation produces a very large impact pressure which correlates with the lung pressure and glottal width. Larger lung pressure and a narrower glottal width increase the impact pressure. The impact pressure was found to be roughly the square root of lung pressure. In the inferior-superior direction, the maximum impact pressure is related to the narrowest glottis. In the anterior-posteriorfirection, the greatest impact pressure appears at the midpoint of the vocal fold. The match between our numerical simulations and clinical observations suggests that this self-oscillating finite-element model might be valuable for predicting mechanical trauma of the vocal folds.     

Application of a finite-element model to low-frequency sound insulation in dwellings

The sound transmission between adjacent rooms has been modeled using a finite-element method. Predicted sound-level difference gave good agreement with experimental data using a full-scale and a quarter-scale model. Results show that the sound insulation characteristics of a party wall at low frequencies strongly depend on the modal characteristics of the sound field of both rooms and of the partition. The effect of three edge conditions of the separating wall on the sound-level difference at low frequencies was examined: simply supported, clamped, and a combination of clamped and simply supported. It is demonstrated that a clamped partition provides greater sound-level difference at low frequencies than a simply supported. It also is confirmed that the sound-pressure level difference is lower in equal room than in unequal room configurations.     

A matrix model of the axle vibration of a piezoelectric motor

In this work, a matrix model of the axle vibration of a piezoelectric motor is proposed. The stator of this motor is composed of a thin piezoelectric membrane and a steel axle fitted at the center of the membrane. The rotor consists of a cylinder-shaped permanent magnet, pressed in contact with the other end of the axle by means of the magnetic forces. A travelling wave is excited in the membrane by using four electrodes and four, properly delayed, driving signals. The rotating flexural displacement of the membrane produces a wide precessional motion of the axle. In this way, a continuous slipping takes place between the axle and the rotor, and therefore, a torque is transmitted to the rotor. In this paper, the precessional motion of the axle is modeled as the composition of two transverse vibrations belonging to two perpendicular planes passing through the axle. The axle, vibrating in its transverse mode, is modeled as a two-port system: the input is the bending moment supplied by the membrane, and the output is the transverse force at the terminal end of the axle. With this model, we have computed the trasmission transfer function as a function of frequency, and the transversal displacement along the axle at its resonance frequency. The computed results are in reasonable agreement with experimental interferometric measurements carried out on a prototype.     

A finite-element model of the aperture method for determining the effective radiating area of physiotherapy treatment heads

This paper describes a theoretical study of the way in which a circular aperture placed in front of a plane-piston modifies the ultrasonic field it generates. Specifically, the radiated acoustic power transmitted by the aperture and the radiation force experienced by an absorbing target placed in the transmitted beam, are evaluated at a distance from the exit-side of the aperture. The calculations used a finite element (FE) method, in conjunction with a surface Helmholtz integral formulation to solve the fluid/structure interaction problem. The PAFEC (Program for Automatic Finite Element Computation) vibroacoustics software was used for the FE calculations and the implementation of the surface Helmholtz integral formulation method. Acoustic pressures and particle velocities were computed at required points, whilst accounting for the presence of the aperture in the medium, together with its dynamic properties when subjected to an incident sound field. This enabled the calculation of the radiation force on the absorber and its variation with the applied aperture diameter was investigated. As part of the validation process for the new FE aperture model, the ratio of radiation force to radiated acoustic power obtained using the FE method in the unapertured case, through the use of the Rayleigh integral, yielded good agreement with results obtained through an analytical solution. The study has been carried out to provide a better understanding of the factors affecting the measurement uncertainty for the aperture method of determining the effective radiating area (A(ER)) of physiotherapy ultrasound treatment heads.     

Compliance profiles derived from a three-dimensional finite-element model of the basilar membrane

A finite-element analysis is used to explore the impact of elastic material properties, boundary conditions, and geometry, including coiling, on the spatial characteristics of the compliance of the unloaded basilar membrane (BM). It is assumed that the arcuate zone is isotropic and the pectinate zone orthotropic, and that the radial component of the effective Young's modulus in the pectinate zone decreases exponentially with distance from base to apex. The results concur with tonotopic characteristics of compliance and neural data. Moreover, whereas the maximum compliance in a radial profile is located close to the boundary between the two zones in the basal region, it shifts to the midpoint of the pectinate zone for the apical BM; the width of the profile also expands. This shift begins near the 1 kHz characteristic place for guinea pig and the 2.4 kHz place for gerbil. Shift and expansion are not observed for linear rather than exponential decrease of the radial component of Young's modulus. This spatial change of the compliance profile leads to the prediction that mechanical excitation in the apical region of the organ of Corti is different to that in the basal region.     

On the damped frequency response of a finite-element model of the cat eardrum

This article presents frequency responses calculated using a three-dimensional finite-element model of the cat eardrum that includes damping. The damping is represented by both mass-proportional and stiffness-proportional terms. With light damping, the frequency responses of points on the eardrum away from the manubrium display numerous narrow minima and maxima, the frequencies and amplitudes of which are different for different positions on the eardrum. The frequency response on the manubrium is smoother than that on the eardrum away from the manubrium. Increasing the degree of damping smooths the frequency responses both on the manubrium and on the eardrum away from the manubrium. The overall displacement magnitudes are not significantly reduced even when the damping is heavy enough to smooth out all but the largest variations. Experimentally observed frequency responses of the cat eardrum are presented for comparison with the model results.     

A non-linear fluid force model for vortex-induced vibration of an elastic cylinder

Even under the assumption of a sinusoidal lift and drag force at a single frequency for a stationary cylinder in a cross flow, higher harmonics that represent non-linearity in the fluid-structure interaction process are present. This fact is considered in the formulation of a non-linear fluid force model for a freely vibrating cylinder in a cross flow. The force model is developed based on an iterative process and the modal analysis approach. The fluid force components in the model can be evaluated from measured vibration data with the help of the auto-regressive moving averaging (ARMA) technique. An example is used to illustrate that non-linear (higher order) force components are present at resonance, even for a case with relatively weak fluid-structure interaction. Further analysis reveals that the fluid force components are dependent on structural damping and mass ratio. The non-linear fluid force model is further modified by taking these considerations into account and is used to predict the dynamic characteristics of a freely vibrating cylinder over a range of Reynolds numbers, mass and structural damping ratios. On comparison with measurements obtained from four different experiments and predictions made by previous single-degree-of-freedom model, good agreement is found over a wide range of these parameters.     

An automated parameter selection procedure for finite-element model updating and its applications

Finite-element model updating is an inverse problem to identify and correct uncertain modeling parameters, which leads to better predictions of the dynamic behavior of a target structure. Unlike other inverse problems, the restrictions on selecting parameters are very high since the updated model should maintain its physical meaning. That is, only the regions with modeling errors should be parameterized and the variations of the parameters should be kept small while the updated results give acceptable correlations with experimental data. To avoid an ill-conditioned numerical problem, the number of parameters should be kept as small as possible. Thus it is very difficult to select an adequate set of updating parameters which meet all these requirements. In this paper, the importance of updating parameter selection is illustrated through a case study, and an automated procedure to guide the parameter selection is suggested based on simple observations. The effectiveness of the suggested procedure is tested with two example problems, one is a simulated case study and the other is a real engineering structure.