Department of Mechanical Engineering, University of Colorado Denver, Denver, CO, USA.
Center of Bioengineering, University of Colorado Anschutz Medical Campus, Aurora, CO, USA.
Ann Biomed Eng. 2022 Dec;50(12):1954-1963. doi: 10.1007/s10439-022-03016-w. Epub 2022 Jul 21.
Cartilage contact pressures are major factors in osteoarthritis etiology and are commonly estimated using finite element analysis (FEA). FEA models often include subject-specific joint geometry, but lack subject-specific joint kinematics and muscle forces. Musculoskeletal models use subject-specific kinematics and muscle forces but often lack methods for estimating cartilage contact pressures. Our objective was to adapt an elastic foundation (EF) contact model within OpenSim software to predict hip cartilage contact pressures and compare results to validated FEA models. EF and FEA models were built for five subjects. In the EF models, kinematics and muscle forces were applied and pressure was calculated as a function of cartilage overlap depth. Cartilage material properties were perturbed to find the best match to pressures from FEA. EF models with elastic modulus = 15 MPa and Poisson's ratio = 0.475 yielded results most comparable to FEA, with peak pressure differences of 4.34 ± 1.98 MPa (% difference = 39.96 ± 24.64) and contact area differences of 3.73 ± 2.92% (% difference = 13.4 ± 11.3). Peak pressure location matched between FEA and EF for 3 of 5 subjects, thus we do not recommend this model if the location of peak contact pressure is critically important to the research question. Contact area magnitudes and patterns matched reasonably between FEA and EF, suggesting that this model may be useful for questions related to those variables, especially if researchers desire inclusion of subject-specific geometry, kinematics, muscle forces, and dynamic motion in a computationally efficient framework.
软骨接触压力是骨关节炎病因的主要因素,通常使用有限元分析 (FEA) 进行估计。FEA 模型通常包括特定于主体的关节几何形状,但缺乏特定于主体的关节运动学和肌肉力量。肌肉骨骼模型使用特定于主体的运动学和肌肉力量,但通常缺乏估计软骨接触压力的方法。我们的目标是在 OpenSim 软件中适应弹性基础 (EF) 接触模型,以预测髋关节软骨接触压力,并将结果与经过验证的 FEA 模型进行比较。EF 和 FEA 模型是为五个主体构建的。在 EF 模型中,应用运动学和肌肉力量,并根据软骨重叠深度计算压力。通过扰动软骨材料特性来找到与 FEA 压力最匹配的模型。EF 模型的弹性模量为 15 MPa,泊松比为 0.475,与 FEA 的结果最匹配,峰值压力差异为 4.34 ± 1.98 MPa(差异百分比= 39.96 ± 24.64),接触面积差异为 3.73 ± 2.92%(差异百分比= 13.4 ± 11.3)。5 个主体中有 3 个主体的 FEA 和 EF 的峰值压力位置匹配,因此,如果峰值接触压力的位置对研究问题至关重要,则不建议使用此模型。FEA 和 EF 之间的接触面积大小和模式匹配合理,表明该模型可能对与这些变量相关的问题有用,特别是如果研究人员希望在计算效率高的框架中包含特定于主体的几何形状、运动学、肌肉力量和动态运动。