Detonating cord output power evaluation method and system

Abstract

The present application relates to the field of performance evaluation of initiating explosive, and particularly relates to a method and system for evaluating the output power of detonating cord. The method comprises the following steps: firstly, obtaining a relative height field matrix of the dent of the steel block, which represents the three-dimensional micro-topography of the dent of the steel block; then, determining the detonation shock center of the relative height field matrix and converting it into a polar coordinate energy flow field, and calculating the gradient data along the radial direction; based on the radiation attenuation law of the shock wave center, constructing and calculating the asymmetry penalty coefficient representing the difference in the circumferential energy distribution, and the calculation includes the weighting processing based on the distance attenuation law; finally, obtaining the physical maximum depth, using the coefficient to nonlinearly correct the physical maximum depth to obtain the effective power value and determine the eligibility. By introducing the physical field characteristic analysis, the present application can effectively identify and eliminate the pseudo-qualified products with the depth meeting the standard but the detonation wave front being inclined, and significantly improves the accuracy and reliability of the acceptance test of the detonating cord.

Description

Technical Field

[0007] ,

[0006] , ,

[0005] , , ,

[0001] The present invention relates to the field of performance evaluation of initiating explosive devices, and particularly to a method and system for evaluating the output power of detonating cords. Background Art

[0002] A detonating cord is a key core component for transmitting high-energy detonation waves in the separation system of aerospace launch vehicles, the canopy ejection system of fighter jets, and the detonating train of weapons. In actual application scenarios, the detonating cord needs to stably and accurately transmit the detonation energy to the next-stage initiating explosive device within an extremely short time. The stability of its output power directly determines the success or failure of the entire stage separation or emergency rescue system. Therefore, during the product development and mass production acceptance stages, its output power must be extremely strictly evaluated.

[0003] Currently, the commonly used acceptance standard in the industry mainly adopts the steel block indentation method. In this method, the detonating cord is detonated on a steel block of a specified material, and then a depth gauge or microscope is used to measure the depth of the deepest point of the indentation on the surface of the steel block. If this physical depth reaches the preset index, the batch of detonating cords is judged to be qualified. Although the operation of this method is simple, its physical essence is based on the measurement of a single-point geometric quantity, defaulting to the linear assumption that depth is power, and ignoring the complexity of the spatial distribution of explosion energy.

[0004] However, in the actual production process, due to the limited charging process, there may be defects in the detonating cord such as uneven density of the core charge, eccentricity of the core, or inconsistent wall thickness of the casing. These defects will cause the detonation wave front to tilt during propagation, resulting in the so-called "crooked handle" phenomenon. In this case, the explosion energy will shift severely to one side. Although the local deepest point it causes on the steel block may still meet the requirement of 0.3 mm, its effective working area is significantly reduced, and the energy coupling efficiency is significantly decreased. The existing single-depth measurement technology cannot sense this three-dimensional spatial energy skew, resulting in such pseudo-qualified products flowing into the subsequent assembly links, which is extremely likely to cause major hidden dangers such as interruption of the energy transfer in the detonation channel or insufficient initiation ability. Therefore, there is an urgent need for an evaluation method that can break through the limitation of single geometric quantity measurement, deeply analyze the three-dimensional microscopic morphology of the indentation, quantify the stability of the explosion process, and correct the depth index in combination with the physical field characteristics, so as to solve the technical problem that the existing technology cannot accurately identify the energy skew defect. Summary of the Invention

[0005] In order to solve the technical problem of how to sense the energy skew in three-dimensional space, the present invention provides a method and system for evaluating the output power of detonating cords.

[0006] In the first aspect, the present invention provides a method for evaluating the output power of detonating cords, adopting the following technical solution:

[0007] A method for evaluating the output power of detonating cords includes the steps:

[0008] Obtain the relative height field matrix of the steel block indentation, which characterizes the three-dimensional micromorphology of the steel block indentation;

[0009] The detonation impact center of the relative height field matrix is ​​determined, a polar coordinate system is established with the detonation impact center as the origin, the relative height field matrix is ​​converted into a polar coordinate energy flow field, and the gradient data along the radial direction is calculated.

[0010] Based on the polar coordinate energy flow field and the gradient data, an asymmetric penalty coefficient is constructed and calculated. The asymmetric penalty coefficient characterizes the energy distribution difference of the explosion shock wave in the circumferential direction, and the calculation process includes weighting based on the distance attenuation law.

[0011] The maximum depth of the indentation in the steel block is obtained, and the maximum depth is nonlinearly corrected using the asymmetric penalty coefficient to obtain the effective power value. The qualification of the detonating cord's output power is then determined based on the effective power value.

[0012] This invention establishes a two-dimensional evaluation system of depth and energy distribution by introducing full-field three-dimensional topographic data, thereby improving the accuracy of detonating cord output power assessment.

[0013] Preferably, one method for obtaining the relative height field matrix of the steel block indentation specifically includes:

[0014] Control the four-directional ring light source to illuminate in a time-division manner, and use a vertical industrial camera to collect multiple grayscale images of the steel block indentation.

[0015] Detect the pixel grayscale values ​​in the grayscale image and remove specular reflection saturation points whose grayscale values ​​exceed the camera's dynamic range threshold;

[0016] Using the remaining unsaturated image data, the normal vector field of the steel block indentation surface is calculated based on the Lambertian reflection model;

[0017] The relative height field matrix is ​​obtained by integrating and reconstructing the normal vector field.

[0018] This invention utilizes multi-angle illumination information to effectively overcome the strong reflective interference of metal surfaces, and can quickly obtain high-precision microstructures without the need for expensive laser equipment, making it suitable for mass production line inspection.

[0019] Preferably, another method for obtaining the relative height field matrix of the steel block indentation specifically includes:

[0020] The single-axis motion module equipped with a line laser profilometer is controlled to sweep across the surface of the steel block indentation at a uniform speed to collect point cloud depth data;

[0021] The point cloud depth data is subjected to gridded interpolation and Gaussian smoothing filtering to remove speckle noise and generate the relative height field matrix.

[0022] This invention directly acquires absolute depth data, is unaffected by changes in ambient light, exhibits high data stability, and is suitable for laboratory environments with extremely high measurement accuracy requirements.

[0023] Preferably, determining the detonation impact center of the relative height field matrix specifically includes:

[0024] Calculate the geometric moments of the relative height field matrix;

[0025] The geometric centroid of the relative height field matrix is ​​determined by the ratio of the first moment to the zeroth moment, and the geometric centroid is marked as the detonation impact center.

[0026] Preferably, the method for constructing and calculating the asymmetric penalty coefficient specifically includes:

[0027] ;

[0028] ;

[0029] in, This indicates that the radius is in polar coordinates. Circumferential energy variation factor at the location; Indicates the radius and angle The radial gradient value at that location; Represents radius The average radial gradient at all angles; Pi; () is an exponential function with the natural constant as its base; The radius of the current calculation point; The attenuation constant characterizing the radiation properties at the center of the shock wave; This represents the asymmetric penalty coefficient; This represents the maximum radius of the effective computational domain of the relative height field matrix; Represents radius The average depth at that location; This indicates the maximum depth of the dent in the steel block.

[0030] This invention introduces an exponential decay weight, which forces the focus to the core region that best reflects the quality of the core, thereby effectively suppressing unrelated deformation noise caused by stress reflection at the edge of the steel block, resulting in higher accuracy of the calculated asymmetric penalty coefficient.

[0031] Preferably, another method for constructing and calculating the asymmetric penalty coefficient specifically includes:

[0032] Perform a fast Fourier transform on the circumferential height sequence at the same radius in the polar coordinate energy flow field;

[0033] Extract the sum of energy of the transformed non-DC harmonic components and multiply it by an attenuation weight that increases with the radius to obtain the frequency domain variation factor;

[0034] The frequency domain variation factors at all radii are weighted, accumulated, and normalized to obtain the asymmetric penalty coefficient.

[0035] This invention utilizes frequency domain features to separate overall depth and shape distortion, enabling more sensitive detection of elliptic or polygonal distortions of indentations that are difficult to detect with the naked eye, thereby improving the detection rate of minute instability defects.

[0036] Preferably, the value of the attenuation constant is set to be between one-third and one-half of the radius of the indentation on the steel block.

[0037] Preferably, the method for nonlinearly correcting the maximum depth using the asymmetric penalty coefficient specifically includes:

[0038]

[0039] in, This indicates the corrected effective power value; This indicates the maximum depth of the dent in the steel block as measured. ( ) represents the hyperbolic tangent function; This represents the instability sensitivity coefficient related to the hardness of the steel block material. This represents the asymmetric penalty coefficient.

[0040] Preferably, the value of the instability sensitivity coefficient is positively correlated with the hardness of the steel block material. For 45# steel, the value of the instability sensitivity coefficient ranges from 3.0 to 4.0.

[0041] Secondly, the present invention provides a detonating cord output power assessment system, which adopts the following technical solution:

[0042] A detonating cord output power assessment system includes a processor and a memory, wherein the memory stores computer program instructions, and when the computer program instructions are executed by the processor, the aforementioned detonating cord output power assessment method is implemented.

[0043] By adopting the above technical solution, a computer program is generated from the above-mentioned method for evaluating the output power of a detonating cord, and stored in a memory for loading and execution by a processor. This allows for the creation of a terminal device based on the memory and processor, making it convenient to use.

[0044] The present invention has the following technical effects:

[0045] This invention constructs a circumferential energy variation factor to map the three-dimensional microscopic morphological features of an indentation into an asymmetric penalty coefficient, and effectively corrects the physical maximum depth. This method can effectively identify and eliminate pseudo-qualified products that meet the physical depth requirements but exhibit severe asymmetry in energy output distribution, thereby significantly improving the detonation reliability and safety of aerospace-grade detonating cord during batch production acceptance.

[0046] Furthermore, this invention incorporates the propagation characteristics of explosive shock waves in a medium and introduces a center radiation attenuation weighting mechanism based on radial distance. This mechanism considers that the deformation of the core region is directly determined by the explosive pressure, while the deformation of the edge region is significantly affected by the reflection of stress waves from the medium. This allows the evaluation algorithm to actively suppress unrelated deformation noise at the edges of the steel block, focusing feature extraction on the core region that best characterizes the quality of the explosive core, thereby obtaining more accurate evaluation results.

[0047] Furthermore, this invention employs high dynamic range photometric stereo vision technology. By setting an adaptive threshold based on a percentage of the dynamic range, it achieves effective removal of specular reflection highlights, reconstructing a high-precision deformation field without the need for expensive laser scanning equipment. This solution boasts low hardware cost and fast imaging speed, meeting the high-throughput full-inspection requirements of industrial production lines and overcoming the shortcomings of traditional contact measurement methods, such as low efficiency and susceptibility to human error. Attached Figure Description

[0048] Figure 1 This is a flowchart of a method for evaluating the output power of a detonating cord according to an embodiment of the present invention;

[0049] Figure 2 A schematic diagram of the polar coordinate energy flow field provided in an embodiment of the present invention;

[0050] Figure 3 This is a schematic diagram of the circumferential energy variation factor provided in an embodiment of the present invention;

[0051] Figure 4 This is a schematic diagram illustrating the final evaluation conclusion provided in an embodiment of the present invention. Detailed Implementation

[0052] This invention discloses a method for evaluating the output power of detonating cord, referring to... Figure 1 This includes steps S1-S4:

[0053] S1: Obtain the relative height field matrix of the steel block indentation, which characterizes the three-dimensional micromorphology of the steel block indentation.

[0054] It is important to note that in explosion mechanics, the steel block indentation is not merely a deep pit; its internal surface microscopic undulations, such as textures and ridge orientations, faithfully record the historical process of the detonation wave impacting the metal lattice. Traditional depth gauge measurements can only obtain a single extreme point on the Z-axis, losing most of the morphological information and thus being unusable for analyzing the uniformity of energy distribution. Therefore, to facilitate subsequent analysis of energy distribution, high-precision 3D reconstruction techniques are needed to convert the physical steel block indentation into a digital matrix, providing comprehensive data support for subsequent analysis.

[0055] Preferably, as an example, the relative height field matrix is ​​obtained, including:

[0056] First, using a data acquisition device consisting of four ring light sources spaced 90 degrees apart and a vertical industrial camera, the light sources are controlled to illuminate in stages to capture images of the steel block indentation, obtaining four grayscale images under different lighting angles.

[0057] Next, specular highlight detection is performed on these four grayscale images. If the grayscale value of a pixel exceeds the camera's dynamic range, it is identified as a specular saturation point and discarded, retaining only the unsaturated pixel data. For example, the camera's dynamic range is taken as... .

[0058] Understandably, this step can effectively eliminate the interference of strong reflections from the metal surface on the calculation of the normal vector.

[0059] Then, based on the Lambertian reflection model, a system of equations was established using the remaining image data to calculate the normal vector field of each pixel on the surface of the steel block indentation.

[0060] Finally, the normal vector field is reconstructed by global integration using the global least squares integral method of the Poisson equation, and the output is the relative height field matrix that can characterize the three-dimensional micro-morphology of the steel block indentation.

[0061] Optional, as an example, obtain the relative height field matrix, including:

[0062] Control the single-axis motion module equipped with a high-precision line laser profilometer to sweep across the surface of the steel block indentation at a constant speed, and directly collect the point cloud depth data of the surface.

[0063] Subsequently, the point cloud depth data is interpolated by gridding to map it onto a regular two-dimensional grid, and a Gaussian smoothing filter algorithm is used to remove high-frequency noise caused by laser speckle effect, finally outputting a smooth and continuous relative height field matrix.

[0064] It should be noted that this method directly measures geometric distance and is suitable for scenarios that are not sensitive to lighting models.

[0065] S2: Determine the detonation impact center of the relative height field matrix, establish a polar coordinate system with the detonation impact center as the origin, convert the relative height field matrix into a polar coordinate energy flow field, and calculate the gradient data along the radial direction.

[0066] It should be noted that in order to analyze whether there is energy skew in the detonating cord, the energy distribution of the steel block indentation needs to be analyzed. The height field matrix describing the surface information of the steel block indentation is constructed in a Cartesian coordinate system, but the row and column structure in the Cartesian coordinate system cannot well describe the symmetry information of the indentation. The polar coordinate system can better reflect and describe the symmetry information of the indentation, so a coordinate system transformation is required for the height field matrix.

[0067] Preferably, as an example, the detonation impact center of the relative height field matrix is ​​determined, a polar coordinate system is established with the detonation impact center as the origin, the relative height field matrix is ​​converted into a polar coordinate energy flow field, and the gradient data along the radial direction is calculated, including:

[0068] First, calculate the geometric moments of the relative height field matrix, and determine the geometric centroid of the matrix by the ratio of the first moment to the zeroth moment, marking it as the detonation impact center.

[0069] Next, taking the detonation impact center as the pole, the cell positions in the relative height field matrix are converted into radii and angles. The relative height field matrix in Cartesian coordinates is resampled into the polar coordinate grid using a bilinear interpolation algorithm to obtain the polar coordinate energy flow field.

[0070] Finally, the polar coordinate energy flow field is differentially calculated along the radial axis to obtain the radial gradient data by calculating the rate of change of each point along the radial direction.

[0071] Figure 2 This is a schematic diagram of the polar coordinate energy flow field. The left side of the diagram represents the polar coordinate energy flow field of a qualified product. This energy flow field appears as uniform, continuous blue concentric circles, indicating that the radial gradient decreases only with the radius, while maintaining a high degree of consistency in the circumferential direction. The left side of the diagram represents the polar coordinate energy flow field of a false qualified product. This energy flow field shows a significant red eccentric distortion region, visually reflecting the accumulation and deflection of explosion energy in a specific direction.

[0072] Optionally, as an example, the detonation impact center of the relative height field matrix is ​​determined, a polar coordinate system is established with the detonation impact center as the origin, the relative height field matrix is ​​converted into a polar coordinate energy flow field, and the radial gradient data is calculated, including:

[0073] First, traverse the entire relative height field matrix, search for the cell with the smallest value, and determine its coordinates as the detonation impact center.

[0074] Subsequently, a polar coordinate mapping relationship was established based on the detonation impact center, and the relative height field matrix was converted into a polar coordinate energy flow field.

[0075] Finally, when calculating the gradient, the Sobel operator is first used to calculate the gradient in the X and Y directions in the original Cartesian coordinate system, and then the gradient is projected onto the radial direction to output the radial gradient data.

[0076] Understandably, this method ensures that the detonation impact center is always located at the point of most severe physical damage, making it more robust to samples with severe eccentricity.

[0077] S3: Based on the polar coordinate energy flow field and the gradient data, construct and calculate the asymmetric penalty coefficient. The asymmetric penalty coefficient characterizes the energy distribution difference of the explosion shock wave in the circumferential direction, and the calculation process includes weighting based on the distance attenuation law.

[0078] It should be noted that if the detonating cord has a uniform energy output, the degree of damage at the same radius on the steel block should be consistent. The deformation at the edge of the steel block is often affected by stress wave reflection, thus the edge of the steel block may exhibit a certain degree of fluctuation. Therefore, this can be used to analyze energy distribution indicators.

[0079] Preferably, as an example, based on the polar coordinate energy flow field and the gradient data, an asymmetric penalty coefficient is constructed and calculated. This asymmetric penalty coefficient characterizes the energy distribution difference of the explosion shock wave in the circumferential direction, and the calculation process includes weighting based on distance attenuation laws, including:

[0080] First, calculate the circumferential energy variation factor, which specifically satisfies the following relationship:

[0081]

[0082] in, The radius of the current calculation point. The attenuation constant is used to characterize the radiation properties at the center of the shock wave. For radius Circumferential energy variation factor at the location; For radius Angle is The radial gradient value at that location; For radius The average radial gradient at all angles; Pi () is an exponential function with the natural constant as its base;

[0083] It is understandable that the square root part in the relation calculates the standard deviation of the gradient, which represents the degree of non-uniformity of energy distribution at the same radius location; For distance decay weights, It is the attenuation constant, for example. Data is taken from one-third to one-half of the maximum radius of the effective computational domain of the relative height field matrix. Distance attenuation weighting can achieve this: as the radius... As the weight increases, the weight value decreases rapidly, thereby suppressing the interference of edge deformation caused by stress wave reflection far from the center.

[0084] Figure 3 This diagram illustrates the circumferential energy variation factor. The horizontal axis represents the distance from the detonation impact center, and the vertical axis represents the circumferential energy variation factor. The green curve shows the time-series data of the circumferential energy variation factor for the qualified sample. This curve shows that the circumferential energy variation factor remains consistently low, close to zero, indicating extremely stable energy output. The red dashed line shows the time-series data of the circumferential energy variation factor for the pseudo-qualified sample. This curve exhibits a large peak near the detonation impact center, and its value gradually decreases with increasing distance from the detonation impact center, thus accurately reflecting the attenuation characteristics weighted by the central radiation attenuation model.

[0085] Then, the asymmetric penalty coefficient is calculated, specifically satisfying the following relation:

[0086]

[0087] in, This is the asymmetric penalty coefficient; The radius of the effective computational domain of the relative height field matrix; It is the circumferential energy variation factor; radius The average depth at that location; This represents the maximum depth of the dent in the steel block.

[0088] It is understandable that the integration operation accumulates the circumferential energy variation factors at all radii, thus obtaining the overall skewness of the energy distribution. As a weighted term, this term further enhances the weight of the deepest depression, thus placing more trust in the calculation results at the center and effectively suppressing the interference of edge deformation caused by stress wave reflection far from the center. The larger the value, the more severe the energy deviation of the detonating cord.

[0089] Optionally, as an example, based on the polar coordinate energy flow field and the gradient data, an asymmetric penalty coefficient is constructed and calculated. This asymmetric penalty coefficient characterizes the energy distribution difference of the explosion shock wave in the circumferential direction, and the calculation process includes weighting based on distance attenuation laws, including:

[0090] Extracting each radius in the polar coordinate energy flow field The corresponding circumferential height sequence.

[0091] A Fast Fourier Transform (FFT) is performed on the circumferential height sequence to obtain the spectral coefficients. The sum of squared magnitudes of the non-DC harmonic components is extracted and denoted as the radius. The distortion energy at this point represents the distortion energy of the waveform, such as ellipticization and triangulation distortion.

[0092] Multiply the distortion energy by the same distance attenuation weight. The calculation results obtained for all radii are summed, and then the sum is normalized to obtain the asymmetric penalty coefficient.

[0093] Understandably, this method reflects the irregularity of the shape from the perspective of the frequency domain, and it has a better recognition ability for scenes with periodic information.

[0094] S4: Obtain the maximum depth of the steel block indentation, use the asymmetric penalty coefficient to perform nonlinear correction on the maximum depth to obtain the effective power value, and determine the qualification of the detonating cord output power based on the effective power value.

[0095] It should be noted that the traditional method measures the detonating cord's output power using depth information; however, energy skew reduces the actual effectiveness of the detonating cord. Therefore, in order to accurately measure the detonating cord's output power, the depth index needs to be corrected using the previously calculated index that reflects energy skew.

[0096] Preferably, as an example, the maximum depth of the steel block indentation is obtained, and the maximum depth is nonlinearly corrected using the asymmetric penalty coefficient to obtain an effective power value. The qualification of the detonating cord's output power is then determined based on the effective power value, including:

[0097] First, the maximum depth of the steel block indentation is corrected using a global asymmetric penalty coefficient to obtain the corrected effective power value. The specific relationship is as follows:

[0098]

[0099] in, This is the corrected effective power value; This represents the maximum depth of the dent in the steel block. This is the asymmetric penalty coefficient; This is the instability sensitivity coefficient, for example, for 45# steel. Take 3.5; It is the hyperbolic tangent function.

[0100] Understandably, when the detonating cord is of good quality, Approaching ,but Approaching The correction factor approaches ,at this time Approaching The product will maintain its original evaluation criteria.

[0101] When the detonating cord has energy deflection Increase Approaching 1, the correction factor approaches 1. At this point, the calculation is... A value that is too small is considered seriously unqualified.

[0102] Then, the calculated corrected effective power value is compared with the standard threshold. If the corrected effective power value is greater than the standard threshold, the output is "qualified"; if the corrected effective power value is not greater than the standard threshold, the output is "unqualified". For example, the standard threshold is 0.3mm.

[0103] Figure 4 This is a schematic diagram of the final evaluation results. The left side of the image shows the evaluation results obtained by the traditional method. The image shows that the physical depth of the counterfeit product is 0.520mm, which is even larger than the 0.499mm of the qualified product, and both exceed the national standard threshold of 0.3mm. This proves that using existing technology can easily misclassify defective products as superior products.

[0104] The left side of the image shows the evaluation results obtained by this invention. The image shows that the effective power value of the qualified sample is 0.451mm, far exceeding the 0.3mm threshold, proving that this invention can accurately detect qualified samples without false positives. The effective power value of the false qualified sample is 0.142mm, far below the 0.3mm threshold. The red arrow in the image clearly indicates the interception effect, proving that this invention, through an asymmetric penalty mechanism, successfully filters out false qualified samples that meet the depth requirement but have skewed energy.

[0105] This invention also discloses a detonating cord output power assessment system, including a processor and a memory. The memory stores computer program instructions, which, when executed by the processor, implement a detonating cord output power assessment method according to the present invention.

[0106] The system also includes other components well known to those skilled in the art, such as communication buses and communication interfaces, the settings and functions of which are known in the art and will not be described in detail here.

[0107] In this invention, the aforementioned memory can be any tangible medium containing or storing a program that can be used or combined with an instruction execution system, apparatus, or device. For example, a computer-readable storage medium can be any suitable magnetic or magneto-optical storage medium, such as resistive random access memory (DRAM), dynamic random access memory (DRAM), static random access memory (SRAM), enhanced dynamic random access memory (DRAM), high-bandwidth memory, hybrid memory cube, etc., or any other medium that can be used to store desired information and can be accessed by an application, module, or both. Any such computer storage medium can be part of a device or accessible to or connected to a device.

Claims

1. A method for evaluating the output power of a detonating cord, characterized in that, The method comprises the steps of: obtaining a relative height field matrix of the steel block indentation, the relative height field matrix representing a three-dimensional micro-topography of the steel block indentation; determining a detonation impact center of the relative height field matrix, establishing a polar coordinate system with the detonation impact center as the origin, converting the relative height field matrix into a polar coordinate energy flow field, and calculating gradient data along the radial direction; based on the polar coordinate energy flow field and the gradient data, constructing and calculating an asymmetry penalty coefficient, the asymmetry penalty coefficient representing the energy distribution difference of the explosion shock wave in the circumferential direction, and the calculation process includes weighted processing based on the distance decay law, specifically including: performing fast Fourier transform on the circumferential height sequence at the same radius in the polar coordinate energy flow field; extracting the energy sum of the transformed non-direct current harmonic components and multiplying it by the decay weight increasing with the radius to obtain the frequency domain heterogeneity factor; performing weighted accumulation and normalization processing on the frequency domain heterogeneity factors at all radii to obtain the asymmetry penalty coefficient; obtaining the maximum depth of the steel block indentation, and nonlinearly correcting the maximum depth by using the asymmetry penalty coefficient to obtain an effective power value, specifically including: ； wherein, represents the corrected effective power value; represents the maximum depth of the dent of the steel block measured; ( ) represents the hyperbolic tangent function; represents the instability sensitivity coefficient related to the hardness of the steel block material; represents the asymmetry penalty coefficient determining the eligibility of the output power of the detonating cord according to the effective power value.

2. The method for evaluating the output power of a detonating cord according to claim 1, wherein, One of the methods for obtaining the relative height field matrix of the steel block indentation specifically includes: controlling the four-azimuth ring light source to light up at different times, and using a vertical industrial camera to collect multiple gray-scale images of the steel block indentation; detecting the pixel gray-scale values in the gray-scale images, and removing the mirror reflection saturation points with gray-scale values exceeding the camera dynamic range threshold; using the remaining unsaturated image data to calculate the normal vector field of the steel block indentation surface based on the Lambertian reflection model; integrating and reconstructing the normal vector field to obtain the relative height field matrix.

3. The method of claim 1, wherein, Another method for obtaining the relative height field matrix of the steel block indentation specifically includes: controlling the single-axis motion module carrying the line laser profiler to uniformly scan the surface of the steel block indentation to collect point cloud depth data; performing gridding interpolation and Gaussian smoothing filter processing on the point cloud depth data to remove speckle noise and generate the relative height field matrix.

4. The method of claim 1, wherein, The method for determining the detonation impact center of the relative height field matrix specifically includes: calculating the geometric moments of the relative height field matrix; determining the geometric centroid of the relative height field matrix by using the ratio of the first moment to the zero moment, and marking the geometric centroid as the detonation impact center.

5. The method of claim 1, wherein, Another method for constructing and calculating the asymmetry penalty coefficient specifically includes: ； ； wherein, represents the circumferential energy anomaly factor at the polar coordinate system under the radius ; represents the radial gradient value at the radius and the angle ; represents the average radial gradient at all angles at the radius ; is the constant pi; () is the exponential function with the natural constant as the base; is the radius of the current calculation point; is the decay constant representing the central radiation characteristic of the shock wave; represents the asymmetry penalty coefficient; represents the maximum radius of the effective calculation domain of the relative height field matrix; represents the average depth at the radius ; represents the maximum depth of the dent of the steel block.

6. The method of claim 5, wherein, The value range of the attenuation constant is set to be between one-third and one-half of the radius of the steel block indentation.

7. The method of claim 1, wherein, The value of the instability sensitivity coefficient is positively correlated with the hardness of the steel block material, and for a 45 steel material, the value range of the instability sensitivity coefficient is 3.0 to 4.

0.

8. A system for assessing the output power of a detonating cord, characterized in that The method comprises: a processor and a memory, the memory storing computer program instructions, when the computer program instructions are executed by the processor, a detonating cord output power evaluation method according to any one of claims 1-7 is realized.

Images