×

用随机参数量化超弹性软组织模型中的不确定性。 (英语) Zbl 1460.74065号

总结:我们提出了一种简单的开源半侵入式计算方法,通过软组织的超弹性模型传播不确定性。该方法比标准蒙特卡罗方法快两个数量级。可以通过调整几行(FEniCS)代码来更改感兴趣的材料模型。该方法能够(1)为用户提供关于实际感兴趣的数量的统计置信区间,例如肿瘤或器官中目标部位的位移;(2) 量化器官对材料模型的相关参数的响应的灵敏度。我们采用这种方法来确定大脑中目标运动的置信区间。我们还表明,对于所考虑的边界条件,与三个最有影响的参数相比,Ogden-Holzapfel-like模型的五个参数对目标区位移的影响可以忽略不计。基准问题和所有相关数据作为补充材料提供。

MSC公司:

74升15 生物力学固体力学
74B20型 非线性弹性
74S60系列 应用于固体力学问题的随机和其他概率方法
92立方厘米 生物力学
PDF格式BibTeX公司 XML格式引用
全文: DOI程序 链接

参考文献:

[1] 巴布斯·卡,I。;Baumann,C。;Oden,J.,扩散问题的间断hp有限元方法:一维分析,计算。数学。申请。,37, 9, 103-122 (1999) ·Zbl 0940.65076号
[2] 博尔达斯,S。;Duflot,M.,扩展有限元的导数恢复和后验误差估计,计算。方法应用。机械。工程,196,35,3381-3399(2007)·Zbl 1173.74401号
[3] Bui,H.P。;托马尔,S。;Bordas,S.P.A.,《手术模拟的实时错误控制》(Brieu,M.;Bordass,S.;Leriche,E.;Sherwin,S.,EUROMECH学术讨论会595,生物力学和计算机辅助手术符合医学现实(2017),欧洲力学学会:欧洲力学学会里尔中心学院,法国维伦纽夫dAscq)
[4] Pham博士。;徐,C。;Prince,J.,《医学图像分割的当前方法》,《生物识别年鉴》。工程,2,1,315-337(2000)
[5] 克莱因,S。;范德海德,美国。;嘴唇,I。;van Vulpen,M。;凝视,M。;Pluim,J.,使用局部互信息通过图谱匹配在三维MR图像中自动分割前列腺,医学物理。,35, 4, 1407-1417 (2008)
[6] Miller,K.,《患者特定应用的计算生物力学》,Ann.Biomed。工程,44,1,1-2(2016)
[7] Viceconti先生。;Bellingeri,L。;克里斯托弗利尼,L。;Toni,A.,《人类股骨不同自动网格生成方法的比较研究》,医学工程物理。,20, 1, 1-10 (1998)
[8] 费伦特,M。;纳巴维,A。;Macq,B。;Jolesz,F.A。;Kikinis,R。;Warfield,S.,使用有限元生物力学模型注册术中大脑的三维Mr图像,IEEE Trans。医学成像,20,12,1384-1397(2001)
[9] Belytschko,T。;帕里米,C。;莫尔斯,N。;北苏库马尔。;Usui,S.,《由隐式曲面定义的实体的结构化扩展有限元方法》,国际期刊数值。方法工程,56,4,609-635(2003)·Zbl 1038.74041号
[10] Nguyen,V。;Rabczuk,T。;博尔达斯,S。;Duflot,M.,《无网格方法:综述和计算机实现方面》,数学。计算。模拟。,79, 3, 763-813 (2008) ·兹比尔1152.74055
[11] 斯特吉奥普洛斯,N。;Young博士。;Rogge,T.,动脉血流的计算机模拟及其在动脉和主动脉狭窄中的应用,J.Biomech。,25, 12, 1477-1488 (1992)
[12] 科廷,S。;Delingette,H。;Ayache,N.,手术模拟软组织的实时弹性变形,IEEE Trans。视觉。计算。组,5,1,62-73(1999)
[13] Holzapfel,G.A。;Ogden,R.W.,《分子、细胞和组织水平的生物力学建模》(2009年),施普林格出版社
[14] Payan,Y.,《计算机辅助手术的软组织生物力学建模》(2012),Springer
[15] Geris,L。;Gomez-Cabrero博士,《生物学的不确定性》(2016),施普林格出版社·Zbl 1325.92003年
[16] Miller,K。;Chinzei,K.,紧张状态下脑组织的力学特性,J.生物医学。,35, 4, 483-490 (2002)
[17] 拉布斯,K.M。;Puttlitz,C.M.,脑白质在双向拉伸和结构-力学关系中的各向异性超弹性本构模型,J.Mech。行为。生物识别。材料。,62, 195-208 (2016)
[18] Chatelin,S。;Constantinesco,A。;Willinger,R.,《脑组织机械测试五十年:从体外到体内研究》,生物流变学,47255-276(2010)
[19] 拉克萨里,K。;沙菲安,M。;Darvish,K.,有限压缩下脑组织的本构模型,J.Biomech。,45, 642-646 (2012)
[20] 安吉拉·米海(Angela Mihai,L.)。;Chin,L。;Janmey,P。;Goriely,A.,《适用于大脑和脂肪组织的超弹性本构模型的比较》,J.R.Soc.Interface,12,110(2015)
[21] de Rooij,R。;Kuhl,E.,《脑组织的本构建模:当前观点》,应用。机械。修订版,68,1(2016)
[22] Goriely,A。;布代,S。;Kuhl,E.,第二章-神经力学:从神经元到大脑,应用力学进展,第48卷,79-139(2015),Elsevier
[23] Joldes,G。;Wittek,A。;Miller,K.,GPU应用于神经外科模拟的实时非线性有限元计算。方法应用。机械。工程,199,49-52,3305-3314(2010)·Zbl 1225.92021号
[24] Caflisch,C.A.,蒙特卡罗和准蒙特卡罗方法,《数值学报》。,7, 1-49 (1998) ·Zbl 0949.65003号
[25] 加尼姆,R。;Spanos,P.,《随机有限元:谱方法》(1991),多佛出版物:多佛出版物,米诺拉,纽约·兹比尔0722.73080
[26] Matthies,G.H.,《随机有限元:随机偏微分方程的计算方法》,J.Appl。数学。机械。,88, 849-873 (2008) ·Zbl 1158.65009号
[27] Giraldi,L。;刘,D。;马蒂斯·H·G。;Nouy,A.,是入侵还是不入侵?参数和随机方程的解——“普通香草”Galerkin案例,SIAM J.Sci。计算。,3622720-2744(2014)·Zbl 1310.65132号
[28] Dopico-González,C。;新,A。;Browne,M.,《股骨特性和植入物设计几何形状的非骨水泥髋关节置换效果的概率有限元分析》,J.Biomech。,43, 3, 512-520 (2010)
[29] 比勒,J。;吉,M.W。;Wall,W.A.,基于贝叶斯多保真度方案的复杂大规模生物力学问题中的有效不确定性量化,生物医学。模型。机械硫醇。,14, 3, 489-513 (2015)
[30] Laz,P.J。;Browne,M.,《矫形生物力学概率分析综述》,《工程医学杂志》,224927-943(2010)
[31] 普兰蒂夫,R。;彼得利克,I。;Courtecuisse,H。;特里维松,R。;Radoux,J。;Cotin,S.,基于Atlas的生物力学模拟边界条件转移,医学图像计算。计算。协助。干预。,17, 33-40 (2014)
[32] Risholm,P。;Fedorov,A。;帕斯利,J。;Tuncali,K。;科马克,R。;Wells,W.,前列腺图像的概率非刚性配准:建模和量化不确定性,IEEE生物医学成像国际研讨会论文集,553-556(2011)
[33] 莫汉,P.S。;Nair,P.B。;Keane,A.J.,随机域上确定性线性椭圆型偏微分方程的随机投影方案,Int.J.Numer。方法。工程,85,874-895(2011)·Zbl 1217.76065号
[34] 林格伦,F。;H街。;Lindström,J.,《高斯场和高斯-马尔可夫随机场之间的明确联系:随机偏微分方程方法》,J.R.Stat.Soc.Ser。B(Stat.Methodol.),第73、4、423-498页(2011年)·Zbl 1274.62360号
[35] 阿德勒·R·J。;Taylor,J.,《随机场与几何》(2007),Springer-Verlag·Zbl 1149.60003号
[36] 曹毅。;侯赛尼,M.Y。;Zhang,T.A.,利用灵敏度衍生物改进采样方法,AIAA J.,42,2,815-822(2004)
[37] 刘,Y。;侯赛尼,M.Y。;Økten,G.,基于灵敏度导数的蒙特卡罗方差减少方法的优化,应用。数字。数学。,72,完成,160-171(2013)·兹比尔1302.65013
[38] P.Hauseux,J.S.Hale,S.P.Bordas,用灵敏度导数驱动的蒙特卡罗方法求解随机burgers方程,无花果10.6084/m9图3561306;P.Hauseux,J.S.Hale,S.P.Bordas,用灵敏度导数驱动的蒙特卡罗方法求解随机burgers方程,无花果10.6084/m9图3561306
[39] Hauseux,P。;Hale,J.S。;Bordas,S.P.,用高级有限元模型的导数加速蒙特卡罗估计,计算。方法应用。机械。工程,318917-936(2017)·Zbl 1439.65006号
[40] 阿尔奈斯,M。;布莱切塔,J。;Hake,J。;A.、J。;Kehlet,B。;Logg,A。;Richardson,C。;Ring,J。;罗杰斯,M。;Wells,G.,The Fenics项目1.5版,Arch。数字。软质。,3, 100 (2015)
[41] Logg,A。;Wells,G.,Dolfin:自动化有限元计算,ACM Trans。数学。软质。,37, 2, 20:1-20:28 (2010) ·Zbl 1364.65254号
[42] 罗伯特·C·P。;Casella,G.,Monte Carlo统计方法,Springer Texts in Statistics(2005),Springer-Verrag New York,Inc.:Springer-Verlag New York,Inc.美国新泽西州塞考克斯
[43] Farrell,体育。;哈姆,D.A。;Funke,S.W。;Rognes,M.E.,《高级瞬态有限元程序伴随的自动推导》,SIAM J.Sci。计算。,35、4、C369-C393(2013)·兹比尔1362.65103
[44] 范伯格,J。;Langtangen,H.,Chaospy:一种用于设计不确定性量化方法的开源工具,J.Compute。科学。,11, 46-57 (2015)
[45] Sobol,I.,非线性数学模型的全局敏感性指数及其蒙特卡罗估计,数学。计算。模拟。,55, 1-3, 271-280 (2001) ·Zbl 1005.65004号
[46] Saltelli,A。;Annoni,P。;阿齐尼,I。;Campolongo,F。;Ratto,M。;Tarantola,S.,基于方差的模型输出敏感性分析。总敏感性指数的设计和估计,计算。物理学。社区。,181, 2, 259-270 (2010) ·Zbl 1219.93116号
[47] Giles,M.,多级蒙特卡罗方法,数字学报。,24, 259-328 (2015) ·兹比尔1316.65010
[48] Sobol’,I.M.,具有额外统一属性的均匀分布序列,USRR Compute。数学。数学。物理。,16, 236-242 (1977) ·Zbl 0391.10033号
[49] 诺瓦克,W。;Tenkleve,S。;Cirpka,O.A.,相互依赖量的线性化互协方差和自协方差矩阵的有效计算,数学。地质。,35, 1, 53-66 (2003) ·Zbl 1302.86028号
[50] A.Logg,K.A.Mardal,G.N.Wells(编辑),《有限元法自动求解微分方程》,计算科学与工程讲义第84卷,施普林格,2012年。doi:10.1007/978-3642-23099-8;A.Logg,K.A.Mardal,G.N.Wells(编辑),《有限元法自动求解微分方程》,计算科学与工程讲义第84卷,施普林格,2012年。doi:10.1007/978-3642-23099-8·兹比尔1247.65105
[51] Holzapfel,G。;Ogden,R.,《被动心肌的本构模型:基于结构的材料表征框架》,Philos。事务处理。R.Soc.伦敦。数学。物理学。工程科学。,367, 1902, 3445-3475 (2009) ·Zbl 1185.74060号
[52] 佩雷斯,F。;Granger,B.,IPython:交互式科学计算系统,Compute。科学。工程,9,3,21-29(2007)
[53] P.Hauseux,J.S.Hale,S.Cotin,S.P.Bordas,用灵敏度导数驱动的蒙特卡罗方法求解随机超弹性方程,无花果10.6084/m9.图片共享4900298;P.Hauseux,J.S.Hale,S.Cotin,S.P.Bordas,用灵敏度导数驱动的蒙特卡罗方法求解随机超弹性方程,无花果10.6084/m9图4900298
[54] 巴伦西亚。;布拉斯,B。;Ortega,J.H.,不同开颅手术和实体模型的脑移位现象建模,J.Appl。数学。,1-21 (2012) ·Zbl 1235.74376号
[55] Pezzuto,S。;Ambrosi,D。;Quarteroni,A.,《心肌力学正交异性激活素模型及其数值近似》,《欧洲医学杂志》。A固体,48,83-96(2014)·Zbl 1406.74502号
[56] 巴拉伊,S。;Abhyankar,S。;亚当斯,M。;Brown,J.等人。;布鲁纳,P。;Buschelman,K。;达尔星。;埃伊霍特,V。;格罗普,W。;考希克,D。;Knepley,M。;柯夫曼,L。;鲁普,K。;B.史密斯。;扎皮尼,S。;张,H。;Zhang,H.,PETSc用户手册,技术报告,ANL-95/11-3.7版(2016),阿贡国家实验室
[57] 法尔古特,R。;Yang,U.,Hypre:一个高性能预处理器库,(Sloot,P.M.a.;Hoekstra,a.G.;Tan,C.J.K.;Dongarra,J.J.,《国际计算科学会议论文集》。《国际计算科学会议论文集》,计算机科学课堂讲稿第2331卷(2002年),《施普林格:施普林格-柏林-海德堡》,第632-641页·Zbl 1056.65046号
[58] 张伯兰,E。;Fortin,A。;Fortin,M.,大变形弹性问题的一些有限元离散性能比较,计算。结构。,88, 1112, 664-673 (2010)
[59] Goriely,A。;Geers,M.G.D。;Holzapfel,G.A。;Jayamohan,J。;Jérusalem,A。;Sivaloganathan,S。;斯奎尔(W.Squier)。;van Dommelen,J。;沃特斯,S。;Kuhl,E.,《大脑力学:观点、挑战和机遇》,《生物技术》。模型。机械硫醇。,14, 5, 931-965 (2015)
[60] Feng,Y。;冈本,R.J。;Namani,R。;杰宁,G.M。;Bayly,P.V.,脑组织机械各向异性的测量及其对白质横向各向同性材料模型的影响,J.Mech。行为。生物识别。材料。,23,补编C,117-132(2013)
[61] Zander,E.K.,随机问题的张量近似方法(2012),Shaker
[62] Sudret,B.,使用多项式混沌展开的全局灵敏度分析。,Reliab公司。工程系统。安全。,93, 7, 964-979 (2008)
[63] Wiener,N.,《齐次混沌》,美国数学杂志。,60, 897-936 (1938)
[64] Malliavin,P.,随机分析,Grundlehren dermathematischen Wissenschaften(1997),Springer:Springer Berlin,纽约·Zbl 0878.60001号
[65] 拉佩尔,H。;洛杉矶Beex。;Bordas,S.P.A.,识别粘弹性参数的贝叶斯推断,机械。时间依赖。材料。,1-32 (2017)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。