[发明专利]一种考虑应力状态的骨结构预测方法

权利要求书

1.一种考虑应力状态的骨结构预测方法，其特征在于，包括：

步骤一：确定要进行预测的骨结构的拓扑范围，并根据力学原理抽象出力学边界条件和载荷；

步骤二：对步骤一确定的骨结构进行建模；

步骤三：对建立的骨结构模型进行有限元网格划分，得到骨结构有限元网格模型；

步骤四：定义骨结构的初始表观密度，计算骨的材料性能参数弹性模量，并对每个单元定义材料属性；

步骤五：对骨结构有限元网格模型加载力学边界条件和载荷；

步骤六：对骨结构进行有限元分析求解；

步骤七：提取骨结构有限元模型各单元的主应力；

步骤八：根据各单元的主应力判断各单元的应力状态；

步骤九：基于各单元的应力状态，根据力学原理选择各单元适用的力学强度理论；

步骤十：根据确定的力学强度理论应用对应的相当应力计算公式进行各单元力学激励的计算；

步骤十一：应用骨重建控制方程计算下一步的表观密度：

步骤十二：根据各单元的新的表观密度计算各单元新的弹性模量材料参数：

步骤十三：对各单元重新定义材料参数；

步骤十四：骨量收敛判断，骨量趋于稳定，退出骨重建迭代循环，结束骨结构预测计算；否则转至步骤六开始下一步的迭代计算；

步骤十五：进入有限元分析后处理，进行各单元骨的表观密度显示，即得到预测的骨结构。

2.根据权利要求1所述一种考虑应力状态的骨结构预测方法，其特征在于，步骤四中的计算骨的材料性能参数弹性模量，具体为：

E＝Cρ^γ------------------------------(1)

式中E为弹性模量，ρ为骨的表观密度，用来表征骨内部结构特性；C、γ为常数。

3.根据权利要求1所述一种考虑应力状态的骨结构预测方法，其特征在于，步骤十中相当应力计算公式为：

式中S为力学激励，σ₁、σ₂、σ₃为主应力，且σ₁≥σ₂≥σ₃，μ为材料泊松比，[σ_t]为材料的抗拉许用应力，[σ_c]为材料的抗压许用应力。

4.根据权利要求1所述一种考虑应力状态的骨结构预测方法，其特征在于，步骤十一中计算下一步的表观密度：

ρ_t＝ρ_t+Δt+B(S-K)Δt-------------------------------(4)

式(3)为微分方程，式(4)为数值计算迭代公式，式中ρ为骨的表观密度，单位体积骨的质量，用来表征骨内部结构特性；B为骨重建速率；S为力学激励，用来表征力学效应的幅度，取相当应力作为力学激励；K为骨重建平衡态的力学激励；ρ_t为t时刻骨的表观密度，ρ_t+Δt为t+Δt时刻骨的表观密度，Δt为积分步长。