CN109800466A - A kind of RBF analogy method of airport asphalt concrete pavement deformation process - Google Patents
A kind of RBF analogy method of airport asphalt concrete pavement deformation process Download PDFInfo
- Publication number
- CN109800466A CN109800466A CN201811572527.6A CN201811572527A CN109800466A CN 109800466 A CN109800466 A CN 109800466A CN 201811572527 A CN201811572527 A CN 201811572527A CN 109800466 A CN109800466 A CN 109800466A
- Authority
- CN
- China
- Prior art keywords
- error
- reflection parameters
- airport
- parameters range
- covering
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 239000011384 asphalt concrete Substances 0.000 title claims abstract description 39
- 230000008569 process Effects 0.000 title claims abstract description 35
- 238000012360 testing method Methods 0.000 claims abstract description 40
- 238000004458 analytical method Methods 0.000 claims abstract description 16
- 230000006870 function Effects 0.000 claims abstract description 16
- 239000010426 asphalt Substances 0.000 claims abstract description 13
- 230000001133 acceleration Effects 0.000 claims abstract description 11
- 238000009825 accumulation Methods 0.000 claims abstract description 10
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000010276 construction Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims description 32
- 230000000873 masking effect Effects 0.000 claims description 26
- 230000004044 response Effects 0.000 claims description 23
- 230000011514 reflex Effects 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 5
- 230000021615 conjugation Effects 0.000 claims description 4
- ONUFESLQCSAYKA-UHFFFAOYSA-N iprodione Chemical compound O=C1N(C(=O)NC(C)C)CC(=O)N1C1=CC(Cl)=CC(Cl)=C1 ONUFESLQCSAYKA-UHFFFAOYSA-N 0.000 claims description 4
- 238000000465 moulding Methods 0.000 claims description 4
- 238000002310 reflectometry Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000008901 benefit Effects 0.000 claims description 3
- 230000001934 delay Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 230000015271 coagulation Effects 0.000 claims 1
- 238000005345 coagulation Methods 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 239000002689 soil Substances 0.000 claims 1
- 238000013507 mapping Methods 0.000 abstract description 6
- 238000012423 maintenance Methods 0.000 description 6
- 238000005457 optimization Methods 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 2
- 239000004567 concrete Substances 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
The invention belongs to Aerodrome Construction technical fields, disclose a kind of RBF analogy method of airport asphalt concrete pavement deformation process, carry out the deformation analysis of airport asphalt concrete pavement using the ANN simulation algorithm of asphalt track plastic deformation accumulation；After the asphalt concrete pavement deformation analysis of airport, test section altitude data is obtained using the multiple vehicle-mounted laser sensors, acceleration transducer and range sensor that are embedded in simulation asphalt track again, by computer to current test point compared with the multiple test points in front, filter out maxima and minima, the very poor value between maximum and minimum is calculated, determines whether road surface deforms based on very poor value.Theoretical proof of the present invention, RBF network is the optimal network for completing mapping function in feedforward network；It is connected to the network weight and output is in a linear relationship；Classification capacity is good；Learning process fast convergence rate.
Description
Technical field
The invention belongs to Aerodrome Construction technical field more particularly to a kind of airport asphalt concrete pavement deformation processes
RBF analogy method.
Background technique
Currently, the prior art commonly used in the trade is such that
The rapid growth of aircraft enlargement and air traffic amount, the presentation for causing the breakage of airfield pavement structure more obvious
In face of us.From the point of view of current airport operating conditions both domestic and external, considerable airfield pavement uses year in not up to design
Before limit, pavement structure with regard to it is premature occur crack, deformation and lose use value, can usually bring huge maintenance pressure and warp
Ji burden.
Boeing in 1992 starts to be dedicated to developing big aircraft, and can airport pavement conditions at that time be faced with and bear
The technical problems such as big aircraft loads, in this context, U.S. FAA and Boeing establish the association of a joint study and exploitation
View, has set up US airports road face test center NAPTF, is substantially carried out the full scale test research of different pavement.
NAPTF is located at the William Hughes technique center of Atlantic Ocean city International airport, and NAPTF is designed as an indoor straight line
Formula load test center.
It is tested around FAA full size airfield pavement, people from different perspectives carry out airfield pavement using different theories method
Effective analysis and research.
Guo Hua has carried out series of statistical analysis to NAPTF airfield runway test data, and has studied airport pavement plate emphatically
Curling problem；Kim et al. has carried out axial symmetry and three-dimensional to NAPTF bituminous concrete airfield runway application ABAQUS respectively to be had
Limit Meta Model analysis；You Qinglong, Ling Jianming etc. consider that with B777-300ER, A380-800 be the new generation large aircraft system of representative
The Airport Asphalt mechanical response analysis of power, by ABAQUS Three-D limited meta software, establish flexible pavement structure and
Semi-rigid road face structural finite element model.
Up to the present, American National road face test center has carried out nine independent pilot projects, and what it is for analysis is
MFC (medium strength flexible in the airport FAA asphalt concrete pavement permanent deformation accumulation test CC1
Conventional) pilot project, i.e., the bituminous concrete track testing project of medium subgrade strength.
The permanent deformation accumulation test of airport flexible pavement: the single wheel loads of CC1 test are 204.116kN (45000 pounds), tire
Pressure is 1303kPa.Aircraft load test carries out under Automatic Control mode, and the reciprocal travel speed that loads is set as definite value
8km/h, this speed are equivalent to the speed of coast period of the aircraft from cabin to before taking off.Laterally offset amount is controlled in 1524mm
Within.US airports road face test center is mounted with 1050 sensors in total on the airfield runway of test, for measuring it
Water content, temperature, wheel load related with stress and deformation.When test vehicle is moved along test runway, counted with one
The data collection system of calculation machine control carrys out the related datas such as response and the performance of automatic collection runway.
In conclusion problem of the existing technology is:
In the detection of prior art airport asphalt concrete pavement deformation process, best approximation properties, and part cannot be obtained
There are problems；It is poor to output and input mapping function；It is connected to the network weight and output is in a linear relationship indefinite；Classification capacity
Difference；Learning process convergence rate is slower.
The data analysis of safe, smooth, fast express agency cannot be provided in the prior art for airport, cause air transportation efficiency
Difference, airfield runway service life is short, wastes a large amount of maintenances, maintenance costs.
Summary of the invention
In view of the problems of the existing technology, the present invention provides a kind of airport asphalt concrete pavement deformation processes
RBF analogy method,
The invention is realized in this way a kind of RBF analogy method of airport asphalt concrete pavement deformation process, the machine
The RBF analogy method of asphalt concrete pavement deformation process includes:
Step 1 carries out the change of airport asphalt concrete pavement using the ANN simulation algorithm of asphalt track plastic deformation accumulation
Conformal analysis；
Step 2, after the asphalt concrete pavement deformation analysis of airport, then using the multiple vehicles for being embedded in simulation asphalt track
It carries laser sensor, acceleration transducer and range sensor and obtains the every 1mm section relative elevation SOURCE data in test section,
Each test point obtains a section relative elevation, by computer to current test point compared with the multiple test points in front
Compared with filtering out maxima and minima, calculate the very poor value between maximum and minimum, whether become based on very poor value judgement road surface
Shape；
When very poor value is less than the standard of 10mm, then there is no deformations on road surface；Road surface is then determined when very poor value is greater than 10mm
There are deformations.
Further, the ANN simulation algorithm of the asphalt track plastic deformation accumulation are as follows:
gi=Rg(fi-1,fi,fi+1；eg,tg)
ki=Rk(fi-1,fi,fi+1；ek,tk)
hi=Rh(fi-1,fi,fi+1；eh,th)
x1< x2< L < xi< L < xn
Wherein, 1 < i < n-2；
Further, every using multiple vehicle-mounted laser sensors, acceleration transducer and range sensor acquisition test section
In 1mm section relative elevation SOURCE data, the detection model of vehicle-mounted laser sensor are as follows:
Wherein, τ is delay skew, and f is Doppler frequency shift, 0 < a, b < α/2, x*(t) conjugation for indicating x (t), as x (t)
When for real signal, x (t)< p >=| x (t) |< p >sgn(x(t))；When x (t) is time multiplexed signal, [x (t)]< p >=| x (t) |p-1x*
(t)；
Acceleration transducer receives signal y (t) and indicates are as follows:
Y (t)=x (t)+n (t)；
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the parsing shape of x (t)
Formula indicates are as follows:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2, Λ, M-1, M are
Order of modulation, an=ej2πε/M, ε=0,1,2, Λ, M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIt indicates
Carrier frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π]；
The signal model of range sensor indicates are as follows:
R (t)=x1(t)+x2(t)+L+xn(t)+v(t)
Wherein, xiIt (t) is each signal component of time-frequency overlapped signal, each component signal is independently uncorrelated, and n is time-frequency weight
The number of folded signal component, θkiIndicate the modulation to each signal component carrier phase, fciFor carrier frequency, AkiBelieve for i-th
Amplitude number at the k moment, TsiFor Baud Length；
Computer compared with the multiple test points in front, filters out maxima and minima to current test point, calculates most
Very poor value greatly between minimum determines whether road surface deforms based on very poor value, specifically includes:
Step 1, using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse and rings
It answers, reflection parameters correspond to Time Domain Reflectometry response, and configured transmission corresponds to time-domain transmission response；
Step 2 constructs four according to the position of the first two pulse on a timeline in reflex response and transmission response respectively
A time domain gating function；
Step 3 selects the reflex response and transmission response of time domain in step 1 using gate function in step 2
It is logical, the first two pulse in reflex response and transmission response is extracted respectively；
Time domain impulse after time domain gating is passed through Fourier transformation respectively, obtains frequency domain gated data by step 4；
Step 5 contains the information of circuit-under-test in frequency domain gated data, using obtained information, according to formula construction
Compensation factor Fcf(i)；Using following formula, compensation factor F is constructedcf(i):
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step；
R (i) is ratio factor；
Fcf(i) compensation factor；
Step 6 is not covered the reflection parameters F of error using occlusion compensation formulaS11(i) and configured transmission FS21
(i)；
Step 7, the reflection parameters range for not covering error are not covered according to the topological order of PI-PO for each
The reflection parameters range of error generates all feasible divisions of k-, and it is feasible that each reflection parameters range without covering error generates k-
Divide specific progress according to the following formula:
Wherein,Represent with ,+represents or, k be LUT input limit, input (v) expression without masking error reflection
The input set of parameter area v, u are without covering one in the reflection parameters range v input set of error without masking error
Reflection parameters range, f (k, v) indicate the feasible division of all k- of the reflection parameters range v without covering error, and f (k, u) is indicated
The feasible division of all k- of the reflection parameters range u of error is not covered.
Further, it needs to follow the steps below before step 1:
The measurement parameter of vector network analyzer is set as needed first, obtains the scattering parameter of tested network entirety,
Including reflection parameters and configured transmission；
Secondly scattering parameter is pre-processed, zero padding is carried out to data sequence, carries out Fast Fourier Transform (FFT)；According to Fu
In after leaf inverse transformation the distribution situation and resolution ratio of time domain impulse requirement, select different window functions to Fourier transformation before
Data are handled.
Further, step 7 specifically includes:
Traversal initialization forward: the first step institute is initialized as either with or without the depth of reflection parameters range for covering error
0, corresponding area stream is initialized as 0；The depth that PI does not cover the reflection parameters range output side of error is initialized as 1, right
The area stream answered is initialized as 0；
Second step judges that all reflection parameters ranges without covering error have been accessed, if so, end step
Third step is gone to, is otherwise continued to execute: being traversed forward according to the topological order of PI-PO, is taken not visited without covering error
Reflection parameters range v, it is no masking error all divisions for meeting following formula of reflection parameters range v in selection area stream most
Small division XvAs optimal dividing:
depth(x)≤Odepth-height(v)；
X indicates a division without the reflection parameters range v for covering error, and depth (x) indicates to divide the depth of x,
Odepth indicates optimal depth, and height (v) indicates the height of the reflection parameters range v without covering error；
Wherein, the area stream calculation for not covering the division of the reflection parameters range v of error is as follows:
Wherein, ξ is arbitrarily small random number, and iedge (v) indicates the input of the reflection parameters range v without covering error
Line set, AvIndicate influence of the reflection parameters range v without covering error to area itself；
Update is respectively depth (X without the depth and area stream for covering the reflection parameters range v of errorv) and af (Xv)；
The depth for updating any one output side e of the reflection parameters range v without masking error is dept (hXv)+
Dela (ye), area stream areWherein, delay (e) indicates the associated time delays of output side e, and oedge (v) expression does not have
There is the output line set of the reflection parameters range v of masking error；
Third step, traversal initialization backward: initialization root collection is combined into all reflection parameters ranges without covering error,
And by be initialized as 1 either with or without the height of reflection parameters range of masking error；
Judge whether either with or without masking error reflection parameters range be accessed, if it is, terminate, otherwise after
It is continuous to execute: according to the inverse topological order of PO-PI, the not visited reflection parameters range v without covering error to be taken out, if v
In set root, calculate:
H=max { height (e): e ∈ oedge (v) }；
Wherein, height (e) is appointing in the output line set oedge (v) without the reflection parameters range v of masking error
The height on meaning side, h is then the maximum of the height on all sides in the output line set for do not cover the reflection parameters range v of error
Value；
Update the optimal dividing X obtained without the reflection parameters range v for covering errorvInterior is any without covering error
The height of reflection parameters range u is height (u)=max { height (u), h }, for XvAny input side e update its height
It spends height (e)=max { height (e), delay (e)+h }, update set root is root ∪ inode (Xv), inode (Xv)
Indicate support without the optimal dividing X of the reflection parameters range v of masking errorvThe tail on input side do not cover the reflection of error
Parameter area.
Another object of the present invention is to provide a kind of RBF moulds for realizing the airport asphalt concrete pavement deformation process
The computer program of quasi- method.
Another object of the present invention is to provide a kind of RBF moulds for realizing the airport asphalt concrete pavement deformation process
The information data processing terminal of quasi- method.
Another object of the present invention is to provide a kind of computer readable storage mediums, including instruction, when it is in computer
When upper operation, so that computer executes the RBF analogy method of the airport asphalt concrete pavement deformation process.
Another object of the present invention is to provide a kind of RBF moulds for realizing the airport asphalt concrete pavement deformation process
The RBF simulation system of the airport asphalt concrete pavement deformation process of quasi- method.
In conclusion advantages of the present invention and good effect are as follows:
A kind of RBF analogy method of airport asphalt concrete pavement deformation process provided by the invention, it has uniquely most
The good characteristic approached, and exist without local minimum problem；RBF neural have it is stronger output and input mapping function, and
And theoretical proof RBF network in feedforward network is the optimal network for completing mapping function；It is connected to the network weight and output is in line
Sexual intercourse；Classification capacity is good；Learning process fast convergence rate.The present invention can support to improve airfield pavement using quality, road face
Shield and maintenance provide necessary reference, provide for airport security, smooth, fast express agency and guarantee data analysis, improve aviation and transport
Movement Capabilities promote economic fast development；Extend airfield runway service life, saves a large amount of maintenances, maintenance costs.
The present invention is based on heuritic approaches, and technology-mapped is divided into logic optimization and structure optimization, structure optimization part
Using DAG model, it is divided into division and generates, divides selection and LUT three steps of mapping, divide the think of for generating and using Dynamic Programming
Think, the reflection parameters range for not covering error for each quickly generates all feasible divisions of k-；Selection is divided based on one kind
The number of iterations can be with the heuristic thought of the iteration of adaptively changing, by the iteration for repeatedly traversing forward with traversing backward, constantly
Optimisation technique mapping as a result, and reduce unnecessary iteration, optimize the efficiency of technology-mapped, final choice go out to be delayed and
Area optimised division set simultaneously, compared to existing delay and area algorithm, the present invention no matter technology-mapped matter
It all improves a lot in the efficiency of amount or technology-mapped；Number can be simulated for the RBF of airport asphalt concrete pavement deformation process
According to acquisition provide necessary condition.
The present invention has modified the reflection parameters range area stream calculation formula for not covering error, improves and divides selection
Randomness.
The detection model of vehicle-mounted laser sensor of the present invention are as follows:
Wherein, τ is delay skew, and f is Doppler frequency shift, 0 < a, b < α/2, x*(t) conjugation for indicating x (t), as x (t)
When for real signal, x (t)< p >=| x (t) |< p >sgn(x(t))；When x (t) is time multiplexed signal, [x (t)]< p >=| x (t) |p-1x*
(t)；
Acceleration transducer receives signal y (t) and indicates are as follows:
Y (t)=x (t)+n (t)；
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the parsing shape of x (t)
Formula indicates are as follows:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2, Λ, M-1, M are
Order of modulation, an=ej2πε/M, ε=0,1,2, Λ, M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIt indicates
Carrier frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π]；
The signal model of range sensor indicates are as follows:
R (t)=x1(t)+x2(t)+L+xn(t)+v(t)
Wherein, xiIt (t) is each signal component of time-frequency overlapped signal, each component signal is independently uncorrelated, and n is time-frequency weight
The number of folded signal component, θkiIndicate the modulation to each signal component carrier phase, fciFor carrier frequency, AkiBelieve for i-th
Amplitude number at the k moment, TsiFor Baud Length；
Computer compared with the multiple test points in front, filters out maxima and minima to current test point, calculates most
Very poor value greatly between minimum determines whether road surface deforms based on very poor value, specifically includes:
Step 1, using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse and rings
It answers, reflection parameters correspond to Time Domain Reflectometry response, and configured transmission corresponds to time-domain transmission response；
Step 2 constructs four according to the position of the first two pulse on a timeline in reflex response and transmission response respectively
A time domain gating function；
Step 3 selects the reflex response and transmission response of time domain in step 1 using gate function in step 2
It is logical, the first two pulse in reflex response and transmission response is extracted respectively；
Time domain impulse after time domain gating is passed through Fourier transformation respectively, obtains frequency domain gated data by step 4；
Step 5 contains the information of circuit-under-test in frequency domain gated data, using obtained information, according to formula construction
Compensation factor Fcf(i)；Using following formula, compensation factor F is constructedcf(i):
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step；
R (i) is ratio factor；
Fcf(i) compensation factor；
Step 6 is not covered the reflection parameters F of error using occlusion compensation formulaS11(i) and configured transmission FS21
(i)；
Step 7, the reflection parameters range for not covering error are not covered according to the topological order of PI-PO for each
The reflection parameters range of error generates all feasible divisions of k-, and it is feasible that each reflection parameters range without covering error generates k-
Division specifically carries out according to the following formula, the RBF analogue data for being embodied as airport asphalt concrete pavement deformation process of above scheme
The acquisition of accuracy provides guarantee.
Detailed description of the invention
Fig. 1 is the RBF analogy method flow chart of asphalt concrete pavement deformation process in airport provided in an embodiment of the present invention.
Fig. 2 is Dynamic Simulation Results analysis schematic diagram provided in an embodiment of the present invention.
Fig. 3 is dynamic thickness and Amplitude Ration analysis schematic diagram provided in an embodiment of the present invention.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to
Limit the present invention.
The RBF analogy method of asphalt concrete pavement deformation process in airport provided in an embodiment of the present invention mainly includes
The ANN simulation algorithm of asphalt track plastic deformation accumulation.
Specifically have:
Such as Fig. 1, the RBF analogy method of asphalt concrete pavement deformation process in airport provided in an embodiment of the present invention is described
The RBF analogy method of airport asphalt concrete pavement deformation process includes:
S101: asphalt concrete pavement deformation in airport is carried out using the ANN simulation algorithm of asphalt track plastic deformation accumulation
Analysis；
S102: it after the asphalt concrete pavement deformation analysis of airport, then uses and is embedded in the multiple vehicle-mounted of simulation asphalt track
Laser sensor, acceleration transducer and range sensor obtain the every 1mm section relative elevation SOURCE data in test section, often
One test point obtains a section relative elevation, by computer to current test point compared with the multiple test points in front,
Maxima and minima is filtered out, the very poor value between maximum and minimum is calculated, determines whether road surface deforms based on very poor value；When
Then deformation is not present in road surface when very poor value is less than the standard of 10mm；Then determine that road surface has deformation when very poor value is greater than 10mm.
The ANN simulation algorithm of the asphalt track plastic deformation accumulation are as follows:
gi=Rg(fi-1,fi,fi+1；eg,tg)
ki=Rk(fi-1,fi,fi+1；ek,tk)
hi=Rh(fi-1,fi,fi+1；eh,th)
x1< x2< L < xi< L < xn
Wherein, 1 < i < n-2；
Further, every using multiple vehicle-mounted laser sensors, acceleration transducer and range sensor acquisition test section
In 1mm section relative elevation SOURCE data, the detection model of vehicle-mounted laser sensor are as follows:
Wherein, τ is delay skew, and f is Doppler frequency shift, 0 < a, b < α/2, x*(t) conjugation for indicating x (t), as x (t)
When for real signal, x (t)< p >=| x (t) |< p >sgn(x(t))；When x (t) is time multiplexed signal, [x (t)]< p >=| x (t) |p-1x*
(t)；
Acceleration transducer receives signal y (t) and indicates are as follows:
Y (t)=x (t)+n (t)；
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the parsing shape of x (t)
Formula indicates are as follows:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2, Λ, M-1, M are
Order of modulation, an=ej2πε/M, ε=0,1,2, Λ, M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIt indicates
Carrier frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π]；
The signal model of range sensor indicates are as follows:
R (t)=x1(t)+x2(t)+L+xn(t)+v(t)
Wherein, xiIt (t) is each signal component of time-frequency overlapped signal, each component signal is independently uncorrelated, and n is time-frequency weight
The number of folded signal component, θkiIndicate the modulation to each signal component carrier phase, fciFor carrier frequency, AkiBelieve for i-th
Amplitude number at the k moment, TsiFor Baud Length；
Computer compared with the multiple test points in front, filters out maxima and minima to current test point, calculates most
Very poor value greatly between minimum determines whether road surface deforms based on very poor value, specifically includes:
Step 1, using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse and rings
It answers, reflection parameters correspond to Time Domain Reflectometry response, and configured transmission corresponds to time-domain transmission response；
Step 2 constructs four according to the position of the first two pulse on a timeline in reflex response and transmission response respectively
A time domain gating function；
Step 3 selects the reflex response and transmission response of time domain in step 1 using gate function in step 2
It is logical, the first two pulse in reflex response and transmission response is extracted respectively；
Time domain impulse after time domain gating is passed through Fourier transformation respectively, obtains frequency domain gated data by step 4；
Step 5 contains the information of circuit-under-test in frequency domain gated data, using obtained information, according to formula construction
Compensation factor Fcf(i)；Using following formula, compensation factor F is constructedcf(i):
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step；
R (i) is ratio factor；
Fcf(i) compensation factor；
Step 6 is not covered the reflection parameters F of error using occlusion compensation formulaS11(i) and configured transmission FS21
(i)；
Step 7, the reflection parameters range for not covering error are not covered according to the topological order of PI-PO for each
The reflection parameters range of error generates all feasible divisions of k-, and it is feasible that each reflection parameters range without covering error generates k-
Divide specific progress according to the following formula:
Wherein,Represent with ,+represents or, k be LUT input limit, input (v) expression without masking error reflection
The input set of parameter area v, u are without covering one in the reflection parameters range v input set of error without masking error
Reflection parameters range, f (k, v) indicate the feasible division of all k- of the reflection parameters range v without covering error, and f (k, u) is indicated
The feasible division of all k- of the reflection parameters range u of error is not covered.
It needs to follow the steps below before step 1:
The measurement parameter of vector network analyzer is set as needed first, obtains the scattering parameter of tested network entirety,
Including reflection parameters and configured transmission；
Secondly scattering parameter is pre-processed, zero padding is carried out to data sequence, carries out Fast Fourier Transform (FFT)；According to Fu
In after leaf inverse transformation the distribution situation and resolution ratio of time domain impulse requirement, select different window functions to Fourier transformation before
Data are handled.
Step 7 specifically includes:
Traversal initialization forward: the first step institute is initialized as either with or without the depth of reflection parameters range for covering error
0, corresponding area stream is initialized as 0；The depth that PI does not cover the reflection parameters range output side of error is initialized as 1, right
The area stream answered is initialized as 0；
Second step judges that all reflection parameters ranges without covering error have been accessed, if so, end step
Third step is gone to, is otherwise continued to execute: being traversed forward according to the topological order of PI-PO, is taken not visited without covering error
Reflection parameters range v, it is no masking error all divisions for meeting following formula of reflection parameters range v in selection area stream most
Small division XvAs optimal dividing:
depth(x)≤Odepth-height(v)；
X indicates a division without the reflection parameters range v for covering error, and depth (x) indicates to divide the depth of x,
Odepth indicates optimal depth, and height (v) indicates the height of the reflection parameters range v without covering error；
Wherein, the area stream calculation for not covering the division of the reflection parameters range v of error is as follows:
Wherein, ξ is arbitrarily small random number, and iedge (v) indicates the input of the reflection parameters range v without covering error
Line set, AvIndicate influence of the reflection parameters range v without covering error to area itself；
Update is respectively depth (X without the depth and area stream for covering the reflection parameters range v of errorv) and af (Xv)；
The depth for updating any one output side e of the reflection parameters range v without masking error is dept (hXv)+
Dela (ye), area stream areWherein, delay (e) indicates the associated time delays of output side e, and oedge (v) expression does not have
Cover the output line set of the reflection parameters range v of error；
Third step, traversal initialization backward: initialization root collection is combined into all reflection parameters ranges without covering error,
And by be initialized as 1 either with or without the height of reflection parameters range of masking error；
Judge whether either with or without masking error reflection parameters range be accessed, if it is, terminate, otherwise after
It is continuous to execute: according to the inverse topological order of PO-PI, the not visited reflection parameters range v without covering error to be taken out, if v
In set root, calculate:
H=max { height (e): e ∈ oedge (v) }；
Wherein, height (e) is appointing in the output line set oedge (v) without the reflection parameters range v of masking error
The height on meaning side, h is then the maximum of the height on all sides in the output line set for do not cover the reflection parameters range v of error
Value；
Update the optimal dividing X obtained without the reflection parameters range v for covering errorvInterior is any without covering error
The height of reflection parameters range u is height (u)=max { height (u), h }, for XvAny input side e update its height
It spends height (e)=max { height (e), delay (e)+h }, update set root is root ∪ inode (Xv), inode (Xv)
Indicate support without the optimal dividing X of the reflection parameters range v of masking errorvThe tail on input side do not cover the reflection of error
Parameter area.
Invention is further explained with reference to the accompanying drawing.
As shown in Fig. 2, Dynamic Simulation Results analysis chart
As shown in figure 3, dynamic thickness and Amplitude Ration analysis chart.
(a) 1,2 layer of ratio and 1, the 3 layer of value of ratio at 0.985-1.025；
(b) 1,2 layer of ratio and 1, the 3 layer of value of ratio at 0.88-1；
(c) 1,2 layer of ratio and 1, the 3 layer of value of ratio at 1-1.09；
(d) 1,2 layer of ratio and 1,3 layer of ratio and 2, the 3 layers of value of ratio at 0-60.
In the above-described embodiments, can come wholly or partly by software, hardware, firmware or any combination thereof real
It is existing.When using entirely or partly realizing in the form of a computer program product, the computer program product include one or
Multiple computer instructions.When loading on computers or executing the computer program instructions, entirely or partly generate according to
Process described in the embodiment of the present invention or function.The computer can be general purpose computer, special purpose computer, computer network
Network or other programmable devices.The computer instruction may be stored in a computer readable storage medium, or from one
Computer readable storage medium is transmitted to another computer readable storage medium, for example, the computer instruction can be from one
A web-site, computer, server or data center pass through wired (such as coaxial cable, optical fiber, Digital Subscriber Line (DSL)
Or wireless (such as infrared, wireless, microwave etc.) mode is carried out to another web-site, computer, server or data center
Transmission).The computer-readable storage medium can be any usable medium or include one that computer can access
The data storage devices such as a or multiple usable mediums integrated server, data center.The usable medium can be magnetic Jie
Matter, (for example, floppy disk, hard disk, tape), optical medium (for example, DVD) or semiconductor medium (such as solid state hard disk Solid
State Disk (SSD)) etc..
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Made any modifications, equivalent replacements, and improvements etc., should be included in X of the present invention within mind and principle1Protection scope within.
Claims (9)
1. a kind of RBF analogy method of airport asphalt concrete pavement deformation process, which is characterized in that the airport pitch coagulation
The RBF analogy method of soil surface deformation process includes:
Step 1 carries out the deformation point of airport asphalt concrete pavement using the ANN simulation algorithm of asphalt track plastic deformation accumulation
Analysis；
Step 2 after the asphalt concrete pavement deformation analysis of airport, then uses and is embedded in the multiple vehicle-mounted sharp of simulation asphalt track
Optical sensor, acceleration transducer and range sensor obtain the every 1mm section relative elevation SOURCE data in test section, each
A test point obtains a section relative elevation, by computer to current test point compared with the multiple test points in front, sieve
Maxima and minima is selected, the very poor value between maximum and minimum is calculated, determines whether road surface deforms based on very poor value；
When very poor value is less than the standard of 10mm, then there is no deformations on road surface；Then determine that road surface exists when very poor value is greater than 10mm
Deformation.
2. the RBF analogy method of asphalt concrete pavement deformation process in airport as described in claim 1, which is characterized in that institute
State the ANN simulation algorithm of asphalt track plastic deformation accumulation are as follows:
gi=Rg(fi-1,fi,fi+1；eg,tg)
ki=Rk(fi-1,fi,fi+1；ek,tk)
hi=Rh(fi-1,fi,fi+1；eh,th)
x1< x2< L < xi< L < xn
Wherein, 1 < i < n-2；
3. the RBF analogy method of asphalt concrete pavement deformation process in airport as described in claim 1, which is characterized in that adopt
The every 1mm section relative elevation in test section is obtained with multiple vehicle-mounted laser sensors, acceleration transducer and range sensor
In SOURCE data, the detection model of vehicle-mounted laser sensor are as follows:
Wherein, τ is delay skew, and f is Doppler frequency shift, 0 < a, b < α/2, x*(t) conjugation for indicating x (t), when x (t) is real
When signal, x (t)< p >=| x (t) |< p >sgn(x(t))；When x (t) is time multiplexed signal, [x (t)]< p >=| x (t) |p-1x*(t)；
Acceleration transducer receives signal y (t) and indicates are as follows:
Y (t)=x (t)+n (t)；
Wherein, x (t) is digital modulation signals, and n (t) is the impulsive noise of obedience standard S α S distribution, the analytical form table of x (t)
It is shown as:
Wherein, N is sampling number, anFor the information symbol of transmission, in MASK signal, an=0,1,2, Λ, M-1, M are modulation
Order, an=ej2πε/M, ε=0,1,2, Λ, M-1, g (t) expression rectangle molding pulse, TbIndicate symbol period, fcIndicate carrier wave
Frequency, carrier wave initial phaseIt is the equally distributed random number in [0,2 π]；
The signal model of range sensor indicates are as follows:
R (t)=x1(t)+x2(t)+L+xn(t)+v(t)
Wherein, xiIt (t) is each signal component of time-frequency overlapped signal, each component signal is independently uncorrelated, and n is time-frequency overlapping letter
The number of number component, θkiIndicate the modulation to each signal component carrier phase, fciFor carrier frequency, AkiExist for i-th of signal
The amplitude at k moment, TsiFor Baud Length；
Computer compared with the multiple test points in front, filters out maxima and minima to current test point, calculate it is maximum with
Very poor value between minimum determines whether road surface deforms based on very poor value, specifically includes:
Step 1, using inverse Fourier transform, the scattering parameter that frequency domain is measured transforms to time domain, obtains time domain impulse response,
Reflection parameters correspond to Time Domain Reflectometry response, and configured transmission corresponds to time-domain transmission response；
Step 2, according to the position of the first two pulse on a timeline in reflex response and transmission response, when constructing four respectively
Domain gate function；
Step 3 gates the reflex response and transmission response of time domain in step 1 using gate function in step 2, point
Indescribably take out the first two pulse in reflex response and transmission response；
Time domain impulse after time domain gating is passed through Fourier transformation respectively, obtains frequency domain gated data by step 4；
Step 5 contains the information of circuit-under-test in frequency domain gated data, using obtained information, is compensated according to formula construction
Factor Fcf(i)；Using following formula, compensation factor F is constructedcf(i):
G1(i)~G4It (i) is the frequency domain gated data obtained in upper step；
R (i) is ratio factor；
Fcf(i) compensation factor；
Step 6 is not covered the reflection parameters F of error using occlusion compensation formulaS11(i) and configured transmission FS21(i)；
Step 7, the reflection parameters range for not covering error do not cover error according to the topological order of PI-PO for each
Reflection parameters range generate all feasible divisions of k-, each reflection parameters range generation feasible division of k- without masking error
Specifically carry out according to the following formula:
Wherein,Represent with ,+represents or, k be LUT input limit, input (v) expression without masking error reflection parameters
The input set of range v, u are without covering one in the reflection parameters range v input set of the error reflection without covering error
Parameter area, f (k, v) indicate that the feasible division of all k- of the reflection parameters range v without covering error, f (k, u) expression do not have
Cover the feasible division of all k- of the reflection parameters range u of error.
4. the RBF analogy method of asphalt concrete pavement deformation process in airport as claimed in claim 3, which is characterized in that
It needs to follow the steps below before step 1:
The measurement parameter of vector network analyzer is set as needed first, obtains the scattering parameter of tested network entirety, including
Reflection parameters and configured transmission；
Secondly scattering parameter is pre-processed, zero padding is carried out to data sequence, carries out Fast Fourier Transform (FFT)；According to Fourier
The requirement of the distribution situation and resolution ratio of time domain impulse, selects different window functions to the data before Fourier transformation after inverse transformation
It is handled.
5. the RBF analogy method of asphalt concrete pavement deformation process in airport as claimed in claim 3, which is characterized in that step
Rapid seven specifically include:
The first step, traversal initialization forward: being initialized as 0 either with or without the depth for the reflection parameters range for covering error for institute, right
The area stream answered is initialized as 0；The depth that PI does not cover the reflection parameters range output side of error is initialized as 1, corresponding
Area stream is initialized as 0；
Second step judges that all reflection parameters ranges without covering error have been accessed, if so, end step is gone to
Otherwise third step continues to execute: traversing forward, take not visited without the anti-of masking error according to the topological order of PI-PO
Parameter area v is penetrated, selection area stream is the smallest in all divisions for meeting following formula of reflection parameters range v of no masking error
Divide XvAs optimal dividing:
depth(x)≤Odepth-height(v)；
X indicates a division without the reflection parameters range v for covering error, and depth (x) indicates to divide the depth of x, Odepth
Indicate optimal depth, height (v) indicates the height of the reflection parameters range v without covering error；
Wherein, the area stream calculation for not covering the division of the reflection parameters range v of error is as follows:
Wherein, ξ is arbitrarily small random number, and iedge (v) indicates the input side collection of the reflection parameters range v without covering error
It closes, AvIndicate influence of the reflection parameters range v without covering error to area itself；
Update is respectively depth (X without the depth and area stream for covering the reflection parameters range v of errorv) and af (Xv)；
The depth for updating any one output side e of the reflection parameters range v without masking error is dept (hXv)+dela
(ye), area stream isWherein, delay (e) indicates the associated time delays of output side e, and oedge (v) expression is not covered
The output line set of the reflection parameters range v of error；
Third step, traversal initialization backward: initialization root collection is combined into all reflection parameters ranges without covering error, and will
Be initialized as 1 either with or without the height of reflection parameters range of masking error；
Judge whether either with or without masking error reflection parameters range be accessed, if it is, terminate, otherwise continue to hold
Row: according to the inverse topological order of PO-PI, the not visited reflection parameters range v without covering error is taken out, if v is collecting
It closes in root, calculates:
H=max { height (e): e ∈ oedge (v) }；
Wherein, height (e) is any side in the output line set oedge (v) without the reflection parameters range v of masking error
Height, h is then the maximum value of the height on all sides in the output line set for do not cover the reflection parameters range v of error；
Update the optimal dividing X obtained without the reflection parameters range v for covering errorvInterior any reflection without covering error
The height of parameter area u is height (u)=max { height (u), h }, for XvAny input side e update its height
Height (e)=max { height (e), delay (e)+h }, update set root are root ∪ inode (Xv), inode (Xv) table
Show support without the optimal dividing X of the reflection parameters range v of masking errorvInput side tail do not cover error reflection ginseng
Number range.
6. a kind of RBF analogy method for realizing asphalt concrete pavement deformation process in airport described in Claims 1 to 5 any one
Computer program.
7. a kind of RBF analogy method for realizing asphalt concrete pavement deformation process in airport described in Claims 1 to 5 any one
Information data processing terminal.
8. a kind of computer readable storage medium, including instruction, when run on a computer, so that computer is executed as weighed
Benefit requires the RBF analogy method of asphalt concrete pavement deformation process in airport described in 1-5 any one.
9. a kind of RBF analogy method for realizing asphalt concrete pavement deformation process in airport described in Claims 1 to 5 any one
Airport asphalt concrete pavement deformation process RBF analog control system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811572527.6A CN109800466A (en) | 2018-12-21 | 2018-12-21 | A kind of RBF analogy method of airport asphalt concrete pavement deformation process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811572527.6A CN109800466A (en) | 2018-12-21 | 2018-12-21 | A kind of RBF analogy method of airport asphalt concrete pavement deformation process |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109800466A true CN109800466A (en) | 2019-05-24 |
Family
ID=66557374
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811572527.6A Pending CN109800466A (en) | 2018-12-21 | 2018-12-21 | A kind of RBF analogy method of airport asphalt concrete pavement deformation process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109800466A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104931912A (en) * | 2015-06-12 | 2015-09-23 | 西安电子科技大学 | Time domain compensation method of vector network analyzer |
| CN105488285A (en) * | 2015-12-04 | 2016-04-13 | 西安电子科技大学 | Efficient FPGA technology mapping algorithm |
| CN107102593A (en) * | 2017-06-21 | 2017-08-29 | 武汉万千无限科技有限公司 | A kind of intelligent remote video monitoring control system based on computer internet technology |
-
2018
- 2018-12-21 CN CN201811572527.6A patent/CN109800466A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104931912A (en) * | 2015-06-12 | 2015-09-23 | 西安电子科技大学 | Time domain compensation method of vector network analyzer |
| CN105488285A (en) * | 2015-12-04 | 2016-04-13 | 西安电子科技大学 | Efficient FPGA technology mapping algorithm |
| CN107102593A (en) * | 2017-06-21 | 2017-08-29 | 武汉万千无限科技有限公司 | A kind of intelligent remote video monitoring control system based on computer internet technology |
Non-Patent Citations (2)
| Title |
|---|
| 尹福成: "机场道面变形与开裂模式人工智能分析方法", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅱ辑》 * |
| 时宁等: "沥青路面变形类病害快速检测方法", 《广西大学学报:自然科学版》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Filippone | Aircraft noise prediction | |
| Stephan et al. | Hybrid simulation of wake-vortex evolution during landing on flat terrain and with plate line | |
| CN103513235B (en) | Clear sky aircraft wake stable section radar scattering characteristic computing method | |
| CN115525978B (en) | Speed inlet numerical wave generation method for hydrodynamic analysis of wave condition of aircraft | |
| Broeren et al. | Low-Reynolds number aerodynamics of an 8.9% scale semispan swept wing for assessment of icing effects | |
| Loubeau et al. | Effects of meteorological variability on sonic boom propagation from hypersonic aircraft | |
| CN115752708A (en) | Airport single-point noise prediction method based on deep time convolution network | |
| Botero et al. | Aerodynamic verification and validation of SUAVE | |
| CN109800466A (en) | A kind of RBF analogy method of airport asphalt concrete pavement deformation process | |
| Hrastovec et al. | Machine learning model for aircraft performances | |
| Busch et al. | Aerodynamic simulation of a horn-ice accretion on a subscale model | |
| Zhou et al. | Development of a Real-Time In-Flight Ice Detection System via Computational Aeroacoustics and Bayesian Neural Networks | |
| Park et al. | Summary of the 2008 NASA fundamental aeronautics program sonic boom prediction workshop | |
| CN105118332A (en) | Air traffic control analog simulation abnormality detection method and device based on clustering analysis method | |
| Danielle et al. | Model for the calculation of the radar cross section of wake vortices of take-off and landing airplanes | |
| CN104698444B (en) | Method for determining influences of array element position errors on three-dimensional imaging quality of array antenna | |
| Lockard et al. | Influence of spanwise boundary conditions on slat noise simulations | |
| Ikuta et al. | Noise prediction for aircraft during approach by machine learning using measured sound source spectra, flight parameters, and aircraft specifications | |
| CN110488237B (en) | Method for classifying number of unmanned aerial vehicle rotors based on frequency modulation continuous wave radar | |
| Jacob et al. | Toward hybrid CAA with ground effects | |
| Wiberg | Large-scale swept-wing ice accretion modeling in the NASA Glenn Icing Research Tunnel using LEWICE3D | |
| Lei et al. | A review of airport runway capacity evaluation model | |
| Shinohara et al. | Study of lateral attenuation under meteorological conditions for airport noise modeling | |
| CN109829190A (en) | A kind of parameters analysis method of airport asphalt concrete pavement deformation simulative | |
| Pruis et al. | Fast-time Wake Vortex Model Predictions Compared with Observations Behind Landing Aircraft Near the Ground |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190524 |
|
| RJ01 | Rejection of invention patent application after publication |