Concrete water pipe cooling temperature field calculation method, system, equipment and medium

Abstract

The invention provides a concrete water pipe cooling temperature field calculation method, system, equipment and medium, and belongs to the field of concrete construction.The concrete water pipe cooling temperature field calculation method comprises the steps that a concrete three-dimensional model containing a water pipe is constructed, and independent grid division is conducted on concrete and the water pipe in the concrete three-dimensional model; grid nodes of the divided concrete units and the water pipe units are in a non-contact relationship; establishing a geometric space mapping relation between the divided concrete units and the water pipe units, and determining geometric relation parameters of the concrete units and the water pipe units according to the mapping relation; calculating the concrete temperature and the water pipe temperature according to the geometrical relationship parameters; and updating the concrete temperature and the water pipe temperature according to the concrete pouring moment to obtain a concrete temperature field and a water pipe temperature field. According to the method, the interdependence of concrete and water pipe grid division is eliminated, heat exchange and temperature calculation of the concrete unit and the water pipe unit are achieved through space mapping, and the calculation cost is effectively reduced while the calculation precision is guaranteed.

Description

Technical Field

[0001] The present invention belongs to the technical field of concrete construction, and particularly relates to a method, system, device and medium for calculating the cooling temperature field of concrete water pipes. Background Art

[0002] During the solidification process of mass concrete, due to the release of cement hydration heat and the inhomogeneity of the concrete itself, unpredictable temperature deformation and temperature stress usually occur. The magnitude and direction of these temperature deformations and stresses are highly uncertain. When the surface temperature of the concrete differs greatly from the external environmental temperature, temperature cracks are more likely to appear, and some of these cracks may develop into through-cracks, thereby weakening the integrity and strength of the mass concrete structure. Therefore, it is very necessary to effectively control the temperature of mass concrete, and water pipe cooling is one of the effective methods to control the temperature of mass concrete at present. By numerically simulating the cooling temperature field of the water pipes in mass concrete, the temperature stress condition of the concrete can be effectively predicted, so as to adjust the construction plan and reduce the possibility of dam cracking.

[0003] However, there are many problems in the existing calculation methods for the temperature field of mass concrete, especially in the calculation of the heat-fluid coupling method: First, heat-fluid calculations usually require complex numerical calculations, including techniques such as finite element analysis, to simulate the temperature distribution of concrete. For mass concrete, due to the excessive number of elements, the computational time cost is immeasurable; Second, the heat-fluid coupling algorithm usually requires complex and fine mesh generation, especially in the area near the water pipes, where more refined meshing is required. To achieve the coupling of the model, the mesh nodes of the water pipe model need to be combined with the mesh nodes of the concrete model, and this process will also increase the complexity and difficulty of the calculation. These requirements limit the application of the heat-fluid coupling algorithm in dealing with actual engineering models with irregular shapes or structures. Summary of the Invention

[0004] In order to solve the problems of high calculation complexity and difficulty in the existing calculation methods for the temperature field of mass concrete, the present invention provides a method for calculating the cooling temperature field of concrete water pipes to balance the calculation efficiency of the cooling temperature field of mass concrete water pipes and the calculation accuracy of the water pipe temperature, which has important theoretical research significance and engineering application value.

[0005] To achieve the above object, the present invention provides the following technical solutions:

[0006] A method for calculating the cooling temperature field of concrete water pipes, comprising:

[0007] Build a 3D concrete model containing water pipes, and separately perform independent mesh generation on the concrete and water pipes in the 3D concrete model. The mesh nodes of the concrete elements and the water pipe elements after division have no contact relationship; and establish the geometric space mapping relationship between the divided concrete elements and the water pipe elements, and determine the geometric relationship parameters between the concrete elements and the water pipe elements according to the mapping relationship. The geometric space mapping relationship characterizes the number pairing relationship of the water pipe elements passing through the concrete elements, as well as the direction and length of the water pipe elements within the concrete elements;

[0008] According to the geometric relationship parameters between the concrete elements and the water pipe elements, for the concrete elements, calculate the heat exchange amount of the concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete elements themselves, and calculate the concrete temperature through the heat exchange amount of the concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete elements themselves; for the water pipe elements, calculate the heat exchange amount of the water pipes to the concrete, and calculate the water pipe temperature through the heat exchange amount of the water pipes to the concrete and the cooling water flow rate;

[0009] Update the concrete temperature and the water pipe temperature according to the concrete pouring time to obtain the concrete temperature field and the water pipe temperature field.

[0010] Preferably, the geometric relationship parameters include the node numbers and node space coordinates corresponding to the divided concrete elements, the space coordinates of the water pipe nodes, the different water pipe numbers corresponding to both ends of the concrete elements, the length of the water pipes in the concrete, the different concrete numbers corresponding to the two nodes of the water pipe elements, and the water pipe node numbers corresponding to the two nodes of the water pipe elements.

[0011] Preferably, the heat exchange amount H of the concrete cooling the water pipes w is:

[0012] H w = A pipe β(T c - T p );

[0013] In the formula: A = πD pip l unit is the heat exchange area between the unit length water pipe and the concrete; D pipe is the water pipe diameter, l unit is the unit length; β is the convective heat transfer coefficient between the water pipe and the concrete; T c is the concrete temperature, T p is the water pipe temperature;

[0014] The heat release amount H of the hydration heat of the concrete elements themselves y is specifically:

[0015]

[0016] Where: Q0 is the final heat of hydration when τ→∞, Q0 = 270 kJ / kg; m is the parameter for controlling the heat release rate; τ is the concrete pouring time, τ0 is the initial concrete pouring time, and τ0 = 0.

[0017] Preferably, the control equation for the concrete temperature T c (x,y,z,τ) is:

[0018]

[0019] Where: (x, y, z) are the spatial coordinates of the concrete node, τ is the concrete pouring time, λ c is the thermal conductivity of the concrete, c c is the specific heat capacity of the concrete, ρ c is the density of the concrete, H w is the heat exchange amount of the concrete for the water pipe cooling, H y is the heat release amount of the hydration heat of the concrete unit itself.

[0020] Preferably, the control equation for the water pipe temperature T p (x,τ) is:

[0021]

[0022] Where: x is the displacement along the flow direction of the cooling water in the water pipe, τ is the concrete pouring time, λ w is the thermal conductivity of the cooling water, H c is the heat exchange amount of the water pipe for the concrete; ρ w is the density of the cooling water, c w is the specific heat capacity of the cooling water, w w is the flow rate of the cooling water.

[0023] Preferably, the cooling water flows from node M to node N, and the control equations for the water pipe unit nodes M and N are:

[0024]

[0025] Where, is the temperature of the cooling water flowing into node M; v w is the flow velocity of the cooling water; is the temperature of the cooling water flowing out of node N; H q1 is the heat carried away by the cooling water from node M of the water pipe unit, H q2 is the heat carried away by the cooling water from node N;

[0026] The heat carried away by the cooling water from nodes M and N of the water pipe unit are respectively:

[0027]

[0028] In the formula: β is the convective heat transfer coefficient between the water pipe and the concrete, is the temperature at node I of the concrete element; is the temperature at node J of the concrete element.

[0029] The present invention also provides a concrete water pipe cooling temperature field calculation system, including:

[0030] A preprocessing module, which is used to construct a three-dimensional concrete model containing water pipes, independently perform mesh division on the concrete and water pipes in the three-dimensional concrete model, and there is no contact relationship between the mesh nodes of the divided concrete elements and water pipe elements; and establish a geometric space mapping relationship between the divided concrete elements and water pipe elements, and determine the geometric relationship parameters of the concrete elements and water pipe elements according to the mapping relationship; the geometric space mapping relationship characterizes the number pairing relationship of the water pipe elements passing through the concrete elements, as well as the direction and length of the water pipe elements within the concrete elements;

[0031] A heat exchange quantity calculation module, which is used to calculate the heat exchange quantity of the concrete cooling the water pipe and the heat release quantity of the hydration heat of the concrete element itself for the concrete element according to the geometric relationship parameters of the concrete element and the water pipe element, and calculate the concrete temperature through the heat exchange quantity of the concrete cooling the water pipe and the heat release quantity of the hydration heat of the concrete element itself; for the water pipe element, calculate the heat exchange quantity of the water pipe for the concrete, and calculate the water pipe temperature through the heat exchange quantity of the water pipe for the concrete and the cooling water flow rate;

[0032] A temperature field calculation module, which is used to update the concrete temperature and the water pipe temperature according to the concrete pouring time, and obtain the concrete temperature field and the water pipe temperature field.

[0033] The present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory, and the processor executes the computer program to implement the steps of any one of the concrete water pipe cooling temperature field calculation methods.

[0034] The present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is loaded by a processor, it can execute the steps of any one of the concrete water pipe cooling temperature field calculation methods.

[0035] The concrete water pipe cooling temperature field calculation method provided by the present invention has the following beneficial effects:

[0036] In the present invention, the concrete and the water pipe are independently meshed. After meshing, the mesh nodes of the concrete elements and the water pipe elements have no contact relationship, avoiding the complexity of node coupling. At the same time, the heat exchange between the two is processed through a mapping relationship. The independent mesh division reduces the number of meshes. Especially around the water pipe, there is no need for particularly fine meshes, nor is it necessary to align the nodes of the water pipe and the concrete, simplifying the mesh generation process. Especially for irregular structures, there is no need to adjust the meshes to match the nodes of the two. Because the mapping relationship can handle the interaction between different meshes, the overall number of elements is reduced and the calculation time is decreased. At the same time, through the mapping relationship, the heat exchange amount between the concrete elements and the water pipe elements can be calculated separately, avoiding the complex calculations of strong coupling and reducing the calculation difficulty. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the embodiments of the present invention and their design schemes, the accompanying drawings required for the present embodiment will be briefly introduced below. The accompanying drawings in the following description are only partial embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other accompanying drawings can be obtained based on these drawings.

[0038] Figure 1 It is a flowchart of the calculation method for the cooling temperature field of the concrete water pipe in Embodiment 1 of the present invention;

[0039] Figure 2 It is a layout diagram of the casting block and the water pipe;

[0040] Figure 3 It is a distribution diagram of the concrete temperature field at the highest temperature;

[0041] Figure 4 It is a temperature curve diagram of a certain concrete characteristic point;

[0042] Figure 5 It is a flowchart of the preprocessing module;

[0043] Figure 6 It is a schematic diagram of the preprocessing program for the water pipe equivalence. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0044] In order to enable those skilled in the art to better understand the technical solutions of the present invention and implement them, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments. The following embodiments are only used to more clearly illustrate the technical solutions of the present invention and cannot be used to limit the protection scope of the present invention.

[0045] In the description of the present invention, it should be understood that the orientation or positional relationship indicated by terms such as "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "axial", "radial", "circumferential", etc. is based on the orientation or positional relationship shown in the drawings. It is only for the convenience of describing the technical solution of the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation. Therefore, it should not be construed as a limitation to the present invention.

[0046] In addition, terms such as "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying relative importance. In the description of the present invention, it should be noted that unless otherwise clearly specified or limited, the terms "connected" and "connected" should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected or indirectly connected through an intermediate medium. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances. In the description of the present invention, unless otherwise stated, the meaning of "a plurality of" is two or more, which will not be elaborated here.

[0047] Embodiment 1

[0048] The present invention provides a method for calculating the cooling temperature field of a concrete water pipe, as Figure 1 shown, which specifically includes the following steps:

[0049] Step 1: Construct a three-dimensional concrete model containing the water pipe, and independently perform mesh division on the concrete and the water pipe in the three-dimensional concrete model. After division, the mesh nodes of the concrete elements and the water pipe elements have no contact relationship; and establish a geometric space mapping relationship between the divided concrete elements and the water pipe elements, and determine the geometric relationship parameters of the concrete elements and the water pipe elements according to the mapping relationship. The geometric space mapping relationship specifically means pairing the unit numbers of the water pipe elements passing through the concrete elements, and giving the direction and length of the water pipe elements in the concrete elements, etc.

[0050] Specifically, the geometric relationship parameters include the node numbers and node space coordinates corresponding to the divided concrete elements, the space coordinates of the water pipe nodes, the different water pipe numbers corresponding to both ends of the concrete elements, the length of the water pipe in the concrete, the different concrete numbers corresponding to the two nodes of the water pipe element, and the water pipe node numbers corresponding to the two nodes of the water pipe element.

[0051] Through an external call method in which the relationship between concrete and water pipes is independent, the dependence of concrete on the unit grid of water pipes is eliminated, and the spatial mapping association between concrete and water pipes is re-established, achieving the purpose of heat transfer calculation between water pipes and concrete, effectively avoiding problems such as difficult mesh division caused by grid dependence between concrete and water pipes and the resulting poor calculation efficiency.

[0052] Step 2: According to the geometric relationship parameters of concrete units and water pipe units, for concrete units, calculate the heat exchange amount of concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete units themselves, and calculate the concrete temperature through the heat exchange amount of concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete units themselves; for water pipe units, calculate the heat exchange amount of the water pipes to the concrete, and calculate the water pipe temperature through the heat exchange amount of the water pipes to the concrete and the cooling water flow rate.

[0053] The water pipe grid is arranged independently of the concrete unit grid, and the unit nodes are respectively associated with the concrete unit nodes through spatial mapping information, and the heat exchange amount between the concrete unit and the corresponding water pipe node is calculated through their correlation to simulate the convective heat exchange of concrete under the cooling of the water pipes.

[0054] (1) Calculation of concrete unit temperature

[0055] The control equation for the concrete temperature T c (x, y, z, τ) is:

[0056]

[0057] In the formula: (x, y, z) are the spatial coordinates of the concrete nodes, τ is the concrete pouring time, λ c is the thermal conductivity of concrete, c c is the specific heat capacity of concrete, ρ c is the density of concrete, H w is the heat exchange amount of concrete cooling the water pipes, H y is the heat release amount of the hydration heat of the concrete unit itself.

[0058] The heat release amount of the hydration heat of the concrete unit H y Specifically:

[0059]

[0060] In the formula: Q0 is the final hydration heat at τ→∞, Q0 = 270 kJ / kg; m is the parameter for controlling the heat release rate; τ is the concrete pouring time, τ0 is the initial concrete pouring time, τ0 = 0.

[0061] The heat exchange amount of concrete cooling the water pipes H w is:

[0062] Hw = A pipe β(T c - T p );

[0063] In the formula: A = πD pipe l unit is the heat transfer area between the unit - length water pipe and the concrete; D pipe is the diameter of the water pipe, l unit is the unit length; β is the convective heat - transfer coefficient between the water pipe and the concrete; T c is the concrete temperature, T p is the water - pipe temperature.

[0064] (2) Calculation of the water - pipe unit temperature

[0065] The control equation of the water - pipe temperature T p (x, τ) is:

[0066]

[0067] In the formula: x is the displacement along the flow direction of the cooling water in the water pipe, τ is the concrete pouring time, λ w is the thermal conductivity of the cooling water, H c is the heat transfer amount from the water pipe to the concrete; ρ w is the density of the cooling water, c w is the specific heat capacity of the cooling water, v w is the flow rate of the cooling water.

[0068] The cooling water flows from node M to node N. The control equations for the water - pipe unit nodes M and N can be written as:

[0069]

[0070] In the formula, is the temperature of the cooling water flowing into node M; v w is the flow velocity of the cooling water; is the temperature of the cooling water flowing out of node N; H q1 is the heat carried away by the cooling water from node M of the water - pipe unit, H q2 is the heat carried away by the cooling water from node N.

[0071] The heat carried away by the cooling water from nodes M and N of the water - pipe unit are respectively:

[0072]

[0073] In the formula: β is the convective heat - transfer coefficient between the water pipe and the concrete, is the temperature of node I of the concrete unit; is the temperature of node J of the concrete unit.

[0074] Step 3: Update the concrete temperature and the water pipe temperature according to the concrete pouring time to obtain the concrete temperature field and the water pipe temperature field.

[0075] Through the innovative strategy of independent grid division plus geometric space mapping, the present invention breaks through the bottleneck of traditional heat flow coupling calculation, significantly improves the efficiency, accuracy and engineering adaptability of the large-volume concrete temperature field simulation, and is especially applicable to practical projects with complex shapes and multiple water pipe intersections.

[0076] Based on the same inventive concept, the present invention also provides a concrete water pipe cooling temperature field calculation system, including:

[0077] A preprocessing module, which is used to construct a three-dimensional concrete model containing water pipes, independently divide the grids of the concrete and the water pipes in the three-dimensional concrete model respectively, and the grid nodes of the divided concrete elements and water pipe elements have no contact relationship; and establish a geometric space mapping relationship between the divided concrete elements and water pipe elements, and determine the geometric relationship parameters of the concrete elements and water pipe elements according to the mapping relationship; the geometric space mapping relationship characterizes the number pairing relationship of the water pipe elements passing through the concrete elements, as well as the direction and length of the water pipe elements within the concrete elements;

[0078] A heat exchange quantity calculation module, which is used to calculate the heat exchange quantity of the concrete cooling the water pipe and the heat release quantity of the hydration heat of the concrete element itself for the concrete element according to the geometric relationship parameters of the concrete element and the water pipe element, and calculate the concrete temperature through the heat exchange quantity of the concrete cooling the water pipe and the heat release quantity of the hydration heat of the concrete element itself; for the water pipe element, calculate the heat exchange quantity of the water pipe to the concrete, and calculate the water pipe temperature through the heat exchange quantity of the water pipe to the concrete and the cooling water flow rate;

[0079] A temperature field calculation module, which is used to update the concrete temperature and the water pipe temperature according to the concrete pouring time to obtain the concrete temperature field and the water pipe temperature field.

[0080] The present invention adopts the core concept of independent water pipe grid and concrete grid, independently creates the water pipe grid and the grid of the concrete element respectively, and the grid nodes of the concrete and the water pipe in the model have no contact relationship.

[0081] The calculation method of the heat exchange quantity calculation module is: on the basis of the hydration heat of the concrete, call the corresponding water pipe temperature and the water pipe length parameter inside the concrete to calculate the heat exchange quantity, and obtain the temperature field.

[0082] The calculation method of the heat exchange quantity calculation module is: the water pipe calls the heat exchange quantity of the corresponding numbered concrete, and then imports the heat flux caused by the cooling water flow to calculate the water pipe temperature.

[0083] The pre - processing module's required pre - data mainly includes: ① Parameters related to the concrete model grid (element - corresponding node numbers, node spatial coordinates); ② Water pipe nodes (in the form of spatial coordinates). The output results of the pre - processing mainly include: ① Different water pipe numbers corresponding to both ends of the concrete element, the length of the water pipe in the concrete; ② Different concrete numbers corresponding to the two nodes of the water pipe element; ③ Water pipe node numbers corresponding to the two nodes of the water pipe element.

[0084] Each module in the above - mentioned concrete water pipe cooling temperature field calculation system can be implemented in whole or in part by software, hardware, and their combinations. Each of the above - mentioned modules can be embedded in the processor of a computer device in hardware form or be independent of it, or can be stored in the memory of the computer device in software form so that the processor can call and execute the operations corresponding to each of the above - mentioned modules.

[0085] The present invention also provides a computer device, including a memory, a processor, and a computer program stored on the memory. The processor executes the computer program to implement the steps in the embodiment of the concrete water pipe cooling temperature field calculation method. For the specific implementation method, reference can be made to the method embodiment and will not be elaborated here.

[0086] Furthermore, the present invention also provides a non - transitory computer - readable storage medium containing instructions, on which a computer program is stored. For example, a memory containing instructions, and the above - mentioned instructions can be executed by the processor of the computer device to complete the above - mentioned method. For example, the non - transitory computer - readable storage medium can be a ROM, a random access memory (RAM), a CD - ROM, a magnetic tape, a floppy disk, and an optical data storage device, etc. When the computer program is executed by the processor, it can implement the steps in the embodiment of the concrete water pipe cooling temperature field calculation method. For the specific implementation method, reference can be made to the method embodiment and will not be elaborated here.

[0087] The following combines the attached Figure 5 and Figure 6 and describes the specific implementation manner of the present invention in detail. The concrete water pipe cooling temperature field calculation method provided by the present invention includes the following steps:

[0088] The first step is the spatial mapping method of the water pipe element and the concrete element. Refer to Figure 6As shown in the figure, the water pipe passes through concrete units 1, 2, and 4 respectively. The preprocessing module can solve the intersection points of each concrete unit and the water pipe unit, calculate the length Lunit of the water pipe inside the concrete unit based on the distance between the intersection points, and calculate the next unit through which the water pipe passes by obtaining the information of the plane where the intersection points are located, so as to calculate the order in which the water pipe passes through the concrete units. A virtual water pipe unit shown as a light blue dotted line is generated according to the order of the concrete units. Concrete unit 1 corresponds to virtual water pipes 1 and 2 and the length Lunit1 of the water pipe inside the concrete; concrete unit 2 corresponds to virtual water pipes 2 and 3 and the length Lunit2 of the water pipe inside the concrete; concrete unit 4 corresponds to virtual water pipes 3 and 4 and the length Lunit4 of the water pipe inside the concrete. The preprocessing module will output the following preprocessing information based on this information: ① The different water pipe numbers corresponding to both ends of the concrete unit and the length of the water pipe inside the concrete; ② The different concrete numbers corresponding to the two nodes of the water pipe unit; ③ The water pipe node numbers corresponding to the two nodes of the water pipe unit.

[0089] Second step, temperature analysis by the finite element method. The heat exchange quantity calculation module will use the preprocessing information generated by the preprocessing module, call the water pipe temperature corresponding to the concrete and the length Lunit of the water pipe inside the unit to calculate the temperature field of convective heat transfer and hydration heat, and record the total amount of convective heat transfer of the current unit. The heat exchange quantity calculation module retrieves the total amount of convective heat transfer of the corresponding concrete unit, then calculates the temperature of the water pipe unit after considering the heat flux, and stores the temperature in the memory. During the calculation of the concrete temperature field and the water pipe temperature field, the grid calculation of the concrete unit only involves the heat release amount of its own unit's hydration heat and an externally imported convective heat transfer amount; the water pipe unit only involves the heat flux of its own unit and the convective heat transfer amount transferred by the concrete. Among them, the concrete unit grid and the water pipe unit grid are independent of each other and have no direct physical contact. Their heat exchange form is realized in the following way: The concrete unit calls the water pipe temperature of the corresponding water pipe unit stored in the memory, calculates the convective heat transfer amount of its own water pipe, and stores the convective heat transfer amount in the memory; the water pipe unit calls the convective heat transfer amount of the corresponding concrete unit stored in the memory to calculate its own temperature, and stores the temperature in the memory to complete the cycle. This is the specific form of the core concept that the water pipe grid and the concrete grid are independent of each other.

[0090] Calculation example model parameters

[0091] Take a cuboid casting block with 3 layers of metal water pipes buried in the concrete, arranged in a rectangular form, with a water pipe spacing of 2.0 m × 2.0 m and a single water pipe length of 125 m. The four sides of the concrete are adiabatic boundaries, as Figure 2 shown.

[0092] Among them, the concrete pouring temperature is 16 °C, and 8 °C chilled water is passed immediately after pouring for 20 days. The flow direction of the cooling water in the water pipe is changed every 1 day. The secondary water supply starts 90 days after pouring and stops when the temperature of the casting block is reduced to 18 °C.

[0093] The adiabatic temperature rise of concrete is calculated by the following formula:

[0094] θ(τ) = 29τ / (τ + 1.604); where τ is the concrete pouring time.

[0095] Other thermodynamic parameters are shown in the following table:

[0096] Table 1 Other parameter table

[0097]

[0098] Example of fine calculation of concrete water pipe cooling with independent grid:

[0099] 1. Model the concrete and perform grid meshing according to the example, and generate a result file containing all geometric relationship parameters.

[0100] 2. Use the preprocessing module to read the following parameters in the result file: ① Parameters related to the concrete model grid (element corresponding node numbers, node spatial coordinates); ② Water pipe nodes (in the form of spatial coordinates); Generate the corresponding information between the concrete and the water pipe: ① Different water pipe numbers corresponding to both ends of the concrete element, length of the water pipe in the concrete; ② Different concrete numbers corresponding to the two nodes of the water pipe element; ③ Water pipe node numbers corresponding to the two nodes of the water pipe element.

[0101] 3. Use the heat transfer calculation module to read the information generated by the preprocessing module, extract the convective heat transfer amount and heat of hydration release amount for calculating the water pipe temperature in the memory, and store the convective heat transfer amount in the memory.

[0102] 4. Use the heat transfer calculation module to read the convective heat transfer amount of the corresponding concrete element in the memory, calculate the heat flux of its own element to obtain the element temperature field, and finally import the temperature into the memory.

[0103] 5. ABAQUS will automatically read this temperature information, output and generate the temperature field, Figure 3 which is the distribution map of the concrete temperature field at the highest temperature. Figure 4 which is the temperature curve graph of a certain concrete characteristic point, and the coordinates of the characteristic point are (6, 12.5, 1); The present invention can effectively calculate the distribution map of the concrete temperature field, and has certain reference value for the improvement of the engineering example plan.

[0104] The concrete water pipe cooling temperature field calculation method proposed by the present invention effectively solves the core problem in the existing large-volume concrete heat flow coupling calculation through the independent grid division and geometric space mapping strategy. The specific advantages are as follows:

[0105] 1. Solve the problem of high computational time cost

[0106] The concrete and the water pipe grid are independently partitioned. The water pipe units do not need to be aligned with the concrete nodes, avoiding the need to densify the grid to match the nodes. The independent grid partitioning reduces the number of units. The direction and length of the water pipe units within the concrete units are pre-determined through the mapping relationship and directly substituted into the heat transfer formula without complex numerical integration. The traditional method takes about 10 days for calculation due to the huge number of grids. The present invention reduces the number of units by 40% through independent grid partitioning, shortens the calculation time to about 3 days, and the error between the temperature field distribution and the actual measurement is <5%.

[0107] 2. Solve the problem of complex grid partitioning

[0108] The concrete grid can be partitioned according to the natural boundaries of the structure (such as by casting layer, crack position, etc.) without accommodating the shape of the water pipes. The water pipe grid can be independently partitioned along a straight line or a serpentine path, avoiding the forced generation of high-quality tetrahedral grids in complex structures. The different grids are directly correlated using geometric relationship parameters (such as water pipe unit number, penetration direction) without physical node connections. Node coupling is replaced by spatial mapping, and the absence of node contact reduces the difficulty of grid generation.

[0109] 3. Solve the problem of poor adaptability to irregular structures

[0110] The independent grid of the present invention adapts to complex geometries. The concrete grid can flexibly handle irregular boundaries (such as curves, holes), and the water pipe grid can be freely arranged (such as spirals, branches). For example: for special-shaped bridge piers or structures with prestressed ducts, the water pipe grid can be independently generated without changing the concrete grid. The mapping relationship compensates for the grid mismatch. Even if the water pipe units penetrate multiple concrete units, the heat transfer amount of each unit can still be calculated separately through the mapping relationship, avoiding errors caused by grid mismatch. For example: when a serpentine water pipe penetrates multiple concrete units, the heat transfer amount of each unit is distributed according to the length ratio of the water pipe within that unit.

[0111] Those skilled in the art should understand that the embodiments of the present invention can provide a method, a system, or a computer program product. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) containing computer-usable program codes.

[0112] The present invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram can be implemented by computer program instructions, as well as the combination of flows and / or blocks in the flowchart and / or block diagram. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing devices generate means for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or multiple blocks.

[0113] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory generate a manufactured article including instruction means, and the instruction means implements the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or multiple blocks.

[0114] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or multiple blocks.

[0115] It should be noted that the above specific embodiments can enable those skilled in the art to understand the present invention more comprehensively, but do not limit the present invention in any way. Therefore, although this specification and the embodiments have described the present invention in detail, those skilled in the art should understand that the present invention can still be modified or equivalently replaced; and all technical solutions and improvements that do not depart from the spirit and scope of the present invention are covered by the protection scope of the patent of the present invention. Any reference signs in the claims should not be regarded as limiting the claimed claims. Any simple changes or equivalent replacements of technical solutions that can be obviously obtained by those skilled in the art within the technical scope disclosed by the present invention all belong to the protection scope of the present invention.

Claims

1. A method for calculating the cooling temperature field of a concrete water pipe, characterized in that, Including: Construct a three-dimensional concrete model including water pipes, independently perform mesh division on the concrete and water pipes in the three-dimensional concrete model, and there is no contact relationship between the mesh nodes of the divided concrete elements and water pipe elements; and establish a geometric space mapping relationship between the divided concrete elements and water pipe elements, and determine the geometric relationship parameters of the concrete elements and water pipe elements according to the mapping relationship; the geometric space mapping relationship characterizes the number pairing relationship of the water pipe elements passing through the concrete elements, as well as the direction and length of the water pipe elements within the concrete elements; According to the geometric relationship parameters of the concrete elements and water pipe elements, for the concrete elements, calculate the heat exchange amount of the concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete elements themselves, and calculate the concrete temperature through the heat exchange amount of the concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete elements themselves; for the water pipe elements, calculate the heat exchange amount of the water pipes to the concrete, and calculate the water pipe temperature through the heat exchange amount of the water pipes to the concrete and the cooling water flow rate; Update the concrete temperature and water pipe temperature according to the concrete pouring time to obtain the concrete temperature field and the water pipe temperature field.

2. The calculation method of the cooling temperature field of the concrete water pipe according to claim 1, characterized in that The geometric relationship parameters include the node numbers and node space coordinates corresponding to the divided concrete elements, the space coordinates of the water pipe nodes, the different water pipe numbers corresponding to both ends of the concrete elements, the length of the water pipes in the concrete, the different concrete numbers corresponding to the two nodes of the water pipe elements, and the water pipe node numbers corresponding to the two nodes of the water pipe elements.

3. The method for calculating the cooling temperature field of the concrete water pipe according to claim 2, wherein The heat exchange quantity H of concrete for cooling the water pipe w is as follows: H w = A pipe β(T c - T p ); Where: A = πD pipe l unit is the heat exchange area between the unit length of the water pipe and the concrete; D pipe is the diameter of the water pipe, l unit is the unit length; β is the convective heat transfer coefficient between the water pipe and the concrete; T c is the concrete temperature, T p is the water pipe temperature; The heat release amount H of the hydration heat of the concrete unit itself y Specifically: In the formula: Q0 is the final hydration heat at τ→∞, Q0 = 270 kJ / kg; m is the parameter controlling the heat release rate; τ is the concrete pouring time, τ0 is the initial concrete pouring time, τ0 = 0.

4. The method for calculating the cooling temperature field of a concrete water pipe according to claim 3, wherein The temperature T of the concrete c (x, y, z, τ) has the following governing equation: Where: (x, y, z) are the spatial coordinates of the concrete joint, τ is the concrete pouring time, and λ c is the concrete thermal conductivity, c c is the specific heat capacity of the concrete, ρ c is the concrete density, H w is the heat exchange amount of the concrete for the water pipe cooling, H y is the heat release amount of the hydration heat of the concrete element itself.

5. The method for calculating the cooling temperature field of a concrete water pipe according to claim 4, characterized in that The temperature T of the water pipe p (x,τ) has the following governing equation: Where: x is the displacement along the flow direction of the cooling water in the water pipe, τ is the concrete pouring time, and λ w is the thermal conductivity of the cooling water, and H c is the heat exchange amount of the water pipe to the concrete; ρ w is the density of the cooling water, c w is the specific heat capacity of the cooling water, and w w is the flow rate of the cooling water.

6. The calculation method of the cooling temperature field of the concrete water pipe according to claim 5, wherein, The cooling water flows from node M to node N, and the control equations for the water pipe element nodes M and N are: In the formula, is the temperature of the cooling water flowing into node M; v w is the cooling water flow rate; is the temperature of the cooling water flowing out of node N; H q1 is the heat taken away by the cooling water from node M of the water pipe unit, H q2 is the heat taken away by the cooling water from node N; The heat carried away by the cooling water from the water pipe element nodes M and N is respectively: where: β is the convective heat transfer coefficient between the water pipe and the concrete; is the temperature at Node I of the concrete element; is the temperature at Node J of the concrete element.

7. A calculation system for the cooling temperature field of a concrete water pipe, characterized in that Including: A preprocessing module for constructing a three-dimensional concrete model including water pipes, independently performing mesh division on the concrete and water pipes in the three-dimensional concrete model, and there is no contact relationship between the mesh nodes of the divided concrete elements and water pipe elements; and establishing a geometric space mapping relationship between the divided concrete elements and water pipe elements, and determining the geometric relationship parameters of the concrete elements and water pipe elements according to the mapping relationship; the geometric space mapping relationship characterizes the number pairing relationship of the water pipe elements passing through the concrete elements, as well as the direction and length of the water pipe elements within the concrete elements; A heat exchange amount calculation module for calculating the heat exchange amount of the concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete elements themselves for the concrete elements according to the geometric relationship parameters of the concrete elements and water pipe elements, and calculating the concrete temperature through the heat exchange amount of the concrete cooling the water pipes and the heat release amount of the hydration heat of the concrete elements themselves; for the water pipe elements, calculating the heat exchange amount of the water pipes to the concrete, and calculating the water pipe temperature through the heat exchange amount of the water pipes to the concrete and the cooling water flow rate; A temperature field calculation module for updating the concrete temperature and water pipe temperature according to the concrete pouring time to obtain the concrete temperature field and the water pipe temperature field.

8. A computer device, comprising a memory, a processor, and a computer program stored on the memory, characterized in that, The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is loaded by the processor, it is capable of executing the steps of the method according to any one of claims 1 to 6.

Images