Concrete structure cracking analysis method and system based on extended finite element method

Abstract

The present invention relates to the technical field of numerical simulation of fracture mechanics, and provides a method and system for analyzing cracking of concrete structure based on extended finite element method, wherein the method comprises: establishing a grid model representing the concrete structure; constructing an extended finite element grid model for the grid model, including drawing geometric parameter conditions of the grid model, drawing geometric parameters of the analysis object and the location of the preset crack; applying boundary conditions to the extended finite element grid model, performing grid division, and then defining the physical and mechanical parameters of the grid model; adding a first single degree of freedom unit and a second single degree of freedom unit of the boundary conditions that can reflect the cracking parameters of the crack to the extended finite element grid model to form an extended finite element analysis model; performing extended finite element calculation through the extended finite element analysis model, observing the cracking process of the grid model, analyzing the stress characteristics of the model, and realizing the cracking analysis of the concrete structure. The present invention can accurately simulate the cracking process of the concrete structure.

Description

Technical Field

[0001] The invention relates to the technical field of fracture mechanics numerical simulation, and in particular to a concrete structure cracking analysis method and system based on an extended finite element method. Background Art

[0002] The traction separation law used in the numerical simulation calculation of concrete cracking will cause the softening part to reduce the local strength of the corresponding cracked unit. When the stress reaches the peak strength, the quasi-static finite element calculation cannot converge to the equilibrium solution.

[0003] The following difficulties exist when using the traditional finite element method to model discontinuous surfaces such as cracks: first, the geometric discontinuity of the mesh model needs to be divided; second, the mesh model needs to be finely divided to describe the singular asymptotic field on the crack propagation path; third, the mesh information needs to be continuously updated to map the crack geometry. The extended finite element method is an extended form of the traditional finite element method based on the concept of unit decomposition. In the extended finite element method, the discontinuity of the crack is guaranteed by enrichment functions and additional degrees of freedom, and the unit decomposition allows the local enrichment function to be included in the finite element approximation. As a very effective method for simulating cracks, the extended finite element method can simulate the initiation and propagation of discrete cracks along any path without the need to re-mesh the damaged material to match the geometry of the discontinuous position. However, the extended finite element method has high requirements for the quality of mesh model division, and the convergence of the model is difficult to guarantee. Summary of the invention

[0004] The purpose of the present invention is to solve at least one technical problem in the background technology and to provide a method and system for analyzing cracking of concrete structures based on an extended finite element method.

[0005] To achieve the above object, the present invention provides a method for analyzing cracking of concrete structure based on extended finite element method, comprising:

[0006] Create a mesh model to represent the concrete structure;

[0007] An extended finite element mesh model is constructed for the mesh model, including drawing geometric parameter conditions of the mesh model, drawing geometric parameters of the analysis object, and the location of the preset cracks;

[0008] Apply boundary conditions to the extended finite element mesh model, perform mesh division, and after the mesh division is completed, define the physical and mechanical parameters of the mesh model;

[0009] Define the boundary conditions and structural load conditions of the extended finite element mesh model;

[0010] Adding a first single-degree-of-freedom unit and a second single-degree-of-freedom unit that can reflect the boundary conditions of cracking parameters of the crack to the extended finite element mesh model to form an extended finite element analysis model;

[0011] The extended finite element analysis model is used to perform extended finite element calculations, observe the cracking process of the grid model, analyze the stress characteristics of the model, and realize the cracking analysis of the concrete structure.

[0012] According to one aspect of the present invention, a first single-degree-of-freedom unit and a second single-degree-of-freedom unit that can reflect the boundary conditions of the cracking parameters of the crack are added to the extended finite element mesh model to form an extended finite element analysis model:

[0013] Based on the extended finite element mesh model, The relational expression defines a first single degree of freedom unit Element1 and a second single degree of freedom unit Element2 that can reflect the boundary conditions of cracking parameters of cracks in the extended finite element mesh model;

[0014] Where f is the vertical force, u is the vertical displacement, N C is the control node, N1 is the crack tip node, and y represents the vertical direction.

[0015] According to one aspect of the present invention, and The mapping relationship between the crack tip node N1 vertical displacement and the control node external force state is established through the transpose of the stiffness matrix defined by the first single degree of freedom unit Element1, as shown in the formula: ;

[0016] Small numbers close to zero are defined in the diagonal of the stiffness matrix to avoid singularities.

[0017] According to one aspect of the present invention, for controlling the node displacement and the outer boundary N L The mapping relationship between the vertical loads of the nodes is established by defining the stiffness matrix of the second single-degree-of-freedom unit Element2, as shown in the formula: ;

[0018] By controlling the vertical force at the node and the external boundary N L The displacement of the node can define the tip of the crack, and by extending the finite element, the relationship between the nodes and the control nodes of the entire potential crack damage area can be defined.

[0019] According to one aspect of the present invention, for the formula The mapping relationship can be expressed as By defining the stiffness matrix in , the average displacement of the crack propagation path is linearly related to the vertical force in the control node.

[0020] To achieve the above object, the present invention also provides a concrete structure cracking analysis system based on extended finite element method, comprising:

[0021] A mesh model building module builds a mesh model representing the concrete structure;

[0022] An extended finite element mesh model construction module is used to construct an extended finite element mesh model for the mesh model, including drawing geometric parameter conditions of the mesh model, drawing geometric parameters of the analysis object, and the location of the preset cracks;

[0023] An extended finite element mesh model parameter definition module applies boundary conditions to the extended finite element mesh model, performs mesh division, and defines the physical and mechanical parameters of the mesh model after the mesh division is completed;

[0024] Define the boundary conditions and structural load conditions of the extended finite element mesh model;

[0025] An extended finite element analysis model forming module is configured to add a first single degree of freedom unit and a second single degree of freedom unit capable of reflecting boundary conditions of cracking parameters of cracks in the extended finite element mesh model to form an extended finite element analysis model;

[0026] The concrete structure cracking analysis module performs extended finite element calculation through the extended finite element analysis model, observes the cracking process of the grid model, analyzes the stress characteristics of the model, and realizes the concrete structure cracking analysis.

[0027] To achieve the above-mentioned purpose, the present invention also provides an electronic device, including a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein when the computer program is executed by the processor, the method for analyzing cracking of a concrete structure based on the extended finite element method as described above is implemented.

[0028] To achieve the above objectives, the present invention also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the method for analyzing cracking of a concrete structure based on the extended finite element method as described above is implemented.

[0029] According to the technical solution of the present invention, the present invention can achieve the following beneficial effects:

[0030] The invention can accurately simulate the cracking process of brittle materials by improving the convergence of the extended finite element mesh model, adopts the large-scale general finite element software Abaqus to perform feasibility verification analysis, and the convergence control relation is compiled through the INP file.

[0031] The present invention can accurately reflect the unstable state after cracking by improving the mesh model (extended finite element analysis model) of convergence, which proves its improvement on convergence. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] Figure 1 A flowchart schematically showing a method for analyzing cracking of a concrete structure based on an extended finite element method according to an embodiment of the present invention;

[0033] Figure 2 This is the main view of the concrete structure mesh model;

[0034] Figure 3 It is the side view of the concrete structure mesh model;

[0035] Figure 4 Schematic diagram of the mesh model to improve convergence;

[0036] Figure 5 Preset structure diagram for 2D model cracks;

[0037] Figure 6 This is a two-dimensional model crack extension structure diagram;

[0038] Figure 7 Preset structural diagrams for 3D model cracks;

[0039] Figure 8 This is a diagram of the crack expansion structure of the three-dimensional model. DETAILED DESCRIPTION

[0040] The present invention will now be discussed with reference to exemplary embodiments. It should be understood that the embodiments discussed are only to enable those skilled in the art to better understand and thus implement the present invention, rather than implying any limitation on the scope of the present invention.

[0041] As used herein, the term “including” and variations thereof are to be interpreted as open-ended terms meaning “including, but not limited to.” The term “based on” is to be interpreted as “based, at least in part, on.” The terms “one embodiment” and “an embodiment” are to be interpreted as “at least one embodiment.”

[0042] Figure 1 The flowchart of the method for analyzing the cracking of concrete structure based on the extended finite element method according to one embodiment of the present invention is schematically shown. Figure 1 As shown, in this embodiment, the concrete structure cracking analysis method based on the extended finite element method includes:

[0043] Create a mesh model to represent the concrete structure;

[0044] An extended finite element mesh model is constructed for the mesh model, including drawing the geometric parameter conditions of the mesh model and drawing the location of the preset cracks;

[0045] Apply boundary conditions to the extended finite element mesh model, perform mesh division, and after the mesh division is completed, define the physical and mechanical parameters of the mesh model;

[0046] Define the boundary conditions and structural load conditions of the extended finite element mesh model;

[0047] Adding a first single-degree-of-freedom unit and a second single-degree-of-freedom unit that can reflect the boundary conditions of cracking parameters of the crack to the extended finite element mesh model to form an extended finite element analysis model;

[0048] By extending the finite element analysis model to perform extended finite element calculations, the cracking process of the grid model is observed, the stress characteristics of the model are analyzed, and the cracking analysis of the concrete structure is realized.

[0049] In this embodiment, a mesh model representing the concrete structure is established, that is, the mesh model of the concrete structure is constructed by using ABAQUS software.

[0050] In this embodiment, the geometric parameter conditions for drawing the grid model are to draw the length, width and height of the concrete component. The location of drawing the preset crack is to draw the initial location of the crack.

[0051] In this embodiment, the physical and mechanical parameters of the grid model are defined to include: elastic modulus, bulk density, compressive strength, tensile strength, and Poisson's ratio.

[0052] In this embodiment, the boundary conditions are constraint conditions of the outer boundaries of the extended finite element mesh model, and the structural force load conditions are external force load conditions applied to the extended finite element mesh model.

[0053] In this embodiment, the extended finite element mesh model is constructed by ABAQUS software, and the constructed extended finite element mesh model is meshed and various parameter conditions are added by ABAQUS software. In addition, after adding the first single degree of freedom unit and the second single degree of freedom unit that can reflect the boundary conditions of the cracking parameters of the crack in the user-defined module of ABAQUS software, the extended finite element analysis model that can simulate the development and evolution of the crack is formed.

[0054] Further, according to an embodiment of the present invention, a first single degree of freedom unit and a second single degree of freedom unit that can reflect the boundary conditions of the cracking parameters of the crack are added to the extended finite element mesh model to form an extended finite element analysis model:

[0055] Based on the extended finite element mesh model, The relational expression defines a first single degree of freedom unit Element1 and a second single degree of freedom unit Element2 that can reflect the boundary conditions of cracking parameters of cracks in the extended finite element mesh model;

[0056] Where f is the vertical force, u is the vertical displacement, N C is the control node, N1 is the crack tip node, and y represents the vertical direction.

[0057] Further, according to one embodiment of the present invention, for and The mapping relationship between the crack tip node N1 vertical displacement and the control node external force state is established through the transpose of the stiffness matrix defined by the first single degree of freedom unit Element1, as shown in the formula: ;

[0058] Small numbers close to 0 are defined in the diagonal of the stiffness matrix to avoid singularities.

[0059] In this embodiment, through the definition of the above matrix, the displacement of the crack tip can be made linearly related to the vertical force in the control node.

[0060] Further, according to an embodiment of the present invention, for controlling the node displacement and the outer boundary N L The mapping relationship between the vertical loads of the nodes is established by defining the stiffness matrix of the second single-degree-of-freedom unit Element2, as shown in the formula: ;

[0061] By controlling the vertical force at the node and the external boundary N L The displacement of the node can define the tip of the crack, and by extending the finite element, the relationship between the nodes and the control nodes of the entire potential crack damage area can be defined.

[0062] Further, according to one embodiment of the present invention, for the formula The mapping relationship can be expressed as By defining the stiffness matrix in , the average displacement of the crack propagation path is linearly related to the vertical force in the control node.

[0063] Furthermore, to achieve the above-mentioned purpose, the present invention also provides a concrete structure cracking analysis system based on extended finite element method, comprising:

[0064] A mesh model building module builds a mesh model representing the concrete structure;

[0065] An extended finite element mesh model construction module is used to construct an extended finite element mesh model for the mesh model, including drawing geometric parameter conditions of the mesh model, drawing geometric parameters of the analysis object, and the location of the preset cracks;

[0066] An extended finite element mesh model parameter definition module applies boundary conditions to the extended finite element mesh model, performs mesh division, and defines physical and mechanical parameters of the mesh model after the mesh division is completed;

[0067] Define the boundary conditions and structural load conditions of the extended finite element mesh model;

[0068] An extended finite element analysis model forming module is configured to add a first single degree of freedom unit and a second single degree of freedom unit capable of reflecting boundary conditions of cracking parameters of cracks in the extended finite element mesh model to form an extended finite element analysis model;

[0069] The concrete structure cracking analysis module performs extended finite element calculation through the extended finite element analysis model, observes the cracking process of the grid model, analyzes the stress characteristics of the model, and realizes the concrete structure cracking analysis.

[0070] In this embodiment, a mesh model representing the concrete structure is established, that is, the mesh model of the concrete structure is constructed by using ABAQUS software.

[0071] In this embodiment, the geometric parameter conditions for drawing the grid model are to draw the length, width and height of the concrete component. The location of drawing the preset crack is to draw the initial location of the crack.

[0072] In this embodiment, the physical and mechanical parameters of the grid model are defined to include: elastic modulus, bulk density, compressive strength, tensile strength, and Poisson's ratio.

[0073] In this embodiment, the boundary conditions are constraint conditions of the outer boundaries of the extended finite element mesh model, and the structural force load conditions are external force load conditions applied to the extended finite element mesh model.

[0074] In this embodiment, the extended finite element mesh model is constructed by ABAQUS software, and the constructed extended finite element mesh model is meshed and various parameter conditions are added by ABAQUS software. In addition, after adding the first single degree of freedom unit and the second single degree of freedom unit that can reflect the boundary conditions of the cracking parameters of the crack in the user-defined module of ABAQUS software, the extended finite element analysis model that can simulate the development and evolution of the crack is formed.

[0075] Further, according to an embodiment of the present invention, a first single degree of freedom unit and a second single degree of freedom unit that can reflect the boundary conditions of the cracking parameters of the crack are added to the extended finite element mesh model to form an extended finite element analysis model:

[0076] Based on the extended finite element mesh model, The relational expression defines a first single degree of freedom unit Element1 and a second single degree of freedom unit Element2 that can reflect the boundary conditions of cracking parameters of cracks in the extended finite element mesh model;

[0077] Where f is the vertical force, u is the vertical displacement, N C is the control node, N1 is the crack tip node, and y represents the vertical direction.

[0078] Further, according to one embodiment of the present invention, for and The mapping relationship between the crack tip node N1 vertical displacement and the control node external force state is established through the transpose of the stiffness matrix defined by the first single degree of freedom unit Element1, as shown in the formula: ;

[0079] Small numbers close to zero are defined in the diagonal of the stiffness matrix to avoid singularities.

[0080] In this embodiment, through the definition of the above matrix, the displacement of the crack tip can be made linearly related to the vertical force in the control node.

[0081] Further, according to an embodiment of the present invention, for controlling the node displacement and the outer boundary N L The mapping relationship between the vertical loads of the nodes is established by defining the stiffness matrix of the second single-degree-of-freedom unit Element2, as shown in the formula: ;

[0082] By controlling the vertical force at the node and the external boundary N L The displacement of the node can define the tip of the crack, and by extending the finite element, the relationship between the nodes and the control nodes of the entire potential crack damage area can be defined.

[0083] Further, according to one embodiment of the present invention, for the formula The mapping relationship can be expressed as By defining the stiffness matrix in , the average displacement of the crack propagation path is linearly related to the vertical force in the control node.

[0084] Furthermore, to achieve the above-mentioned purpose, the present invention also provides an electronic device, including a processor, a memory, and a computer program stored in the memory and executable on the processor, wherein when the computer program is executed by the processor, the method for analyzing cracking of a concrete structure based on the extended finite element method as described above is implemented.

[0085] Furthermore, to achieve the above-mentioned purpose, the present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the concrete structure cracking analysis method based on the extended finite element method as described above is implemented.

[0086] Furthermore, based on the above-mentioned solution of the present invention, the effectiveness of the technical solution of the present invention is verified in the form of a specific embodiment in combination with the accompanying drawings. Example

[0087] In order to verify whether the extended finite element mesh model with improved convergence (i.e., extended finite element analysis model) can accurately simulate the cracking process of brittle materials, a feasibility verification analysis was conducted using the large-scale general finite element software ABAQUS, and the convergence control relationship was compiled through the INP file. The indoor model test studied by ERDOGAN and SIH was selected as the model background, and two 2D and 3D mesh models were established to verify the versatility of the extended finite element network model with improved convergence.

[0088] The numerical experimental specimen model (grid model) is 6m long, 1m thick, and 3m wide. The specific model dimensions are shown in Figure 2 and Figure 3 At the same time, in order to reflect the fracture behavior of the specimen under the mixed loading mode, in the static analysis step, the upper and lower boundaries are set with forced displacement boundary conditions to the left and right (0.01m), and at the same time, the upper and lower boundaries are set with forced displacement boundary conditions to the upward and downward directions (8×10-4m).

[0089] The material parameters used in the mesh model are shown in Table 1. The maximum principal stress failure criterion is selected in the initial stage of crack damage, and the damage evolution criterion based on fracture energy is selected in the damage extension stage.

[0090]

[0091] For model calculation verification, first, define the control node through the INP file input method of "*NODE,NSET" as follows. It can be placed anywhere in the model. Next, define the first reference point according to the equation. According to the strength matrix transpose given by the general settings of the finite element software, add a very small number in the diagonal to avoid singularity, that is, 1.e-10. Figure 4As shown, the INP file is compiled. Specifically, the element (the first single-degree-of-freedom element) is specified as containing 2 nodes (N C and N1), a linear unit with only 2nd degree of freedom, the first reference point N1 and the control node N C Defined by the user. Then, define the second custom element (the second single-degree-of-freedom element). The definition process is the same as the above element, and the second custom element connects the control node to the outer boundary. Use a set to identify each node, and define constraint equations to ensure that the displacements of the nodes in the outer boundary are equal in the vertical direction. At this time, the global matrix will no longer be symmetrical, and concentrated forces can be specified at the control node through the "*CLOAD" statement.

[0092] Figure 5-Figure 6 and Figure 7-Figure 8 These are the calculation results of the XFEM model (i.e., the extended finite element analysis model, i.e., the improved mesh model). Cracks extend in both the two-dimensional model and the three-dimensional model under the mixed loading mode. The cracks extend along the mid-span direction and have an angle of 66.5° with the vertical direction. The two sets of numerical analysis results are consistent. In addition, the mesh models have completed the calculation without any non-convergence phenomenon, and the numerical calculation results are basically consistent with the following formula and the theoretical solution of the crack extension angle of 70° given by SIH. It can be seen that the use of the XFEM mesh model with improved convergence to calculate the crack damage problem of brittle materials can accurately reflect the structural crack damage failure mode and crack extension morphology. Therefore, the extended finite element mesh model with improved convergence can be used to numerically analyze the crack damage mechanism of tunnel linings.

[0093] ;

[0094] In the formula, K I and K II They are the stress intensity factors of type I and II cracks respectively.

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

[0096] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and equipment can refer to the corresponding processes in the aforementioned method implementation methods, and will not be repeated here.

[0097] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the modules is only a logical function division. There may be other division methods in actual implementation, such as multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or modules, which can be electrical, mechanical or other forms.

[0098] The modules described as separate components may or may not be physically separated, and the components shown as modules may or may not be physical modules, that is, they may be located in one place or distributed on multiple network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the embodiments of the present invention.

[0099] In addition, each functional module in the embodiment of the present invention may be integrated into one processing module, or each module may exist physically separately, or two or more modules may be integrated into one module.

[0100] If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention can be essentially or partly embodied in the form of a software product that contributes to the prior art. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the energy-saving signal sending / receiving method of each embodiment of the present invention. The aforementioned storage medium includes: various media that can store program codes, such as USB flash drives, mobile hard disks, ROM, RAM, magnetic disks, or optical disks.

[0101] The above description is only a preferred embodiment of the present application and an explanation of the technical principles used. Those skilled in the art should understand that the scope of the invention involved in the present application is not limited to the technical solution formed by a specific combination of the above technical features, but should also cover other technical solutions formed by any combination of the above technical features or their equivalent features without departing from the inventive concept. For example, the above features are replaced with the technical features with similar functions disclosed in this application (but not limited to) by each other.

[0102] It should be understood that the size of the serial numbers of each step in the content of the invention and the implementation methods of the present invention does not absolutely 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 implementation methods of the present invention.

Claims

1. A concrete structure cracking analysis method based on the extended finite element method, characterized in that: include: Create a mesh model to represent the concrete structure; An extended finite element mesh model is constructed for the mesh model, including drawing geometric parameter conditions of the mesh model, drawing geometric parameters of the analysis object, and the location of the preset cracks; Apply boundary conditions to the extended finite element mesh model, perform mesh division, and after the mesh division is completed, define the physical and mechanical parameters of the mesh model; Define the boundary conditions and structural load conditions of the extended finite element mesh model; Adding a first single-degree-of-freedom unit and a second single-degree-of-freedom unit that can reflect the boundary conditions of cracking parameters of the crack to the extended finite element mesh model to form an extended finite element analysis model; Performing extended finite element calculations through the extended finite element analysis model, observing the cracking process of the grid model, analyzing the stress characteristics of the model, and realizing cracking analysis of the concrete structure; The first single-degree-of-freedom unit and the second single-degree-of-freedom unit that can reflect the boundary conditions of the cracking parameters of the crack are added to the extended finite element mesh model to form an extended finite element analysis model: Based on the extended finite element mesh model, The relational expression defines a first single degree of freedom unit Element1 and a second single degree of freedom unit Element2 that can reflect the boundary conditions of cracking parameters of cracks in the extended finite element mesh model; Where f is the vertical force, u is the vertical displacement, N C is the control node, N1 is the crack tip node, y represents the vertical direction, N L N is the outer boundary L node.

2. The method for analyzing cracking of concrete structure based on extended finite element method according to claim 1, characterized in that: for and The mapping relationship between the crack tip node N1 vertical displacement and the control node external force state is established through the transpose of the stiffness matrix defined by the first single degree of freedom unit Element1, as shown in the formula: ; Small numbers close to 0 are defined in the diagonal of the stiffness matrix to avoid singularities.

3. The method for analyzing cracking of concrete structure based on extended finite element method according to claim 1, characterized in that: For the control node displacement and the outer boundary N L The mapping relationship between the vertical loads of the nodes is established by defining the stiffness matrix of the second single-degree-of-freedom unit Element2, as shown in the formula: ; By controlling the vertical force at the node and the external boundary N L The displacement of the node can define the tip of the crack, and by extending the finite element, the relationship between the nodes and the control nodes of the entire potential crack damage area can be defined.

4. The method for analyzing cracking of concrete structure based on extended finite element method according to claim 1, characterized in that: For the formula The mapping relationship can be expressed as By defining the stiffness matrix in , the average displacement of the crack propagation path is linearly related to the vertical force in the control node.

5. The concrete structure cracking analysis system based on the extended finite element method is characterized by: include: A mesh model building module builds a mesh model representing the concrete structure; An extended finite element mesh model construction module is used to construct an extended finite element mesh model for the mesh model, including drawing geometric parameter conditions of the mesh model, drawing geometric parameters of the analysis object, and the location of the preset cracks; An extended finite element mesh model parameter definition module applies boundary conditions to the extended finite element mesh model, performs mesh division, and defines physical and mechanical parameters of the mesh model after the mesh division is completed; Define the boundary conditions and structural load conditions of the extended finite element mesh model; An extended finite element analysis model forming module is configured to add a first single degree of freedom unit and a second single degree of freedom unit capable of reflecting boundary conditions of cracking parameters of cracks in the extended finite element mesh model to form an extended finite element analysis model; A concrete structure cracking analysis module performs extended finite element calculations through the extended finite element analysis model, observes the cracking process of the grid model, analyzes the stress characteristics of the model, and realizes concrete structure cracking analysis; The first single-degree-of-freedom unit and the second single-degree-of-freedom unit that can reflect the boundary conditions of the cracking parameters of the crack are added to the extended finite element mesh model to form an extended finite element analysis model: Based on the extended finite element mesh model, The relational expression defines a first single degree of freedom unit Element1 and a second single degree of freedom unit Element2 that can reflect the boundary conditions of cracking parameters of cracks in the extended finite element mesh model; Where f is the vertical force, u is the vertical displacement, N C is the control node, N1 is the crack tip node, y represents the vertical direction, N L N is the outer boundary L node.

6. An electronic device, characterized in that The method comprises a processor, a memory and a computer program stored in the memory and executable on the processor, wherein when the computer program is executed by the processor, the method for analyzing cracking of a concrete structure based on the extended finite element method as described in any one of claims 1 to 5 is implemented.

7. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for analyzing cracking of a concrete structure based on the extended finite element method according to any one of claims 1 to 5 is implemented.

Images