Technical Brief

Dynamic Simulation of Viscoelastic Soft Tissue in Acoustic Radiation Force Creep Imaging

[+] Author and Article Information
Xiaodong Zhao

Department of Mechanical and Aerospace Engineering,
Rutgers, the State University of New Jersey,
98 Brett Road,
Piscataway, NJ 08854-8058
e-mail: xiaodong.zhao@rutgers.edu

Assimina A. Pelegri

Fellow ASME
Department of Mechanical and Aerospace Engineering,
Rutgers, the State University of New Jersey,
98 Brett Road,
Piscataway, NJ 08854-8058
e-mail: pelegri@jove.rutgers.edu

1Corresponding author.

Manuscript received October 28, 2013; final manuscript received June 20, 2014; accepted manuscript posted July 1, 2014; published online July 15, 2014. Assoc. Editor: Jeffrey Ruberti.

J Biomech Eng 136(9), 094502 (Jul 15, 2014) (7 pages) Paper No: BIO-13-1510; doi: 10.1115/1.4027934 History: Received October 28, 2013; Revised June 20, 2014; Accepted July 01, 2014

Acoustic radiation force (ARF) creep imaging applies step ARF excitation to induce creep displacement of soft tissue, and the corresponding time-dependent responses are used to estimate soft tissue viscoelasticity or its contrast. Single degree of freedom (SDF) and homogeneous analytical models have been used to characterize soft tissue viscoelasticity in ARF creep imaging. The purpose of this study is to investigate the fundamental limitations of the commonly used SDF and homogeneous assumptions in ARF creep imaging. In this paper, finite element (FE) models are developed to simulate the dynamic behavior of viscoelastic soft tissue subjected to step ARF. Both homogeneous and heterogeneous models are studied with different soft tissue viscoelasticity and ARF configurations. The results indicate that the SDF model can provide good estimations for homogeneous soft tissue with high viscosity, but exhibits poor performance for low viscosity soft tissue. In addition, a smaller focal region of the ARF is desirable to reduce the estimation error with the SDF models. For heterogeneous media, the responses of the focal region are highly affected by the local heterogeneity, which results in deterioration of the effectiveness of the SDF and homogeneous simplifications.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Sugimoto, T., Ueha, S., and Itoh, K., 1990, “Tissue Hardness Measurement Using the Radiation Force of Focused Ultrasound,” Proceedings of the IEEE Ultrasonics Symposium, Honolulu, HI, Dec. 4–7, pp. 1377–1380.
Sarvazyan, A. P., Rudenko, O. V., Swanson, S. D., Fowlkes, J. B., and Emelianov, S. Y., 1998, “Shear Wave Elasticity Imaging: A New Ultrasonic Technology of Medical Diagnostics,” Ultrasound Med. Biol., 24(9), pp. 1419–1435. [CrossRef] [PubMed]
Bercoff, J., Tanter, M., and Fink, M., 2004, “Supersonic Shear Imaging: A New Technique for Soft Tissue Elasticity Mapping,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, 51(4), pp. 396–409. [CrossRef]
Nightingale, K., McAleavey, S., and Trahey, G., 2003, “Shear-Wave Generation Using Acoustic Radiation Force: In Vivo and Ex Vivo Results,” Ultrasound Med. Biol., 29(12), pp. 1715–1723. [CrossRef] [PubMed]
Chen, S., Urban, M., Pislaru, C., Kinnick, R., Zheng, Y., Yao, A., and Greenleaf, J., 2009, “Shearwave Dispersion Ultrasound Vibrometry (SDUV) for Measuring Tissue Elasticity and Viscosity,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, 56(1), pp. 55–62. [CrossRef]
Vappou, J., Maleke, C., and Konofagou, E. E., 2009, “Quantitative Viscoelastic Parameters Measured by Harmonic Motion Imaging,” Phys. Med. Biol., 54(11), pp. 3579–3594. [CrossRef] [PubMed]
Walker, W. F., Fernandez, F. J., and Negron, L. A., 2000, “A Method of Imaging Viscoelastic Parameters With Acoustic Radiation Force,” Phys. Med. Biol., 45(6), pp. 1437–1447. [CrossRef] [PubMed]
Viola, F., and Walker, W. F., 2003, “Radiation Force Imaging of Viscoelastic Properties With Reduced Artifacts,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, 50(6), pp. 736–742. [CrossRef]
Mauldin, F. W., Haider, M. A., Loboa, E. G., Behler, R. H., Euliss, L. E., Pfeiler, T. W., and Gallippi, C. M., 2008, “Monitored Steady-State Excitation and Recovery (MSSR) Radiation Force Imaging Using Viscoelastic Models,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, 55(7), pp. 1597–1610. [CrossRef]
Amador, C., Urban, M. W., Chen, S., and Greenleaf, J. F., 2012, “Loss Tangent and Complex Modulus Estimated by Acoustic Radiation Force Creep and Shear Wave Dispersion,” Phys. Med. Biol., 57(5), pp. 1263–1282. [CrossRef] [PubMed]
Scola, M. R., Baggesen, L. M., and Gallippi, C. M., 2012, “Multi-Push (MP) Acoustic Radiation Force (ARF) Ultrasound for Assessing Tissue Viscoelasticity, In Vivo,” 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), San Diego, CA, Aug. 28–Sept. 1, pp. 2323–2326.
Palmeri, M. L., McAleavey, S. A., Fong, K. L., Trahey, G. E., and Nightingale, K. R., 2006, “Dynamic Mechanical Response of Elastic Spherical Inclusions to Impulsive Acoustic Radiation Force Excitation,” IEEE Trans. Ultrason., Ferroelectr., Freq. Control, 53(11), pp. 2065–2079. [CrossRef]
Zhao, X., and Pelegri, A. A., 2014, “Modelling of Global Boundary Effects on Harmonic Motion Imaging of Soft Tissues,” Comput. Methods Biomech. Biomed. Eng., 17(9), pp. 1021–1031. [CrossRef]
Zhang, Y., Hall, L. O., Goldgof, D. B., and Sarkar, S., 2006, “A Constrained Genetic Approach for Computing Material Property of Elastic Objects,” IEEE Trans. Evol. Comput., 10(3), pp. 341–357. [CrossRef]
Catheline, S., Gennisson, J. L., Delon, G., Fink, M., Sinkus, R., Abouelkaram, S., and Culioli, J., 2004, “Measurement of Viscoelastic Properties of Homogeneous Soft Solid Using Transient Elastography: An Inverse Problem Approach,” J. Acoust. Soc. Am., 116(6), pp. 3734–3741. [CrossRef] [PubMed]
Skovoroda, A. R., and Sarvazyan, A. P., 1999, “Determination of Viscoelastic Shear Characteristics of a Medium From Its Response to Focused Ultrasonic Loading,” Biophysics, 44, pp. 325–329.
Nightingale, K. R., Rouze, N. C., Wang, M. H., Zhai, L., and Palmeri, M. L., 2011, “Comparison of Qualitative and Quantitative Acoustic Radiation Force Based Elasticity Imaging Methods,” 2011 IEEE International Symposium on IEEE, Biomedical Imaging: From Nano to Macro, Chicago, IL, Mar. 30–Apr. 2, pp. 1606–1609.
Bower, A. F., 2010, Applied Mechanics of Solids, CRC Press, Boca Raton, FL.
ABAQUS 6.11 Documentation, 2011, Providence, RI, Dassault Systèmes Simulia Corp.
Sarvazyan, A. P., Skovoroda, A. R., Emelianov, S. Y., Fowlkes, J. B., Pipe, J. G., Adler, R. S., Buxton, R. B., and Carson, P. L., 1995, “Biophysical Bases of Elasticity Imaging,” Acoustical Imaging, Springer, New York, pp. 223–240.
Wells, P. N., and Liang, H. D., 2011, “Medical Ultrasound: Imaging of Soft Tissue Strain and Elasticity,” J. R. Soc., Interface, 8(64), pp. 1521–1549. [CrossRef]
Maleke, C., Luo, J., Gamarnik, V., Lu, X. L., and Konofagou, E. E., 2010, “Simulation Study of Amplitude-Modulated (AM) Harmonic Motion Imaging (HMI) for Stiffness Contrast Quantification With Experimental Validation,” Ultrason. Imaging, 32(3), pp. 154–176. [CrossRef] [PubMed]
Nightingale, K., Palmeri, M., and Trahey, G., 2006, “Analysis of Contrast in Images Generated With Transient Acoustic Radiation Force,” Ultrasound Med. Biol., 32(1), pp. 61–72. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Model diagram and validation: (a) Diagram of the axisymmetric model with a spherical inclusion in the center of the model and (b) axial displacement induced by ARF for the homogeneous case with μ = 3 kPa. The horizontal axis in (b) is the axial distance from the focal center.

Grahic Jump Location
Fig. 2

Creep displacement responses to step ARF for soft tissue with different viscoelasticity. Three time constants are studied: (a) τ = 0.0003 s, (b) τ = 0.0009 s, and (c) τ = 0.0027 s. The corresponding normalized creep displacement responses are shown in (d), (e), and (f).

Grahic Jump Location
Fig. 3

REE of τ by fitting the SDF models with the FE simulated creep displacement responses. (a) SDF model without considering the inertial effect and (b) SDF model with the inertial component included. Each marker represents cases with the same shear viscosity.

Grahic Jump Location
Fig. 4

Creep displacement responses to different ARF configurations: (a) ARFs with different magnitude, but the same distribution, i.e., Fnumber = 0.83; and (b) ARFs with different distribution, but the same magnitude, i.e., fo. The corresponding normalized creep displacement responses are shown in (c) and (d). The model is homogeneous with μ = 3 kPa, and τ = 0.0009 s.

Grahic Jump Location
Fig. 5

Creep displacement responses at the origin of the heterogeneous models with different inclusion sizes: (a) Sphere with diameter 3 mm and (b) Sphere with diameter 6 mm. The corresponding normalized creep displacement responses are shown in (c) and (d). The solid black line denotes the homogeneous case. τ = 0.0003 s is for both background and inclusion.

Grahic Jump Location
Fig. 6

Axial normal strain field near the region of the 3-mm-diameter spherical inclusion after a 10 ms step ARF excitation. The dimension of the region is h = 6 mm (axial length) and r = 3 mm (radial length). The time constant is 0.0009 s and the shear moduli are: (a) μΒ = 3 kPa, and μΙ = 3 kPa; (b) μΒ = 0.3 kPa, and μΙ = 3 kPa; and (c) μΒ = 30 kPa, and μΙ = 3 kPa.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In