CN110108354A - A method of passing through noise in small size semianechoic room interior prediction vehicle - Google Patents
A method of passing through noise in small size semianechoic room interior prediction vehicle Download PDFInfo
- Publication number
- CN110108354A CN110108354A CN201910474575.XA CN201910474575A CN110108354A CN 110108354 A CN110108354 A CN 110108354A CN 201910474575 A CN201910474575 A CN 201910474575A CN 110108354 A CN110108354 A CN 110108354A
- Authority
- CN
- China
- Prior art keywords
- point
- equivalent source
- noise
- near field
- vehicle
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M17/00—Testing of vehicles
- G01M17/007—Wheeled or endless-tracked vehicles
Abstract
The present invention provides a kind of method for passing through noise in small size semianechoic room interior prediction vehicle, this method is based on interpolation time-domain equivalent source method and is realized, it include: to establish coordinate system, equivalent source point and its mirror point position, the arrangement near field measurement point and far field microphone array target point of sound source are set, calculate the equivalent source point and the mirror point to the near field measurement point and the far field microphone array target point distance, and the acoustic pressure measured in conjunction near field measurement point, the source strength of each equivalent source is calculated, the acoustic pressure of far field microphone array target point is rebuild.The present invention reduces requirement when actual vehicle passes through noise testing to anechoic room bulk, realizes that vehicle is carried out in the semianechoic room of small size quickly and accurately to be predicted by noise, and can reduce the investment of measurement process.
Description
Technical field
The present invention relates to noise testing technical fields, logical in small size semianechoic room interior prediction vehicle more particularly to one kind
The method for crossing noise.
Background technique
With the development and continuous improvement of people's living standards of society, requirement of the people to living environment is also constantly
It is promoted, various noise problems are increasingly taken seriously, such as road traffic noise problem, it seriously affects the daily life of people
The problems such as living, causing life worried and sleep disturbance to people.In order to correctly assess vehicle in urban road making an uproar when driving
Sound is horizontal, and International Organization for standardization enacts and enforces ISO 362-1:2007, which is tested using acceleration noise and at the uniform velocity made an uproar
Sound test combines and what the mode that weights assessed vehicle passes through noise level.China is by testing New Draft outside noise room
(GB1495-20xx) it is currently in the stage of consulting on, it is contemplated that the year two thousand twenty is formally issued and implemented.
Automobile accelerates through noise testing international standard (ISO 362-1:2007) in measuring instrument, measurement place, measurement
Environment and measurement vehicle prepare 4 aspects and propose tightened up requirement, and vehicle is increased by noise room external pelivimetry difficulty, measures number
It is reduced according to qualification rate, it is time-consuming.And vehicle is by interference of the noise testing vulnerable to the factors such as climatic environment and driving technology, it is real
It is poor to test reproducibility of results.For economical and efficient implement to measure, it is logical that vehicle is carried out in domestic and foreign scholars' research in semianechoic room
Cross the method and technology of noise prediction.1999, Fry et al. proposed that carrying out vehicle based on neural network numerical analysis method passes through
Noise prediction, 2001, Kim et al. passed letter by near field measurement and sound and carries out through noise prediction and noise source contribution amount point
Analysis, 2006, Genuit et al. using acoustic array technology carries out in anechoic room vehicle pass through noise testing experimental study, 2013
Year, Belgian Univ Louvain is based on quick multistage Element BEM (FMBEM) and assesses vehicle by noise, together by CAD model
When consider vehicle movement caused by Doppler effect influence, Walsh etc. summarize presently, there are interior pass through noise testing side
Method and its advantage and disadvantage study the characteristic of each noise source of vehicle.
These methods can avoid vehicle from passing through the outdoor test of noise, realize to vehicle through the prediction of noise, and
And the characteristics of respectively having oneself, but there is also certain limitations.For example, the base that international standard (ISO 362-3:2016) is advocated
It needs to carry out in large-sized semianechoic room of 30m × 20m by measurement method in noise room in the vehicle of acoustic array technology,
And experiment at least needs 42 microphones, experimental cost is too high.Vehicle based on quick multistage Element BEM (FMBEM) passes through
Although noise evaluation method can be under conditions of not carrying out experiment test, realization is carried out in early stage vehicle design phase by making an uproar
Sound prediction, but the calculating of its numerical value needs to expend long time, and precision of prediction is relatively low.
Summary of the invention
In view of the above deficiencies of the prior art, the purpose of the present invention is to provide one kind in small size semianechoic room
Method of the vehicle by noise is predicted, to can fast and effeciently predict vehicle by making an uproar in the semianechoic room of small size
Sound replaces actual noise source using equivalent source point, by arranging near field measurement point, reduces the requirement to measurement bulk；It is logical
The source strength of each equivalent source of acoustic pressure reverse measured near field measurement point is crossed, to realize to far field microphone array target point sound field
It is rebuild.The invention, which can also be realized, passes through each noise source contribution amount of noise and its noise source specificity analysis to vehicle.
To achieve the above object, the present invention provides a kind of in the side that small size semianechoic room interior prediction vehicle passes through noise
Method, it is described it is a kind of in small size semianechoic room interior prediction vehicle by the method for noise the following steps are included:
A) three-dimensional cartesian coordinate system is defined, the three-dimensional cartesian coordinate system is using room floor of eliminating the noise as X-Y axial plane；
B) the equivalent source point in several vehicle noise sources is established, and arranges several near field measurement points around vehicle, is recorded
The coordinate of the equivalent source point and near field measurement point on the coordinate system, and the near field is measured using measuring device
The acoustic pressure of measurement point；
C) mirror point of the equivalent source point based on noise elimination room floor is established, and records the mirror point in the coordinate
Coordinate in system；
D) position that far field microphone array target point is defined in the vehicle two sides, records the far field microphone array
Coordinate of the target point in the coordinate system；
E) it according to the coordinate of the equivalent source point, the mirror point and the far field microphone array target point, calculates
The distance of the equivalent source point and the mirror point to the far field microphone array target point；
F) according to the coordinate of the equivalent source point, the mirror point and near field measurement point, the equivalent source is calculated
The distance of point and the mirror point to the near field measurement point；
G) according to the distance and the near field measurement of the equivalent source point and the mirror point to the near field measurement point
Point acoustic pressure calculates each equivalent source source strength；
H) according to each equivalent source source strength and the equivalent source point and the mirror point to the far field microphone array
The distance of column target point calculates the acoustic pressure of the far field microphone array target point.
Optionally, the quantity of far field microphone array target point described in each side of vehicle is between 18-24.
Optionally, the step f) includes step f1) it is surveyed according to the equivalent source point and the mirror point to the near field
The distance of amount point rebuilds semi-free field space acoustic pressure transmission function, and the acoustic pressure transmission function is
Wherein, δ () is Dirac function, and c is the velocity of sound, r (s) be the equivalent source point to the near field measurement point away from
It is distance of the mirror point to the near field measurement point from r'(s), t is indicated when the near field measurement point measures acoustic pressure
It carves.
Optionally, near field measurement point acoustic pressure described in the step g) is each equivalent source source strength and acoustic pressure transmission function
Convolution superposition value.
Optionally, the source of each equivalent source is calculated in the step g) by cubic spline interpolation time-domain equivalent source method
By force.
Optionally, the step g) includes step g1) calculate each equivalent source delay time T, τ=ti-rHmn/ c,
In, c is the velocity of sound, tiThe specific moment of acoustic pressure, r are measured for the near field measurement pointHmnIt is the near field measurement point to described equivalent
The distance of source point.
Optionally, the step g) includes step g2) the calculating mirror point delay time T ', τ '=ti-r'Hmn/ c,
In, c is the velocity of sound, tiThe specific moment of acoustic pressure, r' are measured for the near field measurement pointHmnFor the near field measurement point to the mirror
The distance of picture point.
Optionally, all far field microphone array target points are in tiThe acoustic pressure at specific moment and all equivalent source points exist
τjTransfer matrix between the source strength at moment is
Wherein, τjIndicate discrete, the r to the equivalent source point delay timeS11-rSLNIndicate that any equivalent source point arrives
The distance of the far field microphone array target point, r'S11-r'SLNIndicate the mirror point to the far field microphone array mesh
The distance of punctuate, φj(τ) is cubic spline functions,Indicate that any equivalent source point corresponds to tiThe specific moment
Delay time,Indicate that any mirror point corresponds to tiThe delay time at specific moment.
As described above, compared with the prior art, a kind of pass through in small size semianechoic room interior prediction vehicle of the invention is made an uproar
The method of sound, has the advantages that, such as:
Using the present invention, it can be achieved that carrying out avoiding a large amount of vehicle by noise prediction for vehicle to be measured in anechoic room
It is tested by noise room external pelivimetry, makes the interference of measurement not climate environmental factor；When noise elimination chamber size is unsatisfactory for iso standard
It is required that size when, this method may be implemented to be effectively reduced measurement by noise in half noise elimination interior prediction vehicle of small size
Requirement to anechoic room bulk reduces the cost of noise testing；The present invention is by introducing cubic spline interpolation time-domain equivalent
The method in source avoids carrying out the transmission function in acoustic propagation path experiment measurement or Numerical evaluation, realizes that vehicle is each by noise
The assessment of the size of noise source contribution amount.
Detailed description of the invention
Fig. 1 is shown as the step flow chart in small size semianechoic room interior prediction vehicle by the method for noise.
Fig. 2 is shown as the coordinate system schematic diagram established.
Fig. 3 is shown as microphone array target point theory sound pressure level in far field under at the uniform velocity operating condition and compares with sound pressure level result is calculated
Figure.
Fig. 4 be shown as under at the uniform velocity operating condition the theoretical third-octave at far field microphone array target point peak sound pressure moment with
Calculate third-octave result comparison chart.
Fig. 5 is shown as microphone array target point theory sound pressure level in far field under accelerating mode and compares with sound pressure level result is calculated
Figure.
Fig. 6 be shown as under accelerating mode the theoretical third-octave at far field microphone array target point peak sound pressure moment with
Calculate third-octave result comparison chart.
Specific embodiment
Illustrate embodiments of the present invention below by way of specific specific example, those skilled in the art can be by this specification
Other advantages and efficacy of the present invention can be easily understood for disclosed content.The present invention can also pass through in addition different specific realities
The mode of applying is embodied or practiced, the various details in this specification can also based on different viewpoints and application, without departing from
Various modifications or alterations are carried out under spirit of the invention.Meanwhile in this specification it is cited as "upper", "lower", " preceding ",
The terms such as " rear ", "left", "right", " centre " are merely convenient to narration and are illustrated, rather than to limit the enforceable model of the present invention
It encloses, relativeness is altered or modified, under the content of no substantial changes in technology, when being also considered as the enforceable scope of the present invention.
Please refer to Fig. 1-6.It should be noted that only the invention is illustrated in a schematic way for diagram provided in the present embodiment
Basic conception, only shown in schema then with related component in the present invention rather than component count, shape when according to actual implementation
Shape and size are drawn, when actual implementation kenel, quantity and the ratio of each component can arbitrarily change for one kind, and its component cloth
Office's kenel may also be increasingly complex.
The present invention provides a kind of method in small size semianechoic room interior prediction vehicle by noise, referring to Fig. 1, described
In small size semianechoic room interior prediction vehicle by the method for noise the following steps are included:
Step S1 is executed, defines three-dimensional cartesian coordinate system, the three-dimensional cartesian coordinate system is using room floor of eliminating the noise as X-Y axis
Plane.
As an example, referring to Fig. 2, the vehicle is placed in the semianechoic room of small size, with noise elimination room floor
Reflecting surface i.e. as shown in Figure 2 defines three-dimensional cartesian coordinate system as X-Y axial plane, with perpendicular to semianechoic room ground simultaneously
The direction definition for being directed toward vehicle up direction is Z axis positive direction, to be parallel to semianechoic room ground and be directed toward vehicle front as X-axis pros
To establishing rectangular coordinate system to be directed toward copilot side such as right direction as Y-axis positive direction.
Step S2 is executed, establishes the equivalent source point in several vehicle noise sources, and arrange that several near fields are surveyed around vehicle
Point is measured, records the coordinate of the equivalent source point and near field measurement point on the coordinate system, and survey using measuring device
Measure the acoustic pressure at the near field measurement point.
It should be noted that the progress vehicle proposed by the present invention in the semianechoic room of smaller size passes through noise prediction
Method by arrangement near field measurement point, realize in the semianechoic room interior prediction vehicle of small size through noise, it is logical to reduce vehicle
Requirement of the noise testing to anechoic room bulk is crossed, while this method is also applied for other large scale dead room environments.
In one embodiment, by taking vehicle in " traveling " as an example, to the Main Noise Sources in its driving process, such as air inlet
N number of equivalent source point is arranged in the vehicle noises such as mouth, engine, gearbox, exhaust system source position, and at vehicle peripheral disposition M
Near field measurement point, wherein M >=N, the shape of near field measurement point surround vehicle periphery substantially with such as sphere or cuboid, to use up
Vehicle may be obtained in detail passes through noise information.
Step S3 is executed, establishes the mirror point of the equivalent source point based on noise elimination room floor, and record the mirror point
Coordinate in the coordinate system.
Referring to Fig. 2, specifically, having an infinitely great boundary face Z=0, space locating for vehicle being divided into upper half-space
With the lower half-space, the mirror point is in mirror image with the boundary face of the infinity with the equivalent source point.
As an example, recording the coordinate (X of an equivalent source S of the vehicle such as vehicles,Ys,Zs) and near field survey
The coordinate (X, Y, Z) of point is measured, and the coordinate of the mirror point S ' of equivalent source S is (Xs,Ys,-ZS)。
Step S4 is executed, the position of far field microphone array target point is defined in the vehicle two sides, records the far field
Coordinate of the microphone array target point in the coordinate system.
As an example, there be 18-24 microphone in each side of vehicle, they constitute far field microphone array target point, Y direction
Apart from 7.45-7.55 meters of vehicle, 1.18-1.22 meters of far field microphone array target point height, each far field microphone array in X-axis
0.8-1.2 meters of column target point spacing chooses any far field microphone array target point l, and records the far field microphone array mesh
Coordinate (the X of punctuatel,Yl,Zl).In one embodiment, the coordinate of certain far field microphone array target point l be (1,7.5,
1.2)。
It should be noted that the present invention is by being defined on far field microphone array aiming spot by surveying in noise room
Microphone array column position is measured, arrangement far field microphone array is avoided, vehicle is allowed not have to cloth by the indoor measurement of noise
Far field microphone array is set, to realize that quickly and accurately predicts vehicle in smaller size anechoic room passes through noise.
Step S5 is executed, according to the equivalent source point and its mirror point and the seat of the far field microphone array target point
The distance of mark, the calculating equivalent source point and its mirror point to the far field microphone array target point.
Step S6 is executed, according to the equivalent source point and its mirror point and the coordinate of near field measurement point, calculates institute
State equivalent source point and its mirror point to the near field measurement point distance.
As an example, the equivalent source S and its mirror point S ' is denoted as r to the distance between near field measurement point respectivelysWith
rs', wherein
It executes step S6 and executes there are also an execution step S61 between step S7, the acoustic pressure for establishing half space sound field is transmitted
Function
Wherein, in formula (2), δ (x) is Dirac delta function, and c is the velocity of sound, and t is indicated in the near field measurement point measurement sound
At the time of pressure.
Execute step S7, according to the distance of the equivalent source point and its mirror point to the near field measurement point and it is described closely
Field measurement point acoustic pressure calculates each equivalent source source strength.
As an example, utilizing such as time-domain equivalent source thought in executing step S7, then m-th of near field measurement point is in tiWhen
The acoustic pressure at quarter can be superimposed to obtain by the source strength of N number of equivalent source with the convolution of impulse response function.
As an example, executing step S7 includes executing step S71: calculating the equivalent source point delay time T, τ=ti-
rHmn/ c, wherein c is the velocity of sound, tiThe specific moment of acoustic pressure, r are measured for the near field measurement pointHmnIt is arrived for the near field measurement point
The distance of the equivalent source point.
As example processing steps S7 include execute step S72: calculate the mirror point delay time T ', τ '=ti-
r'Hmn/ c, wherein c is the velocity of sound, tiThe specific moment of acoustic pressure, r' are measured for the near field measurement pointHmnFor the near field measurement point
To the distance of the mirror point.According to executing step S7, executing step S71 and executing step S72, i.e.,
In formula (3), qn(ti) indicate n-th of equivalent source in tiThe source strength at moment, rHmnWith r'HmnRespectively m-th of near field is surveyed
Amount o'clock is to the distance between n-th of point source of sound and its mirror point.
Enable τHmn=ti-rHmn/c、τ'Hmn=ti-r'Hmn/ c, then formula (3) can be converted into
Similarly, first of far field microphone array target point is in tiThe acoustic pressure at moment
It can be seen that the source strength being calculated in formula (4), which can not be directly used in formula (5), calculates far field microphone
The sound field of array target point.Therefore, it is the reconstruction for realizing sound field, cubic spline functions is introduced into equation (4) reciprocity
It imitates source node noise intensity and carries out interpolation, before this, discrete obtain need to be carried out to delay time (source signal time) τ
τj=τ0+jΔτ (6)
Wherein, j=1,2, L, J, Δ τ indicate delay time step-length, τ0Indicate the initial delay moment, i.e., as τ≤τ0When, etc.
Imitating source strength is zero.
Interpolation is carried out to source strength equivalent in formula (4) using cubic spline functions, can be obtained
It should be noted that being carried out to vehicle by noise by introducing spline interpolation time-domain equivalent source method
Sound field prediction, is conducive to analyze influence of each noise source to vehicle by noise, such as the shadow of independent analysis vehicle tire noise
It rings, so that measurement is not limited, analysis can be accurately measured for each noise source.
Wherein,Indicate that n-th of equivalent source point postpones the source strength at moment, Φ at j-thj(τ) is cubic spline interpolation letter
Number, expression formula are
Formula (7) substitution formula (4) can be obtained
M near field measurement point can be obtained in t by equation (9)iThe acoustic pressure at specific moment
Wherein:
Therefore, the key of problem is converted to the equivalent source strength for calculating each discrete delay moment, and far field biography can be realized
The acoustic pressure of sound device array target point calculates.It notices
And if only if τj-1≤τi Hmn≤τj+1With τj-1≤τi' Hmn≤τj+1When,Not equal to zero, with formula (11)
First item analysis for, enable
τ1=t1-Rmin/ c Δ t=V τ (12)
Then τj-1≤τi Hmn≤τj+1It can be expressed as
t1-Rmin/c+(j-2)Vτ≤t1-Rmn/c+(i-1)Vτ≤t1-Rmin/c+jVτ (13)
It is available by formula (13) and formula (8):
φj(τi Hmn)=φj-1(τi-1 Hmn)=φj-2(τi-2 Hmn)=L=φ1(τi-j+1 Hmn)
(14)
Similarly, the second part of formula (11) also available identical conclusion, therefore just have
Therefore equation (10) can simplify for
Wherein:
The characteristics of according to interpolating function, if makingIt is not zero, then
τ0≤τi-j+1 Hmn≤τ2 (18)
Wherein:
τ0=τ1-Δτ (19)
It can be obtained by formula (12) and formula (18)
Then have
It enablesThen the nonzero element of G is located at section [1, K], therefore equation (10) can be further
It is expressed as
Therefore, the equivalent source strength of i-th of discrete instants can be expressed as
Step S80 is executed, according to each equivalent source source strength and the equivalent source point and its mirror point to the far field
The distance of microphone array target point calculates the acoustic pressure of the far field microphone array target point.
As an example, by the equivalent source strongly expressed of i-th of discrete instants known to above-mentioned execution step S70, therefore far field
The acoustic pressure of microphone array target point l
Wherein,Indicate all far field microphone array target points in tiThe acoustic pressure at specific moment and all equivalent source points
In τjTransfer matrix between the source strength at moment,
It should be noted that in formula (25), i=1,2, L, floor ((τJ+rSmin/ c)/Δ t), wherein rSminFor
rSmnWith r'SmnIn minimum value, floor (x), which is rounded downwards i.e., to be indicated to take less than or equal to x and the nearest integer of distance x.
As an example, please referring to Fig. 3-4-5-6, the effective of verifying the method for the present invention will be gone from two angles of time-domain and frequency-domain
Property, specifically, in experiment tested with neutral gear (1500rpm) constant rotational speed Work condition analogue vehicle at the uniform velocity operating condition, it is complete with neutral gear
Throttle accelerating mode simulating vehicle accelerating mode is tested, and each near field measurement point, equivalent source point and its mirror point, far field are recorded
The position coordinates of microphone array target point, wherein certain far field microphone array target point l coordinate choose for example (1,7.5,
1.2)。
Fig. 3-4 is please referred to, at the uniform velocity under operating condition, the calculated result and measurement result of far field microphone array target point sound pressure level
Almost the same, the difference of the two is controlled substantially within ± 1dB.And the sound that the method for the present invention is rebuild it can be seen from spectrogram
The third-octave figure and actual measured results for peak value moment of arbitrarily downgrading are almost the same, and the amplitude of each frequency band is essentially identical.
Please refer to Fig. 5-6, under accelerating mode, the calculated result of the sound pressure level of far field microphone array target point and measurement are tied
Fruit is almost the same, and the difference of the two is controlled substantially within ± 1dB.And the method for the present invention reconstruction it can be seen from spectrogram
The third-octave figure and actual measured results of sound pressure level peak value moment are almost the same, and the amplitude of each frequency band is essentially identical.
By above-mentioned test result, it can be seen that, one kind of the present invention passes through in small size semianechoic room interior prediction vehicle
The method of noise, it can be achieved that predict indoors vehicle by noise, avoid a large amount of vehicle by noise room external pelivimetry, reduce
Influence of the climatic environment to measurement result；Meanwhile when noise elimination chamber size is unsatisfactory for the size of iso standard requirement, this method can
To realize in half noise elimination interior prediction vehicle of small size through noise, measurement place problem is efficiently solved, experimental cost is reduced,
Save experimental period.
In conclusion the present invention provides a kind of method in small size semianechoic room interior prediction vehicle by noise, it can be real
Present Indoo r prediction vehicle passes through noise, experimental result favorable repeatability；Meanwhile it solving to be unsatisfactory for ISO when noise elimination chamber size
When measuring size of standard requirements, realization are directed to being predicted by noise for vehicle to be measured in the anechoic room of smaller size,
Cost of the vehicle by noise testing is reduced, a large amount of time of measuring is saved；It, can be more to be measured simultaneously using the present invention
The size of each Main Noise Sources relative contribution amount of vehicle, the influence that analysis each noise source of vehicle passes through noise to vehicle.Therefore, originally
Invention effectively overcomes many disadvantages in the prior art and has high industrial utilization value, is greatly expanded applicable model
It encloses.
The principle of the present invention and its effect is only illustrated in examples detailed above, and is not intended to limit the present invention.It is any to be familiar with
The personage of this technology all without departing from the spirit and scope of the present invention, carries out modifications and changes to above-described embodiment.Therefore,
Such as those of ordinary skill in the art is completed without departing from the spirit and technical ideas disclosed in the present invention
All equivalent modifications or change, should be covered by the claims of the present invention.
Claims (8)
1. a kind of method for passing through noise in small size semianechoic room interior prediction vehicle, which comprises the following steps:
A) three-dimensional cartesian coordinate system is defined, the three-dimensional cartesian coordinate system is using room floor of eliminating the noise as X-Y axial plane；
B) it establishes the equivalent source point in several vehicle noise sources, and arranges several near field measurement points around vehicle, described in record
The coordinate of equivalent source point and near field measurement point on the coordinate system, and the near field measurement is measured using measuring device
Acoustic pressure at point position；
C) mirror point of the equivalent source point based on noise elimination room floor is established, and records the mirror point in the coordinate system
Coordinate；
D) position that far field microphone array target point is defined in the vehicle two sides, records the far field microphone array target
Coordinate of the point in the coordinate system；
E) according to the coordinate of the equivalent source point, the mirror point and the far field microphone array target point, described in calculating
The distance of equivalent source point and the mirror point to the far field microphone array target point；
F) according to the coordinate of the equivalent source point, the mirror point and near field measurement point, calculate the equivalent source point and
Distance of the mirror point to the near field measurement point；
G) according to the distance and near field measurement point sound of the equivalent source point and the mirror point to the near field measurement point
Pressure calculates each equivalent source source strength；
H) according to each equivalent source source strength and the equivalent source point and the mirror point to the far field microphone array mesh
The distance of punctuate calculates the acoustic pressure of the far field microphone array target point.
2. a kind of method in small size semianechoic room interior prediction vehicle by noise according to claim 1, feature
It is, the quantity of far field microphone array target point is between 18-24 described in each side of vehicle.
3. a kind of method in small size semianechoic room interior prediction vehicle by noise according to claim 1, feature
It is, the step f) includes step f1) according to the distance of the equivalent source point and the mirror point to the near field measurement point
Semi-free field space acoustic pressure transmission function is established, the acoustic pressure transmission function is
Wherein, δ () is Dirac function, and c is the velocity of sound, and r (s) is distance of the equivalent source point to the near field measurement point, r'
It (s) is the distance of the mirror point to the near field measurement point, t is indicated at the time of the near field measurement point measures acoustic pressure.
4. a kind of method in small size semianechoic room interior prediction vehicle by noise according to claim 3, feature
It is, near field measurement point acoustic pressure described in the step g) is that the equivalent source point source strength is superimposed with the convolution of acoustic pressure transmission function
Value.
5. a kind of method in small size semianechoic room interior prediction vehicle by noise according to claim 1, feature
It is, the step g) calculates the source strength of each equivalent source by spline interpolation time-domain equivalent source method.
6. a kind of method in small size semianechoic room interior prediction vehicle by noise according to claim 1, feature
It is, the step g) includes step g1) calculate the equivalent source point delay time T, τ=ti-rHmn/ c, wherein c is the velocity of sound,
tiThe specific moment of acoustic pressure, r are measured for the near field measurement pointHmnFor the distance of the near field measurement point to the equivalent source point.
7. a kind of method in small size semianechoic room interior prediction vehicle by noise according to claim 1, feature
It is, the step g) includes step g2) the calculating mirror point delay time T ', τ '=ti-r'Hmn/ c, wherein c is sound
Speed, tiThe specific moment of acoustic pressure, r' are measured for the near field measurement pointHmnFor the near field measurement point to the mirror point away from
From.
8. a kind of described in -7 according to claim 1 pass through the method for noise, spy in small size semianechoic room interior prediction vehicle
Sign is that all far field microphone array target points are in tiThe acoustic pressure at specific moment and all equivalent source points are in τjMoment
Transfer matrix between source strength is
Wherein, τjIndicate discrete, the r to the equivalent source point delay timeS11-rSLNIndicate any equivalent source point described in
The distance of far field microphone array target point, r 'S11-r‘SLNIndicate the mirror point to the far field microphone array target point
Distance, φj(τ) is cubic spline functions,Indicate that any equivalent source point corresponds to tiThe specific moment prolongs
The slow time,Indicate that any mirror point corresponds to tiThe delay time at specific moment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910474575.XA CN110108354B (en) | 2019-06-03 | 2019-06-03 | Method for predicting vehicle passing noise in small-size semi-anechoic chamber |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910474575.XA CN110108354B (en) | 2019-06-03 | 2019-06-03 | Method for predicting vehicle passing noise in small-size semi-anechoic chamber |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110108354A true CN110108354A (en) | 2019-08-09 |
| CN110108354B CN110108354B (en) | 2021-06-04 |
Family
ID=67493560
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910474575.XA Active CN110108354B (en) | 2019-06-03 | 2019-06-03 | Method for predicting vehicle passing noise in small-size semi-anechoic chamber |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110108354B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110907192A (en) * | 2019-11-18 | 2020-03-24 | 东风商用车有限公司 | Method for reducing passing noise of whole vehicle |
| CN112254978A (en) * | 2020-09-09 | 2021-01-22 | 交通运输部公路科学研究所 | Method for measuring and evaluating passing noise source intensity and directivity of vehicle |
| CN113567156A (en) * | 2021-07-22 | 2021-10-29 | 青岛双星轮胎工业有限公司 | Device and method for testing tire noise by using semi-anechoic chamber flat steel belt testing machine |
| CN114936415A (en) * | 2022-04-27 | 2022-08-23 | 浙江大学 | Propeller tone frequency prediction method |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007133484A (en) * | 2005-11-08 | 2007-05-31 | Mitsubishi Electric Corp | Generation noise simulated measurement device and method for electronic substrate |
| CN102790931A (en) * | 2011-05-20 | 2012-11-21 | 中国科学院声学研究所 | Distance sense synthetic method in three-dimensional sound field synthesis |
| CN105929253A (en) * | 2016-06-15 | 2016-09-07 | 北京航空航天大学 | Method using dipole moment model to inverse PCB circuit electromagnetic field |
| CN105975713A (en) * | 2016-05-20 | 2016-09-28 | 北京航空航天大学 | Method for predicting electromagnetic radiation of electronic circuit through source reconstruction |
| CN106777844A (en) * | 2017-03-17 | 2017-05-31 | 合肥工业大学 | A kind of Forecasting Methodology of tyre noise loudspeaker enlarge-effect |
| CN107566970A (en) * | 2017-07-20 | 2018-01-09 | 西北工业大学 | A kind of medium-high frequency Reconstruction of Sound Field method inside enclosed environment |
| CN108627240A (en) * | 2018-03-23 | 2018-10-09 | 南京大学 | A kind of evaluation method of constant volume speed sound source low frequency acoustical power |
| CN109030029A (en) * | 2018-08-09 | 2018-12-18 | 合肥工业大学 | Pass through the controller for automobile throttle pedal and method of noise testing for vehicle |
-
2019
- 2019-06-03 CN CN201910474575.XA patent/CN110108354B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007133484A (en) * | 2005-11-08 | 2007-05-31 | Mitsubishi Electric Corp | Generation noise simulated measurement device and method for electronic substrate |
| CN102790931A (en) * | 2011-05-20 | 2012-11-21 | 中国科学院声学研究所 | Distance sense synthetic method in three-dimensional sound field synthesis |
| CN105975713A (en) * | 2016-05-20 | 2016-09-28 | 北京航空航天大学 | Method for predicting electromagnetic radiation of electronic circuit through source reconstruction |
| CN105929253A (en) * | 2016-06-15 | 2016-09-07 | 北京航空航天大学 | Method using dipole moment model to inverse PCB circuit electromagnetic field |
| CN106777844A (en) * | 2017-03-17 | 2017-05-31 | 合肥工业大学 | A kind of Forecasting Methodology of tyre noise loudspeaker enlarge-effect |
| CN107566970A (en) * | 2017-07-20 | 2018-01-09 | 西北工业大学 | A kind of medium-high frequency Reconstruction of Sound Field method inside enclosed environment |
| CN108627240A (en) * | 2018-03-23 | 2018-10-09 | 南京大学 | A kind of evaluation method of constant volume speed sound source low frequency acoustical power |
| CN109030029A (en) * | 2018-08-09 | 2018-12-18 | 合肥工业大学 | Pass through the controller for automobile throttle pedal and method of noise testing for vehicle |
Non-Patent Citations (5)
| Title |
|---|
| 庞川博 等: ""考虑地面影响的半空间气动噪声数值计算方法"", 《航空学报》 * |
| 张洋: ""非稳态声场预测与重建的时域边界元法及其非稳定性研究"", 《中国优秀硕士学位论文全文数据库》 * |
| 王德威 等: ""高速列车车外噪声预测建模与声源贡献量分析"", 《中南大学学报》 * |
| 邓江华 等: ""基于近场声全息的车内噪声贡献量"", 《吉林大学学报》 * |
| 陈剑 等: ""汽车加速度通过噪声室内预测方法研究"", 《噪声与振动控制》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110907192A (en) * | 2019-11-18 | 2020-03-24 | 东风商用车有限公司 | Method for reducing passing noise of whole vehicle |
| CN112254978A (en) * | 2020-09-09 | 2021-01-22 | 交通运输部公路科学研究所 | Method for measuring and evaluating passing noise source intensity and directivity of vehicle |
| CN112254978B (en) * | 2020-09-09 | 2022-05-17 | 交通运输部公路科学研究所 | Method for measuring and evaluating passing noise source intensity and directivity of vehicle |
| CN113567156A (en) * | 2021-07-22 | 2021-10-29 | 青岛双星轮胎工业有限公司 | Device and method for testing tire noise by using semi-anechoic chamber flat steel belt testing machine |
| CN114936415A (en) * | 2022-04-27 | 2022-08-23 | 浙江大学 | Propeller tone frequency prediction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110108354B (en) | 2021-06-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110108354A (en) | A method of passing through noise in small size semianechoic room interior prediction vehicle | |
| Siltanen et al. | Rays or waves? Understanding the strengths and weaknesses of computational room acoustics modeling techniques | |
| CN109813501B (en) | Method, device and system for measuring leakage position of gas pipeline | |
| CN104574484B (en) | A kind of method and apparatus producing picture dynamic effect based on interactive operation | |
| CN106680376A (en) | Sound insulation measurement system and measurement method based on three-dimensional sound intensity array | |
| CN101957442A (en) | Sound source positioning device | |
| CN111159928A (en) | Transformer noise calculation method and system based on multi-line sound source model | |
| JPH04148840A (en) | Method for simulating tire noise and simulator | |
| CN110049408A (en) | A kind of microphone speaker array formation optimization method | |
| CN104407328A (en) | Method and system for positioning sound source in enclosed space based on spatial pulse response matching | |
| CN104678359B (en) | A kind of porous sound holographic method of sound field identification | |
| CN103091395A (en) | Method and device for three-dimensionally dynamically detecting fluid phase state in porous medium under reservoir conditions | |
| CN107566969A (en) | A kind of enclosed environment internal low-frequency Reconstruction of Sound Field method | |
| EP1682856A1 (en) | A method of determining the sound pressure resulting from a surface element of a sound emitting surface | |
| CN109100685A (en) | A kind of passive acoustic direction blending algorithm of two-sided quaternary cross battle array | |
| CN108627801A (en) | Movable body position estimating system, device and method | |
| CN103712652B (en) | Multi-physical field measuring method and device based on acoustic sensing technology | |
| CN109307608B (en) | Method for plugging and repairing train sound leakage part in whole vehicle semi-anechoic chamber environment | |
| CN206563731U (en) | A kind of sound insulation measurement system arranged based on three dimensional sound powerful formation | |
| CN108154557B (en) | Triangular surface tiling method for silent area based on home environment | |
| CN109783969A (en) | A method of different-thickness homogeneous material acoustic characteristic is calculated using transfer matrix | |
| CN117233266A (en) | Full waveform inversion guided wave tomography method based on cyclic neural network | |
| CN109188378B (en) | Ultrasonic Doppler echo signal simulation method | |
| CN110261472B (en) | Method for measuring and reconstructing hearth section flow field by using wavelet basis function through acoustic wave method | |
| CN108304653B (en) | Quantitative and rapid evaluation technology for environmental vibration influence of crystal cultivation process |
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 |