Plate shape regulation and control efficiency coefficient prediction method in cold-rolled plate shape control

Abstract

The invention provides a strip shape regulation and control efficacy coefficient prediction method in cold-rolled strip shape control, and belongs to the technical field of metallurgical rolling. The method comprises the following steps: firstly, calculating a plate shape regulation efficiency coefficient matrix of different preset working condition points through off-line simulation, improving an optimal weight factor between the preset working condition points and an actual working condition point by utilizing production data, and then establishing an optimal weight factor prediction model of the actual working condition point by adopting a data driving mode. And the plate shape regulation and control efficiency coefficient under any rolling working condition is accurately obtained. The method for predicting the optimal weight factor of the actual working condition point not only has relatively high precision, but also has relatively good stability and relatively high convergence speed, and a new method is provided for optimizing the strip shape regulation and control efficacy coefficient.

Description

Technical Field

[0001] The invention belongs to the technical field of metallurgical rolling and relates to a method for predicting a plate shape regulation efficiency coefficient in cold-rolled plate shape control. Background Art

[0002] The shape control efficiency coefficient can be regarded as a quantitative description of the influence of the shape adjustment mechanism on the strip shape. Only by obtaining an accurate shape control efficiency coefficient can we accurately calculate the optimal adjustment amount of the shape adjustment mechanism such as bending roll, lateral displacement, and tilting required to eliminate the shape deviation. Therefore, the accuracy of the shape control efficiency coefficient directly affects the accuracy of the shape closed-loop control. How to obtain a high-precision shape control efficiency coefficient online is a fundamental issue in shape control.

[0003] Many researchers at home and abroad have used influence function method, finite element analysis and other methods to analyze the adjustment characteristics of plate shape adjustment mechanism and its influence on plate shape, providing a theoretical basis for the calculation of plate shape control efficiency coefficient. Since the factors that affect the calculation accuracy of plate shape control efficiency coefficient include strip width, strip thickness, rolling force, pressure reduction, etc., the plate shape control efficiency coefficient obtained by traditional methods such as single constant parameter mechanism model or offline simulation model has poor accuracy, lacks quantitative analysis of the transition process between two steady states, and cannot meet the plate shape control accuracy requirements under all working conditions. In order to solve this problem, some scholars have proposed a model-free adaptive optimization based entirely on rolling process data to improve the accuracy of plate shape control efficiency coefficient. However, due to the complex on-site environment and the existence of noise in the collected data, the adaptation without initial model often deviates from the correct optimization direction and fails to achieve the purpose of improving accuracy. Summary of the invention

[0004] In order to solve the above technical problems, the purpose of the present invention is to provide a method for predicting the plate shape control efficiency coefficient in cold-rolled plate shape control, and adopt a hybrid modeling method that combines simulation modeling and data-driven modeling into a series structure to establish an optimal weight factor prediction model for the actual working condition point. First, a three-dimensional finite element simulation model of the strip cold rolling process is established to obtain the plate shape control efficiency coefficient matrix of the preset working condition point, and the production data is used to improve the optimal weight factor between the preset working condition point and the actual working condition point. Then, a data-driven method is used to establish the optimal weight factor prediction model for the actual working condition point, and then the plate shape control efficiency coefficient under any rolling condition is accurately obtained.

[0005] The present invention provides a method for predicting a plate shape control efficiency coefficient in cold rolling plate shape control, comprising: Step 1: Set rolling force, strip width and strip thickness as parameters of the working point, and select multiple preset working points in the simulation experiment; Step 2: Obtain the plate shape control efficiency coefficient matrix of each preset operating point through offline simulation, and export the measured rolling data from the rolling mill; Step 3: Use the mean square error to detect and eliminate abnormal values ​​of the plate shape change in the measured rolling data; Step 4: Use Lagrange interpolation method to fill the missing data of the plate shape change after removing the outliers; Step 5: Obtain the actual operating points from the measured rolling data in sequence, and calculate the optimal weight factor array between the actual operating point and the preset operating point by using the DBO algorithm; Step 6: The parameters of the actual operating point and the optimal weight factor array constitute a data set and divide it into a training set and a test set; Step 7: Establish a BP neural network prediction model, train the BP neural network prediction model through the training set, and obtain the optimal number of hidden layer nodes; Step 8: Use the DBO algorithm to optimize the initial weights and thresholds of the BP neural network prediction model; Step 9: Use the trained BP neural network prediction model to predict the plate shape control efficiency coefficient.

[0006] The method for predicting the efficiency coefficient of plate shape regulation in cold rolling plate shape control of the present invention has the following beneficial effects: (1) First, the plate shape control efficiency coefficient matrix under different working conditions is calculated through offline simulation, and then the value range of the optimal weight factor is limited by the Euclidean distance method to prevent deviation from the correct optimization direction during the online adaptive process.

[0007] (2) The weight factor between the actual operating point and the preset operating point is determined by using the three-dimensional characteristic operating condition space representation technology, which significantly improves the calculation accuracy compared with the traditional two-dimensional operating condition division method.

[0008] (3) A data-driven approach is used to establish a prediction model for the optimal weight factor of the actual operating point, and the optimal weight factor is found using the actual operating point data, thereby accurately obtaining the plate shape control efficiency coefficient matrix under any rolling condition. The optimal weight factor prediction method for the actual operating point has not only high accuracy, but also good stability and fast convergence speed. BRIEF DESCRIPTION OF THE DRAWINGS

[0009] Figure 1 It is a flow chart of a method for predicting the efficiency coefficient of plate shape regulation in cold-rolled plate shape control of the present invention; Figure 2 a is the measured plate shape change at a certain actual working point before the mean square error detection of the present invention; Figure 2 b is the plate shape change at a certain actual working point after Lagrangian interpolation filling of the present invention; Figure 3 It is a schematic diagram of the positions of the actual operating point and the 8 related preset operating points in the present invention; Figure 4 This is a flow chart of the present invention using the DBO algorithm to optimize the BP neural network. DETAILED DESCRIPTION

[0010] like Figure 1 As shown, a method for predicting the efficiency coefficient of plate shape control in cold rolling plate shape control of the present invention comprises the following steps: Step 1: Set rolling force, strip width and strip thickness as the parameters of the working point, and select multiple preset working points in the simulation experiment, specifically: Step 1.1: Determine the upper and lower limits of the rolling force, strip width, and strip thickness parameter values ​​during actual rolling, and use them as the upper and lower limits of the parameter values ​​of the preset operating point to ensure that the actual operating point parameter values ​​during actual rolling fall within the preset operating point parameter value range.

[0011] Step 1.2: Divide the parameter value intervals of rolling force, strip width and strip thickness at the preset working condition into V Each parameter is divided into V +1 parameter value, and arrange and combine these three parameter values ​​to get ( V +1) 3 nodes, which serve as preset operating points in the offline simulation experiment.

[0012] Step 2: Obtain the plate shape control efficiency coefficient matrix of each preset working point through offline simulation, and export the measured rolling data from the rolling mill, specifically: Step 2.1: Use the parameter values ​​of each preset operating point to conduct offline simulation experiments, and give different plate shape adjustment mechanisms the adjustment amount to form a plate shape adjustment amount matrix , Δ u j For the j The adjustment amount of the plate shape adjustment mechanism is used to obtain the plate shape change matrix in the plate width direction. , Δ y i The width direction of the strip steel plate i The change of plate shape at each measuring point.

[0013] In specific implementation, the plate shape adjustment mechanism includes working roll bending, intermediate roll bending, roll lateral shifting and roll tilting, etc.

[0014] Step 2.2: Calculate the plate shape control efficiency coefficient matrix for each preset working point through the plate shape change matrix in the plate width direction and the plate shape adjustment matrix. The calculation formula is: ; In the formula, Eff is the plate shape control efficiency coefficient matrix of the preset working point, eff ij For the j The plate shape adjustment mechanism is located in the plate width direction. i The plate shape control efficiency coefficient of each measuring point.

[0015] Step 2.3: After obtaining the plate shape control efficiency coefficient matrix of several preset operating points and the parameter values ​​of the preset operating points, make a table and then save it in the form of a text file.

[0016] Step 2.4: Obtain the actual rolling data from the rolling mill, including the rolling force, strip width, strip thickness, adjustment amount of the plate shape adjustment mechanism, and the plate shape change at each measuring point in the strip width direction. Organize these data and make a table, and save them in the measured rolling data Excel file.

[0017] Step 3: Use the mean square error to detect and eliminate abnormal values ​​of the plate shape change in the measured rolling data, specifically: The mean square error is used to detect abnormal values ​​in each group of plate shape changes under actual working conditions. μ- 3 σ , μ+ 3 σ ] are considered valid data, otherwise they are considered abnormal data. μ and σ The calculation formula is as follows: ; Where: Δ y oi The actual working point O The width direction of the board i The measured plate shape change at each measuring point is M The number of measuring points in the width direction of the strip.

[0018] Figure 2 a It is the measured plate shape change before mean square error detection under the parameters of rolling force 255KN, strip width 1275mm and strip thickness 0.5mm.

[0019] Step 4: Use the Lagrange interpolation method to fill the missing data of the plate shape change after removing the outliers, specifically: Step 4.1: Select a set of plate shape variation from the plate shape variation data set after removing outliers, and select the plate shape variation that has not been removed from this set of plate shape variation. L Measuring points: ,in C i Indicates the width direction of the board i The serial number of the measurement point, Δy oi Indicates the width direction of the board i The change of plate shape at each measuring point.

[0020] Step 4.2: According to the given L measurement points, construct the Lagrange basic polynomial corresponding to each measurement point: ; in, l i ( C ) is the i The Lagrange basic polynomial of the measurement points is: C r Indicates the serial number of the measurement point to be eliminated.

[0021] Step 4.3: Use the Lagrange basic polynomial corresponding to each measurement point obtained in step 4.2 l i ( C ) Calculate the known L The Lagrange interpolation polynomial corresponding to the measurement points L ( C ), which is the first r The change of plate shape at each measuring point Δ y or , complete the filling of missing characteristic parameters, Lagrange interpolation polynomial L ( C )for: .

[0022] Figure 2 b It is the change of plate shape after Lagrangian interpolation filling under the parameters of rolling force 255KN, strip width 1275mm and strip thickness 0.5mm.

[0023] Step 5: Obtain the actual operating points from the measured rolling data in sequence, and calculate the optimal weight factor array between the actual operating point and the preset operating point through the DBO algorithm, specifically: Step 5.1: Figure 3 As shown, let the strip width be the X axis, the strip thickness be the Y axis, and the rolling force be the Z axis. The three parameters form a three-dimensional space structure, and the three coordinate axes have V +1 parameter value, which divides the three-dimensional space structure into V 3The actual working point O is obtained from the measured rolling data in sequence, the small cube where the actual working point O is located is determined, and the 8 vertices of the small cube, i.e., the 8 preset working points related to the actual working point O, are determined.

[0024] Step 5.2: Calculate the weight factor γ between the actual operating point O and the eight related preset operating points by using the Euclidean distance method. 1、 γ 2、 γ 3、 γ 4、 γ 5、 γ 6、 γ 7、 γ 8, , the calculation formula of the weight factor calculated by the Euclidean distance method is as follows: ; ; ; ; ; ; ; ; Where: x 0. y 0. z 0 are the strip width value, strip thickness value and rolling force value at the actual working point O respectively; x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8 are the strip width values ​​of 8 preset working points, y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8 are the strip thickness values ​​of 8 preset working points, z 1. z 2. z 3. z 4. z 5. z 6. z 7. z8 are rolling force values ​​of 8 preset working points respectively.

[0025] Step 5.3: The weight factor calculated by the Euclidean distance method is not the optimal weight factor, but it represents the similarity between the preset operating point and the actual operating point to a certain extent. The optimal weight factor value range is set based on the weight factor calculated by the Euclidean distance method to prevent the model from deviating from the correct direction when searching for the optimal weight factor. The optimal weight factor value range is as follows: ; Where: The actual working point and k The optimal weight factor among the preset operating points, γ k is the distance between the actual operating point calculated by the Euclidean distance method and the k The weight factor between the preset operating points, k =1,2,3,4,5,6,7,8.

[0026] Step 5.4: Use the DBO algorithm to obtain the optimal weight factor array. The mathematical model of the DBO algorithm is as follows: create R =90 arrays, each consisting of H The maximum number of iterations of the DBO algorithm in this implementation is T =100, this step is to find the optimal weight factor array, so in this step H =8, each number represents a weight factor, and the values ​​of the 8 weight factors are randomly generated within the range of the weight factors in step 5.3. Use the evaluation function to evaluate the quality of the array. The lower the evaluation function value, the better the array. The evaluation function calculation method is as follows: ; Where: J q is the evaluation function, g i is the weight coefficient of each measuring point in the board width direction. The weight coefficient of the measuring points within 60% of the board width center is 0.8, and the weight coefficient of the measuring points on the remaining two sides is 1; Δ u oj is the actual working point O j The actual measured adjustment amount of the plate shape adjustment mechanism, N is the number of plate shape adjustment mechanisms; Eff oq The actual operating point O is created by the DBO algorithm. q The weight factor array is the board shape control efficiency coefficient matrix calculated by the weight factor, and the calculation formula is: ; Where: Eff s is the first value related to the actual operating point O calculated in step 2 s The plate shape control efficiency coefficient matrix of the preset working point, X qs Created for q The first of the weight factor arrays s Number, s =1,2,3,4,5,6,7,8.

[0027] Step 5.5: Create R The weight factor arrays are divided into four parts according to 6:6:7:11. The evaluation functions of all weight factor arrays are calculated. The four weight factor arrays are sorted from low to high according to the evaluation function. At the same time, all weight factor arrays are sorted from low to high according to the evaluation function. The array with the highest evaluation function among all weight factor arrays is taken as X worst , the array with the lowest evaluation function is X best The array with the smallest evaluation function in the second part of the weight factor array is X gbest , the array with the smallest evaluation function in the third part of the weight factor array is X lbest , each part of the weight factor array is updated through different iterative methods to find the weight factor array with the highest adaptability.

[0028] (1) There are two iterative update formulas for the weight factor array in the first part, which are: ; ; Where: t Indicates the current iteration number, Indicates t After the iteration n array, 0< c ≤0.2 is a constant value. B is a constant value between 0 and 1. α is a natural coefficient assigned a value of -1 or 1; Indicates t After the iteration, the array with the worst fitness is randomly generated in [0,π]. When θ=0, θ=π / 2, and θ=π, the array will not be updated in this round.

[0029] Each array in the first part will take a random number between 0 and 1 before iterative update. If the random number is less than 0.9, it will be iteratively updated according to the first formula, otherwise it will be iteratively updated according to the second formula.

[0030] (2) The iterative update formula of the second part weight factor array is: ; Where: B 1, B 2 is a 1×8 independent random variable, Lb * , Ub * is the upper and lower bounds of the second part of the weight factor array, and the calculation formula is: ; Where: For the t After iterations, the array with the smallest evaluation function in the second part of the weight factor array is Lb is the lower limit array of the optimal weight factor value range in step 5.3, Ub is the upper limit array of the optimal weight factor value range in step 5.3, R =(1- t ) / T .

[0031] (3) The iterative update formula of the weight factor array in the third part is: ; Where: C 1 is a random number that follows a normal distribution. C 2 is a 1×8 random vector between (0,1); Lb l , Ub l is the upper and lower bounds of the third part of the array, and the calculation formula is: ; Where: For the t The array with the smallest evaluation function in the third part of the weight factor array after iterations.

[0032] (4) The fourth part of the array iteration update formula is: ; Where: G is a 1×D random vector between (0,1), S is a constant.

[0033] Each time the array is iterated and updated, the evaluation function is calculated once, and the new ;Will X worst Arrays and The evaluation function of the array is compared, and the one with the higher evaluation function becomes the newX worst ;Will X best Arrays and The evaluation function of the array is compared, and the one with the lower evaluation function becomes the new X best ;Will X gbest Arrays and The evaluation function of the array is compared, and the one with the lower evaluation function becomes the new X gbest ;Will X lbest Arrays and The evaluation function of the array is compared, and the one with the lower evaluation function becomes the new X lbest .

[0034] Step 5.6: Determine the number of iterations t Is it greater than the maximum number of iterations? T If it is not satisfied, continue to step 5.5. If it is satisfied, stop the iteration and get T The array with the smallest evaluation function after generation X best , and use this set of arrays as the optimal weight factor array of the actual operating point O.

[0035] Step 5.7: Find the optimal weight factor array for all actual operating points in the measured production data, and make a table of the rolling force, strip width, strip thickness and optimal weight factor array of these actual operating points as the data set for BP neural network training. Shuffle the data in the data set, use the first 80% as the training set, and the last 20% as the test set.

[0036] Step 6: The parameters of the actual operating points and the optimal weight factor array constitute a data set and divide it into a training set and a test set.

[0037] Step 7: Establish a BP neural network prediction model. Train the BP neural network prediction model with the training set to obtain the optimal number of hidden layer nodes, specifically: Step 7.1: Determine the BP neural network structure, where the input nodes are 3, which are rolling force, strip width and strip thickness respectively; the output nodes are 8, which are the optimal weight factors between the actual working point and the 8 related preset working points. The range of the number of hidden layer nodes of the BP neural network is determined according to the empirical formula: ; Where: , h is the number of hidden layer nodes of BP neural network, a is the number of input nodes, b is the number of output nodes.

[0038] The calculated range of the number of nodes in the hidden layer of the BP neural network is h =[5,13].

[0039] Step 7.2: Input the training set, set the initial weights and thresholds of the BP neural network, including the connection weights between the input layer and the hidden layer. w ef , the threshold between the input layer and the hidden layer is b f ; The connection weight between the hidden layer and the output layer is w fg , the threshold between the hidden layer and the output layer is b g .

[0040] Neural Network Input P e and output Q The calculation formula is: ; Where: P e For the neural network e input variables, F is the Sigmoid activation function, f Indicates f hidden nodes.

[0041] Step 7.3: Calculate the predicted output value of the BP neural network under different numbers of hidden layer nodes Q dg and target value The mean square error is calculated as: ; Where, MSE is the total mean square error of the BP neural network output. E is the number of samples, Q dg For the g The output of d The predicted output value of sample data, For the g The output of d The target value of the sample data.

[0042] Step 7.4: Take the number of hidden layer nodes when the mean square error value is the smallest as the optimal number of hidden layer nodes, and take the number of hidden layer nodes as 10.

[0043] The experimental results are shown in Table 1. When the number of hidden layer nodes is 10, the minimum mean square error is 0.0038310, so the BP neural network is set to a 3-10-8 structure.

[0044] Table 1 Mean square error of training set under different hidden layer nodes

[0045] Step 8: Use the DBO algorithm to optimize the initial weights and thresholds of the BP neural network prediction model to improve the prediction accuracy and convergence speed of the BP neural network. Figure 4 As shown, specifically: Step 8.1: Create R = 90 arrays, each array consists of the initial weights and thresholds of the BP neural network w ef , b f , w fg , b g The elements of v =121 numbers, let R In the array u The array is a u ,Right now: ; in ; ; ; ; Step 8.2: Use fitness to represent the prediction ability of the BP neural network under different initial weights and thresholds, and calculate the predicted output values ​​of the BP neural network under different initial weights and thresholds. Q dg and target value The root mean square error is used as the fitness value of all arrays. The better the prediction ability of the neural network, the smaller the fitness value.

[0046] Step 8.3: Use the DBO algorithm to iteratively update the array, find the array with the minimum fitness, and determine the optimal initial weights and thresholds of the BP neural network; follow step 5.5 to update the array. R Arrays T =100 rounds of iterative updates, select the array with the smallest fitness after 100 iterative updates a best , as the optimal initial weight and threshold.

[0047] Step 8.4: Use the test set to verify the prediction accuracy of the trained BP neural network prediction model.

[0048] Step 9: Use the trained BP neural network prediction model to predict the plate shape control efficiency coefficient.

[0049] The above description is only a preferred embodiment of the present invention and is not intended to limit the concept of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A method for predicting the efficiency coefficient of plate shape control in cold rolling plate shape control, characterized in that: include: Step 1: Set rolling force, strip width and strip thickness as parameters of the working point, and select multiple preset working points in the simulation experiment; Step 2: Obtain the plate shape control efficiency coefficient matrix of each preset operating point through offline simulation, and export the measured rolling data from the rolling mill; Step 3: Use the mean square error to detect and eliminate abnormal values ​​of the plate shape change in the measured rolling data; Step 4: Use Lagrange interpolation method to fill the missing data of the plate shape change after removing the outliers; Step 5: Obtain the actual operating points from the measured rolling data in sequence, and calculate the optimal weight factor array between the actual operating point and the preset operating point by using the DBO algorithm; Step 6: The parameters of the actual operating point and the optimal weight factor array constitute a data set and divide it into a training set and a test set; Step 7: Establish a BP neural network prediction model, train the BP neural network prediction model through the training set, and obtain the optimal number of hidden layer nodes; Step 8: Use the DBO algorithm to optimize the initial weights and thresholds of the BP neural network prediction model; Step 9: Use the trained BP neural network prediction model to predict the plate shape control efficiency coefficient.

2. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 1, characterized in that: The step 1 is specifically as follows: Step 1.1: Determine the upper and lower limits of the rolling force, strip width, and strip thickness parameter values ​​during actual rolling, and use them as the upper and lower limits of the parameter values ​​of the preset working point, so as to ensure that the actual working point parameter values ​​during actual rolling fall within the preset working point parameter value range; Step 1.2: Divide the parameter value intervals of rolling force, strip width and strip thickness at the preset working condition into V Each parameter is divided into V +1 parameter value, and arrange and combine these three parameter values ​​to get ( V +1) 3 nodes, which serve as preset operating points in the offline simulation experiment.

3. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 1, characterized in that: The step 2 is specifically as follows: Step 2.1: Use the parameter values ​​of each preset operating point to conduct offline simulation experiments, and give different plate shape adjustment mechanisms the adjustment amount to form a plate shape adjustment amount matrix , Δ u j For the j The adjustment amount of the plate shape adjustment mechanism is used to obtain the plate shape change matrix in the plate width direction. , Δ y i The width direction of the strip steel plate i The change of plate shape at each measuring point; Step 2.2: Calculate the plate shape control efficiency coefficient matrix for each preset working point through the plate shape change matrix in the plate width direction and the plate shape adjustment matrix. The calculation formula is: ； In the formula, Eff is the plate shape control efficiency coefficient matrix of the preset working point, eff ij For the j The plate shape adjustment mechanism is located in the plate width direction. i Plate shape control efficiency coefficient of each measuring point; Step 2.3: After obtaining the plate shape control efficiency coefficient matrix of several preset working points and the parameter values ​​of the preset working points, make a table and then save it in the form of a text file; Step 2.4: Obtain the actual rolling data from the rolling mill, including the rolling force, strip width, strip thickness, adjustment amount of the plate shape adjustment mechanism, and the plate shape change at each measuring point in the strip width direction. Organize these data and make a table, and save them in the measured rolling data Excel file.

4. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 3, characterized in that: The step 3 is specifically as follows: The mean square error is used to detect abnormal values ​​in each group of plate shape changes under actual working conditions. μ- 3 σ , μ+ 3 σ ] are considered valid data, otherwise they are considered abnormal data. μ and σ The calculation formula is as follows: ； Where: Δ y oi The actual working point O The width direction of the board i The measured plate shape change at each measuring point is M The number of measuring points in the width direction of the strip.

5. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 4, characterized in that: The step 4 is specifically as follows: Step 4.1: Select a set of plate shape variation from the plate shape variation data set after removing outliers, and select the plate shape variation that has not been removed from this set of plate shape variation. L Measuring points: ,in C i Indicates the width direction of the board i The serial number of the measurement point, Δ y oi Indicates the width direction of the board i The change of plate shape at each measuring point; Step 4.2: According to the given L measurement points, construct the Lagrange basic polynomial corresponding to each measurement point: ； in, l i ( C ) is the i The Lagrange basic polynomial of the measurement points is: C r Indicates the serial number of the measurement point to be eliminated; Step 4.3: Use the Lagrange basic polynomial corresponding to each measurement point obtained in step 4.2 l i ( C ) Calculate the known L The Lagrange interpolation polynomial corresponding to the measurement points L ( C ), which is the first r The change of plate shape at each measuring point Δ y or , complete the filling of missing characteristic parameters, Lagrange interpolation polynomial L ( C )for: 。 6. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 5, characterized in that: The step 5 is specifically as follows: Step 5.1: Let the strip width be the X axis, the strip thickness be the Y axis, and the rolling force be the Z axis. The three parameters form a three-dimensional space structure. The three coordinate axes have V +1 parameter value, which divides the three-dimensional space structure into V 3 Small cubes of the same size, each cube has 8 vertices, i.e. 8 preset working points; The actual working point O is obtained from the measured rolling data in sequence, a small cube where the actual working point O is located is determined, and 8 vertices of the small cube are determined, that is, 8 preset working points related to the actual working point O; Step 5.2: Calculate the weight factor γ between the actual operating point O and the eight related preset operating points by using the Euclidean distance method. 1、 γ 2、 γ 3、 γ 4、 γ 5、 γ 6、 γ 7、 γ 8, , the calculation formula of the weight factor calculated by the Euclidean distance method is as follows: ； ； ； ； ； ； ； ； Where: x 0. y 0. z 0 are the strip width value, strip thickness value and rolling force value at the actual working point O respectively; x 1. x 2. x 3. x 4. x 5. x 6. x 7. x 8 are the strip width values ​​of 8 preset working points, y 1. y 2. y 3. y 4. y 5. y 6. y 7. y 8 are the strip thickness values ​​of 8 preset working points, z 1. z 2. z 3. z 4. z 5. z 6. z 7. z 8 are rolling force values ​​of 8 preset working conditions respectively; Step 5.3: Set the optimal weight factor value range based on the weight factor calculated by the Euclidean distance method to prevent the model from deviating from the correct direction when searching for the optimal weight factor. The optimal weight factor value range is as follows: ； Where: The actual working point and k The optimal weight factor among the preset operating points, γ k is the distance between the actual operating point calculated by the Euclidean distance method and the k The weight factor between the preset operating points, k =1,2,3,4,5,6,7,8; Step 5.4: Use the DBO algorithm to obtain the optimal weight factor array. The mathematical model of the DBO algorithm is as follows: create R = 90 weight factor arrays, each weight factor array consists of 8 weight factors, and the maximum number of iterations of the DBO algorithm is T =100, 8 weight factors are randomly generated within the range of the optimal weight factor values ​​in step 5.3; an evaluation function is used to evaluate the quality of each weight factor array. The lower the evaluation function value, the better the weight factor array. The evaluation function calculation method is as follows: ； Where: J q is the evaluation function, g i is the weight coefficient of each measuring point in the board width direction. The weight coefficient of the measuring points within 60% of the board width center is 0.8, and the weight coefficient of the measuring points on the remaining two sides is 1; Δ u oj is the actual working point O j The actual measured adjustment amount of the plate shape adjustment mechanism, N is the number of plate shape adjustment mechanisms; Eff oq The actual operating point O is created by the DBO algorithm. q The weight factor array is the board shape control efficiency coefficient matrix calculated by the weight factor, and the calculation formula is: ； Where: Eff s is the first value related to the actual operating point O calculated in step 2 s The plate shape control efficiency coefficient matrix of the preset working point, X qs Created for q The first of the weight factor arrays s Number, s =1,2,3,4,5,6,7,8; Step 5.5: Divide the created R weight factor arrays into 4 parts, calculate the evaluation function of all weight factor arrays, sort the four weight factor arrays from low to high according to the evaluation function, and sort all weight factor arrays from low to high according to the evaluation function. The array with the highest evaluation function among all weight factor arrays is taken as X worst , the array with the lowest evaluation function is X best ; The array with the smallest evaluation function in the second part of the weight factor array is X gbest , the array with the smallest evaluation function in the third part of the weight factor array is X lbest , each part of the weight factor array is updated through different iterations to find the weight factor array with the highest adaptability; (1) There are two iterative update formulas for the weight factor array in the first part, which are: ； ； Where: t Indicates the current iteration number, Indicates t After the iteration n arrays, 0< c ≤0.2 is a constant value. B is a constant value between 0 and 1. α is a natural coefficient assigned a value of -1 or 1; Indicates t The array with the worst fitness after iterations, θ is randomly generated in [0,π]. When θ=0, θ=π / 2, θ=π, the array will not be updated in this round; Each array in the first part will take a random number between 0 and 1 before iterative update. If the random number is less than 0.9, it will be iteratively updated according to the first formula, otherwise it will be iteratively updated according to the second formula; (2) The iterative update formula of the second part weight factor array is: ； Where: B 1, B 2 is a 1×8 independent random variable, Lb * , Ub * is the upper and lower bounds of the second part of the weight factor array, and the calculation formula is: ； Where: For the t After iterations, the array with the smallest evaluation function in the second part of the weight factor array is Lb is the lower limit array of the optimal weight factor value range in step 5.3, Ub is the upper limit array of the optimal weight factor value range in step 5.3, R =(1- t ) / T ; (3) The iterative update formula of the weight factor array in the third part is: ； Where: C 1 is a random number that follows a normal distribution. C 2 is a 1×8 random vector between (0,1); Lb l , Ub l is the upper and lower bounds of the third part of the array, and the calculation formula is: ； Where: For the t The array with the smallest evaluation function in the third part of the weight factor array after iterations; (4) The iterative update formula of the fourth part weight factor array is: ； Where: G is a 1×8 random vector between (0,1), S is a constant; Each time the array is iterated and updated, the evaluation function is calculated once, and the new ;Will X worst Arrays and The evaluation function of the array is compared, and the one with the higher evaluation function becomes the new X worst ;Will X best Arrays and The evaluation function of the array is compared, and the one with the lower evaluation function becomes the new X best ;Will X gbest Arrays and The evaluation function of the array is compared, and the one with the lower evaluation function becomes the new X gbest ;Will X lbest Arrays and The evaluation function of the array is compared, and the one with the lower evaluation function becomes the new X lbest ; Step 5.6: Determine the number of iterations t Is it greater than the maximum number of iterations? T If it is not satisfied, continue to step 5.

5. If it is satisfied, stop the iteration and get T The array with the smallest evaluation function after generation X best , take this set of arrays as the optimal weight factor array of the actual operating point O; Step 5.7: Find the optimal weight factor array for all actual operating points in the measured production data, and make a table of the rolling force, strip width, strip thickness and optimal weight factor array of these actual operating points as the data set for BP neural network training. Shuffle the data in the data set, use the first 80% as the training set, and the last 20% as the test set.

7. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 6, characterized in that: The step 7 is specifically as follows: Step 7.1: Determine the BP neural network structure, where the input nodes are 3, which are rolling force, strip width and strip thickness respectively; the output nodes are 8, which are the optimal weight factors between the actual working point and the 8 related preset working points. The range of the number of hidden layer nodes of the BP neural network is determined according to the empirical formula: ； Where: , h is the number of hidden layer nodes of BP neural network, a is the number of input nodes, b is the number of output nodes; The calculated range of the number of nodes in the hidden layer of the BP neural network is h =[5,13]; Step 7.2: Input the training set, set the initial weights and thresholds of the BP neural network, including the connection weights between the input layer and the hidden layer. w ef , the threshold between the input layer and the hidden layer is b f ; The connection weight between the hidden layer and the output layer is w fg , the threshold between the hidden layer and the output layer is b g ; Neural Network Input P e and output Q The calculation formula is: ； Where: P e For the neural network e input variables, F is the Sigmoid activation function, f Indicates f hidden nodes; Step 7.3: Calculate the predicted output value of the BP neural network under different numbers of hidden layer nodes Q dg and target value The mean square error is calculated as: ； Where, MSE is the total mean square error of the BP neural network output. E is the number of samples, Q dg For the g The output of d The predicted output value of sample data, For the g The output of d The target value of sample data; Step 7.4: Take the number of hidden layer nodes when the mean square error value is the smallest as the optimal number of hidden layer nodes, and take the number of hidden layer nodes as 10.

8. The method for predicting the efficiency coefficient of shape control in cold rolling shape control according to claim 7, characterized in that: The step 8 is as follows: Step 8.1: Create R = 90 arrays, each array consists of the initial weights and thresholds of the BP neural network w ef , b f , w fg , b g The elements of v =121 numbers, let R In the array u The array is a u ,Right now: ； in ； ； ； ； Step 8.2: Use fitness to represent the prediction ability of the BP neural network under different initial weights and thresholds, and calculate the predicted output values ​​of the BP neural network under different initial weights and thresholds. Q dg and target value The root mean square error is used as the fitness value of all arrays. The better the prediction ability of the neural network, the smaller the fitness value. Step 8.3: Use the DBO algorithm to iteratively update the array, find the array with the minimum fitness, and determine the optimal initial weights and thresholds of the BP neural network; follow step 5.5 to update the array. R Arrays T =100 rounds of iterative updates, select the array with the smallest fitness after 100 iterative updates a best , as the optimal initial weight and threshold; Step 8.4: Use the test set to verify the prediction accuracy of the trained BP neural network prediction model.

Images