ABSTRACT
Computational biomechanics plays an important role in biomedical engineering: using modeling to understand pathophysiology, treatment and device design. While experimental evidence indicates that the mechanical response of most tissues is viscoelastic, current biomechanical models in the computational community often assume hyperelastic material models. Fractional viscoelastic constitutive models have been successfully used in literature to capture viscoelastic material response; however, the translation of these models into computational platforms remains limited. Many experimentally derived viscoelastic constitutive models are not suitable for three-dimensional simulations. Furthermore, the use of fractional derivatives can be computationally prohibitive, with a number of current numerical approximations having a computational cost that is ðª ( N T 2 ) and a storage cost that is ðª(NT ) (NT denotes the number of time steps). In this paper, we present a novel numerical approximation to the Caputo derivative which exploits a recurrence relation similar to those used to discretize classic temporal derivatives, giving a computational cost that is ðª(NT ) and a storage cost that is fixed over time. The approximation is optimized for numerical applications, and an error estimate is presented to demonstrate the efficacy of the method. The method, integrated into a finite element solid mechanics framework, is shown to be unconditionally stable in the linear viscoelastic case. It was then integrated into a computational biomechanical framework, with several numerical examples verifying the accuracy and computational efficiency of the method, including in an analytic test, in an analytic fractional differential equation, as well as in a computational biomechanical model problem.
ABSTRACT
Solid tumour growth is often associated with the accumulation of mechanical stresses acting on the surrounding host tissue. Due to tissue nonlinearity, the shear modulus of the peri-tumoural region inherits a signature from the tumour expansion which depends on multiple factors, including the soft tissue constitutive behaviour and its stress/strain state. Shear waves used in MR-elastography (MRE) sense the apparent change in shear modulus along their propagation direction, thereby probing the anisotropic stiffness field around the tumour. We developed an analytical framework for a heterogeneous shear modulus distribution using a thick-shelled sphere approximation of the tumour and soft tissue ensemble. A hyperelastic material (plastisol) was identified to validate the proposed theory in a phantom setting. A balloon-catheter connected to a pressure sensor was used to replicate the stress generated from tumour pressure and growth while MRE data were acquired. The shear modulus anisotropy retrieved from the reconstructed elastography data confirmed the analytically predicted patterns at various levels of inflation. An alternative measure, combining the generated deformation and the local wave direction and independent of the reconstruction strategy, was also proposed to correlate the analytical findings with the stretch probed by the waves. Overall, this work demonstrates that MRE in combination with non-linear mechanics, is able to identify the apparent shear modulus variation arising from the strain generated by a growth within tissue, such as an idealised model of tumour. Investigation in real tissue represents the next step to further investigate the implications of endogenous forces in tissue characterisation through MRE.
Subject(s)
Elasticity Imaging Techniques/methods , Magnetic Resonance Imaging/methods , Neoplasms/diagnosis , Anisotropy , Biomimetic Materials , Elasticity Imaging Techniques/instrumentation , Humans , Magnetic Resonance Imaging/instrumentation , Phantoms, Imaging , Stress, MechanicalABSTRACT
Understanding the biomechanics of the heart in health and disease plays an important role in the diagnosis and treatment of heart failure. The use of computational biomechanical models for therapy assessment is paving the way for personalized treatment, and relies on accurate constitutive equations mapping strain to stress. Current state-of-the art constitutive equations account for the nonlinear anisotropic stress-strain response of cardiac muscle using hyperelasticity theory. While providing a solid foundation for understanding the biomechanics of heart tissue, most current laws neglect viscoelastic phenomena observed experimentally. Utilizing experimental data from human myocardium and knowledge of the hierarchical structure of heart muscle, we present a fractional nonlinear anisotropic viscoelastic constitutive model. The model is shown to replicate biaxial stretch, triaxial cyclic shear and triaxial stress relaxation experiments (mean error â¼7.68%), showing improvements compared to its hyperelastic (mean error â¼24%) counterparts. Model sensitivity, fidelity and parameter uniqueness are demonstrated. The model is also compared to rate-dependent biaxial stretch as well as different modes of biaxial stretch, illustrating extensibility of the model to a range of loading phenomena. STATEMENT OF SIGNIFICANCE: The viscoelastic response of human heart tissues has yet to be integrated into common constitutive models describing cardiac mechanics. In this work, a fractional viscoelastic modeling approach is introduced based on the hierarchical structure of heart tissue. From these foundations, the current state-of-the-art biomechanical models of the heart muscle are transformed using fractional viscoelasticity, replicating passive muscle function across multiple experimental tests. Comparisons are drawn with current models to highlight the improvements of this approach and predictive responses show strong qualitative agreement with experimental data. The proposed model presents the first constitutive model aimed at capturing viscoelastic nonlinear response across multiple testing regimes, providing a platform for better understanding the biomechanics of myocardial tissue in health and disease.
Subject(s)
Models, Biological , Myocardium , Anisotropy , Biomechanical Phenomena , Elasticity , Humans , Stress, Mechanical , ViscosityABSTRACT
Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) allow investigations in a human cardiac model system, but disorganized mechanics and immaturity of hPSC-CMs on standard two-dimensional surfaces have been hurdles. Here, we developed a platform of micron-scale cardiac muscle bundles to control biomechanics in arrays of thousands of purified, independently contracting cardiac muscle strips on two-dimensional elastomer substrates with far greater throughput than single cell methods. By defining geometry and workload in this reductionist platform, we show that myofibrillar alignment and auxotonic contractions at physiologic workload drive maturation of contractile function, calcium handling, and electrophysiology. Using transcriptomics, reporter hPSC-CMs, and quantitative immunofluorescence, these cardiac muscle bundles can be used to parse orthogonal cues in early development, including contractile force, calcium load, and metabolic signals. Additionally, the resultant organized biomechanics facilitates automated extraction of contractile kinetics from brightfield microscopy imaging, increasing the accessibility, reproducibility, and throughput of pharmacologic testing and cardiomyopathy disease modeling.
Subject(s)
Heart/growth & development , Myocardium , Myocytes, Cardiac/cytology , Pluripotent Stem Cells/cytology , Biomechanical Phenomena , Calcium/metabolism , Cell Culture Techniques , Dimethylpolysiloxanes , Electrophysiological Phenomena , Gene Expression Profiling , High-Throughput Screening Assays/instrumentation , Humans , Lab-On-A-Chip Devices , Models, Cardiovascular , Myocardial Contraction , Myocardium/cytology , Myocardium/metabolism , Myofibrils/metabolism , Reproducibility of ResultsABSTRACT
Soft tissue mechanical characterisation is important in many areas of medical research. Examples span from surgery training, device design and testing, sudden injury and disease diagnosis. The liver is of particular interest, as it is the most commonly injured organ in frontal and side motor vehicle crashes, and also assessed for inflammation and fibrosis in chronic liver diseases. Hence, an extensive rheological characterisation of liver tissue would contribute to advancements in these areas, which are dependent upon underlying biomechanical models. The aim of this paper is to define a liver constitutive equation that is able to characterise the nonlinear viscoelastic behaviour of liver tissue under a range of deformations and frequencies. The tissue response to large amplitude oscillatory shear (1-50%) under varying preloads (1-20%) and frequencies (0.5-2 Hz) is modelled using viscoelastic-adapted forms of the Mooney-Rivlin, Ogden and exponential models. These models are fit to the data using classical or modified objective norms. The results show that all three models are suitable for capturing the initial nonlinear regime, with the latter two being capable of capturing, simultaneously, the whole deformation range tested. The work presented here provides a comprehensive analysis across several material models and norms, leading to an identifiable constitutive equation that describes the nonlinear viscoelastic behaviour of the liver.
Subject(s)
Elasticity , Liver/physiology , Nonlinear Dynamics , Animals , Biomechanical Phenomena , Cattle , Models, Biological , Viscosity , Weight-BearingABSTRACT
The solid and fluid pressures of tumours are often elevated relative to surrounding tissue. This increased pressure is known to correlate with decreased treatment efficacy and potentially with tumour aggressiveness and therefore, accurate noninvasive estimates of tumour pressure would be of great value. We present a proof-of-concept method to infer the total tumour pressure, that is the sum of the fluid and solid parts, by examining stiffness in the peritumoural tissue with MR elastography and utilising nonlinear biomechanical models. The pressure from the tumour deforms the surrounding tissue leading to changes in stiffness. Understanding and accounting for these biases in stiffness has the potential to enable estimation of total tumour pressure. Simulations are used to validate the method with varying pressure levels, tumour shape, tumour size, and noise levels. Results show excellent matching in low noise cases and still correlate well with higher noise. Percent error remains near or below 10% for higher pressures in all noise level cases. Reconstructed pressures were also calculated from experiments with a catheter balloon embedded in a plastisol phantom at multiple inflation levels. Here the reconstructed pressures generally match the increases in pressure measured during the experiments. Percent errors between average reconstructed and measured pressures at four inflation states are 17.9%, 52%, 23.2%, and 0.9%. Future work will apply this method to in vivo data, potentially providing an important biomarker for cancer diagnosis and treatment.
Subject(s)
Elasticity Imaging Techniques/methods , Neoplasms/pathology , Biomechanical Phenomena , Computer Simulation , Humans , Models, Theoretical , Phantoms, Imaging , PressureABSTRACT
Characterisation of soft tissue mechanical properties is a topic of increasing interest in translational and clinical research. Magnetic resonance elastography (MRE) has been used in this context to assess the mechanical properties of tissues in vivo noninvasively. Typically, these analyses rely on linear viscoelastic wave equations to assess material properties from measured wave dynamics. However, deformations that occur in some tissues (e.g. liver during respiration, heart during the cardiac cycle, or external compression during a breast exam) can yield loading bias, complicating the interpretation of tissue stiffness from MRE measurements. In this paper, it is shown how combined knowledge of a material's rheology and loading state can be used to eliminate loading bias and enable interpretation of intrinsic (unloaded) stiffness properties. Equations are derived utilising perturbation theory and Cauchy's equations of motion to demonstrate the impact of loading state on periodic steady-state wave behaviour in nonlinear viscoelastic materials. These equations demonstrate how loading bias yields apparent material stiffening, softening and anisotropy. MRE sensitivity to deformation is demonstrated in an experimental phantom, showing a loading bias of up to twofold. From an unbiased stiffness of [Formula: see text] Pa in unloaded state, the biased stiffness increases to 9767.5 [Formula: see text]1949.9 Pa under a load of [Formula: see text] 34% uniaxial compression. Integrating knowledge of phantom loading and rheology into a novel MRE reconstruction, it is shown that it is possible to characterise intrinsic material characteristics, eliminating the loading bias from MRE data. The framework introduced and demonstrated in phantoms illustrates a pathway that can be translated and applied to MRE in complex deforming tissues. This would contribute to a better assessment of material properties in soft tissues employing elastography.