X-type six-roller mill plate shape forecasting method

Abstract

The invention provides an X-type six-roller mill plate shape forecasting method, and belongs to the technical field of metallurgical rolling. The method comprises the following steps: collecting equipment parameters of the X-type six-roller mill, strip characteristic parameters and corresponding rolling process parameters; dividing units and solving influence coefficients; solving the plastic deformation of the strip; carrying out iterative calculation on the roll system elastic deformation model and the strip plastic deformation model until convergence; forecasting the transverse distribution value of the pre-tensile stress of the strip during rolling; and strip outlet plate shape transverse distribution during rolling is forecasted. According to the method, on the basis of the strip plastic deformation model and the roll system elastic deformation model, the pre-tensile stress after rolling is calculated, and then strip shape forecasting is achieved.

Description

Technical Field

[0001] The invention relates to the technical field of metallurgical rolling, and in particular to a method for predicting the plate shape of an X-type six-roller rolling mill. Background Art

[0002] With the continuous development of the automotive, aerospace and other industries, the quality requirements for copper strips are also increasing, and the demand for thin strip products is particularly prominent. The traditional four- and six-roll rolling mills have a large working diameter and poor thinning ability due to the limitation of the minimum rollable thickness. [1] In addition, the rigidity of the traditional rolling mill may not meet the requirements when the copper strip is cold rolled. Although the twenty-high rolling mill meets the requirements for producing high-quality copper strip in terms of rolling capacity and rigidity, it has a complex rolling mill structure. [2] The X-type six-roller mill has the disadvantages of simple structure, low cost, high production efficiency, and the rigidity and thinning capacity of the mill can better meet the requirements of producing copper thin strips. [3-4] In the actual production process of the X-type six-high rolling mill, the flatness problem is more obvious due to the thinness of the rolled strip. It is necessary to establish a mathematical model to quantitatively analyze the flatness problem and adjust the parameters to effectively control the flatness. By analyzing the mechanism of the rolling mill and establishing a flatness prediction model, the flatness can be predicted and feedback controlled. Although there are a lot of studies on the flatness prediction model at home and abroad, [5] However, there is no literature to establish a plate shape prediction model for X-type six-high rolling mill.

[0003] References

[0004] [1] Wang Dongcheng, Ren Yinjie, Lü Youchuang, et al. Rapid prediction model of flatness of Sendzimir 20-roll mill [J / OL]. Steel, 1-16.

[0005] [2] Zhang Qingdong, Dai Chang, Wen Jie, et al. Simulation study on the plate shape control performance of 20-high Sendzimir mill [J]. Steel Rolling, 2013, 30(03): 1-6.

[0006] [3] Chen Chuncan. Research on key technologies and equipment development of X-type 6-roll copper foil rolling mill[J]. Nonferrous Metal Processing, 2018, 47(01): 1-4+16.

[0007] [4]Cao Wenjuan. Research and development of copper foil working roll of X-type six-roll mill[J]. Heat Treatment, 2020, 35(06): 32-36.

[0008] [5] Wang Dongcheng, Zhang Wei, Wang Zhijie, et al. Rapid prediction method of strip shape of six-roll cold rolling mill [J]. Steel, 2016, 51(09): 51-56. Summary of the invention

[0009] According to the technical problem that the X-type six-roll mill has poor plate shape (straightness) during the production process, a plate shape prediction method for an X-type six-roll mill is provided. The present invention is a plate shape prediction method for an X-type six-roll mill. The present invention performs iterative calculations between the roll system elastic deformation model and the strip plastic deformation model until convergence, calculates the pre-tension stress after rolling, and then realizes plate shape prediction.

[0010] The technical means adopted by the present invention are as follows:

[0011] A method for predicting the shape of an X-type six-high rolling mill comprises the following steps:

[0012] Step 1: Obtain the basic equipment parameters of the X-type six-high rolling mill and the rolling process parameters of the strip;

[0013] Step 2: Divide the roller and the strip into units, a total of t units; the coordinate of the center point of each unit width is y i , i=1,2,…,t, each unit width is Δy i , i=1,2,…,t;

[0014] Step 3: Solve the support roller deflection influence coefficient α bij , working roll deflection influence coefficient matrix α wij And the working roll bending influence coefficient α Fwi ;

[0015] Step 4: Preset the initial value h of strip exit thickness distribution 1i ;

[0016] Step 5: Calculate the strip plastic deformation by using the strip element variation method in the three-dimensional plastic deformation model;

[0017] Step 6: Calculate the rolling force p i ;

[0018] Step 7: Preset the initial value of the pressure distribution between the rollers to q i ;

[0019] Step 8: Solve the roller flattening influence coefficient;

[0020] Step 9: List the deformation coordination equation between the rolls of the roll system;

[0021] Step 10: List the force and moment balance equations of the working roll;

[0022] Step 11: Establish a linear equation system to solve the roller pressure;

[0023] Step 12: Determine whether the calculated value q′ satisfies the inequality max|q′-q|<0.01; if not, set q=q+0.5(q′-q) and proceed to step 8; if satisfied, proceed to step 13;

[0024] Step 13: Calculate the transverse distribution of the strip thickness after rolling;

[0025] Step 14: Determine whether the calculated value h′ satisfies the inequality max|h′-h|<0.0001; if not, set h=h+0.001(h′-h), and proceed to step 5; if satisfied, proceed to step 15;

[0026] Step 15: Output the transverse distribution value σ of the strip pre-tension stress 1i ;

[0027] Step 16: According to the transverse distribution value σ of the strip pre-tension stress 1i , calculate the strip shape distribution BX at the current moment i ;

[0028]

[0029] Furthermore, in step 1, the basic equipment parameters of the X-type six-high rolling mill include: the roll length, roll radius, crown, elastic modulus and Poisson's ratio of the support roll and the working roll, which are respectively denoted as L b , L w , R b , R w , ΔD b , ΔD w 、E b 、E w ,υ b and w ; The distance between the support rollers on both sides is L s , the distance between the two bending roll hydraulic cylinders of the working roll is L Fw The angle between the line connecting the center points of the cross sections of the support roll and the working roll and the vertical direction is α.

[0030] Furthermore, in step 1, the rolling process parameters of the strip include: the average thickness of the strip material Average thickness of strip after rolling Width B, strip elastic modulus E s , Strip Poisson's ratio υ s , Strip deformation resistance σ s , strip shear deformation resistance k s , average post-tensioning stress of the strip Average pre-tension stress of strip Working roll bending force F w , calculate the allowable error ε q and ε h and the relaxation factor χ q and x h .

[0031] Furthermore, in step 2, the coordinate origin is selected at the left pressure support of the support roller; the roller and the strip are divided into units as a whole, and a total of t units are divided; the coordinate of the center point of the width of each unit is y i , i=1,2,…,t, each unit width is Δy i , i=1,2,…,t; where the starting unit number of the contact part between the backup roll and the working roll is t1, and the ending unit number is t2; the starting unit number of the contact part between the working roll and the strip is t3, and the ending unit number is t4;

[0032] When calculating the plastic deformation of the strip, the outlet lateral displacement of the strip unit node line needs to be taken as an unknown quantity, so the strip unit node lines are renumbered; the two node lines corresponding to unit t3 are numbered 0 and 1, and the two node lines corresponding to unit t4 are numbered n-1 and n, n=t4-t3+1, and the node line numbers of other units can be obtained by analogy.

[0033] Furthermore, in step 5, the strip plastic deformation is calculated by the strip element variation method in the three-dimensional plastic deformation model as follows:

[0034]

[0035] Among them, u0, u1, ...u n They represent the lateral displacement of the strip unit node line outlet, which is an unknown quantity; α1~α n , β1~β n They are related to unit width, strip elastic modulus, strip Poisson's ratio, strip shear deformation resistance, neutral point thickness, work roll elastic flattening radius, friction coefficient, average friction stress of contact surface in deformation zone, pre-rolling thickness, reduction, deformation zone length, and stress state coefficient respectively; k s Indicates deformation resistance; Among them, h 1i ,h 0i Indicates the outlet thickness and inlet thickness of the strip element, l 0i Indicates the length of incoming strip material. Indicates the lateral average length of the incoming material.

[0036] Furthermore, in step 8, the elastic flattening influence coefficient γ between the rollers is i , the equation of the roller flattening amount and roller pressure distribution is:

[0037] δ i =γ i q i ,(i=t1~t2);

[0038] Among them, δ i Indicates the elastic flattening between the rollers, q i Represents the roller pressure of the i-th unit.

[0039] Furthermore, in step 9, the inter-roller deformation coordination equation is:

[0040] f wi =2cosα(f bi +δ i +ΔD i ),(i=t1~t2);

[0041] Among them, f wi represents the deflection of the i-th unit of the working roll in the vertical direction; f bi It represents the deflection of the i-th unit of the support roll along the line connecting the center points of the support roll and the working roll. ΔD represents the no-load gap between the rolls, i.e., the crown. ΔD i =(ΔD wi +ΔD bi ) / 2.

[0042] Furthermore, in step 10, the force and moment balance equation of the working roll is:

[0043]

[0044] Among them, p i Represents the rolling pressure of the i-th unit.

[0045] Furthermore, in step 13, h 1i The calculation formula is:

[0046] h 1i =s0+2f wi +2δ wi +ΔD i ,(i=t3~t4);

[0047] Among them, δ wi Indicates the amount of elastic flattening of the working roll surface in contact with the workpiece.

[0048] Furthermore, in step 15, σ 1i The calculation formula is:

[0049]

[0050] Compared with the prior art, the present invention has the following advantages:

[0051] The present invention establishes a plate shape prediction model by analyzing the mechanism of the rolling mill, and can predict and feedback control the plate shape. Specifically, the plate shape can be predicted according to the rolling parameters, and the rolling parameters can be adjusted in time; the initial rolling parameters can be set according to the process parameters before rolling, and the adjustment amount of the parameters during subsequent rolling can be reduced; the rolling state of the X-type six-roll rolling mill can be simulated and designed, providing a reference for the model design, and also providing a reference for the roll type optimization design and other aspects. BRIEF DESCRIPTION OF THE DRAWINGS

[0052] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0053] Figure 1 It is a schematic diagram of the overall process of the present invention.

[0054] Figure 2 It is a schematic diagram of the segmented roll system of the X-type six-roller rolling mill of the present invention.

[0055] Figure 3 This is a diagram showing the distribution angles and numbering of the roll system of an X-type six-roller rolling mill of the present invention.

[0056] Figure 4 This is a diagram showing the correspondence between the unit node line number and the unit number of the present invention.

[0057] Figure 5 The lateral distribution diagrams of plate shapes of Examples 1 and 2 of the present invention are shown in FIG. 1 , where (a) is the lateral distribution diagram of plate shape of Example 1; (b) is the lateral distribution diagram of plate shape of Example 2. DETAILED DESCRIPTION

[0058] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0059] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0060] like Figure 1 As shown, the present invention provides a method for predicting the shape of an X-type six-high rolling mill, comprising the following steps:

[0061] Step 1: Obtain the basic equipment parameters of the X-type six-high mill, including: the roll length, roll radius, crown, elastic modulus and Poisson's ratio of the support roll and the working roll, respectively denoted as L b , L w , R b , R w , ΔD b , ΔD w 、E b 、E w ,υ b and w , the distance between the two sides of the support point is L s , the distance between the two bending roll hydraulic cylinders of the working roll is L Fw The angle between the line connecting the center points of the cross sections of the support roll and the working roll and the vertical direction is α.

[0062] At the same time, the rolling process parameters of the strip are obtained, including: the average thickness of the incoming strip Average thickness of strip after rolling Width B, strip elastic modulus E s , Strip Poisson's ratio v s , Strip deformation resistance σ s , strip shear deformation resistance k s , average post-tensioning stress of the strip Average pre-tension stress of strip Working roll bending force F w , calculate the allowable error ε q and ε h , relaxation factor χ q and x h .

[0063] Step 2: Unit division, such as Figure 2As shown. The origin of the coordinates is selected at the left support point of the support roller. The roller and the strip are divided into units, a total of t units. The coordinate of the center point of each unit width is y i (i=1,2,…,t), each unit width is Δy i (i=1,2,…,t). Among them, the starting unit number of the contact part between the support roll and the working roll is t1, and the ending unit number is t2; the starting unit number of the contact part between the working roll and the strip is t3, and the ending unit number is t4. When calculating the plastic deformation of the strip, the exit lateral displacement of the strip unit node line needs to be taken as an unknown quantity. Therefore, the strip unit node lines are renumbered, such as Figure 4 As shown. The two node lines corresponding to unit t3 are numbered 0 and 1, and the two node lines corresponding to unit t4 are numbered n-1 and n, (n=t4-t3+1). The node line numbers of other units can be obtained by analogy;

[0064] Step 3: Solve the deflection influence coefficient of each roll and the bending influence coefficient of the working roll, including: Calculate the deflection influence coefficient α of the working roll wij , Working roll bending influence coefficient α Fwi and the support roller deflection influence coefficient α bij ;

[0065] Step 4: Given the initial value h of the strip exit thickness distribution 1i ;

[0066] Step 5: Calculate the strip plastic deformation by the strip element variation method in the three-dimensional plastic deformation model; the formula is as follows:

[0067]

[0068] Among them, u0, u1, ...u n They represent the lateral displacement of the strip unit node line outlet, which is an unknown quantity; α1~α n , β1~β n They are related to unit width, strip elastic modulus, strip Poisson's ratio, strip shear deformation resistance, neutral point thickness, work roll elastic flattening radius, friction coefficient, average friction stress of contact surface in deformation zone, pre-rolling thickness, reduction, deformation zone length, and stress state coefficient respectively; k s Indicates deformation resistance; Among them, h 1i ,h 0i Indicates the outlet thickness and inlet thickness of the strip element, l 0i Indicates the length of incoming strip material. Indicates the lateral average length of the incoming material.

[0069] Step 6: Calculate the rolling force p i ;

[0070]

[0071] Among them, ω o Indicates the influence coefficient of various factors on rolling force except tension; It represents the influence coefficient of tension on rolling force; l represents the length of deformation zone.

[0072] Step 7: Given the initial value of the pressure distribution between the rollers as q i ;

[0073] Step 8: Solve the roller flattening influence coefficient and the roller elastic flattening influence coefficient γ i , the equation of the roller flattening amount and roller pressure distribution is:

[0074] δ i =γ i q i ,(i=t1~t2);

[0075] Among them, δ i Indicates the elastic flattening between the rollers, q i Represents the roller pressure of the i-th unit.

[0076] Step 9: List the deformation coordination equation between the rolls of the roll system. The deformation coordination equation between the rolls is:

[0077] f wi =2cosα(f bi +δ i +ΔD i ),(i=t1~t2);

[0078]

[0079] Among them, f w Indicates the deflection of the working roll in the vertical direction; f b Indicates the deflection of the support roll along the line connecting the center points of the support roll and the working roll cross section; Indicates the rigid displacement of the axis in the vertical direction of the working roll; C1 indicates the displacement of the axis at the left end of the working roll body; C2 indicates the displacement of the axis at the right end of the working roll body; ΔD indicates the no-load clearance (convexity) between the rolls, ΔD i =(ΔD wi +ΔD bi ) / 2.

[0080] Step 10: List the force and moment balance equations of the working roll. The force and moment balance equations of the working roll are:

[0081]

[0082] Among them, p i Represents the rolling pressure of the i-th unit.

[0083] Step 11: Establish a linear equation system to solve the roller pressure;

[0084] Step 12: Determine whether the inequality max|q′-q|<0.01 holds; if not, set q=q+0.5(q′-q) and proceed to step 8; if yes, proceed to step 13;

[0085] Step 13: Solve for the transverse distribution of strip thickness after rolling, h 1i The calculation formula is:

[0086] h 1i =s0+2f wi +2δ wi +ΔD i ,(i=t3~t4).

[0087] Among them, δ wi Indicates the elastic flattening amount of the working roll surface in contact with the workpiece;

[0088] Step 14: Determine whether the inequality max|h′-h|<0.0001 holds; if not, set h=h+0.001(h′-h) and proceed to step 5; if yes, proceed to step 15;

[0089] Step 15: Output the transverse distribution value σ of the strip pre-tension stress 1i , σ 1i The calculation formula is:

[0090]

[0091] Step 16: According to the transverse distribution value σ of the strip pre-tension stress 1i , calculate the strip shape distribution BX at the current moment i ;

[0092]

[0093] Example 1

[0094] Taking a certain X-type six-roller rolling mill as an example, according to Figure 1 The overall calculation flow chart of the X-type six-high rolling mill shape prediction method shown in the figure, first, in step 1, the basic equipment parameters of the X-type six-high rolling mill are collected: the radius R of the working roll body w =62.5mm, support roller radius R b =300mm, working roll elastic modulus E w =236GPa, elastic modulus of support roller E b =236GPa, Poisson's ratio of working roll w =0.3, Poisson's ratio of support roller b=0.3, the distance between the left and right pressing support points is L s =1000mm, distance between two bending roll hydraulic cylinders of working roll L Fw =1000mm, working roll body length L w =800mm, support roller body length L b =800mm, the angle between the line connecting the center points of the working roll and the support roll and the vertical direction is α=50°, and the crown of the working roll is ΔD w =0mm, the convexity of the support roller is ΔD b =0mm.

[0095] At the same time, in step 1, the rolling process parameters of typical specifications of products are collected, mainly including the average thickness of the incoming strip Average thickness of strip after rolling Width B=

[0096] 600mm, strip elastic modulus E s =194GPa, strip Poisson's ratio v s =0.3, strip deformation resistance σ s =300MPa, strip shear deformation resistance k s =173.1MPa, average post-tension stress of strip Average pre-tension stress of strip Working roll bending force F w =0KN, calculate the allowable error ε q =0.01, and ε h =0.0001, relaxation factor χ q =0.5 and χ h =0.001;

[0097] Then, in step 2, unit division is performed, and the strip and roll system unit division is as follows Figure 2 As shown. The origin of the coordinates is selected at the left support point of the support roller. The roller and the strip are divided into units, a total of 100 units. The coordinates of the center point of each unit width are

[0098] y i (i=1,2,…,100), each unit width is Δy i (i=1,2,…,100). The starting unit number of the contact part between the support roll and the working roll is 11, and the ending unit number is 90; the starting unit number of the contact part between the working roll and the strip is 21, and the ending unit number is 80. When calculating the plastic deformation of the strip, the exit lateral displacement of the strip unit node line needs to be taken as an unknown quantity. Therefore, the strip unit node lines are renumbered, such as Figure 3The two node lines corresponding to unit 21 are numbered 0 and 1, the two node lines corresponding to unit 80 are numbered 59 and 60, and the node line numbers of other units can be deduced by analogy;

[0099] Then, in step 3, the deflection influence coefficient of each roll and the bending influence coefficient of the working roll are solved, including: calculating the deflection influence coefficient α of the working roll wij 、Working roll bending influence coefficient α Fwi and the support roller deflection influence coefficient α bij ;

[0100] Then, in step 4, the initial value h of the strip exit thickness distribution is given 1i ;

[0101] Subsequently, in step 5, the strip plastic deformation is calculated by the strip element variation method in the three-dimensional plastic deformation model; the formula is as follows:

[0102]

[0103] Subsequently, in step 6, the rolling force p is calculated i ;

[0104]

[0105] Then, in step 7, the initial value of the pressure distribution between the rollers is given as q i ;

[0106] Then, in step 8, the roller flattening influence coefficient, roller elastic flattening influence coefficient γ i , the equation of the roller flattening amount and roller pressure distribution is:

[0107] δ i =γ i q i ,(i=11~90);

[0108] Then, in step 9, the deformation coordination equation between the rolls of the roll system is listed, and the deformation coordination equation between the rolls is:

[0109] f wi =2cosα(f bi +δ i +ΔD i ),(i=11~90);

[0110]

[0111] Then, in step 10, the force and moment balance equation of the working roll is listed, and the force and moment balance equation of the working roll is:

[0112]

[0113] Then, in step 11, a linear equation system is established to solve the pressure between the rollers;

[0114] Then, in step 12, determine whether the inequality max|q′-q|<0.01 holds; if not, set q=q+0.5(q′-q) and proceed to step 8; if so, proceed to step 13;

[0115] Then, in step 13, the lateral distribution of the strip thickness after rolling, h 1i The calculation formula is:

[0116] h 1i =s0+2f wi +2δ wi +ΔD i ,(i=21~80).

[0117] Then, in step 14, determine whether the inequality max|h′-h|<0.0001 holds; if not, set h=h+0.001(h′-h), and proceed to step 5; if so, proceed to step 15;

[0118] Then, in step 15, the transverse distribution value σ of the strip pre-tension stress is output 1i , σ 1i The calculation formula is:

[0119]

[0120] Then, in step 16, according to the transverse distribution value σ of the strip pre-tension stress 1i , calculate the strip shape distribution BX at the current moment i ,See Figure 5 (a);

[0121]

[0122] Example 2

[0123] Taking a certain X-type six-roller rolling mill as an example, according to Figure 1 The overall calculation flow chart of the X-type six-high rolling mill shape prediction method shown in the figure, first, in step 1, the basic equipment parameters of the X-type six-high rolling mill are collected: the radius R of the working roll body w =62.5mm, support roller radius R b =300mm, working roll elastic modulus E w =236GPa, elastic modulus of support roller E b =236GPa, Poisson's ratio of working roll w =0.3, Poisson's ratio of support roller u b =0.3, the distance between the left and right pressing points is L s=1000mm, distance between two bending roll hydraulic cylinders of working roll L Fw =1000mm, working roll body length L w =800mm, support roller body length L b =800mm, the angle between the line connecting the center points of the working roll and the support roll and the vertical direction is α=50°, and the crown of the working roll is ΔD w =0mm, the convexity of the support roller is ΔD b =0mm.

[0124] At the same time, in step 1, the rolling process parameters of typical specifications of products are collected, mainly including the average thickness of the incoming strip Average thickness of strip after rolling Width B=

[0125] 600mm, strip elastic modulus E s =194GPa, strip Poisson's ratio v s =0.3, strip deformation resistance σ s =300MPa, strip shear deformation resistance k s =173.1MPa, average post-tension stress of strip Average pre-tension stress of strip Working roll bending force F w =0KN, calculate the allowable error ε q =0.01, and ε h =0.0001, relaxation factor χ q =0.5 and χ h =0.001;

[0126] Then, in step 2, unit division is performed, and the strip and roll system unit division is as follows Figure 2 As shown. The origin of the coordinates is selected at the left support point of the support roller. The roller and the strip are divided into units, a total of 100 units. The coordinates of the center point of each unit width are

[0127] y i (i=1,2,…,100), each unit width is Δy i (i=1,2,…,100). The starting unit number of the contact part between the support roll and the working roll is 11, and the ending unit number is 90; the starting unit number of the contact part between the working roll and the strip is 21, and the ending unit number is 80. When calculating the plastic deformation of the strip, the exit lateral displacement of the strip unit node line needs to be taken as an unknown quantity. Therefore, the strip unit node lines are renumbered, such as Figure 3 The two node lines corresponding to unit 21 are numbered 0 and 1, the two node lines corresponding to unit 80 are numbered 59 and 60, and the node line numbers of other units can be deduced by analogy;

[0128] Then, in step 3, the deflection influence coefficient of each roll and the bending influence coefficient of the working roll are solved, including: calculating the deflection influence coefficient α of the working roll wij 、Working roll bending influence coefficient α Fwi and the support roller deflection influence coefficient α bij ;

[0129] Then, in step 4, the initial value h of the strip exit thickness distribution is given 1i ;

[0130] Subsequently, in step 5, the strip plastic deformation is calculated by the strip element variation method in the three-dimensional plastic deformation model; the formula is as follows:

[0131]

[0132] Subsequently, in step 6, the rolling force p is calculated i ;

[0133]

[0134] Then, in step 7, the initial value of the pressure distribution between the rollers is given as q i ;

[0135] Then, in step 8, the roller flattening influence coefficient, roller elastic flattening influence coefficient γ i , the equation of the roller flattening amount and roller pressure distribution is:

[0136] δ i =γ i q i ,(i=11~90);

[0137] Then, in step 9, the deformation coordination equation between the rolls of the roll system is listed, and the deformation coordination equation between the rolls is:

[0138] f wi =2cosα(f bi +δ i +ΔD i ),(i=11~90);

[0139]

[0140] Then, in step 10, the force and moment balance equation of the working roll is listed, and the force and moment balance equation of the working roll is:

[0141]

[0142] Then, in step 11, a linear equation system is established to solve the pressure between the rollers;

[0143] Then, in step 12, determine whether the inequality max|q′-q|<0.01 holds; if not, set q=q+0.5(q′-q) and proceed to step 8; if so, proceed to step 13;

[0144] Then, in step 13, the lateral distribution of the strip thickness after rolling, h 1i The calculation formula is:

[0145] h 1i =s0+2f wi +2δ wi +ΔD i ,(i=21~80).

[0146] Then, in step 14, determine whether the inequality max|h′-h|<0.0001 holds; if not, set h=h+0.001(h′-h), and proceed to step 5; if so, proceed to step 15;

[0147] Then, in step 15, the transverse distribution value σ of the strip pre-tension stress is output 1i , σ 1i The calculation formula is:

[0148]

[0149] Then, in step 16, according to the transverse distribution value σ of the strip pre-tension stress 1i , calculate the strip shape distribution BX at the current moment i ,See Figure 5 (b);

[0150]

[0151] The serial numbers of the above embodiments of the present invention are only for description and do not represent the advantages or disadvantages of the embodiments.

[0152] In the above embodiments of the present invention, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0153] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. Among them, the device embodiments described above are only schematic. For example, the division of the units can be a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of units or modules, which can be electrical or other forms.

[0154] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple units. Some or all of the units may be selected according to actual needs to achieve the purpose of the present embodiment.

[0155] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0156] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.

[0157] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for predicting the shape of an X-type six-high rolling mill, characterized in that: The following steps are involved: Step 1: Obtain the basic equipment parameters of the X-type six-high rolling mill and the rolling process parameters of the strip; Step 2: Divide the roller and the strip into units, a total of t units; the coordinate of the center point of each unit width is y i , i=1,2,…,t, each unit width is Δy i , i=1,2,…,t; Step 3: Solve the support roller deflection influence coefficient α bij , working roll deflection influence coefficient matrix α wij And the working roll bending influence coefficient α Fwi ; Step 4: Preset the initial value h of strip exit thickness distribution 1i ; Step 5: Calculate the strip plastic deformation by using the strip element variation method in the three-dimensional plastic deformation model; Step 6: Calculate the rolling force p i ; Step 7: Preset the initial value of the pressure distribution between the rollers to q i ; Step 8: Solve the roller flattening influence coefficient; Step 9: List the deformation coordination equation between the rolls of the roll system; Step 10: List the force and moment balance equations of the working roll; Step 11: Establish a linear equation system to solve the roller pressure; Step 12: Determine whether the calculated value q′ satisfies the inequality max|q′-q|<0.01; if not, set q=q+0.5(q′-q) and proceed to step 8; if satisfied, proceed to step 13; Step 13: Calculate the transverse distribution of the strip thickness after rolling; Step 14: Determine whether the calculated value h′ satisfies the inequality max|h′-h|<0.0001; if not, set h=h+0.001(h′-h), and proceed to step 5; if satisfied, proceed to step 15; Step 15: Output the transverse distribution value σ of the strip pre-tension stress 1i ; Step 16: According to the transverse distribution value σ of the strip pre-tension stress 1i , calculate the strip shape distribution BX at the current moment i ; 2. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 1, the basic equipment parameters of the X-type six-high rolling mill include: the roll length, roll radius, crown, elastic modulus and Poisson's ratio of the support roll and the working roll, respectively denoted as L b , L w , R b , R w , ΔD b , ΔD w 、E b 、E w ,υ b and w ; The distance between the support rollers on both sides is L s , the distance between the two bending roll hydraulic cylinders of the working roll is L Fw The angle between the line connecting the center points of the cross sections of the support roll and the working roll and the vertical direction is α.

3. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 1, the rolling process parameters of the strip include: the average thickness of the strip material Average thickness of strip after rolling Width B, strip elastic modulus E s , Strip Poisson's ratio v s , strip deformation resistance σ s , strip shear deformation resistance k s , average post-tensioning stress of the strip Average pre-tension stress of strip Working roll bending force F w , calculate the allowable error ε q and ε g and the relaxation factor χ q and x h .

4. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 2, the coordinate origin is selected at the left pressure support of the support roller; the roller and the strip are divided into units as a whole, and a total of t units are divided; the coordinate of the center point of the width of each unit is y i , i=1,2,…,t, each unit width is Δy i , i=1,2,…,t; where the starting unit number of the contact part between the backup roll and the working roll is t1, and the ending unit number is t2; the starting unit number of the contact part between the working roll and the strip is t3, and the ending unit number is t4; When calculating the plastic deformation of the strip, the outlet lateral displacement of the strip unit node line needs to be taken as an unknown quantity, so the strip unit node lines are renumbered; the two node lines corresponding to unit t3 are numbered 0 and 1, and the two node lines corresponding to unit t4 are numbered n-1 and n, n=t4-t3+1, and the node line numbers of other units can be obtained by analogy.

5. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 5, the strip plastic deformation is calculated by the strip element variation method in the three-dimensional plastic deformation model as follows: Among them, u0, u1, ...u n They represent the lateral displacement of the strip unit node line outlet, which is an unknown quantity; α1~α n , β1~β n They are related to unit width, strip elastic modulus, strip Poisson's ratio, strip shear deformation resistance, neutral point thickness, work roll elastic flattening radius, friction coefficient, average friction stress of contact surface in deformation zone, pre-rolling thickness, reduction, deformation zone length, and stress state coefficient respectively; k s Indicates deformation resistance; Among them, h 1i ,h 0i Indicates the outlet thickness and inlet thickness of the strip element, l 0i Indicates the length of incoming strip material. Indicates the lateral average length of the incoming material.

6. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 8, the elastic flattening influence coefficient γ between rollers i , the equation of the roller flattening amount and roller pressure distribution is: d i =c i q i ,(i=t1~t2); Among them, δ i Indicates the elastic flattening between the rollers, q i Represents the roller pressure of the i-th unit.

7. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 9, the inter-roller deformation coordination equation is: f wi =2cosα(f bi +d i +ΔD i ),(i=t1~t2); Among them, f wi represents the deflection of the i-th unit of the working roll in the vertical direction; f bi It represents the deflection of the i-th unit of the support roll along the line connecting the center points of the support roll and the working roll. ΔD represents the no-load gap between the rolls, i.e., the crown. ΔD i =(ΔD wi +ΔD bi ) / 2.

8. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 10, the force and moment balance equation of the working roll is: Among them, p i Represents the rolling pressure of the i-th unit.

9. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 13, h 1i The calculation formula is: h 1i =s0+2f wi +2δ wi +ΔD i ,(i=t3~t4); Among them, δ wi Indicates the amount of elastic flattening of the working roll surface in contact with the workpiece.

10. The method for predicting the shape of an X-type six-high rolling mill according to claim 1, characterized in that: In step 15, σ 1i The calculation formula is:

Images