CN114117679A - 一种厚板弯曲成形过程中中性层偏移量的计算方法 - Google Patents

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

Info

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
Application number
CN202111460909.1A
Other languages
English (en)
Other versions
CN114117679B (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.)
Taiyuan University of Science and Technology
Original Assignee
Taiyuan University of Science and Technology
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 Taiyuan University of Science and Technology filed Critical Taiyuan University of Science and Technology
Priority to CN202111460909.1A priority Critical patent/CN114117679B/zh
Publication of CN114117679A publication Critical patent/CN114117679A/zh
Application granted granted Critical
Publication of CN114117679B publication Critical patent/CN114117679B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force 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:利用分层拉伸断裂实验,获得厚向强度线性分布曲线,根据厚板弯曲理论,建立应力-应变数学模型:
Figure DEST_PATH_IMAGE001
(1)
K--硬化系数,
n--硬化指数,
Q--反映厚向强度分布不均匀程度的系数,
z--板层距离中间层的垂直距离,
e--理论中性层偏移量,
t--厚板的厚度,
E--弹性模量;
步骤2:在拉伸区和压缩区分别取微元体进行受力平衡分析,并在极限状态下得到中性层偏移公式:
厚板弯曲的拉伸区的力学平衡关系式如下:
Figure DEST_PATH_IMAGE002
(2)
厚板弯曲的压缩区的力学平衡关系式如下:
Figure DEST_PATH_IMAGE003
(3)
σx--第一主应力,
σy--第二主应力;其中:
Figure DEST_PATH_IMAGE004
β--中间主应力影响系数;
ρ--中性层的曲率半径,
α——微分单元体的长度,
ρya——微分单元体在压缩区的曲率半径,
ρla——微分单元体在拉伸区的曲率半径;
当中性层发生塑形变形时,拉伸区和压缩区的应力相等,即:
ρla=ρya (4)
根据公式(2)、(3)(4)可推导出理论中性层偏移量的计算公式:
Figure DEST_PATH_IMAGE005
(5)
ρ--中性层的曲率半径,
ρmin--最小弯曲曲率半径,
ρmax--最大弯曲曲率半径。
对于给定材料,工人只需测量弯曲半径、板厚即可利用公式(5)轻松计算出中性层偏移量,方便工人实际操作。
本发明的有益效果是:
1. 本发明提出一种确定中性层偏移量的方法,并给出计算公式。通过中性层偏移量的确定,可以协助制定更为准确的弯曲工艺参数。
2.本发明提出的确定中性层偏移量的方法可以提高厚壁直缝焊管的成形精度,降低生产成本。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为弯曲成形的模型的示意图;
图2为厚板弯曲的应力平衡关系图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分的实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
实施例1
以厚板的JCO成形为例:
如图1-2所示,一种厚板JCO成形过程中中性层偏移量的计算方法,包括以下步骤:
步骤1:利用分层拉伸断裂实验,获得厚向强度线性分布曲线,根据厚板弯曲理论,建立应力-应变数学模型:
Figure 813259DEST_PATH_IMAGE001
(1)
K--硬化系数,
n--硬化指数,
Q--反映厚向强度分布不均匀程度的系数,
z--板层距离中间层的垂直距离,
e--理论中性层偏移量,
t--厚板的厚度,
E--弹性模量;
步骤2:在拉伸区和压缩区分别取微元体进行受力平衡分析,并在极限状态下得到中性层偏移公式:
厚板弯曲的拉伸区的力学平衡关系式如下:
Figure 25673DEST_PATH_IMAGE002
(2)
厚板弯曲的压缩区的力学平衡关系式如下:
Figure 805410DEST_PATH_IMAGE003
(3)
σx--第一主应力,
σy--第二主应力;其中:
Figure 200619DEST_PATH_IMAGE004
β--中间主应力影响系数;
ρ--中性层的曲率半径,
α——微分单元体的长度,
ρya——微分单元体在压缩区的曲率半径,
ρla——微分单元体在拉伸区的曲率半径;
当中性层发生塑形变形时,拉伸区和压缩区的应力相等,即:
ρla=ρya (4)
根据公式(2)、(3)(4)可推导出理论中性层偏移量的计算公式:
Figure 749412DEST_PATH_IMAGE005
(5)
ρ--中性层的曲率半径,
ρ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,通过计算可以得到模拟中性层偏移量=4.79mm
根据模拟结果可以得到板层距离中间层的垂直距离z为7.47mmt=60mmQ设为0.2,ρ=172.53mm,ρ min =150mm,ρ max =210mm,将其代入公式(5)中,可以得到理论中性层偏移量e=3.00mm。在模拟中的应变云图中相隔一个单元体选取两个单元体之间的结点,通过测量以及代入公式(5)中可以得到在厚板的厚度为60mm,压下量为30mm下的模拟中性层偏移数据和理论中性层偏移数据。通过计算模拟中性层偏移量与理论偏移量数据之间误差为1.62%。因此对于给定材料,工人只需测量弯曲半径、板厚即可利用本发明所述公式(5)轻松计算出中性层偏移量,方便工人实际操作。本发明提出的确定中性层偏移量的方法可以提高厚壁直缝焊管的成形精度,降低生产成本。
利用origin软件对所得到数据进行曲线拟合,采用多项式拟合,模拟的中性层偏移量的拟合函数为:
Figure DEST_PATH_IMAGE006
(6)
其中,R=0.97288(R为曲线相关系数)。
理论的中性层偏移量的拟合函数为:
Figure DEST_PATH_IMAGE007
(7)
其中,R=0.90515(R为曲线相关系数)。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进、部件拆分或组合等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种厚板弯曲成形过程中中性层偏移量的计算方法,其特征在于,包括以下步骤:
步骤1:利用分层拉伸断裂实验,获得厚向强度线性分布曲线,根据厚板弯曲理论,建立应力-应变数学模型:
Figure 335797DEST_PATH_IMAGE001
(1)
K--硬化系数,
n--硬化指数,
Q--反映厚向强度分布不均匀程度的系数,
z--板层距离中间层的垂直距离,
e--理论中性层偏移量,
t--厚板的厚度,
E--弹性模量;
步骤2:在拉伸区和压缩区分别取微元体进行受力平衡分析,并在极限状态下得到中性层偏移公式:
厚板弯曲的拉伸区的力学平衡关系式如下:
Figure 580833DEST_PATH_IMAGE002
(2)
厚板弯曲的压缩区的力学平衡关系式如下:
Figure 891729DEST_PATH_IMAGE003
(3)
σ x --第一主应力,
σ y --第二主应力;其中:
Figure 552517DEST_PATH_IMAGE004
β--中间主应力影响系数;
ρ--中性层的曲率半径,
α——微分单元体的长度,
ρ ya ——微分单元体在压缩区的曲率半径,
ρ la ——微分单元体在拉伸区的曲率半径;
当中性层发生塑形变形时,拉伸区和压缩区的应力相等,即:
ρ la= ρ ya (4)
根据公式(2)、(3)(4)可推导出理论中性层偏移量的计算公式:
Figure 898048DEST_PATH_IMAGE005
(5)
ρ--中性层的曲率半径,
ρ min --最小弯曲曲率半径,
ρ max --最大弯曲曲率半径。
CN202111460909.1A 2021-12-03 2021-12-03 一种厚板弯曲成形过程中中性层偏移量的计算方法 Active CN114117679B (zh)

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)

* Cited by examiner, † Cited by third party
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 西北工业大学 一种管材弯曲中性层偏移位置确定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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