A numerical simulation method for complex stress aqueduct structures in industrial atmospheric environment

Abstract

The invention provides a numerical simulation method for a complex stressed aqueduct structure in an industrial atmospheric environment, and relates to the technical field of numerical simulation of aqueduct structures. The method comprises: obtaining each subcomponent of a target aqueduct structure to obtain a subcomponent list A; obtaining a constitutive relationship corresponding to a test piece corresponding to each subcomponent in A in the industrial atmospheric environment where the target aqueduct structure is located to obtain a constitutive relationship list set X corresponding to A; establishing a target aqueduct structure model corresponding to the target aqueduct structure according to an initial constitutive relationship corresponding to each subcomponent in A, and running the target aqueduct structure model to obtain a stress cloud map and a damage cloud map corresponding to the target aqueduct structure model; if the damage cloud map corresponding to the target aqueduct structure has no damage area, then according to X and the stress cloud map corresponding to the target aqueduct structure model, the weak area and the dangerous time of the target aqueduct structure are accurately determined, so as to provide a reference for the detection and maintenance of the aqueduct structure.

Description

Technical Field

[0001] The invention relates to the technical field of numerical simulation of aqueduct structures, and in particular to a numerical simulation method for aqueduct structures subjected to complex forces in industrial atmospheric environments. Background Art

[0002] As the main building material, the performance of aqueduct gradually degrades with the change of time and environmental conditions. When the aqueduct is in an industrial atmosphere, the industrial atmosphere contains a large amount of acidic gas. After the acidic gas mixes with the moisture in the air, it reacts with the cement stone and aggregate inside the aqueduct, reducing the alkalinity of the aqueduct, causing corrosion of the aqueduct structure, and leading to the degradation of the aqueduct strength.

[0003] When the aqueduct is subjected to external forces, cracks or changes in pore structure may occur inside it. These changes may change the permeability of the aqueduct, making it easier or harder for acidic substances in the industrial atmosphere to enter the aqueduct, thereby accelerating or slowing down the erosion process of the aqueduct. For example, under tension, the pores of the aqueduct expand, and acidic substances are more likely to corrode the aqueduct through the gaps. Conversely, under compression, the pores of the aqueduct shrink, which can hinder or delay the corrosion of the aqueduct through the gaps by acidic substances. That is, at the same time, the aqueduct under different stress states is corroded by acidic substances to different degrees, and its actual strength and performance change to different degrees.

[0004] In engineering practice, for an aqueduct structure that is subjected to complex forces in an industrial atmospheric environment, its strength will degrade over time. In order to prevent a part of the aqueduct structure from breaking or collapsing, the aqueduct structure can be monitored in real time. However, when the aqueduct structure is large, real-time monitoring of all its areas is time-consuming and labor-intensive. On the other hand, real-time monitoring is meaningless in the early stages of the construction of the aqueduct structure, because the initial strength is very high and there is no need for monitoring. Therefore, how to accurately determine the weak areas of the aqueduct structure and the time period of damage has become a technical problem that needs to be solved urgently. Summary of the invention

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

[0006] According to a first aspect of the present application, a numerical simulation method for a complex stressed aqueduct structure in an industrial atmospheric environment is provided, the method comprising the following steps:

[0007] S100, obtaining each subcomponent of the target aqueduct structure to obtain a subcomponent list A=(A1, A2, ..., A i , …, A n ), i=1, 2,...,n; among them, A i is the i-th subcomponent of the target aqueduct structure, and n is the number of subcomponents of the target aqueduct structure;

[0008] S200, obtaining the constitutive relations of the specimens corresponding to each subcomponent in A in the industrial atmospheric environment where the target aqueduct structure is located, so as to obtain a constitutive relation list set X=(X1, X2, …, X i , …, X n ); where X i A i The corresponding constitutive relation list; X i =(X i 1，1 , X i 1，2 ,…,X i 1，p , …, X i 1，q , …, X i j，1 , X i j，2 ,…,X i j，p , …, X i j，q , …, X i m,1 , X i m，2 , …, X i m，p , …, X i m，q ),j=1,2,…,m;p=1,2,…,q;X i j，p A i The corresponding constitutive relationship after the corresponding specimen is applied with p times of the preset force F in the industrial atmosphere environment where the target aqueduct structure is located and maintained for j preset time lengths T; m is the number of preset time lengths; q is the preset multiple of the applied preset force;

[0009] S300, establishing a target aqueduct structure model corresponding to the target aqueduct structure according to the initial constitutive relationship corresponding to each subcomponent in A, and running the target aqueduct structure model to obtain a stress cloud map and a damage cloud map corresponding to the target aqueduct structure model;

[0010] S400: If there is no damaged area in the damage cloud map corresponding to the target aqueduct structure, determine the weak area and dangerous time of the target aqueduct structure according to X and the stress cloud map corresponding to the target aqueduct structure model.

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

[0012] The numerical simulation method for a complex stressed aqueduct structure in an industrial atmospheric environment of the present invention obtains each subcomponent of a target aqueduct structure to obtain a subcomponent list A; obtains the constitutive relationship corresponding to a test piece corresponding to each subcomponent in A in the industrial atmospheric environment where the target aqueduct structure is located to obtain a constitutive relationship list set X corresponding to A; establishes a target aqueduct structure model corresponding to the target aqueduct structure according to the initial constitutive relationship corresponding to each subcomponent in A, and runs the target aqueduct structure model to obtain a stress cloud map and a damage cloud map corresponding to the target aqueduct structure model; if the damage cloud map corresponding to the target aqueduct structure has no damage area, then according to X and the stress cloud map corresponding to the target aqueduct structure model, the weak area and the dangerous time of the target aqueduct structure are accurately determined, so as to provide a reference for the detection and maintenance of the aqueduct structure. BRIEF DESCRIPTION OF THE DRAWINGS

[0013] 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.

[0014] Figure 1 A flow chart of a numerical simulation method for a complex stressed aqueduct structure in an industrial atmospheric environment provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0015] The technical solutions in the embodiments of the present application are described clearly and completely below. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by technicians in the relevant technical field without creative work are within the scope of protection of the present application.

[0016] It is noteworthy that in the following description, specific details such as specific system structures, technologies, etc. are proposed for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present application. However, it should be clear to those skilled in the art that the present application can also be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to prevent unnecessary details from hindering the description of the present application.

[0017] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of the described features, integrals, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components and / or their collections. It should also be understood that when the term "and / or" is used in the present specification and the appended claims, it refers to any combination of one or more of the associated listed items and all possible combinations, and includes these combinations.

[0018] In the description of the present specification and the appended claims, if the term "if" is used, it can be interpreted as "in the case of" or "when" or "uponce" or "in response to determining" or "in response to detecting" according to the context. Similarly, the phrase "if it is determined" or "if [the described condition or event] is detected" can be interpreted as meaning "uponce it is determined" or "in response to determining" or "uponce [the described condition or event] is detected" or "in response to detecting [the described condition or event]" according to the context.

[0019] In the description of this application specification and the appended claims, the terms "first", "second", "third", etc., are only used to distinguish the description and are not to be understood as indicating or implying relative importance, nor are they used to describe a specific order or sequence.

[0020] It is worth noting that the "f(i)" used in this application cannot be simply understood as a fixed numerical value. Instead, it should be understood that the value of "f(i)" changes with the value of i. For example, if the value range of i is 1 to 2, "f(1)" may be equal to "f(2)" or may not be equal to "f(2)". That is, the number of XX information in different XX information lists may be the same or different, and is generally determined according to the actual application scenario. Other similar expressions such as "g(j)" should be understood in the same way and will not be elaborated on.

[0021] References to "one embodiment" or "some embodiments" etc. described in the specification of this application mean that one or more embodiments of the present application include specific features, structures or characteristics described in conjunction with the embodiment. Therefore, the phrases "in one embodiment", "in some embodiments", "in some other embodiments", "in some other embodiments", etc. appearing in different places in this specification do not necessarily refer to the same embodiment, but mean "one or more but not all embodiments", unless otherwise specifically emphasized in other ways.

[0022] It should be understood that the order of execution of the steps in the following embodiments does not necessarily mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0023] The following will refer to Figure 1 The flowchart of the numerical simulation method for complex stress-bearing aqueduct structure in industrial atmospheric environment shown in the figure introduces a numerical simulation method for complex stress-bearing aqueduct structure in industrial atmospheric environment.

[0024] The numerical simulation method for a complex stressed aqueduct structure in an industrial atmospheric environment may include the following steps:

[0025] S100, obtaining each subcomponent of the target aqueduct structure to obtain a subcomponent list A=(A1, A2, ..., A i , …, A n ), i=1, 2,...,n; among them, A i is the i-th subcomponent of the target aqueduct structure, and n is the number of subcomponents of the target aqueduct structure.

[0026] In this embodiment, the target aqueduct structure can be an aqueduct structure that has been put into use, or it can be an aqueduct structure that has been completed and is about to be put into use; it can be understood that the aqueduct structure can be a single relatively regular aqueduct structure, for example: a rectangular aqueduct structure; the aqueduct structure can also be an aqueduct structure composed of multiple sub-components; each sub-component of the aqueduct structure can be obtained to obtain A.

[0027] S200, obtaining the constitutive relations of the specimens corresponding to each subcomponent in A in the industrial atmospheric environment where the target aqueduct structure is located, so as to obtain a constitutive relation list set X=(X1, X2, …, X i , …, X n ); where X i A i The corresponding constitutive relation list; X i =(X i 1，1 , X i 1，2 ,…,X i 1，p , …, X i 1，q , …, X i j，1 , X i j，2 ,…,X i j，p , …, X i j，q , …, X i m,1, X i m，2 , …, X i m，p , …, X i m，q ),j=1,2,…,m;p=1,2,…,q;X i j，p A i The corresponding constitutive relationship of the corresponding specimen after applying p times the preset force F in the industrial atmospheric environment where the target aqueduct structure is located and maintaining it for j preset time periods T; m is the number of preset time periods; q is the preset multiple of the applied preset force.

[0028] In this embodiment, a number of test pieces with the same shape as each sub-component can be made, and then the test pieces corresponding to each sub-component are placed in an industrial atmospheric environment that is the same as the industrial atmospheric environment in which the target aqueduct structure is located, and a preset multiple of a preset force is applied to each test piece, for example: 1 times the preset force F is applied to the first test piece corresponding to the first sub-component, and 2 times the preset force F is applied to the second test piece corresponding to the first sub-component; similarly, the retention time of each test piece corresponding to each sub-component in the industrial atmospheric environment can be controlled, for example: the first test piece corresponding to the first sub-component is kept in the industrial atmospheric environment for 1 preset time period, and the second test piece corresponding to the first sub-component is kept in the industrial atmospheric environment for 2 preset time periods; thereby, a list of constitutive relations corresponding to each sub-component can be obtained; it should be noted that those skilled in the art can use existing constitutive relation determination methods to determine the constitutive relations of each test piece according to actual needs, which will not be elaborated here.

[0029] Furthermore, the preset force applied by the subcomponent may be pressure or tension; specifically, if a subcomponent is actually subjected to pressure, pressure is applied to the corresponding specimen; if a subcomponent is actually subjected to tension, tension is applied to the corresponding specimen.

[0030] Furthermore, each subcomponent has the same structural shape as the corresponding test piece, thereby avoiding the influence of the structural shape on the constitutive relationship, making the constitutive relationship of the corresponding subcomponent obtained through the test piece more accurate.

[0031] S300, establishing a target aqueduct structure model corresponding to the target aqueduct structure according to the initial constitutive relationship corresponding to each subcomponent in A, and running the target aqueduct structure model to obtain a stress cloud map and a damage cloud map corresponding to the target aqueduct structure model.

[0032] In this embodiment, the target aqueduct structure model corresponding to the target aqueduct structure is composed of several subcomponents. The more subcomponents there are, the closer a single subcomponent is to being subjected to unidirectional force. The initial constitutive relationship corresponding to each subcomponent in A can be determined by a standard aqueduct specimen under normal conditions. The remaining parameters of the subcomponents can be determined by the actual materials used, and then loads and constraints are added according to the actual situation, and the model is run to obtain the stress cloud map and damage cloud map of the aqueduct structure model.

[0033] S400: If there is no damaged area in the damage cloud map corresponding to the target aqueduct structure, determine the weak area and dangerous time of the target aqueduct structure according to X and the stress cloud map corresponding to the target aqueduct structure model.

[0034] Further, step S400 may include the following steps:

[0035] S410, obtaining a preset value W=1 and an intermediate aqueduct structure model QR.

[0036] S420, setting QR as the target aqueduct structure model after operation.

[0037] S430, according to the stress cloud diagram corresponding to each subcomponent of QR, determine the maximum principal stress of each subcomponent to obtain a subcomponent maximum principal stress list σ max =(σ max，1 , σ max，2 , …, σ max，i , …, σ max，n ); where σ max，i A i The corresponding maximum principal stress; σ max，i =MAX(σ i，x ,σ i，y ,σ i，z );σ i，x A i The stress component in the X-axis direction, σ i，y A i The stress component in the Y-axis direction, σ i，z A i The stress component in the Z-axis direction.

[0038] In this embodiment, σ i，x A i The average stress of all elements on the surface in the X-axis direction; σ i，y A i The average stress of all elements on the surface in the Y-axis direction; σ i，z A iThe average stress of all elements on the surface in the Z-axis direction; the stress of each section of the subcomponent can be directly viewed through the view section manager. The more subcomponents there are, the closer it is to unidirectional force, the more obvious the maximum principal stress is, but the calculation amount is greater.

[0039] S440, traversal σ max , if r×F / D≤σ max，i <(r+1)×F / D, then replace A in QR i The constitutive relation of is changed to X i j，r , to get the Wth updated aqueduct model; D is σ max，i Corresponding cross-sectional area; j=1, 2,…, m; r=1, 2,…, q-1; enter S450.

[0040] In this embodiment, the initial value of i is 1.

[0041] S450, running the Wth update of the aqueduct model.

[0042] S460, set QR to the Wth updated aqueduct model after running, and enter S430; assign W to W+1.

[0043] S470: If the constitutive relations of all subcomponents in QR are not changed, proceed to S480.

[0044] S480, determining the damaged area of ​​the damage cloud map corresponding to the Wth updated aqueduct model after running as a weak area, and obtaining the number NUM' of preset retention time lengths corresponding to the subcomponents corresponding to the weak area.

[0045] S490, determining NUM'×T as the danger time corresponding to the target aqueduct structure.

[0046] In this embodiment, if the Wth updated aqueduct model after running does not need to change the constitutive relationship of any sub-component, it means that the Wth updated aqueduct model after running at this time is closest to the actual situation of the target aqueduct structure remaining in the industrial atmospheric environment for a corresponding time.

[0047] The damaged area of ​​the damage cloud map corresponding to the Wth updated aqueduct model after this operation can be obtained, the damaged area can be determined as a weak area, and the number NUM' of the preset maintenance time corresponding to the subcomponents corresponding to the weak area can be obtained, so that the dangerous time corresponding to the target aqueduct structure can be obtained as NUM'×T.

[0048] Furthermore, the weak areas can be inspected and repaired intensively within the time of NUM'×T to avoid the fracture or collapse of the aqueduct structure.

[0049] Furthermore, after step S400, the method may further include the following steps:

[0050] S500: If a damage area exists in the damage cloud map corresponding to the target aqueduct structure, the damage area in the damage cloud map is determined as a weak area corresponding to the target aqueduct structure.

[0051] In this embodiment, after running the target aqueduct structure model corresponding to the target aqueduct structure, if there is a damaged area in the damage cloud map corresponding to the target aqueduct structure, it means that the target aqueduct structure model corresponding to the target aqueduct structure has been damaged; at this time, the damaged area in the damage cloud map can be directly determined as the weak area corresponding to the target aqueduct structure, and no subsequent simulation operation is required.

[0052] S600: Obtain the number NUM of preset retention time periods corresponding to the subcomponents corresponding to the weak area.

[0053] S700, determining NUM×T as the danger time corresponding to the target aqueduct structure.

[0054] In this embodiment, through the above steps, the number NUM of preset durations corresponding to the destruction areas in the destruction cloud corresponding to the target aqueduct structure can be obtained, so that the danger time corresponding to the target aqueduct structure can be obtained as NUM×T.

[0055] In this embodiment, the aqueduct in the pressure area is less corroded by acidic substances in the industrial atmospheric environment, and the strength degradation of the aqueduct is slower. Therefore, the constitutive relationship of the aqueduct area with different stress states is continuously replaced through iteration to reflect the influence of the stress state on the strength degradation of the aqueduct. The weak areas of the aqueduct structure and the dangerous time of damage are judged through the aqueduct damage cloud map, and the aqueduct structure model with complex stress in the industrial atmospheric environment is simulated. It provides a reference for the monitoring and maintenance of the aqueduct structure.

[0056] The numerical simulation method for a complex stressed aqueduct structure in an industrial atmospheric environment of the present embodiment obtains each subcomponent of the target aqueduct structure to obtain a subcomponent list A; obtains the constitutive relationship corresponding to the test piece corresponding to each subcomponent in A in the industrial atmospheric environment where the target aqueduct structure is located to obtain a constitutive relationship list set X corresponding to A; establishes a target aqueduct structure model corresponding to the target aqueduct structure according to the initial constitutive relationship corresponding to each subcomponent in A, and runs the target aqueduct structure model to obtain a stress cloud map and a damage cloud map corresponding to the target aqueduct structure model; if the damage cloud map corresponding to the target aqueduct structure has no damage area, then according to X and the stress cloud map corresponding to the target aqueduct structure model, the weak area and dangerous time of the target aqueduct structure are accurately determined, thereby providing a reference for the detection and maintenance of the aqueduct structure.

[0057] 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.

[0058] 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.

[0059] 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.

[0060] 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.

[0061] 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.

[0062] 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).

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

[0064] 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.

[0065] 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).

[0066] 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.

[0067] 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).

[0068] 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.

[0069] 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.

[0070] 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.

[0071] 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.

[0072] 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.

[0073] 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 complex stress-bearing aqueduct structures in industrial atmospheric environments, characterized in that: The method comprises the following steps: S100, obtaining each subcomponent of the target aqueduct structure to obtain a subcomponent list A=(A1, A2, ..., A i , …, A n ), i=1, 2,...,n; among them, A i is the i-th subcomponent of the target aqueduct structure, and n is the number of subcomponents of the target aqueduct structure; S200, obtaining the constitutive relations of the specimens corresponding to each subcomponent in A in the industrial atmospheric environment where the target aqueduct structure is located, so as to obtain a constitutive relation list set X=(X1, X2, …, X i , …, X n ); where X i A i The corresponding constitutive relation list; X i =(X i 1，1 , X i 1，2 ,…,X i 1，p , …, X i 1，q , …, X i j，1 , X i j，2 ,…,X i j，p , …, X i j，q , …, X i m,1 , X i m，2 , …, X i m，p , …, X i m，q ),j=1,2,…,m;p=1,2,…,q;X i j，p A i The corresponding constitutive relationship after the corresponding specimen is applied with p times of the preset force F in the industrial atmosphere environment where the target aqueduct structure is located and maintained for j preset time lengths T; m is the number of preset time lengths; q is the preset multiple of the applied preset force; S300, establishing a target aqueduct structure model corresponding to the target aqueduct structure according to the initial constitutive relationship corresponding to each subcomponent in A, and running the target aqueduct structure model to obtain a stress cloud map and a damage cloud map corresponding to the target aqueduct structure model; S400, if the damage cloud map corresponding to the target aqueduct structure has no damage area, then determine the weak area and dangerous time of the target aqueduct structure according to X and the stress cloud map corresponding to the target aqueduct structure model; After step S400, the method further comprises the following steps: S500, if the damage cloud map corresponding to the target aqueduct structure has a damage area, the damage area in the damage cloud map is determined as the weak area corresponding to the target aqueduct structure; S600, obtaining the number NUM of preset retention time lengths corresponding to the subcomponents corresponding to the weak area; S700, determining NUM×T as the danger time corresponding to the target aqueduct structure; Step S400 includes the following steps: S410, obtaining a preset value W=1 and a middle aqueduct structure model QR; S420, setting QR as the target aqueduct structure model after operation; S430, according to the stress cloud diagram corresponding to each subcomponent of QR, determine the maximum principal stress of each subcomponent to obtain a subcomponent maximum principal stress list σ max =(σ max，1 , σ max，2 , …, σ max，i , …, σ max，n ); where σ max，i A i The corresponding maximum principal stress; σ max，i =MAX(σ i，x ,σ i，y ,σ i，z );σ i，x A i The stress component in the X-axis direction, σ i，y A i The stress component in the Y-axis direction, σ i，z A i The stress component in the Z-axis direction; S440, traversal σ max , if r×F / D≤σ max，i <(r+1)×F / D, then replace A in QR i The constitutive relation of is changed to X i j，r , to get the Wth updated aqueduct model; D is σ max，i Corresponding cross-sectional area; j=1, 2, ..., m; r=1, 2, ..., q-1; F is the preset force; enter S450; S450, running the Wth update of the aqueduct model; S460, setting QR to the Wth updated aqueduct model after running, and entering S430; assigning W to W+1; S470, if the constitutive relations of all subcomponents in QR are not changed, proceed to S480; S480, determining the damaged area of ​​the damage cloud map corresponding to the Wth updated aqueduct model after running as a weak area, and obtaining the number NUM' of preset retention time lengths corresponding to the subcomponents corresponding to the weak area; S490, determining NUM'×T as the danger time corresponding to the target aqueduct structure.

2. The numerical simulation method for complex stress-bearing aqueduct structure in industrial atmospheric environment according to claim 1 is characterized in that: The preset force F is tension or compression.

3. The numerical simulation method for complex stress-bearing aqueduct structure in industrial atmospheric environment according to claim 1, characterized in that: Each subcomponent has the same structural shape as the corresponding test piece.

4. The numerical simulation method for complex stress-bearing aqueduct structure in industrial atmospheric environment according to claim 1, characterized in that: σ i，x A i The average stress of all elements on the surface in the X-axis direction; σ i，y A i The average stress of all elements on the surface in the Y-axis direction; σ i，z A i The average stress of all elements on the surface in the Z-axis direction.

5. The numerical simulation method for complex stress-bearing aqueduct structure in industrial atmospheric environment according to claim 1, characterized in that: The operational target aqueduct structure model includes: Add loads and constraints according to the actual usage conditions corresponding to the target aqueduct structure, and run the target aqueduct structure model.

Images