为获得有限元计算模型的机械参数，需要进行生物力学实验，生物力学实验的静态分析适用于瞬态变化或平衡状态，可用于研究软组织在某一时刻的特征变化。而受试者年龄的变化会严重影响组织结构和机械性能，影响膝关节中各个结构材料属性的选择。在有限元建模过程中，按照受试者年龄段选择合适的材料特性和材料参数可以更好地描述膝关节的力学特性，这对提高有限元的精度至关重要。下面将分别描述不同年龄阶段膝关节各结构的模型材料特性。 2.1   膝关节有限元模型研究情况   膝关节有限元模型的关键时间点及重要研究介绍，见表1。STEELE等[10]使用0.5 T MRI单元对没有膝关节创伤史的运动女性右膝进行了成像。FUCHS等[11]提出了一种能够通过CT扫描确定膝关节表面坐标的方法。WAWRO等[12]提出一种新的非线性有限元分析软件。DONG等[13]提出了一个3D膝关节有限元模型，以确定软骨缺损对缺损边缘周围应力分布的影响。í?IGUEZ-MACEDO等[14]分析并比较了人造和健康膝关节有限元模型皮质骨和小梁骨的应力分布情况。CHEN等[15]基于CT图像数据建立动态膝关节有限元模型，并验证该模型的有效性。LAHKAR等[16]报道了一种从双平面X射线片和CT数据得到韧带预应力的新方法。LU等[17]通过模拟方法来预测反向运动跳跃期间膝关节内的应力分布，该模拟方法将肌肉骨骼模型与膝关节的有限元模型相结合。ALONSO等[18]施加恒定的负荷提高矫形器在治疗儿童期和青春期早期膝关节侧偏方面的有效性。"

2.2   骨(股骨、胫骨、腓骨、髌骨)    膝关节中骨(股骨、胫骨、腓骨、髌骨)的解剖学、生理学、生化、矿化和病理学方面随年龄的变化已有大量文献，而关于膝关节骨材料特性随年龄变化的文献很少。 2.2.1   骨密度与年龄之间的关系   皮质骨是一种有机结构，矿物质含量约占人骨骼的80%。影响人类皮质骨质量的病理因素包括成人骨质疏松症[19]、儿童骨质疏松症[20]、克罗恩-S病或骨化症[21-22]。相关研究发现，由于人口老龄化的加剧，预测到2050年骨质疏松性骨折的直接成本将达到767亿欧元[23]。儿童时期的低骨量被认为是晚年骨质疏松的高危因素[24]，相关研究认为评估儿童骨密度非常重要[25-26]。骨密度是评估矿物质状况的黄金标准参数之一，首先需要使用双能X射线吸收仪进来检测。相关研究使用霍洛奇双能X射线吸收仪(霍洛奇公司，Waltham，12 MA，USA)对7-11岁儿童进行骨密度扫描[27]。骨密度与骨微结构无关。多项研究表明，超声波测量能够评估探索骨区域(弹性和结构)的数量和质量[28-29]。KIM等[30]将儿童膝关节有限元计算模型中股骨材料的弹性模量与泊松比分别设置为17 000 MPa、0.3；LIMBERT等[31]与LANCIANESE等[32]将儿童胫骨材料的弹性模量与泊松比分别设置为12 200 MPa、0.3；BERTEAU等[33]将儿童腓骨材料的弹性模量与泊松比分别设置为15 500 MPa、0.24；KIAPOUR等[34]将儿童髌骨材料的弹性模量与泊松比分别设置为15 000 MPa、0.3。而有关成年人膝关节有限元计算模型中股骨、胫骨、腓骨、髌骨的弹性模量和泊松比设置为17 000 MPa、0.3[35-36]。这是由于骨密度随着年龄的增长而增加，在成年后达到峰值[37]。 2.2.2   骨密度与弹性模量之间的关系   王倞等[38]将10例成人新鲜尸体作为样本源，选取胫骨近端、大转子、股 骨颈、肱骨头和椎体5个部位的松质骨，采用CT扫描试样骨密度，通过力学测试分析分析不同部位松质骨的弹性模量，发现骨密度与弹性模量之间存在较高的线性相关性(0.850 > R2 > 0.785)和幂次相关性(0.871> R2 > 0.825)。袁嘉尧等[39]通过Pearson相关性分析发现，中老年人整体骨密度与身高、体质量、骨矿含量(整体、头部、主干、上肢、下肢)、肌肉质量(整体、头部、主干、下肢)、头部脂肪质量呈显著正相关，而与年龄呈显著负相关。从中年到老年阶段，膝关节中股骨、胫骨、腓骨、髌骨的弹性模量会随着年龄的增长而降低，进而回归到儿童时期的骨弹性模量。 2.2.3   年龄与弹性模量之间的关系   有关年轻骨骼强度和弹性模量的数据非常缺乏，多篇论文报道了超声轴向透射数据的年龄依赖性[40-42]，只有2项体外机械研究对儿童皮质骨弹性模量进行了定量研究[43-44]，这2项研究都是对干燥样本进行破坏性测试，儿童骨骼实验值都支持理论优化假设[45]，即骨骼弹性模量会从新生儿值上升到成人值，这一假设目前被用于儿科计算方法。儿童骨骼体外超声波速度和弹性模量值与现有文献中的老年成人值接近[46-47]。 评估儿童皮质骨弹性特性是生物力学工程界及更广泛医疗保健专业人员面临的重大挑战[48]。即使采用X射线断层扫描、MRI成像和超声波成像等经典临床模式，由于缺乏儿童骨骼的参考值，也可能导致诊断不当。相关研究利用超声波和机械测量方法提供了儿童(4-16岁，平均年龄10岁)皮质骨的弹性特性值，并与老年人(75 岁以上)进行了比较，首先通过超声波方法评估动态弹性模量和泊松比，其次通过三点微弯曲测试估算静态弹性模量，结果显示：儿童样本在 10 MHz 频率下测得的纵向和横向波速平均值分别为(3.2±0.5) mm/μs 和(1.8±0.1) mm/μs，老年人样本的波速分别为(3.5±0.2) mm/μs 和(1.9±0.1) mm/μs，儿童样本的平均动态弹性模量和平均静态弹性模量分别为 (15.5±3.4) GPa 和(9.1±3.5) GPa；老年人样本的平均动态弹性模量和平均静态弹性模量分别为(16.7±1.9) GPa和(5.8±2.1) GPa，表明儿童和老年人骨骼的动态弹性模量、泊松比和静态弹性模量在相同范围内，没有统计差异[49]。 UEHIRA[50]发现成年人皮质骨的抗压和抗拉强度都随着年龄的增长而下降。抗压强度的下降归因于骨密度的下降，抗拉强度的下降归因于哈弗管横截面积的增加。根据SUGIYAMA[51]的说法，皮质骨的洛氏硬度随着年龄的增长而降低，这种降低是由骨质疏松症引起的，而不是由于骨骼中钙含量的变化。BLACK等[52]对8例胫骨中段的成活皮质骨动态力学性能进行了研究，发现35.4 Hz处的拉伸弹性模数随年龄的增加而下降，而353.6 Hz处的拉伸弹性模数则没有这种变化。WILSON[53]研究了35-89岁个体24个骨骼标本的压缩强度和单位体积骨矿物质的变化，35岁以后骨骼标本的最大抗压强度每10年下降约7%，单位体积骨矿含量每10年下降约3.3%。 不同年龄阶段骨结构的材料特性，见表2。"

2.3   软骨   在膝关节有限元计算模型中，由于软骨结构光滑、柔软(刚性较小)、厚度较薄，因此在膝关节有限元建模过程中较为困难[54]，也很难追踪软骨形态随年龄的演变。SI等[55]利用3-T超导扫描仪采集膝关节MRI图像，通过深度学习自动、精确地分割软骨，计算出14个解剖区域的软骨厚度，发现股骨软骨厚度始终大于胫骨软骨相应部位，除胫骨软骨内侧和外侧外，所有区域的软骨都随着年龄的增长而变薄。BLAZEK等[56]发现软骨厚度与体质量×身高呈正相关。相关研究发现，兔膝关节软骨弹性模量与关节软骨厚度呈负相关[57]。 由上可知，随着年龄的增加股骨软骨厚度在减少，而随着股骨软骨厚度的减少其弹性模量在增加。但是DéMARTEAU等[58]通过比较老年人股骨软骨与成年牛肱骨软骨发现，老年人股骨软骨具有更大的弹性模量。作者认为这是由于研究不同类型软骨得出的结论，不影响前文的结论。该文对后续学者研究软骨厚度与弹性模量直接的线性或非线性关系具有一定的指导意义。 2.4   韧带   膝关节韧带是由软组织组成的弹性带，具有复杂的微观结构和生物力学特性，对膝关节的运动学特性和应力承受能力至关重要[59-60]。从弹性一维单元到具有解剖逼真形状的复杂超弹性三维结构，已被广泛应用于膝关节韧带研究[61]。在膝关节3D模型中包含韧带的最直观方法是使用实体元素[62]。MRI扫描与3D重建软件相结合，提供了一种方便的方法来创建关节的详细模型，包括特定患者的韧带几何形状及其插入位置[63]。MRI扫描也可以很容易与CT相结合，CT提供了骨骼结构的高分辨率成像，这种方法还便于通过使用面对面接触来模拟韧带包裹，这在有限元程序中都可以实现[64]。韧带的材料特性(弹性模量、泊松比)选择对于计算精度非常重要，相关研究发现，青壮年人体标本中韧带的弹性模量、极限拉伸应力和破坏应变大约是60岁及以上人体标本的两三 倍[65]，青壮年人体标本的主要破坏模式是韧带断裂，而老年人标本的主要破坏模式是韧带插入部位下方的骨骼撕裂，产生的差异和老年人标本韧带连接部位骨量减少有关。实验结果表明，随着年龄的增长，韧带的强度和刚度特性会显著降低，其程度超出预期。将年龄较大人体韧带标本实验结果应用于年轻或运动型个体的韧带破坏机制时，应谨慎解释，人们认为老年人的较差韧带生物力学特性减弱，同时关节固定也会导致韧带出现类似的功能性损伤。KASPERCZYK等[66]通过对6名年轻(平均年龄30岁)和6名年长(平均年龄64.7岁)活动水平相似供者的前交叉韧带和后交叉韧带进行了生物力学测试，发现年轻/年长前交叉韧带的最大应力为24/21 N/mm2)、年轻/年长前交叉韧带的弹性模量为144/129 MPa，年轻/年长前交叉韧带应变为0.25/0.28 mm。 综上所述，随着年龄的增长韧带弹性模量会出现一定程度的降低(出现减低的原因是由于随着年龄的增长韧带成纤维细胞的力学和生物学特性出现一定程度的退化)，但是年轻人与老年人韧带之间的弹性模量无统计学差异，表明老年人并不一定会出现交叉韧带功能性损伤。 2.5   半月板    2.5.1   半月板的纳米力学特征随年龄的变化情况   半月板为人类膝关节提供基本支持并包含一系列功能，例如调节膝关节的负荷传递[67-69]。半月板含有沿圆周排列的胶原纤维，能够分布胫骨平台上的载荷，从而强化关节软骨并预防骨关节炎[70-73]。半月板通过称为“半月板附着物”的韧带结构与胫骨相连，其基质主要由蛋白多糖、弹性蛋白、糖脂、成纤维细胞以及60%-70%的水分组成[74-75]。这些强化的胶原纤维帮助半月板将载荷扩散到附着点，逐渐将软组织过渡到骨骼中，减少应力集中[76-78]。 随着年龄的增长，半月板的细胞外基质发生分子和结构变化，导致基质退化[79-80]，这一退化过程的特征是合成代谢和分解代谢失衡，细胞减少了重要生长因子的产生，增加了基质降解酶和炎症因子的产生[81]。研究表明，衰老会导致半月板的胶原纤维排列紊乱、细胞减少等生化变化[77-78，82-84]。尽管衰老和骨关节炎引起的半月板退化已通过组织学和免疫组织化学分析研究，但尚未在微观结构层面详细研究其生物力学变化[85-87]。基于现有知识，衰老会改变半月板细胞外基质的生物力学特性，进而影响关节负荷时的反应，导致骨关节炎的发生[88]。研究还显示，随着年龄和疾病的进展，半月板的区域性机械性能逐渐模糊，机械完整性下降，这会破坏组织的生物力学功能，增加骨关节炎的风险。随着年龄增长，纳米力学性能的异质性逐渐增大[89-91]。 2.5.2   半月板的纳米力学特征与弹性模量之间的关系   正常、老年及骨关节炎患者的半月板组织在纳米力学特性上表现出显著差异[92]。年轻健康组织的弹性模量呈单峰分布，各区域表现出明确的峰值：外侧130-150 kPa、内侧70-90 kPa。老化组织同样呈现单峰分布，但其范围更宽，通常呈双峰，并且偏向较高的硬度值。整体上，双峰分布由一个较窄的第一峰和一个较宽的第二峰构成，表明退变相关的硬化区域有所增加。 骨关节炎不仅作为一种独立的疾病表现出特异的临床和生化特征，还随着组织的硬化而导致压缩刚度的增加，最终导致抗拉强度的下降，引发生物力学失效[93-95]。基质硬化源于细胞外基质中与年龄相关的分子变化，例如细胞外基质分子的异常沉积及交联，已知这些变化能够影响基质的机械特性，进而导致细胞功能异常[96-100]。半月板细胞对机械应力高度敏感，尤其是在关节受力时。细胞外基质中这些变化会促发或加剧与骨关节炎相关的细胞表型改变[101-103]。然而，基质力学与骨关节炎发展的确切关联仍有待进一步研究。 不同年龄阶段软骨、韧带与半月板结构的材料特性，见表3-5。"

综上所述，人类半月板组织的局部纳米力学特性会随着年龄的增长而发生变化。文献表明，儿童与成人半月板的弹性模量大致相同，但在从成人到老年阶段，弹性模量逐渐降低[35-36，105，108]，这一变化与半月板退行性变和骨关节炎的发生密切相关。因此，未来研究应深入探讨半月板的生物力学和生化特性之间的关联，尤其是在退变和骨关节炎的发展过程中。了解这些年龄相关的退变机制，将为预防和治疗骨关节炎提供重要的理论依据和策略支持。"