Bayesian inverse problems for soft tissue mechanics

软组织力学的贝叶斯反问题

• 批准号: 2596737
• 金额: --
• 依托单位: University of Manchester
• 依托单位国家: 英国
• 项目类别: Studentship
• 财政年份: 2021
• 资助国家: 英国
• 起止时间: 2021 至 无数据
• 项目状态: 未结题
• 来源: https://gtr.ukri.org/projects?ref=studentship-2596737
• 关键词: Bayesian; inverse; problems; soft; tissue

Over the last few decades, a large body of literature has focussed on mathematically modelling the mechanical behaviour of biological soft tissues within a continuum mechanics framework. Particular attention has focused on deriving constitutive equations within the context of large strain (nonlinear) elasticity and viscoelasticity. Many of the most widely used models take a phenomenological approach (for example, by proposing a strain energy which is an exponential function of the strain invariants, with free parameters that can be fitted to experimental data); however, this approach cannot be used to predict the effects of microstructural changes on macroscale mechanics. An alternative approach explicitly models the microstructure of the soft tissue, with measurable parameters to describe the geometry of the collagen fibre network that makes up the tissue, along with simple (linear) constitutive equations for the collagen fibres themselves. These linear constitutive models fit experimental data on individual collagen fibres well; however, due to the difficulty in precisely measuring the geometries of, and forces acting on, fibres that can be as small as tens of nanometres in diameter, the reported values of the constitutive parameters (e.g. the collagen Young's modulus) vary by orders of magnitude. This uncertainty in the values of these parameters makes it difficult to make quantitative predictions with deterministic, microstructural models; therefore, an alternative approach is needed which accounts for the uncertainty. In this project, we will invoke the Bayesian framework, which, by incorporating the models of soft tissues with prior beliefs and macroscale experimental data, will give us a posterior probability distribution of the parameters conditioned on the observations. These distributions contain not only information about the likely values of the parameters, but also allow us to quantify and assess the uncertainty inherent in the estimates that arise from them. In practice this can be achieved by implementing Markov chain Monte Carlo (MCMC) methods, which through construction of an ergodic Markov chain with stationary distribution equal to the posterior distribution, allows us to sample from the distribution in order to characterise it. This might be done through the implementation of existing MCMC methods, or through the design of new methodologies which are able to efficiently target the types of distributions which arise from these inverse problems.

在过去的几十年里，大量的文献集中于在连续介质力学框架内对生物软组织的力学行为进行数学建模。特别关注在大应变(非线性)弹性和粘弹性的背景下推导本构方程。许多最广泛使用的模型采用现象学方法(例如，通过提出应变不变量的指数函数的应变能，其自由参数可以与实验数据拟合)；然而，这种方法不能用于预测微观结构变化对宏观尺度力学的影响。另一种方法明确地模拟软组织的微结构，用可测量的参数来描述构成组织的胶原纤维网络的几何形状，以及胶原纤维本身的简单(线性)本构方程。这些线性本构模型很好地符合单个胶原纤维的实验数据；然而，由于难以精确测量直径可小至数十纳米的纤维的几何形状和作用力，已报道的本构参数(例如，胶原杨氏模数)的值按数量级变化。这些参数值的这种不确定性使得用确定性的微观结构模型进行定量预测变得困难；因此，需要一种解释这种不确定性的替代方法。在这个项目中，我们将引用贝叶斯框架，通过结合具有先验信念的软组织模型和宏观实验数据，我们将给出基于观测条件的参数的后验概率分布。这些分布不仅包含有关参数的可能值的信息，而且还允许我们量化和评估由它们引起的估计中固有的不确定性。在实践中，这可以通过实现马尔可夫链蒙特卡罗(MCMC)方法来实现，该方法通过构造平稳分布等于后验分布的遍历马尔可夫链，允许我们从分布中采样以对其进行表征。这可以通过实施现有的MCMC方法来实现，或者通过设计能够有效地针对由这些逆问题产生的分布类型的新方法来实现。