CN103698220B - 基于抛物线模型的管材成形性能自适应测试方法 - Google Patents
基于抛物线模型的管材成形性能自适应测试方法 Download PDFInfo
- Publication number
- CN103698220B CN103698220B CN201310722598.0A CN201310722598A CN103698220B CN 103698220 B CN103698220 B CN 103698220B CN 201310722598 A CN201310722598 A CN 201310722598A CN 103698220 B CN103698220 B CN 103698220B
- Authority
- CN
- China
- Prior art keywords
- value
- tubing
- bulging
- curve
- pressure
- 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
- 238000010998 test method Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 43
- 238000012360 testing method Methods 0.000 claims abstract description 39
- 230000008569 process Effects 0.000 claims abstract description 25
- 238000005482 strain hardening Methods 0.000 claims description 78
- 238000004088 simulation Methods 0.000 claims description 33
- 239000012530 fluid Substances 0.000 claims description 22
- 239000007788 liquid Substances 0.000 claims description 18
- 238000007789 sealing Methods 0.000 claims description 15
- 238000013459 approach Methods 0.000 claims description 14
- 238000002474 experimental method Methods 0.000 claims description 11
- 229910000838 Al alloy Inorganic materials 0.000 claims description 7
- 229910000975 Carbon steel Inorganic materials 0.000 claims description 7
- 239000010962 carbon steel Substances 0.000 claims description 7
- 239000010935 stainless steel Substances 0.000 claims description 7
- 229910001220 stainless steel Inorganic materials 0.000 claims description 7
- 230000008676 import Effects 0.000 claims description 6
- 238000003556 assay Methods 0.000 claims description 4
- 229910000861 Mg alloy Inorganic materials 0.000 claims description 3
- 229910001069 Ti alloy Inorganic materials 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- 238000005520 cutting process Methods 0.000 abstract description 2
- 238000011056 performance test Methods 0.000 abstract description 2
- 239000007787 solid Substances 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000009864 tensile test Methods 0.000 description 2
- 238000003466 welding Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 208000020442 loss of weight Diseases 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000026676 system process Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明公开了一种基于抛物线模型的管材成形性能自适应测试方法,属于管材成形性能测试领域。本发明首先对管材的胀形高度和内压进行实时测量;其次通过抛物线模型初步算出胀形最高处的实时壁厚,结合管材初始尺寸,进而算出相关曲线和数值;然后代入自适应有限元求解***进行模拟;最后对比试验所得的相同曲线并通过不断迭代修正,最终得到管材真实的力学性能参数。本发明只需测量一次管材壁厚,整个胀形过程连续进行,测试过程简单高效;无需切割测试管材进行壁厚测量,即可得到管材真实双向应力状态下的力学性能,所得管材力学性能参数真实可靠,为后续加工奠定了坚实基础。
Description
技术领域
本发明涉及一种管材成形性能自适应测试方法,具体说是一种基于抛物线模型的管材成形性能自适应测试方法,属于管材成形性能测试领域。
背景技术
目前,管材结构件在航空、航天工程、船舶、汽车等应用领域占有较大的比例,其对于提高航空、航天器、船舶、汽车的整体结构性能及有效减重具有极其重要的实际意义。
通过数值模拟的方法研究管材塑性成形过程,能够大量降低工艺的开发时间和成本,这种方法已经被广泛的应用。其准确性取决于管材的成形性能参数(如加工硬化系数K值、加工硬化指数n值)获得的准确性。
现有的管材力学性能测试方法主要包括:(1)、对于一些焊接管材,一般通过测量其卷焊前的原始板材的拉伸力学性能,代替焊接管材的力学性能。此方法没有考虑轧制方向和加工硬化的影响效果,使得所测得的力学性能与管材的真实力学性能相差较大。(2)、在管材上沿着轴向直接截取拉伸试样,进行单向拉伸实验。这种方法只能测得管材轴向力学性能,无法反映环向或者其他方向的力学性能。如果在环向或者其他方向上截取试样,因为管材呈曲面形状,需要展平试样才能进行单向拉伸实验,但是展平过程中的加工硬化必然使试样力学性能发生变化,管材直径越小影响越大,从而无法得到管材准确的力学性能参数。(3)、在管材上沿着环向截取拉伸试样,进行环向拉伸实验。这种方法可以测得管材的环向力学性能,但是试样与夹具摩擦较大,对测得的管材力学性能影响较大,而且管材塑性变形时,是受双向载荷作用的一种综合状态,无法用单一纵向拉伸实验或者环向拉伸试样进行表征。(4)、中断实验法。当管材发生一定量的胀形变形后中断实验,对管材的壁厚进行测量。测量后,将管材进行重新密封、胀形,然后再中断实验进行壁厚的测量。如此反复,直到管材发生破裂。该方法,虽可以获得不同胀形阶段管材的壁厚,但是由于管材经历了多次的加载、卸载、再加载的过程,管材的力学性能已经发生变化,因此所得结果并不可靠。而且,其实验操作繁琐,效率低下。
2013年4月17日,中国发明专利申请CN2013100258589公开了一种基于壁厚线性模型的管材力学性能液压胀形测试方法,其通过破裂点壁厚的测量,假设出壁厚线性模型,即每个时刻的壁厚等于原始壁厚减去胀形高度乘的常数倍,直接计算管材各个时刻的壁厚,再通过计算获得管材的等效应力应变曲线,确定管材材料的K值、n值和抗拉强度。该发明方法中的关键点在于假设的壁厚线性模型,但是其却只是众多模型中的一个,无法验证该模型的准确性,并且其测试结果没有修正功能,与实际的管材力学性能差距较大。
发明内容
本发明提供了一种基于抛物线模型的管材成形性能自适应测试方法,旨在解决现有的管材力学性能测试方法中无法准确得到管材真实的力学性能,影响了管材塑性成形过程的有限元模拟和实际成形的技术问题。
为了解决上述技术问题。本发明提供了一种基于抛物线模型的管材成形性能自适应测试方法,包括以下步骤:
1)、测量测试管材的初始壁厚th0、外半径RO0、内半径RI0、胀形区长度l0、二分之一的胀形区长度w=l0/2;
2)、将测试管材放置在上模和下模之间,并通过左密封冲头和右密封冲头对测试管材的两端进行密封;
3)、向测试管材内充入高压液体介质,高压液体介质的压力范围是1-350Mpa;
4)、通过传感器实时记录测试管材内各个时刻的高压液体的压力P(t)、最高点的胀形高度BHout(t);
5)、根据抛物线模型的几何关系,计算获得胀形过程中每一时刻的壁厚th(t):
外部轮廓抛物线方程,即Rout(z,t):
内部轮廓抛物线方程,即Rin(z,t):
实时厚度th(t)变化关系式:
th(t)=BHout(t)+tho-BHin(t)(3)
由积不变原则所得关系式,即实时管材体积Vtube等于原始管材体积V0:
将(1)(2)(3)式代入(4),消去Rout、Rin、BHout(t)和w,整理得到:
6)、通过以下计算获得管材的等效应力应变曲线:
根据米塞斯等效应力应变模型:
环向应变:
厚度方向应变:
环向应力:
轴向应力:
其中R1为轴向曲率半径:
R2为环向曲率半径:
7)、根据等效应力应变曲线,确定测试管材的加工硬化系数K值、加工硬化指数n值;
8)、将步骤7)确定的加工硬化系数K值、加工硬化指数n值代入自适应有限元求解***,对在此力学性能参数情况下的测试管材进行自由胀形模拟,得到模拟的液体压力和胀形最高点的高度关系曲线;
9)、利用线性插值的方法将步骤8)得到的模拟的压力和高度关系曲线数值数组补全,使其能够与测试所得的液体压力和胀形高度关系曲线数值数组进行比较;
10)、通过判断两条曲线的边界数值相减后的乘积,进行判断曲线是否相交,并根据相同内压下,胀形高度数值差异的绝对值总和的平均值比较范围的最高边界处真实压力与计算压力的差异Xe、比较范围的最低边界处真实压力与计算压力的差异X0、逼近容差R,对加工硬化系数K值、加工硬化指数n值进行修改,如下列公式;
式中,Ki是第i次模拟的加工硬化系数K值,ni是第i次模拟的加工硬化指数n值,同理,Ki+1是第i+1次模拟的加工硬化系数K值,ni+1是第i次模拟的加工硬化指数n值;
11)、将修改后的加工硬化系数K值、加工硬化指数n值再次代入自适应有限元求解***,重复步骤8)、9)、10),直到模拟出的液体压力和胀形高度关系曲线,无限接近试验所得的液体压力和胀形高度关系曲线,取此时的K值、n值为管材的真实力学性能参数。
本发明中,所述步骤1)中测试管材为不锈钢管材、铝合金管材、碳钢管材、镁合金管材或钛合金管材。
本发明中,所述步骤8)中的自适应有限元求解***流程为:
1)、在自适应有限元求解***中导入管材胀形实验中的压力-位移实验数据,在主界面中输入原始管材的外直径、内直径、壁厚、模具内间距,构建所述步骤1)到步骤7)的解析模型,求得管材真实-应变曲线;
2)、用最小二乘法求得K值、n值,在***界面中选择两端固定胀形或者是自由胀形,并输入加载路径、管材长度、模具圆角半径,进而生成有限元模型文件;
3)、导入有限元模型文件,将解析值作为模型的初始值代入到模型中,激活有限元求解器,获得模拟P-H曲线i、Ki、ni;
4)、将获得模拟P-H曲线i、Ki、ni与试验所得的P-H曲线相比较,是否匹配,如果误差比设定的允许误差大,则修正此次模拟的加工硬化系数K值、加工硬化指数n值,代入有限元模型中,重新模拟计算,直到所得模拟P-H曲线与实验所得P-H曲线无限接近,误差在允许误差范围内,则停止循环,输出此时的加工硬化系数K值、加工硬化指数n值。
本发明的有益效果在于:(1)、本发明只需在开始前测量一次管材壁厚,不需要测量破裂处壁厚,胀形过程中的任意时刻的壁厚通过理论计算获得,整个胀形过程连续进行,任意时刻的液体压力和胀形高度可以连续的测量和记录,简单易行,极大地提高了实验效率。(2)、本发明无需切割测试管材进行壁厚测量,即可得到管材真实双向应力状态下的力学性能;通过K值、n值的修正和反复迭代,使最终得到的模拟液体压力和胀形高度关系曲线无限接近试验所得真实液体压力和胀形高度关系曲线,进而为管材力学性能参数真实性,提供了可靠的参照对象和检验结果,为后续加工奠定了坚实基础。
附图说明
图1为本发明基于抛物线模型的管材力学性能自适应有限元测试方法流程图;
图2为本发明中自适应有限元求解***流程图;
图3为本发明中实验装置示意图及剖面图;
图4为本发明中抛物线模型示意图;
图中,1-测试管材,2-上模,3-下模,4-左密封冲头,5-右密封冲头,6-传感器,th0-初始壁厚、RO0-原始管材的外半径、RI0-原始管材的内半径、l0-胀形区长度、Bout-管材外部实时半径、Bin-管材内部实时半径,w-二分之一的胀形区长度、z-管材轴向方向、r-为管材径向方向、BHout(t)-外壁最高点的胀形高度、BHin(t)-内壁最高点的胀形高度、t-管壁实时厚度。
具体实施方式
下面结合附图1、2、3、4及以分别测试不锈钢管材、铝合金管材、碳钢管材的力学性能为具体实施例详细说明本发明基于抛物线模型的管材力学性能自适应有限元测试方法。
实施例1
第一步,选择直径为100mm,壁厚为1mm,长度为500mm的不锈钢管材1,测量该管材的初始壁厚th0=1mm、外半径RO0=50mm、内半径RI0=49mm、胀形区长度l0=100mm,二分之一的胀形区长度w=l0/2=50mm,z是图4中抛物线模型中X(z)轴的横坐标,在运算中可以抵消掉,不影响运算过程;
第二步,将不锈钢管材1放置在上模2和下模3之间,合模后,通过左密封冲头4和右密封冲头5对其两端进行密封,如图3所示;
第三步、通过液压泵提供高压液体介质,沿着与液压泵相连接的密封冲头内部的导油槽,向不锈钢管材1内充入高压液体介质,高压液体的压力范围是1-350Mpa;
第四步、通过传感器6实时记录不锈钢管材1内各个时刻的高压液体的压力P(t)、最高点的胀形高度BHout(t),例如某一时刻P(t)=12Mpa,BHout(t)=0.5mm,最高点为原点截面处,即z=0;
第五步、根据抛物线模型的几何关系,如图4中所示,计算获得胀形过程中每一时刻的壁厚th(t),根据步骤四中的某一时刻内压、高度值为例:
外部轮廓抛物线方程,即Rout(z,t):
内部轮廓抛物线方程,即Rin(z,t):
实时厚度th(t)变化关系式:
th(t)=BHout(t)+tho-BHin(t)(3)
由积不变原则所得关系式,即实时管材体积Vtube等于原始管材体积V0:
将(1)(2)(3)式代入(4),消去Rout、Rin、BHout(t)和w,整理得到:
由此式可得th(t)=0.9899mm;
第六步、通过以下计算获得管材的等效应力应变曲线:
根据米塞斯等效应力应变模型:
环向应变:
厚度方向应变:
环向应力:
轴向应力:
其中R1为轴向曲率半径:
R2为环向曲率半径:
第七步、根据等效应力应变曲线,采用本行业内公认的计算方法,即将等效应力应变进行对数化处理,选取若干个点,进行最小二乘直线拟合,拟合直线的截距为n,斜率为logK,进而确定待测试管材的加工硬化系数K值、加工硬化指数n值;
第八步、将上述加工硬化系数K值、加工硬化指数n值代入自适应有限元求解***,利用其中的有限元模拟软件,对在此力学性能参数情况下的管材,进行自由胀形模拟,得到模拟的液体压力和胀形最高点的高度关系曲线,如图2所示;
第九步、利用线性插值的方法将模拟的压力和高度关系曲线数值数组补全,使其能够与测试所得的液体压力和胀形高度关系曲线数值数组进行比较;
第十步、通过判断上述两条曲线的边界数值相减后的乘积,进行判断曲线是否相交,并根据相同内压下,胀形高度数值差异的绝对值总和的平均值比较范围的最高边界处真实压力与计算压力的差异Xe、比较范围的最低边界处真实压力与计算压力的差异X0、逼近容差R,对加工硬化系数K值、加工硬化指数n值进行修改,如以公式:
式中,Ki是第i次模拟的加工硬化系数K值,ni是第i次模拟的加工硬化指数n值,同理,Ki+1是第i+1次模拟的加工硬化系数K值,ni+1是第i次模拟的加工硬化指数n值;
第十一步、将修改后的加工硬化系数K值、加工硬化指数n值再次代入自适应有限元求解***,重复步骤八、步骤九、步骤十,直到模拟出的液体压力和胀形高度关系曲线,无限接近试验所得的液体压力和胀形高度关系曲线,取此时的K值、n值为管材的真实力学性能参数。
实施例2
第一步,选择直径为100mm,壁厚为1mm,长度为500mm的铝合金管材1,测量该管材的初始壁厚th0=1mm、外半径RO0=50mm、内半径RI0=49mm、胀形区长度l0=100mm,即w=l0/2=50mm,z是图4中抛物线模型中X(z)轴的横坐标,在运算中可以抵消掉,不影响运算过程;
第二步,将铝合金管材1放置在上模2和下模3之间,合模后,通过左密封冲头4和右密封冲头5对其两端进行密封,如图3所示;
第三步、通过液压泵提供高压液体介质,沿着与液压泵相连接的密封冲头内部的导油槽,向铝合金管材1内充入高压液体介质,高压液体的压力范围是1-350Mpa;
第四步、通过传感器6实时记录铝合金管材1内各个时刻的高压液体的压力P(t)、最高点的胀形高度BHout(t),例如某一时刻P(t)=10Mpa,BHout(t)=0.5mm,最高点为原点截面处,即z=0;
第五步、根据抛物线模型的几何关系,如图4中所示,计算获得胀形过程中每一时刻的壁厚th(t),根据步骤四中的某一时刻内压、高度值为例:
外部轮廓抛物线方程,即Rout(z,t):
内部轮廓抛物线方程,即Rin(z,t):
实时厚度th(t)变化关系式:
th(t)=BHout(t)+tho-BHin(t)(3)
由积不变原则所得关系式,即实时管材体积Vtube等于原始管材体积V0:
将(1)(2)(3)式代入(4),消去Rout、Rin、BHout(t)和w,整理得到:
由此式可得th(t)=0.9899mm;
第六步、通过以下计算获得管材的等效应力应变曲线:
根据米塞斯等效应力应变模型:
环向应变:
厚度方向应变:
环向应力:
轴向应力:
其中R1为轴向曲率半径:
R2为环向曲率半径:
第七步、根据等效应力应变曲线,采用本行业内公认的计算方法,即将等效应力应变进行对数化处理,选取若干个点,进行最小二乘直线拟合,拟合直线的截距为n,斜率为logK,进而确定待测试管材的加工硬化系数K值、加工硬化指数n值;
第八步、将上述加工硬化系数K值、加工硬化指数n值代入自适应有限元求解***,利用其中的有限元模拟软件,对在此力学性能参数情况下的管材,进行自由胀形模拟,得到模拟的液体压力和胀形最高点的高度关系曲线,如图2所示;
第九步、利用线性插值的方法将模拟的压力和高度关系曲线数值数组补全,使其能够与测试所得的液体压力和胀形高度关系曲线数值数组进行比较;
第十步、通过判断上述两条曲线的边界数值相减后的乘积,进行判断曲线是否相交,并根据相同内压下,胀形高度数值差异的绝对值总和的平均值比较范围的最高边界处真实压力与计算压力的差异Xe、比较范围的最低边界处真实压力与计算压力的差异X0、逼近容差R,对加工硬化系数K值、加工硬化指数n值进行修改,如以下公式:
式中,Ki是第i次模拟的加工硬化系数K值,ni是第i次模拟的加工硬化指数n值,同理,Ki+1是第i+1次模拟的加工硬化系数K值,ni+1是第i次模拟的加工硬化指数n值;
第十一步、将修改后的加工硬化系数K值、加工硬化指数n值再次代入自适应有限元求解***,重复步骤八、步骤九、步骤十,直到模拟出的液体压力和胀形高度关系曲线,无限接近试验所得的液体压力和胀形高度关系曲线,取此时的K值、n值为管材的真实力学性能参数。
实施例3
第一步,选择直径为100mm,壁厚为1mm,长度为500mm的碳钢管材1,测量该管材的初始壁厚th0=1mm、外半径RO0=50mm、内半径RI0=49mm、胀形区长度l0=100,即w=l0/2=50mm,z是图4中抛物线模型中X(z)轴的横坐标,在运算中可以抵消掉,不影响运算过程;
第二步,将碳钢管材1放置在上模2和下模3之间,合模后,通过左密封冲头4和右密封冲头5对待测试管材1的两端进行密封,如图3所示;
第三步、通过液压泵提供高压液体介质,沿着与液压泵相连接的密封冲头内部的导油槽,向碳钢管材1内充入高压液体介质,高压液体的压力范围是1-350Mpa;
第四步、通过传感器6实时记录碳钢管材1内各个时刻的高压液体的压力P(t)、最高点的胀形高度BHout(t),例如某一时刻P(t)=14Mpa,BHout(t)=0.5mm,最高点为原点截面处,即z=0;
第五步、根据抛物线模型的几何关系,如图4中所示,计算获得胀形过程中每一时刻的壁厚th(t),根据步骤四中的某一时刻内压、高度值为例:
外部轮廓抛物线方程,即Rout(z,t):
内部轮廓抛物线方程,即Rin(z,t):
实时厚度th(t)变化关系式:
th(t)=BHout(t)+tho-BHin(t)(3)
由积不变原则所得关系式,即实时管材体积Vtube等于原始管材体积V0:
将(1)(2)(3)式代入(4),消去Rout、Rin、BHout(t)和w,整理得到:
由此式可得th(t)=0.9899mm;
第六步、通过以下计算获得管材的等效应力应变曲线:
根据米塞斯等效应力应变模型:
环向应变:
厚度方向应变:
环向应力:
轴向应力:
其中R1为轴向曲率半径:
R2为环向曲率半径:
第七步、根据等效应力应变曲线,采用本行业内公认的计算方法,即将等效应力应变进行对数化处理,选取若干个点,进行最小二乘直线拟合,拟合直线的截距为n,斜率为logK,进而确定待测试管材的加工硬化系数K值、加工硬化指数n值;
第八步、将上述加工硬化系数K值、加工硬化指数n值代入自适应有限元求解***,利用其中的有限元模拟软件,对在此力学性能参数情况下的管材,进行自由胀形模拟,得到模拟的液体压力和胀形最高点的高度关系曲线,如图2所示;
第九步、利用线性插值的方法将模拟的压力和高度关系曲线数值数组补全,使其能够与测试所得的液体压力和胀形高度关系曲线数值数组进行比较;
第十步、通过判断上述两条曲线的边界数值相减后的乘积,进行判断曲线是否相交,并根据相同内压下,胀形高度数值差异的绝对值总和的平均值比较范围的最高边界处真实压力与计算压力的差异Xe、比较范围的最低边界处真实压力与计算压力的差异X0、逼近容差R,对加工硬化系数K值、加工硬化指数n值进行修改,如以下公式:
式中,Ki是第i次模拟的加工硬化系数K值,ni是第i次模拟的加工硬化指数n值,同理,Ki+1是第i+1次模拟的加工硬化系数K值,ni+1是第i次模拟的加工硬化指数n值;
第十一步、将修改后的加工硬化系数K值、加工硬化指数n值再次代入自适应有限元求解***,重复步骤八、步骤九、步骤十,直到模拟出的液体压力和胀形高度关系曲线,无限接近试验所得的液体压力和胀形高度关系曲线,取此时的K值、n值为管材的真实力学性能参数。
上述实施例1至3中自适应有限元求解***流程如下:
第一步、打开自适应有限元求解***软件,进入***主界面,根据导入按键提示,导入管材胀形实验中的压力—位移实验数据,在主界面中输入原始管材的外直径、内直径、壁厚、模具内间距,构建上述实施例中步骤一到步骤七的解析模型,求得管材真实应力—应变曲线;
第二步、用最小二乘法求得K值、n值,在***界面中选择两端固定胀形或者是自由胀形,并输入加载路径、管材长度、模具圆角半径,进而生成有限元模型文件;
第三步、导入有限元模型文件,将解析值作为模型的初始值代入到模型中,激活有限元求解器,获得模拟P-H曲线i、Ki、ni,其中i是第i次模拟结果,Ki是第i次模拟结果的加工硬化系数K值,ni是第i次模拟结果的加工硬化指数n值;
第四步、将获得模拟P-H曲线i、Ki、ni与试验所得的P-H曲线相比较,是否匹配,如果误差比设定的允许误差大,则修正此次模拟的加工硬化系数K值、加工硬化指数n值,代入有限元模型中,重新模拟计算,直到所得模拟P-H曲线与实验所得P-H曲线无限接近,误差在允许误差范围内,则停止循环,输出此时的加工硬化系数K值、加工硬化n值。
上述方法同样适用于镁合金管材、钛合金管材和其他常见管材的力学性能测试,步骤与上实施例相同,在此不作进一步说明。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以做出若干改进,这些改进也应视为本发明的保护范围。
Claims (3)
1.一种基于抛物线模型的管材成形性能自适应测试方法,其特征在于包括以下步骤:
1)、测量测试管材(1)的初始壁厚th0、外半径RO0、内半径RI0、胀形区长度l0、二分之一的胀形区长度w=l0/2;
2)、将测试管材(1)放置在上模(2)和下模(3)之间,并通过左密封冲头(4)和右密封冲头(5)对测试管材(1)的两端进行密封;
3)、向测试管材(1)内充入高压液体介质,高压液体介质的压力范围是1-350Mpa;
4)、通过传感器(6)实时记录测试管材(1)内各个时刻的高压液体的压力P(t)、最高点的胀形高度BHout(t);
5)、根据抛物线模型标注的几何关系,计算获得胀形过程中每一时刻的壁厚th(t):
外部轮廓抛物线方程,即Rout(z,t):
式(1)中,z为管材在抛物线模型中X(z)轴的横坐标;
内部轮廓抛物线方程,即Rin(z,t):
式(2)中,BHin(t)为管材内壁到胀形最高点内壁的垂直距离;
实时厚度th(t)变化关系式:
th(t)=BHout(t)+tho-BHin(t)(3)
由体积不变原则所得关系式,即实时管材体积Vtube等于原始管材体积V0:
将(1)(2)(3)式代入(4),消去Rout、Rin、BHout(t)和w,整理得到:
6)、通过以下计算获得管材的等效应力应变曲线:
根据米塞斯等效应力应变模型:
环向应变:
厚度方向应变:
环向应力:
轴向应力:
其中R1为轴向曲率半径:
R2为环向曲率半径:
7)、根据等效应力应变曲线,确定测试管材(1)的加工硬化系数K值、加工硬化指数n值;
8)、将步骤7)确定的加工硬化系数K值、加工硬化指数n值代入自适应有限元求解***,对在此力学性能参数情况下的测试管材(1)进行自由胀形模拟,得到模拟的液体压力和胀形最高点的高度关系曲线;
9)、利用线性插值的方法将步骤8)得到的模拟的压力和高度关系曲线数值数组补全,使其能够与测试所得的液体压力和胀形高度关系曲线数值数组进行比较;
10)、通过判断两条曲线的边界数值相减后的乘积,进行判断曲线是否相交,并根据相同内压下,胀形高度数值差异的绝对值总和的平均值比较范围的最高边界处真实压力与计算压力的差异Xe、比较范围的最低边界处真实压力与计算压力的差异X0、逼近容差R,对加工硬化系数K值、加工硬化指数n值进行修改,如下列公式;
式中,Ki是第i次模拟的加工硬化系数K值,ni是第i次模拟的加工硬化指数n值,同理,Ki+1是第i+1次模拟的加工硬化系数K值,ni+1是第i次模拟的加工硬化指数n值;
11)、将修改后的加工硬化系数K值、加工硬化指数n值再次代入自适应有限元求解***,重复步骤8)、9)、10),直到模拟出的液体压力和胀形高度关系曲线,无限接近试验所得的液体压力和胀形高度关系曲线,取此时的K值、n值为管材的真实力学性能参数。
2.根据权利要求1所述的基于抛物线模型的管材成形性能自适应测试方法,其特征在于:所述步骤1)中测试管材(1)为不锈钢管材、铝合金管材、碳钢管材、镁合金管材或钛合金管材。
3.根据权利要求1或2所述的基于抛物线模型的管材成形性能自适应测试方法,其特征在于所述步骤8)中的自适应有限元求解***流程为:
1)、在自适应有限元求解***中导入管材胀形实验中的压力-位移实验数据,在主界面中输入原始管材的外直径、内直径、壁厚、模具内间距,构建所述步骤1)到步骤7)的解析模型,求得管材真实-应变曲线;
2)、用最小二乘法求得K值、n值,在***界面中选择两端固定胀形或者是自由胀形,并输入加载路径、管材长度、模具圆角半径,进而生成有限元模型文件;
3)、导入有限元模型文件,将解析值作为模型的初始值代入到模型中,激活有限元求解器,获得模拟P-H曲线i、Ki、ni;
4)、将获得模拟P-H曲线i、Ki、ni与试验所得的P-H曲线相比较,是否匹配,如果误差比设定的允许误差大,则修正此次模拟的加工硬化系数K值、加工硬化指数n值,代入有限元模型中,重新模拟计算,直到所得模拟P-H曲线与实验所得P-H曲线无限接近,误差在允许误差范围内,则停止循环,输出此时的加工硬化系数K值、加工硬化指数n值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310722598.0A CN103698220B (zh) | 2013-12-24 | 2013-12-24 | 基于抛物线模型的管材成形性能自适应测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310722598.0A CN103698220B (zh) | 2013-12-24 | 2013-12-24 | 基于抛物线模型的管材成形性能自适应测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103698220A CN103698220A (zh) | 2014-04-02 |
CN103698220B true CN103698220B (zh) | 2016-01-27 |
Family
ID=50359824
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310722598.0A Active CN103698220B (zh) | 2013-12-24 | 2013-12-24 | 基于抛物线模型的管材成形性能自适应测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103698220B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110553917B (zh) * | 2019-09-11 | 2022-03-18 | 北京东方雨虹防水技术股份有限公司 | 一种防水材料多向拉伸的试验***及试验方法 |
CN111238699A (zh) * | 2019-12-13 | 2020-06-05 | 中国石油化工股份有限公司 | 变形管道应力测试确定方法 |
CN113838396B (zh) * | 2020-06-24 | 2024-04-16 | 深圳市国显科技有限公司 | 一种vcom的烧录方法 |
CN112287486B (zh) * | 2020-10-30 | 2024-03-19 | 江苏科技大学 | 一种潜水器圆柱壳无模自由胀形过程的性能预测方法 |
CN114433705B (zh) * | 2021-12-28 | 2023-12-26 | 台州通禾流体控制股份有限公司 | 一种管材内高压成形的起皱智能控制方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU459705A1 (ru) * | 1972-11-04 | 1975-02-05 | Уральский политехнический институт им.С.М.Кирова | Устройство дл измерени деформации грунта |
CN101793647A (zh) * | 2010-02-04 | 2010-08-04 | 西北工业大学 | 薄壁管件周向力学性能测试方法 |
CN103048203A (zh) * | 2013-01-23 | 2013-04-17 | 哈尔滨工业大学 | 基于壁厚线性模型的管材力学性能液压胀形测试方法 |
-
2013
- 2013-12-24 CN CN201310722598.0A patent/CN103698220B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU459705A1 (ru) * | 1972-11-04 | 1975-02-05 | Уральский политехнический институт им.С.М.Кирова | Устройство дл измерени деформации грунта |
CN101793647A (zh) * | 2010-02-04 | 2010-08-04 | 西北工业大学 | 薄壁管件周向力学性能测试方法 |
CN103048203A (zh) * | 2013-01-23 | 2013-04-17 | 哈尔滨工业大学 | 基于壁厚线性模型的管材力学性能液压胀形测试方法 |
Non-Patent Citations (3)
Title |
---|
管材胀形工艺分类及其变形力学特征;王同海;《锻压技术》;19991231(第4期);第30-32页 * |
管材自由胀形时极限载荷及成形极限的确定;杨连发;《塑性工程学报》;20060228;第13卷(第1期);第13-17页 * |
管材自由胀形时胀形区轮廓形状的影响因素;林艳丽;《金属学报》;20100630;第46卷(第6期);第729-735页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103698220A (zh) | 2014-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103698220B (zh) | 基于抛物线模型的管材成形性能自适应测试方法 | |
CN103048203A (zh) | 基于壁厚线性模型的管材力学性能液压胀形测试方法 | |
Netto et al. | On the effect of corrosion defects on the collapse pressure of pipelines | |
Błachut | Buckling of externally pressurized steel toriconical shells | |
Lianfa et al. | Determination of stress–strain relationship of tubular material with hydraulic bulge test | |
Dodge et al. | Stress indices and flexibility factors for moment loadings on elbows and curved pipe | |
CN104019790B (zh) | 水下耐压结构长时间累积变形的测量装置及方法 | |
CN110763568B (zh) | 一种管材任意方向厚向异性系数的确定方法 | |
Zhang et al. | Buckling properties of bulged barrels under external pressure | |
Guo et al. | Experimental and numerical investigation of bulging behaviour of hyperelastic textured tubes | |
Yu et al. | The cross-over mechanisms of integral buckle arrestors for offshore pipelines | |
He et al. | A novel test method for continuous nonlinear biaxial tensile deformation of sheet metals by bulging with stepped-dies | |
Kiani et al. | Numerical Modeling and Analytical Investigation of Autofrettage Process on the Fluid End Module of Fracture Pumps | |
Osadchii et al. | Shaping of thick sheet in the production of welded large-diameter pipe | |
Gelin et al. | Modelling and control of hydroforming processes for flanges forming | |
Karabegović et al. | Experimental modeling of fluid pressure during hydroforming of welded plates | |
CN208313506U (zh) | 一种气体压力传感器标定装置 | |
Tripathi et al. | Pressure drop analysis for hydraulic valves | |
CN110763566A (zh) | 一种各向异性管材环向厚向异性系数的确定方法 | |
Zhao et al. | Failure Analysis of Large-Diameter Coiled Tubing Based on Diameter Growth | |
Gholipour et al. | Severity of the bend and its effect on the subsequent hydroforming process for aluminum alloy tube | |
CN215448783U (zh) | 一种压力敏感元件形变试验装置 | |
Rasty et al. | Isothermal and thermomechanical finite-element analysis of the tube drawing process using a fixed tapered plug | |
Baharom et al. | Combined hardening parameters of high strength steel under low cycle fatigue | |
CN118111814B (zh) | 一种气液两相混合物压缩性能精确测定的方法 |
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 | ||
CB03 | Change of inventor or designer information |
Inventor after: Tao Jie Inventor after: Guo Xunzhong Inventor after: Jin Xianyu Inventor after: Fan Minyu Inventor after: Tao Jieguo Xunzhong, Jin Xianyu, Fan Minyu, Qitong Baihui Inventor before: Jin Xianyu Inventor before: Fan Minyu Inventor before: Qitong Baihui Inventor before: Guo Xunzhong Inventor before: Tao Jie |
|
CB03 | Change of inventor or designer information |