Macro-micro dual-scale simulation method for concrete structure damage

Abstract

The invention relates to a macro-micro dual-scale simulation method for concrete structure damage. The method is suitable for the technical field of concrete structure multi-scale analysis. According to the technical scheme, the method comprises the steps that macroscale finite element mesh generation is conducted on a concrete structure, and a macroscopic mesh containing all structural analysis areas is obtained; assuming that the analysis scale of each macroscopic unit integral point in the macroscopic grid in the initial state is a macroscopic scale and is classified as a macroscopic scale analysis area; the concrete loading process is divided into N increment steps, and finite element equilibrium iteration is carried out on any increment step; and in any iteration step of any increment step, judging whether each macroscopic unit integration point in the macroscopic scale analysis area meets a preset analysis scale adaptive transformation criterion, and if so, transforming the analysis scale of the corresponding macroscopic unit integration point from a macroscopic scale to a macro-micro dual scale, and classifying the macro-micro dual scale into a macro-micro dual scale analysis area.

Description

Technical Field

[0001] The invention relates to a macroscopic and microscopic dual-scale simulation method for concrete structure damage, which is applicable to the field of multi-scale analysis of concrete structures. Background Art

[0002] Concrete structures are widely used in engineering fields such as water conservancy, civil engineering, and transportation, and are an important part of the national infrastructure system. In complex environments, concrete structures are prone to damage and cracking under various loads, which seriously affects their service performance and may even lead to failure. Therefore, accurately simulating the damage evolution process of concrete structures is of great significance for reasonably evaluating their working properties and ensuring their safety.

[0003] At present, the simulation of the damage evolution process of concrete structures is carried out on a single macroscopic structural scale, that is, a damage constitutive model is established through mathematical abstraction based on the analysis and summary of the macroscopic mechanical properties of concrete obtained from the test, and the macroscopic phenomenological simulation of the damage evolution process of concrete structures is carried out on this basis. The highly random microscopic structure of concrete materials is not considered in the simulation process. The above-mentioned prior art has the following deficiencies:

[0004] 1. Concrete damage and cracking show obvious cross-scale characteristics, and its evolution process is directly related to the microscopic material structure. However, in the existing technology, the simulation of the damage evolution process of concrete structure is only carried out on a single macroscopic scale, which is unable to accurately simulate the damage evolution process of concrete structure with macro-micro cross-scale characteristics.

[0005] 2. In the prior art, it is assumed that the concrete material is "uniform", and the complex microstructure of the concrete material cannot be considered.

[0006] 3. In the existing technology, the degree of damage of concrete in different stress stages can only be quantitatively described on a macro scale, but it is impossible to quantitatively describe the degree of damage of concrete on both macro and micro scales at the same time. Summary of the invention

[0007] The technical problem to be solved by the present invention is: in view of the above-mentioned problems, a macro- and micro-scale dual-scale simulation method for concrete structure damage is provided.

[0008] The technical solution adopted by the present invention is: a macro-micro dual-scale simulation method for concrete structure damage, characterized by comprising:

[0009] Perform macro-scale finite element meshing on the concrete structure to obtain a macro-grid that includes the entire structural analysis area;

[0010] It is assumed that the analysis scale of each macro-unit integration point in the macro-grid at the initial state is macro-scale and is classified as the macro-scale analysis area;

[0011] The concrete loading process is divided into N incremental steps, and finite element equilibrium iteration is carried out for any incremental step;

[0012] In any iteration step of any incremental step, it is determined whether each macro unit integration point in the macro scale analysis area satisfies the preset analysis scale adaptive transformation criterion;

[0013] If it is satisfied, the analysis scale of the corresponding macro unit integration point is transformed from macro scale to macro-micro dual scale and classified as macro-micro dual scale analysis area;

[0014] The finite element analysis model is used to simulate the damage evolution process of the concrete structure, in which the integral points in the macro-scale analysis area adopt the macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with the meso-scale finite element analysis model.

[0015] The judgment basis of the analysis scale adaptive transformation criterion includes the stress of the macro-unit integration point in the macro-scale analysis area.

[0016] The analysis scale adaptive transformation criteria include:

[0017] σ1-[σ]≥0

[0018] Where σ1 is the stress at the macro unit integration point located in the macro scale analysis area. The first principal stress, with tension as positive; [σ] is the maximum allowable tensile stress.

[0019] The macro-scale finite element analysis model includes an elastic constitutive relation.

[0020] The micro-scale finite element analysis model is adaptively generated according to a given side length, concrete aggregate content and particle size range, and includes aggregate units, mortar units and interface transition zone units, wherein the aggregate adopts an elastic constitutive relationship, and the mortar and interface transition zone adopt an elastic-plastic damage constitutive relationship.

[0021] The given side length dimension is determined based on a coupled volume method, which refers to taking the volume or area of ​​the microscopic model associated with a macroscopic unit integration point as the integrated volume or area of ​​the integration point.

[0022] The given side length dimensions include

[0023] In the three-dimensional case, the shape of the mesoscopic model is a cube, and in the two-dimensional case, the shape of the mesoscopic model is a square, and the side length of the mesoscopic model is l m The calculation formula is as follows:

[0024]

[0025] In the formula, S IP and V IP They represent the integrated area and integrated volume of the macro unit integration points respectively.

[0026] The micro-scale finite element analysis model includes:

[0027] The macro-element integration point strain is converted into the displacement boundary conditions of the associated micro-scale finite element analysis model.

[0028] The method of converting the macro unit integration point strain into the displacement boundary condition of the micro-scale finite element analysis model associated therewith includes:

[0029]

[0030] In the formula, is the macro unit integration point strain; is the linear displacement of the boundary of the meso-scale finite element analysis model; x is the coordinate of the boundary node of the meso-scale finite element analysis model; Fluctuating displacements of the model boundary for mesoscale finite element analysis.

[0031] The micro-scale finite element analysis model includes:

[0032] The stress field, stiffness matrix and damage variable field of the micro-scale finite element analysis model are homogenized into the associated macro-element integration point stress, constitutive matrix and damage variable.

[0033] The stress field of the micro-scale finite element analysis model is homogenized into the associated macro-element integration point stress, including:

[0034]

[0035] In the formula, is the stress of the macro-unit integration point located in the macro-micro dual-scale analysis area; σ is the stress of the integration point in the micro-scale finite element analysis model associated with the macro-unit integration point; Ω is the integration range of the micro-scale model associated with the macro-unit integration point.

[0036] The stiffness matrix of the micro-scale finite element analysis model is homogenized into the associated macro-element integration point constitutive matrix, including:

[0037]

[0038] Where D M is the constitutive matrix of the macro-element integration point located in the macro-micro dual-scale analysis area; K is the vertex coordinate of the micro-scale finite element analysis model associated with the macro-element integration point, i, j = 1, 2; M =KAA -K AN K NN - 1 K NA , K NN , K NA , K AN , K AA These are the four sub-matrices of the stiffness matrix of the meso-scale finite element analysis model, which are composed of stiffness matrix elements related to non-vertex node degrees of freedom, non-vertex node degrees of freedom and vertex structural degrees of freedom, vertex node degrees of freedom and non-vertex node degrees of freedom, and vertex node degrees of freedom.

[0039] The damage variable field of the micro-scale finite element analysis model is homogenized into the macro-element integration point damage variables associated with it, including:

[0040]

[0041] Where, d M is the damage variable of the macro unit integration point; d i is the damage variable of the i-th integration point in the damage loading area in the mesoscale finite element analysis model; Ω d,i is the integration range of the i-th integration point in the damage loading region in the mesoscale finite element analysis model; Ω d is the extent of the damaged loading zone in the meso-scale finite element analysis model.

[0042] A macro-micro dual-scale simulation device for concrete structure damage, comprising:

[0043] The meshing module is used to perform macro-scale finite element meshing on the concrete structure to obtain a macro-grid that includes the entire structural analysis area;

[0044] The grid initialization module is used to assume that the analysis scale of each macro unit integration point in the macro grid is macro scale at the initial state and is classified as a macro scale analysis area;

[0045] The regional update module is used to divide the concrete loading process into N incremental steps. For any incremental step, finite element equilibrium iteration is carried out; in any iteration step of any incremental step, it is judged whether each macro-unit integration point in the macro-scale analysis area meets the preset analysis scale adaptive transformation criterion; if it meets the criterion, the analysis scale of the corresponding macro-unit integration point is transformed from the macro scale to the macro-micro dual scale, and it is classified as the macro-micro dual scale analysis area;

[0046] The partition simulation module is used to simulate the damage evolution process of concrete structure using the finite element analysis model, in which the integral points in the macro-scale analysis area adopt the macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with the micro-scale finite element analysis model.

[0047] A storage medium stores a computer program executable by a processor, wherein the steps of the method are implemented when the computer program is executed.

[0048] A data processing device comprises a memory and a processor. The memory stores a computer program executable by the processor. When the computer program is executed, the steps of the method are implemented.

[0049] The beneficial effect of the present invention is that: the present invention determines whether each macro-unit integral point in the macro-scale analysis area meets the preset analysis scale adaptive transformation criterion in any iterative step of any incremental step, so as to dynamically adjust the macro-scale analysis area and the macro-micro dual-scale analysis area during the simulation process to reduce the calculation scale, breaking through the bottleneck that the traditional macro-scale analysis cannot accurately simulate the damage evolution process of the concrete structure with macro-micro cross-scale characteristics, and has a high calculation and analysis efficiency, providing an effective solution for efficiently and accurately simulating the damage evolution process of the concrete structure.

[0050] The present invention performs partitioning based on the adaptive transformation criterion of analysis scale. The integral points in the macro-scale analysis area adopt the macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with the meso-scale finite element analysis model. This can not only consider the influence of the mesostructure of the concrete material in the process of simulating the damage evolution of the concrete structure, but also quantitatively simulate the concrete damage evolution process at the macro-micro dual-scale. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 Schematic diagram of a flow chart of an embodiment of the present invention.

[0052] Figure 2 Schematic diagram of the structural dimensions, loading and boundary conditions of an embodiment of the present invention.

[0053] Figure 3 This is a macro grid diagram of the entire structural analysis area according to an embodiment of the present invention.

[0054] Figure 4 An example diagram of a microscopic model associated with a macroscopic unit integration point that is adaptively generated according to an embodiment of the present invention.

[0055] Figure 5 The following is a diagram of a scale adaptive transformation process according to an embodiment of the present invention.

[0056] Figure 6This is an example diagram of the damage evolution process at a microscopic scale according to an embodiment of the present invention.

[0057] Figure 7 This is a diagram of the macro-scale damage evolution process of an embodiment of the present invention. DETAILED DESCRIPTION

[0058] Example 1: Figure 1 As shown, this embodiment is a macro-micro dual-scale simulation method for concrete structure damage, comprising the following steps:

[0059] S100, performing macro-scale finite element meshing on the concrete structure to obtain a macro-mesh that includes the entire structural analysis area.

[0060] Figure 2 In order to simulate the tensile damage evolution process of the concrete dog-bone specimen in this embodiment, the specimen is 480 mm high and 320 mm wide. A displacement constraint is applied to the bottom of the specimen and a vertical displacement load is applied to the top with a value of 0.15 mm.

[0061] In this embodiment, a four-node reduced integration element is used to perform macro-scale meshing of the concrete dog-bone specimen, and the number of elements is 95. Figure 3 It is a macroscopic mesh obtained by finite element meshing that includes the entire structural analysis area.

[0062] S200, assuming that the analysis scale of each macro unit integration point in the macro grid in the initial state is macro scale, and all are classified as macro scale analysis area.

[0063] In this embodiment, elastic material parameters are assigned to all macro units in the macro scale analysis area, and the macro scale finite element analysis model is established using elastic constitutive relations. The elastic parameters are the elastic modulus E and Poisson's ratio μ of the concrete material, which are taken as 33.19 GPa and 0.2 respectively in this embodiment.

[0064] S300, dynamically updating the macro-scale analysis region and the macro-micro dual-scale analysis region on the concrete structure based on the analysis scale adaptive transformation criterion.

[0065] S310, dividing the concrete loading process into N incremental steps, and performing macro-scale finite element equilibrium iteration for any incremental step;

[0066] S320. In any iteration step of any incremental step, determine whether each macro-unit integration point in the macro-scale analysis area satisfies a preset analysis scale adaptive transformation criterion.

[0067] The analysis scale adaptive transformation criterion in this embodiment refers to the judgment criterion used when the analysis scale of the macro unit integration point is transformed from the macro scale to the macro-micro dual scale. In this example, the judgment is made based on the stress of each macro unit integration point in the macro scale analysis area, and the expression is as follows:

[0068]

[0069] In the formula, is the stress at the macro-element integration point located in the macro-scale analysis area;

[0070] In this embodiment, a specific expression of formula (1) is given based on the maximum tensile stress criterion:

[0071] σ1-[σ]≥0(2)

[0072] Where σ1 is the stress at the macro unit integration point located in the macro scale analysis area. The first principal stress (with tension as positive); [σ] is the maximum allowable tensile stress.

[0073] S330. If the macro-unit integration point in the macro-scale analysis area meets the analysis scale adaptive transformation criterion, the analysis scale of the corresponding macro-unit integration point is transformed from the macro scale to the macro-micro dual scale and is classified as the macro-micro dual scale analysis area; if not, the analysis scale of the corresponding macro-unit integration point remains the macro scale.

[0074] In this embodiment, the macro-scale analysis area refers to the area in the structural analysis area where the analysis scale has not been transformed from the initial macro-scale to the macro-micro dual scale, and the macro-micro dual scale analysis area refers to the area in the structural analysis area where the analysis scale is transformed from the macro-scale to the macro-micro dual scale. The macro-scale analysis area and the macro-micro dual scale analysis area do not overlap and together constitute the entire analysis area of ​​the concrete structure.

[0075] S400, using finite element analysis model to simulate the damage evolution process of concrete structure, where the integral points in the macro-scale analysis area adopt the macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with the micro-scale finite element analysis model.

[0076] In this embodiment, the analysis scale of any integral point of any unit in the macro-scale analysis area is macro-scale, and an elastic constitutive relationship is adopted to establish a macro-scale finite element analysis model.

[0077] In this example, the analysis scale of any integration point of any unit in the macro-micro dual-scale analysis area is macro-micro dual-scale, there is a micro-scale finite element analysis model associated with it, and the constitutive relation of the integration point is determined by the macro-microscale connection.

[0078] In this embodiment, the meso-scale finite element analysis model adopts a rectangular coordinate system with the origin placed at the geometric center of the model. The model is adaptively generated based on a given side length, concrete aggregate content, and particle size range, and has a random meso-structure, including aggregate units, mortar units, and interface transition zone units. The elastic constitutive relationship is adopted for the aggregate, and the elastoplastic damage constitutive relationship is adopted for the mortar and the interface transition zone (ITZ).

[0079] In this embodiment, the aggregate content of the concrete is 50%, and the aggregate particle size ranges from 5 to 8 mm. Table 1 shows the constitutive parameters of the aggregate, mortar and interface transition zone.

[0080] Table 1 Constitutive parameters of aggregate, mortar and interface transition zone (ITZ)

[0081]

[0082] In this example, the fixed side length of the mesoscale finite element analysis model is determined based on the coupled volume method. The coupled volume refers to taking the area of ​​the mesoscale model associated with the macro unit integration point as the integral area of ​​the integration point. The shape of the mesoscale model is a square, and the side length of the mesoscale model is l m The calculation formula is as follows:

[0083]

[0084] In the formula, S IP Represents the integral area of ​​the macro unit integration point. Figure 4 Example of adaptively generated edge lengths l associated with two macroelement integration points m Two mesoscopic models with diameters of 38.9 mm and 41.2 mm respectively.

[0085] In this embodiment, the macro-micro scale connection includes two links: localization and homogenization. Localization refers to converting the strain of the macro unit integration point into the displacement boundary condition of the associated micro model, and homogenization refers to homogenizing the stress field, stiffness matrix and damage variable field of the micro model into the stress, constitutive matrix and damage variable of the macro unit integration point associated with it.

[0086] In this example, the transformation of the macro unit integration point strain into the displacement boundary condition of the associated micro model refers to the transformation of the macro unit integration point strain into the periodic displacement boundary condition of the associated micro model. The specific transformation formula is:

[0087]

[0088] In the formula, is the macro unit integration point strain; is the linear displacement of the boundary of the mesoscopic model; x is the coordinate of the node at the boundary of the mesoscopic model; is the fluctuating displacement of the boundary of the mesoscopic model.

[0089] In this example, the specific formula for homogenizing the stress field of the microscopic model into the stress of the associated macroscopic unit integration point is:

[0090]

[0091] In the formula, is the stress of the macro-unit integration point located in the macro-micro dual-scale analysis area; σ is the stress of the integration point in the micro-model associated with the macro-unit integration point; Ω is the integration range of the micro-model associated with the macro-unit integration point.

[0092] In this embodiment, the specific formula for homogenizing the stiffness matrix of the microscopic model into the associated macroscopic unit integration point constitutive matrix is:

[0093]

[0094] Where D M is the constitutive matrix of the macro-element integration point located in the macro-micro dual-scale analysis area; is the vertex coordinate of the microscopic model associated with the macro unit integration point, i,j = 1,2; K M =K AA -K AN K NN -1 K NA , K NN , K NA , K AN , K AA are the four sub-matrices of the microscopic model stiffness matrix, which are composed of stiffness matrix elements related to non-vertex node degrees of freedom, non-vertex node degrees of freedom and vertex structural degrees of freedom, vertex node degrees of freedom and non-vertex node degrees of freedom, and vertex node degrees of freedom.

[0095] In this embodiment, the specific formula for homogenizing the damage variable field of the microscopic model into the damage variable of the macroscopic unit integration point associated therewith is:

[0096]

[0097] Where, d M is the damage variable of the macro unit integration point; d i is the damage variable of the i-th integration point in the damage loading area in the mesoscopic model; Ω d,i is the integration range of the i-th integration point in the damage loading area in the mesoscopic model; Ω d is the range of the damage loading zone in the mesoscopic model.

[0098] In this embodiment, Figure 5The adaptive transformation process of analysis scale is given; Figure 6 The mesoscale damage evolution of the mesomodel associated with three macrocell integration points is illustrated; Figure 7 The macroscopic damage evolution process is given.

[0099] Embodiment 2: This embodiment is a macro-micro dual-scale simulation device for concrete structure damage, including: a grid generation module, a grid initialization module, a regional update module and a partition simulation module.

[0100] In this example, the mesh generation module is used to perform macro-scale finite element mesh generation on the concrete structure to obtain a macro-scale mesh that includes the entire structural analysis area; the mesh initialization module is used to assume that the analysis scale of each macro-unit integration point in the macro-grid is macro-scale in the initial state and is classified as a macro-scale analysis area.

[0101] In this embodiment, the regional update module is used to divide the concrete loading process into N incremental steps. For any incremental step, finite element equilibrium iteration is carried out; in any iteration step of any incremental step, it is judged whether each macro-unit integration point in the macro-scale analysis area meets the preset analysis scale adaptive transformation criterion; if it meets the criterion, the analysis scale of the corresponding macro-unit integration point is transformed from the macro scale to the macro-micro dual scale, and is classified as the macro-micro dual scale analysis area.

[0102] In this embodiment, the partition simulation module is used to simulate the damage evolution process of the concrete structure using a finite element analysis model, wherein the integral points in the macro-scale analysis area use a macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with a micro-scale finite element analysis model.

[0103] Embodiment 3: This embodiment is a storage medium on which a computer program that can be executed by a processor is stored. When the computer program is executed, the steps of the macro-micro dual-scale simulation method for concrete structure damage in Embodiment 1 are implemented.

[0104] Embodiment 4: This embodiment is a data processing device having a memory and a processor. The memory stores a computer program that can be executed by the processor. When the computer program is executed, the steps of the macro-micro dual-scale simulation method of concrete structure damage in Embodiment 1 are implemented.

Claims

1. A macro-micro dual-scale simulation method for concrete structure damage, characterized in that: include: Perform macro-scale finite element meshing on the concrete structure to obtain a macro-grid that includes the entire structural analysis area; It is assumed that the analysis scale of each macro-unit integration point in the macro-grid at the initial state is macro-scale and is classified as the macro-scale analysis area; The concrete loading process is divided into N incremental steps, and finite element equilibrium iteration is carried out for any incremental step; In any iteration step of any incremental step, it is determined whether each macro unit integration point in the macro scale analysis area satisfies the preset analysis scale adaptive transformation criterion; If it is satisfied, the analysis scale of the corresponding macro unit integration point is transformed from macro scale to macro-micro dual scale and classified as macro-micro dual scale analysis area; The finite element analysis model is used to simulate the damage evolution process of the concrete structure, in which the integral points in the macro-scale analysis area adopt the macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with the meso-scale finite element analysis model.

2. The macro-micro dual-scale simulation method for concrete structure damage according to claim 1 is characterized in that: The judgment basis of the analysis scale adaptive transformation criterion includes the stress of the macro-unit integration point in the macro-scale analysis area.

3. The macro-micro dual-scale simulation method for concrete structure damage according to claim 2 is characterized in that: The analysis scale adaptive transformation criteria include: σ1-[σ]≥0 Where σ1 is the stress at the macro unit integration point located in the macro scale analysis area. The first principal stress, with tension as positive; [σ] is the maximum allowable tensile stress.

4. The macro-micro dual-scale simulation method for concrete structure damage according to claim 1 is characterized in that: The macro-scale finite element analysis model includes an elastic constitutive relation.

5. The macro-micro dual-scale simulation method for concrete structure damage according to claim 1 is characterized in that: The micro-scale finite element analysis model is adaptively generated according to a given side length, concrete aggregate content and particle size range, and includes aggregate units, mortar units and interface transition zone units, wherein the aggregate adopts an elastic constitutive relationship, and the mortar and interface transition zone adopt an elastic-plastic damage constitutive relationship.

6. The macro-micro dual-scale simulation method for concrete structure damage according to claim 5 is characterized in that: The given side length dimension is determined based on a coupled volume method, which refers to taking the volume or area of ​​the microscopic model associated with a macroscopic unit integration point as the integrated volume or area of ​​the integration point.

7. The macro-micro dual-scale simulation method for concrete structure damage according to claim 6 is characterized in that: The given side length dimensions include In the three-dimensional case, the shape of the mesoscopic model is a cube, and in the two-dimensional case, the shape of the mesoscopic model is a square, and the side length of the mesoscopic model is l m The calculation formula is as follows: In the formula, S IP and V IP They represent the integrated area and integrated volume of the macro unit integration points respectively.

8. The macro-micro dual-scale simulation method for concrete structure damage according to claim 5 is characterized in that: The micro-scale finite element analysis model includes: The macro-element integration point strain is converted into the displacement boundary conditions of the associated micro-scale finite element analysis model.

9. The macro-micro dual-scale simulation method for concrete structure damage according to claim 8 is characterized in that: The method of converting the macro unit integration point strain into the displacement boundary condition of the micro-scale finite element analysis model associated therewith includes: In the formula, is the macro unit integration point strain; is the linear displacement of the boundary of the meso-scale finite element analysis model; x is the coordinate of the boundary node of the meso-scale finite element analysis model; Fluctuating displacements of the model boundary for mesoscale finite element analysis.

10. The macro-micro dual-scale simulation method for concrete structure damage according to claim 5, characterized in that: The micro-scale finite element analysis model includes: The stress field, stiffness matrix and damage variable field of the micro-scale finite element analysis model are homogenized into the associated macro-element integration point stress, constitutive matrix and damage variable.

11. The macro-micro dual-scale simulation method for concrete structure damage according to claim 10, characterized in that: The stress field of the micro-scale finite element analysis model is homogenized into the associated macro-element integration point stress, including: In the formula, is the stress of the macro-unit integration point located in the macro-micro dual-scale analysis area; σ is the stress of the integration point in the micro-scale finite element analysis model associated with the macro-unit integration point; Ω is the integration range of the micro-scale model associated with the macro-unit integration point.

12. The macro-micro dual-scale simulation method for concrete structure damage according to claim 10, characterized in that: The stiffness matrix of the micro-scale finite element analysis model is homogenized into the associated macro-element integration point constitutive matrix, including: Where D M is the constitutive matrix of the macro-element integration point located in the macro-micro dual-scale analysis area; K is the vertex coordinate of the micro-scale finite element analysis model associated with the macro-element integration point, i, j = 1, 2; M =K AA -K AN K NN -1 K NA , K NN , K NA , K AN , K AA These are the four sub-matrices of the stiffness matrix of the meso-scale finite element analysis model, which are composed of stiffness matrix elements related to non-vertex node degrees of freedom, non-vertex node degrees of freedom and vertex structural degrees of freedom, vertex node degrees of freedom and non-vertex node degrees of freedom, and vertex node degrees of freedom.

13. The macro-micro dual-scale simulation method for concrete structure damage according to claim 10, characterized in that: The damage variable field of the micro-scale finite element analysis model is homogenized into the damage variables of the macro-element integration points associated with it, including: Where, d M is the damage variable of the macro unit integration point; d i is the damage variable of the i-th integration point in the damage loading area in the mesoscale finite element analysis model; Ω d,i is the integration range of the i-th integration point in the damage loading region in the mesoscale finite element analysis model; Ω d is the extent of the damage loading zone in the meso-scale finite element analysis model.

14. A dual-scale simulation device for concrete structure damage at macro and micro scales, characterized in that: include: The meshing module is used to perform macro-scale finite element meshing on the concrete structure to obtain a macro-grid that includes the entire structural analysis area; The grid initialization module is used to assume that the analysis scale of each macro unit integration point in the macro grid is macro scale at the initial state and is classified as a macro scale analysis area; The regional update module is used to divide the concrete loading process into N incremental steps. For any incremental step, finite element equilibrium iteration is carried out; In any iteration step of any incremental step, it is determined whether each macro-unit integration point in the macro-scale analysis area satisfies the preset analysis scale adaptive transformation criterion; if so, the analysis scale of the corresponding macro-unit integration point is transformed from the macro scale to the macro-micro dual scale and classified as the macro-micro dual scale analysis area; The partition simulation module is used to simulate the damage evolution process of concrete structure using the finite element analysis model, in which the integral points in the macro-scale analysis area adopt the macro-scale finite element analysis model, and the integral points in the macro-micro dual-scale analysis area are associated with the micro-scale finite element analysis model.

15. A storage medium having stored thereon a computer program executable by a processor, characterized in that: When the computer program is executed, the steps of the method according to any one of claims 1 to 13 are implemented.

16. A data processing device comprising a memory and a processor, wherein the memory stores a computer program executable by the processor, wherein: When the computer program is executed, the steps of the method according to any one of claims 1 to 13 are implemented.

Images