CN113500100A - 基于轧制接触界面分段模型上力学参数的辊缝控制方法 - Google Patents
基于轧制接触界面分段模型上力学参数的辊缝控制方法 Download PDFInfo
- Publication number
- CN113500100A CN113500100A CN202110815211.0A CN202110815211A CN113500100A CN 113500100 A CN113500100 A CN 113500100A CN 202110815211 A CN202110815211 A CN 202110815211A CN 113500100 A CN113500100 A CN 113500100A
- Authority
- CN
- China
- Prior art keywords
- rolling
- rolled piece
- thickness
- contact interface
- force
- 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
Images
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
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Control Of Metal Rolling (AREA)
Abstract
本发明提供一种基于轧制接触界面分段模型上力学参数的辊缝控制方法,具体实施步骤为:S1、根据轧件厚度和变形区长度,将轧制变形区按照轧制接触界面摩擦力特征划分为不同变形区微元体;S2、根据轧制变形区的微元体,建立轧制接触界面内轧件金属流动的体积表达式;S3、根据轧制接触界面上接触弧长度l和轧件的平均厚度的比值大小,对轧制接触界面进行分段划分并计算各段内的单位轧制压力;S4、根据S3计算的单位轧制力,分别求得对应的总轧制力Pi;S5、将S4计算的总轧制力Pi的值与设定的允许偏差范围e值进行比较,进行动态调整。本发明可以对轧制中的轧制力进行准确计算,并及时发现现场生产中实际轧制力与设定轧制力的偏差,并做出调整。
Description
技术领域
本发明涉及钢铁冶金技术领域,特别涉及一种基于轧制界面分段模型上力学参数的辊缝控制方法。
背景技术
随着我国经济的飞速发展,传统产业向新型产业转型的进程加速,汽车制造、航空航天、精密仪器等行业的技术水平日益提高,各行各业对于高品质热轧板带材的需求也明显增加。在实际生产中,还存在着许多问题。比如因轧机振动而引起的板带材产品表面质量问题。连续、高速、重载和稳定轧制是轧钢生产追求的终极目标,但几乎所有的轧机设备在轧制生产的过程中都不能避免轧机稳定性和轧机振动问题。
在生产中,影响板带材表面质量的几种振动可以主要分为三类:轧机***垂直振动、水平振动和扭转振动。其中轧制界面受到全部三种振动的影响,是实现稳定轧制的关键所在。
轧制界面是轧辊和轧件的工作界面,轧制界面的动态行为对于轧制过程的振动有着至关重要的影响。轧制界面的动态行为与轧制界面的摩擦状态和轧制力密切相关。根据摩擦状态的不同,找到轧制力与轧机振动之间的关系,对于及时调整轧制参数并更好地生产高品质的板带材具有重要意义。
根据变形区长度与轧件的平均高度的比值,将轧制界面分为三类:第一类包含两段滑动区,两段制动区和一段停滞区;第二类包含两段滑动区和一段停滞区;第三类包含一段停滞区。
在轧制力、停滞区长度、动态速度模型中,存在互相关联耦合的参数,这使实现稳定轧制变得更加复杂。轧制变形区与轧制界面汇集了轧制过程动态力学参数信息。在轧制过程非稳态时,轧机振动会引起轧制变形区力学参数变化,力学参数的变化又将引起轧机运动状态发生动态变化。轧机动态运动和轧制变形区动态力学参数的相互耦合作用,决定了轧机***的动态特性行为特性,直接影响产品品质。因此,建立一种方法来准确计算相关力学参数并针对轧辊控制参数进行调整,对于轧制稳定性的提升有着重要意义和实用价值。
发明内容
针对现有技术存在的问题,本发明提供一种基于轧制界面分段模型上力学参数的辊缝控制方法,主要根据轧制界面的力学参数将其划分不同类型,并针对不同类型计算,从而可以进行更准确的计算,并针对计算结果提出轧辊压下量的控制变化量,从而提升轧制稳定性和产品质量。
本发明提供了一种基于轧制接触界面分段模型上力学参数的辊缝控制方法,具体实施步骤如下:
S1、根据轧件厚度和变形区长度,将轧制变形区按照轧制接触界面摩擦力特征划分为不同变形区微元体;
S2、根据轧制变形区的微元体,建立轧制接触界面内轧件金属流动的体积表达式,得到轧件在轧制接触界面上沿轧辊接触弧各点的水平速度;
S21、根据轧制理论,求得轧件沿轧制接触界面的中性角,具体表达式如下:
其中,μ为摩擦系数,Δh为轧件的压下量,R为轧辊的半径;
S22、根据轧制接触界面的几何分布状态特征,求得轧件的接触弧函数,具体表达式如下:
hx=2zx
其中,hB为轧件出口的厚度,R为轧辊半径,x为轧件和轧辊的接触界面上任意位置的横坐标;
S23、根据轧件的几何参数特征,求得轧制接触界面中性点处轧件的厚度,具体表达式如下:
其中,hB为轧后轧件的厚度,R为轧辊的半径,γ为轧制接触界面的动态中性角;
S31、根据轧制理论,将接触弧长度l和轧件的平均厚度的比值大于5 的设为第一类型;将接触弧长度l和轧件的平均厚度的比值在2至5之间的设为第二类型;将接触弧长度l和轧件的平均厚度的比值小于2的设为第三类型;
S32、计算第一种类型中轧件变形区各段的单位轧制压力;
S33、计算第二种类型中轧件变形区各段的单位轧制压力;
S34、计算第三种类型中轧件变形区各段的单位轧制压力;
S4、根据步骤S32至S34计算的单位轧制力,分别求得第一种类型、第二种类型和第三种类型对应的总轧制力Pi(i=1,2,3);
S41、根据步骤S32求得的单位轧制力,得到第一种类型的总轧制力P1;
S411、计算各分段区域点A-F对应的坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xF表示第一种类型中各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数;
S412、将步骤S32求得的各段轧制力进行求和,得到第一种轧制接触界面类型下的总轧制力P1,具体表达式如下:
其中,px表示任意段处的单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCE,PDF是制动区单位轧制力,PEF是停滞区的单位轧制力;
S42、根据步骤S33求得的单位轧制力,得到第二种轧制接触界面类型的总轧制力P2;
S421、计算各分段区域点A-F对应坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xD表示第一种类型中各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数;
S422、将步骤S33求得的各段轧制力进行求和,得到第二种轧制接触界面类型下的总轧制力P2,具体表达式如下:
其中,px表示任意段处的单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCD是停滞区单位轧制力;
S43、根据步骤S34求得的单位轧制力,由此得到第三种轧制接触界面类型的总轧制力P3;
S431、计算各分段区域点AB对应坐标,具体表达式如下:
xA=l
xB=0;
S432、将步骤S34求得的各段轧制力进行求和,得到第三种轧制接触界面类型下的总轧制力P3,具体表达式如下:
其中,hx表示任意位置处板带厚度,dhx表示任意位置出的厚度变化率,θ表示对应弦与轧件中心线的夹角,μ表示轧制界面摩擦系数,hnd为板带中性面厚度,k表示摩擦力达到最大值时的常值;
S5、将步骤S4计算的总轧制力Pi(i=1,2,3)的值与设定的允许偏差范围e 值进行比较,进行动态调整;
S51、若|P设-P|<e时,则不需要进行轧辊的控制;
S52、若|P设-P|>e且P设>P时,则调整轧制接触界面的辊缝值,如果设定的总的轧制力更大,则将轧件新的出口厚度值为h'B=95%hB,即控制轧辊压下,将h′B替换hB代入到上述步骤S3至S4中,重新计算,直至|P设-P|<e为止;如果计算的总轧制力更大,则将轧件新的出口厚度值增加5%,即控制轧辊抬升,将h′B替换hB代入到上述步骤S3至S4中,重新计算,直至|P设-P|<e为止。此时的hB值为最终轧辊应调整到的最终值。
可优选的是,所述步骤S32具体包括以下步骤:
S322、根据步骤S321中变形区的长度,求得轧件变形区各段的单位轧制力:
前滑区BD区间单位轧制力的具体表达式如下:
后滑区AC区间单位轧制力的具体表达式如下:
制动区CE和DF区间单位轧制力的具体表达式分别如下:
其中,k表示摩擦力达到最大值时的常值,hC,hD对应制动区分段C、D位置处的轧件厚度,hx表示对应分段区间内d任意位置处的板带厚度,δ表示相对应段的相关系数,μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角;
停滞区EF区间单位轧制力的具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,θ表示对应弦与轧件中心线的夹角,hnd为板带中性面厚度,hx为对应停滞区任意位置轧件厚度,pE表示E点应力值,hE表示E点对应板带厚度。
可优选的是,所述步骤S33具体包括以下步骤:
S332、根据步骤S331中变形区的长度,求得轧件变形区各段的单位轧制力:
前滑区BD区间单位轧制压力具体表达式如下:
后滑区AC区间单位轧制压力具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角,hx为对应后滑区任意位置轧件厚度,hA表示A点对应轧制界面出口处的板带厚度;
停滞区CD区间单位轧制压力具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度,hC表示C点对应板带厚度,pC表示C点应力值。
可优选的是,所述第一种类型中停滞区EF区间接触弧长度lEF的具体表达式为:
其中,hE为x等于l/2时轧件的厚度,θEN为停滞区EF区间内任意位置接触弧与水平面的夹角。
可优选的是,所述第二种类型中停滞区CD区间接触弧长度lCD的具体表达式为:
其中,l表示变形区长度,k表示摩擦力达到最大值时的常值;δ表示相对应段的相关系数μ表示摩擦系数;θ表示对应弦与轧件中心线的夹角,hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度,hA表示A点对应板带厚度,Δh表示轧件压下量。
可优选的是,所述步骤S412中各段轧制力的具体表达式为:
x∈(xF,xE)
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度,hA~hE表示A~E点对应的板带厚度,l表示变形区长度,
可优选的是,所述步骤S422中各段轧制力的具体表达式为:
其中,xA~xD表示第一种类型中各点的坐标计算值,hA,hB,hD表示A,B, D点对应的板带厚度,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角,hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度。
本发明与现有技术相比,具有如下优点:
1.本发明可以对轧制中的轧制力进行准确计算,并及时发现现场生产中实际轧制力与设定轧制力发生偏差情况,做出轧辊控制调整。
2.本发明通过轧制力稳定控制,改善因轧制力不当而产生的轧制过程板带振动问题,从而提高板带材质量,延长轧机设备寿命,节约成本并相应地减少不必要的产能浪费。
3.本发明将理论研究和实际生产相结合,是产学研体系的一种代表形式,使生产更高效,研究更具针对性和实用性。
附图说明
图1为本发明基于轧制界面分段模型上力学参数的辊缝控制方法中第一种类型下单位摩擦力沿接触弧的分布图;
图2为本发明基于轧制界面分段模型上力学参数的辊缝控制方法中第二种类型下单位摩擦力沿接触弧的分布图;
图3为本发明基于轧制界面分段模型上力学参数的辊缝控制方法中第三种类型下单位摩擦力沿接触弧的分布图;
图4为本发明基于轧制界面分段模型上力学参数的辊缝控制方法的流程图。
具体实施方式
为详尽本发明之技术内容、所达成目的及功效,以下将结合说明书附图进行详细说明。
如图4所示,本发明主要是基于轧制接触界面分段模型上的力学参数进行计算,求解动态轧制接触界面动态接触参数和单位轧制压力,并针对计算结果提出轧辊压下量的控制变化量,从而提升轧制稳定性和产品质量。在本发明的一个优选实施方式中,本发明的基于轧制接触界面分段模型上力学参数的辊缝控制方法,包括以下步骤:
S1、根据轧件厚度和变形区长度,将轧制变形区按照轧制接触界面摩擦力特征划分为不同变形区微元体。
S2、根据轧制变形区的微元体,建立轧制接触界面内轧件金属流动的体积表达式,得到轧件在轧制接触界面上沿轧辊接触弧各点的水平速度。
S4、根据步骤S32至S34计算的单位轧制力,分别求得第一种类型、第二种类型和第三种类型对应的总轧制力Pi(i=1,2,3)。
S5、将步骤S4计算的总轧制力Pi(i=1,2,3)的值与设定的允许偏差范围 e值进行比较,进行动态调整,e根据使用环境与不同用途确定。
进一步的,步骤S2中的获得轧件在轧制接触界面上沿轧辊接触弧各点的水平速度的方法包括:
S21、根据轧制理论,求得轧件沿轧制接触界面的中性角,具体表达式如下:
其中,μ为摩擦系数,Δh为轧件的压下量,R为轧辊的半径。
S22、根据轧制接触界面的几何分布状态特征,求得轧件的接触弧函数,具体表达式如下:
hx=2zx
其中,hB为轧件出口的厚度,R为轧辊半径,x为轧件和轧辊的接触界面上任意位置的横坐标。
S23、根据轧件的几何参数特征,求得轧制接触界面中性点处轧件的厚度,具体表达式如下:
其中,hB为轧后轧件的厚度,R为轧辊的半径,γ为轧制接触界面的动态中性角。
进一步的,步骤S3中的获得轧制接触界面各段内的单位轧制压力的方法包括,
S31、根据轧制理论,将接触弧长度l和轧件的平均厚度的比值大于5 的设为第一类型;将接触弧长度l和轧件的平均厚度的比值在2至5之间的设为第二类型;将接触弧长度l和轧件的平均厚度的比值小于2的设为第三类型。
进一步的,为了更好的描述第一种类型,如图1所示,设A点和B点处的厚度分别为轧件的入口厚度和轧件的出口厚度;AC段和BD段为滑动区, CE段和DF段为制动区,EF段为停滞区。
如图2所示,设A点和B点处的厚度分别为轧件的入口厚度和轧件的出口厚度;AC段和BD段为滑动区,CD段为停滞区。
如图3所示,设A点和B点处的厚度分别为轧件的入口厚度和轧件的出口厚度;AB段为停滞区。
S32、计算第一种类型中轧件变形区各段的单位轧制压力:
ˉ
S321、根据轧制接触界面上,接触弧长度l和轧件的平均厚度h的比值大于5时轧件与轧辊接触摩擦力的状态,将轧件变形区分为前滑区、后滑区、制动区和停滞区。
S322、根据步骤S321中变形区的长度,求得轧件变形区各段的单位轧制力:
前滑区BD区间单位轧制力的具体表达式如下:
后滑区AC区间单位轧制力的具体表达式如下:
制动区CE和DF区间单位轧制力的具体表达式分别如下:
停滞区EF区间单位轧制力的具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx为对应停滞区任意位置轧件厚度,pE表示轧制界面E点应力值,hE表示E点对应板带厚度,hnd表示轧制接触界面中性点处的轧件厚度。
S33、计算第二种类型中轧件变形区各段的单位轧制压力:
S332、根据步骤S331中变形区的长度,求得轧件变形区各段的单位轧制力:
前滑区BD区间单位轧制压力具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx为对应前滑区任意位置轧件厚度,hB表示轧制界面入口位置B点处对应板带厚度;
后滑区AC区间单位轧制压力具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx为对应后滑区任意位置轧件厚度,hA表示轧制界面出口位置处A点对应板带厚度;
停滞区CD区间单位轧制压力具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx对应停滞区任意位置轧件厚度,pC表示C点应力值,hnd表示轧制接触界面板带中性面厚度,hC表示C点对应板带厚度,pC表示轧制界面C点应力值。
S34、计算第三种类型中轧件变形区各段的单位轧制压力:
其中,l表示变形区长度,k表示摩擦力达到最大值时的常值;δ表示相对应段的相关系数μ表示摩擦系数;θ表示对应弦与轧件中心线的夹角,hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度,hA表示轧制界面出口位置处A点对应板带厚度。
具体而言,第一种类型中停滞区EF区间接触弧长度lEF的具体表达式为:
其中,hE为x等于l/2时轧件的厚度,θEN为停滞区任意位置接触弧与水平面的夹角。
具体而言,第二种类型中停滞区CD区间接触弧长度lCD的具体表达式为:
进一步的,步骤S4中的求得每种类型对应的总轧制力Pi(i=1,2,3)的具体实施过程如下:
S41、根据步骤S32求得的单位轧制力,得到第一种类型的总轧制力P1;
S411、计算各分段区域点A-F对应的坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xF表示第一种类型中各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数。
S412、将步骤S32求得的各段轧制力进行求和,得到第一种轧制接触界面类型下的总轧制力P1,具体表达式如下:
其中,px表示任意分段单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCE,PDF是制动区单位轧制力,PEF是停滞区的单位轧制力。
进一步的,各段轧制力的具体表达式为:
x∈(xF,xE)
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度,hA~ hE表示A~E点对应的板带厚度,l表示变形区长度,
S42、根据步骤S33求得的单位轧制力,得到第二种轧制接触界面类型的总轧制力P2:
S421、计算各分段区域点A-D对应坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xD表示第一种类型中各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数。
S422、将步骤S33求得的各段轧制力进行求和,得到第二种轧制接触界面类型下的总轧制力P2,具体表达式如下:
其中,px表示任意分段单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCD是停滞区单位轧制力,PEF是停滞区的单位轧制力。
进一步的,S422中各段轧制力的具体表达式为:
其中,xA~xD表示第一种类型中各点的坐标计算值,hA,hB,hD表示A, B,D点对应板带厚度,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角;hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度。
S43、根据步骤S34求得的单位轧制力,由此得到第三种轧制接触界面类型的总轧制力P3:
S431、计算各分段区域点AB对应坐标,具体表达式如下:
xA=l
xB=0
S432、将步骤S34求得的各段轧制力进行求和,得到第三种轧制接触界面类型下的总轧制力P3,具体表达式如下:
其中,hx表示任意位置处板带厚度,dhx表示任意位置出的厚度变化率,θ表示对应弦与轧件中心线的夹角,μ表示轧制界面摩擦系数,hnd为板带中性面厚度,k表示摩擦力达到最大值时的常值。
进一步的,步骤S5中的动态调整,具体过程如下:
S51、若|P设-P|<e时,则不需要进行轧辊的控制;
S52、若|P设-P|>e且P设>P时,则调整轧制接触界面的辊缝值,如果设定的总的轧制力更大,则将轧件新的出口厚度值为h'B=95%hB,即控制轧辊压下,将h′B替换hB代入到上述步骤S3至S4中,重新计算,直至|P设-P|<e为止;如果计算的总轧制力更大,则将轧件新的出口厚度值增加5%,即控制轧辊抬升,将h′B替换hB代入到上述步骤S3至S4中,重新计算,直至|P设-P|<e为止。此时的hB值为最终轧辊应调整到的最终值。
以下结合实施例对本发明一种基于轧制界面分段模型上力学参数的辊缝控制方法做进一步描述:
确定输入的轧制参数,如表1所示。
表1轧制参数值
具体实施过程如下:
S1、根据本实施例中轧件的厚度和变形区的长度,将轧制变形区按照轧制接触界面摩擦力特征划分为不同变形区微元体。
S2、根据轧制变形区的微元体,建立轧制接触界面内轧件金属流动的体积表达式,得到轧件在轧制接触界面上沿轧辊接触弧各点的水平速度。
S3、根据出口厚度和入口厚度的值,取两者的平均值得到轧件的平均厚度,即(hA+hB)/2=11.5mm;轧件的变形区长度为:根据求得的轧制接触界面上接触弧长度l和轧件的平均厚度的比值大小,对轧制接触界面进行分段划分并计算各段内的单位轧制压力。
S33、计算第二种类型中轧件变形区各段的单位轧制压力:
S332、根据步骤S331中变形区的长度,求得轧件变形区各段的单位轧制力:
前滑区BD区间单位轧制压力具体表达式如下:
后滑区AC区间单位轧制压力具体表达式如下:
停滞区CD区间单位轧制压力具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角, hx对应停滞区任意位置轧件厚度,pC表示C点应力值,hnd表示轧制接触界面板带中性面厚度,hC表示C点对应板带厚度。
S4、根据步骤S32至S34计算的单位轧制力,分别求得第一种类型、第二种类型和第三种类型对应的总轧制力Pi(i=1,2,3)。
S42、根据步骤S33求得的单位轧制力,得到第二种轧制接触界面类型的总轧制力P2;
S421、计算各分段区域点A-D对应坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xD表示第一种类型中各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数。
S422、将步骤S33求得的各段轧制力进行求和,得到第二种轧制接触界面类型下的总轧制力P2,具体表达式如下:
其中,px表示任意分段单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCD是停滞区单位轧制力,PEF是停滞区的单位轧制力。
其中,xA~xD表示第一种类型中各点的坐标计算值,hA,hB,hD表示A, B,D点对应板带厚度,k表示轧制接触界面中摩擦力达到最大值时的常值,δ表示相对应段的相关系数μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角,hx对应停滞区任意位置轧件厚度,hnd表示轧制接触界面板带中性面厚度。
S5、进行判断,总轧制力的设定值与计算值差值是否在允许范围e(本实施例子中e取5%P)内。判断|P设-P|<5%P设是否成立,若成立,则流程结束,轧辊无需调整;若不成立,判断P设>P是否成立,若成立,将hB2=95%hB作为新的hB代入步骤一并重新进行计算直至结束;若P设>P不成立,将hB2=105%hB 作为新的hB代入步骤一并重新进行计算直至结束。结束时hB2的值为轧件出口厚度值,(hB2-hB1)为轧辊的调整量。
通过上述实施例,本发明通过对轧制力求解模型进行优化求解,得到了精度更高的轧制力;基于轧制力反馈控制理论,对轧制力进行修正,保证轧制过程轧制界面受载稳定性。本发明可以保证轧制过程工艺参数的稳定性,进而提高板带生产精度和生产效率。
以上所述的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (8)
1.一种基于轧制接触界面分段模型上力学参数的辊缝控制方法,其特征在于,具体实施步骤如下:
S1、根据轧件厚度和变形区长度,将轧制变形区按照轧制接触界面摩擦力特征划分为不同变形区微元体;
S2、根据轧制变形区的微元体,建立轧制接触界面内轧件金属流动的体积表达式,得到轧件在轧制接触界面上沿轧辊接触弧各点的水平速度;
S21、求得轧件沿轧制接触界面的中性角,具体表达式如下:
其中,μ为轧制界面摩擦系数,Δh为轧件的压下量,R为轧辊的半径;
S22、根据轧制接触界面的几何分布状态特征,获得轧件的接触弧函数,具体表达式如下:
hx=2zx
其中,zx表示接触弧函数,dhx表示任意位置出的厚度变化率,hB为轧件出口的厚度,R为轧辊半径,x为轧件和轧辊的接触界面上任意位置的横坐标;
S23、根据轧件的几何参数特征,求得轧制接触界面中性点处轧件的厚度,具体表达式如下:
其中,hB为轧后轧件的厚度,R为轧辊的半径,γ为轧制接触界面的动态中性角;
S31、根据轧制理论,将接触弧长度l和轧件的平均厚度的比值大于5的设为第一类型;将接触弧长度l和轧件的平均厚度的比值在2至5之间的设为第二类型;将接触弧长度l和轧件的平均厚度的比值小于2的设为第三类型;
S32、计算第一种类型中轧件变形区各段的单位轧制压力;
S33、计算第二种类型中轧件变形区各段的单位轧制压力;
S34、计算第三种类型中轧件变形区各段的单位轧制压力;
S4、根据步骤S32至S34计算的单位轧制力,分别求得第一种类型、第二种类型和第三种类型对应的总轧制力Pi(i=1,2,3);
S41、根据步骤S32求得的单位轧制力,得到第一种类型的总轧制力P1;
S411、计算各分段区域点A-F对应的坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xF表示第一种类型中各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数;
S412、将步骤S32求得的各段轧制力进行求和,得到第一种轧制接触界面类型下的总轧制力P1,具体表达式如下:
其中,px表示任意分段内的单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCE,PDF是制动区单位轧制力,PEF是停滞区的单位轧制力;
S42、根据步骤S33求得的单位轧制力,得到第二种轧制接触界面类型的总轧制力P2;
S421、计算各分段区域点A-F对应坐标,具体表达式如下:
xB=0
xA=l
其中,xA~xD表示第一种类型中A~D各点的坐标计算值,l为接触弧长,R为工作辊半径,Δh表示板带压下量,μ表示轧制界面摩擦系数;
S422、将步骤S33求得的各段轧制力进行求和,得到第二种轧制接触界面类型下的总轧制力P2,具体表达式如下:
其中,px表示任意分段单位轧制力,PAC是前滑区的单位轧制力,PBD是后滑区的单位轧制力,PCD是停滞区单位轧制力,PEF是停滞区的单位轧制力;
S43、根据步骤S34求得的单位轧制力,由此得到第三种轧制接触界面类型的总轧制力P3;
S431、计算各分段区域点AB对应坐标,具体表达式如下:
xA=l
xB=0;
S432、将步骤S34求得的各段轧制力进行求和,得到第三种轧制接触界面类型下的总轧制力P3,具体表达式如下:
其中,hx表示任意位置处板带厚度,dhx表示任意位置出的厚度变化率,θ表示对应弦与轧件中心线的夹角,μ表示轧制界面摩擦系数,hnd为板带中性面厚度,k表示摩擦力达到最大值时的常值;
S5、将步骤S4计算的总轧制力Pi(i=1,2,3)的值与设定的允许偏差范围e值进行比较,进行动态调整;
S51、若|P设-P|<e时,则不需要进行轧辊的控制;
S52、若|P设-P|>e且P设>P时,则调整轧制接触界面的辊缝值,如果设定的总的轧制力更大,则将轧件新的出口厚度值为h'B=95%hB,即控制轧辊压下,将h′B替换hB代入到上述步骤S3至S4中,重新计算,直至|P设-P|<e为止;如果计算的总轧制力更大,则将轧件新的出口厚度值增加5%,即控制轧辊抬升,将h′B替换hB代入到上述步骤S3至S4中,重新计算,直至|P设-P|<e为止,此时的hB值为最终轧辊应调整到的最终值。
2.根据权利要求1所述的基于轧制接触界面分段模型上力学参数的辊缝控制方法,其特征在于,所述步骤S32具体包括以下步骤:
S322、根据步骤S321中变形区的长度,求得轧件变形区各段的单位轧制力:
前滑区BD区间单位轧制力的具体表达式如下:
后滑区AC区间单位轧制力的具体表达式如下:
制动区CE和DF区间单位轧制力的具体表达式分别如下:
其中,k表示摩擦力达到最大值时的常值,hC,hD对应制动区分界点C、D位置处的轧件厚度,hx表示制动区任意位置处板带厚度,δ表示相对应段的相关系数,μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角,hnd表示轧制接触界面板带中性面厚度;
停滞区EF区间单位轧制力的具体表达式如下:
其中,k表示轧制接触界面中摩擦力达到最大值时的常值,μ表示摩擦系数,θ表示对应弦与轧件中心线的夹角,hx为对应停滞区任意位置轧件厚度,pE表示E点对应轧制应力值,hE表示E点对应板带厚度,hnd表示轧制接触界面板带中性面厚度。
3.根据权利要求1所述的基于轧制接触界面分段模型上力学参数的辊缝控制方法,其特征在于,所述步骤S33具体包括以下步骤:
S332、根据步骤S331中变形区的长度,求得轧件变形区各段的单位轧制力;
前滑区BD区间单位轧制压力具体表达式如下:
后滑区AC区间单位轧制压力具体表达式如下:
停滞区CD区间单位轧制压力具体表达式如下:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110815211.0A CN113500100B (zh) | 2021-07-19 | 2021-07-19 | 基于轧制接触界面分段模型上力学参数的辊缝控制方法 |
JP2022091125A JP7212976B1 (ja) | 2021-07-19 | 2022-06-03 | 圧延接触界面セグメント化モデルの力学パラメータに基づくロールギャップ制御方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110815211.0A CN113500100B (zh) | 2021-07-19 | 2021-07-19 | 基于轧制接触界面分段模型上力学参数的辊缝控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113500100A true CN113500100A (zh) | 2021-10-15 |
CN113500100B CN113500100B (zh) | 2022-04-26 |
Family
ID=78013850
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110815211.0A Active CN113500100B (zh) | 2021-07-19 | 2021-07-19 | 基于轧制接触界面分段模型上力学参数的辊缝控制方法 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7212976B1 (zh) |
CN (1) | CN113500100B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115990624B (zh) * | 2023-03-22 | 2023-05-30 | 太原理工大学 | 一种差厚板增厚制备过程中轧制力的预测方法 |
CN116984385B (zh) * | 2023-09-26 | 2024-01-05 | 太原理工大学 | 冷连轧带钢由薄到厚规格轧制中轧制力确定方法及装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5609054A (en) * | 1991-09-10 | 1997-03-11 | Nippon Steel Corporation | Rolling mill for flat products |
JPH11156412A (ja) * | 1997-12-01 | 1999-06-15 | Nkk Corp | 液圧圧下式圧延機 |
CN101566538A (zh) * | 2009-05-27 | 2009-10-28 | 东北大学 | 一种中厚板轧制过程中轧件塑性系数在线获取方法 |
CN101633003A (zh) * | 2009-07-07 | 2010-01-27 | 东北大学 | 周期性变厚度带材轧制过程中厚度的控制方法及控制*** |
CN101797588A (zh) * | 2010-04-01 | 2010-08-11 | 中色科技股份有限公司 | 一种热轧机二级控制轧制力预报方法 |
CN102039317A (zh) * | 2009-10-10 | 2011-05-04 | 北京理工大学 | 一种中厚板轧机液压压下***压力超限保护方法 |
CN102581035A (zh) * | 2012-01-30 | 2012-07-18 | 中冶南方工程技术有限公司 | 一种冷轧带钢板形前馈控制*** |
CN103008360A (zh) * | 2012-10-30 | 2013-04-03 | 中冶南方(武汉)信息技术工程有限公司 | 一种确定冷轧机工作辊温度场及热膨胀的方法 |
CN105195524A (zh) * | 2015-10-21 | 2015-12-30 | 东北大学 | 一种冷轧轧制升降速过程中板带厚度补偿控制方法 |
CN107138537A (zh) * | 2017-06-08 | 2017-09-08 | 北京科技大学 | 一种采用电动和液压压下的铝板热轧机厚度控制方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS609511A (ja) * | 1983-06-29 | 1985-01-18 | Hitachi Ltd | 圧延機における寸法制御方法 |
JPH05337527A (ja) * | 1992-06-08 | 1993-12-21 | Nkk Corp | 鋼ストリップの圧延方法 |
FR2783444B1 (fr) * | 1998-09-21 | 2000-12-15 | Kvaerner Metals Clecim | Procede de laminage d'un produit metallique |
CN104001739B (zh) * | 2014-05-14 | 2015-12-09 | 中冶南方工程技术有限公司 | 一种冷轧平整轧制压力的获取方法及装置 |
CN104998913B (zh) * | 2015-06-29 | 2017-03-15 | 东北大学 | 一种冷轧轧制过程中电机功率的预测方法 |
CN106623443B (zh) * | 2016-12-28 | 2019-01-22 | 苏州大学 | 一种厚板粗轧阶段轧制力预报方法 |
-
2021
- 2021-07-19 CN CN202110815211.0A patent/CN113500100B/zh active Active
-
2022
- 2022-06-03 JP JP2022091125A patent/JP7212976B1/ja active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5609054A (en) * | 1991-09-10 | 1997-03-11 | Nippon Steel Corporation | Rolling mill for flat products |
JPH11156412A (ja) * | 1997-12-01 | 1999-06-15 | Nkk Corp | 液圧圧下式圧延機 |
CN101566538A (zh) * | 2009-05-27 | 2009-10-28 | 东北大学 | 一种中厚板轧制过程中轧件塑性系数在线获取方法 |
CN101633003A (zh) * | 2009-07-07 | 2010-01-27 | 东北大学 | 周期性变厚度带材轧制过程中厚度的控制方法及控制*** |
CN102039317A (zh) * | 2009-10-10 | 2011-05-04 | 北京理工大学 | 一种中厚板轧机液压压下***压力超限保护方法 |
CN101797588A (zh) * | 2010-04-01 | 2010-08-11 | 中色科技股份有限公司 | 一种热轧机二级控制轧制力预报方法 |
CN102581035A (zh) * | 2012-01-30 | 2012-07-18 | 中冶南方工程技术有限公司 | 一种冷轧带钢板形前馈控制*** |
CN103008360A (zh) * | 2012-10-30 | 2013-04-03 | 中冶南方(武汉)信息技术工程有限公司 | 一种确定冷轧机工作辊温度场及热膨胀的方法 |
CN105195524A (zh) * | 2015-10-21 | 2015-12-30 | 东北大学 | 一种冷轧轧制升降速过程中板带厚度补偿控制方法 |
CN107138537A (zh) * | 2017-06-08 | 2017-09-08 | 北京科技大学 | 一种采用电动和液压压下的铝板热轧机厚度控制方法 |
Non-Patent Citations (2)
Title |
---|
李毅: "柔性轧制成形区对三维曲面件成形结果的影响", 《柔性轧制成形区对三维曲面件成形结果的影响》 * |
黄红乾: "半固态合金强流变轧制单位轧制压力分布的数学模型和理论研究", 《半固态合金强流变轧制单位轧制压力分布的数学模型和理论研究》 * |
Also Published As
Publication number | Publication date |
---|---|
JP2023017105A (ja) | 2023-02-03 |
JP7212976B1 (ja) | 2023-01-26 |
CN113500100B (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113500100B (zh) | 基于轧制接触界面分段模型上力学参数的辊缝控制方法 | |
CN101537431A (zh) | 冷轧薄带钢表面粗糙度在线预测和控制方法 | |
CN102451838B (zh) | 一种消除钢板热轧过程镰刀弯缺陷的方法 | |
CN101003063A (zh) | 在平整轧制过程中带钢的延伸率、板形及综合控制方法 | |
CN110202001A (zh) | 超薄超宽不锈钢产品轧制的辊系配备方法及轧制方法 | |
CN106077098B (zh) | 一种双锥度工作辊及其辊形设计方法 | |
CN111832155B (zh) | 一种考虑中间坯凸度的热轧粗轧负荷分配方法 | |
KR890003802B1 (ko) | 형상제어 조작량 설정값의 결정장치 | |
CN101422785B (zh) | 一种调整二次冷轧荫罩带钢机组中间辊窜动的方法 | |
CN110860561A (zh) | 短流程冷、热轧协调控制电工钢平直截面的方法 | |
CN113751510A (zh) | 一种带钢的轧制控制方法 | |
CN113361179A (zh) | 一种中厚板轧机角轧满负荷规程分配方法 | |
JP4948301B2 (ja) | 冷間圧延における形状制御方法 | |
CN110227721B (zh) | 一种用于变速过程中轧机辊缝的控制方法 | |
CN112958633A (zh) | 基于来料镰刀弯的精轧带钢头部预摆调平控制方法 | |
JP4671544B2 (ja) | 冷間タンデム圧延における板幅制御方法 | |
CN111906152B (zh) | 一种控制中厚板板形的精轧阶段负荷调整方法 | |
JPH06254611A (ja) | 形状クラウン制御熱間圧延機および圧延方法 | |
Wang et al. | Modeling and optimization of threading process for shape control in tandem cold rolling | |
CN114918258A (zh) | 一种冷连轧机基于大数据的来料断面形状窜辊补偿方法 | |
JP4813014B2 (ja) | 冷間タンデム圧延機の形状制御方法 | |
CN114505352B (zh) | 一种同板差窜辊前馈动态控制方法 | |
CN107999546A (zh) | 基于辊端压靠治理极薄带钢的边部碎浪的方法 | |
Ogarkov et al. | Analysis and mathematical modeling for microgeometry reproduction of a textured roll surface on a cold-rolled galvanized strip during temper rolling | |
CN116371943A (zh) | 一种ds四辊轧机的有载辊缝计算方法 |
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 |