A numerical simulation method for aqueduct structure

Abstract

The present invention provides a numerical simulation method for aqueduct structures, which relates to the field of numerical simulation of aqueduct structures. The method includes: determining the highest temperature T corresponding to the aqueduct structure max and the lowest temperature T min ; obtaining a specified temperature list T; acquiring a number of preset concrete parameters of the concrete cube corresponding to the aqueduct structure at each specified temperature in T to obtain a concrete parameter list set A corresponding to T; performing a temperature field simulation on the finite element model WR corresponding to the aqueduct structure according to T to obtain a temperature field list B corresponding to the aqueduct structure; obtaining the loading duration corresponding to each temperature field in B to obtain a loading duration list t corresponding to B; re-performing a simulation on WR and updating the concrete parameters corresponding to the temperature field during the simulation to determine the simulation value of the aqueduct structure. The method of the present invention can improve the accuracy of the simulation value of the aqueduct structure.

Description

Technical Field

[0001] The invention relates to the field of numerical simulation of aqueduct structures, and in particular to a numerical simulation method for aqueduct structures. Background Art

[0002] Aqueducts, especially large aqueducts, are all spatial plate-beam structures with complex structures and thin-walled and irregular shapes. The inner surface is the water boundary, and the other surfaces are the boundaries in contact with the air. Due to the difference in water temperature and air temperature in the same season, and the different deformation constraint boundaries of each surface, this will inevitably lead to certain temperature stress in the aqueduct structure. Looking back at the aqueducts built at home and abroad, almost all of them are cracked, and temperature stress is one of the main reasons. For thin-walled structures such as aqueducts, the temperature stress generated will inevitably cause microcracks on the surface of the aqueduct concrete, and it is very easy to develop into through cracks, threatening the normal operation of the aqueduct and reducing the service life of the aqueduct.

[0003] Existing numerical simulation methods generally use local meteorological data, combined with actual initial conditions and boundary conditions, to calculate the two-dimensional temperature field of the aqueduct structure and couple it with the stress field. When the above method is used to simulate the numerical value of the aqueduct, the various parameters of the concrete material do not change with changes in temperature and other conditions, and the cracking state of the concrete is not considered. Therefore, the simulation value of the aqueduct obtained by the above method is of low accuracy. Summary of the invention

[0004] In view of the above technical problems, the technical solution adopted by the present invention is:

[0005] According to a numerical simulation method for aqueduct structure provided in the present application, the method comprises the following steps:

[0006] S100, determining the maximum temperature T corresponding to the aqueduct structure according to the historical meteorological data of the area where the aqueduct structure is located within a preset historical time period max and the minimum temperature T min ;

[0007] S200, for T min and T max The corresponding temperature intervals are evenly divided to obtain the specified temperature list T = (T1, T2, ..., T i , …, T n ), i=1, 2, …, n; where T i For T min and T max The corresponding temperature interval is divided into the i-th specified temperature, n is the T min and T max The number of specified temperatures obtained by dividing the corresponding temperature interval; T1=T min , Tn =T max ;

[0008] S300, obtaining a number of preset concrete parameters of a concrete cube corresponding to the aqueduct structure at each specified temperature in T, so as to obtain a concrete parameter list set A corresponding to T = (A1, A2, ..., A i , …, A n ); where A i The concrete cube corresponding to the aqueduct structure is at T i List of concrete parameters below;

[0009] S400, based on T, a temperature field simulation is performed on the finite element model WR corresponding to the aqueduct structure to obtain a temperature field list B=(B1, B2, …, B p , …, B q ), p = 1, 2, …, q; where B p is the pth temperature field obtained by simulating the temperature field of WR, q is the number of temperature fields obtained by simulating the temperature field of WR; TB r+1 -TB r =α×(T n -T1) / (n-1); TB r To simulate B r The maximum temperature corresponding to the aqueduct structure at 1 0000 ℃; r = 1, 2, ..., q-1; α is the preset weight;

[0010] S500, obtaining the loading duration corresponding to each temperature field in B, so as to obtain a loading duration list t=(t1, t2, ..., t p ,…,t q ); where t p To get B p The corresponding loading time;

[0011] S600, re-simulate WR according to T, A and t, and update the concrete parameters corresponding to the temperature field during the simulation process to determine the simulation values ​​of the aqueduct structure.

[0012] Further, step S600 includes the following steps:

[0013] S610, obtaining a preset value H=1 and an intermediate model QR;

[0014] S620, obtain QR=WR;

[0015] S630, simulate QR and set the loading time history. When the loading time reaches t H When , stop the simulation to get the Hth initial model W H ;

[0016] S640, according to B H , W H The upper temperature is T i -(T n -T1) / (2×(n-1)) to T i + (T n -T1) / (2×(n-1)) is divided into the i-th sub-region X i ;

[0017] S650, X i The corresponding concrete parameter is set to A i , to get W H The corresponding Hth target model U H ;

[0018] S660, if H < q and the maximum tensile stress of QR is less than or equal to the concrete tensile stress corresponding to the aqueduct structure, then obtain QR = U H , and H=H+1, proceed to S630;

[0019] S670, if H<q and the maximum tensile stress of QR is greater than the tensile stress of the concrete corresponding to the aqueduct structure, then select a unit whose tensile stress is greater than the tensile stress of the concrete, cut the unit along the tensile stress direction with a preset plane, so that the unit generates two separated surfaces, and set a hard contact and frictionless contact relationship on the newly generated contact surface; obtain H=H+1, and enter S630;

[0020] S680, if H=q, stop the simulation and determine the crack position of the aqueduct structure based on the simulation results.

[0021] Furthermore, the preset concrete parameters include: elastic modulus, Young's stiffness, Poisson's ratio, shear modulus, thermal conductivity, specific heat capacity, thermal expansion coefficient and density.

[0022] Furthermore, α<1.

[0023] Furthermore, the initial parameters of WR are set according to the influence range of the water flow temperature in the aqueduct structure; wherein, the concrete parameters within the influence range of the water flow temperature are assigned to the concrete parameters at the corresponding temperature, and the concrete parameters outside the influence range of the water flow temperature are assigned to the concrete parameters under the initial conditions.

[0024] Furthermore, the water flow temperature influence range is measured by an infrared thermal imager.

[0025] The present invention has at least the following beneficial effects:

[0026] The numerical simulation method for aqueduct structure of the present invention obtains several preset concrete parameters of the concrete cube corresponding to the aqueduct structure at each specified temperature to obtain a concrete parameter list set while only considering the temperature influence; according to each specified temperature, the temperature field analysis is performed on the finite element model corresponding to the aqueduct structure to obtain a temperature field list corresponding to the aqueduct structure and a loading time corresponding to each temperature field; then, the finite element model of the aqueduct structure is re-simulated, and the concrete parameters of the corresponding area of ​​the temperature field are updated for each loading time corresponding to the temperature field, so as to obtain the simulation value of the aqueduct structure; in the process of simulating the finite element model of the aqueduct structure, the changes in the temperature field and the corresponding changes in the concrete parameters are simultaneously considered, so that the simulation process is closer to the real environment, thereby improving the accuracy of the simulation value of the aqueduct structure. BRIEF DESCRIPTION OF THE DRAWINGS

[0027] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0028] Figure 1 A flow chart of a numerical simulation method for an aqueduct structure provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0029] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0030] It should be noted that, based on the present disclosure, those skilled in the art should understand that an aspect described herein can be implemented independently of any other aspect, and two or more of these aspects can be combined in various ways. For example, any number of aspects described herein can be used to implement a device and / or practice a method. In addition, other structures and / or functionalities other than one or more of the aspects described herein can be used to implement this device and / or practice this method.

[0031] The following will refer to Figure 1 The flowchart of the numerical simulation method for aqueduct structure is shown, and a numerical simulation method for aqueduct structure is introduced.

[0032] The principle of the method of this embodiment is as follows:

[0033] It is understandable that the goal of thermal field solution is temperature field. After obtaining the temperature field result, the temperature change of each node of the unit will be introduced as thermal load into the analysis of stress field. The temperature field and stress field are coupled and calculated through the thermal expansion coefficient of the material to obtain the stress component of each unit node. It should be noted that the various parameters of concrete materials change with temperature and other conditions. Therefore, if the concrete parameters at the corresponding temperature can be given to the concrete of each part of the aqueduct structure in real time, the influence of temperature effect can be better reflected.

[0034] Based on the above principles, the numerical simulation method for aqueduct structure may include the following steps:

[0035] S100, determining the maximum temperature T corresponding to the aqueduct structure according to the historical meteorological data of the area where the aqueduct structure is located within a preset historical time period max and the minimum temperature T min .

[0036] In this embodiment, the aqueduct structure is usually built in a fixed area. The historical meteorological data of the area in a preset historical period can be obtained through historical records, so that the maximum temperature and the minimum temperature of the area where the aqueduct structure is located can be obtained; the functional relationship between the temperature of the aqueduct structure and the temperature of the area can be obtained by fitting a large amount of data. Therefore, the maximum temperature T corresponding to the aqueduct structure can be further obtained. max and the minimum temperature T min ; The preset historical time period can be set according to actual needs, for example: the preset historical time period is 5 years or 10 years.

[0037] S200, for T min and T max The corresponding temperature intervals are evenly divided to obtain the specified temperature list T = (T1, T2, ..., T i ,…,T n ), i=1, 2, …, n; where T i For T min and T max The corresponding temperature interval is divided into the i-th specified temperature, n is the T min and T max The number of specified temperatures obtained by dividing the corresponding temperature interval; T1=T min , T n =T max .

[0038] In this embodiment, after obtaining T min and T max After that, you can min and T maxThe corresponding temperature interval is evenly divided to obtain T; that is, the temperature difference between any two adjacent specified temperatures in T is the same.

[0039] S300, obtaining a number of preset concrete parameters of a concrete cube corresponding to the aqueduct structure at each specified temperature in T, so as to obtain a concrete parameter list set A corresponding to T = (A1, A2, ..., A i , …, A n ); where A i The concrete cube corresponding to the aqueduct structure is at T i List of concrete parameters below.

[0040] In this embodiment, the specific parameters of the concrete used in the construction of the aqueduct structure are known, and therefore, the concrete cube corresponding to the aqueduct structure can be obtained; further, a number of preset concrete parameters of the concrete cube corresponding to the aqueduct structure at each specified temperature in T can be measured in a laboratory environment; it should be noted that those skilled in the art can use the existing concrete parameter measurement method to measure a number of preset concrete parameters of the concrete cube corresponding to the aqueduct structure at each specified temperature in T according to actual needs, which will not be elaborated here.

[0041] Furthermore, the preset concrete parameters include: elastic modulus, Young's stiffness, Poisson's ratio, shear modulus, thermal conductivity, specific heat capacity, thermal expansion coefficient and density.

[0042] S400, based on T, a temperature field simulation is performed on the finite element model WR corresponding to the aqueduct structure to obtain a temperature field list B=(B1, B2, …, B p , …, B q ), p = 1, 2, …, q; where B p is the pth temperature field obtained by simulating the temperature field of WR, q is the number of temperature fields obtained by simulating the temperature field of WR; TB r+1 -TB r =α×(T n -T1) / (n-1); TB r To simulate B r The maximum temperature corresponding to the aqueduct structure at q-1; r=1, 2, …, q-1; α is the preset weight.

[0043] In this embodiment, for the aqueduct structure, a corresponding finite element model can be established, and the heat conduction equation can be determined, and the initial conditions and boundary conditions are substituted to consider the temperature field of the aqueduct structure; when water flows through the aqueduct structure, the temperature of the concrete in contact with it is equal to the water temperature when water flows through it, and the influence range of the water temperature is determined by tools such as infrared thermal imagers, and the concrete parameters within the influence range are assigned to the material parameters at the corresponding temperature, and the rest are assigned to the material parameters under the initial conditions; for example, the concrete parameters of the range affected by the water flow temperature on the inner side of the aqueduct structure can be specifically set according to the water temperature, and the water temperature can be obtained according to the relationship between the water temperature and the ambient temperature, and the range affected by the ambient temperature on the outer side can be set according to the ambient temperature.

[0044] The temperature field of WR is simulated under the condition that only the temperature effect is considered. When the maximum temperature increases by α×(T n -T1) / (n-1), record the corresponding loading time, set the loading time from the beginning to the first maximum temperature rise α×(T n The loading time is t1, and the second maximum temperature rises by α×(T n -T1) / (n-1) has a loading time of t2, and so on, to get B.

[0045] Furthermore, α<1, the value of α can be set according to actual needs, and the smaller α is, the higher the simulation accuracy is.

[0046] S500, obtaining the loading duration corresponding to each temperature field in B, so as to obtain a loading duration list t=(t1, t2, ..., t p ,…,t q ); where t p To get B p The corresponding loading time.

[0047] In this embodiment, in the above step S400, the loading duration corresponding to each temperature field is recorded. Therefore, the loading duration corresponding to each temperature field in B can be obtained to obtain the loading duration list t corresponding to B.

[0048] S600, re-simulate WR according to T, A and t, and update the concrete parameters corresponding to the temperature field during the simulation process to determine the simulation values ​​of the aqueduct structure.

[0049] Further, step S600 may include the following steps:

[0050] S610, obtaining a preset value H=1 and a preset intermediate model QR.

[0051] S620, obtain QR=WR.

[0052] S630, simulate QR and set the loading time history. When the loading time reaches t H When , stop the simulation to get the Hth initial model W H .

[0053] S640, according to B H , W H The upper temperature is T i -(T n -T1) / (2×(n-1)) to T i + (T n -T1) / (2×(n-1)) is divided into the i-th sub-region X i .

[0054] S650, X i The corresponding concrete parameter is set to A i , to get W H The corresponding Hth target model U H .

[0055] S660, if H < q and the maximum tensile stress of QR is less than or equal to the concrete tensile stress corresponding to the aqueduct structure, then obtain QR = U H , and H=H+1, enter S630.

[0056] S670, if H<q and the maximum tensile stress of QR is greater than the tensile stress of the concrete corresponding to the aqueduct structure, then select the unit whose tensile stress is greater than the tensile stress of the concrete, cut the unit along the tensile stress direction with a preset plane, so that the unit produces two separated surfaces, and set hard contact and frictionless contact relationship on the newly generated contact surface; obtain H=H+1, and enter S630.

[0057] S680, if H=q, stop the simulation and determine the crack position of the aqueduct structure based on the simulation results.

[0058] In this embodiment, the loading time history is set. When the loading time reaches t1, the calculation is stopped and the model at this time is derived. According to the calculated temperature field, combined with T, the model is divided into regions in the cross-sectional direction. i -(T n -T1) / (2×(n-1)) to T i + (T n -T1) / (2×(n-1)) is divided into a region and named as the i-th sub-region X i , give X i The concrete gives A i The corresponding concrete parameters are recalculated. When the loading time reaches t2, the model is exported, the area is re-divided, the concrete parameters are re-assigned, and the calculation is repeated again.

[0059] Assume that after k cycles, the maximum tensile stress in the aqueduct is greater than the concrete tensile stress. The model at this time is exported, and the unit with tensile stress greater than the concrete tensile stress is selected. The unit is cut along the stress direction with each tiny plane to produce two separated surfaces. Hard contact and frictionless contact relationship are set on the newly generated contact surface, and then the model is imported and recalculated. In this way, the unit at this location is prevented from transmitting tensile stress and cracking is approximately simulated.

[0060] Furthermore, a calculation stop condition can be set, such as when the concrete crack grows to a certain length, the calculation can be stopped.

[0061] Conventional numerical simulation methods are difficult to consider the effects of temperature stress and concrete cracking. The method in this embodiment takes into account the effect of carbonization on concrete and can more accurately simulate the performance of reinforced concrete. At the same time, as the crack width of the aqueduct structure continues to increase, when the crack reaches a certain length, the structure can be considered to be damaged, which has a certain predictive effect on the service life of the structure.

[0062] In this embodiment, when only the temperature influence is considered, several preset concrete parameters of the concrete cube corresponding to the aqueduct structure at each specified temperature are obtained to obtain a concrete parameter list set; according to each specified temperature, the temperature field analysis is performed on the finite element model corresponding to the aqueduct structure to obtain a temperature field list corresponding to the aqueduct structure and a loading time corresponding to each temperature field; then, the finite element model of the aqueduct structure is re-simulated, and the concrete parameters of the corresponding area of ​​the temperature field are updated for each loading time corresponding to the temperature field, so as to obtain the simulated numerical value of the aqueduct structure; in the process of simulating the finite element model of the aqueduct structure, the changes in the temperature field and the corresponding changes in the concrete parameters are simultaneously considered, so that the simulation process is closer to the real environment, thereby improving the accuracy of the simulated numerical value of the aqueduct structure.

[0063] In addition, although the steps of the method in the present disclosure are described in a specific order in the drawings, this does not require or imply that the steps must be performed in this specific order, or that all the steps shown must be performed to achieve the desired results. Additionally or alternatively, some steps may be omitted, multiple steps may be combined into one step, and / or one step may be decomposed into multiple steps, etc.

[0064] An embodiment of the present invention also provides a non-transitory computer-readable storage medium, which can be set in an electronic device to store at least one instruction or at least one program related to implementing a method in a method embodiment. The at least one instruction or the at least one program is loaded and executed by the processor to implement the method provided in the above embodiment.

[0065] The program product may adopt any combination of one or more readable media. The readable medium may be a readable signal medium or a readable storage medium. The readable storage medium may be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples (non-exhaustive list) of readable storage media include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0066] Computer readable signal media may include data signals propagated in baseband or as part of a carrier wave, in which readable program code is carried. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. Readable signal media may also be any readable medium other than a readable storage medium, which may send, propagate, or transmit a program for use by or in conjunction with an instruction execution system, apparatus, or device.

[0067] The program code embodied on the readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wired, optical cable, RF, etc., or any suitable combination of the foregoing.

[0068] Program code for performing the operations of the present application may be written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Java, C++, etc., and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user computing device, partially on the user device, as a separate software package, partially on the user computing device and partially on a remote computing device, or entirely on a remote computing device or server. In the case of a remote computing device, the remote computing device may be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computing device (e.g., using an Internet service provider to connect through the Internet).

[0069] An embodiment of the present invention further provides an electronic device, comprising a processor and the aforementioned non-transitory computer-readable storage medium.

[0070] The electronic device is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0071] The electronic device is presented in the form of a general-purpose computing device. The components of the electronic device may include, but are not limited to: the at least one processor mentioned above, the at least one memory mentioned above, and a bus connecting different system components (including the memory and the processor).

[0072] The memory stores program codes, which can be executed by the processor, so that the processor executes the steps in various embodiments described in this specification.

[0073] The memory may include readable media in the form of volatile memory, such as random access memory (RAM) and / or cache memory, and may further include read only memory (ROM).

[0074] The memory may also include a program / utility having a set (at least one) of program modules, such program modules including but not limited to: an operating system, one or more application programs, other program modules and program data, each of which or some combination may include an implementation of a network environment.

[0075] The bus may represent one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, or a local bus using any of a variety of bus architectures.

[0076] The electronic device may also communicate with one or more external devices (e.g., keyboards, pointing devices, Bluetooth devices, etc.), may communicate with one or more devices that enable a user to interact with the electronic device, and / or may communicate with any device (e.g., routers, modems, etc.) that enables the electronic device to communicate with one or more other computing devices. Such communication may be performed through an input / output (I / O) interface. Furthermore, the electronic device may also communicate with one or more networks (e.g., local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) through a network adapter. The network adapter communicates with other modules of the electronic device through a bus. It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in conjunction with the electronic device, including but not limited to: microcode, device drivers, redundant processors, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, etc.

[0077] Through the description of the above implementation, it is easy for those skilled in the art to understand that the example implementation described here can be implemented by software, or by software combined with necessary hardware. Therefore, the technical solution according to the implementation of the present disclosure can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a USB flash drive, a mobile hard disk, etc.) or on a network, including several instructions to enable a computing device (which can be a personal computer, a server, a terminal device, or a network device, etc.) to execute the method according to the implementation of the present disclosure.

[0078] An embodiment of the present invention further provides a computer program product, which includes program code. When the program product is run on an electronic device, the program code is used to enable the electronic device to execute the steps of the method according to various exemplary embodiments of the present invention described above in this specification.

[0079] Although some specific embodiments of the present invention have been described in detail by way of example, it should be understood by those skilled in the art that the above examples are only for illustration, not for limiting the scope of the present invention. It should also be understood by those skilled in the art that various modifications may be made to the embodiments without departing from the scope and spirit of the present invention.

Claims

1. A numerical simulation method for aqueduct structure, characterized in that: The method comprises the following steps: S100, determining the maximum temperature T corresponding to the aqueduct structure according to the historical meteorological data of the area where the aqueduct structure is located within a preset historical time period max and the minimum temperature T min ; S200, for T min and T max The corresponding temperature intervals are evenly divided to obtain the specified temperature list T = (T1, T2, ..., T i , …, T n ), i=1, 2, …, n; where T i For T min and T max The corresponding temperature interval is divided into the i-th specified temperature, n is the T min and T max The number of specified temperatures obtained by dividing the corresponding temperature interval; T1=T min , T n =T max ; S300, obtaining a number of preset concrete parameters of a concrete cube corresponding to the aqueduct structure at each specified temperature in T, so as to obtain a concrete parameter list set A corresponding to T = (A1, A2, ..., A i , …, A n ); where A i The concrete cube corresponding to the aqueduct structure is at T i List of concrete parameters below; S400, based on T, a temperature field simulation is performed on the finite element model WR corresponding to the aqueduct structure to obtain a temperature field list B=(B1, B2, …, B p , …, B q ), p = 1, 2, …, q; where B p is the pth temperature field obtained by simulating the temperature field of WR, q is the number of temperature fields obtained by simulating the temperature field of WR; TB r+1 -TB r =α×(T n -T1) / (n-1); TB r To simulate B r The maximum temperature corresponding to the aqueduct structure at 1 0000 ℃; r = 1, 2, ..., q-1; α is the preset weight; S500, obtaining the loading duration corresponding to each temperature field in B, so as to obtain a loading duration list t=(t1, t2, ..., t p ,…,t q ); where t p To get B p The corresponding loading time; S600, re-simulate WR according to T, A and t, and update the concrete parameters corresponding to the temperature field during the simulation process to determine the simulation values ​​of the aqueduct structure.

2. The numerical simulation method for aqueduct structure according to claim 1, characterized in that: Step S600 includes the following steps: S610, obtaining a preset value H=1 and a preset intermediate model QR; S620, obtain QR=WR; S630, simulate QR and set the loading time history. When the loading time reaches t H When , stop the simulation to get the Hth initial model W H ; S640, according to B H , W H The upper temperature is T i -(T n -T1) / (2×(n-1)) to T i + (T n -T1) / (2×(n-1)) is divided into the i-th sub-region X i ; S650, X i The corresponding concrete parameter is set to A i , to get W H The corresponding Hth target model U H ; S660, if H < q and the maximum tensile stress of QR is less than or equal to the concrete tensile stress corresponding to the aqueduct structure, then obtain QR = U H , and H=H+1, proceed to S630; S670, if H<q and the maximum tensile stress of QR is greater than the tensile stress of the concrete corresponding to the aqueduct structure, then select a unit whose tensile stress is greater than the tensile stress of the concrete, cut the unit along the tensile stress direction with a preset plane, so that the unit generates two separated surfaces, and set a hard contact and frictionless contact relationship on the newly generated contact surface; obtain H=H+1, and enter S630; S680, if H=q, stop the simulation and determine the crack position of the aqueduct structure based on the simulation results.

3. The numerical simulation method for aqueduct structure according to claim 1, characterized in that: The preset concrete parameters include: elastic modulus, Young's stiffness, Poisson's ratio, shear modulus, thermal conductivity, specific heat, thermal expansion coefficient and density.

4. The numerical simulation method for aqueduct structure according to claim 1, characterized in that: α＜1。 5. The numerical simulation method for aqueduct structure according to claim 1, characterized in that: The initial parameters of WR are set according to the influence range of the water flow temperature in the aqueduct structure; among them, the concrete parameters within the influence range of the water flow temperature are assigned to the concrete parameters at the corresponding temperature, and the concrete parameters outside the influence range of the water flow temperature are assigned to the concrete parameters under the initial conditions.

6. The numerical simulation method for aqueduct structure according to claim 5, characterized in that: The water flow temperature influence range is measured by an infrared thermal imager.

Images