CN114912248A - 热镀锌机组运行过程中沉没辊系受力预报方法 - Google Patents
热镀锌机组运行过程中沉没辊系受力预报方法 Download PDFInfo
- Publication number
- CN114912248A CN114912248A CN202210359339.5A CN202210359339A CN114912248A CN 114912248 A CN114912248 A CN 114912248A CN 202210359339 A CN202210359339 A CN 202210359339A CN 114912248 A CN114912248 A CN 114912248A
- Authority
- CN
- China
- Prior art keywords
- roller
- roll
- strip steel
- sink
- stress
- 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
Links
- 238000005246 galvanizing Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013277 forecasting method Methods 0.000 title claims abstract description 13
- 229910000831 Steel Inorganic materials 0.000 claims description 84
- 239000010959 steel Substances 0.000 claims description 84
- 230000000087 stabilizing effect Effects 0.000 claims description 44
- 239000011701 zinc Substances 0.000 claims description 37
- HCHKCACWOHOZIP-UHFFFAOYSA-N Zinc Chemical compound [Zn] HCHKCACWOHOZIP-UHFFFAOYSA-N 0.000 claims description 29
- 229910052725 zinc Inorganic materials 0.000 claims description 29
- 238000004364 calculation method Methods 0.000 claims description 28
- 230000005484 gravity Effects 0.000 claims description 18
- 238000005096 rolling process Methods 0.000 claims description 14
- 238000004519 manufacturing process Methods 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 10
- 239000007788 liquid Substances 0.000 claims description 7
- 239000002893 slag Substances 0.000 claims description 5
- 238000003780 insertion Methods 0.000 claims description 4
- 230000037431 insertion Effects 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000008021 deposition Effects 0.000 claims description 2
- 238000004148 unit process Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 238000005457 optimization Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
-
- C—CHEMISTRY; METALLURGY
- C23—COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; CHEMICAL SURFACE TREATMENT; DIFFUSION TREATMENT OF METALLIC MATERIAL; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL; INHIBITING CORROSION OF METALLIC MATERIAL OR INCRUSTATION IN GENERAL
- C23C—COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; SURFACE TREATMENT OF METALLIC MATERIAL BY DIFFUSION INTO THE SURFACE, BY CHEMICAL CONVERSION OR SUBSTITUTION; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL
- C23C2/00—Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor
- C23C2/04—Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor characterised by the coating material
- C23C2/06—Zinc or cadmium or alloys based thereon
-
- C—CHEMISTRY; METALLURGY
- C23—COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; CHEMICAL SURFACE TREATMENT; DIFFUSION TREATMENT OF METALLIC MATERIAL; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL; INHIBITING CORROSION OF METALLIC MATERIAL OR INCRUSTATION IN GENERAL
- C23C—COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; SURFACE TREATMENT OF METALLIC MATERIAL BY DIFFUSION INTO THE SURFACE, BY CHEMICAL CONVERSION OR SUBSTITUTION; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL
- C23C2/00—Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor
- C23C2/34—Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor characterised by the shape of the material to be treated
- C23C2/36—Elongated material
- C23C2/40—Plates; Strips
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Organic Chemistry (AREA)
- Materials Engineering (AREA)
- Mechanical Engineering (AREA)
- Metallurgy (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Coating With Molten Metal (AREA)
Abstract
本发明提供了热镀锌机组运行过程中沉没辊系受力预报方法，涉及沉没辊系受力预报技术领域，该方法是根据锌锅内浸没在锌液内沉没辊系工作状态，实现对沉没辊、稳定辊、矫正辊的受力预报，本发明能够根据冷轧带钢的现场生产情况，充分结合热镀锌机组的设备特点，通过对机组沉没辊系相关参数收集，有效解决了在正常工况下沉没辊系受力计算，实现对沉没辊系动态载荷的预报，为后续沉没辊系受力动态控制及沉没辊系工艺参数的优化提供了理论基础及指导作用。
Description
技术领域
本发明涉及热镀锌技术领域，特别涉及热镀锌机组运行过程中沉没辊系受力预报方法。
背景技术
热镀锌是板带产品的最后一道工序，板带表面质量，如镀锌厚度均匀性、上下表面镀层厚度偏差量是热镀板带表面质量的主要衡量指标。板带在热镀锌工艺时，首先经过连续退火工艺过程，而后进入锌锅内的沉没辊系(沉没辊、稳定辊及矫正辊构成)，随后出锌锅进入气刀处，刮去板带表面的多余的锌层，最后冷却成形。
由于锌锅内含有很多锌渣，如浮渣和沉渣，同时锌锅是处于高温、锌液场等多物理场作用下沉没辊系各辊受力比较复杂。为研究实现沉没辊系在锌锅内处于最佳工作状态，处于振动幅度最小，需要掌握辊系的受力状态，实现对其精确化控制，从而对沉没辊系工艺参数优化提供理论基础。
现有技术中，如发明CN101201276-沉没辊两端轴颈受力的检测方法和装置中，通过检测带钢左右两边部的振动信号，来检测沉没辊两端轴颈受力，但不能在线实现，即不能热镀锌机组运行过程中沉没辊系受力预报。
可见，目前现场难以实现沉没辊系各辊的受力状态预报，也就无法实现相关工艺参数优化。因此，为了保证热镀锌机组的正常生产，提高带钢产品镀锌质量，减小沉没辊系在运行过程中振动及磨损，亟需热镀锌机组运行过程中沉没辊系受力预报方法。
发明内容
有鉴于此，本发明提供了热镀锌机组运行过程中沉没辊系受力预报方法，该方法是根据锌锅内浸没在锌液内沉没辊系工作状态，实现对沉没辊、稳定辊、矫正辊的受力预报，能够在热镀锌机组运行过程中进行沉没辊系受力预报，进而为相关工艺参数优化奠定了基础。
为此，本发明提供了以下技术方案：
本发明提供了一种热镀锌机组运行过程中沉没辊系受力预报方法，所述方法包括以下步骤：
2、一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，所述方法包括以下步骤：
步骤A、设定沉没辊系中辊的辊号λ，初值λ＝1；
步骤B、获取镀锌机组沉没辊系辊的设备参数以及带钢生产工艺参数；
步骤C、根据几何法推导的公式，并结合实际工况的矫正辊插入量确定沉没辊系各辊在工作状态下的包角θ；
步骤D、将辊接触区带钢沿辊周向划分m个横向条元；再沿任意辊接触区带钢横向条元的横向划分n个子单元，并计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θi，zj)；
步骤E、计算辊号为λ的辊对带钢单元的接触压力Ni,j和滚动摩擦力fi,j，计算公式为：
fi,j＝τk·Ni,j；
其中，Ti,j为带钢单元旋入端的张力，Ti,j＝T0；T0为带钢初始设定张力；T′i,j为带钢单元旋出端的张力，
Δgi,j为带钢单元受到重力与浮力的合力，
H、ρFe、ρZn、v、B为带钢生产工艺参数，分别表示带材的厚度、带钢密度、锌锅内高温下锌液密度、带钢运行速度、带材宽度；R为辊的设备参数，表示辊半径；τk为辊面子单元所在辊的全周向区域摩擦系数；
步骤F、计算带钢对各辊的辊面滚动摩擦力矩M以及各辊受到的总接触压力N：
步骤G、计算各辊的轴端支撑力F以及摩擦阻力矩Md，计算公式为：
Md＝μFd/2；
其中，ΔG为沉没辊系重力与浮力的合力，ΔG＝G-ρZngV；其中，g为重力加速常量；V、G、L、d、μ为辊的设备参数，分别表示辊的有效体积、辊的重力、辊的辊身长、辊轴颈直径、辊的轴颈与轴瓦间初始摩擦系数；B为带钢生产工艺参数，表示带材宽度；
步骤H、判断λ≤3是否成立，若λ≤3成立，则令λ＝λ+1，并返回执行步骤B，进行下一个辊的受力预报；若λ≤3不成立，则输出沉没辊系各辊受到的总接触压力N、带钢的摩擦驱动力矩M、轴端支持力F及其摩擦阻力矩Md。
进一步地，所述沉没辊系包括：沉没辊、稳定辊及矫正辊。
进一步地，确定沉没辊系各辊在工作状态下的包角，包括：
判断λ＝1是否成立，若λ＝1成立，则：
并令R*＝R，θ*＝θ，
辊包角θ的正向角θ+、负向角θ-为：
若不成立，则
判断λ＝2是否成立，若λ＝2成立，则
若λ＝2不成立，则
其中，θ0为带钢与辊包角初始值，Δs为矫正辊插入量，R*、θ*、
是过程变量，初值均为0。
进一步地，分别计算与辊接触带钢的任意(i,j)子单元坐标值(θi,zj)包括：
进一步地，还包括：计算带钢对各辊的接触压力N的方向角φ，计算公式为：
进一步地，辊面子单元所在辊的全周向区域摩擦系数τk的计算方式如下：
在各辊沿其辊面轴向均匀选取若干个子单元n*，并通过测量传感器测量得出在各辊表面上对应这些点的锌渣沉积厚度hk，其中k＝1,2,...,n*，n*≤n；
计算辊面子单元所在辊的全周向区域摩擦系数τk，计算公式为：
其中，
为辊的辊面平均锌渣沉积厚度，
τ为带钢生产工艺参数，表示带钢与辊的辊面间的初始摩擦系数。
本发明的优点和积极效果：本发明能够根据冷轧成品带钢热镀锌的现场生产情况，充分结合热镀锌机组的设备特点，通过对机组沉没辊系相关参数收集，基于沉没辊系在线锌渣测量传感器、辊系动态包角求解算法，利用本发明在线辊系离散条元法计算方法和技术，可有效解决了在正常工况下沉没辊系受力计算，实现对沉没辊系动态载荷的预报，为后续沉没辊系受力动态控制及沉没辊系工艺参数的优化提供了理论基础及指导作用。而目前热镀锌机组无法实现辊系在线锌渣测量，也无法对辊系包角动态计算预测，更无法实现沉没辊系动态载荷预报。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中沉没辊系组成示意图；
图2(a)为本发明实施例中对与沉没辊接触带钢划分横向条元；
图2(b)为本发明实施例中对带钢横向条元划分子单元；
图3(a)为本发明实施例中对与稳定辊接触带钢划分横向条元；
图3(b)为本发明实施例中对带钢横向条元划分子单元；
图4(a)为本发明实施例中对与矫正辊接触带钢划分横向条元；
图4(b)为本发明实施例中对带钢横向条元划分子单元；
图5为本发明实施例中沉没辊系各辊受力预报流程图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
以如图1所示的某镀锌机组沉没辊系各辊(沉没辊系包括沉没辊、稳定辊、矫正辊)为例，结合图5示出的沉没辊系各辊受力预报流程图，对本发明提供的热镀锌机组运行过程中沉没辊系受力预报方法进行详细说明。
实施例1：
步骤A1，设定沉没辊系中辊的辊号λ初值λ＝1，矫正辊、沉没辊和稳定辊的辊号依次为1、2、3；并依次取λ＝2、λ＝3。
步骤B1，收集镀锌机组沉没辊系设备参数，包括辊系半径R，辊重力G，辊的有效体积V，辊的辊身长为L，辊轴颈直径d，辊的轴颈与轴瓦间初始摩擦系数μ，相关设备参数如表1(镀锌机组沉没辊、稳定辊设备参数)、表2(镀锌机组沉没辊、稳定辊轴端零件参数)所示。
表1
| 名称 | 半径/mm | 重力/N | 浮力/N | 有效体积/m3 | 辊身长/m |
| 沉没辊 | 300 | 9000 | 8760 | 0.1253 | 1.6 |
| 稳定辊 | 100 | 2000 | 2500 | 0.0326 | 1.6 |
| 矫正辊 | 100 | 2000 | 2500 | 0.0326 | 1.6 |
表2
步骤C1，收集带钢生产工艺参数，包括带钢运行速度v＝120m/min，带钢初始设定张力T0＝15kN，带材的宽度B＝1000mm及厚度H＝1.2mm，带钢密度ρFe＝7850kg/m3，带钢与沉没辊、稳定辊、矫正辊包角初始值分别为130°、0°、5°，锌锅内高温下锌液密度ρZn＝6510kg/m3，带钢与沉没辊、稳定辊、矫正辊的辊面间的初始摩擦系数都为0.04，矫正辊插入量Δs＝5mm，设定相关过程变量初值R*＝0，θ*＝0，
步骤D1，根据辊的重力G、有效体积V和锌锅内高温下锌液密度ρZn计算辊的重力与浮力的合力ΔG；计算公式为：ΔG＝G-ρZngV；其中，g为重力加速常量，取9.8N/kg；可计算出沉没辊、稳定辊、矫正辊重力与浮力的合力分别为994.7N、258.9N、258.9N；
步骤E1，确定沉没辊系各辊在工作状态下的包角：
E11)判断λ＝1是否成立；若λ＝1成立，则：
并令R*＝R，θ*＝θ，
辊包角的正向角θ+、负向角θ-为：
若不成立，则转入E12)；
E12)辊包角θ为：
判断λ＝2是否成立；若λ＝2成立，则
若不成立，则
沉没辊系各辊的计算结果如表3所示。
表3
| 沉没辊 | 稳定辊 | 矫正辊 | |
| 包角/° | 135 | 25.2 | 65.3 |
| 正向角/° | 40 | 0 | 32.65 |
| 负向角/° | 95 | -25.2 | -32.65 |
步骤F1，将各辊包角范围内接触带钢划分别沿辊周向划分条元、沿轴向划分子单元，并分别计算与辊接触带钢的任意(i,j)子单元坐标值(θi,zj)，计算公式为：
其中，θ为辊的包角，θ-为辊包角的负向角；m为将辊接触区带钢沿辊周向划分的横向条元数目，n为沿任意横向条元横向划分的子单元数目。沉没辊、稳定辊及矫正辊所划分的数目分别为5,5,5；再沿任意沉没辊、稳定辊及矫正辊接触区带钢横向条元的横向所划分分别为5,5,5。沉没辊、稳定辊、矫正辊上任意(i,j)子单元坐标值的计算结果如表4、5、6所示。
表4
| 单元 | i＝1 | i＝2 | i＝3 | i＝4 | i＝5 |
| j＝1 | (26°，0.32) | (52°，0.32) | (78°，0.32) | (104°，0.32) | (130°，0.32) |
| j＝2 | (26°，0.64) | (52°，0.64) | (78°，0.64) | (104°，0.64) | (130°，0.64) |
| j＝3 | (26°，0.96) | (52°，0.96) | (78°，0.96) | (104°，0.96) | (130°，0.96) |
| j＝4 | (26°，1.28) | (52°，1.28) | (78°，1.28) | (104°，1.28) | (130°，1.28) |
| j＝5 | (26°，1.6) | (52°，1.6) | (78°，1.6) | (104°，1.6) | (130°，1.6) |
表5
| 单元 | i＝1 | i＝2 | i＝3 | i＝4 | i＝5 |
| j＝1 | (5°，0.32) | (10°，0.32) | (15°，0.32) | (20°，0.32) | (25°，0.32) |
| j＝2 | (5°，0.64) | (10°，0.64) | (15°，0.64) | (20°，0.64) | (25°，0.64) |
| j＝3 | (5°，0.96) | (10°，0.96) | (15°，0.96) | (20°，0.96) | (25°，0.96) |
| j＝4 | (5°，1.28) | (10°，1.28) | (15°，1.28) | (20°，1.28) | (25°，1.28) |
| j＝5 | (5°，1.6) | (10°，1.6) | (15°，1.6) | (20°，1.6) | (25°，1.6) |
表6
| 单元 | i＝1 | i＝2 | i＝3 | i＝4 | i＝5 |
| j＝1 | (13°，0.32) | (26°，0.32) | (39°，0.32) | (52°，0.32) | (65°，0.32) |
| j＝2 | (13°，0.64) | (26°，0.64) | (39°，0.64) | (52°，0.64) | (65°，0.64) |
| j＝3 | (13°，0.96) | (26°，0.96) | (39°，0.96) | (52°，0.96) | (65°，0.96) |
| j＝4 | (13°，1.28) | (26°，1.28) | (39°，1.28) | (52°，1.28) | (65°，1.28) |
| j＝5 | (13°，1.6) | (26°，1.6) | (39°，1.6) | (52°，1.6) | (65°，1.6) |
步骤G1，在各辊沿其辊面轴向均匀选取若干个子单元n*，通过测量传感器测量得出在各辊表面上对应这些点的锌渣沉积厚度hk，其中k＝1,2,..n.*,；
运行过程线上测量，测量传感器沿辊轴向均匀移动，等间隔在辊面采样出若干点，每个采样点对应带钢子单元；
计算辊面子单元所在辊的全周向区域摩擦系数τk，计算公式为：
其中，
为辊的辊面平均锌渣沉积厚度，
沉没辊、稳定辊及矫正辊的辊面子单元所在辊的全周向区域摩擦系数
的计算结果如表7、表8、表9所示。
表7
表8
表9
步骤H1，当θi在辊的包角范围内，计算该辊对带钢单元的接触压力Ni,j和滚动摩擦力fi,j，计算公式为：
fi,j＝τk·Ni,j；
其中，
Ti,j＝T0；
步骤I1，计算带钢对各辊的辊面滚动摩擦力矩M以及各辊受到的总接触压力N及其接触角度φ计算公式为：
所计算带钢对沉没辊系各辊的辊面滚动摩擦力矩、接触压力及方向角如表10所示。
表10
| 滚动摩擦力矩/(kN.m) | 接触压力(kN) | 接触压力方向角 | |
| 沉没辊 | 8.3 | 13.2 | 0.64 |
| 稳定辊 | 5.4 | 6.5 | 0.27 |
| 矫正辊 | 5.4 | 6.5 | 0.27 |
步骤J1，计算各辊的轴端支撑力F以及摩擦阻力矩Md，计算公式为：
Md＝μFd/2；
计算得到的沉没辊系各辊的轴端支撑力及摩擦阻力矩如表11所示。
表11
| 摩擦阻力矩/(kN.m) | 轴端支撑力(kN) | |
| 沉没辊 | 0.86 | 0.91 |
| 稳定辊 | 0.49 | 0.83 |
| 矫正辊 | 0.49 | 0.83 |
步骤K1，判断λ≤3是否成立；若成立，则令λ＝λ+1，并转入步骤B)；若不成立，则输出沉没辊系各辊受到的总接触压力N、带钢的摩擦驱动力矩M、轴端支持力及其摩擦阻力矩Md。
最后输出镀锌机组沉没辊系受力的预报值与实测值(实测值是通过布置在沉没辊系各辊轴端的力传感器测量获得的)对比，如表12所示。
表12
对于沉没辊、稳定辊及矫正辊的滚动摩擦力矩的实测值与预报值相对误差分别为10.7％、12.9％、6.9％，对于沉没辊、稳定辊及矫正辊的接触压力的实测值与预报值相对误差分别为3.9％、10.1％、4.8％，对于沉没辊、稳定辊及矫正辊接触压力方向角的实测值与预报值相对误差分别为11.1％、12.9％、12.9％，对于沉没辊、稳定辊及矫正辊摩擦阻力矩的实测值与预报值相对误差分别为11.3％、2.1％、12.2％，对于沉没辊、稳定辊及矫正辊轴端支撑力的实测值与预报值相对误差分别为2.2％、12.1％、12.1％，以上预测值与实测值都在机组工况要求15％以内，这证实本发明有很好实际工程价值。
实施例2：
步骤A2，设定沉没辊系中辊的辊号λ初值λ＝1，矫正辊、沉没辊和稳定辊的辊号依次为1、2、3；并依次取λ＝2、λ＝3。
步骤B2，收集镀锌机组沉没辊系设备参数，包括辊系半径R，辊重力G，辊的有效体积V，辊的辊身长为L，辊轴颈直径d，辊的轴颈与轴瓦间初始摩擦系数μ，相关设备参数如表13(镀锌机组沉没辊、稳定辊、矫正辊设备参数)、表14(镀锌机组沉没辊、稳定辊、矫正辊轴端零件参数)所示。
表13
| 名称 | 半径/mm | 重力/N | 浮力/N | 有效体积/m<sup>3</sup> | 辊身长/m |
| 沉没辊 | 350 | 9780 | 9260 | 0.1361 | 1.6 |
| 稳定辊 | 150 | 2400 | 2340 | 0.0435 | 1.6 |
| 矫正辊 | 150 | 2400 | 2340 | 0.0435 | 1.6 |
表14
步骤C2，收集带钢生产工艺参数，包括带钢运行速度v＝120m/min，带钢初始设定张力T0＝17kN，带材的宽度B＝900mm及厚度H＝0.75mm，带钢密度ρFe＝7850kg/m3，带钢与沉没辊、稳定辊、矫正辊包角初始值分别为132°、0°、7°，锌锅内高温下锌液密度ρZn＝6510kg/m3，带钢与沉没辊、稳定辊、矫正辊的辊面间的初始摩擦系数都为0.05，矫正辊插入量Δs＝7mm，设定相关过程变量初值R*＝0，θ*＝0，
步骤D2，按照与实施例1相同的方式计算辊的重力与浮力的合力ΔG；计算出的沉没辊、稳定辊、矫正辊重力与浮力的合力分别为520N、60N、60N；
步骤E2，按照与实施例1相同的方式确定沉没辊系各辊在工作状态下的包角；沉没辊系各辊的包角计算结果如表15所示。
表15
| 沉没辊 | 稳定辊 | 矫正辊 | |
| 包角/° | 142 | 26.9 | 78.3 |
| 正向角/° | 40 | 0 | 39.16 |
| 负向角/° | 102 | -26.9 | -39.16 |
步骤F2，将各辊包角范围内接触带钢划分别沿辊周向划分条元、沿轴向划分子单元，并分别计算与辊接触带钢的任意(i,j)子单元坐标值(θi,zj)，计算公式为：
其中，θ为辊的包角，θ-为辊包角的负向角；m为将辊接触区带钢沿辊周向划分的横向条元数目，n为沿任意横向条元横向划分的子单元数目。沉没辊、稳定辊及矫正辊所划分的数目分别为6,6,6；再沿任意沉没辊、稳定辊及矫正辊接触区带钢横向条元的横向所划分分别为6,6,6。沉没辊、稳定辊、矫正辊上任意(i,j)子单元坐标值的计算结果如表16、17、18所示。
表16
表17
表18
步骤G2、在各辊沿其辊面轴向均匀选取若干个子单元n*，测量得出在各辊表面上对应这些点的锌渣沉积厚度hk，其中k＝1,2,...,n*；
计算辊面子单元所在辊的全周向区域摩擦系数τk，计算公式为：
其中，
为辊的辊面平均锌渣沉积厚度，
沉没辊、稳定辊及矫正辊的辊面子单元所在辊的全周向区域摩擦系数
的计算结果如表19、表20、表21所示。
表19
表20
表21
步骤H2，当θi在辊的包角范围内，计算该辊对带钢单元的接触压力Ni,j和滚动摩擦力fi,j，计算公式为：
fi,j＝τk·Ni,j；
其中，
Ti,j＝T0；
步骤I2，计算带钢对各辊的辊面滚动摩擦力矩M以及各辊受到的总接触压力N及其接触角度φ计算公式为：
所计算带钢对沉没辊系各辊的辊面滚动摩擦力矩、接触压力及方向角如表22所示。
表22
| 滚动摩擦力矩/(kN.m) | 接触压力(kN) | 接触压力方向角 | |
| 沉没辊 | 12.3 | 25.3 | 0.79 |
| 稳定辊 | 8.7 | 7.0 | 0.36 |
| 矫正辊 | 9.1 | 8.3 | 0.51 |
步骤J2，计算各辊的轴端支撑力F以及摩擦阻力矩Md，计算公式为：
Md＝μFd/2；
其计算沉没辊系各辊的轴端支撑力及摩擦阻力矩如表23所示。
表23
| 摩擦阻力矩/(kN.m) | 轴端支撑力(kN) | |
| 沉没辊 | 1.58 | 0.79 |
| 稳定辊 | 1.02 | 0.48 |
| 矫正辊 | 1.17 | 0.70 |
步骤K2，判断λ≤3是否成立；若成立，则令λ＝λ+1，并转入步骤B)；若不成立，则输出沉没辊系各辊受到的总接触压力N、带钢的摩擦驱动力矩M、轴端支持力及其摩擦阻力矩Md。
最后输出镀锌机组沉没辊系受力的预报值与实测值(实测值是通过布置在沉没辊系各辊轴端的力传感器测量获得的)对比，如表24所示。
表24
对于沉没辊、稳定辊及矫正辊的滚动摩擦力矩的实测值与预报值相对误差分别为5.1％、10.1％、4.2％，对于沉没辊、稳定辊及矫正辊的接触压力的实测值与预报值相对误差分别为4.5％、6.7％、9.2％，对于沉没辊、稳定辊及矫正辊接触压力方向角的实测值与预报值相对误差分别为3.6％、5.8％、6.3％，对于沉没辊、稳定辊及矫正辊摩擦阻力矩的实测值与预报值相对误差分别为11.3％、10.9％、3.3％，对于沉没辊、稳定辊及矫正辊轴端支撑力的实测值与预报值相对误差分别为11.2％、9.4％、2.9％，以上预测值与实测值都在机组工况要求15％以内，这证实本发明具有很好的实际工程价值。
实施例1、2分别针对两种规格带钢，实施例1针对适用于生产厚度大于1mm带钢的热镀锌机组在工作条件下沉没辊系受力状态的实时监测，实施例2针对厚度小于1mm带钢的热镀锌机组在工作条件下沉没辊系受力状态的实时监测，这两种实施例1、2具有代表性，可作为机组生产相关规格带钢的参考。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (6)
1.一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，所述方法包括以下步骤：
步骤A、设定沉没辊系中辊的辊号λ，初值λ＝1；
步骤B、获取镀锌机组沉没辊系辊的设备参数以及带钢生产工艺参数；
步骤C、根据几何法推导的公式，并结合实际工况的矫正辊插入量确定沉没辊系各辊在工作状态下的包角θ；
步骤D、将辊接触区带钢沿辊周向划分m个横向条元；再沿任意辊接触区带钢横向条元的横向划分n个子单元，并计算与沉没辊接触带钢的任意(i,j)子单元坐标值(θi，zj)；
步骤E、计算辊号为λ的辊对带钢单元的接触压力Ni,j和滚动摩擦力fi,j，计算公式为：
fi,j＝τk·Ni,j；
其中，Ti,j为带钢单元旋入端的张力，Ti,j＝T0；T0为带钢初始设定张力；T′i,j为带钢单元旋出端的张力，
Δgi,j为带钢单元受到重力与浮力的合力，
H、ρFe、ρZn、v、B为带钢生产工艺参数，分别表示带材的厚度、带钢密度、锌锅内高温下锌液密度、带钢运行速度、带材宽度；R为辊的设备参数，表示辊半径；τk为辊面子单元所在辊的全周向区域摩擦系数；
步骤F、计算带钢对各辊的辊面滚动摩擦力矩M以及各辊受到的总接触压力N：
步骤G、计算各辊的轴端支撑力F以及摩擦阻力矩Md，计算公式为：
Md＝μFd/2；
其中，ΔG为沉没辊系重力与浮力的合力，ΔG＝G-ρZngV；其中，g为重力加速常量；V、G、L、d、μ为辊的设备参数，分别表示辊的有效体积、辊的重力、辊的辊身长、辊轴颈直径、辊的轴颈与轴瓦间初始摩擦系数；B为带钢生产工艺参数，表示带材宽度；
步骤H、判断λ≤3是否成立，若λ≤3成立，则令λ＝λ+1，并返回执行步骤B，进行下一个辊的受力预报；若λ≤3不成立，则输出沉没辊系各辊受到的总接触压力N、带钢的摩擦驱动力矩M、轴端支持力F及其摩擦阻力矩Md。
2.根据权利要求1所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，所述沉没辊系包括：沉没辊、稳定辊及矫正辊。
3.根据权利要求2所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，确定沉没辊系各辊在工作状态下的包角，包括：
判断λ＝1是否成立，若λ＝1成立，则：
并令R*＝R，θ*＝θ，
辊包角θ的正向角θ+、负向角θ-为：
若不成立，则
判断λ＝2是否成立，若λ＝2成立，则
若λ＝2不成立，则
其中，θ0为带钢与辊包角初始值，Δs为矫正辊插入量，R*、θ*、
是过程变量，初值均为0。
4.根据权利要求3所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，分别计算与辊接触带钢的任意(i,j)子单元坐标值(θi,zj)包括：
5.根据权利要求1所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，还包括：计算带钢对各辊的接触压力N的方向角φ，计算公式为：
6.根据权利要求1所述的一种热镀锌机组运行过程中沉没辊系受力预报方法，其特征在于，辊面子单元所在辊的全周向区域摩擦系数τk的计算方式如下：
在各辊沿其辊面轴向均匀选取若干个子单元n*，并通过测量传感器测量得出在各辊表面上对应这些点的锌渣沉积厚度hk，其中k＝1,2,...,n*，n*≤n；
计算辊面子单元所在辊的全周向区域摩擦系数τk，计算公式为：
其中，
为辊的辊面平均锌渣沉积厚度，
τ为带钢生产工艺参数，表示带钢与辊的辊面间的初始摩擦系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210359339.5A CN114912248B (zh) | 2022-04-06 | 热镀锌机组运行过程中沉没辊系受力预报方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210359339.5A CN114912248B (zh) | 2022-04-06 | 热镀锌机组运行过程中沉没辊系受力预报方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114912248A true CN114912248A (zh) | 2022-08-16 |
| CN114912248B CN114912248B (zh) | 2024-06-25 |
Family
ID=
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116103597A (zh) * | 2023-03-15 | 2023-05-12 | 燕山大学 | 一种热镀铝锌过程锌锅内温度变化预报方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004043858A (ja) * | 2002-07-10 | 2004-02-12 | Jfe Steel Kk | 連続溶融金属めっき装置のシンクロール |
| CN101201276A (zh) * | 2006-12-14 | 2008-06-18 | 宝山钢铁股份有限公司 | 沉没辊两端轴颈受力的检测方法和装置 |
| CN106568371A (zh) * | 2016-11-07 | 2017-04-19 | 北京航空航天大学 | 一种基于分形的轴向柱塞泵转子‑配流盘动态混合润滑建模方法 |
| CN111570534A (zh) * | 2020-06-18 | 2020-08-25 | 燕山大学 | 酸轧机组基于大数据理论的热轧特性对摩擦系数影响预测方法 |
| CN112974544A (zh) * | 2021-03-09 | 2021-06-18 | 燕山大学 | 一种y型轧机从动工作辊力矩动态调节装置及其方法 |
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004043858A (ja) * | 2002-07-10 | 2004-02-12 | Jfe Steel Kk | 連続溶融金属めっき装置のシンクロール |
| CN101201276A (zh) * | 2006-12-14 | 2008-06-18 | 宝山钢铁股份有限公司 | 沉没辊两端轴颈受力的检测方法和装置 |
| CN106568371A (zh) * | 2016-11-07 | 2017-04-19 | 北京航空航天大学 | 一种基于分形的轴向柱塞泵转子‑配流盘动态混合润滑建模方法 |
| CN111570534A (zh) * | 2020-06-18 | 2020-08-25 | 燕山大学 | 酸轧机组基于大数据理论的热轧特性对摩擦系数影响预测方法 |
| CN112974544A (zh) * | 2021-03-09 | 2021-06-18 | 燕山大学 | 一种y型轧机从动工作辊力矩动态调节装置及其方法 |
Non-Patent Citations (1)
| Title |
|---|
| 魏宝民等: "热轧卷取后的冷却工艺对低碳铝镇静钢热镀锌板耐时效性能的影响", 《燕山大学学报》, 25 May 2021 (2021-05-25) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116103597A (zh) * | 2023-03-15 | 2023-05-12 | 燕山大学 | 一种热镀铝锌过程锌锅内温度变化预报方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2017054427A1 (zh) | 连续变厚度带材的热镀锌层厚度控制系统及方法 | |
| CN106984651B (zh) | 一种提高轧件厚度控制精度的在线控制系统 | |
| CN101107378B (zh) | 热浸镀的热轧钢带的受控厚度减薄及其中所用的设备 | |
| CN100357953C (zh) | 金属板带自动矫直工艺参数的优化方法 | |
| CN114672753A (zh) | 一种热镀锌沉没辊系转动失效预报方法 | |
| CN110318013B (zh) | 一种热镀锌带钢的生产方法 | |
| CN110064667B (zh) | 一种钢板层流冷却方法 | |
| CN107876565A (zh) | 一种确定轧辊磨损量的方法及装置 | |
| CN102363159A (zh) | 一种单张板精密冷轧板厚测量系统的厚度控制方法 | |
| JP4919857B2 (ja) | 圧延方向に板厚がテーパ状に変化する片面テーパ鋼板の製造方法 | |
| CN108202085A (zh) | 冷连轧机的板厚控制装置 | |
| CN114912248A (zh) | 热镀锌机组运行过程中沉没辊系受力预报方法 | |
| JP2012171008A (ja) | 調質圧延方法および調質圧延装置、ならびに鋼帯の製造方法 | |
| CN114912248B (zh) | 热镀锌机组运行过程中沉没辊系受力预报方法 | |
| CN112052610A (zh) | 一种整辊式板形检测辊张力与包角设定的计算方法 | |
| CN108655188B (zh) | 一种板形测量误差补偿方法 | |
| CN113732063A (zh) | 用于压花辊轧带材的方法和装置 | |
| EP3369836B1 (en) | Method for manufacturing hot-dip galvanized steel sheet | |
| CN111382521A (zh) | 一种适用于连退机组炉内带钢残余应力的预报方法 | |
| CN110227722B (zh) | 基于轧辊换辊周期内轧辊辊型的优化方法 | |
| FI84791B (fi) | Foerfarande foer valsning av band i en valsanlaeggning samt dess styrsystem. | |
| CN108213088A (zh) | 一种冷轧钢板湿平整轧制过程中平整液流量的控制方法 | |
| CN107686958B (zh) | 一种连续热浸镀锌钢板成品厚度的精准控制方法 | |
| CN112916625A (zh) | 冷轧带钢板形控制的控制方法及控制系统 | |
| CN115480597B (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 |