一种厚板弯曲成形过程中中性层偏移量的计算方法与流程

1.本发明属于弯曲成形的技术领域，具体涉及一种厚板在弯曲成形过程中考虑中性层偏移的计算新方法。


背景技术：

2.厚板金属弯曲成形是制备厚壁直缝焊管的重要变形过程，例如uoe、jcoe、rbe等成形方式中厚板金属弯曲成形是制定成形工艺的重点研究部分。但弯曲成形方式在成形过程中的控制比较粗糙，对生产效率、钢管焊接及外观质量都产生了一定的影响。实际生产过程中更多仍然依靠经验进行模具设计和工艺参数调整。而且对焊管生产过程中复杂的变形认识相对较浅，理论分析相对不足，生产上还存在很多问题。同时管坯成形的多道次压弯过程是一个复杂的非线性过程，包含弹塑性材料的非线性、大挠度的几何非线性、边界条件的非线性，而且在成形过程中发生的一些物理现象，如金属的流动、接触面上的摩擦、变形过程中应力应变变化情况等。因此考虑中性层偏移对厚板弯曲成形过程中的影响，推导出理论中性层偏移量的计算公式。从而提高弯曲成形工艺，降低生产成本，提高产品的质量。


技术实现要素：

3.鉴于此，本发明的目的在于，提供一种厚板弯曲成形过程中中性层偏移量的计算方法，能够降低生产成本，提高厚壁直缝焊管的产品精度。
4.为了达到上述发明目的，进而采取的技术方案如下：一种厚板弯曲成形过程中中性层偏移量的计算方法，包括以下步骤：步骤1：利用分层拉伸断裂实验，获得厚向强度线性分布曲线，根据厚板弯曲理论，建立应力-应变数学模型：
ꢀꢀꢀꢀ
（1）k
‑‑
硬化系数，n
‑‑
硬化指数，q
‑‑
反映厚向强度分布不均匀程度的系数，z
‑‑
板层距离中间层的垂直距离，e
‑‑
理论中性层偏移量，t
‑‑
厚板的厚度，e
‑‑
弹性模量；步骤2：在拉伸区和压缩区分别取微元体进行受力平衡分析，并在极限状态下得到中性层偏移公式：厚板弯曲的拉伸区的力学平衡关系式如下：
ꢀꢀꢀ
（2）厚板弯曲的压缩区的力学平衡关系式如下：
ꢀꢀꢀꢀ
（3） σx
‑‑
第一主应力，σy
‑‑
第二主应力；其中：β
‑‑
中间主应力影响系数；ρ
‑‑
中性层的曲率半径，α——微分单元体的长度，ρ
ya
——微分单元体在压缩区的曲率半径，ρ
la
——微分单元体在拉伸区的曲率半径；当中性层发生塑形变形时，拉伸区和压缩区的应力相等，即：ρ
la=
ρ
ya
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）根据公式（2）、（3）（4）可推导出理论中性层偏移量的计算公式：
ꢀꢀꢀꢀꢀꢀꢀꢀ
（5）ρ
‑‑
中性层的曲率半径，ρ
min
‑‑
最小弯曲曲率半径，ρ
max
‑‑
最大弯曲曲率半径。
5.对于给定材料，工人只需测量弯曲半径、板厚即可利用公式（5）轻松计算出中性层偏移量，方便工人实际操作。
6.本发明的有益效果是：1. 本发明提出一种确定中性层偏移量的方法，并给出计算公式。通过中性层偏移量的确定，可以协助制定更为准确的弯曲工艺参数。
7.2.本发明提出的确定中性层偏移量的方法可以提高厚壁直缝焊管的成形精度，降低生产成本。
附图说明
8.构成本技术的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：图1为弯曲成形的模型的示意图；图2为厚板弯曲的应力平衡关系图。
具体实施方式
9.需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
10.为了使本技术领域的人员更好地理解本技术方案，下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分的实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本技术保护的范围。
11.实施例1以厚板的jco成形为例：如图1-2所示，一种厚板jco成形过程中中性层偏移量的计算方法，包括以下步骤：步骤1：利用分层拉伸断裂实验，获得厚向强度线性分布曲线，根据厚板弯曲理论，建立应力-应变数学模型：
ꢀꢀꢀꢀ
（1）k
‑‑
硬化系数，n
‑‑
硬化指数，q
‑‑
反映厚向强度分布不均匀程度的系数，z
‑‑
板层距离中间层的垂直距离，e
‑‑
理论中性层偏移量，t
‑‑
厚板的厚度，e
‑‑
弹性模量；步骤2：在拉伸区和压缩区分别取微元体进行受力平衡分析，并在极限状态下得到中性层偏移公式：厚板弯曲的拉伸区的力学平衡关系式如下：
ꢀꢀꢀ
（2）厚板弯曲的压缩区的力学平衡关系式如下：
ꢀꢀꢀꢀ
（3） σx
‑‑
第一主应力，σy
‑‑
第二主应力；其中：β
‑‑
中间主应力影响系数；ρ
‑‑
中性层的曲率半径，α——微分单元体的长度，ρ
ya
——微分单元体在压缩区的曲率半径，ρ
la
——微分单元体在拉伸区的曲率半径；当中性层发生塑形变形时，拉伸区和压缩区的应力相等，即：
ρ
la=
ρ
ya
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
（4）根据公式（2）、（3）（4）可推导出理论中性层偏移量的计算公式：
ꢀꢀꢀꢀꢀꢀꢀꢀ
（5）ρ
‑‑
中性层的曲率半径，ρ
min
‑‑
最小弯曲曲率半径，ρ
max
‑‑
最大弯曲曲率半径。
12.实施例2为了验证本发明，选用厚板材料为镍基合金为例。
13.厚板的基本参数为：厚板的宽度w=1000mm，厚板的长度l=2000mm，厚板的厚度t=60mm，杨氏模量e=205gpa，泊松比为0.3，密度为7.82 g/cm3。上模曲率半径为150mm，下模曲率半径为150mm，下模跨距为770mm。上模和下模设置为刚体。其压下量设为30mm。
14.从模拟结果中的应变云图选取厚板中间两个单元体之间的结点，通过测量可以的得到：模型在压下量为30mm时拉伸区的弹性高度y1为2.10mm，压缩区的弹性高度y2为7.47mm，通过计算可以得到模拟中性层偏移量e
ˊ
=4.79mm。
15.根据模拟结果可以得到板层距离中间层的垂直距离z为7.47mm，t=60mm，q设为0.2，ρ=172.53mm，ρ
min
=150mm，ρ
max =210mm，将其代入公式（5）中，可以得到理论中性层偏移量e=3.00mm。在模拟中的应变云图中相隔一个单元体选取两个单元体之间的结点，通过测量以及代入公式（5）中可以得到在厚板的厚度为60mm，压下量为30mm下的模拟中性层偏移数据和理论中性层偏移数据。通过计算模拟中性层偏移量与理论偏移量数据之间误差为1.62%。因此对于给定材料，工人只需测量弯曲半径、板厚即可利用本发明所述公式（5）轻松计算出中性层偏移量，方便工人实际操作。本发明提出的确定中性层偏移量的方法可以提高厚壁直缝焊管的成形精度，降低生产成本。
16.利用origin软件对所得到数据进行曲线拟合，采用多项式拟合，模拟的中性层偏移量的拟合函数为：
ꢀꢀ
（6）其中，r=0.97288（r为曲线相关系数）。
17.理论的中性层偏移量的拟合函数为：
ꢀꢀꢀꢀꢀꢀꢀ
（7）其中，r=0.90515（r为曲线相关系数）。
18.以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进、部件拆分或组合等，均应包含在本发明的保护范围之内。