An Interval Parameter Calculation Method of Asphalt Pavement Structure Design Based on Point Numerical Algorithm

In the asphalt pavement structure design method, the structural analysis and design are generally performed in the form of point values. However, determining the point value form of design parameters based on the statistical analysis theory cannot fully reflect the complex properties such as variability and uncertainty of parameters. In order to further improve the reliability and practicability of pavement design parameters, in this article, we have introduced the interval number representation that can better reflect the complex nature of parameters; but the interval number algorithm is too complicated and common calculation tools and software are difficult to adopt, which limit the wide application of interval analysis to some extent. The article analyzes the algorithm of interval numbers, focusing on the analysis of interval numbers of unary and binary functions. In this way, the point number operation can be used to obtain the interval number result of the function consistent with the interval number algorithm, which avoids the complicated interval number operation process and the interval expansion. The point numerical function algorithm of interval numbers is verified by design parameters and the calculation of asphalt pavement structure such as axle load conversion, cumulative equivalent axis calculation, calculation of foundation layer tensile stress of each structure layer, calculation of mixture penetration strength, fatigue cracking check of asphalt mixture layer, permanent deformation check, and vertical pressure strain test of roadbed top surface. In conclusion, this research provides a simple and easy way to implement the application of mathematical tools for interval analysis, which is suitable for direct use for existing point numerical calculation tools and software.


Introduction
In the current theory and method of pavement structure design, the pavement structure analysis and design are carried out with representative values, i.e., point value mode [1]. Pavement structure design parameters are mostly expressed in the form of point values, such as traffic load parameters (axial load spectrum), fatigue equation parameters, permanent deformation of the asphalt layer, and so on. Pavement structure design parameters are generally obtained by means of testing. Due to the representative problem of sampling of materials, the accuracy and stability of the test equipment, and the objective existence of test data measurement error, the point value of parameters is determined according to the statistical analysis theory and method. It cannot fully reflect the complex nature of parameters such as volatility and uncertainty.
From the quantum effect and the uncertainty principle, the "point true value" of the material parameters of the pavement structure cannot be finalized according to its nature; however, pavement design specifications, including asphalt and cement pavement, and design parameters are expressed in terms of point values. From the three stages of pavement engineering design, construction, and acceptance, some indicators are given a range of intervals in acceptance, but in the design and construction process, point values are given. e link between the two is the use of the statistical analysis theory. In fact, the design and construction processes are generally erratic. erefore, the introduction to the interval number representation instead of the point value is an objective requirement for characterizing the complex nature of the material parameters of the pavement structure and is also consistent with the actual situation where the pavement design deviates and the construction has errors.
Since the emergence of the interval analysis theory, it has received extensive attention and a large number of books have been published on it [2][3][4].
e articles on interval analysis are increasing, and interval analysis begins to move from theory to practice and plays its role in more and more fields. e important applications of interval analysis in the world include material mechanics, structural mechanics, and so on, which began in the 1980s [5][6][7].
In the field of engineering, Su et al. [8] studied the reliability analysis method of response surface based on interval variables; Tang et al. [9] proposed the fatigue crack checking process of cement-stabilized macadam is analyzed by using measurement uncertainty and interval analysis theory in order to improve the efficiency of fatigue test data; Su et al. [10] incorporated a new uncertainty analysis method for engineering structures, the interval analysis method, into the sensitivity analysis of engineering parameters, obtained a new sensitivity analysis method for engineering parameters, and further broadened the theory of interval analysis methods. Application areas: Huang [11] considered the basic parameters that affect the particle size improvement of the roadbed thickness as interval variables. Using interval analysis to study the thickness can better overcome the complexity of the environment and related parameters and the problem of limited information; in the study by Wang et al. [12], starting from data processing, the concept of interval number is given and a multiattribute decision-making model for interval numbers is proposed. e model provides a standardized processing method for the interval number decision matrix. Wang and Tang [13] discussed the influence of different subgrade soil parameters on the subgrade elastic modulus and analyzed the subgrade elastic modulus under equilibrium humidity; Zhang et al. [14] proposed the interval fuzzy evaluation analysis method of highway subgrade stability in a karst area; Liu and Han [15] used fuzzy reliability in the reliability design of the pavement structure, which can further consolidate the performance of the pavement and extend the service life of the pavement; Yu et al. [16] applied the interval analysis model to the stability analysis of slopes. By citing the idea of interval mathematics, the interval limit equilibrium method was used to derive the minimum safety factor interval of the slope, and on this basis, the slope was nonprobabilistically reliable. Degree analysis: Tang and Zheng pointed out that there are two kinds of ill-posed problems in the parameter inversion of the Duncan-Zhang model of the soil and proposed the concepts of interval well-defined and interval ill-posed in parameter inversion; and given the definitions and analysis of interval, the basic theoretical issues of geotechnical engineering were studied and the basic theoretical framework of interval analysis soil mechanics was focused on [17,18]; Xie et al. [19] introduced the interval analysis method to analyze the fatigue test and data of semirigid base materials and established the fatigue strength equation of semirigid materials in the form of interval parameters; Impollonia and Muscolino [20] proposed a method of evaluating the static response of structures of spaced axial stiffness. An interval analysis method of fatigue crack propagation (FCG) life prediction was proposed by Long et al. [21]. Liu et al. [22] proposed a new static response interval uncertainty analysis method. Sofa et al. [23] proposed an interval finite element analysis method developed by additional unit intervals in the improved interval analysis framework. Dimarogonas [24] proposed an interval analysis of the vibration system. Donald and Chen explored the method of slope stability analysis [25]. Tonon et al. [26] proposed a stochastic set theory about the range of parameters in rock engineering. Giasi et al. [27] studied the fuzzy reliability of the Aliano slope. Schweiger and Peschl [28] proposed the reliability analysis of the stochastic finite element method in geotechnical engineering. For the most part, human factors and environmental conditions affect the measurement results of test parameters. When these digital fuzzy parameters are analyzed in the form of intervals, the accuracy of the measurement results can be more accurately reflected from the aspects of influence analysis, model optimization calculation, and the improvement of the system's predictive ability, so as to draw a new solution [29][30][31][32][33][34][35][36].
However, due to the complexity of the interval algorithm, special interval analysis calculation tools and software are generally used in the calculation of expressions in the form of interval parameters. Because such computing tools and software are not popular, the types are too small and most of them are written by foreign scholars, which greatly limits the learning and use of domestic scholars. Based on this, the article deeply analyzes the algorithm of interval numbers and focuses on the interval number operation of unary and binary functions. By transforming the function expression and other methods, the interval number operation of the function is changed into the upper and lower endpoint value of the independent interval number. In this way, the point numerical operation can be used to obtain the function interval number result consistent with the interval number algorithm, which avoids the complicated interval number operation process and the interval expansion. It can also provide direct application of interval calculation and analysis for various calculation tools and software that currently use point numerical algorithms.

Interval Number Algorithm
Interval mathematics is a mathematical theory defined on the set of intervals. A continuous subset X � [X, X] on the set of real numbers R is called a real interval. X is the lower endpoint value of the interval number X, and X is the upper 2 Advances in Civil Engineering endpoint value of the interval number X. An important rule of interval number is X ≤ X.
e following are the four arithmetic rules of interval number: Interval addition Interval subtraction Interval multiplication Equation (4) shows that for basic interval operations, interval division with 0 is not allowed.
Interval arithmetic has its own characteristics. For example, to find the trigonometric function Sin In fact, the value of a function of a certain interval variable is the function value field corresponding to the interval variable and is not the interval value obtained by the combination of the endpoint values. For example, the interval variable of Sin[0, π] is [0, π], and the corresponding range is [0, 1] instead of [0, 0], as shown in Figure 1. erefore, it is necessary to deeply analyze the algorithm of interval numbers in interval analysis, and especially the interval number algorithm for general functions, because it has a strong practical significance of the correct use of interval analysis. us,

Point Numerical Operation Method for the Interval Number of Unary Function
Here, whether or not the interval number   [27,36].

Calculation Example
It is noted that Based on the interval number algorithm, the interval number result of f(x) is obtained, which is consistent with the number of intervals obtained by combining the upper and lower endpoint values.

General Conclusion
(1) A general conclusion is now given. Referring to the nature of higher mathematics for function derivatives, the monotonicity of functions is defined as follows: If for any two points x, x ∈ D, when x < x, there is a constant: , the function f(x) is said to be monotonically reduced in D (strictly).
Definition 2 (see Ref. [37]): (2) Set f(x) � ax n + b to avoid the interval width. Related literature [17] has proved that when the independent variable x appears only once in the expression f(x) on the right side, the interval expansion of the function f(x) result is the interval , the argument x in the right expression only appears twice, and the interval number of the function result obtained by taking the upper and lower endpoint values of the argument interval is not applicable.
, the result of interval number function calculation is [6,10], which is consistent with the result of using the interval number function operation rule directly.  [5,6], and this result is consistent with the operation result of the interval  (1)- (4) are the binary-time functions addition, subtraction, multiplication, and division, the specific method of calculating the upper and lower endpoint values is adopted. Investigating equations (1)-(4), the addition rule is easy to understand and remember, but the subtraction method is more complicated.
We can calculate the interval number of interval subtraction according to the following steps: (1) e sign of the subtracted number is unchanged (2) Add a negative sign to the subtraction, according to the number of intervals, and change the upper and lower endpoints of the interval of the subtraction, (1) First case: at this time, no matter if X Y and X Y are big or small; in the four values X Y, X Y, X Y, XY, the smallest is X Y and the largest is XY. In this case, when all the interval parameters are greater than 0, it is very convenient to use the interval endpoint value to calculate the interval number.
Evidently, when the number of the two multiplied intervals is a negative interval, XY ≤ X Y, X Y ≤ X Y, no matter which is X Y and X Y and which is small; in the four values X Y, X Y, X Y, XY, the smallest is XY and the largest is X Y. We can directly use the lower endpoint value of X to multiply the lower endpoint value of Y in order to obtain the upper endpoint value of the interval number; multiply the upper endpoint value of X by the upper endpoint value of Y to obtain the lower endpoint value of the interval number.
e value of the four endpoints X, X, Y, and Y is either positive or negative. Among the four values X Y, X Y, X Y, XY, the judgment of the size relationship between them is complicated by the existence of the sign, such as

Engineering Application Examples.
In general, in engineering applications, some of the parameters including the equations used to fit the series of parameters are mostly positive intervals. For example, in pavement engineering, the uniaxial compressive strength, four-point bending strength, and splitting strength test results of the cement-stabilized macadam base material can be fitted by an equation containing a natural logarithmic function, such as the following equation: In equation (5), S c is strength, t is age, and a, b, and c are all fitting parameters. According to the interval four algorithm, the two interval numbers are added together, and the upper and lower endpoint values can be directly added to obtain the interval number result. e interval number calculation formula of equation (6) is e key is to check whether [b ln(t + c), b ln(t + c)] can be the number of intervals in which the upper and lower endpoint values are calculated to obtain the result.
If Consider a more general case of a binary function: assume f(x, y) � x a y b , where x, y > 0 and a, b > 0. Obviously, f(x, y) in this case-x a and y b are both increasing functions-satisfies the first case of interval multiplication. erefore, the upper and lower endpoint values of x, y and a, b can be directly used to calculate the interval number of f(x, y), which is In equation (7), it is assumed that x, y > 0 and a,  If f(x, y) � x a + y b and x, y > 0, a, b > 0, the number of intervals of f(x, y) can also be calculated directly by using the upper and lower endpoint values of x, y and a, b, which is e division of the binary function can be changed into the multiplication of binary function by changing the formula. Referring to the above, the interval number results can be directly calculated by using the upper and lower endpoint values. For example, in the early fatigue damage research, the Miner linear model is an earlier theory, and its basic structure is as follows (9): where N f is the fatigue life and N is the number of loads. Equation (9) is rewritten into equation (10): where [N f , N f ] of equation (9) is rewritten as

Point Number Operating Method of the Interval Number of Multivariate Functions
For ternary function expressions, let f(x, y, z, · · ·) � x a y b z c , · · · , x, y, z, · · · > 0, a, b, c · · · > 0. Obviously, in this case, f(x, y, z, · · ·) is an increasing function, which satisfies the first case of interval multiplication, so the upper and lower endpoint values of x a , y b , z c , · · · can be directly used to calculate the interval number of f(x, y, z, · · ·), which is f(x, y, z, · · ·), f(x, y, z, · · ·) � x a y b z c , · · · , x a y b z c , · · · . (11)

If
f(x, y, z, · · ·) � x a + y b + z c + · · ·, and x, y, z, · · · > 0, a, b, c, · · · > 0, you can also directly use the upper and lower endpoint values of x, y, z, · · · , a, b, c, · · · to calculate the interval number of f(x, y, z, · · ·), which is f(x, y, z, · · ·), f(x, y, z, · · ·) � x a + y b + z c + · · · , x a + y b + z c + · · · . (12) e division of multivariate functions can be transformed into the multiplication of multivariate functions by changing their formulas. Referring to the above, the results of the interval numbers can be directly calculated by using the upper and lower endpoints.

Examples of Asphalt Pavement Structure
Design and Checking Calculation e basic data of asphalt pavement structure design, namely traffic volume data, are very different from the theoretical prediction. It is unreasonable to use the point numerical traffic volume to design the asphalt pavement structure. It is possible to consider the number of traffic intervals to design the pavement structure.

Axle Load Conversion and Accumulated Equivalent Axis Interval Number Calculation.
e axle load conversion at each level uses the following equation: N 1 is the standard axle load equivalent axle number, times/day; n i is the converted vehicle axle load action number at all levels, times/day; p is the standard axle load, kN; p i is the converted vehicle axle load at all levels, kN; k is the converted vehicle type; c 1 is the axle number coefficient, c 1 � 1 + 1.2(m − 1), and m is axle number. When the shaft spacing is greater than 3 m, the calculation is based on a single axle load. When the shaft spacing is less than 3 m, the axle coefficient should be considered; c 2 is the the wheel group coefficient, and the single-wheel group is 6.4, the twowheel group is 1.0, and the four-wheel group is 0.38.
As we all know, n i , the number of axle loads converted into vehicles at all levels, is obtained directly from the traffic flow observation data. Due to the observation of traffic flow, there are various statistical methods, such as manual, electronic, and toll stations. ere are more or less statistical errors in these statistical methods. e traditional calculation of n i in the form of point values cannot fully reflect the complex characteristics of traffic flow. If the statistical error of a certain class of axle load traffic (vehicle/day) is ±1%, then the traffic volume interval parameter is[n i − n i × 1%, n i + n i × 1%].
According to the point numerical operation method of the interval number of the unary function in the second section, equation (13) directly calculates the lower endpoint value of the equivalent axis sub-N 1 interval number of the standard axle load using the lower endpoint value n i − n i × 6 Advances in Civil Engineering 1% and calculates the standard using the upper endpoint value n i + n i × 1%. e upper endpoint value of the number of equivalent N 1 subintervals of the axle load, assuming the axle number coefficient c 1 is 1, the wheel group coefficient c 2 is 1, p i is 55.1 kN, p is 100 kN of the standard axle load, vehicle type k is 1, and n i is 1000 times a day, is calculated as follows: us, the formula for calculating the number of intervals for the cumulative equivalent axis N e is N e is the cumulative equivalent axle number; η is the lane coefficient value according to specifications; t is the design years; c is the design years, the annual average growth rate of traffic volume forecast; and N 1 is the standard axle load equivalent axle number, next/day. Now, if

Interval Number Calculation of Structural Coefficient of Tensile Strength.
e road grade coefficient A c is 1.0, the surface layer is asphalt concrete, the surface type coefficient A s is 1.0, the semirigid base course is 1.0, and the pavement structure type coefficient A b is 1.0. For the asphalt concrete surface, the interval number of the structural coefficient of tensile strength is calculated as follows: e stable aggregates for inorganic binders: For inorganic binder stabilized fine-grained soils: In equations (16)- (18), N e , N e ′ -cumulative equivalent axes. Among them, N e ′ is the cumulative equivalent axis for calculating the bending and tensile stresses at the bottom of the semirigid base course. From Section 2, Example 3, we can see that equations (16)-(18) directly calculate the lower endpoint value K s of the equivalent axle number K s interval of the standard axle load using the lower endpoint value N e , N e ′ and calculate the upper endpoint value K s of the equivalent axle number K s interval of the standard axle load using the upper endpoint value N e , N ′ .

Calculation of the Number of Tensile Stress Intervals in Each Layer of the Material.
e formula for calculating the allowable tensile stress of the pavement structural layer material is σ R is the allowable tensile stress of pavement structural layer materials, MPa; ultimate tensile strength of σ s is the structure layer material, MPa, determined by experiment; and structural coefficient of K s is the tensile strength.
For the asphalt concrete surface layer, if the tensile strength structural coefficient is the interval number, then For the calculation of equation (20), one thing to note is since [K s , K s ] is a divisor, it needs to be rewritten as For the calculation of equation (21), the lower endpoint value σ R of the allowable tensile stress σ R interval number can be directly calculated by using the lower endpoint value 1/K s and the upper endpoint value σ R of the allowable tensile stress σ R interval number can be calculated by the upper endpoint value 1/K s .  Table 1 shows the calculation results of the allowable tensile stress design interval for each structural layer of a road surface.

Interval Number Calculation of Asphalt Mixture Penetration Strength.
e penetration strength of asphalt mixture of inorganic binder stabilized base asphalt pavement, asphalt binder stabilized subbase asphalt pavement, and cement concrete base asphalt pavement should meet the following requirements: Specific explanations of the parameters in equation (23) can be found in formula 5.5.8-1 of Highway Asphalt Pavement Design Code (JTG D50-2017). If the cumulative equivalent axis N e5 is considered as interval number [N e5 , N e5 ], the formula for calculating the interval number of penetration strength is as follows: According to Section 2, equation (24) can directly calculate the lower endpoint value R τs of the penetration strength R τs interval number using the lower endpoint value N e5 and the upper endpoint value R τs of the R τs interval number using the upper endpoint value N e5 .

Fatigue Cracking Check of Asphalt Mixture Layer.
e fatigue cracking life of the asphalt mixture should be calculated according to the pavement structure, and the bottom strain of the asphalt mixture layer is calculated according to the following equation: For a detailed explanation of the parameters in equation (25) According to the interval division calculation rule, combined with the third section of the article and equation (10), equation (26) can be rewritten as At this time, the lower endpoint value N f1 of the fatigue cracking life interval number can be directly calculated by using the lower endpoint value 1/ε a and the upper endpoint value N f1 of the fatigue cracking life interval number can be calculated by using the upper endpoint value 1/ε a .

Verification of Permanent Deformation of Asphalt
According to the second section, equation (29) can directly calculate the lower endpoint value R ai of the i-layer permanent deformation amount R ai interval using the lower endpoint value N e3 and the upper endpoint value R ai of the R ai interval number using the upper endpoint value N e3 . Since the index of N e4 , that is, − 0.21, is negative, it needs to be rewritten as

Vertical Compressive Strain
After this, according to the second section, equation (32) can directly calculate the lower endpoint value ε z of the allowable vertical compressive strain ε z interval number of the top surface of the subgrade using the lower endpoint value 1/N e4 and the upper endpoint value of the number ε z of the ε z interval using the upper endpoint value 1/N e4 .
In addition to the penetration strength of the asphalt mixture, fatigue cracking check, permanent deformation check, vertical pressure, and strain test of the roadbed top surface, the new road asphalt pavement design specification also stipulates subgrade elastic modulus, asphalt surface low temperature cracking index, antifreeze, design parameters, and so on and performs indicator checks such as thickness check and acceptance deflection. e specific analysis can be carried out according to the formula. Referring to the methods described in Sections 2-4, the expression is rewritten or some parameters are simplified, and the design parameters and the calculation formulas for the calculation index suitable for directly calculating the interval number of the upper and lower endpoints of the independent variable are established.

Conclusions
is paper discusses the transformation of the complex interval number function algorithm into the point numerical algorithm and provides a simple and easy interval analysis process and method for the various calculation tools and software that generally use the point numerical algorithm. Taking the calculation and verification of interval parameters of asphalt pavement structure as an example, the theory and method discussed in the article are verified, and the following conclusions are obtained: (2) In the unary function, if the argument in the right expression only appears once, and the monotonicity of the function is consistent within the interval of the independent variable, the upper and lower endpoint values of the argument interval can be used to obtain the interval number of the function. (3) For a binary function, especially when multiplying two independent variables, the positive and negative values of the independent variable interval must be considered. When the first and second cases in Section 3 of this paper are satisfied, the number of intervals can be obtained by calculating the upper and lower endpoint values of the two independent variable intervals. (4) For the multivariate function, the situation is more complicated, and there are many factors to consider. We can learn this from the example in Section 4 of the article. (5) e division of unary, binary, and multivariate functions can be multiplied by the change in the expression, and then the upper and lower endpoint values are used to calculate the interval number as described in this paper. (6) When calculating the interval number of the function using the upper and lower endpoint values of the independent variable, if the coefficient or the index in front of the independent variable is greater than 0, and the function expression only adds and multiplies the two operations, the interval number calculation is relatively simple. However, when the coefficient or exponent of the independent variable is less than 0, and the function expression has two operations of subtraction and division, the interval number calculation is more complicated and needs careful analysis.
In conclusion, this research provides a convenient and easy-to-implement process and method for the popularization of interval analysis theory. In the current analysis and calculation of road mechanics, most of the parameters involved are positive numbers, that is, greater than zero. In addition, the mechanical parameter measurement uncertainty report provided by the national standard "Measurement Uncertainty Evaluation and Representation" (GB/T 27418-2017) is mostly a positive interval value. In practical applications, the data processing software and various other point numerical calculation tools can be directly used to quickly obtain the interval number of parameters, thereby avoiding complicated interval algorithm programming. However, further analysis of the mechanical mechanics of the pavement involves complex operations such as differentiation, integration, and even matrix, especially numerical simulation. In future, how to simplify the calculation of interval numbers in such complex operations is an important issue to be solved.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.