CN111896373A - 一种测定等效塑性应变成形极限图的试验与计算方法 - Google Patents

一种测定等效塑性应变成形极限图的试验与计算方法 Download PDF

Info

Publication number
CN111896373A
CN111896373A CN202010618096.3A CN202010618096A CN111896373A CN 111896373 A CN111896373 A CN 111896373A CN 202010618096 A CN202010618096 A CN 202010618096A CN 111896373 A CN111896373 A CN 111896373A
Authority
CN
China
Prior art keywords
strain
forming limit
true
uniform
equivalent plastic
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.)
Granted
Application number
CN202010618096.3A
Other languages
English (en)
Other versions
CN111896373B (zh
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.)
Wuhan Shangshan Simulation Technology Co ltd
Original Assignee
Wuhan Shangshan Simulation Technology 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 Wuhan Shangshan Simulation Technology Co ltd filed Critical Wuhan Shangshan Simulation Technology Co ltd
Priority to CN202010618096.3A priority Critical patent/CN111896373B/zh
Publication of CN111896373A publication Critical patent/CN111896373A/zh
Application granted granted Critical
Publication of CN111896373B publication Critical patent/CN111896373B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • 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/006—Crack, flaws, fracture or rupture
    • G01N2203/0067—Fracture or rupture
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General 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)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种测定等效塑性应变成形极限图的试验与计算方法，其步骤包括：测定材料的真实均匀应变和真实断裂应变；计算不同简单应变路径下的主应变和次应变；绘制并修正等效塑性应变成形极限图(EPS‑FLD)。本发明通过理论创新避免了现有技术的各种不足，为新形势下的汽车用钢提供了技术解决方案，同时，将材料均匀成形极限和断裂成形极限的测定集成在单轴拉伸DIC试验中，在无需额外增加试验并可达到更好的技术目标的条件下，无需额外增加试验，节省了宝贵的人力成本和设备成本，因此，本发明具有重要的理论与工程实践意义。

Description

一种测定等效塑性应变成形极限图的试验与计算方法
技术领域
本发明涉及一种测定成形极限图的方法，尤其是涉及金属材料的等效塑性应变成形极限图的试验与计算及其修正方法。
背景技术
成形极限图(Forming Limit Diagram，简称FLD)于上个世纪60年代由Keeler(SAEtechnical paper 650535)和Goodwin(SAE technical paper 680093)共同提出。Keeler提出了右半边区域用于板料受双拉变形的简单应变状态区域的分析，Goodwin在Keeler的基础上扩展了左半边区域用于板料受拉-压变形的简单应变状态区域的分析。
FLD描述了在不同简单应变路径下的主应变(纵坐标)和次应变(横坐标)临近局部颈缩状态的极限组合，其简单应变路径覆盖从纯单轴拉伸到纯双轴拉伸状态，对应的简单应变路径β的取值范围为[-0.5,1]。
目前，成形极限图已经被工业界认为是评价板料成形性的有效工具，形成了成熟的行业标准，如GB-T/15825.8-2008(《金属薄板成形性能与试验方法第8部分:成形极限图(FLD)测定指南》)，该标准针对如何测定板料的FLD制定了详细的试验方法和流程。简单说，首先，制作几种简单应变路径下的试验试样并进行试验；其次，试验完成之后，选取并测定临界网格圆的主应变和次应变，在成形极限图中标注测量点；最后，根据测量点的分布状态，拟合得到板料的成形极限曲线(Forming Limit Curve，简称FLC)，其形状一般为“V”字型。
在标准的成形极限图中，测量主应变和次应变的临界网格圆的选取是根据是否处于发生局部颈缩的临界状态为判断标准，具体在工程上的处理方法是：将位于缩颈区但未破裂的网格圆作为临界网格圆，或者，将紧靠缩颈或裂纹的网格圆作为临界网格圆，或者，将与缩颈或裂纹横贯其中部之网格圆相邻的网格圆作为临界网格圆，因此，FLD不能用于在冲压分析中判断材料失效的标准。
然而，在标准的成形极限图中，若测量主应变和次应变的临界网格圆的选取根据是否处于发生断裂的临界状态为判断标准，则拟合得到板料的断裂成形极限(FractureForming Limit，简称FFL)，若简单应变路径β的取值范围为[0,1]，根据体积不变条件：
ε1+ε2+ε3＝0 (1)
在断裂临界状态，厚度方向的应变近似等于常数ε3≈const.，则(1)式可变为：
ε1＝-ε2+const. (2)
因此，在成形极限图中，FFL近似为一条斜率为“-1”的直线。
正是由于FLD为方便研究板料成形极限、评价板料成形性能及解决板料成形领域中众多难题提供了技术基础和实用判据，因此，FLD被集成到了商用冲压仿真软件中用于板料成形严重度的判断标准。
但是，本申请的发明人在实现本发明的过程中，经过研究发现：传统的FLD判据标准在测试方法上存在技术缺陷：
第一，不同简单应变路径下测定的主应变和次应变，对临界网格圆的变形程度判断，难以保证在同一个变形水平，FLC曲线在形式上不符合理论预测结果。
第二，试验测量的大变形条件下的主应变和次应变，是在特定尺寸下测量的，两者之比并不一定等于该试验测试条件下对应的β，其值的准确性严重依赖于网格尺寸，一般而言，只有在网格尺寸等于1.0mm时，主应变与次应变之比，才能与β值严格对应。
第三，需要通过复杂的试验得到，且试验结果受许多因素的影响，如测试工具、试样的料厚、n值等。
此外，随着先进高强钢(AHSS)的大量应用，传统的FLD已经被反复证明不能预测AHSS的冲压成形(SAE Technical Paper 2014-01-0985)，特别是针对AHSS常见的边缘开裂和局部失效问题，传统的FLD判据标准完全无能为力。
因此，由于现有技术在测试方法和工程应用中存在的各种不足，使得传统的FLD受到了质疑，一直被不断的改进和发展。本发明正是对传统FLD的补充和完善，所依据的理论基础如下：
针对临界网格圆，根据体积不变条件，则有：
ε1+ε2+ε3＝0 (3)
其中，ε1、ε2和ε3分别为三个方向的主应变，同时，r值的定义如下：
Figure BDA0002561774170000021
则ε2和ε3采用ε1、r值表示如下：
Figure BDA0002561774170000022
Figure BDA0002561774170000023
在单轴拉伸状态下，等效应变
Figure BDA0002561774170000031
表示如下：
Figure BDA0002561774170000032
将(5)和(6)式带入(7)式，则有：
Figure BDA0002561774170000033
采用
Figure BDA0002561774170000034
和r值表示ε1、ε2和ε3如下：
Figure BDA0002561774170000035
Figure BDA0002561774170000036
Figure BDA0002561774170000037
根据(9)、(10)和(11)式，临界网格圆的应变状态(ε1、ε2和ε3)的平方和为常数，如下：
Figure BDA0002561774170000038
其物理意义表示：基于应变状态表征的材料屈服面或后继屈服面，其应变强度或等效应变是一不变的定值，该定值只取决于临界网格圆的变形程度，而与应变路径或应变状态无关，根据(12)式，材料屈服面或后继屈服面是半径为
Figure BDA0002561774170000039
的标准球面。根据(3)、(7)和(12)式，表征材料采用了等向硬化模型。
为了便于工程应用，可将(12)式表征的三维屈服面降为二维表达，将(3)式带入(7)式，则有：
Figure BDA00025617741700000310
其物理意义表示：以主应变ε1为纵坐标，次应变ε2为横坐标，材料屈服面或后继屈服面的形状为一标准椭圆，并将(9)和(10)式带入(13)式中，(13)式同样得到满足。
根据(9)和(10)式，式中r值表示简单应变状态的取值范围，
Figure BDA00025617741700000311
表示特定的变形程度，只要测量任意简单应变路径下的等效应变，则可以计算出具有相同变形水平的其它简单应变路径下的主应变和次应变，因此，在理论上可绘制出材料的成形极限图FLD。比如，从单轴拉伸试验中，根据试样均匀颈缩极限确定等效均匀应变
Figure BDA0002561774170000041
针对冲压成形，r值的取值范围一般为[-0.5,1]，从而，绘制出均匀成形极限曲线，该曲线是椭圆的一部分，因而，均匀成形极限曲线是外凸的，同样，若从单轴拉伸试验中测定等效断裂应变
Figure BDA0002561774170000042
则可绘制出断裂成形极限曲线。
发明内容
1、本发明解决的技术问题
针对现有技术中的不足，本发明所要解决的技术问题是：在单轴拉伸试验中，测定材料的等效均匀应变
Figure BDA0002561774170000043
和等效断裂应变
Figure BDA0002561774170000044
基于材料的等向硬化模型，材料的均匀成形极限和断裂成形极限下的等效塑性应变与临界网格圆所处的应力应变状态无关，并根据单轴拉伸下的
Figure BDA0002561774170000045
和
Figure BDA0002561774170000046
即可从理论上计算出处于同一变形程度，如均匀成形极限或断裂成形极限下的任意简单应变状态下的真实主应变εmajor和真实次应变εminor，从而，无需通过现有的FLD测定方法，即可绘制出等效塑性应变成形极限图(Equivalent Plastic StrainForming Limit Diagram，简称“EPS-FLD”)，直接避免了现有技术中的不足。
2、本发明的技术方案
为了实现本发明所要解决的技术问题的目的，本发明提供了一种测定等效塑性应变成形极限图的试验与计算方法，其步骤包括：
步骤一、制作被测试材料的单轴拉伸试验试样，并进行单轴拉伸试验，测定被测试材料的真实均匀应变
Figure BDA0002561774170000047
和真实断裂应变
Figure BDA0002561774170000048
所述的真实均匀应变
Figure BDA0002561774170000049
是指：在均匀颈缩变形阶段，试样上初始长宽高分别为l0·w0·h0的有限体A始终保持均匀变形直到判断试样均匀颈缩结束，其中，l0为初始标距，w0为试样初始宽度，h0为试样厚度，有限体A保持极限均匀变形的真实均匀应变
Figure BDA00025617741700000410
采用下式计算：
Figure BDA00025617741700000411
其中，lu是指初始标距l0在均匀颈缩结束时的变形长度；
所述的真实断裂应变
Figure BDA00025617741700000412
是指：在局部颈缩变形阶段，在试样断裂中心区域，始终存在初始长宽高分别为l0'·l0'·h0的有限体始终保持均匀变形直到判断材料失效，取满足该条件的最大初始标距为lmax的有限体B保持极限均匀变形的真实断裂应变
Figure BDA0002561774170000051
采用下式计算：
Figure BDA0002561774170000052
其中，lf是指有限体B的最大初始标距lmax在断裂时刻的变形长度；
所述的真实均匀应变
Figure BDA0002561774170000053
亦等于等效均匀应变
Figure BDA0002561774170000054
所述的真实断裂应变
Figure BDA0002561774170000055
亦等于等效断裂应变
Figure BDA0002561774170000056
步骤二、按以下两式分别计算不同简单应变路径下的真实主应变εmajor和真实次应变εminor：
Figure BDA0002561774170000057
Figure BDA0002561774170000058
其中，β是真实次应变εminor与真实主应变εmajor之比，其某一特定的取值表示某一特定的简单应变路径状态，取值范围为[-∞,+∞]，
Figure BDA0002561774170000059
为等效塑性应变，其取值与应变状态无关；
或者，按以下两式分别计算不同简单应变路径下的真实主应变εmajor和真实次应变εminor：
Figure BDA00025617741700000510
Figure BDA00025617741700000511
其中，r是塑性应变比，其某一特定的取值表示某一特定的简单应变路径状态，取值范围为[-∞,+∞]，
Figure BDA00025617741700000512
为等效塑性应变，其取值与应变状态无关。
步骤三、以真实主应变εmajor为纵坐标，真实次应变εminor为横坐标，绘制等效塑性应变成形极限曲线，若取等效塑性应变
Figure BDA00025617741700000513
等于真实均匀应变
Figure BDA00025617741700000514
绘制的曲线称之为均匀成形极限曲线；若取等效塑性应变
Figure BDA00025617741700000515
等于真实断裂应变
Figure BDA00025617741700000516
绘制的曲线称之为断裂成形极限曲线；两条曲线的形状均为椭圆，将两条曲线绘制于同一张图片中，则得到被测材料的等效塑性应变成形极限图。
进一步的，所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：在步骤一中，所述的真实均匀应变
Figure BDA00025617741700000517
是指初始标距l0为10mm，试样上的有限体A始终处于均匀颈缩变形下的最大真实应变；所述的真实断裂应变
Figure BDA00025617741700000518
是指有限体B的最大初始标距为lmax为1.0mm，则所取的长宽高为1.0mm·1.0mm·h0的有限体B始终处于均匀颈缩变形下的最大真实应变。
进一步的，所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：在步骤三中，取成形极限图所覆盖的简单应变路径状态，其范围是从纯双轴拉伸状态逐步变化到纯单轴拉伸状态，若在步骤二中采用(3)和(4)式计算，则β对应的取值范围是[1，-0.5,]，若在步骤二中采用(5)和(6)式计算，则r对应的取值范围是[-0.5,1]，得到与传统成形极限图具有相同简单应变路径范围的被测试材料的等效塑性应变成形极限图。
进一步的，所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：测定材料1.0mm标距下的塑性应变比r1.0，将r值的取值范围由[-0.5,1]更新为[-0.5,r1.0]，则对应的等效塑性应变成形极限图随之更新。
进一步的，所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：测定材料在单轴拉伸状态下的边缘损伤系数η，在局部等效塑性应变成形极限图上，将两条成形极限曲线上真实次应变εminor的横坐标值小于
Figure BDA0002561774170000061
的数据删除，得到进一步修正后的等效塑性应变成形极限图，所述的边缘损伤系数η按下式计算：
Figure BDA0002561774170000062
其中，真实断裂应变
Figure BDA0002561774170000063
由步骤一中的单轴拉伸试验测定，扩孔率λ由材料扩孔试验测定。
进一步的，所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：在步骤一中，在不同的应变速率下进行单轴拉伸试验，则得到各应变速率下的等效塑性应变成形极限图。
3、本发明的有益效果
本发明提供了一种测定等效塑性应变成形极限图的试验与计算方法，避免了现有技术的不足，相比现有的技术方案，其产生的有益效果主要体现在以下两个方面：
第一，由于现有的FLD，在理论上存在缺陷——FLC是内凹的，实际上是外凸的，在方法上不严谨——主应变和次应变测定没有经过体积不变条件验证，正是由于传统FLD在理论和方法上的不足，导致在实际工程应用中存在各种问题，已经不能满足新形势下汽车用钢的轻量化技术需求，因此，本发明的EPS-FLD通过理论创新避免了现有技术的各种不足，为新形势下的汽车用钢提供了技术解决方案。
第二，由于现有的FLD，需要制作不同简单应变路径下的测量试样，并采用专门的试验设备进行测量，因此，测定材料的FLD费时费力，本发明将材料的均匀成形极限和断裂成形极限的测定集成在单轴拉伸试验中，无需额外增加试验，且可以实现更准确的技术目标，节省了宝贵的人力成本和设备成本。
综上所述，通过本发明的实施，可以以极低的成本获取重要的材料成形参数，因此，本发明具有重要的理论与工程实践意义及非常广阔的应用前景。
附图说明
下面结合附图对本发明作进一步详细说明。
图1是本实施例流程示意图；
图2是材料DP780的均匀成形极限曲线；
图3是材料DP780的完整均匀成形极限曲线；
图4是材料DP780的断裂成形极限曲线；
图5是材料DP780的完整断裂成形极限曲线；
图6是材料DP780的完整等效塑性应变成形极限图；
图7是材料DP780的等效塑性应变成形极限图；
图8是材料DP780的均匀变形和变形的简单应变路径区域；
图9是考虑塑性应变比影响，修正后的等效塑性应变成形极限图；
图10是考虑边缘损伤系数影响，修正后的等效塑性应变成形极限图；
图11是考虑塑性应变比和边缘损伤系数影响，修正后的等效塑性应变成形极限图。
具体实施方式
下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述。
实施例
本实施例以材料DP780/1.4mm，采用的拉伸试验机型号为Zwick/Roell Z050、DIC测试***为德国GOM公司的ARAMIS***，详细说明本发明的测定材料等效塑性应变成形极限图的试验与计算方法，如图1所示，其包括如下步骤：
步骤一、按照“一种测定材料真实应力应变曲线的试验与计算方法，申请号201910801077.1”中所描述的方法，测定被测试材料的真实均匀应变(True UniformStrain)
Figure BDA0002561774170000071
和真实断裂应变(True Fracture Strain)
Figure BDA0002561774170000072
单轴拉伸试验在准静态拉伸速率下进行。
所述的真实均匀应变
Figure BDA0002561774170000073
是指：在单轴拉伸应变状态下的均匀颈缩变形阶段，初始长宽高分别为10mm·12mm·1.4mm的有限体A始终保持均匀变形直到判断试样均匀颈缩结束，以初始标距l0＝10mm计算的最大真实应变。由DIC技术测量的初始标距l0在均匀颈缩结束时的变形长度lu＝12.254mm，代入下式计算真实均匀应变
Figure BDA0002561774170000074
Figure BDA0002561774170000081
所述的真实断裂应变
Figure BDA0002561774170000082
是指：在单轴拉伸应变状态下的局部颈缩变形阶段，在试样断裂中心区域，有限体B始终保持均匀变形直到判断材料失效，满足该条件的有限体B的最大初始标距lmax＝1.0mm时计算的最大真实应变，有限体B的初始长宽高分别为1.0mm·1.0mm·1.4mm。由DIC技术测量的最大初始标距lmax在断裂时刻的变形长度lf＝1.537mm，采用下式计算：
Figure BDA0002561774170000083
所述的真实均匀应变
Figure BDA0002561774170000084
亦等于等效均匀应变
Figure BDA0002561774170000085
所述的真实断裂应变
Figure BDA0002561774170000086
亦等于等效断裂应变
Figure BDA0002561774170000087
步骤二、按以下两式分别计算不同简单应变路径下的真实主应变εmajor和真实次应变εminor：
Figure BDA0002561774170000088
Figure BDA0002561774170000089
其中，r是塑性应变比，取值范围为[-∞,+∞]，其某一特定的取值表示某一特定的简单应变路径状态，
Figure BDA00025617741700000810
为等效塑性应变，其值与应变状态无关。
若等效塑性应变
Figure BDA00025617741700000811
分别取真实均匀应变
Figure BDA00025617741700000812
和真实断裂应变
Figure BDA00025617741700000813
r取值范围为[-∞,+∞]，按照(5)式和(6)式计算的真实主应变(Major TrueStrain)εmajor和真实次应变(Minor True Strain)εminor如表1所示。
步骤三、以真实主应变εmajor为纵坐标，真实次应变εminor为横坐标，根据表1的计算结果绘制等效塑性应变成形极限图(EPS-FLD)。
若取等效塑性应变
Figure BDA00025617741700000814
等于真实均匀应变
Figure BDA00025617741700000815
绘制的曲线称之为均匀成形极限曲线，如图2所示；将表1中均匀变形极限下的真实主应变εmajor和真实次应变εminor均乘以-1，保持两者的对应关系不变，得到新的计算数据，根据表1和新的计算数据，绘制出完整的均匀成形极限曲线，如图3所示，其形状为椭圆。
若取等效塑性应变
Figure BDA00025617741700000816
等于真实断裂应变
Figure BDA00025617741700000817
绘制的曲线称之为断裂成形极限曲线，如图4所示；将表1中断裂变形极限下的真实主应变εmajor和真实次应变εminor均乘以-1，保持两者的对应关系不变，得到新的计算数据，根据表1和新的计算数据，绘制出完整的断裂成形极限曲线，如图5所示，其形状为椭圆。
表1不同简单应变状态下的真实主应变和真实次应变
Figure BDA0002561774170000091
将图3和图5中形状为标准椭圆的两条曲线的等效塑性应变成形极限曲线，绘制于同一张图中，则得到被测材料的等效塑性应变成形极限图，如图6所示，其本质是：基于等向硬化模型，由三维主应变空间(球形空间)初始屈服面及后继屈服面退化或降维的二维主应变平面；r值的实质是对简单应变路径的表征参数。
在表1中，仅取r的取值范围[-0.5,1]对应的真实主应变εmajor和真实次应变εminor数据，在同一张图中绘制均匀成形极限曲线和断裂成形极限曲线，得到工程上应用的等效塑性应变成形极限图，结果如图7所示，其物理意义是简单应变路径状态的范围是从纯双轴拉伸状态逐步变化到纯单轴拉伸状态，与传统的成形极限图表示的简单应变路径范围一致，可用于冲压成形仿真材料DP780的成形极限图的参数输入。
在图7中，简单应变路径r＝1、均匀成形极限曲线(Uniform_EPS-FLC)和简单应变路径r＝-0.5所围成的扇形区域，如图8中的Zone-1区域，表示应变状态处于此区域的材料是安全的，同时，表示零件处于全局或整体的均匀变形状态。若均匀成形极限曲线(Uniform_EPS-FLC)向外移动表示成形广度越大，Zone-1越大表示相同的变形可以分布或扩散到更大的区域，且应变状态分布范围也更大，因此，可成形更复杂的零件，即均匀成形能力越强。
在图7中，简单应变路径r＝1、断裂成形极限曲线(Fracture_EPS-FLC)、简单应变路径r＝-0.5和均匀成形极限曲线(Uniform_EPS-FLC)所围成的扇形区域，如图8中的Zone-2区域，表示应变状态处于此区域的材料存在不同程度的风险，同时，表示零件处于局部变形状态——增加的变形由零件上的局部区域承担，增加的变形越大，开裂风险越大。若断裂成形极限曲线(Fracture_EPS-FLC)向外移动表示成形深度越大，Zone-2越大表示变形只发生在局部区域的应变状态范围，其区域越大表示局部成形能力越好。
按照“一种测定金属材料塑性应变比的试验与计算方法，申请号201910801202.9”中所描述的方法，测定材料DP780在1.0mm标距下的塑性应变比r1.0＝0.78，将r值的取值范围为由[-0.5,1]更新为[-0.5,0.78]，则对应的如图7所示的等效塑性应变成形极限图随之更新，结果如图9所示，图中阴影部分的数据将从成形极限图中删除(如图9中的Zone-3区域)。
材料DP780的真实断裂应变
Figure BDA0002561774170000101
由扩孔试验测量的扩孔率λ＝23，材料DP780在单轴拉伸状态下的边缘损伤系数η按下式计算：
Figure BDA0002561774170000102
则按
Figure BDA0002561774170000103
计算的带损伤的真实次应变εminor＝-0.104。
在如图7所示的等效塑性应变成形极限图上，将两条成形极限曲线上横坐标值小于真实次应变-0.104的数据删除，得到修正后的等效塑性应变成形极限图，结果如图10所示，图中阴影部分(如图10中的Zone-4区域)的数据将从成形极限图中删除。
在如图9所示的等效塑性应变成形极限图上，将两条成形极限曲线上横坐标值小于真实次应变-0.104的数据删除，得到修正后的等效塑性应变成形极限图，结果如图11所示，图中阴影部分(如图11中的Zone-5区域)的数据将从成形极限图中删除，修正后的成形极限图仅用于零件上存在边缘损伤的自由边区域，非自由边区域依然适应于未修正的成形极限图。
以上实施方式仅为本发明的示例性实施方式，不用于限制本发明，本发明的保护范围由权利要求书限定。本领域技术人员在本发明的实质性保护范围内，对本发明做出的各种修改或等同替换也落在本发明的保护范围内。

Claims (6)

1.一种测定等效塑性应变成形极限图的试验与计算方法，其步骤包括：
步骤一、制作被测试材料的单轴拉伸试验试样，并进行单轴拉伸试验，测定被测试材料的真实均匀应变
Figure FDA0002561774160000011
和真实断裂应变
Figure FDA0002561774160000012
所述的真实均匀应变
Figure FDA0002561774160000013
是指：在均匀颈缩变形阶段，试样上初始长宽高分别为l0·w0·h0的有限体A始终保持均匀变形直到判断试样均匀颈缩结束，其中，l0为初始标距，w0为试样初始宽度，h0为试样厚度，有限体A保持极限均匀变形的真实均匀应变
Figure FDA0002561774160000014
采用下式计算：
Figure FDA0002561774160000015
其中，lu是指初始标距l0在均匀颈缩结束时的变形长度；
所述的真实断裂应变
Figure FDA0002561774160000016
是指：在局部颈缩变形阶段，在试样断裂中心区域，始终存在初始长宽高分别为l0'·l0'·h0的有限体始终保持均匀变形直到判断材料失效，取满足该条件的最大初始标距为lmax的有限体B保持极限均匀变形的真实断裂应变
Figure FDA0002561774160000017
采用下式计算：
Figure FDA0002561774160000018
其中，lf是指有限体B的最大初始标距lmax在断裂时刻的变形长度；
所述的真实均匀应变
Figure FDA0002561774160000019
亦等于等效均匀应变
Figure FDA00025617741600000110
所述的真实断裂应变
Figure FDA00025617741600000111
亦等于等效断裂应变
Figure FDA00025617741600000112
步骤二、按以下两式分别计算不同简单应变路径下的真实主应变εmajor和真实次应变εminor：
Figure FDA00025617741600000113
Figure FDA00025617741600000114
其中，β是真实次应变εminor与真实主应变εmajor之比，其某一特定的取值表示某一特定的简单应变路径状态，取值范围为[-∞,+∞]，
Figure FDA00025617741600000115
为等效塑性应变，其取值与应变状态无关；
或者，按以下两式分别计算不同简单应变路径下的真实主应变εmajor和真实次应变εminor：
Figure FDA00025617741600000116
Figure FDA0002561774160000021
其中，r是塑性应变比，其某一特定的取值表示某一特定的简单应变路径状态，取值范围为[-∞,+∞]，
Figure FDA0002561774160000022
为等效塑性应变，其取值与应变状态无关。
步骤三、以真实主应变εmajor为纵坐标，真实次应变εminor为横坐标，绘制等效塑性应变成形极限曲线，若取等效塑性应变
Figure FDA0002561774160000023
等于真实均匀应变
Figure FDA0002561774160000024
绘制的曲线称之为均匀成形极限曲线；若取等效塑性应变
Figure FDA0002561774160000025
等于真实断裂应变
Figure FDA0002561774160000026
绘制的曲线称之为断裂成形极限曲线；两条曲线的形状均为椭圆，将两条曲线绘制于同一张图片中，则得到被测材料的等效塑性应变成形极限图。
2.根据权利要求1所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：在步骤一中，所述的真实均匀应变
Figure FDA0002561774160000027
是指初始标距l0为10mm，试样上的有限体A始终处于均匀颈缩变形下的最大真实应变；所述的真实断裂应变
Figure FDA0002561774160000028
是指有限体B的最大初始标距为lmax为1.0mm，则所取的长宽高为1.0mm·1.0mm·h0的有限体B始终处于均匀颈缩变形下的最大真实应变。
3.根据权利要求2所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：在步骤三中，取成形极限图所覆盖的简单应变路径状态，其范围是从纯双轴拉伸状态逐步变化到纯单轴拉伸状态，若在步骤二中采用(3)和(4)式计算，则β对应的取值范围是[1，-0.5,]，若在步骤二中采用(5)和(6)式计算，则r对应的取值范围是[-0.5,1]，得到与传统成形极限图具有相同简单应变路径范围的被测试材料的等效塑性应变成形极限图。
4.根据权利要求3所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：测定材料1.0mm标距下的塑性应变比r1.0，将r值的取值范围由[-0.5,1]更新为[-0.5,r1.0]，则对应的等效塑性应变成形极限图随之更新。
5.根据权利要求4所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：测定材料在单轴拉伸状态下的边缘损伤系数η，在局部等效塑性应变成形极限图上，将两条成形极限曲线上真实次应变εminor的横坐标值小于
Figure FDA0002561774160000029
的数据删除，得到进一步修正后的等效塑性应变成形极限图，所述的边缘损伤系数η按下式计算：
Figure FDA00025617741600000210
其中，真实断裂应变
Figure FDA00025617741600000211
由步骤一中的单轴拉伸试验测定，扩孔率λ由材料扩孔试验测定。
6.根据权利要求1或2或5所述的测定等效塑性应变成形极限图的试验与计算方法，其特征在于：在步骤一中，在不同的应变速率下进行单轴拉伸试验，则得到各应变速率下的等效塑性应变成形极限图。
CN202010618096.3A 2020-06-30 2020-06-30 一种测定等效塑性应变成形极限图的试验与计算方法 Active CN111896373B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010618096.3A CN111896373B (zh) 2020-06-30 2020-06-30 一种测定等效塑性应变成形极限图的试验与计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010618096.3A CN111896373B (zh) 2020-06-30 2020-06-30 一种测定等效塑性应变成形极限图的试验与计算方法

Publications (2)

Publication Number Publication Date
CN111896373A true CN111896373A (zh) 2020-11-06
CN111896373B CN111896373B (zh) 2023-03-24

Family

ID=73191134

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010618096.3A Active CN111896373B (zh) 2020-06-30 2020-06-30 一种测定等效塑性应变成形极限图的试验与计算方法

Country Status (1)

Country Link
CN (1) CN111896373B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112964555A (zh) * 2021-02-10 2021-06-15 肖锋 一种测定等效塑性应变成形极限图的试验与计算方法
CN113010932A (zh) * 2021-02-10 2021-06-22 肖锋 一种成形极限应变云图的后处理方法
CN113176738A (zh) * 2021-03-23 2021-07-27 首钢集团有限公司 控制冲压成型性的方法及装置
CN113624602A (zh) * 2021-07-29 2021-11-09 中国科学院金属研究所 管材成形极限图右侧区域曲线的实验装置及构建方法
CN114322915A (zh) * 2021-12-06 2022-04-12 武汉钢铁有限公司 一种材料成形极限的测定方法及装置
CN114509341A (zh) * 2022-02-23 2022-05-17 北京理工大学重庆创新中心 一种测量材料试样拉伸断裂试验过程中应力三轴度的方法
CN114544383A (zh) * 2021-12-23 2022-05-27 武汉上善仿真科技有限责任公司 一种测定极限尖冷弯角和等效断裂应变的试验方法
CN115326504A (zh) * 2022-07-26 2022-11-11 燕山大学 用于薄钢板边缘拉伸开裂极限的评价方法及其评价装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101379381A (zh) * 2006-02-01 2009-03-04 新日本制铁株式会社 断裂预测方法
CN101598720A (zh) * 2009-06-26 2009-12-09 东北大学 一种建立相变诱发塑性钢板成形极限图的实验方法
CN101599094A (zh) * 2009-06-26 2009-12-09 东北大学 一种建立相变诱发塑性钢板成形极限图的模型方法
CN109870362A (zh) * 2019-03-04 2019-06-11 燕山大学 一种高强铝合金板材的断裂成形极限图建立方法及***
CN110411837A (zh) * 2019-08-28 2019-11-05 武汉上善仿真科技有限责任公司 一种测定材料标距效应曲线的试验与计算方法
CN110411836A (zh) * 2019-08-28 2019-11-05 武汉上善仿真科技有限责任公司 一种测定金属材料塑性应变比的试验与计算方法
CN110501224A (zh) * 2019-08-28 2019-11-26 肖锋 一种测定材料真实应力应变曲线的试验与计算方法
CN111292813A (zh) * 2020-01-17 2020-06-16 中国石油天然气集团有限公司 一种基于断裂准则绘制各向异性材料断裂成形极限图的方法及其使用方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101379381A (zh) * 2006-02-01 2009-03-04 新日本制铁株式会社 断裂预测方法
US20090177417A1 (en) * 2006-02-01 2009-07-09 Shigeru Yonemura Fracture prediction method, device, a program arrangement and computer-accessible medium therefor
CN101598720A (zh) * 2009-06-26 2009-12-09 东北大学 一种建立相变诱发塑性钢板成形极限图的实验方法
CN101599094A (zh) * 2009-06-26 2009-12-09 东北大学 一种建立相变诱发塑性钢板成形极限图的模型方法
CN109870362A (zh) * 2019-03-04 2019-06-11 燕山大学 一种高强铝合金板材的断裂成形极限图建立方法及***
CN110411837A (zh) * 2019-08-28 2019-11-05 武汉上善仿真科技有限责任公司 一种测定材料标距效应曲线的试验与计算方法
CN110411836A (zh) * 2019-08-28 2019-11-05 武汉上善仿真科技有限责任公司 一种测定金属材料塑性应变比的试验与计算方法
CN110501224A (zh) * 2019-08-28 2019-11-26 肖锋 一种测定材料真实应力应变曲线的试验与计算方法
CN111292813A (zh) * 2020-01-17 2020-06-16 中国石油天然气集团有限公司 一种基于断裂准则绘制各向异性材料断裂成形极限图的方法及其使用方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
曹宏深: "复合变形路径下的板材冲压成形极限应变(Ⅱ)――冲压成形极限应变的计算与试验", 《机械工程学报》 *
魏毅静等: "金属板材成形性能的力学分析", 《宝钢技术》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113010932A (zh) * 2021-02-10 2021-06-22 肖锋 一种成形极限应变云图的后处理方法
CN112964555A (zh) * 2021-02-10 2021-06-15 肖锋 一种测定等效塑性应变成形极限图的试验与计算方法
CN113176738A (zh) * 2021-03-23 2021-07-27 首钢集团有限公司 控制冲压成型性的方法及装置
CN113176738B (zh) * 2021-03-23 2022-06-21 首钢集团有限公司 控制冲压成型性的方法及装置
CN113624602B (zh) * 2021-07-29 2022-07-15 中国科学院金属研究所 管材成形极限图右侧区域曲线的实验装置及构建方法
CN113624602A (zh) * 2021-07-29 2021-11-09 中国科学院金属研究所 管材成形极限图右侧区域曲线的实验装置及构建方法
CN114322915A (zh) * 2021-12-06 2022-04-12 武汉钢铁有限公司 一种材料成形极限的测定方法及装置
CN114322915B (zh) * 2021-12-06 2023-10-03 武汉钢铁有限公司 一种材料成形极限的测定方法及装置
CN114544383A (zh) * 2021-12-23 2022-05-27 武汉上善仿真科技有限责任公司 一种测定极限尖冷弯角和等效断裂应变的试验方法
CN114544383B (zh) * 2021-12-23 2024-03-15 武汉上善仿真科技有限责任公司 一种测定极限尖冷弯角和等效断裂应变的试验方法
CN114509341B (zh) * 2022-02-23 2023-08-25 北京理工大学重庆创新中心 一种测量材料试样拉伸断裂试验过程中应力三轴度的方法
CN114509341A (zh) * 2022-02-23 2022-05-17 北京理工大学重庆创新中心 一种测量材料试样拉伸断裂试验过程中应力三轴度的方法
CN115326504A (zh) * 2022-07-26 2022-11-11 燕山大学 用于薄钢板边缘拉伸开裂极限的评价方法及其评价装置
CN115326504B (zh) * 2022-07-26 2024-06-04 燕山大学 用于薄钢板边缘拉伸开裂极限的评价方法及其评价装置

Also Published As

Publication number Publication date
CN111896373B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
CN111896373B (zh) 一种测定等效塑性应变成形极限图的试验与计算方法
Hussaini et al. Development of experimental and theoretical forming limit diagrams for warm forming of austenitic stainless steel 316
Butuc et al. An experimental and theoretical analysis on the application of stress-based forming limit criterion
CN109933925B (zh) 一种金属板材的冲压成形性能预测方法
Choi et al. Effect of punch speed on amount of springback in U-bending process of auto-body steel sheets
CN111950098A (zh) 基于有限元仿真的金属板料冲孔工艺参数优化方法
CN113420391A (zh) 一种获得复杂应力状态下材料高精度硬化模型参数的方法
CN113764056B (zh) 一种获得材料多应变率下高精度硬化模型参数的方法
CN112257302B (zh) 一种采用极限厚度构建板料成形极限曲线的方法
Lemoine et al. Flow curve determination at large plastic strain levels: limitations of the membrane theory in the analysis of the hydraulic bulge test
CN106769439A (zh) 一种管线钢热轧卷板屈服强度的测试方法
Kuwabara et al. Hole expansion simulation considering the differential hardening of a sheet metal
CN112964555A (zh) 一种测定等效塑性应变成形极限图的试验与计算方法
JP4696893B2 (ja) 鋼管の局部座屈性能評価方法、鋼管の材質設計方法、鋼管の製造方法
CN114386216B (zh) 一种基于剪切应变路径构建高强钢加工硬化的方法
Tang et al. Analysis of shear stress wrinkling of asymmetric sheet specimen under offset loading
KR100402923B1 (ko) 금속재료,특히판금의2축변형거동을결정하기위한측정방법
CN112345379A (zh) 一种板材弯曲回弹特性的测试装置及其测试方法
CN111272551A (zh) 一种高速拉伸试验数据曲线处理方法
CN111475896B (zh) 一种棒材矫直应力中性层位置确定方法
JP2015047605A (ja) 曲げ内側割れの評価方法および予測方法
Raupach et al. Mesh refinement study and experimental validation for stretch bending of sheet metals
JP7464176B1 (ja) プレス成形割れ判定方法、装置及びプログラム、並びにプレス成形品の製造方法
CN102269677A (zh) 塑性应变比r值异常值的处理方法
CN113010932A (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