法律状态公告日
法律状态信息
法律状态
2020-07-10
未缴年费专利权终止 IPC(主分类):G01N3/42 授权公告日:20170201 终止日期:20190723 申请日:20140723
专利权的终止
2017-02-01
授权
授权
2015-01-14
实质审查的生效 IPC(主分类):G01N3/42 申请日:20140723
实质审查的生效
2014-12-24
公开
公开
技术领域
本发明属于材料力学性能测试领域。具体涉及一种利用仪器化压入仪和Berkovich压头测试金属材料应变硬化指数、弹性模量、条件屈服强度σ0.2及强度极限σb的方法。
背景技术
仪器化压入测试技术通过实时同步测量作用于金刚石压头上的压入载荷与金刚石压头压入被测材料表面的压入深度获得压入载荷-位移曲线,根据仪器化压入响应与被测材料力学性能参数间的无量纲函数关系式,可识别被测材料的诸多力学性能参数。
材料弹性模量的仪器化压入测试主要有W.C.Oliver和G.M.Pharr提出的“Oliver-Pharr方法”或“斜率方法”和马德军提出的“马德军方法”或“纯能量方法”。“斜率方法”的理论基础为小变形弹性理论,由于未考虑被测材料在压头作用下的塑性行为和几何大变形,使得“斜率方法”在应用于低应变硬化指数的被测材料时,测试结果严重偏离弹性模量真值。“纯能量方法”考虑了材料、几何和接触边界条件的非线性,其弹性模量的测试精度因此高于“斜率方法”。尽管如此,“纯能量方法”依然存在一定的理论测试误差,该误差源于被测材料的应变硬化指数未知,因此设法识别被测试材料的应变硬化指数是提高材料弹性模量仪器化压入测试精度的唯一有效途径。
材料应变硬化指数与屈服强度的仪器化压入测试目前存在基于球形压头的单一球压头压入法和基于不同锥顶角的多个锥压头压入法,其中应用单一球压头压入法遇到的困难是制造半径为几个或几十微米的球形压头其几何加工精度难以满足测试要求,因此,基于球形压头的材料应变硬化指数与屈服强度的仪器化压入测试方法在实际应用或工程化方面难有作为。应用多个锥压头压入法不存在压头制造方面的问题,但测试过程需要更换不同锥顶角的棱锥压头,同时需要对仪器柔度进行重新标定,而精确标定仪器柔度既耗时又困难,因此应用多锥压头压入法进行测试其效率较低。
针对目前金属材料弹塑性参数仪器化压入测试中存在的问题,本发明提出了一种基于Berkovich压痕的金属材料应变硬化指数、弹性模量、条件屈服强度σ0.2及强度极限σb的仪器化压入测试方法。
发明内容
本发明的目的是提供一种基于Berkovich压痕的金属材料弹塑性参数仪器化压入测试方法,利用该方法可以确定金属材料的弹塑性参数包括应变硬化指数、弹性模量、条件屈服强度σ0.2及强度极限σb。与使用两个或两个以上不同锥顶角的棱锥压头仪器化压入测试方法相比,该方法仅使用单一Berkovich压头对金属材料实施仪器化压入测试并辅以Berkovich压痕几何参数测试即可确定金属材料的应变硬化指数n、弹性模量E、条件屈服强度σ0.2及强度极限σb,避免了测试前需要单独设计加工不同于标准凌锥压头锥顶角的非标准棱锥压头问题,以及测试过程中需要更换压头及由此导致的需要对仪器柔度进行重新标定的问题,提高了测试效率。
为了实现上述目的,本发明采用如下的技术方案:
一种基于Berkovich压痕的金属材料弹塑性参数仪器化压入测试方法,该方法利用单一Berkovich压头仪器化压入金属材料所得载荷-位移曲线及压痕确定金属材料的应变硬化指数、弹性模量、条件屈服强度σ0.2及强度极限σb;首先,利用Berkovich压痕中边距与名义中边距之比和仪器化压入比功确定金属材料的应变硬化指数;其次,利用仪器化压入比功、仪器化压入名义硬度及测试所得应变硬化指数确定金属材料的弹性模量;最后,利用仪器化压入比功、仪器化压入名义硬度及测试所得弹性模量和应变硬化指数确定金属材料的条件屈服强度σ0.2与强度极限σb。具体包括以下步骤:
1)利用仪器化压入仪和金刚石Berkovich压头对被测材料实施某一最大压入载荷为Pm的仪器化压入测试,获得压入载荷-位移曲线,同时利用该曲线确定金刚石Berkovich压头的最大压入深度hm、名义硬度Hn=Pm/A(hm),其中,A(hm)为对应最大压入深度时的金刚石Berkovich压头横截面积,当不考虑金刚石Berkovich压头尖端钝化时
2)通过分别积分载荷-位移曲线关系中的加载曲线和卸载曲线计算压入加载功Wt、卸载功We,并在此基础上计算压入比功We/Wt;
3)借助显微镜分别量取Berkovich压痕中心至三个压痕边界的距离:d1、d2和d3,并确定中边距d=(d1+d2+d3)/3及其与名义中边距dn=hmtan65.274°之比d/dn;
4)根据4个不同硬化指数(n1=0,n2=0.15,n3=0.30,n4=0.45)下 的仪器化压入比功We/Wt与d/dn的关系
>
进一步根据非负原则确定被测试材料的应变硬化指数n:
n=max{n′,0}
表1.多项式系数aij(i=1,...,4;j=0,1,2)的取值
5)根据4个不同硬化指数ni(i=1,2,3,4)下的仪器化压入比功We/Wt与比值Hn/Ec的关系
>
进一步根据仪器化压入名义硬度Hn及比值Hn/Ec确定被测试材料与金刚石Berkovich压头的联合弹性模量Ec:
Ec=Hn/(Hn/Ec)
及被测试材料的弹性模量E:
>
其中,金刚石Berkovich压头的弹性模量Ei=1141GPa,泊松比νi=0.07,被测试材料的泊松比ν可根据材料手册确定;
表2.多项式系数bij(i=1,...,4;j=0,...,6)的取值
6)根据4个不同硬化指数ni(i=1,2,3,4)及3个不同被测试材料与金刚石压头平面应变弹性模量之比ηj(j=1,2,3)(η1=0.0671,η2=0.1917,η3=0.3834)下的仪器化压入比功We/Wt与屈服强度同名义硬度的比值关系
>
进一步根据仪器化压入名义硬度Hn及比值σy/Hn确定被测试材料的屈服强度σy:
σy=Hn(σy/Hn)
及由关系式σ0.2=σy1-n[σ0.2+0.002E]n确定被测试材料的条件屈服强度σ0.2;
表3.多项式系数cijk(i=1,...,4;j=1,2,3;k=0,...,6)的取值
7)计算εy=σy/E,并由关系式
>
其中,步骤5)中,如果被测材料的泊松比不能由材料手册确定,则取值为0.3。
与使用两个或两个以上不同锥顶角的棱锥压头仪器化压入测试方法相比,本发明仅使用单一Berkovich压头对金属材料实施仪器化压入测试并辅以Berkovich压痕几何参数测试即可确定金属材料的应变硬化指数n、弹性模量E、条件屈服强度σ0.2及强度极限σb,避免了测试前需要单独设计加工不同于标准凌锥压头锥顶角的非标准棱锥压头问题,以及测试过程中需要更换压头及由此导致的需要对仪器柔度进行重新标定的问题,提高了测试效率。
附图说明:
图1a是鼓凸情况下的Berkovich压痕示意图;
图1b是沉陷情况下的Berkovich压痕示意图;
图2是Berkovich压头示意图;
图3是仪器化压入加、卸载曲线及加、卸载功示意图;
图4a是对应n=0,η分别取0.0671、0.1917和0.3834三个数值的Hn/Er-We/Wt关系图;
图4b是对应n=0.15,η分别取0.0671、0.1917和0.3834三个数值的Hn/Er-We/Wt关系图;
图4c是对应n=0.30,η分别取0.0671、0.1917和0.3834三个数值的Hn/Er-We/Wt关系图;
图4d是对应n=0.45,η分别取0.0671、0.1917和0.3834三个数值的Hn/Er-We/Wt关系图;
图5a是对应n=0,η分别取0.0671、0.1917和0.3834三个数值的Hn/Ec-We/Wt关系图;
图5b是对应n=0.15,η分别取0.0671、0.1917和0.3834三个数值的Hn/Ec-We/Wt关系图;
图5c是对应n=0.30,η分别取0.0671、0.1917和0.3834三个数值的Hn/Ec-We/Wt关系图;
图5d是对应n=0.45,η分别取0.0671、0.1917和0.3834三个数值的Hn/Ec-We/Wt关系图;
图6是式(16)所代表的n分别取0、0.15、0.30和0.45时的Hn/Ec-We/Wt关系图;
图7是对应不同n和η的d/dn-We/Wt关系图;
图8a是对应η=0.0671,n分别取0、0.15、0.30和0.45四个数值的σy/Hn-We/Wt关系图;
图8b是对应η=0.1917,n分别取0、0.15、0.30和0.45四个数值的σy/Hn-We/Wt关系图;
图8c是对应η=0.3834,n分别取0、0.15、0.30和0.45四个数值的σy/Hn-We/Wt关系图;
图9是6061铝合金的仪器化压入载荷-位移曲线;
图10是S45C碳钢的仪器化压入载荷-位移曲线;
图11是SS316不锈钢的仪器化压入载荷-位移曲线;
图12是黄铜的仪器化压入载荷-位移曲线;
图13是分别采用仪器化压入测试和标准单轴拉伸测试所得6061铝合金的真实应力-应变关系的比较;
图14是分别采用仪器化压入测试和标准单轴拉伸测试所得S45C碳钢的真实应力-应变关系的比较;
图15是分别采用仪器化压入测试和标准单轴拉伸测试所得SS316不锈钢的真实应力-应变关系的比较;
图16是分别采用仪器化压入测试和标准单轴拉伸测试所得黄铜的真实应力-应变关系的比较。
具体实施方式
以下通过结合附图对本发明的方法进行详细说明,但这些实施例仅仅是例示的目的,并不旨在对本发明的范围进行任何限定。本申请提出了一种基于Berkovich压痕的金属材料弹塑性参数仪器化压入测试方法,该方法利用单一Berkovich压头仪器化压入金属材料所得载荷-位移曲线及压痕确定金属材料的应变硬化指数、弹性模量、条件屈服强度σ0.2及强度极限σb;首先,利用Berkovich压痕中边距与名义中边距之比和仪器化压入比功确定金属材料的应变硬化指数;其次,利用仪器化压入比功、仪器化压入名义硬度及测试所得应变硬化指数确定金属材料的弹性模量;最后,利用仪器化压入比功、仪器化压入名义硬度及测试所得弹性模量和应变硬化指数确定金属材料的条件屈服强度σ0.2与强度极限σb。具体包括以下步骤:
1)利用仪器化压入仪和金刚石Berkovich压头对被测材料实施某一最大压入载荷为Pm的仪器化压入测试,获得压入载荷-位移曲线,同时利用该曲线确定金刚石Berkovich压头的最大压入深度hm、名义硬度Hn=Pm/A(hm),其中,A(hm)为对应最大压入深度时的金刚石Berkovich压 头横截面积,当不考虑金刚石Berkovich压头尖端钝化时
2)通过分别积分载荷-位移曲线关系中的加载曲线和卸载曲线计算压入加载功Wt、卸载功We,并在此基础上计算压入比功We/Wt;
3)借助显微镜分别量取Berkovich压痕中心至三个压痕边界的距离:d1、d2和d3,并确定中边距d=(d1+d2+d3)/3及其与名义中边距dn=hmtan65.274°之比d/dn;
4)根据4个不同硬化指数(n1=0,n2=0.15,n3=0.30,n4=0.45)下的仪器化压入比功We/Wt与d/dn的关系
>
进一步根据非负原则确定被测试材料的应变硬化指数n:
n=max{n′,0}
表1.多项式系数aij(i=1,...,4;j=0,1,2)的取值
5)根据4个不同硬化指数ni(i=1,2,3,4)下的仪器化压入比功We/Wt与比值Hn/Ec的关系
>
进一步根据仪器化压入名义硬度Hn及比值Hn/Ec确定被测试材料与金刚石Berkovich压头的联合弹性模量Ec:
Ec=Hn/(Hn/Ec)
及被测试材料的弹性模量E:
>
其中,金刚石Berkovich压头的弹性模量Ei=1141GPa,泊松比νi=0.07,被测试材料的泊松比ν可根据材料手册确定;
表2.多项式系数bij(i=1,...,4;j=0,...,6)的取值
6)根据4个不同硬化指数ni(i=1,2,3,4)及3个不同被测试材料与金刚石压头平面应变弹性模量之比ηj(j=1,2,3)(η1=0.0671,η2=0.1917,η3=0.3834)下的仪器化压入比功We/Wt与屈服强度同名义硬度的比值关系
>
进一步根据仪器化压入名义硬度Hn及比值σy/Hn确定被测试材料的屈服强度σy:
σy=Hn(σy/Hn)
及由关系式σ0.2=σy1-n[σ0.2+0.002E]n确定被测试材料的条件屈服强度σ0.2;
表3.多项式系数cijk(i=1,...,4;j=1,2,3;k=0,...,6)的取值
7)计算εy=σy/E,并由关系式
>
以下详细说明本发明的形成过程。Berkovich压痕示意图如附图1a和附图1b所示,定义Berkovich压痕中边距d为Berkovich压痕中心至三个压痕边界距离d1、d2及d3的平均值,即d=(d1+d2+d3)/3。金刚石Berkovich压头示意图如附图2所示,根据最大压入深度hm定义Berkovich名义压痕中边距dn=hmtan65.274°。仪器化压入载荷-位移曲线示意图如附图3所示,纵轴表示压入载荷P,横轴表示压入深度h,加载曲线为1,卸载曲线为2,加载功Wt区域为3,卸载功We区域为4。仪器化压入所设定的最大压入载荷为Pm,与之相对应的最大压入深度为hm。用A(hm)表示金刚石Berkovich压头在最大压入深度hm位置处的金刚石Berkovich压头横截面积,则名义硬度Hn被定义为最大压入载荷Pm与金刚石Berkovich压头横截面积A(hm)之比,即Hn=Pm/A(hm)。进一步定义仪器化压入加载功Wt和卸载功We分别为实施仪器化压入时金刚石Berkovich压头在加载阶段和卸载阶段所做的功,其值分别等于加载曲线和卸载曲线与仪器化压入载荷-位移曲线横坐标所围面积。仪器化压入比功We/Wt为卸载功We与加载功Wt的比值。
将金刚石Berkovich压头视为弹性体,其弹性模量与泊松比分别用Ei和vi表示;被测材料视为弹塑性体,其单轴真实应力-应变关系由线弹性和Hollomon幂硬化函数组成,同时其弹性模量与泊松比分别用E和ν表示,屈服强度与应变硬化指数分别用σy和n表示。基于上述设定及忽略金刚石Berkovich压头与被测试材料间的摩擦,则仪器化压入名义硬度Hn、仪器化压入比功We/Wt及Berkovich压痕中边距与名义中边距的比值d/dn可以分别表示为被测材料的屈服强度σy、应变硬化指数n、弹性模量E、泊松比ν与金刚石Berkovich压头的弹性模量Ei、泊松比vi以及最大压入 深度hm的函数:
>
>
>
其中E/(1-ν2)和
E/(1-ν2)=(η+1)Er (4)
>
于是,式(1)、(2)和(3)可以被改写为:
Hn=ΓH2(σy,n,Er,η,hm) (6)
We/Wt=ΓW2(σy,n,Er,η,hm) (7)
d/dn=ΓD2(σy,n,Er,η,hm) (8)
应用量纲∏定理,式(6)、(7)和(8)可简化为:
Hn/Er=ΓH3(σy/Er,n,η) (9)
We/Wt=ΓW3(σy/Er,n,η) (10)
d/dn=ΓD3(σy/Er,n,η) (11)
由式(10)可得:
>
将式(12)代入式(9)和式(11)得:
Hn/Er=ΓH4(We/Wt,n,η) (13)
d/dn=ΓD4(We/Wt,n,η) (14)
由式(12)和式(13)可得:
σy/Hn=Γ5(We/Wt,n,η) (15)
通过有限元数值模拟可获得式(13)、式(14)和式(15)的显式解。模拟中金刚石Berkovich压头的弹性模量取值为Ei=1141GPa,泊松比取值为νi=0.07。被测材料弹性模量E的取值分别设为70GPa、200GPa和400GPa;屈服强度σy的取值范围为0.7~160000MPa;应变硬化指数n的取值为0、0.15、0.3和0.45;泊松比ν取固定值0.3。被测材料与金刚石Berkovich压头的平面应变弹性模量之比η分别为0.0671、0.1917和0.3834;被测材料与金刚石Berkovich压头间的接触摩擦系数取值为零。
附图4a、附图4b、附图4c和附图4d为对应不同n和η的Hn/Er-We/Wt关系图,从图中可以看出,对于确定的应变硬化指数n,η对Hn/Er-We/Wt关系有一定的影响,这表明折合弹性模量Er不能准确反映被测材料和金刚石Berkovich压头之间的综合弹性效应。为此,定义联合弹性模量 >并将其替代折合弹性模量Er可得Hn/Ec-We/Wt关系,结果如附图5a、附图5b、附图5c和附图5d所示,从图中可以看出,对于确定的应变硬化指数n,Hn/Ec-We/Wt关系几乎不受η的影响。于是,可以利用多项式函数对应变硬化指数n的4个不同取值情况下的Hn/Ec-WE/Wt关系进行曲线拟合,结果表示为:
>
其中,i=1,...,4分别对应应变硬化指数n的4个不同取值:0,0.15,0.3,0.45;系数bij(j=0,...,6)的取值见表2。式(16)所代表的n分别取0、0.15、0.30和0.45时的Hn/Ec-We/Wt关系如附图6所示。
表2.系数bij(i=1,...,4;j=0,...,6)的取值
附图7为对应不同n和η的d/dn-We/Wt关系图,从图中可以看出,对于确定的应变硬化指数n,η对d/dn-We/Wt关系的影响可以忽略。因此,可以利用多项式函数对应变硬化指数n的4个不同取值情况下的d/dn-We/Wt关系进行曲线拟合,结果表示为:
>
其中,i=1,...,4分别对应应变硬化指数n的4个不同取值:0,0.15,0.3,0.45;系数aii(j=0,1,2)的取值见表1。
表1.系数aij(i=1,...,4;j=0,1,2)的取值
附图8a、附图8b和附图8c为对应不同n和η的σy/Hn-We/Wt关系图。利用多项式函数对σy/Hn-We/Wt关系进行拟合,结果可表示为:
>
其中,i=1,...,4对应n的取值为0,0.15,0.3,0.45;j=1,2,3对应η的取值为0.0671,0.1917,0.3834;系数cijk(k=0,...,6)的取值见表3。
表3.系数cijk(i=1,...,4;j=1,2,3;k=0,...,6)的取值
应用实施例
选择6061铝合金、S45C碳钢、SS316不锈钢和黄铜进行仪器化压入实验。根据发明人所提实验步骤,应用自行研制且已获得国家发明专利授权的高精度仪器化压入仪[马德军,宋仲康,郭俊宏,陈伟.一种高精度压入仪及金刚石压头压入试样深度的计算方法.专利号:ZL201110118464.9]和金刚石Berkovich压头对6061铝合金、S45C碳钢、SS316不锈钢和黄铜不同区域重复进行5次仪器化压入实验。图9、图10、图11和图12分别为6061铝合金、S45C碳钢、SS316不锈钢和黄铜的仪器化压入载荷-位移曲线。应用光学显微镜可分别观测6061铝合金、S45C碳钢、SS316不锈钢以及黄铜的Berkovich压痕中边距。
根据仪器化压入载荷-位移曲线及Berkovich压痕,可以分别确定被测材料的仪器化压入名义硬度Hn、仪器化压入比功We/Wt及Berkovich压痕中边距与名义中边距的比值d/dn,结果如表4,在此基础上应用发明人 所提方法便可确定被测试材料的应变硬化指数n、弹性模量E、条件屈服强度σ0.2及强度极限σb。为了与标准单轴拉伸试验结果进行比较,将仪器化压入实验所用6061铝合金、S45C碳钢、SS316不锈钢和黄铜的相同材料分别制成标准单轴拉伸试样,并对其分别实施2次标准单轴拉伸试验,以2次试验的平均值作为材料单轴拉伸试验的测试结果,则由标准单轴拉伸试验测定的6061铝合金的弹性模量、应变硬化指数、条件屈服强度及强度极限分别为E单轴=71GPa、n单轴=0.052、σ0.2单轴=299.37MPa及σb单轴=366.25MPa;由标准单轴拉伸试验测定的S45C碳钢的弹性模量、应变硬化指数、条件屈服强度及强度极限分别为E单轴=201GPa、n单轴=0.176、σ0.2单轴=431.08MPa及σb单轴=612.84MPa;由标准单轴拉伸试验测定的SS316不锈钢的弹性模量、应变硬化指数、条件屈服强度及强度极限分别为E单轴=184GPa、n单轴=0.134、σ0.2单轴=610.11MPa及σb单轴=827.51MPa;由标准单轴拉伸试验测定的黄铜的弹性模量、应变硬化指数、条件屈服强度及强度极限分别为E单轴=83GPa、n单轴=0.125、σ0.2单轴=346.67MPa及σb单轴=421.23MPa。将6061铝合金、S45C碳钢、SS316不锈钢和黄铜的弹性模量、应变硬化指数、条件屈服强度和强度极限的仪器化压入测试结果与单轴拉伸试验结果进行比较,可以确定仪器化压入测试结果的测试误差:EErr=(E-E单轴)/E单轴、Δn=n-n单轴、σ0.2Err=(σ0.2-σ0.2单轴)/σ0.2单轴及σbErr=(σb-σb单轴)/σb单轴,结果见表4。从表中可以看出,6061铝合金、S45C碳钢、SS316不锈钢和黄铜的弹性模量相对测试误差分别为-4.52%、-1.66%、-2.03%和12.79%,应变硬化指数的绝对测试误差分别为0.015、0.011、0.020和-0.013,条件屈服强度σ0.2的相对测试误差分别为0.98%、-11.74%、3.79%和-1.42%,强度极限σb的相对测试误差分别为-9.64%、7.70%、9.88%和1.45%。进一步根据仪器化压入实验测得的6061铝合金、S45C碳钢、SS316不锈钢和黄铜的应变硬化指数n、弹性模量E和条件屈服强度σ0.2的平均值可以绘制其真实应力-应变关系,该关系与标准单轴拉伸试验测得的真实应力-应变关系的比较如附图13、图14、图15和图16所示,在附图13、图14、图15和图16中,横轴为真实应变ε,纵轴为真实应力σ,虚线为仪器化压入测试结果,粗实线为单轴拉伸试验一,细实线为单轴拉伸试验二。从图中可以看出两者具有较好的一致性。纵观以上实验结果表明,发明人所提基于Berkovich压痕的金属材料弹塑性参数仪器化压入测试方法是可行和非常有效的。
表4.6061铝合金、S45C碳钢、SS316不锈钢和黄铜弹塑性参数仪器化压入测试结果与测试误差
尽管上文对本发明的具体实施方式给予了详细描述和说明,但是应该指明的是,我们可以依据本发明的构想对上述实施方式进行各种等效改变和修改,其所产生的功能作用仍未超出说明书及附图所涵盖的精神时,均应在本发明的保护范围之内。
机译: 自动化音乐创作和生成系统,自动化音乐创作和生成过程,自动化音乐创作和生成,玩具乐器,音乐伴奏和音乐创作玩具乐器,自动化创作玩具乐器系统和音乐生成,电子信息处理和显示系统,企业基于互联网的一流音乐创作和生成系统,用于自动生成和传送数字复合音乐的网络系统,用于音乐环境的基于独立音乐的音乐创作和表演系统人工智能,基于音乐的自主创作过程音乐的生成和表演人工智能,自主分析仪器系统,用于建立自动音乐创作和生成引擎的网络,几何方法音乐理论系统操作参数映射,以自动方式构成和生成数字音乐的方法,参数转换
机译: 压痕装置,仪器化的测量系统以及通过该压痕方法确定材料的机械性能的方法
机译: 压头压痕参数的仪器化测量系统