Roll shaft temperature field acquisition method and device, roll shaft equipment and medium

Abstract

The invention relates to the technical field of industrial thermal analysis, and discloses a roll shaft temperature field acquisition method and device, roll shaft equipment and a medium, and the roll shaft temperature field acquisition method comprises the steps: carrying out the discrete division of a roll shaft based on the spatial structure information of a target roll shaft, and determining nodes and a node spatial distribution relation; based on the nodes and the node space distribution relation, constructing a heat conduction topological structure, and identifying the nodes at the boundary positions to set corresponding heat boundary conditions; establishing a heat balance equation according to the heat conduction topological structure and a heat boundary condition; performing temperature setting on nodes forming the heat balance equation to form an initial temperature distribution state; iterative solution is carried out in the temperature distribution state, a node temperature distribution result is obtained, and then a steady-state temperature field distribution result corresponding to the target roller shaft is generated. According to the method, rapid modeling and visual output of the temperature field are achieved, rapid reuse of multiple working conditions is supported, and the calculation efficiency of engineering analysis and the flexibility of design optimization are remarkably improved.

Description

Technical Field

[0001] This application relates to the field of industrial thermal analysis technology, and in particular to a method, apparatus, roller equipment and medium for obtaining the temperature field of a roller. Background Technology

[0002] In industrial manufacturing processes such as battery electrode rolling, the roller shaft, as a key heat transfer component, directly affects the forming quality and process stability of the material. To achieve accurate analysis of the roller shaft temperature distribution, current engineering practices primarily employ finite element method (FEM) or computational fluid dynamics (CFD) software for simulation calculations. These methods typically establish a three-dimensional geometric model of the roller shaft, define the material's thermal properties (such as thermal conductivity and specific heat capacity), heat source conditions (such as heating rod power), and boundary conditions (such as ambient temperature and convective heat transfer coefficient), and then use a numerical solver to simulate the heat conduction process, thereby obtaining the temperature field distribution results. Theoretically, this method, based on Fourier's law of thermal conductivity and the principle of energy conservation, can model and solve heat conduction processes with high accuracy under complex structures and multi-physics coupling conditions, making it a widely used technique in the field of thermal analysis engineering. However, in practical engineering applications, existing simulation methods suffer from problems such as complex modeling, low computational efficiency, and reliance on high-performance computing resources, making it difficult to meet the practical needs of rapid design and multi-condition optimization. Furthermore, simulation methods require high precision and appropriateness in mesh generation; improper mesh generation can affect the accuracy of calculation results. When dealing with dynamic boundary conditions such as natural convection and heat source distribution, simulation results may also deviate due to model simplification or inaccurate parameter settings, thus limiting their applicability in practical engineering. Summary of the Invention

[0003] In view of this, the embodiments of this application provide a method, apparatus, roller equipment and medium for obtaining roller temperature field, which can effectively solve the problems of complex modeling, low calculation efficiency, strong dependence on high-performance computing resources and inaccurate boundary condition handling in existing roller temperature field calculation methods.

[0004] In a first aspect, embodiments of this application provide a method for obtaining the temperature field of a roller, comprising: The spatial structure information of the target roller shaft is obtained, and the target roller shaft is discretized based on the spatial structure information to determine the nodes used for thermal calculation and the spatial distribution relationship of the nodes. A heat conduction topology is constructed based on the nodes and their spatial distribution, and nodes at the boundary positions are identified to set corresponding thermal boundary conditions. Based on the aforementioned thermal conduction topology and thermal boundary conditions, a discretized thermal equilibrium equation is established. An initial temperature is set for each node constituting the heat balance equation to form an initial temperature distribution state. Iterative solutions are performed under the initial temperature distribution state to obtain the node temperature distribution results that satisfy the preset convergence conditions; Based on the node temperature distribution results, the steady-state temperature field distribution results corresponding to the target roller are generated.

[0005] In some embodiments, obtaining the spatial structure information of the target roller shaft, and discretizing the target roller shaft based on the spatial structure information to determine the nodes used for thermal calculation and the spatial distribution relationship of the nodes, includes: Extract the geometric dimensional parameters of the target roller shaft; The target roller shaft is divided radially and axially based on the geometric parameters to obtain a grid cell with axisymmetric characteristics; The center position of the grid cell is used as the node coordinates to determine the nodes used to participate in the thermal calculation; Based on the coordinate information of the nodes, the distribution relationship of the nodes in the radial and axial directions is calculated to obtain the spatial distribution relationship of the nodes.

[0006] In some embodiments, the step of constructing a heat conduction topology based on the nodes and their spatial distribution relationships, and identifying nodes at boundary locations to set corresponding thermal boundary conditions, includes: Based on the distribution relationship of each node in the radial and axial directions, the thermal connection relationship between the node and its adjacent nodes is determined, and a thermal conduction topology is constructed to describe the heat conduction path. Based on the heat conduction topology, a set of heat flow conduction paths is generated, and the heat flow conduction path between each node and its adjacent nodes is determined. Based on the set of heat conduction paths, identify the set of nodes located at the boundary. Based on the spatial location of the node set, set the corresponding thermal boundary conditions.

[0007] In some embodiments, establishing a discretized thermal equilibrium equation based on the thermal conduction topology and the thermal boundary conditions includes: Based on the thermal boundary conditions, determine the thermal input or thermal dissipation terms for each node; By combining the heat flow conduction path and the heat input or heat dissipation terms, the heat balance equations for the corresponding nodes are established according to Fourier's law of thermal conductivity.

[0008] In some embodiments, setting an initial temperature for each node constituting the heat balance equation to form an initial temperature distribution includes: An initial temperature value is set for each node that constitutes the heat balance equation; Substituting the initial temperature value into the heat balance equation generates the corresponding initial temperature distribution state.

[0009] In some embodiments, the iterative solution performed under the initial temperature distribution state to obtain the node temperature distribution result that satisfies the preset convergence condition includes: Based on the initial temperature distribution state, the heat balance equation is iteratively calculated to obtain the temperature distribution result of the current iteration. Calculate the temperature change based on the temperature distribution results of the current iteration and the previous iteration; Determine whether the temperature change satisfies the preset convergence condition. If not, continue the next iteration calculation based on the temperature distribution result of the current iteration until the preset convergence condition is met. If it is satisfied, use the temperature distribution result of the current iteration as the node temperature distribution result that satisfies the convergence condition.

[0010] In some embodiments, generating the steady-state temperature field distribution result corresponding to the target roller shaft based on the node temperature distribution result includes: The temperature distribution results of nodes that meet the preset convergence conditions are mapped to the spatial coordinates of the corresponding nodes to construct the temperature distribution relationship. Based on the temperature distribution relationship, the steady-state temperature field of the target roller shaft in the discrete calculation region is reconstructed, and the steady-state temperature field distribution result corresponding to the target roller shaft is output.

[0011] Secondly, embodiments of this application provide a device for acquiring the temperature field of a roller, comprising: The partitioning module is used to acquire the spatial structure information of the target roller shaft, and to discretize the target roller shaft based on the spatial structure information to determine the nodes used for thermal calculation and the spatial distribution relationship of the nodes. The condition setting module is used to construct a heat conduction topology based on the nodes and their spatial distribution relationship, and to identify nodes at the boundary positions in order to set the corresponding thermal boundary conditions. The equation construction module is used to establish thermal balance equations characterizing the conservation of nodal energy based on the thermal conduction topology and the thermal boundary conditions. The temperature setting module is used to set the initial temperature for each node that constitutes the heat balance equation, thereby forming an initial temperature distribution state. The solution module is used to perform iterative solutions under the initial temperature distribution state to obtain the node temperature distribution results that satisfy the preset convergence conditions; The result acquisition module is used to generate the steady-state temperature field distribution result corresponding to the target roller shaft based on the node temperature distribution result.

[0012] Thirdly, embodiments of this application provide a roller device, the roller device including a processor and a memory, the memory storing a computer program, and the processor executing the computer program to implement the roller temperature field acquisition method of the first aspect described above.

[0013] Fourthly, embodiments of this application provide a computer-readable storage medium, wherein when the computer program is executed on a processor, it implements the roller temperature field acquisition method described in the first aspect.

[0014] The embodiments of this application have the following beneficial effects: By acquiring the spatial structure information of the target roller shaft and discretizing the roller shaft structure based on this information, the nodes involved in the thermal calculation and their spatial distribution relationships are determined. A heat conduction topology is constructed based on the node distribution, and boundary nodes are identified to set corresponding thermal boundary conditions. A heat balance equation characterizing the energy conservation of the nodes is established based on the heat conduction topology and boundary conditions. An initial temperature distribution is set for the nodes, and iterative solutions are performed. The timing of iteration termination is determined using preset convergence conditions. Finally, the steady-state temperature field distribution of the roller shaft is generated based on the node temperature distribution results obtained from the iteration. This achieves visualized output of the temperature field and rapid reuse across multiple operating conditions, improving the efficiency and flexibility of engineering analysis. The method of this application does not require complex modeling and high-density mesh generation, has high computational efficiency, low resource consumption, and supports rapid iteration across multiple operating conditions, making it suitable for roller shaft heating design and optimization in engineering practice. Attached Figure Description

[0015] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0016] Figure 1 A flowchart of a method for obtaining the roller temperature field according to an embodiment of this application is shown; Figure 2 A schematic diagram of the roller structure in the roller temperature field acquisition method of this application is shown; Figure 3 A schematic diagram of the roller dimensional parameters in the method for obtaining the roller temperature field according to an embodiment of this application is shown; Figure 4 A schematic diagram of the discrete structure configuration of the roller shaft in the method for obtaining the roller shaft temperature field according to an embodiment of this application is shown; Figure 5A schematic diagram of the thickness of the equally divided structure is shown in the method for obtaining the roller temperature field according to an embodiment of this application; Figure 6 The diagram shows a typical boundary region of the roller structure in an axisymmetric section in the roller temperature field acquisition method of the present application embodiment; Figure 7 A schematic diagram of the nodal structure of the roller structure in the method for obtaining the roller temperature field according to an embodiment of this application is shown; Figure 8 A schematic diagram of the boundary nodes of the roller structure in the roller temperature field acquisition method of this application is shown; Figure 9 A schematic diagram of nodes A to F in the roller temperature field acquisition method of an embodiment of this application is shown; Figure 10 This illustration shows an example of calculating the roller temperature field using Excel in the method for obtaining the roller temperature field according to an embodiment of this application. Figure 11 A schematic diagram of temperature result output in the roller temperature field acquisition method of this application embodiment is shown; Figure 12 A schematic diagram of a two-dimensional temperature cloud map is shown in the method for obtaining the roller temperature field according to an embodiment of this application; Figure 13 A schematic diagram of steady-state thermal field simulation results is shown in the method for obtaining the roller temperature field according to an embodiment of this application; Figure 14 A comparative schematic diagram of the distribution of roller surface temperature along the length direction in the roller temperature field acquisition method of the embodiments of this application is shown; Figure 15 A schematic diagram of a method for obtaining the temperature field of a roller in an embodiment of this application is shown. Detailed Implementation

[0017] The technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments.

[0018] The components of the embodiments of this application described and illustrated in the accompanying drawings can be arranged and designed in a variety of different configurations. Therefore, the following detailed description of the embodiments of this application provided in the drawings is not intended to limit the scope of the claimed application, but merely to illustrate selected embodiments of the application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application.

[0019] In the following text, the terms "comprising," "having," and their cognates, which may be used in various embodiments of this application, are intended only to indicate a particular feature, number, step, operation, element, component, or combination thereof, and should not be construed as primarily excluding the presence of one or more other features, numbers, steps, operations, elements, components, or combinations thereof, or adding the possibility of one or more combinations thereof. Furthermore, the terms "first," "second," "third," etc., are used only for distinguishing descriptions and should not be construed as indicating or implying relative importance.

[0020] Unless otherwise specified, 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 the various embodiments of this application pertain. Terms (such as those defined in commonly used dictionaries) shall be interpreted as having the same meaning as in their contextual meaning in the relevant technical field and shall not be construed as having an idealized or overly formal meaning, unless clearly defined in the various embodiments of this application.

[0021] The following detailed description of some embodiments of this application is provided in conjunction with the accompanying drawings. Unless otherwise specified, the following embodiments and features can be combined with each other.

[0022] Considering that existing methods for calculating roller temperature fields suffer from problems such as complex modeling, low computational efficiency, strong dependence on high-performance computing resources, and limited accuracy when dealing with dynamic boundary conditions, this application proposes a method for obtaining roller temperature fields. By discretizing the roller structure based on its axisymmetric structural characteristics, constructing the heat conduction relationship between nodes using Fourier's law of thermal conductivity, and establishing a heat balance equation in the form of energy conservation, an iterative solution method is used to efficiently calculate the temperature field. This method achieves fast, stable, and high-precision solutions for roller temperature fields without the need for complex 3D modeling and mesh generation.

[0023] The following describes the method for obtaining the temperature field of the roller shaft using some specific embodiments.

[0024] Figure 1 A flowchart illustrating a method for obtaining the roller temperature field according to an embodiment of this application is shown. Exemplarily, the method for obtaining the roller temperature field includes the following steps: Step S100: Obtain the spatial structure information of the target roller shaft, and discretize the target roller shaft based on the spatial structure information to determine the nodes used for thermal calculation and the spatial distribution relationship of the nodes.

[0025] The spatial structure information refers to the geometric parameters and axisymmetric physical structure of the target roller shaft, including inner and outer radii, length, number of segments, and symmetry conditions. This information supports mesh generation and node modeling operations in subsequent thermal field calculations. Demonstratively, an axisymmetric simplified model is used, uniformly dividing the target roller shaft into n parts along its circumference. Calculating the thermal field distribution of any one part is sufficient to represent the overall temperature distribution. Each part is discretized using multi-layer discretization in both the radial and axial directions to form a spatial mesh structure. The node coordinates are determined by calculating the center positions of each mesh, thus constructing a complete discretized node system and the spatial distribution relationship of the nodes.

[0026] In an optional embodiment, step S100 includes the following sub-steps: S101, extract the geometric parameters of the target roller.

[0027] Among them, geometric dimension parameters refer to the variables used to characterize the structural dimensions of the roller shaft, including the inner radius r and the center radius of the roller shaft. Outer radius and axial length distribution parameters , This is used to define the boundaries of the thermal computation region and the dimensions of the discrete structural elements. For example, as shown... Figure 2 As shown, the target roller shaft is a symmetrical hollow structure with a central hole for placing an electric heating rod. The sidewalls are in contact with air, creating a typical heating-dissipation environment. Further, as... Figure 3 As shown, based on the requirements of heat conduction modeling, key variables related to thermal boundaries and heat conduction paths in the structural dimensions are extracted, including radius r and median radius r. Medium radius ,length ,length Five key dimensions are used for subsequent steps to perform segmentation and temperature distribution calculations.

[0028] S102, the target roller shaft is divided radially and axially based on geometric dimension parameters to obtain a mesh element with axisymmetric characteristics.

[0029] The radial and axial divisions refer to dividing the roller shaft into several discrete regions along its radius and length directions, based on the axisymmetric structural assumption, to form structural units used to establish thermal equilibrium relationships at different locations. Exemplarily, the radial division employs a two-segment structure, extending from the inner radius r to the middle radius r. , and from the mid-radius To outer radius It consists of two segments, and the corresponding number of divisions is: , ; The lengths of the two ends of the roller shaft are respectively measured in the axial direction Intermediate length Discretized , Therefore, the expression for the inter-element thickness can be obtained: ; To further clarify the grid cell arrangement relationship of the above-mentioned partitioned structure, Figure 4 The discrete structural configuration of the roller after being divided in the radial and axial directions is shown. The grid regions correspond to different structural segments, and a corresponding number of equally divided units are set in each region. This structural division ensures the spatial continuity of the heat conduction path of each part of the model and matches the actual thermal behavior of the structure.

[0030] S103 uses the center position of the grid cell as the node coordinate to determine the nodes used for thermal calculation.

[0031] In this context, a node is a discrete point in a discrete structure used to establish the expression of heat conduction and energy conservation. It represents the temperature state of each tiny grid cell, and its center coordinates determine its heat transfer relationship with adjacent nodes in the model. Demonstratively, in each pre-divided grid cell, its center is taken as the node location, and thermally conductive connections are established between nodes. Once the node locations are determined, they can be used as calculation points.

[0032] S104. Based on the coordinate information of the nodes, calculate the distribution relationship of the nodes in the radial and axial directions to obtain the spatial distribution relationship of the nodes.

[0033] The spatial distribution of nodes refers to the relative position of any node in the entire heat conduction topology and its physical connection with adjacent nodes (top, bottom, left, right, inside, outside), forming the basis for constructing the heat balance equation. Exemplarily, in the axial direction, each equal segment is... (For example The radial direction is divided according to unequal thickness. Note that, see [link to relevant documentation]. Figure 4 The thickness of the above equally divided structure is not consistent in the z-direction (circumferential direction), that is, the thickness of each column of the discretized grid along the y-axis is equal, but the thickness of each row along the x-axis is different. The formula for calculating the thickness of each discretized grid is as follows: ; This leads to the expression for the radial thickness difference between the two segments: ; ; Where a and b represent the thickness difference between adjacent nodes within the radius intervals r-R1 and R1-R2, respectively. By progressively superimposing a or b, the specific position coordinates of each node in the radial direction can be obtained, such as... Figure 5As shown, once the spatial distribution relationship between nodes is determined, the heat conduction direction of each node and the heat conduction path it participates in can be identified, providing a clear physical basis for the next step of constructing the heat conduction topology.

[0034] Step S200: Construct a heat conduction topology based on nodes and their spatial distribution, and identify nodes at the boundary to set corresponding thermal boundary conditions.

[0035] In this context, the heat conduction topology refers to the topological network that describes the heat conduction path by connecting adjacent grid nodes; boundary nodes are nodes located at the spatial boundaries of the roller shaft, requiring specific boundary conditions to be applied during heat conduction. Exemplarily, the radial and axial coordinate information of each node is obtained, and the connection relationships between nodes are determined to construct a heat connection topology. Based on this, nodes located at the structural edges in the topology are further identified, classified as boundary nodes, and assigned corresponding convective heat transfer or heat source boundary conditions according to the region where the node is located.

[0036] In an optional embodiment, step S200 includes the following sub-steps: S201, based on the radial and axial distribution relationship of each node, determines the thermal connection relationship between the node and its adjacent nodes, and constructs a thermal conduction topology to describe the heat conduction path.

[0037] Thermal connectivity refers to the topological relationship between a node and its adjacent nodes that have a heat transfer path in space. For example, based on the known node positions and distribution sequence in the grid, adjacent nodes in the radial and axial directions (up, down, left, and right) are found for each node, and connecting edges are established between them. The resulting node-edge structure forms a heat conduction topology graph, which is used for subsequent heat flow path extraction and heat balance calculation.

[0038] S202 generates a set of heat flow conduction paths based on the heat conduction topology, and determines the heat flow conduction path between each node and its adjacent nodes.

[0039] Here, a heat conduction path refers to a heat conduction path represented by the connecting edges in the topology, and includes parameters such as thermal conductivity, connection width, and distance; the set of heat conduction paths refers to a complete set of path information covering all thermal connections between nodes. Exemplarily, all connecting edges in the heat conduction topology are traversed, and the expression for the heat flux on each path is calculated based on the node number, geometric parameters (such as cross-sectional area and distance), and thermal conductivity k connected to each edge.

[0040] S203 identifies the set of nodes at the boundary position based on the set of heat flow conduction paths.

[0041] The set of nodes at boundary locations refers to the peripheral nodes that have thermal connections with nodes in only some directions along the heat flow path. These typically correspond to geometric boundaries, heated surfaces, or convective heat transfer surfaces. For example, all nodes in the heat flow path set are scanned. If a node has no corresponding thermal connection edge in a specific direction (such as the outer diameter or end), it is determined to be a boundary node, and its number is recorded in the boundary node set for subsequent boundary condition configuration.

[0042] S204, set the corresponding thermal boundary conditions according to the spatial location of the node set.

[0043] Thermal boundary conditions refer to the heat input or heat dissipation terms applied to the boundary nodes, including electric heating flux q or natural convection heat transfer terms. Setting boundary conditions refers to determining the physical region where a node is located based on its spatial position and assigning a corresponding boundary heat flux form. For example, positional analysis is performed on each node in the boundary node set: if it is located on the inner wall of the roller shaft, it is determined to be within the area of ​​the electric heating rod, and a fixed heating intensity q is specified as the heat input; if it is located on the outer surface of the roller shaft and in contact with air, the heat transfer coefficient h is set based on convection and the ambient temperature. And construct the corresponding heat dissipation terms; if it is located on the axisymmetric plane, it is assumed that there is no heat transport behavior in the region, and an adiabatic boundary condition is applied, that is, the heat flux is zero.

[0044] like Figure 6 As shown, the roller structure contains three typical boundary regions in an axisymmetric section: an electric heating rod is installed in the central cavity, forming an "electric heating surface"; two "symmetric planes" exist in the axial direction, where heat does not undergo lateral transfer and are designated as heat flux-free surfaces; the outer region is exposed to the air environment, forming an "external heat exchange surface," where nodes dissipate heat outward through natural convection, simulating the heat exchange process between the roller and the environment. That is, based on the position of the boundary nodes in these typical physical surfaces, corresponding thermal boundary conditions are applied to facilitate the subsequent construction of the thermal balance equations.

[0045] Step S300: Based on the heat conduction topology and thermal boundary conditions, establish the discretized heat balance equation.

[0046] The discretized heat balance equation refers to a set of energy conservation equations constructed using nodes as basic units, based on the heat conduction relationship between the node and its neighboring nodes, and the boundary heat input or dissipation conditions. Exemplarily, the heat conduction topology and boundary condition configuration results output from the above steps are received. For each grid node, its connection path with neighboring nodes and the corresponding heat conduction flux expression are extracted, and whether it is a boundary node is determined whether it includes a heating power term or a convective heat transfer term. Based on this, heat balance expressions are written for all nodes individually, and the discretized heat balance equations for the entire domain are constructed by summarizing them.

[0047] In an optional embodiment, step S300 includes the following sub-steps: S301, based on the thermal boundary conditions, determine the thermal input or thermal dissipation terms for each node.

[0048] Here, the heat input term refers to the heat flux per unit area applied to the heated surface node, often expressed as q; the heat dissipation term refers to the natural convection heat transfer between the boundary node and the environment, usually expressed as... ,in The convective heat transfer coefficient is... For ambient temperature, The current temperature of the node. As an example, the thermal conductivity of the roller is set to [value missing]. The heat transfer coefficient with air is The ambient temperature is The output power of the central heating rod is , The radius of the inner wall of the central axis. , For different lengths of the roller shaft, the heating power per unit area of ​​the inner wall is calculated as follows: ; In addition, such as Figure 7 The discretized geometric mesh of the roller shown extracts the geometric center coordinates of each mesh cell, which serves as the heat conduction node represented by that cell. Each node has unique horizontal (x) and vertical (y) coordinates in a two-dimensional plane, forming a data reference point set for thermal topology construction and thermal balance calculation. This node set includes both internal nodes and boundary nodes located in the boundary region; that is, boundary heat input and heat conduction paths are applied based on this node set. In other words, Figure 7 This illustrates the heat conduction path between any node (m,n) and its four neighboring nodes. Taking the left neighboring node (m−1,n) as an example, the expression for the heat flow transferred from it to (m,n) is: ; Where: k is the thermal conductivity; The horizontal length of the node; The longitudinal distance between node centers; The thickness of the heat transfer surface in this direction; Let (m, n) represent the temperatures of the left adjacent node (m−1, n) and the current node (m, n), respectively. Ultimately, all boundary heat inputs and heat transfer terms will work together to construct the heat balance expression for each node.

[0049] Furthermore, in addition to the heat flow corresponding to the left node In addition, a similar method can be used to derive the heat conduction on the upper, right, and lower sides of a node, denoted as […]. , , Applying the principle of energy conservation to the element (m,n), its heat balance equation can be obtained as follows: ; In one alternative implementation, the type of heat input to the node needs to be defined before establishing the heat balance equation. Figure 8 This paper illustrates six typical boundary node types that may exist in a roller mesh structure, including different heat input, heat dissipation, and transition connection types, denoted as A to F. Type A nodes: inner surface nodes of the heating cavity, receiving heat flow input q; Type B nodes: outer surface of the roller, exchanging heat with the outside air; Type C nodes: structural turning points; Types D, E, and F nodes: intermediate transition areas or corner points, also exhibiting the combined effects of heat flow input and multi-directional conduction. It can be understood that... Figure 8 The positions marked at each point are the boundary nodes for which a separate heat balance expression will be established in the subsequent derivation.

[0050] S302 combines the heat flow conduction path with the heat input or heat dissipation terms, and establishes the energy conservation expression for the corresponding node based on Fourier's law of thermal conductivity.

[0051] The energy conservation expression refers to the heat input and output balance formula established at each node, used to describe the state where heat does not accumulate at that node under steady state. Fourier's law of thermal conductivity is a fundamental law characterizing heat transfer behavior and is applicable to heat flow calculations between nodes. Exemplarily, based on the set of heat flow conduction paths, the connection relationships between each node and its adjacent nodes, along with the corresponding spatial distance, connection width, and thermal conductivity parameters, are extracted. Combined with boundary heat input or heat dissipation terms, all heat flow input and output terms are summed to establish the node's heat balance formula. For internal nodes, the heat transfer term is entirely calculated from the heat conduction paths between their adjacent nodes; for boundary nodes, it includes both heating power from the heat source and convective heat transfer with the environment. Finally, an energy conservation expression is established for each node, and all expressions are recorded in a sparse structure to form the overall discretized heat balance equation. Figure 9 As shown, the following describes the process of constructing the boundary thermal equilibrium expressions from node A to node F (wherein, the heat source power density is...). Thermal conductivity k, convective heat transfer coefficient h): Node A: Electrically heated boundary + natural convection; where node A is located at the junction of the heated surface and the natural heat transfer boundary, and contains both heat source input and convective heat dissipation terms. Exemplarily, boundary parameters are extracted based on the schematic boundary structure of node A shown in the figure, where: The horizontal length between node centers; The longitudinal distance between node centers; These represent the temperatures of the upper adjacent node (m, n+1), the right adjacent node (m+1, n), and the current node (m, n), respectively. w1 and w2 represent the thicknesses of the heat transfer surfaces at nodes (m,n) and (m,n+1), respectively. The following heat balance expression is constructed: ; Its explicit solution can be transformed into: ; Node B: Electrically heated boundary (lower + right); Node B is located in the middle of the electrically heated surface, with heat conduction paths on all sides, and the heat input mainly comes from the electrical heating power term. Exemplarily, boundary parameters are extracted based on the schematic boundary structure of node B shown in the figure, where: The horizontal length between node centers; These are the vertical distances between the left and right centers of nodes (m,n), respectively. These represent the temperatures of the upper adjacent node (m, n+1), the right adjacent node (m+1, n), the left adjacent node (m+1, n), and the current node (m, n), respectively. w1 and w2 represent the thicknesses of the heat transfer surfaces at nodes (m,n) and (m,n+1), respectively. The heat balance equation is constructed based on the thermal conductivity and geometric parameters in each direction along the heat flow path: ; Its solution is expressed as: ; Node C: Natural heat transfer boundary + bidirectional heat conduction; Node C is located above the natural heat transfer surface, exchanging heat through convection. The boundary term is mainly... As an example, boundary parameters are extracted based on the schematic boundary structure of node C shown in the figure, where: These represent the horizontal lengths between the lower and upper centers of nodes (m,n), respectively. The longitudinal distance between node centers; These represent the temperatures of the upper adjacent node (m, n+1), the right adjacent node (m+1, n), the lower adjacent node (m, n-1), and the current node (m, n), respectively. w1 and w2 are the thicknesses of the heat transfer surfaces at nodes (m,n), (m,n+1), and (m,n-1), respectively. Considering the effects of top-to-bottom heat conduction and natural convection, this node is constructed as follows: ; The expression for the solution is: ; Node D: Double natural heat transfer boundary; wherein node D is in contact with two natural convection surfaces simultaneously. Exemplarily, boundary parameters are extracted based on the schematic boundary structure of node D shown in the figure, wherein: The horizontal length between node centers; The longitudinal distance between node centers; These represent the temperatures of the right adjacent node (m+1,n), the lower adjacent node (m,n-1), and the current node (m,n), respectively. w2 and w2 are the thicknesses of the heat transfer surfaces at nodes (m,n) and (m,n-1), respectively. A dual heat transfer term is established: ; The converted expression is: ; Node E: Natural heat transfer + multi-directional heat conduction boundary; the boundary of node E is relatively complex, containing three-way heat conduction and natural convection terms. Exemplarily, boundary parameters are extracted based on the schematic boundary structure of node E shown in the figure, where: These represent the horizontal lengths between the lower and upper centers of nodes (m,n), respectively. These are the vertical distances between the left and right centers of nodes (m,n), respectively. These represent the temperatures of the upper adjacent node (m, n+1), the right adjacent node (m+1, n), the lower adjacent node (m, n-1), the left adjacent node (m-1, n), and the current node (m, n), respectively. w1 and w2 are the thicknesses of the heat transfer surfaces at nodes (m,n), (m,n+1), and (m,n-1), respectively. The following heat balance expression is constructed: ; The expression for the solution is: ; Node F: A four-way heat-conducting closed boundary; Node F is located at the inner closed boundary, with no heat transfer, and only the heat conduction path is considered. Exemplarily, the boundary parameters are extracted based on the schematic boundary structure of node F shown in the figure, where: These represent the horizontal lengths between the lower and upper centers of nodes (m,n), respectively. These are the vertical distances between the left and right centers of nodes (m,n), respectively. These represent the temperatures of the upper adjacent node (m, n+1), the right adjacent node (m+1, n), the lower adjacent node (m, n-1), the left adjacent node (m-1, n), and the current node (m, n), respectively. Let w1 and w2 be the thicknesses of the heat transfer surfaces at nodes (m,n), (m,n+1), and (m,n-1), respectively. Construct the following heat balance equation: ; Its conversion form is: ; Step S400: Set the initial temperature for each node that constitutes the heat balance equation to form the initial temperature distribution state.

[0052] The initial temperature setting refers to assigning a starting temperature value to each computational node before iteratively solving the heat balance equations. The initial temperature distribution refers to the two-dimensional temperature field structure formed by the initial temperatures of all nodes, used for iterative calculations. For example, first, all nodes participating in the solution of the heat balance equations are identified, and based on physical field assumptions or empirical temperature values, initial temperature values ​​are assigned to each node uniformly or regionally. For instance, the initial temperature of all nodes is set to the ambient temperature. Alternatively, different initial values ​​can be set for different regions based on boundary characteristics. After initialization, these temperature values ​​are substituted into the thermal balance expression of the corresponding node to form the input temperature field state for the first iteration.

[0053] In an optional embodiment, step S400 includes the following sub-steps: S401 sets the initial temperature value for the nodes that constitute the heat balance equation.

[0054] Setting the initial temperature value refers to assigning an initial temperature to each discrete node before the iteration to facilitate the initial energy conservation calculation. This value can be a globally uniform constant or different initial values ​​can be set based on spatial location. For example, in scenarios with symmetrical structures and a relatively uniform expected temperature distribution, the initial temperature of all nodes can be set to [value missing]. If a temperature difference is predicted near the boundary, initial high or low temperature conditions can be specified for each boundary node. After setting, a data structure corresponding to the node index and its initial value will be generated, providing a starting benchmark for the subsequent temperature iteration process.

[0055] S402, based on the initial temperature value, is substituted into the heat balance equation to generate the corresponding initial temperature distribution state.

[0056] The initial temperature distribution state refers to the numerical state of the heat flow between nodes and the heat balance equations before they are iteratively solved under the initial temperature conditions. This state is used to trigger the subsequent solver startup. For example, the initial temperatures of all nodes are substituted into their corresponding heat balance expressions to form the reference state used to calculate the temperature in the first iteration. This initial state includes the current temperature value of each node, its connection relationships with adjacent nodes, and boundary heat sources or convection conditions. As the starting point for the thermal iteration solution, it will be used to perform temperature updates in the next step.

[0057] Step S500: Iteratively solve the problem under the initial temperature distribution state to obtain the node temperature distribution result that satisfies the preset convergence condition.

[0058] Iterative solution refers to the process of approximating the steady-state temperature distribution by continuously updating the heat balance equations under a given initial temperature condition. Preset convergence criteria refer to the temperature change used to determine whether the iteration is complete; absolute or relative convergence criteria are commonly used. For example, in the default case, the convergence criterion is whether the temperature difference between two iterations for all nodes is less than a set absolute threshold (e.g., ...). This criterion can be used as a convergence criterion. However, when there are regions with high heating power or low thermal conductivity, and steep local temperature gradients, this criterion may not be sufficient to reflect the convergence status. Therefore, a relative convergence criterion can also be used, namely: ; in, For example, the relative error threshold. or The relative criterion can avoid the risk of division by zero in the low-temperature region and has a relatively fair tolerance for errors in different gradient regions, thus improving the overall convergence efficiency.

[0059] In an optional embodiment, step S500 includes the following sub-steps: S501, based on the initial temperature distribution state, iteratively calculates the heat balance equation to obtain the temperature distribution result of the current iteration.

[0060] In this process, iterative calculation refers to updating the heat balance equation node by node based on the current node's temperature field, generating a new round of temperature distribution results. The current iteration round refers to the iteration state numbered k-th round. For example, the temperature result of round k-1 is used as input, and substituted into the node heat balance equations one by one. Combining the connection relationships between adjacent nodes, thermal conductivity, and boundary heat source terms, the new temperature of round k is obtained. After all nodes have been updated, a new round of temperature distribution output is generated.

[0061] S502, calculate the temperature change based on the temperature distribution results of the current iteration and the previous iteration.

[0062] Among them, temperature change is an indicator used to measure the degree of temperature difference between two iterations, and can be expressed as absolute change. or relative change To account for the temperature range differences in different regions, the above calculation is performed for each node as an example, obtaining the current change value of that node, and the maximum change among all nodes is counted as the input for convergence judgment.

[0063] S503, determine whether the temperature change meets the preset convergence condition. If not, continue the next iteration calculation based on the temperature distribution result of the current iteration until the preset convergence condition is met. If it is met, use the temperature distribution result of the current iteration as the node temperature distribution result that meets the convergence condition.

[0064] The convergence condition determination refers to the process of terminating or continuing the iteration by comparing whether the maximum temperature change is lower than a threshold. Exemplarily, an absolute threshold is preset. or relative threshold The maximum change value obtained in step S502 is compared with this value for judgment. If no criterion is met, the iteration continues to return to step S501; if one of them is met, it is determined that the current temperature field has converged and stabilized.

[0065] Step S600: Based on the node temperature distribution results, generate the steady-state temperature field distribution results corresponding to the target roller.

[0066] The node temperature distribution result refers to the set of steady-state temperature values ​​obtained by iteratively solving the thermal balance equations and converging at each discrete node. The steady-state temperature field distribution result refers to the temperature distribution in two-dimensional or three-dimensional space constructed by mapping the node temperature values ​​to their spatial coordinates in the target roller. This temperature field is used to reflect the temperature distribution in various regions of the target roller under steady-state conditions.

[0067] In an optional embodiment, step S600 includes the following sub-steps: S601 maps the node temperature distribution results that meet the preset convergence conditions to the spatial coordinates of the corresponding nodes to construct the temperature distribution relationship.

[0068] The temperature distribution relationship refers to the correspondence between the node temperature value and its spatial position in the axial and radial directions, reflecting the spatial trend of temperature variation. Exemplarily, based on the node coordinate information and corresponding steady-state temperature established in the aforementioned thermal topology modeling process, the triplet set of all nodes is transformed into a spatial temperature distribution relationship. Through this mapping, the temperature variation trend from the roller center to the outer edge and from the center to the end can be intuitively displayed.

[0069] S602 reconstructs the steady-state temperature field of the target roller in the discrete calculation region based on the temperature distribution relationship, and outputs the steady-state temperature field distribution result corresponding to the target roller.

[0070] Reconstructing the steady-state temperature field refers to filling the unsampled spatial points within the entire computational region with temperature data from discrete nodes through interpolation or fitting, generating a continuous temperature distribution map. The steady-state temperature field distribution result refers to saving or displaying the temperature field in the form of an image, matrix, or numerical grid. For example, based on the point set generated by S601, a bilinear interpolation algorithm or a discrete point cloud fitting algorithm is used to complete the temperature values ​​in the intermediate region, ultimately obtaining a two-dimensional temperature distribution map or a three-dimensional temperature field map covering the entire computational region for visualization and post-processing analysis.

[0071] Based on the above-mentioned steady-state temperature field solution process, this application further constructs a complete engineering calculation example, and uses Excel to complete the iterative solution of the steady-state temperature field and the input processing of boundary conditions. Figure 10 The diagram shows the key parameter settings and roller geometry used in the simulation calculation, including: total roller length 1000mm, end diameter 600mm (corresponding to radius 0.3m), middle diameter 1100mm (corresponding to radius 0.55m), middle electric heating through-hole diameter 100mm, end and middle lengths 250mm and 500mm respectively, and thermal conductivity of [missing information]. The ambient temperature is 25℃, and the heat transfer coefficient is... The internal cavity surface temperature is set to The step lengths, tangent parameters, and ratio coefficients for each discrete step are set uniformly in Excel. The figure also illustrates the relationships between geometric parameters such as R1, R2, L1, and L2.

[0072] After setting the parameters, construct the energy conservation equations for the two-dimensional heat conduction problem in Excel, and perform initial temperature setting and numerical iteration solution according to the set step size and boundary conditions. Figure 11 The typical output results of this numerical solution process are shown: Figure 11 The image shows the main interface for temperature distribution, which visualizes the overall trend of temperature field changes. Figure 11 The upper right corner shows the temperature distribution curve of the roller surface along the longitudinal direction, which reflects the gradient distribution of heat from the middle to both ends. Figure 11 The bottom right corner shows the temperature values ​​of the key sections at the right, middle, and left ends, along with the temperature differences and thermal power data between them. To further verify the solution accuracy, Figure 12 The results from the Excel calculation are displayed as a two-dimensional temperature contour plot, reflecting the gradient structure of temperature change from bottom to top. Figure 13 To obtain steady-state thermal field simulation results under the same boundary conditions using finite element simulation software (such as ANSYS), it can be observed that... Figure 12 Similar heat distribution patterns were observed. Finally, to verify the accuracy of the numerical method results, a quantitative comparison was made between the theoretical calculations and simulation results. Figure 14The diagram shows a comparison of the temperature distribution along the length of the roller surface, demonstrating a high degree of consistency between the two methods in terms of temperature change trends and numerical values, thus confirming the effectiveness and engineering applicability of the method presented in this application.

[0073] Figure 15 A schematic diagram of a roller temperature field acquisition device according to an embodiment of this application is shown. Exemplarily, the roller temperature field acquisition device 100 includes: The partitioning module 110 is used to acquire the spatial structure information of the target roller shaft, and to discretize the target roller shaft based on the spatial structure information to determine the nodes used for thermal calculation and the spatial distribution relationship of the nodes. The condition setting module 120 is used to construct a heat conduction topology based on the nodes and their spatial distribution relationship, and to identify nodes at the boundary positions in order to set corresponding thermal boundary conditions. The equation construction module 130 is used to establish a thermal balance equation characterizing the conservation of nodal energy based on the thermal conduction topology and the thermal boundary conditions. Temperature setting module 140 is used to set the initial temperature of the nodes constituting the heat balance equation to form an initial temperature distribution state; The solver module 150 is used to perform iterative solutions under the initial temperature distribution state to obtain the node temperature distribution results that satisfy the preset convergence conditions. The result acquisition module 160 is used to generate the steady-state temperature field distribution result corresponding to the target roller shaft based on the node temperature distribution result.

[0074] It is understood that the apparatus of this embodiment corresponds to the method of the above embodiments, and the options in the above embodiments are also applicable to this embodiment, so they will not be described again here.

[0075] This application also provides a roller device, which, by way of example, includes a processor and a memory, wherein the memory stores a computer program, and the processor, by running the computer program, causes the roller device to perform the functions of the various modules in the above-described method or apparatus.

[0076] The processor can be an integrated circuit chip with signal processing capabilities. The processor can be a general-purpose processor, including at least one of a Central Processing Unit (CPU), Graphics Processing Unit (GPU), Network Processor (NP), Digital Signal Processor (DSP), Application-Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. The general-purpose processor can be a microprocessor or any conventional processor, capable of implementing or executing the methods, steps, and logic block diagrams disclosed in the embodiments of this application.

[0077] The memory can be, but is not limited to, Random Access Memory (RAM), Read Only Memory (ROM), Programmable Read-Only Memory (PROM), Erasable Programmable Read-Only Memory (EPROM), Electrically Erasable Programmable Read-Only Memory (EEPROM), etc. The memory is used to store computer programs, and the processor can execute the computer programs accordingly after receiving execution instructions.

[0078] This application also provides a computer-readable storage medium for storing the computer program used in the aforementioned roller device. For example, the computer-readable storage medium may include, but is not limited to, various media capable of storing program code, such as a USB flash drive, a portable hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0079] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can also be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the flowcharts and block diagrams in the accompanying drawings show the architecture, functionality, and operation of possible implementations of apparatus, methods, and computer program products according to various embodiments of this application. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that, in alternative implementations, the functions marked in the blocks may occur in a different order than those marked in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagram and / or flowchart, and combinations of blocks in the block diagram and / or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or action, or using a combination of dedicated hardware and computer instructions.

[0080] In addition, the functional modules or units in the various embodiments of this application can be integrated together to form an independent part, or each module can exist independently, or two or more modules can be integrated to form an independent part.

[0081] If the aforementioned functions are implemented as software functional modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a smartphone, personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application.

[0082] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any changes 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.

Claims

1. A method for obtaining a temperature field of a roll, characterized in that, The method comprises: acquiring spatial structure information of a target roller shaft, and discretely dividing the target roller shaft based on the spatial structure information to determine nodes participating in heat calculation and node spatial distribution relationship; constructing a heat conduction topology structure based on the nodes and the node spatial distribution relationship, and identifying nodes at boundary positions to set corresponding heat boundary conditions; establishing a discretized heat balance equation according to the heat conduction topology structure and the heat boundary conditions; setting initial temperatures for each node constituting the heat balance equation to form an initial temperature distribution state; iteratively solving under the initial temperature distribution state to obtain a node temperature distribution result satisfying a preset convergence condition; generating a steady-state temperature field distribution result corresponding to the target roller shaft based on the node temperature distribution result.

2. The method of claim 1, wherein The acquiring spatial structure information of a target roller shaft, and discretely dividing the target roller shaft based on the spatial structure information to determine nodes participating in heat calculation and node spatial distribution relationship comprises: extracting geometric size parameters of the target roller shaft; dividing the target roller shaft in the radial and axial directions based on the geometric size parameters to obtain grid cells with axial symmetry characteristics; determining nodes participating in the heat calculation by taking the center positions of the grid cells as node coordinates; calculating the distribution relationship of the nodes in the radial and axial directions according to the coordinate information of the nodes to obtain the node spatial distribution relationship.

3. The method of claim 2, wherein The constructing a heat conduction topology structure based on the nodes and the node spatial distribution relationship, and identifying nodes at boundary positions to set corresponding heat boundary conditions comprises: determining the heat connection relationship between the nodes and adjacent nodes based on the distribution relationship of each node in the radial and axial directions to construct a heat conduction topology structure for describing heat conduction paths; generating a heat flow conduction path set based on the heat conduction topology structure to determine the heat flow conduction path between each node and the adjacent nodes; identifying a node set at boundary positions based on the heat flow conduction path set; setting corresponding heat boundary conditions according to the spatial positions of the node set.

4. The method of claim 3, wherein The establishing a discretized heat balance equation according to the heat conduction topology structure and the heat boundary conditions comprises: determining a heat input term or a heat dissipation term of each node according to the heat boundary conditions; combining the heat flow conduction path and the heat input term or the heat dissipation term to establish a heat balance equation of the corresponding node according to the Fourier heat conduction law.

5. The method of claim 1, wherein The setting initial temperatures for each node constituting the heat balance equation to form an initial temperature distribution state comprises: setting an initial temperature value for each node constituting the heat balance equation; generating a corresponding initial temperature distribution state by substituting the initial temperature value into the heat balance equation.

6. The method of claim 1, wherein The iteratively solving under the initial temperature distribution state to obtain a node temperature distribution result satisfying a preset convergence condition comprises: iteratively calculating the heat balance equation based on the initial temperature distribution state to obtain a temperature distribution result of the current iteration round; According to the temperature distribution result of the current iteration round and the last iteration round, a temperature change amount is calculated; It is judged whether the temperature change amount meets a preset convergence condition, if not, the next iteration calculation is continued based on the temperature distribution result of the current iteration round until the preset convergence condition is reached; if yes, the temperature distribution result of the current iteration round is taken as the node temperature distribution result meeting the convergence condition.

7. The method of claim 6, wherein The node temperature distribution result is used to generate a steady-state temperature field distribution result corresponding to the target roller shaft, including: The node temperature distribution result meeting the preset convergence condition is mapped to the spatial coordinates of the corresponding node to construct a temperature distribution relationship; According to the temperature distribution relationship, the steady-state temperature field of the target roller shaft in the discrete calculation region is reconstructed, and the steady-state temperature field distribution result corresponding to the target roller shaft is output.

8. A device for acquiring the temperature field of a roller, characterized in that, Including: The division module is used to obtain the spatial structure information of the target roller shaft, and discretely divide the target roller shaft based on the spatial structure information to determine the nodes participating in heat calculation and the node spatial distribution relationship; The condition setting module is used to construct a heat conduction topological structure based on the node and node spatial distribution relationship, and identify the nodes at the boundary position to set the corresponding heat boundary condition; The equation set construction module is used to establish a heat balance equation representing the energy conservation of the nodes according to the heat conduction topological structure and the heat boundary condition; The temperature setting module is used to set the initial temperature of each node constituting the heat balance equation to form an initial temperature distribution state; The solving module is used to perform iterative solving under the initial temperature distribution state to obtain the node temperature distribution result meeting the preset convergence condition; The result acquisition module is used to generate the steady-state temperature field distribution result corresponding to the target roller shaft based on the node temperature distribution result.

9. A roller arrangement, characterized by The roller shaft device includes a processor and a memory, the memory stores a computer program, and the processor is used to execute the computer program to implement the roller shaft temperature field acquisition method of any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, It stores a computer program, which implements the roller shaft temperature field acquisition method according to any one of claims 1-7 when executed on a processor.

Images