CN101599094B - 一种建立相变诱发塑性钢板成形极限图的模型方法 - Google Patents

一种建立相变诱发塑性钢板成形极限图的模型方法 Download PDF

Info

Publication number
CN101599094B
CN101599094B CN2009100122578A CN200910012257A CN101599094B CN 101599094 B CN101599094 B CN 101599094B CN 2009100122578 A CN2009100122578 A CN 2009100122578A CN 200910012257 A CN200910012257 A CN 200910012257A CN 101599094 B CN101599094 B CN 101599094B
Authority
CN
China
Prior art keywords
strain
steel plate
induced plasticity
transformation
plasticity steel
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.)
Expired - Fee Related
Application number
CN2009100122578A
Other languages
English (en)
Other versions
CN101599094A (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.)
Northeastern University China
Original Assignee
Northeastern University China
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 Northeastern University China filed Critical Northeastern University China
Priority to CN2009100122578A priority Critical patent/CN101599094B/zh
Publication of CN101599094A publication Critical patent/CN101599094A/zh
Application granted granted Critical
Publication of CN101599094B publication Critical patent/CN101599094B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

一种建立相变诱发塑性钢板成形极限图的模型方法,涉及成形极限图的建立方法,按以下步骤进行:(1)测定TRIP钢板的厚度和应变硬化指数;(2)计算平面应变状态下相变诱发塑性钢板的极限应变值;(3)建立应变坐标系,将数据合成在坐标系上,建立该相变诱发塑性钢板的成形极限图。本发明是解决TRIP钢板料冲压问题的一个有效的工具,能成功预测TRIP钢板料失效,有效评价其成形性能,减少了现有建立成形极限图的测试工作量。

Description

一种建立相变诱发塑性钢板成形极限图的模型方法
技术领域
本发明涉及钢板成形极限图的建立方法,特别涉及一种建立相变诱发塑性钢板成形极限图的模型方法。
背景技术
TRIP(Transformation Induced Plasticity)钢又称相变诱发塑性钢,是近几年发展起来的一种高强度高塑性钢。用作汽车钢板可以减轻车重,降低油耗,同时有较强的能量吸收能力,能够抵御撞击时的塑性变形,显著提高汽车的安全等级,具有明显的优越性。板料成形是金属成形的一个重要组成部分,它在航空航天、汽车、轻工及家电等制造业中有着非常广泛的应用。在板料成形过程中,人们迫切需要知道板料能够承受的成形极限及在成形过程中发生破坏的时间和位置,以便采取措施加以避免。板料成形中的基本实验只能对板料的成形性能做出定性综合的一般评价,模拟实验又只能对少数典型工序,在较为单纯的条件下进行,故所得结果很难对复杂零件的成形性能做出确切的判断,也不能很好地处理生产中遇到的具体问题。而成形极限图(forming limit diagram,FLD)的出现为人们研究板料成形极限以及评价板料成形性能提供了基础,是解决板料冲压问题的一个有效的工具,能成功预测板料失效,有效评价其成形性能。成形极限图的获得通常采用标准的试验装置,通过改变试件的宽度和润滑条件,基于网格应变分析技术来获得极限应变数据。从关于成形极限图的测试方法可知,其测试方法虽然在不断的完善,但是总的来说还存在费工、费料和费时等缺点。
发明内容
本发明的目的是针对现有建立板带材成形极限图需要标准试验装置、制备不同宽度尺寸试件和改变润滑条件等费工、费料和费时等缺点,本发明提供一种只需测定钢板的厚度和应变硬化指数就可以建立相变诱发塑性钢板成形极限图的模型方法,减少建立相变诱发塑性钢板成形极限图的测试工作量,成功预测相变诱发塑性钢板料失效,有效评价其成形性能。
实现本发明目的的技术解决方案如下:
1、测定相变诱发塑性钢板的厚度和应变硬化指数。
2、计算该相变诱发塑性钢板平面应变状态下的极限应变值,计算公式为:
FLD 0 = ( 23.3 + 360 t 25.4 ) × n 21 - - - ( 1 )
其中t为相变诱发塑性钢板的厚度,单位为mm;n为相变诱发塑性钢板的应变硬化指数;FLD0为相变诱发塑性钢板平面应变状态下的极限应变值。
工程主应变和工程次应变的测试方法按国家标准GB/T15825.8-1995《金属薄板成形性能与试验方法-成形极限图(FLD)试验》进行,工程主应变为e1,工程次应变为e2,当工程次应变e2=0时,成形极限图中平面应变区域的工程主应变e1便是FLD0;测得FLD0与钢板厚度以及钢板硬化指数的关系如(1)式所示。
3、建立相变诱发塑性钢板成形极限图建立的工程主应变e1与工程次应变e2的关系式为:
以工程主应变e1为纵坐标,工程次应变e2为横坐标,建立应变坐标系,将根据公式(1)和(2)计算获得的数据合成在坐标系上,建立该相变诱发塑性钢板的成形极限图。
成形极限图中工程次应变e2<0的部分为左半区域,工程次应变e2>0的部分为右半区域。
成形极限图的左半区域中,工程主应变(e1)与工程次应变(e2=0)的关系等式如下:
e2<0时,e1=-e2+FLD0                (3)
即由平面应变点FLD0在左侧引一条与横轴负方向夹角为45°的直线,构成成形极限图的左板区域。
成形极限图的右半区域中,在平面变形条件下,对于胀形区,在一定的应变比下,宽度方向的真应变与厚度方向的真实应变存在的相关用Cr表示为:
Cr = ϵ w ϵ t - - - ( 4 )
εw为宽度方向的真应变,εt为厚度方向的应变;根据体积不变原理,长度方向的真应变εl与εw和εt之间,存在以下关系:
ϵ l + ϵ w + ϵ t = 0 ⇔ ϵ l = - ϵ w ( 1 + C r ) / C r - - - ( 5 )
根据工程应变和真实应变的关系ε=ln(1+e),可以得到以下关系式:
ln ( 1 + e 1 ) = ln ( 1 + e 2 ) - ( 1 + C r ) / C r - - - ( 6 )
上式中设 ( 1 + e 2 ) - ( 1 + C r ) / C r = C 1 ( 1 + e 2 ) C , 由上式推导可得工程主应变e1与工程次应变e2的关系式为:
e1=C1×(1+e2)C-1                    (7)
其中C为常数。从成形极限图可知,当e2=0时,e1=FLD0,从而可得常数C1=1+FLD0。于是由公式(7)可得:
e1=(1+FLD0)×(1+e2)C-1              (8)
选用典型TRIP钢的FLD0数据,通过最小二乘法原理进行拟合分析得到C≈0.5。
从而成形极限图右半区域的工程主应变e1与工程次应变e2的关系式为:
e2>0时,e1=(1+FLD0)×(1+e2)0.5-1               (9)
工程主应变和工程次应变的测试方法按国家标准GB/T15825.8-1995《金属薄板成形性能与试验方法-成形极限图(FLD)试验》进行,工程主应变为e1,工程次应变为e2,当工程次应变e2=0时,成形极限图中平面应变区域的工程主应变e1便是FLD0;测得FLD0与钢板厚度以及钢板硬化指数的关系如(1)式所示。
上述的相变诱发塑性钢板应变硬化指数的测试方法为:通过单轴拉伸实验测得其工程应力应变曲线,在5%至35%的均匀塑性变形范围内将其转变为真应力应变曲线,真应力与真应变关系式为:
σ=k·εn                 (10)
上式中,σ为真应力、ε为真应变、n为相变诱发塑性钢板的应变硬化指数、k为系数,将式(10)转变成对数方程:
lnσ=lnk+nlnε            (11)
在双对数坐标平面上的直线的斜率即为相变诱发塑性钢板的应变硬化指数n。
本发明通过模型方法建立相变诱发塑性钢板的成形极限图,该方法是解决TRIP钢板料冲压问题的一个有效的工具,能成功预测相变诱发塑性钢板料失效,有效评价其成形性能。模型方法只需测定钢板的厚度和应变硬化指数就可以建立相变诱发塑性钢板成形极限图,减少了现有建立成形极限图的测试工作量。
附图说明
图1为本发明实施例1中的成形极限图的模型计算结果和实验结果对比图;
图2为本发明实施例2中的成形极限图的模型计算结果和实验结果对比图;
图3为本发明实施例3中的成形极限图的模型计算结果和实验结果对比图;
图4为本发明实施例4中的成形极限图的模型计算结果和实验结果对比图;
图中a为实验数据,b为建立的模型。
具体实施方式
本发明实施例中TRIP钢板的制备方法分为两种:1、将低碳低硅钢冶炼成钢锭,成分按重量百分比为含C0.18~0.20%,Mn1.30~1.45%,Si0.60~0.70%,Al0.50~0.60%,P≤0.009%,S≤0.007%,Nb≤0.040%,余量为Fe。钢锭经1200℃保温1.5~2hr后,锻造成30mm厚的锻坯方料。热轧至3.5mm厚,再冷轧至1.0mm厚,获得冷轧钢板,将冷轧钢板在800℃的氯化盐中保温5min,然后淬入温度为400℃的硝酸盐中保温5min,空冷至室温,制备成TRIP钢板。2、将低碳高锰钢冶炼成钢锭,其化学成分C600ppm,Mn23.8%,Si3.0%,Al2.7%,余量为Fe。钢锭经1150℃保温1.5~2hr后,锻造成30mm厚的锻坯方料。热轧至5.0mm厚后再冷轧至1.0mm厚。将冷轧钢板加热到1100℃保温15min,然后水冷至室温,制备成TRIP钢板。
本发明实施例中工程主应变和工程次应变的测试方法为:按国家标准GB/T15825.8-1995《金属薄板成形性能与试验方法-成形极限图(FLD)试验》进行,实验设备为BCS-30D型通用板料成形试验机,试样为长度180mm,宽度分别为180mm、160mm、120mm、100mm、80mm、40mm和20mm的矩形试样,每组取3块,在试样上印制,网格形式为直径2.5mm相切圆。
通过成形试验机对试样进行胀形实验,至试样发生破裂,选择靠近裂纹区的网格作为临界网格点,拍摄至少一个二维影像图,采用ASAME自动应变测量系统将二维影像图处理成三维影像图,对临界网格圆进行分析,临界网格圆为试样上紧靠颈缩或裂纹的网格圆。
测量临界网格圆的长轴直径和短轴直径,然后根据以下公式计算工程主应变和工程次应变:
e1=(d1-d0)/d0×100(%)              (12)
e2=(d2-d0)/d0×100(%)              (13)
式中e1为工程主应变,e2为工程次应变,d1为临界网格圆长轴直径,d2为临界网格圆短轴直径,d0为网格圆初始直径。
实施例1
首先制取相变诱发塑性钢板,方法是将低碳硅钢冶炼成钢锭,其成分按重量百分比为含C0.19%,Mn1.32%,Si 0.70%,Al 0.51%,P 0.0068%,S 0.0068%,余量为Fe。钢锭经1200℃保温1.5~2hr后,锻造成30mm厚的锻坯方料,热轧至3.5mm厚,再冷轧至1.0mm厚,获得冷轧钢板,将冷轧钢板切割成合理尺寸的小块,在800℃的氯化盐中保温5min,然后淬入温度为400℃的硝酸盐中保温5min,空冷至室温,制备成TRIP钢板。
测定该钢板的应变硬化指数n为0.323,t=1.0mm,根据公式(1)计算得FLD0=0.576。建立相变诱发塑性钢板成形极限图的工程主应变e1与工程次应变e2的关系式为:
e2<0时,e1=-e2+0.576
e2=0时,e1=FLD0=0.576
e2>0时,e1=(1+0.576)×(1+e2)0.5-1=1.576(1+e2)0.5-1
以工程主应变e1为纵坐标,工程次应变e2为横坐标,建立应变坐标系,将根据公式计算获得的数据合成在坐标系上,建立该相变诱发塑性钢板的的成形极限图。
对制备的TRIP钢板进行成形实验,实验结果与成形极限图的模型计算结果如图1所示。从此图可以看出,按模型关系式(2)计算的成形极限的左半支曲线(e2<0)、平面应变状态下的极限应变的FLD0值(e2=0)和右半支曲线(e2>0)与实验数据点的符合性都很好。从而表明本方法提供的建立成形极限图的模型方法能够准确地建立相变诱发塑性钢板的成形极限图。
从图1还可以看出,该TRIP钢板具有良好的成形性,其平面应变状态下的极限应变的FLD0值分别为0.576。在成形曲线以上的区域是TRIP钢板将发生失效的区域,而在成形曲线以下的区域属于TRIP钢板的安全区域。
实施例2
首先制取相变诱发塑性钢板,方法是将低碳硅钢冶炼成钢锭,其成分按重量百分比为含C0.18%,Mn1.38%,Si 0.67%,Al 0.56%,P 0.0075%,S 0.0048%,Nb0.014%,余量为Fe。钢锭经1200℃保温1.5~2hr后,锻造成30mm厚的锻坯方料,热轧至3.5mm厚,再冷轧至1.0mm厚,获得冷轧钢板,将冷轧钢板切割成合理尺寸的小块,在800℃的氯化盐中保温5min,然后淬入温度为400℃的硝酸盐中保温5min,空冷至室温,制备成TRIP钢板。
测定该钢板的应变硬化指数n为0.325,t=1.0mm,根据公式(1)计算得FLD0=0.580。建立相变诱发塑性钢板成形极限图建立的工程主应变e1与工程次应变e2的关系式为:
e2<0时,e1=-e2+0.580
e2=0时,e1=FLD0=0.580
e2>0时,e1=(1+0.580)×(1+e2)0.5-1=1.580(1+e2)0.5-1
以工程主应变e1为纵坐标,工程次应变e2为横坐标,建立应变坐标系,将根据公式计算获得的数据合成在坐标系上,建立该相变诱发塑性钢板的的成形极限图。
对制备的TRIP钢板进行成形实验,实验结果与成形极限图的模型计算结果如图2所示。
实施例3
首先制取相变诱发塑性钢板,方法是将低碳硅钢冶炼成钢锭,其成分按重量百分比为含C0.20%,Mn1.45%,Si 0.61%,Al 0.58%,P 0.0083%,S 0.0050%,Nb0.036%,余量为Fe。钢锭经1200℃保温1.5~2hr后,锻造成30mm厚的锻坯方料,热轧至3.5mm厚,再冷轧至1.0mm厚,获得冷轧钢板,将冷轧钢板切割成合理尺寸的小块,在800℃的氯化盐中保温5min,然后淬入温度为400℃的硝酸盐中保温5min,空冷至室温,制备成TRIP钢板。
测定该钢板的应变硬化指数n为0.310,t=1.0mm,根据公式(1)计算得FLD0=0.553。建立相变诱发塑性钢板成形极限图建立的工程主应变e1与工程次应变e2的关系式为:
e2<0时,e1=-e2+0.553
e2=0时,e1=FLD0=0.553
e2>0时,e1=(1+0.553)×(1+e2)0.5-1=1.553(1+e2)0.5-1
根据该公式建立该TRIP钢板的成形极限图。
对制备的TRIP钢板进行成形实验,实验结果与成形极限图的模型计算结果如图3所示。
实施例4
首先制取相变诱发塑性钢板,方法是将低碳高锰钢冶炼成钢锭,其化学成分C600ppm,Mn23.8%,Si3.0%,Al2.7%,余量为Fe。钢锭经1150℃保温1.5~2hr后,锻造成30mm厚的锻坯方料。热轧至5.0mm厚后再冷轧至1.0mm厚。将冷轧钢板加热到1100℃保温15min,然后水冷至室温,制备成TRIP钢板。
测定该钢板的应变硬化指数n为0.242,t=1.0mm,根据公式(1)计算得FLD0=0.432。建立相变诱发塑性钢板成形极限图建立的工程主应变e1与工程次应变e2的关系式为:
e2<0时,e1=-e2+0.432
e2=0时,e1=FLD0=0.432
e2>0时,e1=(1+0.432)×(1+e2)0.5-1=1.432(1+e2)0.5-1
以工程主应变e1为纵坐标,工程次应变e2为横坐标,建立应变坐标系,将根据公式计算获得的数据合成在坐标系上,建立该相变诱发塑性钢板的成形极限图。
对制备的TRIP钢板进行成形实验,实验结果与成形极限图的模型计算结果如图4所示。

Claims (1)

1.一种建立相变诱发塑性钢板成形极限图的模型方法,其特征在于按以下步骤进行:(1)测定相变诱发塑性钢板的厚度和应变硬化指数;(2)计算该相变诱发塑性钢板平面应变状态下的极限应变值,计算公式为:
FLD 0 = ( 23.3 + 360 t 25.4 ) × n 21 ;
其中t为相变诱发塑性钢板的厚度,单位为mm;n为相变诱发塑性钢板的应变硬化指数,FLD0为极限应变值;
(3)建立相变诱发塑性钢板成形极限图的工程主应变e1与工程次应变e2的关系式为:
Figure FSB00000316689400012
以工程主应变e1为纵坐标,工程次应变e2为横坐标,建立应变坐标系,将根据上述两个公式计算获得的数据合成在坐标系上,建立该相变诱发塑性钢板的成形极限图;所述的应变硬化指数的测试方法为:通过单轴拉伸实验测取工程应力应变曲线,在5%至35%的均匀塑性变形范围内将其转变为真应力应变曲线,真应力与真应变关系式为:
σ=k·εn
上式中,σ为真应力、ε为真应变、n为相变诱发塑性钢板的应变硬化指数、k为系数,将上式转变成对数方程:
lnσ=lnk+nlnε
在双对数坐标平面上的直线的斜率即为相变诱发塑性钢板的应变硬化指数n;
所述的相变诱发塑性钢板成分按重量百分比为含C0.19%,Mn1.32%,Si 0.70%,Al 0.51%,P 0.0068%,S 0.0068%,余量为Fe;或者为含C0.18%,Mn1.38%,Si 0.67%,Al 0.56%,P 0.0075%,S 0.0048%,Nb0.014%,余量为Fe;或者为含C0.20%,Mn1.45%,Si 0.61%,Al 0.58%,P 0.0083%,S 0.0050%,Nb0.036%,余量为Fe;或者为含C600ppm,Mn23.8%,Si3.0%,Al2.7%,余量为Fe;上述相变诱发塑性钢板的应变硬化指数分别为0.323、0.325、0.310、0.242,厚度为1mm,FLD0分别为0.576、0.580、0.553、0.432。
CN2009100122578A 2009-06-26 2009-06-26 一种建立相变诱发塑性钢板成形极限图的模型方法 Expired - Fee Related CN101599094B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100122578A CN101599094B (zh) 2009-06-26 2009-06-26 一种建立相变诱发塑性钢板成形极限图的模型方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100122578A CN101599094B (zh) 2009-06-26 2009-06-26 一种建立相变诱发塑性钢板成形极限图的模型方法

Publications (2)

Publication Number Publication Date
CN101599094A CN101599094A (zh) 2009-12-09
CN101599094B true CN101599094B (zh) 2011-04-06

Family

ID=41420537

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100122578A Expired - Fee Related CN101599094B (zh) 2009-06-26 2009-06-26 一种建立相变诱发塑性钢板成形极限图的模型方法

Country Status (1)

Country Link
CN (1) CN101599094B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102661899A (zh) * 2012-05-07 2012-09-12 同济大学 金属板材成形极限图的建立和使用方法
CN108398339A (zh) * 2017-02-05 2018-08-14 鞍钢股份有限公司 一种基于成形安全裕度的汽车钢板冲压性能判别方法
CN112926173B (zh) * 2019-12-06 2024-03-01 上海梅山钢铁股份有限公司 一种热轧高强钢板成形极限图的计算方法
CN111639405B (zh) * 2020-05-29 2022-07-29 燕山大学 一种板壳起皱失稳极限图的数值模拟求解及绘制方法
CN111896373B (zh) * 2020-06-30 2023-03-24 武汉上善仿真科技有限责任公司 一种测定等效塑性应变成形极限图的试验与计算方法

Also Published As

Publication number Publication date
CN101599094A (zh) 2009-12-09

Similar Documents

Publication Publication Date Title
Lee et al. Numerical modeling of dynamic recrystallization during nonisothermal hot compression by cellular automata and finite element analysis
CN101599094B (zh) 一种建立相变诱发塑性钢板成形极限图的模型方法
Meng et al. Testing of hot-finished high strength steel SHS and RHS under combined compression and bending
Shan et al. Prediction of martensitic transformation and deformation behavior in the TRIP steel sheet forming
CN102661899A (zh) 金属板材成形极限图的建立和使用方法
Dai et al. Investigation on the influence of damage to springback of U-shape HSLA steel plates
Pack et al. Numerical failure analysis of three-point bending on martensitic hat assembly using advanced plasticity and fracture models for complex loading
Jiang et al. A novel model to predict U-bending springback and time-dependent springback for a HSLA steel plate
CN102818763B (zh) 一种适合生产现场的热轧钢板残余应力计算方法
Choi et al. Determining the coefficients of a homogeneous anisotropic hardening model for ultrathin steel sheets
CN101598720A (zh) 一种建立相变诱发塑性钢板成形极限图的实验方法
CN103871098B (zh) 一种变厚板(vrb)三维热成形极限图的构建和使用方法
Feng et al. High temperature deformation behavior and constitutive modeling for 20CrMnTiH steel
Mei et al. High strain rate behaviour of stainless-clad bimetallic steel
Liu et al. A novel approach to determine plastic hardening curves of AA7075 sheet utilizing hydraulic bulging test at elevated temperature
CN112926173A (zh) 一种热轧高强钢板成形极限图的计算方法
Singh et al. Effect of process parameters on roll separating force, driving torque and end crop length during grooved hot rolling of SAE 1020 steel
Singh et al. Effect of process parameters on performance of grooved hot rolling of SAE 4340 steel bars
Wang et al. Investigation on forming methods in rubber pad forming process used for fabricating Cu/Ni clad foil microchannel
Liu et al. Local heating-aided hot blanking of quenched ultra-high-strength steel BR1500HS
Kusumi et al. Formability and FEM simulation of steel sheets in the hot stamping process
Lu et al. Springback behavior of tailor rolled blank in U-shape forming
Chen et al. Hot deformation behaviour and Hansel-Spittel constitutive model of Cr5 alloy for heavy backup roll
Luo et al. Numerical Analysis of AHSS Fracture in a Stretch‐bending Test
Sui et al. Analysis on shear deformation for high manganese austenite steel during hot asymmetrical rolling process using finite element method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110406

Termination date: 20120626