Mohammadi H, Pietruszczak S, Quenneville C E
Department of Civil Engineering, McMaster University, Hamilton, Ontario, Canada.
Department of Civil Engineering, McMaster University, Hamilton, Ontario, Canada.
J Mech Behav Biomed Mater. 2021 Mar;115:104283. doi: 10.1016/j.jmbbm.2020.104283. Epub 2020 Dec 16.
The primary purpose of this paper is to outline a methodology for evaluating the likelihood of cortical bone fracture in the proximal femur in the event of a sideways fall. The approach includes conducting finite element (FE) analysis in which the cortical bone is treated as an anisotropic material, and the admissibility of the stress field is validated both in tension and compression regime. In assessing the onset of fracture, two methodologies are used, namely the Critical Plane approach and the Microstructure Tensor approach. The former is employed in the tension regime, while the latter governs the conditions at failure in compression. The propagation of localized damage is modeled using a constitutive law with embedded discontinuity (CLED). In this approach, the localized deformation is described by a homogenization procedure in which the average properties of cortical tissue intercepted by a macrocrack are established. The key material properties governing the conditions at failure are specified from a series of independent material tests conducted on cortical bone samples tested at different orientations relative to the loading direction. The numerical analysis deals with simulations of experiments involving the sideways fall, and the results are compared with the experimental data. This includes both the evolution of fracture pattern and the local load-displacement characteristics. The proposed approach is numerically efficient, and the results do not display a pathological mesh-dependency. Also, in contrast to the XFEM approach, the analysis does not require any extra degrees of freedom.
本文的主要目的是概述一种方法,用于评估在发生侧方跌倒时股骨近端皮质骨骨折的可能性。该方法包括进行有限元(FE)分析,其中将皮质骨视为各向异性材料,并在拉伸和压缩状态下对应力场的可接受性进行验证。在评估骨折的发生时,使用了两种方法,即临界平面方法和微观结构张量方法。前者用于拉伸状态,而后者则控制压缩失效条件。使用具有嵌入式不连续性的本构定律(CLED)对局部损伤的扩展进行建模。在这种方法中,局部变形通过均匀化过程来描述,在该过程中确定宏观裂纹截获的皮质组织的平均特性。控制失效条件的关键材料特性是根据对相对于加载方向以不同取向测试的皮质骨样本进行的一系列独立材料测试确定的。数值分析处理涉及侧方跌倒的实验模拟,并将结果与实验数据进行比较。这包括骨折模式的演变和局部载荷 - 位移特性。所提出的方法在数值上是有效的,并且结果不会显示出病态的网格依赖性。此外,与扩展有限元法(XFEM)方法不同,该分析不需要任何额外的自由度。