Dynamic simulation of viscoelastic soft tissues in harmonic motion imaging application

A finite element model was built to simulate the dynamic behavior of soft tissues subjected to sinusoidal excitation during harmonic motion imaging. In this study, soft tissues and tissue-like phantoms were modeled as isotropic, viscoelastic, and nearly incompressible media. A 3D incompressible mixed u–p element of eight nodes, S1P0, was developed to accurately calculate the stiffness matrix for soft tissues. The finite element equations of motion were solved using the Newmark method. The Voigt description for tissue viscosity was applied to estimate the relative viscous coefficient from the phase shift between the response and excitation in a harmonic case. After validating our model via ANSYS simulation and experiments, a MATLAB finite element program was then employed to explore the effect of excitation location, viscosity, and multiple frequencies on the dynamic displacement at the frequency of interest.