A ferroelectret nanogenerator (FENG) was recently developed as a flexible energy harvesting device with bi-directional capability between electrical and mechanical energy domains, and its use as a loudspeaker/microphone was demonstrated. Dependencies of Sound Pressure Levels (SPLs) generated by FENG due to an AC voltage stimulus, surface area, geometric shape, and addition of layers are presented here. Also, the relation between the sound output to the electrical input is studied and shown to be linear, which demonstrates that these flexible loudspeakers have low distortion within the human audible range of 20 Hz to 20 kHz. A study for ultrasonic frequencies up to 40 kHz is also presented. A theoretical model relating the electrical and acoustical domain of the FENG is developed based on the experimental observations made and using Boundary Element Methods (BEM) to accurately mimic the testing environment for simulation purposes. The comparison between this model and the actual behavior is presented under several cases and observed to be closely correlated.

While simulations of the measured biodynamic responses of the whole human body or body segments to vibration are conventionally interpreted as summaries of biodynamic measurements, and the resulting models are considered quantitative, this study looked at these simulations from a different angle: model calibration. The specific aims of this study are to review and clarify the theoretical basis for model calibration, to help formulate the criteria for calibration validation, and to help appropriately select and apply calibration methods. In addition to established vibration theory, a novel theorem of mechanical vibration is also used to enhance the understanding of the mathematical and physical principles of the calibration. Based on this enhanced understanding, a set of criteria was proposed and used to systematically examine the calibration methods. Besides theoretical analyses, a numerical testing method is also used in the examination. This study identified the basic requirements for each calibration method to obtain a unique calibration solution. This study also confirmed that the solution becomes more robust if more than sufficient calibration references are provided. Practically, however, as more references are used, more inconsistencies can arise among the measured data for representing the biodynamic properties. To help account for the relative reliabilities of the references, a baseline weighting scheme is proposed. The analyses suggest that the best choice of calibration method depends on the modeling purpose, the model structure, and the availability and reliability of representative reference data.

This paper provides new formulations to derive the impulse response matrix, which is then used in the problem of load identification with application to wind induced vibration. The applied loads are inversely identified based on the measured structural responses by solving the associated discrete ill-posed problem. To this end - based on an existing parametric structural model - the impulse response functions of acceleration, velocity and displacement have been computed. Time discretization of convolution integral has been implemented according to an existing and a newly proposed procedure, which differ in the numerical integration methods. The former was evaluated based on a constant rectangular approximation of the sampled data and impulse response function in a number of steps corresponding to the sampling rate, while the latter interpolates the sampled data in an arbitrary number of sub-steps and then integrates over the sub-steps and steps. The identification procedure was implemented for a simulation example as well as an experimental laboratory case. The ill-conditioning of the impulse response matrix made it necessary to use Tikhonov regularization to recover the applied force from noise polluted measured response. The optimal regularization parameter has been obtained by L-curve and GCV method. The results of simulation represent good agreement between identified and measured force. In the experiments the identification results based on the measured displacement as well as acceleration are provided. Further it is shown that the accuracy of experimentally identified load depends on the sensitivity of measurement instruments over the different frequency ranges.

Breath sounds are often used to aid in the diagnosis of pulmonary disease. Mechanical and numerical models could be used to enhance our understanding of relevant sound transmission phenomena. Sound transmission in an airway mimicking phantom was investigated using a mechanical model with a branching airway network embedded in a compliant viscoelastic medium. The Horsfield self-consistent model for the bronchial tree was adopted to topologically couple the individual airway segments into the branching airway network. The acoustics of the bifurcating airway segments were measured by microphones and calculated analytically. Airway phantom surface motion was measured using scanning laser Doppler vibrometry. Finite element simulations of sound transmission in the airway phantom were performed. Good agreement was achieved between experiments and simulations. The validated computational approach can provide insight into sound transmission simulations in real lungs.

The relationship between the vibration transmissibility and driving-point response functions (DPRFs) of the human body is important for understanding vibration exposures of the system and for developing valid models. This study identified their theoretical relationship and demonstrated that the sum of the DPRFs can be expressed as a linear combination of the transmissibility functions of the individual mass elements distributed throughout the system. The relationship is verified using several human vibration models. This study also clarified the requirements for reliably quantifying transmissibility values used as references for calibrating the system models. As an example application, this study used the developed theory to perform a preliminary analysis of the method for calibrating models using both vibration transmissibility and DPRFs. The results of the analysis show that the combined method can theoretically result in a unique and valid solution of the model parameters, at least for linear systems. However, the validation of the method itself does not guarantee the validation of the calibrated model, because the validation of the calibration also depends on the model structure and the reliability and appropriate representation of the reference functions. The basic theory developed in this study is also applicable to the vibration analyses of other structures.

The method of tailored Green's functions advocated by Doak ( A254 (1960) 129 - 145.) for the solution of aeroacoustic problems is used to analyse the contribution of the mucosal wave to self-sustained modulation of air flow through the glottis during the production of voiced speech. The amplitude and phase of the aerodynamic surface force that maintains vocal fold vibration are governed by flow separation from the region of minimum cross-sectional area of the glottis, which moves back and forth along its effective length accompanying the mucosal wave peak. The correct phasing is achieved by asymmetric motion of this peak during the opening and closing phases of the glottis. Limit cycle calculations using experimental data of Berry ( 110 (2001) 2539 - 2547.) obtained using an excised canine hemilarynx indicate that the mechanism is robust enough to sustain oscillations over a wide range of voicing conditions.

The objectives of this study were to develop models of the hand-arm system in the three orthogonal directions ( , and ) and to enhance the understanding of the hand vibration dynamics. A four-degrees-of-freedom (DOF) model and 5-DOF model were used in the simulation for each direction. The driving-point mechanical impedances distributed at the fingers and palm of the hand reported in a previous study were used to determine the parameters of the models. The 5-DOF models were generally superior to the 4-DOF models for the simulation. Hence, as examples of applications, the 5-DOF models were used to predict the transmissibility of a vibration-reducing glove and the vibration transmissibility on the major substructures of the hand-arm system. The model-predicted results were also compared with the experimental data reported in two other recent studies. Some reasonable agreements were observed in the comparisons, which provided some validation of the developed models. This study concluded that the 5-DOF models are acceptable for helping to design and analyze vibrating tools and anti-vibration devices. This study also confirmed that the 5-DOF model in the direction is acceptable for a coarse estimation of the biodynamic responses distributed throughout the major substructures of the hand-arm system. Some interesting phenomena observed in the experimental study of the biodynamic responses in the three directions were also explained in this study.

Model updating procedures are applied in order to improve the matching between experimental data and corresponding model output. The updated, i.e. improved, finite element (FE) model can be used for more reliable predictions of the structural performance in the target mechanical environment. The discrepancies between the output of the FE-model and the results of tests are due to the uncertainties that are involved in the modeling process. These uncertainties concern the structural parameters, measurement errors, the incompleteness of the test data and also the FE-model itself. The latter type of uncertainties is often referred to as model uncertainties and is caused by simplifications of the real structure that are made in order to reduce the complexity of reality. Several approaches have been proposed for taking model uncertainties into consideration, where the focus of this manuscript will be set on the updating procedure within the Bayesian statistical framework. A numerical example involving different degrees of nonlinearity will be used for demonstrating how this type of uncertainty is considered within the Bayesian updating procedure.

An analysis is made of the fluid-structure interactions involved in the production of voiced speech. It is usual to avoid time consuming numerical simulations of the aeroacoustics of the vocal tract and glottis by the introduction of Fant's 'reduced complexity' equation for the glottis volume velocity Q (G. Fant, Acoustic Theory of Speech Production, Mouton, The Hague 1960). A systematic derivation is given of Fant's equation based on the nominally exact equations of aerodynamic sound. This can be done with a degree of approximation that depends only on the accuracy with which the time-varying flow geometry and surface-acoustic boundary conditions can be specified, and replaces Fant's original 'lumped element' heuristic approach. The method determines all of the effective 'source terms' governing Q. It is illustrated by consideration of a simplified model of the vocal system involving a self-sustaining single-mass model of the vocal folds, that uses free streamline theory to account for surface friction and flow separation within the glottis. Identification is made of a new source term associated with the unsteady vocal fold drag produced by their oscillatory motion transverse to the mean flow.

The steady state response motion of a base excited cantilever beam with circular cross-section excited by a unidirectional displacement will fall along a straight line. However, achieving straight-line motion with a real cantilever beam of circular cross-section is difficult to accomplish. This is due to the fact that nonlinear effects, small deviations from circularity, asymmetric boundary conditions, and actuator cross coupling can induce whirling. The vast majority of previous work on cantilever beam whirling has focused on the effects of system nonlinearities. We show that whirling is a much broader problem in the design of resonant beam scanners in that the onset of whirling does not depend on large amplitude of motion. Rather, whirling is the norm in real systems due to small system asymmetries and actuator cross coupling. It is therefore necessary to control the growth of the whirling motion when a unidirectional beam motion is desired. We have developed a novel technique to identify the two eigen directions of the beam. Base excitation generated by virtual electrodes along these orthogonal eigen axes of the cantilever beam system generates tip vibration without whirl. This leads to accurate open loop control of the motion of the beam through the combined actuation of two pairs of orthogonally placed actuator electrodes.

Quadratic matrix polynomials are fundamental to vibration analysis. Because of the predetermined interconnectivity among the constituent elements and the mandatory nonnegativity of the physical parameters, most given vibration systems will impose some inherent structure on the coefficients of the corresponding quadratic matrix polynomials. In the inverse problem of reconstructing a vibration system from its observed or desirable dynamical behavior, respecting the intrinsic structure becomes important and challenging both theoretically and practically. The issue of whether a structured inverse eigenvalue problem is solvable is problem dependent and has to be addressed structure by structure. In an earlier work, physical systems that can be modeled under the paradigm of a serially linked mass-spring system have been considered via specifically formulated inequality systems. In this paper, the framework is generalized to arbitrary generally linked systems. In particular, given any configuration of interconnectivity in a mass-spring system, this paper presents a mechanism that systematically and automatically generates a corresponding inequality system. A numerical approach is proposed to determine whether the inverse problem is solvable and, if it is so, computes the coefficient matrices while providing an estimate of the residual error. The most important feature of this approach is that it is problem independent, that is, the approach is general and robust for any kind of physical configuration. The ideas discussed in this paper have been implemented into a software package by which some numerical experiments are reported.

The flow-induced vibrations of a single-layer vocal fold model were investigated as a function of vocal fold stiffness, and subglottal and supraglottal acoustic loading. Previously, it was reported that the single-layer vocal fold model failed to vibrate when short, clinically-relevant tracheal tubes were used. Moreover, it was reported that the model had a propensity to be acoustically driven, and aerodynamically driven vibration was observed only when a vertical restraint was applied superiorly to the vocal folds. However, in this study involving a wider range of source/tract conditions, the previous conclusions were shown to apply only for the special case of a stiff vocal fold model, for which self-oscillation occurred only when the vocal fold vibration synchronized to either a subglottal or supraglottal resonance. For a more general case, when vocal fold stiffness was decreased, the model did exhibit self-oscillation at short tracheal tubes, and no vertical restraint was needed to induce aerodynamically driven phonation. Nevertheless, the vocal fold vibration transitioned from aerodynamically-driven to acoustically-driven vibration when one of the subglottal resonance frequencies approximated one of the natural frequencies of the vocal folds. In this region, strong superior-inferior vibrations were observed, the phonation threshold pressure was significantly reduced, and the phonation onset frequency was heavily influenced by the dominant acoustic resonance. For acoustically-driven phonation, a compliant subglottal system always lowered phonation threshold. However, an inertive vocal tract could either increase or decrease phonation threshold pressure, depending on the phonation frequency.

Nonlinear dynamic analysis and model simulations are used to study the nonlinear dynamic characteristics of vocal folds with vocal tremor, which can typically be characterized by low frequency modulation and aperiodicity. Tremor voices from patients with disorders such as paresis, Parkinson's disease, hyperfunction, and adductor spasmodic dysphonia show low-dimensional characteristics, differing from random noise. Correlation dimension analysis statistically distinguishes tremor voices from normal voices. Furthermore, a nonlinear tremor model is proposed to study the vibrations of the vocal folds with vocal tremor. Fractal dimensions and positive Lyapunov exponents demonstrate the evidence of chaos in the tremor model, where amplitude and frequency play important roles in governing vocal fold dynamics. Nonlinear dynamic voice analysis and vocal fold modeling may provide a useful set of tools for understanding the dynamic mechanism of vocal tremor in patients with laryngeal diseases.