Strip steel cold rolling non-uniform speed field rolling force prediction method and system

Abstract

The application discloses a strip steel cold rolling non-uniform speed field rolling force prediction method and system, comprising: obtaining the process parameters and strip parameters of cold continuous rolling; according to the process parameters and strip parameters, establishing the roll contact arc curve and its first derivative equation; according to the roll contact arc curve, establishing the non-uniform speed field and strain speed field of the rolling deformation zone; according to the non-uniform speed field and strain speed field, calculating multiple power items of the plastic deformation zone of the cold rolling slab; according to the multiple power items, constructing the total power functional and performing the minimum processing on the total power functional to obtain the neutral angle; and according to the minimum total power functional and the neutral angle, calculating the rolling torque and rolling force. The application can accurately capture the speed gradient difference in the thickness and width directions, effectively reduces the load calculation error under the complex working condition, simultaneously provides a reliable theoretical basis for the shape control and process optimization, and takes into account the high-precision calculation and good engineering applicability.

Description

Technical Field

[0001] This invention belongs to the field of cold-rolled strip rolling, and particularly relates to a method and system for predicting rolling force in non-uniform speed fields during cold rolling of strip steel. Background Technology

[0002] Rolling, as one of the core processes in the field of metal plastic forming, plays an irreplaceable and crucial role in the continuous and large-scale production of sheet and strip products. By applying controllable pressure to metal sheets and strips through rolls, continuous plastic deformation is achieved, thereby precisely obtaining the target geometry, dimensional accuracy, and excellent internal microstructure. Among these parameters, rolling force, as a core process parameter characterizing the essence of the rolling process, directly determines the quality and process stability of the sheet and strip products through accurate prediction. Especially against the backdrop of the continuously growing demand for high-quality sheet metal in high-end manufacturing, developing high-precision rolling force models has become a key scientific issue driving the advancement of rolling technology. In previous idealized theoretical models, researchers often simplified the steady-state rolling process to uniform deformation, assuming that the velocity field of the metal in the deformation zone is uniform along the thickness and width directions, with a continuous velocity gradient change only along the rolling direction.

[0003] However, numerous experimental observations and numerical simulations clearly demonstrate that the material flow behavior during actual steady-state rolling is far from the idealized uniform state. A significant non-uniform velocity field exists within the deformation zone, and this spatial non-uniformity is an inherent physical characteristic of the rolling process. Traditional rolling force models based on the assumption of uniform deformation cannot accurately describe the true velocity field characteristics of metal flow within the deformation zone, and struggle to capture velocity gradient differences along the thickness and width directions. This results in insufficient accuracy in predicting rolling forces under complex conditions, large load calculation errors, and consequently, affects the accuracy of shape control and process optimization. Therefore, how to construct a mathematical model that can accurately describe the spatial distribution characteristics of the non-uniform velocity field within the deformation zone, thereby improving the accuracy of rolling force models in representing complex deformation processes, has become a pressing technical problem to be solved in this field. Summary of the Invention

[0004] To address the aforementioned technical problems, the present invention aims to provide a model prediction method capable of accurately describing the spatial distribution characteristics of the non-uniform velocity field within the rolling deformation zone. This invention aims to overcome the limitations of traditional uniform deformation theory by deeply revealing the true velocity field characteristics of metal flow, establishing a more accurate rolling force model, and thus more accurately predicting the magnitude of the rolling force generated during the rolling process, thereby improving the characterization accuracy of complex deformation processes.

[0005] To achieve the above objectives, this invention provides a method and system for predicting rolling force in a non-uniform velocity field during cold rolling of strip steel. Specifically, the method for predicting rolling force in a non-uniform velocity field during cold rolling of strip steel includes: Obtain the process parameters and strip parameters for cold continuous rolling; Based on the process parameters and strip parameters, establish the roll contact arc curve and its first derivative equation; Based on the roll contact arc curve, establish the non-uniform velocity field and strain velocity field of the rolling deformation zone; Based on the non-uniform velocity field and strain velocity field, calculate various power terms in the plastic deformation zone of the cold-rolled slab. Based on the various power terms, a total power functional is constructed and the total power functional is minimized to obtain the neutral angle; The rolling torque and rolling force are calculated based on the minimum total power functional and the neutral angle.

[0006] Preferably, the process of obtaining the process parameters and strip parameters for cold continuous rolling includes: Obtain the inlet thickness, outlet thickness, inlet width, initial roll radius, roll linear speed, front tension, and back tension.

[0007] Preferably, the process of establishing the roll contact arc curve and its first derivative equation includes: Based on the boundary conditions and geometric characteristics of the strip deformation zone, the roll contact arc curve and its first derivative equation are established.

[0008] Preferably, the process of establishing the non-uniform velocity field and strain velocity field in the rolling deformation zone includes: Based on the velocity boundary conditions of the strip rolling deformation zone, the principle of constant volume, and the roll contact arc curve, the non-uniform velocity field and strain velocity field of the rolling deformation zone are established.

[0009] Preferably, the process of calculating multiple power terms in the plastic deformation zone of cold-rolled slab includes: Calculate the internal deformation power based on the non-uniform velocity field and strain velocity field; The frictional power is calculated based on the non-uniform velocity field and the frictional conditions of the rolled piece surface. Calculate the shear power based on the shear force at the velocity discontinuity at the entry point of the workpiece; Calculate the tension power based on the influence of front tension and back tension on the rolling process.

[0010] Preferably, after obtaining the rolling force, the method further includes iterative calculation of the roll radius; The iterative calculation process includes: The actual radius of the roll is calculated based on the rolling force, Poisson's ratio of the roll, and the elastic modulus of the roll. Substitute the updated actual roll radius into the steps of establishing the roll contact arc curve, re-establish the roll contact arc curve and iteratively calculate the rolling force until the convergence condition is met.

[0011] This invention also provides a rolling force prediction system for non-uniform velocity fields in cold rolling of strip steel, comprising: The parameter input module is used to obtain the process parameters and strip parameters of cold continuous rolling; The geometric modeling module, connected to the parameter input module, is used to establish the roll contact arc curve and its first derivative equation based on the process parameters and strip parameters. A velocity field construction module, connected to the geometric modeling module, is used to establish a non-uniform velocity field and strain velocity field in the rolling deformation zone based on the roll contact arc curve. The power calculation module, connected to the velocity field construction module, is used to calculate various power terms in the plastic deformation zone of the cold-rolled slab based on the non-uniform velocity field and strain velocity field. An optimization solution module, connected to the power calculation module, is used to construct a total power functional based on the various power terms, and to minimize the total power functional to obtain the neutral angle. The rolling force calculation module, connected to the optimization solution module, is used to calculate the rolling torque and rolling force based on the minimum total power functional and the neutral angle.

[0012] Preferably, the process parameters and strip parameters acquired by the parameter input module include inlet thickness, outlet thickness, inlet width, initial radius of the roll, linear speed of the roll, front tension, and back tension.

[0013] Preferably, the power calculation module includes: An internal deformation power calculation unit is used to calculate the internal deformation power based on the non-uniform velocity field and strain velocity field. Friction power calculation unit, used to calculate friction power based on the non-uniform velocity field and friction conditions of the workpiece surface; The shear power calculation unit is used to calculate the shear power based on the shear force at the velocity discontinuity at the entry point of the workpiece. The tension power calculation unit is used to calculate the tension power based on the influence of the front tension and the back tension on the rolling process.

[0014] Preferably, the system further includes an iterative correction module, which is connected to the parameter input module and the rolling force calculation module respectively. The iterative correction module is used to calculate the actual radius of the roll based on the rolling force, the roll Poisson's ratio and the roll elastic modulus, and feed the updated roll radius back to the geometric modeling module. The geometric modeling module re-establishes the roll contact arc curve based on the updated roll radius and iteratively calculates the rolling force until the convergence condition is met.

[0015] Compared with the prior art, the present invention has the following advantages and technical effects: This invention overcomes the limitations of the uniform deformation assumption in traditional rolling theory. By constructing a non-uniform velocity field model, it accurately describes the spatial non-uniform distribution characteristics of the velocity field within the deformation zone, significantly improving the prediction accuracy of the rolling force model. Verified by actual production data, the maximum prediction error of this invention's model is 4.78%, and the average error is 3.92%, fully meeting the accuracy requirements for industrial applications. This technical solution not only achieves a refined characterization of metal flow behavior, accurately capturing velocity gradient differences in the thickness and width directions and effectively reducing load calculation errors under complex working conditions, but also provides a reliable theoretical basis for shape control and process optimization, balancing high-precision calculation with good engineering applicability. Attached Figure Description

[0016] The accompanying drawings, which form part of this application, are used to provide a further understanding of this application. The illustrative embodiments and descriptions of this application are used to explain this application and do not constitute an undue limitation of this application. In the drawings: Figure 1 This is a schematic diagram of the cold-rolled deformation zone according to an embodiment of the present invention; Figure 2 This is a data graph showing the measured rolling force values ​​in an embodiment of the present invention. Detailed Implementation

[0017] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0018] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0019] Example 1 This embodiment provides a method for predicting rolling force in a non-uniform velocity field during cold rolling of strip steel, including: Obtain the process parameters and strip parameters for cold continuous rolling; Based on the process parameters and strip parameters, establish the roll contact arc curve and its first derivative equation; Based on the roll contact arc curve, establish the non-uniform velocity field and strain velocity field of the rolling deformation zone; Based on the non-uniform velocity field and strain velocity field, calculate various power terms in the plastic deformation zone of cold-rolled slab; Based on multiple power terms, a total power functional is constructed and minimized to obtain the neutral angle; Calculate the rolling torque and rolling force based on the minimum total power functional and the neutral angle.

[0020] Furthermore, the process of obtaining the process parameters and strip parameters for cold continuous rolling includes: Obtain the inlet thickness, outlet thickness, inlet width, initial roll radius, roll linear speed, front tension, and back tension.

[0021] Furthermore, this embodiment collects and inputs the process and strip parameters of cold continuous rolling, including the inlet thickness 2h1, the outlet thickness 2h0, the inlet width 2B, the initial roll radius R0, and the roll linear speed v. r , pretension σ f and after-tension σ b wait.

[0022] Furthermore, the process of establishing the roll contact arc curve and its first derivative equation includes: Based on the boundary conditions and geometric characteristics of the strip deformation zone, the roll contact arc curve and its first derivative equation are established.

[0023] Furthermore, in this embodiment, the thickness and width directions of the sheet metal within the rolling deformation zone exhibit strict geometric symmetry. To simplify calculations and improve efficiency, only one-quarter of the symmetrical body is considered. Based on the boundary conditions and geometric characteristics of the strip deformation zone, the roll contact arc curve and its first-order derivative equation are established: MERGEFORMAT (1) MERGEFORMAT (2) MERGEFORMAT (3) Where x is any position within the deformation region, h x It is half the thickness of the rolled piece at any position x in the deformation zone.

[0024] Furthermore, the process of establishing the non-uniform velocity field and strain velocity field in the rolling deformation zone includes: Based on the velocity boundary conditions of the strip rolling deformation zone, the principle of constant volume, and the roll contact arc curve, the non-uniform velocity field and strain velocity field of the rolling deformation zone are established.

[0025] Furthermore, this embodiment establishes a non-uniform velocity field and strain velocity field in the rolling deformation zone that satisfy the motion permitting conditions, based on the velocity boundary conditions of the strip rolling deformation zone, the principle of constant volume, and the roll contact arc curve.

[0026] MERGEFORMAT (4) MERGEFORMAT (5) MERGEFORMAT (6) MERGEFORMAT (7) MERGEFORMAT (8) Where y is the vertical coordinate within the deformation zone during rolling, and v x It is the horizontal flow velocity within the deformation zone during the rolling process, v y H(x) is the vertical flow velocity within the deformation zone during rolling, H(x) is the non-uniformity coefficient, U is the volumetric flow rate per second, and h is the flow rate in the vertical direction. n It is half the thickness at the neutral plane. θ n It is the angle between the neutral plane and the y-axis. It is the strain component in the horizontal direction. It is the strain component in the vertical direction.

[0027] Furthermore, the process of calculating various power terms in the plastic deformation zone of cold-rolled slab includes: Calculate the internal deformation power based on the non-uniform velocity field and strain velocity field; Calculate the friction power based on the non-uniform velocity field and the friction conditions of the rolled workpiece surface; Calculate the shear power based on the shear force at the velocity discontinuity at the entry point of the workpiece; Calculate the tension power based on the influence of front tension and back tension on the rolling process.

[0028] Furthermore, based on the established non-uniform velocity field and strain velocity field, this embodiment calculates the internal deformation power, friction power, shear power, and tension power of the plastic deformation zone of the cold-rolled slab; Specifically, the internal plastic deformation power W of the rolled piece p : MERGEFORMAT (9) MERGEFORMAT (10) MERGEFORMAT (11) MERGEFORMAT (12) Where D1, D2, and D3 are integration constants, and σ s The deformation resistance of the material; Surface friction power of rolled workpiece W f : MERGEFORMAT(13) MERGEFORMAT (14) Where A1 to A3 and D4 to D6 are integration constants, and k is the shear stress of the material; Shearing power of rolled product (W) s : According to slip line theory, when a velocity discontinuity exists, the power of the shear force on that surface must be considered. For the rolling process, the flow velocity of the workpiece within the entire deformation zone changes continuously in both the vertical and horizontal directions. Only in the vertical direction at the inlet does the flow velocity abruptly change from zero to a certain value. Therefore, the velocity discontinuity exists only at the inlet of the workpiece, where the shear power is: MERGEFORMAT (15) Tensile stress power W t : In cold rolling, the transverse stress distribution of the strip between the rolls is balanced by front and rear tension to suppress shape defects such as edge waviness and center waviness. Changes in tension affect all rolling parameters. Therefore, to make the rolling model more accurate, the influence of tensile stress needs to be considered. The tensile stress power during rolling can be expressed as: MERGEFORMAT (16) Furthermore, this embodiment minimizes the total power functional to obtain the neutral angle; Specifically, this involves summing the internal deformation power, friction power, shear power, and tension power to obtain the total power functional W. z .

[0029] MERGEFORMAT (17) Minimizing the total power functional yields the neutral angle. The total power is a function of the neutral angle. According to the first variational principle for rigid-plastic materials, the true solution minimizes the total power function, and the minimum total power can be obtained by differentiating the total power with respect to the neutral angle and setting it equal to zero.

[0030] MERGEFORMAT (18) Obtain the rolling force model; When the total power is minimized, the calculated neutral angle is also the true neutral angle during the rolling process. After obtaining the minimum total power and the actual neutral angle, W p W f and W s The sum of these values ​​is the minimum. Therefore, the rolling torque M can be obtained: MERGEFORMAT (19) The rolling force is: MERGEFORMAT (20) Furthermore, after obtaining the rolling force, the calculation of the roll radius is also included. The iterative calculation process includes: The actual radius of the roll is calculated based on the rolling force, the roll Poisson's ratio, and the roll's elastic modulus. Substitute the updated actual roll radius into the steps of establishing the roll contact arc curve, re-establish the roll contact arc curve and iteratively calculate the rolling force until the convergence condition is met.

[0031] Furthermore, regarding the iterative calculation of the roll radius in this embodiment, during the rolling process, the roll undergoes significant elastic flattening, requiring iterative calculation. The actual radius of the roll during the rolling process is... : MERGEFORMAT (21) In the formula v roll E is the Poisson's ratio of the rolls. roll This refers to the elastic modulus of the roll.

[0032] Because the roll radius and rolling force are coupled, the roll radius needs to be iteratively calculated using the rolling force model until the required accuracy is met. In this embodiment, the iteration stops when the absolute value of the difference between the roll radius values ​​in two consecutive iterations is less than a certain value.

[0033] The convergence condition is: MERGEFORMAT (22) In the formula For the first The roll radius value at the next iteration; For the first The roll radius value at the next iteration.

[0034] This embodiment overcomes the limitations of the uniform deformation assumption in traditional rolling theory. By accurately describing the spatial non-uniform distribution characteristics of the velocity field within the deformation zone, it significantly improves the prediction accuracy of the rolling force model and effectively reduces the calculation error of loads under complex working conditions. This technical solution not only achieves a refined characterization of metal flow behavior and accurately captures velocity gradient differences in the thickness and width directions, but also provides a reliable theoretical basis for shape control and process optimization, balancing high-precision calculation with good engineering applicability.

[0035] To further optimize the scheme, as a preferred implementation method, the rolling force model in this embodiment is illustrated using PDA data from a 2180mm cold rolling production line (5-stand continuous rolling) at a steel plant. Taking cold-rolled commercial steel sheet (SPCC) as an example, the strip width is 1313mm, the material thickness at the S1 stand entrance is 2.45mm, and after passing through 5 mills, the final thickness at the exit is 0.59mm. Specific parameters for each stand are shown in Table 1.

[0036] Table 1 Taking the process parameters of the first rack as an example, the detailed calculation steps are as follows: Based on the process parameters of the first stand in cold rolling mill, the strip's inlet thickness 2h1 = 2.45 mm, outlet thickness 2h0 = 1.63 mm, inlet width B = 1313 mm, and pretension σ are determined. f =91MPa and back tension σ b =55.00MPa, initial roll radius R0=258mm, roll linear velocity v r =4.66m / s.

[0037] A two-dimensional schematic diagram of one-quarter of the cold-rolled deformation zone in this embodiment is shown below. Figure 1 As shown, let the x and y axes represent the length and thickness directions of the cold-rolled strip, respectively. The initial radius of the roll is R0, the slab inlet thickness is 2h1, the outlet thickness is 2h0, and θ is the angle between the line connecting any point in the deformation zone and the center of the roll. Based on the boundary conditions and geometric characteristics of the strip deformation zone, the roll contact arc curve and its first derivative equation are established: MERGEFORMAT (23) MERGEFORMAT (24) MERGEFORMAT (25) Establish a non-uniform velocity field and strain velocity field in the rolling deformation zone that meet the motion permitting conditions.

[0038] MERGEFORMAT (26) MERGEFORMAT (27) MERGEFORMAT (28) MERGEFORMAT (29) MERGEFORMAT (30) Based on the non-uniform velocity field and strain rate, the total power functional is obtained by calculating the internal deformation power, friction power, shear power, and tension power. MERGEFORMAT (31) According to different neutral angles θ n The corresponding total power functional is then used to obtain its minimum value: MERGEFORMAT (32) θ was calculated n =0.012.

[0039] Based on the coupling between rolling force and roll flattening radius, the rolling force meeting the convergence condition is calculated through iterative calculation and compared with the measured value, as shown in the figure. Figure 2 As shown in the figure. A comparison of the calculated and measured rolling force values ​​for the five stands is shown in Table 2.

[0040] Table 2 As can be seen from Table 2, the maximum error of the model in this invention is 4.78%, and the average error is 3.92%. In summary, the rolling force model in this embodiment fully meets the accuracy requirements of industrial applications. This result fully verifies that the rolling force model based on a non-uniform velocity field proposed in this paper has higher accuracy and reliability, and can more realistically reflect the physical nature of metal flow within the rolling deformation zone. It has a positive promoting effect on the optimization and control of the rolling process.

[0041] Example 2 Based on the same inventive concept, this embodiment also provides a rolling force prediction system for non-uniform speed fields in cold rolling of strip steel, including: The parameter input module is used to obtain the process parameters and strip parameters of cold continuous rolling; The geometric modeling module, connected to the parameter input module, is used to establish the roll contact arc curve and its first derivative equation based on process parameters and strip parameters. The velocity field construction module, connected to the geometric modeling module, is used to establish the non-uniform velocity field and strain velocity field of the rolling deformation zone based on the roll contact arc curve. The power calculation module, connected to the velocity field construction module, is used to calculate various power terms in the plastic deformation zone of cold-rolled slabs based on non-uniform velocity fields and strain velocity fields. The optimization solution module is connected to the power calculation module. It is used to construct the total power functional based on multiple power terms and to minimize the total power functional to obtain the neutral angle. The rolling force calculation module, connected to the optimization solution module, is used to calculate the rolling torque and rolling force based on the minimum total power functional and the neutral angle.

[0042] Furthermore, the process parameters and strip parameters acquired by the parameter input module include inlet thickness, outlet thickness, inlet width, initial radius of the roll, linear speed of the roll, front tension, and back tension.

[0043] Furthermore, the power calculation module includes: Internal deformation power calculation unit, used to calculate internal deformation power based on non-uniform velocity field and strain velocity field; Friction power calculation unit, used to calculate friction power based on non-uniform velocity field and friction conditions of workpiece surface; The shear power calculation unit is used to calculate the shear power based on the shear force at the velocity discontinuity at the entry point of the workpiece. The tension power calculation unit is used to calculate the tension power based on the influence of the front tension and the back tension on the rolling process.

[0044] Furthermore, the system also includes an iterative correction module, which is connected to the parameter input module and the rolling force calculation module respectively. The iterative correction module is used to calculate the actual roll radius based on the rolling force, the roll Poisson's ratio and the roll elastic modulus, and feed the updated roll radius back to the geometric modeling module. The geometric modeling module re-establishes the roll contact arc curve based on the updated roll radius and iteratively calculates the rolling force until the convergence condition is met.

[0045] The strip cold rolling non-uniform velocity field rolling force prediction system provided in this embodiment has all the advantages of the strip cold rolling non-uniform velocity field rolling force prediction method provided in Embodiment 1.

[0046] Example 3 This embodiment also discloses a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method described in Embodiment 1.

[0047] Example 4 This embodiment also discloses a computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the steps of the method described in Embodiment 1.

[0048] Example 5 This embodiment also discloses a computer program product, including a computer program that, when executed by a processor, implements the steps of the method described in Embodiment 1.

[0049] The above are merely preferred embodiments of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for predicting rolling force in a non-uniform velocity field during cold rolling of strip steel, characterized in that, include: Obtain the process parameters and strip parameters for cold continuous rolling; Based on the process parameters and strip parameters, establish the roll contact arc curve and its first derivative equation; Based on the roll contact arc curve, establish the non-uniform velocity field and strain velocity field of the rolling deformation zone; Based on the non-uniform velocity field and strain velocity field, calculate various power terms in the plastic deformation zone of the cold-rolled slab. Based on the various power terms, a total power functional is constructed and the total power functional is minimized to obtain the neutral angle; The rolling torque and rolling force are calculated based on the minimum total power functional and the neutral angle.

2. The method according to claim 1, characterized in that, The process of obtaining the process parameters and strip parameters for cold continuous rolling includes: Obtain the inlet thickness, outlet thickness, inlet width, initial roll radius, roll linear speed, front tension, and back tension.

3. The method according to claim 1, characterized in that, The process of establishing the roll contact arc curve and its first derivative equation includes: Based on the boundary conditions and geometric characteristics of the strip deformation zone, the roll contact arc curve and its first derivative equation are established.

4. The method according to claim 1, characterized in that, The process of establishing the non-uniform velocity field and strain velocity field in the rolling deformation zone includes: Based on the velocity boundary conditions of the strip rolling deformation zone, the principle of constant volume, and the roll contact arc curve, the non-uniform velocity field and strain velocity field of the rolling deformation zone are established.

5. The method according to claim 1, characterized in that, The process of calculating the various power terms in the plastic deformation zone of cold-rolled slab includes: Calculate the internal deformation power based on the non-uniform velocity field and strain velocity field; The frictional power is calculated based on the non-uniform velocity field and the frictional conditions of the rolled piece surface. Calculate the shear power based on the shear force at the velocity discontinuity at the entry point of the workpiece; Calculate the tension power based on the influence of front tension and back tension on the rolling process.

6. The method according to claim 1, characterized in that, After obtaining the rolling force, the method also includes iterative calculation of the roll radius; The iterative calculation process includes: The actual radius of the roll is calculated based on the rolling force, Poisson's ratio of the roll, and the elastic modulus of the roll. Substitute the updated actual roll radius into the steps of establishing the roll contact arc curve, re-establish the roll contact arc curve and iteratively calculate the rolling force until the convergence condition is met.

7. A rolling force prediction system for non-uniform velocity fields in cold rolling of strip steel, characterized in that, include: The parameter input module is used to obtain the process parameters and strip parameters of cold continuous rolling; The geometric modeling module, connected to the parameter input module, is used to establish the roll contact arc curve and its first derivative equation based on the process parameters and strip parameters. A velocity field construction module, connected to the geometric modeling module, is used to establish a non-uniform velocity field and strain velocity field in the rolling deformation zone based on the roll contact arc curve. The power calculation module, connected to the velocity field construction module, is used to calculate various power terms in the plastic deformation zone of the cold-rolled slab based on the non-uniform velocity field and strain velocity field. An optimization solution module, connected to the power calculation module, is used to construct a total power functional based on the various power terms, and to minimize the total power functional to obtain the neutral angle. The rolling force calculation module, connected to the optimization solution module, is used to calculate the rolling torque and rolling force based on the minimum total power functional and the neutral angle.

8. The system according to claim 7, characterized in that, The process parameters and strip parameters acquired by the parameter input module include inlet thickness, outlet thickness, inlet width, initial radius of the roll, linear speed of the roll, front tension, and back tension.

9. The system according to claim 7, characterized in that, The power calculation module includes: An internal deformation power calculation unit is used to calculate the internal deformation power based on the non-uniform velocity field and strain velocity field. Friction power calculation unit, used to calculate friction power based on the non-uniform velocity field and friction conditions of the workpiece surface; The shear power calculation unit is used to calculate the shear power based on the shear force at the velocity discontinuity at the entry point of the workpiece. The tension power calculation unit is used to calculate the tension power based on the influence of the front tension and the back tension on the rolling process.

10. The system according to claim 7, characterized in that, The system also includes an iterative correction module, which is connected to the parameter input module and the rolling force calculation module respectively. The iterative correction module is used to calculate the actual roll radius based on the rolling force, the roll Poisson's ratio and the roll elastic modulus, and feed the updated roll radius back to the geometric modeling module. The geometric modeling module re-establishes the roll contact arc curve based on the updated roll radius and iteratively calculates the rolling force until the convergence condition is met.

Images