CN115791467A - 基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法 - Google Patents

基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法 Download PDF

Info

Publication number
CN115791467A
CN115791467A CN202111054494.8A CN202111054494A CN115791467A CN 115791467 A CN115791467 A CN 115791467A CN 202111054494 A CN202111054494 A CN 202111054494A CN 115791467 A CN115791467 A CN 115791467A
Authority
CN
China
Prior art keywords
unloading
press
calculating
stage
indentation
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.)
Pending
Application number
CN202111054494.8A
Other languages
English (en)
Inventor
吴海利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Electric Power Generation Equipment Co Ltd
Original Assignee
Shanghai Electric Power Generation Equipment Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Electric Power Generation Equipment Co Ltd filed Critical Shanghai Electric Power Generation Equipment Co Ltd
Priority to CN202111054494.8A priority Critical patent/CN115791467A/zh
Publication of CN115791467A publication Critical patent/CN115791467A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,包括如下步骤:使硬质球形压头对齐于金属试样的待测试压入区域;使用硬质球形压头在所述待测试压入区域实施多级压入卸载试验,记录压入载荷与压入深度之间的关系曲线;在压入载荷与压入深度之间的关系曲线上提取每一级压入卸载试验所对应的最大压入载荷Fmi、最大压入深度hmi、加载斜率Sli和初始卸载斜率Sui,并进一步计算得出包括卸载弹性模量E、塑性延伸强度Rpx、抗拉强度Rm在内的拉伸力学参数并利用幂函数拟合获得的反映应力与应变关系的本构曲线,解决了现有压入测试方法获取的强度预测值与真实值偏差较大和有限元模拟方程过于繁琐难以工程应用的技术问题。

Description

基于球形多级压入卸载法的材料力学参数与本构曲线的获取 方法
技术领域
本发明涉及材料测试领域,特别是一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法。
背景技术
拉伸强度和应力应变关系曲线是金属材料的基本力学表征,传统力学试验方法需要从原材料或结构上切取拉伸试样,在试验机上进行拉伸试验获取材料力学性能。然而对于在役设备、焊接件、微小结构件和表面处理如涂层等零构件无法加工试样进行传统拉伸试验,大量学者希望通过压入试验来获取材料的力学性能。压入试验无需取样,通过对表面微区施加压入载荷的简单操作获取压入载荷与压入深度之间的关系曲线,并由计算公式获得材料的力学性能。
目前压入法大致可分为经验物理法和有限元模拟数值分析法。经验物理法主要是定义代表应力和代表应变,分别近似为材料的应力和应变,从而实现从压入载荷与压入深度之间的关系曲线转化为应力与应变关系的本构曲线。大部分经验物理模型均基于纯弹性材料和理想线性变形的假设,材料表面的加工状态、压痕边缘鼓凸或沉陷现象以及非线性塑性变形没有被充分考虑和修正,导致结果与真值严重偏离,应用性较差。有的推算方法是以压入部分曲线为内凹形状为前提对压入曲线进行幂律拟合作为力学经验推算公式的基础,对于压入曲线呈直线或外凸形状的材料并不适用。
基于有限元模拟的数值分析方法考虑了材料加工硬化行为对表征应力-表征应变的影响,针对不同本构方程建立材料单轴力学性能与压入性能的关联方程。此类表达式完全依赖有限元计算,模拟方程中存在大量需要通过有限元结果标定的系数,方程式冗长繁琐,而且对于不同塑性特征的材料,需要重新标定系数。数值分析法并未从理论上描述材料的力学行为,缺乏明确的物理意义,参数的选择和迭代过程较为复杂,难以进行工程应用。
因此需要一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,用于无法加工试样进行传统拉伸试验获取材料力学参数的测试环境,解决现有压入测试方法获取的强度预测值与真实值偏差较大和有限元模拟方程过于繁琐难以工程应用的技术问题。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,解决现有压入测试方法获取的强度预测值与真实值偏差较大和有限元模拟方程过于繁琐难以工程应用的技术问题。
为实现上述目的及其他相关目的,本发明提供一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,包括如下步骤:
使硬质球形压头对齐于金属试样的待测试压入区域;
使用硬质球形压头在所述待测试压入区域实施多级压入卸载试验,记录压入载荷与压入深度之间的关系曲线;
在压入载荷与压入深度之间的关系曲线上提取每一级压入卸载试验所对应的最大压入载荷Fmi、最大压入深度hmi、加载斜率Sli和初始卸载斜率Sui,测量金属试样的表面压痕直径df
计算第一级卸载对应的修正残余深度hc1和等效弹性模量E1 *,获得卸载弹性模量E;
计算第i级卸载对应的修正残余深度hci,获得等效弹性模量Ei *
计算第i级卸载前最大压入载荷Fmi对应的应力σmi
通过第i级卸载对应的应力σmi和等效弹性模量Ei *计算第i级卸载的总弹性应变εtei
计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,表征总弹性应变εtei与总应变εti的比值,分解得到第i级卸载的塑性应变量εpi
将塑性应变量εpi和应力σmi进行幂函数拟合,获得压入强化系数KI和压入硬化指数nI;由压入强化系数KI和压入硬化指数nI获得金属试样所用材料的本构曲线;
计算金属试样所用材料的塑性延伸强度Rpx和抗拉强度Rm
优选的,所述获取方法还包括:对金属试样的待测试压入区域进行表面处理,待测试压入区域的表面的粗糙度小于0.8μm,通过这种方式进一步提高了通过试验获取的数据的准确性。
优选的,所述硬质球形压头的直径为D;将最后一级压入卸载试验所对应的最大压入载荷Fmi记为Fmf;满足以下条件:对于黑色金属,最大压入载荷Fmf满足Fmf=294D2;对于有色金属,最大压入载荷Fmf满足Fmf=98D2,保证了得到的压入载荷与压入深度之间的关系曲线的准确性。
优选的,所述使用硬质球形压头在所述待测试压入区域实施多级压入卸载试验的的步骤包括:始终保持所述硬质球形压头垂直于金属试样的待测试压入区域,确保产生的压入深度和压痕直径是载荷直接作用的结果,而不是其法向的分解载荷。
优选的,将最后一级压入卸载试验所对应的最大压入载荷Fmi记为Fmf;在0~Fmf的载荷范围内实施压入卸载试验的级数至少8级;保证实施多级压入卸载的级数大于等于8,能够保证通过试验获得的压入载荷与压入深度之间的关系曲线更加准确,且获得数据足够丰富。
优选的,所述修正残余深度hci的计算方法为:
根据关系式
Figure BDA0003254045900000031
计算第i级卸载对应的理论残余深度h* ci,其中α=0.75;
再根据关系式
Figure BDA0003254045900000032
计算塑性系数β,其中h* cf为关系式
Figure BDA0003254045900000033
取值i=f得到的最后一级卸载后的理论残余深度;hcf为从压入载荷与压入深度之间的关系曲线上获取的最后一级卸载的真实残余深度;Fmf为从压入载荷与压入深度之间的关系曲线上获取的最后一级卸载前的最大压入载荷;
将塑性系数β代入关系式
Figure BDA0003254045900000034
计算第i级卸载对应的修正残余深度hci;本方法根据从压入载荷与压入深度之间的关系曲线上获取的参数,并通过计算得到修正残余深度hci,为后续计算其他拉伸力学参数提供了基础。
优选的,所述等效弹性模量Ei *的计算方法为:首先根据关系式Aci=π(Dhci-hci 2)计算得出第i级卸载的接触投影面积Aci;再代入关系式
Figure BDA0003254045900000035
计算等效弹性模量Ei *;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到包括等效弹性模量Ei *,方便快捷。
优选的,所述计算第一级卸载对应的修正残余深度hc1和等效弹性模量E1 *,获得卸载弹性模量E的步骤包括:取i=1,得出第一级卸载的修正残余深度hc1,计算得到第一级卸载对应的等效弹性模量E1 *并代入关系式
Figure BDA0003254045900000036
计算得到金属试样的卸载弹性模量E;其中E压头为硬质球形压头的卸载弹性模量,v压头为硬质球形压头的泊松比;v为金属试样材料的泊松比,对于一般合金材料,取v=0.3;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到卸载弹性模量E,方便快捷。
优选的,所述计算第i级卸载前最大压入载荷Fmi对应的应力σmi的步骤包括:根据关系式Asmi=πDhmi计算第i级卸载前最大压入深度对应的压痕表面积Asmi,再代入关系式
Figure BDA0003254045900000041
计算第i级卸载前最大压入载荷Fmi对应的应力σmi;其中M为材料常数,对于金属材料,M的取值范围为2.8~3;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到最大压入载荷Fmi对应的应力σmi,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
优选的,所述表征总弹性应变εtei的计算步骤包括:根据关系式
Figure BDA0003254045900000042
计算第i级卸载对应的总弹性应变εtei;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到第i级卸载的总弹性应变εtei,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
优选的,所述计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,表征总弹性应变εtei与总应变εti的比值,分解得到第i级卸载的塑性应变量εpi的步骤包括:根据关系式
Figure BDA0003254045900000043
计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,即总弹性应变εtei与总应变εti的比值;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
优选的,所述分解得到第i级卸载的塑性应变量εpi的步骤包括:根据关系式
Figure BDA0003254045900000044
计算得到金属试样的第i级卸载的塑性应变量εpi;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到第i级卸载的塑性应变量εpi,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
优选的,所述将塑性应变量εpi和应力σmi进行幂函数拟合,获得压入强化系数KI和压入硬化指数nI的步骤包括:将每一级卸载对应的塑性应变量εpi和应力σmi按式
Figure BDA0003254045900000045
进行幂函数拟合,获取压入强化系数KI和压入硬化指数nI;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过对获得的试验数据进行幂函数拟合得到压入强化系数KI和压入硬化指数nI,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
优选的,所述由压入强化系数KI和压入硬化指数nI获得反映应力与应变关系的本构曲线的步骤包括:将压入强化系数KI、压入硬化指数nI以及卸载弹性模量E,代入关系式
Figure BDA0003254045900000051
获得用于反映应力与应变关系的本构曲线;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过幂函数拟合以及计算获得了反映应力与应变关系的本构曲线。
优选的,所述塑性延伸强度Rpx的计算步骤包括:将压入强化系数KI、压入硬化指数nI以及塑性延伸量的百分数x代入关系式
Figure BDA0003254045900000052
计算塑性延伸强度Rpx;本方法根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到塑性延伸强度Rpx,方便快捷。
优选的,所述抗拉强度Rm的计算步骤包括:
根据关系式Asmi=πDhmi计算获得最后一级卸载前最大压入深度hm对应的压痕表面积Asm
再根据关系式
Figure BDA0003254045900000053
计算获得压痕直径df对应的压痕表面积Asf
将Asm以及Asf代入关系式
Figure BDA0003254045900000054
得出比例系数γ;
根据关系式
Figure BDA0003254045900000055
计算得到最后一级卸载前的最大应力σmf
将比例系数γ以及最后一级卸载前的最大应力σmf代入关系式Rm=γσmf计算得到抗拉强度Rm;本方法根据测量获得的试验数据,进而通过计算得到抗拉强度Rm,方便快捷。
如上所述,本发明的一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,具有以下有益效果:本发明通过硬质球形压头和对待测的金属试样进行多级压入卸载试验,并根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到包括卸载弹性模量E、塑性延伸强度Rpx、抗拉强度Rm在内的拉伸力学参数并利用幂函数拟合获得的反映应力与应变关系的本构曲线;进一步解决了现有压入测试方法获取的强度预测值与真实值偏差较大和有限元模拟方程过于繁琐难以工程应用的技术问题。
附图说明
图1显示为本发明一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法的示意图;
图2显示为本发明一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法的压入载荷与压入深度之间的关系曲线示意图;
图3显示为本发明一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法获得的反映应力与应变关系的本构曲线与拉伸试验获得的反映应力与应变关系的本构曲线的对比示意图。
元件标号说明
1 硬质球形压头
2 金属试样
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效。
须知,本说明书所附图中所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容所能涵盖的范围内。同时,本说明书中所引用的如“上”、“下”、“左”、“右”、“中间”及“一”等的用语,亦仅为便于叙述的明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
如图1所示,本发明提供一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其步骤包括:
使硬质球形压头1对齐于金属试样2的待测试压入区域;
使用硬质球形压头1在所述待测试压入区域实施多级压入卸载试验,记录压入载荷与压入深度之间的关系曲线;
在压入载荷与压入深度之间的关系曲线上提取每一级压入卸载试验所对应的最大压入载荷Fmi、最大压入深度hmi、加载斜率Sli和初始卸载斜率Sui,测量金属试样2的表面压痕直径df
计算第一级卸载对应的修正残余深度hc1和等效弹性模量E1 *,获得卸载弹性模量E;
计算第i级卸载对应的修正残余深度hci,获得等效弹性模量Ei *
计算第i级卸载前最大压入载荷Fmi对应的应力σmi
通过第i级卸载对应的应力σmi和等效弹性模量Ei *计算第i级卸载的总弹性应变εtei
计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,表征总弹性应变εtei与总应变εti的比值,分解得到第i级卸载的塑性应变量εpi
将塑性应变量εpi和应力σmi进行幂函数拟合,获得压入强化系数KI和压入硬化指数nI;由压入强化系数KI和压入硬化指数nI获得金属试样2所用材料的本构曲线;
计算金属试样2所用材料的塑性延伸强度Rpx和抗拉强度Rm
在本实施例中,上述获取方法还包括:
对金属试样2的待测试压入区域进行表面处理,待测试压入区域的表面的粗糙度小于0.8μm,进一步提高了通过试验获取的数据的准确性。
在本实施例中,如图1所示,硬质球形压头1的直径为D;将最后一级压入卸载试验所对应的最大压入载荷Fmi记为Fmf;满足以下条件:对于黑色金属,最大压入载荷Fmf满足Fmf=294D2;对于有色金属,最大压入载荷Fmf满足Fmf=98D2,保证了得到的压入载荷与压入深度之间的关系曲线的准确性。
进一步的,在本实施例中,上述使用硬质球形压头1在待测试压入区域实施多级压入卸载试验的的步骤包括:始终保持硬质球形压头1垂直于金属试样2的待测试压入区域,确保产生的压入深度和压痕直径是载荷直接作用的结果,而不是其法向的分解载荷。
在本实施例中,将最后一级压入卸载试验所对应的最大压入载荷Fmi记为Fmf;在0~Fmf的载荷范围内实施压入卸载试验的级数至少8级;保证实施多级压入卸载的级数大于等于8,通过足够多的压入卸载试验的级数控制,能够保证通过试验获得的压入载荷与压入深度之间的关系曲线更加准确,且获得数据足够丰富。
进一步的,在本实施例中,如图1所示,使用直径为2.5mm的碳化钨合金材质的硬质球形压头1对10Cr材质的金属试样2表面进行单点法多级压入卸载试验,待测试压入区域表面应平整,粗糙度小于0.8μm;压入试验采用载荷控制,加载速率30N/s,卸载速率为50N/s,在0~1838N压向载荷范围内进行9级等间隔载荷压入卸载操作,记录压入卸载过程中的压入载荷与压入深度之间的关系曲线,如图2所示;使用光学显微镜测量10Cr材料表面的压痕直径df=0.90mm。
在本实施例中,多级压入卸载试验采用单点法,单点法是在一个压入点以预设压入速度和预设卸载速度进行多次多级压入卸载操作;单点法实施简单,速度更快,适用于材料试验表面空间有限的情况;进一步的,多级压入卸载试验还可以采用多点法;多点法是在不同压入点分别以预设压入速度和预设卸载速度施加不同压入载荷进行单次压入卸载操作,有效测试次数不少于8次,各压入点之间的间距大于10D;多点法适用于无法实施多级压入卸载操作的压入设备和材料试验表面空间大的情况。
进一步的,在本实施例中,所述修正残余深度hci的计算方法为:
根据关系式
Figure BDA0003254045900000081
计算第i级卸载对应的理论残余深度h* ci,其中α=0.75;
再根据关系式
Figure BDA0003254045900000082
计算塑性系数β,其中h* cf为关系式
Figure BDA0003254045900000083
取值i=f得到的最后一级卸载后的理论残余深度;hcf为从压入载荷与压入深度之间的关系曲线上获取的最后一级卸载的真实残余深度;Fmf为从压入载荷与压入深度之间的关系曲线上获取的最后一级卸载前的最大压入载荷;
将塑性系数β代入关系式
Figure BDA0003254045900000084
计算第i级卸载对应的修正残余深度hci;本实施例根据从如图2所示的压入载荷与压入深度之间的关系曲线上获取的参数,并通过计算得到修正残余深度hci,为后续计算其他拉伸力学参数提供了基础。
进一步的,在本实施例中,等效弹性模量Ei *的计算方法为:首先根据关系式Aci=π(Dhci-hci 2)计算得出第i级卸载的接触投影面积Aci;再代入关系式
Figure BDA0003254045900000085
计算等效弹性模量Ei *;本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到等效弹性模量Ei *,方便快捷。
在本实施例中,计算第一级卸载对应的修正残余深度hc1和等效弹性模量E1 *,获得卸载弹性模量E的步骤包括:取i=1,得出第一级卸载的修正残余深度hc1,计算得到第一级卸载对应的等效弹性模量E1 *并代入关系式
Figure BDA0003254045900000086
计算得到金属试样2的卸载弹性模量E;其中E压头为硬质球形压头1的弹性模量,v压头为硬质球形压头1的泊松比,在本实施例中,碳化钨合金球形压头1的E压头=710Gpa,v压头=0.21;v为金属试样2的材料泊松比,取v=0.3;进一步计算得到金属试样2的弹性模量E=219.91GPa;本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量并获得的试验数据,进而通过计算得到卸载弹性模量E,方便快捷。
在本实施例中,计算第i级卸载前最大压入载荷Fmi对应的应力σmi的步骤包括:根据关系式Asmi=πDhmi计算第i级卸载前最大压入深度对应的压痕表面积Asmi,再代入关系式
Figure BDA0003254045900000091
计算第i级卸载前最大压入载荷Fmi对应的应力σmi;其中M为材料常数,对于金属材料,M的取值范围为2.8~3,在本实施例中,取值M=2.9;本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到最大压入载荷Fmi对应的应力σmi,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
进一步的,在本实施例中,总弹性应变εtei的计算步骤包括:
根据关系式
Figure BDA0003254045900000092
计算第i级卸载对应的总弹性应变εtei;本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量并获得的试验数据,进而通过计算得到第i级卸载的总弹性应变εtei,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
在本实施例中,计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,表征总弹性应变εtei与总应变εti的比值,分解得到第i级卸载的塑性应变量εpi的步骤包括:
根据关系式
Figure BDA0003254045900000093
计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,即总弹性应变εtei与总应变εti的比值;
根据关系式
Figure BDA0003254045900000094
计算得到金属试样2的第i级卸载的塑性应变量εpi;本实施例根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到第i级卸载的塑性应变量εpi,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
在本实施例中,将上述将塑性应变量εpi和应力σmi进行幂函数拟合,获得压入强化系数KI和压入硬化指数nI的步骤包括:将每一级卸载对应的塑性应变量εpi和应力σmi按式
Figure BDA0003254045900000095
进行幂函数拟合,获取压入强化系数KI和压入硬化指数nI;在本实施例中获取压入强化系数KI=1023.7MPa和压入硬化指数nI=0.0477,本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过对获得的试验数据进行幂函数拟合得到压入强化系数KI和压入硬化指数nI,为后续获得反映应力与应变关系的本构曲线提供了数据支持。
在本实施例中,上述由压入强化系数KI和压入硬化指数nI获得反映应力与应变关系的本构曲线的步骤包括:将压入强化系数KI、压入硬化指数nI以及卸载弹性模量E,代入关系式
Figure BDA0003254045900000101
获得用于反映应力与应变关系的本构曲线,本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过幂函数拟合以及计算获得了反映应力与应变关系的本构曲线。
将本实施例获取的反映应力与应变关系的本构曲线与拉伸试验获取的反映应力与应变关系的本构曲线对比如图3所示,可知两条曲线基本重合。
在本实施例中,上述塑性延伸强度Rpx的计算步骤包括:将塑性延伸量的百分数x代入关系式
Figure BDA0003254045900000102
计算塑性延伸强度Rpx;在本实施例中采用x=0.2,代入数值计算得到Rp0.2=761MPa,本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到塑性延伸强度Rpx,方便快捷。
在本实施例中,上述抗拉强度Rm的计算步骤包括:
根据关系式Asmi=πDhmi计算获得最后一级卸载前最大压入深度hm对应的压痕表面积Asm
再根据关系式
Figure BDA0003254045900000103
计算获得压痕直径df对应的压痕表面积Asf
将Asm以及Asf代入关系式
Figure BDA0003254045900000104
得出比例系数γ;
根据关系式
Figure BDA0003254045900000105
计算得到最后一级卸载前的最大应力σmf
将比例系数γ以及最后一级卸载前的最大应力σmf代入关系式Rm=γσmf计算得到抗拉强度Rm;代入数值计算得到Rm=933MPa;本实施例根据测试得到的如图2所示的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到抗拉强度Rm,方便快捷。
将本发明实施例获取的力学参数和本构参数与材料拉伸试验结果对比如下表,可知本发明的实施例获得的材料卸载弹性模量E、屈服强度Rp0.2、抗拉强度Rm、强化系数K和强化指数n与拉伸试验结果相比偏差很小,证明了本发明技术的有效性和准确性。
Figure BDA0003254045900000111
综上所述,本发明的一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,通过硬质球形压头1和对待测的金属试样2进行多级压入卸载试验,并根据测试得到的压入载荷与压入深度之间的关系曲线以及测量获得的试验数据,进而通过计算得到材料力学参数并利用幂函数拟合获得的反映应力与应变关系的本构曲线;解决了现有压入测试方法获取的强度预测值与真实值偏差较大和有限元模拟方程过于繁琐难以工程应用的技术问题。所以,本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。

Claims (10)

1.本发明提供一种基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于,包括以下步骤:
使硬质球形压头(1)对齐于金属试样(2)的待测试压入区域;
使用硬质球形压头(1)在所述待测试压入区域实施多级压入卸载试验,记录压入载荷与压入深度之间的关系曲线;
在压入载荷与压入深度之间的关系曲线上提取每一级压入卸载试验所对应的最大压入载荷Fmi、最大压入深度hmi、加载斜率Sli和初始卸载斜率Sui,测量金属试样(2)的表面压痕直径df
计算第一级卸载对应的修正残余深度hc1和等效弹性模量E1 *,获得卸载弹性模量E;
计算第i级卸载对应的修正残余深度hci,获得等效弹性模量Ei *
计算第i级卸载前最大压入载荷Fmi对应的应力σmi
通过第i级卸载对应的应力σmi和等效弹性模量Ei *计算第i级卸载的总弹性应变εtei
计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,表征总弹性应变εtei与总应变εti的比值,分解得到第i级卸载的塑性应变量εpi
将塑性应变量εpi和应力σmi进行幂函数拟合,获得压入强化系数KI和压入硬化指数nI
由压入强化系数KI和压入硬化指数nI获得金属试样(2)所用材料的本构曲线;
计算金属试样(2)所用材料的塑性延伸强度Rpx和抗拉强度Rm
2.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:
所述硬质球形压头(1)的直径为D,将最后一级压入卸载试验所对应的最大压入载荷Fmi记为Fmf,满足以下条件:对于黑色金属,最大压入载荷Fmf满足Fmf=294D2;对于有色金属,最大压入载荷Fmf满足Fmf=98D2
待测试压入区域表面应平整,粗糙度小于0.8μm;始终保持所述硬质球形压头(1)垂直于金属试样(2)的待测试压入区域,在0~Fmf的载荷范围内实施压入卸载试验的级数至少8级。
3.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述计算等效弹性模量Ei *的步骤包括:
根据关系式
Figure FDA0003254045890000021
计算第i级卸载对应的理论残余深度h* ci,其中α=0.75;
再根据关系式
Figure FDA0003254045890000022
计算塑性系数β,其中h* cf为关系式
Figure FDA0003254045890000023
取值i=f得到的最后一级卸载后的理论残余深度;hcf为从压入载荷与压入深度之间的关系曲线上获取的最后一级卸载的真实残余深度;Fmf为从压入载荷与压入深度之间的关系曲线上获取的最后一级卸载前的最大压入载荷;
将塑性系数β代入关系式
Figure FDA0003254045890000024
计算第i级卸载对应的修正残余深度hci
根据关系式Aci=π(Dhci-hci 2)计算得出第i级卸载的接触投影面积Aci;再代入关系式
Figure FDA0003254045890000025
计算第i级卸载对应的等效弹性模量Ei *
4.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述计算第一级卸载对应的修正残余深度hc1和等效弹性模量E1 *,获得卸载弹性模量E的步骤包括:
取i=1,得出第一级卸载的修正残余深度hc1,计算得到第一级卸载对应的等效弹性模量E1 *并代入关系式
Figure FDA0003254045890000026
计算得到金属试样(2)的卸载弹性模量E;其中E压头为硬质球形压头(1)的卸载弹性模量,v压头为硬质球形压头(1)的泊松比;v为金属试样(2)材料的泊松比。
5.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述计算第i级卸载前最大压入载荷Fmi对应的应力σmi的步骤包括:
根据关系式Asmi=πDhmi计算第i级卸载前最大压入深度对应的压痕表面积Asmi,代入关系式
Figure FDA0003254045890000027
计算第i级卸载前最大压入载荷Fmi对应的应力σmi;其中M为材料常数,对于金属材料,M的取值范围为2.8~3。
6.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述表征总弹性应变εtei的计算步骤包括:
根据关系式
Figure FDA0003254045890000031
计算第i级卸载对应的总弹性应变εtei
7.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi和所述分解得到第i级卸载的塑性应变量εpi的步骤包括:
根据关系式
Figure FDA0003254045890000032
计算第i级卸载前后的加载斜率Sli和卸载斜率Sui的比值pi,即总弹性应变εtei与总应变εti的比值。
再根据关系式
Figure FDA0003254045890000033
计算得到金属试样(2)的第i级卸载的塑性应变量εpi
8.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述将塑性应变量εpi和应力σmi进行幂函数拟合,获得压入强化系数KI和压入硬化指数nI以及本构曲线的步骤包括:
将每一级卸载对应的塑性应变量εpi和应力σmi按式
Figure FDA0003254045890000034
进行幂函数拟合,获取压入强化系数KI和压入硬化指数nI
再将压入强化系数KI、压入硬化指数nI以及卸载弹性模量E,代入关系式
Figure FDA0003254045890000035
获得用于反映应力与应变关系的本构曲线。
9.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述塑性延伸强度Rpx的计算步骤包括:
将压入强化系数KI、压入硬化指数nI和塑性延伸量的百分数x代入关系式
Figure FDA0003254045890000036
计算塑性延伸强度Rpx
10.根据权利要求1所述的基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法,其特征在于:所述抗拉强度Rm的计算步骤包括:
根据关系式Asmi=πDhmi计算获得最后一级卸载前最大压入深度hm对应的压痕表面积Asm
再根据关系式
Figure FDA0003254045890000041
计算获得压痕直径df对应的压痕表面积Asf
将Asm和Asf代入关系式
Figure FDA0003254045890000042
得出比例系数γ;
根据关系式
Figure FDA0003254045890000043
计算得到最后一级卸载前的最大应力σmf
将比例系数γ以及最后一级卸载前的最大应力σmf代入关系式Rm=γσmf计算得到抗拉强度Rm
CN202111054494.8A 2021-09-09 2021-09-09 基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法 Pending CN115791467A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111054494.8A CN115791467A (zh) 2021-09-09 2021-09-09 基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111054494.8A CN115791467A (zh) 2021-09-09 2021-09-09 基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法

Publications (1)

Publication Number Publication Date
CN115791467A true CN115791467A (zh) 2023-03-14

Family

ID=85416917

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111054494.8A Pending CN115791467A (zh) 2021-09-09 2021-09-09 基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法

Country Status (1)

Country Link
CN (1) CN115791467A (zh)

Similar Documents

Publication Publication Date Title
CN108645704B (zh) 基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法
CN109299568B (zh) 基于纳米压痕试验的焊接接头本构模型反推方法
WO2014141794A1 (ja) 応力-ひずみ関係シミュレート方法、スプリングバック量予測方法およびスプリングバック解析装置
CN106644711B (zh) 一种延性材料单轴本构关系测试方法
CN103886125B (zh) 一种钛合金热复合成形数值模拟方法
CN111896373B (zh) 一种测定等效塑性应变成形极限图的试验与计算方法
CN109933925B (zh) 一种金属板材的冲压成形性能预测方法
CN108645706B (zh) 一种通过硬化强度和抗拉强度预测金属材料疲劳强度的方法
CN111339703A (zh) 大应变条件下材料硬化行为的虚拟预测方法
CN113420391B (zh) 一种获得复杂应力状态下材料高精度硬化模型参数的方法
CN110702513A (zh) 一种金属棒材大应变范围硬化曲线的试验测量方法
CN110749510B (zh) 基于有限元仿真检测金属材料弯曲性能的方法
CN110837675B (zh) 一种优化的断裂准则预测差厚板断裂的方法、装置和系统
CN113764056B (zh) 一种获得材料多应变率下高精度硬化模型参数的方法
Strzelecki et al. Analysis of axial load and bending load effects on the fatigue life
CN115791467A (zh) 基于球形多级压入卸载法的材料力学参数与本构曲线的获取方法
RU2324918C1 (ru) Способ оценки предельной деформации при локальной листовой штамповке
CN113358678B (zh) α钛变形过程介观应力和织构的半定量预测及可视化方法
Li et al. Identification of material parameters from punch stretch test
Zhang et al. A macro-pillar compression technique for determining true stress-strain curves of steels
RU2617798C1 (ru) Способ определения пластичности металлов и сплавов
CN114386216A (zh) 一种基于剪切应变路径构建高强钢加工硬化的方法
CN112924278B (zh) 一种用于高能重离子辐照样品的小冲杆测试装置和方法
CN115512792A (zh) 一种钛合金材料应变硬化模型的参数识别方法
CN115114824B (zh) 一种弹、塑性应变能比例叠加的延性材料压入试验方法

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