CN106964654A - 一种求解轧制力的建模和仿真方法 - Google Patents
一种求解轧制力的建模和仿真方法 Download PDFInfo
- Publication number
- CN106964654A CN106964654A CN201610977304.2A CN201610977304A CN106964654A CN 106964654 A CN106964654 A CN 106964654A CN 201610977304 A CN201610977304 A CN 201610977304A CN 106964654 A CN106964654 A CN 106964654A
- Authority
- CN
- China
- Prior art keywords
- modeling
- zone
- slippage
- force
- simulation
- 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.)
- Pending
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21B—ROLLING OF METAL
- B21B37/00—Control devices or methods specially adapted for metal-rolling mills or the work produced thereby
- B21B37/58—Roll-force control; Roll-gap control
-
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mechanical Engineering (AREA)
- Control Of Metal Rolling (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种求解轧制力的建模和仿真方法,利用Simulink建立联合仿真总***来求解轧制力计算方程——卡尔曼微分方程,利用Simulink建立的联合仿真总***包括五个子***:后滑区子***、前滑区子***、停止子***、动态变形抗力子***和最大接触角子***。本发明采用Simulink强大的运算功能实现了轧制塑性区变形的求解,通过建模、仿真一目了然的可以看出单位轧制力在变形区的变化,也从另个侧面反映了中性角的存在,并可以计算其值βF。
Description
技术领域
本发明涉轧钢领域中轧制力的计算,具体涉及一种求解轧制力的建模和仿真方法。
背景技术
平辊轧制时轧制力的理论计算,从1925年卡尔曼建立微分平衡方程以来,已经有九十多年历史。现有的轧制力计算公式,大部分是在卡尔曼方程基础上,加上各种假设条件解出来的。而且为了更加容易求解轧制力,变形抗力K一般根据钢种取固定值(K=1.15σs,其中σs:钢的屈服强度),但往往由于轧制过程中,由于钢的屈服强度是变化的,因此会大大影响轧制力的计算精度。
采用卡尔曼微分方程计算轧制力,由于该方程是一个带有非常量系数的普通一次线性微分方程。由于前滑和后滑的符号变化我们没有一个解析解,这一过程反复重复直到前滑应力与后滑应力相等。此时的β为中性角βF,βF所对应的点为中性点F。
卡尔曼微分方程的建立过程:卡尔曼微分方程及其数值积分推导所用假设条件为:
一、卡尔曼微分方程及其数值积分推导所用假设条件为:
1)变形区内沿轧件横断面高度上的各点的金属流动速度、应力及变形均匀分布;
2)当(或)的比值很大时,认为宽展很小,可以忽略,Δb=0。这样把三个方向都有变形的空间问题变成了只有两个方向变形的平面变形问题;
3)轧制时,轧件高向、纵向和横向的变形都与主应力方向一致,忽略了切应力的影响;
4)认为金属质点在变形过程中,性质处处相同;
5)上下轧辊辊径相等,并做匀速运动,不产生惯性力,轧辊和机架为刚体,即不产生弹性变形。
二、卡尔曼单位压力微分方程推导:
在变形区内的后滑区任取一微分体积abcd,其厚度为dx,如图1所示。高度在微分体的右侧为hx,左侧为hx+dhx;弧长近似为弦长,弦长轧件宽度为单位1,如说明书附图中图1所示。
在接触表面上轧辊对轧件的作用力有径向单位压力p及单位摩擦力t,在后滑区,接触面上金属质点向着轧辊转动相反方向滑动,它们在接触弧ab上的合力水平投影为:
式中“+”代表后滑区摩擦力水平方向水平投影与轧制方向一致,“-”代表前滑区摩擦力水平方向水平投影与轧制方向相反。θ—ab弧切线与水平面所夹成的夹角,亦即相对应的圆心角。
根据假设条件1),则作用在微分体积两侧的应力各为σx及σx+dσx,设2y=hx,则其合力为:σx2y和(σx+dσx)2(y+dy)。因为微分体积abcd处于平衡,故作用在微分体上的水平力和为零(∑x=0),即平衡微分方程为:
Sx=2Bsxy-2B(sx+dsx)(y+dy)-2ptanqdxB±2tdxB=0
原设没有宽展,并取tanθ=dy/dx,忽略高阶项,对上式进行简化,可以得到:
为了对方程求解,必须找出单位压力p与水平压应力σx之间的关系。根据假设,设水平压应力σx和垂直压应力σy为主应力,则可得出:
忽略第二项,则同时s1=-sx。带入屈服条件式则:-sx-(-p)=K,或p-sx=K。式中K平面变形抗力,K=1.15σs。由以上条件可以得到卡尔曼单位压力微分方程的一般式:
轧制过程中采用乳液润滑因而产生滑动摩擦,摩擦系数μ取值在0.01到0.1之间,所以库仑摩擦定律是有效的,即:t(θ)=μp(θ)。
本发明所述方案,采用Simulink强大的运算功能实现了轧制塑性区变形的求解,通过建模、仿真一目了然的可以看出单位轧制力在变形区的变化,也从另个侧面反映了中性角的存在,并可以计算其值βF。
发明内容
为实现上述目的,本发明采用的技术方案如下:
一种求解轧制力的建模和仿真方法,所述求解轧制力的建模和仿真方法是利用Simulink建立联合仿真总***求解轧制力计算方程——卡尔曼微分方程,所述利用Simulink建立的联合仿真总***包括五个子***:后滑区子***、前滑区子***、停止子***、动态变形抗力子***和最大接触角子***。
所述求解轧制力计算方程——卡尔曼微分方程中,仿真函数是关于单位轧制压力和接触角θ的函数,方程中取dx=Rcosqdq,dy=Rsinqdq则卡尔曼方程变为关于变量θ的微分方程:
对于所述变形抗力kf***来说,在冷连轧过程中,带钢变形抗力是决定模型计算精度关键因素,冷连轧机一般采用变形抗力kf与累计厚度压缩比εΣ的关系建立如下公式:
为断面收缩率为零时屈服极限、为加工硬化曲线指数、εΣ为断面收缩率即忽略宽展时的累计厚度压缩比,这里:εΣ=[he-h(β)]/he。
进一步的,建模仿真后,单位轧制力最小值对应的角度为中性角
所述求解轧制力的建模和仿真方法是利用Simulink对前滑区***和后滑区***进行建模,并将前滑区***和后滑区***计算结果用加法模块求和进行联合建模。
前滑区与后滑区建模不同在于摩擦系数μ=artg(sita)转换为弧度时乘上了值为-1的增益模块;仿真时以角度β作为时间变量,即横坐标,单位采用弧度制。两次建模采用向同一方向进行积分,避开了前滑区计算由最大β角向中性角βF积分;利用两次计算时摩擦力的相反方向可以同为从0弧度向最大弧度的积分。
采用Simulink强大的运算功能实现了轧制塑性区变形的求解,通过建模、仿真一目了然的可以看出单位轧制力在变形区的变化,也从另个侧面反映了中性角的存在,并可以计算其值βF。
附图说明
图1为卡尔曼模型示意图;
图2为变形抗力kf***图;
图3为后滑区轧制过程基本方程的建模图;
图4为前滑区轧制过程基本方程的建模图;
图5为后滑区仿真结果图;
图6为前滑区仿真结果图;
图7为联合建模图;
图8为联合仿真结果图。
具体实施方式
为了更好地理解本发明,下面附图与实施例进一步清楚阐述本发明的内容,但本发明的保护内容不仅仅局限于下面的实施例,具体实施方式如下:
一种求解轧制力的建模和仿真方法,所述求解轧制力的建模和仿真方法是利用Simulink建立联合仿真总***求解轧制力计算方程——卡尔曼微分方程,所述利用Simulink建立的联合仿真总***包括三个主要子***:后滑区子***、前滑区子***和停止子***以及两个动态变形抗力子***和最大接触角子***。
所述求解轧制力计算方程——卡尔曼微分方程中,仿真函数是关于单位轧制压力和接触角θ的函数,方程中取dx=R cos q d q, dy=R sin q d q则卡尔曼方程变为关于变量θ的微分方程:
变形抗力kf***见图2在冷连轧过程中,带钢变形抗力是决定模型计算精度关键因素,冷连轧机一般采用变形抗力kf与累计厚度压缩比εΣ的关系建立如下公式:
为断面收缩率为零时屈服极限、为加工硬化曲线指数、εΣ为断面收缩率即忽略宽展时的累计厚度压缩比,这里:εΣ=[he-h(β)]/he。
后滑区轧制过程基本方程的建模见图3。
前滑区轧制过程基本方程的建模见图4。
仿真结果
前滑区与后滑区建模不同在于摩擦系数μ=artg(sita)转换为弧度时乘上了值为-1的增益模块;仿真时以角度β作为时间变量,即横坐标,单位采用弧度制。两次建模采用向同一方向进行积分,避开了前滑区计算由最大β角向中性角βF积分;利用两次计算时摩擦力的相反方向可以同为从0弧度向最大弧度的积分。如图5、图6所示。
Simulink建立联合仿真总***,将前滑区***(forward_slip)和后滑区***(backward_slip)计算结果用加法模块求和建模,如图7。仿真结果类似抛物线如图8显示,右侧单位轧制力为零点处的横坐标即为中性角βF的值。
这里采用联合仿真的数据是某冷轧厂宽度为1米冷轧板的生产实际数据,如表1所示。
表1 生产数据
Claims (6)
1.一种求解轧制力的建模和仿真方法,其特征在于:所述求解轧制力的建模和仿真方法是利用Simulink建立联合仿真总***来求解轧制力计算方程——卡尔曼微分方程,所述利用Simulink建立的联合仿真总***包括五个子***:后滑区子***、前滑区子***、停止子***、动态变形抗力子***和最大接触角子***。
2.根据权利要求1所述的一种求解轧制力的建模和仿真方法,其特征在于:所述求解轧制力计算方程——卡尔曼微分方程中,仿真函数是关于单位轧制压力和接触角θ的函数,方程中取dx=Rcosqdq,dy=Rsinqdq则卡尔曼方程变为关于变量θ的微分方程:
3.根据权利要求1所述的一种求解轧制力的建模和仿真方法,其特征在于:对于所述变形抗力***来说,在冷连轧过程中,带钢变形抗力是决定模型计算精度关键因素,冷连轧机一般采用变形抗力与累计厚度压缩比的关系建立如下公式:
为断面收缩率为零时屈服极限、为加工硬化曲线指数、εΣ为断面收缩率即忽略宽展时的累计厚度压缩比,这里:εΣ=[he-h(β)]/he。
4.根据权利要求2所述的一种求解轧制力的建模和仿真方法,其特征在于:建模仿真后单位轧制力最小值对应的角度为中性角。
5.根据权利要求1所述的一种求解轧制力的建模和仿真方法,其特征在于:所述求解轧制力的建模和仿真方法是利用Simulink对前滑区***和后滑区***进行建模,并将前滑区***和后滑区***计算结果用加法模块求和,进行联合建模。
6.根据权利要求5所述的一种求解轧制力的建模和仿真方法,其特征在于:前滑区与后滑区建模不同在于摩擦系数μ=artg(sita)转换为弧度时,乘上了值为-1的增益模块;仿真时,以角度β作为时间变量,即横坐标,单位采用弧度制;两次建模采用向同一方向进行积分,避开了前滑区计算由最大β角向中性角βF积分;利用两次计算时摩擦力的相反方向可以同为从0弧度向最大弧度的积分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610977304.2A CN106964654A (zh) | 2016-11-07 | 2016-11-07 | 一种求解轧制力的建模和仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610977304.2A CN106964654A (zh) | 2016-11-07 | 2016-11-07 | 一种求解轧制力的建模和仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106964654A true CN106964654A (zh) | 2017-07-21 |
Family
ID=59334349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610977304.2A Pending CN106964654A (zh) | 2016-11-07 | 2016-11-07 | 一种求解轧制力的建模和仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106964654A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107983772A (zh) * | 2017-11-27 | 2018-05-04 | 河北工业大学 | 一种板材冷轧过程中可靠性有限元模型的建立方法 |
CN108655179A (zh) * | 2018-05-31 | 2018-10-16 | 太原晋西春雷铜业有限公司 | 金属板带材轧制前滑控制方法 |
CN109359429A (zh) * | 2018-11-29 | 2019-02-19 | 中冶南方工程技术有限公司 | 一种板带轧制前滑值计算方法及装置 |
CN109711101A (zh) * | 2019-01-24 | 2019-05-03 | 太原科技大学 | 一种求解厚钢板同速异径蛇形轧制力能参数的方法 |
CN109840373A (zh) * | 2019-01-24 | 2019-06-04 | 太原科技大学 | 一种求解厚钢板同径异速蛇形轧制力能参数的方法 |
CN109885972A (zh) * | 2019-03-22 | 2019-06-14 | 东北大学 | 异步冷轧板带的变形区划分及各分区占比的计算方法 |
CN109926456A (zh) * | 2019-03-25 | 2019-06-25 | 燕山大学 | 一种混合润滑状态下的轧制力预报方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030055604A (ko) * | 2001-12-27 | 2003-07-04 | 주식회사 포스코 | 스테인레스강 열간압연시 압연하중 예측방법 |
CN101604144A (zh) * | 2008-06-13 | 2009-12-16 | 中国科学院金属研究所 | 一种板材轧制在线控制模型的建模方法 |
CN101869914A (zh) * | 2009-04-22 | 2010-10-27 | 宝山钢铁股份有限公司 | 一种精轧带钢的厚度控制方法及装置 |
CN102294364A (zh) * | 2010-06-22 | 2011-12-28 | 宝山钢铁股份有限公司 | 极薄板平整机的轧制力预设定方法 |
CN103447308A (zh) * | 2013-07-02 | 2013-12-18 | 西安理工大学 | 一种轧机半实物仿真***及仿真方法 |
CN104624669A (zh) * | 2015-02-10 | 2015-05-20 | 北京科技大学 | 一种电火花毛化工作辊生产带钢的粗糙度预测和控制方法 |
-
2016
- 2016-11-07 CN CN201610977304.2A patent/CN106964654A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20030055604A (ko) * | 2001-12-27 | 2003-07-04 | 주식회사 포스코 | 스테인레스강 열간압연시 압연하중 예측방법 |
CN101604144A (zh) * | 2008-06-13 | 2009-12-16 | 中国科学院金属研究所 | 一种板材轧制在线控制模型的建模方法 |
CN101869914A (zh) * | 2009-04-22 | 2010-10-27 | 宝山钢铁股份有限公司 | 一种精轧带钢的厚度控制方法及装置 |
CN102294364A (zh) * | 2010-06-22 | 2011-12-28 | 宝山钢铁股份有限公司 | 极薄板平整机的轧制力预设定方法 |
CN103447308A (zh) * | 2013-07-02 | 2013-12-18 | 西安理工大学 | 一种轧机半实物仿真***及仿真方法 |
CN104624669A (zh) * | 2015-02-10 | 2015-05-20 | 北京科技大学 | 一种电火花毛化工作辊生产带钢的粗糙度预测和控制方法 |
Non-Patent Citations (4)
Title |
---|
张贞贞等: "铝热轧机二级自动化***中轧制力模型的探究", 《物理测试》 * |
王平等: "《金属塑性成形力学》", 31 August 2006, 冶金工业出版社 * |
赵志业: "《金属塑性加工力学》", 30 November 1987, 冶金工业出版社 * |
邹家祥: "《轧钢机械》", 31 May 1995, 冶金工业出版社 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107983772A (zh) * | 2017-11-27 | 2018-05-04 | 河北工业大学 | 一种板材冷轧过程中可靠性有限元模型的建立方法 |
CN107983772B (zh) * | 2017-11-27 | 2019-06-04 | 河北工业大学 | 一种板材冷轧过程中可靠性有限元模型的建立方法 |
CN108655179A (zh) * | 2018-05-31 | 2018-10-16 | 太原晋西春雷铜业有限公司 | 金属板带材轧制前滑控制方法 |
CN108655179B (zh) * | 2018-05-31 | 2019-08-23 | 太原晋西春雷铜业有限公司 | 金属板带材轧制前滑控制方法 |
CN109359429A (zh) * | 2018-11-29 | 2019-02-19 | 中冶南方工程技术有限公司 | 一种板带轧制前滑值计算方法及装置 |
CN109711101A (zh) * | 2019-01-24 | 2019-05-03 | 太原科技大学 | 一种求解厚钢板同速异径蛇形轧制力能参数的方法 |
CN109840373A (zh) * | 2019-01-24 | 2019-06-04 | 太原科技大学 | 一种求解厚钢板同径异速蛇形轧制力能参数的方法 |
CN109885972A (zh) * | 2019-03-22 | 2019-06-14 | 东北大学 | 异步冷轧板带的变形区划分及各分区占比的计算方法 |
CN109926456A (zh) * | 2019-03-25 | 2019-06-25 | 燕山大学 | 一种混合润滑状态下的轧制力预报方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106964654A (zh) | 一种求解轧制力的建模和仿真方法 | |
Guo et al. | Modification of equation of motion of fluid-conveying pipe for laminar and turbulent flow profiles | |
Chen et al. | Calculation of rolling pressure distribution and force based on improved Karman equation for hot strip mill | |
Hong et al. | A parametric study on forming length in roll forming | |
Mousavi et al. | Three dimensional numerical analyses of asymmetric rolling | |
Xie et al. | Rigid–viscoplastic dynamic explicit FEA of the ring rolling process | |
Bel’skii et al. | Influence of the cross section of hot-rolled steel on the flatness of cold-rolled strip | |
CN106623443B (zh) | 一种厚板粗轧阶段轧制力预报方法 | |
CN109359429B (zh) | 一种板带轧制前滑值计算方法及装置 | |
Sorochan et al. | Mathematical model of plates and strips rolling for calculation of energy power parameters and dynamic loads | |
Wang et al. | Analysis of asymmetrical rolling of unbonded clad sheet by slab method considering vertical shear stress | |
Zhang et al. | A fast rigid-plastic finite element method for online application in strip rolling | |
Kim et al. | A study of the application of upper bound method to the CONFORM process | |
Wang et al. | Mathematical and numerical analysis of cross-directional control for SmartCrown rolls in strip mill | |
Ma et al. | Design of continuous variable curvature roll shape and straightening process research for two-roll straightener of bar | |
Knight et al. | Investigations into the influence of asymmetric factors and rolling parameters on strip curvature during hot rolling | |
Meybodi et al. | A general methodology for bearing design in non-symmetric T-shaped sections in extrusion process | |
Pater | Numerical simulation of the cross wedge rolling process including upsetting | |
Xu et al. | Research and application on slab camber control model in hot rolling | |
CN105013835A (zh) | 冷连轧机组极薄带轧制中基于热凸度的原始辊缝设定方法 | |
CN109926456B (zh) | 一种混合润滑状态下的轧制力预报方法 | |
Hsu et al. | Two analytical models of double-layer clad sheet compression forming based on the upper bound and the slab methods | |
Khrapov et al. | Lagrange-Eulerian method for numerical integration of the gas dynamics equations: parallel implementation on GPUs | |
Shi et al. | A comparison of methods to estimate the roll torque in thin strip rolling | |
Minutolo et al. | Dimensional analysis of a new type of groove for steel rebar rolling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170721 |