Search

Dynamic mechanical analysis (DMA) is a common way to measure the mechanical properties of materials as functions of frequency. Traditionally, a viscoelastic mechanical model is applied and current DMA techniques fit an analytical approximation to measured dynamic motion data by neglecting inertial forces and adding empirical correction factors to account for transverse boundary displacements. Here, a finite-element (FE) approach to processing DMA data was developed to estimate poroelastic material properties. Frequency-dependent inertial forces, which are significant in soft media and often neglected in DMA, were included in the FE model. The technique applies a constitutive relation to the DMA measurements and exploits a nonlinear inversion to estimate the material properties in the model that best fit the model response to the DMA data. A viscoelastic version of this approach was developed to validate the approach by comparing complex modulus estimates to the direct DMA results. Both analytical and FE poroelastic models were also developed to explore their behavior in the DMA testing environment. All of the models were applied to tofu as a representative soft poroelastic material that is a common phantom in elastography imaging studies. Five samples of three different stiffnesses were tested from 1-14 Hz with rough platens placed on the top and bottom surfaces of the material specimen under test to restrict transverse displacements and promote fluid-solid interaction. The viscoelastic models were identical in the static case, and nearly the same at frequency with inertial forces accounting for some of the discrepancy. The poroelastic analytical method was not sufficient when the relevant physical boundary constraints were applied, whereas the poroelastic FE approach produced high quality estimates of shear modulus and hydraulic conductivity. These results illustrated appropriate shear modulus contrast between tofu samples and yielded a consistent contrast in hydraulic conductivity as well.

Purpose: Descriptions of the structure of brain tissue as a porous cellular matrix support application of a poroelastic (PE) mechanical model which includes both solid and fluid phases. However, the majority of brain magnetic resonance elastography (MRE) studies use a single phase viscoelastic (VE) model to describe brain tissue behavior, in part due to availability of relatively simple direct inversion strategies for mechanical property estimation. A notable exception is low frequency intrinsic actuation MRE, where PE mechanical properties are imaged with a nonlinear inversion algorithm., Methods: This paper investigates the effect of model choice at each end of the spectrum of in vivo human brain actuation frequencies. Repeat MRE examinations of the brains of healthy volunteers were used to compare image quality and repeatability for each inversion model for both 50 Hz externally produced motion and ≈1 Hz intrinsic motions. Additionally, realistic simulated MRE data were generated with both VE and PE finite element solvers to investigate the effect of inappropriate model choice for ideal VE and PE materials., Results: In vivo, MRE data revealed that VE inversions appear more representative of anatomical structure and quantitatively repeatable for 50 Hz induced motions, whereas PE inversion produces better results at 1 Hz. Reasonable VE approximations of PE materials can be derived by equating the equivalent wave velocities for the two models, provided that the timescale of fluid equilibration is not similar to the period of actuation. An approximation of the equilibration time for human brain reveals that this condition is violated at 1 Hz but not at 50 Hz. Additionally, simulation experiments when using the "wrong" model for the inversion demonstrated reasonable shear modulus reconstructions at 50 Hz, whereas cross-model inversions at 1 Hz were poor quality. Attenuation parameters were sensitive to changes in the forward model at both frequencies, however, no spatial information was recovered because the mechanisms of VE and PE attenuation are different., Conclusions: VE inversions are simpler with fewer unknown properties and may be sufficient to capture the mechanical behavior of PE brain tissue at higher actuation frequencies. However, accurate modeling of the fluid phase is required to produce useful mechanical property images at the lower frequencies of intrinsic brain motions.

Purpose: Recently, the attenuating behavior of soft tissue has been addressed in magnetic resonance elastography by the inclusion of a damping mechanism in the methods used to reconstruct the resulting mechanical property image. To date, this mechanism has been based on a viscoelastic model for material behavior. Rayleigh, or proportional, damping provides a more generalized model for elastic energy attenuation that uses two parameters to characterize contributions proportional to elastic and inertial forces. In the case of time-harmonic vibration, these two parameters lead to both the elastic modulus and the density being complex valued (as opposed to the case of pure viscoelasticity, where only the elastic modulus is complex valued)., Methods: This article presents a description of Rayleigh damping in the time-harmonic case, discussing the differences between this model and the viscoelastic damping models. In addition, the results from a subzone based Rayleigh damped elastography study of gelatin and tofu phantoms are discussed, along with preliminary results from in vivo breast data., Results: Both the phantom and the tissue studies presented here indicate a change in the Rayleigh damping structure, described as Rayleigh composition, between different material types, with tofu and healthy tissue showing lower Rayleigh composition values than gelatin or cancerous tissue., Conclusions: It is possible that Rayleigh damping elastography and the concomitant Rayleigh composition images provide a mechanism for differentiating tissue structure in addition to measuring elastic stiffness and attenuation. Such information could be valuable in the use of Rayleigh damped magnetic resonance elastography as a diagnostic imaging tool.

Elastography is an emerging imaging technique that focuses on assessing the resistance to deformation of soft biological tissues in vivo. Magnetic resonance elastography (MRE) uses measured displacement fields resulting from low-amplitude, low-frequency (10 Hz-1 kHz) time-harmonic vibration to recover images of the elastic property distribution of tissues including breast, liver, muscle, prostate, and brain. While many soft tissues display complex time-dependent behavior not described by linear elasticity, the models most commonly employed in MRE parameter reconstructions are based on elastic assumptions. Further, elasticity models fail to include the interstitial fluid phase present in vivo. Alternative continuum models, such as consolidation theory, are able to represent tissue and other materials comprising two distinct phases, generally consisting of a porous elastic solid and penetrating fluid. MRE reconstructions of simulated elastic and poroelastic phantoms were performed to investigate the limitations of current-elasticity-based methods in producing accurate elastic parameter estimates in poroelastic media. The results indicate that linearly elastic reconstructions of fluid-saturated porous media at amplitudes and frequencies relevant to steady-state MRE can yield misleading effective property distributions resulting from the complex interaction between their solid and fluid phases.

A reconstruction process featuring full parameterization of the three dimensional, time-harmonic equations of linear elasticity is developed and reconstructed property images are presented from simulation-based investigation. While interesting in its own right through the potential for increased adaptability of these reconstructive elastic imaging techniques, this study also presents a set of analysis tools used to study the poor convergence behavior found in the case of tissue like conditions (i.e. nearly incompressible materials). The choice of elastic properties for imaging in elastography research remains an open question at this point; the use of the stability and sensitivity-based analytical methods described here will help to predict and understand the value and reliability of different parameterizations of elasticity imaging. Additionally, though results indicate significant work needs to be done to achieve effective multiparameter reconstructive imaging, the methods detailed here offer the promise of increased flexibility and sophistication in elastographic imaging techniques.