CN108281193B - 一种股骨头力学仿真模拟模型处理方法 - Google Patents
一种股骨头力学仿真模拟模型处理方法 Download PDFInfo
- Publication number
- CN108281193B CN108281193B CN201810062242.1A CN201810062242A CN108281193B CN 108281193 B CN108281193 B CN 108281193B CN 201810062242 A CN201810062242 A CN 201810062242A CN 108281193 B CN108281193 B CN 108281193B
- Authority
- CN
- China
- Prior art keywords
- femoral head
- simulation
- unit
- elastic
- extracting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供一种股骨头力学仿真模型处理方法,方法包括:根据灰度值对计算机断层扫描CT图像提取股骨头几何模型;将股骨头几何模型进行网格划分,得到与股骨头几何模型对应股骨头仿真模拟单元,将灰度值转化成密度值后,进而得到仿真模型单元数量百分比‑密度值关系曲线;将股骨头仿真模拟单元划分为多个部分;确定各个部分的股骨头仿真模拟单元的弹塑性本构参数;对添加弹塑性本构参数的股骨头仿真模拟单元进行预处理,输出仿真结果;对仿真结果进行提取和分析。将股骨头根据密度值等分成多个部分,很好的处理了股骨头非均质材料特性的问题。
Description
技术领域
本发明涉及骨力学仿真模拟技术领域,尤其涉及一种股骨头力学仿真模型处理方法。
背景技术
骨科临床治疗方案的制定和治疗效果严重依赖于骨的力学属性。然而,由于可用于实验的人骨样本稀缺,可以取到的样品尺寸过小,加之骨是具有多级网状复杂结构的材料,这些原因都加大了实验的难度,因此采用实验的手段去研究骨的力学属性受到限制。取而代之,有限元仿真模拟就成为了分析骨力学行为的一个简便快捷的方法。
股骨头骨折是老年常见病,其原因在于骨质疏松导致股骨头力学属性改变,遇到难以承受的外部载荷而发生骨折。现有针对股骨头的力学仿真模拟方法存在的问题主要有:1、股骨头是一种多级网状、非各向同性的材料,力学参数在不同位置、不同方向都存在差异,而目前模拟多简化为各向同性的均匀材料模型,忽略了这些要素;2、股骨头材料在受到外部载荷时,实际上会呈现出先弹性后塑性(屈服)的行为,然而目前的模拟方法多是把股骨头当作弹性材料处理,这和实际情况不符;3、用应力值作为判断股骨头损坏的参数不够准确,不同样品松质骨试件在屈服时对应的应力会有较大的随机性,但对应的应变相对稳定,这是由股骨头本身的网状结构决定的。
因此,现有仿真技术,模拟结果不尽人意,和临床应用尚有较大差距。
发明内容
本发明实施例提供一种股骨头力学仿真模型处理方法,用以解决现有仿真技术,模拟结果不尽人意,和临床应用尚有较大差距的问题。
本发明实施例提供一种股骨头力学仿真模型处理方法,包括:
根据灰度值对计算机断层扫描CT图像提取股骨头几何模型;
将所述股骨头几何模型进行网格划分,得到与所述股骨头几何模型对应的股骨头仿真模拟单元,将灰度值转化成密度值后,进而得到仿真模型单元数量百分比-密度值关系曲线;
将所述股骨头仿真模拟单元划分为多个部分,其中,划分规则为所述股骨头仿真模型单元按密度值均等划分后的单元数量百分比-密度值所形成的柱体总面积与所述曲线所形成的总面积的偏差不大于5%。
确定各个部分股骨头仿真模拟单元的弹塑性本构参数;
对添加所述弹塑性本构参数的所述股骨头仿真模拟单元进行预处理,输出仿真结果;
对所述仿真结果进行提取和分析。
在一可能的实施例中,所述对添加所述弹塑性本构参数的所述股骨头仿真模拟单元进行预处理,输出仿真结果,包括:
对所述添加所述弹塑性本构参数的所述股骨头仿真模拟单元再次添加边界条件,经加载后输出仿真结果。
在一可能的实施例中,所述将所述股骨头仿真模拟单元划分为多个部分,包括:
采用密度值将所述股骨头仿真模拟单元均等划分为多个部分。
在一可能的实施例中,所述确定各部分股骨头仿真模拟单元的弹塑性本构参数,包括:
对所述股骨头仿真模拟单元采用压缩线性强化弹塑性本构和拉伸线性强化弹塑性本构。
在一可能的实施例中,所述划分规则还包括:所述股骨头仿真模拟单元的总数量为N,划分部分数量为M,每个部分的单元数为n,每个部分等分区间的数量取2Mn/N的向上整数。
在一可能的实施例中,所述对所述仿真结果进行提取和分析,包括:
以应变超过屈服点的应变作为股骨头松质骨损坏的判断参数,对所述输出结果进行提取和分析;
其中,所述压缩屈服点的应变为0.0085,所述拉伸屈服点的应变为0.006,所述拉伸断裂点的应变为0.009。
将股骨头根据密度值等分成多个部分,很好的处理了股骨头非均质材料特性的问题。由于每个部分都采用线性强化弹塑性模型,且各自本构唯一性由该部分股骨头本身决定,使仿真模拟结果和真实情况更加吻合。在结果分析时,用应变值作为参数来判断松质骨的受力响应,更能反应股骨头内部的多孔网状结构,比用应力值作为参数来判断更加符合实际情况,结果也更加稳定。
附图说明
图1为本发明实施例提供的一种股骨头力学仿真模型处理方法流程示意图;
图2为本发明按密度值等分后的单元数量百分比-密度值柱状示意图;
图3为本发明压缩线性强化弹塑性力学模型示意图;
图4为本发明拉伸线性强化弹塑性力学模型示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图和实施例,对本发明实施例中的技术方案进行清楚地描述。
图1为本发明实施例提供的一种股骨头力学仿真模拟方法流程示意图,参照图1,该方法具体包括:
S101、根据灰度值对计算机断层扫描CT图像提取股骨头几何模型。
S102、将所述股骨头几何模型进行网格划分,得到与所述股骨头几何模型对应股骨头仿真模拟单元,将灰度值转化成密度值后,进而可以得到仿真模型单元数量百分比-密度值关系曲线。
S103、将所述股骨头仿真模拟单元划分为多个部分。
具体地,采用密度值将所述股骨头仿真模型单元划分为多个部分;其中,划分规则为所述股骨头仿真模型单元按密度值均等划分后的单元数量百分比-密度值所形成的柱体总面积与所述曲线所形成的总面积的偏差不大于5%。
以密度值相差50kg/m3做为一个区间进行划分,分成10个部分后,得到如图2所示的柱状图,柱体总面积与所述曲线所形成的总面积的偏差不大于5%;股骨头仿真模拟单元的总数量为N,每个部分的单元数为n,每个部分等分区间的数量取20n/N的向上整数可选地,在本实施例中,还可以根据实际需要按照划分规则,将股骨头仿真模型单元划分为多个部分,如划分为八个部分,对此本实施例不作具体限定。
通过以下公式赋予股骨头力学参数(弹性模量E),
灰度值与密度之间的公式:
ρ=131+1.067CT
式中ρ为密度值,单位为kg/m3;CT为CT值;
密度与弹性模量之间的公式:
E=0.01ρ1.86
式中ρ为密度值,单位为kg/m3,E为弹性模量。
每个密度值可根据股骨头的灰度值得到,在根据密度值赋予力学参数时,取该密度值区间的平均值进行计算。
图2所示为根据股骨头密度值区间,将密度值每隔50kg/m3作为一个区间进行划分,得到的股骨头单元数量百分比-密度值柱状示意图,整个股骨头由这10个不同的部分组成,每个部分等分区间的数量为1。
S104、确定所述至少各个部分股骨头仿真模拟单元的弹塑性本构参数。
将S103所得的各个部分的所述股骨头仿真模拟单元再导入到有限元分析软件中,步骤三所得弹性模量E为本构曲线弹性部分,如图三所示,该模型取压缩线性强化弹塑性本构和拉伸线性强化弹塑性本构关系(应变值压缩为负,拉伸为正),线性强化塑性部分的斜率取0.05E。
图3所示为该模型计算的本构方程曲线,其中εC为压缩屈服点的应变,取值为0.0085。
图4为本发明拉伸线性强化弹塑性力学模型示意图,其中,εT为压缩屈服点的应变,取值为0.006,εb为断裂点的应变,取值为0.009。
S105、对添加所述弹塑性本构参数的所述股骨头仿真模拟单元进行预处理,输出仿真结果。
具体地,对所述添加所述弹塑性本构参数的所述股骨头仿真模拟单元再次添加边界条件,经加载后输出仿真结果。
S106、对所述仿真结果进行提取和分析。
具体地,以应变超过屈服点作为股骨头松质骨损坏的判断参数,对所述输出结果进行提取和分析;
其中,所述压缩屈服点的应变为0.0085;所述拉伸屈服点的应变为0.006,拉伸断裂点的应变为0.009。
将股骨头根据密度值等分成多个部分,很好的处理了股骨头非均质材料特性的问题。由于每个部分都采用线性强化弹塑性模型,且各自本构唯一性由该部分股骨头本身决定,使仿真模拟结果和真实情况更加吻合。在结果分析时,用应变值作为参数来判断松质骨的受力响应,更能反应股骨头内部的多孔网状结构,比用应力值作为参数来判断更加符合实际情况,结果也更加稳定。
至此,本领域技术人员应认识到,虽然本文已详尽示出和描述了本发明的多个示例性实施例,但是,在不脱离本发明精神和范围的情况下,仍可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其他变型或修改。因此,本发明的范围应被理解和认定为覆盖了所有这些其他变型或修改。
Claims (6)
1.一种股骨头力学仿真模型处理方法,其特征在于,包括:
根据灰度值对计算机断层扫描CT图像提取股骨头几何模型;
将所述股骨头几何模型进行网格划分,得到与所述股骨头几何模型对应股骨头仿真模拟单元,将灰度值转化成密度值后,进而得到仿真模型单元数量百分比-密度值关系曲线;
将所述股骨头仿真模拟单元划分为多个部分,其中,划分规则为股骨头仿真模型单元按密度值均等划分后的单元数量百分比-密度值所形成的柱体总面积与所述曲线所形成的总面积的偏差不大于5%;
确定各个部分股骨头仿真模拟单元的弹塑性本构参数;
对添加所述弹塑性本构参数的所述股骨头仿真模拟单元进行预处理,输出仿真结果;
对所述仿真结果进行提取和分析。
2.根据权利要求1所述的方法,其特征在于,所述对添加所述弹塑性本构参数的所述股骨头仿真模拟单元进行预处理,输出仿真结果,包括:
对所述添加所述弹塑性本构参数的所述股骨头仿真模拟单元再次添加边界条件,经加载后输出仿真结果。
3.根据权利要求1所述的方法,其特征在于,所述将所述股骨头仿真模拟单元根据所述划分规则划分为多个部分,包括:
采用密度值将所述股骨头仿真模拟单元均等划分为多个部分。
4.根据权利要求3所述的方法,其特征在于,所述划分规则还包括:所述股骨头仿真模拟单元的总数量为N,划分部分数量为M,每个部分的单元数为n,每个部分等分区间的数量取2Mn/N的向上整数。
5.根据权利要求1所述的方法,其特征在于,所述确定各部分股骨头仿真模拟单元的弹塑性本构参数,包括:
对所述股骨头仿真模拟单元采用压缩线性强化弹塑性本构和拉伸线性强化弹塑性本构。
6.根据权利要求1所述的方法,其特征在于,所述对所述仿真结果进行提取和分析,包括:
以应变超过屈服点的应变作为股骨头松质骨损坏的判断参数,对所述输出结果进行提取和分析;
其中,压缩屈服点的应变为0.0085,拉伸屈服点的应变为0.006,拉伸断裂点的应变为0.009。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810062242.1A CN108281193B (zh) | 2018-01-23 | 2018-01-23 | 一种股骨头力学仿真模拟模型处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810062242.1A CN108281193B (zh) | 2018-01-23 | 2018-01-23 | 一种股骨头力学仿真模拟模型处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108281193A CN108281193A (zh) | 2018-07-13 |
CN108281193B true CN108281193B (zh) | 2021-06-25 |
Family
ID=62804578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810062242.1A Active CN108281193B (zh) | 2018-01-23 | 2018-01-23 | 一种股骨头力学仿真模拟模型处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108281193B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110833472B (zh) * | 2019-11-22 | 2020-12-11 | 吉林大学 | 基于3d打印的个体化定制膝关节仿生假体的制作方法 |
CN112137719B (zh) * | 2020-09-16 | 2021-07-06 | 中国科学院力学研究所 | 一种基于ct影像的3d打印仿真骨的制造方法 |
CN115131276B (zh) * | 2022-03-03 | 2023-07-25 | 中国人民解放军总医院第四医学中心 | 骨骼受力分布的获取方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104462636A (zh) * | 2014-09-30 | 2015-03-25 | 西安科技大学 | 基于伞状股骨头支撑器的坏死股骨头修复模型的建模方法 |
CN106372349A (zh) * | 2016-09-09 | 2017-02-01 | 吉林大学 | 构建乘员腿部快速参数化模型设计方法及系统 |
-
2018
- 2018-01-23 CN CN201810062242.1A patent/CN108281193B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104462636A (zh) * | 2014-09-30 | 2015-03-25 | 西安科技大学 | 基于伞状股骨头支撑器的坏死股骨头修复模型的建模方法 |
CN106372349A (zh) * | 2016-09-09 | 2017-02-01 | 吉林大学 | 构建乘员腿部快速参数化模型设计方法及系统 |
Non-Patent Citations (3)
Title |
---|
Development of Femoral Head Interior Supporting Device and 3D Finite Element Analysis of its Application in the Treatment of Femoral Head Avascular Necrosis;Dongmin Xiao et al.;《Medical Science Monitor》;20150526;全文 * |
Mechanics of femoral head osteonecrosis using three-dimensional finite element method;Jin-Won Yang et al.;《Archives of Orthopaedic and Trauma Surgery》;20010905;全文 * |
Prosthetic Modeling of the Femoral Head Based on Mesh Segmentation;Bin Liu et al.;《2011 International Coriference on Electronics and Optoelectronics》;20111231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108281193A (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108281193B (zh) | 一种股骨头力学仿真模拟模型处理方法 | |
Chevalier et al. | A patient-specific finite element methodology to predict damage accumulation in vertebral bodies under axial compression, sagittal flexion and combined loads | |
Hambli | A quasi-brittle continuum damage finite element model of the human proximal femur based on element deletion | |
DeSilva et al. | A comparative study of the trabecular bony architecture of the talus in humans, non-human primates, and Australopithecus | |
Gocha et al. | Spatial variation in osteon population density at the human femoral midshaft: histomorphometric adaptations to habitual load environment | |
Charlebois et al. | A nonlocal constitutive model for trabecular bone softening in compression | |
Fitton et al. | The Impact of Simplifications on the Performance of a Finite Element Model of a M acaca fascicularis Cranium | |
Young et al. | Cranial biomechanics of Diplodocus (Dinosauria, Sauropoda): testing hypotheses of feeding behaviour in an extinct megaherbivore | |
Miszkiewicz | Investigating histomorphometric relationships at the human femoral midshaft in a biomechanical context | |
Fortuny et al. | Skull mechanics and the evolutionary patterns of the otic notch closure in capitosaurs (Amphibia: Temnospondyli) | |
Hambli et al. | Integrated remodeling-to-fracture finite element model of human proximal femur behavior | |
DeSantis et al. | Assessing niche conservatism using a multiproxy approach: dietary ecology of extinct and extant spotted hyenas | |
Liu et al. | Decoupled effects of bone mass, microarchitecture and tissue property on the mechanical deterioration of osteoporotic bones | |
Pitfield et al. | Cortical histomorphometry of the human humerus during ontogeny | |
Bettamer et al. | Using visual image measurements to validate a novel finite element model of crack propagation and fracture patterns of proximal femur | |
Chamrad et al. | On the level of computational model of a human skull: a comparative study | |
Be'ery-Lipperman et al. | A method of quantification of stress shielding in the proximal femur using hierarchical computational modeling | |
Lee et al. | Modeling and simulation of osteoporosis and fracture of trabecular bone by meshless method | |
Mirzaei et al. | Finite element prediction and experimental verification of the failure pattern of proximal femur using quantitative computed tomography images | |
Roque et al. | Tortuosity and elasticity study of distal radius trabecular bone | |
Hosseinitabatabaei et al. | Introduction of maximum stress parameter for the evaluation of stress shielding around orthopedic screws in the presence of bone remodeling process | |
Stefanek et al. | Comparison of linear and nonlinear stepwise μFE displacement predictions to digital volume correlation measurements of trabecular bone biopsies | |
Akhlaghi et al. | Investigation on primary stability of dental implants through considering peri-implant bone damage, caused by small and large deformations: A validated non-linear micro finite element study | |
Coiner‐Collier et al. | Trabecular anisotropy in the primate mandibular condyle is associated with dietary toughness | |
Sahli et al. | Structural parameters determining the strength of the porcine vertebral body affected by tumours |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |