[发明专利]一种基于纤维方向的关节软骨的建模计算方法

说明书

技术领域

本发明涉及一种关节软骨有限元仿真技术，特别涉及一种基于纤维方向的关节软骨的建模计算方法

背景技术

软骨组织因其力学传递的功能，软骨没有血液供应，损伤后很难实现自我修复，损伤程度过大还会使关节伤残成为人体的支架不可或缺的重要部分，其生长、发育及对正常生理功能的维持都需要载荷的作用。建立真实的软骨模型，从而探索软骨损伤状态下的应力应变关系，进而提出预防和治疗机制。

在关节软骨研究建模初期，研究学者将关节软骨视为各向同性单相线弹性材料，然而，单相模型得出的结果与实际情况和实验情况的对比数据出现很大的差异。Mow和Lai等人基于混合物理论而提出了软骨的两相介质模型，使两相理论得到发展，即将软骨视为由线弹性固相和非粘性液相组成的不可压缩的混合物，由于与单相模型相比，两相模型更加精确，所以大部分的研究学者都已两相模型来对关节软骨进行研究，但是,两相模型只能很好的模拟围限压缩试验，在非围限实验中的模拟计算中与实际情况存在很大差异，不能够完整的表现出软骨的力学性能。

发明内容

本发明为了解决固液两相模型中在受载荷时无法准确表现其横向变形以及软骨各层次力学性质不同的问题。

一种基于纤维方向的关节软骨的建模计算方法，其特征是，具体步骤如下：

(1)几何模型的建立,在Abaqus中建立关节软骨几何模型，该软骨几何模型采用8节点轴对称两相多孔CAX8PR单元，在每个节点处计算固体和孔隙压力的总应力；

(2)建立关节软骨力学模型，并在其中考虑胶原纤维的分布以及不同层次的渗透率对关节软骨的影响，在fortran中进行材料子程序的编写，并嵌入到 abaqus中；

(3)在步骤(2)后，在CAE界面上的节点之间添加弹簧单元，在inp文件中修改弹簧单元参数，提交Abaqus中分析，用newton法求解控制方程，进行仿真分析。

根据权利要求1所述的一种基于纤维方向的关节软骨建模计算方法，其特征是，所述的步骤(1)的建立软骨几何模型，该几何模型分为两个层区，浅表层代表模型整个模型的最上方百分之十。

根据权利要求1所述的一种基于纤维方向的关节软骨建模计算方法，其特征是，所述的步骤(2)的建立关节软骨力学模型，其本构模型的基本框架为非线性超弹性两相多孔介质模型；

将关节软骨视为由基质固体，理想流体，和胶原纤维组成的三相混合物，且各个组相具有独立的运动规律，用S表示固体相，F表示液体相，由关节软骨的不可压缩性和力学特性得到关节软骨的力学本构方程为：

σ^F＝-n^FpI (2) 式中有效柯西应力张量，n^s为固体相体积分数，p代表液体压力，I是单位矩阵，ρ^S是固体密度，Ψ^S是固体Helmhooltz自由能方程，C_S为右柯西-格林张量，为固体相的材料变形梯度。

根据权利要求1所述的一种基于纤维方向的关节软骨建模计算方法，其特征是，所述步骤(2)中的考虑胶原纤维的分布对关节软骨的影响，在此非线性各向异性超弹性本构模型框架基础上考虑关节软骨的纤维特性，胶原纤维被假设为两个规整的纤维族，这些纤维族对称性的分布在软骨的轴向线两边，并且具有相同的刚度。这些纤维族被假定认为在关节软骨轴线具有一个参考方向，分别用和来表达平均纤维角度，在ABAQUS中定义了两类纤维，两类胶原纤维随着其对称轴发散，纤维方向的发散水平反映出纤维在软骨中的密度，其纤维优选方向由纤维密度函数来确定；

U代表应变能密度函数，C₁₀，k₁，k₂为材料参数，D为温度相关材料参数，是第一Green应变不变量，J是弹性体积比，是右Cauchy-Green张量偏差C与纤维族的平均纤维角度A的伪不变量。k代表纤维的分散度，N为纤维族数量，Ea为拟应变量；

在所建立的纤维模型中，每一族胶原纤维的方向有一个平均优先角度，平均优先角度沿该方向是分散的，定义为纤维族的松散度k，其表达式为：

ρ(θ)是方向密度函数，代表了在平均方向的范围[θ,θ+Δθ]中的规范化纤维的数量分布情况，当k＝0时，代表纤维都是沿同一方向排列的，没有分散，当 k＝1/3时，表示纤维是随机分布的，因此材料就变成了各向同性，纤维方向角度对模型的应变能函数没有影响。