CN107657097A - The method for numerical simulation of mine road roughness - Google Patents
The method for numerical simulation of mine road roughness Download PDFInfo
- Publication number
- CN107657097A CN107657097A CN201710832754.7A CN201710832754A CN107657097A CN 107657097 A CN107657097 A CN 107657097A CN 201710832754 A CN201710832754 A CN 201710832754A CN 107657097 A CN107657097 A CN 107657097A
- Authority
- CN
- China
- Prior art keywords
- road roughness
- mine
- numerical simulation
- transporting equipment
- mine road
- 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
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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Discrete Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of method for numerical simulation of mine road roughness, including：Determine the upper limit, the lower limit of the intrinsic frequency of transporting equipment, and the common speed scope of transporting equipment；Mine road surface sample length, sampling number are obtained according to intrinsic frequency, common speed；The discrete power spectral density of mine road roughness is obtained according to sample length, sampling number, Uneven road degrees of data is obtained according to power spectral density.The present invention can rapidly and accurately obtain Uneven road degrees of data, and can consider the influence of transporting equipment eigentone and the speed of service.
Description
Technical field
The present invention relates to a kind of mine transportation equipment designing technique, specifically, is related to a kind of number of mine road roughness
It is worth analogy method.
Background technology
Mine road surface is the most important running environment of various mine transportation equipments, and road roughness be it is various transport set
Standby main exciting source.Mine Uneven road degrees of data is accurately grasped for improving transporting equipment smooth-going in product design
Property, analysis parts fatigue properties, improve product life-span have great significance.
Mine Uneven road degrees of data can be obtained by following two methods：
(1) method of experimental test；
The cycle of this method experimental test is long, cost is high, and the later stage needs substantial amounts of data processing, in addition mine road surface
It is varied, it is impossible to which that each road surface is tested.
(2) given mine road roughness power spectral density is transformed to road roughness.
There are trigonometric series method, wave filter shaped white noise method, integration unit white noise method etc. at present, obtained by these methods
Road roughness power spectral density there is certain error compared with given power spectral density.
And above two method does not consider the influence of transporting equipment eigentone and the speed of service.
The content of the invention
Technical problem solved by the invention is to provide a kind of method for numerical simulation of mine road roughness, can be quick
Uneven road degrees of data is obtained exactly, and can consider the shadow of transporting equipment eigentone and the speed of service
Ring.
Technical scheme is as follows：
A kind of method for numerical simulation of mine road roughness, including：
Determine the upper limit, the lower limit of the intrinsic frequency of transporting equipment, and the common speed scope of transporting equipment；
Mine road surface sample length, sampling number are obtained according to intrinsic frequency, common speed；
The discrete power spectral density of mine road roughness is obtained according to sample length, sampling number, according to power spectrum
Density obtains Uneven road degrees of data.
Further：The upper limit n of spatial frequency1=f1/v1, wherein f1It is the upper limit of transporting equipment intrinsic frequency, v1It is conventional
The lower limit of speed；The lower limit n of spatial frequency2=f2/v2, wherein f2It is the lower limit of transporting equipment intrinsic frequency, v2It is common speed
The upper limit.
Further：Sampling number N=L/ Δsl, distance samples interval deltal≤1/2n1, sample length L >=1/n2。
Further：Power spectral density Gq(nk)=Gq(n0)(n/n0)-w, wherein Gq(n0) it is irregularity coefficient, n is space frequency
Rate, n0Value is 0.1m-1, frequency index W=2.
Further：According to formulaRoad roughness discrete Fourier transform modulus value can be obtained；
According to formulaRoad roughness discrete Fourier transform plural form can be obtained；According to formulaInverse Fourier transform obtains Uneven road degrees of data.
Compared with prior art, the technology of the present invention effect includes：
The present invention can rapidly and accurately obtain Uneven road degrees of data, and can consider transporting equipment and inherently shake
The influence of dynamic frequency and the speed of service.
Relative to traditional testing method or method for numerical simulation, the characteristics of having quick and precisely of the invention, and
The influence of transporting equipment eigentone and the speed of service is above all considered, for accurately being commented at the beginning of design
Estimate equipment ride comfort, analyze parts fatigue properties, improve product life-span have great significance.
The present invention has considered two factors of transporting equipment intrinsic frequency and travelling speed, using Fourier Tranform and Fu
The method of sharp leaf inverse transformation obtains the uneven degrees of data on road surface, and can be realized by software programming.
Brief description of the drawings
Fig. 1 is the flow chart of the method for numerical simulation of chats mountain of the present invention road roughness；
Fig. 2 is road roughness curve map in the present invention.
Embodiment
Technical solution of the present invention is elaborated below with reference to example embodiment.However, example embodiment can
Implement in a variety of forms, and be not understood as limited to embodiment set forth herein；On the contrary, these embodiments are provided so that
The design of example embodiment more comprehensively and completely, and is comprehensively communicated to those skilled in the art by the present invention.
As shown in figure 1, it is the flow chart of the method for numerical simulation of chats mountain of the present invention road roughness.
The method for numerical simulation of mine road roughness, the shadow of transporting equipment intrinsic frequency and travelling speed is considered
Ring, specifically include following steps：
In order to reach above-mentioned technical purpose, technical scheme is as follows：
Step 1：Determine the upper limit, the lower limit of the intrinsic frequency of transporting equipment, and transporting equipment common speed scope；
Step 2：Determine the upper limit, the lower limit of spatial frequency；
Spatial frequency upper limit n1=f1/v1, wherein f1It is the upper limit of transporting equipment intrinsic frequency, v1It is under common speed
Limit；The lower limit n of spatial frequency2=f2/v2, wherein f2It is the lower limit of transporting equipment intrinsic frequency, v2It is the upper limit of common speed.
Step 3：The sample length on mine road surface is determined, when calculating Road Surface Power Spectrum Density to avoid frequency alias, away from
From sampling interval Δl≤1/2n1, sample length L >=1/n2。
Step 4：Determine the sampling number on mine road surface, sampling number N=L/ Δsl。
Step 5：Calculate the discrete power spectral density of mine road roughness；
According to requirements of the national standard, Road Surface Power Spectrum Density Gq(n)=Gq(n0)(n/n0)-w, wherein Gq(n0) it is unevenness
Coefficient, different values is taken according to the difference of pavement grade；N is spatial frequency, and it is the inverse of wavelength X, is represented in every meter of length
Periodicity comprising ripple；n0It is reference frequency, value 0.1m-1；W is frequency index, and it determines Road Surface Power Spectrum Density
Structure, take W=2.Gq(n0)、n、n0, W be, it is known that can thus draw discrete Road Surface Power Spectrum Density Gq(nk)。
Step 6：According to formulaRoad roughness discrete Fourier transform mould can be obtained
Value；
Step 7：According to formulaRoad roughness discrete Fourier transform plural number shape can be obtained
Formula；
Step 8：According to formulaInverse Fourier transform obtains Uneven road degrees of data.
Embodiment 1
1) by taking the prevailing traffic equipment quarry tipper of surface mine as an example, its intrinsic frequency scope is 1-30Hz, often
It is 10km/h~35km/h (2.8m/s~9.7m/s) with velocity interval.
2) according to formula n=f/v, it may be determined that the upper limit of spatial frequency is 10.7m-1, lower limit 0.1m-1, illustrate only to need
In 0.1m-1~10.7m-1In the range of simulated roadway unevenness, can to cover the intrinsic frequency of quarry tipper vibrational system,
When thereby may be ensured that the performance for utilizing obtained road roughness to study quarry tipper, tally with the actual situation.
3) according to formula：Distance samples interval delta l≤1/2n1=1/ (2*10.7)=0.047m, sample length L >=1/n2
=1/0.1=10m takes Δ l=0.01, L=256m in this example.
4) according to formula N=L/ Δ l, it is N=256/0.01=25600 to determine sampling number.
5) further, according to formula Gq(n)=Gq(n0)(n/n0)-w, discrete Road Surface Power Spectrum Density G can be obtainedq
(n0)、Gq(n1)、Gq(n2)、Gq(n3)…….G in this exampleq(n0) irregularity coefficient value is 256mm2/m-1, equivalent to country
The requirement on C levels road surface in standard.
6) according to formulaRoad roughness discrete Fourier transform modulus value can be obtained | X0|、|
X1|、|X2|、|X3|……。
7) according to formulaRoad roughness discrete Fourier transform plural form X can be obtained0、X1、
X2、X3……。
8) according to formulaUneven road degrees of data x is obtained to obtain inverse Fourier transform0、
x1、x2、x3……。
Road roughness curve can be drawn by being drawn using MATLAB softwares, as shown in Fig. 2 being the Uneven road of example 1
Write music line.
Term used herein is explanation and exemplary and nonrestrictive term.Because the present invention can be with a variety of
Form specific implementation without departing from the spiritual or substantive of invention, it should therefore be appreciated that above-described embodiment be not limited to it is any foregoing
Details, and should widely being explained in the spirit and scope that appended claims are limited, thus fall into claim or its etc.
Whole changes and remodeling in the range of effect all should be appended claims and covered.
Claims (5)
1. a kind of method for numerical simulation of mine road roughness, including：
Determine the upper limit, the lower limit of the intrinsic frequency of transporting equipment, and the common speed scope of transporting equipment；
Mine road surface sample length, sampling number are obtained according to intrinsic frequency, common speed；
The discrete power spectral density of mine road roughness is obtained according to sample length, sampling number, according to power spectral density
Obtain Uneven road degrees of data.
2. the method for numerical simulation of mine road roughness as claimed in claim 1, it is characterised in that：The upper limit n of spatial frequency1
=f1/v1, wherein f1It is the upper limit of transporting equipment intrinsic frequency, v1It is the lower limit of common speed；The lower limit n of spatial frequency2=f2/
v2, wherein f2It is the lower limit of transporting equipment intrinsic frequency, v2It is the upper limit of common speed.
3. the method for numerical simulation of mine road roughness as claimed in claim 2, it is characterised in that：Sampling number N=L/ Δsl,
Distance samples interval deltal≤1/2n1, sample length L >=1/n2。
4. the method for numerical simulation of mine road roughness as claimed in claim 3, it is characterised in that：Power spectral density Gq(nk)
=Gq(n0)(n/n0)-w, wherein Gq(n0) it is irregularity coefficient, n is spatial frequency, n0Value is 0.1m-1, frequency index W=2.
5. the method for numerical simulation of mine road roughness as claimed in claim 3, it is characterised in that：According to formulaRoad roughness discrete Fourier transform modulus value can be obtained；According to formula
Road roughness discrete Fourier transform plural form can be obtained；According to formulaFourier inversion
Get Uneven road degrees of data in return.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710832754.7A CN107657097A (en) | 2017-09-15 | 2017-09-15 | The method for numerical simulation of mine road roughness |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710832754.7A CN107657097A (en) | 2017-09-15 | 2017-09-15 | The method for numerical simulation of mine road roughness |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107657097A true CN107657097A (en) | 2018-02-02 |
Family
ID=61129858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710832754.7A Pending CN107657097A (en) | 2017-09-15 | 2017-09-15 | The method for numerical simulation of mine road roughness |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107657097A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108804814A (en) * | 2018-06-08 | 2018-11-13 | 吉林大学 | A method of extraction spectrum of road surface roughness parameter |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103318180A (en) * | 2013-05-28 | 2013-09-25 | 万向钱潮股份有限公司 | Vehicular road surface irregularity automatic identification system and method |
| CN105427375A (en) * | 2015-12-16 | 2016-03-23 | 北京汽车研究总院有限公司 | Method and system for generating three-dimensional road profile |
-
2017
- 2017-09-15 CN CN201710832754.7A patent/CN107657097A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103318180A (en) * | 2013-05-28 | 2013-09-25 | 万向钱潮股份有限公司 | Vehicular road surface irregularity automatic identification system and method |
| CN105427375A (en) * | 2015-12-16 | 2016-03-23 | 北京汽车研究总院有限公司 | Method and system for generating three-dimensional road profile |
Non-Patent Citations (1)
| Title |
|---|
| 任思润 等: "基于功率谱密度的路面不平整度数值模拟", 《黑龙江科技信息》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108804814A (en) * | 2018-06-08 | 2018-11-13 | 吉林大学 | A method of extraction spectrum of road surface roughness parameter |
| CN108804814B (en) * | 2018-06-08 | 2020-05-15 | 吉林大学 | Method for extracting road surface spectrum parameters |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Acoustic emission detection of rail defect based on wavelet transform and Shannon entropy | |
| CN106770015A (en) | A kind of oil property detection method based on the similar differentiation of principal component analysis | |
| Liu et al. | Construction of road roughness in left and right wheel paths based on PSD and coherence function | |
| CN107657097A (en) | The method for numerical simulation of mine road roughness | |
| CN103527172B (en) | Variable rock electric coupling index water saturation computational methods | |
| Liu et al. | An Introduction of a Robust OMA Method: CoS‐SSI and Its Performance Evaluation through the Simulation and a Case Study | |
| Li et al. | Influence of road surface characteristics on tire–road noise for thin-layer surfacings | |
| Lai et al. | Noise reduction characteristics of macroporous asphalt pavement based on a weighted sound pressure level sensor | |
| KR20160113231A (en) | Roadway infrastructure monitoring based on aggregated mobile vehicle communication parameters | |
| Li et al. | International roughness index and a new solution for its calculation | |
| Sun et al. | Nonstationary dynamic pavement loads generated by vehicles traveling at varying speed | |
| CN103926215A (en) | Method for quickly detecting asphalt penetration | |
| CN102095399B (en) | Road spectrum based method for extracting quantified road characteristic parameters | |
| CN111221038B (en) | Method and device for quantitatively predicting thickness of thin reservoir | |
| Moeinikia et al. | Experimental Design in reservoir simulation: an integrated solution for uncertainty analysis, a case study | |
| Adu-Gyamfi et al. | Critical analysis of different Hilbert-Huang algorithms for pavement profile evaluation | |
| CN110569608B (en) | Method and device for determining vibration severity of vehicle-mounted platform | |
| CN103207932A (en) | Method for determining automobile wheelbases to conduct load analysis | |
| Zhang et al. | Updating fatigue damage coefficient in railway bridge design code in China | |
| Praticò et al. | On the dependence of acoustic pore shape factors on porous asphalt volumetrics | |
| Wu et al. | Fast calibration for vibration-based pavement roughness measurement based on model updating of vehicle dynamics | |
| Ren et al. | Piecewise polynomial fitting with trend item removal and its application in a cab vibration test | |
| Neunzert | Progress in Industrial Mathematics at ECMI 94: Progress in Industrial Mathematics at ECMI 94 | |
| Lin | Development of interharmonics identification using enhanced‐fft algorithm | |
| CN206593873U (en) | Acoustic experiment rotary drum and acoustic experiment device |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180202 |