首页 / 专利库 / 物理 / 应力分布 / 一种骨表面残余应力分布测试方法

一种骨表面残余应分布测试方法

阅读:385发布:2020-05-11

专利汇可以提供一种骨表面残余应分布测试方法专利检索,专利查询,专利分析的服务。并且本 发明 公开一种骨表面残余应 力 分布测试方法:步骤1建立了骨不同深度 应力 与应变的关系;步骤2建立了骨的应力与深度的关系;步骤3通过等式变换消除了实验最终要测的未知量,引入了柔度系数这一量,建立了应变与柔度系数和待求系数的关系并通过有限元模型计算出了柔度系数;步骤4通过改进的实验测得了应变值;步骤5通过实验测得了骨这一 各向异性 材料横向和纵向的 弹性模量 ;步骤6通过求出的柔度系数和实验测得的应变解出了待求系数,从而最终求出了骨表面残余应力与深度的关系。本发明的有益效果是:通过独创的实验手段和建模结合的方式,提出了一种新的测量骨表面残余 应力分布 的测试方法,能够较为准确地测量骨表面附近残余应力分布。,下面是一种骨表面残余应分布测试方法专利的具体信息内容。

1.一种骨表面残余应分布测试方法,其特征在于:该方法包含以下六个步骤:
步骤1、建立应力与应变的关系
假设残余应力在z方向是均匀分布的,残余应力能够通过以下公式与应变联系起来式中,ε(ai)是槽深度为ai时测得的应变值,函数G(x,ai)给出了在槽深度为ai时深度x处单位应力的应变响应;
步骤2、建立残余应力与深度的关系
残余应力与深度的关系可以近似为
式中,Aj是待求系数,n是槽深度方向上被划分的单元数量,Uj(x)是单位脉冲,当深度x满足aj-1≤x≤aj时,Uj(x)=1;当x满足x≤aj-1,x≥aj时,Uj(x)=0;
步骤3、求解柔度系数Cij
将式(2)代入式(1),得到
式中Cij表示柔度矩阵公式,它可以表示为
式中,E为弹性模量
从式(1)和式(4)可知,矩阵中的元素Cij代表在aj-1≤x≤aj时,不同增量下,单位应力引起的表面上的应变,即当槽在深度方向有i个增量单元时,在第j个增量单元上施加单位应力,此时表面上的应变响应就是柔度系数Cij;
Cij的值可以通过建立槽的有限元模型,然后逐步施加单位应力计算出来,通过删除单元格来模拟槽深度的扩展,但是这种方法只能获得同一厚度平面的平均应力,而无法得到在同一个厚度沿切口方向的应力分布;因此将每个厚度平面分解为多个相互平行的单元,逐个分别施加单位应力,求解每个单元应力对于表面位移的贡献;通过反解位移传递函数矩阵的形式就能获得残余应力沿切口方向分布的情况,再以逐步递进加深切口的方式,得到残余应力沿厚度方向变化情况,从而组合得到整个切割面的三维残余应力分布;
步骤4、用股骨进行切槽实验测量骨头不同深度处的应变
将牛股骨先进行处理,除去末端和软组织,将其浸泡在5%磷酸盐缓冲溶液中48小时;
然后将其固定在夹具上,将试件连同固定它的夹具在室温下浸入箱中,在试件上方安装切槽装置;对试件的表面进行检查,寻找一个平坦的区域并绘制该区域的平坦度示意图;接下来确定预期要开的槽大致所处的表面区域,这个区域在槽长度内的平均平坦度必须在
0.013毫米内,然后标记这个区域;暂时取出试件,用氰基丙烯酸酯胶粘剂将应变片附着在预期要开的槽附近的表面上,应变片的边缘距槽的边缘1毫米;等待涂层干燥后,将试件放回水箱,调整预期要开的槽的方向以确保其水平,将试件浸没在水箱中24小时以进一步水合,使得试件温度与周围大量水的温度保持平衡;在切槽之前,检查应变读数的热稳定性,在实验持续时间内,应变的波动不超过2με;
确定所要开的槽的长度、位置和深度;使用精确的平移台设置槽的深度;以小于每分钟
10转的速度手动转动钻杆并平移铣刀来切槽;槽经过十步最终被加深到0.61毫米,前六次每次加深0.051毫米,后四次每次加深0.076毫米;在每个切槽步骤后,在沿着槽长度的四个位置测量实际深度,并取这些值的平均值作为槽的实际深度,槽的深度需达到0.0025毫米的精确度;对应于有限元模型中的这个深度的应变值通过测得的应变与深度关系插值得到;
步骤5、计算纵向弹性模量EL和横向弹性模量ET
将来自牛股骨的不同深度层的试件仔细研磨,然后室温下在5%磷酸盐缓冲盐溶液中储存48h;接下来将这些试件放在跨度为8毫米的微型三点弯装置中进行弯曲后测量挠度;
用一个100倍显微镜观察每个试件的中点挠度,中点挠度与中点弯矩有如下关系式中ω表示试件中点的位移,M表示试件中点的弯矩,l表示微型三点弯装置的跨度,E表示弹性模量,I表示试件截面的惯性矩;E的值可以通过式(5)求出来,在牛股骨的横向和纵向分别截取试件进行测试,可以分别求出横向和纵向弹性模量ET和EL;
步骤6、求解骨表面不同深度层的残余应力
将式(3)表示为{{ε}=[C]{A},根据步骤3测得的柔度系数Cij和步骤4测得的应变可以解得式(2)中的系数Aj
{A}=([C]T[C])-1CT{εmeas}   (6)
将Aj代入式(2)即可求出骨表面残余应力与深度的关系。

说明书全文

一种骨表面残余应分布测试方法

【技术领域】

[0001] 本发明属于生物实验力学领域,具体涉及一种骨表面残余应力分布测试方法。【背景技术】
[0002] 在人体的组织和器官当中,残余应力普遍存在并扮演着非常重要的色。组织和器官在其生长、重构和形态变化的过程中均会产生相应的残余应力。而反过来,这些残余应力在组织器官生长、重构和形态变化过程中将起着非常重要的作用。同样,残余应力在骨结构的生长、成形和日常修复中起着至关重要的作用。骨表面残余应力的分布研究对于骨组织成长与重建机理研究对于治疗相应的骨科疾病,研制仿生骨结构至关重要。但是由于目前缺乏表面残余应力的有效测量手段,相关研究才刚刚起步。
[0003] 随着残余应力在部分生物结构中的存在、形成和作用逐渐被人们了解和认识,部分生物组织中残余应力的测量手段和方法也在不断发展。但是,由于生物组织往往具有各向异性,非线性,粘弹性,同时在不同层间性质还会发生变化,几何结构也比工程构件更为复杂;并且生物组织离体后残余应力可能会发生变化,材料预处理方式,选取的材料体积等都可能极大地影响实验结果,这就给传统的残余应力的测试方法带来了极大的挑战。到目前为止,骨头里的残余应力主要通过X射线衍射来进行研究。但是研究表明,辐射剂量和样品中量都会影响这种方法的实验精度。据此,我们提出了一种骨表面残余应力分布测试方法,该方法用一个逐渐加深的槽来研究股骨表面附近不同深度层的残余应力,克服了X射线衍射法的缺陷。【发明内容】
[0004] 本发明的目的是提供一种骨表面残余应力分布测试方法,即在骨表面引入一个逐渐加深的槽,在骨头中,残余应力σ垂直于槽并随着深度x以某种关系变化。这个槽的产生释放了残余应力,使得物体表面产生了垂直于槽的应变ε,这通常可以通过安装在槽位置附近(根据需要也可以安装在对面)的应变片测量出来。测得的应变随深度变化的数值可以和计算模型一起,用于计算残余应力σ随深度的变化规律。
[0005] 本发明为一种骨表面残余应力分布测试方法,它包含以下六个步骤:
[0006] 步骤1、建立应力与应变的关系
[0007] 假设残余应力在z方向是均匀分布的,残余应力能够通过以下公式与应变联系起来
[0008]
[0009] 式中,ε(ai)是槽深度为ai时测得的应变值,函数G(x,ai)给出了在槽深度为ai时深度x处单位应力的应变响应。
[0010] 步骤2、建立残余应力与深度的关系
[0011] 残余应力与深度的关系可以近似为
[0012]
[0013] 式中,Aj是待求系数,n是槽深度方向上被划分的单元数量,Uj(x)是单位脉冲,当深度x满足aj-1≤x≤aj时,Uj(x)=1;当x满足x≤aj-1,x≥aj时,Uj(x)=0。
[0014] 步骤3、求解柔度系数Cij
[0015] 将式(2)代入式(1),得到
[0016]
[0017] 式中Cij表示柔度矩阵公式,它可以表示为
[0018]
[0019] 式中,E为弹性模量
[0020] 从式(1)和式(4)可知,矩阵中的元素Cij代表在aj-1≤x≤aj时,不同增量下,单位应力引起的表面上的应变,即当槽在深度方向有i个增量单元时,在第j个增量单元上施加单位应力,此时表面上的应变响应就是柔度系数Cij。
[0021] 如图2所示,Cij的值可以通过建立槽的有限元模型,然后逐步施加单位应力计算出来,通过删除单元格来模拟槽深度的扩展。但是这种方法只能获得同一厚度平面的平均应力,而无法得到在同一个厚度沿切口方向的应力分布。因此本方法将每个厚度平面分解为多个相互平行的单元,逐个分别施加单位应力,求解每个单元应力对于表面位移的贡献。如图3所示,通过反解位移传递函数矩阵的形式就能获得残余应力沿切口方向分布的情况,再以逐步递进加深切口的方式,得到残余应力沿厚度方向变化情况,从而组合得到整个切割面的三维残余应力分布。
[0022] 步骤4、用骨头进行切槽实验测量骨头不同深度处的应变
[0023] 将牛股骨先进行处理,除去末端和软组织,将其浸泡在5%磷酸盐缓冲溶液(PBS)中48小时;然后将其固定在夹具上,将试件连同固定它的夹具在室温下浸入水箱中,在试件上方安装切槽装置;对试件的表面进行检查,寻找一个平坦的区域并绘制该区域的平坦度示意图;接下来确定预期要开的槽大致所处的表面区域,这个区域在槽长度内的平均平坦度必须在0.013毫米内,然后标记这个区域;暂时取出试件,用氰基丙烯酸酯胶粘剂将应变片附着在预期要开的槽附近的表面上,应变片的边缘距槽的边缘1毫米;等待涂层干燥后,将试件放回水箱,调整预期要开的槽的方向以确保其水平,将试件浸没在水箱中24小时以进一步水合,使得试件温度与周围大量水的温度保持平衡;在切槽之前,检查应变读数的热稳定性,在实验持续时间内,应变的波动不超过2με。
[0024] 确定所要开的槽的长度、位置和深度。使用精确的平移台设置槽的深度。以小于每分钟10转的速度手动转动钻杆并平移铣刀来切槽。槽经过十步最终被加深到0.61毫米,前六次每次加深0.051毫米,后四次每次加深0.076毫米。在每个切槽步骤后,在沿着槽长度的四个位置测量实际深度,并取这些值的平均值作为槽的实际深度,槽的深度需达到0.0025毫米的精确度。对应于有限元模型中的这个深度的应变值通过测得的应变与深度关系插值得到。
[0025] 步骤5、计算纵向弹性模量EL和横向弹性模量ET
[0026] 将来自骨头的不同深度层的试件仔细研磨,然后室温下在5%磷酸盐缓冲盐溶液中储存48h。接下来将这些试件放在跨度为8毫米的微型三点弯装置中进行弯曲后测量挠度。用一个100倍显微镜观察每个试件的中点挠度,中点挠度与中点弯矩有如下关系[0027]
[0028] 式中ω表示试件中点的位移,M表示试件中点的弯矩,l表示微型三点弯装置的跨度,E表示弹性模量,I表示试件截面的惯性矩。E的值可以通过式(5)求出来,在骨头的横向和纵向分别截取试件进行测试,可以分别求出横向和纵向弹性模量ET和EL。
[0029] 步骤6、求解骨表面不同深度层的残余应力
[0030] 将式(3)表示为{{ε}=[C]{A},根据步骤3测得的柔度系数Cij和步骤4测得的应变可以解得式(2)中的系数Aj
[0031] {A}=([C]T[C])-1CT{εmeas}   (6)
[0032] 将Aj代入式(2)即可求出骨表面残余应力与深度的关系。
[0033] 通过上述六个步骤,最终达到了骨表面残余应力分布测试的目的。其中步骤1建立了骨不同深度应力与应变的关系;步骤2建立了骨的应力与深度的关系;步骤3通过等式变换消除了实验最终要测的未知量,引入了柔度系数这一量,建立了应变与柔度系数和待求系数的关系并通过有限元模型计算出了柔度系数;步骤4通过改进的实验测得了应变值;步骤5通过实验测得了骨这一各向异性材料横向和纵向的弹性模量;步骤6通过求出的柔度系数和实验测得的应变解出了待求系数,从而最终求出了骨表面残余应力与深度的关系。
[0034] 本发明一种骨表面残余应力分布测试方法的有益效果是:通过独创的实验手段和建模结合的方式,提出了一种新的测量骨表面残余应力分布的测试方法,能够较为准确地测量骨表面附近残余应力分布。【附图说明】
[0035] 图1为本方法流程图
[0036] 图2为在不同深度施加单位应力以求出柔度矩阵Cij示意图。
[0037] 图3为在深度和切口方向同时分割单元然后施加单位应力以求出柔度矩阵Cij示意图。
[0038] 图4为进行切槽实验的装置示意图。【具体实施方式】
[0039] 下面结合附图对本发明作详细说明。
[0040] 本发明提出了一种骨表面残余应力分布测试方法,如图1所示,它包括以下六个步骤:
[0041] 步骤一:建立应力与应变的关系
[0042] 假设残余应力在z方向是均匀分布的,残余应力与应变的关系一般表示为[0043]
[0044] 式中,ε(ai)是槽深度为ai时测得的应变值,函数G(x,ai)给出了在槽深度为ai时深度x处单位应力的应变响应。
[0045] 步骤二:建立应力与深度的关系
[0046] 残余应力与深度的关系可以近似为
[0047]
[0048] 式中,Aj是待求系数,Uj(x)是单位脉冲,当深度x满足aj-1≤x≤aj时,Uj(x)=1;当x满足x≤aj-1,x≥aj时,Uj(x)=0。
[0049] 步骤三:用有限元模型求解柔度系数Cij
[0050] 将式(8)代入式(7),得到
[0051]
[0052] 式中Cij表示柔度矩阵公式,它可以表示为
[0053]
[0054] 式中,E为弹性模量。
[0055] 采用有限元程序建立一个有限元模型来测定逐渐加深的狭槽的柔度系数Cij。槽经过十步最终被加深到最终深度。如图3所示,将切口沿着槽的深度方向和切口方向分割成多个相互平行的单元,柔度系数Cij的值通过在每个单元上一步一步施加一个单位应力计算出来。通过删除单元格来模拟槽深度的扩展,用位移数据来计算应变。
[0056] 步骤四:用骨头进行切槽实验测量应变
[0057] 将牛股骨进行处理,除去末端和软组织,然后将其在5%磷酸盐缓冲溶液(PBS)中浸泡48小时。每个试件长度在125到150毫米之间,横截面宽度为50毫米到75毫米。接下来,将试件固定在定制的夹具上,使其稳定夹持试件并且使切槽能够在两个方向精确地调整倾斜角度。然后,将试件连同固定它的夹具在室温下浸入水箱,在试件上方安装一个定制的切槽装置。这个装置是由精密移动平台和导管组成,精密移动平台用来控制导管在x、y、z三个方向上精密移动,导管用来安装显微镜和钻头等其他装置。将显微镜管以紧密配合的方式滑入装置的导管中,并对试件的表面进行检查以找到一个平坦的区域。接下来,取出显微镜管,将位移探头滑入导管中,用来绘制该区域的平坦度示意图。通过调整试件的位置和倾斜度,可以确定预期要开的槽大致所处的表面区域,使得这个区域在槽长度内的平均平坦度在0.013毫米内,然后用防水墨水标记这个区域。暂时从水箱中取出试件,用氰基丙烯酸酯胶粘剂将应变片附着在预期要开的槽附近的表面上,应变片的边缘距槽的边缘1毫米。将应变片及其末端覆盖一层聚酯涂层用以防水,并将切槽的区域暂时遮挡以防止被涂层覆盖。等待20分钟涂层干燥后,将试件放回水箱,调整预期要开的槽的方向以确保其水平。将试件浸没在水箱中24小时以进一步水合,使得试件温度与周围大量水的温度保持平衡。在切槽之前,检查应变读数的热稳定性。在实验持续时间内,应变的波动不超过2με。
[0058] 确定所要开的槽的长度、位置和深度。切槽是通过将带有立铣刀的钻杆滑动到图4所示的导管中来进行的。使用精确的平移台(y方向)来设置槽的深度。通过以小于每分钟10转的速度手动转动钻杆并使用另一个平移台在z方向逐渐平移立铣刀来切槽。槽经过十步最终被加深到0.61毫米,前六次每次加深0.051毫米,后四次每次加深0.076毫米。经过预实验验证,采用这个数值进行实验效果最好,当然,采用其他数值也能得出结果。实验中每一次加深过后槽的实际深度并不与有限元模型中一模一样。因此在每个切槽步骤后,在沿着槽长度的四个位置测量实际深度,并取这些值的平均值作为槽的实际深度,槽的深度需达到0.0025毫米的精确度。对应于有限元模型中的这个深度的应变值通过测得的应变与深度关系插值得到。
[0059] 步骤五:计算纵向弹性模量EL和横向弹性模量ET
[0060] 骨材料被认为是各向异性的,本方法通过实验来求与深度相关的EL和ET的值而不是采用其他人得到的骨头的弹性模量值。将来自不同骨头的不同深度层的长度为9毫米左右、矩形截面宽为1毫米左右、厚度为0.4到0.5毫米的试件仔细研磨,然后在室温下储存在5%磷酸盐缓冲盐溶液48小时后中。将这些试件放在跨度为8毫米的微型三点弯装置中进行测试。用一个100倍显微镜观察每个试件的中点挠度,由式(5)可得
[0061]
[0062] 式中ω表示试件中点的位移,M表示试件中点的弯矩,l表示微型三点弯装置的跨度,E表示弹性模量,I表示试件截面的惯性矩。E的值可以通过式(11)求出来,在骨头的横向和纵向分别截取试件进行测试,可以分别求出横向和纵向弹性模量ET和EL。
[0063] 步骤六:求解残余应力
[0064] 将式(9)表示为{{ε}=[C]{A},根据步骤三测得的柔度系数Cij和步骤四测得的应变可以解得式(8)中的系数Aj
[0065] {A}=([C]T[C])-1CT{εmeas}   (12)
[0066] 将Aj代入式(8)即可求出骨表面不同深度的残余应力。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈