Method and system for calculating tensile cracking damage of concrete panel

Abstract

The invention discloses a concrete panel tensile crack damage calculation method and a concrete panel tensile crack damage calculation system, which combine extended finite element calculation with a concrete tensile crack damage calculation model, accurately simulate crack behaviors in concrete, and are particularly suitable for numerical simulation requirements under the conditions of complex crack networks and large-scale damage. According to the method, the expansion finite element method is introduced, subtle changes of crack starting, expansion paths and final forms are captured on the premise that the calculation cost is not increased, the common problem of grid dependence in a traditional finite element method is avoided, fine division or frequent updating of grids is not needed, and the calculation efficiency is improved. The method reduces the preprocessing time and iteration times, reduces the consumption of overall computing resources, can accurately describe the expansion path of the crack and the nonlinear mechanical behavior of the concrete, and is suitable for the simulation of a complex crack network.

Description

Technical Field

[0001] This invention belongs to the field of concrete tensile cracking damage calculation, and relates to a method and system for calculating concrete panel tensile cracking damage. Background Technology

[0002] Concrete, as a primary material widely used in engineering structures, exhibits cracking behavior under stress conditions, which is a key issue affecting structural safety. Traditional finite element methods (FEMs) require frequent mesh remapping when dealing with concrete cracking, resulting in high computational costs and insufficient accuracy. Concrete is a quasi-brittle material, and its fracture characteristics differ from brittle materials and general metallic materials. Linear elastic fracture mechanics and nonlinear fracture mechanics for metallic materials are not entirely applicable to concrete. While plastic damage models can effectively describe the nonlinear behavior and damage evolution process of concrete, they still have limitations in handling complex cracking problems. Therefore, there is an urgent need for a numerical simulation method that combines the extended finite element method with plastic damage models to achieve accurate and efficient simulation of concrete cracking behavior. Summary of the Invention

[0003] The purpose of this invention is to solve the problems that traditional methods in the prior art require frequent mesh re-drilling, resulting in high computational costs, and that existing linear elastic fracture mechanics and nonlinear fracture mechanics for metallic materials are not entirely applicable to concrete and cannot describe the nonlinear behavior and damage evolution process of concrete. This invention provides a method and system for calculating tensile cracking damage of concrete panels.

[0004] To achieve the above objectives, the present invention employs the following technical solution: A method for calculating tensile cracking damage of a concrete panel, comprising the following steps: Obtain the basic parameters of concrete and construct a calculation model for tensile cracking damage of concrete based on the basic parameters of concrete. Construct a simulation equation for tensile damage of concrete based on a calculation model of tensile cracking damage in concrete; Extended finite element method (EFB) calculations were performed on tensile cracking damage in concrete. Damage variables were calculated based on the concrete tensile damage simulation equations, and the damage variables were compared with the set damage variable thresholds. When the damage variable is greater than the damage variable threshold, the current element is determined to be in a cracked state. The XFEM crack propagation model is then introduced, and the XFEM crack propagation model calculates the crack morphology information of the cracked concrete based on the cracked state.

[0005] A further improvement of the present invention is that: The process of obtaining basic concrete parameters and constructing a tensile cracking damage calculation model for concrete based on these parameters includes: By introducing the Lee-Fenves model framework, adopting the uniaxial tensile-compression constitutive relation of fracture energy, and introducing a damage factor, the calculation model for tensile damage in concrete is improved.

[0006]

[0007]

[0008]

[0009] In the formula, and These represent compressive stress and compressive strain, respectively. Indicates the initial cutting modulus. and This represents the standard value of axial compressive strength and its corresponding peak strain. Indicates the compressive fracture energy. denoted by , b represents the characteristic length of the element, and b represents the ratio of equivalent plastic strain to inelastic strain; This represents the elastic stiffness under undamaged conditions.

[0010] The simulation equation for tensile damage of concrete based on the concrete tensile cracking damage calculation model includes:

[0011] in, The stress it is subjected to; For the corresponding strain; For the corresponding elastic modulus; For tensile damage variables, in finite element analysis, the damage variable values ​​of different elements at different times are in a process of change. If the damage variable of the current element... Greater than That is, the state of the previous unit is determined to be cracked.

[0012] The XFEM crack propagation model is introduced, which calculates the crack morphology information of cracked concrete based on the cracking state, including: The XFEM crack propagation model assumes a scaling factor for the external load. Loading is performed, and the system's equilibrium equations are expressed as follows:

[0013] in, Indicates the loading mode. Represents the external force vector. ; For internal forces of the system, The equivalent nodal force of the force acting on the virtual crack surface;

[0014]

[0015] In the formula, The normal vector of the crack surface; For crack opening displacement; The surface force is the force on the crack surface. The softening rate is used to reflect the relationship between the surface force on the crack and the crack opening displacement; Based on the system's equilibrium equations, the crack morphology information calculated by the XFEM crack propagation model includes crack location, displacement field of the crack region, and crack opening displacement.

[0016] The calculation of the crack location includes: The XFEM crack propagation model is based on the level set function. The location of the crack is represented by the level set function, expressed as a signed distance function:

[0017]

[0018] The equation for crack propagation and evolution is:

[0019] when When known, This represents the velocity of a point on the moving interface along the normal direction.

[0020] The calculation of the displacement field in the crack region includes: The XFEM crack propagation model describes the displacement jump and tip singularity of cracks using enrichment functions, which include standard shape functions. Jump shape function H and singular functions The displacement field in the crack region can then be described by the following equation:

[0021] In the formula, , , , These are the nodal degrees of freedom of traditional finite element, nodal degrees of freedom of crack penetration element, nodal degrees of freedom of crack tip element, and nodal degrees of freedom of crack intersection element; , , , These are, respectively, the set of nodes for conventional finite elements, the set of nodes for through elements, the set of nodes for elements with crack tips, and the set of nodes for crack intersections; , , respectively, are the node shape functions of the conventional element and the node shape function of the crack cutting element; q and nx are the number of intersecting nodes and the total number of intersecting nodes, respectively; , , These are the enrichment functions for nodes of through elements, nodes of elements with crack tips, and nodes of intersections, respectively.

[0022] The calculation of the crack opening displacement includes: The virtual crack opening displacement is described by the following formula:

[0023] In the formula, and These represent the displacement fields of the solids on either side of the crack. and For the enrichment function of the nodes of the unit that runs through both sides, , These represent the degrees of freedom at the crack penetration element nodes and the degrees of freedom at the crack tip element nodes, respectively. This is the normal vector of the crack surface.

[0024] A system for calculating tensile cracking damage of concrete panels, comprising: The tensile cracking damage model construction module is used to obtain basic concrete parameters and construct a tensile cracking damage calculation model for concrete based on these parameters. The tensile damage simulation equation construction module is used to construct the tensile damage simulation equation of concrete based on the concrete tensile cracking damage calculation model. The cracking judgment module is used to perform finite element calculations on tensile cracking damage of concrete. It calculates damage variables based on the concrete tensile damage simulation equation, compares the damage variables with the set damage variable threshold, and judges the current element to be in a cracked state when the damage variable is greater than the damage variable threshold. The crack morphology calculation module is used to import the XFEM crack propagation model. The XFEM crack propagation model calculates the crack morphology information of cracked concrete based on the cracking state.

[0025] A terminal device includes a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of any of the methods described in this invention.

[0026] A computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of any of the methods described in this invention.

[0027] Compared with the prior art, the present invention has the following beneficial effects: This invention discloses a method for calculating tensile cracking damage in concrete panels. It combines extended finite element method (EFM) calculation with a concrete tensile cracking damage calculation model to accurately simulate crack behavior in concrete. It is particularly suitable for numerical simulation needs under complex crack networks and large-scale damage conditions. By introducing the extended finite element method, it captures subtle changes in crack initiation, propagation path, and final morphology without increasing computational costs, avoiding the mesh dependency problem common in traditional finite element methods. The method of this invention does not require fine meshing or frequent updates, reducing preprocessing time and iteration count, and reducing overall computational resource consumption. It can accurately describe crack propagation paths and nonlinear mechanical behavior of concrete, and is suitable for simulating complex crack networks. Attached Figure Description

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

[0029] Figure 1 This is a flowchart of a method for calculating tensile cracking damage in concrete panels according to the present invention. Figure 2 This is a uniaxial stress-strain curve of concrete according to the present invention; Figure 3 This is an XFEM crack propagation model diagram of the present invention; Figure 4 This is a diagram of the XFEM virtual crack model of the present invention. Detailed Implementation

[0030] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. The components of the embodiments of the present invention described and shown in the accompanying drawings can generally be arranged and designed in various different configurations.

[0031] Therefore, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0032] It should be noted that similar labels and letters in the following figures indicate similar items. Therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures.

[0033] In the description of the embodiments of the present invention, it should be noted that if terms such as "upper," "lower," "horizontal," or "inner" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings, or the orientation or positional relationship commonly used when the product of the invention is in use, they are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of the present invention. Furthermore, terms such as "first" and "second" are only used to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0034] Furthermore, the use of the term "horizontal" does not imply that the component must be absolutely horizontal, but rather that it can be slightly tilted. For example, "horizontal" simply means that its direction is more horizontal than "vertical," and does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0035] In the description of the embodiments of the present invention, it should also be noted that, unless otherwise explicitly specified and limited, the terms "set," "install," "connect," and "link" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral connection; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; and they can refer to the internal connection of two components. Those skilled in the art can understand the specific meaning of the above terms in the present invention according to the specific circumstances.

[0036] The present invention will now be described in further detail with reference to the accompanying drawings: See Figure 1 This invention discloses a method for calculating tensile cracking damage in concrete panels. This concrete cracking simulation method, which couples the extended finite element method with a plastic damage model, can accurately describe the crack propagation path and the nonlinear mechanical behavior of concrete, and is suitable for simulating complex crack networks.

[0037] Specifically, the following steps are included: Step 1: Determine a suitable CDP model that ensures computational convergence. Improve the plastic damage model based on the selected model framework, and establish a calculation model for tensile cracking damage of concrete based on the plastic damage model; See Figure 2 Within the Lee-Fenves model framework, a uniaxial tensile-compression constitutive relation based on fracture energy is adopted, and a damage factor is introduced to improve the plastic damage model: Lee-Fenves' plasticity theory states that the total strain of concrete should consist of the following two parts:

[0038] in, For the strain of the elastic component, For the strain of the plastic portion, the stress relationship can be expressed as:

[0039] In the formula, denoted as concrete stress; D represents the degradation damage variable due to decreased elastic stiffness (0-1 indicates no damage to complete damage); Effective stress; The elastic stiffness under undamaged conditions can be expressed as:

[0040]

[0041]

[0042]

[0043]

[0044] In the formula, and For compressive stress and compressive strain, The initial cutting modulus, and The standard value of axial compressive strength and its corresponding peak strain are given. The compressive fracture energy, denoted as the element characteristic length, and b as the ratio of equivalent plastic strain to inelastic strain; the uniaxial compressive stress-strain relationship of concrete can be divided into the quasi-elastic stage, the strengthening stage, and the softening stage.

[0045] Step 2: Based on the determined calculation model, establish the simulation equation for tensile damage in concrete; The simulation equation for tensile damage in concrete is described as follows:

[0046] In the formula, The stress it is subjected to; For the corresponding strain; For the corresponding elastic modulus; Tensile damage variable Step 3: Determine the relationship between tensile damage and the strength value of concrete cracking to determine whether cracking has occurred; A finite element calculation model of ABAQUS concrete structure was established based on actual engineering projects. The concrete finite element model is refined by extending the finite element mesh model, including drawing the geometric parameters of the concrete structure being analyzed, selecting appropriate geometric parameter conditions for the model, determining the location of cracks, and refining the mesh around the cracks. Mesh the established extended finite element model; After the mesh generation is completed, the mechanical parameters of the concrete model are defined, including elastic parameters and plastic parameters, etc. Define the boundary conditions, constraints, and structural load conditions for the extended finite element model; The damage evolution equation of the concrete plastic damage constitutive model is defined in the calculation. The calculation method of tensile cumulative damage of concrete panel is written into ABAQUS software through a subroutine. The three-dimensional model of concrete panel is built through ABAQUS software so that it can output the tensile cumulative damage variables of concrete panel.

[0047] S4, in XFEM calculations, for cracked concrete, calculates the crack location, crack area displacement field, and crack opening displacement.

[0048] Using the XFEM extended finite element function in ABAQUS software, the cumulative tensile damage of a concrete panel was calculated. The damage variable values ​​of different elements at different times are in a changing process. Greater than When the concrete cracks, XFEM begins to calculate the crack parameters.

[0049] XFEM builds upon the traditional finite element method by adding enhancement functions for the crack and crack tip to describe crack discontinuities. See also Figure 3 This represents the propagation of any crack within an element in the XFEM model (where the nodes of the elements through which the crack propagates are labeled with hollow squares, and the element where the crack tip is located is represented by a hollow circle).

[0050] In specific calculations: In ABAQUS, XFEM scales the assumed external loads using a scaling factor. Loading is performed, external force vector Represented as , Indicating the loading mode, the equilibrium equations calculated by the XFEM system can be expressed as:

[0051] In the formula, For internal forces of the system, The equivalent nodal force is the force acting on the virtual crack surface.

[0052]

[0053]

[0054] In the formula, The normal vector of the crack surface; For crack opening displacement; The surface force is the force on the crack surface. The softening rate is used to reflect the relationship between the surface force on the crack and the crack opening displacement.

[0055] Furthermore, describe the location of the crack: XFEM uses the Level Set Method to represent the geometry of cracks, where level set functions are used. To indicate the location of a crack, the horizontal set function is usually represented by the signed distance function:

[0056]

[0057] Furthermore, the crack propagation evolution equation is:

[0058] when When known, This represents the velocity of a point on the moving interface along the normal direction, and is used to describe the evolution of the crack.

[0059] Furthermore, the displacement field in the crack region is described: XFEM describes the displacement jump and tip singularity of cracks using enrichment functions, which include standard shape functions. Jump shape function H and singular functions In XFEM, these functions are combined to describe the displacement field in the crack region:

[0060] In the formula, , , , These are the nodal degrees of freedom of traditional finite element, nodal degrees of freedom of crack penetration element, nodal degrees of freedom of crack tip element, and nodal degrees of freedom of crack intersection element; , , , These are, respectively, the set of nodes for conventional finite elements, the set of nodes for through elements, the set of nodes for elements with crack tips, and the set of nodes for crack intersections; , These are the nodal shape functions of a conventional element and the nodal shape functions of a crack-cutting element, respectively. q , nx These are the number of cross nodes and the total number of cross nodes, respectively. , , These are the enrichment functions for nodes of through elements, nodes of elements with crack tips, and nodes of intersections, respectively.

[0061] Furthermore, the virtual crack opening displacement is described: See Figure 4 Virtual crack opening displacement It can be represented as:

[0062] In the formula, and These represent the displacement fields of the solids on either side of the crack. and For the enrichment function of the nodes of the unit that runs through both sides, , These represent the degrees of freedom at the crack penetration element nodes and the degrees of freedom at the crack tip element nodes, respectively. This is the normal vector of the crack surface.

[0063] The plastic damage calculation method for this concrete panel can calculate the tensile cracking damage of the concrete panel and determine the crack opening. Through XFEM calculation, the crack location, crack area displacement field, and crack opening displacement of the concrete panel rockfill dam panel can be obtained.

[0064] This invention discloses a system for calculating tensile cracking damage of concrete panels, comprising: The tensile cracking damage model construction module is used to obtain basic concrete parameters and construct a tensile cracking damage calculation model for concrete based on these parameters. The tensile damage simulation equation construction module is used to construct the tensile damage simulation equation of concrete based on the concrete tensile cracking damage calculation model. The cracking judgment module is used to perform finite element calculations on tensile cracking damage of concrete. It calculates damage variables based on the concrete tensile damage simulation equation, compares the damage variables with the set damage variable threshold, and judges the current element to be in a cracked state when the damage variable is greater than the damage variable threshold. The crack morphology calculation module is used to import the XFEM crack propagation model. The XFEM crack propagation model calculates the crack morphology information of cracked concrete based on the cracking state.

[0065] This invention creatively combines the Extended Finite Element Method (XFEM) with a plastic damage model, significantly improving the efficiency and accuracy of concrete crack simulation. This innovative technology can not only efficiently and accurately simulate crack behavior in concrete, but is also particularly suitable for numerical simulation needs under complex crack networks and large-scale damage conditions, providing a solid theoretical foundation and technical guarantee for the design optimization and safety assessment of concrete structures.

[0066] By introducing the extended finite element method, this invention can capture subtle changes in crack initiation, propagation path, and final morphology without increasing computational costs, avoiding the mesh dependency problem common in traditional finite element methods, thereby greatly improving the accuracy and reliability of simulation results.

[0067] The method of this invention eliminates the need for fine mesh subdivision or frequent updates, reducing preprocessing time and iteration count, and decreasing overall computational resource consumption. Furthermore, it supports adaptive mesh refinement and other advanced algorithms, further enhancing the flexibility and efficiency of the simulation.

[0068] A schematic diagram of a terminal device according to an embodiment of the present invention. The terminal device of this embodiment includes: a processor, a memory, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the steps in the various method embodiments described above. Alternatively, when the processor executes the computer program, it implements the functions of each module / unit in the various device embodiments described above.

[0069] The computer program can be divided into one or more modules / units, which are stored in the memory and executed by the processor to complete the present invention.

[0070] The terminal device can be a desktop computer, laptop computer, cloud server, or other device with strong computing power. The terminal device may include, but is not limited to, a processor and memory.

[0071] The optimal choice for the processor is a multi-core high-speed central processing unit (CPU).

[0072] The memory can be used to store the computer program and / or module. The processor implements various functions of the terminal device by running or executing the computer program and / or module stored in the memory and calling the data stored in the memory.

[0073] If the modules / units integrated into the terminal device are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, all or part of the processes in the methods of the above embodiments can also be implemented by a computer program instructing related hardware. The computer program can be stored in a computer-readable storage medium, and when executed by a processor, it can implement the steps of the various method embodiments described above. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or certain intermediate forms. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording media, USB flash drives, portable hard drives, magnetic disks, optical disks, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signals, telecommunication signals, and software distribution media, etc. It should be noted that the content included in the computer-readable medium can be appropriately added or removed according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, computer-readable media do not include electrical carrier signals and telecommunication signals.

[0074] The above are merely preferred embodiments of the present invention and are not intended to limit the present invention. Various modifications and variations can be made to the present invention by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for calculating tensile cracking damage in concrete panels, characterized in that, Includes the following steps: Obtain the basic parameters of concrete and construct a calculation model for tensile cracking damage of concrete based on the basic parameters of concrete. Construct a simulation equation for tensile damage of concrete based on a calculation model of tensile cracking damage in concrete; Extended finite element method (EFB) calculations were performed on tensile cracking damage in concrete. Damage variables were calculated based on the concrete tensile damage simulation equations, and the damage variables were compared with the set damage variable thresholds. When the damage variable is greater than the damage variable threshold, the current element is determined to be in a cracked state. The XFEM crack propagation model is then introduced, and the XFEM crack propagation model calculates the crack morphology information of the cracked concrete based on the cracked state.

2. The method for calculating tensile cracking damage of a concrete panel according to claim 1, characterized in that, The process of obtaining basic concrete parameters and constructing a tensile cracking damage calculation model for concrete based on these parameters includes: By introducing the Lee-Fenves model framework, adopting the uniaxial tensile-compression constitutive relation of fracture energy, and introducing a damage factor, the calculation model for tensile damage in concrete is improved. In the formula, and These represent compressive stress and compressive strain, respectively. Indicates the initial cutting modulus. and This represents the standard value of axial compressive strength and its corresponding peak strain. Indicates the compressive fracture energy. denoted by , b represents the characteristic length of the element, and b represents the ratio of equivalent plastic strain to inelastic strain; This represents the elastic stiffness under undamaged conditions.

3. The method for calculating tensile cracking damage of a concrete panel according to claim 1, characterized in that, The simulation equation for tensile damage of concrete based on the concrete tensile cracking damage calculation model includes: in, The stress it is subjected to; For the corresponding strain; For the corresponding elastic modulus; For tensile damage variables, in finite element analysis, the damage variable values ​​of different elements at different times are in a process of change. If the damage variable of the current element... Greater than That is, the state of the previous unit is determined to be cracked.

4. The method for calculating tensile cracking damage of a concrete panel according to claim 1, characterized in that, The XFEM crack propagation model is introduced, which calculates the crack morphology information of cracked concrete based on the cracking state, including: The XFEM crack propagation model assumes a scaling factor for the external load. Loading is performed, and the system's equilibrium equations are expressed as follows: in, Indicates the loading mode. Represents the external force vector. ; For internal forces of the system, The equivalent nodal force of the force acting on the virtual crack surface; In the formula, The normal vector of the crack surface; For crack opening displacement; For the surface force of the crack; The softening rate is used to reflect the relationship between the surface force on the crack and the crack opening displacement; Based on the system's equilibrium equations, the crack morphology information calculated by the XFEM crack propagation model includes crack location, displacement field of the crack region, and crack opening displacement.

5. The method for calculating tensile cracking damage of a concrete panel according to claim 4, characterized in that, The calculation of the crack location includes: The XFEM crack propagation model is based on the level set function. The location of the crack is represented by the level set function, expressed as a signed distance function: The equation for crack propagation and evolution is: when When known, This represents the velocity of a point on the moving interface along the normal direction.

6. The method for calculating tensile cracking damage of a concrete panel according to claim 4, characterized in that, The calculation of the displacement field in the crack region includes: The XFEM crack propagation model describes the displacement jump and tip singularity of cracks using enrichment functions, which include standard shape functions. Jump shape function H and singular functions The displacement field in the crack region can then be described by the following equation: In the formula, , , , These are the nodal degrees of freedom of traditional finite element, nodal degrees of freedom of crack penetration element, nodal degrees of freedom of crack tip element, and nodal degrees of freedom of crack intersection element; , , , These are, respectively, the set of nodes for conventional finite elements, the set of nodes for through elements, the set of nodes for elements with crack tips, and the set of nodes for crack intersections; , , respectively, are the node shape functions of the conventional element and the node shape function of the crack cutting element; q and nx are the number of intersecting nodes and the total number of intersecting nodes, respectively; , , These are the enrichment functions for nodes of through elements, nodes of elements with crack tips, and nodes of intersections, respectively.

7. The method for calculating tensile cracking damage of a concrete panel according to claim 4, characterized in that, The calculation of the crack opening displacement includes: The virtual crack opening displacement is described by the following formula: In the formula, and These represent the displacement fields of the solids on either side of the crack. and For the enrichment function of the nodes of the unit that runs through both sides, , These represent the degrees of freedom at the crack penetration element nodes and the degrees of freedom at the crack tip element nodes, respectively. This is the normal vector of the crack surface.

8. A system for calculating tensile cracking damage of concrete panels, characterized in that, include: The tensile cracking damage model construction module is used to obtain basic concrete parameters and construct a tensile cracking damage calculation model for concrete based on these parameters. The tensile damage simulation equation construction module is used to construct the tensile damage simulation equation of concrete based on the concrete tensile cracking damage calculation model. The cracking judgment module is used to perform finite element calculations on tensile cracking damage of concrete. It calculates damage variables based on the concrete tensile damage simulation equation, compares the damage variables with the set damage variable threshold, and judges the current element to be in a cracked state when the damage variable is greater than the damage variable threshold. The crack morphology calculation module is used to import the XFEM crack propagation model. The XFEM crack propagation model calculates the crack morphology information of cracked concrete based on the cracking state.

9. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method as described in any one of claims 1-7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method as described in any one of claims 1-7.

Images