CN105373685B - 一种钢管连轧数值模拟过程高精度壁厚的确定方法 - Google Patents

一种钢管连轧数值模拟过程高精度壁厚的确定方法 Download PDF

Info

Publication number
CN105373685B
CN105373685B CN201510961177.2A CN201510961177A CN105373685B CN 105373685 B CN105373685 B CN 105373685B CN 201510961177 A CN201510961177 A CN 201510961177A CN 105373685 B CN105373685 B CN 105373685B
Authority
CN
China
Prior art keywords
steel pipe
wall thickness
node
interpolation
cubic spline
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
Application number
CN201510961177.2A
Other languages
English (en)
Other versions
CN105373685A (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.)
Northeastern University China
Original Assignee
Northeastern University China
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 Northeastern University China filed Critical Northeastern University China
Priority to CN201510961177.2A priority Critical patent/CN105373685B/zh
Publication of CN105373685A publication Critical patent/CN105373685A/zh
Application granted granted Critical
Publication of CN105373685B publication Critical patent/CN105373685B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Control Of Metal Rolling (AREA)

Abstract

本发明一种钢管连轧数值模拟过程高精度壁厚的确定方法，属于轧制技术领域，本发明基于三次样条插值方法对变形后钢管内外表面节点坐标进行插值，获得内外表面平滑的插值曲线，再求得与计算壁厚处的内表面节点处切线相垂直的直线方程，求其与外表面插值曲线的交点，进而计算钢管变形后的壁厚值；由于三次样条曲线具有在插值点上的一阶、二阶导数连续，具有连续的、曲率变化均匀的特点，可以在最大程度上保证了钢管壁厚的计算精度；这种方法对壁厚的求解精度远高于直接使用节点坐标计算或对节点坐标使用线性插值、抛物线插值等计算壁厚的方法，对于无缝钢管连轧过程的工艺调整和参数优化具有重要的实际意义。

Description

一种钢管连轧数值模拟过程高精度壁厚的确定方法
技术领域
本发明属于轧制技术领域，具体涉及一种钢管连轧数值模拟过程高精度壁厚的确定方法。
背景技术
无缝钢管连轧过程的金属流动规律异常复杂，使用有限元软件对其连轧过程进行数值模拟是常见的孔型优化和缺陷分析的手段；对钢管连轧过程使用有限元软件模拟得到最终结果后，初始划分的、排列规整的单元和节点坐标随着钢管的变形也逐渐发生变化，作为评价钢管连轧孔型参数设计是否合理的重要参数，钢管壁厚值随着单元的变形难以直接从模拟结果中获得精确值，如图1中图(a)和图(b)所示，变形前钢管内外表面对应节点间距离为初始钢管壁厚，变形后钢管内外表面节点位置均发生了变化，对应节点间距离并不能代表壁厚。
在得到钢管成形数值模拟结果后，常见的壁厚的求解方法是在内外壁上取两个对应的节点，然后直接求它们之间的距离，这样求解出来的结果无法保证精度，因为钢管三维轧制是大塑性变形过程，金属流动性大，变形前所取的代表壁厚的两个节点在变形之后相对位置已经发生了变化，即使提高划分单元的数量，用这种方式计算的钢管壁厚也是不准确的。
发明内容
针对现有技术的不足，本发明提出一种钢管连轧数值模拟过程高精度壁厚的确定方法，以达到提高钢管壁厚计算精度的目的。
一种钢管连轧数值模拟过程高精度壁厚的确定方法，包括以下步骤：
步骤1、采用有限元分析方法对钢管连轧过程进行数值模拟；
步骤2、根据钢管变形后所要计算壁厚的位置点，提取钢管横断面该位置点附近内外表面节点坐标；
步骤3、基于所提取的内外表面节点坐标值，采用三次样条插值方法重建钢管内外表面曲线；
步骤4、通过计算壁厚的位置点，做钢管横断面内表面样条曲线的切线；
步骤5、通过计算壁厚的位置点，做所述切线的垂线，获得该垂线与钢管横断面外表面样条曲线的交点；
步骤6、获得上述交点与计算壁厚的位置点之间的距离，所述的距离即为钢管壁厚。
步骤2所述的提取钢管横断面该位置点附近内外表面节点坐标，所提取的节点的个数根据实际需求而定。
本发明优点：
本发明提出一种钢管连轧数值模拟过程高精度壁厚的确定方法，基于三次样条插值方法对变形后钢管内外表面节点坐标进行插值，获得内外表面平滑的插值曲线，再求得与计算壁厚处的内表面节点处切线相垂直的直线方程，求其与外表面插值曲线的交点，进而计算钢管变形后的壁厚值；由于三次样条曲线具有在插值点上的一阶、二阶导数连续，具有连续的、曲率变化均匀的特点，可以在最大程度上保证了钢管壁厚的计算精度；本发明采用的方法可求出任意变形时刻钢管变形后的高精度壁厚分布曲线，由于三次样条曲线具有连续、曲率变化均匀的特点，这种方法对壁厚的求解精度远高于直接使用节点坐标计算或对节点坐标使用线性插值、抛物线插值等计算壁厚的方法，保证了根据有限元软件数值模拟结果对孔型参数优化评价的准确性，对于无缝钢管连轧过程的工艺调整和参数优化具有重要的实际意义。
附图说明
图1为本发明一种实施例的钢管横断面上变形前后内外表面节点坐标示意图，其中，图(a)为变形前内外表面节点坐标示意图，图(b)为变形后内外表面节点坐标示意图；1-内壁节点，2-外壁节点；
图2为本发明一种实施例的钢管连轧数值模拟过程高精度壁厚的确定方法流程图；
图3为本发明一种实施例的钢管内表面节点3处壁厚求解示意图。
具体实施方式
下面结合附图对本发明一种实施例做进一步说明。
本发明实施例中，钢管连轧数值模拟过程高精度壁厚的确定方法，方法流程图如图2所示，包括以下步骤：
步骤1、采用有限元分析方法对钢管连轧过程进行数值模拟；
本发明实施例中，壁厚的计算数据来源于采用有限元软件对某5机架钢管连轧过程进行数值模拟后的结果，末机架的壁厚代表整个连轧机组轧制结束厚度，选择第5机架的模拟结果中的内外表面节点实施计算；
本发明实施例中，轧制产品的工艺参数如下：
●钢种：Q235
●连轧入口温度：1050℃
●芯棒直径：181.2mm
●连轧入口荒管尺寸：218mm×14.38mm
●连轧出口荒管尺寸：190mm×5.5mm
步骤2、根据钢管变形后所要计算壁厚的位置点，提取钢管横断面该位置点附近内外表面节点坐标；
本发明实施例中，如图3所示，在钢管连轧数值模拟结果中，要想求出钢管内表面节点3处的壁厚，首先以内表面节点3为中心，向两个相邻方向对称提取节点坐标，本实施例在内外表面各提取了5个节点。坐标值如下：
表1提取的节点坐标集合
步骤3、基于所提取的内外表面节点坐标值，采用三次样条插值方法重建钢管内外表面曲线；
本发明实施例中，对内、外表面的节点坐标分别进行三次样条插值，插值的附加的边界条件取自然边界条件，即两端的二次导数为0，建立三次样条曲线参数求解方程组，通过追赶法求解方程组获得内、外表面共两条三次样条插值曲线的描述参数。
本发明实施例中，样条插值的节点数共5个，编号为j＝0，1，...，4，则在每个区间[xj，xj+1]内的三次样条函数表达形式S(x)为：
其中，xj表示第j个节点的横坐标，xj+1表示第j+1个节点的横坐标，yj表示第j个节点的纵坐标，yj+1表示第j+1个节点的纵坐标，hj＝xj+1-xj，Mj和Mj+1表示三次样条曲线参数，是需要求解的未知数；
根据三次样条函数特性：
S′(xj+0)＝S′(xj-0) (2)
得到如下方程组：
μjMj-1+2Mj+λjMj+1＝dj (3)
其中，
在节点集合中的两个端点处选择为自然边界条件的情况：Mo＝f″0＝0，Mn＝f″n＝0，这里n＝4，方程组(3)表示为如下矩阵形式：
此矩阵严格对角占优，采用追赶法可求解出三次样条插值方程系数Mi，即得到了钢管内、外表面的三次样条插值平滑曲线S1(x)、S2(x)；
本发明实施例中，将内外节点坐标值分别带入公式(3)中，得到未知数的矩阵表示形式，通过追赶法求解三次样条曲线参数Mj，编号为j＝0，1，...，4，最终得到内外表面的分段插值样条函数，由于要求解内表面节点3位置处的壁厚，本发明实施例中列出内外表面节点2、3和节点3、4处的插值函数：
内表面节点2与节点3之间的三次样条插值函数S1_23(x)为：
S1-23〔x)＝23.6×〔90.85-x)3-21.6×(90.85-x)-6.608×(x-91.5)3-16.67×(x-91.5)
内表面节点3与节点4之间的三次样条插值函数S1_34(x)为：
S1_34(x)＝-7.58×(90.29-x)3-19.87×(90.29-x)-30.33×(x-90.85)
外表面节点2与节点3之间的三次样条插值函数S2_23(x)为：
S2-23(x)＝103.9×(96.51-x)3-34.12×〔96.51-x)+1.7×〔x-96.8)3-46.5×〔x-96.8)
外表面节点3与节点4之间的三次样条插值函数S2_34(x)为：
S2_34〔x)＝0.97×〔96.0-x)3-26.78×〔96.0-x)-38.0×〔x-96.51)
步骤4、通过计算壁厚的位置点，做钢管横断面内表面样条曲线的切线；
本发明实施例中，为了能够精确的计算钢管壁厚，首先以钢管内表面节点为基础，在内表面三次样条插值曲线S1(x)上，求解钢管内表面需要计算壁厚的节点3位置处与三次样条插值曲线S1-23(x)(或S1-34(x))相切的直线方程Y1i；
Y1i＝-3.3x+312.3
步骤5、通过计算壁厚的位置点，做所述切线的垂线，获得该垂线与钢管横断面外表面样条曲线的交点；
本发明实施例中，求解通过内表面节点3，并且与直线方程Y1i相垂直的直线方程Y2i；
Y2i＝0.304x-15.1
本发明实施例中，直线方程Y2i与外表面样条插值曲线S2(x)的交点可能落在外表面节点2、3之间，或者3、4之间，需要根据实际情况单独判断；采用牛顿迭代方法求解直线方程Y2i与钢管外表面节点三次样条插值曲线S2(x)的交点，得到的结果坐标(x，y)如下：
x＝96.345，y＝14.195
步骤6、获得上述交点与计算壁厚的位置点之间的距离，所述的距离即为钢管壁厚。
本发明实施例中，内表面节点3的坐标位置为已知(x3，y3)，与计算得到的交点(x，y)之间的距离即为所要求解的壁厚Di：
按照以上方法，可以求解钢管塑性变形后任意节点处的壁厚。

Claims (2)

1.一种钢管连轧数值模拟过程高精度壁厚的确定方法，其特征在于，包括以下步骤：
步骤1、采用有限元分析方法对钢管连轧过程进行数值模拟；
步骤2、根据钢管变形后所要计算壁厚的位置点，提取钢管横断面该位置点附近内外表面节点坐标；
步骤3、基于所提取的内外表面节点坐标值，采用三次样条插值方法重建钢管内外表面曲线；
所述对内、外表面的节点坐标分别进行三次样条插值，插值的附加的边界条件取自然边界条件，即两端的二次导数为0，建立三次样条曲线参数求解方程组，通过追赶法求解方程组获得内、外表面共两条三次样条插值曲线的描述参数；
步骤4、通过计算壁厚的位置点，做钢管横断面内表面样条曲线的切线；
步骤5、通过计算壁厚的位置点，做所述切线的垂线，获得该垂线与钢管横断面外表面样条曲线的交点；
步骤6、获得上述交点与计算壁厚的位置点之间的距离，所述的距离即为钢管壁厚。
2.根据权利要求1所述的钢管连轧数值模拟过程高精度壁厚的确定方法，其特征在于，步骤2所述的提取钢管横断面该位置点附近内外表面节点坐标，所提取的节点的个数根据实际需求而定。
CN201510961177.2A 2015-12-17 2015-12-17 一种钢管连轧数值模拟过程高精度壁厚的确定方法 Active CN105373685B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510961177.2A CN105373685B (zh) 2015-12-17 2015-12-17 一种钢管连轧数值模拟过程高精度壁厚的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510961177.2A CN105373685B (zh) 2015-12-17 2015-12-17 一种钢管连轧数值模拟过程高精度壁厚的确定方法

Publications (2)

Publication Number Publication Date
CN105373685A CN105373685A (zh) 2016-03-02
CN105373685B true CN105373685B (zh) 2018-05-04

Family

ID=55375880

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510961177.2A Active CN105373685B (zh) 2015-12-17 2015-12-17 一种钢管连轧数值模拟过程高精度壁厚的确定方法

Country Status (1)

Country Link
CN (1) CN105373685B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110852007B (zh) * 2019-10-28 2021-05-14 北京科技大学 一种考虑不均匀变形热的大方坯轧制温度场计算方法
CN111215592B (zh) * 2020-02-21 2020-11-17 东北大学 一种重轨钢连铸坯偏析缺陷遗传性研究方法及其应用

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101097439A (zh) * 2006-06-30 2008-01-02 宝山钢铁股份有限公司 张减机壁厚前馈结合反馈控制方法
CN101347795A (zh) * 2008-09-03 2009-01-21 天津钢管集团股份有限公司 连轧机更换轧制规格的调整方法
CN101826117A (zh) * 2009-03-04 2010-09-08 中国核电工程有限公司 一种管道系统有限单元法力学计算模型制作方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101097439A (zh) * 2006-06-30 2008-01-02 宝山钢铁股份有限公司 张减机壁厚前馈结合反馈控制方法
CN101347795A (zh) * 2008-09-03 2009-01-21 天津钢管集团股份有限公司 连轧机更换轧制规格的调整方法
CN101826117A (zh) * 2009-03-04 2010-09-08 中国核电工程有限公司 一种管道系统有限单元法力学计算模型制作方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周期变厚度带材轧制控制系统开发研究;吴志强;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20150715;第B022-161页 *

Also Published As

Publication number Publication date
CN105373685A (zh) 2016-03-02

Similar Documents

Publication Publication Date Title
CN103440382B (zh) 一种基于逆向反求技术的3d打印塑件产品的方法
CN105373685B (zh) 一种钢管连轧数值模拟过程高精度壁厚的确定方法
CN101672637A (zh) 一种复杂曲面的数字化检测方法
CN108682043A (zh) 一种基于参数映射的复杂曲面测量规划方法
CN109013717B (zh) 一种热连轧中间坯心部温度计算方法
CN103729505A (zh) 一种基于cfd的阀门当量长度计算方法
CN103035029A (zh) 用离散裂缝端点变形消除油藏数值模型极小化网格的方法
CN111046601A (zh) 基于工件变形代理模型的薄壁曲面加工进给速度优化方法
CN103810326A (zh) 一种鱼池流态数字化分析与构建方法
CN103464475B (zh) 一种基于关联神经网络的热轧卷取温度预报方法
CN201514207U (zh) 一种复杂曲面的数字化检测系统
CN101949680B (zh) 一种型面测量用叶片标准件的设计方法
CN113935209A (zh) 一种基于fluent的辊式连续淬火过程温度场仿真方法
CN106202686B (zh) 一种涡轮盘等温模锻预成形坯料的多目标设计方法
CN105855297B (zh) 一种提高热轧首块无取向硅钢头部厚度精度的控制方法
CN114417674A (zh) 一种圆角方坯连铸凝固传热的有限元计算方法
Jiang et al. Method for improving transverse wall thickness precision of seamless steel tube based on tube rotation
CN103484642B (zh) 一种研究制管变形对管线钢氢致开裂性能影响的方法
CN107030121B (zh) 一种连铸坯感应加热快速自适应温控方法
CN102722661B (zh) 一种高速列车通过长大隧道时隧道出口微压波的预测方法
CN116147724B (zh) 一种适用于超声水表的计量方法
CN115659875B (zh) 一种棒束燃料组件试验装置及其试验参数获取方法
CN105855298B (zh) 一种高精度含硼钢热轧轧制力计算方法
CN103884306A (zh) 一种研究大口径热挤压成型三通壁厚的测试方法
CN114491855A (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
GR01 Patent grant
GR01 Patent grant