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
offset
curvature
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

    • 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：利用分层拉伸断裂实验，获得厚向强度线性分布曲线，根据厚板弯曲理论，建立应力-应变数学模型：
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，通过计算可以得到模拟中性层偏移量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软件对所得到数据进行曲线拟合，采用多项式拟合，模拟的中性层偏移量的拟合函数为：
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
CN102313677B (zh) 一种识别焊缝不同区域材料动态力学性能参数的检测方法
CN112949108B (zh) 热轧高强钢板形缺陷全流程预测方法及图形用户界面装置
CN116078831B (zh) 一种板带冷轧动态减薄生产过程中轧制力的计算方法
CN118211349B (zh) 一种考虑初始残余应力的板材拉伸弯曲矫直通用解析方法
JP6172401B2 (ja) 金属板の圧延制御方法、圧延制御装置及び圧延金属板の製造方法
CN111487869B (zh) 一种矫直控制方法及装置
CN105956270A (zh) 端部接触式少片端部加强型主副簧各片应力的计算方法
Abathun et al. Effects of manufacturing methods and production routes on residual stresses of rectangular and square hollow steel sections: a review
CN115186523A (zh) 一种基于仿真的高强钢焊后热影响区的本构模型建立方法
CN119416590B (zh) 一种涉及力学多参数变化的带钢拉伸弯曲矫直计算方法
Larrañaga et al. Development and validation of a numerical model for sheet metal roll forming
CN114117679A (zh) 一种厚板弯曲成形过程中中性层偏移量的计算方法
Liu et al. Shape change simulation analysis of wheel steel in a four-high hot rolling mill
CN112207135A (zh) 一种冷轧复合双金属板平直控制的方法
CN116371942A (zh) 一种基于带钢横向强度不均的横向厚度分布预测方法
CN115017450A (zh) 一种厚壁管材制备过程中回弹量的计算方法
CN105840702A (zh) 非端部接触式少片抛物线型变截面副簧长度的设计方法
CN110617933B (zh) 一种双对称截面开口薄壁梁冲击载荷测量方法
CN111842547B (zh) 一种具有拉压不对称性金属板材矫直工艺新方法
CN105414179A (zh) 一种改善铁素体不锈钢板带组织性能的冷轧方法及其分析方法
CN112926146B (zh) 一种针对冷轧轧辊的在线疲劳预测方法及系统
CN110688787B (zh) 局部受横向弯矩作用的长条形基座的螺栓群受力确定方法
CN111241640B (zh) 一种计算棒材弯曲过程中径向截面中性层偏移量的方法
CN114386234B (zh) 一种厚板jco成形工艺压下量的计算新方法
CN109977457A (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