Nakamachi Eiji, Noma Tomohiro, Nakahara Kaito, Tomita Yoshihiro, Morita Yusuke
Department of Biomedical Engineering, Doshisha University, Kyotanabe, Kyoto, 610-0394, Japan.
Int J Numer Method Biomed Eng. 2017 Nov;33(11). doi: 10.1002/cnm.2864. Epub 2017 Mar 3.
The articular cartilage of a knee joint has a variety of functions including dispersing stress and absorbing shock in the tissue and lubricating the surface region of cartilage. The metabolic activity of chondrocytes under the cyclic mechanical stimulations regenerates the morphology and function of tissues. Hence, the stress evaluation of the chondrocyte is a vital subject to assess the regeneration cycle in the normal walking condition and predict the injury occurrence in the accidents. Further, the threshold determination of stress for the chondrocytes activation is valuable for development of regenerative bioreactor of articular cartilage. In this study, in both macroscale and microscale analyses, the dynamic explicit finite element (FE) method was used for the solid phase and the smoothed particle hydrodynamics (SPH) method was used for the fluid phase. In the homogenization procedure, the representative volume element for the microscale finite element model was derived by using the multiphoton microscope measured 3D structure comprising 3 different layers: surface, middle, and deep layers. The layers had different anisotropic structural and rigidity characteristics because of the collagen fiber orientation. In both macroscale and microscale FE analyses, the visco-anisotropic hyperelastic constitutive law was used. Material properties were identified by experimentally determined stress-strain relationships of 3 layers. With respect to the macroscale and microscale SPH models for non-Newtonian viscous fluid, the previous observation results of interstitial fluid and proteoglycan were used to perform parameter identifications. Biphasic multiscale FE and SPH analyses were conducted under normal walking conditions. Therefore, the hydrostatic and shear stresses occurring in the chondrocytes caused by the compressive load and shear viscous flow were evaluated. These stresses will be used to design an ex-vivo bioreactor to regenerate the damaged articular cartilage, where chondrocytes are seeded in the culture chamber. To know the stress occurred on and in the chondrocytes is vitally important not only to understand the normal metabolic activity of the chondrocyte but also to develop a bioreactor of articular cartilage regeneration as the knee joint disease treatment.
膝关节的关节软骨具有多种功能,包括分散组织中的应力和吸收冲击,以及润滑软骨表面区域。在周期性机械刺激下,软骨细胞的代谢活动可使组织的形态和功能得以再生。因此,软骨细胞的应力评估是评估正常行走状态下再生周期以及预测事故中损伤发生情况的关键课题。此外,确定软骨细胞激活的应力阈值对于关节软骨再生生物反应器的开发具有重要价值。在本研究中,在宏观和微观尺度分析中,固相采用动态显式有限元(FE)方法,流体相采用光滑粒子流体动力学(SPH)方法。在均匀化过程中,通过使用多光子显微镜测量的包含表面、中间和深层3个不同层的三维结构,推导出微观尺度有限元模型的代表性体积单元。由于胶原纤维的取向,这些层具有不同的各向异性结构和刚度特性。在宏观和微观有限元分析中,均采用了粘弹性各向异性超弹性本构定律。通过实验确定的3层应力-应变关系来识别材料特性。对于非牛顿粘性流体的宏观和微观SPH模型,利用先前关于组织液和蛋白聚糖的观察结果进行参数识别。在正常行走条件下进行了双相多尺度有限元和SPH分析。因此,评估了由压缩载荷和剪切粘性流在软骨细胞中产生的静水压力和剪应力。这些应力将用于设计体外生物反应器,以再生受损的关节软骨,其中软骨细胞接种在培养室中。了解软骨细胞上和内部产生的应力不仅对于理解软骨细胞的正常代谢活动至关重要,而且对于开发作为膝关节疾病治疗手段的关节软骨再生生物反应器也至关重要。