CN114117679A - 一种厚板弯曲成形过程中中性层偏移量的计算方法 - Google Patents
一种厚板弯曲成形过程中中性层偏移量的计算方法 Download PDFInfo
- Publication number
- CN114117679A CN114117679A CN202111460909.1A CN202111460909A CN114117679A CN 114117679 A CN114117679 A CN 114117679A CN 202111460909 A CN202111460909 A CN 202111460909A CN 114117679 A CN114117679 A CN 114117679A
- Authority
- CN
- China
- Prior art keywords
- neutral layer
- thick plate
- bending forming
- curvature
- radius
- 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
Links
- 230000007935 neutral effect Effects 0.000 title claims abstract description 45
- 238000005452 bending Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000008569 process Effects 0.000 title claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000009826 distribution Methods 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- 239000010410 layer Substances 0.000 claims description 41
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 229910000831 Steel Inorganic materials 0.000 claims description 5
- 239000010959 steel Substances 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 3
- 239000011229 interlayer Substances 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000013508 migration Methods 0.000 claims description 3
- 230000005012 migration Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 10
- 239000000463 material Substances 0.000 abstract description 5
- 230000008859 change Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract description 2
- 239000002184 metal Substances 0.000 description 3
- 229910052751 metal Inorganic materials 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- PXHVJJICTQNCMI-UHFFFAOYSA-N Nickel Chemical compound [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- BSYNRYMUTXBXSQ-UHFFFAOYSA-N Aspirin Chemical compound CC(=O)OC1=CC=CC=C1C(O)=O BSYNRYMUTXBXSQ-UHFFFAOYSA-N 0.000 description 1
- 241000405070 Percophidae Species 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052759 nickel Inorganic materials 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000003466 welding 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/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Bending Of Plates, Rods, And Pipes (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明公开了一种厚板弯曲成形过程中中性层偏移量的计算方法,厚板弯曲成形过程复杂多变,包含弹塑性材料的非线性,大挠度的几何非线性等以及变形过程中应力应变的变化情况。目前厚板弯曲成形中并未考虑中性层偏移的因素,其原因在于无法精确给定中性层偏移量,从而导致弯曲成形质量较差。为提高弯曲成形质量,降低生产成本,提高生产率,本发明提出一种厚板在弯曲成形过程中中性层偏移量的计算方法。利用实验测量得到厚向强化线性分布曲线,再根据弹塑性力学理论和中厚板弯曲理论,通过对厚板在弯曲时取微分单元体进行力学分析,讨论一种极限状态,推导出理论中性层偏移的计算公式。
Description
技术领域
本发明属于弯曲成形的技术领域,具体涉及一种厚板在弯曲成形过程中考虑中性层偏移的计算新方法。
背景技术
厚板金属弯曲成形是制备厚壁直缝焊管的重要变形过程,例如UOE、JCOE、RBE等成形方式中厚板金属弯曲成形是制定成形工艺的重点研究部分。但弯曲成形方式在成形过程中的控制比较粗糙,对生产效率、钢管焊接及外观质量都产生了一定的影响。实际生产过程中更多仍然依靠经验进行模具设计和工艺参数调整。而且对焊管生产过程中复杂的变形认识相对较浅,理论分析相对不足,生产上还存在很多问题。同时管坯成形的多道次压弯过程是一个复杂的非线性过程,包含弹塑性材料的非线性、大挠度的几何非线性、边界条件的非线性,而且在成形过程中发生的一些物理现象,如金属的流动、接触面上的摩擦、变形过程中应力应变变化情况等。因此考虑中性层偏移对厚板弯曲成形过程中的影响,推导出理论中性层偏移量的计算公式。从而提高弯曲成形工艺,降低生产成本,提高产品的质量。
发明内容
鉴于此,本发明的目的在于,提供一种厚板弯曲成形过程中中性层偏移量的计算方法,能够降低生产成本,提高厚壁直缝焊管的产品精度。
为了达到上述发明目的,进而采取的技术方案如下:
一种厚板弯曲成形过程中中性层偏移量的计算方法,包括以下步骤:
步骤1:利用分层拉伸断裂实验,获得厚向强度线性分布曲线,根据厚板弯曲理论,建立应力-应变数学模型:
K--硬化系数,
n--硬化指数,
Q--反映厚向强度分布不均匀程度的系数,
z--板层距离中间层的垂直距离,
e--理论中性层偏移量,
t--厚板的厚度,
E--弹性模量;
步骤2:在拉伸区和压缩区分别取微元体进行受力平衡分析,并在极限状态下得到中性层偏移公式:
厚板弯曲的拉伸区的力学平衡关系式如下:
厚板弯曲的压缩区的力学平衡关系式如下:
σx--第一主应力,
β--中间主应力影响系数;
ρ--中性层的曲率半径,
α——微分单元体的长度,
ρya——微分单元体在压缩区的曲率半径,
ρla——微分单元体在拉伸区的曲率半径;
当中性层发生塑形变形时,拉伸区和压缩区的应力相等,即:
ρla=ρya (4)
根据公式(2)、(3)(4)可推导出理论中性层偏移量的计算公式:
ρ--中性层的曲率半径,
ρmin--最小弯曲曲率半径,
ρmax--最大弯曲曲率半径。
对于给定材料,工人只需测量弯曲半径、板厚即可利用公式(5)轻松计算出中性层偏移量,方便工人实际操作。
本发明的有益效果是:
1. 本发明提出一种确定中性层偏移量的方法,并给出计算公式。通过中性层偏移量的确定,可以协助制定更为准确的弯曲工艺参数。
2.本发明提出的确定中性层偏移量的方法可以提高厚壁直缝焊管的成形精度,降低生产成本。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为弯曲成形的模型的示意图;
图2为厚板弯曲的应力平衡关系图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分的实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
实施例1
以厚板的JCO成形为例:
如图1-2所示,一种厚板JCO成形过程中中性层偏移量的计算方法,包括以下步骤:
步骤1:利用分层拉伸断裂实验,获得厚向强度线性分布曲线,根据厚板弯曲理论,建立应力-应变数学模型:
K--硬化系数,
n--硬化指数,
Q--反映厚向强度分布不均匀程度的系数,
z--板层距离中间层的垂直距离,
e--理论中性层偏移量,
t--厚板的厚度,
E--弹性模量;
步骤2:在拉伸区和压缩区分别取微元体进行受力平衡分析,并在极限状态下得到中性层偏移公式:
厚板弯曲的拉伸区的力学平衡关系式如下:
厚板弯曲的压缩区的力学平衡关系式如下:
σx--第一主应力,
β--中间主应力影响系数;
ρ--中性层的曲率半径,
α——微分单元体的长度,
ρya——微分单元体在压缩区的曲率半径,
ρla——微分单元体在拉伸区的曲率半径;
当中性层发生塑形变形时,拉伸区和压缩区的应力相等,即:
ρla=ρya (4)
根据公式(2)、(3)(4)可推导出理论中性层偏移量的计算公式:
ρ--中性层的曲率半径,
ρmin--最小弯曲曲率半径,
ρmax--最大弯曲曲率半径。
实施例2
为了验证本发明,选用厚板材料为镍基合金为例。
厚板的基本参数为:厚板的宽度W=1000mm,厚板的长度L=2000mm,厚板的厚度t=60mm,杨氏模量E=205GPa,泊松比为0.3,密度为7.82 g/cm 3。上模曲率半径为150mm,下模曲率半径为150mm,下模跨距为770mm。上模和下模设置为刚体。其压下量设为30mm。
从模拟结果中的应变云图选取厚板中间两个单元体之间的结点,通过测量可以的得到:模型在压下量为30mm时拉伸区的弹性高度y 1为2.10mm,压缩区的弹性高度y 2为7.47mm,通过计算可以得到模拟中性层偏移量eˊ=4.79mm。
根据模拟结果可以得到板层距离中间层的垂直距离z为7.47mm,t=60mm,Q设为0.2,ρ=172.53mm,ρ min =150mm,ρ max =210mm,将其代入公式(5)中,可以得到理论中性层偏移量e=3.00mm。在模拟中的应变云图中相隔一个单元体选取两个单元体之间的结点,通过测量以及代入公式(5)中可以得到在厚板的厚度为60mm,压下量为30mm下的模拟中性层偏移数据和理论中性层偏移数据。通过计算模拟中性层偏移量与理论偏移量数据之间误差为1.62%。因此对于给定材料,工人只需测量弯曲半径、板厚即可利用本发明所述公式(5)轻松计算出中性层偏移量,方便工人实际操作。本发明提出的确定中性层偏移量的方法可以提高厚壁直缝焊管的成形精度,降低生产成本。
利用origin软件对所得到数据进行曲线拟合,采用多项式拟合,模拟的中性层偏移量的拟合函数为:
其中,R=0.97288(R为曲线相关系数)。
理论的中性层偏移量的拟合函数为:
其中,R=0.90515(R为曲线相关系数)。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进、部件拆分或组合等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种厚板弯曲成形过程中中性层偏移量的计算方法,其特征在于,包括以下步骤:
步骤1:利用分层拉伸断裂实验,获得厚向强度线性分布曲线,根据厚板弯曲理论,建立应力-应变数学模型:
K--硬化系数,
n--硬化指数,
Q--反映厚向强度分布不均匀程度的系数,
z--板层距离中间层的垂直距离,
e--理论中性层偏移量,
t--厚板的厚度,
E--弹性模量;
步骤2:在拉伸区和压缩区分别取微元体进行受力平衡分析,并在极限状态下得到中性层偏移公式:
厚板弯曲的拉伸区的力学平衡关系式如下:
厚板弯曲的压缩区的力学平衡关系式如下:
σ x --第一主应力,
β--中间主应力影响系数;
ρ--中性层的曲率半径,
α——微分单元体的长度,
ρ ya ——微分单元体在压缩区的曲率半径,
ρ la ——微分单元体在拉伸区的曲率半径;
当中性层发生塑形变形时,拉伸区和压缩区的应力相等,即:
ρ la= ρ ya (4)
根据公式(2)、(3)(4)可推导出理论中性层偏移量的计算公式:
ρ--中性层的曲率半径,
ρ min --最小弯曲曲率半径,
ρ max --最大弯曲曲率半径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111460909.1A CN114117679B (zh) | 2021-12-03 | 2021-12-03 | 一种厚板弯曲成形过程中中性层偏移量的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111460909.1A CN114117679B (zh) | 2021-12-03 | 2021-12-03 | 一种厚板弯曲成形过程中中性层偏移量的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114117679A true CN114117679A (zh) | 2022-03-01 |
CN114117679B CN114117679B (zh) | 2024-03-29 |
Family
ID=80365683
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111460909.1A Active CN114117679B (zh) | 2021-12-03 | 2021-12-03 | 一种厚板弯曲成形过程中中性层偏移量的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114117679B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012033039A (ja) * | 2010-07-30 | 2012-02-16 | Nippon Steel Corp | 材料の曲げ破断予測方法および装置、ならびにプログラムおよび記録媒体 |
CN103678932A (zh) * | 2013-12-24 | 2014-03-26 | 太原科技大学 | 用多极边界元法计算辊式矫直机矫直力的方法 |
CN104866641A (zh) * | 2014-12-10 | 2015-08-26 | 太原科技大学 | 一种棒材二辊矫直的回弹预测模型 |
CN109396226A (zh) * | 2018-12-28 | 2019-03-01 | 西北工业大学 | 一种管材弯曲中性层偏移位置确定方法 |
-
2021
- 2021-12-03 CN CN202111460909.1A patent/CN114117679B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012033039A (ja) * | 2010-07-30 | 2012-02-16 | Nippon Steel Corp | 材料の曲げ破断予測方法および装置、ならびにプログラムおよび記録媒体 |
CN103678932A (zh) * | 2013-12-24 | 2014-03-26 | 太原科技大学 | 用多极边界元法计算辊式矫直机矫直力的方法 |
CN104866641A (zh) * | 2014-12-10 | 2015-08-26 | 太原科技大学 | 一种棒材二辊矫直的回弹预测模型 |
CN109396226A (zh) * | 2018-12-28 | 2019-03-01 | 西北工业大学 | 一种管材弯曲中性层偏移位置确定方法 |
Non-Patent Citations (1)
Title |
---|
殷;梁荣;: "非对称截面梁纯弯曲中性层偏移理论解析", 燕山大学学报, no. 05, 30 September 2016 (2016-09-30) * |
Also Published As
Publication number | Publication date |
---|---|
CN114117679B (zh) | 2024-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112949108B (zh) | 热轧高强钢板形缺陷全流程预测方法及图形用户界面装置 | |
CN100357953C (zh) | 金属板带自动矫直工艺参数的优化方法 | |
CN105930559A (zh) | 一种确定管材数控弯曲回弹角的方法 | |
CN105956270B (zh) | 端部接触式少片端部加强型主副簧各片应力的计算方法 | |
Larrañaga et al. | Development and validation of a numerical model for sheet metal roll forming | |
JP6172401B2 (ja) | 金属板の圧延制御方法、圧延制御装置及び圧延金属板の製造方法 | |
CN118211349B (zh) | 一种考虑初始残余应力的板材拉伸弯曲矫直通用解析方法 | |
CN114117679A (zh) | 一种厚板弯曲成形过程中中性层偏移量的计算方法 | |
Sun et al. | Experimental investigation and prediction of the maximum edge longitudinal membrane strain and springback of Chain-die-formed AHSS U-channels using response surface methodology | |
CN110605531A (zh) | 双螺旋箱型钢结构旋转楼梯的加工方法 | |
CN110303070B (zh) | 一种双金属复合板矫直工艺新方法 | |
Görtan et al. | Roll forming of branched profiles | |
Elyasi et al. | Experimental investigation and numerical simulation of the effect of type of bending die on the quality of tube forming in rotary draw bending process | |
CN111842547B (zh) | 一种具有拉压不对称性金属板材矫直工艺新方法 | |
CN109877168A (zh) | 一种建立特厚板轧制力模型的方法 | |
CN105843989B (zh) | 非端部接触双加强型少片主副簧的副簧起作用载荷验算方法 | |
Liu et al. | Shape error prediction and compensation of three-dimensional surface in flexibly-reconfigurable roll forming | |
CN114386234B (zh) | 一种厚板jco成形工艺压下量的计算新方法 | |
Groche et al. | Benefits of stress superposition in combined bending-linear flow splitting process | |
Shim et al. | Improvement strategy for edge waviness in roll bending process of corrugated sheet metals | |
CN111842508A (zh) | 一种冷轧带材的轧后板形表示方法 | |
CN105956311A (zh) | 非端部接触式少片端部加强型副簧根部厚度的设计方法 | |
CN109977457A (zh) | 一种考虑温卷影响的加钒钢制筒节极限载荷预测方法 | |
CN111475896B (zh) | 一种棒材矫直应力中性层位置确定方法 | |
CN113664039B (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 |