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.     

Finite element dynamic analysis of soft tissues using state-space model

A finite element (FE) model is employed to investigate the dynamic response of soft tissues under external excitations, particularly corresponding to the case of harmonic motion imaging. A solid 3D mixed ‘u–p’ element S8P0 is implemented to capture the near-incompressibility inherent in soft tissues. Two important aspects in structural modelling of these tissues are studied; these are the influence of viscous damping on the dynamic response and, following FE-modelling, a developed state-space formulation that valuates the efficiency of several order reduction methods. It is illustrated that the order of the mathematical model can be significantly reduced, while preserving the accuracy of the observed system dynamics. Thus, the reduced-order state-space representation of soft tissues for general dynamic analysis significantly reduces the computational cost and provides a unitary framework for the ‘forward’ simulation and ‘inverse’ estimation of soft tissues. Moreover, the results suggest that damping in soft-tissue is significant, effectively cancelling the contribution of all but the first few vibration modes.     

Effects of friction on the unconfined compressive response of articular cartilage: a finite element analysis

A finite element analysis is used to study a previously unresolved issue of the effects of platen-specimen friction on the response of the unconfined compression test; effects of platen permeability are also determined. The finite element formulation is based on the linear KLM biphasic model for articular cartilage and other hydrated soft tissues. A Galerkin weighted residual method is applied to both the solid phase and the fluid phase, and the continuity equation for the intrinsically incompressible binary mixture is introduced via a penalty method. The solid phase displacements and fluid phase velocities are interpolated for each element in terms of unknown nodal values, producing a system of first order differential equations which are solved using a standard numerical finite difference technique. An axisymmetric element of quadrilateral cross-section is developed and applied to the mechanical test problem of a cylindrical specimen of soft tissue in unconfined compression. These studies show that interfacial friction plays a major role in the unconfined compression response of articular cartilage specimens with small thickness to diameter ratios.     

A novel two-layer, coupled finite element approach for modeling the nonlinear elastic and viscoelastic behavior of human erythrocytes

A novel finite element approach is presented to simulate the mechanical behavior of human red blood cells (RBC, erythrocytes). As the RBC membrane comprises a phospholipid bilayer with an intervening protein network, we propose to model the membrane with two distinct layers. The fairly complex characteristics of the very thin lipid bilayer are represented by special incompressible solid shell elements and an anisotropic viscoelastic constitutive model. Properties of the protein network are modeled with an isotropic hyperelastic third-order material. The elastic behavior of the model is validated with existing optical tweezers studies with quasi-static deformations. Employing material parameters consistent with literature, simulation results are in excellent agreement with experimental data. Available models in literature neglect either the surface area conservation of the RBC membrane or realistic loading conditions of the optical tweezers experiments. The importance of these modeling assumptions, that are both included in this study, are discussed and their influence quantified. For the simulation of the dynamic motion of RBC, the model is extended to incorporate the cytoplasm. This is realized with a monolithic fully coupled fluid-structure interaction simulation, where the fluid is described by the incompressible Navier–Stokes equations in an arbitrary Lagrangian Eulerian framework. It is shown that both membrane viscosity and cytoplasm viscosity have significant influence on simulation results. Characteristic recovery times and energy dissipation for varying strain rates in dynamic laser trap experiments are calculated for the first time and are found to be comparable with experimental data.     

Finite element modeling of impulsive excitation and shear wave propagation in an incompressible,transversely isotropic medium

Elastic properties of materials can be measured by observing shear wave propagation following localized, impulsive excitations and relating the propagation velocity to a model of the material. However, characterization of anisotropic materials is difficult because of the number of elasticity constants in the material model and the complex dependence of propagation velocity relative to the excitation axis, material symmetries, and propagation directions. In this study, we develop a model of wave propagation following impulsive excitation in an incompressible, transversely isotropic (TI) material such as muscle. Wave motion is described in terms of three propagation modes identified by their polarization relative to the material symmetry axis and propagation direction. Phase velocities for these propagation modes are expressed in terms of five elasticity constants needed to describe a general TI material, and also in terms of three constants after the application of two constraints that hold in the limit of an incompressible material. Group propagation velocities are derived from the phase velocities to describe the propagation of wave packets away from the excitation region following localized excitation. The theoretical model is compared to the results of finite element (FE) simulations performed using a nearly incompressible material model with the five elasticity constants chosen to preserve the essential properties of the material in the incompressible limit. Propagation velocities calculated from the FE displacement data show complex structure that agrees quantitatively with the theoretical model and demonstrates the possibility of measuring all three elasticity constants needed to characterize an incompressible, TI material.     

Microplane constitutive model and computational framework for blood vessel tissue

This paper presents a nonlinearly elastic anisotropic microplane formulation in 3D for computational constitutive modeling of arterial soft tissue in the passive regime. The constitutive modeling of arterial (and other biological) soft tissue is crucial for accurate finite element calculations, which in turn are essential for design of implants, surgical procedures, bioartificial tissue, as well as determination of effect of progressive diseases on tissues and implants. The model presented is defined at a lower scale (mesoscale) than the conventional macroscale and it incorporates the effect of all the (collagen) fibers which are anisotropic structural components distributed in all directions within the tissue material in addition to that of isotropic bulk tissue. It is shown that the proposed model not only reproduces Holzapfel's recent model but also improves on it by accounting for the actual three-dimensional distribution of fiber orientation in the arterial wall, which endows the model with advanced capabilities in simulation of remodeling of soft tissue. The formulation is flexible so that its parameters could be adjusted to represent the arterial wall either as a single material or a material composed of several layers in finite element analyses of arteries. Explicit algorithms for both the material subroutine and the explicit integration with dynamic relaxation of equations of motion using finite element method are given. To circumvent the slow convergence of the standard dynamic relaxation and small time steps dictated by the stability of the explicit integrator, an adaptive dynamic relaxation technique that ensures stability and fastest possible convergence rates is developed. Incompressibility is enforced using penalty method with an updated penalty parameter. The model is used to simulate experimental data from the literature demonstrating that the model response is in excellent agreement with the data. An experimental procedure to determine the distribution of fiber directions in 3D for biological soft tissue is suggested in accordance with the microplane concept. It is also argued that this microplane formulation could be modified or extended to model many other phenomena of interest in biomechanics.     

Foot arch deformation and plantar fascia loading during running with rearfoot strike and forefoot strike: A dynamic finite element analysis

Forefoot strike becomes popular among runners because it facilitates better impact attenuation. However, forefoot strike may overload the plantar fascia and impose risk of plantar fasciitis. This study aimed to examine and compare the foot arch deformation and plantar fascia tension between different foot strike techniques in running using a computational modelling approach. A three-dimensional finite element foot model was reconstructed from the MRI of a healthy runner. The foot model included twenty bones, bulk soft tissue, ligaments, tendons, and plantar fascia. The time-series data of segmental kinematics, foot muscle force, and ankle joint reaction force were derived from a musculoskeletal model of the same participant based on the motion capture analysis and input as the boundary conditions for the finite element analysis. Rearfoot strike and forefoot strike running were simulated using a dynamic explicit solver. The results showed that, compared to rearfoot strike, forefoot strike reduced the foot arch height by 9.12% and increased the medial longitudinal arch angle by 2.06%. Forefoot strike also increased the plantar connective tissues stress by 18.28–200.11% and increased the plantar fascia tensile force by 18.71–109.10%. Although it is currently difficult to estimate the threshold value of stress or force that results in injury, forefoot strike runners appeared to be more vulnerable to plantar fasciitis.     

A large-strain finite element formulation for biological tissues with application to mitral valve leaflet tissue mechanics

This paper presents a finite element formulation suitable for large-strain modeling of biological tissues and uses this formulation to implement an accurate finite element model for mitral valve leaflet tissue. First, an experimentally derived strain energy function is obtained from literature. This function is implemented in finite elements using the mixed pressure-displacement formulation. A modification is made to aid in maintaining positive definiteness of the stiffness matrix at low strains. The numerical implementation is shown to be accurate in representing the analytical model of material behavior. The mixed formulation is useful for modeling of soft biological tissues in general, and the model presented here is applicable to finite element simulation of mitral valve mechanics.     

Analysis of the dynamic strains in a fingertip exposed to vibrations: Correlation to the mechanical stimuli on mechanoreceptors

The reduction in vibrotactile sensitivity in the fingertip is assumed to be associated with the exposure of the tissues to repetitive, non-physiological strains during dynamic loading. Experimental results demonstrated that the magnitude of a vibration-induced temporary threshold shift is dependent upon the vibration frequency of both the exposure and testing stimuli. In the present study, the frequency-dependent strain imposed on cutaneous and subcutaneous tissues of the fingertip is analyzed theoretically using a finite element model. The proposed fingertip model is two-dimensional and includes major anatomical substructures: skin, subcutaneous tissue, bone, and nail. The soft tissues (skin and subcutaneous tissues) were assumed to be nonlinearly elastic and viscoelastic, while the bone and nail were considered as linearly elastic. Simulations were performed for the contact between the fingertip and a flat surface for four different pre-compressions (0.5, 1.0, 1.5, and 2.0 mm). The frequency-dependent distributions of the dynamic strain magnitudes in the soft tissues were investigated. The model predictions indicated that the vibration exposure at a frequency range from 63 to 250 Hz will induce excessive dynamic strain in the deep zone of the finger tissues, effectively inhibiting the high-frequency mechanoreceptors; while the vibration exposure at low frequency (less than 31.5 Hz) tends to induce excessive dynamic strain in superficial layer in the tissues, inhibiting the low-frequency mechanoreceptors. These theoretical predictions are consistent with the experimental observations in literature. The proposed model can be used to predict the responses of the soft tissues in different depths to vibration exposures, providing valuable information and data that are essential for improving vibrotactile perception tests.     

OSAHS患者上气道和软腭的流固耦合有限元模型的建立

目的:重建OSAHS患者上气道和软腭的流固耦合有限元模型,研究OSAHS患者上气道及软腭气流动力学特征,为进一步探讨OSAHS的的发病机制奠定基础。方法:对一名中度OSAHS患者的上气道及周围组织进行MRI扫描,将以DICOM格式存储的扫描数据导入Mimics15.0软件中进行预处理,得到上气道和软腭的模型;再利用逆向工程软件Geomagic Studio 2013建立了2 mm气道壁;然后在3-D重建软件NX中,生成气道壁和气道以及软腭之间的组合模型;最后将该组合模型导入ANSYS Workbench13.0软件中,通过网格划分、定义材料属性、设定模型的边界条件操作建立了上气道和软腭的流固耦合有限元模型。结果:利用Mimics、Ansys等软件建立了完整的上气道和软腭的流固耦合有限元模型。共得到气道:2806835单元和529281个节点;气道壁:2304348单元和3487609个节点;软腭:131855单元和204784个节点。结论:本研究建立的上气道及软腭的流固耦合有限元模型符合人体的生物力学特点,为下一步的数值模拟实验提供了一个更真实、可靠的模型。     

Implementation of controlling strategy in a biomechanical lower limb model with active muscles for coupling multibody dynamics and finite element analysis

Computational biomechanics for human body modeling has generally been categorized into two separated domains: finite element analysis and multibody dynamics. Combining the advantages of both domains is necessary when tissue stress and physical body motion are both of interest. However, the method for this topic is still in exploration. The aim of this study is to implement unique controlling strategies in finite element model for simultaneously simulating musculoskeletal body dynamics and in vivo stress inside human tissues. A finite element lower limb model with 3D active muscles was selected for the implementation of controlling strategies, which was further validated against in-vivo human motion experiments. A unique feedback control strategy that couples together a basic Proportion-Integration-Differentiation (PID) controller and generic active signals from Computed Muscle Control (CMC) method of the musculoskeletal model or normalized EMG singles was proposed and applied in the present model. The results show that the new proposed controlling strategy show a good correlation with experimental test data of the normal gait considering joint kinematics, while stress distribution of local lower limb tissue can be also detected in real-time with lower limb motion. In summary, the present work is the first step for the application of active controlling strategy in the finite element model for concurrent simulation of both body dynamics and tissue stress. In the future, the present method can be further developed to apply it in various fields for human biomechanical analysis to monitor local stress and strain distribution by simultaneously simulating human locomotion.     

3D finite element modeling of pelvic organ prolapse

Objectives: The purpose of this study is to develop a validated 3D finite element model of the pelvic floor system which can offer insights into the mechanics of anterior vaginal wall prolapse and have the ability to assess biomedical device treatment methods. The finite element results should accurately mimic the clinical findings of prolapse due to intra-abdominal pressure (IAP) and soft tissues impairment conditions. Methods: A 3D model of pelvic system was created in Creo Parametric 2.0 based on MRI Images, which included uterus, cervix, vagina, cardinal ligaments, uterosacral ligaments, and a simplified levator plate and rectum. The geometrical model was imported into ANSYS Workbench 14.5. Mechanical properties of soft tissues were based on experimental data of tensile test results from current literature. Studies were conducted for IAP loadings on the vaginal wall and uterus, increasing from lowest to extreme values. Results: Anterior vaginal wall collapse occurred at an IAP value corresponding to maximal valsalva and showed similar collapsed shape as clinical findings. Prolapse conditions exhibited high sensitivity to vaginal wall stiffness, whereas healthy tissues was found to support the vagina against prolapse. Ligament impairment was found to have only a secondary effect on prolapse.     

Smooth muscle stiffness variation due to external longitudinal oscillations

This research focuses on an in vitro investigation of the stiffness changes of contracted airway smooth muscles (ASM) subjected to external longitudinal oscillations. ASM tissues were dissected from excised pig tracheas and stimulated by a chemical stimulus (acetylcholine, 10(-3) M) to produce maximum contractions. The tissues were then systematically excited with external oscillations. Various frequencies, amplitudes and durations were used in the experiments to determine stiffness changes in response to these variations. Force changes were recorded to reflect the muscle stiffness changes. Two stiffness definitions were used to quantify the results, dynamic stiffness to reflect variations during oscillation and static stiffness to reflect the net effect of oscillation. Under isometric contractions, these two stiffnesses were determined before, during and after oscillations. Incorporating an empirical stiffness equation, a two-dimensional finite element model (FEM) was developed to generalize the tissue responses to oscillation. The main outcomes from this work are: the dynamic stiffness has the tendency to decrease as the frequency and/or amplitude of external oscillation increases; the static stiffness has the tendency of decreasing with an increase in the frequency and/or amplitude of excitation until it reaches almost a constant value for frequencies at and above 25 Hz. The difference in the behavior of the dynamic and static stiffness changes may be attributed to the effect of elasticity and mass inertia that are involved in the dynamic motion. The findings of this research are in agreement with the hypothesis that oscillations exert a direct action on the contractile processes by causing an increased rate of actin-myosin detachments.     

An accurate finite element model of the cervical spine under quasi-static loading

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.