1
|
Brown AL, Sexton ZA, Hu Z, Yang W, Marsden AL. Computational approaches for mechanobiology in cardiovascular development and diseases. Curr Top Dev Biol 2024; 156:19-50. [PMID: 38556423 DOI: 10.1016/bs.ctdb.2024.01.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/02/2024]
Abstract
The cardiovascular development in vertebrates evolves in response to genetic and mechanical cues. The dynamic interplay among mechanics, cell biology, and anatomy continually shapes the hydraulic networks, characterized by complex, non-linear changes in anatomical structure and blood flow dynamics. To better understand this interplay, a diverse set of molecular and computational tools has been used to comprehensively study cardiovascular mechanobiology. With the continual advancement of computational capacity and numerical techniques, cardiovascular simulation is increasingly vital in both basic science research for understanding developmental mechanisms and disease etiologies, as well as in clinical studies aimed at enhancing treatment outcomes. This review provides an overview of computational cardiovascular modeling. Beginning with the fundamental concepts of computational cardiovascular modeling, it navigates through the applications of computational modeling in investigating mechanobiology during cardiac development. Second, the article illustrates the utility of computational hemodynamic modeling in the context of treatment planning for congenital heart diseases. It then delves into the predictive potential of computational models for elucidating tissue growth and remodeling processes. In closing, we outline prevailing challenges and future prospects, underscoring the transformative impact of computational cardiovascular modeling in reshaping cardiovascular science and clinical practice.
Collapse
Affiliation(s)
- Aaron L Brown
- Department of Mechanical Engineering, Stanford University, Stanford, CA, United States
| | - Zachary A Sexton
- Department of Bioengineering, Stanford University, Stanford, CA, United States
| | - Zinan Hu
- Department of Mechanical Engineering, Stanford University, Stanford, CA, United States
| | - Weiguang Yang
- Department of Pediatrics, Stanford University, Stanford, CA, United States
| | - Alison L Marsden
- Department of Bioengineering, Stanford University, Stanford, CA, United States; Department of Pediatrics, Stanford University, Stanford, CA, United States.
| |
Collapse
|
2
|
Schwarz EL, Pfaller MR, Szafron JM, Latorre M, Lindsey SE, Breuer CK, Humphrey JD, Marsden AL. A Fluid-Solid-Growth Solver for Cardiovascular Modeling. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2023; 417:116312. [PMID: 38044957 PMCID: PMC10691594 DOI: 10.1016/j.cma.2023.116312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2023]
Abstract
We implement full, three-dimensional constrained mixture theory for vascular growth and remodeling into a finite element fluid-structure interaction (FSI) solver. The resulting "fluid-solid-growth" (FSG) solver allows long term, patient-specific predictions of changing hemodynamics, vessel wall morphology, tissue composition, and material properties. This extension from short term (FSI) to long term (FSG) simulations increases clinical relevance by enabling mechanobioloigcally-dependent studies of disease progression in complex domains.
Collapse
Affiliation(s)
- Erica L Schwarz
- Department of Bioengineering, Stanford Univeristy, Stanford, CA 94306, USA
| | - Martin R Pfaller
- Department of Pediatrics - Cardiology, Stanford Univeristy, Stanford, CA 94306, USA
| | - Jason M Szafron
- Department of Pediatrics - Cardiology, Stanford Univeristy, Stanford, CA 94306, USA
| | - Marcos Latorre
- Center for Research and Innovation in Bioengineering, Universitat Politècnica de València, València 46022, Spain
| | - Stephanie E Lindsey
- Department of Pediatrics - Cardiology, Stanford Univeristy, Stanford, CA 94306, USA
| | - Christopher K Breuer
- Department of Surgery, Nationwide Children's Hospital, Columbus, OH 43210, USA
- Center for Regenerative Medicine, Abigail Wexner Research Institute, Nationwide Children's Hospital, Columbus, OH 43215, USA
| | - Jay D Humphrey
- Department of Biomedical Engineering, Yale Univeristy, New Haven, CT 06520, USA
| | - Alison L Marsden
- Department of Bioengineering, Stanford Univeristy, Stanford, CA 94306, USA
- Department of Pediatrics - Cardiology, Stanford Univeristy, Stanford, CA 94306, USA
| |
Collapse
|
3
|
Cheng F, Watton PN, Pederzani G, Kurobe M, Takaoka EI, Chapple C, Birder L, Yoshimura N, Robertson AM. A constrained mixture-micturition-growth (CMMG) model of the urinary bladder: Application to partial bladder outlet obstruction (BOO). J Mech Behav Biomed Mater 2022; 134:105337. [PMID: 35863296 PMCID: PMC9835014 DOI: 10.1016/j.jmbbm.2022.105337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 04/13/2022] [Accepted: 06/24/2022] [Indexed: 01/14/2023]
Abstract
We present a constrained mixture-micturition-growth (CMMG) model for the bladder. It simulates bladder mechanics, voiding function (micturition) and tissue adaptations in response to altered biomechanical conditions. The CMMG model is calibrated with both in vivo and in vitro data from healthy male rat urinary bladders (cystometry, bioimaging of wall structure, mechanical testing) and applied to simulate the growth and remodeling (G&R) response to partial bladder outlet obstruction (BOO). The bladder wall is represented as a multi-layered, anisotropic, nonlinear constrained mixture. A short time scale micturition component of the CMMG model accounts for the active and passive mechanics of voiding. Over a second, longer time scale, G&R algorithms for the evolution of both cellular and extracellular constituents act to maintain/restore bladder (homeostatic) functionality. The CMMG model is applied to a spherical membrane model of the BOO bladder utilizing temporal data from an experimental male rodent model to parameterize and then verify the model. Consistent with the experimental studies of BOO, the model predicts: an initial loss of voiding capacity followed by hypertrophy of SMC to restore voiding function; bladder enlargement; collagen remodeling to maintain its role as a protective sheath; and increased voiding duration with lower average flow rate. This CMMG model enables a mechanistic approach for investigating the bladder's structure-function relationship and its adaption in pathological conditions. While the approach is illustrated with a conceptual spherical bladder model, it provides the basis for application of the CMMG model to anatomical geometries. Such a mechanistic approach has promise as an in silico tool for the rational development of new surgical and pharmacological treatments for bladder diseases such as BOO.
Collapse
Affiliation(s)
- Fangzhou Cheng
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, United States
| | - Paul N Watton
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, United States; Department of Computer Science & Insigneo Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom.
| | - Giulia Pederzani
- Department of Computer Science & Insigneo Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom
| | - Masahiro Kurobe
- Department of Urology, University of Pittsburgh, Pittsburgh, United States
| | - Ei-Ichiro Takaoka
- Department of Urology, University of Pittsburgh, Pittsburgh, United States
| | - Chris Chapple
- Sheffield Teaching Hospitals NHS Foundation Trust, Sheffield, United Kingdom
| | - Lori Birder
- Sheffield Teaching Hospitals NHS Foundation Trust, Sheffield, United Kingdom; Department of Medicine, University of Pittsburgh, United States
| | - Naoki Yoshimura
- Department of Urology, University of Pittsburgh, Pittsburgh, United States
| | - Anne M Robertson
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, United States
| |
Collapse
|
4
|
Goswami S, Li DS, Rego BV, Latorre M, Humphrey JD, Karniadakis GE. Neural operator learning of heterogeneous mechanobiological insults contributing to aortic aneurysms. J R Soc Interface 2022; 19:20220410. [PMID: 36043289 PMCID: PMC9428523 DOI: 10.1098/rsif.2022.0410] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2022] [Accepted: 08/05/2022] [Indexed: 11/12/2022] Open
Abstract
Thoracic aortic aneurysm (TAA) is a localized dilatation of the aorta that can lead to life-threatening dissection or rupture. In vivo assessments of TAA progression are largely limited to measurements of aneurysm size and growth rate. There is promise, however, that computational modelling of the evolving biomechanics of the aorta could predict future geometry and properties from initiating mechanobiological insults. We present an integrated framework to train a deep operator network (DeepONet)-based surrogate model to identify TAA contributing factors using synthetic finite-element-based datasets. For training, we employ a constrained mixture model of aortic growth and remodelling to generate maps of local aortic dilatation and distensibility for multiple TAA risk factors. We evaluate the performance of the surrogate model for insult distributions varying from fusiform (analytically defined) to complex (randomly generated). We propose two frameworks, one trained on sparse information and one on full-field greyscale images, to gain insight into a preferred neural operator-based approach. We show that this continuous learning approach can predict the patient-specific insult profile associated with any given dilatation and distensibility map with high accuracy, particularly when based on full-field images. Our findings demonstrate the feasibility of applying DeepONet to support transfer learning of patient-specific inputs to predict TAA progression.
Collapse
Affiliation(s)
- Somdatta Goswami
- Division of Applied Mathematics, Brown University, Providence, RI, USA
| | - David S. Li
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA
| | - Bruno V. Rego
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA
| | - Marcos Latorre
- Centre for Research and Innovation in Bioengineering, Universitat Politècnica de València, València, Spain
| | - Jay D. Humphrey
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA
| | - George Em Karniadakis
- Division of Applied Mathematics, Brown University, Providence, RI, USA
- School of Engineering, Brown University, Providence, RI, USA
| |
Collapse
|
5
|
Ng E, Looi LJC. Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:10213-10251. [PMID: 36031992 DOI: 10.3934/mbe.2022479] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Abdominal aortic aneurysms are serious and difficult to detect, conditions can be deadly if they rupture. In this study, the heat transfer and flow physics of Abdominal Aortic Aneurysm (AAA) were discussed and associated with cardiac cycle to illustrate the cardiac thermal pulse (CTP) of AAA. A CTP and infrared thermography (IRT) evaluation-based on AAA and abdomen skin surface detection method was proposed, respectively. Infrared thermography (IRT) is a promising imaging technique that may detect AAA quicker and cheaper than other imaging techniques (as biomarker). From CFD rigid-wall and FSI Analysis, the transient bioheat transfer effect resulted in a distinct thermal signature (circular thermal elevation) on the temperature profile of midriff skin surface, at both regular body temperature and supine position, under normal clinical temperature. However, it is important to note that thermography is not a perfect technology, and it does have some limitations, such as lack of clinical trials. There is still work to be done to improve this imaging technique and make it a more viable and accurate method for detecting abdominal aortic aneurysms. However, thermography is currently one of the most convenient technologies in this field, and it has the potential to detect abdominal aortic aneurysms earlier than other techniques. CTP, on the other hand, was used to examine the thermal physics of AAA. In CFD rigid-wall Analysis, AAA had a CTP that only responded to systolic phase at regular body temperature. In contrast, a healthy abdominal aorta displayed a CTP that responded to the full cardiac cycle, including diastolic phase at all simulated cases. Besides, the findings from FSI Analysis suggest the influence of numerical simulation techniques on the prediction of thermal physics behaviours of AAA and abdominal skin surface. Lastly, this study correlated the relationship between natural convective heat transfer coefficient with AAA and provided reference for potential clinical diagnostic using IRT in clinical implications.
Collapse
Affiliation(s)
- Eyk Ng
- School of Mechanical and Aerospace Engineering, College of Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798
| | - Leonard Jun Cong Looi
- School of Mechanical and Aerospace Engineering, College of Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798
| |
Collapse
|
6
|
Haughton J, Cotter SL, Parnell WJ, Shearer T. Bayesian inference on a microstructural, hyperelastic model of tendon deformation. J R Soc Interface 2022; 19:20220031. [PMID: 35582809 PMCID: PMC9114946 DOI: 10.1098/rsif.2022.0031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2022] [Accepted: 04/25/2022] [Indexed: 11/30/2022] Open
Abstract
Microstructural models of soft-tissue deformation are important in applications including artificial tissue design and surgical planning. The basis of these models, and their advantage over their phenomenological counterparts, is that they incorporate parameters that are directly linked to the tissue's microscale structure and constitutive behaviour and can therefore be used to predict the effects of structural changes to the tissue. Although studies have attempted to determine such parameters using diverse, state-of-the-art, experimental techniques, values ranging over several orders of magnitude have been reported, leading to uncertainty in the true parameter values and creating a need for models that can handle such uncertainty. We derive a new microstructural, hyperelastic model for transversely isotropic soft tissues and use it to model the mechanical behaviour of tendons. To account for parameter uncertainty, we employ a Bayesian approach and apply an adaptive Markov chain Monte Carlo algorithm to determine posterior probability distributions for the model parameters. The obtained posterior distributions are consistent with parameter measurements previously reported and enable us to quantify the uncertainty in their values for each tendon sample that was modelled. This approach could serve as a prototype for quantifying parameter uncertainty in other soft tissues.
Collapse
Affiliation(s)
- James Haughton
- Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
| | - Simon L. Cotter
- Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
| | - William J. Parnell
- Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
| | - Tom Shearer
- Department of Mathematics, University of Manchester, Manchester M13 9PL, UK
- Department of Materials, University of Manchester, Manchester M13 9PL, UK
| |
Collapse
|
7
|
Nabaei M. Cerebral aneurysm evolution modeling from microstructural computational models to machine learning: A review. Comput Biol Chem 2022; 98:107676. [DOI: 10.1016/j.compbiolchem.2022.107676] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2021] [Revised: 02/13/2022] [Accepted: 03/30/2022] [Indexed: 11/03/2022]
|
8
|
Laubrie JD, Mousavi SJ, Avril S. About prestretch in homogenized constrained mixture models simulating growth and remodeling in patient-specific aortic geometries. Biomech Model Mechanobiol 2022; 21:455-469. [PMID: 35067825 PMCID: PMC8940846 DOI: 10.1007/s10237-021-01544-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Accepted: 10/01/2021] [Indexed: 11/10/2022]
Abstract
Evolution of mechanical and structural properties in the Ascending Thoracic Aorta (ATA) is the results of complex mechanobiological processes. In this work, we address some numerical challenges in order to elaborate computational models of these processes. For that, we extend the state of the art of homogenized constrained mixture (hCM) models. In these models, prestretches are assigned to the mixed constituents in order to ensure local mechanical equilibrium macroscopically, and to maintain a homeostatic level of tension in collagen fibers microscopically. Although the initial prestretches were assumed as homogeneous in idealized straight tubes, more elaborate prestretch distributions need to be considered for curved geometrical models such as patient-specific ATA. Therefore, we introduce prestretches having a three-dimensional gradient across the ATA geometry in the homeostatic reference state. We test different schemes with the objective to ensure stable growth and remodeling (G&R) simulations on patient-specific curved vessels. In these simulations, aneurysm progression is triggered by tissue changes in the constituents such as mass degradation of intramural elastin. The results show that the initial prestretches are not only critical for the stability of numerical simulations, but they also affect the G&R response. Eventually, we submit that initial conditions required for G&R simulations need to be identified regionally for ensuring realistic patient-specific predictions of aneurysm progression.
Collapse
|
9
|
Volumetric growth of soft tissues evaluated in the current configuration. Biomech Model Mechanobiol 2022; 21:569-588. [PMID: 35044527 PMCID: PMC8940838 DOI: 10.1007/s10237-021-01549-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2021] [Accepted: 12/17/2021] [Indexed: 11/02/2022]
Abstract
AbstractThe growth and remodelling of soft tissues plays a significant role in many physiological applications, particularly in understanding and managing many diseases. A commonly used approach for soft tissue growth and remodelling is volumetric growth theory, introduced in the framework of finite elasticity. In such an approach, the total deformation gradient tensor is decomposed so that the elastic and growth tensors can be studied separately. A critical element in this approach is to determine the growth tensor and its evolution with time. Most existing volumetric growth theories define the growth tensor in the reference (natural) configuration, which does not reflect the continuous adaptation processes of soft tissues under the current configuration. In a few studies where growth from a loaded configuration was considered, simplifying assumptions, such as compatible deformation or geometric symmetries, were introduced. In this work, we propose a new volumetric growth law that depends on fields evaluated in the current configuration, which is residually stressed and loaded, without any geometrical restrictions. We illustrate our idea using a simplified left ventricle model, which admits inhomogeneous growth in the current configuration. We compare the residual stress distribution of our approach with the traditional volumetric growth theory, that assumes growth occurring from the natural reference configuration. We show that the proposed framework leads to qualitative agreements with experimental measurements. Furthermore, using a cylindrical model, we find an incompatibility index that explains the differences between the two approaches in more depth. We also demonstrate that results from both approaches reach the same steady solution published previously at the limit of a saturated growth. Although we used a left ventricle model as an example, our theory is applicable in modelling the volumetric growth of general soft tissues.
Collapse
|
10
|
Vastmans J, Maes L, Peirlinck M, Vanderveken E, Rega F, Kuhl E, Famaey N. Growth and remodeling in the pulmonary autograft: Computational evaluation using kinematic growth models and constrained mixture theory. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2022; 38:e3545. [PMID: 34724357 DOI: 10.1002/cnm.3545] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 10/26/2021] [Accepted: 10/26/2021] [Indexed: 06/13/2023]
Abstract
Computational investigations of how soft tissues grow and remodel are gaining more and more interest and several growth and remodeling theories have been developed. Roughly, two main groups of theories for soft tissues can be distinguished: kinematic-based growth theory and theories based on constrained mixture theory. Our goal was to apply these two theories on the same experimental data. Within the experiment, a pulmonary artery was exposed to systemic conditions. The change in diameter was followed-up over time. A mechanical and microstructural analysis of native pulmonary artery and pulmonary autograft was conducted. Whereas the kinematic-based growth theory is able to accurately capture the growth of the tissue, it does not account for the mechanobiological processes causing this growth. The constrained mixture theory takes into account the mechanobiological processes including removal, deposition and adaptation of all structural constituents, allowing us to simulate a changing microstructure and mechanical behavior.
Collapse
Affiliation(s)
- Julie Vastmans
- Biomechanics Section, Mechanical Engineering Department, KU Leuven, Leuven, Belgium
| | - Lauranne Maes
- Biomechanics Section, Mechanical Engineering Department, KU Leuven, Leuven, Belgium
| | - Mathias Peirlinck
- Department of Mechanical Engineering, Stanford University, Stanford, California, USA
- IBiTech-bioMMeda, Department of Electronics and Information Systems, Ghent University, Ghent, Belgium
| | - Emma Vanderveken
- Department of Cardiovascular Sciences, KU Leuven, Leuven, Belgium
| | - Filip Rega
- Department of Cardiovascular Sciences, KU Leuven, Leuven, Belgium
| | - Ellen Kuhl
- Department of Mechanical Engineering, Stanford University, Stanford, California, USA
| | - Nele Famaey
- Biomechanics Section, Mechanical Engineering Department, KU Leuven, Leuven, Belgium
| |
Collapse
|
11
|
Humphrey JD. Constrained Mixture Models of Soft Tissue Growth and Remodeling - Twenty Years After. JOURNAL OF ELASTICITY 2021; 145:49-75. [PMID: 34483462 PMCID: PMC8415366 DOI: 10.1007/s10659-020-09809-1] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Accepted: 12/29/2020] [Indexed: 05/06/2023]
Abstract
Soft biological tissues compromise diverse cell types and extracellular matrix constituents, each of which can possess individual natural configurations, material properties, and rates of turnover. For this reason, mixture-based models of growth (changes in mass) and remodeling (change in microstructure) are well-suited for studying tissue adaptations, disease progression, and responses to injury or clinical intervention. Such approaches also can be used to design improved tissue engineered constructs to repair, replace, or regenerate tissues. Focusing on blood vessels as archetypes of soft tissues, this paper reviews a constrained mixture theory introduced twenty years ago and explores its usage since by contrasting simulations of diverse vascular conditions. The discussion is framed within the concept of mechanical homeostasis, with consideration of solid-fluid interactions, inflammation, and cell signaling highlighting both past accomplishments and future opportunities as we seek to understand better the evolving composition, geometry, and material behaviors of soft tissues under complex conditions.
Collapse
Affiliation(s)
- J D Humphrey
- Department of Biomedical Engineering, Yale University, New Haven, CT 06520 USA
| |
Collapse
|
12
|
Eichinger JF, Haeusel LJ, Paukner D, Aydin RC, Humphrey JD, Cyron CJ. Mechanical homeostasis in tissue equivalents: a review. Biomech Model Mechanobiol 2021; 20:833-850. [PMID: 33683513 PMCID: PMC8154823 DOI: 10.1007/s10237-021-01433-9] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 02/04/2021] [Indexed: 12/20/2022]
Abstract
There is substantial evidence that growth and remodeling of load bearing soft biological tissues is to a large extent controlled by mechanical factors. Mechanical homeostasis, which describes the natural tendency of such tissues to establish, maintain, or restore a preferred mechanical state, is thought to be one mechanism by which such control is achieved across multiple scales. Yet, many questions remain regarding what promotes or prevents homeostasis. Tissue equivalents, such as collagen gels seeded with living cells, have become an important tool to address these open questions under well-defined, though limited, conditions. This article briefly reviews the current state of research in this area. It summarizes, categorizes, and compares experimental observations from the literature that focus on the development of tension in tissue equivalents. It focuses primarily on uniaxial and biaxial experimental studies, which are well-suited for quantifying interactions between mechanics and biology. The article concludes with a brief discussion of key questions for future research in this field.
Collapse
Affiliation(s)
- Jonas F Eichinger
- Institute for Computational Mechanics, Technical University of Munich, 85748, Munich, Germany.,Institute of Continuum and Materials Mechanics, Hamburg University of Technology, 21073, Hamburg, Germany
| | - Lea J Haeusel
- Institute for Computational Mechanics, Technical University of Munich, 85748, Munich, Germany
| | - Daniel Paukner
- Institute of Continuum and Materials Mechanics, Hamburg University of Technology, 21073, Hamburg, Germany.,Institute of Material Systems Modeling, Helmholtz-Zentrum Geesthacht, 21502, Geesthacht, Germany
| | - Roland C Aydin
- Institute of Material Systems Modeling, Helmholtz-Zentrum Geesthacht, 21502, Geesthacht, Germany
| | - Jay D Humphrey
- Department of Biomedical Engineering, Yale University, New Haven, CT, 06520, USA
| | - Christian J Cyron
- Institute of Continuum and Materials Mechanics, Hamburg University of Technology, 21073, Hamburg, Germany. .,Institute of Material Systems Modeling, Helmholtz-Zentrum Geesthacht, 21502, Geesthacht, Germany.
| |
Collapse
|
13
|
Teixeira FS, Neufeld E, Kuster N, Watton PN. Modeling intracranial aneurysm stability and growth: an integrative mechanobiological framework for clinical cases. Biomech Model Mechanobiol 2020; 19:2413-2431. [PMID: 32533497 PMCID: PMC7603456 DOI: 10.1007/s10237-020-01351-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Accepted: 05/12/2020] [Indexed: 11/03/2022]
Abstract
We present a novel patient-specific fluid-solid-growth framework to model the mechanobiological state of clinically detected intracranial aneurysms (IAs) and their evolution. The artery and IA sac are modeled as thick-walled, non-linear elastic fiber-reinforced composites. We represent the undulation distribution of collagen fibers: the adventitia of the healthy artery is modeled as a protective sheath whereas the aneurysm sac is modeled to bear load within physiological range of pressures. Initially, we assume the detected IA is stable and then consider two flow-related mechanisms to drive enlargement: (1) low wall shear stress; (2) dysfunctional endothelium which is associated with regions of high oscillatory flow. Localized collagen degradation and remodelling gives rise to formation of secondary blebs on the aneurysm dome. Restabilization of blebs is achieved by remodelling of the homeostatic collagen fiber stretch distribution. This integrative mechanobiological modelling workflow provides a step towards a personalized risk-assessment and treatment of clinically detected IAs.
Collapse
Affiliation(s)
| | - Esra Neufeld
- IT’IS Foundation & ETH Zürich, Zürich, Switzerland
| | - Niels Kuster
- IT’IS Foundation & ETH Zürich, Zürich, Switzerland
| | - Paul N. Watton
- Department of Computer Science, Insigneo Institute for in silico Medicine, University of Sheffield, Sheffield, UK
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, USA
| |
Collapse
|
14
|
Ghavamian A, Mousavi SJ, Avril S. Computational Study of Growth and Remodeling in Ascending Thoracic Aortic Aneurysms Considering Variations of Smooth Muscle Cell Basal Tone. Front Bioeng Biotechnol 2020; 8:587376. [PMID: 33224937 PMCID: PMC7670047 DOI: 10.3389/fbioe.2020.587376] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Accepted: 09/28/2020] [Indexed: 11/16/2022] Open
Abstract
In this paper, we investigate the progression of Ascending Thoracic Aortic Aneurysms (ATAA) using a computational model of Growth and Remodeling (G&R) taking into account the composite (elastin, four collagen fiber families and Smooth Muscle Cells—SMCs) and multi-layered (media and adventitia) nature of the aorta. The G&R model, which is based on the homogenized Constrained Mixture theory, is implemented as a UMAT in the Abaqus finite-element package. Each component of the mixture is assigned a strain energy density function: nearly-incompressible neo-Hookean for elastin and Fung-type for collagen and SMCs. Active SMCs tension is additionally considered, through a length-tension relationship having a classic inverted parabola shape, in order to investigate its effects on the progression of ATAA in a patient-specific model. A sensitivity analysis is performed to evaluate the potential impact of variations in the parameters of the length-tension relationships. These variations reflect in variations of SMCs normal tone during ATAA progression, with active stress contributions ranging between 30% (best case scenario) and 0% (worst case scenario) of the total wall circumferential stress. Low SMCs active stress in the worst case scenarios, in fact, affect the rates of collagen deposition by which the elastin loss is gradually compensated by collagen deposition in the simulated ATAA progression, resulting eventually in larger aneurysm diameters. The types of length-tension relationships leading to a drop of SMCs active stress in our simulations reveal a critical condition which could also result in SMCs apoptosis.
Collapse
Affiliation(s)
- Ataollah Ghavamian
- Mines Saint-Etienne, Université Lyon, Université Jean Monnet, INSERM, U 1059 Sainbiose, Centre CIS, Saint-Étienne, France
| | - S Jamaleddin Mousavi
- Mines Saint-Etienne, Université Lyon, Université Jean Monnet, INSERM, U 1059 Sainbiose, Centre CIS, Saint-Étienne, France
| | - Stéphane Avril
- Mines Saint-Etienne, Université Lyon, Université Jean Monnet, INSERM, U 1059 Sainbiose, Centre CIS, Saint-Étienne, France
| |
Collapse
|
15
|
Latorre M, Humphrey JD. Fast, Rate-Independent, Finite Element Implementation of a 3D Constrained Mixture Model of Soft Tissue Growth and Remodeling. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2020; 368:113156. [PMID: 32655195 PMCID: PMC7351114 DOI: 10.1016/j.cma.2020.113156] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Constrained mixture models of soft tissue growth and remodeling can simulate many evolving conditions in health as well as in disease and its treatment, but they can be computationally expensive. In this paper, we derive a new fast, robust finite element implementation based on a concept of mechanobiological equilibrium that yields fully resolved solutions and allows computation of quasi-equilibrated evolutions when imposed perturbations are slow relative to the adaptive process. We demonstrate quadratic convergence and verify the model via comparisons with semi-analytical solutions for arterial mechanics. We further examine the enlargement of aortic aneurysms for which we identify new mechanobiological insights into factors that affect the nearby non-aneurysmal segment as it responds to the changing mechanics within the diseased segment. Because this new 3D approach can be implemented within many existing finite element solvers, constrained mixture models of growth and remodeling can now be used more widely.
Collapse
Affiliation(s)
- Marcos Latorre
- Department of Biomedical Engineering Yale University, New Haven, CT, USA
| | - Jay D. Humphrey
- Department of Biomedical Engineering Yale University, New Haven, CT, USA
- Vascular Biology and Therapeutics Program Yale School of Medicine, New Haven, CT, USA
| |
Collapse
|
16
|
Guglietta F, Behr M, Biferale L, Falcucci G, Sbragaglia M. On the effects of membrane viscosity on transient red blood cell dynamics. SOFT MATTER 2020; 16:6191-6205. [PMID: 32567630 DOI: 10.1039/d0sm00587h] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Computational Fluid Dynamics (CFD) is currently used to design and improve the hydraulic properties of biomedical devices, wherein the large scale blood circulation needs to be simulated by accounting for the mechanical response of red blood cells (RBCs) at the mesoscale. In many practical instances, biomedical devices work on time-scales comparable to the intrinsic relaxation time of RBCs: thus, a systematic understanding of the time-dependent response of erythrocyte membranes is crucial for the effective design of such devices. So far, this information has been deduced from experimental data, which do not necessarily adapt to the broad variety of fluid dynamic conditions that can be encountered in practice. This work explores the novel possibility of studying the time-dependent response of an erythrocyte membrane to external mechanical loads via mesoscale numerical simulations, with a primary focus on the detailed characterisation of the RBC relaxation time tc following the arrest of the external mechanical load. The adopted mesoscale model exploits a hybrid Immersed Boundary-Lattice Boltzmann Method (IB-LBM), coupled with the Standard Linear Solid (SLS) model to account for the RBC membrane viscosity. We underscore the key importance of the 2D membrane viscosity μm to correctly reproduce the relaxation time of the RBC membrane. A detailed assessment of the dependencies on the typology and strength of the applied mechanical loads is also provided. Overall, our findings open interesting future perspectives for the study of the non-linear response of RBCs immersed in time-dependent strain fields.
Collapse
Affiliation(s)
- Fabio Guglietta
- Department of Physics & INFN, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy. and Chair for Computational Analysis of Technical Systems (CATS), RWTH Aachen University, 52056 Aachen, Germany and Computation-Based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Str., 2121 Nicosia, Cyprus
| | - Marek Behr
- Chair for Computational Analysis of Technical Systems (CATS), RWTH Aachen University, 52056 Aachen, Germany
| | - Luca Biferale
- Department of Physics & INFN, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy.
| | - Giacomo Falcucci
- Department of Enterprise Engineering "Mario Lucertini", University of Rome "Tor Vergata", Via del Politecnico 1, 00133 Rome, Italy and John A. Paulson School of Engineering and Applied Physics, Harvard University, 33 Oxford Street, 02138 Cambridge, Massachusetts, USA
| | - Mauro Sbragaglia
- Department of Physics & INFN, University of Rome "Tor Vergata", Via della Ricerca Scientifica 1, 00133, Rome, Italy.
| |
Collapse
|
17
|
Niestrawska JA, Augustin CM, Plank G. Computational modeling of cardiac growth and remodeling in pressure overloaded hearts-Linking microstructure to organ phenotype. Acta Biomater 2020; 106:34-53. [PMID: 32058078 PMCID: PMC7311197 DOI: 10.1016/j.actbio.2020.02.010] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2019] [Revised: 02/06/2020] [Accepted: 02/07/2020] [Indexed: 12/25/2022]
Abstract
Cardiac growth and remodeling (G&R) refers to structural changes in myocardial tissue in response to chronic alterations in loading conditions. One such condition is pressure overload where elevated wall stresses stimulate the growth in cardiomyocyte thickness, associated with a phenotype of concentric hypertrophy at the organ scale, and promote fibrosis. The initial hypertrophic response can be considered adaptive and beneficial by favoring myocyte survival, but over time if pressure overload conditions persist, maladaptive mechanisms favoring cell death and fibrosis start to dominate, ultimately mediating the transition towards an overt heart failure phenotype. The underlying mechanisms linking biological factors at the myocyte level to biomechanical factors at the systemic and organ level remain poorly understood. Computational models of G&R show high promise as a unique framework for providing a quantitative link between myocardial stresses and strains at the organ scale to biological regulatory processes at the cellular level which govern the hypertrophic response. However, microstructurally motivated, rigorously validated computational models of G&R are still in their infancy. This article provides an overview of the current state-of-the-art of computational models to study cardiac G&R. The microstructure and mechanosensing/mechanotransduction within cells of the myocardium is discussed and quantitative data from previous experimental and clinical studies is summarized. We conclude with a discussion of major challenges and possible directions of future research that can advance the current state of cardiac G&R computational modeling. STATEMENT OF SIGNIFICANCE: The mechanistic links between organ-scale biomechanics and biological factors at the cellular size scale remain poorly understood as these are largely elusive to investigations using experimental methodology alone. Computational G&R models show high promise to establish quantitative links which allow more mechanistic insight into adaptation mechanisms and may be used as a tool for stratifying the state and predict the progression of disease in the clinic. This review provides a comprehensive overview of research in this domain including a summary of experimental data. Thus, this study may serve as a basis for the further development of more advanced G&R models which are suitable for making clinical predictions on disease progression or for testing hypotheses on pathogenic mechanisms using in-silico models.
Collapse
Affiliation(s)
- Justyna A Niestrawska
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz 8010, Austria
| | - Christoph M Augustin
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz 8010, Austria.
| | - Gernot Plank
- Gottfried Schatz Research Center: Division of Biophysics, Medical University of Graz, Graz 8010, Austria; BioTechMed-Graz, Austria
| |
Collapse
|
18
|
Laubrie JD, Mousavi JS, Avril S. A new finite-element shell model for arterial growth and remodeling after stent implantation. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2020; 36:e3282. [PMID: 31773919 DOI: 10.1002/cnm.3282] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 10/05/2019] [Accepted: 10/27/2019] [Indexed: 06/10/2023]
Abstract
The goal of this paper is to study computationally how blood vessels adapt when they are exposed to a mechanobiological insult, namely, a sudden change of their biomechanical conditions such as proteolytic injuries or implantation. Adaptation occurs through growth and remodeling (G&R), consisting of mass production or removal of structural proteins, such as collagen, until restoring the initial homeostatic biomechanical conditions. In some circumstances, the initial conditions can never be recovered, and arteries evolve towards unstable pathological conditions, such as aneurysms, which are responsible for significant morbidity and mortality. Therefore, computational predictions of G&R under different circumstances can be helpful in understanding fundamentally how arterial pathologies progress. For that, we have developed a low-cost open-source finite-element 2D axisymmetric shell model (FEM) of the arterial wall. The constitutive equations for static equilibrium used to model the stress-strain behavior and the G&R response are expressed within the homogenized constrained mixture theory. The originality is to integrate the layer-specific behavior of both arterial layers (media and adventitia) into the model. Considering different mechanobiological insults, our results show that the resulting arterial dilatation is strongly correlated with the media thickness. The adaptation to stent implantation is particularly interesting. For large stent oversizing ratios, the artery cannot recover from the mechanobiological insult and dilates forever, whereas dilatation stabilizes after a transient period for more moderate oversizing ratios. We also show that stent implantation induces a different response in an aneurysm or in a healthy artery, the latter yielding more unstable G&R. Finally, our G&R model can efficiently predict, with very low computational cost, fundamental aspects of arterial adaptation induced by clinical procedures.
Collapse
Affiliation(s)
- Joan D Laubrie
- Mines Saint-Étienne, Univ Lyon, Univ Jean Monnet, INSERM, U 1059 Sainbiose, Centre CIS, F - 42023, Saint-Étienne, France
| | - Jamaleddin S Mousavi
- Mines Saint-Étienne, Univ Lyon, Univ Jean Monnet, INSERM, U 1059 Sainbiose, Centre CIS, F - 42023, Saint-Étienne, France
| | - Stéphane Avril
- Mines Saint-Étienne, Univ Lyon, Univ Jean Monnet, INSERM, U 1059 Sainbiose, Centre CIS, F - 42023, Saint-Étienne, France
| |
Collapse
|
19
|
Patient-specific predictions of aneurysm growth and remodeling in the ascending thoracic aorta using the homogenized constrained mixture model. Biomech Model Mechanobiol 2019; 18:1895-1913. [DOI: 10.1007/s10237-019-01184-8] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Accepted: 06/05/2019] [Indexed: 12/19/2022]
|
20
|
Salman HE, Ramazanli B, Yavuz MM, Yalcin HC. Biomechanical Investigation of Disturbed Hemodynamics-Induced Tissue Degeneration in Abdominal Aortic Aneurysms Using Computational and Experimental Techniques. Front Bioeng Biotechnol 2019; 7:111. [PMID: 31214581 PMCID: PMC6555197 DOI: 10.3389/fbioe.2019.00111] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2019] [Accepted: 05/02/2019] [Indexed: 11/13/2022] Open
Abstract
Abdominal aortic aneurysm (AAA) is the dilatation of the aorta beyond 50% of the normal vessel diameter. It is reported that 4-8% of men and 0.5-1% of women above 50 years of age bear an AAA and it accounts for ~15,000 deaths per year in the United States alone. If left untreated, AAA might gradually expand until rupture; the most catastrophic complication of the aneurysmal disease that is accompanied by a striking overall mortality of 80%. The precise mechanisms leading to AAA rupture remains unclear. Therefore, characterization of disturbed hemodynamics within AAAs will help to understand the mechanobiological development of the condition which will contribute to novel therapies for the condition. Due to geometrical complexities, it is challenging to directly quantify disturbed flows for AAAs clinically. Two other approaches for this investigation are computational modeling and experimental flow measurement. In computational modeling, the problem is first defined mathematically, and the solution is approximated with numerical techniques to get characteristics of flow. In experimental flow measurement, once the setup providing physiological flow pattern in a phantom geometry is constructed, velocity measurement system such as particle image velocimetry (PIV) enables characterization of the flow. We witness increasing number of applications of these complimentary approaches for AAA investigations in recent years. In this paper, we outline the details of computational modeling procedures and experimental settings and summarize important findings from recent studies, which will help researchers for AAA investigations and rupture mechanics.
Collapse
Affiliation(s)
| | - Burcu Ramazanli
- Department of Mechanical Engineering, Middle East Technical University, Ankara, Turkey
| | - Mehmet Metin Yavuz
- Department of Mechanical Engineering, Middle East Technical University, Ankara, Turkey
| | | |
Collapse
|
21
|
Latorre M, Humphrey JD. Mechanobiological Stability of Biological Soft Tissues. JOURNAL OF THE MECHANICS AND PHYSICS OF SOLIDS 2019; 125:298-325. [PMID: 31543547 PMCID: PMC6754118 DOI: 10.1016/j.jmps.2018.12.013] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Like all other materials, biological soft tissues are subject to general laws of physics, including those governing mechanical equilibrium and stability. In addition, however, these tissues are able to respond actively to changes in their mechanical and chemical environment. There is, therefore, a pressing need to understand such processes theoretically. In this paper, we present a new rate-based constrained mixture formulation suitable for studying mechanobiological equilibrium and stability of soft tissues exposed to transient or sustained changes in material composition or applied loading. These concepts are illustrated for canonical problems in arterial mechanics, which distinguish possible stable versus unstable mechanobiological responses. Such analyses promise to yield insight into biological processes that govern both health and disease progression.
Collapse
Affiliation(s)
- Marcos Latorre
- Department of Biomedical Engineering Yale University, New Haven, CT 06520, USA
- Corresponding author: (Marcos Latorre), (Jay D. Humphrey)
| | - Jay D. Humphrey
- Department of Biomedical Engineering Yale University, New Haven, CT 06520, USA
- Vascular Biology and Therapeutics Program Yale School of Medicine, New Haven, CT 06520, USA
| |
Collapse
|
22
|
Bhogal P, Pederzani G, Grytsan A, Loh Y, Brouwer PA, Andersson T, Gundiah N, Robertson AM, Watton PN, Söderman M. The unexplained success of stentplasty vasospasm treatment : Insights using Mechanistic Mathematical Modeling. Clin Neuroradiol 2019; 29:763-774. [PMID: 30915482 DOI: 10.1007/s00062-019-00776-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2018] [Accepted: 03/08/2019] [Indexed: 11/29/2022]
Abstract
BACKGROUND Cerebral vasospasm (CVS) following subarachnoid hemorrhage occurs in up to 70% of patients. Recently, stents have been used to successfully treat CVS. This implies that the force required to expand spastic vessels and resolve vasospasm is lower than previously thought. OBJECTIVE We develop a mechanistic model of the spastic arterial wall to provide insight into CVS and predict the forces required to treat it. MATERIAL AND METHODS The arterial wall is modelled as a cylindrical membrane using a constrained mixture theory that accounts for the mechanical roles of elastin, collagen and vascular smooth muscle cells (VSMC). We model the pressure diameter curve prior to CVS and predict how it changes following CVS. We propose a stretch-based damage criterion for VSMC and evaluate if several commercially available stents are able to resolve vasospasm. RESULTS The model predicts that dilatation of VSMCs beyond a threshold of mechanical failure is sufficient to resolve CVS without damage to the underlying extracellular matrix. Consistent with recent clinical observations, our model predicts that existing stents have the potential to provide sufficient outward force to successfully treat CVS and that success will be dependent on an appropriate match between stent and vessel. CONCLUSION Mathematical models of CVS can provide insights into biological mechanisms and explore treatment approaches. Improved understanding of the underlying mechanistic processes governing CVS and its mechanical treatment may assist in the development of dedicated stents.
Collapse
Affiliation(s)
- P Bhogal
- Department of Interventional Neuroradiology, The Royal London Hospital, Whitechapel Road, E1 1BB, London, UK.
| | - G Pederzani
- Department of Computer Science, University of Sheffield, Sheffield, UK.,Insigneo Institute for in silico Medicine, University of Sheffield, Sheffield, UK
| | - A Grytsan
- Insigneo Institute for in silico Medicine, University of Sheffield, Sheffield, UK
| | - Y Loh
- Uniformed Services University, University of California, Los Angeles, USA.,Swedish Neuroscience Institute, 550 17th Avenue Seattle, 98122, Washington, USA
| | - P A Brouwer
- The Karolinska University Hospital, 171 76, Stockholm, Sweden
| | - T Andersson
- The Karolinska University Hospital, 171 76, Stockholm, Sweden.,AZ Groeninge, Kortrijk, Belgium
| | - Namrata Gundiah
- Department of Mechanical Engineering, Indian Institute of Science, Bangalore, India
| | - Anne M Robertson
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, PA, USA
| | - Paul N Watton
- Department of Computer Science, University of Sheffield, Sheffield, UK.,Insigneo Institute for in silico Medicine, University of Sheffield, Sheffield, UK.,Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, PA, USA
| | | |
Collapse
|
23
|
Braeu FA, Aydin RC, Cyron CJ. Anisotropic stiffness and tensional homeostasis induce a natural anisotropy of volumetric growth and remodeling in soft biological tissues. Biomech Model Mechanobiol 2018; 18:327-345. [PMID: 30413985 DOI: 10.1007/s10237-018-1084-x] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 10/16/2018] [Indexed: 12/20/2022]
Abstract
Growth in soft biological tissues in general results in anisotropic changes of the tissue geometry. It remains a key challenge in biomechanics to understand, quantify, and predict this anisotropy. In this paper, we demonstrate that anisotropic tissue stiffness and the well-known mechanism of tensional homeostasis induce a natural anisotropy of the geometric changes resulting from volumetric growth in soft biological tissues. As a rule of thumb, this natural anisotropy makes differential tissue volume elements dilate mainly in the direction(s) of lowest stiffness. This simple principle is shown to explain the experimentally observed growth behavior in a host of different soft biological tissues without relying on any additional heuristic assumptions or quantities (such as ad hoc defined growth tensors).
Collapse
Affiliation(s)
- F A Braeu
- Institute for Computational Mechanics, Technical University of Munich, Munich, Germany
| | - R C Aydin
- Institute for Computational Mechanics, Technical University of Munich, Munich, Germany
| | - Christian J Cyron
- Institute of Continuum Mechanics and Materials Mechanics, Hamburg University of Technology, Eissendorfer Strasse 42, 21073, Hamburg, Germany. .,Institute of Materials Research, Materials Mechanics, Helmholtz-Zentrum Geesthacht, Geesthacht, Germany.
| |
Collapse
|
24
|
Owen B, Bojdo N, Jivkov A, Keavney B, Revell A. Structural modelling of the cardiovascular system. Biomech Model Mechanobiol 2018; 17:1217-1242. [PMID: 29911296 PMCID: PMC6154127 DOI: 10.1007/s10237-018-1024-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Accepted: 04/25/2018] [Indexed: 02/02/2023]
Abstract
Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.
Collapse
Affiliation(s)
- Benjamin Owen
- School of Mechanical, Aerospace and Civil Engineering, University of Manchester, George Begg Building, Manchester, M1 3BB, UK.
| | - Nicholas Bojdo
- School of Mechanical, Aerospace and Civil Engineering, University of Manchester, George Begg Building, Manchester, M1 3BB, UK
| | - Andrey Jivkov
- School of Mechanical, Aerospace and Civil Engineering, University of Manchester, George Begg Building, Manchester, M1 3BB, UK
| | - Bernard Keavney
- Division of Cardiovascular Sciences, University of Manchester, AV Hill Building, Manchester, M13 9PT, UK
| | - Alistair Revell
- School of Mechanical, Aerospace and Civil Engineering, University of Manchester, George Begg Building, Manchester, M1 3BB, UK
| |
Collapse
|
25
|
Bevan T, Merabet N, Hornsby J, Watton PN, Thompson MS. A biomechanical model for fibril recruitment: Evaluation in tendons and arteries. J Biomech 2018; 74:192-196. [PMID: 29636179 DOI: 10.1016/j.jbiomech.2018.03.047] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Revised: 02/19/2018] [Accepted: 03/24/2018] [Indexed: 11/20/2022]
Abstract
Simulations of soft tissue mechanobiological behaviour are increasingly important for clinical prediction of aneurysm, tendinopathy and other disorders. Mechanical behaviour at low stretches is governed by fibril straightening, transitioning into load-bearing at recruitment stretch, resulting in a tissue stiffening effect. Previous investigations have suggested theoretical relationships between stress-stretch measurements and recruitment probability density function (PDF) but not derived these rigorously nor evaluated these experimentally. Other work has proposed image-based methods for measurement of recruitment but made use of arbitrary fibril critical straightness parameters. The aim of this work was to provide a sound theoretical basis for estimating recruitment PDF from stress-stretch measurements and to evaluate this relationship using image-based methods, clearly motivating the choice of fibril critical straightness parameter in rat tail tendon and porcine artery. Rigorous derivation showed that the recruitment PDF may be estimated from the second stretch derivative of the first Piola-Kirchoff tissue stress. Image-based fibril recruitment identified the fibril straightness parameter that maximised Pearson correlation coefficients (PCC) with estimated PDFs. Using these critical straightness parameters the new method for estimating recruitment PDF showed a PCC with image-based measures of 0.915 and 0.933 for tendons and arteries respectively. This method may be used for accurate estimation of fibril recruitment PDF in mechanobiological simulation where fibril-level mechanical parameters are important for predicting cell behaviour.
Collapse
Affiliation(s)
- Tim Bevan
- Institute of Biomedical Engineering, University of Oxford, United Kingdom
| | - Nadege Merabet
- Institute of Biomedical Engineering, University of Oxford, United Kingdom
| | - Jack Hornsby
- Institute of Biomedical Engineering, University of Oxford, United Kingdom
| | - Paul N Watton
- Department of Computer Science & INSIGNEO Institute for In Silico Medicine, University of Sheffield, United Kingdom; Mechanical Engineering and Materials Science, University of Pittsburgh, USA
| | - Mark S Thompson
- Institute of Biomedical Engineering, University of Oxford, United Kingdom.
| |
Collapse
|
26
|
Kemmerling EMC, Peattie RA. Abdominal Aortic Aneurysm Pathomechanics: Current Understanding and Future Directions. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2018; 1097:157-179. [DOI: 10.1007/978-3-319-96445-4_8] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
27
|
Growth Description for Vessel Wall Adaptation: A Thick-Walled Mixture Model of Abdominal Aortic Aneurysm Evolution. MATERIALS 2017; 10:ma10090994. [PMID: 28841196 PMCID: PMC5615649 DOI: 10.3390/ma10090994] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Revised: 08/21/2017] [Accepted: 08/23/2017] [Indexed: 12/20/2022]
Abstract
(1) Background: Vascular tissue seems to adapt towards stable homeostatic mechanical conditions, however, failure of reaching homeostasis may result in pathologies. Current vascular tissue adaptation models use many ad hoc assumptions, the implications of which are far from being fully understood; (2) Methods: The present study investigates the plausibility of different growth kinematics in modeling Abdominal Aortic Aneurysm (AAA) evolution in time. A structurally motivated constitutive description for the vessel wall is coupled to multi-constituent tissue growth descriptions; Constituent deposition preserved either the constituent’s density or its volume, and Isotropic Volume Growth (IVG), in-Plane Volume Growth (PVG), in-Thickness Volume Growth (TVG) and No Volume Growth (NVG) describe the kinematics of the growing vessel wall. The sensitivity of key modeling parameters is explored, and predictions are assessed for their plausibility; (3) Results: AAA development based on TVG and NVG kinematics provided not only quantitatively, but also qualitatively different results compared to IVG and PVG kinematics. Specifically, for IVG and PVG kinematics, increasing collagen mass production accelerated AAA expansion which seems counterintuitive. In addition, TVG and NVG kinematics showed less sensitivity to the initial constituent volume fractions, than predictions based on IVG and PVG; (4) Conclusions: The choice of tissue growth kinematics is of crucial importance when modeling AAA growth. Much more interdisciplinary experimental work is required to develop and validate vascular tissue adaption models, before such models can be of any practical use.
Collapse
|
28
|
Thompson MS, Bajuri MN, Khayyeri H, Isaksson H. Mechanobiological modelling of tendons: Review and future opportunities. Proc Inst Mech Eng H 2017; 231:369-377. [DOI: 10.1177/0954411917692010] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Tendons are adapted to carry large, repeated loads and are clinically important for the maintenance of musculoskeletal health in an increasing, actively ageing population, as well as in elite athletes. Tendons are known to adapt to mechanical loading. Also, their healing and disease processes are highly sensitive to mechanical load. Computational modelling approaches developed to capture this mechanobiological adaptation in tendons and other tissues have successfully addressed many important scientific and clinical issues. The aim of this review is to identify techniques and approaches that could be further developed to address tendon-related problems. Biomechanical models are identified that capture the multi-level aspects of tendon mechanics. Continuum whole tendon models, both phenomenological and microstructurally motivated, are important to estimate forces during locomotion activities. Fibril-level microstructural models are documented that can use these estimated forces to detail local mechanical parameters relevant to cell mechanotransduction. Cell-level models able to predict the response to such parameters are also described. A selection of updatable mechanobiological models is presented. These use mechanical signals, often continuum tissue level, along with rules for tissue change and have been applied successfully in many tissues to predict in vivo and in vitro outcomes. Signals may include scalars derived from the stress or strain tensors, or in poroelasticity also fluid velocity, while adaptation may be represented by changes to elastic modulus, permeability, fibril density or orientation. So far, only simple analytical approaches have been applied to tendon mechanobiology. With the development of sophisticated computational mechanobiological models in parallel with reporting more quantitative data from in vivo or clinical mechanobiological studies, for example, appropriate imaging, biochemical and histological data, this field offers huge potential for future development towards clinical applications.
Collapse
Affiliation(s)
- Mark S Thompson
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, UK
| | - M Nazri Bajuri
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, UK
- Faculty of Biosciences & Medical Engineering, Universiti Teknologi Malaysia, Johor Bahru, Malaysia
| | - Hanifeh Khayyeri
- Department of Biomedical Engineering, Lund University, Lund, Sweden
| | - Hanna Isaksson
- Department of Biomedical Engineering, Lund University, Lund, Sweden
| |
Collapse
|
29
|
Gasser TC, Grytsan A. Biomechanical modeling the adaptation of soft biological tissue. CURRENT OPINION IN BIOMEDICAL ENGINEERING 2017. [DOI: 10.1016/j.cobme.2017.03.004] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
30
|
Abstract
Abdominal aortic aneurysm (AAA) is a localized enlargement of the abdominal aorta, such that the diameter exceeds 3 cm. The natural history of AAA is progressive growth leading to rupture, an event that carries up to 90% risk of mortality. Hence there is a need to predict the growth of the diameter of the aorta based on the diameter of a patient’s aneurysm at initial screening and aided by non-invasive biomarkers. IL-6 is overexpressed in AAA and was suggested as a prognostic marker for the risk in AAA. The present paper develops a mathematical model which relates the growth of the abdominal aorta to the serum concentration of IL-6. Given the initial diameter of the aorta and the serum concentration of IL-6, the model predicts the growth of the diameter at subsequent times. Such a prediction can provide guidance to how closely the patient’s abdominal aorta should be monitored. The mathematical model is represented by a system of partial differential equations taking place in the aortic wall, where the media is assumed to have the constituency of an hyperelastic material.
Collapse
|
31
|
Cyron CJ, Humphrey JD. Growth and Remodeling of Load-Bearing Biological Soft Tissues. MECCANICA 2017; 52:645-664. [PMID: 28286348 PMCID: PMC5342900 DOI: 10.1007/s11012-016-0472-5] [Citation(s) in RCA: 71] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
The past two decades reveal a growing role of continuum biomechanics in understanding homeostasis, adaptation, and disease progression in soft tissues. In this paper, we briefly review the two primary theoretical approaches for describing mechano-regulated soft tissue growth and remodeling on the continuum level as well as hybrid approaches that attempt to combine the advantages of these two approaches while avoiding their disadvantages. We also discuss emerging concepts, including that of mechanobiological stability. Moreover, to motivate and put into context the different theoretical approaches, we briefly review findings from mechanobiology that show the importance of mass turnover and the prestressing of both extant and new extracellular matrix in most cases of growth and remodeling. For illustrative purposes, these concepts and findings are discussed, in large part, within the context of two load-bearing, collagen dominated soft tissues - tendons/ligaments and blood vessels. We conclude by emphasizing further examples, needs, and opportunities in this exciting field of modeling soft tissues.
Collapse
Affiliation(s)
- C J Cyron
- Institute for Computational Mechanics, Technische Universität München, Garching, Germany
| | - J D Humphrey
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA; Vascular Biology and Therapeutics Program, Yale School of Medicine, New Haven, CT, USA
| |
Collapse
|
32
|
Comellas E, Gasser TC, Bellomo FJ, Oller S. A homeostatic-driven turnover remodelling constitutive model for healing in soft tissues. J R Soc Interface 2016; 13:rsif.2015.1081. [PMID: 27009177 DOI: 10.1098/rsif.2015.1081] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2015] [Accepted: 03/01/2016] [Indexed: 01/08/2023] Open
Abstract
Remodelling of soft biological tissue is characterized by interacting biochemical and biomechanical events, which change the tissue's microstructure, and, consequently, its macroscopic mechanical properties. Remodelling is a well-defined stage of the healing process, and aims at recovering or repairing the injured extracellular matrix. Like other physiological processes, remodelling is thought to be driven by homeostasis, i.e. it tends to re-establish the properties of the uninjured tissue. However, homeostasis may never be reached, such that remodelling may also appear as a continuous pathological transformation of diseased tissues during aneurysm expansion, for example. A simple constitutive model for soft biological tissues that regards remodelling as homeostatic-driven turnover is developed. Specifically, the recoverable effective tissue damage, whose rate is the sum of a mechanical damage rate and a healing rate, serves as a scalar internal thermodynamic variable. In order to integrate the biochemical and biomechanical aspects of remodelling, the healing rate is, on the one hand, driven by mechanical stimuli, but, on the other hand, subjected to simple metabolic constraints. The proposed model is formulated in accordance with continuum damage mechanics within an open-system thermodynamics framework. The numerical implementation in an in-house finite-element code is described, particularized for Ogden hyperelasticity. Numerical examples illustrate the basic constitutive characteristics of the model and demonstrate its potential in representing aspects of remodelling of soft tissues. Simulation results are verified for their plausibility, but also validated against reported experimental data.
Collapse
Affiliation(s)
- Ester Comellas
- International Center for Numerical Methods in Engineering (CIMNE), Campus Nord UPC, Building C1, c/Gran Capita s/n, 08034 Barcelona, Spain Department of Strength of Materials and Structural Engineering, ETSECCPB, Universitat Politcnica de Catalunya, Barcelona Tech (UPC), Campus Nord, Building C1, c/Jordi Girona 1-3, 08034 Barcelona, Spain
| | - T Christian Gasser
- Department of Solid Mechanics, School of Engineering Sciences, KTH Royal Institute of Technology, Teknikringen 8, 100 44 Stockholm, Sweden
| | - Facundo J Bellomo
- INIQUI (CONICET), Faculty of Engineering, National University of Salta, Av. Bolivia 5150, 4400 Salta, Argentina
| | - Sergio Oller
- International Center for Numerical Methods in Engineering (CIMNE), Campus Nord UPC, Building C1, c/Gran Capita s/n, 08034 Barcelona, Spain Department of Strength of Materials and Structural Engineering, ETSECCPB, Universitat Politcnica de Catalunya, Barcelona Tech (UPC), Campus Nord, Building C1, c/Jordi Girona 1-3, 08034 Barcelona, Spain
| |
Collapse
|
33
|
Homogenized constrained mixture models for anisotropic volumetric growth and remodeling. Biomech Model Mechanobiol 2016; 16:889-906. [DOI: 10.1007/s10237-016-0859-1] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2016] [Accepted: 11/18/2016] [Indexed: 10/20/2022]
|
34
|
Hiorns JE, Jensen OE, Brook BS. Static and dynamic stress heterogeneity in a multiscale model of the asthmatic airway wall. J Appl Physiol (1985) 2016; 121:233-47. [PMID: 27197860 DOI: 10.1152/japplphysiol.00715.2015] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2015] [Accepted: 05/15/2016] [Indexed: 12/13/2022] Open
Abstract
Airway hyperresponsiveness (AHR) is a key characteristic of asthma that remains poorly understood. Tidal breathing and deep inspiration ordinarily cause rapid relaxation of airway smooth muscle (ASM) (as demonstrated via application of length fluctuations to tissue strips) and are therefore implicated in modulation of AHR, but in some cases (such as application of transmural pressure oscillations to isolated intact airways) this mechanism fails. Here we use a multiscale biomechanical model for intact airways that incorporates strain stiffening due to collagen recruitment and dynamic force generation by ASM cells to show that the geometry of the airway, together with interplay between dynamic active and passive forces, gives rise to large stress and compliance heterogeneities across the airway wall that are absent in tissue strips. We show further that these stress heterogeneities result in auxotonic loading conditions that are currently not replicated in tissue-strip experiments; stresses in the strip are similar to hoop stress only at the outer airway wall and are under- or overestimates of stresses at the lumen. Taken together these results suggest that a previously underappreciated factor, stress heterogeneities within the airway wall and consequent ASM cellular response to this micromechanical environment, could contribute to AHR and should be explored further both theoretically and experimentally.
Collapse
Affiliation(s)
- J E Hiorns
- School of Mathematical Sciences, University of Nottingham, Nottingham, United Kingdom; and
| | - O E Jensen
- School of Mathematics, University of Manchester, Manchester, United Kingdom
| | - B S Brook
- School of Mathematical Sciences, University of Nottingham, Nottingham, United Kingdom; and
| |
Collapse
|
35
|
Seyedsalehi S, Zhang L, Choi J, Baek S. Prior Distributions of Material Parameters for Bayesian Calibration of Growth and Remodeling Computational Model of Abdominal Aortic Wall. J Biomech Eng 2016. [PMID: 26201289 DOI: 10.1115/1.4031116] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
For the accurate prediction of the vascular disease progression, there is a crucial need for developing a systematic tool aimed toward patient-specific modeling. Considering the interpatient variations, a prior distribution of model parameters has a strong influence on computational results for arterial mechanics. One crucial step toward patient-specific computational modeling is to identify parameters of prior distributions that reflect existing knowledge. In this paper, we present a new systematic method to estimate the prior distribution for the parameters of a constrained mixture model using previous biaxial tests of healthy abdominal aortas (AAs). We investigate the correlation between the estimated parameters for each constituent and the patient's age and gender; however, the results indicate that the parameters are correlated with age only. The parameters are classified into two groups: Group-I in which the parameters ce, ck1, ck2, cm2,Ghc, and ϕe are correlated with age, and Group-II in which the parameters cm1, Ghm, G1e, G2e, and α are not correlated with age. For the parameters in Group-I, we used regression associated with age via linear or inverse relations, in which their prior distributions provide conditional distributions with confidence intervals. For Group-II, the parameter estimated values were subjected to multiple transformations and chosen if the transformed data had a better fit to the normal distribution than the original. This information improves the prior distribution of a subject-specific model by specifying parameters that are correlated with age and their transformed distributions. Therefore, this study is a necessary first step in our group's approach toward a Bayesian calibration of an aortic model. The results from this study will be used as the prior information necessary for the initialization of Bayesian calibration of a computational model for future applications.
Collapse
|
36
|
Aparício P, Thompson MS, Watton PN. A novel chemo-mechano-biological model of arterial tissue growth and remodelling. J Biomech 2016; 49:2321-30. [PMID: 27184922 DOI: 10.1016/j.jbiomech.2016.04.037] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2016] [Accepted: 04/18/2016] [Indexed: 02/08/2023]
Abstract
Arterial growth and remodelling (G&R) is mediated by vascular cells in response to their chemical and mechanical environment. To date, mechanical and biochemical stimuli tend to be modelled separately, however this ignores their complex interplay. Here, we present a novel mathematical model of arterial chemo-mechano-biology. We illustrate its application to the development of an inflammatory aneurysm in the descending human aorta. The arterial wall is modelled as a bilayer cylindrical non-linear elastic membrane, which is internally pressurised and axially stretched. The medial degradation that accompanies aneurysm development is driven by an inflammatory response. Collagen remodelling is simulated by adaption of the natural reference configuration of constituents; growth is simulated by changes in normalised mass-densities. We account for the distribution of attachment stretches that collagen fibres are configured to the matrix and, innovatively, allow this distribution to remodel. This enables the changing functional role of the adventitia to be simulated. Fibroblast-mediated collagen growth is represented using a biochemical pathway model: a system of coupled non-linear ODEs governs the evolution of fibroblast properties and levels of key biomolecules under the regulation of Transforming Growth Factor (TGF)-β, a key promoter of matrix deposition. Given physiologically realistic targets, different modes of aneurysm development can be captured, while the predicted evolution of biochemical variables is qualitatively consistent with trends observed experimentally. Interestingly, we observe that increasing the levels of collagen-promoting TGF-β results in arrest of aneurysm growth, which seems to be consistent with experimental evidence. We conclude that this novel Chemo-Mechano-Biological (CMB) mathematical model has the potential to provide new mechanobiological insight into vascular disease progression and therapy.
Collapse
Affiliation(s)
- Pedro Aparício
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, UK.
| | - Mark S Thompson
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, UK.
| | - Paul N Watton
- Department of Computer science, University of Sheffield, Sheffield, UK; INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, UK.
| |
Collapse
|
37
|
Volokh KY, Aboudi J. Aneurysm strength can decrease under calcification. J Mech Behav Biomed Mater 2015; 57:164-74. [PMID: 26717251 DOI: 10.1016/j.jmbbm.2015.11.012] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2015] [Revised: 11/04/2015] [Accepted: 11/11/2015] [Indexed: 12/18/2022]
Abstract
Aneurysms are abnormal dilatations of vessels in the vascular system that are prone to rupture. Prediction of the aneurysm rupture is a challenging and unsolved problem. Various factors can lead to the aneurysm rupture and, in the present study, we examine the effect of calcification on the aneurysm strength by using micromechanical modeling. The calcified tissue is considered as a composite material in which hard calcium particles are embedded in a hyperelastic soft matrix. Three experimentally calibrated constitutive models incorporating a failure description are used for the matrix representation. Two constitutive models describe the aneurysmal arterial wall and the third one - the intraluminal thrombus. The stiffness and strength of the calcified tissue are simulated in uniaxial tension under the varying amount of calcification, i.e. the relative volume of the hard inclusion within the periodic unit cell. In addition, the triaxiality of the stress state, which can be a trigger for the cavitation instability, is tracked. Results of the micromechanical simulation show an increase of the stiffness and a possible decrease of the strength of the calcified tissue as compared to the non-calcified one. The obtained results suggest that calcification (i.e. the presence of hard particles) can significantly affect the stiffness and strength of soft tissue. The development of refined experimental techniques that will allow for the accurate quantitative assessment of calcification is desirable.
Collapse
Affiliation(s)
| | - Jacob Aboudi
- Faculty of Engineering, Tel Aviv University, Israel.
| |
Collapse
|
38
|
Farsad M, Zeinali-Davarani S, Choi J, Baek S. Computational Growth and Remodeling of Abdominal Aortic Aneurysms Constrained by the Spine. J Biomech Eng 2015; 137:2397298. [PMID: 26158885 PMCID: PMC4574855 DOI: 10.1115/1.4031019] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2014] [Revised: 06/27/2015] [Indexed: 01/01/2023]
Abstract
Abdominal aortic aneurysms (AAAs) evolve over time, and the vertebral column, which acts as an external barrier, affects their biomechanical properties. Mechanical interaction between AAAs and the spine is believed to alter the geometry, wall stress distribution, and blood flow, although the degree of this interaction may depend on AAAs specific configurations. In this study, we use a growth and remodeling (G&R) model, which is able to trace alterations of the geometry, thus allowing us to computationally investigate the effect of the spine for progression of the AAA. Medical image-based geometry of an aorta is constructed along with the spine surface, which is incorporated into the computational model as a cloud of points. The G&R simulation is initiated by local elastin degradation with different spatial distributions. The AAA-spine interaction is accounted for using a penalty method when the AAA surface meets the spine surface. The simulation results show that, while the radial growth of the AAA wall is prevented on the posterior side due to the spine acting as a constraint, the AAA expands faster on the anterior side, leading to higher curvature and asymmetry in the AAA configuration compared to the simulation excluding the spine. Accordingly, the AAA wall stress increases on the lateral, posterolateral, and the shoulder regions of the anterior side due to the AAA-spine contact. In addition, more collagen is deposited on the regions with a maximum diameter. We show that an image-based computational G&R model not only enhances the prediction of the geometry, wall stress, and strength distributions of AAAs but also provides a framework to account for the interactions between an enlarging AAA and the spine for a better rupture potential assessment and management of AAA patients.
Collapse
Affiliation(s)
- Mehdi Farsad
- Department of Mechanical Engineering,
Michigan State University,
East Lansing, MI 48824
e-mail:
| | | | - Jongeun Choi
- Associate Professor
Department of Mechanical Engineering,
Michigan State University,
East Lansing, MI 48824
- Department of Electrical and
Computer Engineering,
Michigan State University,
East Lansing, MI 48824
e-mail:
| | - Seungik Baek
- Associate Professor
Department of Mechanical Engineering,
Michigan State University,
East Lansing, MI 48824
e-mail:
| |
Collapse
|
39
|
Wu J, Shadden SC. Coupled Simulation of Hemodynamics and Vascular Growth and Remodeling in a Subject-Specific Geometry. Ann Biomed Eng 2015; 43:1543-54. [PMID: 25731141 PMCID: PMC4497867 DOI: 10.1007/s10439-015-1287-6] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2014] [Accepted: 02/19/2015] [Indexed: 01/20/2023]
Abstract
A computational framework to couple vascular growth and remodeling (G&R) with blood flow simulation in a 3D patient-specific geometry is presented. Hyperelastic and anisotropic properties are considered for the vessel wall material and a constrained mixture model is used to represent multiple constituents in the vessel wall, which was modeled as a membrane. The coupled simulation is divided into two time scales-a longer time scale for G&R and a shorter time scale for fluid dynamics simulation. G&R is simulated to evolve the boundary of the fluid domain, and fluid simulation is in turn used to generate wall shear stress and transmural pressure data that regulates G&R. To minimize required computation cost, the fluid dynamics are only simulated when G&R causes significant vascular geometric change. For demonstration, this coupled model was used to study the influence of stress-mediated growth parameters, and blood flow mechanics, on the behavior of the vascular tissue growth in a model of the infrarenal aorta derived from medical image data.
Collapse
Affiliation(s)
- Jiacheng Wu
- Mechanical Engineering, University of California, 5126 Etcheverry Hall, Berkeley, CA, 94720-1740, USA
| | | |
Collapse
|
40
|
NABAEI MALIKEH, FATOURAEE NASSER. A 3D MODEL FOR MURAL-CELL-MEDIATED DESTRUCTIVE REMODELING DURING EARLY DEVELOPMENT OF A CEREBRAL ANEURYSM. J MECH MED BIOL 2015. [DOI: 10.1142/s0219519415500347] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Development of a diagnostic tool for predicting the behavior of cerebral aneurysms was the inspiration of many research groups in recent years. In the present study a fluid–solid-growth (FSG) model for the early development of a cerebral aneurysm was presented in a 3D model of the internal carotid artery (ICA). This model is the result of two parallel mechanisms: first, defining arterial wall as a living tissue with the ability of degradation, growth and remodeling and second, full coupling of the wall and the blood flow. Taking into account the shear dependent nature of elastin degradation and mural-cell-mediated destructive activities, here, the degradation process has been linked to high effective stress of the vascular wall. The evolving properties of the elastinous and collagenous constituents have been predicted during the early development of the aneurysm and the code is applicable to more complicated aneurismal growth models.
Collapse
Affiliation(s)
- MALIKEH NABAEI
- Biological Fluid Mechanics Research Laboratory, Faculty of Biomedical Engineering, Amirkabir University of Technology, (Tehran Polytechnic), Tehran, P. O. Box 15875-3413, I.R. Iran, Postal Code 15914, Iran
| | - NASSER FATOURAEE
- Biological Fluid Mechanics Research Laboratory, Faculty of Biomedical Engineering, Amirkabir University of Technology, (Tehran Polytechnic), Tehran, P. O. Box 15875-3413, I.R. Iran, Postal Code 15914, Iran
| |
Collapse
|
41
|
Grytsan A, Watton PN, Holzapfel GA. A Thick-Walled Fluid–Solid-Growth Model of Abdominal Aortic Aneurysm Evolution: Application to a Patient-Specific Geometry. J Biomech Eng 2015; 137:2020812. [DOI: 10.1115/1.4029279] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2014] [Indexed: 11/08/2022]
Abstract
We propose a novel thick-walled fluid–solid-growth (FSG) computational framework for modeling vascular disease evolution. The arterial wall is modeled as a thick-walled nonlinearly elastic cylindrical tube consisting of two layers corresponding to the media-intima and adventitia, where each layer is treated as a fiber-reinforced material with the fibers corresponding to the collagenous component. Blood is modeled as a Newtonian fluid with constant density and viscosity; no slip and no-flux conditions are applied at the arterial wall. Disease progression is simulated by growth and remodeling (G&R) of the load bearing constituents of the wall. Adaptions of the natural reference configurations and mass densities of constituents are driven by deviations of mechanical stimuli from homeostatic levels. We apply the novel framework to model abdominal aortic aneurysm (AAA) evolution. Elastin degradation is initially prescribed to create a perturbation to the geometry which results in a local decrease in wall shear stress (WSS). Subsequent degradation of elastin is driven by low WSS and an aneurysm evolves as the elastin degrades and the collagen adapts. The influence of transmural G&R of constituents on the aneurysm development is analyzed. We observe that elastin and collagen strains evolve to be transmurally heterogeneous and this may facilitate the development of tortuosity. This multiphysics framework provides the basis for exploring the influence of transmural metabolic activity on the progression of vascular disease.
Collapse
Affiliation(s)
- Andrii Grytsan
- Department of Solid Mechanics, Royal Institute of Technology (KTH), Teknikringen 8d, Stockholm 10044, Sweden
| | - Paul N. Watton
- Department of Computer Science, University of Sheffield, Sheffield, UK
- INSIGNEO Institute of In Silico Medicine, University of Sheffield, Sheffield, UK
| | - Gerhard A. Holzapfel
- Institute of Biomechanics, Graz University of Technology, Kronesgasse 5-I, Graz 8010, Austria e-mail:
| |
Collapse
|
42
|
Volokh KY. Cavitation instability as a trigger of aneurysm rupture. Biomech Model Mechanobiol 2015; 14:1071-9. [PMID: 25637515 DOI: 10.1007/s10237-015-0655-3] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2014] [Accepted: 01/22/2015] [Indexed: 11/26/2022]
Abstract
Aneurysm formation and growth is accompanied by microstructural alterations in the arterial wall. Particularly, the loss of elastin may lead to tissue disintegration and appearance of voids or cavities at the micron scale. Unstable growth and coalescence of voids may be a predecessor and trigger for the onset of macroscopic cracks. In the present work, we analyze the instability of membrane (2D) and bulk (3D) voids under hydrostatic tension by using two experimentally calibrated constitutive models of abdominal aortic aneurysm enhanced with energy limiters. The limiters provide the saturation value for the strain energy, which indicates the maximum energy that can be stored and dissipated by an infinitesimal material volume. We find that the unstable growth of voids can start when the critical stress is considerably less than the aneurysm strength. Moreover, this critical stress may even approach the arterial wall stress in the physiological range. This finding suggests that cavitation instability can be a rational indicator of the aneurysm rupture.
Collapse
Affiliation(s)
- K Y Volokh
- Faculty of Civil and Environmental Engineering, Technion - I.I.T., Haifa, Israel,
| |
Collapse
|
43
|
Hald ES, Steucke KE, Reeves JA, Win Z, Alford PW. Long-term vascular contractility assay using genipin-modified muscular thin films. Biofabrication 2014; 6:045005. [PMID: 25245868 DOI: 10.1088/1758-5082/6/4/045005] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Vascular disease is a leading cause of death globally and typically manifests chronically due to long-term maladaptive arterial growth and remodeling. To date, there is no in vitro technique for studying vascular function over relevant disease time courses that both mimics in vivo-like tissue structure and provides a simple readout of tissue stress. We aimed to extend tissue viability in our muscular thin film contractility assay by modifying the polydimethylsiloxane (PDMS) substrate with micropatterned genipin, allowing extracellular matrix turnover without cell loss. To achieve this, we developed a microfluidic delivery system to pattern genipin and extracellular matrix proteins on PDMS prior to cell seeding. Tissues constructed using this method showed improved viability and maintenance of in vivo-like lamellar structure. Functional contractility of tissues fabricated on genipin-modified substrates remained consistent throughout two weeks in culture. These results suggest that muscular thin films with genipin-modified PDMS substrates are a viable method for conducting functional studies of arterial growth and remodeling in vascular diseases.
Collapse
Affiliation(s)
- Eric S Hald
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455, USA
| | | | | | | | | |
Collapse
|
44
|
Nabaei M, Fatouraee N. Microstructural modelling of cerebral aneurysm evolution through effective stress mediated destructive remodelling. J Theor Biol 2014; 354:60-71. [PMID: 24657628 DOI: 10.1016/j.jtbi.2014.03.020] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2013] [Revised: 03/10/2014] [Accepted: 03/11/2014] [Indexed: 11/25/2022]
Abstract
Recently, researchers have shown an increased interest in the biomechanical modelling of cerebral aneurysm development. In the present study a fluid-solid-growth model for the formation of a fusiform aneurysm has been presented in an axi-symmetric geometry of the internal carotid artery. This model is the result of two parallel mechanisms: first, defining arterial wall as a living tissue with the ability of degradation, growth and remodelling and second, full coupling of the wall and the blood flow. Here for the first time the degradation of elastin has been defined as a function of vascular wall effective stress to take into account the shear dependent nature of degradation and the mural-cell-mediated destructive activities. The model has been stabilized in size and mechanical properties and is consistent with other computational or clinical studies. Furthermore, the evolving microstructural properties of the wall during the evolution process have been predicted.
Collapse
|
45
|
Wilson JS, Humphrey JD. Evolving anisotropy and degree of elastolytic insult in abdominal aortic aneurysms: potential clinical relevance? J Biomech 2014; 47:2995-3002. [PMID: 25086482 DOI: 10.1016/j.jbiomech.2014.07.003] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Revised: 05/05/2014] [Accepted: 07/02/2014] [Indexed: 12/16/2022]
Abstract
Accurately estimating patient-specific rupture risk remains a primary challenge in timing interventions for abdominal aortic aneurysms (AAAs). By re-analyzing published biaxial mechanical testing data from surgically repaired human AAAs, material anisotropy emerged as a potentially important determinant of patient-specific lesion progression. That is, based on a new classification scheme, we discovered that anisotropic aneurysmal specimens correlated with increased patient age at surgery when compared with more isotropic specimens (79.7 vs. 70.9 years, p<0.002), despite no significant difference in maximum diameter. Furthermore, using an idealized axisymmetric, finite-element growth and remodeling model of AAA progression, we found that both the initial axial extent of elastin loss and ongoing damage to elastin in the shoulder region of the AAA directly affected the degree of anisotropy as the lesion evolved, with more extensive insults increasing the anisotropy. This effect appeared to be mediated by alterations in axial loading and subsequent differences in orientation of deposited collagen fibers. While the observed increased age before surgical intervention may suggest a potential benefit of anisotropic remodeling, future biaxial tests coupled with pre-surgical data on expansion rates and detailed theoretical analyses of the biostability of a lesion as a function of anisotropy will be required to verify its clinical relevance to patient-specific rupture risk.
Collapse
Affiliation(s)
- John S Wilson
- Department of Biomedical Engineering, Yale University, New Haven, CT 06520, USA
| | - J D Humphrey
- Department of Biomedical Engineering, Yale University, New Haven, CT 06520, USA; Vascular Biology and Therapeutics Program, Yale School of Medicine, New Haven, CT 06520, USA.
| |
Collapse
|
46
|
Lanir Y. Mechanistic micro-structural theory of soft tissues growth and remodeling: tissues with unidirectional fibers. Biomech Model Mechanobiol 2014; 14:245-66. [DOI: 10.1007/s10237-014-0600-x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2013] [Accepted: 05/23/2014] [Indexed: 10/25/2022]
|
47
|
Aparício P, Mandaltsi A, Boamah J, Chen H, Selimovic A, Bratby M, Uberoi R, Ventikos Y, Watton PN. Modelling the influence of endothelial heterogeneity on the progression of arterial disease: application to abdominal aortic aneurysm evolution. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2014; 30:563-586. [PMID: 24424963 DOI: 10.1002/cnm.2620] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/27/2013] [Revised: 09/19/2013] [Accepted: 11/18/2013] [Indexed: 06/03/2023]
Abstract
We sophisticate a fluid-solid growth computational framework for modelling aneurysm evolution. A realistic structural model of the arterial wall is integrated into a patient-specific geometry of the vasculature. This enables physiologically representative distributions of haemodynamic stimuli, obtained from a rigid-wall computational fluid dynamics analysis, to be linked to growth and remodelling algorithms. Additionally, a quasistatic structural analysis quantifies the cyclic deformation of the arterial wall so that collagen growth and remodelling can be explicitly linked to the cyclic deformation of vascular cells. To simulate aneurysm evolution, degradation of elastin is driven by reductions in wall shear stress (WSS) below homeostatic thresholds. Given that the endothelium exhibits spatial and temporal heterogeneity, we propose a novel approach to define the homeostatic WSS thresholds: We allow them to be spatially and temporally heterogeneous. We illustrate the application of this novel fluid-solid growth framework to model abdominal aortic aneurysm (AAA) evolution and to examine how the influence of the definition of the WSS homeostatic threshold influences AAA progression. We conclude that improved understanding and modelling of the endothelial heterogeneity is important for modelling aneurysm evolution and, more generally, other vascular diseases where haemodynamic stimuli play an important role.
Collapse
Affiliation(s)
- P Aparício
- Systems Biology Doctoral Training Centre, University of Oxford, Oxford, UK
| | | | | | | | | | | | | | | | | |
Collapse
|
48
|
Balakhovsky K, Jabareen M, Volokh KY. Modeling rupture of growing aneurysms. J Biomech 2013; 47:653-8. [PMID: 24359675 DOI: 10.1016/j.jbiomech.2013.11.049] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2013] [Revised: 11/19/2013] [Accepted: 11/27/2013] [Indexed: 11/27/2022]
Abstract
Growth and rupture of aneurysms are driven by micro-structural alterations of the arterial wall yet precise mechanisms underlying the process remain to be uncovered. In the present work we examine a scenario when the aneurysm evolution is dominated by turnover of collagen fibers. In the latter case it is natural to hypothesize that rupture of individual fibers (or their bonds) causes the overall aneurysm rupture. We examine this hypothesis in computer simulations of growing aneurysms in which constitutive equations describe both collagen evolution and failure. Failure is enforced in constitutive equations by limiting strain energy that can be accumulated in a fiber. Within the proposed theoretical framework we find a range of parameters that lead to the aneurysm rupture. We conclude in a qualitative agreement with clinical observations that some aneurysms will rupture while others will not.
Collapse
Affiliation(s)
- K Balakhovsky
- Faculty of Civil and Environmental Engineering, Technion - I.I.T., Israel
| | - M Jabareen
- Faculty of Civil and Environmental Engineering, Technion - I.I.T., Israel
| | - K Y Volokh
- Faculty of Civil and Environmental Engineering, Technion - I.I.T., Israel; Department of Structural Engineering, Ben-Gurion University of the Negev, Israel.
| |
Collapse
|
49
|
Wilson JS, Virag L, Di Achille P, Karsaj I, Humphrey JD. Biochemomechanics of intraluminal thrombus in abdominal aortic aneurysms. J Biomech Eng 2013; 135:021011. [PMID: 23445056 DOI: 10.1115/1.4023437] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Most computational models of abdominal aortic aneurysms address either the hemodynamics within the lesion or the mechanics of the wall. More recently, however, some models have appropriately begun to account for the evolving mechanics of the wall in response to the changing hemodynamic loads. Collectively, this large body of work has provided tremendous insight into this life-threatening condition and has provided important guidance for current research. Nevertheless, there has yet to be a comprehensive model that addresses the mechanobiology, biochemistry, and biomechanics of thrombus-laden abdominal aortic aneurysms. That is, there is a pressing need to include effects of the hemodynamics on both the development of the nearly ubiquitous intraluminal thrombus and the evolving mechanics of the wall, which depends in part on biochemical effects of the adjacent thrombus. Indeed, there is increasing evidence that intraluminal thrombus in abdominal aortic aneurysms is biologically active and should not be treated as homogeneous inert material. In this review paper, we bring together diverse findings from the literature to encourage next generation models that account for the biochemomechanics of growth and remodeling in patient-specific, thrombus-laden abdominal aortic aneurysms.
Collapse
Affiliation(s)
- J S Wilson
- Department of Biomedical Engineering, Yale University, New Haven, CT 06520, USA
| | | | | | | | | |
Collapse
|
50
|
Martufi G, Christian Gasser T. Review: the role of biomechanical modeling in the rupture risk assessment for abdominal aortic aneurysms. J Biomech Eng 2013; 135:021010. [PMID: 23445055 DOI: 10.1115/1.4023254] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
AAA disease is a serious condition and a multidisciplinary approach including biomechanics is needed to better understand and more effectively treat this disease. A rupture risk assessment is central to the management of AAA patients, and biomechanical simulation is a powerful tool to assist clinical decisions. Central to such a simulation approach is a need for robust and physiologically relevant models. Vascular tissue senses and responds actively to changes in its mechanical environment, a crucial tissue property that might also improve the biomechanical AAA rupture risk assessment. Specifically, constitutive modeling should not only focus on the (passive) interaction of structural components within the vascular wall, but also how cells dynamically maintain such a structure. In this article, after specifying the objectives of an AAA rupture risk assessment, the histology and mechanical properties of AAA tissue, with emphasis on the wall, are reviewed. Then a histomechanical constitutive description of the AAA wall is introduced that specifically accounts for collagen turnover. A test case simulation clearly emphasizes the need for constitutive descriptions that remodels with respect to the mechanical loading state. Finally, remarks regarding modeling of realistic clinical problems and possible future trends conclude the article.
Collapse
Affiliation(s)
- Giampaolo Martufi
- Department of Solid Mechanics, School of Engineering Sciences, Royal Institute of Technology (KTH), Osquars Backe 1, SE-100 44 Stockholm, Sweden.
| | | |
Collapse
|