Steel enterprise panoramic intelligent collaborative management platform based on industrial big data

Abstract

The application is steel enterprise panoramic intelligent collaborative management platform based on industrial big data, relates to the technical field of industrial big data, and comprises the following modules: a multi-dimensional physical quantity acquisition and initial boundary matrix construction module, a transient heat conduction calculation engine module based on a partial differential equation, a stacking topology optimization and total enthalpy value evaluation module, and a physical scheduling instruction conversion and equipment collaboration module.In the application, array type multi-physical quantity acquisition, millisecond level time synchronization, steel grade temperature segmented thermal property interpolation, and unified three-dimensional space coordinate system are used to construct a global initial boundary matrix that can be directly used for numerical solution; relying on a three-dimensional transient heat conduction partial differential equation and a finite difference method, a multi-physical field boundary of contact thermal resistance, natural convection, and thermal radiation is coupled to realize high-precision quantitative calculation of the three-dimensional temperature field, temperature drop rate, and temperature gradient of the steel billet at the target extraction moment.

Description

Technical Field

[0001] This invention relates to the field of industrial big data technology, and in particular to a panoramic intelligent collaborative management platform for steel enterprises based on industrial big data. Background Technology

[0002] Currently, the continuous casting billet storage process has significant shortcomings in terms of comprehensive physical quantity perception and the construction of initial thermodynamic boundaries. Traditional storage methods only employ single-point, discontinuous temperature detection and lack an array-based multi-dimensional physical acquisition system, making it impossible to obtain the temperature field and refined spatial coordinates of the entire billet surface. Thermophysical parameters are mostly fixed constants, without segmented modeling and interpolation correction according to steel grade and temperature. The spatial coordinate system is not consistent with the storage equipment benchmark, making it difficult to construct a standardized initial boundary matrix that can be directly used for numerical solutions. This results in a lack of accurate and unified input foundation for subsequent thermodynamic calculations.

[0003] Existing technologies generally lack the ability to accurately solve transient heat conduction based on partial differential equations. A three-dimensional transient heat conduction numerical model for high-temperature steel billets has not been established, and multi-physics boundary coupling of contact thermal resistance, natural convection, and thermal radiation has not been achieved. The finite difference discretization scheme, time-stepping iteration, and sparse matrix solution mechanism are imperfect, resulting in low accuracy and poor stability of temperature field calculations. It is impossible to quantitatively predict the three-dimensional temperature distribution, temperature drop rate, and temperature gradient of the steel billet at any extraction time, making it difficult to support the accurate assessment of the thermal state of the steel billet.

[0004] Therefore, a panoramic intelligent collaborative management platform for steel enterprises based on industrial big data is proposed to address the aforementioned issues. Summary of the Invention

[0005] The purpose of this invention is to propose a panoramic intelligent collaborative management platform for steel enterprises based on industrial big data in order to solve the above-mentioned problems.

[0006] To achieve the above objectives, the present invention adopts the following technical solution: A comprehensive intelligent collaborative management platform for steel enterprises based on industrial big data includes: The multidimensional physical quantity acquisition and initial boundary matrix construction module is configured to construct a global initial boundary matrix for thermodynamic numerical solutions by acquiring billet temperature, thermal properties, spatial coordinates and storage environment parameters. The transient heat conduction calculation engine module based on partial differential equations is configured to solve the three-dimensional temperature field of the billet at the time of target extraction by coupling contact thermal resistance and convection-radiation boundary based on the three-dimensional transient heat conduction equation. The stacking topology optimization and total enthalpy evaluation module is configured to maximize the total sensible enthalpy as the optimization objective, screen the optimal stacking topology, and prevent billet cooling cracking by using the phase transformation temperature threshold. The physical scheduling instruction conversion and equipment coordination module is configured to convert the optimal stacking topology into the three-dimensional coordinates and timing action instructions of the unmanned crane, and realize the equipment execution through industrial communication.

[0007] Preferably, the multidimensional physical quantity acquisition and initial boundary matrix construction module specifically includes: Real-time acquisition of multi-point temperature values ​​from the upper surface, lower surface, and four sides of the steel billet, synchronously latching the timestamps of the PLC system, and recording the moment the steel billet cutting is completed. Timing of billet unloading , When entering the slab warehouse ; The system has a built-in database and high storage density. Specific heat capacity at constant pressure Thermal conductivity ; Based on the billet grade and the current average temperature, linear interpolation is used to calculate the temperature at that point. , , value; Using one corner of the slab warehouse as the origin of the coordinate system The length direction is Axial and width directions are Axis and height directions are Establish an absolute three-dimensional Cartesian coordinate system based on the axes; Upon entering the warehouse, steel billets are automatically matched and assigned to corresponding grid cells based on their actual length, width, and height. The set of grid IDs and center coordinates of each steel billet are recorded. Boundary coordinate range, total number of grid cells occupied; The initial temperature field of a single steel billet , , , Spatial grid range, timestamp encapsulated as The initial matrix is ​​a 6-dimensional matrix, and the matrices of all steel billets are summed to form the global initial boundary matrix.

[0008] Preferably, the method further includes: Temperature and humidity sensors and wind speed sensors were evenly distributed within the slab storage area, positioned at the four corners of the storage area, to collect ambient temperature data. air velocity ; convective heat transfer coefficient The calculations were performed using natural convection correlations. Radiative heat transfer The calculations are performed using the Stefan-Boltzmann law; The final output is a third type of boundary condition that couples convective and radiative heat dissipation. The numerical input is used to the computation engine, and the boundary condition is directly substituted into the discrete equation in the form of nodal heat flux density.

[0009] Preferably, the transient heat conduction calculation engine module based on partial differential equations specifically includes: Three-dimensional transient heat conduction equation: ; in, For temperature; For time; Finite Difference Method (FDM) Discretization: Spatial discretization: using a central difference scheme for... Discretization of directional second-order partial derivatives; Time discretization: Implicit Euler scheme is used, with a time step size of [missing information]. fixed; Grid node numbering: Assigning grid nodes to each steel billet internally Three-dimensional indexing, establishing discrete algebraic equations node by node; Calculation of contact thermal resistance and multi-boundary coupling: Contact thermal resistance The value is determined based on the surface roughness of the steel billet, the thickness of the oxide scale, and the stacking pressure; It automatically matches the corresponding contact thermal resistance based on the number of stacked layers. Numerical value.

[0010] Preferably, the process further includes boundary partitioning: Steel billet-steel billet contact surface: thermally conductive boundary, considering contact thermal resistance, and substituted into the three-dimensional transient heat conduction finite difference discretization equation at the boundary grid node in the form of heat flux density; Steel billet-air exposed surface: convection and radiation coupled boundary, direct calculation of nodal heat flux; Steel billet bottom-pad beam contact surface: fixed contact thermal resistance boundary, and the thermal resistance value is uniformly fixed; All boundary conditions are transformed into nodal temperature algebraic equations and substituted into the global discrete equation system. Time step iteration and target time solution: Inputs: Initial temperature field, property matrix, boundary conditions, planned extraction time. ; Iterative process: from the moment of entering the slab warehouse Beginning, with Using a time step size, solve the linear equations step by step and update the temperature of all nodes. Convergence criterion: The current step is considered converged when the maximum temperature change in a single step meets the preset requirements; if it does not converge, the iteration margin is automatically reduced and the calculation is repeated until the conditions are met. Output: Arrived At time t, output the three-dimensional temperature field of each steel billet. Core temperature, surface temperature, maximum temperature difference, temperature gradient, and cooling rate.

[0011] Preferably, the stacking topology optimization and total enthalpy evaluation module specifically includes: Total enthalpy formula: ; Numerical integration method: The three-dimensional trapezoidal integration method is used to integrate the volume of each grid cell occupied by the steel billet. ; For a single grid volume ( ); Multi-constraint deterministic topology optimization: Topology scheme generation: Enumerate all stacking methods that satisfy physical constraints, and the system generates candidate schemes according to the thermal insulation logic; Constraints include equipment constraints, safety constraints, and process constraints; A heuristic traversal optimization method is used to calculate the enthalpy value of all feasible solutions and directly select the appropriate one. The largest possible solution is the optimal topology.

[0012] Preferably, the method further includes phase change threshold interception and security determination: extract Lowest local temperature of steel billet at any time ,like If the cooling rate exceeds the material's critical cooling rate, it is also considered a risky solution. Automatically reject the proposed solution and recalculate the temperature field until... And the cooling rate meets the standard; Output the billet number, target coordinates, stacking relationship, insulation method, expected enthalpy value, temperature index, and safety judgment result of the optimal topology to form a unique executable solution.

[0013] Preferably, the physical scheduling instruction conversion and device coordination module specifically includes: The theoretical coordinates of thermodynamic topology are converted into absolute pulse coordinates that can be recognized by the overhead crane PLC using rigid body coordinate transformation formulas. The action commands are broken down into transfer commands, clamping commands, and positioning commands. Stacking action sequence generation: Generate an action sequence with absolute timestamps according to preset principles.

[0014] Preferably, the method further includes: The system continuously compares the crane's return to position signal, action completion signal, fault signal, and current coordinates with the optimal topology in real time to determine whether the execution result is consistent with the optimal topology. If the deviation exceeds the preset allowable range, the execution is judged to be abnormal. After the exception handling action is completed, the new coordinates and temperature field of the billet are updated, and the initial boundary matrix is ​​written synchronously as the initial condition for the next scheduling.

[0015] In summary, due to the adoption of the above technical solution, the beneficial effects of the present invention are: 1. This invention constructs a global initial boundary matrix that can be directly used for numerical solution by array-type multi-physical quantity acquisition, millisecond-level time synchronization, segmented thermal property interpolation of steel temperature, and a unified three-dimensional spatial coordinate system; relying on the three-dimensional transient heat conduction partial differential equation and the finite difference method, it couples the multi-physical field boundaries of contact thermal resistance, natural convection, and thermal radiation to achieve high-precision quantitative calculation of the three-dimensional temperature field, temperature drop rate, and temperature gradient of the billet at the target extraction time.

[0016] 2. This invention, under multiple constraints of equipment, safety, and process, employs a combination of heuristic traversal and genetic algorithms to quickly output the optimal stacking scheme; it incorporates a built-in martensitic phase transformation temperature threshold protection mechanism to determine the risk of cooling cracking in real time and automatically iteratively optimize the insulation topology, eliminating billet quality defects from the source; through precise coordinate mapping and standardized instruction decomposition, the thermodynamically optimal topology is transformed into a time-sequence action that can be executed by the unmanned overhead crane, relying on industrial communication to achieve a complete closed loop of instruction issuance, trajectory planning, collision avoidance, and execution feedback; the solution eliminates manual experience-based scheduling, improves scheduling efficiency and equipment automation level, reduces safety hazards and quality fluctuations, and helps steel enterprises achieve panoramic intelligent collaboration and quality improvement in the continuous casting billet storage process. Attached Figure Description

[0017] Further details, features, and advantages of this application are disclosed in the following description of exemplary embodiments in conjunction with the accompanying drawings, in which: Figure 1 This is a system structure diagram of the present invention. Detailed Implementation

[0018] Several embodiments of this application will now be described in more detail with reference to the accompanying drawings to enable those skilled in the art to implement this application. This application may be embodied in many different forms and for various purposes and should not be limited to the embodiments set forth herein. These embodiments are provided to make this application thorough and complete, and to fully convey the scope of this application to those skilled in the art. The embodiments described do not limit this application.

[0019] Unless otherwise defined, all terms used herein (including technical and scientific terms) shall have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains. It will be further understood that terms such as those defined in commonly used dictionaries shall be interpreted as having a meaning consistent with their meaning in the relevant field and / or the context of this specification, and shall not be interpreted in an idealized or overly formal sense unless expressly defined herein.

[0020] Example 1

[0021] Its specific implementation method is combined with the appendix Figure 1 Please provide a detailed explanation.

[0022] Appendix Figure 1 The structural block diagram of the panoramic intelligent collaborative management platform for steel enterprises based on industrial big data provided in the embodiments of the present invention shows the connection relationship between the multi-dimensional physical quantity acquisition and initial boundary matrix construction module and the physical scheduling instruction conversion and equipment coordination module, and marks the main functional interaction flow of each module.

[0023] In this embodiment, it includes: The multidimensional physical quantity acquisition and initial boundary matrix construction module is configured to construct a global initial boundary matrix for thermodynamic numerical solutions by acquiring billet temperature, thermal properties, spatial coordinates and storage environment parameters. Specifically, it includes: This module is responsible for converting the billet and storage environment into standardized physical parameters that can be directly input into partial differential equations.

[0024] Initial temperature field acquisition and mapping: Hardware data acquisition: An infrared temperature array sensor (non-contact, temperature accuracy ±1℃, temperature range 300℃~1200℃) is connected to the output roller conveyor of the continuous casting machine. ≥8 temperature measurement points are evenly arranged along the length of the billet and ≥4 temperature measurement points are evenly arranged along the width. The acquisition frequency is 10Hz, acquiring multi-point temperature values ​​of the upper surface, lower surface, and four sides of the billet in real time. The sensor and roller conveyor encoder are linked in real time to ensure automatic triggering of data acquisition when the billet reaches the temperature measurement area, avoiding blank sampling and missed data acquisition. A complete surface temperature matrix is ​​generated for each billet.

[0025] Time synchronization: When acquiring temperature data, the PLC system timestamp is simultaneously latched to accurately record the moment when the billet cutting is completed. Timing of billet unloading , When entering the slab warehouse Each timestamp is accurate to the millisecond level and is used for subsequent initial calculations of temperature decay, ensuring that the time base for thermodynamic calculations is completely consistent.

[0026] Thermophysical property parameter acquisition: The system has a built-in thermophysical property database segmented by steel grade and temperature, storing information including carbon equivalent, alloy composition, and density for different grades such as Q235, Q355, 45#, high-strength structural steel, and alloy steel. Specific heat capacity at constant pressure Thermal conductivity Piecewise functions / discrete data tables that vary with temperature; the database is segmented and stored in 50°C intervals, covering the entire range of 300°C to 1200°C.

[0027] The system uses linear interpolation based on the billet grade and current average temperature to calculate the accurate temperature at a given point. , , The values ​​are not fixed constants, ensuring that the calculation benchmark is unique and reproducible. If the temperature falls between two segmented intervals, bilinear interpolation is automatically performed to ensure that the physical property parameters remain continuous without abrupt changes.

[0028] 3D spatial mesh generation and property matrix construction: Definition of spatial coordinate system: The origin is a corner of the slab warehouse. The length direction is Axial and width directions are Axis and height directions are A three-dimensional Cartesian coordinate system is established, with coordinate units in mm and a positioning accuracy of ±5 mm. This coordinate system shares the same reference with the unmanned overhead crane and warehouse positioning system, eliminating coordinate transformation bias.

[0029] Mesh discretization rule: Distribute the entire slab library according to... 500mm axis step size 500mm axis step size The axis step size is 300mm, which is divided into uniform three-dimensional mesh units. Each mesh unit is assigned a unique ID number. For steel billets with irregular cross sections, the system automatically and adaptively adjusts the local mesh density to ensure that the mesh completely fits the outline of the steel billet.

[0030] Upon entering the warehouse, steel billets are automatically matched and assigned to corresponding grid cells based on their actual length, width, and height. The set of grid IDs and center coordinates of each steel billet are recorded. By defining the boundary coordinate range and the total number of grid cells occupied, the spatial location of the billet can be accurately expressed digitally.

[0031] Initial property matrix generation: The initial temperature field of a single steel billet , , , Spatial grid range, timestamp encapsulated as A 6-dimensional initial matrix is ​​generated, and the matrices of all steel billets are aggregated into a global initial boundary matrix, stored in a floating-point array format, and directly used as input to the computation engine. Each row of the matrix corresponds to a grid node, and each column corresponds to a physical parameter, which can be directly read and called by the numerical solution module.

[0032] Real-time collection and calculation of warehouse environment boundary conditions: Environmental data collection: ≥4 sets of temperature and humidity sensors and wind speed sensors are evenly distributed within the slab warehouse, located at the four corners of the warehouse area. Data collection range: Ambient temperature (-10℃~60℃, accuracy ±0.5℃), airflow rate (0~5m / s, accuracy ±0.1m / s), sampling frequency 1Hz.

[0033] The system performs arithmetic averaging on multiple sets of sensor data to eliminate errors caused by local environmental fluctuations and obtain globally unified environmental boundary parameters.

[0034] convective heat transfer coefficient The calculations were performed using natural convection correlations: ; in: It is a constant (0.54 for horizontal surfaces and 0.59 for vertical surfaces). The thermal conductivity of air, For Grashof numbers, For Prandtl numbers, The characteristic length of the steel billet; The system has a built-in air property parameter table that automatically retrieves parameters based on the ambient temperature. , The numerical value is automatically calculated in real time based on the billet's placement posture using the corresponding formula. The value is calculated and the result is rounded to three decimal places.

[0035] Radiative heat transfer Calculation: Using the Stefan-Boltzmann law: ; in: The surface emissivity of the steel billet is 0.8~0.9 (for oxidized surfaces). For Boltzmann constant, The surface temperature of the steel billet is used as the final output. The third type of boundary condition, which is coupled with convective and radiative heat dissipation, is numerically input into the calculation engine. The boundary condition is directly substituted into the discrete equation in the form of nodal heat flux density (the discrete equation corresponds to the three-dimensional transient heat conduction equation of the transient heat conduction calculation engine module based on partial differential equations).

[0036] Data validation and fault tolerance: If the sensor data acquisition is abnormal (overheating, disconnection, data jump), the system will automatically use the valid data from the previous moment to interpolate and complete the data. If the acquisition is abnormal for 3 consecutive times, an alarm will be triggered and the system will automatically switch to redundant sensor data. If all sensors fail, the system will automatically call the historical average environmental value for the same period to ensure that the boundary conditions are continuous and uninterrupted, and will not affect the subsequent calculation process.

[0037] The transient heat conduction calculation engine module based on partial differential equations is configured to use the three-dimensional transient heat conduction equations, coupled with contact thermal resistance and convection-radiation boundary, and iteratively solve the three-dimensional temperature field of the billet at the target extraction time using the finite difference method. Specifically, it includes: Three-dimensional transient heat conduction equation (isotropic, no internal heat source): ; in, Temperature, the temperature field inside the material, is represented by spatial coordinates. and time function ; For time, the time variable in the heat conduction process; Physical meaning: The rate of change of internal energy of a steel billet per unit volume = heat flux density divergence in three directions, which fully describes the internal heat conduction process of the steel billet; Applicable conditions: solid steel billet, no phase change (before phase change), no flow, no internal heat source, applicable to the full temperature range cooling process of high-temperature steel billets after continuous casting.

[0038] The rate of change of internal energy per unit volume of material over time (i.e., heat storage rate). Within a unit volume, through , , Net heat transfer rate in three directions (inflow / outflow); Finite Difference Method (FDM) Discretization: Spatial discretization: using a central difference scheme for... Directional second-order partial derivative discretization: ; in , The average thermal conductivity between nodes is calculated using an arithmetic mean, ensuring the accuracy of thermal conductivity calculations for heterogeneous materials. For steel billets Grid step size in direction (spatial distance between adjacent nodes); Time discretization: Implicit Euler scheme (unconditionally stable), time step size A fixed step size of 1 second is used to avoid calculation divergence; for long-term prediction scenarios, the step size can be automatically increased but not exceeding 5 seconds, while ensuring calculation accuracy.

[0039] Grid node numbering: Assigning grid nodes to each steel billet internally A three-dimensional index is used to establish discrete algebraic equations for each node. All node equations are aggregated to form a global sparse linear equation system, and the size of the equation system corresponds exactly to the number of grids.

[0040] Calculation of contact thermal resistance and multi-boundary coupling: Contact thermal resistance Values: The system has a built-in database of billet-to-billet contact thermal resistance, and the value range is determined based on billet surface roughness, oxide scale thickness, and stacking pressure. The greater the stacking pressure, the larger the contact area. The smaller; The system automatically matches the corresponding contact thermal resistance based on the number of stacked layers. Numerical value, interfacial heat flux density: , , These represent the temperatures at the interface of the two contacting steel billets; this achieves heat flow coupling between the two billets, ensuring continuous heat transfer at the interface.

[0041] Boundary partitioning: Billet-to-billet contact surface: thermally conductive boundary. Considering contact thermal resistance, the heat flux density is substituted into the three-dimensional transient heat conduction finite difference discretization equation at the boundary mesh nodes. The heat flux density generated by the contact thermal resistance is then used to solve the problem. As a known term, the internal heat conduction term in the corresponding direction in the original equation is replaced, thereby completing the numerical application of the heat conduction boundary condition. The process is as follows: by Taking direction as an example, internal nodes The discrete equation (derived from the central difference) is: ; This equation describes the rate of temperature change of the internal nodes, equal to the rate of temperature change from... , , Net heat flow in three directions.

[0042] Then, substitute the boundary heat flux density into the discrete equations of the boundary nodes: When node When located on the billet-to-billet contact surface, it is no longer an internal node, but a boundary node. At this point, The directional thermal conductivity will be affected by the heat flux density generated by the contact thermal resistance. The boundary node discrete equations, after substitution, become: ; Among them, heat flux density The contact thermal resistance formula gives: ; Billet-air exposed surface: convection and radiation coupling boundary. , The natural convection heat transfer coefficient is calculated from the ambient temperature, air velocity, and billet posture. The temperature of the exposed surface of the steel billet. The temperature of the air in the storage environment. The emissivity of the steel billet surface (typically 0.8~0.9 for anodized surfaces) is given. Using the Stefan-Boltzmann constant, the heat flux of nodal heat dissipation is directly calculated; The contact surface between the bottom of the billet and the support beam has a fixed thermal resistance boundary, and the thermal resistance value is uniformly set to a fixed value. ; All boundary conditions are transformed into algebraic equations of nodal temperatures, which are then substituted into the global discrete equation system without altering the structure of the equation system, thus ensuring the stability of the solution.

[0043] Time step iteration and target time solution: Inputs: Initial temperature field, property matrix, boundary conditions, planned extraction time. (Issued by the steel rolling process, accurate to the second); Iterative process: from the moment of entering the slab warehouse Beginning, with The linear equations are solved step by step, with a step size of 1, and the temperature of all nodes is updated. The solver adopts the catch-up method and is optimized for sparse matrices to improve computational efficiency. Convergence criterion: The current step is considered converged when the maximum temperature change in a single step meets the preset requirement (<0.1℃); if it does not converge, the iteration margin is automatically reduced and the calculation is repeated until the condition is met. Output: Arrived At time t, output the three-dimensional temperature field of each steel billet. Core temperature, surface temperature, maximum temperature difference, temperature gradient, and cooling rate are all stored as grid node data for subsequent optimization modules to access.

[0044] Calculate fault tolerance: If the number of iterations exceeds the limit or the system of equations is singular, the system will automatically adjust the time step slightly. (Recalculate); if multiple calculations fail, automatically simplify the boundary conditions (ignore radiation and retain only heat conduction and convection) to ensure that the results are output, and record the calculation mode for subsequent verification.

[0045] The stacking topology optimization and total enthalpy evaluation module is configured to maximize the total sensible enthalpy as the optimization objective, select the optimal stacking topology under equipment, safety and process constraints, and prevent billet cooling cracking by using the phase transformation temperature threshold. Specifically, it includes: Calculation of total sensible enthalpy: Total enthalpy formula: ; For the entire volume occupied by the steel billet Perform a triple integral to accumulate the enthalpy values ​​of all infinitesimal elements; Numerical integration method: The three-dimensional trapezoidal integration method is used to integrate the volume of each grid cell occupied by the steel billet. ; For a single grid volume ( During integration, the results are accumulated grid by grid, and empty grids not occupied by steel billets are automatically excluded to ensure accurate integration results.

[0046] Calculation range: enthalpy values ​​for a single stack or the entire warehouse can be selected. The larger the value, the less heat is lost, making it the optimal core indicator; the system also calculates the average temperature and temperature drop as auxiliary evaluation indicators.

[0047] Multi-constraint deterministic topology optimization: Topology scheme generation: Enumerate all stacking methods that meet physical constraints, including stacking layers of 1 to 6, side by side, sandwich, triangular, and high and low temperature combination stacking; the system prioritizes generating candidate schemes based on the insulation logic of high temperature inside, low temperature outside, high temperature on top, and low temperature on bottom, reducing the calculation of invalid schemes.

[0048] The constraints include: Equipment constraints: The maximum load capacity of the overhead crane shall be less than or equal to the rated load (e.g., 30t). The operating distance does not exceed the warehouse body, and a single hook can only lift one steel billet at a time; Safety constraints: stacking height-to-width ratio ≤ 2.0, stacking center of gravity offset ≤ 100mm, interlayer misalignment not exceeding 1 / 5 of the billet width; Process constraints: first-in, first-out; rolling time window error ≤ 10 min; steel billets of the same order should be stacked together. Optimization Algorithm: A heuristic traversal optimization is used to calculate the enthalpy value of all feasible solutions and directly select the optimal one. The largest solution is the optimal topology; if the number of solutions is greater than 100, a genetic algorithm is used to accelerate the search for the optimal solution. The population size is 50, the number of iterations is 20, the crossover probability is 0.7, the mutation probability is 0.1, and the best individual is retained in each generation to ensure that the solution converges quickly to the global optimum.

[0049] Scheme sorting: Sort all feasible schemes from highest to lowest enthalpy value to generate alternative schemes, which can be directly switched when the main scheme fails.

[0050] It also includes phase change threshold interception and security determination: Phase transformation database: The system contains the critical points of martensitic phase transformation for various steel grades. The thermometer is calibrated based on material test data, such as high-strength steel. carbon steel The database allows for precise searching by steel grade.

[0051] Judgment rule: Extract Lowest local temperature of steel billet at any time ,like If the cooling rate exceeds the material's critical cooling rate, it is also considered a risky solution. Processing logic: Automatically reject the proposed solution, forcibly replace it with a thermally insulated topology that uses high-temperature steel billets to cover low-temperature steel billets, reduces exposed area, increases stacking density, and ensures airtight stacking, and recalculates the temperature field until... And the cooling rate meets the standard; Solution output: The system outputs the billet number, target coordinates, stacking relationship, insulation method, expected enthalpy, temperature index, and safety judgment result of the optimal topology, forming a unique executable solution. The solution comes with complete calculation basis and is traceable and verifiable.

[0052] The physical scheduling instruction conversion and equipment coordination module is configured to convert the optimal stacking topology into the three-dimensional coordinates and timing action instructions of the unmanned crane, and realize precise execution and full-process control of the equipment through industrial communication; Specifically, it includes: 3D coordinate transformation and command generation: Coordinate mapping: The theoretical coordinates of the thermodynamic topology are converted into absolute pulse coordinates that can be recognized by the overhead crane PLC using rigid body coordinate transformation formulas. ; in The calibration coefficients for the overhead crane (determined by equipment installation and commissioning, and fixed in the system) are as follows: the origin pulse corresponds to the physical coordinate 0, the extreme position pulse corresponds to the physical maximum value, and the coefficients are obtained by linear fitting.

[0053] Action instructions are broken down as follows: Transfer command: Target coordinates and running speed (alignment at high speed 50m / min, low speed 5m / min, automatic deceleration when 500mm from the target position); Clamping commands: clamping pressure, opening degree, lifting height; clamping closing accuracy ±2mm. Positioning command: Positioning accuracy ±10mm, timeout judgment 30s, automatic retry if not in place within timeout.

[0054] Stacking operation timing locking and sequence output: Stacking action timing generation: According to the preset principle (far first then near, top first then bottom, no collision principle), generate an action sequence with absolute timestamp, with an action interval of ≥2s to avoid interference from multiple overhead cranes; the system automatically performs spatial collision detection, predicts the intersection of motion trajectories, and adjusts the timing to avoid interference areas.

[0055] Standard instruction example: Timestamp: 2026-03-02 10:00:00

[0056] Control the overhead crane to move from the origin to The clamps grip the steel billet #A (800℃) and lift it to... Run to Place it above steel billet #B (500℃), release the clamps, and return to the standby position; Instruction format: Standard industrial string instructions are used, including: action code, coordinates, speed, time, billet number, and check digit. The instruction length is fixed, which facilitates PLC parsing and execution.

[0057] Communication and closed-loop feedback: Communication protocol: Communicates with the overhead crane PLC via Profinet / ModbusTCP, with a data transmission frequency of 50Hz. Commands are retransmitted 3 times after a timeout, with each retransmission spaced 500ms apart. The communication message includes a command header, data area, and checksum to ensure error-free transmission.

[0058] Execution feedback: The system compares the overhead crane return to position signal, action completion signal, fault signal, and current coordinates in real time with the optimal topology to see if the execution result is consistent with the optimal topology. If the deviation exceeds the preset allowable range (±10mm), the execution is judged to be abnormal. Anomaly Handling: If execution fails (collision, timeout, fault), the system will automatically recalculate the current temperature field, generate a corrected topology and emergency instructions, prioritize the safety of billet insulation, and ensure that the process is not interrupted.

[0059] Execution closed loop: After the anomaly handling action is completed, the new coordinates and temperature field of the billet are updated and synchronously written into the initial boundary matrix as the initial conditions for the next scheduling, forming a complete and controllable closed loop of physical perception, numerical calculation, topology optimization, equipment execution, and data update. Each scheduling is based on the real-time physical state and there is no cumulative error.

[0060] The above formulas are all dimensionless calculations. The formulas are derived from software simulations based on a large amount of collected data to obtain the most recent real-world results. The preset parameters in the formulas are set by those skilled in the art according to the actual situation.

[0061] The foregoing has only described certain exemplary embodiments of the present invention by way of illustration. Undoubtedly, those skilled in the art can modify the described embodiments in various ways without departing from the spirit and scope of the present invention. Therefore, the foregoing drawings and descriptions are illustrative in nature and should not be construed as limiting the scope of protection of the claims of the present invention.

[0062] It should be noted that, in this document, the use of relational terms such as "first" and "second" is merely for distinguishing one entity or operation from another, and does not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes the element.

[0063] It should be understood that in the various embodiments of this application, the order of the above-mentioned processes does not imply the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application.

[0064] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0065] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0066] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0067] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0068] The above description is merely a specific embodiment 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.

[0069] The foregoing has only described certain exemplary embodiments of the present invention by way of illustration. Undoubtedly, those skilled in the art can modify the described embodiments in various ways without departing from the spirit and scope of the present invention. Therefore, the foregoing drawings and descriptions are illustrative in nature and should not be construed as limiting the scope of protection of the claims of the present invention.

Claims

1. Steel enterprise panoramic intelligent collaborative management platform based on industrial big data, characterized in that, include: The multidimensional physical quantity acquisition and initial boundary matrix construction module is configured to acquire billet temperature, thermal properties, spatial coordinates and storage environment parameters to construct a global initial boundary matrix for thermodynamic numerical solution; the global initial boundary matrix serves as a data carrier and is output to the transient heat conduction calculation engine module based on partial differential equations. The transient heat conduction calculation engine module based on partial differential equations is configured to receive the global initial boundary matrix and substitute it into the three-dimensional transient heat conduction equation, coupled with the contact thermal resistance and the convection-radiation boundary, to solve for the three-dimensional temperature field of the billet at the target extraction time; the three-dimensional temperature field is used to calculate the total sensible enthalpy of each billet, as well as to extract the lowest local temperature and cooling rate for safety determination. The stacking topology optimization and total enthalpy evaluation module is configured to maximize the total sensible enthalpy as the optimization objective, select the optimal stacking topology, and prevent the billet from cooling and cracking by using the phase transformation temperature threshold. The theoretical coordinates in the optimal stacking topology are converted into absolute pulse coordinates that can be recognized by the crane PLC through rigid body coordinate transformation, and used to generate transfer instructions, clamping instructions and positioning instructions. The physical scheduling instruction conversion and equipment coordination module is configured to convert the optimal stacking topology into the three-dimensional coordinates and timing action instructions of the unmanned crane, and realize the equipment execution through industrial communication. After the equipment execution is completed, the new coordinates and temperature field of the billet are updated and synchronously written into the initial boundary matrix as the initial conditions for the next scheduling.

2. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 1, characterized in that, The multidimensional physical quantity acquisition and initial boundary matrix construction module specifically includes: Real-time acquisition of multi-point temperature values ​​from the upper surface, lower surface, and four sides of the steel billet, synchronously latching the timestamps of the PLC system, and recording the moment the steel billet cutting is completed. Timing of billet unloading , When entering the slab warehouse ; The system has a built-in database and high storage density. Specific heat capacity at constant pressure thermal conductivity ; Based on the billet grade and the current average temperature, linear interpolation is used to calculate the temperature at that point. , , value; Using one corner of the slab warehouse as the origin of the coordinate system The length direction is Axial and width directions are Axis and height directions are Establish an absolute three-dimensional Cartesian coordinate system based on the axes; Upon entering the warehouse, steel billets are automatically matched and assigned to corresponding grid cells based on their actual length, width, and height. The set of grid IDs and center coordinates of each steel billet are recorded. Boundary coordinate range, total number of grid cells occupied; The initial temperature field of a single steel billet , , , Spatial grid range, timestamp encapsulated as The initial matrix is ​​a 6-dimensional matrix, and the matrices of all steel billets are summed to form the global initial boundary matrix.

3. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 2, characterized in that, Also includes: Temperature and humidity sensors and wind speed sensors were evenly distributed within the slab storage area, positioned at the four corners of the storage area, to collect ambient temperature data. air velocity ; convective heat transfer coefficient The calculations were performed using natural convection correlations. Radiative heat transfer The calculations are performed using the Stefan-Boltzmann law; The final output is a third type of boundary condition that couples convective and radiative heat dissipation. The numerical input is used to the computation engine, and the boundary condition is directly substituted into the discrete equation in the form of nodal heat flux density.

4. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 1, characterized in that, The transient heat conduction calculation engine module based on partial differential equations specifically includes: Three-dimensional transient heat conduction equation: ； in, For temperature; For time; Finite Difference Method (FDM) Discretization: Spatial discretization: using a central difference scheme for... Discretization of the second-order partial derivatives in the direction; Time discretization: Implicit Euler scheme is used, time step size fixed; Grid node numbering: Assigning grid nodes to each steel billet internally Three-dimensional indexing, establishing discrete algebraic equations node by node; Calculation of contact thermal resistance and multi-boundary coupling: Contact thermal resistance The value is determined based on the surface roughness of the steel billet, the thickness of the oxide scale, and the stacking pressure; It automatically matches the corresponding contact thermal resistance based on the number of stacked layers. Numerical value.

5. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 4, characterized in that, It also includes boundary partitioning: Steel billet-steel billet contact surface: thermally conductive boundary, considering contact thermal resistance, and substituted into the three-dimensional transient heat conduction finite difference discretization equation at the boundary grid node in the form of heat flux density; Steel billet-air exposed surface: convection and radiation coupled boundary, direct calculation of nodal heat flux; Steel billet bottom-pad beam contact surface: fixed contact thermal resistance boundary, and the thermal resistance value is uniformly fixed; All boundary conditions are transformed into nodal temperature algebraic equations and substituted into the global discrete equation system. Time step iteration and target time solution: Inputs: Initial temperature field, property matrix, boundary conditions, planned extraction time. ; Iterative process: from the moment of entering the slab warehouse Beginning, with Using a time step size, solve the linear equations step by step and update the temperature of all nodes; Convergence criterion: The current step is considered converged when the maximum temperature change in a single step meets the preset requirements; if it does not converge, the iteration margin is automatically reduced and the calculation is repeated until the conditions are met. Output: Arrived At time t, output the three-dimensional temperature field of each steel billet. Core temperature, surface temperature, maximum temperature difference, temperature gradient, and cooling rate.

6. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 1, characterized in that, The stacking topology optimization and total enthalpy evaluation module specifically includes: Total enthalpy formula: ； Numerical integration method: The three-dimensional trapezoidal integration method is used to integrate the volume of each grid cell occupied by the steel billet. ;in, For a single grid volume ( ); Multi-constraint deterministic topology optimization: Topology scheme generation: Enumerate all stacking methods that satisfy physical constraints, and the system generates candidate schemes according to the insulation logic; Constraints include equipment constraints, safety constraints, and process constraints; A heuristic traversal optimization method is used to calculate the enthalpy value of all feasible solutions and directly select the appropriate one. The largest possible solution is the optimal topology.

7. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 6, characterized in that, It also includes phase change threshold interception and security determination: extract Lowest local temperature of steel billet at any time ,like If the cooling rate exceeds the material's critical cooling rate, it is also considered a risky solution. Automatically reject the proposed solution and recalculate the temperature field until... And the cooling rate meets the standard; Output the billet number, target coordinates, stacking relationship, insulation method, expected enthalpy value, temperature index, and safety judgment result of the optimal topology to form a unique executable solution.

8. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 1, characterized in that, The physical scheduling instruction conversion and device coordination module specifically includes: The theoretical coordinates of thermodynamic topology are converted into absolute pulse coordinates that can be recognized by the overhead crane PLC using rigid body coordinate transformation formulas. The action commands are broken down into transfer commands, clamping commands, and positioning commands. Stacking action sequence generation: Generate an action sequence with absolute timestamps according to preset principles.

9. The panoramic intelligent collaborative management platform for steel enterprises based on industrial big data as described in claim 8, characterized in that, Also includes: The system continuously compares the crane's return to position signal, action completion signal, fault signal, and current coordinates with the optimal topology in real time. If the deviation exceeds the preset allowable range, the system determines that the execution is abnormal. After the exception handling action is completed, the new coordinates and temperature field of the billet are updated, and the initial boundary matrix is ​​written synchronously as the initial condition for the next scheduling.

Images