Dynamic Inversion Analysis of Structural Layer Modulus of Semirigid Base Pavement considering the Influence of Temperature and Humidity

This paper is aimed to solve the overlearning problem of the neural network algorithm used to calculate the asphalt concrete pavement structural modulus in reverse. The ﬁreﬂy algorithm was adapted to optimize the selection of support vector machine (SVM) parameters. Based on the optimized SVM model, a new method for dynamic inversion of the semirigid base asphalt concrete pavement structural layer modulus was presented. The results show that the absolute value of relative error of each layer modulus is not more than 3.73% by using the proposed method. Then, the inﬂuences of temperature and humidity on the inversion modulus of semirigid base asphalt concrete pavement in the seasonal frozen area were analyzed, and the correction formula of the inversion modulus was established. The paper is of practical signiﬁcance for improving the safety performance of semirigid base pavement in the seasonal frozen area in China.


Introduction
In recent years, China's highway construction has developed rapidly in terms of construction scale and quantity. e highway has reached a total mileage of 5 million kilometers by 2019, including 143,000 km of expressways. e semirigid base asphalt concrete pavement is most widely used in expressways [1]. In seasonal frozen areas such as Northeast China, Inner Mongolia, Qinghai, Gansu, Ningxia, and Tibet, the asphalt concrete pavement structure bears the frostthawing damage under low temperature and damp environment, and its bearing capacity and stability of material will be greatly affected. ese factors result in pavement cracks, cavities, subsidence, rutting, and other diseases, which seriously affect the traffic operation and normal use of expressways and even threaten the development of society and economy [2]. At present, in the design methods of asphalt pavement in the world, static load or approximate equivalent static load is generally used, but the static load mode cannot reflect the actual stress state of the pavement, so the influence of dynamic load on the design of the pavement structure has become the focus of research. Golewski studies the design of the concrete structure under dynamic load and puts forward the reference criteria for the design, construction, and maintenance of the dynamic loading structure, which provides a reference method for the design of pavement concrete structure [3,4]. e bearing capacity of semirigid base asphalt concrete pavement is important to evaluate pavement performance. Falling weight deflectometer (FWD) is widely used as a kind of nondestructive testing equipment for the asphalt concrete pavement structure. Its principle is to carry out inversion analysis on the modulus of the asphalt concrete pavement structural layer according to the deflection value measured by FWD under impact load [4]. Scrivener et al. first proposed the inverse modulus calculation method of FWD deflection basin by solving the Burmister double layer system and compiled the inversion nomograph [5]. However, the regression analysis method cannot be used to process largescale FWD-measured data with poor inversion accuracy and versatility. Bush developed BISDEF and CHEVDEF software for inversion of the pavement structural layer modulus by using the iterative method [6]. Based on the system identification, Wang used this method to develop a software SIDMOD for inversion of the pavement structural layer modulus [7]. Nevertheless, the inversion results obtained by the iterative method are greatly affected by initial value and iterative method, and the iterative method converges slowly. Based on the database search method, Utahj et al. developed the MODULUS software for inversion of the pavement structural layer modulus [8]. Nonetheless, the database search method takes a long time to calculate and needs to select a more strict range of modulus. Meier et al. used the BP neural network to inverse the modulus and developed ANOVA software for inversion of pavement structural layer modulus [9]. ere are some problems in using the neural network method such as oscillation and divergence out of improper initial value. Fwa et al. studied the inverse calculation of modulus by adopting the genetic algorithm and developed the inverse software NUS-GABACK [10]. Due to the slow calculation speed, of the genetic algorithm, NUS-GABACK software cannot realize the large-scale inverse calculation needed in practical engineering.
e above studies are all static inversion methods. Wang developed a dynamic inversion program integrated with the artificial neural network and genetic algorithm (ANN-GA) to calculate the modulus of the structural layer. According to the asphalt concrete pavement deflection measured by FWD, the service condition of asphalt pavement can be evaluated [11]. You et al. developed a dynamic inversion model based on ANN to predict the interlayer contact conditions and structural layer modulus of the pavement structure in view of the asphalt concrete pavement interlayer contact conditions [12]. Considering the influence of FWD dynamic load characteristics, viscoelastic material parameters, and nonlinear material parameters, Li et al. developed a dynamic inversion program based on ANN in combination with the genetic algorithm [13]. However, the neural network may cause the problem of overlearning.
Asphalt concrete pavement material is a kind of viscoelastic material whose strength and stability are easily affected by temperature changes, so the bearing capacity of the asphalt concrete pavement structure will change with temperature. Subgrade is a linear structure built in the natural environment, and its structural performance will change with the freeze-thaw cycle in the seasonal frozen area. e specific performance is that the change of seasonal climate causes the moisture in the subgrade to freeze-thaw.
is affects the structural performance of the subgrade and ultimately leads to the destruction of the asphalt concrete pavement structure. erefore, the influence of temperature and humidity should be considered when calculating the modulus of the pavement structure layer. Scholars have carried out a lot of research on the relationship between temperature and humidity in seasonally frozen areas and the modulus of pavement material. Chamberlain et al. found that the pore ratio of fine-grained soil changes with freezethaw cycle and modulus, and the permeability increases with the continuous freeze-thaw cycle [14]. Observing the change law of soil material properties after the first freeze-thaw change, Simonsen et al. found that the modulus value of subgrade decreased after the first freeze-thaw [15]. Taking remolded soil as the research object, Mao studied the change law of soil modulus of resilience with water content and freeze-thaw cycles. It was found that water content in subgrade is the most important factor affecting the modulus of resilience, and the freeze-thaw cycle is a secondary one under the conditions of optimum water content and degree of compaction [16]. Chen used the self-developed multifunctional subgrade freeze-thaw cycle test device to paper the repeated freeze-thaw characteristics of subgrade soil in a seasonal frozen area under different water compensating state. Analyzing the influence of initial degree of compaction, freeze-thaw cycle times, and water content on modulus of resilience, he revealed the change law of modulus of resilience and gave the modulus reduction coefficient of subgrade in a seasonal frozen area [17]. Wang used the selfdeveloped remote temperature and humidity acquisition system for subgrade in a seasonal frozen area to monitor the temperature and humidity status of subgrade, revealing the temperature and humidity change law of pavement structure in the seasonal frozen area. Meanwhile, he also analyzed the influence of water content, freeze-thaw cycle times, and degree of compaction on subgrade strength in a seasonal frozen area through the indoor freeze-thaw cycle test [18]. To sum up, although there are many research studies on the asphalt concrete pavement in the seasonal frozen area, most of them are about the influence of environmental factors on pavement materials and subgrade physical characteristics in the frozen soil area. ere is no relevant research on the influence of temperature and humidity in the seasonal frozen area on the inversion modulus of the semirigid base pavement structure layer. erefore, this paper introduces the SVM method, optimizes the inversion of the asphalt concrete pavement structural layer modulus by the firefly algorithm considering the effect of temperature and humidity, and analyzes the influence of temperature and humidity on the inversion modulus of semirigid base asphalt pavement in the seasonal frozen area. e research results have guiding significance to ensure the safety performance of the asphalt concrete pavement structure in the seasonal frozen area and extend its service cycle. is paper takes FWD load as the research object and selects the depth direction 6m and radial direction 6m as calculation elements. According to the theory of pavement structure mechanics, the displacement, stress, and strain of the axisymmetric structure have nothing to do with the value of θ, but only with the value of R and Z. erefore, the rectangular ring element is used as the research object of the axisymmetric structure, and only the section on the coordinate R and z plane needs to be studied. e finite element model is shown in Figure 1.

Calculation of Pavement Dynamic Response.
According to Hamilton's variation principle, the dynamic balance equation considering damping effect is [19] [ where [M] is the global mass matrix; € δ is the global ac- { }is the equivalent node load vector; α and β are the damping coefficients, α � ξω 1 and β � ξ/ω 1 ; ω 1 is the structural fundamental frequency; ξ is the mode damping ratio. e calculation range of the model, boundary treatment, and characteristics of the rectangular ring element are known from [20,21].
In this paper, the Wilson-θ direct integration method [20,21] is selected to calculate the dynamic balance equation, and θ � 1.4 is taken in the calculation. e specific steps are as follows.

Initial Calculation.
Calculate the following constants: e effective stiffness matrix is established: e effective stiffness matrix is decomposed:

Calculation of Each Time
Step. Calculate the effective load vector: e equation of motion is solved by Calculate acceleration, velocity, and displacement of t 1 + Δt time: Based on the above solution method, the dynamic response calculation program of the semirigid base pavement structure is compiled in this paper.

Example Analysis
(1) In a three-layer semirigid base pavement system, the elastic modulus of each layer is, respectively, 1200 MPa, 1400 Mpa, and 60 MPa, the thickness 18 cm, 32 cm, and 550 cm, the Poisson's ratio 0.35, 0.30, and 0.40, and the density 2400 kg/m 3 , 2300 kg/ m 3 , and 1800 kg/m 3 . In this paper, Dynatest8000 FWD is used as the applied load. In order to simplify the calculation process, the load application curve of FWD is simplified to a sine semiperiodic function with the peak value of 721 kpa, the action time of 32 ms, and the loading radius of 0.15 m. e load time history curve and transducer distribution position of the simplified FWD are shown in Figure 2 and Table 1.
e calculated deflection results in this paper are basically consistent with the measured deflection basin data at different measuring points with small error and high accuracy from Figure 5, which can be applied to actual asphalt concrete pavement engineering.
Eight-node linear hexahedron element (C3D8R) is selected as the calculation element of the finite element model. e boundary conditions are set as follows: there is no X-direction displacement in the left-right direction and no Y-direction displacement in the front-back direction. All nodes on the bottom are completely fixed, and the layers are completely continuous. e finite element model is shown in Figure 3. It can be seen from Figure 4 that the calculation results at different times obtained by applying the Advances in Civil Engineering method in this paper are basically consistent with those of ABAQUS, which proves that the dynamic response calculation program has high accuracy.
(2) In order to further verify the accuracy of the program, based on FWD-measured deflection basin data and the position of the measuring points in [22], the calculation program of this paper is used to calculate the corresponding data. e calculated results are compared with the measured data, as shown in Figure 5.

Modulus Inversion of Asphalt Concrete Pavement Structural Layer.
is paper introduces the SVM model in the modulus inversion of the asphalt concrete pavement structure layer and uses the SVM model instead of the numerical model to carry out the inversion of the modulus, considering the disadvantages of the previous inversion methods of pavement modulus. Based on the structural risk minimization criterion, SVM can balance the relationship between training error and complexity of learning machine. It can improve inversion efficiency in the application of engineering practice. SVM can solve the "overlearning" problem faced by the neural network well and make the calculation results have better generalization ability. It has been widely used in the fields of automatic control, pattern recognition, and rock mechanic parameter inversion [23]. e inversion process of the pavement structural modulus using the SVM model is as follows [24].
whereX � (x 1 , x 2 , . . . , x m )is the parameter of the pavement structural layer to be inverted and Y is the deflection value of a measuring point measured in the field. e number of SVM(x 1 , x 2 , . . . , x m ) is the same as the number of deflection transducer, and n SVM models correspond to n deflection values. e SVM model y(X) which reflects the mapping relationship between structural layer modulus and pavement deflection is established by enough learning samples: where y(X) is a set of displacements of a line corresponding to parameter X � (x 1 , x 2 , . . . , x m ); l is the number of learning samples; k(x, x i ) is the kernel function; α, α * is the Lagrange multiplier of the SVM algorithm. From the SVM model, it can be seen that the parameter α, α * and the offset b are important parameters for creating the SVM model, and their values can be obtained by solving the constrained quadratic programming problem: An efficient optimization algorithm is used to solve the parameters α, α * and the offset b to establish the standard SVM model. e optimal solution is constantly searched in the range of modulus. When the error between the deflection value calculated by the SVM model and the measured deflection reaches the minimum, a set of modulus parameters corresponding to the calculated deflection can be determined as the real modulus of the asphalt concrete pavement structure layer. e inversion process is completed.

Basic Steps of SVM Inversion Optimization Based on Firefly Algorithm.
e selection and optimization of model parameters have a great influence on the establishment of SVM model. e reasonable model parameters are one of the important factors to ensure the accuracy and reliability of the inversion results of the SVM model. e SVM model has two parameters, the basic parameter of the model (such as penalty factor C and insensitive loss parameter ε) and the parameter related to kernel function (such as the kernel width σ in radial basis function, the constant c in polynomial kernel function, and the exponential function d). Previous methods of parameter selection mainly rely on subjective experience and grid search, so it is difficult to find the best parameters [25,26]. In recent years, scholars use the heuristic intelligent search algorithm with global optimization ability to select parameters of SVM and get good results. erefore, the firefly algorithm is used to obtain the optimal parameters of the SVM model in this paper.

Introduction of Firefly
Algorithm. Suppose n fireflies are randomly and evenly distributed in the solution space and each firefly is given a unique initial brightness value. e function value of the firefly's position will determine the brightness of the firefly. e relationship between the position and the brightness of the firefly is the better the position, the higher the brightness. Each firefly has its own dynamic decision domain r i d , in which the firefly moves towards the brighter firefly, forming a domain set. en, the position, brightness, and decision domain of the firefly are updated [27]. e main phases of the firefly algorithm are as follows: (1) e updating phrase of the dynamic decision domain: where r i d (t + 1) is the t + 1 generation decision range of the t generation firefly i; r s is the perception range of the firefly; n t is the domain threshold which can control the number of fireflies; β is the constant  In this phase, each firefly first selects the firefly whose brightness is higher than itself within the range of its decision domain to form the domain set. e formula is where ‖x j (t) − x i (t)‖ is the Euclidean distance between two fireflies; x j (t)is the firefly position; l i (t) is the fluorescein value of the t iteration. e formula of the moving probability is as follows: where p ij is the probability of firefly i moving to firefly j.
After the composition of the domain set, the position is updated according to the following formula: (3) e updating phrase of updating fluorescein value: where l i (t) represents the fluorescein value of the firefly i of the t generation; p is the vanishing rate; c is the parameter to measure the fitness function value; J(x i (t)) is the fitness function value. e main steps of the firefly algorithm are as follows: (1) Random initialization of firefliesi(i � 1, 2, . . . , n)

Optimization of SVM Model Parameters Based on
Firefly Algorithm. is paper adopts the firefly algorithm to optimize the SVM model and then solves the problem of "overlearning" and "underlearning" when the SVM model is not the optimal parameter. e specific steps of the firefly algorithm to optimize SVM model parameters are as follows: (1) Giving initial values to the parameters of each firefly, and setting the initial position of the firefly with the range of initial position parameters (C min ∼ C max , σ min ∼ σ max ) and the maximum number of iterations; (2) Calculating brightness value as brightness update; (3) Obtaining the domain set of each firefly and randomly selecting the firefly i in the set as its moving direction; (4) Calculating the training value of SVM after position update. If the training value of the updated SVM is better than the previous one, the position will be replaced, otherwise it will not change; (5) Updating the dynamic decision domain; (6) Calculating the above process in cycle. If the search times reach its maximum, the cycle process will be interrupted. e SVM will be trained according to the optimal solution in the process until the end condition is satisfied.  Table 2. e range of modulus parameters of each layer is 1000-4000 MPa for surface course, 400-3000 MPa for base course, and 30-150 MPa for subgrade. e modulus parameters of each layer are substituted into the dynamic response calculation program to obtain the corresponding pavement deflection. ese data form a sample pair, and thus 12 test samples and 100 learning samples are obtained.

Example Analysis
It can be seen from Figure 6 that the modulus of asphalt concrete pavement surface course, base layer, and subgrade obtained by this method is between the regression model algorithm and the BP model algorithm. Furthermore, the inversion results have good stability by this method. e three methods are basically consistent with each other, which proves that the method in this paper has high accuracy and completely meets the practical engineering calculation requirements.
is paper adopts standard SVM in which the loss function is an insensitive loss function and the kernel function is a radial basis kernel function (RBF) [28].
rough the training of learning samples and test 6 Advances in Civil Engineering

Advances in Civil Engineering
samples, the optimal model parameters and the offset b values of SVM corresponding to each deflection transducer are obtained, as shown in Table 3.
After obtaining the optimal SVM model, the SVM model instead of the numerical model combines with the firefly algorithm to carry out the inversion. e inversion results are shown in Table 4. e results reveal that the average absolute value of relative error of surface course modulus inversion results is 0.97%, and the maximum value is of 2.16%. e average absolute value of relative error absolute value of the base course modulus inversion result is 0.82%, and the maximum value is of 2.31%. e average absolute value of relative error absolute value of the subgrade structure layer modulus prediction result is 1.91%, and the maximum value is of 3.73%. Based on the results, the SVM model has high inversion efficiency and accuracy, which can fully meet the practical engineering demands.
(2) In order to further verify the accuracy of the inversion program, the inversion method proposed in this paper is used to carry out the modulus inversion of each structural layer of pavement with the measured deflection data [22]. e results are compared with those of the two methods in [22]. e results are presented in Table 5.

Temperature Correction of Inversion Surface Modulus in Seasonal Frozen
Area. For asphalt concrete pavement, the asphalt mixture is a kind of temperature sensitive material. Its strength, stiffness, and deformation resistance will change with temperature. When the asphalt concrete pavement temperature increases, the strength and stiffness of the asphalt mixture will decrease, causing the asphalt concrete pavement to easily undergo high-temperature shear failure under the vehicular load. When the asphalt concrete pavement temperature decreases, the deformation capacity of the asphalt mixture decreases, causing the asphalt concrete pavement to easily undergo low-temperature shrinkage cracking. Pavement modulus which can reflect the strength and stiffness of pavement is the main parameter of pavement structure design. erefore, the analysis of the relationship between asphalt concrete pavement surface modulus and temperature is of great significance to the evaluation of pavement structure bearing capacity in the seasonal frozen area.
In this paper, the temperature field model of Zhang is selected as the research object [29]. e model is as follows.
Pavement temperature along the depth direction: where H is the distance between a certain point within the asphalt concrete pavement and the surface; T H is the temperature at the depth of h; T 60 is the temperature at 60 cm of subgrade; T S is the asphalt concrete pavement temperature or actual pavement temperature; a ∼ e is the undetermined regression coefficient. According to the hot linear elastic layer system theory and in consideration of the structural characteristics and stress state of asphalt pavement, the following assumptions are set: (1) Each structural layer of the asphalt concrete pavement is isotropic and homogeneous thermoelasticity material (2) It is assumed that only the material parameters of the surface course are affected by temperature and the material parameters of the base course and subgrade are regarded as constant (3) e structural layers are completely continuous Based on the above hypothesis analysis, the following calculation examples are constructed for analysis: a threelayer semirigid base pavement system in the seasonal frozen area is selected, including 18 cm asphalt concrete surface course and 32 cm cement-stabilized macadam base and subgrade. Poisson's ratio of each layer in the asphalt concrete pavement is, respectively, 0.35, 0.30, and 0.40 and the density is 2400 kg/m 3 , 2300 kg/m 3 , and 1800 kg/m 3 . e relationship between the inversion modulus of asphalt concrete pavement surface course and temperature can be obtained by embedding the temperature field prediction model into the inversion program. e calculation results are shown in Figure 7 and Table 6. e inverse modulus of asphalt pavement has a good exponential relationship with the asphalt concrete pavement temperature from Figure 7. e inverse modulus of the surface course will decrease with the increase in temperature. e correlation formula of E-T is obtained by regression analysis: e correlation coefficient R 2 between them is equal to 0.9896. e temperature correction coefficient K is defined as the ratio of the modulus E of asphalt pavement under T temperature and the modulus E S under standard temperature T S . T S is taken as 20°C, that is,

Advances in Civil Engineering
where K is the correction coefficient of the inversion modulus of the surface course and T is the expressed temperature.

Humidity Correction of Inversion Subgrade Modulus in
Seasonal Frozen Area. Temperature in different seasons fluctuates greatly in the seasonal frozen area. Water in the subgrade turns into ice at the temperature below 0°C, which changes the properties of subgrade structure materials as well as the position and volume of soil particles. In contrast, ice turns into water at the temperature above 0°C, which causes the volume of soil particles to become smaller. Both factors will result in repeated frost heaving and thawing settlement deformation of pavement. erefore, considering    1  2342  2198  2486  2599  2466  2732  147  155  138  2  4300  4177  4423  3850  3737  3962  136  130  141  3  5187  5225  5149  4495  4463  4526  222  210  234  4  3938  4003  3872  3819  3679  3957  150  156  144  5  3009  2916  3102  2157  2075  2239  178  185  170  6  2694  2638  2749  3159  3077  3241  154  146  161  7  3427  3481  3372  2629  2595  2663 e correlation between subgrade modulus and humidity can be obtained by embedding the model into the inversion program. e calculation results are shown in Figure 8 and Table 7. Figure 8 shows that the inversion modulus of subgrade decreases with the increase of water content, which is an exponential relationship. e correlation formula of inversion modulus of subgrade and water content is obtained from regression analysis as follows: e correlation coefficient R 2 between them is equal to 0.9833. e humidity correction coefficient K E is defined as the ratio between the modulus E of subgrade with water content w and the modulus E S under the condition of optimal water content w S . e optimal water content w S is 16%, that is,  where K E is the correction coefficient of inversion modulus of subgrade and w is the water content.

Conclusion
In this paper, in view of the existing problems in the inversion of the asphalt concrete pavement structure layer modulus, the support vector machine model is introduced into the asphalt concrete pavement modulus inverse analysis to study the influences of the temperature and humidity in the seasonal frozen area on the inversion modulus. e research conclusions are as follows: (1) After inversion by using the support vector machine model optimized by the firefly algorithm, the average of relative error absolute value of the pavement structural layer modulus is within 2% and the maximum value is within 4%, which meets the accuracy requirements. (2) Considering the influence of temperature and humidity on the inversion modulus, the inversion modulus of pavement surface layer will decrease with the increase of temperature, and the inversion modulus of subgrade will decrease with the increase of water content.
(3) According to the inversion results, the modified formulas for the inversion modulus of the semirigid base pavement surface layer and subgrade varying with temperature and humidity are established: K � e −0.033(T− 20) and K E � e −0.049(w−16) .

Data Availability
e data used to support the findings of this study are not available.

Conflicts of Interest
e authors declare that they have no conflicts of interest.  Advances in Civil Engineering 11