Concrete damage analysis method and related equipment

Abstract

The invention relates to the technical field of concrete damage analysis, and provides a concrete damage analysis method and related equipment, and the method comprises the steps: obtaining acoustic emission data of a plurality of damage areas of a target concrete area; determining the damage type of the damage area based on the acoustic emission data of each damage area; for each damage area, calculating an acoustic emission b value of the damage area according to the acoustic emission data of the damage area, and calculating spatial-temporal characteristics of the damage area; determining an associated damage area corresponding to each damage area from all the damage areas; and performing damage analysis on the target concrete area based on the damage types, the acoustic emission b values, the spatial-temporal characteristics and the corresponding associated damage areas of all the damage areas. According to the method, the accuracy of concrete damage analysis can be improved.

Description

Technical Field

[0001] This application relates to the field of concrete damage analysis technology, and in particular to a concrete damage analysis method and related equipment. Background Technology

[0002] Concrete damage identification is a crucial research area in structural health monitoring, stemming from the inherent properties of concrete and the demands of engineering safety. For concrete, tensile and shear damage are the dominant types of cracking damage under external loads. Under different loading conditions, the formation and evolution of cracking damage are guided by different types of damage. These damages not only affect the structural appearance but can also accelerate steel corrosion, reduce structural load-bearing capacity and durability, and even lead to major safety accidents. Traditional manual inspection methods are inefficient and struggle to detect hidden damage. Furthermore, the increasing size and complexity of modern engineering structures place higher demands on the early and accurate identification of damage. Therefore, exploring the classification of concrete structural damage and the laws governing crack evolution is of significant theoretical and engineering importance for accurately assessing the structural damage state.

[0003] In damage identification methods, acoustic emission (AE) technology is widely used in concrete fracture monitoring due to its high sensitivity. The nearest neighbor distance method in seismology provides a new approach for analyzing the correlation of AE events during concrete fracture. This method effectively characterizes the spatiotemporal features of crack evolution by measuring the spatial distance between events. However, the traditional nearest neighbor distance method is sensitive to data scale; if the features are not properly processed, certain dimensions may dominate the distance calculation, thus affecting the accuracy of concrete damage analysis. Summary of the Invention

[0004] This application provides a concrete damage analysis method and related equipment, which can solve the problem of low accuracy in concrete damage analysis.

[0005] In a first aspect, embodiments of this application provide a concrete damage analysis method, which includes:

[0006] Acquire acoustic emission data from multiple damaged areas within the target concrete region;

[0007] The damage type of each damaged area is determined based on the acoustic emission data of that area.

[0008] For each damaged area, the acoustic emission b-value of the damaged area is calculated based on the acoustic emission data of the damaged area, and the spatiotemporal characteristics of the damaged area are also calculated.

[0009] Identify the associated damage regions for each damage region from all damage regions; the associated damage regions are the other damage regions with the greatest spatiotemporal correlation to the damage region.

[0010] Damage analysis is performed on the target concrete area based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas of all damaged areas.

[0011] Optionally, the damage type of each damaged region can be determined based on the acoustic emission data of that region, including:

[0012] For each damaged area, perform the following steps:

[0013] Multiple probability densities of the damaged region are calculated based on acoustic emission data of the damaged region; each probability density corresponds one-to-one with a different damage type.

[0014] The attribution probability of each damage type is calculated based on the probability density corresponding to each damage type; the attribution probability is used to describe the probability that the damaged area belongs to the corresponding damage type.

[0015] The damage type with the highest probability of attribution is taken as the damage type of the damage region.

[0016] Optionally, multiple probability densities of the damaged region can be calculated based on the acoustic emission data of the damaged region, including:

[0017] Through the formula:

[0018] ;

[0019] Calculate the first probability density ;

[0020] in, Acoustic emission data representing the damaged area, Indicates the first A vector of mean values ​​from a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution. Representing dimension, , Indicates the number of damage types.

[0021] Optionally, the attribution probability of each damage type is calculated based on the probability density corresponding to each damage type, including:

[0022] Through the formula:

[0023] ;

[0024] Calculate the first Probability of each damage type ;

[0025] in, Represents latent variables. Indicates the first The mixing coefficients of a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution. Indicates the first A vector of mean values ​​from a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution.

[0026] Optionally, the acoustic emission b-value of the damaged area is calculated based on the acoustic emission data of the damaged area, including:

[0027] Through the formula:

[0028] ;

[0029] Calculate the acoustic emission b-value of the damaged area ;

[0030] in, This represents the amplitude gradation statistical matrix. Represents an empirical constant. This represents the acoustic emission amplitude in the acoustic emission data of the damaged area.

[0031] Optionally, the spatiotemporal characteristics of the damaged region are calculated, including:

[0032] Construct a sphere with the center point of the target concrete area as its center, covering all damaged areas;

[0033] Calculate the Euclidean distance from each damaged area to the center of the ball;

[0034] The maximum radius is defined by the largest Euclidean distance, and multiple radius intervals are divided within the range of 0 to the maximum radius.

[0035] Identify the radius range of the damaged area and calculate the spatiotemporal characteristics of the damaged area based on the number of damaged areas within the radius range.

[0036] Optionally, the spatiotemporal characteristics of the damaged region are calculated based on the number of damaged regions within the radius range of the damaged region, including:

[0037] Through the formula:

[0038] ;

[0039] Calculate the spatiotemporal characteristics of the damaged region ;

[0040] in, This indicates the number of damaged areas within the radius range of the damaged area. Represents a constant. This represents the maximum of the two boundary values ​​within the radius range of the damaged area.

[0041] Optionally, the associated damage regions corresponding to each damage region are determined from all damage regions, including:

[0042] For each damaged area, perform the following steps:

[0043] Calculate the regional distance between the damaged area and each other damaged area;

[0044] For each other damaged area, the spatiotemporal correlation parameters between the damaged area and other damaged areas are calculated based on the area distances corresponding to the other damaged areas.

[0045] The other damaged region with the smallest spatiotemporal correlation parameter is taken as the associated damaged region.

[0046] Optionally, based on the regional distances corresponding to other damaged regions, spatiotemporal correlation parameters between the damaged region and other damaged regions are calculated, including:

[0047] Through the formula:

[0048] ;

[0049] Calculate the first The damage area and the first Spatiotemporal correlation parameters between damaged regions ;

[0050] in, Indicates the first The damage area and the first The temporal proximity value between the damaged areas Indicates the first The damage area and the first The regional distance between each damaged area , , A set of numbers representing the damaged areas.

[0051] Secondly, embodiments of this application provide a concrete damage analysis device, comprising:

[0052] The acquisition module is used to acquire acoustic emission data from multiple damaged areas of the target concrete area;

[0053] The first determining module is used to determine the damage type of each damaged area based on the acoustic emission data of each damaged area.

[0054] The calculation module is used to calculate the acoustic emission b-value of each damaged area based on the acoustic emission data of the damaged area, and to calculate the spatiotemporal characteristics of the damaged area.

[0055] The second determination module is used to determine the associated damage region corresponding to each damage region from all damage regions; the associated damage region is the other damage region with the greatest spatiotemporal correlation with the damage region.

[0056] The damage analysis module is used to perform damage analysis on the target concrete area based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas of all damaged areas.

[0057] Thirdly, embodiments of this application provide a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the concrete damage analysis method described above.

[0058] Fourthly, embodiments of this application provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described concrete damage analysis method.

[0059] The above-mentioned solution in this application has the following beneficial effects:

[0060] In the embodiments of this application, acoustic emission data of multiple damaged areas of the target concrete region are acquired. Then, the damage type of each damaged area is determined based on its acoustic emission data. For each damaged area, its acoustic emission b-value is calculated based on the acoustic emission data, and its spatiotemporal characteristics are also calculated. Next, the associated damaged areas corresponding to each damaged area are identified from all damaged areas. Finally, damage analysis is performed on the target concrete region based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damaged areas of all damaged areas. Determining the damage type of each damaged area and calculating its acoustic emission b-value and spatiotemporal characteristics allows for the analysis of the characteristics of damaged areas of different damage types in the fracture process and in space and time. Identifying associated damaged areas allows for the analysis of the correlation between damaged areas in the target concrete region. Combining damage type, acoustic emission b-value, spatiotemporal characteristics, and associated damaged areas improves the comprehensiveness of the damage analysis, allowing for multi-dimensional analysis of concrete damage and thus enhancing the accuracy of concrete damage analysis.

[0061] Other beneficial effects of this application will be described in detail in the following detailed description section. Attached Figure Description

[0062] To more clearly illustrate the technical solutions in the embodiments of this application, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0063] Figure 1 A flowchart of a concrete damage analysis method provided in an embodiment of this application;

[0064] Figure 2 A schematic diagram of damage categories provided for an embodiment of this application;

[0065] Figure 3 This is a schematic diagram illustrating the change in acoustic emission b-value according to an embodiment of this application;

[0066] Figure 4 This is a schematic diagram illustrating the change of the d value according to an embodiment of this application;

[0067] Figure 5 This is a schematic diagram of the spatiotemporal correlation parameter distribution provided in an embodiment of this application;

[0068] Figure 6 This is a schematic diagram of acoustic emission event distribution provided in an embodiment of this application;

[0069] Figure 7 This is a schematic diagram of tensile damage distribution provided in an embodiment of this application;

[0070] Figure 8 This is a schematic diagram of shear damage distribution provided in an embodiment of this application;

[0071] Figure 9 This is a schematic diagram of the structure of a concrete damage analysis device provided in an embodiment of this application;

[0072] Figure 10 This is a schematic diagram of the structure of a terminal device provided in an embodiment of this application. Detailed Implementation

[0073] In the following description, specific details such as particular system architectures and techniques are set forth for illustrative purposes and not for limitation, in order to provide a thorough understanding of the embodiments of this application. However, those skilled in the art will understand that this application may also be implemented in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, apparatuses, circuits, and methods have been omitted so as not to obscure the description of this application with unnecessary detail.

[0074] It should be understood that, when used in this application specification and the appended claims, the term "comprising" indicates the presence of the described features, integrals, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components and / or a collection thereof.

[0075] It should also be understood that the term “and / or” as used in this application specification and the appended claims means any combination of one or more of the associated listed items and all possible combinations, and includes such combinations.

[0076] As used in this application specification and the appended claims, the term "if" may be interpreted, depending on the context, as "when," "once," "in response to determination," or "in response to detection." Similarly, the phrase "if determined" or "if detected [the described condition or event]" may be interpreted, depending on the context, as meaning "once determined," "in response to determination," "once detected [the described condition or event]," or "in response to detection [the described condition or event]."

[0077] Furthermore, in the description of this application and the appended claims, the terms "first," "second," "third," etc., are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0078] References to "one embodiment" or "some embodiments" as described in this specification mean that one or more embodiments of this application include a specific feature, structure, or characteristic described in connection with that embodiment. Therefore, the phrases "in one embodiment," "in some embodiments," "in other embodiments," "in still other embodiments," etc., appearing in different parts of this specification do not necessarily refer to the same embodiment, but rather mean "one or more, but not all, embodiments," unless otherwise specifically emphasized. The terms "comprising," "including," "having," and variations thereof mean "including but not limited to," unless otherwise specifically emphasized.

[0079] To address the low accuracy of existing concrete damage analysis methods, this application provides a concrete damage analysis method. This method acquires acoustic emission data from multiple damaged areas within a target concrete region. Then, based on the acoustic emission data of each damaged area, it determines the damage type. For each damaged area, it calculates the acoustic emission b-value and its spatiotemporal characteristics. Next, it identifies the associated damaged areas corresponding to each damaged area. Finally, based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damaged areas of all damaged areas, it performs damage analysis on the target concrete region. Specifically, determining the damage type of each damaged area and calculating its acoustic emission b-value and spatiotemporal characteristics allows for the analysis of the characteristics of different damage types in the fracture process and in space and time. Identifying associated damaged areas allows for the analysis of the relationships between damaged areas within the target concrete region. Combining damage type, acoustic emission b-value, spatiotemporal characteristics, and associated damaged areas improves the comprehensiveness of the damage analysis, allowing for multi-dimensional analysis of concrete damage and thus enhancing the accuracy of concrete damage analysis.

[0080] The concrete damage analysis method provided in this application will be described exemplarily below.

[0081] like Figure 1 As shown, the concrete damage analysis method provided in this application includes the following steps:

[0082] Step 11: Obtain acoustic emission data from multiple damaged areas of the target concrete area.

[0083] The target concrete area mentioned above is the concrete area that needs to be damaged (such as building foundation, bridge pier, etc.). The damaged area is the area in the target concrete area where damage occurs. The acoustic emission data includes the average frequency (AF), rise time to amplitude ratio (RA), amplitude, and other data of the acoustic emission signal generated when the concrete undergoes deformation.

[0084] In some embodiments of this application, acoustic emission data of the damaged area can be obtained by deploying devices such as acoustic emission sensors.

[0085] For example, after acquiring acoustic emission data, it needs to be standardized, and the expression is:

[0086] ;

[0087] in, This represents the standardized acoustic emission data. Represents acoustic emission data. This represents the mean of acoustic emission data. This represents the standard deviation of acoustic emission data.

[0088] Step 12: Determine the damage type of each damaged area based on the acoustic emission data of that area.

[0089] The above-mentioned damage types are tensile damage or shear damage.

[0090] In some embodiments of this application, the step of determining the damage type of a damaged region based on the acoustic emission data of each damaged region includes:

[0091] For each damaged area, perform the following steps:

[0092] The first step is to calculate multiple probability densities of the damaged area based on the acoustic emission data of the damaged area.

[0093] Multiple probability densities correspond one-to-one with multiple damage types.

[0094] Specifically, through the formula:

[0095] ;

[0096] Calculate the first probability density .

[0097] in, Acoustic emission data representing the damaged area, Indicates the first A vector of mean values ​​from a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution. Representing dimension, , Indicates the number of damage types.

[0098] The Gaussian distribution described above is the probability distribution of the corresponding damage type.

[0099] The second step is to calculate the attribution probability of each damage type based on the probability density corresponding to each damage type.

[0100] Attribution probability describes the probability that a damaged area belongs to a corresponding damage type.

[0101] Specifically, through the formula:

[0102] ;

[0103] Calculate the first Probability of each damage type .

[0104] in, Represents latent variables. Indicates the first The mixing coefficients of a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution. Indicates the first A vector of mean values ​​from a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution.

[0105] The third step is to select the damage type with the highest probability of attribution as the damage type of the damage region.

[0106] For example, a Gaussian Mixture Model (GMM) can be used to determine the damage type of a damaged region based on its acoustic emission data. Specifically, the GMM parameters are first initialized: weight α initialization: equal weights α1 = α2 = 0.5. mean μ initialization: using the initial cluster centers of K-means as the mean. covariance Σ initialization: initialized as an identity matrix [[1,0],[0,1]], indicating the assumption that there is no correlation between features. Weight α: initial mixing coefficients for each Gaussian distribution; mean μ: center point of each Gaussian distribution; covariance Σ: describes the shape of the distribution.

[0107] GMM assumes that the data is generated by a linear combination of K Gaussian distributions, with the following probability density function:

[0108] ;

[0109] in, Represents a probability distribution. Indicates model parameters, Indicates the first The mixing coefficients of a Gaussian distribution satisfy the following conditions: . :No. The probability density function of a Gaussian distribution:

[0110]

[0111] During the model training phase, the Expectation Maximization Algorithm (EM) is used to estimate the GMM parameters:

[0112] E-step (expected step):

[0113] formula:

[0114] ;

[0115] ;

[0116] in, These are latent variables, representing samples. Does it belong to the first A Gaussian distribution, Indicates sample For the The probability of belonging to a Gaussian distribution.

[0117] M-step (Maximization Step) formula:

[0118] (1) Update weights:

[0119] ;

[0120] (2) Update the mean:

[0121] ;

[0122] (3) Update covariance:

[0123] ;

[0124] Regularization terms (e.g.) This prevents the covariance matrix from becoming singular (non-invertible).

[0125] Convergence criterion: Log-likelihood change < threshold value (e.g., e) -5 ).

[0126] Log-likelihood formula:

[0127] ;

[0128] The following example illustrates this step.

[0129] Damage types such as Figure 2 As shown, Figure 2 This is a schematic diagram of loss types in a three-dimensional Cartesian coordinate system constructed with the vertex of the target concrete region as the origin. In the figure, the horizontal axis represents the X coordinate of the damaged region, the vertical axis represents the Y coordinate, and the vertical axis represents the Z coordinate. Cluster1 represents the damaged region corresponding to tensile damage, and Cluster2 represents the damaged region corresponding to shear damage.

[0130] Step 13: For each damaged area, calculate the acoustic emission b-value of the damaged area based on the acoustic emission data of the damaged area, and calculate the spatiotemporal characteristics of the damaged area.

[0131] The aforementioned acoustic emission b-value describes the ratio of small to large acoustic emission events during concrete cracking, reflecting the concrete cracking process. The aforementioned spatiotemporal characteristics describe the degree of damage in the damaged area over time.

[0132] In some embodiments of this application, the steps of calculating the acoustic emission b-value of the damaged region based on the acoustic emission data of the damaged region and calculating the spatiotemporal characteristics of the damaged region include:

[0133] The first step is to calculate the acoustic emission b-value of the damaged area based on the acoustic emission data of the damaged area.

[0134] Specifically, through the formula:

[0135] ;

[0136] Calculate the acoustic emission b-value of the damaged area .

[0137] in, This represents the amplitude gradation statistical matrix. Represents an empirical constant. This represents the acoustic emission amplitude in the acoustic emission data of the damaged area.

[0138] It should be noted that in the above amplitude grading statistical matrix, each row corresponds to the amplitude of the acoustic emission signal in one acoustic emission data set, and each column corresponds to an amplitude threshold (the amplitude threshold is a preset value). In practical applications, in order to consider the changes in concrete damage over time, a window size is set, and the data is slid according to the time of acoustic emission data acquisition, and all acoustic emission data under each window are statistically analyzed into an amplitude grading statistical matrix.

[0139] For example, the acoustic emission b value is as follows: Figure 3 As shown, Figure 3 The horizontal axis represents the acquisition time of the acoustic emission data, in seconds (s), and the vertical axis represents the acoustic emission b-value. GC242 shows the curve of the acoustic emission b-value changing over time. Figure 3 'a' represents the acoustic emission value 'b' of the damaged area during tensile damage. Figure 3 b represents the acoustic emission b-value of the damaged area due to shear damage.

[0140] It can be seen that the tensile parameter b gradually decreases with fluctuation over time, indicating that tensile cracks are constantly being generated; while in the shear diagram, the most fine shear cracks are generated at 180s, and macroscopic shear cracks are formed at 280s.

[0141] The second step is to calculate the spatiotemporal characteristics of the damaged area.

[0142] First, construct a sphere with the center point of the target concrete area as its center, covering all damaged areas.

[0143] Then, calculate the Euclidean distance from each damaged area to the center of the ball.

[0144] Then, taking the largest Euclidean distance as the maximum radius, multiple radius intervals are divided within the numerical range from 0 to the maximum radius.

[0145] For example, if the maximum radius is 5 meters, then the range of 0 to 5 meters can be divided into 5 radius intervals: 0-1 meters, 1-2 meters, 2-3 meters, 3-4 meters, and 4-5 meters.

[0146] Finally, the radius range of the damaged area is identified, and the spatiotemporal characteristics of the damaged area are calculated based on the number of damaged areas within the radius range.

[0147] Specifically, through the formula:

[0148] ;

[0149] Calculate the spatiotemporal characteristics of the damaged region .

[0150] in, This indicates the number of damaged areas within the radius range of the damaged area. Represents a constant. This represents the maximum of the two boundary values ​​within the radius range of the damaged area.

[0151] For example, the Euclidean distance between the center point of the damaged area and the center of the sphere can be calculated to obtain the Euclidean distance between the damaged area and the center of the sphere. The step size is equal for all radius intervals. If the Euclidean distance between the damaged area and the center of the sphere falls within a radius interval, then the damaged area is considered to be within that radius interval. The formula for calculating the Euclidean distance is:

[0152] ;

[0153] in, Represents Euclidean distance. The x-coordinate of the sphere's center. The ordinate represents the center of the sphere. The vertical coordinate of the sphere's center. The x-coordinate represents the center point of the damaged area. The ordinate represents the center point of the damaged area. The vertical coordinates of the center point of the damaged area.

[0154] For example, the values ​​of spatiotemporal features are as follows: Figure 4As shown in the figure, the horizontal axis represents the acquisition time of acoustic emission data in seconds (s), and the vertical axis represents the d-value (i.e., the value of the spatiotemporal feature). GC242 represents the curve of the d-value changing over time. Figure 4 'a' represents the 'd' value of tensile damage. Figure 4 b represents the d value of shear damage.

[0155] It can be seen that the tensile parameter d gradually increases with the passage of time, and the longer time taken indicates that the specimen undergoes a longer tensile damage accumulation process; the shear parameter d also gradually increases with the passage of time, reaching its lowest point at 250s, which is when a large crack is generated.

[0156] Step 14: Identify the associated damage regions corresponding to each damage region from all damage regions.

[0157] The aforementioned associated damage areas are the other damage areas with the greatest spatiotemporal correlation to the damage area.

[0158] In some embodiments of this application, the step of determining the associated damage region corresponding to each damage region from all damage regions includes:

[0159] For each damaged area, perform the following steps:

[0160] The first step is to calculate the regional distance between the damaged area and each other damaged area.

[0161] For example, the regional distance between two damaged regions can be obtained by calculating the Euclidean distance between the center points of the two damaged regions.

[0162] The second step is to calculate the spatiotemporal correlation parameters between each other damaged area based on the regional distances between the other damaged areas.

[0163] Specifically, through the formula:

[0164] ;

[0165] Calculate the first The damage area and the first Spatiotemporal correlation parameters between damaged regions .

[0166] in, Indicates the first The damage area and the first The temporal proximity value between the damaged areas Indicates the first The damage area and the first The regional distance between each damaged area , , A set of numbers representing the damaged areas.

[0167] For example, the above parameters It can be obtained by calculating the difference in the acquisition time of acoustic emission data from the two damaged areas.

[0168] The third step is to identify the other damaged regions with the smallest spatiotemporal correlation parameters as the associated damaged regions corresponding to the damaged regions.

[0169] For example, the distribution statistics of spatiotemporal correlation parameters are performed, and the distribution results are as follows: Figure 5 As shown in the figure, the horizontal axis represents the logarithm of the spatiotemporal correlation parameter. The vertical axis represents the frequency, and the curve represents the Gaussian fitted function curve. Figure 5 'a' represents the spatiotemporal correlation parameters of tensile damage, and the function curve is... Coefficient of determination , Figure 5 b represents the spatiotemporal correlation parameters of shear damage, and the function curve is... Coefficient of determination .

[0170] Two-dimensional joint density distribution of time-space components, as shown Figure 6 As shown, Figure 6 The horizontal axis represents The vertical axis is Density represents density, and the area inside the ellipse represents background acoustic emission events. It can be seen that background events appear in all areas of the density map. Shear damage mainly occurs near macroscopic cracks, with fewer edge events and no obvious damage in other areas.

[0171] Step 15: Based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas of all damaged areas, perform damage analysis on the target concrete area.

[0172] For example, the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas of each damaged region in the target concrete area are integrated into a single tag data set. Then, the damage status of the target concrete is analyzed based on all the tag data. Three-dimensional visualization technology can be used to display all the tag data, revealing the relationships between regions of different damage types in the target concrete across three dimensions: acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas. This allows for a comprehensive assessment of damage types (including tensile and shear damage), evolution stages, and spatial distribution. Based on this visualization, analysis of damage in the target concrete area can be achieved in terms of damage type, spatiotemporal characteristics, and associated damage, such as analysis of the distribution of damage types and the correlation evolution of damage within the target concrete area.

[0173] For example, the frequency of acoustic emission events (i.e., the areas where damage occurs) of different damage types can be statistically analyzed to assess the damage status of the target concrete area.

[0174] The distribution of acoustic emission events under tensile damage at different thresholds is as follows: Figure 7 As shown in the figure, the horizontal axis represents the Y-axis coordinate, with the unit being millimeters (mm), the vertical axis represents the frequency, and the curve is a Gaussian fitted function curve. Figure 7 'a' represents the original distribution, and the function curve is... Coefficient of determination , Figure 7 b represents The distribution over time, the function curve is Coefficient of determination , Figure 7 c is The distribution over time, the function curve is Coefficient of determination , Figure 7 d is The distribution over time, the function curve is Coefficient of determination .

[0175] The distribution of acoustic emission events under different thresholds for shear damage is as follows: Figure 8 As shown in the figure, the horizontal axis represents the Y-axis coordinate, with the unit being millimeters (mm), the vertical axis represents the frequency, and the curve is a Gaussian fitted function curve. Figure 8 'a' represents the original distribution, and the function curve is... Coefficient of determination , Figure 8 b represents The distribution over time, the function curve is Coefficient of determination , Figure 8 c is The distribution over time, the function curve is Coefficient of determination , Figure 8 d is The distribution over time, the function curve is Coefficient of determination .

[0176] It is evident that acoustic emission (AE) events are more densely distributed in the spatiotemporal region near macroscopic cracks, exhibiting shorter intervals between adjacent events. The analysis results indicate that selecting an appropriate nearest neighbor distance threshold can not only distinguish between background and clustered events but also filter out more closely related events within the clustered events.

[0177] It is worth mentioning that by determining the damage type of the damaged area and calculating the acoustic emission b-value and spatiotemporal characteristics of the damaged area, it is possible to analyze the characteristic information of the damaged area of ​​different damage types in the fracture process and in the spatiotemporal space. It is also possible to identify the associated damaged areas of the damaged area and analyze the correlation between the damaged areas of the target concrete area. By combining the damage type, acoustic emission b-value, spatiotemporal characteristics, and associated damaged areas, damage analysis of the target concrete area can be performed, which improves the comprehensiveness of the damage analysis, analyzes the damage of concrete from multiple dimensions, and thus improves the accuracy of concrete damage analysis.

[0178] The concrete damage analysis device provided in this application is described below as an example.

[0179] like Figure 9 As shown, this application embodiment provides a concrete damage analysis device, the concrete damage analysis device 900 including:

[0180] Acquisition module 901 is used to acquire acoustic emission data of multiple damaged areas in the target concrete area;

[0181] The first determining module 902 is used to determine the damage type of the damaged area based on the acoustic emission data of each damaged area;

[0182] The calculation module 903 is used to calculate the acoustic emission b-value of each damaged area based on the acoustic emission data of the damaged area, and to calculate the spatiotemporal characteristics of the damaged area.

[0183] The second determining module 904 is used to determine the associated damage region corresponding to each damage region from all damage regions; the associated damage region is the other damage region with the greatest spatiotemporal correlation with the damage region.

[0184] The damage analysis module 905 is used to perform damage analysis on the target concrete area based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas of all damaged areas.

[0185] It should be noted that the information interaction and execution process between the above-mentioned devices / units are based on the same concept as the method embodiments of this application. For details on their specific functions and technical effects, please refer to the method embodiments section, and they will not be repeated here.

[0186] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is merely an example. In practical applications, the above functions can be assigned to different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit. Furthermore, the specific names of the functional units and modules are only for easy differentiation and are not intended to limit the scope of protection of this application. The specific working process of the units and modules in the above system can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.

[0187] like Figure 10 As shown, an embodiment of this application provides a terminal device, wherein the terminal device D10 of this embodiment includes: at least one processor D100 ( Figure 10 The diagram shows only one processor, a memory D101, and a computer program D102 stored in the memory D101 and executable on the at least one processor D100, wherein the processor D100 executes the computer program D102 to implement the steps in any of the above method embodiments.

[0188] Specifically, when the processor D100 executes the computer program D102, it acquires acoustic emission data from multiple damaged areas of the target concrete region. Then, based on the acoustic emission data of each damaged area, it determines the damage type of that area. Next, for each damaged area, it calculates the acoustic emission b-value and its spatiotemporal characteristics based on the acoustic emission data. Then, it identifies the associated damaged areas corresponding to each damaged area from all the damaged areas. Finally, based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damaged areas of all damaged areas, it performs damage analysis on the target concrete region. Determining the damage type of a damaged area and calculating its acoustic emission b-value and spatiotemporal characteristics allows for the analysis of the characteristics of damaged areas of different damage types in the fracture process and in space and time. Identifying associated damaged areas allows for the analysis of the correlation between damaged areas in the target concrete region. Combining damage type, acoustic emission b-value, spatiotemporal characteristics, and associated damaged areas to perform damage analysis on the target concrete region improves the comprehensiveness of the damage analysis, allowing for multi-dimensional analysis of concrete damage and thus improving the accuracy of concrete damage analysis.

[0189] The processor D100 can be a central processing unit (CPU), or it can be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor can be a microprocessor or any conventional processor.

[0190] In some embodiments, the memory D101 may be an internal storage unit of the terminal device D10, such as a hard disk or memory of the terminal device D10. In other embodiments, the memory D101 may be an external storage device of the terminal device D10, such as a plug-in hard disk, smart media card (SMC), secure digital card (SD), flash card, etc., equipped on the terminal device D10. Furthermore, the memory D101 may include both internal and external storage units of the terminal device D10. The memory D101 is used to store the operating system, applications, bootloader, data, and other programs, such as the program code of the computer program. The memory D101 can also be used to temporarily store data that has been output or will be output.

[0191] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps described in the various method embodiments above.

[0192] This application provides a computer program product that, when run on a terminal device, enables the terminal device to implement the steps described in the various method embodiments above.

[0193] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it 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 of this application can 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 at least: any entity or device capable of carrying the computer program code to the concrete damage analysis method apparatus / terminal device, a recording medium, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electrical carrier signal, a telecommunication signal, and a software distribution medium. Examples include USB flash drives, portable hard drives, magnetic disks, or optical disks. In some jurisdictions, according to legislation and patent practice, computer-readable media cannot be electrical carrier signals or telecommunication signals.

[0194] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail or recorded in a certain embodiment, please refer to the relevant descriptions of other embodiments.

[0195] Those skilled in the art will recognize that the units and algorithm steps of the various examples 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 implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0196] The above description is the preferred embodiment of this application. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of this invention, and these improvements and modifications should also be considered within the scope of protection of this application.

Claims

1. A method for analyzing concrete damage, characterized in that, include: Acquire acoustic emission data from multiple damaged areas within the target concrete region; The damage type of each damaged region is determined based on the acoustic emission data of that region. For each of the damaged regions, the acoustic emission b-value of the damaged region is calculated based on the acoustic emission data of the damaged region, and the spatiotemporal characteristics of the damaged region are also calculated. Identify the associated damage region corresponding to each of the damage regions from all the damage regions; The associated damage region is the other damage region with the greatest spatiotemporal correlation to the damage region; Damage analysis is performed on the target concrete area based on the damage type, acoustic emission b-value, spatiotemporal characteristics, and corresponding associated damage areas of all damaged areas.

2. The concrete damage analysis method according to claim 1, characterized in that, Determining the damage type of each damaged region based on acoustic emission data includes: For each of the aforementioned damaged areas, the following steps are performed: Based on the acoustic emission data of the damaged region, multiple probability densities of the damaged region are calculated; each probability density corresponds one-to-one with a different damage type. The attribution probability of each damage type is calculated based on the probability density corresponding to each damage type; the attribution probability is used to describe the probability that the damaged area belongs to the corresponding damage type. The damage type with the highest probability of attribution is taken as the damage type of the damage region.

3. The concrete damage analysis method according to claim 2, characterized in that, The calculation of multiple probability densities of the damaged region based on the acoustic emission data of the damaged region includes: Through the formula: ； Calculate the first probability density ; in, Acoustic emission data representing the damaged area, Indicates the first A vector of mean values ​​from a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution. Representing dimension, , Indicates the number of damage types.

4. The concrete damage analysis method according to claim 3, characterized in that, The calculation of the attribution probability of each damage type based on the probability density corresponding to each damage type includes: Through the formula: ； Calculate the first Probability of each damage type ; in, Represents latent variables. Indicates the first The mixing coefficients of a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution. Indicates the first A vector of mean values ​​from a Gaussian distribution. Indicates the first The mixing coefficients of a Gaussian distribution.

5. The concrete damage analysis method according to claim 1, characterized in that, The step of calculating the acoustic emission b-value of the damaged region based on the acoustic emission data of the damaged region includes: Through the formula: ； Calculate the acoustic emission b-value of the damaged area ; in, This represents the amplitude gradation statistical matrix. Represents an empirical constant. This represents the acoustic emission amplitude in the acoustic emission data of the damaged area.

6. The concrete damage analysis method according to claim 1, characterized in that, The calculation of the spatiotemporal characteristics of the damaged region includes: Construct a sphere with the center point of the target concrete area as its center, covering all damaged areas; Calculate the Euclidean distance from each damaged area to the center of the sphere; The maximum radius is defined by the largest Euclidean distance, and multiple radius intervals are divided within the range of 0 to the maximum radius. Identify the radius range in which the damaged area is located, and calculate the spatiotemporal characteristics of the damaged area based on the number of damaged areas within the radius range in which the damaged area is located.

7. The concrete damage analysis method according to claim 6, characterized in that, The step of calculating the spatiotemporal characteristics of the damaged region based on the number of damaged regions within the radius interval of the damaged region includes: Through the formula: ； Calculate the spatiotemporal characteristics of the damaged region ; in, This indicates the number of damaged areas within the radius range of the damaged area. Represents a constant. This represents the maximum of the two boundary values ​​within the radius range of the damaged area.

8. The concrete damage analysis method according to claim 1, characterized in that, The step of determining the associated damage region corresponding to each damage region from all damage regions includes: For each of the aforementioned damaged areas, the following steps are performed: Calculate the regional distance between the damaged area and each other damaged area; For each of the other damaged regions, the spatiotemporal correlation parameters between the damaged region and the other damaged regions are calculated based on the region distance corresponding to the other damaged regions. The other damaged region with the smallest spatiotemporal correlation parameter is taken as the associated damaged region corresponding to the damaged region.

9. The concrete damage analysis method according to claim 8, characterized in that, The calculation of the spatiotemporal correlation parameters between the damaged area and the other damaged areas based on the regional distances corresponding to the other damaged areas includes: Through the formula: ； Calculate the first The damage area and the first Spatiotemporal correlation parameters between damaged regions ; in, Indicates the first The damage area and the first The temporal proximity value between the damaged areas Indicates the first The damage area and the first The regional distance between each damaged area , , A set of numbers representing the damaged areas.

10. 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 concrete damage analysis method as described in any one of claims 1 to 9.

Images