1. Introduction
Aeolian sand is a widely used material in engineering that can be used as a fine aggregate to prepare concrete to ensure that it meets the requirements of general engineering applications [
1]. Research status of aeolian sand concrete has been studied and applied as engineering material for a long time, and many scholars regarded aeolian sand as fine aggregate. The use of aeolian sand to replace all or part of river sand in the preparation of concrete had a certain scope and application prospect in the area where river sand resources are scarce [
2]. As for the workability of aeolian sand concrete, the research shows that the partial replacement of river sand by aeolian sand can improve the workability of concrete [
3,
4]. Both the collapse and water absorption of concrete increased with the increase in the aeolian sand replacement rate [
5]. For the mechanical properties of aeolian sand concrete, the research shows that the strength of aeolian sand concrete was inversely proportional to amount of aeolian sand [
6]. The mechanical properties of aeolian sand concrete can be improved effectively by using fiber, chemical admixtures and mineral admixtures [
7,
8]. As for the frost resistance of aeolian sand concrete, the research shows that aeolian sand within 30% can effectively inhibit the freeze–thaw damage of light aggregate concrete [
9]. The introduction of aeolian sand can affect the spatial arrangement of hydration products [
10].
Regarding the chosen concrete reliability analysis method, a considerable amount of research has been carried out by domestic and international scholars. Zhao Gaosheng et al. used the finite volume method and the Monte Carlo simulation method to study the reliability of marine concrete with different sections of the material being subjected to chloride ion corrosion. The results of this study demonstrate that the finite volume method yields high accuracy in reliability analysis and highlights the significant impact of section shape on the durability of concrete structures [
11]. Jiang et al. proposed a specific reliability evaluation method based on multiple degradation factors. The copula function was applied to combine two edge distribution functions in order to derive the joint distribution function, and the residual durability reliability model was subsequently established. The validity of the model was later verified using experimental data [
12]. Lydia et al. established a mechanical reliability model of reliability index and failure probability for square composite columns filled with ordinary concrete and high-performance concrete under axial compression by using the response surface method. The findings of their study indicate that the properties and dimensions of the column play a crucial role in determining its strength and reliability. This approach was applied to assess the responsiveness of random parameters in terms of structural reliability [
13]. Qiao Hongxia et al. established a freeze–thaw failure reliability calculation model based on Palmgren’s theory to predict the remaining life of ceramic powder reclaimed concrete. The reliability calculation model of freeze–thaw failure based on the Palmgren model has been proven to be highly reliable through calculation and verification. The model accurately depicts the correlation between the reliability of ceramic powder recycled concrete and freeze–thaw cycles, making it a practical and applicable tool in various scenarios [
14]. The concrete life analysis of probability theory can be divided into two categories: (1) the reliability function of concrete can be established by using the Wiener distribution probability method to reflect the remaining life of specimens with the optimal ratio [
15,
16], and (2), based on the Weibull probability method, a failure model can be established with concrete mass loss and relative dynamic elastic modulus as indexes for freeze–thaw reliability analysis [
17,
18]. According to the latest domestic and international research findings in this field, a significant portion of studies are primarily focused on qualitative research. The utilized analysis methodology predominantly delineates the degradation patterns through experimental data, with minimal utilization of artificial intelligence models for advanced life prediction and reliability analysis. Simultaneously, the traditional approaches exhibit certain limitations in that they rely exclusively on linear regression, making it challenging to ensure precise fitting accuracy for reliability models characterized by strong nonlinearity.
A neural network can accurately fit nonlinear functions; as a result, it possesses a multitude of applications in material constitutive relations and parameter fitting. Gasperlin M. et al. proposed the basic mathematical model of a neural network to predict the viscoelastic behavior of an emulsion system with a certain degree of accuracy [
19]. Xue J. et al. introduced a highly effective approach leveraging artificial neural networks. Taking typical concrete materials as an example, the authors established a macroscopic analytical strength criterion from three steps of the nonlinear homogenization process. This model can effectively predict the friction coefficient and cohesion of porous cement slurry on a microscopic scale with good accuracy [
20]. Al-Haik M. S. et al. studied the stress–relaxation behavior of polymer composites based on artificial neural networks and predicted a broader range of nonlinear models more accurately [
21]. Logsig function and Tansig function are commonly used as activation functions of neurons in traditional neural networks. As general function approximators, they have the following disadvantages: the calculation process of the neural network is separated from the fitting equation, which leads to poor generalization ability and low accuracy of the fitting calculation [
22].
Based on the recent results of domestic and international research, most studies remain in the qualitative research stage, and the analysis method utilized primarily describes the performance degradation law with experimental data [
23]. When probability theory is used for analysis, life prediction and reliability are usually analyzed via linear regression. Because the correlation between damage parameters cannot be accurately described and the failure surface is highly nonlinear, the accuracy and efficiency of reliability analysis are not high. In the present study, aeolian sand with low-quality natural resources was used as concrete fine aggregate. Firstly, we carried out a frost resistance durability test on the aeolian sand concrete, and the microstructural characterization of the concrete was used to determine the durability damage and deterioration mechanism of the concrete through the use of scanning electron microscopy (SEM) and X-ray diffraction (XRD). Secondly, the joint probability density function of damage parameters was constructed using a statistical method and the Copula function to accurately describe the correlation between damage parameters. Finally, due to the highly nonlinear failure surface, the dual neural network method was used to calculate the reliability and failure probability of wind-accumulated sand concrete under different mix ratios, and the reliability function was established to describe the relationship between the reliability of concrete and the freeze–thaw cycle. This method was used to reflect the concrete’s remaining life, providing theoretical guidance for engineering practice and relevant evaluation and identification.
4. Construction of Aeolian Sand Concrete Damage Parameter Joint Probability Density Function
The sample size of concrete damage parameter data obtained in the experimental section of the present study is small, and the distribution law of data samples cannot be known due to many uncertain factors. Therefore, the key to accurately describing the parameter joint probability density function is determining whether reliability can be accurately calculated.
4.1. Damage Parameter Distribution Law
The statistical analysis software SPSS (version 26) was used to draw the distribution histogram of the relative dynamic elastic modulus
P and mass loss rate
W of OC–20% under 200 freeze–thaw cycles, as shown in
Figure 7. It can be seen from the figure that the distribution of
P and
W is basically symmetric.
Figure 8 is the P–P figure of
P and
W, and it can be seen that each observation point is roughly around a straight line.
In summary, the relative dynamic elastic modulus
P and mass loss rate
W distributions can be considered normal distributions. Thus, it is necessary to write the two–parameter edge probability density function
and
as
where
a1 = 60.6598 and
b1 = 7.3395;
a2 = 0.8499 and
b2 = 0.0101.
4.2. Copula Method for Series Structure Failure Mode
For a series structure system, a single component failure will lead to structural system failure, and then, the series system reliability
Pr can be expressed as follows:
The safety domain of series system
is obtained, and the failure domain of two variables series structure failure modes is shown in
Figure 9.
If we have the joint probability density function
fX(
x) of random variable
x = [
x1,
x2, …,
xn]
T, then the theoretical reliability of the series system can be determined.
As shown in Equation (6), accurately determining the joint probability density function of concrete damage parameters is crucial for determining the reliability of concrete structures. However, constructing this joint probability density function has consistently posed a significant challenge. Based on this challenge, a construction method of performance parameters’ joint probability density function based on the Copula function is proposed in the present paper.
According to Sklar’s theorem [
27], the joint distribution function of variable
x1,
x2,
…,
xn is expressed as
where
is the edge distribution function of the variable
;
is
a Copula function; and
is
the related parameters of the Copula function.
Taking the derivative of Equation (5), the joint probability density function of variables
x1,
x2, …,
xn is written as
where
is
the edge probability density function of variables
x1,
x2,
…,
xn and
is
the density function of the Copula function.
Table 2 shows the common Copula function types.
4.3. Construction of Concrete Structures Damage Parameter Joint Probability Density Function
The joint probability density function fX(P,W) construction of the relative dynamic elastic modulus P and mass loss rate W under 200 freeze–thaw cycles of OC–20% is taken as an example to illustrate the method.
We first drew the test scatter plot for
P and
W, with the black dots shown in
Figure 10. As can be seen from the figure, there is a positive correlation between the two parameters; therefore, the function that can express the positive correlation was selected as the Copula function. To construct a two–parameter joint probability density function, relevant parameter values of each Copula function were obtained, as shown in
Table 3.
Figure 10 shows the scatter plots simulated by each Copula function when the number of simulations
n = 200. It can be seen from the figure that the scatter plot generated by each Copula function can suitably cover the original observed data of relative dynamic elastic modulus and mass loss rate, demonstrating the ability of the Copula function to fit the measured data.
The joint probability density function
fX(
P,
W) of relative dynamic elastic modulus
P and mass loss rate
W can be written as:
5. Analysis of Reliability Results Based on the Dual Neural Network
Considering that the performance function of the concrete structure is Z = G(
P,
W), the joint probability distribution function of basic random variable
X = (
P,
W) is
FX(
P,
W), and the joint probability density function is
fX(
P,
W), then the reliable probability
Pr of the concrete structure can be written as follows:
It is widely known that a concrete structural system consists of two series of performance function components. According to the Standard Test Method for Long–Term Performance and Durability of Ordinary Concrete (GB/T 50082-2009), once the relative dynamic elastic modulus reaches 60% or the mass loss rate reaches 5%, the concrete specimen meets the failure criteria. Its performance function expression is shown in Equation (11):
5.1. Dual Neural Network Method
The reliability of a concrete structure can be calculated based on the defined series structural reliability.
where the indicated function
takes
on 1 when min (
)
> 0; otherwise, it is equal to 0.
It can be observed that neural network A, which comprises three layers, can represent the original function of multiple integrals, as stated by a multi–layer neural network capable of approximating any nonlinear function with high precision. The network structure is shown in
Figure 11. The input and output relation in scalar form is:
It can be respectively determined via the partial derivation of variables
If
, Equation (14) can be written as a function of network
output and input variables:
If
,
Equation (15) can be simplified as
The neural network B structure is illustrated in
Figure 12.
Evidentially, the original function of integrand
in Equation (16) is
. Therefore, neural network A and neural network B are considered dual neural networks [
28].
By comparing Equations (13) and (16), it can be seen that the original function network and integrand function network have the same network structure and are both three–layer neural networks with n input, single output, and m hidden layer elements. The above two networks are related in connection weight, threshold, and activation function type. In the integrated function network B, the connection weight from the input layer to the hidden layer element is ,
the connection weight from the hidden layer to the output layer element is . Original function network
A the connection weight from the input layer to the hidden layer unit is ,
the connection weight of the hidden layer to the output layer unit is ,
the threshold value of the output layer unit is c. Based on the above
structure, it can be observed that when network B in a dual neural network
satisfies the given relation and approximates the integrand function in the
integral, network A approximates the original function.
It can be seen from the theory of multivariate integration that the integral value can be expressed as the weighted algebraic sum of the original function
at
the hypercube vertices
,
specifically expressed as follows:
The dual neural network method was used to integrate Equation (12), and the reliability of OC–20% under 200 freeze–thaw cycles was determined.
The variables
P and
W are divided into equal parts in the interval 100, and the integrand function network B is crossed in pairs to form the input sample point, and the network output value of the corresponding sample point is calculated using the integrand function. The training sample is shown in
Table 4.
Using the dual neural network relation, the original function network A is obtained. The vertices of the hypercube in the input sample of network B are simulated through network A, and the sample points are shown in
Table 5. The numerical solution of the integral, namely the structural reliability, can be obtained by substituting the simulation values into Equation (17).
Through the above calculation, the structural reliability and failure probability of OC–40%, OC–60%, LAC–20%, LAC–40% and LAC–60% under 200 freeze–thaw cycles can be similarly obtained, as outlined in
Table 5 and
Table 6. To further illustrate the performance of the proposed method, the Weibull probability distribution model (WPDM) was used for evaluation.
5.2. Weibull Probability Distribution Model (WPDM)
The decay model expression of aeolian sand concrete total failure energy
U under freeze–thaw cycles was established using the Weibull probability distribution as follows:
where
U0 represents the total failure energy of concrete in the absence of freeze–thaw action;
Un denotes the total failure energy of concrete after
n freeze–thaw cycles; and
Dn is the damage variable after
n times freeze–thaw.
Assume that
fn represents the probability density function following
n freeze–thaw cycles.
where
a is the shape parameter and
b is the size parameter.
The failure probability distribution function is derived by integrating Equation (19).
The relationship between failure probability and the number of freeze–thaw cycles is shown in Equation (20). The freeze–thaw failure of aeolian sand concrete is the result of the long–term accumulation of freeze–thaw damage; thus, the damage variable and failure probability are simultaneously superimposed. After
n times freeze–thaw cycles of aeolian sand concrete,
Dn represents the damage degree and
Fn represents failure probability. When the concrete reaches freeze–thaw failure,
Dn =
Fn = 1, the failure probability is equivalent to the damage degree:
5.3. Analysis of Reliability Results
In addition, Monte Carlo simulation (MCS) with 1 million sample points was used to evaluate the reliability of the series structure system, and the results are listed in
Table 6 and
Table 7.
Similarly, according to the above method, the reliability and failure probability can be determined under 25, 50, …, and 175 freeze–thaw cycles; a custom neural network was constructed with reliability and failure probability as inputs and the number of freeze–thaw cycles as outputs. According to the empirical equation [
29], the number of hidden layer elements was taken as the hidden layer
and
activation function
.
After training 1000 steps, the relative dynamic modulus and mass loss rate
under the same freeze–thaw cycle were taken as the
x– and
y–axis
coordinates, and the fitting curve was drawn, as shown in
Figure 13. The maximum relative error of the
fitting results compared with the training output sample is listed in
Table 8.
If either the relative dynamic elastic modulus loss or mass loss meets the failure standard, the material is considered to have failed. Following training of the custom neural network, the predicted input samples that have reached failure are taken as the input of the neural network, and the number of freeze–thaw cycles when failure is obtained, as shown in
Table 9.
From the reliability calculation results of different methods shown in
Table 6 and
Table 7, it can be seen that with MCS as the exact solution, the calculation accuracy of the proposed method depicted in this paper is better than that of the WPDM, and the maximum relative error is less than 2%.
It can be seen from
Table 8 that the fitting accuracy of the custom neural network model for the freeze–thaw cycle failure of aeolian sand concrete is higher. Among the various types of aeolian sand concrete, the highest relative error is 4.41 × 10
−4, with the remaining errors all falling below this value, with some found to be even as low as 1.0 × 10
−6 orders of magnitude, demonstrating excellent calculation accuracy.
Based on the proposed method described in this paper, the reliability analysis model was obtained to accurately calculate the reliability of various types of aeolian sand concrete to accurately predict the number of freeze–thaw cycles of various types of aeolian sand concrete during failure.
As can be seen from
Table 9, OC–60% and LAC–60% fail first, and the number of freeze–thaw cycles for these samples is 162 and 127, respectively, which is consistent with the test results. This result verifies the effectiveness of the custom neural network model. The freeze–thaw cycle durability of the aeolian sand concrete with equal replacement rates of 20% and 40% is evidently better than that of the aeolian sand concrete with 60% aeolian sand concrete, and the optimal ratio is the aeolian sand concrete with a replacement rate of 40%. The reliability of light aggregate concrete with aeolian sand replacement rates of 20% and 40% is slightly better than that of ordinary concrete with aeolian sand, with the best concrete being LAC–40%, and the number of freeze–thaw cycles reaching 378.