CN110146249A - A kind of pressure-measuring pipe road optimum design method for wind tunnel test - Google Patents

A kind of pressure-measuring pipe road optimum design method for wind tunnel test Download PDF

Info

Publication number
CN110146249A
CN110146249A CN201910342726.6A CN201910342726A CN110146249A CN 110146249 A CN110146249 A CN 110146249A CN 201910342726 A CN201910342726 A CN 201910342726A CN 110146249 A CN110146249 A CN 110146249A
Authority
CN
China
Prior art keywords
pressure
pipe road
measuring pipe
pressure signal
signal
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
Application number
CN201910342726.6A
Other languages
Chinese (zh)
Inventor
伞冰冰
邱冶
李剑鸿
马云龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201910342726.6A priority Critical patent/CN110146249A/en
Publication of CN110146249A publication Critical patent/CN110146249A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/02Wind tunnels
    • G01M9/04Details
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Abstract

The invention discloses a kind of pressure-measuring pipe road optimum design methods for wind tunnel test, the specific steps are as follows: and it is step 1, theoretical based on circuited transmission, obtain the theoretical frequency response function of pressure-measuring pipe road system;Step 2, the Optimized model of multistage pressure-measuring pipe road system is established;Step 3, using pressure-measuring pipe road inside radius or duct length as optimized variable, the optimal design parameters of pipe-line system are determined using genetic Optimization Algorithm;Step 4, according to two evaluation indexes of coefficient of wind pres time-histories and Power spectral density, effect of optimization is evaluated;Step 5, final wind tunnel test pressure-measuring pipe road system and its high-precision modification method are determined.The signal distortion problem that effective solution of the present invention wind tunnel test pressure-measuring pipe road generates when longer, it is simple and effective, and can carry out deeper research based on this.

Description

A kind of pressure-measuring pipe road optimum design method for wind tunnel test
Technical field
The present invention relates to a kind of pressure-measuring pipe road optimum design methods for wind tunnel test.
Background technique
Recent decades, electronics science field is fast-developing, and the appearance and innovation of the equipment such as sensor, scanning valve are made us Pulsating wind pressure can be measured more accurately to meet the Research Requirements in relation to Wind Engineering field.
In wind tunnel test, there are certain length for pressure-measuring pipe road system, this will certainly be to the pulsating pressure in pipe-line system Force signal has an impact, mainly acted on air damping based on.For heavy damping pipe-line system, it will usually cause fluctuation pressure Decaying on signal amplitude, and can then signal amplitude be caused to amplify underdamping pipe-line system.In turn, fluctuation pressure passes through one Signal distortion can be generated after measured length pressure-measuring pipe road system to sensor.Therefore, it is opened for the pressure-measuring pipe road system of wind tunnel test Opening up optimization design research has important theory significance and engineering practical value.
Summary of the invention
In order to solve the above technical problems, the present invention provides a kind of pressure-measuring pipe road optimum design method for wind tunnel test.
The present invention uses following technical scheme to solve above-mentioned technical problem:
The present invention provides a kind of pressure-measuring pipe road optimum design method for wind tunnel test, and steps are as follows for specific method:
Step 1: it is theoretical based on circuited transmission, obtain the frequency response function in pressure-measuring pipe road;
Step 2: establishing the Optimized model in multistage pressure-measuring pipe road;
Step 3: using pressure-measuring pipe road inside radius or duct length as optimized variable, being determined using genetic Optimization Algorithm and survey pressure The optimal design parameters of pipeline.
As a further optimization solution of the present invention, the frequency response function of pressure-measuring pipe road described in step 1 are as follows:
Wherein, m11It (f) is the top left hand element of the transfer matrix M of pressure-measuring pipe road whole story end parameters relationship, M=M1,2… Mn,n+1, n is the pipeline section number for constituting pressure-measuring pipe road, Mi,i+1For the transfer matrix of i-th of pipeline section, i=1,2 ..., n.
As further technical solution of the present invention, the Optimized model in the multistage pressure-measuring pipe road established in step 2 specifically:
Objective function min f (l1,l2,...,ln,r1,r2,...rn);
Constraint condition
Wherein, i=1,2 ..., n, liFor the length of i-th of pipeline section, L is pressure-measuring pipe road total length, and n is to constitute pressure-measuring pipe The pipeline section number in road, riFor the radius of i-th of pipeline section, ri u、ri lRespectively riUpper and lower limit,H (f) is the frequency response function in pressure-measuring pipe road, f0For pulsating wind pressure letter Number maximum frequency.
The present invention also provides a kind of modification methods of the pressure signal of pressure-measuring pipe road acquisition, and the pressure-measuring pipe road is using as above The optimum design method optimizes, the modification method the following steps are included:
Step A: frequency f=[0, fs/ 2] in range, if frequency response function H (f) ∝ 1 in pressure-measuring pipe road, pressure-measuring pipe The pressure signal of road acquisition is no to then follow the steps 2 without amendment;Wherein, fsFor the sample frequency of pressure scanning valve;
Step B: the pressure signal that pressure-measuring pipe road acquires is carried out according to coefficient of wind pres time-histories or Power spectral density Amendment, obtains true pressure signal.
As further technical solution of the present invention, pressure-measuring pipe road is acquired according to coefficient of wind pres time-histories in step B pressure Force signal is modified, and obtains true pressure signal, specifically:
1. being Fourier transform FFT to the pressure signal y (t) of pressure-measuring pipe road acquisition obtains fourier coefficient sequence in frequency domain Y(f);
2. obtaining Fourier transform coefficient sequence X (f)=Y (f)/H (f) of the true pressure signal x (t) in frequency domain;
3. carrying out inverse fourier transform to X (f) obtained in 2., true pressure signal x (t) is finally obtained.
As further technical solution of the present invention, pressure-measuring pipe road is acquired according to Power spectral density in step B Pressure signal be modified, obtain true pressure signal, specifically:
In formula,For pressure-measuring pipe road acquisition pressure signal power spectral density function,For true pressure The power spectral density function of signal.
The invention adopts the above technical scheme compared with prior art, has following technical effect that
1) present invention is simple clear, using convenient;Exploitativeness is strong;
2) pressure-measuring pipe road optimum design method proposed by the invention is able to solve multistage pressure-measuring pipe in wind tunnel pressure measuring test Signal distortion problem caused by road system improves the measurement accuracy of pulsating wind pressure signal.
Detailed description of the invention
Fig. 1 is a kind of flow chart of pressure-measuring pipe road optimization design and modification method for wind tunnel test;
Fig. 2 is the simple cross-section straight pipeline schematic diagram that length is l;
Fig. 3 is the common typical pressure-measuring pipe road schematic diagram of wind-tunnel;
Piping frequency response function curve when Fig. 4 is is calculating variable with pipeline inside radius before corresponding optimization and after optimization;
Piping frequency response function curve when Fig. 5 is is calculating variable with length before corresponding optimization and after optimization;
Fig. 6 is that the pressure-measuring pipe road for being 0.8m to length inputs four groups of different frequency sine pressure waves, utilizes the frequency response of pipeline The correction result comparison diagram that function is modified the pressure signal of output end, wherein (a) is input frequency f=20Hz, (b) It is (c) input frequency f=100Hz to input frequency f=60Hz, is (d) input frequency f=120Hz;
Fig. 7 is compared for the collected pulsating wind pressure in different length pressure-measuring pipe road under log-log coordinate from spectrum correction result Figure, wherein (a) is 0.5m length pressure-measuring pipe road, (b) is 0.8m length pressure-measuring pipe road, (c) is 1m length pressure-measuring pipe road, (d) For 1.2m length pressure-measuring pipe road;
Fig. 8 is length l=0.8m pressure-measuring pipe road spectrum error distribution schematic diagram certainly.
Specific embodiment
Embodiments of the present invention are described below in detail, the example of the embodiment is shown in the accompanying drawings, wherein from beginning Same or similar element or element with the same or similar functions are indicated to same or similar label eventually.Below by ginseng The embodiment for examining attached drawing description is exemplary, and for explaining only the invention, and is not construed as limiting the claims.
Those skilled in the art can understand that unless otherwise defined, all terms used herein (including skill Art term and scientific term) there is meaning identical with the general understanding of those of ordinary skill in fields of the present invention.Also It should be understood that those terms such as defined in the general dictionary should be understood that have in the context of the prior art The consistent meaning of meaning will not be explained in an idealized or overly formal meaning and unless defined as here.
Technical solution of the present invention is described in further detail with reference to the accompanying drawing:
A kind of pressure-measuring pipe road optimum design method for wind tunnel test of the present invention, as shown in Figure 1, specific method step is such as Under:
Step 1: it is theoretical based on circuited transmission, obtain the frequency response function in pressure-measuring pipe road.
Pressure-measuring pipe road system is subjected to blocking first, the stream of a unit is indicated with flow resistance, fluid capacitance, conductance and influenza Bulk properties.Fluid pipe road carries out dynamic analysis in frequency domain again, by continuity equation, the equation of momentum, energy equation and State equation obtains propagation operator and characteristic impedance under different flox conditions, and combines (the pipeline output of specific load impedance The ratio of end pressure and flow) and the source impedance ratio of end pressure and flow (pipeline input), it can will be at pressure-measuring pipe road system Manage into a pipeline fluid mechanics problem.
For simple cross-section pressure-measuring pipe road, the length is l for order, as shown in Figure 2, it is assumed that following three conditions are set up:
(1) air makees the Laminar Flow of microvariations in round tube;
(2) due to bore be much smaller than air-flow wavelength, pressure is propagated in pipe can be considered plane wave, axial heat transfer item and Radial heat transfer item is compared to negligible;
(3) for gas in pipe, pressure-measuring pipe is isothermal rigid walls pipe.
According to these it is assumed that the met differential equation of gas flowing can simplify in managing are as follows:
Continuity equation:
The equation of momentum:
Energy equation:
State equation:
Wherein, x is the distance at distance signal beginning, and r is pipe radius, and ρ is fluid density, and u is fluid velocity, and T is temperature Degree, p is Fluid pressure, and γ is adiabatic exponent, and v is kinematic viscosity, and σ is Pu Langte number;Subscript " 0 " indicates time average, ρ0 It is the time average of fluid density, T0It is the time average of temperature, v0It is the time average of kinematic viscosity, σ0It is general youth The time average of special number.
For cross-section round tube, since above-mentioned differential equation group considers heat transfer effect and viscosity simultaneously, therefore it is considered It is the accurate model for carrying out pipeline fluid dynamic analysis.It is transformed to above-mentioned equation group to solve:
Wherein:
Wherein, P (x, s) is the Laplace transform of pressure flow p (x, t);Q (x, s) is volume flow's Laplace transform;Subscript " 0 " indicates time average,0It is the time average of fluid density, A is circulation area of pipeline (r0 2), r0It is the time average of pipe inside radius, I1、I0It is single order and zeroth order Bessel function of imaginary argument respectively.v0It is kinematic viscosity Time average,0It is the time average of Prandtl number, is adiabatic exponent, a0It is the time average of the medium velocity of sound.Single order and Zeroth order Bessel function of imaginary argument are as follows:
By the met wave equation of the available flowing of formula (1-5):
In (1-9) formula:Title is propagation constant, the solution of equation (1-9) are as follows:
P (x, s)=C1e-χ(s)x+C2eχ(s)x (1-10)
In formula, right end two respectively indicate incidence wave and transmitted wave, are obtained by (1-9) formula:
C1,C2It for integral constant, is determined by boundary condition, simple cross-section straight tube such as Fig. 4 for being l to length, in input terminal Have:
X=0, P (x, s)=P1(s), Q (x, s)=Q1(s)
In formula, PiIt (s) is the Laplace transformation of pressure at two ends, QiIt (s) is Laplace transformation (the wherein i of both ends cross-sectional volume flow =1,2 respectively indicate the left and right ends of pipeline)
And then it obtains:
C1=[P1(s)+Zc(s)Q1(s)]/2,C1=[P1(s)-Zc(s)Q1(s)]/2
It is located at output end: x=l, P (x, s)=P2(s), Q (x, s)=Q2(s), then it is solved by formula (1-10), (1-11):
Above formula provides the matrix algebraic eqation relationship of cross-section straight tube inlet and outlet parameter.Wherein:Referred to as Characteristic impedance,Referred to as propagation operator.Notice that transfer function matrix is the function of s complex variable, takes S=2 π if, f are the maximum frequency for testing pulsating wind pressure signal of interest.
Common pressure-measuring pipe road is usually by 4 cross-section pipeline section groups of different length and different tube diameters in wind tunnel test At as shown in Figure 3.Connecting pipe then is pressed to 4 sections of common survey, is closed according to the transmitting that (1-12) formula obtains whole story end two o'clock parameter System are as follows:
M in formulai,i+1For the transfer matrix of i-th of pipeline section, M=M is enabled1,2...M4,5, indicate that pressure-measuring pipe road beginning terminal Parameters close The transfer matrix of system is 2 × 2 rank matrixes.End 5 connects modular pressure, it is believed that be closed end its flow is zero, therefore has Q5=0, Then the pressure dependence at tube body system whole story end is obtained are as follows:
P1=m11P5 (1-14)
In formula, m11For the top left hand element of Metzler matrix, it is the letter of pipeline geometric parameter and dielectric property and complex variable s Number, takes s=2 π if, and enable
It is spreaded to from above and constitutes the pipeline section number in pressure-measuring pipe road when being n, the end of n-th pipeline section connects modular pressure, there is Qn+1= 0, andThen P1=m11Pn+1,
Above formula, that is, pressure-measuring pipe road frequency response function.
Step 2: the Optimized model in multistage pressure-measuring pipe road is established, even model theory frequency response function described in step 1 is 1。
Actually this ideal pressure-measuring pipe road is since technical conditions are not present, but we can be to pressure-measuring pipe road In each section of pipe range optimization is adjusted with internal diameter, to obtain relatively flat frequency response function curve in a certain frequency range. At this point, we can think that the transmitting of pressure signal is not distorted, amendment can be not added in the frequency range and obtain Pressure signal at measuring point, therefore the time consumed by corrected Calculation can be saved.The frequency response function and composition in pressure-measuring pipe road are surveyed The length of each section of pipe of pressure pipeline is related to radius, therefore can construct objective function by the frequency response function of theory analysis:
Wherein, liFor the length of i-th of pipeline section, riFor the radius of i-th of pipeline section, i=1,2 ..., n, n are to constitute to survey pressure The pipeline section number of pipeline, f0It is the maximum frequency of test pulsating wind pressure signal of interest.
Optimization problem can be attributed to constrained nonlinear programming problem:
min f(l1,l2,...,ln,r1,r2,...rn) (2-2)
Constraint condition are as follows:
Wherein, L is pressure-measuring pipe road total length;ri u、ri lRespectively riUpper and lower limit.
Step 3: using pressure-measuring pipe road inside radius or duct length as optimized variable, being determined using genetic Optimization Algorithm and survey pressure The optimal design parameters of pipeline.
Now by taking pressure-measuring pipe road inside radius is optimized variable as an example, specific step is as follows for genetic algorithm:
It (1) is design variable by identified each pipe radius, it is random to generate initially generation individual;
(2) enabling pipeline total length is L, a length of l of metal tube1For constraint condition, calculating target function;
(3) judge fitness, execute genetic algorithm selection, intersect, four variation, migration operations, optimize, obtain new Individual;
(4) judge whether the evolutionary generation of initial individuals has reached preset value.If having reached preset value, stop excellent Change, exports the individual information in (3) step;If not up to preset value, (2) step is returned, 2-3 step is continued to execute, carries out Optimization.
The invention also discloses a kind of modification methods of the pressure signal of pressure-measuring pipe road acquisition, according to the pressure-measuring pipe road Above method optimization design obtains.
Judge whether to need further to correct the pressure signal of acquisition, the frequency in the pressure-measuring pipe road mainly obtained according to optimization Function H (f) is rung to be differentiated.The sample frequency that assumed stress scans valve is fs, in frequency f=[0, fs/ 2] in range, when frequency When ringing function H (f) ∝ 1, illustrate the obtained pressure-measuring pipe road system of optimization to fluctuation pressure effect of signals very little, can approximation ignore. Otherwise, the amplitude attenuation of fluctuation pressure signal be can not ignore, and to obtain high-precision pressure measurement signal, then be needed according to coefficient of wind pres Time-histories or Power spectral density are modified, to obtain true pressure signal.
Step 4: according to two evaluation indexes of coefficient of wind pres time-histories and Power spectral density, effect of optimization being commented Valence;
Step 5: determining final wind tunnel test pressure-measuring pipe road system and its high-precision modification method.
(1) coefficient of wind pres time-histories
It is anti-to release the true pressure of input according to the frequency response function H (f) and output pressure signal y (t) of known pressure-measuring pipe road system Force signal x (t), specific makeover process sum up as follows:
1. doing Fourier transform (FFT) to output signal y (t) obtains fourier coefficient sequence Y (f) in frequency domain;
2. finding out Fourier transform coefficient sequence X (f)=Y (f)/H (f) of the input signal in frequency domain;
3. carrying out inverse fourier transform to X (f), the time series x (t) of input signal is finally obtained.
(2) Power spectral density is corrected
If model surface measuring point pressure signal CPPower spectral density function beThe frequency response function in pressure-measuring pipe road is H (f), then real sensor measures fluctuation pressure signal CPxPower spectral density function are as follows:
By the pulsating wind pressure power spectral density function of the available pressure signal when not being distorted of above formula are as follows:
Directly had with the actual frequency response function of pipeline if not considering signal correction by the error caused by signal distortion It closes, theoretically mean square deviation evaluated error and Power estimation error caused by pipeline distortion can be corrected completely by formula (5-2).
Corresponding Power estimation error calculation formula are as follows:
Its mean square deviation evaluated error are as follows:
Specific embodiment
In the following with reference to the drawings and specific embodiments, the present invention is furture elucidated.It should be understood that following specific embodiments are only used In illustrating the present invention rather than limit the scope of the invention.
Following two operating condition is now taken to optimize pipe-line system:
(1) it when each one timing of length of tube, is optimized using each segment pipe inside radius as variable
Signal sampling frequencies needed for general wind tunnel test should take pulsating wind pressure signal in 100Hz~300Hz, this operating condition Frequency f0=150Hz (corresponding sample frequency is 300Hz).Pipeline total length L=1.2m is enabled, pipeline section number is 4, in addition to connecting mould Metal tube between type and PVC hose is regular length (l1), the pipe range perseverance of remaining each pipe takes 0.4m (l2=l3=l4= 0.4m).Using the inside radius of each segment pipe of pressure-measuring pipe road system as variable is calculated, with r0=[0.5,0.5,0.5,0.5] makees It is calculated for iterative initial value, as metal pipe range l1When=1cm, each pipe radius optimal solution for being calculated are as follows:
R=[0.47834,0.74267,0.34104,0.24562]
It is piping frequency response function curve when calculating variable before corresponding optimization and after optimization that Fig. 4, which is with pipeline inside radius,. Corresponding objective function maximum value is 0.79753%.
(2) it when each one timing of pipe inside radius, is optimized using each segment pipe length as variable
F is taken in this operating condition0=150Hz.Although using pipe radius as the available more satisfied optimum results of variable, The requirement of processing technology is allowed for, so that each pipe radius reaches optimal value and cannot achieve, each section of caliber can be enabled thus It is to calculate variable to optimize with duct length for fixed value.The caliber of each pipe takes r respectively1=0.35mm, r2=0.7mm, r3 =0.35mm, r4=0.25mm.As metal pipe range l1It is to calculate variable with its excess-three segment pipe length, with l when=1cm0= [0.4,0.4,0.4] is calculated as iterative initial value.Each length of tube optimal solution being calculated are as follows:
L=[l2,l3,l4]=[0.33674,0.25412,0.43996]
It is piping frequency response function curve when calculating variable before corresponding optimization and after optimization that Fig. 5, which is with length,.Accordingly Objective function maximum value is 0.56128%, it can be seen that the frequency width curve after pipe-line system is optimized is [0, f0] in close to flat.
(3) pressure signal time-histories is corrected
Four groups of different frequency sine pressure waves are inputted to the pressure-measuring pipe road system that length is 0.8m, utilize the frequency response letter of pipeline The pressure signal of several pairs of output ends is modified, and correction result is respectively as shown in (a) to (d) in Fig. 6, and wherein original signal is It is obtained by the collected pressure signal of very short duct institute.From figure it can be found that with frequency input signal increase, it is defeated The frequency response function that the attenuation degree of signal is significantly increased, and obtained using theoretical calculation out is modified distorted signal, so that it may To obtain with original signal very close to revise signal.
(a) to (d) in Fig. 7 be set forth different length pressure-measuring pipe road under log-log coordinate (0.5m, 0.8m, 1m and 1.2m) the collected pulsating wind pressure of institute uses Dimensionless Form from correction result, power spectrum is composed, and wherein dotted line is the pipeline | H(f)|2, it can be found that being modified available more true autopower spectral density letter to distorted signal using frequency response function Number curve.
Fig. 8 is length l=0.8m pressure-measuring pipe road spectrum error distribution schematic diagram certainly.It can be seen that being greater than the letter of 60Hz to frequency Number, Power estimation error can be greater than 50%, only signal of the frequency less than 20Hz, and Power estimation error just can be less than 10%.And it uses Pipeline frequency response function, which is directly corrected, can eliminate this error.
The technical means disclosed in the embodiments of the present invention is not limited to the technical means disclosed in the above technical means, and further includes Technical solution consisting of any combination of the above technical features.
The above, the only specific embodiment in the present invention, but scope of protection of the present invention is not limited thereto, appoints What is familiar with the people of the technology within the technical scope disclosed by the invention, it will be appreciated that expects transforms or replaces, and should all cover Within scope of the invention, therefore, the scope of protection of the invention shall be subject to the scope of protection specified in the patent claim.

Claims (6)

1. a kind of pressure-measuring pipe road optimum design method for wind tunnel test, which is characterized in that steps are as follows for specific method:
Step 1: it is theoretical based on circuited transmission, obtain the frequency response function in pressure-measuring pipe road;
Step 2: establishing the Optimized model in multistage pressure-measuring pipe road;
Step 3: using pressure-measuring pipe road inside radius or duct length as optimized variable, pressure-measuring pipe road being determined using genetic Optimization Algorithm Optimal design parameters.
2. a kind of pressure-measuring pipe road optimum design method for wind tunnel test according to claim 1, which is characterized in that step The frequency response function in pressure-measuring pipe road described in rapid 1 are as follows:
Wherein, m11It (f) is the top left hand element of the transfer matrix M of pressure-measuring pipe road whole story end parameters relationship, M=M1,2…Mn,n+1, n For the pipeline section number for constituting pressure-measuring pipe road, Mi,i+1For the transfer matrix of i-th of pipeline section, i=1,2 ..., n.
3. a kind of pressure-measuring pipe road optimum design method for wind tunnel test according to claim 1, which is characterized in that step The Optimized model in the multistage pressure-measuring pipe road established in rapid 2 specifically:
Objective function minf (l1,l2,...,ln,r1,r2,...rn);
Constraint condition
Wherein, i=1,2 ..., n, liFor the length of i-th of pipeline section, L is pressure-measuring pipe road total length, and n is to constitute pressure-measuring pipe road Pipeline section number, riFor the radius of i-th of pipeline section, ri u、ri lRespectively riUpper and lower limit,H (f) is the frequency response function in pressure-measuring pipe road, f0For pulsating wind pressure letter Number maximum frequency.
4. a kind of modification method of the pressure signal of pressure-measuring pipe road acquisition, the pressure-measuring pipe road is used as appointed in claims 1 to 3 Optimum design method described in one optimizes, which comprises the following steps:
Step A: frequency f=[0, fs/ 2] in range, if frequency response function H (f) ∝ 1 in pressure-measuring pipe road, the acquisition of pressure-measuring pipe road Pressure signal without amendment, it is no to then follow the steps 2;Wherein, fsFor the sample frequency of pressure scanning valve;
Step B: being modified according to the pressure signal that coefficient of wind pres time-histories or Power spectral density acquire pressure-measuring pipe road, Obtain true pressure signal.
5. a kind of modification method of the pressure signal of pressure-measuring pipe road acquisition according to claim 4, which is characterized in that step It is modified in B according to the pressure signal that coefficient of wind pres time-histories acquires pressure-measuring pipe road, obtains true pressure signal, specifically:
1. being Fourier transform FFT to the pressure signal y (t) of pressure-measuring pipe road acquisition obtains fourier coefficient sequence Y in frequency domain (f);
2. obtaining Fourier transform coefficient sequence X (f)=Y (f)/H (f) of the true pressure signal x (t) in frequency domain;
3. carrying out inverse fourier transform to X (f) obtained in 2., true pressure signal x (t) is finally obtained.
6. a kind of modification method of the pressure signal of pressure-measuring pipe road acquisition according to claim 4, which is characterized in that step It is modified in B according to the pressure signal that Power spectral density acquires pressure-measuring pipe road, obtains true pressure signal, specifically Are as follows:
In formula,For pressure-measuring pipe road acquisition pressure signal power spectral density function,For true pressure signal Power spectral density function.
CN201910342726.6A 2019-04-26 2019-04-26 A kind of pressure-measuring pipe road optimum design method for wind tunnel test Pending CN110146249A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910342726.6A CN110146249A (en) 2019-04-26 2019-04-26 A kind of pressure-measuring pipe road optimum design method for wind tunnel test

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910342726.6A CN110146249A (en) 2019-04-26 2019-04-26 A kind of pressure-measuring pipe road optimum design method for wind tunnel test

Publications (1)

Publication Number Publication Date
CN110146249A true CN110146249A (en) 2019-08-20

Family

ID=67594548

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910342726.6A Pending CN110146249A (en) 2019-04-26 2019-04-26 A kind of pressure-measuring pipe road optimum design method for wind tunnel test

Country Status (1)

Country Link
CN (1) CN110146249A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110826155A (en) * 2019-09-27 2020-02-21 哈尔滨工程大学 Multi-branch pipeline design method based on genetic algorithm
CN111737886A (en) * 2020-08-07 2020-10-02 中国空气动力研究与发展中心低速空气动力研究所 Wind tunnel test scheduling method and system
CN112197935A (en) * 2020-10-19 2021-01-08 重庆大学 Method for acquiring frequency response value at any frequency under any pipe length, pressure measuring pipeline correcting method and storage medium
CN112763181A (en) * 2020-12-29 2021-05-07 中国航空工业集团公司西安飞机设计研究所 Method for determining sampling parameters of pulsating pressure wind tunnel test signals
CN114910242A (en) * 2022-07-18 2022-08-16 中国空气动力研究与发展中心低速空气动力研究所 Pressure measurement device, pressure measurement method, pressure stability judgment method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104374468A (en) * 2014-10-24 2015-02-25 哈尔滨工业大学深圳研究生院 Intrinsic frequency measuring method for long and thin object
WO2015113844A1 (en) * 2014-01-31 2015-08-06 Andrew Wireless Systems Gmbh Method for compensating a temperature drift of a microwave filter
CN105488282A (en) * 2015-11-30 2016-04-13 南京航空航天大学 Cutting parameter segmentation and variable cutting depth optimizing method based on dynamic machining feature
CN107887903A (en) * 2017-10-31 2018-04-06 深圳供电局有限公司 Consider the micro-capacitance sensor robust Optimization Scheduling of element frequency characteristic

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015113844A1 (en) * 2014-01-31 2015-08-06 Andrew Wireless Systems Gmbh Method for compensating a temperature drift of a microwave filter
CN104374468A (en) * 2014-10-24 2015-02-25 哈尔滨工业大学深圳研究生院 Intrinsic frequency measuring method for long and thin object
CN105488282A (en) * 2015-11-30 2016-04-13 南京航空航天大学 Cutting parameter segmentation and variable cutting depth optimizing method based on dynamic machining feature
CN107887903A (en) * 2017-10-31 2018-04-06 深圳供电局有限公司 Consider the micro-capacitance sensor robust Optimization Scheduling of element frequency characteristic

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110826155A (en) * 2019-09-27 2020-02-21 哈尔滨工程大学 Multi-branch pipeline design method based on genetic algorithm
CN110826155B (en) * 2019-09-27 2022-05-20 哈尔滨工程大学 Multi-branch pipeline design method based on genetic algorithm
CN111737886A (en) * 2020-08-07 2020-10-02 中国空气动力研究与发展中心低速空气动力研究所 Wind tunnel test scheduling method and system
CN112197935A (en) * 2020-10-19 2021-01-08 重庆大学 Method for acquiring frequency response value at any frequency under any pipe length, pressure measuring pipeline correcting method and storage medium
CN112763181A (en) * 2020-12-29 2021-05-07 中国航空工业集团公司西安飞机设计研究所 Method for determining sampling parameters of pulsating pressure wind tunnel test signals
CN114910242A (en) * 2022-07-18 2022-08-16 中国空气动力研究与发展中心低速空气动力研究所 Pressure measurement device, pressure measurement method, pressure stability judgment method and system

Similar Documents

Publication Publication Date Title
CN110146249A (en) A kind of pressure-measuring pipe road optimum design method for wind tunnel test
US8910530B2 (en) Method and apparatus for the measurement of flow in gas or oil pipes
CN110197049B (en) Non-metal pipeline leakage positioning method based on transient inverse problem
CN110513603B (en) Non-metal pipeline leakage positioning method based on inverse transient analysis method
CN110645483B (en) Urban buried pipeline early leakage diagnosis method based on spectrum analysis
CN103742794A (en) Simulating device and simulating method for pipeline leakage acoustic emission signals
CN107085123A (en) A kind of multipath pipe flow speed measurement apparatus and its method based on acoustic technique
Li et al. Propagation of monopole source excited acoustic waves in a cylindrical high-density polyethylene pipeline
CN108931270B (en) Two-phase flow parameter detection method based on porous throttling and acoustic emission technology
CN110231408A (en) A kind of method and device measuring material acoustics constant
CN106018561A (en) System and method for measuring sound wave amplitude attenuation coefficients in different pipeline structures
Lopes et al. Development of 3-D ultrasonic anemometer with nonorthogonal geometry for the determination of high-intensity winds
CN103438931A (en) Integrated moist steam flow dryness measuring device and measuring method
Wang et al. A new method for processing ultrasonic gas flowmeter signal in wet gas
Kupnik et al. Numerical simulation of ultrasonic transit-time flowmeter performance in high temperature gas flows
Ferrante et al. Experimental characterization of PVC-O pipes for transient modeling
Kim Advanced numerical and experimental transient modelling of water and gas pipeline flows incorporating distributed and local effects.
CN105868491B (en) A kind of foreseeable method of plate sound insulation property based on Least squares approach theory
CN110631958A (en) Gas-liquid two-phase flow parameter measuring method
Gao et al. Application of state vector formalism and Legendre polynomial hybrid method in the longitudinal guided wave propagation analysis of composite multi-layered pipes
CN114459674B (en) Dynamic calibration method and system for amplitude-frequency characteristics of high-frequency pressure sensor
Yu et al. The inverse of material properties of functionally graded pipes using the dispersion of guided waves and an artificial neural network
Kawell et al. A theoretical investigation of the nonlinear acoustic response of abrupt area contractions
CN110260948B (en) Liquid level measuring method based on acoustic resonance frequency nonlinear filtering
CN113865787A (en) Time domain calibration method for pressure measurement pipeline signal

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: 20190820