CN108133082B - 一种基于有限元模拟确定压痕应变法中应力测量常数的方法 - Google Patents
一种基于有限元模拟确定压痕应变法中应力测量常数的方法 Download PDFInfo
- Publication number
- CN108133082B CN108133082B CN201711277423.8A CN201711277423A CN108133082B CN 108133082 B CN108133082 B CN 108133082B CN 201711277423 A CN201711277423 A CN 201711277423A CN 108133082 B CN108133082 B CN 108133082B
- Authority
- CN
- China
- Prior art keywords
- strain
- calibration
- stress
- finite element
- stress measurement
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000007373 indentation Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 36
- 238000005259 measurement Methods 0.000 title claims abstract description 35
- 239000000463 material Substances 0.000 claims abstract description 26
- 238000002474 experimental method Methods 0.000 claims abstract description 16
- 238000011088 calibration curve Methods 0.000 claims abstract description 13
- 239000007769 metal material Substances 0.000 claims abstract description 9
- 238000005516 engineering process Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 abstract description 14
- 238000009662 stress testing Methods 0.000 abstract description 3
- 229910000838 Al alloy Inorganic materials 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 229910052751 metal Inorganic materials 0.000 description 3
- 229910000851 Alloy steel Inorganic materials 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 240000005373 Panax quinquefolius Species 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000005483 Hooke's law Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
-
- 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/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0676—Force, weight, load, energy, speed or acceleration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于有限元模拟确定压痕应变法应力测量常数的方法,属于残余应力测试技术领域。应力测量常数是压痕应变增量与弹性应变的关系系数,一般通过实验标定获得,要求标定试板为无应力状态并具有足够大的尺寸,这给标定实验带来很大难度。该方法首先通过实验获得无应力材料中的压痕应变增量;然后建立有限元模拟标定模型,以零应力下的应变增量为依据,来确定相应材料的压头载荷;最后计算各弹性应变下的应变增量;然后以弹性应变为横坐标、应变增量为纵坐标绘制标定曲线,通过标定曲线实现应力测量常数的非实验标定。本发明能够获得更高准确度和更好规律性的标定曲线,可以作为大多数金属材料的非实验标定方法。
Description
技术领域
本发明涉及残余应力测试技术领域,具体涉及一种基于有限元模拟确定压痕应变法中应力测量常数的方法。
背景技术
压痕应变法是根据压痕诱导获得的应变变量信息(称为压痕叠加应变增量,简称应变增量),通过理论分析和数据归纳,最终利用胡克定律而求出构件中残余应力的一种几乎无损的应力测试方法。从测试的准确性出发,它本身只是压痕引起的应变变化,是一个宏观的不依赖于材料微观组织的变化量,只与构件中的应力水平有关,所以准确度只和测试者的水平和应力测量系数的标定结果有关。
应力测量常数(即压痕应变增量与弹性应变的关系系数)的确定是决定应力测量准确与否的关键。应力测量常数的确定都是采用实验标定方法获得,其规律见图1所示。所用的标定试板须要无应力状态,即标定前要对试板进行消除应力处理,试板的应力状态直接影响标定结果的准确性和稳定性。同时,为满足加载条件,试板尺寸也有相应要求,对于无法提供符合条件的标定试板,标定试验无法完成。这就给一些材料是否能够采用压痕应变法测量残余应力带来了很大的局限性。
引用有限元数值模拟技术,可以实现应力测量常数的非实验标定。
现有技术中曾采用模拟计算的方法对特定材料的标定结果进行过研究[文献1:孙渊,王庆明,夏风芳,等.残余应力测量法中压痕标定实验的分析[J].机械制造,2006:44:70-72;文献2:曲鹏程.屈服强度对压痕应变法测量残余应力的影响[D].沈阳:中国科学院金属研究所,2006;文献3:孟宪陆.不同应力场中压痕应变法的数值模拟[D].沈阳:中国科学院金属研究所,2007;文献4:刘生.材料性能和方向性对压痕应变法测量残余应力的影响[D].沈阳:中国科学院金属研究所,2013.],并获得了一些成果。但是,鉴于模型的局限性,有关结论并不能用于其它类型的材料。其它关于压痕应变法应力测量常数的确定未见报道,且目前没有与本发明相关的专利。本发明旨在建立可以准确获得大部分金属材料标定曲线的通用模型,解决一些金属材料不能实施实验标定和实验标定中误差较大的问题,以提高压痕应变法测试残余应力的精确性和扩大其应用范围。
发明内容
本发明的目的在于提供一种基于有限元模拟确定压痕应变法中应力测量常数的方法,本发明采用有限元数值模拟方法,进行压痕应变法应力测量常数的非实验标定,能够获得大多数金属材料的应力测量常数,解决实验标定中的各种困难和可能误差,同时该方法较之实验标定结果,规律性更强,准确度更高。
为实现上述目的,本发明所采用的技术方案如下:
一种基于有限元模拟确定压痕应变法应力测量常数的方法,该方法是采用数值模拟技术实现金属材料应力测量常数的非实验标定,具体包括如下步骤:
(1)通过实验获得无应力材料中的压痕应变增量;
(2)建立有限元模拟标定模型,以步骤(1)实验中获得的零应力下的应变增量为依据,并根据有限元模拟标定模型,来确定相应材料的压头载荷;
(3)基于该材料的拉伸曲线,计算各弹性应变下的应变增量;然后以弹性应变εe为横坐标、应变增量Δε为纵坐标绘制标定曲线,通过标定曲线实现应力测量常数的非实验标定。
上述步骤(3)中,通过该材料的拉伸曲线,获得初始屈服强度σy和屈服应变εy,并进一步计算各弹性应变下的应变增量。
上述步骤(3)中,应力测量常数的非实验标定过程为:在标定曲线中,将施加单向拉伸或压缩应变方向的各点按3次方进行拟合,在高值单向拉伸应力情况下,横向应变输出偏离3次方拟和曲线的各点按直线拟合,由拟合后得到的两个方程计算出的方程系数即为压痕应变法的应力测量常数。
上述步骤(2)中,所建立的有限元模拟标定模型尺寸为48mm×48mm×4mm,最小单元尺寸为0.04mm×0.03mm×0.004mm。
本发明具有如下优点和有益效果:
1、本发明方法能够准确确定大多数金属材料应力测量常数的非实验标定方法,该方法解决了实验标定中的各种困难和可能误差(包括实验误差、标定试板准备、试板零应力处理、标定装置的可靠性等)。
2、本发明通过引入有限元数值模拟技术,通过对动态压入材料的特性分析,以无应力下的应变增量实验数据为依据,基于传统拉伸曲线,针对不同材料动态确定合理的压入载荷,建立可以准确确定应力测量常数的标定模型。模拟标定结果和实验标定结果进行对比验证,证明本发明具有更好的规律性和更高的准确性,可以作为大多数金属材料的非实验标定方法,能够获得更高准确度和更好规律性的标定曲线。
附图说明
图1为实验标定曲线示意图。
图2为模拟标定模型。
图3为铝合金7N01的模拟和标定结果对比。
图4为低合金钢Q345的模拟和标定结果对比。
具体实施方式
以下结合实施例和附图详述本发明,其中实验标定部分按照GB/T24179-2009《金属材料残余应力测定压痕应变法》进行。
本发明为基于有限元模拟确定压痕应变法应力测量常数的方法,该方法是采用数值模拟技术实现金属材料应力测量常数的非实验标定,具体过程为:首先,通过实验获得材料零应力下的应变增量;然后采用图2所示的数值模拟标定模型,以零应力下的实验数据为依据,来确定相应材料的压头载荷;最后,计算各弹性应变下的应变增量,获得标定曲线,实现应力测量常数的非实验确定。
本发明方法中,模型尺寸和单元大小应适当才能保证计算精度,并降低运算成本。本模型采用48×48×4mm,最小单元尺寸为0.04×0.03×0.004mm。
实施例1:
首先获得铝合金7N01材料的拉伸性能曲线,得到初始屈服强度σy和屈服应变εy。通过实验获得无应力材料中的压痕应变增量。建立有限元模拟标定模型,压头直径φ1.588mm。以实验得到的无应力下的应变增量为依据,确定压头静载压入的载荷大小。相应的,获得在-0.9εy~0.9εy区间单向拉伸或压缩条件下任意弹性应变和应变增量的关系。将施加单向拉伸或压缩应变方向的结果按3次方进行拟合,与其垂直方向的结果若偏离3次方拟和曲线可按直线拟合,获得的方程系数即为压痕应变法的应力测量常数。
图3为铝合金7N01材料的有限元模拟结果和实验结果的对比图(仅以加载方向的模拟曲线和实验数据示例)。对比结果显示,模拟标定结果和实验标定结果呈现很好的一致性。由于实验时标定试板状态的影响和操作时不可避免的误差,导致实验标定数据较为分散,而模拟标定数据表现出了更好的规律性。
实施例2:
首先获得Q345低合金钢材料的拉伸性能曲线,得到初始屈服强度σy和屈服应变εy。通过实验获得无应力材料中的压痕应变增量。建立有限元模拟标定模型,压头直径1.588mm。以实验得到的无应力下的应变增量为依据,确定压头静载压入的载荷大小。相应的,获得在-0.9εy~0.9εy区间单向拉伸或压缩条件下任意弹性应变和应变增量的关系。将施加单向拉伸或压缩应变方向的结果按3次方进行拟合,与其垂直方向的结果若偏离3次方拟和曲线可按直线拟合,获得的方程系数即为压痕应变法的应力测量常数。
图4为Q345材料的有限元模拟结果和实验结果的对比图(仅以加载方向的模拟曲线和实验数据示例)。对比结果显示,模拟标定结果和实验标定结果呈现很好的一致性。由于实验时标定试板状态的影响和操作时不可避免的误差,导致实验标定数据较为分散,而模拟标定数据表现出了更好的规律性。
Claims (3)
1.一种基于有限元模拟确定压痕应变法应力测量常数的方法,其特征在于:该方法是采用数值模拟技术实现金属材料应力测量常数的非实验标定,具体包括如下步骤:
(1)通过实验获得无应力材料中的压痕应变增量;
(2)建立有限元模拟标定模型,以步骤(1)实验中获得的零应力下
的应变增量为依据,并根据有限元模拟标定模型,来确定相应材料的压头载荷;
(3)基于该材料的拉伸曲线,计算各弹性应变下的应变增量;然后以弹性应变为横坐标、应变增量为纵坐标绘制标定曲线,通过标定曲线实现应力测量常数的非实验标定;应力测量常数的非实验标定过程为:在标定曲线中,将施加单向拉伸或压缩应变方向的各点按3次方进行拟合;在高值单向拉伸应力情况下,横向应变输出偏离3 次方拟和曲线的各点按直线拟合,由拟合后得到的两个方程计算出的方程系数即为压痕应变法的应力测量常数。
2.根据权利要求1所述的基于有限元模拟确定压痕应变法应力测量常数的方法,其特征在于:步骤(3)中,通过该材料的拉伸曲线,获得初始屈服强度和屈服应变,并进一步计算各弹性应变下的应变增量。
3.根据权利要求1所述的基于有限元模拟确定压痕应变法应力测量常数的方法,其特征在于:步骤(2)中,所建立的有限元模拟标定模型尺寸为48mm×48mm×4mm,最小单元尺寸为0.04mm×0.03mm×0.004mm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711277423.8A CN108133082B (zh) | 2017-12-06 | 2017-12-06 | 一种基于有限元模拟确定压痕应变法中应力测量常数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711277423.8A CN108133082B (zh) | 2017-12-06 | 2017-12-06 | 一种基于有限元模拟确定压痕应变法中应力测量常数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108133082A CN108133082A (zh) | 2018-06-08 |
CN108133082B true CN108133082B (zh) | 2021-04-20 |
Family
ID=62389999
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711277423.8A Active CN108133082B (zh) | 2017-12-06 | 2017-12-06 | 一种基于有限元模拟确定压痕应变法中应力测量常数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108133082B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086529B (zh) * | 2018-08-07 | 2023-05-30 | 中国科学院金属研究所 | 一种基于零压力下应变增量确定压痕应变法中应力计算函数的方法 |
CN109145417B (zh) * | 2018-08-07 | 2023-01-17 | 中国科学院金属研究所 | 一种基于材料力学性能直接确定压痕应变法应力计算函数的方法 |
CN109470564B (zh) * | 2018-11-28 | 2021-08-06 | 南京航空航天大学 | 夹持段打滑后复合材料拉伸试验夹紧力测试装置及测试方法 |
CN111141437B (zh) * | 2020-01-02 | 2021-12-21 | 中车青岛四方机车车辆股份有限公司 | 电阻点焊接头内部残余应力测量方法 |
CN111366283A (zh) * | 2020-03-30 | 2020-07-03 | 中国科学院金属研究所 | 一种辅助压痕应变法检测焊接结构件残余应力的方法 |
CN111735728B (zh) * | 2020-05-22 | 2023-03-31 | 扬州市舜意机械有限公司 | 一种减振涂层力学性能识别方法及装置 |
CN113297772B (zh) * | 2021-06-07 | 2022-04-01 | 南昌航空大学 | 一种使用有限元软件模拟复合材料残余应力分布的方法 |
CN116577190B (zh) * | 2023-07-13 | 2023-11-28 | 南京华建检测技术有限公司 | 一种t型实验试块的智能化检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7992449B1 (en) * | 2008-02-15 | 2011-08-09 | Mahmoud Khaled M | Method for assessment of cable strength and residual life |
CN102288499A (zh) * | 2011-08-31 | 2011-12-21 | 湖南大学 | 一种识别焊缝不同区域材料静态力学性能参数的检测方法 |
KR20130123699A (ko) * | 2012-05-03 | 2013-11-13 | 한국표준과학연구원 | 압흔영상분석을 기반으로 하는 구형 압흔의 유동곡선 획득방법 |
CN103528887A (zh) * | 2013-10-24 | 2014-01-22 | 吉林大学 | 原位拉/压-扭转复合载荷材料微观力学测试平台 |
CN106202631A (zh) * | 2016-06-28 | 2016-12-07 | 湖南大学 | 一种冲压成形高强钢的材料参数获取方法 |
CN106326567A (zh) * | 2016-08-25 | 2017-01-11 | 武汉大学 | Abaqus中钢‑聚丙烯混杂纤维混凝土塑性本构参数的取值方法 |
CN107024401A (zh) * | 2017-04-20 | 2017-08-08 | 西北工业大学 | 获取金属材料各向异性和拉压非对称性的方法及*** |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
US20130192384A1 (en) * | 2011-08-15 | 2013-08-01 | E I Du Pont De Nemours And Company | Method for testing tensile strength of an electrically nonconductive material |
CN105512400A (zh) * | 2015-12-09 | 2016-04-20 | 大连理工大学 | 一种脆性材料切削过程仿真方法 |
-
2017
- 2017-12-06 CN CN201711277423.8A patent/CN108133082B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7992449B1 (en) * | 2008-02-15 | 2011-08-09 | Mahmoud Khaled M | Method for assessment of cable strength and residual life |
CN102288499A (zh) * | 2011-08-31 | 2011-12-21 | 湖南大学 | 一种识别焊缝不同区域材料静态力学性能参数的检测方法 |
KR20130123699A (ko) * | 2012-05-03 | 2013-11-13 | 한국표준과학연구원 | 압흔영상분석을 기반으로 하는 구형 압흔의 유동곡선 획득방법 |
CN103528887A (zh) * | 2013-10-24 | 2014-01-22 | 吉林大学 | 原位拉/压-扭转复合载荷材料微观力学测试平台 |
CN106202631A (zh) * | 2016-06-28 | 2016-12-07 | 湖南大学 | 一种冲压成形高强钢的材料参数获取方法 |
CN106326567A (zh) * | 2016-08-25 | 2017-01-11 | 武汉大学 | Abaqus中钢‑聚丙烯混杂纤维混凝土塑性本构参数的取值方法 |
CN107024401A (zh) * | 2017-04-20 | 2017-08-08 | 西北工业大学 | 获取金属材料各向异性和拉压非对称性的方法及*** |
Non-Patent Citations (3)
Title |
---|
Determination of equi-biaxial residual stress and plastic properties in structural steel using instrumented indentation;Pham TH 等;《Materials Science & Engineering A》;20170314;第688卷;第352-363页 * |
三点弯曲标定试验中双向主应变之间的关系;陈静 等;《兵器材料科学与工程》;20081130;第31卷(第6期);第37-40页 * |
金属材料残余应力的测定方法;巴发海 等;《理化检验-物理分册》;20171108;第53卷(第11期);第771-777页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108133082A (zh) | 2018-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108133082B (zh) | 一种基于有限元模拟确定压痕应变法中应力测量常数的方法 | |
CN108982223B (zh) | 金属圆棒试样单轴拉伸试验大应变范围应力应变的测量方法 | |
CN108982222B (zh) | 金属板试样单轴拉伸试验大应变范围应力应变的测量方法 | |
Wang et al. | An Experimental‐Numerical Combined Method to Determine the True Constitutive Relation of Tensile Specimens after Necking | |
CN110717298B (zh) | 一种积分型静力渐进损伤破坏模型的试验数据处理方法 | |
CN108108530A (zh) | 一种适用于结构连接件的疲劳寿命校准方法 | |
CN108844824B (zh) | 一种基于圆锥形压头的已知材料残余应力测定方法 | |
CN110927054B (zh) | 金属孔结构的冷挤压强化摩擦系数测试方法 | |
CN108645706A (zh) | 一种通过硬化强度和抗拉强度预测金属材料疲劳强度的方法 | |
Zafer et al. | Design and development of strain gauge pressure transducer working in high pressure range of 500 MPa using autofrettage and finite element method | |
CN108548720B (zh) | I型裂纹弹塑性理论公式获取延性材料j阻力曲线的方法 | |
CN109446658B (zh) | 一种钛合金锥形压痕约束因子和代表应变的确定方法 | |
CN104237037B (zh) | 基于Berkovich压痕的材料弹塑性参数仪器化压入测试方法 | |
CN103047939A (zh) | 光纤光栅应变传感器工程适用性评价方法 | |
CN109145417B (zh) | 一种基于材料力学性能直接确定压痕应变法应力计算函数的方法 | |
CN110618029A (zh) | 极限载荷极短时间工况下滚珠丝杠副疲劳弹性寿命测试方法 | |
CN109490334B (zh) | 一种运用残余应力预测模型的t字型锻件无损测试方法 | |
CN115479853A (zh) | 一种砂土循环动力响应预测方法 | |
CN104077444B (zh) | 一种压痕数据分析方法 | |
CN111623909B (zh) | 一种双叉臂二力杆的受力测试方法 | |
Chen et al. | A novel method to determine full-range hardening curve for metal bar using hyperbolic shaped compression specimen | |
KR20230006308A (ko) | 금형수명 예측 시스템 및 방법 | |
CN105808846A (zh) | 基于弹簧模型的干涉配合定量计算方法 | |
CN109086529B (zh) | 一种基于零压力下应变增量确定压痕应变法中应力计算函数的方法 | |
CN113218571A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |