CN105022923B - A kind of computational methods of the mutual iteration of roll-force and rolling temperature - Google Patents
A kind of computational methods of the mutual iteration of roll-force and rolling temperature Download PDFInfo
- Publication number
- CN105022923B CN105022923B CN201510423680.2A CN201510423680A CN105022923B CN 105022923 B CN105022923 B CN 105022923B CN 201510423680 A CN201510423680 A CN 201510423680A CN 105022923 B CN105022923 B CN 105022923B
- Authority
- CN
- China
- Prior art keywords
- roll
- force
- rolling temperature
- rolling
- coefficient
- 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
Links
Landscapes
- Control Of Metal Rolling (AREA)
Abstract
The present invention relates to a kind of roll-force and the computational methods of the mutual iteration of rolling temperature, after primary Calculation roll-force, according to the corresponding roll flattening radius of gained tube rolling simulation and rolling temperature, rolling temperature according to obtained by calculating calculates resistance of deformation, coefficient of friction to reach the purpose according to rolling temperature calculating roll-force, then whether the relatively more front and rear roll-force calculated is of substantially equal, if the two difference very little, namely reached the condition of convergence, then iteration terminates, and takes and finally calculates obtained roll-force, rolling temperature as result of calculation.If being unsatisfactory for the condition of convergence, roll flattening radius and rolling temperature are calculated according to the roll-force calculated again, roll-force is then calculated according to rolling temperature again, and judges whether to reach the condition of convergence again.The present invention take into account roll-force, both rolling temperatures in real production process and influence each other, the fact that there is Thermal-mechanical Coupling, pass through the method for iterative numerical, by the calculating process of roll-force and rolling temperature as an entirety, reduction rolling power and the interactional objective circumstances of rolling temperature, to improve the precision of prediction of roll-force and rolling temperature model simultaneously.
Description
Technical field
The present invention relates to a kind of roll-force and the computational methods of the mutual iteration of rolling temperature.
Background technology
In the Hot Rolled Steel Producing operation of rolling, rolling procedure is to formulate one of most basic content of hot-rolled production process.And
The quality of rolling procedure, directly affects the difficulty of subsequent control and determines the quality of final products, so formulating rational rolling
Code is significant to improving the steel plate strip quality of production.And rational rolling procedure is formulated, accurate mathematical modeling
It is crucial.Wherein, prediction of rolling force model, rolling temperature forecast model is the key model during rolling procedure making.It is accurate
The size of true prediction roll-force and rolling temperature during steel plate rolling is the most basic guarantor of rational rolling procedure
Barrier.In Hot Strip Rolling, rolling temperature directly affects resistance of deformation size of the rolled piece in the operation of rolling and then affected
The size of roll-force, and on the other hand, rolled piece temperature change caused by friction, plastic deformation is how many again straight with roll-force size
Correlation is connect, so, Thermal-mechanical Coupling relation is there is between rolling temperature and roll-force.And in the prior art, calculating rolling
When power and rolling temperature, be all by two model separate computations, it is artificial to have separated influencing each other for roll-force and rolling temperature,
So as to cause predicted roll-force and rolling temperature and measured value to there is certain error.
The content of the invention
The invention aims to overcome the deficiencies in the prior art mutually to be changed there is provided a kind of roll-force and rolling temperature
The computational methods in generation.Roll-force and temperature computation will be worked as integral, reduction rolling power and the interactional reality of rolling temperature
Border situation, to reach while improving roll-force and the effect of rolling temperature precision of prediction.
The present invention is to be achieved through the following technical solutions：
The computational methods of the mutual iteration of roll-force and rolling temperature involved in the present invention, realize that step is as follows：
1st, a kind of computational methods of the mutual iteration of roll-force and rolling temperature, it is characterised in that comprise the following steps：
1) it is positive integer and i >=1 to take counter i=1, wherein i, calculates the water for considering the deformed area contact arc length after flattening
Flat projectionmm；Roller radius R after presetting flatteningi', resistance to plane deformation KiAnd friction coefficient μiInitially
Value, i.e. Ri'=R, Ki=K0, μi=μ0, R is working roll original radius, mm, Δ h=h1-h2For drafts, h1、h2Respectively plate
Tape entry, exit thickness, mm, K0For resistance to plane deformation initial value, MPa, μ0For coefficient of friction initial value；
2) roll-force P is determinedi=BmldiKiQp, kN；In formula,For the mean breadth of rolled piece in deformed area, B, b
Before respectively rolling, wide, the mm that rolls back plate；QpFor plane stress state coefficient；
3) according to step 2) determine roll-force Pi, it is determined that roller radius after flatteningmm；According to
Roller radius R after flatteningi+1' calculate contact arc length
4) according to step 2) determine roll-force Pi, calculate rolling temperature ti, DEG C；Wherein, the temperature produced by plastic deformation
Rise℃；Wherein, ρ is rolled piece proportion, kg/m3；C is rolled piece specific heat, J/ (kgm3), ηPFor plasticity
Deform coefficient of thermal efficiency；Friction Temperature Rise℃；Wherein, ηfFor friction thermal efficiency system
Number, f is advancing slip value, b sliding values after being；
5) according to step 4) obtained by rolling temperature ti, calculate resistance of deformationMPa, root
According to resistance of deformation σsiCalculation Plane resistance of deformation Ki+1=1.155 σsi, wherein, ε is true strain；For rate of deformation, s-1；For thermodynamic temperature, K；a1~a4For undetermined coefficient；
6) according to step 4) obtained by rolling temperature ti, calculate friction coefficient μi+1=K1K2K3K4(0.877-
0.0039ti), wherein, K1For correction factor relevant with roll material and surface state, K2For amendment relevant with rolling stock species
Coefficient；K3For correction factor relevant with mill speed, K4Coefficient is influenceed for tension force；
7) roll-force P is determinedi+1, wherein, Pi+1=Ki+1QpBmld(i+1)；
8) judge whether to meet the condition of convergence | (Pi+1-Pi)/Pi| ＜ 10-5；
If 9) be unsatisfactory for the condition of convergence, following steps are carried out：
1. calculateWith
2. rolling temperature t is calculatedi+1, wherein, calculating rolling temperature ti+1When,
3. roller radius, i.e. P after roll-force, rolling temperature and flattening are updatedi=Pi+1, ti=ti+1, Ri+1'=Ri+2'；Meter
Device i is calculated from increasing 1, i.e. i=i+1
It is then back to step 5)；
If 10) meet the condition of convergence, roll-force Pi+1As final roll-force, and according to this roll-force Pi+1Calculating is obtained
Final rolling temperature ti+1。
Compared with prior art, the present invention is beneficial has the technical effect that：Based on roll-force in actual production process and rolling
The calculating process of rolling temperature and roll-force is regarded one by the interactional objective fact of temperature, the method circulated using numerical value
It is overall, to reach while the technique effect of roll-force and rolling temperature is improved, to reduce calculation error, raising model prediction essence
Degree.
Brief description of the drawings
Fig. 1 is structural representation of the invention.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
Referring to Fig. 1, calculation flow charts of the Fig. 1 for the present invention.
The computational methods of the mutual iteration of roll-force and rolling temperature involved in the present invention, realize that step is as follows：
1st, a kind of computational methods of the mutual iteration of roll-force and rolling temperature, it is characterised in that comprise the following steps：
1) it is positive integer and i >=1 to take counter i=1, wherein i, calculates the water for considering the deformed area contact arc length after flattening
Flat projectionmm；Roller radius R after presetting flatteningi', resistance to plane deformation KiAnd friction coefficient μiInitial value,
That is Ri'=R, Ki=K0, μi=μ0, R is working roll original radius, mm, Δ h=h1-h2For drafts, h1、h2Respectively strip enters
Mouth, exit thickness, mm, K0For resistance to plane deformation initial value, MPa, μ0For coefficient of friction initial value；
2) roll-force P is determinedi=BmldiKiQp, kN；In formula,For the mean breadth of rolled piece in deformed area, B, b
Before respectively rolling, wide, the mm that rolls back plate；QpFor plane stress state coefficient；
3) according to step 2) determine roll-force Pi, it is determined that roller radius after flatteningmm；According to
Roller radius R after flatteningi+1' calculate contact arc length
4) according to step 2) determine roll-force Pi, calculate rolling temperature ti, DEG C；Wherein, the temperature produced by plastic deformation
Rise℃；Wherein, ρ is rolled piece proportion, kg/m3；C is rolled piece specific heat, J/ (kgm3), ηPFor plasticity
Deform coefficient of thermal efficiency；Friction Temperature Rise℃；Wherein, ηfFor friction thermal efficiency system
Number, f is advancing slip value, b sliding values after being；
5) according to step 4) obtained by rolling temperature ti, calculate resistance of deformationMPa, root
According to resistance of deformation σsiCalculation Plane resistance of deformation Ki+1=1.155 σsi, wherein, ε is true strain；For rate of deformation, s-1；For thermodynamic temperature, K；a1~a4For undetermined coefficient；
6) according to step 4) obtained by rolling temperature ti, calculate friction coefficient μi+1=K1K2K3K4(0.877-
0.0039ti), wherein, K1For correction factor relevant with roll material and surface state, K2For amendment relevant with rolling stock species
Coefficient；K3For correction factor relevant with mill speed, K4Coefficient is influenceed for tension force；
7) roll-force P is determinedi+1, wherein, Pi+1=Ki+1QpBmld(i+1)；
8) judge whether to meet the condition of convergence | (Pi+1-Pi)/Pi| ＜ 10-5；
If 9) be unsatisfactory for the condition of convergence, following steps are carried out：
1. calculateWith
2. rolling temperature t is calculatedi+1, wherein, calculating rolling temperature ti+1When,
3. roller radius, i.e. P after roll-force, rolling temperature and flattening are updatedi=Pi+1, ti=ti+1, Ri+1'=Ri+2'；Meter
Device i is calculated from increasing 1, i.e. i=i+1
It is then back to step 5)；
If 10) meet the condition of convergence, roll-force Pi+1As final roll-force, and according to this roll-force Pi+1Calculating is obtained
Final rolling temperature ti+1。
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
Any modifications, equivalent substitutions and improvements made within refreshing and principle etc., should be included in the scope of the protection.
Claims (1)
1. the computational methods of roll-force and the mutual iteration of rolling temperature, it is characterised in that comprise the following steps：
1) it is positive integer and i >=1 to take counter i=1, wherein i, calculates and considers that the level of the deformed area contact arc length after flattening is thrown
Shadowmm；Roller radius R after presetting flatteningi', resistance to plane deformation KiAnd friction coefficient μiInitial value, i.e.,
Ri'=R, Ki=K0, μi=μ0, R is working roll original radius, mm, Δ h=h1-h2For drafts, h1、h2Respectively strip enters
Mouth, exit thickness, mm, K0For resistance to plane deformation initial value, MPa, μ0For coefficient of friction initial value；
2) roll-force P is determinedi=BmldiKiQp, kN；In formula,For the mean breadth of rolled piece in deformed area, B, b difference
For before rolling, wide, the mm that rolls back plate；QpFor plane stress state coefficient；
3) according to step 2) determine roll-force Pi, it is determined that roller radius after flatteningmm；According to flattening
Roller radius R afterwardsi+1' calculate contact arc length
4) according to step 2) determine roll-force Pi, calculate rolling temperature ti, DEG C；Wherein, the temperature rise produced by plastic deformation℃；Wherein, ρ is rolled piece proportion, kg/m3；C is rolled piece specific heat, J/ (kgm3), ηPBecome for plasticity
Shape coefficient of thermal efficiency；Friction Temperature Rise℃；Wherein, ηfFor friction coefficient of thermal efficiency,
F is advancing slip value, b sliding values after being；
5) according to step 4) obtained by rolling temperature ti, calculate resistance of deformationMPa, according to change
Shape drag σsiCalculation Plane resistance of deformation Ki+1=1.155 σsi, wherein, ε is true strain；For rate of deformation, s-1；For thermodynamic temperature, K；a1~a4For undetermined coefficient；
6) according to step 4) obtained by rolling temperature ti, calculate friction coefficient μi+1=K1K2K3K4(0.877-0.0039ti), its
In, K1For correction factor relevant with roll material and surface state, K2For correction factor relevant with rolling stock species；K3For with
The relevant correction factor of mill speed, K4Coefficient is influenceed for tension force；
7) roll-force P is determinedi+1, wherein, Pi+1=Ki+1QpBmld(i+1)；
8) judge whether to meet the condition of convergence | (Pi+1-Pi)/Pi| ＜ 10-5；
If 9) be unsatisfactory for the condition of convergence, following steps are carried out：
1. calculate
2. rolling temperature t is calculatedi+1, wherein, calculating rolling temperature ti+1When,
3. roller radius, i.e. P after roll-force, rolling temperature and flattening are updatedi=Pi+1, ti=ti+1, Ri+1'=Ri+2'；Calculator i
From increasing 1, i.e. i=i+1
It is then back to step 5)；
If 10) meet the condition of convergence, roll-force obtains final rolling temperature as final roll-force, and according to this tube rolling simulation
Degree.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510423680.2A CN105022923B (en) | 2015-07-19 | 2015-07-19 | A kind of computational methods of the mutual iteration of roll-force and rolling temperature |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510423680.2A CN105022923B (en) | 2015-07-19 | 2015-07-19 | A kind of computational methods of the mutual iteration of roll-force and rolling temperature |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105022923A CN105022923A (en) | 2015-11-04 |
| CN105022923B true CN105022923B (en) | 2017-10-31 |
Family
ID=54412889
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510423680.2A Active CN105022923B (en) | 2015-07-19 | 2015-07-19 | A kind of computational methods of the mutual iteration of roll-force and rolling temperature |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105022923B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106557621A (en) * | 2016-11-07 | 2017-04-05 | 南阳师范学院 | A kind of modeling and simulation method for calculating roll elastic flattening radius |
| CN110802114B (en) * | 2019-10-30 | 2021-08-03 | 中冶陕压重工设备有限公司 | Method for calculating rolling force of cold-rolled strip |
| CN111036693B (en) * | 2019-11-15 | 2022-02-22 | 中冶华天工程技术有限公司 | Rolling force energy checking and calculating method for hot-rolled high-speed wire rod |
| CN112453071B (en) * | 2020-11-17 | 2022-07-01 | 太原理工大学 | Method for predicting rolling force and thickness of each layer of cold-rolled metal composite plate |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3628358A (en) * | 1969-10-07 | 1971-12-21 | Gen Electric | Method of revising workpiece temperature estimates or measurements using workpiece deformation behavior |
| CN101733289A (en) * | 2009-12-23 | 2010-06-16 | 北京理工大学 | Method for dynamically setting tandem rolling schedule of hot rolling strip |
| CN102941232A (en) * | 2012-11-12 | 2013-02-27 | 东北大学 | Method for controlling hot continuous rolling finish rolling process |
| CN103028603A (en) * | 2012-11-29 | 2013-04-10 | 一重集团大连设计研究院有限公司 | Method for optimizing rolling scope of hot strip rolling |
-
2015
- 2015-07-19 CN CN201510423680.2A patent/CN105022923B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3628358A (en) * | 1969-10-07 | 1971-12-21 | Gen Electric | Method of revising workpiece temperature estimates or measurements using workpiece deformation behavior |
| CN101733289A (en) * | 2009-12-23 | 2010-06-16 | 北京理工大学 | Method for dynamically setting tandem rolling schedule of hot rolling strip |
| CN102941232A (en) * | 2012-11-12 | 2013-02-27 | 东北大学 | Method for controlling hot continuous rolling finish rolling process |
| CN103028603A (en) * | 2012-11-29 | 2013-04-10 | 一重集团大连设计研究院有限公司 | Method for optimizing rolling scope of hot strip rolling |
Non-Patent Citations (2)
| Title |
|---|
| 6016铝合金热变形行为及热连轧过程数值模拟研究;刁金鹏;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20120716(第8期);第B022-137页 * |
| 板带热连轧轧制力及其设定的研究;刘洋;《中国优秀硕士学位论文全文数据库 工程科技I辑》;20071215(第6期);第B022-138页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105022923A (en) | 2015-11-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105022923B (en) | A kind of computational methods of the mutual iteration of roll-force and rolling temperature | |
| CN102941232B (en) | Method for controlling hot continuous rolling finish rolling process | |
| CN103302104A (en) | Method for manufacturing hot rolled silicon steel | |
| CN106483845B (en) | A kind of method and system improving the forecast of straightener force-power parameter | |
| CN103761370A (en) | Method for predicting strip surface coefficients of heat transfer in hot rolling procedures | |
| CN101507978B (en) | Worker roller abrasion forecasting method of cold continuous rolling device based on mechanism and operation condition | |
| CN102601128B (en) | Temperature difference control method for rolling pieces for steckel mill | |
| CN104298884A (en) | Finite element and finite difference coupling method for fast calculating rolled piece section temperature | |
| CN103028615A (en) | Method for predicting temperature evolution in hot continuous rolling process of strip steel | |
| CN105363794A (en) | Finish rolling energy-saving control method based on mechanical performance prediction model and rolling energy consumption model | |
| CN103008360B (en) | A kind of method determining cold rolling mill work roller temperature field and thermal expansion | |
| CN111570534B (en) | Method for predicting influence of hot rolling characteristics of acid rolling mill group on friction coefficient | |
| CN105234186B (en) | Cold continuous rolling process is controlled to the rolling schedule optimization method of target with electric power consumption per ton steel | |
| CN103878186A (en) | Method for determining hot rolled strip steel laminar cooling temperature | |
| CN103100564A (en) | Novel rolling process self-adaptive control method | |
| CN106475421A (en) | A kind of control method of middle base supplied materials convexity and device | |
| CN105234188B (en) | Mill speed optimization method in cold continuous rolling process with Result coutrolling as target | |
| CN105013835A (en) | Original roller seam setting method based on thermal crown in ultra-thin strip rolling conducted by cold continuous rolling unit | |
| CN103586289B (en) | The draught pressure establishing method of rougher of hot strip mill district edge rolling | |
| CN102553945A (en) | Abnormal shape forecasting method suitable for four-high rolling mill | |
| Tao et al. | Improvement of 3-D FEM coupled model on strip crown in hot rolling | |
| CN107234135B (en) | One kind being suitable for hot tandem and exports belt steel surface roughness control method | |
| CN110227721A (en) | A kind of control method for mill roll-gap in speed-change process | |
| CN104537136A (en) | A method for compensating defect of diameter reducer of roller of six-roller mill | |
| CN105537284A (en) | Temperature forecasting method for finish rolling inlet during hot continuous rolling |
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 |