不同运动强度对骨质疏松影响的数值模拟

1、不同运动强度对骨质疏松影响的数值模拟         11-02-09 10:05:00     编辑：studa20                 作者：朱东，闫华，朱伟民，赵双龙【摘要】  探讨数值模拟不同运动强度对因废用引起的骨质疏松的影响。依据骨的功能适应性和力学调控系统理论，遵循骨重建的生理过程，

2、建立了一个带有时间历程的描述松质骨骨质疏松过程的计算模型，用C语言编制程序。得到了不同运动负荷对因废用引起的骨质疏松（孔隙率）的影响曲线。结果表明：较大的运动强度可以减缓因废用引起的骨质疏松。这对于运动防治骨质疏松有一定的参考价值。 【关键词】  基本多细胞单位；孔隙率；运动强度；骨质疏松；松质骨；模拟Computer Simulation on the Influence of    Different Exercise Intensity on OsteoporosisZHU Dong1，YAN Hua2，ZHU Weimin1，ZHAO Shuan

3、glong1    (1.The First Hospital， Jilin University， Changchun 130025，China;    2.Jilin Communications Polytechnic College，Changchun 130022；    3.The Mechanical Department of Jilin University,Changchun,130022)    Abstract：The influence of com

4、puter simulation of different exercise intensity of osteoporosis caused by disuse.We established a cancellous model of osteoporosis bone with time history based on bone functional adaptation and Frosts mechanical theory and physiological process of bone remodeling was programmed by C language. It pr

5、oduced the curve of osteoporosis influence caused by disuse under different exercise load conditions. The results show that porosity can be reduced with greater exercise intensity and slow osteoporosis of disuse.These results are valuable for the prevention and treatment of osteoporosis by exercise.

6、    Key words：Basic multicellular unit；Porosity ；Different exercise intensity；Osteoporosis；Cancellous；Computer simulation    1  引  言    随着社会的发展,世界已进入了老龄化社会,骨质疏松症(osteoporosis,OP) 的发生率也逐年增加，已成为威胁中老年人身体健康的疾病之一。运动疗法是一种非药物、非侵入性、无毒副作用的治疗骨质疏松症的方法，已成为骨

7、力学研究和骨质疏松研究的一个新的课题。    骨再造（包括骨建造和骨重建）贯穿人的一生，对其进行模拟是骨生物力学的一个研究热点。很多研究者已经用计算模型来验证数学规律控制骨的力学适应性的假设1，骨再造是生物学和力学因素耦合共同作用下发生的，因此只有综合考虑这两方面的耦联关系才有可能正确反映骨再造的过程。    2  材料和方法    2.1  材料    我们选择一个取自股骨近端的松质骨体积单元作为施加载荷的模型，见图1(a)，在此骨单元的上下两侧作用着均

8、布变化力。图1(b)是从图1(a)的体积单元中选取的具有代表性的横截面，里面分布着一些不规则的骨小梁，其中一些激活的基本多细胞单位(BMU)正在活动，有的进行骨吸收，有的进行骨填充,方框中显示的半圆形缺损是激活的BMU正在骨小梁表面进行骨重建，在松质骨中，骨小梁表面BMU吸收和填充的截面形状一般认为是沟状或槽状2，见图2。    2.2  方法    松质骨的骨质疏松模型是从现在已知的骨重建动态生理过程出发，细致描述BMU的破骨、成骨过程，融入合理的力学因素形成一个计算模型。    =(1-P)

9、×0(1)    其中，0取为2.0 g/cm3。    弹性模量E（MPa）是表示骨力学性能的一个重要指标，表示材料抵抗弹性变形能力的大小，孔隙率越大，表观密度和弹性模量就越小。我们通过表观密度与弹性模量的关系来表现孔隙率变化对骨组织力学性能（骨质疏松）的影响。朱兴华3等提出了一个较为科学合理的弹性模量和表观密度的分段关系式：    E=1.007×2  0.25    0.255×  0.250.4    2.9722-0.933  0.41.2    =ni=1siqRli（3）    是力学激励，我们考虑了一种载荷工况，s是图1(a)中所示工况体积单元产生的峰值应变，经计算为2.5×10-4。按正常中老年人每天行走2 km，每步0.75 m，则载荷频率Rl1300 cycles per day（cpd）。q是应变权重指数5，一般取值为4。由式（3）可得初始力学激励为：0=0.508×10-11cpd。    参加不同体育运