A Multiple Step Active Stiffness Integration Scheme to Couple a Stochastic Cross-Bridge Model and Continuum Mechanics for Uses in Both Basic Research and Clinical Applications of Heart Simulation.
A Multiple Step Active Stiffness Integration Scheme to Couple a Stochastic Cross-Bridge Model and Continuum Mechanics for Uses in Both Basic Research and Clinical Applications of Heart Simulation.
复制标题
多个步骤的主动刚度整合方案，将随机跨桥模型和连续力学融入心脏模拟的基础研究和临床应用中。
DOI：
10.3389/fphys.2021.712816
复制
发表时间：
2021
影响因子：
4
通讯作者：
Washio T
中科院分区：
文献类型：
--
作者：
Yoneda K;Okada JI;Watanabe M;Sugiura S;Hisada T;Washio T
在跳动的心脏的多尺度模拟中，收缩蛋白的快速随机构象变化和确定性的宏观结果（如心室压力和容积）之间的时间尺度上的非常大的差异阻碍了两个尺度的有效耦合算法的实施。此外，在连续介质力学中还没有很好地建立考虑由马达蛋白的跨桥活动引起的肌肉刚度的动态变化。为了克服这些问题，我们提出了一个多时间步长的计划称为多步主动刚度积分计划（MusAsi）的Monte Carlo（MC）多个步骤和隐式有限元（FE）时间积分步骤的耦合。该方法的重点是主动张力刚度矩阵，其中的主动张力导数关于当前的位移在有限元模型中被正确地集成到总刚度矩阵，以避免不稳定。对MC模型中使用的样本数量和时间步长组合的敏感性分析证实了MusAsi的准确性和鲁棒性，我们得出结论，在每个有限元中使用几百个马达蛋白质的1.25 ms FE时间步长和0.005 ms MC多个步长的组合在准确性和计算时间之间的权衡中是适当的。此外，对于由45，000个四面体单元组成的双心室FE模型，使用320个核心的传统并行计算机系统可以在1.5小时内计算一次心跳。这些结果支持MusAsi用于分子机制和心输出量之间关系的基础研究和围手术期预测的临床应用的实用性。
In a multiscale simulation of a beating heart, the very large difference in the time scales between rapid stochastic conformational changes of contractile proteins and deterministic macroscopic outcomes, such as the ventricular pressure and volume, have hampered the implementation of an efficient coupling algorithm for the two scales. Furthermore, the consideration of dynamic changes of muscle stiffness caused by the cross-bridge activity of motor proteins have not been well established in continuum mechanics. To overcome these issues, we propose a multiple time step scheme called the multiple step active stiffness integration scheme (MusAsi) for the coupling of Monte Carlo (MC) multiple steps and an implicit finite element (FE) time integration step. The method focuses on the active tension stiffness matrix, where the active tension derivatives concerning the current displacements in the FE model are correctly integrated into the total stiffness matrix to avoid instability. A sensitivity analysis of the number of samples used in the MC model and the combination of time step sizes confirmed the accuracy and robustness of MusAsi, and we concluded that the combination of a 1.25 ms FE time step and 0.005 ms MC multiple steps using a few hundred motor proteins in each finite element was appropriate in the tradeoff between accuracy and computational time. Furthermore, for a biventricular FE model consisting of 45,000 tetrahedral elements, one heartbeat could be computed within 1.5 h using 320 cores of a conventional parallel computer system. These results support the practicality of MusAsi for uses in both the basic research of the relationship between molecular mechanisms and cardiac outputs, and clinical applications of perioperative prediction.
登录
查看更多内容
影响因子：
3.1
作者：
Dabiri, Yaghoub;Van der Velden, Alex;Guccione, Julius M.
通讯作者：
Guccione, Julius M.
DOI：
10.1615/intjmultcompeng.2011002360
发表时间：
2012-01-01
影响因子：
1.4
作者：
Chapelle, D.;Le Tallec, P.;Sorine, M.
通讯作者：
Sorine, M.
DOI：
10.1016/j.cma.2020.113506
发表时间：
2021-01-01
影响因子：
7.2
作者：
Regazzoni, F.;Quarteroni, A.
通讯作者：
Quarteroni, A.
影响因子：
3.4
作者：
Rice, JJ;Stolovitzky, G;de Tombe, PP
通讯作者：
de Tombe, PP
影响因子：
1
作者：
Azzolin, Luca;Dede, Luca;Quarteroni, Alfio
通讯作者：
Quarteroni, Alfio