Concrete structure array ultrasonic sound field numerical calculation method and terminal equipment

Abstract

The invention discloses a concrete structure array ultrasonic field numerical calculation method, which comprises the following steps: S1, establishing a two-dimensional geometric model of a concrete structure in COMSOLMultiphysics, the two-dimensional geometric model comprising a concrete base material and defects, and setting a sound source, boundary conditions and a selected physical field on the two-dimensional geometric model to form a finite element model; s2, sound field simulation parameters are set in the finite element model; s3, based on the simulation parameters, COMSOLMultiphysics is operated to carry out sound field calculation and analysis so as to obtain sound field data in the concrete structure, the sound field data comprise time and pressure data sets, based on material physical attributes of mortar and aggregate in the concrete structure, the sound velocity of the mortar and the sound velocity of the aggregate are calculated respectively, and the sound velocity of the mortar and the sound velocity of the aggregate are calculated respectively; calculating the equivalent sound velocity of the ultrasonic wave propagating in the concrete structure through a weighted average method; and S4, presetting a precision requirement, comparing the defect shape, position and geometric dimension of the concrete structure based on a two-stage verification mode, and judging whether the error precision meets the preset precision requirement or not.

Description

Technical Field

[0001] The present invention relates to the technical field of acoustic detection of concrete structures, and in particular to a numerical calculation method for an array ultrasonic sound field of a concrete structure and terminal equipment. Background Art

[0002] Concrete, a commonly used structural material in construction, is widely used in various structures such as bridges, buildings, and roads. Due to the complexity of concrete structures and their potential defects (such as cracks and pores), accurately calculating the acoustic field in concrete structures is crucial for structural health monitoring, design optimization, and defect detection.

[0003] Acoustic imaging technology is an important nondestructive testing method for concrete structural defects. Traditional sound field calculation methods are typically based on a homogeneous medium model. However, concrete is a composite material composed of aggregate and mortar. Sound waves propagate at different speeds in different media, and its internal sound velocity distribution is complex. This makes it difficult for traditional sound velocity measurement methods to accurately capture the actual sound velocity characteristics of concrete, leading to errors in sound field calculations and affecting accuracy in practical applications.

[0004] This method may lead to errors in the calculation of the sound field, affecting the accuracy in practical applications. To solve this problem, the present invention proposes a numerical calculation method for the ultrasonic sound field of a concrete structure array based on equivalent sound velocity. This method combines COMSOL finite element simulation technology to obtain sound field data and MATLAB full matrix method full focus imaging technology to reconstruct the sound field image. The method is then cross-verified with the defect location and size of the actual concrete test block to accurately simulate the sound field distribution in the concrete structure, thereby improving the accuracy and efficiency of non-destructive detection of concrete structure defects. Summary of the Invention

[0005] The purpose of the present invention is to solve the problems of the above-mentioned background technology and to propose a numerical calculation method for the ultrasonic sound field of a concrete structure array and a terminal device.

[0006] The purpose of the present invention can be achieved through the following technical solutions:

[0007] A numerical calculation method for ultrasonic sound field of concrete structure array, comprising:

[0008] S1. Establishing a two-dimensional geometric model of a concrete structure in COMSOL Multiphysics, wherein the two-dimensional geometric model includes a concrete substrate and defects, and setting a sound source, boundary conditions, and selecting a physical field on the two-dimensional geometric model to form a finite element model;

[0009] S2. Setting sound field simulation parameters in the finite element model, wherein the simulation parameters include sound source type, frequency, excitation signal, and material physical properties, wherein the material physical properties include density, Young's modulus, and Poisson's ratio;

[0010] S3. Based on the simulation parameters, running COMSOL Multiphysics to calculate and analyze the sound field to obtain sound field data in the concrete structure, wherein the sound field data includes a time and pressure data group;

[0011] Based on the material physical properties of mortar and aggregate in concrete structure, the sound velocity of mortar and aggregate are calculated respectively, and the equivalent sound velocity of ultrasonic wave propagation in concrete structure is calculated by weighted average method.

[0012] S4. Preset an accuracy requirement and compare the shape, location and geometric dimensions of the defects in the concrete structure based on a two-level verification method to determine whether the error accuracy meets the preset accuracy requirement;

[0013] The two-level verification method includes: full-matrix full-focus imaging verification using MATLAB and array ultrasonic detection verification.

[0014] As a further solution of the present invention: in the step S1, the aggregate in the concrete base material selects a two-level continuous grading scheme. In the two-level continuous grading, the aggregate is added using fine aggregate with a particle size of 5mm to 20mm and coarse aggregate with a particle size of 20mm to 40mm, and the mass ratio of fine aggregate to coarse aggregate is 5.5:4.5. In the concrete base material, the aggregate accounts for 40% and the mortar accounts for 60%.

[0015] As a further solution of the present invention: in step S2, the specific operation of setting the sound source type is: establishing an array element model on the upper surface of the two-dimensional geometric model, the array element model parameters are set to the number of array elements, array element width and array element spacing, the number of array elements is set to 8, the array element width is set to 20 mm, and the array element spacing is set to 25 mm, establishing line segments on the two-dimensional geometric model to replace array element probes, and assigning sound source attributes to each array element.

[0016] As a further solution of the present invention: in step S2, the frequency f0 is set to 50kHz, and the excitation signal is set to a three-cycle Hanning window function, whose expression is: f(t) = 0.5*(1-cos(2π*f0*t) / 3))*sin(2π*f0*t), where f0 is frequency, t is time, and t≤3 / f0.

[0017] As a further solution of the present invention: in the step S3, when obtaining the sound field data in the concrete structure, the two-dimensional geometric model needs to be meshed, and the mesh is a free triangle mesh with a mesh size of V / f0 / 6 to V / f0 / 10, where V is the sound velocity of the aggregate or mortar.

[0018] As a further solution of the present invention: in step S3, the calculation formulas of the sound velocity of the mortar and the sound velocity of the aggregate are as follows:

[0019]

[0020] Among them, V 骨料 is the aggregate sound velocity, V 砂浆 is the mortar sound velocity, E 骨料 is the Young's modulus of the aggregate, E 砂浆 is the Young's modulus of the mortar, ρ 骨料 is the density of aggregate, ρ 砂浆 is the density of mortar, σ 骨料 is the Poisson's ratio of aggregate, σ 砂浆 is the Poisson's ratio of the mortar.

[0021] As a further embodiment of the present invention, in step S3, the sound field data is collected and calculated using a parametric array element scanning method, with each element independently calculated. The output time step is range(0, T0 / 20, 25*T0), where T0 = 1 / f0, the calculated propagation time is 25 times T0, and the sound field data is collected every T0 / 20.

[0022] As a further solution of the present invention: in step S3, the calculation formula of the equivalent sound speed is as follows: Among them, V 骨料 is the aggregate sound velocity, V 砂浆 is the mortar sound speed, S 骨料 is the aggregate area, S 砂浆 is the mortar area, S 骨料 With S 砂浆 Both are two-dimensional target area.

[0023] As a further solution of the present invention: in the step S4, the operation mode of the MATLAB full-matrix method full-focus imaging verification is: comparing the sound field image generated by the MATLAB full-focus imaging with the actual produced concrete test block containing defects, and verifying whether the position and size of the MATLAB defect imaging at the equivalent sound speed meet the preset accuracy requirements; the operation mode of the array ultrasonic detection verification is: inputting the equivalent sound speed into the array ultrasonic detection equipment, detecting the defects of the actual produced concrete test block containing defects, obtaining the test data file and importing it into the concrete design and analysis software, generating a defect sound field image, comparing the output defect sound field image with the defect position and size of the actual concrete test block, and verifying whether the position and size of the array ultrasonic detection defect imaging at the equivalent sound speed meet the preset accuracy requirements; if either the MATLAB full-matrix method full-focus imaging verification or the array ultrasonic detection verification verification does not meet the preset accuracy requirements, the Young's modulus of the material is adjusted, and steps S2 to S4 are repeated until both verification methods meet the preset accuracy requirements.

[0024] A terminal device is applicable to a numerical calculation method of ultrasonic sound field of concrete structure array, comprising: a processor, wherein the processor is used to run COMSOL Multiphysics and MATLAB.

[0025] Beneficial effects of the present invention:

[0026] (1) In the present invention, the influence of concrete defects on the sound field is reflected by calculating the equivalent sound velocity, thereby improving the accuracy of the numerical calculation of the sound field.

[0027] (2) In the present invention, by combining simulation and imaging technology, the parameter settings of the sound field calculation are effectively optimized, thereby improving the efficiency and accuracy of numerical calculation.

[0028] (3) In the present invention, the MATLAB full-matrix method and full-focus imaging technology are used to verify each other with array ultrasonic detection, which increases the mutual verification with the actual concrete structure defect detection imaging test, helps to verify the accuracy of the finite element model and optimize the simulation results, and improves the reliability of the results. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] The present invention will be further described below with reference to the accompanying drawings.

[0030] Figure 1 It is a flow chart of the numerical calculation method of the ultrasonic sound field of the concrete structure array of the present invention;

[0031] Figure 2 This is a flow chart of the numerical calculation method for the ultrasonic sound field of a concrete structure array according to the present invention when it is applied;

[0032] Figure 3This is a schematic diagram of the finite element model of the concrete defect during COMSOL Multiphysics acoustic field simulation of the present invention;

[0033] Figure 4 yes Figure 3 Enlarged view of point A in the middle;

[0034] Figure 5 is a schematic diagram of the array ultrasonic detection device of the present invention when in operation;

[0035] Figure 6 This is the full-focus imaging diagram of the MATLAB full-matrix method of the present invention;

[0036] Figure 7 This is an imaging diagram of the concrete structure under detection by array ultrasonic detection equipment of the present invention. DETAILED DESCRIPTION

[0037] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making any creative efforts shall fall within the scope of protection of the present invention.

[0038] The present invention provides a terminal device comprising a processor for running COMSOL Multiphysics and MATLAB, so that the terminal device can process large-scale finite element models and complex sound field data.

[0039] Specifically, COMSOL Multiphysics and MATLAB software need to be installed and configured in the terminal device, and ensure that the software versions are compatible and have the required functional modules. COMSOL Multiphysics is configured with acoustic modules and finite element modules, and MATLAB is configured with imaging processing modules.

[0040] See also Figure 1-Figure 2 As shown, the present invention is a numerical calculation method for the ultrasonic sound field of a concrete structure array, comprising:

[0041] S1. Establish a two-dimensional geometric model of the concrete structure in COMSOL Multiphysics. The two-dimensional geometric model includes the concrete substrate and defects.

[0042] See also Figure 3-Figure 4As shown, the two-dimensional geometric model of the concrete structure can be created by first setting a shape and size, and then creating it according to the preset shape and size. In this embodiment, the shape and size are a rectangular target area of ​​400mm x 265mm. Prefabricated defects include but are not limited to holes and cracks. In this embodiment, the prefabricated defect is a hole with a diameter of 75mm, its center is on the centerline of the target area, and the defect is 110mm away from the measurement surface.

[0043] A two-stage continuous gradation scheme was used for the aggregate in the concrete matrix. Fine aggregate with a particle size of 5mm to 20mm and coarse aggregate with a particle size of 20mm to 40mm were used, with a mass ratio of 5.5:4.5 for fine aggregate to coarse aggregate. This was used to create different aggregate shapes and evenly distribute them within the two-dimensional target area. Aggregate with a particle size below 5mm was considered the cement mortar matrix and was not subjected to aggregate generation. In the concrete matrix, aggregate accounted for 40% and mortar accounted for 60%.

[0044] Furthermore, sound sources, boundary conditions, and physical fields are set on the two-dimensional geometric model to form a finite element model.

[0045] S2. Set the sound field simulation parameters in the finite element model. The simulation parameters include the sound source type, frequency, excitation signal, and material physical properties. The material physical properties include density, Young's modulus, and Poisson's ratio.

[0046] The specific operation of setting the sound source type is: establishing an array element model on the upper surface of the two-dimensional geometric model, and the array element model parameters include the number of array elements, array element width and array element spacing.

[0047] In this embodiment, the number of array elements is set to 8, the array element width is set to 20 mm, the array element spacing is set to 25 mm, line segments are established on the two-dimensional geometric model to replace the array element probes, and sound source attributes are assigned to each array element. The sound source attributes include but are not limited to sound source intensity and radiation pattern.

[0048] The selected frequency f0 is set to 50kHz, and the excitation signal is set to a three-cycle Hanning window function, whose expression is:

[0049] f(t)=0.5*(1-cos(2π*f0*t) / 3))*sin(2π*f0*t);

[0050] Where f0 is frequency, t is time, and t≤3 / f0.

[0051] The physical field is solid mechanics, and the boundary conditions are that the left and right segments of the concrete model are set as low-reflection boundaries, and the segments of the array element model are all set to specified displacements so that their longitudinal propagation displacement is f(t).

[0052] Then determine the material physical properties of the concrete structure, as shown in the following table:

[0053]

[0054] The table above shows the density, Young's modulus and Poisson's ratio of mortar and aggregate after multiple optimizations at the preset accuracy.

[0055] S3. Based on the simulation parameters, run COMSOL Multiphysics to calculate and analyze the acoustic field to obtain acoustic field data in the concrete structure. This data includes time and pressure data sets. Based on the physical properties of the mortar and aggregate in the concrete structure, the sound velocity of the mortar and aggregate is calculated, respectively. The equivalent sound velocity of the ultrasonic wave propagating in the concrete structure is calculated using a weighted average method.

[0056] Specifically, to obtain acoustic field data in concrete structures, the two-dimensional geometric model needs to be meshed using a free triangle mesh with a mesh size of V / f0 / 6 to V / f0 / 10, where V is the velocity of sound in the aggregate or mortar.

[0057] The average mesh quality of aggregate and mortar is 0.8 ± 0.2 to ensure the accuracy of simulation results, and the mesh is refined to accurately capture the changes in the acoustic field around the defects.

[0058] Furthermore, the calculation formulas for the sound velocity of mortar and aggregate are as follows:

[0059]

[0060] Among them, V 骨料 is the aggregate sound velocity, V 砂浆 is the mortar sound velocity, E 骨料 is the Young's modulus of the aggregate, E 砂浆 is the Young's modulus of the mortar, ρ 骨料 is the density of aggregate, ρ 砂浆 is the density of mortar, σ 骨料 is the Poisson's ratio of aggregate, σ 砂浆 is the Poisson's ratio of the mortar.

[0061] Keeping one decimal place, the calculated aggregate sound velocity is 3116.8m / s and the mortar sound velocity is 2690.5m / s.

[0062] The formula for calculating the equivalent sound velocity of ultrasonic waves propagating in concrete using the weighted average method is as follows:

[0063]

[0064] Among them, V 骨料 is the aggregate sound velocity, V 砂浆 is the mortar sound speed, S 骨料 is the aggregate area, S砂浆 is the mortar area, S 骨料 With S 砂浆 Both are two-dimensional target area.

[0065] The calculated equivalent sound speed is rounded to 2861 m / s.

[0066] In this embodiment, the sound field data is collected and calculated using a data acquisition and calculation method of an array element parametric scan, and each array element is calculated independently, with the output time step being range(0, T0 / 20, 25*T0).

[0067] Where T0=1 / f0, the calculated propagation time is 25 times T0, and the sound field data is collected every T0 / 20. The sound field data includes a full matrix data set of time and pressure.

[0068] S4. Preset an accuracy requirement and compare the shape, location, and geometric dimensions of the concrete structure defects using a two-stage verification method to determine whether the error accuracy meets the preset accuracy requirement. The two-stage verification method includes: MATLAB full-matrix method full-focus imaging verification and array ultrasonic testing verification.

[0069] See also Figure 6 As shown in the figure, the operation method of the MATLAB full matrix method for full focusing imaging verification is as follows: the sound field data obtained in step S3 is imported into the MATLAB software, and the number of array elements N, the sound wave frequency f, the distance d0 between the first probe and the leftmost end, the probe spacing d, the sampling interval t0, the workpiece size (L×W) and the equivalent sound velocity V of the full focusing imaging are input into MATLAB. 等效 , where the i-th probe is the transmitting probe and the j-th is the receiving probe, where i = 1, 2, 3...N, j = 1, 2, 3...N.

[0070] Create an empty matrix P of the same size as the original data. For each pixel (xn, zm), use the following formula to calculate the distance between the transmitting probe and the receiving probe to the focal point:

[0071]

[0072] Calculate the signal arrival time t based on the distance and speed of sound, where And extract the corresponding signal value. Use the following formula to calculate the diffusion angle θ between the probes:

[0073]

[0074] If the diffusion angle θ between the probes is within a constant range and the signal arrival time is within the data length, the signal value is modulated and accumulated into the image matrix P.

[0075] The modulation process is multiplied by the function cos(tanθ) 4 , where tanθ represents the rate of change of the signal's phase, and cos(tanθ) represents the phase-modulated signal. Raising the result of this transformation to the fourth power can be used in filter design or signal analysis. Specifically, raising it to the fourth power amplifies certain signal characteristics (such as high-frequency components or phase characteristics), which can be used to enhance or suppress specific frequency components in filter design.

[0076] The accumulated image matrix P (including the time and pressure matrix data sets) is then processed using the Hilbert transform function to enhance edges and improve imaging quality, thereby generating a high-resolution sound field image.

[0077] In the MATLAB all-focus imaging algorithm, the imagesc function is used to display image data, including the image's color map, color bar, axis labels, and title. The function maps the values ​​of the image matrix P to a colormap, normalizes the image matrix, and displays the result.

[0078] The reconstructed sound field image is analyzed to identify the defect location and size, and the defect location and size are compared with those of the actual defective concrete test blocks to verify whether the location and size of the MATLAB defect imaging at the equivalent sound speed meet the preset accuracy requirements.

[0079] In this embodiment, the preset accuracy requirements for defect location and size are: ≤5% for both. The acoustic field image obtained using MATLAB full-focus imaging technology shows that the defect signal center is 113.64 mm from the measurement surface, with an accuracy error of 3.31%. The defect size is 72.74 mm, with an accuracy error of 3.01%, both meeting the preset accuracy requirements.

[0080] See also Figure 5 、 Figure 7 As shown, the array ultrasonic testing verification operation method is: input the equivalent sound velocity value into the array ultrasonic testing equipment, and detect the defects of the actual produced concrete test blocks with defects to obtain the test data file. The ultrasonic testing equipment model used in this embodiment is: A 1040MIRA 3D array ultrasonic testing instrument.

[0081] The acquired test data file is then saved and imported into the concrete design and analysis software (IntroViewConcrete 7.0) to generate a defect acoustic field image. The generated acoustic field image is analyzed to identify the defect location and size. This image is then compared with the defect location and size of the actual concrete test block to verify that the location and size of the defect imaged by array ultrasonic testing at equivalent sound speed meet the preset accuracy requirements.

[0082] The sound field image obtained using the array ultrasonic testing equipment shows that the center of the defect signal is 112.5mm away from the measuring surface, with an accuracy error of 2.27%, and the defect size is 78mm, with an accuracy error of 4%, both of which meet the preset accuracy requirements.

[0083] It can be understood that if either the MATLAB full-matrix method full-focus imaging verification or the array ultrasonic detection verification verification method does not meet the preset accuracy requirements, it is necessary to adjust and optimize the Young's modulus of the concrete structure material, and repeat steps S2 to S4 until both verification methods meet the preset accuracy requirements and complete the numerical calculation of the concrete structure array ultrasonic sound field.

[0084] There are generally three ways to adjust and optimize the Young's modulus of concrete structural materials, namely: optimizing the Young's modulus of aggregate, optimizing the Young's modulus of mortar, and adjusting the ratio of aggregate and mortar to achieve the optimization of the Young's modulus of concrete structure. The optimization method is an existing technology, so it will not be described in detail here.

[0085] The inventors found that the acoustic imaging technology in the existing technology cannot fully take into account the non-uniformity and internal defects of the concrete structure, which leads to large errors in the sound field calculation, affecting the accuracy, and lacks mutual verification with actual concrete structure defect detection imaging tests.

[0086] In this application, by combining COMSOL finite element simulation, MATLAB full matrix method full focus imaging technology and actual concrete test block array ultrasonic detection technology, the calculation accuracy of the concrete structure sound field is improved and comparative verification is added.

[0087] Compared with traditional acoustic imaging technology and finite element calculation, the present invention reflects the influence of concrete defects on the sound field by calculating the equivalent sound velocity, thereby improving the accuracy of the sound field numerical calculation.

[0088] MATLAB full-matrix method and full-focus imaging technology provide high-resolution sound field images, which help verify the accuracy of finite element models and optimize simulation results.

[0089] By combining simulation and imaging technology, the parameter settings of sound field calculation are effectively optimized, and the efficiency and accuracy of numerical calculation are improved.

[0090] The mutual verification between MATLAB full-matrix method and full-focus imaging technology and array ultrasonic testing is used to increase the mutual verification with the actual concrete structure defect detection imaging test, which helps to verify the accuracy of the finite element model and optimize the simulation results, thereby improving the reliability of the results.

[0091] The above is a detailed description of an embodiment of the present invention, but the content is only a preferred embodiment of the present invention and should not be considered to limit the scope of the present invention. All equivalent changes and improvements made within the scope of the present invention should still fall within the scope of the claims of the present invention.

Claims

1. A numerical calculation method for ultrasonic sound field of concrete structure array, characterized in that: include: S1. Establishing a two-dimensional geometric model of a concrete structure in COMSOL Multiphysics, wherein the two-dimensional geometric model includes a concrete substrate and defects, and setting a sound source, boundary conditions, and selecting a physical field on the two-dimensional geometric model to form a finite element model; S2. Setting sound field simulation parameters in the finite element model, wherein the simulation parameters include sound source type, frequency, excitation signal and material physical properties, wherein the material physical properties include density, Young's modulus and Poisson's ratio; S3, based on the simulation parameters, running COMSOL Multiphysics to calculate and analyze the sound field to obtain the sound field data in the concrete structure, wherein the sound field data includes a time and pressure data group; Based on the material physical properties of mortar and aggregate in concrete structure, the sound velocity of mortar and aggregate are calculated respectively, and the equivalent sound velocity of ultrasonic wave propagating in concrete structure is calculated by weighted average method. S4. Preset an accuracy requirement, compare the defect shape, position and geometric dimensions of the concrete structure based on a two-level verification method, and determine whether the error accuracy meets the preset accuracy requirement; The two-level verification method includes: full-matrix method full-focus imaging verification using MATLAB and array ultrasonic detection verification.

2. The method for numerically calculating the ultrasonic sound field of a concrete structure array according to claim 1 is characterized in that: In the step S1, a two-level continuous grading scheme is selected for the aggregate in the concrete base material. In the two-level continuous grading, the aggregate is placed using fine aggregate with a particle size of 5mm to 20mm and coarse aggregate with a particle size of 20mm to 40mm, and the mass ratio of fine aggregate to coarse aggregate is 5.5:4.

5. In the concrete base material, the aggregate accounts for 40% and the mortar accounts for 60%.

3. The method for numerically calculating the ultrasonic sound field of a concrete structure array according to claim 1 is characterized in that: In the step S2, the specific operation of setting the sound source type is: establishing an array element model on the upper surface of the two-dimensional geometric model, and the array element model parameters are set to include the number of array elements, array element width and array element spacing; The number of array elements is set to 8, the array element width is set to 20 mm, the array element spacing is set to 25 mm, line segments are established on the two-dimensional geometric model to replace array element probes, and sound source attributes are assigned to each array element.

4. A numerical calculation method for ultrasonic sound field of concrete structure array according to claim 3, characterized in that: In step S2, the frequency f0 is set to 50 kHz, and the excitation signal is set to a three-cycle Hanning window function, which is expressed as: f(t)=0.5*(1-cos(2π*f0*t) / 3))*sin(2π*f0*t); Where f0 is frequency, t is time, and t≤3 / f0.

5. A numerical calculation method for ultrasonic sound field of concrete structure array according to claim 4, characterized in that: In step S3, when obtaining the sound field data in the concrete structure, the two-dimensional geometric model needs to be meshed, and the mesh is a free triangle mesh with a mesh size of V / f0 / 6 to V / f0 / 10, where V is the sound velocity of aggregate or mortar.

6. A numerical calculation method for ultrasonic sound field of concrete structure array according to claim 5, characterized in that: In step S3, the calculation formulas of the sound velocity of the mortar and the sound velocity of the aggregate are as follows: Among them, V 骨料 is the aggregate sound velocity, V 砂浆 is the mortar sound velocity, E 骨料 is the Young's modulus of aggregate, E 砂浆 is the Young's modulus of the mortar, ρ 骨料 is the density of aggregate, ρ 砂浆 is the density of mortar, σ 骨料 is the Poisson's ratio of aggregate, σ 砂浆 is the Poisson's ratio of the mortar.

7. A numerical calculation method for ultrasonic sound field of concrete structure array according to claim 6, characterized in that: In step S3, the sound field data is collected and calculated using a data acquisition method of an array element parametric scan, and each array element is independently calculated, and the output time step is range (0, T0 / 20, 25*T0); Where T0=1 / f0, the calculated propagation time is 25 times T0, and the sound field data is collected every T0 / 20.

8. The method for numerically calculating the ultrasonic sound field of a concrete structure array according to claim 6 is characterized in that: In step S3, the calculation formula of the equivalent sound speed is as follows: Among them, V 骨料 is the aggregate sound velocity, V 砂浆 is the mortar sound velocity, S 骨料 is the aggregate area, S 砂浆 is the mortar area, S 骨料 With S 砂浆 Both are two-dimensional target area.

9. The method for numerically calculating the ultrasonic sound field of a concrete structure array according to claim 1 is characterized in that: In step S4, the operation mode of the MATLAB full-matrix method full-focus imaging verification is: comparing the sound field image reconstructed by the MATLAB full-focus imaging with the actually produced concrete test block containing defects, and verifying whether the position and size of the MATLAB defect imaging at the equivalent sound speed meet the preset accuracy requirements; The operation mode of array ultrasonic detection verification is as follows: input the equivalent sound velocity into the array ultrasonic detection equipment, detect the defects of the concrete test block with defects actually produced, obtain the test data file and import it into the concrete design and analysis software, generate the defect sound field image, compare the output defect sound field image with the defect position and size of the actual concrete test block, and verify whether the position and size of the defect imaging of the array ultrasonic detection at the equivalent sound velocity meet the preset accuracy requirements; If either the MATLAB full-matrix method full-focus imaging verification or the array ultrasonic detection verification does not meet the preset accuracy requirements, the Young's modulus of the material is adjusted, and steps S2 to S4 are repeated until both verification methods meet the preset accuracy requirements.

10. A terminal device, applicable to the numerical calculation method of ultrasonic sound field of concrete structure array according to any one of claims 1 to 9, characterized in that: include: A processor for running COMSOL Multiphysics and MATLAB.

Images