CN112765825A - Roadbed surface displacement response determination method considering transverse isotropy and modulus nonlinear distribution, application, equipment and storage medium - Google Patents

Roadbed surface displacement response determination method considering transverse isotropy and modulus nonlinear distribution, application, equipment and storage medium Download PDF

Info

Publication number
CN112765825A
CN112765825A CN202110108452.1A CN202110108452A CN112765825A CN 112765825 A CN112765825 A CN 112765825A CN 202110108452 A CN202110108452 A CN 202110108452A CN 112765825 A CN112765825 A CN 112765825A
Authority
CN
China
Prior art keywords
roadbed
modulus
displacement response
laplace
hankel
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110108452.1A
Other languages
Chinese (zh)
Other versions
CN112765825B (en
Inventor
张军辉
范海山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changsha University of Science and Technology
Original Assignee
Changsha University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changsha University of Science and Technology filed Critical Changsha University of Science and Technology
Priority to CN202110108452.1A priority Critical patent/CN112765825B/en
Publication of CN112765825A publication Critical patent/CN112765825A/en
Application granted granted Critical
Publication of CN112765825B publication Critical patent/CN112765825B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

The invention discloses a roadbed surface displacement response determining method considering transverse isotropy and modulus nonlinear distribution, and application, equipment and a storage medium thereof, wherein the method comprises the following steps: obtaining roadbed material parameters, applying dynamic load on the roadbed surface, and constructing a roadbed modulus function which is distributed in a vertical nonlinear manner; calculating a surface displacement response solution of the roadbed in a Laplace-Hankel domain, wherein the surface displacement response solution simultaneously considers the transverse isotropy and the nonlinear distribution of the roadbed modulus along the vertical direction; and calculating the surface displacement response of the current roadbed model in a time domain range based on a displacement response solution of the current roadbed mechanics model in a Laplace-Hankel domain. The invention fully reflects the non-uniform distribution of the modulus of the actual roadbed along the depth direction, simultaneously considers the transverse isotropy of the roadbed, improves the accuracy of the roadbed surface displacement response, has high precision and high speed, and can obtain the distribution condition of the modulus of the roadbed along the depth direction by combining with the roadbed nondestructive detection technology.

Description

Roadbed surface displacement response determination method considering transverse isotropy and modulus nonlinear distribution, application, equipment and storage medium
Technical Field
The invention belongs to the technical field of road engineering, and relates to a roadbed surface displacement response determination method considering transverse isotropy and modulus nonlinear distribution, application, equipment and a storage medium.
Background
With the development of nondestructive testing technology, some soil-based rapid nondestructive testing instruments are in succession, and the application of portable hammer type deflectometer (PFWD) is the most widespread. The PFWD carries out roadbed parameter inversion calculation by applying impact load on the top surface of the roadbed and combining a certain mechanical model and a parameter adjustment algorithm. However, the current mainstream roadbed parameter inversion calculation is still established on the basis of an isotropic linear elastic system, and the actual characteristics of the roadbed cannot be reflected. At present, a great deal of research results show that the roadbed has stronger stress dependence characteristics, and the modulus of the roadbed generally shows an ascending trend along with the depth along with the increase of the depth. In addition, the transverse isotropy of the subgrade due to natural settlement cannot be ignored in practical calculations.
In part of the prior art, for a pavement structure, the modulus of a material in each layer is considered to be uniform and does not change with space, but the modulus of a roadbed is changed in the vertical direction, the modulus has a nonlinear distribution characteristic in the vertical direction, and the displacement response of the roadbed surface is difficult to accurately determine by adopting the existing method. In part of the prior art, iterative computation is performed through structural layer strain and modulus nonlinearity on the basis of considering roadbed modulus homogeneity, the method is essentially a mathematical iterative method, the nonlinear distribution characteristic of the modulus is not considered on a mechanical model, the iterative method can only approach to a real result, the computation precision is low, repeated iteration is needed, and the computation is slow.
Disclosure of Invention
In order to solve the problems, the invention provides a roadbed surface displacement response determining method considering transverse isotropy and modulus nonlinear distribution and application thereof, which fully reflects the non-uniform distribution of the modulus of an actual roadbed along the depth direction, simultaneously considers the transverse isotropy of the roadbed, improves the accuracy of determining the roadbed surface displacement response, has the characteristics of high precision and high speed, can obtain the distribution condition of the modulus of the roadbed along the depth direction by combining with a roadbed nondestructive detection technology, and solves the problems in the prior art.
The invention adopts the technical scheme that a roadbed surface displacement response determining method considering transverse isotropy and modulus nonlinear distribution is specifically carried out according to the following steps:
step S1, obtaining roadbed material parameters, wherein the roadbed material parameters comprise: vertical modulus E of roadbed top surfacev0Modulus E of road bed infinite depthv∞Roadbed modulus non-uniform coefficient alpha and roadbed modulus ratio neHorizontal poisson's ratio muhVertical poisson's ratio muvSubgrade density rho and subgrade thickness H;
s2, applying dynamic load p (r, t) to the roadbed surface, establishing a cylindrical coordinate system by taking the dynamic load center as the coordinate center, wherein r represents a radial coordinate, z represents a vertical coordinate, phi represents an annular coordinate, and constructing a roadbed modulus function E which is distributed along the vertical non-linearityv(z)＝Ev0+(Ev∞-Ev0)(1-e-αz) (ii) a Calculating a surface displacement response solution of the roadbed in the Laplace-Hankel domain, which simultaneously considers the transverse isotropy and the nonlinear distribution of the roadbed modulus along the vertical direction
Figure BDA0002918403690000021
Step S3, based on the displacement response solution of the current roadbed mechanics model in the Laplace-Hankel domain
Figure BDA0002918403690000022
And calculating the surface displacement response of the current roadbed model in the time domain range.
The method is used for obtaining the vertical modulus distribution state of the on-site roadbed.
An electronic device, comprising:
a memory for storing instructions executable by the processor; and
a processor for executing the instructions to implement the above-described method.
A computer-readable storage medium for storing a computer program which, when executed, is capable of carrying out the above-mentioned method.
The invention has the beneficial effects that:
embodiment of the invention through the construction of EvThe (z) function describes the modulus distribution condition of the roadbed along the depth direction, and meanwhile, a physical equation with transverse isotropy is introduced to derive a variable coefficient partial differential equation set considering the characteristics of the roadbed. And a roadbed dynamic response analytic solution considering roadbed transverse isotropy and modulus along vertical nonlinear distribution and a corresponding determination method are provided through Hankel-Laplace integral transformation and a Frobenius method. The transverse isotropy of the roadbed and the nonlinear distribution characteristics of the modulus along the depth direction are fully considered during roadbed dynamic response calculation, and compared with a traditional elastic system, the method can truly reflect the actual state of the roadbed. Compared with traditional methods such as a finite element method and the like, the method has the characteristics of high precision, high speed and the like.
According to the determining method provided by the embodiment of the invention, the modulus distribution condition of the roadbed can be described, so that the rigidity of the roadbed can be evaluated more accurately and quickly compared with the traditional method; the method is combined with roadbed nondestructive detection, a roadbed vertical modulus curve can be obtained, and the roadbed humidification can be evaluated through the vertical modulus change curve monitored for a long time; and a more reliable theoretical basis is provided for roadbed design, detection, evaluation, research and the like.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram of a roadbed mechanics system with transverse isotropy and vertically nonlinear distribution of modulus under the action of a rigid bearing plate according to an embodiment of the invention.
FIG. 2 is a flow chart of a dynamic response determination method of an embodiment of the present invention.
FIG. 3 is a diagram of an ABAQUS axisymmetric model in an embodiment of the present invention.
Figures 4a-4b are graphs comparing a dynamic response determination method of an embodiment of the present invention to a method using a finite element model, ABAQUS.
Fig. 5 is a comparison graph of a numerical simulation calculation result obtained by a humidity-stress dependent finite element method and an actual real value of the modulus distribution in the depth direction of the roadbed structure obtained by fitting in the embodiment of the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the invention provides a roadbed dynamic response determination method with transverse isotropy and modulus distributed along vertical nonlinearity, which is specifically carried out according to the following steps as shown in fig. 2:
step S1, obtaining roadbed material parameters, wherein the roadbed material parameters comprise: vertical modulus E of roadbed top surfacev0Modulus E of road bed infinite depthv∞Roadbed modulus non-uniform coefficient alpha and roadbed modulus ratio neHorizontal poisson's ratio muhVertical poisson's ratio muvSubgrade density rho and subgrade thickness H; the roadbed modulus non-uniformity coefficient alpha represents the change speed of the roadbed along with the modulus in the depth direction; the subgrade modulus ratio neThe roadbed is uniformly distributed in the whole roadbed space; the thickness H may be finite or infinite.
Step S2, as shown in FIG. 1, applying dynamic load p (r, t) to the roadbed surface, establishing a cylindrical coordinate system with the dynamic load center as the coordinate center, wherein r represents a radial coordinate, z represents a vertical coordinate, and phi represents an annular coordinate, and establishing a roadbed modulus function E which is distributed along the vertical non-linear directionv(z)＝Ev0+(Ev∞-Ev0)(1-e-αz) In FIG. 1, R represents the radius of the rigid carrier plate; calculating a surface displacement response solution of the roadbed in the Laplace-Hankel domain, considering transverse isotropy and modulus along vertical nonlinear distribution
Figure BDA0002918403690000031
Roadbed surface displacement response solution of roadbed mechanics model considering transverse isotropy and modulus in vertical nonlinear distribution in Laplace-Hankel domain
Figure BDA0002918403690000032
Is calculated as:
Figure BDA0002918403690000041
in the formula (1), A1,A2,A3,A4For the boundary condition parameters, by solving the following matrix equation [ equation (2) ]]And (4) obtaining. Wherein, if the roadbed thickness is H → + ∞, there is A3＝A4If 0, then equation (2) may be degenerated into a two-dimensional matrix equation;
Figure BDA0002918403690000042
in the formula (2)
Figure BDA0002918403690000043
The result of Hankel-Laplace transformation for the dynamic load p (r, t),
Figure BDA0002918403690000044
xi and s are Hankel integral transformation and Laplace integral transformation coefficients respectively; ev0Representing the modulus of the top surface of the subgrade;
Figure BDA0002918403690000045
Figure BDA0002918403690000046
Figure BDA0002918403690000047
Figure BDA0002918403690000048
in the formula, alpha is a roadbed modulus non-uniform coefficient; e is constant, xi is modulus E through the roadbed top surfacev0Modulus E at infinite depth of roadbedv∞Obtaining by calculation according to the formula (7); xi represents a dimensionless coefficient of the relative magnitude of modulus of the top surface and the bottom of the roadbed, and the larger the numerical value is, the larger the difference of the upper modulus and the lower modulus of the roadbed is;
Figure BDA0002918403690000049
C11，C12，C13，C33，C44a parameter characterizing the transverse isotropy of the road base material in the problem of axial symmetry, C11＝kne(1-neμv)，
Figure BDA00029184036900000410
C13＝kneμv(1+μh)，
Figure BDA00029184036900000411
C44＝0.5(1+μv)-1；
Figure BDA00029184036900000412
k represents an intermediate variable used for calculating each parameter of transverse isotropy and simplifying an expression; n iseIs the subgrade modulus ratio; mu.sv，μhN is the serial number of the series, namely the Poisson ratio of the roadbed in the vertical direction and the horizontal direction.
On the basis of the roadbed mechanics model shown in fig. 1, partial differential equations shown in formulas (32) to (33) can be obtained by combining dynamic balance equations [ formulas (26) to (27) ] of the axial symmetry problem, physical equations representing axial symmetry transverse isotropy and geometric equations [ formulas (28) to (31) ].
Figure BDA0002918403690000051
Figure BDA0002918403690000052
Figure BDA0002918403690000053
Figure BDA0002918403690000054
Figure BDA0002918403690000055
Figure BDA0002918403690000056
Wherein u isr＝ur(r,z,t)，uz＝uz(r, z, t) are displacement components in the r-direction and the z-direction, respectively; s2, obtaining a displacement analysis solution in a Hankel-Laplace domain, and S3 obtaining a displacement response result in a time domain through integral inverse transformation; ρ is the subgrade density; t is time; sigmar＝σr(r,z,t),σz＝σz(r,z,t)，σφ＝σφ(r, z, t) are stress components in the r direction, the z direction and the phi direction of the cylindrical coordinate system, respectively; tau iszt＝τzt(r, z, t) is the shear stress in the z-r plane.
Figure BDA0002918403690000057
Figure BDA0002918403690000061
In order to solve the partial differential equation set shown in the equations (32) to (33), firstly, Laplace transformation about t is carried out on the equations (32) to (33); performing first-order Hankel transformation on r on the Laplace transformed expression (32); performing zeroth-order Hankel transformation on r for the Laplace transformed expression (33); obtaining a variable coefficient ordinary differential equation set of the formulas (34) to (35), wherein the independent variable z is omitted. Wherein the Laplace transform for x is defined:
Figure BDA0002918403690000062
define the zeroth order Hankel transform for x:
Figure BDA0002918403690000063
define a first order Hankel transform for x:
Figure BDA0002918403690000064
ξ and s represent the Hankel transform and Laplace integral transform coefficients, respectively.
Figure BDA0002918403690000065
Figure BDA0002918403690000066
Figure BDA0002918403690000067
Represents a pair urThe result after the first-order transformation of Hankel-Laplace is carried out,
Figure BDA0002918403690000068
represents a pair uzTo carry out HankAnd (5) performing the result of the el-Laplace zero-order transformation.
To solve the system of constant differential equations with variable coefficients as shown in equations (34) to (35), the following variables are introduced:
ψ＝Ξe-αz (36)
xi calculated according to equation (7), and in subsequent derivations, E is obtained by replacing the independent variable z in equation by the variable ψvThe expression of (ψ) is shown in formula (38). Variable substitution of z by Ψ, Ev(psi) and Ev(z) has the same meaning except that the argument is replaced and its value is unchanged. The independent variable is z in (34) to (35) and ψ in (40) to (43); the purpose is to eliminate Ev(z) an exponential term present, which is converted to algebraic form as shown in (39); only the variable substitutions are made and the expressions themselves are not changed.
Meanwhile, the relation of the formula (39) can be obtained by applying a derivative theory; where f (z) denotes any function with z as argument:
Figure BDA0002918403690000071
Figure BDA0002918403690000072
by bringing expressions (38) to (39) into expressions (34) to (35), a system of constant differential equations with a variable coefficient having a ψ argument can be obtained as shown in expressions (40) to (43).
Wherein formulae (40) to (41) correspond to Ev0<Ev∞In the case of (1), expressions (42) to (43) correspond to Ev0>Ev∞The case (1). The argument ψ is omitted in the formula; where δ is an intermediate calculation parameter, δ ═ Ev0/Ev∞。
Figure BDA0002918403690000073
Figure BDA0002918403690000074
Figure BDA0002918403690000075
Figure BDA0002918403690000076
Based on an infinite series theory, approximation is carried out by using infinite series, and if the ordinary differential equations shown in the formulas (40) to (43) exist, displacement solutions shown in the formulas (44) to (45) exist; wherein
Figure BDA0002918403690000077
To represent
Figure BDA0002918403690000078
A first derivative of (a) is obtained,
Figure BDA0002918403690000079
to represent
Figure BDA00029184036900000710
A second derivative of;
Figure BDA00029184036900000711
to represent
Figure BDA00029184036900000712
A first derivative of (a) is obtained,
Figure BDA00029184036900000713
to represent
Figure BDA00029184036900000714
A second derivative of;
Figure BDA00029184036900000715
Figure BDA00029184036900000716
in the formula, m is a parameter, and n is a serial number of a series sequence;
Figure BDA0002918403690000081
corresponding m in infinite series representing z-direction displacement componentiThe coefficient of the nth order of terms of (a),
Figure BDA0002918403690000082
corresponding m in infinite series expressing r-direction displacement componentiThe coefficient of the nth order of terms of (a); a isnIs composed of
Figure BDA0002918403690000083
Abbreviation of (a), (b)nIs composed of
Figure BDA0002918403690000084
The abbreviation of (1); the series form shown in the formulas (44) to (45) is subjected to displacement solution
Figure BDA0002918403690000085
Figure BDA0002918403690000086
By substituting the equations (40) to (43), the equations shown in the equations (46) to (49) can be obtained. The formulae (46) to (47) correspond to Ev0<Ev∞In the case of (1), expressions (48) to (49) correspond to Ev0>Ev∞The case (1). Where θ is an intermediate calculation parameter, and θ is ρ s2/Ev∞；
Figure BDA0002918403690000087
Figure BDA0002918403690000088
Figure BDA0002918403690000089
Figure BDA00029184036900000810
Let n be 0, and compare equations (46) to (49) on both sides, the following relationship can be obtained:
[α2C44m2-C11ξ2-θ]a0+αξ(C13+C44)mb0＝0 (50)
αξ(C13+C44)ma0-[α2C33m2-C44ξ2-θ]b0＝0 (51)
if equations (50) to (51) are satisfied, m must satisfy equation (52).
Figure BDA00029184036900000811
In the complex domain, there must be four solutions to the equation for m [ equation (52) ] followed by:
Figure BDA00029184036900000812
Figure BDA00029184036900000813
Figure BDA0002918403690000091
Figure BDA0002918403690000092
in formulae (8) to (11), miFor intermediate calculation of the parameters, theta1，Θ2，Θ3Respectively calculating parameters for the intermediate；Θ1＝α2C33C44；
Figure BDA0002918403690000093
Θ3＝(C44ξ2+θ)(C11ξ2+ θ); where θ is an intermediate calculation parameter, θ ═ ρ s2/Ev∞(ii) a Rho is roadbed density, s in the theta expression is Hankel integral transformation and Laplace integral transformation coefficients, and after integral transformation, the coefficients become s equivalent to time t and can be regarded as a variable.
m has four solutions, and it can be seen from equation (44) that the final solution should be in the form of a combination of four infinite series numbers, in order to distinguish the coefficients a in the four infinite series numbersnTherefore, in the equation (12), the superscripts (1), (2), (3), and (4) correspond to 4 solutions of m, respectively. By bringing formula (12) into formula (50), the compound represented by formula (13) can be obtained
Figure BDA0002918403690000094
And (5) expressing.
Figure BDA0002918403690000095
Respectively representing the corresponding m in infinite seriesiThe coefficient of the nth order of the terms of (1) is a specific coefficient sequence;
Figure BDA0002918403690000096
Figure BDA0002918403690000097
the coefficients of the corresponding stages in (46) to (49) are equal, that is, the stages of the same power should be equal on both sides of the equation, thereby obtaining
Figure BDA0002918403690000098
See the expressions (14) and (15).
Figure BDA0002918403690000099
Figure BDA00029184036900000910
Calculation according to the recurrence relation shown in equations (12) to (15)
Figure BDA00029184036900000911
The expanded expressions of the parameters in the formulas (14) and (15) are shown in (16) to (21):
Figure BDA00029184036900000912
Figure BDA00029184036900000913
Figure BDA00029184036900000914
Figure BDA00029184036900000915
Figure BDA0002918403690000101
Figure BDA0002918403690000102
in the formulae (12) to (15),
Figure BDA0002918403690000103
for the intermediate calculation parameters, the calculation is performed in accordance with equations (16) to (21). Wherein the content of the first and second substances,
Figure BDA0002918403690000104
is calculated with Ev0And Ev∞Are different, and are considered in two cases, when Ev0<Ev∞When calculated according to equation (20), when Ev0>Ev∞Then, the calculation is performed according to equation (21), where δ is an intermediate calculation parameter, and δ is equal to E in equation (21)v0/Ev∞。
In formula (20)
Figure BDA0002918403690000105
Introducing intermediate parameters
Figure BDA0002918403690000106
The recursion formula can be obtained by the corresponding number of steps of (46) - (49) being equal, because the result that n is 0 is known, and the result that n is 1 can be obtained by this, and the recursion calculation is carried out in turn, and n is 2, n is 3, and … can be calculated.
According to the ordinary differential solution theory, the equation is solved into a series form, m has corresponding four results, so that the equation has four infinite series solutions in total, and the final general solution formula can be represented by linear combination of the four infinite series solutions, namely the final displacement general solution formula (53) and the final displacement general solution formula (1):
Figure BDA0002918403690000107
Figure BDA0002918403690000108
four undetermined coefficients A exist in the general solution1，A2，A3，A4. To solve for its value, certain boundary conditions need to be brought in. According to fig. 1, boundary conditions as shown in formulas (54) to (57) are introduced; wherein p (r, t) is a roadbed top surface load expression.
Figure BDA0002918403690000109
Figure BDA00029184036900001010
ur|z＝H＝0 (56)
uz|z＝H＝0 (57)
Laplace transform on t and Hankel transform on 0 th order of r are executed on the equations (54) and (57); for equation (55), equation (56) performs Laplace transform on t and Hankel transform of order 1 on r. Then, by introducing the conversion relationships shown in the expressions (7) and (38), the boundary conditions shown in the expressions (58) to (61) can be obtained.
Figure BDA0002918403690000111
Figure BDA0002918403690000112
Figure BDA0002918403690000113
Figure BDA0002918403690000114
The general solution expressions of the formula (53) and the formula (1) are brought into the boundary conditions shown in the formulas (58) to (61), so that a fourth-order matrix equation shown in the formula (2) can be obtained, and the matrix equation shown in the formula (2) is solved, so that the coefficient A to be determined can be obtained1，A2，A3，A4(ii) a In particular, in the case of roadbed thickness H → + ∞, if the expressions (60) to (61) need to be satisfied, a is always present3＝A4The matrix equation shown in equation (2) can be further degenerated into a second-order matrix equation. Determined A1，A2，A3，A4The analytical solution results in the Hankel-Laplace domain can be obtained by substituting the formula (53) and the formula (1). Further, a Hankel inverse transform and a Laplace inverse transform are performed on the result, and a time domain solution of the time domain can be obtained. The calculation process is shown in fig. 2.
In practical calculation, the vertical modulus E of the top surface of the roadbed is firstly calculatedv0Modulus E of road bed infinite depthv∞Roadbed modulus non-uniform coefficient alpha and roadbed modulus ratio neHorizontal poisson's ratio muhVertical poisson's ratio muvM is calculated according to the equations (8) to (11) based on the parameters such as the roadbed density rho and the roadbed thickness H1，m2，m3，m4(ii) a Further, coefficient series were calculated from equations (12) to (21)
Figure BDA0002918403690000115
Further establishing a fourth-order matrix equation according to the formulas (2) to (6), and solving the coefficient A on the basis of the fourth-order matrix equation1,A2,A3,A4. Finally, calculating a roadbed surface dynamic response solution in a Hankel-Laplace domain according to the formula (1)
Figure BDA0002918403690000116
Step S3, based on the displacement response solution of the current roadbed mechanics model in the Laplace-Hankel domain
Figure BDA0002918403690000117
And calculating the surface displacement response of the current roadbed model in the time domain range.
Displacement response solution of current roadbed mechanics model in Laplace-Hankel domain through formula (23)
Figure BDA0002918403690000118
And carrying out Hankel inverse transformation, and then carrying out Laplace inverse transformation on the result of the Hankel inverse transformation through a formula (22) to obtain the surface displacement response of the current roadbed mechanics system in a time domain range:
Figure BDA0002918403690000121
Figure BDA0002918403690000122
wherein the content of the first and second substances,
Figure BDA0002918403690000123
is composed of
Figure BDA0002918403690000124
F (r,0, t)l) Is uz(r,0,tl)；uz(r,0,tl) For the surface displacement response of the current roadbed mechanical model in the time domain range, the corresponding time step is tl，tl＝l·T/NaT is the total length of solution, NaIs the total time increment step number of the solution, l is an increment step sequence, and l is less than or equal to N; r is the horizontal distance from the calculated point position to the load center; Δ t is the time increment step; the formula (22) shows that inverse Lpalace transformation is carried out by Durbin method, L and gamma are calculation parameters of Durbin method, and L.Na50-5000, and 5-10 of gamma and T; j is a unit imaginary number, j2-1; for the Durbin method, the final time solution of 1/4 can generate a large system error, and the solution is usually solved by prolonging the total time length T of solution, through testing, the invention can meet the precision requirement by taking 2 times of time length (when inverse transformation is performed on the solution in the Hankel-Laplace domain, the time length of calculation needs to be given, for example, I needs to calculate 60ms dynamic response, because the Durbin method has certain defects, the time length parameter of calculation is generally 2-3 times in actual calculation, the embodiment of the invention takes 2 times of time length, namely the time length parameter T of calculation takes 120ms, but the final result only takes 60 ms). χ is the upper limit of the numerical integration, AkThe weight of the kth integral node of the 20-node Gaussian interpolation; x is the number ofikIs the kth integral node value, x, of the i integral sections corresponding to the Gaussian integralik＝(i-1)ΔL+xkΔ L is the Gaussian integration segment division length, xkA kth integrator node value that is a 20 node gaussian integral; n is a radical of0Is the total number of segments of the Gaussian integral, N0＝ceil(χ/ΔL) obtained by calculating (χ/Δ L) and rounding up;
Figure BDA0002918403690000125
representing the result of the arbitrary function f (-) after the Hankel-Laplace transform;
Figure BDA0002918403690000126
representing the result of the arbitrary function f (-) after Laplace transformation; re (-) is the operation of taking a complex real part; j. the design is a square0(xikr) represents an independent variable of (x)ikr) function result corresponding to the zero order Bessel function. The Laplace inverse transformation is carried out by Durbin method (first-order autocorrelation inspection method) in the formula (1) and the formula (22), the Hankel inverse transformation is carried out by Gaussian interpolation integration in the formula (23), the calculation processes are all realized by adopting a programming mode, and only a numerical value determination method is introduced here. The numerical value determination methods are all realized by programming, programming languages are not limited, and data calculated in the embodiment of the invention are all calculated by MATLAB 2016 a.
After verification, when the upper limit n of the infinite series sequence in the formula (1) is 100, the upper limit χ of the gaussian integral in the formula (23) and the length Δ L of the integral interval are calculated according to the formulas (24) to (25), the result can be converged to an accurate solution.
Figure BDA0002918403690000127
Figure BDA0002918403690000131
Effect verification:
the comparative calculation was performed using finite element software ABAQUS. The subgrade parameters are shown in table 1, and the defined subprogram UMAT is adopted in ABAQUS to realize the nonlinear distribution of the subgrade along the depth direction. p (r, t) is the load corresponding to the rigid bearing plate with the diameter of 30cm, the load is in a half sine wave form, the peak value is 100kPa, and the action time is 20 ms; and setting the calculation time length to be 40ms and the time step increment to be 1ms during calculation. When the invention is applied to calculation, according to the parameter recommended values of the equations (24) to (25), the calculation parameters are set as: n is 100, χ is 300, Δ L is 10.
Table 1 roadbed structure calculation parameter table
Name (R) Cell type modulus/MPa ne μv(＝μh) ρ/kg·m-3 H/m
Bearing plate CAX4 2.1×105 1.0 0.25 7500 -
Road bed CAX4/CINAX4 100+200(1-e-2z) 1.0 0.35 1816 1.0/2.0
As shown in fig. 3, the entire ABAQUS finite element model is modeled in an axisymmetric manner, and in order to simulate the horizontal direction infinity of the layered system, infinite CINAX4 cells are used for the boundary portion, and CAX4 cells are used for the finite size portion. The calculation result is shown in fig. 4, MATLAB in fig. 4 is a calculation result obtained by using the method of the embodiment of the present invention, ABAQUS is a calculation result using a finite element model, and the ABAQUS finite element can be used as a method for widely checking the calculation correctness. In the calculation, ABAQUS obtained by adopting an Intel i 78750H/8G notebook computer to test needs 5 minutes for completing the calculation, but the method only needs 2 minutes, thereby greatly reducing the calculation cost.
When the non-uniform distribution of the roadbed is calculated according to the prior method, local meteorological data need to be acquired, and the meteorological data mainly comprise: and calculating the non-uniform humidity field distribution of the roadbed by adopting calculation software according to the data of wind speed, temperature, precipitation and the like. On the basis, the non-uniform modulus field of the roadbed is calculated by considering the roadbed stress state and the non-uniform humidity field data and combining the roadbed soil dynamic resilience modulus and humidity coupling equation. The calculation consumption is large, and the method is not suitable for popularization.
The invention employs Ev(z)＝Ev0+(Ev∞-Ev0)(1-e-αz) Approximately describing the vertical modulus distribution of the roadbed, wherein E in the formulav0，Ev∞And α is a fitting parameter. For verifying the rationality of the formula, the roadbed modulus along the depth direction of the center of the driveway in the initial state, the operation 1 year, the operation 2 years and the balance humidity state is calculated by adopting meteorological data of nearly 16 years in Nanchang city and combining a soil-based modulus-humidity coupling equation, and E is adoptedv(z) fitting in conjunction with least squaresAs can be seen from FIG. 5, the fitting effect is good, indicating EvAnd (z) the modulus distribution rule of the roadbed along the depth direction in the operation process can be represented more ideally.
E can be applied no matter how long the roadbed is operatedvThe formula (z) fits the test sample, the fitting effect is good, and the result shows thatvAnd (z) the modulus distribution rule of the roadbed along the depth direction in the operation process can be represented more ideally.
In pair EvOn the basis of (z) verification of rationality, the embodiment of the invention further provides a roadbed surface displacement response determining method based on the roadbed nonlinear distribution, and the determining method can be used for back-calculating the modulus distribution of the roadbed, and specifically comprises the following steps:
actually measuring displacement response on site, and recording the displacement response and load data on site; randomly assuming a series of E within a certain rangev(z) parameter Ev0、Ev∞Alpha and subgrade modulus ratio neWherein E isv0In the range of 20MPa to 100MPa, Ev∞The range of (a) is 100MPa to 300MPa, the range of the roadbed modulus non-uniform coefficient alpha is 0.05 to 5.0, neThe range of (1) to (4); while the remaining parameters generally take the following values: the poisson ratio is 0.35, and the density is 1.8-2.0 g/cm3In a hypothetical series EvOn the basis of the (z) parameter, the obtained surface load data is taken as p (r, t) in S2, and a series of theoretical displacement response calculation results can be obtained by combining and adopting the calculation method provided by the invention and are compared and matched with an actual displacement curve. According to the matching result, the parameter E is adjusted through a parameter adjusting algorithm (such as a genetic algorithm, a group intelligent algorithm and the like)v0、Ev∞、α、neAnd adjusting to enable the theoretical displacement response to be close to the field actual measurement displacement response to the maximum extent. And taking the parameter value corresponding to the most matched theoretical displacement response calculation result as an inversion result of the actual roadbed, thereby obtaining the modulus distribution rule of the roadbed along the depth direction and realizing roadbed rigidity evaluation.
In addition, a large number of theories can be obtained in advance through the calculation method provided by the inventionAnd (3) constructing a database by the aid of the theoretical displacement response, and quickly matching the field actual measurement displacement curve by means of combining the conventional BP neural network and a database search method, so that the vertical modulus distribution state of the field roadbed can be quickly obtained. Obtaining E corresponding to the theoretical displacement response with the closest actual displacement responsev(z) parameter and neAnd thus obtaining a modulus distribution function in the depth direction of the actual roadbed structure. Compared with the traditional method which can only obtain a single structural modulus, the method can more accurately reflect the vertical modulus distribution state of the roadbed by calculation.
Obtaining the result of the distribution of the vertical modulus of the roadbed by iterative calculation of climate data and a series of finite elements, and then using the method E of the inventionvAnd (z) fitting to obtain theoretically true roadbed vertical modulus distribution, and assuming a series of loads p (r, t) to obtain a series of theoretical displacement response calculation results for theoretical analysis and calculation. Such as the impact of overload, calculating subgrade operating areas, etc.
The method is combined with roadbed nondestructive detection, a roadbed vertical modulus curve can be obtained, and roadbed humidification can be evaluated by monitoring the vertical modulus change curve for a long time. The method for evaluating the wetting of the roadbed comprises the following steps: wetting the roadbed to reduce the vertical modulus curve of the roadbed as a whole; the humidity state of the roadbed can be qualitatively evaluated by regularly detecting the roadbed and comparing with roadbed vertical modulus distribution curves in different periods; if the modulus of the top of the roadbed is reduced, the humidity of the roadbed is integrally increased, and if the modulus of the bottom is basically unchanged, the bottom is proved to reach the equilibrium humidity, and the upper part does not reach the equilibrium humidity; the method is that the roadbed nondestructive test using a homogeneous elastomer system as a mechanical model can not be realized at present, and the traditional mechanical model can only obtain one modulus which approximately represents the modulus of the whole roadbed, so that the depth distribution rule can not be obtained.
The advantages of the invention are as follows:
1. at present, the nonuniform modulus distribution of the roadbed is acknowledged by scholars in the field, but in practical application, when the traditional method is used for solving and calculating the dynamic response considering the nonuniform modulus distribution of the roadbed along the vertical direction, a finite element method is mostly needed, although the method can obtain a more accurate result, the calculation speed is slower, the modeling is required to be carried out again every time, and particularly, for a certain specific situation, the method needs to be independently modeled, so that the time consumption is longer, the cost is high, and the calculation is inconvenient to carry out on site. The invention introduces the nonlinear distribution state of the modulus along the depth direction, combines the prior ordinary differential and partial differential solution theories to carry out the relevant derivation of the analytical solution, can be compiled into any program code to be embedded into small equipment in the derivation and calculation process, can be realized by programming, is simple and convenient to calculate, avoids the complicated process of finite elements, can realize the rapid calculation on site, and meets the relevant requirements of actual calculation.
2. Prior art 1 (multilayer displacement response determination method considering interlayer contact conditions and transverse isotropy) is to divide a pavement structure into N layers, but the modulus in each layer is uniform. The invention only considers the roadbed, but the modulus distribution of the roadbed is uneven, and the fundamental differences of the modulus in the mechanical model are larger. The invention provides a method for determining pavement displacement of a roadbed with isotropy in transverse view and nonlinear distribution of modulus under the impact load of nondestructive testing equipment such as PFWD (frequency-modulated wavelength distribution), and the like.
The roadbed surface displacement response determining method considering transverse isotropy and modulus nonlinear distribution can be stored in a computer readable storage medium if the roadbed surface displacement response determining method is realized in the form of a software functional module and is sold or used as an independent product. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the roadbed displacement response determining method considering transverse isotropy and modulus nonlinear distribution according to the embodiment of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a U disk, a removable hard disk, a ROM, a RAM, a magnetic disk, or an optical disk.
The above description is only for the preferred embodiment of the present invention, and is not intended to limit the scope of the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention shall fall within the protection scope of the present invention.

Claims (10)

1. A roadbed surface displacement response determining method considering transverse isotropy and modulus nonlinear distribution is characterized by comprising the following steps:
step S1, obtaining roadbed material parameters, wherein the roadbed material parameters comprise: vertical modulus E of roadbed top surfacev0Modulus E of road bed infinite depthv∞Roadbed modulus non-uniform coefficient alpha and roadbed modulus ratio neHorizontal poisson's ratio muhVertical poisson's ratio muvSubgrade density rho and subgrade thickness H;
s2, applying dynamic load p (r, t) to the roadbed surface, establishing a cylindrical coordinate system by taking the dynamic load center as the coordinate center, wherein r represents a radial coordinate, z represents a vertical coordinate, phi represents an annular coordinate, and constructing a roadbed modulus function E which is distributed along the vertical non-linearityv(z)＝Ev0+(Ev∞-Ev0)(1-e-αz) (ii) a Calculating a surface displacement response solution of the roadbed in the Laplace-Hankel domain, which simultaneously considers the transverse isotropy and the nonlinear distribution of the roadbed modulus along the vertical direction
Figure FDA0002918403680000011
Step S3, based on the displacement response solution of the current roadbed mechanics model in the Laplace-Hankel domain
Figure FDA0002918403680000012
And calculating the surface displacement response of the current roadbed model in the time domain range.
2. The roadbed surface displacement response determination method considering transverse isotropy and modulus nonlinear distribution as claimed in claim 1, wherein in the step S2,
Figure FDA0002918403680000013
Figure FDA0002918403680000014
in the formula (1-1), A1,A2,A3,A4In order to be a parameter of the boundary condition,
Figure FDA0002918403680000019
corresponding m in infinite series expressing r-direction displacement componentiI is 1,2,3, 4; n is the number of the series of stages, m1～m4Is an intermediate parameter; xi denotes a dimensionless coefficient of the relative magnitude of the modulus of the top and bottom subgrade surfaces;
Figure FDA0002918403680000015
Figure FDA0002918403680000016
Figure FDA0002918403680000017
Figure FDA0002918403680000018
in the formulae (1-3) to (1-6), Θ1，Θ2，Θ3Respectively intermediate calculation parameters; theta1＝α2C33C44；
Figure FDA0002918403680000021
Θ3＝(C44ξ2+θ)(C11ξ2+ θ); where θ is an intermediate calculation parameter, θ ═ ρ s2/Ev∞(ii) a Xi and s are Hankel integral transformation coefficients and Laplace integral transformation coefficients respectively, and time t becomes s after integral transformation; c11，C12，C13，C33，C44Parameters for characterizing the transverse isotropy characteristics of the road base material in the axial symmetry problem; c11＝kne(1-neμv)，
Figure FDA0002918403680000022
C13＝kneμv(1+μh)，
Figure FDA0002918403680000023
C44＝0.5(1+μv)-1；
Figure FDA0002918403680000024
k represents an intermediate variable.
3. The method for determining the displacement response of the roadbed surface considering the transverse isotropy and the non-linear distribution of the modulus as claimed in claim 2, wherein in the step S2, the boundary condition parameter A1,A2,A3,A4The method is specifically determined according to the following method:
Figure FDA0002918403680000025
Figure FDA0002918403680000026
Figure FDA0002918403680000027
Figure FDA0002918403680000028
Figure FDA0002918403680000029
Figure FDA00029184036800000210
the result of the Hankel-Laplace transform for the load p (r, t), Ψ1i、Ψ2i、Ψ3i、Ψ4iIn order to calculate the parameters in the middle of the run,
Figure FDA00029184036800000211
corresponding m in infinite series representing z-direction displacement componentiThe coefficient of the nth order of terms of (1).
4. The roadbed surface displacement response determination method considering transverse isotropy and modulus nonlinear distribution as claimed in claim 3, wherein in the step S2,
Figure FDA00029184036800000212
the method is specifically determined according to the following method:
Figure FDA00029184036800000213
Figure FDA00029184036800000214
Figure FDA0002918403680000031
Figure FDA0002918403680000032
in the formulae (1-12) to (1-15),
Figure FDA0002918403680000033
calculating parameters for intermediate calculation according to the formulas (1-16) to (1-21); when E isv0<Ev∞When E is calculated according to the formula (1-20)v0>Ev∞Then, the calculation is performed according to the formula (1-21), where δ is an intermediate calculation parameter, and δ is equal to Ev0/Ev∞；
Figure FDA0002918403680000034
Figure FDA0002918403680000035
Figure FDA0002918403680000036
Figure FDA0002918403680000037
Figure FDA0002918403680000038
Figure FDA0002918403680000039
5. The method for determining the displacement response of the roadbed surface by considering the transverse isotropy and the nonlinear distribution of the modulus as claimed in claim 2, wherein the intermediate parameter m1～m4The determination method comprises the following steps:
on the basis of a roadbed mechanics model, a dynamic balance equation of an axisymmetric problem, a physical equation representing axisymmetric transverse isotropy and a geometric equation are combined to obtain a partial differential equation set shown in formulas (1-22) to (1-23):
Figure FDA00029184036800000310
Figure FDA00029184036800000311
wherein u isr＝ur(r,z,t)，uz＝uz(r, z, t) are displacement components along the radius r direction of the cylindrical coordinate system and the roadbed depth z direction respectively; t is time; sigmar＝σr(r,z,t),σz＝σz(r,z,t),σφ＝σφ(r, z, t) are stress components in the r direction, the z direction and the phi direction, respectively; tau iszt＝τzt(r, z, t) is the shear stress in the z-r plane;
firstly, Laplace conversion of t is carried out on the formulas (1-22) to (1-23); performing first-order Hankel transformation on r on the Laplace transformed expression (1-22); performing zeroth-order Hankel transformation on r on the Laplace transformed (1-23); obtaining a variable coefficient ordinary differential equation set of formulas (1-24) to (1-25), wherein an independent variable z is omitted;
Figure FDA0002918403680000041
Figure FDA0002918403680000042
wherein the content of the first and second substances,
Figure FDA0002918403680000043
represents a pair urThe result after the first-order transformation of Hankel-Laplace is carried out,
Figure FDA0002918403680000044
represents a pair uzPerforming Hankel-Laplace zeroth order transformation;
in order to solve the variable coefficient ordinary differential equation set shown in equations (1-24) to (1-25), a variable ψ is introduced:
ψ＝Ξe-αz (1-26)
replacement of the subgrade modulus function E with the variable psiv(z)＝Ev0+(Ev∞-Ev0)(1-e-αz) Medium independent variable z, to EvThe expression of (ψ) is as shown in the formula (1-27):
Figure FDA0002918403680000045
meanwhile, through derivative theory, the formula (1-28) can be obtained; where f (z) denotes any function with z as argument:
Figure FDA0002918403680000046
the equations (1-27) to (1-28) are taken into the equations (1-24) to (1-25) to obtain a system of constant differential equations with a variable ψ as shown in the equations (1-29) to (1-32); wherein formulae (1-29) to (1-30) correspond to Ev0<Ev∞In the case where the formulae (1-31) to (1-32) correspond to Ev0>Ev∞The case (1); the argument ψ is omitted in the formula; where δ is an intermediate calculation parameter, δ ═ Ev0/Ev∞；
Figure FDA0002918403680000051
Figure FDA0002918403680000052
Figure FDA0002918403680000053
Figure FDA0002918403680000054
In order to solve the ordinary differential equation sets shown in the formulas (1-29) to (1-32), approximation is carried out by using infinite series based on infinite series theory, and the ordinary differential equation sets shown in the formulas (1-29) to (1-32) are assumed to have displacement solutions shown in the formulas (1-33) to (1-34);
Figure FDA0002918403680000055
Figure FDA0002918403680000056
in the formula, m is a parameter, and n is a serial number of a series sequence; the series form shown in the formulas (1-33) to (1-34) is subjected to displacement solution
Figure FDA0002918403680000057
The equations are introduced into the equations (1-29) to (1-32) to obtain the equation relations shown in the equations (1-35) to (1-38); formulae (1-35) to (1-36) correspond to Ev0<Ev∞In the case where the formulae (1-37) to (1-38) correspond to Ev0>Ev∞The case (1);
Figure FDA0002918403680000058
Figure FDA0002918403680000059
Figure FDA0002918403680000061
Figure FDA0002918403680000062
when n is 0, the following relationships can be obtained by comparing equations (1-35) to (1-38):
[α2C44m2-C11ξ2-θ]a0+αξ(C13+C44)mb0＝0 (1-39)
αξ(C13+C44)ma0-[α2C33m2-C44ξ2-θ]b0＝0 (1-40)
if the expressions (1-39) to (1-40) are satisfied, m must satisfy the expression (1-41);
Figure FDA0002918403680000063
in the complex domain, there must be four solutions for equation (1-41) for m, see equations (1-3) through (1-6).
6. The method for determining the displacement response of the roadbed surface considering the transverse isotropy and the nonlinear distribution of the modulus as claimed in claim 3, wherein the boundary conditions are as shown in the formulas (1-42) to (1-45):
Figure FDA0002918403680000064
Figure FDA0002918403680000065
ur|z＝H＝0 (1-44)
uz|z＝H＝0 (1-45)
performing Laplace transformation on t and Hankel transformation of 0 order on r on the equations (1-42) and (1-45); performing Laplace transformation on t and 1-order Hankel transformation on r for equations (1-43) and equations (1-44); on the basis, the transformation relations shown in the formulas (1-2) and (1-27) are introduced, so that the boundary conditions shown in the formulas (1-46) to (1-49) can be obtained;
Figure FDA0002918403680000066
Figure FDA0002918403680000067
Figure FDA0002918403680000071
Figure FDA0002918403680000072
Figure FDA0002918403680000073
is composed of
Figure FDA0002918403680000074
The first derivative of (a) is,
Figure FDA0002918403680000075
is composed of
Figure FDA0002918403680000076
The first derivative of (1) can be obtained by bringing the general solution expression into the boundary conditions shown in (1-46) to (1-49) to obtain a fourth-order matrix equation shown in (1-7), and the coefficient A to be determined can be obtained by solving the matrix equation shown in (1-7)1，A2，A3，A4。
7. The method for determining a subgrade surface displacement response considering transverse isotropy and modulus nonlinear distribution according to claim 1, wherein said step S3 is specifically: displacement response solution of current roadbed mechanics model in Laplace-Hankel domain through formula (1-51)
Figure FDA0002918403680000077
Namely, the formula (1-1) is used for carrying out Hankel inverse transformation, then the Laplace inverse transformation is carried out on the result of the Hankel inverse transformation after the formula (1-50) is used for obtaining the surface displacement response u of the current mechanical model in the time domainz(r,0,tl)；
Figure FDA0002918403680000078
Figure FDA0002918403680000079
Wherein the content of the first and second substances,
Figure FDA00029184036800000710
to represent
Figure FDA00029184036800000711
f(r,0,tl) Represents uz(r,0,tl) (ii) a Time step of tl，tl＝l·T/NaT is the total length of solution, NaIs the total time increment step number of the solution, l is an increment step sequence, and l is less than or equal to N; r is the horizontal distance from the calculated point position to the load center; Δ t is the time increment step; the expression (1-50) shows that inverse Lpalace transform is performed by Durbin method, L and gamma are Durbin method calculation parameters, j is unit imaginary number2-1; χ is the upper limit of the numerical integration, AkThe weight of the kth integral node of the 20-node Gaussian interpolation; x is the number ofikIs the kth integral node value, x, of the i integral sections corresponding to the Gaussian integralik＝(i-1)ΔL+xkΔ L is the Gaussian integration segment division length, xkA kth integrator node value that is a 20 node gaussian integral; n is a radical of0Is the total number of segments of the Gaussian integral, N0Ceil (χ/Δ L), obtained by pair calculation (χ/Δ L) and rounding up; re (-) is the operation of taking a complex real part; j. the design is a square0(xikr) represents an independent variable of (x)ikr) function result corresponding to the zero order Bessel function.
8. A method according to any one of claims 1 to 7 for obtaining a vertical modulus profile of a subgrade in a site.
9. An electronic device, comprising:
a memory for storing instructions executable by the processor; and
a processor for executing the instructions to implement the method of any one of claims 1-7.
10. A computer-readable storage medium for storing a computer program which, when executed, is capable of carrying out the method of any one of claims 1-7.
CN202110108452.1A 2021-01-27 2021-01-27 Roadbed surface displacement response determination method considering modulus nonlinear distribution Active CN112765825B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110108452.1A CN112765825B (en) 2021-01-27 2021-01-27 Roadbed surface displacement response determination method considering modulus nonlinear distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110108452.1A CN112765825B (en) 2021-01-27 2021-01-27 Roadbed surface displacement response determination method considering modulus nonlinear distribution

Publications (2)

Publication Number Publication Date
CN112765825A true CN112765825A (en) 2021-05-07
CN112765825B CN112765825B (en) 2022-07-08

Family

ID=75705990

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110108452.1A Active CN112765825B (en) 2021-01-27 2021-01-27 Roadbed surface displacement response determination method considering modulus nonlinear distribution

Country Status (1)

Country Link
CN (1) CN112765825B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070074565A1 (en) * 2005-09-30 2007-04-05 Paramsothy Jayakumar System for virtual prediction of road loads
US20090070042A1 (en) * 2007-09-11 2009-03-12 Richard Birchwood Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state
CN108517735A (en) * 2018-04-12 2018-09-11 长沙理工大学 A kind of durable asphalt pavement design method and its pavement structure based on bimodular theory
CN110502845A (en) * 2019-08-27 2019-11-26 长沙理工大学 Method is determined based on the PFWD subgrade modulus of genetic algorithm and viscoelasticity theory
CN110658086A (en) * 2019-08-31 2020-01-07 长沙理工大学 Asphalt pavement load response analysis method considering tension-compression modulus difference
CN111475876A (en) * 2020-03-04 2020-07-31 长沙理工大学 Method for obtaining dynamic resilience mechanical characteristic parameters of granules
CN112214817A (en) * 2020-09-29 2021-01-12 长沙理工大学 Multilayer displacement response determination method considering interlayer condition and transverse isotropy

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070074565A1 (en) * 2005-09-30 2007-04-05 Paramsothy Jayakumar System for virtual prediction of road loads
US20090070042A1 (en) * 2007-09-11 2009-03-12 Richard Birchwood Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state
CN108517735A (en) * 2018-04-12 2018-09-11 长沙理工大学 A kind of durable asphalt pavement design method and its pavement structure based on bimodular theory
CN110502845A (en) * 2019-08-27 2019-11-26 长沙理工大学 Method is determined based on the PFWD subgrade modulus of genetic algorithm and viscoelasticity theory
CN110658086A (en) * 2019-08-31 2020-01-07 长沙理工大学 Asphalt pavement load response analysis method considering tension-compression modulus difference
CN111475876A (en) * 2020-03-04 2020-07-31 长沙理工大学 Method for obtaining dynamic resilience mechanical characteristic parameters of granules
CN112214817A (en) * 2020-09-29 2021-01-12 长沙理工大学 Multilayer displacement response determination method considering interlayer condition and transverse isotropy

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JUNHUI ZHANG等: "Back-Calculation of Elastic Modulus of High Liquid Limit Clay Subgrades Based on Viscoelastic Theory and Multipopulation Genetic Algorithm", 《INTERNATIONAL JOURNAL OF GEOMECHANICS》, 31 October 2020 (2020-10-31), pages 1 - 11 *
单景松等: "路面动载响应模型研究", 《振动与冲击》 *
单景松等: "路面动载响应模型研究", 《振动与冲击》, no. 09, 25 September 2007 (2007-09-25), pages 89 - 93 *
曹宽: "荷载和变温综合作用下沥青路面非线性疲劳损伤研究", 《北方交通》, no. 10, 20 October 2018 (2018-10-20), pages 53 - 60 *
颜可珍等: "横观各向同性沥青路面结构动力响应分析", 《公路交通科技》, no. 05, 15 May 2017 (2017-05-15), pages 1 - 9 *

Also Published As

Publication number Publication date
CN112765825B (en) 2022-07-08

Similar Documents

Publication Publication Date Title
Babuska et al. Verification and validation in computational engineering and science: basic concepts
Ryckelynck et al. Multi-level a priori hyper-reduction of mechanical models involving internal variables
Ziyadi et al. Model uncertainty analysis using data analytics for life-cycle assessment (LCA) applications
Young Data-based mechanistic modelling, generalised sensitivity and dominant mode analysis
Stojanovic et al. Adaptive system for dam behavior modeling based on linear regression and genetic algorithms
Jensen et al. The use of updated robust reliability measures in stochastic dynamical systems
CN111324949B (en) Engineering structure flexibility recognition method considering noise influence
Zhang et al. Sobol's sensitivity analysis for a fuel cell stack assembly model with the aid of structure-selection techniques
Levenberg Inverse analysis of viscoelastic pavement properties using data from embedded instrumentation
Pelayo et al. Modal analysis based stress estimation for structural elements subjected to operational dynamic loadings
Orobio et al. Sampling-based sensitivity analysis of the mechanistic–empirical pavement design guide applied to material inputs
CN107341332A (en) A kind of determination method of the evaluation criterion weight of automobile motor drive system
Zou et al. Modal analysis of trough solar collector
CN112163669A (en) Pavement subsidence prediction method based on BP neural network
Shen et al. An optimized discrete grey multi-variable convolution model and its applications
Kassem et al. Quantification of the inherent uncertainty in the relaxation modulus and creep compliance of asphalt mixes
CN112765825B (en) Roadbed surface displacement response determination method considering modulus nonlinear distribution
CN101609411B (en) Embedded software power consumption BP neural network modeling method based on complexity
Coleri et al. Evaluation of laboratory, construction, and performance variability by bootstrapping and monte carlo methods for rutting performance prediction of heavy vehicle simulator test sections
He et al. A multi-scale wavelet finite element model for damage detection of beams under a moving load
Saleeb et al. Numerical simulation techniques for HMA rutting under loaded wheel tester
Hemez et al. The good, the bad and the ugly of predictive science
Song et al. Long-term deformation safety evaluation method of concrete dams based on the time-varying stability of concrete material
CN106709167A (en) Similarity theory-based assembling validity assessment method
Yun et al. Improved selfsim for inverse extraction of nonuniform, nonlinear, and inelastic material behavior under cyclic loadings

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant