EVALUATION METHOD FOR CORROSION DAMAGE EVOLUTION OF UNDERWATER CONCRETE STRUCTURES

Abstract

An evaluation method for corrosion damage evolution of underwater concrete structures includes performing the time reversal test on the concrete beam specimen placed in the water, performing the uniaxial compression test on the concrete cube specimens; immersing the concrete beam specimen and the concrete cube specimens in a hydrochloric acid solution, and performing the time reversal test on the concrete beam specimen on the 10th, 20th and 30th days respectively. At the same time, a concrete cube specimen is taken out to perform the uniaxial compression test on the 10th, 20th and 30th days respectively; and using the above calculation results to evaluate the corrosion evolution process thereof without damaging the underwater concrete structure.

Claims

1. An evaluation method for corrosion damage evolution of underwater concrete structures, comprising: (a) immersing a concrete beam specimen and a plurality of concrete cube specimens in water to obtain a first immersed concrete beam specimen and a plurality of first immersed concrete cube specimens; (b) performing a time reversal test on the first immersed concrete beam specimen in the water to obtain a damage index DI.sup.0 of a concrete in a water-immersed state, performing a uniaxial compression test on the plurality of first immersed concrete cube specimens to obtain a compressive strength F.sub.cp.sup.0 and an elastic modulus E.sub.c.sup.0 of the concrete in the water-immersed state; (c) immersing the first immersed concrete beam specimen and the plurality of first immersed concrete cube specimens in a hydrochloric acid solution to mimic a long-term corrosion damage of the concrete in the water-immersed state to obtain a second immersed concrete beam specimen and a plurality of second immersed concrete cube specimens, performing the time reversal test on the second immersed concrete beam specimen in the hydrochloric acid solution on 10th, 20th, and 30th days respectively to obtain a damage index DI′ of a concrete in a corrosion state, and at the same time, performing the uniaxial compression test on the plurality of second immersed concrete cube specimens on 10th, 20th, and 30th days respectively to obtain a compressive strength F′.sub.cp and an elastic modulus E′.sub.c of the concrete in the corrosion state; (d) calculating a corrosion index CI with formula (1), calculating a concrete compressive strength loss rate LR.sub.F with formula (2), calculating a concrete elastic modulus loss rate LR.sub.E with formula (3), and $\begin{array}{cc}\mathrm{CI}={\mathrm{DI}}^{\prime }-{\mathrm{DI}}^0& \left(1\right)\end{array}$ $\begin{array}{cc}{\mathrm{LR}}_F=\frac{F_{cp}^0-F_{cp}^{\prime }}{F_{\mathrm{cp}}^0}& \left(2\right)\end{array}$ $\begin{array}{cc}{\mathrm{LR}}_E=\frac{E_c^0-E_c^{\prime }}{E_c^0}& \left(3\right)\end{array}$ in the formulas, CI is the corrosion index, DI.sup.0 is the damage index of the concrete in the water-immersed state, DI′ is the damage index of the concrete in the corrosion state, LR.sub.F is the concrete compressive strength loss rate, F.sub.cp.sup.0 is the compressive strength of the concrete in the water-immersed state, F′.sub.cp is the compressive strength of the concrete in the corrosion state, LR.sub.E is the concrete elastic modulus loss rate, E.sub.c.sup.0 is the elastic modulus of the concrete in the water-immersed state, and E′.sub.c is the elastic modulus of the concrete in the corrosion state; and (e) calculating an absolute error δ.sub.CF between the corrosion index and the concrete compressive strength loss rate with formula (4), and calculating the absolute error δ.sub.CE between the corrosion index and the concrete elastic modulus loss rate with formula (5) to realize an evaluation of the corrosion damage evolution of the underwater concrete structures;
δ.sub.CE=CI−LR.sub.F  (4)
δ.sub.CE=CI−LR.sub.E  (5) in the formulas, δ.sub.CF is the absolute error between the corrosion index and the concrete compressive strength loss rate, CI is the corrosion index, LR.sub.F is the concrete compressive strength loss rate, δ.sub.CE is the absolute error between the corrosion index and the concrete elastic modulus loss rate, and LR.sub.E is the concrete elastic modulus loss rate.

2. The evaluation method for the corrosion damage evolution of the underwater concrete structures according to claim 1, wherein a method for obtaining a damage index of a concrete by performing the time reversal test on the concrete beam specimen comprises: (21) determining an excitation signal V.sub.A(t); (22) inputting the excitation signal V.sub.A(t) to a sensor at position A on the concrete beam specimen, an expression of response signals received by a sensor at position B in a frequency domain is formula (6),
V.sub.B(r,ω)=k.sub.A(ω)k.sub.B(ω)G(r,ω)V.sub.A(W)  (6) an expression of the response signals in a time domain is formula (7), $\begin{array}{cc}{V}_B\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{V}_B\left(r,\omega \right)e^{i\omega t}d\omega & \left(7\right)\end{array}$ in the formula, r is a distance from the position A to the position B, k.sub.A is an electromechanical coupling coefficient of the sensor at the position A, k.sub.B is an electromechanical coupling coefficient of the sensor at the position B, and G(r, ω) is a transfer function from the sensor at the position A to the sensor at the position B; V.sub.B(t) is the response signals received by the sensor at the position B in the time domain; (23) applying a rectangular window, as shown in formula (8), to a response signal to remove a crosstalk, and then performing a time reversal to obtain the reversed signal, $\begin{array}{cc}w\left(t\right)=\{\begin{array}{cc}1,& 0\le t\le T\\ 0,& \mathrm{else}\end{array}& \left(8\right)\end{array}$ an expression of a reversed signal in the frequency domain is formula (9),
{circumflex over (V)}.sub.B(r,ω)=k*.sub.A(ω)k*.sub.B(ω)G*(r,ω)V*.sub.A(ω)e.sup.iωt  (9) an expression of a reversed signal in the time domain is formula (10), $\begin{array}{cc}{V}_B\left(-t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{V}_B^{*}\left(r,\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega & \left(10\right)\end{array}$ in the formulas, V*.sub.B is a phase conjugation of V.sub.B, * is a complex conjugation operator, T is a sampling duration, and r is the distance from the A to the B; (24) inputting the reversed signal to the sensor at the position B, then the sensor at the position A receives a focused signal, an expression of the focused signal in the frequency domain is formula (11), and
{circumflex over (V)}.sub.A(r,ω)=k.sub.A(ω)k*.sub.A(ω)k.sub.B(ω)k*.sub.B(ω)G(r,ω)G*(r,ω)V*.sub.A(ω)e.sup.iωT  (11) an expression of the focused signal in the time domain is formula (12); $\begin{array}{cc}\begin{array}{c}{\hat{V}}_A\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\hat{V}}_A\left(r,\omega \right)e^{i\omega T}d\omega \\ =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{k}_A\left(\omega \right)k_A^{*}\left(\omega \right)k_B\left(\omega \right)k_B^{*}\end{array}& \left(12\right)\end{array}$ $\left(\omega \right)G\left(r,\omega \right)G^{*}\left(r,\omega \right)V_A^{*}\left(\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega$ (25) performing a time reversal reconstruction on a focused signal to obtain a reconstructed signal, and an expression of the reconstructed signal in the frequency domain is formula (13),
{tilde over (V)}.sub.A(r,ω)={circumflex over (V)}*.sub.A(r,ω)e.sup.iωT  (13) in the formula, {circumflex over (V)}*.sub.A is a phase conjugation of {circumflex over (V)}.sub.A, and * is the complex conjugation operator; and an expression of a reconstructed signal in the time domain is formula (14); $\begin{array}{cc}\begin{array}{c}{\tilde{V}}_A\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\hat{V}}_A^{*}\left(r,\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega \\ =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{k}_A\left(\omega \right)k_B\left(\omega \right)k_A^{*}\left(\omega \right)k_B^{*}\end{array}& \left(14\right)\end{array}$ $\left(\omega \right)G\left(r,\omega \right)G^{*}\left(r,\omega \right)V_A\left(\omega \right)e^{i\omega T}d\omega$ (26) performing a normalization on the excitation signal V.sub.A(t) and the reconstructed signal with formula (15),
N.sub.A(t)=V.sub.A(t)/max(V.sub.A(t))
Ñ.sub.A(t)={tilde over (V)}.sub.A(t)/max({tilde over (V)}.sub.A(t))  (15) in the formulas, V.sub.A(t) is the excitation signal, {tilde over (V)}.sub.A is the reconstructed signal, N.sub.A is a normalized excitation signal, and Ñ.sub.A is a normalized reconstructed signal; and (27) substituting the normalized excitation signal and the normalized reconstructed signal into formula (16) to calculate a damage index; $\begin{array}{cc}\mathrm{DI}=1-\sqrt{{\left[{\int }_{t_0}^{t_1}{N}_A\left(t\right){\tilde{N}}_A\left(t\right)dt\right]}^2/{\int }_{t_0}^{t_1}{\left[N_A\left(t\right)\right]}^{2}dt{\int }_{t_0}^{t_1}{\left[{\tilde{N}}_A\left(t\right)\right]}^2\mathrm{dt}}\times 100\%& \left(16\right)\end{array}$ in the formula, DI is the damage index, t.sub.0 and t.sub.1 are a start time and an end time of a signal comparison interval separately, N.sub.A(t) is the normalized excitation signal, and N.sub.A(t) is the normalized reconstructed signal.

3. The evaluation method for the corrosion damage evolution of the underwater concrete structures according to claim 2, wherein a method for determining the excitation signal V.sub.A(t) comprises: an expression of a 5-cycle sine function modulated by a Hanning window is formula (17), $\begin{array}{cc}x\left(t\right)=\sin \left(2\pi f_{c}t\right)\left[1-\cos \frac{2\pi f_c^t}{5}\right]& \left(17\right)\end{array}$ in the formula, f.sub.c is a central frequency of modulated signals; a method for determining the central frequency of the modulated signals is as follows: modulating modulated signals with different central frequencies from 0-10 MHz at an interval of 10 kHz, inputting the modulated signals with different central frequencies to the sensor at the position A on the concrete beam specimen separately, response signals are received by the sensor at the position B, comparing the response signals and selecting a modulated signal with a largest amplitude of the response signals as the excitation signal V.sub.A(t).

4. The evaluation method for the corrosion damage evolution of the underwater concrete structures according to claim 1, wherein a compressive strength F.sub.cp of the concrete i is obtained with formula (18), and an elastic modulus E.sub.c of the concrete is obtained with formula (19); $\begin{array}{cc}{F}_{\mathrm{cp}}=\frac{F_{\max }}{A}& \left(18\right)\end{array}$ $\begin{array}{cc}{E}_c=\frac{\left(F_a-F_0\right)}{A}\times \frac{L}{\Delta }& \left(19\right)\end{array}$ in the formulas, F.sub.cp is the compressive strength of the concrete cube specimen, F.sub.max is a failure load, A is a loading area of the concrete cube specimen, E.sub.c is the elastic modulus of the concrete cube specimen, F.sub.a is a load when a stress is ⅔F.sub.cp, F.sub.0 is a load when the stress is ⅓ F.sub.cp, L is a measuring scale distance of the concrete cube specimen, and Δ is a deformation of the concrete cube specimen loaded from F.sub.0 to F.sub.a.

5. The evaluation method for the corrosion damage evolution of the underwater concrete structures according to claim 1, wherein a method for preparing the concrete beam specimen comprises: performing an anti-corrosion, an insulation, and a waterproof treatment on a pair of sensors; building a formwork of the concrete beam specimen; placing the pairs of sensors at predetermined positions in the formwork of the concrete beam specimen; and pouring the concrete, and completing a maintenance according to standards.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to describe the examples and the technical solutions of the present invention, accompanying drawings required by the examples are briefly introduced below. Obviously, the accompanying drawings in the following description are only partial examples of the present invention, and a person of ordinary skill in the art can be able to derive other accompanying drawings from these accompanying drawings without creative efforts. According to the examples of the present invention, corrosion evolution of underwater concrete structures is evaluated by taking corrosion damage caused by immersing C30 concrete specimens in a hydrochloric acid solution with pH=1 as an example.

[0048] FIG. 1 is a flow chart of a method in the present invention;

[0049] FIG. 2 is a schematic diagram of a time reversal test system for a concrete beam specimen in the present invention;

[0050] FIG. 3 is a schematic diagram of the concrete beam specimen in the present invention;

[0051] FIG. 4 is a schematic diagram of sensors in the present invention;

[0052] FIG. 5 is a flow chart of the time reversal test in the present invention;

[0053] FIG. 6 is a diagram of modulated signals with different central frequencies (60-140 kHz) and corresponding response signals in the present invention;

[0054] FIG. 7 is a diagram of amplitudes of response signals with different central frequencies (60-140 kHz) in the present invention;

[0055] FIGS. 8A-8B are the excitation signals with a central frequency of 100 kHz selected in the present invention: (a) time domain diagram; (b) frequency domain diagram;

[0056] FIGS. 9A-9D are diagrams of the response signal obtained by performing the time reversal test on the concrete beam specimen at different corrosion durations in the present invention;

[0057] FIGS. 10A-10D are diagrams of the reversed signal obtained by performing the time reversal test on the concrete beam specimen at different corrosion durations in the present invention;

[0058] FIGS. 11A-11D are diagrams of the focused signal obtained by performing the time reversal test on the concrete beam specimen at different corrosion durations in the present invention;

[0059] FIGS. 12A-12D are diagrams of the excitation signal and the reconstructed signal obtained by performing the time reversal test on the concrete beam specimen at different corrosion durations and after normalization in the present invention; and

[0060] FIG. 13 is stress-strain curves of concrete cube specimens with different corrosion durations obtained from the uniaxial compression test.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0061] In order to make the technical solution of the present invention better understood and implemented by those skilled in the art, the present invention will be described in detail below with reference to FIG. 1-FIG. 13. The following examples are only illustrative of the technical solution of the present invention more clearly and are not intended to limit the scope of protection of the present invention.

Example 1

[0062] According to the evaluation method for corrosion damage evolution of underwater concrete structures, the basic steps are shown in the flow chart of FIG. 1. The time reversal test system of the concrete beam specimen in the present example is shown in FIG. 2, the time reversal test system is composed of a computer, a signal generator, a signal amplifier and an oscilloscope, before the time reversal test, the computer is connected to the signal generator, the signal generator is connected to the oscilloscope, the oscilloscope is connected to the signal amplifier, clamps are led out from the signal amplifier and the signal generator respectively, and are connected to the concrete beam specimen in the experiment to perform the time reversal test. An anti-corrosion tank is filled with water/hydrochloric acid, which is capable to immerse the concrete specimen.

[0063] The present invention will be further described with reference to FIG. 1 and a test process: [0064] (1) the method for preparing the concrete specimens includes:

[0065] performing anti-corrosion, insulation and waterproof treatment on a pair of sensors, with a specific operation method as shown in FIG. 4: selecting d.sub.33-type PZT patches {circle around (1)} having a diameter of 14 mm and a thickness of 1 mm for the sensors, and encapsulating the d.sub.33-type PZT patches with an acrylic pipe {circle around (4)} having an external diameter of 25 mm, a wall thickness of 2 mm and a height of 20 mm and epoxy resin adhesive {circle around (2)}, and leading out leads {circle around (3)} to achieve the objectives of anti-corrosion, insulation and waterproof;

[0066] after that, building a formwork of the concrete beam specimen, placing the sensors at predetermined positions A and B in the formwork, after leading out the leads {circle around (3)}, pouring concrete, and making C30 concrete specimen by selecting ordinary Portland cement with a grade of 32.5, fine aggregates are sands with particle sizes of 0.25-0.5 mm, and coarse aggregates are stones with particle sizes of 5-30 mm according to a mass ratio of 1:0.958:2.462 and a water-cement ratio of 0.38. As shown in FIG. 3, the size of the concrete beam specimen is 500 mm×100 mm×100 mm (length×width×height), the size of concrete cube specimens is 100 mm×100 mm×100 mm, the two sensors are located on the longitudinal axis line 20 mm away from the left and right ends respectively, and maintenance is completed according to standards. [0067] (2) after being maintained, immerse the concrete beam specimen and the cube specimens in water for 10 days to achieve a state of full immersion.

[0068] next, the time reversal test of stress waves is performed on the concrete beam specimen, before the test, as shown in FIG. 3, two clamps are connected to the two leads {circle around (3)}, the two leads {circle around (3)} are led out from the sensor {circle around (5)} at A and the sensor {circle around (6)} at B on the concrete beam specimen separately, and the sensor {circle around (5)} at A and the sensor {circle around (6)} at B are preset on the concrete beam specimen.

[0069] a process of the time reversal test is shown in FIG. 5, the excitation signal V.sub.A(t) is selected at first, and the method for determining the excitation signal V.sub.A(t) includes: [0070] (3) generating a modulated signal with an amplitude of 8 Vp-p, a center frequency of 0-10 mHz and an interval of 10 kHz of the 5-cycle sine function modulated by a Hanning window by the signal generator, with expression as follows:

[00010]$x\left(t\right)=\sin \left(2\pi f_{c}t\right)\left[1-\cos \frac{2\pi f_{c}t}{5}\right]$

[0071] in the formula, f.sub.c is the center frequency of the modulated signal, the modulated signal is input to the sensor {circle around (5)} at A on the concrete beam specimen separately, the response signal is received by the sensor {circle around (6)} at B, a waveform diagram of the response signal in a frequency band of 60-140 kHz is shown in FIG. 6, and the maximum amplitude of the response signal is shown in FIG. 7. Selecting the corresponding signal with the maximum amplitude of the response signal, that is, the modulated signal of the 5-cycle sine function modulated by a Hanning window with the central frequency of 100 kHz as the excitation signal V.sub.A(t) for the follow-up test, with time domain and frequency domain as shown in FIGS. 8A-8B. [0072] (4) inputting the excitation signal V.sub.A(t) modulated by the 5-cycle sine function modulated by a Hanning window with an amplitude of 8 Vp-p and a central frequency of 100 kHz is to the sensor {circle around (5)} at A on the concrete beam specimen through the signal generator, propagating stress waves from A to B, and receiving the response signal V.sub.B(t) by the sensor {circle around (6)} at B, with the expression of the response signal in frequency domain as follows:

V.sub.B(r,ω)=k.sub.A(ω)k.sub.B(ω)G(r,ω)V.sub.A(ω)

[0073] the expression of the response signal in time domain follows:

[00011]$V_B\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{V}_B\left(r,\omega \right)e^{i\omega t}d\omega$

[0074] In the formula, r is the distance from A to B (that is, the distance of propagation of stress waves), k.sub.A is an electromechanical coupling coefficient of the sensor {circle around (5)} at A, k.sub.B is an electromechanical coupling coefficient of the sensor {circle around (6)} at B, and G(r, ω) is the transfer function from the sensor {circle around (5)} at A to the sensor {circle around (6)} at B. The response signal obtained in the test is shown in FIGS. 9A-9D, and in this step, the sensor {circle around (5)} at A and the sensor {circle around (6)} at B are used as an exciter and a receiver respectively; [0075] (5) applying a rectangular window to remove the crosstalk from the response signal V.sub.B(t), and then performing the time reversal to obtain the reversed signal V.sub.B(−t), as shown in FIGS. 10A-10D, the expression of the rectangular window function is as follows:

[00012]$\begin{array}{c}{\hat{V}}_A\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\hat{V}}_A\left(r,\omega \right)e^{i\omega t}d\omega \\ =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{k}_A\left(\omega \right)k_A^{*}\left(\omega \right)k_B\left(\omega \right)k_B^{*}\left(\omega \right)G\left(r,\omega \right)G^{*}\left(r,\omega \right)\\ V_A^{*}\left(\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega \end{array}$

[0076] the expression of the reversed signal in frequency domain is:

{circumflex over (V)}.sub.B(r,ω)=k*.sub.A(ω)k*.sub.B(ω)G*(r,ω)V*.sub.A(ω)e.sup.iωt

[0077] the expression of the reversed signal in time domain is:

[00013]$V_B\left(-t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{V}_B^{*}\left(r,\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega$

[0078] In the formula, V*.sub.B is the phase conjugation of V.sub.B, * is the complex conjugation operator, T is the sampling duration, and r is the distance from A to B; [0079] (6) inputting the reversed signal in step (5) to the sensor {circle around (6)} at B on the concrete beam specimen, and receiving the focused signal by the sensor {circle around (5)} at A, as shown in FIGS. 11A-11D, the expression of the focused signal in frequency domain is as follows:

{circumflex over (V)}.sub.A(r,ω)=k.sub.A(ω)k*.sub.A(ω)k.sub.B(ω)k*.sub.B(ω)G(r,ω)G*(r,ω)V*.sub.A(ω)e.sup.iωT

[0080] the expression of the focused signal in time domain is:

[00014]$\begin{array}{c}{\hat{V}}_A\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\hat{V}}_A\left(r,\omega \right)e^{i\omega t}d\omega \\ =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{k}_A\left(\omega \right)k_A^{*}\left(\omega \right)k_B\left(\omega \right)k_B^{*}\end{array}$$\left(\omega \right)G\left(r,\omega \right)G^{*}\left(r,\omega \right)V_A^{*}\left(\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega$ [0081] (7) performing time reversal reconstruction on the focused signal in step (6) to obtain the reconstructed signal, with the expression of the reconstructed signal in frequency domain as follows:

{tilde over (V)}.sub.A(r,ω)={circumflex over (V)}*.sub.A(r,ω)e.sup.iωT

[0082] In the formula, {circumflex over (V)}*.sub.A is the phase conjugation of {circumflex over (V)}.sub.A, and * is the complex conjugation operator.

[0083] the expression of the reconstructed signal in time domain is:

[00015]$\begin{array}{c}{\tilde{V}}_A\left(t\right)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\widehat{V_A^{*}}\left(r,\omega \right)e^{i\omega T}{e}^{i\omega t}d\omega \\ =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{k}_A\left(\omega \right)k_B\left(\omega \right)k_A^{*}\left(\omega \right)k_B^{*}\end{array}$$\left(\omega \right)G\left(r,\omega \right)G^{*}\left(r,\omega \right)V_A\left(\omega \right)e^{i\omega t}d\omega$ [0084] (8) performing normalization on the excitation signal in step (3) and the reconstructed signal in step (7), with the formula of normalization as follows:

N.sub.A(t)=V.sub.A(t)/max(V.sub.A(t))

Ñ.sub.A(t)={tilde over (V)}.sub.A(t)/max({tilde over (V)}.sub.A(t))

[0085] In the formula, V.sub.A(t) is the excitation signal, {tilde over (V)}.sub.A(t) is the reconstructed signal, N.sub.A(t) is the normalized excitation signal, and Ñ.sub.A(t) is the normalized reconstructed signal; and

[0086] the normalized excitation signal and the normalized reconstructed signal are shown in FIGS. 12A-12D. The damage index can be calculated by substituting the normalized excitation signal and the normalized reconstruction signal into the following formula:

[00016]$\mathrm{DI}=1-\sqrt{{\left[{\int }_{t_0}^{t_1}{N}_A\left(t\right){\tilde{N}}_A\left(t\right)dt\right]}^2/{\int }_{t_0}^{t_1}{\left[N_A\left(t\right)\right]}^{2}dt{\int }_{t_0}^{t_1}{\left[{\tilde{N}}_A\left(t\right)\right]}^2\mathrm{dt}}\times 100\%$

[0087] In the formula, DI is the damage index, t.sub.0 and t.sub.1 are start time and end time of the signal comparison interval separately, N.sub.A(t) is the normalized excitation signal, and Ñ.sub.A(t) is the normalized reconstructed signal.

[0088] Next, the uniaxial compression test is performed on the concrete cube specimens; [0089] (9) performing the uniaxial compression test on concrete cube specimens to obtain the stress-strain curves to calculate the compressive strength and the elastic modulus of concrete with the following formulas:

[00017]$F_{\mathrm{cp}}=\frac{F_{\max }}{A}$$E_c=\frac{\left(F_a-F_0\right)}{A}\times \frac{L}{\Delta }$

[0090] in the formulas, F.sub.cp is the compressive strength of the concrete cube specimen, F.sub.max is the failure load, A is the loading area of the specimen, E.sub.c is the elastic modulus of the concrete cube specimen, F.sub.a is the load when stress is ⅔F.sub.cp, F.sub.0 is the load when stress is ⅓ F.sub.cp, L is the measuring scale distance of the concrete cube specimen, and Δ is the deformation of the concrete cube specimen loaded from F.sub.0 to F.sub.a; [0091] (10) replacing the water in the tank with the hydrochloric acid solution with pH=1, immersing the concrete beam specimen and the cube specimens in the hydrochloric acid solution to introduce corrosion damage, and repeating steps (4) to (9) on the 10th, 20th and 30th days after immersing the concrete beam specimen and the cube specimens in the hydrochloric acid solution, separately, and calculating the damage index, the compressive strength and the elastic modulus of the concrete after corrosion; [0092] (11) using the damage index, the compressive strength and the elastic modulus to calculate the corrosion index, the loss rate of concrete compressive strength and the loss rate of concrete elastic modulus, with the formulas as follows:

[00018]$\mathrm{CI}={\mathrm{DI}}^{\prime }-{\mathrm{DI}}^0$${\mathrm{LR}}_F=\frac{F_{cp}^0-F_{cp}^{\prime }}{F_{\mathrm{cp}}^0}$${\mathrm{LR}}_E=\frac{E_c^0-E_c^{\prime }}{E_c^0}$

[0093] in the formulas, CI is the corrosion index, DI.sup.0 is the damage index of concrete in the water-immersed state, DI′ is the damage index of concrete in the corrosion state, LR.sub.F is the loss rate of concrete compressive strength, F.sub.cp.sup.0 is the compressive strength of concrete in the water-immersed state, F′.sub.cp is the compressive strength of concrete in the corrosion state, LR.sub.E is the loss rate of concrete elastic modulus, E.sub.c.sup.0 is the elastic modulus of concrete in the water-immersed state, and E′.sub.c is the elastic modulus of concrete in the corrosion state. [0094] (12) calculating the absolute error between the corrosion index and the loss rate of concrete compressive strength, and calculating the absolute error between the corrosion index and the loss rate of concrete elastic modulus, with the formulas as follows:

δ.sub.CE=CI−LR.sub.F

δ.sub.CE=CI−LR.sub.E

[0095] in the formulas, δ.sub.CF is the absolute error between the corrosion index and the loss rate of concrete compressive strength, CI is the corrosion index, LR.sub.F is the loss rate of concrete compressive strength, δ.sub.CE is the absolute error between the corrosion index and the loss rate of concrete elastic modulus, and LR.sub.E is the loss rate of concrete elastic modulus.

[0096] FIGS. 12A-12D show the normalized excitation signal and the normalized reconstructed signal, FIG. 13 shows the stress-strain curves of concrete cube specimens with different corrosion durations obtained from the uniaxial compression test, and the above test data are summarized in Table 1 and Table 2 as follows:

TABLE-US-00001 TABLE 1 Corrosion duration t 0 10 20 30 The loss rate of concrete 0 4.46% 9.39% 23.29% compressive strength LR.sub.F Corrosion index CI 0 7.35% 13.39% 23.06% The absolute error δ.sub.CF 0 2.89% 4.00% 0.23%

TABLE-US-00002 TABLE 2 Corrosion duration t 0 10 20 30 The loss rate of concrete elastic 0 12.11% 18.03% 24.08% modulus LR.sub.F Corrosion index CI 0 7.35% 13.39% 23.06% The absolute error δ.sub.CE 0 4.76% 4.64% 1.02%

[0097] It can be seen from Table 1 and Table 2 that the absolute error δ.sub.CF between the corrosion index and the loss rate of concrete compressive strength and the absolute error δ.sub.CE between the corrosion index and the loss rate of concrete elastic modulus are less than 5%, and the corrosion index can reasonably evaluate the corrosion state of concrete. The above results show that the provided method is feasible and effective for evaluating the corrosion damage evolution of underwater concrete structures.

[0098] The evaluation method for corrosion damage evolution of underwater concrete structures provided by the present invention timely and effectively represents degradation degree of mechanical properties of the underwater concrete structures caused by corrosion, provides a feasible method for realizing the corrosion damage evolution evaluation of the underwater concrete structure, and the method is practical and worth popularizing.

[0099] The examples disclosed above are only preferable specific examples of the present invention, but the examples of the present invention are not limited to the above examples, and variations readily conceivable to anyone skilled in the art all fall within the scope of protection of the present invention.