一种骨小梁-骨水泥复合体有限元模型构建方法

1.本发明涉及有限元模型构建技术领域，具体涉及一种骨小梁-骨水泥复合体有限元模型构建方法。

背景技术：

2.现有技术中，建立骨小梁-骨水泥界面作用的有限元模型主要包括2种方法：
3.方法一：1、获取经骨水泥强化固定的人尸体骨，钻取骨小梁-骨水泥复合体样本；2、将骨小梁-骨水泥复合体样本进行micro-ct扫描，获取影像数据；3、计算机逆向工程，重建骨小梁-骨水泥复合体样本的计算机三维模型；4、划分体网格、赋予材料属性、参数，设置边界条件、载荷，进行有限元力学分析。
4.方法二：1、获取人尸体骨，钻取骨小梁样本；2、对骨小梁进行骨水泥灌注；3、待骨水泥完全固化后，对骨小梁-骨水泥复合体样本进行micro-ct扫描；其余步骤与方法1相同。
5.如图1所示，为在对骨小梁-骨水泥复合体样本进行micro-ct扫描，得到扫描影像后，构建骨小梁-骨水泥复合体有限元模型的过程示意图。
6.采用现有技术中的上述方法构建骨小梁-骨水泥复合体有限元模型存在以下缺点：1、获取的经骨水泥强化固定的人尸体骨，骨小梁往往存在不同程度的骨溶解、骨吸收情况，并导致大量孔隙形成，需要在计算机软件进行反向填充，以获取骨小梁-骨水泥完整接触的三维模型，步骤繁杂、耗时巨大。2、骨小梁与骨水泥灰度阈值部分重叠，在计算机软件中分离两者容易造成人为的错误分配。3、经骨水泥强化固定的人尸体骨来源非常有限，且人尸体骨成本高昂，难以推广应用。4、人工灌注骨水泥仍然存在一定数量的孔隙，影响有限元分析的精确度，三维模型仍然需要人工反向填充。

技术实现要素：

7.针对上述问题，本发明提供一种骨小梁-骨水泥复合体有限元模型构建方法。
8.本发明提供的骨小梁-骨水泥复合体有限元模型构建方法包括：
9.获取骨小梁微观影像数据；
10.在骨小梁微观影像数据中选取感兴趣区域，在感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像；
11.对骨小梁三维结构图像中的骨小梁进行打磨光滑处理；
12.根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体；
13.设置网格尺寸，选中流体域，对流体域进行二维面网格细分，并在细分后对流体域进行三维体网格划分；
14.选中骨小梁，对骨小梁进行三维体网格划分；
15.对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给骨小梁-骨水泥复合体中的骨水泥、骨小梁网格赋予
对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。
16.进一步的，获取骨小梁微观影像数据包括：
17.对人体进行高分辨率外周定量ct扫描；
18.或者，对动物骨骼标本或骨科手术过程中产生的废弃骨标本进行micro-ct扫描。
19.进一步的，若骨小梁微观影像数据通过高分辨率外周定量ct扫描得到，在骨小梁微观影像数据中选取感兴趣区域，在感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像包括：
20.将骨小梁微观影像数据导入mimics软件，调整窗宽和窗位以选取感兴趣区域，在感兴趣区域进行阈值分割和区域增长操作，去除孤立结构单元，得到边界清晰的骨小梁三维结构图像并保存为stl格式文件；
21.若骨小梁微观影像数据通过micro-ct扫描得到，在骨小梁微观影像数据中选取感兴趣区域，在感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像包括：
22.将骨小梁微观影像数据导入ctan软件，选取骨小梁微观影像数据的上、下边界，选取感兴趣区域，对感兴趣区域二值化，得到骨小梁三维结构图像并保存为stl格式文件。
23.进一步的，对骨小梁三维结构图像中的骨小梁进行打磨光滑处理为：
24.将骨小梁三维结构图像导入geomagic studio软件，利用geomagic studio软件对骨小梁三维结构图像中的骨小梁进行打磨光滑处理，并将打磨光滑处理后的骨小梁三维结构图像保存为stl格式文件；
25.其中，打磨光滑处理包括：自动分析和修复骨小梁模型、填充孔、删除钉状物并重画骨小梁的二维面网格。
26.进一步的，根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体为：
27.将打磨光滑处理后的骨小梁三维结构图像导入至ansys 19.0软件，在ansys19.0软件中设置流体域参数，根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体，得到包裹后的骨小梁三维结构图并保存为stl格式文件。
28.进一步的，设置网格尺寸，选中流体域，对流体域进行二维面网格细分，并在细分后对流体域进行三维体网格划分为：
29.将包裹后的骨小梁三维结构图导入hypermesh 14.0软件，利用hypermesh14.0软件设置网格尺寸，选中流体域，对流体域进行二维面网格细分，并在细分后对流体域进行三维体网格划分。
30.进一步的，选中骨小梁，对骨小梁进行三维体网格划分为：
31.在hypermesh 14.0软件中，选中骨小梁，对骨小梁进行三维体网格划分。
32.进一步的，对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给骨小梁-骨水泥复合体中的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型为：
33.在hypermesh 14.0软件中，对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给骨小梁-骨水泥复合体中
的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。
34.进一步的，流体域参数包括流体域长度、流体域宽度以及流体域高度。
35.进一步的，材料属性及材料参数包括弹性模量、泊松比。
36.本发明提供的一种骨小梁-骨水泥复合体有限元模型构建方法，至少包括以下有益效果：
37.(1)本发明提供的构建方法，无需预先获取/制备骨小梁-骨水泥复合体样本，只需获取骨小梁微观影像数据即可进行余下的骨小梁-骨水泥复合体有限元模型构建步骤，操作快速、简便，成本低，易推广。
38.(2)本发明提供的构建方法，由于不需要利用骨水泥对尸骨进行强化固定，得到骨小梁-骨水泥复合体样本，所建立的有限元模型不存在孔隙问题，骨小梁与骨水泥完全贴合，提高有限元分析的精准度，分析结果精度高。
附图说明
39.为了更清楚的说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见的，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它附图。
40.图1为现有技术中骨小梁-骨水泥复合体有限元模型的构建过程示意图；
41.图2为本发明一种实施例中的骨小梁-骨水泥复合体有限元模型构建方法流程图；
42.图3为本发明一种实施例中的骨小梁-骨水泥复合体有限元模型构建流程示意图。
具体实施方式
43.下面将结合本发明中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通的技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例，都属于本发明的保护范围。
44.在本发明的一种实施例中，提供一种骨小梁-骨水泥复合体有限元模型构建方法，如图2所示，方法包括以下步骤：
45.步骤s1：获取骨小梁微观影像数据；
46.其中，骨小梁微观影像数据，可以通过对人体进行高分辨率外周定量ct(hr-qct)扫描得到。具体的，人体仰卧位并伸直下肢，在腰部底下加垫固体体模，对腰1-腰5椎体由近端向远端进行螺旋扫描，扫描过程中使扫描线与椎体上下缘平行，扫描分辨率40～60μm，将扫描得到的骨小梁影像数据用标准dicom格式保存。
47.骨小梁微观影像数据，还可以通过对动物骨骼标本或骨科手术过程中产生的废弃骨标本进行micro-ct扫描得到。
48.具体的，对动物骨骼标本进行micro-ct扫描的具体步骤为：对sd大鼠经麻醉法或颈椎脱位法处死后快速离体骨骼标本(如股骨、胫骨、腰椎骨等)，4％多聚甲醛固定48-72小时后进行micro-ct扫描，扫描方向为由骨骼近端向远端垂直于骨骼纵轴进行扫描，扫描分
辨率15μm，扫描得到的骨小梁微观影像数据以bmp格式保存。
49.对骨科手术过程中产生的废弃骨标本(如股骨头)进行micro-ct扫描的具体步骤为：在髋部骨折患者进行髋关节置换术过程中收集股骨头标本，4％多聚甲醛固定48-72小时后进行micro-ct扫描，扫描方向为由股骨头近端向远端垂直于股骨颈纵轴进行扫描，扫描分辨率40μm，以bmp格式保存骨小梁微观影像数据。
50.其中，骨小梁微观影像数据包括多张骨小梁微观影像图像。
51.步骤s2：在骨小梁微观影像数据中选取感兴趣区域，在感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像。
52.其中，若骨小梁微观影像数据是通过对人体进行高分辨率外周定量ct(hr-qct)扫描得到，则步骤s2具体为：将骨小梁微观影像数据导入mimics软件，调整窗宽和窗位(也即调整感兴趣区域的大小和位置)，以在骨小梁微观影像上选取感兴趣区域。在感兴趣区域中进行阈值分割和区域增长操作，去除孤立结构单元，得到边界清晰的骨小梁三维结构图像，保存为stl格式文件。
53.其中，在感兴趣区域中进行阈值分割是根据设定的预设阈值，将感兴趣区域中的骨小梁与背景进行分割。在进行阈值分割后，再利用mimics软件中的区域增长功能去除孤立结构单元，进而得到边界清晰的骨小梁三维结构图像。
54.若骨小梁微观影像数据是通过对动物骨骼标本或骨科手术过程中产生的废弃骨标本进行micro-ct扫描得到，则步骤s2具体为：将骨小梁微观影像数据导入ctan软件，选取骨小梁微观影像集的上、下边界，选取感兴趣区域，对感兴趣区域二值化，得到骨小梁三维结构图像并将骨小梁三维结构图像保存为stl格式文件。
55.其中，骨小梁微观影像集，指的是骨小梁微观影像数据中包含的多张骨小梁微观影像图像的集合。选取上下边界的作用是确定模型的高度(也即感兴趣区域的高度)，如图3中a上的黄矩形的高度就是通过选取上下边界确定。
56.进一步的，在选取感兴趣区域，对感兴趣区域二值化以后，还可以导出感兴趣区域中骨小梁的形态学参数。
57.其中，感兴趣区域中骨小梁的形态学参数包括骨小梁体积分数bv/tv、骨小梁厚度trabecular thickness，骨小梁数量trabecular number等，用来定量描述骨小梁结构质量。保存感兴趣区域中骨小梁的形态学参数，后续可以比较不同感兴趣区域之间的结构差异，并进行统计分析。
58.如图3中的a中方框中选取的即为感兴趣区域。
59.步骤s3：对骨小梁三维结构图像中的骨小梁进行打磨光滑处理。
60.其中，对骨小梁三维结构图像中的骨小梁进行打磨光滑处理为：
61.将骨小梁三维结构图像导入geomagic studio软件，利用geomagic studio软件对骨小梁三维结构图像中的骨小梁进行打磨光滑处理，并将打磨光滑处理后的骨小梁三维结构图像保存为stl格式文件。打磨光滑处理后，即可获取到光整的外表面，在geomagic studio软件中保存为stl格式文件时可自动在外表面生成二维面网格。
62.其中，打磨光滑处理包括：使用网格医生命令自动分析和修复骨小梁模型、填充孔、删除钉状物并重画骨小梁的二维面网格(细化二维面网格尺寸并使其均匀分布，尺寸大小为0.04mm)等。其中，当保存为stl格式文件时会在骨小梁的外表面自动生成二维面网格
(也即重画骨小梁的二维面网格)。
63.如图3中的b即为打磨光滑处理后得到的骨小梁三维结构图像。
64.步骤s4：根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体。
65.其中，步骤s4具体包括：将打磨光滑处理后的骨小梁三维结构图像导入至ansys 19.0软件，在ansys 19.0软件中设置流体域参数，根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体，得到包裹后的骨小梁三维结构图并保存为stl格式文件。
66.其中，在ansys 19.0软件中，以spaceclaim模块打开骨小梁三维结构图像，选择“enclosure”功能对骨小梁进行包裹，设置流体域参数后将其转换为实体，得到包裹后的骨小梁三维结构图并保存为stl格式文件。
67.其中，流体域参数包括流体域长度、流体域宽度以及流体域高度。
68.如图3中的c即为包裹后的骨小梁三维结构图。
69.步骤s5：设置网格尺寸，选中流体域，对流体域进行二维面网格细分，并在细分后对流体域进行三维体网格划分。
70.其中，步骤s5具体包括：将包裹后的骨小梁三维结构图导入hypermesh 14.0软件，利用hypermesh 14.0软件设置网格尺寸，选中流体域，对流体域进行二维面网格细分，并在细分后对流体域进行三维体网格划分。
71.其中，网格尺寸设置依据为网格收敛性分析，即网格无关性验证。在本实施例中，基于网格收敛性分析，可将网格尺寸设置为0.04mm。同时，在本实施例中，网格尺寸设置是在hypermesh 14.0软件菜单栏“2d”选项卡中选中“automesh”功能，从而设置网格尺寸。
72.设置网格尺寸之后，仍在“automesh”功能下选中流体域，对其进行二维面网格细分，从而获得均匀分布的二维面网格并为下一步的三维体网格划分提供模板。如图3中的d所示。
73.在细分后对流体域进行三维体网格划分，是在hypermesh 14.0软件菜单栏“3d”选项卡中选中“tetramesh”功能，选中“tetra mesh”，选中流体域，从而对其进行三维体网格划分。
74.步骤s6：选中骨小梁，对骨小梁进行三维体网格划分；
75.同样的，本步骤中对骨小梁进行三维体网格划分，也为在hypermesh 14.0软件菜单栏“3d”选项卡中选中“tetramesh”功能，选中“tetra mesh”，选中骨小梁，从而对其进行三维体网格划分。
76.步骤s5、步骤s6中的体网格划分均是由hypermesh 14.0软件根据布尔逻辑运算自动区分流体域及骨小梁并自动填充骨小梁间隙。如图3中的e所示。
77.步骤s7：对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给骨小梁-骨水泥复合体中的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。
78.其中，步骤s7具体为：在hypermesh 14.0软件中，对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给骨小
梁-骨水泥复合体中的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。如图3中的f所示。
79.进一步的，预设填充深度由技术人员根据实际需要自行设定，本发明对此不作限制。
80.材料属性包括：固体、管道、壳等。
81.材料参数包括弹性模量、泊松比、密度。骨水泥与骨小梁之间的接触类型包括摩擦、绑定等。
82.其中，材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，均由技术人员根据实际需求自行设定，本发明对此不做限制。
83.更具体的，在hypermesh 14.0软件中，在菜单栏“tool”选项卡中选中“delete”功能，选中骨水泥网格(即流体域网格)进行部分裁剪、删除，从而获得预设填充深度的骨小梁-骨水泥复合体。
84.进一步的，当获得上述的骨小梁-骨水泥复合体有限元模型后，即可对模型运算求解，获取相应的实验数据，如位移、应力。如图3中的g。
85.本发明提供的骨小梁-骨水泥复合体有限元模型构建方法，由于构建有限元模型时，无需获取或者制备骨小梁-骨水泥复合体样本，只需要获取骨小梁微观影像数据，后续利用骨小梁微观影像数据重构骨小梁三维模型，并在骨小梁外部包裹流体域，将流体域作为骨水泥，最终得到骨小梁-骨水泥复合体有限元模型。无需用到人尸体骨，大大降低了成本；并且生成骨小梁-骨水泥复合体有限元模型过程中不存在骨溶解、骨吸收问题，也不存在因人工灌注骨水泥而造成的大量孔隙的生成问题，无需反向填充，操作快速、简便，易于推广。同时，不存在骨小梁与骨水泥灰度阈值部分重叠，建模对两者进行分离时出现错误分配的问题。所建立的有限元模型不存在孔隙问题，骨小梁与骨水泥完全贴合，提高有限元分析的精准度，分析结果精度高。
86.本发明说明书中使用的术语和措辞仅仅为了举例说明，并不意味构成限定。本领域技术人员应当理解，在不脱离所公开的实施方式的基本原理的前提下，对上述实施方式中的各细节可进行各种变化。因此，本发明的范围只由权利要求确定，在权利要求中，除非另有说明，所有的术语应按最宽泛合理的意思进行理解。

技术特征：

1.一种骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述方法包括：获取骨小梁微观影像数据；在所述骨小梁微观影像数据中选取感兴趣区域，在所述感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像；对所述骨小梁三维结构图像中的骨小梁进行打磨光滑处理；根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体；设置网格尺寸，选中流体域，对所述流体域进行二维面网格细分，并在细分后对所述流体域进行三维体网格划分；选中骨小梁，对所述骨小梁进行三维体网格划分；对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给所述骨小梁-骨水泥复合体中的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。2.根据权利要求1所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述获取骨小梁微观影像数据包括：对人体进行高分辨率外周定量ct扫描；或者，对动物骨骼标本或骨科手术过程中产生的废弃骨标本进行micro-ct扫描。3.根据权利要求2所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，若所述骨小梁微观影像数据通过高分辨率外周定量ct扫描得到，所述在所述骨小梁微观影像数据中选取感兴趣区域，在所述感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像包括：将所述骨小梁微观影像数据导入mimics软件，调整窗宽和窗位以选取感兴趣区域，在所述感兴趣区域进行阈值分割和区域增长操作，去除孤立结构单元，得到骨小梁三维结构图像并保存为stl格式文件；若所述骨小梁微观影像数据通过micro-ct扫描得到，所述在所述骨小梁微观影像数据中选取感兴趣区域，在所述感兴趣区域内重建骨小梁三维结构，得到骨小梁三维结构图像包括：将所述骨小梁微观影像数据导入ctan软件，选取骨小梁微观影像数据的上、下边界，选取感兴趣区域，对感兴趣区域二值化，得到骨小梁三维结构图像并保存为stl格式文件。4.根据权利要求3所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，对所述骨小梁三维结构图像中的骨小梁进行打磨光滑处理为：将所述骨小梁三维结构图像导入geomagic studio软件，利用geomagic studio软件对所述骨小梁三维结构图像中的骨小梁进行打磨光滑处理，并将打磨光滑处理后的骨小梁三维结构图像保存为stl格式文件；其中，所述打磨光滑处理包括：自动分析和修复骨小梁模型、填充孔、删除钉状物并重画骨小梁的二维面网格。5.根据权利要求4所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体为：将打磨光滑处理后的骨小梁三维结构图像导入至ansys 19.0软件，在ansys19.0软件
中设置流体域参数，根据设置的流体域参数，将打磨光滑后的骨小梁利用流体域进行包裹后转换为实体，得到包裹后的骨小梁三维结构图并保存为stl格式文件。6.根据权利要求5所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述设置网格尺寸，选中流体域，对所述流体域进行二维面网格细分，并在细分后对所述流体域进行三维体网格划分为：将包裹后的骨小梁三维结构图导入hypermesh 14.0软件，利用hypermesh14.0软件设置网格尺寸，选中流体域，对所述流体域进行二维面网格细分，并在细分后对所述流体域进行三维体网格划分。7.根据权利要求6所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述选中骨小梁，对所述骨小梁进行三维体网格划分为：在所述hypermesh 14.0软件中，选中骨小梁，对所述骨小梁进行三维体网格划分。8.根据权利要求7所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给所述骨小梁-骨水泥复合体中的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型为：在hypermesh 14.0软件中，对划分后的流体域网格进行部分剪裁、删除，以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给所述骨小梁-骨水泥复合体中的骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。9.根据权利要求1所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，所述流体域参数包括流体域长度、流体域宽度以及流体域高度。10.根据权利要求1所述的骨小梁-骨水泥复合体有限元模型构建方法，其特征在于，材料参数包括弹性模量、泊松比。

技术总结

本发明涉及一种骨小梁-骨水泥复合体有限元模型构建方法，方法包括获取骨小梁微观影像数据；在感兴趣区域内重建骨小梁三维结构；对骨小梁三维结构图像中的骨小梁进行打磨光滑处理；利用流体域进行包裹后转换为实体；对流体域进行二维面网格细分，并在细分后对流体域进行三维体网格划分；对骨小梁进行三维体网格划分；以流体域作为骨水泥，获得预设填充深度的骨小梁-骨水泥复合体，分别给骨水泥、骨小梁网格赋予对应的材料属性及材料参数，设置边界条件、载荷及骨水泥与骨小梁之间的接触类型，得到骨小梁-骨水泥复合体有限元模型。本发明的方法省略了传统的构建骨小梁-骨水泥复合体样本的过程，建立的有限元模型精度高，成本低，容易推广。容易推广。容易推广。