The effect of Glisson's capsule on the superficial elasticity measurements of the liver

In the past decade, novel tools for surgical planning and disease diagnosis have been developed to detect the liver's mechanical properties. Some tools utilize superficial indentation type measurements to determine the elasticity of the liver parenchyma and to assume material homogeneity. In fact, the liver is a soft tissue covered with a connective sheathing that is called Glisson's capsule. This article quantifies the effect of this capsule on the measured or "effective" elastic modulus obtained by indentation with a spherical geometry. Two sets of parametric computational studies in which the Glisson capsule thickness and elasticity were varied, demonstrated the relationship between the measured elastic response and the underlying parenchymal elastic response. Previously reported in vivo indentation data on the human liver were utilized to determine the elasticity of its parenchyma. The results indicated a linear relationship between the effective (measured) elastic response and the underlying parenchyma for the Glisson capsule thicknesses considered. When previously published human liver indentation data were analyzed, the measured elastic modulus was nearly 6.9% greater than the underlying parenchyma elastic modulus. Although the analyzed data were obtained from a single liver and yet displayed a significant variation, the Glisson capsule elasticity induced a significant but systematic error as well. The Glisson capsule thickness error was negligible for capsule parameters associated with a normal liver. Based on this work, an emphasis on the Glisson capsule's contribution to the mechanical response of the liver would enhance the clinical potential of indentation-based novel tools for liver care.     

In vivo liver tissue mechanical properties by Transient Elastography: comparison with Dynamic Mechanical Analysis

Understanding the mechanical properties of human liver is one of the most critical aspects of its numerical modeling for medical applications or impact biomechanics. Generally, model constitutive laws come from in vitro data. However, the elastic properties of liver may change significantly after death and with time. Furthermore, in vitro liver elastic properties reported in the literature have often not been compared quantitatively with in vivo liver mechanical properties on the same organ. In this study, both steps are investigated on porcine liver. The elastic property of the porcine liver, given by the shear modulus G, was measured by both Transient Elastography (TE) and Dynamic Mechanical Analysis (DMA). Shear modulus measurements were realized on in vivo and in vitro liver to compare the TE and DMA methods and to study the influence of testing conditions on the liver viscoelastic properties. In vitro results show that elastic properties obtained by TE and DMA are in agreement. Liver tissue in the frequency range from 0.1 to 4 Hz can be modeled by a two-mode relaxation model. Furthermore, results show that the liver is homogeneous, isotropic and more elastic than viscous. Finally, it is shown in this study that viscoelastic properties obtained by TE and DMA change significantly with post mortem time and with the boundary conditions.     

Modelling of global boundary effects on harmonic motion imaging of soft tissues

Biomechanical imaging techniques have been developed for soft tissue characterisation and detection of breast tumours. Harmonic motion imaging (HMI) uses a focused ultrasound technology to generate a harmonic radiation force in a localised region inside a soft tissue. The resulting dynamic response is used to map the local distribution of the mechanical properties of the tissue. In this study, a finite element (FE) model is developed to investigate the effect of global boundary conditions on the dynamic response of a soft tissue during HMI. The direct-solution steady-state dynamic analysis procedure is used to compute the harmonic displacement amplitude in FE simulations. The model is parameterised in terms of boundary conditions and viscoelastic properties, and the corresponding raster-scan displacement amplitudes are captured to examine its response. The effect of the model's global dimensions on the harmonic response is also investigated. It is observed that the dynamic response of soft tissue with high viscosity is independent of the global boundary conditions for regions remote to the boundary; thus, it can be subjected to local analysis to estimate the underlying mechanical properties. However, the dynamic response is sensitive to global boundary conditions for tissue with low viscosity or regions located near to the boundary.     

Mechano-acoustic determination of Young's modulus of articular cartilage

The compressive stiffness of an elastic material is traditionally characterized by its Young's modulus. Young's modulus of articular cartilage can be directly measured using unconfined compression geometry by assuming the cartilage to be homogeneous and isotropic. In isotropic materials, Young's modulus can also be determined acoustically by the measurement of sound speed and density of the material. In the present study, acoustic and mechanical techniques, feasible for in vivo measurements, were investigated to quantify the static and dynamic compressive stiffness of bovine articular cartilage in situ. Ultrasound reflection from the cartilage surface, as well as the dynamic modulus were determined with the recently developed ultrasound indentation instrument and compared with the reference mechanical and ultrasound speed measurements in unconfined compression (n=72). In addition, the applicability of manual creep measurements with the ultrasound indentation instrument was evaluated both experimentally and numerically. Our experimental results indicated that the sound speed could predict 47% and 53% of the variation in the Young's modulus and dynamic modulus of cartilage, respectively. The dynamic modulus, as determined manually with the ultrasound indentation instrument, showed significant linear correlations with the reference Young's modulus (r(2)=0.445, p<0.01, n=70) and dynamic modulus (r(2)=0.779, p<0.01, n=70) of the cartilage. Numerical analyses indicated that the creep measurements, conducted manually with the ultrasound indentation instrument, were sensitive to changes in Young's modulus and permeability of the tissue, and were significantly influenced by the tissue thickness. We conclude that acoustic parameters, i.e. ultrasound speed and reflection, are indicative to the intrinsic mechanical properties of the articular cartilage. Ultrasound indentation instrument, when further developed, provides an applicable tool for the in vivo detection of cartilage mechano-acoustic properties. These techniques could promote the diagnostics of osteoarthrosis.     

Measurement of the dynamic shear modulus of mouse brain tissue in vivo by magnetic resonance elastography

In this study, the magnetic resonance (MR) elastography technique was used to estimate the dynamic shear modulus of mouse brain tissue in vivo. The technique allows visualization and measurement of mechanical shear waves excited by lateral vibration of the skull. Quantitative measurements of displacement in three dimensions during vibration at 1200 Hz were obtained by applying oscillatory magnetic field gradients at the same frequency during a MR imaging sequence. Contrast in the resulting phase images of the mouse brain is proportional to displacement. To obtain estimates of shear modulus, measured displacement fields were fitted to the shear wave equation. Validation of the procedure was performed on gel characterized by independent rheometry tests and on data from finite element simulations. Brain tissue is, in reality, viscoelastic and nonlinear. The current estimates of dynamic shear modulus are strictly relevant only to small oscillations at a specific frequency, but these estimates may be obtained at high frequencies (and thus high deformation rates), noninvasively throughout the brain. These data complement measurements of nonlinear viscoelastic properties obtained by others at slower rates, either ex vivo or invasively.     

Estimation of material elastic moduli in elastography: a local method, and an investigation of Poisson's ratio sensitivity

The local material stiffness of tissues is a well-known indicator of pathology, with locally stiffer tissue related to the possible presence of an abnormal growth in otherwise compliant tissue. Elastography is a non-invasive technique for measuring displacement distributions in loaded tissues within a medical imaging context. From these measured displacement fields, estimated for local strain have been made using well-studied techniques, but the calculation of elastic modulus has been difficult. In this study we show a method for estimating local tissue elastic modulus that gives numerically stable and robust results in test cases, and that is numerically efficient. The method assumes the tissue is isotropic and it requires an independent estimate of tissue Poisson's ratio, but the method reaches a stable result when the estimated Poisson's ratio is in error, and the resulting estimates are not very sensitive to the assumed value.     

High-resolution finite element models with tissue strength asymmetry accurately predict failure of trabecular bone

The ability to predict trabecular failure using microstructure-based computational models would greatly facilitate study of trabecular structure–function relations, multiaxial strength, and tissue remodeling. We hypothesized that high-resolution finite element models of trabecular bone that include cortical-like strength asymmetry at the tissue level, could predict apparent level failure of trabecular bone for multiple loading modes. A bilinear constitutive model with asymmetric tissue yield strains in tension and compression was applied to simulate failure in high-resolution finite element models of seven bovine tibial specimens. Tissue modulus was reduced by 95% when tissue principal strains exceeded the tissue yield strains. Linear models were first calibrated for effective tissue modulus against specimen-specific experimental measures of apparent modulus, producing effective tissue moduli of (mean±S.D.) 18.7±3.4 GPa. Next, a parameter study was performed on a single specimen to estimate the tissue level tensile and compressive yield strains. These values, 0.60% strain in tension and 1.01% strain in compression, were then used in non-linear analyses of all seven specimens to predict failure for apparent tensile, compressive, and shear loading. When compared to apparent yield properties previously measured for the same type of bone, the model predictions of both the stresses and strains at failure were not statistically different for any loading case (p>0.15). Use of symmetric tissue strengths could not match the experimental data. These findings establish that, once effective tissue modulus is calibrated and uniform but asymmetric tissue failure strains are used, the resulting models can capture the apparent strength behavior to an outstanding level of accuracy. As such, these computational models have reached a level of fidelity that qualifies them as surrogates for destructive mechanical testing of real specimens.     

Effect of preservation period on the viscoelastic material properties of soft tissues with implications for liver transplantation

The liver harvested from a donor must be preserved and transported to a suitable recipient immediately for a successful liver transplantation. In this process, the preservation period is the most critical, since it is the longest and most tissue damage occurs during this period due to the reduced blood supply to the harvested liver and the change in its temperature. We investigate the effect of preservation period on the dynamic material properties of bovine liver using a viscoelastic model derived from both impact and ramp and hold experiments. First, we measure the storage and loss moduli of bovine liver as a function of excitation frequency using an impact hammer. Second, its time-dependent relaxation modulus is measured separately through ramp and hold experiments performed by a compression device. Third, a Maxwell solid model that successfully imitates the frequency- and time-dependent dynamic responses of bovine liver is developed to estimate the optimum viscoelastic material coefficients by minimizing the error between the experimental data and the corresponding values generated by the model. Finally, the variation in the viscoelastic material coefficients of bovine liver are investigated as a function of preservation period for the liver samples tested 1 h, 2 h, 4 h, 8 h, 12 h, 24 h, 36 h, and 48 h after harvesting. The results of our experiments performed with three animals show that the liver tissue becomes stiffer and more viscous as it spends more time in the preservation cycle.     

Large deformation shear properties of liver tissue

Characterization of the mechanical properties of soft biological tissues is important for establishing the mechanical tolerances of the tissues, and for input to computational models. In this work, the viscoelastic properties of bovine liver tissue in shear loading have been measured using relaxation and constant shear rate loading. The tissue is nonlinearly viscoelastic for strains greater than 0.2%, has a yield strain of approximately 10, and shows moderate strain-rate sensitivity. The response can be modelled using a nonlinear viscoelastic differential model previously developed for brain tissue.     

A two-layered mechanical model of the rat esophagus. Experiment and theory

Background  

The function of esophagus is to move food by peristaltic motion which is the result of the interaction of the tissue forces in the esophageal wall and the hydrodynamic forces in the food bolus. The structure of the esophagus is layered. In this paper, the esophagus is treated as a two-layered structure consisting of an inner collagen-rich submucosa layer and an outer muscle layer. We developed a model and experimental setup for determination of elastic moduli in the two layers in circumferential direction and related the measured elastic modulus of the intact esophagus to the elastic modulus computed from the elastic moduli of the two layers.     

Characterization of vascular strain during in-vitro angioplasty with high-resolution ultrasound speckle tracking

Background

Ultrasound elasticity imaging provides biomechanical and elastic properties of vascular tissue, with the potential to distinguish between tissue motion and tissue strain. To validate the ability of ultrasound elasticity imaging to predict structurally defined physical changes in tissue, strain measurement patterns during angioplasty in four bovine carotid artery pathology samples were compared to the measured physical characteristics of the tissue specimens.

Methods

Using computational image-processing techniques, the circumferences of each bovine artery specimen were obtained from ultrasound and pathologic data.

Results

Ultrasound-strain-based and pathology-based arterial circumference measurements were correlated with an R² value of 0.94 (p = 0.03). The experimental elasticity imaging results confirmed the onset of deformation of an angioplasty procedure by indicating a consistent inflection point where vessel fibers were fully unfolded and vessel wall strain initiated.

Conclusion

These results validate the ability of ultrasound elasticity imaging to measure localized mechanical changes in vascular tissue.     

An information-based machine learning approach to elasticity imaging

An information-based technique is described for applications in mechanical property imaging of soft biological media under quasi-static loads. We adapted the Autoprogressive method that was originally developed for civil engineering applications for this purpose. The Autoprogressive method is a computational technique that combines knowledge of object shape and a sparse distribution of force and displacement measurements with finite-element analyses and artificial neural networks to estimate a complete set of stress and strain vectors. Elasticity imaging parameters are then computed from estimated stresses and strains. We introduce the technique using ultrasonic pulse-echo measurements in simple gelatin imaging phantoms having linear-elastic properties so that conventional finite-element modeling can be used to validate results. The Autoprogressive algorithm does not require any assumptions about the material properties and can, in principle, be used to image media with arbitrary properties. We show that by selecting a few well-chosen force–displacement measurements that are appropriately applied during training and establish convergence, we can estimate all nontrivial stress and strain vectors throughout an object and accurately estimate an elastic modulus at high spatial resolution. This new method of modeling the mechanical properties of tissue-like materials introduces a unique method of solving the inverse problem and is the first technique for imaging stress without assuming the underlying constitutive model.     

Characterization of bone-implant fixation using modal analysis: Application to a press-fit implant model

Orthopaedic implant fixation is strongly dependant upon the effective mechanical properties of newly formed tissue. In this study, we evaluated the potential of modal analysis to derive viscoelastic properties of periprosthetic tissue. We hypothesized that Young's modulus and loss factor could be obtained by a combined theoretical, computational and experimental modal analysis approach. This procedure was applied to ex vivo specimens from a cylindrical experimental implant placed in cancellous bone in an unloaded press-fit configuration, obtained after a four week observation period. Four sections each from seven textured titanium implants were investigated. The first resonant frequency and loss factor were measured. Average experimentally determined loss factor was 2% (SD 0.4%) and average first resonant frequency was 2.1 KHz (SD: 50). A 2D axisymmetric finite element (FE) model identified effective Young's modulus of tissue using experimental resonant frequencies as input. Average value was 42 MPa (SD: 2.4) and no significant difference between specimens was observed.In this pilot study, the non-destructive method allowed accurate measure of dynamic loss factor and resonant frequency and derivation of effective Young's modulus. Prior to implementing this dynamic protocol for broader mechanical evaluation of experimental implant fixation, further work is needed to determine if this affects results from subsequent destructive shear push-out tests.     

Improved measurement of brain deformation during mild head acceleration using a novel tagged MRI sequence

In vivo measurements of human brain deformation during mild acceleration are needed to help validate computational models of traumatic brain injury and to understand the factors that govern the mechanical response of the brain. Tagged magnetic resonance imaging is a powerful, noninvasive technique to track tissue motion in vivo which has been used to quantify brain deformation in live human subjects. However, these prior studies required from 72 to 144 head rotations to generate deformation data for a single image slice, precluding its use to investigate the entire brain in a single subject. Here, a novel method is introduced that significantly reduces temporal variability in the acquisition and improves the accuracy of displacement estimates. Optimization of the acquisition parameters in a gelatin phantom and three human subjects leads to a reduction in the number of rotations from 72 to 144 to as few as 8 for a single image slice. The ability to estimate accurate, well-resolved, fields of displacement and strain in far fewer repetitions will enable comprehensive studies of acceleration-induced deformation throughout the human brain in vivo.     

Systematic review and meta-analysis of the biomechanical properties of the human dura mater applicable in computational human head models

Accurate biomechanical properties of the human dura mater are required for computational models and to fabricate artificial substitutes for transplantation and surgical training purposes. Here, a systematic literature review was performed to summarize the biomechanical properties of the human dura mater that are reported in the literature. Furthermore, anthropometric data, information regarding the mechanically tested samples, and specifications with respect to the used mechanical testing setup were extracted. A meta-analysis was performed to obtain the pooled mean estimate for the elastic modulus, ultimate tensile strength, and strain at maximum force. A total of 17 studies were deemed eligible, which focused on human cranial and spinal dura mater in 13 and 4 cases, respectively. Pooled mean estimates for the elastic modulus (n?=?448), the ultimate tensile strength (n?=?448), and the strain at maximum force (n?=?431) of 68.1 MPa, 7.3 MPa and 14.4% were observed for native cranial dura mater. Gaps in the literature related to the extracted data were identified and future directions for mechanical characterizations of human dura mater were formulated. The main conclusion is that the most commonly used elastic modulus value of 31.5 MPa for the simulation of the human cranial dura mater in computational head models is likely an underestimation and an oversimplification given the morphological diversity of the tissue in different brain regions. Based on the here provided meta-analysis, a stiffer linear elastic modulus of 68 MPa was observed instead. However, further experimental data are essential to confirm its validity.

     

Mechanical properties of brain tissue in-vivo: experiment and computer simulation

Realistic computer simulation of neurosurgical procedures requires incorporation of the mechanical properties of brain tissue in the mathematical model. Possible applications of computer simulation of neurosurgery include non-rigid registration, virtual reality training and operation planning systems and robotic devices to perform minimally invasive brain surgery. A number of constitutive models of brain tissue, both single-phase and bi-phasic, have been proposed in recent years. The major deficiency of most of them, however, is the fact that they were identified using experimental data obtained in vitro and there is no certainty whether they can be applied in the realistic in vivo setting. In this paper we attempt to show that previously proposed by us hyper-viscoelastic constitutive model of brain tissue can be applied to simulating surgical procedures. An in vivo indentation experiment is described. The force-displacement curve for the loading speed typical for surgical procedures is concave upward containing no linear portion from which a meaningful elastic modulus might be determined. In order to properly analyse experimental data, a three-dimensional, non-linear finite element model of the brain was developed. Magnetic resonance imaging techniques were used to obtain geometric information needed for the model. The shape of the force-displacement curve obtained using the numerical solution was very similar to the experimental one. The predicted forces were about 31% lower than those recorded during the experiment. Having in mind that the coefficients in the model had been identified based on experimental data obtained in vitro, and large variability of mechanical properties of biological tissues, such agreement can be considered as very good. By appropriately increasing material parameters describing instantaneous stiffness of the tissue one is able, without changing the structure of the model, to reproduce experimental curve almost perfectly. Numerical studies showed also that the linear, viscoelastic model of brain tissue is not appropriate for the modelling brain tissue deformation even for moderate strains.     

Vascular elasticity measurement of the great saphenous vein based on optical coherence elastography

Vascular elasticity is important in physiological and clinical problems. The mechanical properties of the great saphenous vein (GSV) deserve attention. This research aims to measure the radial elasticity of ex vivo GSV using the optical coherence elasticity (OCE). The finite element model of the phantom is established, the displacement field is calculated, the radial mechanical characteristics of the simulation body are obtained. Furthermore, we performed OCE on seven isolated GSVs. The strain field is obtained by combining the relationship between strain and displacement to obtain the radial elastic modulus of GSVs. In the phantom experiment, the strain of the experimental region of interest is mainly between 0.1 and 0.4, while the simulation result is between 0.06 and 0.40. The radial elastic modulus of GSVs ranged from 3.83 kPa to 7.74 kPa. This study verifies the feasibility of the OCE method for measuring the radial elastic modulus of blood vessels.     

Rapid identification of elastic modulus of the interface tissue on dental implants surfaces using reduced-basis method and a neural network

This paper proposes a rapid inverse analysis approach based on the reduced-basis method (RBM) and neural network (NN) to identify the “unknown” elastic modulus (Young's modulus) of the interfacial tissue between a dental implant and the surrounding bones. In the present RBM–NN approach, a RBM model is first built to compute displacement responses of dental implant-bone structures subjected to a harmonic loading for a set of “assumed” Young's moduli. The RBM model is then used to train a NN model that is used for actual inverse analysis in real-time. Actual experimental measurements of displacement responses are fed into the trained NN model to inversely determine the “true” elastic modulus of the interfacial tissue. As an example, a physical model of dental implant-bone structure is built and inverse analysis is conducted to verify the present RBM–NN approach. Based on numerical simulation and actual experiments, it is confirmed that the identified results are very accurate, reliable, and the computational saving is very significant. The present RBM–NN approach is found robust and efficient for inverse material characterizations in noninvasive and/or nondestructive evaluations.     

Diffusing wave spectroscopy microrheology of actin filament networks.

Filamentous actin (F-actin), one of the constituents of the cytoskeleton, is believed to be the most important participant in the motion and mechanical integrity of eukaryotic cells. Traditionally, the viscoelastic moduli of F-actin networks have been measured by imposing a small mechanical strain and quantifying the resulting stress. The magnitude of the viscoelastic moduli, their concentration dependence and strain dependence, as well as the viscoelastic nature (solid-like or liquid-like) of networks of uncross-linked F-actin, have been the subjects of debate. Although this paper helps to resolve the debate and establishes the extent of the linear regime of F-actin networks' rheology, we report novel measurements of the high-frequency behavior of networks of F-actin, using a noninvasive light-scattering based technique, diffusing wave spectroscopy (DWS). Because no external strain is applied, our optical assay generates measurements of the mechanical properties of F-actin networks that avoid many ambiguities inherent in mechanical measurements. We observe that the elastic modulus has a small magnitude, no strain dependence, and a weak concentration dependence. Therefore, F-actin alone is not sufficient to generate the elastic modulus necessary to sustain the structural rigidity of most cells or support new cellular protrusions. Unlike previous studies, our measurements show that the mechanical properties of F-actin are highly dependent on the frequency content of the deformation. We show that the loss modulus unexpectedly dominates the elastic modulus at high frequencies, which are key for fast transitions. Finally, the measured mean square displacement of the optical probes, which is also generated by DWS measurements, offers new insight into the local bending fluctuations of the individual actin filaments and shows how they generate enhanced dissipation at short time scales.     

Stochastic multi-scale prediction on the apparent elastic moduli of trabecular bone considering uncertainties of biological apatite (BAp) crystallite orientation and image-based modelling

An assessment of the mechanical properties of trabecular bone is important in determining the fracture risk of human bones. Many uncertainty factors contribute to the dispersion of the estimated mechanical properties of trabecular bone. This study was undertaken in order to propose a computational scheme that will be able to predict the effective apparent elastic moduli of trabecular bone considering the uncertainties that are primarily caused by image-based modelling and trabecular stiffness orientation. The effect of image-based modelling which focused on the connectivity was also investigated. A stochastic multi-scale method using a first-order perturbation-based and asymptotic homogenisation theory was applied to formulate the stochastically apparent elastic properties of trabecular bone. The effective apparent elastic modulus was predicted with the introduction of a coefficient factor to represent the variation of bone characteristics due to inter-individual differences. The mean value of the predicted effective apparent Young's modulus in principal axis was found at approximately 460 MPa for respective 15.24% of bone volume fraction, and this is in good agreement with other experimental results. The proposed method may provide a reference for the reliable evaluation of the prediction of the apparent elastic properties of trabecular bone.