CN111189699B - 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 - Google Patents
一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 Download PDFInfo
- Publication number
- CN111189699B CN111189699B CN201911271059.3A CN201911271059A CN111189699B CN 111189699 B CN111189699 B CN 111189699B CN 201911271059 A CN201911271059 A CN 201911271059A CN 111189699 B CN111189699 B CN 111189699B
- Authority
- CN
- China
- Prior art keywords
- parameters
- model
- calculating
- load
- nanoindentation
- 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
- 239000000463 material Substances 0.000 title claims abstract description 127
- 239000013078 crystal Substances 0.000 title claims abstract description 57
- 229920003023 plastic Polymers 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000002474 experimental method Methods 0.000 title claims abstract description 16
- 238000009864 tensile test Methods 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 230000002068 genetic effect Effects 0.000 claims abstract description 13
- 239000003795 chemical substances by application Substances 0.000 claims abstract description 6
- 238000007373 indentation Methods 0.000 claims description 55
- 238000006073 displacement reaction Methods 0.000 claims description 54
- 238000012360 testing method Methods 0.000 claims description 24
- 238000005070 sampling Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 14
- 238000005457 optimization Methods 0.000 claims description 13
- 239000007769 metal material Substances 0.000 claims description 10
- 238000005482 strain hardening Methods 0.000 claims description 8
- 238000009825 accumulation Methods 0.000 claims description 7
- 238000005498 polishing Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000012512 characterization method Methods 0.000 abstract description 2
- 239000010963 304 stainless steel Substances 0.000 description 3
- 229910000589 SAE 304 stainless steel Inorganic materials 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000002178 crystalline material Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229920001887 crystalline plastic Polymers 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/40—Investigating hardness or rebound hardness
- G01N3/42—Investigating hardness or rebound hardness by performing impressions under a steady load by indentors, e.g. sphere, pyramid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/32—Polishing; Etching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0218—Calculations based on experimental data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/026—Specifications of the specimen
- G01N2203/0286—Miniature specimen; Testing on microregions of a specimen
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于材料力学性能表征技术领域,提供一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法。该方法,首先,运用Oliver‑Pharr方法获得材料弹性模量;其次,使用分段线性/幂律硬化材料模型,联合MATLAB和ABAQUS建立纳米压痕宏观参数反演模型,利用堆积/沉陷参数对实际纳米压痕实验数据进行修正,结合Kriging代理模型和遗传算法计算出压痕材料宏观本构参数;最后,再联合MATLAB和ABAQUS建立拉伸试件基于晶体塑性有限元的多晶有限元模型,根据所求材料本构参数结合Kriging代理模型和遗传算法计算出晶体塑性材料参数。本发明相比于现有技术能够提高计算数据的精确度,显著减少计算量,并提高计算收敛性,在晶体塑性材料参数反演识别中都具有很高的实用价值和参考意义。
Description
技术领域
本发明属于材料力学性能表征技术领域,涉及一种晶体塑性材料参数的基于纳米压痕实验的反演标定金属材料微观本构参数的方法。
背景技术
材料的细观力学行为直接影响材料的强度及其它宏观力学性能,从细观尺度研究材料的力学行为,有助于深化对材料变形与损伤的机理的认识,对材料的使用及性能的改进有重要意义。在细观力学研究中,晶体材料的强化是晶体材料的弹塑性本构描述的重要部分。Peirce 等人在期刊《Acta Metallurgica》1982年第6期中提出了晶体滑移硬化模量的简单形式:
hαβ=qh(γ)(α≠β)
其中hαβ为滑移硬化模量,分为自硬化模量hαα和潜硬化模量hαβ(α≠β),γ为所有滑移系上泰勒累积剪切应变,h0为初始硬化模量,τs为阶段Ⅰ的应力,τ0为初始屈服应力,q为常数。这些晶体塑性本构参数的准确识别是研究材料晶体塑性力学行为的基础。
随着高分辨率测试设备的出现,压痕测试已成为表征各种材料,尤其是体积或尺寸微小的材料的力学性能最常用的技术之一。运用纳米压痕技术,可以对晶体塑性本构参数的进行准确的识别。这种测试方法要求准确了解压痕试样上接触力和接触深度之间的关系。在压痕试验过程中,被压材料的流动可能因其力学性能的不同而不同。压痕接触区域周围的材料可以沿施加荷载的方向向上(堆积)或向下(沉陷)变形。这种表面变形模式会影响压头与试样的真实接触面积,进而会影响测量精度。现有基于纳米压痕的晶体塑性本构参数的识别方法中,需要建立纳米压痕的晶体塑性有限元模型,其计算量大,收敛性差,而且在计算中没有考虑标定参数下压痕模型堆积/沉陷变形与实际测试中堆积/沉陷变形的关系。由于反演问题解的唯一性很难保证,需要反演的参数越少,越容易得到准确解。而运用纳米压痕求解弹性参数(弹性模量)的技术已经非常成熟,因此先求解出弹性参数,然后再反演出材料的塑性参数不仅可以减少计算量,而且更容易得到准确解。然而在没有准确的载荷-位移(压入深度) 曲线时,计算结果可能存在较大的误差,因而难以对材料本构参数进行快速准确的识别。
发明内容
为克服上述现有技术的缺陷,本发明的目的在于提供一种晶体塑性材料参数的基于纳米压痕实验的反演标定金属材料微观本构参数的方法。通过金属材料的纳米压痕测试获得载荷- 位移(压入深度)曲线和接触刚度等压痕响应。通过ABAQUS软件建立纳米压痕的常规有限元模型,运用分段线性/幂律硬化材料模型模拟金属材料上的纳米压痕过程;联合MATLAB 和ABAQUS建立参数反演模型,运用拉丁超立方抽样抽取分段线性/幂律硬化材料模型本构参数作为输入变量,以常规有限元压痕的载荷-位移曲线和压痕堆积/沉陷参数为输出变量,并计算出模拟数据与实验数据的残差,利用MATLAB建立本构参数与残差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出纳米压痕实验材料的分段线性/幂律硬化材料模型本构参数。运用压痕堆积/沉陷参数修正实验中的载荷-位移曲线,重复上述过程,直至前后两次的计算误差在允许范围。然后运用晶体塑性有限元建立拉伸试件的多晶体有限元模型。将晶体塑性本构参数作为设计变量,使用拉丁超立方抽样抽取不同的晶体塑性材料参数计算得到多晶体材料的应力应变曲线,并将其与分段线性/幂律强化模型的应力应变曲线对比,计算两组数据的均方差。以两组数据的均方差最小为目标,采用基于遗传算法的优化方法对Kriging代理模型的相关参数进行单目标寻优,计算出晶体塑性材料参数。
为了达到上述目的,本发明采用的技术方案如下:
一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法,包括以下方法:
步骤1:待测金属材料的纳米压痕实验
1-1:将待测金属材料切割,运用机械抛光和振动抛光,得到符合要求的纳米压痕试件;
1-2:运用纳米压痕***对1-1中的压痕试件进行压痕测试,获得包括载荷-位移曲线、最大载荷、接触刚度和接触硬度的实验压痕响应。运用Oliver-Pharr方法获得材料的弹性模量E。
步骤2:联合MATLAB和ABAQUS建立纳米压痕的基于分段线性/幂律硬化材料模型的常规有限元模型,反求出材料的宏观本构参数(屈服应力σy和应***化指数n)。其中分段线性/幂律硬化材料模型描述为:
其中ε为总应变,σ为应力。
2-1:运用ABAQUS建立纳米压痕的二维轴对称有限元模型。采用位移控制加载的方式计算得到压头沿压入方向的接触反力和位移,输出接触力、试件接触面的接触压力和位移和压头最下方节点的位移,生成输入文件;
2-2:在MATLAB中运用拉丁超立方抽样抽取分段线性/幂律硬化材料模型本构参数,修改2-1中输入文件中的材料参数,计算得到各组抽样参数下的压痕的载荷-位移曲线和压痕堆积/沉陷参数s/h(s为堆积或沉陷高度,出现堆积现象时为s正,出现沉陷现象时为s负;h 为压入深度),并计算出模拟载荷-位移曲线与实验载荷-位移曲线的均方差,利用MATLAB 建立本构参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出实验材料的分段线性/幂律硬化材料模型本构参数,与运用 Oliver-Pharr方法计算材料的弹性模量一起记为C0;
2-3:用2-2中的堆积/沉陷参数修正实验载荷-位移曲线,得到修正的载荷-位移曲线,然后再利用2-2中的模拟载荷-位移曲线计算其与修正的载荷-位移曲线的均方差,重复2-2中的单目标寻优过程(利用MATLAB建立本构参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优),计算出修正后材料的分段线性/幂律硬化材料模型本构参数,与运用Oliver-Pharr方法计算材料的弹性模量一起,记为C1;
2-4:计算2-2中计算出的材料本构参数与2-3中修正的本构参数之间的误差,如果误差在允许范围之内,则以2-3中修正的本构参数为材料的宏观本构参数;如果误差超过允许范围,则将2-3中修正的载荷-位移曲线作为实验载荷-位移曲线,重复2-2,2-3和2-4中的步骤,直至误差在允许范围之内。
步骤3:联合MATLAB和ABAQUS,运用晶体塑性有限元建立拉伸试件的多晶体有限元模型,求解晶体塑性材料参数与分段线性/幂律硬化材料模型本构参数的对应关系,进而反求待测材料的晶体塑性材料参数。
3-1:在ABAQUS中建立标准拉伸试件的晶体塑性有限元模型,运用ABAQUS材料子程序赋予模型晶体塑性材料参数h0、τs和τ0。采用载荷控制加载的方式计算得到材料的应力应变曲线。生成输入文件;
3-2:在MATLAB中运用拉丁超立方抽样抽取晶体塑性材料参数,修改3-1中输入文件中的材料参数,计算得到各组抽样参数下的应力应变曲线,并计算出模拟的应力应变曲线与2-4中宏观材料参数下的应力应变曲线的均方差,利用MATLAB建立晶体塑性材料参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出实验材料的晶体塑性材料参数。
进一步的,所述的步骤2-4中,误差的允许范围为0~2%。
进一步的,所述的步骤2中,材料参数反求的过程中将弹性参数和塑性参数分开,弹性参数运用成熟的理论方法求解,塑性参数运用有限元反演求解。塑性参数运用有限元反演求解过程中,运用拉丁超立方抽样抽取分段线性/幂律硬化材料模型本构参数,修改输入文件中的材料参数计算得到各组抽样参数下的压痕的载荷-位移曲线和压痕堆积/沉陷参数s/h,并计算出模拟载荷-位移曲线与实验载荷-位移曲线的均方差。在所述的步骤3晶体塑性参数反求过程中将压痕的晶体塑性有限元模型转化为压痕的常规有限元模型和拉伸试件的晶体塑性有限元模型。
与现有技术相比,本发明的技术效果在于:
(1)本发明在参数反求过程中将弹性参数和塑性参数分开,在考虑了压痕堆积/沉陷现象对压入深度的影响,并用有限元的方法修正了载荷-位移曲线之后,弹性参数运用成熟的理论方法求解,塑性参数运用有限元反演求解,使材料参数的反求过程中的原始数据的精确度得到了提高,计算量减小;
(2)本发明将晶体塑性参数反求过程中的压痕的晶体塑性有限元模型转化为压痕的常规有限元模型和拉伸试件的晶体塑性有限元模型。由于压痕在模拟的过程中涉及到接触非线性问题,因而压痕的晶体塑性有限元模型计算量大大增加而且收敛性差,压痕的常规有限元模型和拉伸试件的晶体塑性有限元模型计算量小而收敛性好。因此本发明的方法计算量小,计算速度高,计算收敛性高,在晶体塑性材料参数反演识别中都具有很高的实用价值和参考意义。
附图说明
图1为本发明流程示意图;
图2为测试载荷-位移曲线;
图3为纳米压痕的二维轴对称有限元模型;
图4为拉伸试件的有限元模型;
图5为纳米压痕反演的应力应变曲线与拉伸试验的应力应变曲线。
具体实施方式
以下结合具体实施例对本发明作进一步说明。
参见附图1所示,晶体塑性材料参数的基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法的具体实施步骤包括:
步骤1:待测金属材料的纳米压痕实验
1-1:选取304不锈钢材料为试件,将材料切割,运用机械抛光和振动抛光,得到符合要求的纳米压痕试件;
1-2:运用纳米压痕xp***对压痕试件进行压痕测试。测试中,压入深度设置为2微米,获得包括载荷-位移曲线、最大载荷、接触刚度和接触硬度的实验压痕响应。多次重复测试,获得5个以上的有效测试点。测试载荷-位移曲线如图2。同时运用Oliver-Pharr方法计算材料的弹性模量E为196.08GPa。
步骤2:联合MATLAB和ABAQUS建立纳米压痕的基于分段线性/幂律硬化材料模型的常规有限元模型,反求出分段线性/幂律硬化材料模型本构参数(屈服应力σy和应***化指数n)。
2-1:以半锥角为70.3°的圆锥形压头等效为步骤1-2中Berkovich三棱锥压头,用ABAQUS 建立纳米压痕的二维轴对称有限元模型。将压头正下方的材料网格进行局部细化,模型如图 3所示。采用位移控制加载的方式计算得到压头沿压入方向的接触反力和位移,输出接触力、试件接触面的接触压力和位移和压头最下方节点的位移,生成输入文件;
2-2:在MATLAB中运用拉丁超立方抽样抽取60组弹性模量E(在1-2计算出的数值附近选取),屈服应力σy和应***化指数n,计算得到各组抽样参数下的压痕的载荷-位移曲线和压痕堆积/沉陷参数s/h,并计算出模拟载荷-位移曲线与实验载荷-位移曲线的均方差,利用 MATLAB建立本构参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出实验材料的分段线性/幂律硬化材料模型本构参数 (屈服应力σy和应***化指数n),与运用Oliver-Pharr方法计算材料的弹性模量一起,记为C0;
2-3:将实验位移与堆积深度s相加,得到修正的的接触深度,进而得到修正的载荷-位移曲线,然后再利用2-2中的模拟载荷-位移曲线计算其与修正的载荷-位移曲线的均方差,重复 2-2中的单目标寻优过程,计算出修正后材料的分段线性/幂律硬化材料模型本构参数(屈服应力σy和应***化指数n),同时运用Oliver-Pharr方法计算材料的弹性模量,并记为C1;
2-4:计算2-2中计算出的材料参数C0与2-3中的C1之间的误差,如果误差在2%以内,则以2-3中计算的本构参数C1为材料的宏观本构参数;如果误差超过2%,则将2-3中修正的载荷-位移曲线作为实验载荷-位移曲线,重复2-2,2-3和2-4中的步骤,直至误差小于2%。此时获得304不锈钢的弹性模量E为196.12GPa,屈服应力σy为196MPa,应***化指数n为 0.251。
步骤3:联合MATLAB和ABAQUS,运用晶体塑性有限元建立拉伸试件的多晶体有限元模型,求解晶体塑性材料参数(初始屈服应力τ0,初始硬化模量h0,和阶段Ⅰ的应力τs) 与分段线性/幂律硬化材料模型本构参数(屈服应力σy和应***化指数n)的对应关系,进而反求待测材料的晶体塑性材料参数。
3-1:在ABAQUS中建立标准拉伸试件的有限元模型如图4,运用ABAQUS材料子程序赋予模型晶体塑性材料参数,进而建立标准拉伸试件的晶体塑性有限元模型。采用载荷控制加载的方式计算得到材料的应力应变曲线。生成输入文件;
3-2:在MATLAB中运用拉丁超立方抽样抽取60个初始硬化模量和饱和屈服应力的抽样点,修改3-1中输入文件中的材料参数,并计算得到各组抽样参数下的应力应变曲线,并计算出其与2-4中宏观材料参数下的应力应变曲线的均方差,利用MATLAB建立晶体塑性材料参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出待测材料的初始屈服应力τ0为86.11MPa,初始硬化模量h0为220.52MPa,和阶段Ⅰ的应力τs为256.35MPa。
步骤4:为了对反演得到的材料参数进行验证,对相同的304不锈钢材料进行了拉伸试验,得到的应力应变曲线与步骤2中获得的应力应变曲线对比如图5 。由此曲线计算得到的初始屈服应力τ0为85.09MPa,初始硬化模量h0为218.77MPa,和阶段Ⅰ的应力τs为260.52MPa。从对比结果可以看出,两种方法计算出的应力应变曲线和晶体塑性参数相差不大,该反演识别方法合理、有效、精度高,整个反演识别过程正确。
以上所述实施例仅表达本发明的实施方式,但并不能因此而理解为对本发明专利的范围的限制,应当指出,对于本领域的技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些均属于本发明的保护范围。
Claims (4)
1.一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法,其特征在于,所述的实现方法,首先,采用Oliver-Pharr方法获得材料弹性模量使得宏观本构参数反求模型简化;其次,使用分段线性/幂律硬化材料模型,联合MATLAB和ABAQUS建立纳米压痕宏观参数反演模型,利用堆积/沉陷参数对实际纳米压痕实验数据进行修正,结合Kriging代理模型和遗传算法计算出待测材料宏观本构参数;最后,建立拉伸试件基于晶体塑性的多晶体有限元模型,结合Kriging代理模型和遗传算法计算出实验材料晶体塑性材料参数。
2.根据权利要求1所述的一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法,其特征在于,具体包括以下步骤:
步骤1:待测金属材料的纳米压痕实验
1-1:将待测金属材料切割,运用机械抛光和振动抛光,得到符合要求的纳米压痕试件;
1-2:运用纳米压痕***对步骤1-1中的压痕试件进行压痕测试,获得包括载荷-位移曲线、最大载荷、接触刚度和接触硬度的实验压痕响应;运用Oliver-Pharr方法获得材料的弹性模量E;
步骤2:联合MATLAB和ABAQUS建立纳米压痕的基于分段线性/幂律硬化材料模型的常规有限元模型,反求出材料的宏观本构参数:屈服应力σy和应***化指数n;其中分段线性/幂律硬化材料模型描述为:
其中ε为总应变,σ为应力;
2-1:运用ABAQUS建立纳米压痕的二维轴对称有限元模型;采用位移控制加载的方式计算得到压头沿压入方向的接触反力和位移,输出接触力、试件接触面的接触压力和位移和压头最下方节点的位移,生成输入文件;
2-2:在MATLAB中运用拉丁超立方抽样抽取分段线性/幂律硬化材料模型本构参数,修改步骤2-1中输入文件中的材料参数,计算得到各组抽样参数下的压痕的载荷-位移曲线和压痕堆积/沉陷参数s/h,其中,s为堆积或沉陷高度,出现堆积现象时为s正,出现沉陷现象时为s负,h为压入深度;并计算出模拟载荷-位移曲线与实验载荷-位移曲线的均方差,建立本构参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出实验材料的分段线性/幂律硬化材料模型本构参数;
2-3:采用步骤2-2中的压痕堆积/沉陷参数修正实验载荷-位移曲线,得到修正的载荷-位移曲线,然后再利用步骤2-2中的模拟载荷-位移曲线计算其与修正的载荷-位移曲线的均方差,重复步骤2-2中的单目标寻优过程,计算出修正后材料的分段线性/幂律硬化材料模型本构参数;
2-4:计算步骤2-2中计算出的分段线性/幂律硬化材料模型本构参数与步骤2-3中修正的本构参数之间的误差,如果误差在允许范围之内,则以步骤2-3中修正的本构参数为材料的宏观本构参数;如果误差超过允许范围,则将步骤2-3中修正的载荷-位移曲线作为实验载荷-位移曲线,重复步骤2-2、步骤2-3和步骤2-4,直至误差在允许范围之内;
步骤3:联合MATLAB和ABAQUS,运用晶体塑性有限元建立拉伸试件的多晶体有限元模型,求解晶体塑性材料参数与分段线性/幂律硬化材料模型本构参数的对应关系,进而反求待测材料的晶体塑性材料参数;
3-1:在ABAQUS中建立标准拉伸试件的晶体塑性有限元模型;采用载荷控制加载的方式计算得到材料的应力应变曲线;生成输入文件;
3-2:在MATLAB中运用拉丁超立方抽样抽取晶体塑性材料参数,修改步骤3-1中输入文件中的材料参数,计算得到各组抽样参数下的应力应变曲线,并计算出模拟的应力应变曲线与步骤2-4中宏观材料参数下的应力应变曲线的均方差,利用MATLAB建立晶体塑性材料参数与均方差的Kriging代理模型,然后运用遗传算法,以两组数据的均方差最小为目标进行单目标寻优,计算出实验材料的晶体塑性材料参数。
3.根据权利要求2所述的一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法,其特征在于,步骤2中,材料参数反求的过程中将弹性参数和塑性参数分开,弹性参数运用成熟的理论方法求解,塑性参数运用有限元反演求解。
4.根据权利要求2所述的一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法,其特征在于,所述的步骤2-4中,误差的允许范围为0~2%。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911271059.3A CN111189699B (zh) | 2019-12-12 | 2019-12-12 | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 |
PCT/CN2020/127686 WO2021114994A1 (zh) | 2019-12-12 | 2020-11-10 | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 |
US17/281,032 US20210310917A1 (en) | 2019-12-12 | 2020-11-10 | Inversion identification method of crystal plasticity material parameters based on nanoindentation experiments |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911271059.3A CN111189699B (zh) | 2019-12-12 | 2019-12-12 | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111189699A CN111189699A (zh) | 2020-05-22 |
CN111189699B true CN111189699B (zh) | 2021-05-18 |
Family
ID=70705773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911271059.3A Active CN111189699B (zh) | 2019-12-12 | 2019-12-12 | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20210310917A1 (zh) |
CN (1) | CN111189699B (zh) |
WO (1) | WO2021114994A1 (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111189699B (zh) * | 2019-12-12 | 2021-05-18 | 大连理工大学 | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 |
CN112651153B (zh) * | 2020-12-09 | 2024-04-09 | 中南大学 | 一种确定晶体塑性有限元模型材料参数的方法 |
CN113049372B (zh) * | 2020-12-23 | 2024-04-02 | 河南工业大学 | 一种计算塑性混凝土弹性模量的方法 |
CN113435077B (zh) * | 2021-05-24 | 2023-08-22 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种紧固件连接结构中强度参数确定方法 |
CN113670751A (zh) * | 2021-07-30 | 2021-11-19 | 中国航发北京航空材料研究院 | 一种测量多相多孔隙材料弹性模量的方法 |
CN113987763B (zh) * | 2021-10-12 | 2024-04-12 | 广州特种设备检测研究院(广州市特种设备事故调查技术中心、广州市电梯安全运行监控中心) | 一种聚乙烯双曲线本构模型的构建方法及应用 |
CN114397210B (zh) * | 2021-12-28 | 2023-07-07 | 西安电子科技大学 | 一种基于球形压痕测量材料各向异性塑性的贝叶斯方法 |
CN114323951B (zh) * | 2021-12-31 | 2023-07-28 | 华侨大学 | 基于纳米压痕技术的材料本构预模型建立方法 |
CN114492122B (zh) * | 2022-01-12 | 2024-02-23 | 西北工业大学 | 一种数值表征喷丸处理材料表层应***化梯度的方法 |
CN114444367B (zh) * | 2022-01-26 | 2024-05-28 | 中国人民解放军63653部队 | 干沙动力学仿真的离散元参数自动标定方法 |
GB2615332B (en) * | 2022-02-03 | 2024-04-03 | Plastometrex Ltd | Indentation plastometry |
CN114894647B (zh) * | 2022-04-30 | 2024-05-10 | 上海核工程研究设计院股份有限公司 | 一种合金预生膜弹性模量测试方法及其应用 |
CN115472248B (zh) * | 2022-09-23 | 2023-06-13 | 哈尔滨工业大学 | 一种CuZrAl非晶合金纳米压痕测试的分子动力学模拟计算方法 |
CN115831289A (zh) * | 2022-11-29 | 2023-03-21 | 东风华神汽车有限公司 | 一种材料防护结构碰撞性能判断方法、装置及设备 |
CN115879331B (zh) * | 2023-03-01 | 2023-05-19 | 山东科技大学 | 基于开尔文模型的弹簧-阻尼减振结构参数优化分析方法 |
CN116629010B (zh) * | 2023-06-02 | 2024-01-23 | 江苏科技大学 | 一种基于随机过程的退化模型确认及其试验设计方法 |
CN117497069B (zh) * | 2023-10-23 | 2024-05-24 | 华中科技大学 | 一种高聚物材料的超弹性本构模型的构建方法和装置 |
CN117409898B (zh) * | 2023-10-24 | 2024-03-19 | 延边大学 | 一种复合多晶超硬材料及其制备方法 |
CN117877642B (zh) * | 2024-01-16 | 2024-07-23 | 华中科技大学 | 一种焊缝晶体塑性模型本构参数标定方法及*** |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100702422B1 (ko) * | 2005-05-24 | 2007-04-03 | 부산대학교 산학협력단 | 전기아연도금강판의 코팅층에 대한 기계적 특성 평가 방법 |
CN102313677A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊缝不同区域材料动态力学性能参数的检测方法 |
CN102313676A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊点不同区域材料静态力学性能参数的检测方法 |
CN103091185A (zh) * | 2011-10-31 | 2013-05-08 | 湘潭大学 | 一种利用压痕法表征金属材料弹塑性力学性能的方法 |
CN104915478A (zh) * | 2015-05-24 | 2015-09-16 | 浙江大学 | 基于多参数不确定性分析的产品设计模型等效简化方法 |
CN105675419A (zh) * | 2016-01-14 | 2016-06-15 | 西南交通大学 | 双锥形压入预测材料单轴本构关系测定方法 |
CN106407614A (zh) * | 2016-11-08 | 2017-02-15 | 武汉大学 | 一种结合纳米压痕与有限元获取焊缝组织力学参数的方法 |
CN108414379A (zh) * | 2018-03-16 | 2018-08-17 | 太原理工大学 | 一种原位压入测试提取金属弹塑性参数的方法 |
CN108645704A (zh) * | 2018-03-05 | 2018-10-12 | 上海交通大学 | 基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法 |
CN109870376A (zh) * | 2019-02-03 | 2019-06-11 | 浙江大学 | 一种基于纳米压痕和数值模拟反演岩石矿物参数的方法 |
CN110261247A (zh) * | 2019-06-20 | 2019-09-20 | 西北工业大学 | 金属材料各向异性屈服及硬化本构参数同步表征方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109900570B (zh) * | 2019-03-28 | 2020-10-27 | 西安交通大学 | 一种基于纳米压入技术的非晶薄膜塑性变形表征方法 |
CN111189699B (zh) * | 2019-12-12 | 2021-05-18 | 大连理工大学 | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 |
-
2019
- 2019-12-12 CN CN201911271059.3A patent/CN111189699B/zh active Active
-
2020
- 2020-11-10 US US17/281,032 patent/US20210310917A1/en not_active Abandoned
- 2020-11-10 WO PCT/CN2020/127686 patent/WO2021114994A1/zh active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100702422B1 (ko) * | 2005-05-24 | 2007-04-03 | 부산대학교 산학협력단 | 전기아연도금강판의 코팅층에 대한 기계적 특성 평가 방법 |
CN102313677A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊缝不同区域材料动态力学性能参数的检测方法 |
CN102313676A (zh) * | 2011-08-31 | 2012-01-11 | 湖南大学 | 一种识别焊点不同区域材料静态力学性能参数的检测方法 |
CN103091185A (zh) * | 2011-10-31 | 2013-05-08 | 湘潭大学 | 一种利用压痕法表征金属材料弹塑性力学性能的方法 |
CN104915478A (zh) * | 2015-05-24 | 2015-09-16 | 浙江大学 | 基于多参数不确定性分析的产品设计模型等效简化方法 |
CN105675419A (zh) * | 2016-01-14 | 2016-06-15 | 西南交通大学 | 双锥形压入预测材料单轴本构关系测定方法 |
CN106407614A (zh) * | 2016-11-08 | 2017-02-15 | 武汉大学 | 一种结合纳米压痕与有限元获取焊缝组织力学参数的方法 |
CN108645704A (zh) * | 2018-03-05 | 2018-10-12 | 上海交通大学 | 基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法 |
CN108414379A (zh) * | 2018-03-16 | 2018-08-17 | 太原理工大学 | 一种原位压入测试提取金属弹塑性参数的方法 |
CN109870376A (zh) * | 2019-02-03 | 2019-06-11 | 浙江大学 | 一种基于纳米压痕和数值模拟反演岩石矿物参数的方法 |
CN110261247A (zh) * | 2019-06-20 | 2019-09-20 | 西北工业大学 | 金属材料各向异性屈服及硬化本构参数同步表征方法 |
Non-Patent Citations (4)
Title |
---|
A procedure to prevent pile up effects on the analysis of spherical indentation data in elastic–plastic materials.;J. Rodriguez等;《Mechanics of Materials》;20070308;第39卷(第11期);987-997 * |
Development of pile-up during spherical indentation of elastic–plastic solids;B. Taljat等;《International Journal of Solids and Structures》;20040324;第41卷(第14期);3891-3904 * |
Kriging模型及代理优化算法研究进展;韩忠华;《航空学报》;20160329;第37卷(第11期);3197-3225 * |
基于纳米压痕技术和有限元仿真的材料力学性能分析;刘扬等;《武汉理工大学学报(交通科学与工程版)》;20031031;第27卷(第5期);690-693 * |
Also Published As
Publication number | Publication date |
---|---|
WO2021114994A1 (zh) | 2021-06-17 |
CN111189699A (zh) | 2020-05-22 |
US20210310917A1 (en) | 2021-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111189699B (zh) | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 | |
CN103792143B (zh) | 一种单轴拉伸全程真应力应变曲线的快速获取方法 | |
KR101065502B1 (ko) | 파단 예측 방법 | |
CN108645704B (zh) | 基于纳米压痕和有限元模拟反演标定金属材料微观本构参数的方法 | |
CN109900554B (zh) | 一种利用压痕法计算断裂韧性的方法 | |
CN111339703A (zh) | 大应变条件下材料硬化行为的虚拟预测方法 | |
CN105260574A (zh) | 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法 | |
CN111666699B (zh) | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 | |
CN111339594B (zh) | 基于dic技术的近场动力学参数实验反演***及使用方法 | |
CN111855458B (zh) | 一种基于纳米压痕理论的多孔材料本构关系求解方法 | |
Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
CN106407614A (zh) | 一种结合纳米压痕与有限元获取焊缝组织力学参数的方法 | |
CN109948216B (zh) | 总应变能密度修正的缺口件低周疲劳预测方法 | |
CN112651153B (zh) | 一种确定晶体塑性有限元模型材料参数的方法 | |
CN109446658A (zh) | 一种钛合金锥形压痕约束因子和代表应变的确定方法 | |
CN113611377A (zh) | 一种利用晶体塑性模型模拟混合控制蠕变疲劳变形的方法 | |
CN108732048A (zh) | 一种级配碎石重复荷载试验的弹塑性屈服点应力确定方法 | |
Ramos et al. | Cracking study of a reinforced concrete beam | |
KR20160088756A (ko) | 유한요소해에 기초한 경화재의 비커스 압입파괴인성 평가방법 | |
Karpenko et al. | Incremental approach to the nonlinear analysis of reinforcement concrete with cracks at plane stress state | |
CN109992853B (zh) | 一种表面强化金属部件的残余应力场数值重建方法 | |
CN110501177B (zh) | 基于自由端倾角影响线曲率的悬臂梁损伤识别方法 | |
CN109459329B (zh) | 一种预加载条件下钛合金锥形压痕代表应变的确定方法 | |
Giusca et al. | Digital representation of a load cell | |
CN109214081B (zh) | 一种基于数值模拟的7005铝合金微观组织分析***及方法 |
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 |