Concrete crack detection method based on quantum detection and Hamiltonian modeling

Abstract

The invention discloses a concrete crack detection method based on quantum detection and Hamiltonian modeling. The method comprises the following steps: collecting detection data in a quantum beam scanning detection process; calculating to obtain a crack projection size and a crack inclination angle according to the collected detection data; substituting the crack inclination angle into a quantum-geometry joint formula to calculate the depth from the crack to the detection reference surface; wherein the quantum-geometry joint formula is represented as follows: the depth from the crack to the detection reference surface is equal to the product of the crack depth in the normal vector direction of the plane where the crack is located and the cosine value of the crack inclination angle; the crack depth is obtained through inversion calculation by introducing a quantum model of quantum mechanics Hamiltonian. By introducing quantum mechanics Hamiltonian quantity modeling, accurate inversion driven by a physical mechanism is realized, so that the accuracy and reliability of concrete crack detection are improved.

Description

Technical Field

[0001] The present application belongs to the technical field of non-destructive testing of concrete structures, and specifically relates to a concrete crack detection method based on quantum detection and Hamiltonian modeling. Background Art

[0002] Traditional concrete testing technologies (such as ultrasonic method and surface wave method) have the following defects:

[0003] (1) Limited detection depth: The maximum detection depth of the ultrasonic method usually does not exceed 0.5m, which cannot meet the needs of large-volume concrete structures.

[0004] (2) Destructive testing: The surface wave method requires an impact source, which may cause secondary damage to the structure.

[0005] (3) Manual interpretation error: Relying on experience to judge the fracture morphology makes it difficult to quantify complex fracture parameters.

[0006] Although existing quantum detection technology introduces quantum particle beams, it has the following shortcomings:

[0007] (1) The crack depth is calculated using an empirical formula, and the error exceeds 20% when the concrete density fluctuation is greater than 5%.

[0008] (2) The coupling effect between the magnetic field gradient and the fracture inclination angle is not considered, resulting in the failure of depth inversion of inclined fractures (θ>30°). Summary of the Invention

[0009] In order to overcome the shortcomings of the existing technology, this application proposes a concrete crack detection method based on quantum detection and Hamiltonian modeling. This application introduces quantum mechanics Hamiltonian modeling to achieve precise inversion driven by physical mechanisms, thereby improving the accuracy and reliability of concrete crack detection.

[0010] This application is implemented through the following technical solutions:

[0011] A concrete crack detection method based on quantum detection and Hamiltonian modeling, comprising:

[0012] Collect detection data during the quantum beam scanning detection process;

[0013] Calculating the fracture projection size and fracture inclination angle based on the collected detection data;

[0014] The crack inclination angle is substituted into the quantum-geometric joint formula to calculate the depth of the crack from the detection reference plane; wherein the quantum-geometric joint formula is expressed as: the depth of the crack from the detection reference plane is equal to the product of the crack depth in the direction of the normal vector of the plane where the crack is located and the cosine value of the crack inclination angle; the crack depth is obtained by solving a quantum model that introduces the quantum mechanical Hamiltonian.

[0015] In some embodiments, the process of obtaining the crack depth in the direction of the normal vector of the plane where the crack is located includes:

[0016] Calculate the energy difference of quantum beams;

[0017] Calculating the attenuation coefficients of the quantum beam in the crack region and the intact region based on the energy difference of the quantum beam and the material parameters;

[0018] The quantum wave beam in the crack region and the attenuation coefficient of the intact region are substituted into the energy eigenvalue equation for solution to obtain the crack depth; the crack depth is obtained by inverse calculation of a quantum model introducing the quantum mechanics Hamiltonian.

[0019] In some embodiments, the process of constructing the energy eigenvalue equation includes:

[0020] Constructing a total Hamiltonian, wherein the total Hamiltonian is the sum of an electromagnetic field kinetic energy term, a potential field correction term, a spin-magnetic field term, and a spin-orbit term;

[0021] Expanding the total Hamiltonian into the sum of a kinetic energy term and a potential energy term;

[0022] Simplifying the kinetic energy term in the total Hamiltonian to obtain a first simplified total Hamiltonian;

[0023] Simplifying the potential energy term in the total Hamiltonian to obtain a second simplified total Hamiltonian;

[0024] The second simplified total Hamiltonian is converted from three dimensions to one dimension by coordinate system transformation to obtain the final simplified total Hamiltonian;

[0025] Based on the final simplified total Hamiltonian, the one-dimensional steady-state Schrödinger equation is constructed and solved in different regions, finally obtaining the energy eigenvalue equation and the regional analytical solution form.

[0026] In some embodiments, the energy eigenvalue equation is expressed as:

[0027]

[0028] The analytical solution is in the form of:

[0029]

[0030] Where d is the vertical distance between the crack plane and the detection reference plane; L z' is the crack depth in the direction of the normal vector z' of the crack plane; κ is the attenuation coefficient of the complete area; k z' is the quantum wave number in the crack region;

[0031] in,

[0032]

[0033] Where V0 is the initial potential energy of the crack-free concrete region in the external potential field; ΔV is the volume change caused by the crack region in the external potential field; E is the wave number energy; m is the particle mass; is the reduced Planck constant; ΔV' is the potential energy change caused by the crack area in the new coordinate system.

[0034] In some embodiments, the coordinate system transformation process includes:

[0035] The original coordinate system xoy is the detection reference plane, the z axis is the beam emission direction; the angle between the crack plane and the detection reference plane is the crack inclination; the azimuth of the crack plane within the detection reference plane is the crack azimuth;

[0036] The x′oy′ plane of the new coordinate system is the crack plane, and the z′ axis is the direction of the normal vector of the crack plane;

[0037] According to the trigonometric function relationship between the fracture dip angle and the fracture azimuth angle, the transformation matrix is ​​obtained;

[0038] According to the transformation matrix, the three-dimensional kinetic energy term is converted into a one-dimensional kinetic energy term.

[0039] In some embodiments, further comprising:

[0040] The crack width and length are calculated based on the crack projection size, potential energy correction term, and calibration parameters.

[0041] In some embodiments, further comprising:

[0042] Dynamically correct the detection range boundaries.

[0043] On the other hand, this application also proposes a concrete crack detection device based on quantum detection and Hamiltonian modeling, comprising:

[0044] Data acquisition module, used to collect detection data during the quantum beam scanning detection process;

[0045] A calculation module, for calculating the crack projection size and crack inclination angle based on the collected detection data;

[0046] And, a quantum geometric inversion module is used to substitute the crack inclination angle into the quantum-geometric joint formula to calculate the depth of the crack to the detection reference plane; wherein, the quantum-geometric joint formula is expressed as: the depth of the crack to the detection reference plane is equal to the product of the crack depth in the direction of the normal vector of the plane where the crack is located and the cosine value of the crack inclination angle; the crack depth is obtained by inversion calculation of a quantum model introducing the quantum mechanical Hamiltonian.

[0047] In some embodiments, the quantum geometric inversion module includes:

[0048] A crack depth solving unit is used to solve the crack depth in the normal vector direction of the crack plane;

[0049] And, an actual depth solving unit is used to substitute the crack depth and crack inclination into the quantum-geometric joint formula to calculate the depth from the crack to the detection reference plane.

[0050] In some embodiments, further comprising:

[0051] The crack parameter calculation module calculates the crack width and length based on the crack projection size, potential energy correction term, and calibration parameters;

[0052] And / or, a dynamic correction module, used to dynamically correct the detection range boundary.

[0053] This application proposes a concrete crack detection method based on quantum detection and Hamiltonian modeling. This method introduces the quantum mechanical Hamiltonian to implement quantum modeling and combines the quantum model with a geometric model to achieve precise detection of concrete cracks. Furthermore, this application utilizes dynamic boundary correction and multi-beam collaboration to adapt to complex field environments, significantly improving detection efficiency. This application has a wide range of applications, including in hydropower, transportation, and other fields. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] The drawings described herein are used to provide a further understanding of the embodiments of the present application, constitute a part of the present application, and do not constitute a limitation of the embodiments of the present application. In the drawings:

[0055] Figure 1 A schematic flow chart of the concrete crack detection method proposed in an embodiment of the present application;

[0056] Figure 2 Schematic diagram of the quantum beam scanning detection plane;

[0057] Figure 3 This is a schematic diagram for calculating the geometric inclination of the crack;

[0058] Figure 4 Schematic diagram of crack geometry;

[0059] Figure 5 This is a block diagram of the principle of the concrete crack detection device proposed in the embodiment of the present application;

[0060] Reference numerals and corresponding component names:

[0061] 200 - concrete crack detection device, 201 - data acquisition module, 202 - calculation module, 203 - quantum geometry inversion module, 204 - crack parameter calculation module, 205 - dynamic correction module, 301 - crack depth solution unit, 302 - actual depth solution unit. DETAILED DESCRIPTION

[0062] Hereinafter, the terms "include" or "may include" as used in various embodiments of the present application indicate the presence of an invented function, operation, or element, and do not limit the addition of one or more functions, operations, or elements. In addition, as used in various embodiments of the present application, the terms "include," "have," and their cognates are intended only to indicate specific features, numbers, steps, operations, elements, components, or combinations of the foregoing, and should not be understood as excluding the presence of one or more other features, numbers, steps, operations, elements, components, or combinations of the foregoing or the possibility of adding one or more features, numbers, steps, operations, elements, components, or combinations of the foregoing.

[0063] In various embodiments of the present application, the expression "or" or "at least one of A or / and B" includes any or all combinations of the words listed simultaneously. For example, the expression "A or B" or "at least one of A or / and B" may include A, may include B, or may include both A and B.

[0064] The expressions (such as "first", "second", etc.) used in the various embodiments of the present application may modify the various constituent elements in the various embodiments, but may not limit the corresponding constituent elements. For example, the above expressions do not limit the order and / or importance of the elements. The above expressions are only used to distinguish one element from other elements. For example, a first user device and a second user device indicate different user devices, although both are user devices. For example, without departing from the scope of the various embodiments of the present application, a first element may be referred to as a second element, and similarly, a second element may also be referred to as a first element.

[0065] It should be noted that when a component is described as being “connected” to another component, the first component may be directly connected to the second component, and a third component may be “connected” between the first and second components. Conversely, when a component is described as being “directly connected” to another component, it can be understood that there is no third component between the first and second components.

[0066] The terms used in the various embodiments of the application are only used to describe the purpose of specific embodiments and are not intended to limit the various embodiments of the application. As used herein, the singular form is intended to also include the plural form, unless the context clearly indicates otherwise. Unless otherwise limited, all terms used here (including technical terms and scientific terms) have the same meaning as the meaning generally understood by those of ordinary skill in the art of the application. The terms (such as the terms defined in the dictionary generally used) will be interpreted as having the same meaning as the contextual meaning in the relevant technical field and will not be interpreted as having idealized meaning or too formal meaning, unless clearly defined in the various embodiments of the application.

[0067] In order to make the objectives, technical solutions and advantages of this application more clear, the present application is further described in detail below in conjunction with examples and drawings. The schematic implementation methods of this application and their descriptions are only used to explain this application and are not intended to limit this application.

[0068] Example:

[0069] Traditional concrete crack detection technology suffers from limitations such as limited detection depth, destructive testing, and manual interpretation errors. Existing quantum detection technologies use empirical formulas to calculate crack depth, resulting in large errors. They also fail to consider the coupling effect between magnetic field gradient and crack inclination, leading to ineffective crack depth inversion. To address these issues, this embodiment proposes a concrete crack detection method based on quantum detection and Hamiltonian modeling.

[0070] like Figure 1 As shown, the concrete crack detection method proposed in this embodiment includes the following steps:

[0071] Step 110: collecting detection data during the quantum beam scanning detection process.

[0072] In this step 110, multiple cylindrical quantum beams are generated by the device and emitted perpendicularly to the concrete surface; scanning is performed along the first direction (X axis) and the second direction (Y axis), and the electrical signal mutation positions S1, S2, S3, and S4 are recorded, such as Figure 2 and Figure 3 shown.

[0073] The quantum beam energy E is obtained by measuring the quantum host, and the potential energy correction V is obtained con .

[0074] like Figure 3As shown, the crack depths measured by the quantum host are h1=AA', h2=BF, h3=CI, and h4=DE, respectively. Among them, quadrilateral A'FIE is the crack plane, quadrilateral ABCD is the projection of the crack on the detection reference plane, ∠IGC'=θ, the crack width is W=A'E, and the crack length is L=A'F.

[0075] Step 120: Calculate the crack projection size and crack inclination angle based on the detection data.

[0076] Calculate the projection length L in two directions according to the position of the electrical signal mutation x and L y , as shown in formula (1):

[0077] L x =|S3-S1|, L y =|S4-S2| (1)

[0078] It should be noted that the projection length is calculated by taking the average of at least three sets of data, that is, according to the above calculation formula, multiple projection length values ​​are calculated using at least three sets of data, and then the average is taken as the final projection length value.

[0079] The crack inclination angle θ is calculated from the multi-directional projection data, and the calculation formula (Equation (2)) is as follows:

[0080]

[0081] Step 130: Substitute the crack inclination angle into the quantum-geometric joint formula to calculate the distance from the crack to the detection reference plane.

[0082] The crack inclination angle is substituted into the quantum-geometric joint formula (Equation (3)) to calculate the depth h from the crack to the detection reference plane:

[0083]

[0084] Furthermore, the derivation process of the quantum-geometric joint formula is as follows:

[0085] First, combine the quantum and geometric equations and transform the quantum depth L z' Equivalent to the geometric depth d (L z' =d), substitute into the geometric equation h = L z' cosθ, we get the equation shown in formula (4):

[0086]

[0087] Among them, L z' It represents the crack depth (i.e. quantum depth) on the z′ axis in the direction of the crack plane normal vector.

[0088] The quantum wave number k in the crack region z' The attenuation coefficient κ of the complete area is replaced by the function of energy E, as shown in formula (5):

[0089]

[0090] Through algebraic operations, the joint equation can be simplified to the equation shown in formula (6):

[0091]

[0092] Furthermore, the crack depth solution process is:

[0093] Calculate the energy difference of the quantum beam ΔE=E-(V0-ΔV').

[0094] Calculate the quantum wave number k in the crack region based on the energy difference of the quantum beam and the material parameters (m, V0, △V) z' and the attenuation coefficient κ of the complete area.

[0095] Solve the energy eigenvalue equation by numerical methods (such as Newton iteration method):

[0096] Get the crack depth L z’ .

[0097] Furthermore, the energy eigenvalue equation is constructed as follows:

[0098] This embodiment simplifies the Hamiltonian to contain only L by using the principle of gauge invariance. z' The observable form of the potential energy gradient in the crack direction (i.e., the direction of the normal vector of the plane where the crack is located).

[0099] (1) Construct the total Hamiltonian, as shown in formula (7):

[0100]

[0101] in, is the momentum operator, and the expression is is the reduced Planck constant, is the gradient operator; q is the charge of the particle, q = -e (e = 1.6 × 10 -19 C); is the scalar potential; is the vector potential (magnetic vector potential), satisfying ( is the magnetic induction intensity); m is the particle mass; is the spatial position vector; is the Pauli spin matrix; is the electric field strength; C is the speed of light in vacuum (3×10 8 m / s).

[0102] The potential energy correction term in the external field of view (such as gravity field, electric field) is shown in formula (8):

[0103]

[0104] in, is the potential energy from the coordinate zero point to the location of the concrete in the external field of view; V0 is the initial state potential energy of the crack-free concrete area in the external potential field; ΔV is the volume change caused by the crack area in the external potential field; ρ v , ρ v0 are the potential energy per unit volume and the initial value respectively; is the crack surface energy, which is related to the crack width, is the surface energy attenuation function, σ0 is the surface energy of intact concrete, W0 is 1 / 2 of the average diameter of aggregate, and W is the crack length; is the density gradient term around the crack, which represents the material heterogeneity, L c is the average crack width; α, β, γ are the correlation coefficients with the concrete material properties (such as porosity, elastic modulus, etc.), which are obtained through experimental calibration; in addition, according to the electric field The spin-orbit coupling term can be expressed as shown in equation (9):

[0105]

[0106] (2) Based on the total Hamiltonian, construct the time-dependent Schrödinger equation and solve it.

[0107] The time-dependent Schrödinger equation is expressed as shown in equation (10):

[0108]

[0109] The solution process includes:

[0110] First, the total Hamiltonian is expanded into the sum of the kinetic energy term and the potential energy term, as shown in Equation (11):

[0111]

[0112] Convert equation (11) to equation (12):

[0113]

[0114] The three-dimensional momentum operator is the equation shown in formula (13):

[0115]

[0116] In order to simplify the solution process, if the potential field does not explicitly contain time, the variables can be separated The time-dependent Schrödinger equation is converted into the steady-state Schrödinger equation, as shown in Equation (14):

[0117]

[0118] Where E is the energy eigenvalue of the system.

[0119] The steady-state solution consists of two parts: the time part and the spatial part (the eigenvalue equation needs to be solved according to the specific form of the potential field).

[0120] Consider the spatial morphology of the crack: When the crack plane not only forms an inclination angle θ with the xoy plane (i.e., the detection reference plane), but also has an azimuth angle When , quantum mechanics formulas need to adapt to three-dimensional geometric relationships through coordinate system rotation and parameter projection correction. The following is a detailed analysis:

[0121] Define the crack plane direction: inclination angle θ: the angle between the crack plane and the xoy plane (the angle between the normal vector and the z axis); azimuth The azimuth angle of the crack plane in the xoy plane (the angle between the projection of the normal vector in the xoy plane and the x-axis).

[0122] Coordinate system transformation: establish a new coordinate system (x′, y′, z′), where the z′ axis is along the direction of the normal vector of the crack plane; the x′ axis is along the projection direction of the crack in the xoy plane (i.e. direction); 3) the y′ axis is perpendicular to the x′ and z′ axes, forming a right-handed coordinate system.

[0123] There is an angle θ between the actual detection direction (z axis of the original coordinate system) and the normal vector of the crack plane (z′ axis of the new coordinate system). The calculated result L z' Project it to the z-axis to obtain the true crack depth h.

[0124] After establishing a new coordinate system, such as Figure 4 As shown, the crack depth h is related to L in the new coordinate system z' The relationship is shown in formula (15):

[0125] h=L z' cosθ (15)

[0126] By using coordinate transformation, the three-dimensional calculation of the crack with inclined emission direction of the beam instrument is converted into a simple three-dimensional model of the crack with vertical emission direction, thereby simplifying the calculation.

[0127] Simplification of kinetic energy terms: When the vector potential amplitude of the quantum beam satisfies When , the squared vector potential term can be ignored. That is, the squared vector potential term

[0128] The vector potential slowly varying approximation requires that the second-order derivative of the vector potential be negligible That is, the spatial variation scale is much larger than the particle wavelength (λ 粒子 ≤L 矢势 ). At this time, the product of the momentum operator and the vector potential approximately satisfies

[0129] Scalar potential electrostatic field approximation ( Not obviously time-varying), when the crack propagation speed (such as v<10 3 m / s), the scalar potential can be regarded as a quasi-static field, and the time derivative term Can be ignored.

[0130]

[0131] Based on the above simplification, the simplified total Hamiltonian is obtained, as shown in formula (16):

[0132]

[0133] Establishment of potential energy model:

[0134] Assume that the potential energy is uniformly distributed along the x′oy′ plane of the new coordinate system and ignore the lateral variation. x′ , L y′ : They represent the projection lengths of the cracks along the two detection plane directions x′ and y′ in the x′oy′ plane, respectively. The normal vector direction of the crack plane x′oy′ is the z′ axis. L z′ Indicates the length on the z' axis.

[0135] When the crack lateral size (L x′ ,L y′ )>>λ beam, transverse quantum effect can be ignored, potential energy correction value of concrete crack It can be simplified to a one-dimensional distribution in the z′ direction:

[0136] (1) Crack area d≤z′≤d+L z′ :

[0137] make Represents the potential energy correction due to cracks.

[0138] The potential energy is reduced to the equation shown in formula (17):

[0139]

[0140] (2) Complete area (other areas): the potential energy remains at V0.

[0141] (3) The total potential field is expressed as the equation shown in formula (18):

[0142]

[0143] (4) Concrete is a non-magnetic material (relative magnetic permeability μ r ≈1), the magnetic field strength B in the crack area is mainly provided by the external superconducting magnet, taking the typical parameter B=1T, then μ B B≈5.8×10 -5 eV; if the potential energy is corrected V con ~1eV, ​​satisfying (The spin-magnetic field energy is much smaller than the potential energy correction term), so the spin-magnetic field term is ignored.

[0144] (5) The speed of quantum beam particles (such as electrons or neutral atoms) v<<c (the speed of light), meeting the non-relativistic limit condition (β=v / c~10 -4 ); spin-orbit coupling term The magnitude is 10 -8 eV, kinetic energy term The order of magnitude is (~1eV), and the order of magnitude of the spin-orbit coupling term is much smaller than that of the kinetic energy term, so it is ignored.

[0145] (6) The spin polarization rate χ of the quantum beam remains spatially uniform during the detection process And the spin relaxation time is much longer than the wave function evolution time. Decoupled from the spatial coordinates, the spin-related terms in the Hamiltonian degenerate into constants and do not affect the spatial distribution of the wave function, so only the potential energy-dominant term is retained.

[0146] Therefore, based on the simplification of the potential energy term above, the simplified total Hatt density quantity is obtained, as shown in formula (19):

[0147]

[0148] Then the three-dimensional stationary Schrödinger equation is converted to one dimension. The conversion process is as follows:

[0149] (1) Assume that the cracks are uniformly distributed along the x′ and y′ directions, the crack ends are continuous in space, and only the momentum and potential energy change in the z′ direction.

[0150] (2) Transformation of momentum operator: The momentum operator must satisfy after coordinate transformation in, is the gradient operator in the new coordinate system.

[0151] (3) In the rectangular coordinate system, the momentum operator is derived from its differential form in the new coordinate system, or the commutation relation between the momentum operator and the position operator is used for correction.

[0152] The transformation matrix is ​​the equation shown in formula (20):

[0153]

[0154] Under the transformation matrix R, the momentum operator is transformed into For rigid rotation (R is a constant matrix), therefore

[0155] It should be noted that: in the new coordinate system, the vector potential and scalar potential It is necessary to transform synchronously to maintain the invariance of the electromagnetic field. For a uniform magnetic field, the vector potential Scalar potential Transformation, where is the gradient operator in the new coordinate system.

[0156] (4) Assumptions and σ s (z') are all constants.

[0157] (5) Assuming electron beam, low field approximation, three-dimensional electromagnetic field terms Simplify and retain only the kinetic energy dominant term

[0158] (6) Potential energy term It changes only in the z′ direction (crack depth direction) and is uniformly distributed in the transverse direction, that is, there is no change in the x′ and y′ directions.

[0159] The total Hamiltonian is finally simplified to the equation shown in formula (21):

[0160]

[0161] The steady-state equation is shown in equation (22):

[0162]

[0163] Then the one-dimensional steady-state Schrödinger equation is solved. The solution process includes:

[0164] Assume that the crack plane is located at a vertical distance d below the detection reference plane (concrete surface), and the crack depth itself is L z′ The coordinate system is defined as follows:

[0165] z′=0: detection plane (concrete surface); crack area: d≤z′≤d+L z′ ; The detection reference plane, that is, the concrete surface is the xoy plane, and the actual extension plane of the length and width of the concrete crack is x′oy′.

[0166] Wave function is defined by region:

[0167] Region I (z′ < d): Intact concrete, potential energy V con = V0;

[0168] Region II (d ≤ z′ ≤ d + L z′ ): Crack region, potential energy V con = V0 - ΔV′;

[0169] Region III (z′ > d + L z′ ): Intact concrete, potential energy V con = V0.

[0170] For the intact regions (z′ < d or z′ > d + L z′ ):

[0171] Potential energy V con = V0, and the stationary - state equation simplifies to the equation shown in Eq. (23):

[0172]

[0173] When z′ < d, the wave function is in an exponentially - decaying form, as shown in Eq. (24):

[0174]

[0175] When z′ > d + L z′ , the wave function is in an exponentially - decaying form, as shown in Eq. (25):

[0176] <​​​​​​​​​​​​​​​​​​​​​​​​​​​​​​At the crack tip, z′ = d + L z′ , the wave function shown in Equation (29) is obtained:

[0185]

[0186] The derivation process of the analytical solution formula is as follows:

[0187] The wave function in the complete region (z′ < d), as shown in Equation (30):

[0188]

[0189] The wave function in the crack region (d ≤ z' ≤ d + L z' ), as shown in Equation (31):

[0190]

[0191] Where

[0192] The wave function in the complete region (z′ > d + L z′ ), as shown in Equation (32):

[0193]

[0194] Applying the boundary conditions:

[0195] At z′ = d: A = D, Aκ = Ck z' .

[0196] At z′ = d + L z′ :

[0197] From the starting point boundary conditions, we can obtain and A = D, and substitute them into the end boundary conditions to obtain the system of equations shown in Equations (33) and (34):

[0198]

[0199] Transform the above system of equations shown in Equations (33) and (34) (i.e., Equation (33) × (-κk z' ), Equation (34) × k z' ) to obtain the system of equations shown in Equations (35) and (36):

[0200]

[0201] 联立 the transformed system of equations shown in Equations (35) and (36) to obtain the equation shown in Equation (37):

[0202]

[0203] By shifting the terms in the above equation, we can get the equation shown in formula (38):

[0204] (k z' 2 -κ 2 )sin(k z' L z' )=2κk z' cos(k z' L z' ) (38)

[0205] Rewriting the above equation can obtain the energy eigenvalue equation, as shown in equation (39):

[0206]

[0207] in, represents the attenuation coefficient of the complete area;

[0208] represents the quantum wave number in the crack region;

[0209] The energy eigenvalue equation determines the energy E and the crack depth L z' The relationship between them.

[0210] Among them, the attenuation coefficient of the complete area is The quantum wave number in the crack region is

[0211] Combine The energy E can be solved numerically.

[0212] E can be a measured value or derived through the above formula, and is usually a measurable value.

[0213] Finally, L is solved by the energy eigenvalue equation z'

[0214] And the analytical solution is shown in formula (40):

[0215]

[0216] It can be seen from this that the crack depth L z' The potential energy correction value directly affects the energy level splitting; by measuring the energy ΔE=E-(V0-ΔV'), the crack depth L is solved by combining the energy eigenvalue equation. z' .

[0217] Then the distance h from the crack to the detection reference plane is calculated using the quantum-geometric joint formula.

[0218] Furthermore, the detection method proposed in this embodiment also includes:

[0219] The crack width and length are calculated based on the crack projection size, potential energy correction term, and calibration parameters.

[0220] The crack width directly affects the surface energy attenuation function, which is expressed as the equation shown in Equation (41):

[0221]

[0222] In the crack region, the potential energy correction term is expressed as shown in equation (42):

[0223] V con (z')=V0-ΔV' (42)

[0224] Substituting the surface energy attenuation function and sorting it out, the approximate formula is shown in formula (43):

[0225]

[0226] The crack width W is iteratively calculated through experiments and detection potential energy correction terms, combined with calibration parameters α, β, γ, and σ0.

[0227] The crack length is calculated based on the crack width: L=kW, where k represents the proportional coefficient of length to width, which can be determined through experimental calibration.

[0228] For example, crack samples with different aspect ratios (e.g., rectangular cracks with k = 1, 2, 5, and 10) can be prepared. The actual length and width of the cracks are then measured, and the k value for each sample is calculated. For each sample, the crack width W is inverted using quantum mechanics formulas. Assuming k = 1 in the inversion, the theoretical crack length is calculated. The error function is then established: k is optimized by nonlinear fitting (such as the least squares method) to minimize the error.

[0229] Furthermore, the actual length and width of the crack can be calculated by geometric correction with the actual length and width of the crack being the projected length to verify the crack length and width obtained by inversion calculation of the quantum mechanics formula, as shown in the following example: Figure 4 As shown, the calculation formula is shown in formula (44):

[0230]

[0231] The concrete crack detection method proposed in this embodiment introduces quantum mechanics Hamiltonian modeling into concrete detection for the first time, achieving accurate inversion driven by physical mechanisms.

[0232] Furthermore, the concrete crack detection method proposed in this embodiment also includes:

[0233] Dynamically correct the detection range boundaries.

[0234] The specific correction process is as follows:

[0235] First, obtain the boundary correction value S, as shown in formula (45):

[0236]

[0237] Where χ is the dimensionless spin polarization (0≤χ≤1); is the reduced Planck constant; is the electric field gradient; k is the proportionality constant, which can be determined by calibration experiments on standard concrete test blocks (crack width 0.01-1 mm) and satisfies k = 0.68 ± 0.05 (95% confidence interval); m is the electron mass; is the fundamental frequency of the quantum beam (in rad / s).

[0238] Then, based on the boundary correction amount S, the corrected detection range (S1-S, S2+S, S3-S, S4+S) is obtained.

[0239] This embodiment adapts to complex on-site environments through dynamic detection boundary correction and multi-beam collaboration, greatly improving detection efficiency.

[0240] In another embodiment, this embodiment also proposes a concrete crack detection device 200 based on quantum detection and Hamiltonian modeling, such as Figure 5 As shown, the detection device 200 includes:

[0241] The data acquisition module 201 is used to acquire detection data during the quantum beam scanning detection process. The details are as described in the above method and will not be repeated here.

[0242] The calculation module 202 calculates the fracture projection size and fracture inclination angle based on the collected data. The specific calculation process is as described in the above method and will not be repeated here.

[0243] Furthermore, the quantum geometry inversion module 203 is used to substitute the crack inclination angle into the quantum-geometry joint formula to calculate the depth of the crack to the detection reference plane. The specific calculation process is as described in the above method and will not be repeated here.

[0244] Furthermore, the detection device proposed in this embodiment also includes:

[0245] The crack parameter calculation module 204 calculates the crack width and crack length based on the crack projection size, potential energy correction term, and calibration parameters. The specific calculation process is as described in the above method and will not be repeated here.

[0246] Furthermore, the detection device proposed in this embodiment also includes:

[0247] The dynamic correction module 205 is used to dynamically correct the detection range boundary. The specific correction method is as described in the above method and will not be repeated here.

[0248] Furthermore, the quantum geometry inversion module 203 of this embodiment further includes:

[0249] The crack depth solving unit 301 is used to solve the crack depth in the normal vector direction of the plane where the crack is located. The specific solving process is as described in the above method and will not be repeated here.

[0250] Furthermore, the actual depth solving unit 302 substitutes the crack depth and inclination output by the crack depth solving unit 301 into the quantum-geometric joint formula to calculate the depth from the crack to the detection reference plane, which is the actual depth of the crack. The specific calculation process is as described in the above method and will not be repeated here.

[0251] It should be noted that, in the above embodiments, the description of each embodiment has its own focus. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant description of other embodiments.

[0252] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0253] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the steps in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0254] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1a process or multiple processes and / or boxes Figure 1 The function specified in one or more boxes.

[0255] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0256] The specific implementation methods described above further illustrate the purpose, technical solutions and beneficial effects of this application. It should be understood that the above description is only the specific implementation methods of this application and is not intended to limit the scope of protection of this application. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of this application should be included in the scope of protection of this application.

Claims

1. A concrete crack detection method based on quantum detection and Hamiltonian modeling, characterized in that: Collect detection data during the quantum beam scanning detection process; Calculating the fracture projection size and fracture inclination angle based on the collected detection data; The crack inclination angle is substituted into the quantum-geometric joint formula to calculate the depth of the crack from the detection reference plane; wherein the quantum-geometric joint formula is expressed as: the depth of the crack from the detection reference plane is equal to the product of the crack depth in the direction of the normal vector of the plane where the crack is located and the cosine value of the crack inclination angle; the crack depth is calculated by inversion of the quantum model introducing the quantum mechanical Hamiltonian.

2. A concrete crack detection method based on quantum detection and Hamiltonian modeling according to claim 1, characterized in that: The process of obtaining the crack depth in the direction of the normal vector of the plane where the crack is located includes: Calculate the energy difference of quantum beams; Calculating the attenuation coefficients of the quantum beam in the crack region and the intact region based on the energy difference of the quantum beam and the material parameters; The quantum wave beam in the crack region and the attenuation coefficient of the intact region are substituted into the energy eigenvalue equation and solved to obtain the crack depth.

3. A concrete crack detection method based on quantum detection and Hamiltonian modeling according to claim 2, characterized in that: The process of constructing the energy eigenvalue equation includes: Constructing a total Hamiltonian, wherein the total Hamiltonian is the sum of an electromagnetic field kinetic energy term, a potential field correction term, a spin-magnetic field term, and a spin-orbit term; Expanding the total Hamiltonian into the sum of a kinetic energy term and a potential energy term; Simplifying the kinetic energy term in the total Hamiltonian to obtain a first simplified total Hamiltonian; Simplifying the potential energy term in the total Hamiltonian to obtain a second simplified total Hamiltonian; The second simplified total Hamiltonian is converted from three dimensions to one dimension by coordinate system transformation to obtain the final simplified total Hamiltonian; Based on the final simplified total Hamiltonian, the one-dimensional steady-state Schrödinger equation is constructed and solved in different regions, finally obtaining the energy eigenvalue equation and the regional analytical solution form.

4. A method for detecting concrete cracks based on quantum detection and Hamiltonian modeling according to claim 3, characterized in that: The energy eigenvalue equation is expressed as: The analytical solution is in the form of: Where d is the vertical distance between the crack plane and the detection reference plane; L z' is the crack depth in the direction of the crack plane normal vector z'; κ is the attenuation coefficient of the complete area; k z' is the quantum wave number in the crack region; in, Where V0 is the initial potential energy of the crack-free concrete region in the external potential field; ΔV is the volume change caused by the crack region in the external potential field; E is the wave number energy; m is the particle mass; is the reduced Planck constant; ΔV' is the potential energy change caused by the crack area in the new coordinate system.

5. The method for detecting concrete cracks based on quantum detection and Hamiltonian modeling according to claim 3, characterized in that: The coordinate system transformation process includes: The original coordinate system xoy is the detection reference plane, the z axis is the beam emission direction; the angle between the crack plane and the detection reference plane is the crack inclination; the azimuth of the crack plane within the detection reference plane is the crack azimuth; The x′oy′ plane of the new coordinate system is the crack plane, and the z′ axis is the direction of the normal vector of the crack plane; According to the trigonometric function relationship between the fracture dip angle and the fracture azimuth angle, the transformation matrix is ​​obtained; According to the transformation matrix, the three-dimensional kinetic energy term is converted into a one-dimensional kinetic energy term.

6. A method for detecting concrete cracks based on quantum detection and Hamiltonian modeling according to any one of claims 1 to 5, characterized in that: Also includes: The crack width and length are calculated based on the crack projection size, potential energy correction term, and calibration parameters.

7. A method for detecting concrete cracks based on quantum detection and Hamiltonian modeling according to any one of claims 1 to 5, characterized in that: Also includes: Dynamically correct the detection range boundaries.

8. A concrete crack detection device based on quantum detection and Hamiltonian modeling, characterized in that: include: Data acquisition module, used to collect detection data during the quantum beam scanning detection process; A calculation module, for calculating the crack projection size and crack inclination angle based on the collected detection data; And, a quantum geometric inversion module is used to substitute the crack inclination angle into the quantum-geometric joint formula to calculate the depth of the crack to the detection reference plane; wherein, the quantum-geometric joint formula is expressed as: the depth of the crack to the detection reference plane is equal to the product of the crack depth in the direction of the normal vector of the plane where the crack is located and the cosine value of the crack inclination angle; the crack depth is obtained by inversion calculation of a quantum model introducing the quantum mechanical Hamiltonian.

9. The concrete crack detection device based on quantum detection and Hamiltonian modeling according to claim 8, characterized in that: The quantum geometry inversion module includes: A crack depth solving unit is used to solve the crack depth in the normal vector direction of the crack plane; And, an actual depth solving unit is used to substitute the crack depth and crack inclination into the quantum-geometric joint formula to calculate the depth from the crack to the detection reference plane.

10. A concrete crack detection device based on quantum detection and Hamiltonian modeling according to claim 8 or 9, characterized in that: Also includes: The crack parameter calculation module calculates the crack width and length based on the crack projection size, potential energy correction term, and calibration parameters; And / or, a dynamic correction module, used to dynamically correct the detection range boundary.

Images