Method, device and medium for constructing master curve of dynamic viscoelastic behavior of asphalt-based materials

Abstract

The present invention discloses a method, device and medium for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials. The method includes: drawing dynamic modulus data at different test temperatures and frequencies under double logarithmic coordinates and selecting a certain test temperature as a reference temperature, determining the interpolation interval of adjacent test temperatures, calculating the horizontal movement range of dynamic modulus test results at adjacent temperatures, determining the shift value of dynamic modulus test results at adjacent temperatures, calculating the shift factor at all test temperatures, drawing a temperature-shift factor curve, determining the shift factor at any target reference temperature by interpolating the shift factor curve and moving it vertically, and after obtaining the shift factor at the target reference temperature, translating the test results under logarithmic coordinates to obtain the master curve and finally fitting the master curve. The present invention can separate the acquisition of the master curve from the fitting of the master curve, has the characteristics of accuracy and simplicity, and can be widely used in the field of road engineering technology.

Description

Technical Field

[0001] The present invention relates to the technical field of road engineering, and in particular to a method, device and medium for constructing a master curve of dynamic viscoelastic behavior of an asphalt-based material. Background Art

[0002] The mechanical behavior of asphalt and asphalt-based materials such as asphalt mortar and asphalt mixture with asphalt as a binder has typical viscoelastic properties. The time-temperature characteristics have an important influence on their mechanical behavior. In order to study their mechanical behavior in a wide time domain and a wide frequency domain, for asphalt-based materials such as materials with simple thermorheological behavior, the time-temperature equivalence principle can usually be used to translate the data in the narrow time domain and narrow frequency domain at different test temperatures to obtain a smooth curve, which is called the "master curve." The existing specification AASHTO R62 standardizes the drawing of the dynamic modulus master curve of asphalt mixture, that is, the shift factor equation and the relationship between the reduced frequency and the frequency are substituted into the master curve model, and then the function is fitted using the standard solution tool or the custom function fitting method. This method is also commonly used to draw the dynamic modulus master curve of asphalt, but the fitting effect of this method is greatly affected by the selection of the master curve model, the shift factor equation and the selection of the initial value of the fitting. In addition, existing research results show that the interlayer shear strength of asphalt mixture and the mechanical damage behavior of asphalt mixture also have time-temperature equivalent characteristics. For this type of unknown mechanical behavior of the master curve model, the standard method cannot be used to draw the master curve, and the manual translation method is usually used, but subjective factors will seriously affect the accuracy of the translation. Summary of the invention

[0003] In order to solve at least one of the technical problems existing in the prior art to a certain extent, the purpose of the present invention is to provide a method, device and medium for constructing a master curve of the dynamic viscoelastic behavior of asphalt-based materials.

[0004] The technical solution adopted by the present invention is:

[0005] A method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials comprises the following steps:

[0006] Plot the dynamic modulus data at different test temperatures and frequencies on a double logarithmic coordinate system, and select the preset test temperature as the reference temperature;

[0007] Determine the interpolation interval of adjacent test temperatures;

[0008] Calculate the horizontal movement range of the dynamic modulus test results at adjacent temperatures based on the interpolation interval;

[0009] Determine the shift value of the dynamic modulus test results at adjacent temperatures based on the horizontal shift range;

[0010] Calculate the shift factors at all test temperatures based on the shift values;

[0011] A temperature-shift factor curve is plotted in a semi-logarithmic coordinate according to the shift factor, and the shift factor at any target reference temperature is determined by interpolating the shift factor curve and moving it vertically;

[0012] After obtaining the shift factor at the target reference temperature, the master curve is obtained, and the master curve or the shift factor is fitted using the master curve model or the shift factor equation to determine the parameters of the master curve model or the shift factor equation.

[0013] Further, the determining of the interpolation interval of adjacent test temperatures includes:

[0014] The dynamic modulus test results of two adjacent test temperatures T1 and T2 are extracted, where the temperature of T1 is lower than that of T2. The interval [D1, D2] is determined as the interpolation interval by the dynamic modulus D1 of the minimum test frequency at temperature T1 and the dynamic modulus D2 of the maximum test frequency at temperature T2.

[0015] Further, the step of calculating the horizontal movement range of the dynamic modulus test results at adjacent temperatures according to the interpolation interval includes:

[0016] Determine the frequency at which the dynamic modulus is D2 at temperature T1 by interpolation The frequency when the dynamic modulus is D1 at temperature T2 make in is the minimum test frequency at temperature T1, is the maximum test frequency at temperature T2; when the data at temperature T1 is fixed, the horizontal movement range of the data at temperature T2 is: [-max(Δf1,Δf2),-min(Δf1,Δf2)].

[0017] Further, determining the shift value of the dynamic modulus test result at adjacent temperatures according to the horizontal movement range includes:

[0018] Select n interpolation points at logarithmic equal intervals in the interpolation interval [D1, D2], divide the determined horizontal movement range into m parts at logarithmic equal intervals, horizontally shift the dynamic modulus data at temperature T2 according to the divided horizontal movement range, and calculate the norm Norm of the angular frequency difference between the n interpolation points of the two curves in the interval [D1, D2] after the translation, and find the horizontal movement value j* that minimizes Norm; where j* is the horizontal shift value of the data at temperature T2 when the data at temperature T1 is fixed.

[0019] Furthermore, the calculation process of the horizontal shift value j* is shown in the following formula:

[0020]

[0021] In the formula, is the frequency of the i-th interpolation point in [D1, D2] at temperature T2 when the norm value is minimum, is the frequency of the i-th interpolation point in [D1, D2] when the horizontal shift value is the j-th value in [-max(Δf1, Δf2), -min(Δf1, Δf2)] at temperature T2.

[0022] Further, the calculation of the shift factors at all test temperatures according to the shift values ​​includes:

[0023] Calculate the shift values ​​of all adjacent temperatures, sum all the shift values ​​obtained between the remaining test temperatures and the reference temperature, and determine the sign of the value according to the difference between the remaining test temperatures and the reference temperature; when the test temperature is lower than the reference temperature, the sign is positive, and when the test temperature is higher than the reference temperature, the sign is negative.

[0024] Furthermore, the method of drawing a temperature-shift factor curve in a semi-logarithmic coordinate according to the shift factor, and determining the shift factor at any target reference temperature by interpolating the shift factor curve and moving it vertically includes:

[0025] After obtaining the shift factor curve at the test temperature as the reference temperature, the shift factor at the target reference temperature is determined using the interpolation method, and the shift factor at the target reference temperature is subtracted from the shift factors at all test temperatures to obtain the final shift factor curve.

[0026] Furthermore, the method further comprises the following steps:

[0027] When the test data have no overlapping area, a linear equation or a polynomial is used to fit the test data respectively, and the overlapping area is obtained by extrapolation.

[0028] Another technical solution adopted by the present invention is:

[0029] A device for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials, comprising:

[0030] at least one processor;

[0031] at least one memory for storing at least one program;

[0032] When the at least one program is executed by the at least one processor, the at least one processor implements the above method.

[0033] Another technical solution adopted by the present invention is:

[0034] A computer-readable storage medium stores a program executable by a processor, wherein the program executable by the processor is used to execute the method described above when executed by the processor.

[0035] The beneficial effect of the present invention is that the present invention can separate the acquisition of the master curve from the fitting of the master curve, thus overcoming the defect that the traditional standard method needs to preset the master curve model and the shift factor equation. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the embodiments of the present invention or the drawings of related technical solutions in the prior art are introduced below. It should be understood that the drawings introduced below are only for the convenience of clearly describing some embodiments of the technical solutions of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0037] Figure 1 is a schematic diagram of shift factor calculation in an embodiment of the present invention;

[0038] Figure 2 is a schematic diagram of vertical movement of the shift factor in an embodiment of the present invention;

[0039] Figure 3 is a master curve of dynamic modulus and phase angle in an embodiment of the present invention;

[0040] Figure 4 is the fitting of the dynamic modulus and phase angle master curve in the embodiment of the present invention;

[0041] Figure 5 It is a flowchart of the steps of a method for constructing a master curve of dynamic viscoelastic behavior of an asphalt-based material in an embodiment of the present invention. DETAILED DESCRIPTION

[0042] The embodiments of the present invention are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present invention, and are not to be construed as limitations of the present invention. For the step numbers in the following embodiments, they are only provided for the convenience of explanation, and the order between the steps is not limited in any way, and the execution order of each step in the embodiment can be adaptively adjusted according to the understanding of those skilled in the art.

[0043] In the description of the present invention, it should be understood that descriptions involving orientations, such as up, down, front, back, left, right, etc., and orientations or positional relationships indicated are based on the orientations or positional relationships shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as a limitation on the present invention.

[0044] In the description of the present invention, "several" means one or more, "more" means more than two, "greater than", "less than", "exceed" etc. are understood as not including the number itself, and "above", "below", "within" etc. are understood as including the number itself. If there is a description of "first" or "second", it is only used for the purpose of distinguishing the technical features, and cannot be understood as indicating or implying the relative importance or implicitly indicating the number of the indicated technical features or implicitly indicating the order of the indicated technical features.

[0045] In the description of the present invention, unless otherwise clearly defined, terms such as setting, installing, connecting, etc. should be understood in a broad sense, and technicians in the relevant technical field can reasonably determine the specific meanings of the above terms in the present invention based on the specific content of the technical solution.

[0046] In view of the problems existing in the prior art, the present invention proposes a method for separating master curve acquisition from master curve fitting based on the definition of the master curve, which can accurately construct the master curve of the mechanical behavior of asphalt-based materials.

[0047] like Figure 5 As shown, this embodiment provides a method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials, comprising the following steps:

[0048] S1. Plot the dynamic modulus data at different test temperatures and frequencies in a double logarithmic coordinate system, and select the preset test temperature as the reference temperature.

[0049] The dynamic modulus data at different test temperatures and frequencies are plotted on a double logarithmic coordinate system and a certain test temperature is selected as the reference temperature.

[0050] S2. Determine the interpolation interval of adjacent test temperatures.

[0051] The dynamic modulus test results of two adjacent test temperatures T1 and T2 are extracted, where the temperature of T1 is lower than that of T2. The interval [D1, D2] can be determined as the interpolation area based on the dynamic modulus D1 of the minimum test frequency at temperature T1 and the dynamic modulus D2 of the maximum test frequency at temperature T2.

[0052] S3. Calculate the horizontal movement range of the dynamic modulus test results at adjacent temperatures according to the interpolation interval.

[0053] Determine the frequency at which the dynamic modulus is D2 at temperature T1 by interpolation The frequency when the dynamic modulus is D1 at temperature T2 make in is the minimum test frequency at temperature T1, is the maximum test frequency at temperature T2. When the data at temperature T1 is fixed, the horizontal movement range of the data at temperature T2 is [-max(Δf1,Δf2),-min(Δf1,Δf2)].

[0054] S4. Determine the shift value of the dynamic modulus test result at adjacent temperatures according to the horizontal movement range.

[0055] Select n interpolation points in the interval [D1, D2] at logarithmic equal intervals, divide the horizontal movement range determined in step 2 into m parts at logarithmic equal intervals, then horizontally move the dynamic modulus data at temperature T2 according to the divided horizontal movement range and calculate the norm Norm of the angular frequency difference of the n interpolation points of the two curves in the interval [D1, D2] after translation, and find the horizontal movement value j* that minimizes Norm. j* is the horizontal shift value of the data at temperature T2 when the data at temperature T1 is fixed. The calculation process can be expressed by the following formula:

[0056]

[0057] As an optional embodiment, the number n of interpolation points and the number m of logarithmically equally spaced divisions of the horizontal movement range are greater than or equal to 15 and 50, respectively.

[0058] S5. Calculate the shift factors at all test temperatures based on the shift values.

[0059] Repeat steps S2 to S4, calculate the shift values ​​of all adjacent temperatures and sum all the shift values ​​calculated between the remaining test temperatures and the reference temperature, and determine the sign of the value according to the difference between the remaining test temperatures and the reference temperature. When the test temperature is lower than the reference temperature, a positive sign is used, and when it is higher than the test temperature, a negative sign is used.

[0060] S6. Draw a temperature-shift factor curve in a semi-logarithmic coordinate system according to the shift factor, and determine the shift factor at any target reference temperature by interpolating the shift factor curve and moving it vertically.

[0061] A temperature-shift factor curve is drawn, and the shift factor at any target reference temperature is determined by interpolating the shift factor curve and moving it vertically. The specific process is as follows: after obtaining the shift factor curve at the test temperature as the reference temperature, the interpolation method is used to determine the shift factor at the target reference temperature. Finally, the shift factor at the target reference temperature is subtracted from the shift factor at all test temperatures to obtain the final shift factor curve.

[0062] S7. After obtaining the shift factor at the target reference temperature, obtain a master curve, fit the master curve or the shift factor using a master curve model or a shift factor equation, and determine the parameters of the master curve model or the shift factor equation.

[0063] After obtaining the shift factor at the target reference temperature, the test results (dynamic modulus and phase angle) can be translated in logarithmic coordinates to obtain the master curve. After obtaining the master curve, the master curve or shift factor is fitted using the master curve model or shift factor equation, and the parameters of the master curve model or shift factor equation are finally determined.

[0064] The above steps are for the case where the test data of adjacent test temperatures have overlapping areas on the y-axis. When the test data have no overlapping areas, linear equations or polynomials are used to fit the test data respectively, and the overlapping areas are obtained by extrapolation, and then the above steps are used for calculation.

[0065] The above method is explained in detail below with reference to the accompanying drawings and specific embodiments.

[0066] The data used in this example are frequency scanning data of Shell 70# base asphalt at different temperatures, where the test temperatures are 4, 16, 28, 40, 52, 64 and 76°C, the test frequency range is 0.1 to 30 Hz, and 16 test points are taken at equal intervals in logarithms. The master curve is obtained by the following steps.

[0067] (1) The dynamic modulus data at different test temperatures and frequencies were plotted on a double logarithmic coordinate system and 4°C was selected as the reference temperature.

[0068] (2) Determine the interpolation interval of adjacent test temperatures. Extract the dynamic modulus test results of two adjacent test temperatures T1 and T2, where the temperature of T1 is lower than that of T2. The interval [D1, D2] can be determined as the interpolation area based on the dynamic modulus D1 at the minimum test frequency at temperature T1 and the dynamic modulus D2 at the maximum test frequency at temperature T2. Calculate the horizontal movement range of the dynamic modulus test results at adjacent temperatures. Determine the frequency at which the dynamic modulus at temperature T1 is D2 by interpolation. The frequency when the dynamic modulus is D1 at temperature T2 make in is the minimum test frequency at temperature T1, is the maximum test frequency at temperature T2. When the data at temperature T1 is fixed, the horizontal movement range of the data at temperature T2 is [-max(Δf1,Δf2),-min(Δf1,Δf2)]. The optimization algorithm is used to determine the shift value of the dynamic modulus test results at adjacent temperatures. Select 25 interpolation points at logarithmic intervals in the interval [D1,D2], divide the determined horizontal movement range into 100 parts at logarithmic intervals, and then horizontally shift the dynamic modulus data at temperature T2 according to the divided horizontal movement range and calculate the norm Norm of the angular frequency difference of the two curves at the n interpolation points in the interval [D1,D2] after the translation, and find the horizontal movement value j* that minimizes Norm. j* is the horizontal shift value of the data at temperature T2 when the data at temperature T1 is fixed. The calculation process is expressed by the following formula, and the specific calculation process is as follows: Figure 1 As shown, Figure 1 (a) is the dynamic modulus test data before translation, Figure 1 (b) is a schematic diagram of the determined translation range.

[0069]

[0070] (3) Calculate the shift factor at all test temperatures. The specific process is to repeat the above process to obtain the shift values ​​of all adjacent temperatures and sum all the shift values ​​obtained between the remaining test temperatures and the reference temperature and determine the sign of the value according to the difference between the remaining test temperatures and the reference temperature. When the test temperature is lower than the reference temperature, the sign is positive, and when it is higher than the test temperature, the sign is negative.

[0071] (4) Plotting a temperature-shift factor curve in a single logarithmic coordinate system, and determining the shift factor at any target reference temperature by interpolating the shift factor curve and moving it vertically. The specific process is as follows: after obtaining the shift factor curve at the test temperature as the reference temperature, use the interpolation method to determine the shift factor at the target reference temperature. Finally, subtract the shift factor at the target reference temperature from the shift factors at all test temperatures to obtain the final shift factor curve. The calculation process is as follows: Figure 2 shown.

[0072] (5) After obtaining the shift factor at the target reference temperature, the test results (dynamic modulus and phase angle) can be translated in logarithmic coordinates to obtain the master curve. The dynamic modulus master curve and phase angle master curve for the target reference temperatures of 5°C, 15°C and 25°C are shown in Figure 2. Figure 3 As shown, Figure 3 (a) is the master curve of dynamic modulus at 5℃, 15℃ and 25℃ obtained by translation. Figure 3 (b) is the phase angle master curve at 5℃, 15℃ and 25℃ obtained by translation.

[0073] (6) After obtaining the master curve, the CAM model is used to fit the dynamic modulus master curve and phase angle master curve with a target reference temperature of 25°C. The dynamic modulus master curve model and phase angle master curve model of the CAM model are shown in the following formulas. The fitting effect is shown in Figure 4 shown.

[0074] Dynamic modulus:

[0075] Phase Angle:

[0076] After data fitting is completed, it can be used to analyze the rheological behavior of asphalt-based materials or obtain the Prony series for numerical simulation.

[0077] In summary, the method of this embodiment realizes the separation of obtaining the master curve and fitting the master curve, overcomes the defect that the traditional standard method needs to preset the master curve model and the shift factor equation, and has the characteristics of accuracy and simplicity.

[0078] This embodiment also provides a device for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials, comprising:

[0079] at least one processor;

[0080] at least one memory for storing at least one program;

[0081] When the at least one program is executed by the at least one processor, the at least one processor implements Figure 5 The method shown.

[0082] A device for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials in this embodiment can execute a method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials provided by an embodiment of the method of the present invention, can execute any combination of implementation steps of the method embodiment, and has the corresponding functions and beneficial effects of the method.

[0083] The present application also discloses a computer program product or a computer program, which includes a computer instruction stored in a computer-readable storage medium. A processor of a computer device can read the computer instruction from the computer-readable storage medium, and the processor executes the computer instruction, so that the computer device executes Figure 5 The method shown.

[0084] This embodiment also provides a storage medium storing instructions or programs that can execute a method for constructing a main curve of dynamic viscoelastic behavior of asphalt-based materials provided by an embodiment of the method of the present invention. When the instructions or programs are run, any combination of implementation steps of the method embodiment can be executed, and the corresponding functions and beneficial effects of the method can be obtained.

[0085] In some selectable embodiments, the function / operation mentioned in the block diagram may not occur in the order mentioned in the operation diagram. For example, depending on the function / operation involved, the two boxes shown in succession can actually be executed substantially simultaneously or the boxes can sometimes be executed in reverse order. In addition, the embodiment presented and described in the flow chart of the present invention is provided by way of example, for the purpose of providing a more comprehensive understanding of technology. The disclosed method is not limited to the operation and logic flow presented herein. Selectable embodiments are expected, wherein the order of various operations is changed and the sub-operation of a part for which is described as a larger operation is performed independently.

[0086] In addition, although the present invention is described in the context of functional modules, it should be understood that, unless otherwise specified, one or more of the functions and / or features described may be integrated into a single physical device and / or software module, or one or more functions and / or features may be implemented in separate physical devices or software modules. It is also understood that a detailed discussion of the actual implementation of each module is unnecessary for understanding the present invention. More specifically, in view of the properties, functions, and internal relationships of the various functional modules in the device disclosed herein, the actual implementation of the module will be understood within the conventional skills of the engineer. Therefore, those skilled in the art can implement the present invention set forth in the claims without excessive experimentation using ordinary techniques. It is also understood that the specific concepts disclosed are merely illustrative and are not intended to limit the scope of the present invention, which is determined by the full scope of the appended claims and their equivalents.

[0087] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc., which can store program codes.

[0088] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in conjunction with such instruction execution systems, devices or apparatuses. For the purposes of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in conjunction with such instruction execution systems, devices or apparatuses.

[0089] More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk case (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be a paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium, followed by editing, deciphering or, if necessary, processing in another suitable manner, and then stored in a computer memory.

[0090] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0091] In the above description of this specification, the description with reference to the terms "one embodiment / example", "another embodiment / example" or "certain embodiments / examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner.

[0092] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the claims and their equivalents.

[0093] The above is a specific description of the preferred implementation of the present invention, but the present invention is not limited to the above embodiments. Those skilled in the art may make various equivalent modifications or substitutions without violating the spirit of the present invention. These equivalent modifications or substitutions are all included in the scope defined by the claims of this application.

Claims

1. A method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials, characterized in that: The following steps are involved: Plot the dynamic modulus data at different test temperatures and frequencies on a double logarithmic coordinate system, and select the preset test temperature as the reference temperature; Determine the interpolation interval of adjacent test temperatures; Calculate the horizontal movement range of the dynamic modulus test results at adjacent temperatures based on the interpolation interval; Determine the shift value of the dynamic modulus test results at adjacent temperatures based on the horizontal shift range; Calculate the shift factors at all test temperatures based on the shift values; A temperature-shift factor curve is plotted in a semi-logarithmic coordinate according to the shift factor, and the shift factor at any target reference temperature is determined by interpolating the shift factor curve and moving it vertically; After obtaining the shift factor at the target reference temperature, a master curve is obtained, and the master curve or the shift factor is fitted using a master curve model or a shift factor equation to determine the parameters of the master curve model or the shift factor equation; The step of determining the interpolation interval of adjacent test temperatures includes: Extract the dynamic modulus test results of two adjacent test temperatures T1 and T2, where the temperature of T1 is lower than that of T2, and determine the interval [D1, D2] as the interpolation interval by the dynamic modulus D1 of the minimum test frequency at temperature T1 and the dynamic modulus D2 of the maximum test frequency at temperature T2; The step of calculating the horizontal movement range of the dynamic modulus test results at adjacent temperatures according to the interpolation interval includes: Determine the frequency at which the dynamic modulus is D2 at temperature T1 by interpolation The frequency when the dynamic modulus is D1 at temperature T2 make in is the minimum test frequency at temperature T1, is the maximum test frequency at T2 temperature; when the data at T1 temperature is fixed, the horizontal movement range of the data at T2 temperature is: [-max(Δf1,Δf2),-min(Δf1,Δf2)]; Determining the shift value of the dynamic modulus test result at adjacent temperatures according to the horizontal movement range includes: Select n interpolation points at logarithmic equal intervals in the interpolation interval [D1, D2], divide the determined horizontal shift range into m parts at logarithmic equal intervals, horizontally shift the dynamic modulus data at temperature T2 according to the divided horizontal shift range, and calculate the norm Norm of the angular frequency difference of the n interpolation points of the two curves in the interval [D1, D2] after the shift, and find the horizontal shift value j* that minimizes Norm; where j* is the horizontal shift value of the data at temperature T2 when the data at temperature T1 is fixed; The calculation process of the horizontal shift value j* is as follows: In the formula, is the frequency of the i-th interpolation point in [D1, D2] at temperature T2 when the norm value is minimum, is the frequency of the i-th interpolation point in [D1, D2] when the horizontal shift value is the j-th value in [-max(Δf1, Δf2), -min(Δf1, Δf2)] at temperature T2.

2. The method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials according to claim 1, characterized in that: The step of calculating the shift factors at all test temperatures according to the shift values ​​includes: Calculate the shift values ​​of all adjacent temperatures, sum all the shift values ​​obtained between the remaining test temperatures and the reference temperature, and determine the sign of the value according to the difference between the remaining test temperatures and the reference temperature; when the test temperature is lower than the reference temperature, the sign is positive, and when the test temperature is higher than the reference temperature, the sign is negative.

3. The method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials according to claim 1, characterized in that: The method of drawing a temperature-shift factor curve in a semi-logarithmic coordinate according to the shift factor and determining the shift factor at any target reference temperature by interpolating the shift factor curve and moving it vertically includes: After obtaining the shift factor curve at the test temperature as the reference temperature, the shift factor at the target reference temperature is determined using the interpolation method, and the shift factor at the target reference temperature is subtracted from the shift factors at all test temperatures to obtain the final shift factor curve.

4. The method for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials according to claim 1, characterized in that: The following steps are also included: When the test data have no overlapping area, a linear equation or a polynomial is used to fit the test data respectively, and the overlapping area is obtained by extrapolation.

5. A device for constructing a master curve of dynamic viscoelastic behavior of asphalt-based materials, characterized in that: include: at least one processor; at least one memory for storing at least one program; When the at least one program is executed by the at least one processor, the at least one processor implements the method according to any one of claims 1 to 4.

6. A computer-readable storage medium storing a program executable by a processor, characterized in that: The processor-executable program is used to perform the method according to any one of claims 1 to 4 when executed by the processor.

Images