CN107172534A - Build the method that multiple rows of manual simulation rotates sound source generators - Google Patents
Build the method that multiple rows of manual simulation rotates sound source generators Download PDFInfo
- Publication number
- CN107172534A CN107172534A CN201710472621.3A CN201710472621A CN107172534A CN 107172534 A CN107172534 A CN 107172534A CN 201710472621 A CN201710472621 A CN 201710472621A CN 107172534 A CN107172534 A CN 107172534A
- Authority
- CN
- China
- Prior art keywords
- loudspeaker
- row
- multiple rows
- phase
- sound source
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2430/00—Signal processing covered by H04R, not provided for in its groups
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
Abstract
The invention provides a kind of method for building multiple rows of manual simulation's rotation sound source generators, including：Based on target modalities, often row's quantity of loudspeaker, the audible frequency of loudspeaker and for arranging the physical parameter needed for multiple rows of loudspeaker, calculate the respective phase of each loudspeaker of all rows；By the function generator program in computer program, audible frequency, sounding amplitude and the sounding phase of each loudspeaker are controlled respectively, wherein the audible frequency and sounding amplitude of each loudspeaker are identicals；Power amplifier is output to via output card, then by the signal output by amplification to each loudspeaker；The sound sent by all loudspeakers of microphone array synchronous acquisition being arranged circumferentially installed in downstream；And via signal amplifier and capture card, the voice signal gathered is input to computer, the corresponding sound pressure level of each glottis neoplasms is calculated by the computer program, sound source generators are rotated to evaluate constructed multiple rows of manual simulation.Sound source generators are rotated present invention also offers the multiple rows of manual simulation implemented according to this method.
Description
Technical field
Sound source generators technical field is rotated the present invention relates to manual simulation, in particular it relates to multiple rows of artificial for building
The method of simulation rotation sound source generators.
Background technology
Since 21 century, due to the requirement more and more higher to aircraft noise compulsory index, acoustic target is increasingly
It is difficult to realize by passive acoustic means.In the case, by active noise controlling realize aero-engine noise reduction by
Increasing concern.Specifically, one of the Main Noise Sources of aircraft during take-off and landing are the rotors of fan blade
Stator is interfered.The sixties in last century Tyler and Sofrin point out that blade row interactions noise is harmonic wave, and one can be inspired
Individual certain moduli state as blade function.These harmonic waves are all blade passing frequency (BPF) integral multiples, as shown in figure 1, simultaneously
And wider frequency noise has 30dB sound pressure level difference.Active noise controlling target is exactly to reduce the noise of these harmonic waves.
Substantial amounts of tested shows that active noise controlling can significantly reduce noise.And then can be more by analogue simulation
Further analyze experimental result and predict noise reduction.The sound field being dominant in order to manual simulation's different modalities, artificial mould
Intend rotation sound source generators to arise at the historic moment.The instrument is dominant by producing modality-specific installed in the loudspeaker array of pipeline wall
Sound field, main sound field has mode that is dominant.The different mode that are dominant are realized by changing loudspeaker control rule, and can
Further to simulate aero-engine, i.e. entrance, fan, the noise at stator blade and nozzle.
At present, existing manual simulation rotation sound source generators by pipeline wall last week to uniformly installing single raise one's voice
Device array, changes loudspeaker and occurs rule to realize sound field that artificial simulation different modalities are dominant.Specifically, single artificial mould
The operation principle for intending rotation sound source generators is as follows：
Consider that fluid is applied the situation of the sounding of unsteady load by blade.First simplify pipeline and leaf model, expire it
Foot three-dimensional lifting surface theory basic assumption：
1) inner-walls of duct is rigid plane.It is without glutinous, compressible, constant entropy, subacoustic Uniform Flow in pipe.
2) plate blade, and ignore thickness effect.
3) flowing incidence angle is the situation that 0 °, the i.e. angle of attack are 0 °, only considers the unsteady load on blade.
Can be what dipole source was produced by problem reduction to be studied according to assumed above and Lighthill acoustic analogy theories
Sound field, its acoustic pressure is：
Wherein：x(x1, x2, x3) it is point of observation position under rest frame；y(y1, y2, y3) it is sound source position under rest frame
Put；T is observing time；τ is phonation time；S (y) is object boundary surface；fyBy object its surfactant fluid is applied it is non-
Permanent power；G is Green's function, and for the circular ring pipe of compressor, Green's function is represented by
In above formulaFor point of observation position column coordinate representation form under rest frame；Sat to be static
The lower sound source position cylindrical coordinates representation of mark system；α is the wave number in z directions；ω is audible frequency；k0For sound wave wave number；M is equal with n
For integer, the circumference and radial mode shape exponent number of sound field in pipeline are represented respectively.
Assuming that the rotor blade number included in rotor blade row is B, then total acoustic pressure that all blades are produced can must be superimposed
For：
Consider above formula Leaf number summation part.More in view of rotor blade number, it is integer typically to take σ B/2 π.
In summary derive and understand, the lobe numbers of the asynchronous concussion sounding of single rotor blade, rotating speed, blade vibration
There is following relation in frequency, blade phase difference, circumferential mode：
ω=Ωs+mΩ
Wherein s, B, m are integer, and be can be seen that for the unique circumferential mode of each sound field frequency correspondence.According to
Upper result, replaces with loudspeaker, i.e. B by blade and represents loudspeaker number, speaker sound frequency can simulate vibration frequency of blade,
Phase difference σ between adjacent loudspeakers can then be obtained manual simulation and rotated sound source with phase angle between model rotor blade vibration
Control rule of the device under the conditions of single loudspeaker.
Analyzed more than, the quantity of loudspeaker is directly related with the highest mode that is dominant that can be produced, for example：If
Number of loudspeakers is 2N, then this manual simulation rotation sound source generators can only produce the sound field that-N~+N mode is dominant, i.e., most
High mode number is N.However, the mode that aero-engine is produced is general higher, and existing single manual simulation's rotation sound source hair
The quantity for the single loudspeaker that the highest mode that raw device can be generated can be arranged by limited inner-walls of duct circumferential size
Limitation.
For example, on experimental bench, inner-walls of duct uniformly installs 9 loudspeakers in one week, if it is desired to realize that+1 mode is dominant
Sound field, then according to below equation：
It can be calculated that the phase difference of adjacent loudspeakers should be set to -40 degree.
However, calculating m=according to above-mentioned condition ..., -17, -8,1,10,19 ..., although the cut-off effect of pipeline can handle
High order mode is eliminated, but when performed at a sufficiently high frequency, inevitably the multiple mode that are dominant, -5 moulds in such as Fig. 9
State.In this case, the sound field that single loudspeaker the is produced situation that is dominant is also undesirable.
The content of the invention
In order to overcome the influence that existing single manual simulation's rotation sound source generators are limited by number of loudspeakers, this hair
Bright to propose a kind of method for building multiple rows of manual simulation's rotation sound source generators, it to first row loudspeaker by mending
Repay, mode can not only be simulated and be dominant the more significant sound field of situation, but also the limitation of number of loudspeakers, mould can be broken through
Draw up the sound field of higher mode.
According to an aspect of the present invention the method that multiple rows of manual simulation rotates sound source generators, bag are built there is provided a kind of
Include：
Based on target modalities, often row's quantity of loudspeaker, the audible frequency of loudspeaker and for arranging multiple rows of loudspeaker
Required physical parameter, calculates the respective phase of each loudspeaker of all rows；
By the function generator program in computer program, audible frequency, the sounding width of each loudspeaker are controlled respectively
Value and sounding phase, wherein the audible frequency and sounding amplitude of each loudspeaker are identicals；
Power amplifier is output to via output card, then by the signal output by amplification to each loudspeaker；
The sound sent by all loudspeakers of microphone array synchronous acquisition being arranged circumferentially installed in downstream；
And
Via signal amplifier and capture card, the voice signal gathered is input to computer, passes through the computer
Program calculates the corresponding sound pressure level of each glottis neoplasms, and sound source generators are rotated to evaluate constructed multiple rows of manual simulation.
According to embodiment, the respective phase of each loudspeaker of all rows is calculated, can be included：
Based on the target modalities, the quantity for often arranging loudspeaker, calculate adjacent in often row using below equation
Phase difference σ between loudspeaker,
Wherein, m is the target modalities, and S is integer, and B is the quantity for often arranging loudspeaker；
Audible frequency and the physical parameter for arranging needed for multiple rows of loudspeaker based on the loudspeaker, are used
Below equation arranges compensation phase of the loudspeaker relative to the 1st row's loudspeaker to calculate jth
Wherein, the physical parameter being used for needed for arranging multiple rows of loudspeaker includes：kzFor axial wave number；Z is loudspeaker
The distance of Propagation；ω is angular frequency, and ω=2 π f, f are the audible frequency of loudspeaker；c0For local velocity of sound；Ma is flowing horse
Conspicuous number；krFor radial direction wave number, the Size calculation of pipeline that can be by the target modalities and for arranging multiple rows of loudspeaker is obtained；L
For the distance between adjacent two rows loudspeaker；N is total row of multiple rows of loudspeaker；And
Phase in the often row for arrange the phase of the 1st loudspeaker based on the given the 1st, being calculated between adjacent loudspeaker
Poor σ and compensation phase of the jth row's loudspeaker relative to the 1st row's loudspeaker calculatedRaised one's voice to calculate each of all rows
The respective phase of device.
According to embodiment, the respective phase of each loudspeaker of all rows is calculated, can also be included：
Based on the quantity for often arranging loudspeaker, jth row's loudspeaker is calculated using below equation and is raised one's voice relative to the 1st row
The angle, θ of device rotationj,
Based on the target modalities, the quantity for often arranging loudspeaker, two adjacent rows are calculated using below equation and are raised
Phase difference σ between sound device,
Wherein, m is the target modalities, and S is integer, and B is the quantity for often arranging loudspeaker；
Based on the phase difference σ between the two adjacent row's loudspeakers calculated, calculate between loudspeaker adjacent in often arranging
Phase difference n* σ；
Audible frequency and the physical parameter for arranging needed for multiple rows of loudspeaker based on the loudspeaker, are used
Below equation arranges compensation phase of the loudspeaker relative to the 1st row's loudspeaker to calculate jth
Wherein, the physical parameter being used for needed for arranging multiple rows of loudspeaker includes：kzFor axial wave number；Z is loudspeaker
The distance of Propagation；ω is angular frequency, and ω=2 π f, f are the audible frequency of loudspeaker；c0For local velocity of sound；Ma is flowing horse
Conspicuous number；krFor radial direction wave number, the Size calculation of pipeline that can be by the target modalities and for arranging multiple rows of loudspeaker is obtained；L
For the distance between adjacent two rows loudspeaker；N is total row of multiple rows of loudspeaker；And
Based on the given the 1st arrange the phase of the 1st loudspeaker, calculate phase difference σ between two adjacent row's loudspeakers,
Phase difference n* σ in the often row calculated between adjacent loudspeaker raise one's voice with the jth row's loudspeaker calculated relative to the 1st row
The compensation phase of deviceTo calculate the respective phase of each loudspeaker of all rows.
According to embodiment, multiple rows of loudspeaker arrangement be axial direction of each row's loudspeaker along the pipeline in the plane
Projection form the equally distributed array of circumference of B loudspeaker.
According to embodiment, multiple rows of loudspeaker arrangement be axial direction of each row's loudspeaker along the pipeline in the plane
Projection form the equally distributed array of circumference of n × B loudspeaker.
According to embodiment, the 1st given phase for arranging the 1st loudspeaker is preferably zero.
According to embodiment, the computer program is Labview computer programs.
According to another aspect of the present invention there is provided one kind according to any above-mentioned multiple rows of manual simulation's rotation sound source hair of structure
The rotation sound source generating means that the method for raw device is implemented.
Brief description of the drawings
Fig. 1 is aircraft engine fan blade blade row interactions noise pattern.
Fig. 2 is the flow chart for the method that sound source generators are rotated according to the multiple rows of manual simulation of structure of the embodiment of the present invention.
Fig. 3 is the multiple rows of of the method implementation of the multiple rows of manual simulation's rotation sound source generators of structure according to the embodiment of the present invention
Manual simulation rotates the profile of sound source generators.
Fig. 4 is the multiple rows of loudspeaker installation position horizontalization of more preferable sound field of being dominant according to the formation mode of one embodiment of the invention
Face expanded schematic diagram.
Fig. 5 be according to the is formationed mode of one embodiment of the invention be dominant more preferable sound field example 3 row's loudspeakers installation
Position plane expanded schematic diagram.
Fig. 6 is multiple rows of loudspeaker installation site plane exhibition of the higher mode sound field of formation according to another embodiment of the present invention
Open schematic diagram.
Fig. 7 is 2 row's loudspeaker installation sites of the example of the higher mode sound field of formation according to another embodiment of the present invention
Planar development schematic diagram.
Fig. 8 for Fig. 72 row's loudspeakers project axially along the pipe after schematic diagram.
Existing single manual simulation's rotation sound source generators and according to embodiments of the present invention many are shown respectively in Fig. 9 and 10
Row's manual simulation's rotation sound source generators produce the mode of target+4 and are dominant the design sketch of sound field.
Existing single manual simulation rotation sound source generators are shown respectively and according to embodiments of the present invention in Figure 11 and 12
Multiple rows of manual simulation rotation sound source generators produce 7 mode and are dominant the design sketch of sound field.
Embodiment
The multiple rows of manual simulation's rotation sound source generators of structure according to embodiments of the present invention are specifically described referring to the drawings
Method.
Fig. 2 shows the flow of the method for the multiple rows of manual simulation's rotation sound source generators of structure according to embodiments of the present invention
Figure.As shown in Fig. 2 the method that multiple rows of manual simulation rotates sound source generators is built, including：
Based on target modalities, often row's quantity of loudspeaker, the audible frequency of loudspeaker and for arranging multiple rows of loudspeaker
Required physical parameter, calculates the respective phase of each loudspeaker of all rows；
By the function generator program in computer program, audible frequency, the sounding width of each loudspeaker are controlled respectively
Value and sounding phase, wherein the audible frequency and sounding amplitude of each loudspeaker are identicals；
Power amplifier is output to via output card, then by the signal output by amplification to each loudspeaker；
The sound sent by all loudspeakers of microphone array synchronous acquisition being arranged circumferentially installed in downstream；
And
Via signal amplifier and capture card, the voice signal gathered is input to computer, passes through the computer
Program calculates the corresponding sound pressure level of each glottis neoplasms, and sound source generators are rotated to evaluate constructed multiple rows of manual simulation.
It is many that Fig. 3 shows that the method for the multiple rows of manual simulation's rotation sound source generators of structure according to embodiments of the present invention is implemented
Arrange the profile that manual simulation rotates sound source generators.Installed as shown in figure 3, being circumferentially provided with multiple rows of loudspeaker in inner-walls of duct
Hole 301, for installing multiple rows of loudspeaker；And conveyer array 302 is circumferentially provided with the downstream position of inner-walls of duct, is used for
Receive the voice signal that multiple rows of loudspeaker is sent.
Formation mode according to an embodiment of the invention is described below in conjunction with Fig. 4 and Fig. 5 to be dominant the multiple rows of of more preferable sound field
Loudspeaker arrangement.
Specifically, Fig. 4 shows concrete arrangement when n arranges loudspeaker, each row has B loudspeaker altogether.
As shown in figure 4, be dominant more preferably sound field to form mode, each row's loudspeaker in multiple rows of loudspeaker is along pipe
The projection of road axial direction should be completely superposed, and the sounding rule of each row's loudspeaker should all be met：
Wherein, m is target modalities, and S is integer, and B is the quantity of the loudspeaker in each row, and σ is adjacent in each row raises
The phase difference of sound device.
As shown in figure 4, for n row's loudspeakers, No. 1 loudspeaker of each row's loudspeaker, i.e. numbering are 11,21,31,
41st ..., n1 loudspeaker (wherein, first of numbering is row number, and second is when front-seat sequence number), i.e. difference row but
Between the adjacent loudspeaker of numbering (11 and 21,21 and 31 ...) there should be phase difference, the difference for compensating Propagation distance
It is different, make it that there is identical phase with the acoustic pressure of the 1st row's loudspeaker here when traveling to microphone array, play enhancing and make
With.
The method for calculating the phase of each loudspeaker of all rows is as follows：
The acoustic pressure function of each loudspeaker can be expressed as：
Wherein, A is sound pressure amplitude, and ω is angular frequency, and ω=2 π f, f are the audible frequency of loudspeaker, kzFor axial wave number,
Z is the distance of sound Propagation.
Wherein, k0For wave number；krFor radial direction wave number, it can be calculated and obtained by target modalities and line size；Ma is flowing Mach
Number；c0For local velocity of sound.
In embodiments of the present invention, loudspeaker array is synchronous cant, and microphone array is synchronous acquisition.In order to protect
Card at a time, (sound transmission z during jth row's loudspeaker pass to microphone arrayjDistance), the 1st row's loudspeaker is existed
Acoustic pressure herein there is enhancing effect, it is necessary to jth arrange loudspeaker phase be modified, just its sound can be made to pass to microphone
During array position, the sound pressure phase of two row's loudspeakers is identical, if orrection phase place isAlso referred to as compensate phase.
If at the time of microphone pick being t0, then after No. 1 loudspeaker sounding of o.11 and jth, at microphone array
Acoustic pressure is：
Wherein, z1=S, zj=S+ (j-1) * L, j=2,3 ..., n
Ensure that phase is identical, i.e.,
As using the orrection phase place of jth row during n-th row's loudspeaker.It can be seen that orrection phase place is unrelated with S.
When given No. 11 loudspeaker phases are 0 °, based on the phase difference between adjacent loudspeakers in each row, adjacent two
The orrection phase place between No. 1 loudspeaker in the orrection phase place between loudspeaker, such as adjacent two row is arranged, all rows can be obtained
In each loudspeaker phase.
Fig. 5 show to be formed mode be dominant more preferable sound field 3 row's loudspeaker arrangements example.Specifically, as shown in figure 5,
In three row manual simulations rotation sound source generators, experiment condition is：Often row has 9 loudspeakers, no flowing, 20 degree of temperature, pipe
Road internal diameter 600mm, target modalities are+2, and audible frequency f is 1000Hz, adjacent two rows loudspeaker spacing L=0.1m.
According toWherein m=+2, B=9, draw σ=- 80 °, if setting No. 11 loudspeaker phases
For 0 °, then the phase of 12 to No. 19 loudspeakers is -80 ° respectively, -160 °, -240 °, -320 °, -400 °, -480 °, -560 °, -
640°。
According to formula
As j=2, c0=343m/s (mainly relevant with temperature), kr=8.726, ω=6283, Ma=0, then phase difference
Thus, based on the phase difference σ between loudspeaker adjacent in 0 ° given of No. 11 loudspeaker phases, each row=-
Phase difference between 80 ° and adjacent two rows loudspeakerRaising one's voice for numbering 21 can be calculated
Device phase is 0+92.25 °=92.25 °, and No. 22 loudspeaker phases are 92.25 ° -80 °=12.25 °, 23 to No. 29 loudspeaker phases
Position is respectively -67.75 °, -147.75 °, -227.75 °, -307.75 °, -387.75 °, -467.75 °, -547.75 °.
Similarly,I.e. the loudspeaker phase of numbering 31 is 0 °+184.5 °=184.5 °,
No. 32 loudspeaker phases are 184.5 ° -80 °=104.5 °, and 33 to No. 39 loudspeaker phases are respectively 24.5 °, -55.5 °, -
135.5 °, -215.5 °, -295.5 °, -375.5 °, -455.5 °.
Higher the multiple rows of of mode sound field of formation according to another embodiment of the present invention is described below in conjunction with Fig. 6 and Fig. 7 to raise one's voice
Device is arranged.
Fig. 6 shows n row's loudspeakers and the often concrete arrangement of B loudspeaker of row, wherein, in order to form more Gao Mo
State sound field, by loudspeaker array relative to the 1st row array row rotate to an angle, as shown in fig. 6, for jth row (j=2,
3 ..., n) loudspeaker array, the angle that should be rotated relative to the 1st row is θj：
Thus, n, which arranges loudspeaker, can regard every one formation battle array of n selection from single loudspeaker (n*B loudspeaker) as
Arrange and do the axial translation of different distance and obtain.That is, n arranges loudspeaker, along pipeline, axially projection in the plane can be with
Form the circumferentially equally distributed array of n*B loudspeaker.
In order to form Target Sound Field, it is necessary to which n × B loudspeaker arrays after projecting are met
Now the phase difference between adjacent two rows loudspeaker isThat is, numbering be 11,21,31,41 ... n1
Two neighboring loudspeaker between phase difference be σ.And the phase difference between two loudspeakers adjacent in each row, that is, number
For n1, n2, n3 ..., the phase difference between nB two neighboring loudspeaker be n* σ.
Because phase difference now is after projection is translated, when the state being returned to before projection, it is necessary to be raised to jth row
The phase of sound device is modified, it is ensured that it is consistent with the phase that acoustic pressure after projection is passed at microphone array before projecting.
The acoustic pressure function of each loudspeaker can be expressed as：
Wherein, A is sound pressure amplitude, and ω is angular frequency, and ω=2 π f, f are loudspeaker audible frequency, kzFor axial wave number, z
For the distance of sound Propagation.
Wherein, k0For wave number；krFor radial direction wave number, it can be calculated by target modalities and line size；Ma is flowing Mach
Number；c0For local velocity of sound.
In the present embodiment, loudspeaker array is synchronous cant and microphone array is synchronous acquisition, in order to ensure
A certain moment, (sound transmission z during jth row loudspeaker pass to microphone array before projectionjDistance), with after projection (to the 1st
Row's loudspeaker plane) acoustic pressure of the loudspeaker at microphone array there is same phase, it is necessary to arrange jth the phase of loudspeaker
It is modified, if orrection phase place is
If at the time of microphone array is gathered being t0, then jth 1 is with projecting preceding No. 1 loudspeaker of jth in microphone after projecting
Acoustic pressure at array is：
Wherein, z1=S, zj=S+ (j-1) * L, j=2,3 ..., n
Ensure that phase is identical, i.e.,
As using the orrection phase place of jth row during n-th row's loudspeaker.It can be seen that orrection phase place is unrelated with S.
Thus, based on the phase difference between adjacent loudspeakers in No. 11 given loudspeaker phases, each row, adjacent two
The orrection phase place of the phase difference and each row's loudspeaker between loudspeaker is arranged, the phase of all loudspeakers can be obtained.
Fig. 7 shows the example for the 2 row's loudspeaker arrangements to form the higher mode sound field that is dominant.As shown in fig. 7, the 2nd row raises one's voice
Device rotates in a circumferential direction relative to the 1st rowBy the 2nd row to the 1st
The circumferentially equally distributed array of 2 × 9 loudspeakers as shown in Figure 8 is formed after the projection of row section.
Assuming that experiment condition is：Often arrange 9 loudspeakers, no flowing, 20 degree of temperature, internal diameter of the pipeline (diameter) 600mm, target
Mode is+7, and frequency is 2000Hz, adjacent two rows loudspeaker spacing L=0.1m.
According toWherein m=+7, B=9, obtain the phase difference between adjacent two rows loudspeakerIf for example setting No. 11 loudspeaker phases as 0 °, the phase of No. 21 is 0 ° of+σ=- 140 °；
And based on the phase difference n* σ between two loudspeakers adjacent in each row, it is 0 ° of+n* σ to obtain the phase of No. 12 loudspeakers
=-280 °, No. 13 phases to No. 19 loudspeakers are -560 ° respectively, -840 °, -1120 °, -1400 °, -1680 °, -1960 °, -
2240°；22 to No. 29 loudspeaker phases are respectively：- 420 °, -700 °, -980 °, -1260 °, -1540 °, -1820 °, -
2100 °, -2380 °
Meanwhile, second row loudspeaker should also have orrection phase place：
As m=+7, c0=343m/s, kr=28.593, ω=2* π * f=12566, Ma=0, L=0.1m
ThenThus the loudspeaker phase for obtaining numbering 21 is -140 °+131.24 °
=-8.76 °, No. 22 loudspeaker phases are -8.76 ° -280 °=- 288.76 °, 23 to No. 29 loudspeaker phases are changed into respectively -
568.76 °, -848.76 °, -1128.76 °, -1408.76 °, -1688.76 °, -1968.76 °, -2248.76 °.
Existing single manual simulation's rotation sound source generators and according to embodiments of the present invention many are shown respectively in Fig. 9 and 10
Row's manual simulation's rotation sound source generators produce 4 mode and are dominant the design sketch of sound field.Figure 11 and 12 is shown respectively existing single
Manual simulation rotates sound source generators and multiple rows of manual simulation rotation sound source generators according to embodiments of the present invention produce 7 mode
The schematic diagram for the sound field that is dominant.The purpose tested twice it can be seen from Fig. 9 and Figure 10 is to produce the sound field that+4 mode are dominant.
Fig. 9 is the sound field that single 9 loudspeakers are produced, and -5 mode correspondence sound pressure level is very high, but the amount of being dominant only 5dB, and Figure 10 is double
The sound field that 18 loudspeakers of row are produced, it can be seen that it restrained effectively the sound of -5 mode, while enhancing+4 mode
Sound, makes the amount of being dominant reach 20dB.
The purpose tested twice it can be seen from Figure 11 and Figure 12 is to produce the sound field that 7 mode are dominant, and single 9 are raised
The sound field that sound device is produced ,+7 mode correspondence sound pressure level is not highest, -2 mode highest on the contrary, it was demonstrated that the structure of this sound field is
Failure, and it can be recognized from fig. 12 that double loudspeaker restrained effectively the sound of -2 mode, while enhancing+7 mode
Sound, makes the amount of being dominant reach 12dB.
In summary, the invention provides a kind of method for building multiple rows of manual simulation's rotation sound source generators and according to this
Multiple rows of manual simulation rotation sound source generators that method is implemented, it, can not only mould by being compensated to first row loudspeaker
Draw up mode to be dominant the more significant sound field of situation, but also the limitation of number of loudspeakers can be broken through, simulate higher mode
Sound field.
In addition, rotating multiple rows of artificial mould that the method for sound source generators is implemented according to the multiple rows of manual simulation of structure of the present invention
Intending rotation sound source generators has a variety of functions, can carry out kinds of experiments, such as performance test of sound lining, absorbing property test
Deng.And it is possible to realize the simulation of turbomachine internal acoustic field using relatively simple structure, and complete related based on the sound field
Research.More simple single mode or multi-modal sound field is generated to a certain extent, is easy to complete some specific mechanisms researchs, for example
Research of blade row's reflectance factor etc..
Claims (8)
1. a kind of method for building multiple rows of manual simulation's rotation sound source generators, including：
Based on target modalities, often row's quantity of loudspeaker, the audible frequency of loudspeaker and for arranging needed for multiple rows of loudspeaker
Physical parameter, calculate the respective phase of each loudspeaker of all rows；
By the function generator program in computer program, control respectively the audible frequency of each loudspeaker, sounding amplitude and
Sounding phase, wherein the audible frequency and sounding amplitude of each loudspeaker are identicals；
Power amplifier is output to via output card, then by the signal output by amplification to each loudspeaker；
The sound sent by all loudspeakers of microphone array synchronous acquisition being arranged circumferentially installed in downstream；And
Via signal amplifier and capture card, the voice signal gathered is input to computer, passes through the computer program
The corresponding sound pressure level of each glottis neoplasms is calculated, sound source generators are rotated to evaluate constructed multiple rows of manual simulation.
2. according to claim 1 build the method that multiple rows of manual simulation rotates sound source generators, wherein, calculate all rows
The respective phase of each loudspeaker, including：
Based on the target modalities, the quantity for often arranging loudspeaker, adjacent in often row raise one's voice is calculated using below equation
Phase difference σ between device,
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>S</mi>
<mi>B</mi>
<mo>-</mo>
<mi>&sigma;</mi>
<mfrac>
<mi>B</mi>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
</mrow>
Wherein, m is the target modalities, and S is integer, and B is the quantity for often arranging loudspeaker；
Audible frequency and the physical parameter for arranging needed for multiple rows of loudspeaker based on the loudspeaker, using following
Equation arranges compensation phase of the loudspeaker relative to the 1st row's loudspeaker to calculate jth
Wherein, the physical parameter being used for needed for arranging multiple rows of loudspeaker includes：kzFor axial wave number；Z is sound Propagation
Distance；ω is angular frequency, and ω=2 π f, f are the audible frequency of loudspeaker；c0For local velocity of sound；Ma is flowing Mach number；krFor
Radial direction wave number, the Size calculation of pipeline that can be by the target modalities and for arranging multiple rows of loudspeaker is obtained；L is adjacent two row
The distance between loudspeaker；N is total row of multiple rows of loudspeaker；And
Phase difference σ in the often row for arrange the phase of the 1st loudspeaker based on the given the 1st, being calculated between adjacent loudspeaker
With compensation phase of the jth row's loudspeaker relative to the 1st row's loudspeaker calculatedIt is each come each loudspeaker for calculating all rows
From phase.
3. according to claim 1 build the method that multiple rows of manual simulation rotates sound source generators, wherein, calculate all rows
The respective phase of each loudspeaker, including：
Based on the quantity for often arranging loudspeaker, jth row's loudspeaker is calculated using below equation and is revolved relative to the 1st row's loudspeaker
The angle, θ turnedj,
Based on the target modalities, the quantity for often arranging loudspeaker, two adjacent row's loudspeakers are calculated using below equation
Between phase difference σ,
<mrow>
<mi>m</mi>
<mo>=</mo>
<mi>n</mi>
<mi>S</mi>
<mi>B</mi>
<mo>-</mo>
<mi>&sigma;</mi>
<mfrac>
<mrow>
<mi>n</mi>
<mi>B</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
</mrow>
1
Wherein, m is the target modalities, and S is integer, and B is the quantity for often arranging loudspeaker；
Based on the phase difference σ between the two adjacent row's loudspeakers calculated, the phase between loudspeaker adjacent in often arranging is calculated
Potential difference n* σ；
Audible frequency and the physical parameter for arranging needed for multiple rows of loudspeaker based on the loudspeaker, using following
Equation arranges compensation phase of the loudspeaker relative to the 1st row's loudspeaker to calculate jth
Wherein, the physical parameter being used for needed for arranging multiple rows of loudspeaker includes：kzFor axial wave number；Z is sound Propagation
Distance；ω is angular frequency, and ω=2 π f, f are the audible frequency of loudspeaker；c0For local velocity of sound；Ma is flowing Mach number；krFor
Radial direction wave number, the Size calculation of pipeline that can be by the target modalities and for arranging multiple rows of loudspeaker is obtained；L is adjacent two row
The distance between loudspeaker；N is total row of multiple rows of loudspeaker；And
The phase of the 1st loudspeaker is arranged based on the given the 1st, phase difference σ between two adjacent row's loudspeakers is calculated, is counted
Phase difference n* σ in the often row calculated between adjacent loudspeaker arrange loudspeaker relative to the 1st row's loudspeaker with the jth calculated
Compensate phaseTo calculate the respective phase of each loudspeaker of all rows.
4. according to claim 2 build the method that multiple rows of manual simulation rotates sound source generators, wherein, it is described multiple rows of to raise
Sound device is arranged as the projection of axial direction of each row's loudspeaker along the pipeline in the plane and forms the circumference of B loudspeaker uniformly
The array of distribution.
5. according to claim 3 build the method that multiple rows of manual simulation rotates sound source generators, wherein, it is described multiple rows of to raise
Sound device is arranged as the circumference that the projection of axial direction of each row's loudspeaker along the pipeline in the plane forms n × B loudspeaker
Equally distributed array.
6. the method that the multiple rows of manual simulation of structure according to claim 4 or 5 rotates sound source generators, wherein, it is described to give
The 1st fixed phase for arranging the 1st loudspeaker is preferably zero.
7. according to claim 6 build the method that multiple rows of manual simulation rotates sound source generators, wherein, the computer
Program is Labview computer programs.
8. a kind of multiple rows of manual simulation of structure according to any one of claim 1 to 7 rotates the side of sound source generators
Multiple rows of manual simulation rotation sound source generating means that method is implemented.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710472621.3A CN107172534B (en) | 2017-06-20 | 2017-06-20 | Construct the method that multiple rows of manual simulation rotates sound source generators |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710472621.3A CN107172534B (en) | 2017-06-20 | 2017-06-20 | Construct the method that multiple rows of manual simulation rotates sound source generators |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107172534A true CN107172534A (en) | 2017-09-15 |
| CN107172534B CN107172534B (en) | 2019-07-09 |
Family
ID=59819496
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710472621.3A Active CN107172534B (en) | 2017-06-20 | 2017-06-20 | Construct the method that multiple rows of manual simulation rotates sound source generators |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107172534B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111523220A (en) * | 2020-04-17 | 2020-08-11 | 中国空气动力研究与发展中心高速空气动力研究所 | Fan and compressor rotating-static interference broadband noise rapid prediction method considering flow influence |
| CN113179465A (en) * | 2021-04-15 | 2021-07-27 | 上海交通大学 | Pipeline high-order acoustic modal generating device and control method |
| CN114743562A (en) * | 2022-06-09 | 2022-07-12 | 成都凯天电子股份有限公司 | Method and system for recognizing airplane voiceprint, electronic equipment and storage medium |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5221185A (en) * | 1991-08-05 | 1993-06-22 | General Electric Company | Method and apparatus for synchronizing rotating machinery to reduce noise |
| CN101242683A (en) * | 2007-02-09 | 2008-08-13 | 中国科学院声学研究所 | Digital coding and modulation airflow speaker |
| CN101607594A (en) * | 2008-01-30 | 2009-12-23 | 尤洛考普特公司 | The letter road reactive torque blade of optimizing autogiro is with method that minimizes noise jamming and the letter road reactive torque blade that obtains thus |
| CN104406756A (en) * | 2014-11-11 | 2015-03-11 | 沈阳黎明航空发动机(集团)有限责任公司 | System and method for simulating single/dual-sound wave stimulated vibration test of working environment of engine blade |
| CN104751835A (en) * | 2015-03-31 | 2015-07-01 | 中国飞机强度研究所 | Rotary sound source generating device |
| CN106289505A (en) * | 2016-07-21 | 2017-01-04 | 合肥工业大学 | A kind of method separating static sound source radiation sound field and rotation sound source radiated sound field |
-
2017
- 2017-06-20 CN CN201710472621.3A patent/CN107172534B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5221185A (en) * | 1991-08-05 | 1993-06-22 | General Electric Company | Method and apparatus for synchronizing rotating machinery to reduce noise |
| CN101242683A (en) * | 2007-02-09 | 2008-08-13 | 中国科学院声学研究所 | Digital coding and modulation airflow speaker |
| CN101607594A (en) * | 2008-01-30 | 2009-12-23 | 尤洛考普特公司 | The letter road reactive torque blade of optimizing autogiro is with method that minimizes noise jamming and the letter road reactive torque blade that obtains thus |
| CN104406756A (en) * | 2014-11-11 | 2015-03-11 | 沈阳黎明航空发动机(集团)有限责任公司 | System and method for simulating single/dual-sound wave stimulated vibration test of working environment of engine blade |
| CN104751835A (en) * | 2015-03-31 | 2015-07-01 | 中国飞机强度研究所 | Rotary sound source generating device |
| CN106289505A (en) * | 2016-07-21 | 2017-01-04 | 合肥工业大学 | A kind of method separating static sound source radiation sound field and rotation sound source radiated sound field |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111523220A (en) * | 2020-04-17 | 2020-08-11 | 中国空气动力研究与发展中心高速空气动力研究所 | Fan and compressor rotating-static interference broadband noise rapid prediction method considering flow influence |
| CN113179465A (en) * | 2021-04-15 | 2021-07-27 | 上海交通大学 | Pipeline high-order acoustic modal generating device and control method |
| CN113179465B (en) * | 2021-04-15 | 2022-08-19 | 上海交通大学 | Pipeline high-order acoustic modal generating device and control method |
| CN114743562A (en) * | 2022-06-09 | 2022-07-12 | 成都凯天电子股份有限公司 | Method and system for recognizing airplane voiceprint, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107172534B (en) | 2019-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rienstra et al. | A numerical comparison between the multiple-scales and finite-element solution for sound propagation in lined flow ducts | |
| Holste et al. | Noise source identification in a propfan model by means of acoustical near field measurements | |
| Posson et al. | On a uniformly valid analytical rectilinear cascade response function | |
| CN107172534A (en) | Build the method that multiple rows of manual simulation rotates sound source generators | |
| McALPINE et al. | “Buzz-saw” noise: Prediction of the rotor-alone pressure field | |
| Doherty et al. | Impact of turbofan intake distortion on fan noise propagation and generation | |
| US9245089B1 (en) | Statistically based approach to broadband liner design and assessment | |
| Mathews et al. | The acoustic Green's function for swirling flow in a lined duct | |
| Casalino et al. | Lattice-Boltzmann simulation of coaxial jet noise generation | |
| Shur et al. | Experimental/numerical study of ducted-fan noise: effect of duct inlet shape | |
| Weckmüller et al. | Cfd/caa coupling applied to dlr uhbr-fan: Comparison to experimantal data | |
| Sutliff | A 20 year retrospective of the advanced noise control fan–contributions to turbofan noise research | |
| Khemiri et al. | Asymptotic description of damping mistuning effects on the forced response of turbomachinery bladed disks | |
| Moreau et al. | Optimization of tonal noise control with flow obstruction | |
| Rossikhin et al. | Numerical investigation of high bypass ratio fan tone noise | |
| Shur et al. | Effect of inlet distortions on ducted fan noise | |
| Jaron et al. | RANS-informed fan noise prediction: separation and extrapolation of rotor wake and potential field | |
| Zhang et al. | Low-frequency fluctuations and their suppression in the 5.5 m× 4 m aeroacoustic wind tunnel at CARDC | |
| Nark et al. | Broadband liner optimization for the source diagnostic test fan | |
| Vorontsov et al. | Comparison of the noise directivity pattern of the main rotor of a helicopter for flight and hover modes | |
| Roger et al. | Aerodynamic and Sound-Scattering Effects in Rotor-Strut Interaction Noise of Small-Size Drones | |
| Choi et al. | Development of the numerical method for calculating sound radiation from a rotating dipole source in an opened thin duct | |
| Zhang et al. | Acoustic radiation from a semi-infinite duct with a subsonic jet | |
| Rynell et al. | Inclusion of upstream turbulent inflow statistics to numerically acquire proper fan noise characteristics | |
| Chen et al. | Aeroacoustic imaging experiments of airframe noise in lined wall closed-section wind tunnel |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Xiaoyu Inventor after: Li Zhibin Inventor after: Sun Xiaofeng Inventor after: Jing Xiaodong Inventor after: Sun Dakun Inventor after: Du Lin Inventor after: Wang Xuli Inventor before: Wang Xiaoyu Inventor before: Sun Xiaofeng Inventor before: Li Zhibin Inventor before: Jing Xiaodong Inventor before: Sun Dakun Inventor before: Du Lin Inventor before: Wang Xuli |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |