CN103440423A - 一种用于测量生理载荷作用下关节软骨力学性能的模型 - Google Patents
一种用于测量生理载荷作用下关节软骨力学性能的模型 Download PDFInfo
- Publication number
- CN103440423A CN103440423A CN2013103948938A CN201310394893A CN103440423A CN 103440423 A CN103440423 A CN 103440423A CN 2013103948938 A CN2013103948938 A CN 2013103948938A CN 201310394893 A CN201310394893 A CN 201310394893A CN 103440423 A CN103440423 A CN 103440423A
- Authority
- CN
- China
- Prior art keywords
- cartilage
- model
- spring
- rolling
- pressure head
- 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.)
- Granted
Links
Images
Landscapes
- Instructional Devices (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种用于测量生理载荷作用下关节软骨力学性能的模型,包括滑动压缩模型和滚动压缩模型,该模型将软骨分为不同层区进行滑动压缩和滚动压缩仿真,考虑了软骨的各向异性特征及实际的生理载荷条件。同时结合有限元软件ABAQUS中的弹簧单元模拟胶原纤维的拉伸性能,考虑了组成成分的作用;该模型用于测量生理载荷下关节软骨的力学性能,包括滑动压缩测量和滚动压缩测量。本发明的优点是:该有限元仿真方法的运用更加接近关节软骨真实受力环境,更加接近软骨组成成分的功能,因此能准确模拟生理载荷下软骨的微观力学性能。
Description
技术领域
本发明主要涉及生物医学工程领域,特别是一种用于测量生理载荷作用下关节软骨力学性能的模型。
背景技术
关节软骨是人体膝关节活动的重要组成部分,起着缓冲震荡、传递载荷等不可替代的作用。目前由于实验条件的限制,对生理载荷作用下软骨的微观结构和组分变化及真实受力状态的测试很难实现,因此许多学者采用有限元仿真的方法来探究软骨的力学性能。Mow 等人最早在混合物理论的两相多孔介质理论基础上建立了关节软骨的两相多孔模型。该模型由两种内在不可压缩、不可混容的固体基质和液体组成,其中固体基质视为线弹性固体,液体视为理想流体。该模型在一定程度上比较准确的描述了软骨的力学性能。Wilson等人建立了纤维增强、多孔粘弹性关节软骨力学模型,包括各向异性粘弹性胶原蛋白结构和依据蛋白多糖密度变化的膨胀性能,并对软骨进行了应力松弛分析。翟文杰等人考虑关节软骨基质的固相和孔隙液相及胶原纤维增强相,对软骨进行了应力松弛模拟分析。这些研究者基本上是将软骨看作杨氏模量相同的均匀质体,进行一些单纯的压缩模拟或者压痕仿真,没有考虑软骨的各向异性特征及实际的生理载荷条件。正如我们所知,关节软骨精细的结构按照胶原纤维分布排列方式可以分为三层:浅表层、中间层和深层,关节软骨的非线性和各向异性特性不仅体现在它的结构和组成,而且体现于它的生物力学性质中。在人类的日常生活活动中,软骨所承受的主要生理载荷是滚压和滑动载荷,即滚压和滑动载荷是关节软骨生长、发育和维持正常生理功能合适的生理载荷。因此,研究此载荷下软骨不同层区的力学行为是非常重要的,基于此,我们提出一种将滚压(包括滑动)载荷作用于功能化关节软骨构建的有限元方法,并对该载荷下软骨的微观力学性能进行了模拟仿真。
发明内容
本发明的目的是针对上述存在问题,提供一种用于功能化关节软骨的滚压和滑动载荷模型及其测量方法,该测量方法建立的软骨模型更加接近关节软骨真实受力环境,进而能准确模拟和测量生理载荷下软骨的微观力学性能。
本发明的技术方案:
一种用于测量生理载荷作用下关节软骨力学性能的模型,包括滑动压缩模型和滚动压缩模型,其中滑动压缩模型包括滑动压头、软骨模型、分层结构,滑动压头为半圆形,可对软骨施加一定的压缩量并在软骨表面连续滑动,软骨模型由不同层分层结构组成;滚动压缩模型包括滚动压头、软骨模型、分层结构,滚动压头为圆形,可对软骨施加一定的压缩量并在软骨表面连续滚动,软骨模型由不同层分层结构组成;软骨为二维的软骨模型,通过使用ABAQUS中的Shell建立并分割为3-9层,同时根据浅表层、中间层及深层各占15%、55%、30%的比例自上而下定义三组不同层区,每层赋予不同的材料属性,每层的杨氏模量和泊松比取不同的值,其中泊松比为0-0.499,杨氏模量自浅表层到深层递增;采用八节点孔压单元模拟固体基质,节点之间连接有弹簧单元用于模拟胶原纤维的拉伸性能;滑动压头采用ABAQUS中Part下的二维线(Wire)选项建立,类型为解析刚体(Analytical rigid);滚动压缩模型的建立与滑动压缩模型相同,不同之处是滚动压头采用ABAQUS中Part下的二维壳(Shell)选项建立,类型为可变形(Deformable),定义其刚度为210GPa或定义为一种圆形离散刚体线(discrete rigid wire)。
一种所述用于测量生理载荷作用下关节软骨力学性能的模型的应用,用于测量生理载荷下关节软骨的力学性能,包括滑动压缩测量和滚动压缩测量,方法如下:
1)滑动压缩测量
在ABAQUS软件中建立滑动压头和软骨模型,设定滑动压头的下移距离即压缩软骨的量为0-30%的任意值,设定其向右滑动的速度为0-35cm/s的任意值;定义软骨中各层的杨氏模量、泊松比,定义软骨的边界条件,软骨模型左侧界面和右侧界面的液体可自由渗透,底部垂直方向的位移为零;运用子程序定义软骨渗透率,划分网格,径向网格划分比例为0.5:0.25:0.25,也可根据实际需要自定义;单元为八节点孔压单元;在CAE界面上的节点之间建立弹簧单元,所用弹簧单元类型为SpringA,在INP文件中找到Spring, elset=SpringA-spring字段,将其修改为*Spring, nonlinear, elset=SpringA-spring,并在下一行输入力与位移的一一对应数值来确定其非线性,数值范围为(0.04712,0.0005)到(9.651,0.008),最后在ABAQUS/Command中提交分析;
2)滚动压缩测量
在ABAQUS软件中建立滚动压头和软骨模型,定义滚动压头为圆形实体面或圆形离散刚体线,设定其压缩软骨的量为0-30%的任意值,及向右滚动的速度为0-20cm/s的任意值,并根据公式 计算获得绕圆心的角速度;定义软骨中各层的杨氏模量、泊松比,定义软骨的边界条件,软骨模型左侧界面和右侧界面的液体可自由渗透,底部垂直方向的位移为零;运用子程序定义软骨渗透率,划分网格,单元为八节点孔压单元;在CAE界面上的节点之间建立弹簧单元,所用弹簧单元类型为SpringA,在INP文件中找到Spring, elset=SpringA-spring字段,将其修改为*Spring, nonlinear, elset=SpringA-spring,并在下一行输入力与位移的一一对应数值来确定其非线性,数值范围为(0.04712,0.0005)到(9.651,0.008),最后在ABAQUS/Command中提交分析。
本发明的优点是:该方法依据软骨的各向异性结构和液固两相特征,建立软骨模型,并模拟关节软骨在真实生理状态下所受的滑动压缩载荷和滚动压缩载荷,获得其生理载荷作用下的微观力学性能,采用该方法得到的仿真结果更加接近软骨的真实力学性能。
【附图说明】
图1为关节软骨滑动压缩模型示意图。
图2为关节软骨滚动压缩模型示意图。
图中:1a.滑动压头 1b.滚动压头 2.软骨模型 3.分层结构
【具体实施方式】
实施例:
一种用于测量生理载荷作用下关节软骨力学性能的模型,包括滑动压缩模型和滚动压缩模型,其中滑动压缩模型,如图1所示,包括滑动压头1a、软骨模型2、分层结构3,滑动压头1a为半圆形,可施加一定的压缩量压缩软骨2,软骨模型2由不同层分层结构3组成;滚动压缩模型,如图2所示,包括滚动压头1b、软骨模型2、分层结构3,滚动压头1b为圆形,可施加一定的压缩量压缩软骨2,软骨模型2由不同层分层结构3组成;软骨模型为二维的软骨模型,通过使用ABAQUS中的shell建立并分割为3-9层,同时根据浅表层、中间层及深层各占15%、55%、30%的比例自上而下定义三组不同层区,每层赋予不同的材料属性,每层的杨氏模量和泊松比取不同的值,其中泊松比为0-0.499,杨氏模量自浅表层到深层递增;采用八节点孔压单元模拟固体基质,节点之间连接有弹簧单元用于模拟胶原纤维的拉伸性能;滑动压头采用ABAQUS中Part下的二维线(Wire)选项建立,类型为解析刚体(Analytical rigid);滚动压缩模型的建立与滑动压缩模型相同,不同之处是滚动压头采用ABAQUS中Part下的二维壳(Shell)选项建立,类型为可变形(Deformable),定义其刚度为210GPa或定义为一种圆形离散刚体线(discrete rigid wire)。
所述关节软骨生理载荷模型用于测量生理载荷下关节软骨的力学性能,包括滑动压缩测量和滚动压缩测量,方法如下:
1)滑动压缩仿真:在ABAQUS软件中建立滑动压头和软骨模型。设定滑动压头的下移距离,即定义压缩软骨的量10%,并设定其向右滑动的速度2cm/s;建立7层模型,其中第一层为浅表层,第二层到第五层为中间层,第六到第七层为深层,杨氏模量自上而下分别取0.079MPa、0.17MPa、0.27 MPa、0.55 MPa、0.58 MPa、0.73 MPa、1.14 MPa。定义泊松比均为0.499;定义软骨的边界条件:模型左侧边界液体可自由渗透,底部垂直方向位移为零,右侧边界液体可自由渗透;运用子程序定义软骨渗透率;划分网格,单元为八节点孔压单元;在CAE界面上的节点之间建立弹簧单元,所用弹簧单元类型为SpringA,在INP文件中找到Spring, elset=SpringA-spring字段,将其修改为*Spring, nonlinear, elset=SpringA-spring,并在下一行输入力与位移的一一对应数值来确定其非线性,数值范围为(0.04712,0.0005)到(9.651,0.008);最后在ABAQUS/Command中提交分析。
2)滚动压缩仿真:在ABAQUS软件中建立滚动压头和软骨模型,定义滚动压头为圆形实体面,半径为1.5cm,设定压缩软骨的量10%,设定其向右滚动的速度10cm/s,角速度为6.67rad/s;建立7层模型,其中第一层为浅表层,第二层到第五层为中间层,第六到第七层为深层,杨氏模量自上而下分别取0.079MPa、0.17MPa、0.27 MPa、0.55 MPa、0.58 MPa、0.73MPa、1.14 MPa。定义泊松比均为0.499;定义软骨的边界条件:模型左侧边界液体可自由渗透,底部垂直方向位移为零,右侧边界液体可自由渗透;运用子程序定义软骨渗透率;划分网格,单元为八节点孔压单元;在CAE界面上的节点之间建立弹簧单元,所用弹簧单元类型为SpringA,在INP文件中找到Spring, elset=SpringA-spring字段,将其修改为*Spring, nonlinear, elset=SpringA-spring,并在下一行输入力与位移的一一对应数值来确定其非线性,数值范围为(0.04712,0.0005)到(9.651,0.008),最后在ABAQUS/Command中提交分析。
实践表明:该仿真方法效果良好,能真实的反映软骨在滑动和滚压载荷作用下的力学性能。
Claims (2)
1.一种用于测量生理载荷作用下关节软骨力学性能的模型,其特征在于:包括滑动压缩模型和滚动压缩模型,其中滑动压缩模型包括滑动压头、软骨模型、分层结构,滑动压头为半圆形,可对软骨施加一定的压缩量并在软骨表面连续滑动,软骨模型由不同层分层结构组成;滚动压缩模型包括滚动压头、软骨模型、分层结构,滚动压头为圆形,可对软骨施加一定的压缩量并在软骨表面连续滚动,软骨模型由不同层分层结构组成;软骨为二维的软骨模型,通过使用ABAQUS中的Shell建立并分割为3-9层,同时根据浅表层、中间层及深层各占15%、55%、30%的比例自上而下定义三组不同层区,每层赋予不同的材料属性,每层的杨氏模量和泊松比取不同的值,其中泊松比为0-0.499,杨氏模量自浅表层到深层递增;采用八节点孔压单元模拟固体基质,节点之间连接有弹簧单元用于模拟胶原纤维的拉伸性能;滑动压头采用ABAQUS中Part下的二维线选项建立,类型为解析刚体;滚动压缩模型的建立与滑动压缩模型相同,不同之处是滚动压头采用ABAQUS中Part下的二维壳选项建立,类型为可变形,定义其刚度为210GPa或定义为一种圆形离散刚体线。
2.一种权利要求1所述用于测量生理载荷作用下关节软骨力学性能的模型的应用,其特征在于:用于测量生理载荷下关节软骨的力学性能,包括滑动压缩测量和滚动压缩测量,方法如下:
1)滑动压缩测量
在ABAQUS软件中建立滑动压头和软骨模型,设定滑动压头的下移距离即压缩软骨的量为0-30%的任意值,设定其向右滑动的速度为0-35cm/s的任意值;定义软骨中各层的杨氏模量、泊松比,定义软骨的边界条件,软骨模型左侧界面和右侧界面的液体可自由渗透,底部垂直方向的位移为零;运用子程序定义软骨渗透率,划分网格,径向网格划分比例为0.5:0.25:0.25,也可根据实际需要自定义;单元为八节点孔压单元;在CAE界面上的节点之间建立弹簧单元,所用弹簧单元类型为SpringA,在INP文件中找到Spring, elset=SpringA-spring字段,将其修改为*Spring, nonlinear, elset=SpringA-spring,并在下一行输入力与位移的一一对应数值来确定其非线性,数值范围为(0.04712,0.0005)到(9.651,0.008),最后在ABAQUS/Command中提交分析;
2)滚动压缩测量
在ABAQUS软件中建立滚动压头和软骨模型,定义滚动压头为圆形实体面或圆形离散刚体线,设定其压缩软骨的量为0-30%的任意值,及向右滚动的速度为0-20cm/s的任意值,并根据公式 计算获得绕圆心的角速度;定义软骨中各层的杨氏模量、泊松比,定义软骨的边界条件,软骨模型左侧界面和右侧界面的液体可自由渗透,底部垂直方向的位移为零;运用子程序定义软骨渗透率,划分网格,单元为八节点孔压单元;在CAE界面上的节点之间建立弹簧单元,所用弹簧单元类型为SpringA,在INP文件中找到Spring, elset=SpringA-spring字段,将其修改为*Spring, nonlinear, elset=SpringA-spring,并在下一行输入力与位移的一一对应数值来确定其非线性,数值范围为(0.04712,0.0005)到(9.651,0.008),最后在ABAQUS/Command中提交分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310394893.8A CN103440423B (zh) | 2013-09-04 | 2013-09-04 | 一种用于测量生理载荷作用下关节软骨力学性能的模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310394893.8A CN103440423B (zh) | 2013-09-04 | 2013-09-04 | 一种用于测量生理载荷作用下关节软骨力学性能的模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103440423A true CN103440423A (zh) | 2013-12-11 |
CN103440423B CN103440423B (zh) | 2016-05-11 |
Family
ID=49694117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310394893.8A Expired - Fee Related CN103440423B (zh) | 2013-09-04 | 2013-09-04 | 一种用于测量生理载荷作用下关节软骨力学性能的模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103440423B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055848A (zh) * | 2016-07-14 | 2016-10-26 | 哈尔滨理工大学 | 一种基于微观结构参数个体化关节软骨仿真方法 |
CN107808037A (zh) * | 2017-10-10 | 2018-03-16 | 哈尔滨理工大学 | 一种基于纤维方向的关节软骨的建模计算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070254041A1 (en) * | 2006-05-01 | 2007-11-01 | Drapeau Susan J | Demineralized bone matrix devices |
US20080026362A1 (en) * | 2006-07-20 | 2008-01-31 | Kaohsiung Medical University | Three-dimensional culture containing human articular chondrocytes with induced terminal differentiation changes, and preparation process and uses of the same |
CN101250501A (zh) * | 2008-04-11 | 2008-08-27 | 北京航空航天大学 | 模拟膝关节物理环境培养软骨的方法及生物反应器 |
CN102353589A (zh) * | 2011-07-11 | 2012-02-15 | 天津理工大学 | 一种在软骨力学性能测试中运用数字图像相关技术的方法 |
-
2013
- 2013-09-04 CN CN201310394893.8A patent/CN103440423B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070254041A1 (en) * | 2006-05-01 | 2007-11-01 | Drapeau Susan J | Demineralized bone matrix devices |
US20080026362A1 (en) * | 2006-07-20 | 2008-01-31 | Kaohsiung Medical University | Three-dimensional culture containing human articular chondrocytes with induced terminal differentiation changes, and preparation process and uses of the same |
CN101250501A (zh) * | 2008-04-11 | 2008-08-27 | 北京航空航天大学 | 模拟膝关节物理环境培养软骨的方法及生物反应器 |
CN102353589A (zh) * | 2011-07-11 | 2012-02-15 | 天津理工大学 | 一种在软骨力学性能测试中运用数字图像相关技术的方法 |
Non-Patent Citations (3)
Title |
---|
张玉峰 等: "《有限元分析系统ABAQUS中的特征技术》", 《工程图学学报》, no. 5, 31 October 2006 (2006-10-31) * |
王玉主: "《关节软骨的力学性能测试及有限元仿真》", 《中国优秀硕士学位论文全文数据库(医药卫生科技辑)》, no. 10, 15 October 2012 (2012-10-15) * |
翟文杰,翟中勇: "《软骨组织无约束压缩的有限元分析》", 《医用生物力学》, vol. 27, no. 6, 31 December 2012 (2012-12-31) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055848A (zh) * | 2016-07-14 | 2016-10-26 | 哈尔滨理工大学 | 一种基于微观结构参数个体化关节软骨仿真方法 |
CN106055848B (zh) * | 2016-07-14 | 2017-10-13 | 哈尔滨理工大学 | 一种基于微观结构参数个体化关节软骨仿真方法 |
CN107808037A (zh) * | 2017-10-10 | 2018-03-16 | 哈尔滨理工大学 | 一种基于纤维方向的关节软骨的建模计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103440423B (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Hybrid porous micro structured finger skin inspired self-powered electronic skin system for pressure sensing and sliding detection | |
Jin et al. | On the contact behavior of micro-/nano-structured interface used in vertical-contact-mode triboelectric nanogenerators | |
Ehlers et al. | A linear viscoelastic biphasic model for soft tissues based on the theory of porous media | |
Mow et al. | Biphasic creep and stress relaxation of articular cartilage in compression: theory and experiments | |
CN102262699B (zh) | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 | |
Mow et al. | The nonlinear interaction between cartilage deformation and interstitial fluid flow | |
Lin | Alternative form of standard linear solid model for characterizing stress relaxation and creep: Including a novel parameter for quantifying the ratio of fluids to solids of a viscoelastic solid | |
Aernouts et al. | Quantification of tympanic membrane elasticity parameters from in situ point indentation measurements: validation and preliminary study | |
Müller et al. | A hyperelastic model for simulating cells in flow | |
Farukh et al. | Numerical modelling of damage initiation in low-density thermally bonded nonwovens | |
Crawford et al. | Compression-dependent permeability measurement for random soft porous media and its implications to lift generation | |
Recuerda et al. | Influence of experimental protocols on the mechanical properties of the intervertebral disc in unconfined compression | |
Barabadi et al. | On the characterization of lifting forces during the rapid compaction of deformable porous media | |
CN103440423A (zh) | 一种用于测量生理载荷作用下关节软骨力学性能的模型 | |
Dallon et al. | A force based model of individual cell migration with discrete attachment sites and random switching terms | |
Pagano et al. | A microscale-based model for small-strain stiffness in unsaturated granular geomaterials | |
Maleki et al. | Effect of structural distortions on articular cartilage permeability under large deformations | |
Taffetani et al. | A quantitative interpretation of the response of articular cartilage to atomic force microscopy-based dynamic nanoindentation tests | |
Benalla et al. | Analytical basis for the determination of the lacunar–canalicular permeability of bone using cyclic loading | |
CN103206203B (zh) | 油井单一射孔出砂的分析方法 | |
CN107808037A (zh) | 一种基于纤维方向的关节软骨的建模计算方法 | |
Sundara-Rajan et al. | Capacitive sensing of interfacial forces in prosthesis | |
Shegaf et al. | Cartilage biomechanical response differs under physiological biaxial loads and uniaxial cyclic compression | |
Yerrabelli et al. | Modeling the mechanics of fibrous-porous scaffolds for skeletal muscle regeneration | |
Tsai | Viscosity effect on consolidation of poroelastic soil due to groundwater table depression |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160511 Termination date: 20160904 |