Simulating contact using the elastic foundation algorithm in OpenSim
Simulating contact using the elastic foundation algorithm in OpenSim
复制标题
DOI：
10.1016/j.jbiomech.2018.11.025
复制
发表时间：
2019-01-03
影响因子：
2.4
通讯作者：
Baxter, Josh R.
中科院分区：
文献类型：
--
作者：
Hast, Michael W.;Hanson, Brett G.;Baxter, Josh R.
在OpenSim中对关节接触建模还不是很好理解。本研究系统地研究了OpenSim中的弹性地基接触模型的相关变量，通过进行一系列的控制台架实验和伴随的模拟。对四种金属-塑料相互作用进行了建模，包括全膝关节置换术（TKR）模型。在100 - 750 N的循环载荷下记录载荷-位移曲线。将几何形状导入OpenSim，并使用板载弹性基础算法对接触力学进行建模。混合优化算法确定TKR植入物的刚度和耗散系数分别为1.52 x 10（10）N/m和57.7 Ns/m。贴靠力的估计值为设盲实验数据的10.2%（平均均方根误差：76.82 +/- 11.47 N）。在本研究的第二部分中，使用免费提供的eTibia TKR渲染图来测试调谐参数的普遍性。它们还被用来执行材料刚度和网格密度的穿透深度和计算时间方面的敏感性分析。当对具有5000个面的eTibia模型施加1 x 10（10）的刚度时，100 kg载荷导致0.259 mm的穿透。在相同条件下，调整后的模型经历了0.300 mm的穿透。1 x 10（13）和1 x 10（15）N/m之间的材料刚度使计算时间增加12-23倍。这项研究提供了关于OpenSim EF算法的使用非常需要的清晰度。它还展示了OpenSim的实用性，可变形的材料和复杂的几何形状建模，这种方法可以适应自然和手术修改关节作出合理的估计。(C)2018爱思唯尔有限公司版权所有
Modeling joint contact in OpenSim is not well understood. This study systematically investigated the variables associated with the elastic foundation contact model within OpenSim by performing a series of controlled benchtop experiment and concomitant simulations. Four metal-on-plastic interactions were modeled, including a model of a total knee replacement (TKR). Load-displacement curves were recorded during cyclic loading between 100 and 750 N. Geometries were imported and into OpenSim and contact mechanics were modeled with the on-board elastic foundation algorithm. A hybrid optimization algorithm determined that stiffness and dissipation coefficients for TKR implants were 1.52 x 10(10) N/m and 57.7 Ns/m, respectively. Estimations of contact forces were 10.2% of blinded experimental data (average root mean square error: 76.82 +/- 11.47 N). In the second portion of this study, freely available eTibia TKR renderings were used to test the ubiquity of the tuning parameters. They were also used to perform a sensitivity analysis of material stiffness and mesh density with regard to penetration depth and computational time. When a stiffness of 1 x 10(10) was applied to an eTibia model with 5000 faces, a 100 kg load caused 0.259 mm of penetration. Under the same conditions, the tuned model experienced 0.300 mm of penetration. Material stiffnesses between 1 x 10(13) and 1 x 10(15) N/m increased computation time by factors of 12-23. This study provides much needed clarity regarding the use of the OpenSim EF algorithm. It also demonstrates the utility of OpenSim to model deformable materials and complex geometries, and this approach can be adapted to make reasonable estimations for both natural and surgically modified joints. (C) 2018 Elsevier Ltd. All rights reserved.