Non-linear elastic properties of the lingual and facial tissues assessed by indentation technique. Application to the biomechanics of speech production

Applied use of biomechanical measurements from human tissues for the development of medical skills trainers: a scoping review

OBJECTIVE

The objective of this review was to identify quantitative biomechanical measurements of human tissues, the methods for obtaining these measurements, and the primary motivations for conducting biomechanical research.

INTRODUCTION

Medical skills trainers are a safe and useful tool for clinicians to use when learning or practicing medical procedures. The haptic fidelity of these devices is often poor, which may be because the synthetic materials chosen for these devices do not have the same mechanical properties as human tissues. This review investigates a heterogeneous body of literature to identify which biomechanical properties are available for human tissues, the methods for obtaining these values, and the primary motivations behind conducting biomechanical tests.

INCLUSION CRITERIA

Studies containing quantitative measurements of the biomechanical properties of human tissues were included. Studies that primarily focused on dynamic and fluid mechanical properties were excluded. Additionally, studies only containing animal, in silico , or synthetic materials were excluded from this review.

METHODS

This scoping review followed the JBI methodology for scoping reviews and the Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews (PRISMA-ScR). Sources of evidence were extracted from CINAHL (EBSCO), IEEE Xplore, MEDLINE (PubMed), Scopus, and engineering conference proceedings. The search was limited to the English language. Two independent reviewers screened titles and abstracts as well as full-text reviews. Any conflicts that arose during screening and full-text review were mediated by a third reviewer. Data extraction was conducted by 2 independent reviewers and discrepancies were mediated through discussion. The results are presented in tabular, figure, and narrative formats.

RESULTS

Data were extracted from a total of 186 full-text publications. All of the studies, except for 1, were experimental. Included studies came from 33 countries, with the majority coming from the United States. Ex vivo methods were the predominant approach for extracting human tissue samples, and the most commonly studied tissue type was musculoskeletal. In this study, nearly 200 unique biomechanical values were reported, and the most commonly reported value was Young's (elastic) modulus. The most common type of mechanical test performed was tensile testing, and the most common reason for testing human tissues was to characterize biomechanical properties. Although the number of published studies on biomechanical properties of human tissues has increased over the past 20 years, there are many gaps in the literature. Of the 186 included studies, only 7 used human tissues for the design or validation of medical skills training devices. Furthermore, in studies where biomechanical values for human tissues have been obtained, a lack of standardization in engineering assumptions, methodologies, and tissue preparation may implicate the usefulness of these values.

CONCLUSIONS

This review is the first of its kind to give a broad overview of the biomechanics of human tissues in the published literature. With respect to high-fidelity haptics, there is a large gap in the published literature. Even in instances where biomechanical values are available, comparing or using these values is difficult. This is likely due to the lack of standardization in engineering assumptions, testing methodology, and reporting of the results. It is recommended that journals and experts in engineering fields conduct further research to investigate the feasibility of implementing reporting standards.

REVIEW REGISTRATION

Open Science Framework https://osf.io/fgb34.

MR-compatible loading device for assessment of heel pad internal tissue displacements under shearing load

In the last decade, the role of shearing loads has been increasingly suspected to play a determinant impact in the formation of deep pressure ulcers. In vivo observations of such deformations are complex to obtain. Previous studies only provide global measurements of such deformations without getting the quantitative values of the loads that generate these deformations. To study the role that shearing loads have in the etiology of heel pressure ulcers, an MR-compatible device for the application of shearing and normal loads was designed. Magnetic resonance imaging is a key feature that allows to monitor deformations of soft tissues after loading in a non-invasive way. Measuring applied forces in an MR-environment is challenging due to the impossibility to use magnetic materials. In our device, forces are applied through the compression of springs made of polylactide. Shearing and normal loads were applied on the plantar skin of the human heel through a flat plate while acquiring MR images. The device materials did not introduce any imaging artifact and allowed for high quality MR deformation measurements of the internal components of the heel. The obtained subject-specific results are an original data set that can be used in validations for Finite Element analysis and therefore contribute to a better understanding of the factors involved in pressure ulcer development.

Mechanical properties of whole-body soft human tissues: a review

The mechanical properties of soft tissues play a key role in studying human injuries and their mitigation strategies. While such properties are indispensable for computational modelling of biological systems, they serve as important references in loading and failure experiments, and also for the development of tissue simulants. To date, experimental studies have measured the mechanical properties of peripheral tissues (e.g. skin)in-vivoand limited internal tissuesex-vivoin cadavers (e.g. brain and the heart). The lack of knowledge on a majority of human tissues inhibit their study for applications ranging from surgical planning, ballistic testing, implantable medical device development, and the assessment of traumatic injuries. The purpose of this work is to overcome such challenges through an extensive review of the literature reporting the mechanical properties of whole-body soft tissues from head to toe. Specifically, the available linear mechanical properties of all human tissues were compiled. Non-linear biomechanical models were also introduced, and the soft human tissues characterized using such models were summarized. The literature gaps identified from this work will help future biomechanical studies on soft human tissue characterization and the development of accurate medical models for the study and mitigation of injuries.

A computational model of upper airway respiratory function with muscular coupling

A two dimensional finite element model of upper airway respiratory function was developed emphasizing the effects of dilator muscular activation on the human retro-lingual airway. The model utilized an upright mid-sagittal computed tomography of the human head and neck to reconstruct relevant structures of the tongue, mandible, and the hyoid-related soft tissues, along with the retro-lingual airway. The reconstructed geometry was divided into fluid and solid domains and discretized into finite element (FE) meshes used for the computational model. Three cases were investigated: standing position; supine position; and supine position coupled with dilator muscle activation. Computations were performed for the inspiration stage of the breathing cycle, utilizing a fluid-structure interaction (FSI) method to couple structural deformation with airflow dynamics. The spatio-temporal deformation of the structures surrounding the airway wall were predicted to be in general agreement with known changes from upright to supine posture on luminal opening, as well as the distribution of airflow. The model effectively captured the effects of muscular stimulation on the upper airway anatomical changes, the flow characteristics relevant to airway reduction in the supine position and airway enlargement with muscle activation. The smallest airway opening in the retro-lingual section is predicted to occur at the epiglottic region in all the three cases considered, an unexpected vulnerable location of airway obstruction. The model also predicted that hyoid displacement would be associated with recovery from airway collapse. This information may be useful for building more complex models relevant to mechanisms and clinical interventions for obstructive sleep apnea.

Personalized biomechanical tongue models based on diffusion-weighted MRI and validated using optical tracking of range of motion

For advanced tongue cancer, the choice between surgery and organ-sparing treatment is often dependent on the expected loss of tongue functionality after treatment. Biomechanical models might assist in this choice by simulating the post-treatment function loss. However, this function loss varies between patients and should, therefore, be predicted for each patient individually. In the present study, the goal was to better predict the postoperative range of motion (ROM) of the tongue by personalizing biomechanical models using diffusion-weighted MRI and constrained spherical deconvolution reconstructions of tongue muscle architecture. Diffusion-weighted MRI scans of ten healthy volunteers were obtained to reconstruct their tongue musculature, which were subsequently registered to a previously described population average or atlas. Using the displacement fields obtained from the registration, the segmented muscle fiber tracks from the atlas were morphed back to create personalized muscle fiber tracks. Finite element models were created from the fiber tracks of the atlas and those of the individual tongues. Via inverse simulation of a protruding, downward, left and right movement, the ROM of the tongue was predicted. This prediction was compared to the ROM measured with a 3D camera. It was demonstrated that biomechanical models with personalized muscles bundles are better in approaching the measured ROM than a generic model. However, to achieve this result a correction factor was needed to compensate for the small magnitude of motion of the model. Future versions of these models may have the potential to improve the estimation of function loss after treatment for advanced tongue cancer.

Machine-Learning based model order reduction of a biomechanical model of the human tongue

BACKGROUND AND OBJECTIVES

This paper presents the results of a Machine-Learning based Model Order Reduction (MOR) method applied to a complex 3D Finite Element (FE) biomechanical model of the human tongue, in order to create a Digital Twin Model (DTM) that enables real-time simulations. The DTM is designed for future inclusion in a computer assisted protocol for tongue surgery planning.

METHODS

The proposed method uses an "a posteriori" MOR that allows, from a limited number of simulations with the FE model, to predict in real time mechanical responses of the human tongue to muscle activations.

RESULTS

The MOR method is evaluated for simulations associated with separate single tongue muscle activations. It is shown to be able to account with a sub-millimetric spatial accuracy for the non-linear dynamical behavior of the tongue model observed in these simulations.

CONCLUSION

Further evaluations of the MOR method will include tongue movements induced by multiple muscle activations. At this stage our MOR method offers promising perspectives for the use of the tongue model in a clinical context to predict the impact of tongue surgery on tongue mobility. As a long term application, this DTM of the tongue could be used to predict the functional consequences of the surgery in terms of speech production and swallowing.

In-vivo tongue stiffness measured by aspiration: Resting vs general anesthesia

Tongue cancer treatment often results in impaired speech, swallowing, or mastication. Simulating the effect of treatments can help the patient and the treating physician to understand the effects and impact of the intervention. To simulate deformations of the tongue, identifying accurate mechanical properties of tissue is essential. However, not many succeeded in characterizing in-vivo tongue stiffness. Those who did, measured the tongue At Rest (AR), in which muscle tone subsides even if muscles are not willingly activated. We expected to find an absolute rest state in participants 'under General Anesthesia' (GA). We elaborated on previous work by measuring the mechanical behavior of the in-vivo tongue under aspiration using an improved volume-based method. Using this technique, 5 to 7 measurements were performed on 10 participants both AR and under GA. The obtained Pressure-Shape curves were first analyzed using the initial slope and its variations. Hereafter, an inverse Finite Element Analysis (FEA) was applied to identify the mechanical parameters using the Yeoh, Gent, and Ogden hyperelastic models. The measurements AR provided a mean Young's Modulus of 1638 Pa (min 1035 - max 2019) using the Yeoh constitutive model, which is in line with previous ex-vivo measurements. However, while hoping to find a rest state under GA, the tongue unexpectedly appeared to be approximately 2 to 2.5 times stiffer under GA than AR. Explanations for this were sought by examining drugs administered during GA, blood flow, perfusion, and upper airway reflexes, but neither of these explanations could be confirmed.

Quantal biomechanical effects in speech postures of the lips

The unique biomechanical and functional constraints on human speech make it a promising area for research investigating modular control of movement. The present article illustrates how a modular control approach to speech can provide insights relevant to understanding both motor control and observed variation across languages. We specifically explore the robust typological finding that languages produce different degrees of labial constriction using distinct muscle groupings and concomitantly distinct lip postures. Research has suggested that these lip postures exploit biomechanical regions of nonlinearity between neural activation and movement, also known as quantal regions, to allow movement goals to be realized despite variable activation signals. We present two sets of computer simulations showing that these labial postures can be generated under the assumption of modular control and that the corresponding modules are biomechanically robust: first to variation in the activation levels of participating muscles, and second to interference from surrounding muscles. These results provide support for the hypothesis that biomechanical robustness is an important factor in selecting the muscle groupings used for speech movements and provide insight into the neurological control of speech movements and how biomechanical and functional constraints govern the emergence of speech motor modules. We anticipate that future experimental work guided by biomechanical simulation results will provide new insights into the neural organization of speech movements.NEW & NOTEWORTHY This article provides additional evidence that speech motor control is organized in a modular fashion and that biomechanics constrain the kinds of motor modules that may emerge. It also suggests that speech can be a fruitful domain for the study of modularity and that a better understanding of speech motor modules will be useful for speech research. Finally, it suggests that biomechanical modeling can serve as a useful complement to experimental work when studying modularity.

Estimation of the hyperelastic parameters of fresh human oropharyngeal soft tissues using indentation testing

Patient-specific finite element (FE) modeling of the upper airway is an effective tool for accurate assessment of obstructive sleep apnea (OSA) syndrome. It is also useful for planning minimally invasive surgical procedures under severe OSA conditions. A major requirement of FE modeling is having reliable data characterizing the biomechanical properties of the upper airway tissues, particularly oropharyngeal soft tissue. While some data characterizing this tissue's linear elastic regime is available, reliable data characterizing its hyperelasticity is scarce. The aim of the current study is to estimate the hyperelastic mechanical properties of the oropharyngeal soft tissues, including the palatine tonsil, soft palate, uvula, and tongue base. Fresh tissue specimens of human oropharyngeal tissue were acquired from 13 OSA patients who underwent standard surgical procedures. Indentation testing was performed on the specimens to obtain their force-displacement data. To determine the specimens' hyperelastic parameters using these data, an inverse FE framework was utilized. In this work, the hyperelastic parameters corresponding to the commonly used Yeoh and 2nd order Ogden models were obtained. Both models captured the experimental force-displacement data of the tissue specimens reasonably accurately with mean errors of 11.65% or smaller. This study has provided estimates of the hyperelastic parameters of all upper airway soft tissues using fresh human tissue specimens for the first time.

Characterization of Material Properties Based on Inverse Finite Element Modelling

This paper describes a new approach that can be used to determine the mechanical properties of unknown materials and complex material systems. The approach uses inverse finite element modelling (FEM) accompanied with a designed algorithm to obtain the modulus of elasticity, yield stress and strain hardening material constants of an isotropic hardening material model, as well as the material constants of the Drucker–Prager material model (modulus of elasticity, cap yield stress and angle of friction). The algorithm automatically feeds the input material properties data to finite element software and automatically runs simulations to establish a convergence between the numerical loading–unloading curve and the target data obtained from continuous indentation tests using common indenter geometries. A further module was developed to optimise convergence using an inverse FEM analysis interfaced with a non-linear MATLAB algorithm. A sensitivity analysis determined that the dual spherical and Berkovich (S&B) approach delivered better results than other dual indentation methods such as Berkovich and Vickers (B&V) and Vickers and spherical (V&S). It was found that better convergence values can be achieved despite a large variation in the starting parameter values and/or material constitutive model and such behaviour reflects the uniqueness of the dual S&B indentation in predicting complex material systems. The study has shown that a robust optimization method based on a non-linear least-squares curve fitting function (LSQNONLIN) within MATLAB and ABAQUS can be used to accurately predict a unique set of elastic plastic material properties and Drucker–Prager material properties. This is of benefit to the scientific investigation of properties of new materials or obtaining the material properties at different locations of a part which may be not be similar because of manufacturing processes (e.g., different heating and cooling rates at different locations).

An interactive surgical simulation tool to assess the consequences of a partial glossectomy on a biomechanical model of the tongue

Oral cancer surgery has a negative influence on the quality of life (QOL). As a result of the complex physiology involved in oral functions, estimation of surgical effects on functionality remains difficult. We present a user-friendly biomechanical simulation of tongue surgery, including closure with suturing and scar formation, followed by an automated adaptation of a finite element (FE) model to the shape of the tongue. Different configurations of our FE model were evaluated and compared to a well-established FE model. We showed that the post-operative impairment as predicted by our model was qualitatively comparable to a patient case for five different tongue maneuvers.

Sagittal Measurement of Tongue Movement During Respiration: Comparison Between Ultrasonography and Magnetic Resonance Imaging

The tongue makes up the anterior pharyngeal wall and is critical for airway patency. Magnetic resonance imaging (MRI) is commonly used to study pharyngeal muscle function in pharyngeal disorders such as obstructive sleep apnoea. Tagged MRI and ultrasound studies have separately revealed ∼1 mm of anterior tongue movement during inspiration in healthy patients, but these modalities have not been directly compared. In the study described here, agreement between ultrasound and MRI in measuring regional tongue displacement in 21 healthy patients and 21 patients with obstructive sleep apnoea was evaluated. We found good consistency and agreement between the two techniques, with an intra-class correlation coefficient of 0.79 (95% confidence interval: 0.75-0.82) for anteroposterior tongue motion during inspiration. Ultrasound measurements of posterior tongue displacement were 0.24 ± 0.64 mm greater than MRI measurements (95% limits of agreement: 1.03 to -1.49). This may reflect the higher spatial and temporal resolution of the ultrasound technique. This study confirms that ultrasound is a suitable method for quantifying inspiratory tongue movement.

Mechanical behaviors of tension and relaxation of tongue and soft palate: Experimental and analytical modeling

This study is to characterize mechanical properties of uniaxial tension and stress relaxation responses of muscle tissues of tongue and soft palate. Uniaxial tension test and stress relaxation test on 39 fresh tissue samples from four five-month-old boars (65 ± 15 kg) are conducted. Firstly, the rationality of the samples' dimension design and experimenal data measurement is validated by one-way ANOVA F-type test. Mechanical properties, including stress-strain relationship and stress relaxation characteristic, are then investigated in details to show the nonlinear behaviors of the tissue samples clearly. Finally, a constitutive model of representing the mechanical properties is formulated within the nonlinear integral representation framework proposed by Pinkin and Rogers, and corresponding material parameters are fitted to the experimental data based on the Levenberg-Marquardt minimization algorithm. The results of the fitting comparison prove that the formulated constitutive model can capture the observed nonlinear behaviors of the muscle tissue samples in both the axial tension and stress relaxation experiments.

Estimation of the Young's moduli of fresh human oropharyngeal soft tissues using indentation testing

Finite element (FE)-based biomechanical simulations of the upper airway are promising computational tools to study abnormal upper airway deformations under obstructive sleep apnea (OSA) conditions and to help guide minimally invasive surgical interventions in case of upper airway collapse. To this end, passive biomechanical properties of the upper airway tissues, especially oropharyngeal soft tissues, are indispensable. This research aimed at characterizing the linear elastic mechanical properties of the oropharyngeal soft tissues including palatine tonsil, soft palate, uvula, and tongue base. For this purpose, precise indentation experiments were conducted on freshly harvested human tissue samples accompanied by FE-based inversion schemes. To minimize the impact of the probable nonlinearities of the tested tissue samples, only the first quarter of the measured force-displacement data corresponding to the linear elastic regime was utilized in the FE-based inversion scheme to improve the accuracy of the tissue samples' Young's modulus calculations. Measured Young's moduli of the oropharyngeal soft tissues obtained in this study are presented. They include first estimates for palatine tonsil tissue samples while measured Young's moduli of other upper airway tissues were obtained for the first time using fresh human tissue samples.

General theory of skin reinforcement

Macroscopic mechanical properties of human skin in vivo cannot be considered independent of adjacent subcutaneous white adipose tissue (sWAT). The layered system skin/sWAT appears as the hierarchical structural composite in which single layers behave as fiber-reinforced structures. Effective macroscopic mechanical properties of such composites are mainly determined either by the properties of the skin or by those of the sWAT, dependent on the conditions of mechanical loading. Mechanical interactions between the skin and the adjacent sWAT associated with a mismatch in the mechanical moduli of these two layers can lead to production of the skin wrinkles. Reinforcement of the composite skin/sWAT can take place in different ways. It can be provided through reorientation of collagen fibers under applied loading, through production of new bonds between existing collagen fibers and through induction of additional collagen structures. Effectiveness of this type of reinforcement is strongly dependent on the type of mechanical loading. Different physical interventions induce the reinforcement of at least one of these two layers, thus increasing the effective macroscopic stiffness of the total composite. At the same time, the standalone reinforcement of the skin appears to be less effective to achieve a delay or a reduction of the apparent signs of skin aging relative to the reinforcement of the sWAT.

Numerical simulation of interaction between organs and food bolus during swallowing and aspiration

The mechanism of swallowing is still not fully understood, because the process of swallowing is a rapid and complex interaction among several involved organs and the food bolus. In this work, with the aim of studying swallowing and aspiration processes noninvasively and systematically, a computer simulation method for analyzing the involved organs and water (considered as the food bolus) is proposed. The shape and motion of the organs involved in swallowing are modeled in the same way as in our previous study, by using the Hamiltonian moving particle simulation (MPS) method and forced displacements on the basis of motion in a healthy volunteer. The bolus flow is simulated using the explicit MPS method for fluid analysis. The interaction between the organs and the bolus is analyzed using a fluid-structure coupling scheme. To validate the proposed method, the behavior of the simulated bolus flow is compared qualitatively and quantitatively with corresponding medical images. In addition to the healthy motion model, disorder motion models are constructed for reproducing the aspiration phenomenon by computer simulation. The behaviors of the organs and the bolus considered as the food bolus in the healthy and disorder motion models are compared for evaluating the mechanism of aspiration.

Investigating skin penetration depth and shape following needle-free injection at different pressures: A cadaveric study

BACKGROUND AND OBJECTIVES

The effectiveness of needle-free injection devices in neocollagenesis for treating extended skin planes is an area of active research. It is anticipated that needle-free injection systems will not only be used to inject vaccines or insulin, but will also greatly aid skin rejuvenation when used to inject aesthetic materials such as hyaluronic acid, botulinum toxin, and placental extracts. There has not been any specific research to date examining how materials penetrate the skin when a needle-free injection device is used. In this study, we investigated how material infiltrates the skin when it is injected into a cadaver using a needle-free device.

STUDY DESIGN/MATERIALS AND METHODS

Using a needle-free injector (INNOJECTOR™; Amore Pacific, Seoul, Korea), 0.2 ml of 5% methylene blue (MB) or latex was injected into cheeks of human cadavers. The device has a nozzle diameter of 100 µm and produces a jet with velocity of 180 m/s. This jet penetrates the skin and delivers medicine intradermally via liquid propelled by compressed gasses. Materials were injected at pressures of 6 or 8.5 bars, and the injection areas were excised after the procedure. The excised areas were observed visually and with a phototrichogram to investigate the size, infiltration depth, and shape of the hole created on the skin. A small part of the area that was excised was magnified and stained with H&E (×40) for histological examination.

RESULTS

We characterized the shape, size, and depth of skin infiltration following injection of 5% MB or latex into cadaver cheeks using a needle-free injection device at various pressure settings. Under visual inspection, the injection at 6 bars created semi-circle-shaped hole that penetrated half the depth of the excised tissue, while injection at 8.5 bars created a cylinder-shaped hole that spanned the entire depth of the excised tissue. More specific measurements were collected using phototrichogram imaging. The shape of the injection entry point was consistently spherical regardless of the amount of pressure used. When injecting 5% MB at 6 bars, the depth of infiltration reached 2.323 mm, while that at 8.5 bars reached 8.906 mm. The area of the hole created by the 5% MB injection was 0.797 mm(2) at 6 bars and 0.242 mm(2) at 8.5 bars. Latex injections reached a depth of 3.480 mm at 6 bars and 7.558 mm at 8.5 bars, and the areas were measured at 1.043 mm(2) (6 bars) and 0.355 mm(2) (8.5 bars). Histological examination showed that the injection penetrated as deep as the superficial musculoaponeurotic system at 6 bars and the masseter muscle at 8.5 bars.

CONCLUSION

When injecting material into the skin using a pneumatic needle-free injector, higher-pressure injections result in a hole with smaller area than lower-pressure injections. The depth and shape of skin penetration vary according to the amount of pressure applied. For materials of low density and viscosity, there is a greater difference in penetration depth according to the degree of pressure. Lasers Surg. Med. 48:624-628, 2016. © 2016 Wiley Periodicals, Inc.

Development and validation of a computational finite element model of the rabbit upper airway: simulations of mandibular advancement and tracheal displacement

The mechanisms leading to upper airway (UA) collapse during sleep are complex and poorly understood. We previously developed an anesthetized rabbit model for studying UA physiology. On the basis of this body of physiological data, we aimed to develop and validate a two-dimensional (2D) computational finite element model (FEM) of the passive rabbit UA and peripharyngeal tissues. Model geometry was reconstructed from a midsagittal computed tomographic image of a representative New Zealand White rabbit, which included major soft (tongue, soft palate, constrictor muscles), cartilaginous (epiglottis, thyroid cartilage), and bony pharyngeal tissues (mandible, hard palate, hyoid bone). Other UA muscles were modeled as linear elastic connections. Initial boundary and contact definitions were defined from anatomy and material properties derived from the literature. Model parameters were optimized to physiological data sets associated with mandibular advancement (MA) and caudal tracheal displacement (TD), including hyoid displacement, which featured with both applied loads. The model was then validated against independent data sets involving combined MA and TD. Model outputs included UA lumen geometry, peripharyngeal tissue displacement, and stress and strain distributions. Simulated MA and TD resulted in UA enlargement and nonuniform increases in tissue displacement, and stress and strain. Model predictions closely agreed with experimental data for individually applied MA, TD, and their combination. We have developed and validated an FEM of the rabbit UA that predicts UA geometry and peripharyngeal tissue mechanical changes associated with interventions known to improve UA patency. The model has the potential to advance our understanding of UA physiology and peripharyngeal tissue mechanics.

Upper Airway Elasticity Estimation in Pediatric Down Syndrome Sleep Apnea Patients Using Collapsible Tube Theory

Elasticity of the soft tissues surrounding the upper airway lumen is one of the important factors contributing to upper airway disorders such as snoring and obstructive sleep apnea. The objective of this study is to calculate patient specific elasticity of the pharynx from magnetic resonance (MR) images using a 'tube law', i.e., the relationship between airway cross-sectional area and transmural pressure difference. MR imaging was performed under anesthesia in children with Down syndrome (DS) and obstructive sleep apnea (OSA). An airway segmentation algorithm was employed to evaluate changes in airway cross-sectional area dilated by continuous positive airway pressure (CPAP). A pressure-area relation was used to make localized estimates of airway wall stiffness for each patient. Optimized values of patient specific Young's modulus for tissue in the velopharynx and oropharynx, were estimated from finite element simulations of airway collapse. Patient specific deformation of the airway wall under CPAP was found to exhibit either a non-linear 'hardening' or 'softening' behavior. The localized airway and tissue elasticity were found to increase with increasing severity of OSA. Elasticity based patient phenotyping can potentially assist clinicians in decision making on CPAP and airway or tissue elasticity can supplement well-known clinical measures of OSA severity.

Full-surface deformation measurement of anisotropic tissues under indentation

Inverse finite element-based analysis of soft biological tissues is an important tool to investigate their complex mechanical behavior and to develop physical models for medical simulations. Although there have recently been advances in dealing with the computational complexities of modeling biological materials, the collection of a sufficiently dense set of experimental data to properly capture their typically regionally varying properties still remains a critical issue. The aim of this work was to develop and test an optical system that combines 2D-Digital Image Correlation (DIC) and a novel Fringe Projection method with radial sensitivity (RFP) to test soft biological tissues under in vitro indentation. This system has the distinctive capability of using a single camera to retrieve the shape and 3D deformation of the whole upper surface of the indented sample without any blind measurement areas (with exception of the area under the indenter), with nominal depth and in-plane resolution of 0.05 mm and 0.004 mm, respectively. To test and illustrate the capabilities of the developed DIC/RFP system, the in vitro response to indentation of a homogeneous and isotropic latex foam is presented against the response of a slab of porcine ventricular myocardium, a highly in-homogeneous and anisotropic tissue. Our results illustrate the enhanced capabilities of the developed method to capture asymmetry in deformation with respect to standard indentation tests. This feature, together with the possibility of miniaturizing the system into a hand-held probe, makes this method potentially extendable to in vivo settings, alone or in combination with ultrasound measurements.

Linear elastic properties of the facial soft tissues using an aspiration device: towards patient specific characterization

Biomechanical modeling of the facial soft tissue behavior is needed in aesthetic or maxillo-facial surgeries where the simulation of the bone displacements cannot accurately predict the visible outcome on the patient's face. Because these tissues have different nature and elastic properties across the face, depending on their thickness, and their content in fat or muscle, individualizing their mechanical parameters could increase the simulation accuracy. Using a specifically designed aspiration device, the facial soft tissues deformation is measured at four different locations (cheek, cheekbone, forehead, and lower lip) on 16 young subjects. The stiffness is estimated from the deformations generated by a set of negative pressures using an inverse analysis based on a Neo Hookean model. The initial Young's modulus of the cheek, cheekbone, forehead, and lower lip are respectively estimated to be 31.0 kPa±4.6, 34.9 kPa±6.6, 17.3 kPa±4.1, and 33.7 kPa±7.3. Significant intra-subject differences in tissue stiffness are highlighted by these estimations. They also show important inter-subject variability for some locations even when mean stiffness values show no statistical difference. This study stresses the importance of using a measurement device capable of evaluating the patient specific tissue stiffness during an intervention.

Biomechanical properties of the human upper airway and their effect on its behavior during breathing and in obstructive sleep apnea

The upper airway is a complex, multifunctional, dynamic neuromechanical system. Its patency during breathing requires moment-to-moment coordination of neural and mechanical behavior and varies with posture. Failure to continuously recruit and coordinate dilator muscles to counterbalance the forces that act to close the airway results in hypopneas or apneas. Repeated failures lead to obstructive sleep apnea (OSA). Obesity and anatomical variations, such as retrognathia, increase the likelihood of upper airway collapse by altering the passive mechanical behavior of the upper airway. This behavior depends on the mechanical properties of each upper airway tissue in isolation, their geometrical arrangements, and their physiological interactions. Recent measurements of respiratory-related deformation of the airway wall have shown that there are different patterns of airway soft tissue movement during the respiratory cycle. In OSA patients, airway dilation appears less coordinated compared with that in healthy subjects (matched for body mass index). Intrinsic mechanical properties of airway tissues are altered in OSA patients, but the factors underlying these changes have yet to be elucidated. How neural drive to the airway dilators relates to the biomechanical behavior of the upper airway (movement and stiffness) is still poorly understood. Recent studies have highlighted that the biomechanical behavior of the upper airway cannot be simply predicted from electromyographic activity (electromyogram) of its muscles.

The distributed lambda (λ) model (DLM): a 3-D, finite-element muscle model based on Feldman's λ model; assessment of orofacial gestures

PURPOSE

The authors aimed to design a distributed lambda model (DLM), which is well adapted to implement three-dimensional (3-D), finite-element descriptions of muscles.

METHOD

A muscle element model was designed. Its stress-strain relationships included the active force-length characteristics of the λ model along the muscle fibers, together with the passive properties of muscle tissues in the 3-D space. The muscle element was first assessed using simple geometrical representations of muscles in the form of rectangular bars. It was then included in a 3-D face model, and its impact on lip protrusion was compared with the impact of a Hill-type muscle model.

RESULTS

The force-length characteristic associated with the muscle elements matched well with the invariant characteristics of the λ model. The impact of the passive properties was assessed. Isometric force variation and isotonic displacements were modeled. The comparison with a Hill-type model revealed strong similarities in terms of global stress and strain.

CONCLUSION

The DLM accounted for the characteristics of the λ model. Biomechanically, no clear differences were found between the DLM and a Hill-type model. Accurate evaluations of the λ model, based on the comparison between data and simulations, are now possible with 3-D biomechanical descriptions of the speech articulators because of the DLM.

Automatic prediction of tongue muscle activations using a finite element model

Computational modeling has improved our understanding of how muscle forces are coordinated to generate movement in musculoskeletal systems. Muscular-hydrostat systems, such as the human tongue, involve very different biomechanics than musculoskeletal systems, and modeling efforts to date have been limited by the high computational complexity of representing continuum-mechanics. In this study, we developed a computationally efficient tracking-based algorithm for prediction of muscle activations during dynamic 3D finite element simulations. The formulation uses a local quadratic-programming problem at each simulation time-step to find a set of muscle activations that generated target deformations and movements in finite element muscular-hydrostat models. We applied the technique to a 3D finite element tongue model for protrusive and bending movements. Predicted muscle activations were consistent with experimental recordings of tongue strain and electromyography. Upward tongue bending was achieved by recruitment of the superior longitudinal sheath muscle, which is consistent with muscular-hydrostat theory. Lateral tongue bending, however, required recruitment of contralateral transverse and vertical muscles in addition to the ipsilateral margins of the superior longitudinal muscle, which is a new proposition for tongue muscle coordination. Our simulation framework provides a new computational tool for systematic analysis of muscle forces in continuum-mechanics models that is complementary to experimental data and shows promise for eliciting a deeper understanding of human tongue function.

A critical evaluation of gestural stiffness estimations in speech production based on a linear second-order model

PURPOSE

Linear second-order models have often been used to investigate properties of speech production. However, these models are inaccurate approximations of the speech apparatus. This study aims at assessing how reliably stiffness can be estimated from kinematics with these models.

METHOD

Articulatory movements were collected for 9 speakers of German during the production of reiterant CVCV words at varying speech rates. Velocity peaks, movement amplitudes, and gesture durations were measured. In the context of an undamped model, 2 stiffness estimations were compared that should theoretically yield the same result. In the context of a damped model, gestural stiffness and damping were calculated for each gesture.

RESULTS

Numerous cases were found in which stiffness estimations based on the undamped model contradicted each other. Less than 80% of the data were found to be compatible with the properties of the damped model. Stiffness tends to decrease with gestural duration. However, it is associated with a large, unrealistic damping dispersion, making stiffness estimations from kinematic data to a large extent unreliable.

CONCLUSION

Any conclusions about speech control based on stiffness estimations using linear second-order models should therefore be considered with caution.

Assessment of nose protector for sport activities: finite element analysis

There has been a significant increase in the number of facial fractures stemming from sport activities in recent years, with the nasal bone one of the most affected structures. Researchers recommend the use of a nose protector, but there is no standardization regarding the material employed. Clinical experience has demonstrated that a combination of a flexible and rigid layer of ethylene vinyl acetate (EVA) offers both comfort and safety to practitioners of sports. The aim of the present study was the investigation into the stresses generated by the impact of a rigid body on the nasal bone on models with and without an EVA protector. For such, finite element analysis was employed. A craniofacial model was constructed from images obtained through computed tomography. The nose protector was modeled with two layers of EVA (1 mm of rigid EVA over 2 mm of flexible EVA), following the geometry of the soft tissue. Finite element analysis was performed using the LS Dyna program. The bone and rigid EVA were represented as elastic linear material, whereas the soft tissues and flexible EVA were represented as hyperelastic material. The impact from a rigid sphere on the frontal region of the face was simulated with a constant velocity of 20 m s(-1) for 9.1 μs. The model without the protector served as the control. The distribution of maximal stress of the facial bones was recorded. The maximal stress on the nasal bone surpassed the breaking limit of 0.13-0.34 MPa on the model without a protector, while remaining below this limit on the model with the protector. Thus, the nose protector made from both flexible and rigid EVA proved effective at protecting the nasal bones under high-impact conditions.

Simulation of dynamic orofacial movements using a constitutive law varying with muscle activation

This paper presents a biomechanical model of the face to simulate orofacial movements in speech and non-verbal communication. A 3D finite element model, based on medical images of a subject, is presented. A hyperelastic Mooney-Rivlin constitutive law accounts for the non-linear behaviour of facial tissue. Muscle fibres are represented by piece-wise uniaxial tensile element that generate force. The stress stiffening effect, an increase in the stiffness of the muscles when activated, is modelled by varying the constitutive law of the tissue with the level of activation of the muscle. A large number of facial movements occurring during speech and facial mimics are simulated. Results show that our modelling approach provides a realistic account of facial mimics. The differences between dynamic vs. quasi-static simulations are also discussed, proving that dynamic trajectories better fit experimental data.

A biomechanical model of cardinal vowel production: muscle activations and the impact of gravity on tongue positioning

A three-dimensional (3D) biomechanical model of the tongue and the oral cavity, controlled by a functional model of muscle force generation (lambda-model of the equilibrium point hypothesis) and coupled with an acoustic model, was exploited to study the activation of the tongue and mouth floor muscles during the production of French cardinal vowels. The selection of the motor commands to control the tongue and the mouth floor muscles was based on literature data, such as electromyographic, electropalatographic, and cineradiographic data. The tongue shapes were also compared to data obtained from the speaker used to build the model. 3D modeling offered the opportunity to investigate the role of the transversalis, in particular, its involvement in the production of high front vowels. It was found, with this model, to be indirect via reflex mechanisms due to the activation of surrounding muscles, not voluntary. For vowel /i/, local motor command variations for the main tongue muscles revealed a non-negligible modification of the alveolar groove in contradiction to the saturation effect hypothesis, due to the role of the anterior genioglossus. Finally, the impact of subject position (supine or upright) on the production of French cardinal vowels was explored and found to be negligible.

Development and validation of a three-dimensional finite element model of the face

A detailed three-dimensional finite element model of the face is presented in this paper. Bones, muscles, skin, fat, and superficial muscoloaponeurotic system were reconstructed from magnetic resonance images and modeled according to anatomical, plastic, and reconstructive surgery literature. The finite element mesh, composed of hexahedron elements, was generated through a semi-automatic procedure with an effective compromise between the detailed representation of anatomical parts and the limitation of the computational time. Nonlinear constitutive equations are implemented in the finite element model. The corresponding model parameters were selected according to previous work with mechanical measurements on soft facial tissue, or based on reasonable assumptions. Model assumptions concerning tissue geometry, interactions, mechanical properties, and the boundary conditions were validated through comparison with experiments. The calculated response of facial tissues to gravity loads, to the application of a pressure inside the oral cavity and to the application of an imposed displacement was shown to be in good agreement with the data from corresponding magnetic resonance images and holographic measurements. As a first application, gravimetric soft tissue descent was calculated from the long time action of gravity on the face in the erect position, with tissue aging leading to a loss of stiffness. Aging predictions are compared with the observations from an "aging database" with frontal photos of volunteers at different age ranges (i.e., 20-40 years and 50-70 years).

Polyvinyl alcohol cryogel: optimizing the parameters of cryogenic treatment using hyperelastic models

The PVA gels obtained by freezing/thawing cycles of PVA solutions, also called cryogels, exhibit non-linear elastic behavior and can mimic, within certain limits, the behavior of biological soft tissues such as arterial tissue. Several authors have investigated the effects of cryogenic processing parameters on the Young's modulus. However, an elastic modulus does not describe the non-linearity of the cryogel's stress-strain response. This study examines the non-linear elastic response of PVA cryogel under uniaxial tension and investigates how processing parameters such as the concentration, the number of thermal cycles, and the thawing rate affect this response. The relationship between the coefficients of the material model and the processing parameters was interpolated to find the set of parameters that would best approximate the elastic response of healthy porcine coronary arteries under uniaxial tension.

In situ measurement and modeling of biomechanical response of human cadaveric soft tissues for physics-based surgical simulation

BACKGROUND

Development of a laparoscopic surgery simulator that delivers high-fidelity visual and haptic (force) feedback, based on the physical models of soft tissues, requires the use of empirical data on the mechanical behavior of intra-abdominal organs under the action of external forces. As experiments on live human patients present significant risks, the use of cadavers presents an alternative. We present techniques of measuring and modeling the mechanical response of human cadaveric tissue for the purpose of developing a realistic model. The major contribution of this paper is the development of physics-based models of soft tissues that range from linear elastic models to nonlinear viscoelastic models which are efficient for application within the framework of a real-time surgery simulator.

METHODS

To investigate the in situ mechanical, static, and dynamic properties of intra-abdominal organs, we have developed a high-precision instrument by retrofitting a robotic device from Sensable Technologies (position resolution of 0.03 mm) with a six-axis Nano 17 force-torque sensor from ATI Industrial Automation (force resolution of 1/1,280 N along each axis), and used it to apply precise displacement stimuli and record the force response of liver and stomach of ten fresh human cadavers.

RESULTS

The mean elastic modulus of liver and stomach is estimated as 5.9359 kPa and 1.9119 kPa, respectively over the range of indentation depths tested. We have also obtained the parameters of a quasilinear viscoelastic (QLV) model to represent the nonlinear viscoelastic behavior of the cadaver stomach and liver over a range of indentation depths and speeds. The models are found to have an excellent goodness of fit (with R (2) > 0.99).

CONCLUSIONS

The data and models presented in this paper together with additional ones based on the principles presented in this paper would result in realistic physics-based surgical simulators.

Speed–Curvature Relations in Speech Production Challenge the 1/3 Power Law

Relations between tangential velocity and trajectory curvature are analyzed for tongue movements during speech production in the framework of the 1/3 power law, discovered by Viviani and colleagues for arm movements. In 2004, Tasko and Westbury found for American English that the power function provides a good account of speech kinematics, but with an exponent that varies across articulators. The present work aims at broadening Tasko and Westbury's study 1) by analyzing speed–curvature relations for various languages (French, German, Mandarin) and for a biomechanical tongue model simulating speech gestures at various speaking rates and 2) by providing for each speaker or each simulated speaking rate a comparison of results found for the complete set of movements with those found for each movement separately. It is found that the 1/3 power law offers a fair description of the global speed–curvature relations for all speakers and all languages, when articulatory speech data are considered in their whole. This is also observed in the simulations, where the motor control model does not specify any kinematic property of the articulatory paths. However, the refined analysis for individual movements reveals numerous exceptions to this law: the velocity always decreases when curvature increases, but the slope in the log–log representation is variable. It is concluded that the speed–curvature relation is not controlled in speech movements and that it accounts only for general properties of the articulatory movements, which could arise from vocal tract dynamics or/and from stochastic characteristics of the measured signals.

A light sterilizable pipette device for the in vivo estimation of human soft tissues constitutive laws

This paper introduces a new light device for the in vivo estimation of human soft tissues constitutive laws. It consists of an aspiration pipette able to meet the very severe sterilization and handling issues imposed during surgery. The simplicity of the device, free of any electronic circuitry, allows using it as an ancillary instrument. The deformation of the aspired tissue is imaged via a mirror using an external camera. The paper describes the experimental setup as well as the protocol that should be used during surgery. First feasibility measurements are shown for human tongue and forearm skin.

Simulations of the consequences of tongue surgery on tongue mobility: implications for speech production in post-surgery conditions

BACKGROUND

We studied the ability of a three-dimensional (3D) biomechanical model of the oral cavity to predict the consequences of tongue surgery on tongue movements, according to the size and location of the tissue loss and the nature of the flap used by the surgeon.

METHOD

The core of our model consists of a 3D biomechanical model representing the tongue as a finite element structure with hexahedral elements and hyperelastic properties, in which muscles are represented by specific subsets of elements. This model is inserted in the oral cavity including jaw, palate and pharyngeal walls. Hemiglossectomy and large resection of the mouth floor are simulated by removing the elements corresponding to the tissue losses. Three kinds of reconstruction are modelled, assuming flaps with low, medium or high stiffness.

RESULTS

The consequences of these different surgical treatments during the activation of some of the main tongue muscles are shown. Differences in global 3D tongue shape and in velocity patterns are evaluated and interpreted in terms of their potential impact on speech articulation. These simulations have been shown to be efficient in accounting for some of the clinically observed consequences of tongue surgery.

CONCLUSION

Further improvements still need to be done before being able to generate patient-specific models easily and to decrease the computation time significantly. However, this approach should represent a significant improvement in planning tongue surgery systems and should be a very useful means of improving the understanding of muscle behaviour after partial resection.

Intra-operative quantification of the surgical gesture in orbital surgery: application to the proptosis reduction

BACKGROUND

Proptosis is characterized by a protrusion of the eyeball due to an increase of the orbital tissue volume. To recover a normal eyeball positioning, the most frequent surgical technique consists in the osteotomy of orbital walls combined with a loading on the eyeball to initiate tissue decompression. The first biomechanical models dealing with proptosis reduction, validated in one patient, have been previously proposed by the authors.

METHODS

This paper proposed an experimental method to quantify the intra-operative clinical gesture in proptosis reduction, and the pilot study concerned one clinical case. The eyeball's backward displacement was measured by an optical 3D localizer and the load applied by the surgeon was simultaneously measured by a custom-made force gauge. Quasi-static stiffness of the intra-orbital content was evaluated.

FINDINGS

The average values for the whole experiment was 16 N (SD: 3N) for the force exerted by the surgeon and 9 mm (SD: 4mm) for the eyeball backward displacement. The averaged quasi-static stiffness of the orbital content was evaluated to 2.4N/mm (SD: 1.2) and showed a global decrease of 45% post-operatively.

INTERPRETATION

The protocol and the associated custom-designed devices allowed loads, induced displacements and macroscopic stiffness of the orbital content to be measured intra-operatively. The clinical relevance has been demonstrated in a pilot study. To our knowledge, no study has been published allowing the clinical gesture in proptosis reduction to be quantified intra-operatively. Associating an enlarged database and validated patient-related predictive models will reinforce the surgical efficiency and patient comfort contributing to diagnosis and intra-operative guidance.