CN115493800B - Synchronous parallel acquisition system for steady-state pressure and pulsating pressure data and application method - Google Patents
Synchronous parallel acquisition system for steady-state pressure and pulsating pressure data and application method Download PDFInfo
- Publication number
- CN115493800B CN115493800B CN202211437857.0A CN202211437857A CN115493800B CN 115493800 B CN115493800 B CN 115493800B CN 202211437857 A CN202211437857 A CN 202211437857A CN 115493800 B CN115493800 B CN 115493800B
- Authority
- CN
- China
- Prior art keywords
- pressure
- acquisition
- steady
- pulsating
- formula
- 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.)
- Active
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000009530 blood pressure measurement Methods 0.000 claims abstract description 74
- 238000012545 processing Methods 0.000 claims abstract description 30
- 230000035485 pulse pressure Effects 0.000 claims description 20
- 230000006870 function Effects 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 16
- 238000005311 autocorrelation function Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000005314 correlation function Methods 0.000 claims description 3
- 230000010349 pulsation Effects 0.000 claims description 3
- 230000001960 triggered effect Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 abstract description 14
- 238000013461 design Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 239000000306 component Substances 0.000 description 6
- 230000003750 conditioning effect Effects 0.000 description 4
- 238000007726 management method Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000712 assembly Effects 0.000 description 2
- 238000000429 assembly Methods 0.000 description 2
- 238000010219 correlation analysis Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241000135164 Timea Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L11/00—Measuring steady or quasi-steady pressure of a fluid or a fluent solid material by means not provided for in group G01L7/00 or G01L9/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
The invention discloses a synchronous parallel acquisition system and an application method for steady-state pressure and pulsating pressure data, which relate to the field of wind tunnel tests.A separate and same-type AD (analog-digital) digitizer is used by each acquisition unit in each acquisition terminal in a steady-state pressure measurement subsystem and a pulsating pressure measurement subsystem, phase locking is carried out by a phase-locked loop, and a time stamp counter is integrated in a processor of the AD digitizer; the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem provide synchronous clocks through an external synchronous control module; and the acquisition terminals of the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem are communicated with the central processing unit through the matched synchronous triggers. The invention provides a synchronous parallel acquisition system and an application method for steady-state pressure and pulsating pressure data, which realize synchronous parallel acquisition of the steady-state pressure data and the pulsating pressure data in a wind tunnel test so as to ensure the accuracy and time domain correlation of subsequent data processing.
Description
Technical Field
The invention relates to the field of wind tunnel tests. More specifically, the invention relates to a system for synchronously and parallelly acquiring steady-state pressure (data) and pulsating pressure (dynamic) data in a wind tunnel test and an application method thereof.
Background
The types of data required to be collected in the wind tunnel test comprise a plurality of types, and the model surface pressure and the internal cavity pulsating pressure are two very important data.
The measurement of the pressure distribution on the surface of the model is an important means for knowing the flow characteristics of the surrounding model, determining the positions of the maximum pressure point, the minimum pressure point, the separation point and the shock wave position, and determining the positions of the normal force, the pressure difference resistance and the pressure center point of the model, so that the aerodynamic characteristics of the aircraft are researched, the aerodynamic load of the aircraft is obtained, and load data are provided for the structural strength calculation of the aircraft. The conventional test method for the current pressure distribution is to open a plurality of pressure measurement spaces on the surface of a model, wherein the pressure at the pressure measurement holes is different from dozens of pressure measurement spaces to hundreds of pressure measurement spaces, the pressure at the pressure measurement holes is connected to an electronic pressure scanning valve through a ventilation conduit, and the pressure distribution data of the surface of the model is obtained through a pressure sensor and an analog-to-digital (A/D) conversion module of the electronic pressure scanning valve, as shown in figure 1, the electronic pressure scanning valve used for measuring the pressure data of the surface of the model is a data acquisition mode which simultaneously acquires dozens of pressure points to hundreds of pressure points and has a low sampling rate, and the limit sampling rate is generally within 1 KHz.
The pressure that changes rapidly and randomly over time is called the pulsating pressure. In the wind tunnel test, many complex flow phenomena, such as vortex, separation flow, turbulent flow and the like belong to unsteady flow, quite strong pulsating pressure is generated, and unsteady characteristics of the complex flow phenomena can be revealed through measurement of the pulsating pressure. It cannot be accurately described by mathematical expressions, has no instantaneous rules, and can be processed and analyzed only according to the statistical rules. The strong pulsating pressure can excite the structure of the aircraft to vibrate or cause a severe noise environment in the aircraft, and in addition, the surface thermal conductivity of the hypersonic aircraft can be influenced, so that the difficulty is brought to the heat-proof design of the high-speed aircraft. As shown in fig. 2, the method of the pulse pressure test is to mount a pulse pressure sensor on the model, the sensor converts the sensed pulse pressure signal into an electrical signal, and the electrical signal is filtered and amplified, and recorded, stored and processed by a data acquisition system. The dynamic data acquisition system used for pulse pressure data acquisition is an acquisition mode which aims at a small number of sensors, generally does not exceed dozens of sensors, but has a high sampling rate, and the sampling rate can reach 1MHz.
In the current wind tunnel test, the measurement of the surface pressure of the model and the measurement of the pulsating pressure are realized by utilizing respective special equipment, namely an electronic pressure scanning valve and a dynamic data acquisition system, to independently complete data acquisition. The acquisition mode has the following four problems:
the electronic pressure scanning valve system for measuring the surface pressure of the model is provided with a plurality of acquisition terminals, each acquisition terminal is provided with a dozen to dozens of acquisition units, and pressure measurement can be simultaneously carried out on nearly a thousand pressure measurement points. All pressure acquisition units on one terminal share one AD (analog-digital) instrument, and when pressure data are acquired, only a serial mode is adopted, and parallel acquisition of multiple paths of pressure data cannot be realized;
secondly, as shown in fig. 3, because tens of pressure acquisition units of the electronic pressure scanning valve share one AD digitizer, the data acquisition mode can only adopt a polling mode, that is, all the pressure acquisition units of the acquisition terminal record pressure data one by one according to the sequence, the data between two adjacent acquisition units has a delay Δ t in time, the more the pressure acquisition units of the terminal are, the greater the delay of the pressure recording time between the first acquisition unit and the last acquisition unit is (n-1) × Δ t, and n is the number of the acquisition units on the terminal. Meanwhile, in order to reduce the random variation error of the pressure data, a plurality of pressure values of the same acquisition unit are averaged to be used as the finally acquired pressure value, and the delay of the pressure recording time between the first acquisition unit and the last acquisition unit is (nk-1) × Δ t. In addition, in order to ensure the real-time performance of data reading and processing, the sampling rate of the electronic pressure measurement is set to be small, and is generally about 100 Hz. At the moment, the data recording time between the pressure acquisition units has an error of a second level, and because the measured pressure in the wind tunnel test is a fluctuation value, the finally acquired pressure data is not a pressure value at the same moment, and a larger error is caused by the delay of the acquisition time.
Thirdly, as shown in fig. 4, the collection trigger signals among the pressure measurement terminals have a network transmission mode and also have a cable transmission mode. Because the delay of the transmission time and the different phases of the clock signals between different acquisition terminals exist, a random acquisition time delay exists between different terminals, and errors are brought to final data.
Fourthly, as shown in fig. 5, because the model surface pressure data acquisition device and the pulsating pressure data acquisition device are independent from each other, the respective acquired data lack a corresponding relationship in time, and the final data analysis only can use the model surface pressure data as a whole reference, and cannot perform time domain correlation analysis on the model surface pressure change and the pulsating pressure change within one wind tunnel test time, and the time domain correlation analysis between the two is a key for extracting the interactive coupling characteristics of the two.
Disclosure of Invention
An object of the present invention is to solve at least the above problems and/or disadvantages and to provide at least the advantages described hereinafter.
In order to achieve the purposes and other advantages of the invention, the invention provides a synchronous parallel acquisition system of steady-state pressure and pulse pressure data, which comprises a steady-state pressure measurement subsystem and a pulse pressure measurement subsystem which respectively acquire steady-state pressure and pulse pressure, and is characterized in that at least one acquisition terminal is respectively arranged in the steady-state pressure measurement subsystem and the pulse pressure measurement subsystem, each acquisition unit in each acquisition terminal uses an independent and same-type AD (analog-to-digital) digital instrument and carries out phase locking through a phase-locked loop, and a processor of the AD digital instrument is integrated with a timestamp counter;
the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem provide a phase-locked synchronous clock for each acquisition terminal through matched external synchronous control modules;
and the acquisition terminals of the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem are communicated with the central processing unit through the matched synchronous triggers.
An application method of a synchronous parallel acquisition system for steady-state pressure data and pulsating pressure data applied to a wind tunnel test comprises the following steps:
s10, the synchronous control module provides a phase-locked common clock for each acquisition terminal of the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem so as to ensure that each acquisition unit AD digital instrument is in the same state to be triggered;
s11, the central processing unit sends a trigger instruction to each acquisition terminal through the synchronous trigger, each acquisition terminal sets an acquisition starting point for the AD digitizer of each acquisition unit after receiving the trigger instruction, pressure signal acquisition is started, and the clock frequency of each AD digitizer is locked through a phase-locked loop in the acquisition process;
pressure signals acquired by the AD digital instruments are recorded with the sampling period number of each AD digital instrument before the current synchronous trigger event occurs through a timestamp counter, and trigger time errors among different acquisition terminals are identified and corrected, so that synchronous acquisition of signals in each subsystem is realized;
and S12, the acquisition terminals of the subsystems transmit the acquired pressure signals back to the central processing unit through the concentrator, and the two data are aligned at the same time mark by comparing the time marks of the steady-state pressure data and the pulsating pressure data among the subsystems.
Preferably, the method further comprises the following steps:
s13, the central processing unit carries out flow field structure change and background noise removal on data measured by the steady-state pressure measurement subsystem, and carries out pretreatment on pulsating pressure data measured by the pulsating pressure measurement subsystem;
wherein the preprocessing is configured to analyze and process in three aspects of amplitude domain, time domain and frequency domain.
Preferably, in the amplitude domain processing, the root mean square value P is determined by the pulsating pressure rms Sound pressure level SPL, coefficient of pulsating pressure C prms Describing the total intensity of the pulsating pressure in the measured frequency range, and describing the probability of the instantaneous value of the pulsating pressure falling in a certain specified pressure range through a probability density function PDF;
wherein, the P rms The corresponding formula one is:
in formula one, P (T) represents a pressure signal that varies randomly with time, and T represents a measurement time;
the formula two corresponding to the SPL is as follows:
in formula two, P ref Is a reference pressure;
said C is prms The corresponding formula three is:
in the third formula, q represents the incoming flow velocity pressure;
the formula four corresponding to the PDF is as follows:
in equation four, Δ p represents a pressure increase; t is p Representing the total time during which the instantaneous value of the pulsating pressure falls within the range p to p + deltap during the measurement time T.
Preferably, in the time domain processing, the autocorrelation function R (b) is obtained by) Or r: () Waveform describing a pulsating pressure signal and its time delayThe degree of similarity of the waveforms in time by the cross-correlation function R (,) Or r: (,) Describing the similarity degree between the pulse pressure signals of two different measuring points;
in the formula six, R (0) represents the value of the autocorrelation function at a time delay of 0;
in formula seven, p x (t) represents the pulsating pressure signal at point x,indicating the distance from the point xAt another measuring point of the station is delayed in timeA time-dependent pulsating pressure signal;
in the formula eight, R x (0) Andrespectively represent measuring points x and x +Delay in time of autocorrelation function of the pulsating pressure signalValue at = 0.
Preferably, in the frequency domain processing, the distribution of the pulsating pressure energy with frequency is described by a power spectral density function G (f);
wherein, the formula nine corresponding to G (f) is:
in equation nine, f represents frequency and Δ f represents a frequency interval.
The invention at least comprises the following beneficial effects: the system can realize synchronous acquisition of steady-state pressure data and pulsating pressure data of different acquisition systems in a wind tunnel test so as to ensure the accuracy and time domain correlation of subsequent data processing.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention.
Drawings
FIG. 1 is a schematic diagram of a prior art model surface pressure data acquisition approach;
FIG. 2 is a diagram illustrating a method for acquiring pulsating pressure data according to the prior art;
FIG. 3 is a schematic diagram of the error of the electronic pressure measurement data acquisition delay in the prior art
FIG. 4 is a schematic diagram of the triggering time difference between different acquisition terminals for pressure measurement in the prior art
FIG. 5 is a schematic diagram showing the lack of time-domain correlation between the surface pressure of a model and the collected data of the pulsating pressure in the prior art;
FIG. 6 is a schematic diagram of synchronous data acquisition and preprocessing steps of a steady-state flow field and pulsating pressure in the preprocessing process of the present invention;
FIG. 7 is a block diagram of the system of the present invention;
FIG. 8 is a block diagram of the components of one of the subsystems of the present invention;
FIG. 9 is a schematic diagram of the present invention for synchronously collecting data of different frequencies;
the system comprises a steady-state pressure measurement subsystem-1, a pulsating pressure measurement subsystem-2, an acquisition terminal-10, an acquisition unit-11, an AD (analog-to-digital) instrument-12, a phase-locked loop-13, an external synchronous control module-3 and a central processing unit-4.
Detailed Description
The present invention is further described in detail below with reference to the attached drawings so that those skilled in the art can implement the invention by referring to the description text.
It will be understood that terms such as "having," "including," and "comprising," as used herein, do not preclude the presence or addition of one or more other elements or groups thereof.
As shown in fig. 7-8, a system for synchronously acquiring steady-state pressure and pulse pressure data comprises a steady-state pressure measurement subsystem 1 and a pulse pressure measurement subsystem 2 which are used for respectively acquiring steady-state pressure on the surface of a model and pulse pressure at a special part, wherein at least one acquisition terminal 10 is respectively arranged in the steady-state pressure measurement subsystem and the pulse pressure measurement subsystem, each acquisition unit 11 in each acquisition terminal uses an independent and same-type AD digital instrument 12 and is phase-locked by a phase-locked loop 13, and a time stamp counter is integrated in a processor of the AD digital instrument;
the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem provide a phase-locked synchronous clock for each acquisition terminal through the matched external synchronous control module 3;
the steady state pressure measurement subsystem, each acquisition terminal of the pulsating pressure measurement subsystem communicate with central processing unit 4 through the synchronous trigger cooperating, the effect of the synchronous trigger, in the multisystem, provide the triggering signal and insert the time stamp mark through the external synchronous control module, realize the phase synchronization of each acquisition terminal among the systems, and to the independent single system inside, each measurement unit in each acquisition terminal all adopts independent, the A/D digital instrument of the same type, realize the phase synchronization among each measurement terminal and measurement unit through synchronous clock and synchronous bus, namely the realization mode that the system of the invention collects synchronously is:
firstly, an independent AD (analog-to-digital) digital instrument with the same type is used by each acquisition unit in each acquisition terminal in a model surface pressure measurement system and a special part pulsating pressure measurement system, and phase locking is carried out through a phase-locked loop, so that the clock frequencies of all the AD digital instruments of the acquisition units are completely consistent, and the phases do not deviate along with time;
secondly, providing a phase-locked synchronous clock for each acquisition terminal of the model surface pressure measurement system and the special part pulsating pressure measurement system through an external synchronous control module to ensure that the AD digital instrument of each acquisition unit in each acquisition terminal of the two systems is in the same state to be triggered;
thirdly, providing a synchronous trigger for each acquisition terminal of the two systems through an external synchronous control module to determine the starting point of the data recorded by the AD digitizer of each acquisition unit in each acquisition terminal;
and fourthly, introducing a time stamp counter into each AD digital instrument, wherein the time stamp counter records the sampling period number of the AD digital instrument in each acquisition unit before the synchronous trigger event occurs, and is used for identifying and correcting the trigger time error among different acquisition units.
And fifthly, comparing the time marks of the model surface pressure data and the pulse pressure data, and aligning the two data at the same time mark. As shown in fig. 9, by adopting the scheme of the present invention, synchronous data acquisition of different frequencies can be effectively realized, and measurement errors are reduced, so that the measurement structure meets the actual test requirements.
In a specific implementation of the steady-state pressure measurement subsystem, each acquisition unit is configured to include:
the pressure sensors are matched with the positions of the points to be measured on the gas path connecting assembly;
the signal conditioning circuit is in communication connection with each pressure sensor to construct a corresponding measuring channel;
the pressure measurement assembly is in communication connection with each first signal conditioning circuit so as to realize synchronous acquisition of each pressure sensor, and the modular design is adopted, so that the interface is expandable and the overall technical index of the system is not influenced;
each acquisition terminal is configured to include:
the pressure measurement machine box (not shown) is used for integrating pressure measurement components, a first controller with an operating system is arranged in the pressure measurement machine box, the pressure measurement components are inserted into the pressure measurement machine box to complete modular integration of pressure measurement, the machine box is a 4U standard rack machine box and is provided with a core controller slot, at most 8 acquisition board cards can be inserted, the first machine box plug-in card adopts a quick-insertion structure, a panel is provided with a locking screw, the board cards can be reinforced to be installed after being screwed, the core controller is a core component of a steady-state pressure measurement system and is equivalent to a CPU and a mainboard of acquisition equipment, and functions of data acquisition management, acquisition instruction sending, data transmission and the like can be completed. The operating system of the core controller is an RT real-time operating system. When the real-time operating system interacts with external data, the real-time operating system can process the data at a high enough speed, the processing result can control the production process or make a quick response to the processing system within a specified time, all available resources are scheduled to complete real-time tasks, and all real-time tasks are controlled to run in a coordinated and consistent mode. Compared with a time-sharing operating system, the system has the advantages of timely response and high reliability. For the data acquisition processing server, the synchronous acquisition component is like a black box, a preset program in the core controller automatically calls and operates a synchronous acquisition function, an acquisition instruction sent by the server is received by adopting a driving program, and acquired data are sent to the data acquisition processing server;
the portable pressure controller mainly outputs a standard pressure to provide a standard pressure source so as to conveniently finish the calibration and troubleshooting work of the pressure measurement module. In the scheme, the parameter indexes of the pressure measurement components need to meet the requirements of a table 1, and the parameter indexes of the pulsating pressure measurement subsystem need to meet the requirements of a table 2, so that the capability of phase synchronization and high-precision acquisition between measurement points and between the pressure measurement components and pulsating pressure data acquisition is achieved;
TABLE 1
TABLE 2
Each pressure measurement assembly is configured to include:
the first FPGA is used for processing and outputting the signals input by the signal conditioning circuit;
the first phase-locked loop PLL is in communication connection with the first FPGA, the first clock management module and the CPU;
the external communication interface, the synchronous clock signal and the synchronous trigger signal are in communication connection with the first clock management module and the CPU through a control bus;
the central processing unit (also referred to as a control unit) is configured to include:
the rack is provided with a front panel with a signal interface;
a second FPGA30 disposed inside the chassis;
the second phase-locked loop PLL is in communication connection with the second FPGA, the second clock management module and the synchronous signal generator;
the control unit is in communication connection with the pressure measurement assemblies and the dynamic data acquisition assemblies through a shunt, the overall structural design of the control unit adopts a 1U rack type structure, a signal interface is configured on a front panel and can be accessed to input and output signals, a synchronization module is in communication with an upper computer through a USB port, hardware equipment in the system adopts an embedded type structure based on FPGA, and the hardware structure based on FPGA and AD has high real-time performance meeting test conditions, so that the acquisition and storage of signals can be realized, and the advantages of strong synchronism, low power consumption, high reliability and the like are achieved; the control unit is used as an external trigger source and synchronously accessed to a plurality of acquisition terminals through a splitter to realize the layout of synchronous triggering of a plurality of devices;
each measuring unit of the scheme adopts a modular design, and the modules have independent functions and are not coupled. The modularized design is easy to expand the system capacity and improve the maintenance efficiency.
The FGPA realizes the operations of data acquisition, data processing, data transmission and the like in the FPGA through a standard PCIe interface. The method is divided into the following steps according to functions: the data exchange system comprises an interface unit, a processing unit, an output unit and the like, wherein a data flow driving mode is adopted among the units, and after the data processing of the unit is completed, the data are packaged and sent to the next level, so that the data exchange is realized inside the framework.
The FPGA design elements comprise interface design, clock design, reset design, function design and the like. The interface design criteria are: excessive logic is not added, and the influence on the interface time sequence after the added logic is jammed is avoided; clocks include logic clocks, interface clocks, memory clocks, and the like. The logic clock depends on the key path of the logic, and the product performance is improved. When the FPGA realizes the synchronous signal time sequence, the synchronous adoption of the interface is realized by adopting a fixed and accurate interface clock. The memory clock realizes the data synchronous cache, and the data loss or instability caused by the refresh frequency is avoided during the design; the internal reset of the FPGA comprises hard reset, logic reset, soft reset and the like. The hard reset is introduced into the external pin and is fed in when the power is on, so that the whole FPGA logic configuration reaches a stable state after being completed. The logic reset is generated by the internal logic of the FPGA and is used for setting the preparation state of signals such as synchronous control and the like. The soft reset is used in the debugging stage, and the soft reset is inserted into the point to be tested, the fault locating point and the like, so that the problem of segmentation can be quickly located, and the debugging speed is accelerated. The function design comprises the functions of signal acquisition, signal conditioning, synchronous control, data transmission and the like through the FPGA. The principle of functional design is stable and efficient, and the required circuit realization function and the constraint conditions for realizing the circuit, such as speed, power consumption, circuit type and the like, are comprehensively considered.
And the further pulsation pressure measurement subsystem is used for measuring dynamically changed data, and the acquisition terminal and the corresponding acquisition unit of the further pulsation pressure measurement subsystem are similar to those of the steady-state pressure measurement subsystem except for the difference of the types of the sensors.
The pressure which changes randomly along with time is called pulsating pressure, many complex flow phenomena such as vortex, separation flow, turbulent flow and the like belong to unsteady flow, quite strong pulsating pressure is generated, and unsteady characteristics of the complex flow phenomena can be revealed through pulsating pressure measurement. It cannot be accurately described by mathematical expressions, has no instantaneous rules, and is only processed and analyzed according to the statistical rules. As shown in FIG. 6, the invention performs the pre-processing of the pulsating pressure data from three aspects of amplitude domain, time domain and frequency domain. The processed data provides basis for extracting the characteristics of the pulsating pressure change brought by the state change of the model and the state change of the flow field, and particularly provides data support for the correlation characteristic analysis of the static pressure change and the pulsating pressure change in the time domain.
1. In the amplitude domain processing method, the root mean square value P is obtained through the pulsating pressure in the amplitude domain processing rms Sound pressure level SPL, coefficient of pulsating pressure C prms Describing the total intensity of the pulsating pressure in the measured frequency range, and describing the probability of the instantaneous value of the pulsating pressure falling in a certain specified pressure range through a probability density function PDF;
wherein, the P rms The corresponding formula one is:
in formula one, P (T) represents a pressure signal that varies randomly with time, and T represents a measurement time;
the formula two corresponding to the SPL is as follows:
in formula two, P ref Is a reference pressure;
said C is prms The corresponding formula three is:
in the third formula, q represents the incoming flow velocity pressure;
the formula four corresponding to the PDF is as follows:
in equation four, Δ p represents a pressure increase; t is p Representing the total time during which the instantaneous value of the pulsating pressure falls within the range p to p + deltap during the measurement time T.
2. The time domain processing method comprises the following main statistical functions of the time domain:
(1) An autocorrelation function R: () Describing the waveform and time delay of a pulsating pressure signalThe degree of similarity of the time waveforms, sometimes by the autocorrelation coefficient r (m:)) And (4) showing.
In the formula:represents the time delay, and R (0) represents the value of the autocorrelation function at a time delay of 0.
In the formula:a pulsating pressure signal representing the measuring point x;indicating the distance from the point xAt another measuring point of the station is delayed in timeThe pulsating pressure signal of time.
The cross-correlation coefficient describes the similarity degree between pulse pressure signals of two different measuring points or is called the correlation degree, and sometimes the cross-correlation coefficient r (R: (R)), ) Represents:
R x (0) Andrespectively represent measuring points x and x +Delay in time of autocorrelation function of the pulsating pressure signalValue at = 0.
3. In the frequency domain processing method, the main statistical function of the frequency domain is mainly to analyze the pulsating pressure signal in the frequency domain, which is generally called spectral analysis, and the analysis result is spectral lines and curves of various physical quantities with the abscissa as the frequency.
The power spectral density function G (f), describes the distribution of the pulsating pressure energy with frequency.
In the formula: f denotes frequency, and Δ f denotes a frequency interval.
The above scheme is merely illustrative of a preferred example, and is not limiting. When the invention is implemented, appropriate replacement and/or modification can be carried out according to the requirements of users.
The number of apparatuses and the scale of the process described herein are intended to simplify the description of the present invention. Applications, modifications and variations of the present invention will be apparent to those skilled in the art.
While embodiments of the invention have been disclosed above, it is not intended to be limited to the uses set forth in the specification and examples. It can be applied to all kinds of fields suitable for the present invention. Additional modifications will readily occur to those skilled in the art. Therefore, the invention is not to be limited to the specific details and illustrations shown and described herein, without departing from the general concept as defined by the claims and their equivalents.
Claims (3)
1. A synchronous parallel acquisition system of steady-state pressure and pulsating pressure data comprises a steady-state pressure measurement subsystem and a pulsating pressure measurement subsystem which are used for respectively acquiring steady-state pressure and pulsating pressure, and is characterized in that the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem are respectively provided with at least one acquisition terminal, each acquisition unit in each acquisition terminal uses an independent and same-type AD (analog-to-digital) digitizer and is phase-locked by a phase-locked loop, and a processor of the AD digitizer is integrated with a timestamp counter;
the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem provide a phase-locked synchronous clock for each acquisition terminal through matched external synchronous control modules;
and the acquisition terminals of the steady-state pressure measurement subsystem and the pulsation pressure measurement subsystem are communicated with the central processing unit through the matched synchronous triggers.
2. A method for using the synchronous and parallel acquisition system of the steady-state pressure and the pulsating pressure data as claimed in claim 1, which comprises:
s10, the external synchronous control module provides a phase-locked common clock for each acquisition terminal of the steady-state pressure measurement subsystem and the pulsating pressure measurement subsystem so as to ensure that each acquisition unit AD digital instrument is in the same state to be triggered;
s11, the central processing unit sends a trigger instruction to each acquisition terminal through the synchronous trigger, each acquisition terminal sets an acquisition starting point for the AD digitizer of each acquisition unit after receiving the trigger instruction, pressure signal acquisition is started, and the clock frequency of each AD digitizer is locked through a phase-locked loop in the acquisition process;
the pressure signals acquired by each AD digital instrument record the sampling period number of each AD digital instrument before the current synchronous trigger event occurs through a timestamp counter, and trigger time errors among different acquisition terminals are identified and corrected, so that synchronous acquisition of signals in each subsystem is realized;
and S12, the acquisition terminals of the subsystems transmit the acquired pressure signals back to the central processing unit through the concentrator, and the two data are aligned at the same time mark position by comparing the time marks of the steady-state pressure data and the pulsating pressure data among the subsystems.
3. The method for applying the synchronous and parallel acquisition system of the steady-state pressure and the pulse pressure data as claimed in claim 2, further comprising:
s13, the central processing unit carries out flow field structure change and background noise removal on data measured by the steady-state pressure measurement subsystem, and carries out pretreatment on pulsating pressure data measured by the pulsating pressure measurement subsystem;
wherein the preprocessing is configured to analyze and process in three aspects of amplitude domain, time domain and frequency domain;
in amplitude domain processing, the root mean square value P is obtained by the pulsating pressure rms Sound pressure level SPL, coefficient of pulsating pressure C prms Describing the total intensity of the pulsating pressure in the measured frequency range, and describing the probability of the instantaneous value of the pulsating pressure falling in a certain specified pressure range through a probability density function PDF;
wherein, the P rms The corresponding formula one is:
in formula one, P (T) represents a pressure signal that varies randomly with time, and T represents a measurement time;
the formula two corresponding to the SPL is as follows:
in formula two, P ref Is a reference pressure;
said C is prms The corresponding formula three is:
in the formula III, q represents the incoming flow velocity pressure;
the formula four corresponding to the PDF is as follows:
in equation four, Δ p represents a pressure increase; t is p Represents the total time during which the instantaneous value of the pulsating pressure falls within the range of p to p + Δ p during the measurement time T;
in the time domain processing, the similarity degree of the waveform of a pulse pressure signal and the waveform of the pulse pressure signal at the time delay tau is described through an autocorrelation function R (tau) or R (tau), and the similarity degree between pulse pressure signals at two different measuring points is described through a cross-correlation function R (epsilon, tau) or R (epsilon, tau);
in formula five, τ represents a time delay;
the formula six corresponding to r (τ) is:
in formula six, R (0) represents the value of the autocorrelation function at a time delay of 0;
the formula seven corresponding to the R (epsilon, tau) is as follows:
in formula seven, p x (t) represents the pulsating pressure signal at point x,representing the pulsating pressure signal at a time delay tau of another measuring point which is at a distance epsilon from the measuring point x;
the formula eight corresponding to r (epsilon, tau) is:
in the formula eight, R x (0) Andvalues of autocorrelation functions of the pulsating pressure signals at the measuring points x and x + τ, respectively, at a time delay τ = 0;
in the frequency domain processing, the distribution of the pulse pressure energy along with the frequency is described by a power spectral density function G (f);
wherein, the formula nine corresponding to G (f) is:
in equation nine, f represents frequency and Δ f represents a frequency interval.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211437857.0A CN115493800B (en) | 2022-11-17 | 2022-11-17 | Synchronous parallel acquisition system for steady-state pressure and pulsating pressure data and application method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211437857.0A CN115493800B (en) | 2022-11-17 | 2022-11-17 | Synchronous parallel acquisition system for steady-state pressure and pulsating pressure data and application method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115493800A CN115493800A (en) | 2022-12-20 |
CN115493800B true CN115493800B (en) | 2023-02-28 |
Family
ID=85115982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211437857.0A Active CN115493800B (en) | 2022-11-17 | 2022-11-17 | Synchronous parallel acquisition system for steady-state pressure and pulsating pressure data and application method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115493800B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118193317A (en) * | 2024-03-13 | 2024-06-14 | 深圳市蓝色亿点科技有限公司 | Method and system for collecting server environment data |
CN118095143A (en) * | 2024-04-23 | 2024-05-28 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Rapid optimization method for combustion stability of combustion chamber |
CN118644052B (en) * | 2024-08-14 | 2024-10-29 | 克拉玛依创业有限公司 | Oil reservoir drilling real-time operation monitoring system and method |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5756892A (en) * | 1997-02-27 | 1998-05-26 | The United States Of America As Represented By The United States National Aeronautics And Space Administration | Apparatus for measuring ambient pressure within a gaseous flow field |
JP2005354568A (en) * | 2004-06-14 | 2005-12-22 | Sony Corp | Physical information acquisition method, physical information acquisition device and semiconductor device for detecting physical value distribution |
CN101534021A (en) * | 2009-04-27 | 2009-09-16 | 北京四方继保自动化股份有限公司 | Multimode data acquisitions and processing method applied to power automation system |
CN104113342A (en) * | 2013-11-28 | 2014-10-22 | 西安电子科技大学 | High-speed data synchronous circuit used for high-speed digital-to-analog converter |
CN107402597A (en) * | 2017-07-31 | 2017-11-28 | 上海联影医疗科技有限公司 | A kind of data and method, apparatus of clock alignment, medium and magnetic resonance equipment |
CN107656151A (en) * | 2017-08-10 | 2018-02-02 | 国家电网公司 | Realize full-time sequence grid measurement data synchronized sampling and the investigating method of processing |
CN107966235A (en) * | 2017-11-27 | 2018-04-27 | 中国科学院工程热物理研究所 | A kind of high-precision pressure measuring system of variable reference pressure |
WO2019047606A1 (en) * | 2017-09-05 | 2019-03-14 | 百度在线网络技术(北京)有限公司 | Method and fpga for realizing synchronization of video frame data and gps time position |
CN113533815A (en) * | 2021-06-29 | 2021-10-22 | 电子科技大学 | Multi-channel sampling synchronization method based on time stamps |
CN113535620A (en) * | 2021-06-29 | 2021-10-22 | 电子科技大学 | Multichannel synchronous high-speed data acquisition device |
CN114785445A (en) * | 2022-06-20 | 2022-07-22 | 中国空气动力研究与发展中心低速空气动力研究所 | Clock synchronization implementation method of wind tunnel pressure measurement system |
CN114964701A (en) * | 2022-05-07 | 2022-08-30 | 中国航空工业集团公司哈尔滨空气动力研究所 | Model surface pressure continuous measurement test system based on pressure scanning valve |
CN115014696A (en) * | 2022-08-08 | 2022-09-06 | 中国空气动力研究与发展中心高速空气动力研究所 | Method for synchronous acquisition and integrated processing of wind tunnel multi-signal source data |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6948910B2 (en) * | 2002-07-12 | 2005-09-27 | Polacsek Ronald R | Spiral-based axial flow devices |
CN101498757B (en) * | 2009-03-06 | 2012-05-30 | 深圳市双合电气股份有限公司 | Real-time dynamic monitoring and recording equipment of power system |
CN112073460B (en) * | 2020-08-03 | 2022-08-09 | 国网山东省电力公司电力科学研究院 | Data acquisition unit, data transmission method and comprehensive energy monitoring system |
-
2022
- 2022-11-17 CN CN202211437857.0A patent/CN115493800B/en active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5756892A (en) * | 1997-02-27 | 1998-05-26 | The United States Of America As Represented By The United States National Aeronautics And Space Administration | Apparatus for measuring ambient pressure within a gaseous flow field |
JP2005354568A (en) * | 2004-06-14 | 2005-12-22 | Sony Corp | Physical information acquisition method, physical information acquisition device and semiconductor device for detecting physical value distribution |
CN101534021A (en) * | 2009-04-27 | 2009-09-16 | 北京四方继保自动化股份有限公司 | Multimode data acquisitions and processing method applied to power automation system |
CN104113342A (en) * | 2013-11-28 | 2014-10-22 | 西安电子科技大学 | High-speed data synchronous circuit used for high-speed digital-to-analog converter |
CN107402597A (en) * | 2017-07-31 | 2017-11-28 | 上海联影医疗科技有限公司 | A kind of data and method, apparatus of clock alignment, medium and magnetic resonance equipment |
CN107656151A (en) * | 2017-08-10 | 2018-02-02 | 国家电网公司 | Realize full-time sequence grid measurement data synchronized sampling and the investigating method of processing |
WO2019047606A1 (en) * | 2017-09-05 | 2019-03-14 | 百度在线网络技术(北京)有限公司 | Method and fpga for realizing synchronization of video frame data and gps time position |
CN107966235A (en) * | 2017-11-27 | 2018-04-27 | 中国科学院工程热物理研究所 | A kind of high-precision pressure measuring system of variable reference pressure |
CN113533815A (en) * | 2021-06-29 | 2021-10-22 | 电子科技大学 | Multi-channel sampling synchronization method based on time stamps |
CN113535620A (en) * | 2021-06-29 | 2021-10-22 | 电子科技大学 | Multichannel synchronous high-speed data acquisition device |
CN114964701A (en) * | 2022-05-07 | 2022-08-30 | 中国航空工业集团公司哈尔滨空气动力研究所 | Model surface pressure continuous measurement test system based on pressure scanning valve |
CN114785445A (en) * | 2022-06-20 | 2022-07-22 | 中国空气动力研究与发展中心低速空气动力研究所 | Clock synchronization implementation method of wind tunnel pressure measurement system |
CN115014696A (en) * | 2022-08-08 | 2022-09-06 | 中国空气动力研究与发展中心高速空气动力研究所 | Method for synchronous acquisition and integrated processing of wind tunnel multi-signal source data |
Non-Patent Citations (2)
Title |
---|
大型旋转机械振动信号数据采集器的研究;徐峰等;《清华大学学报(自然科学版)》;19980610(第06期);全文 * |
脉动压力分布同步测量的探讨;郭明晻等;《实验力学》;20050330(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115493800A (en) | 2022-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115493800B (en) | Synchronous parallel acquisition system for steady-state pressure and pulsating pressure data and application method | |
CN104749498B (en) | A kind of detection of Portable PD On-Line and diagnostic device | |
CN115508040B (en) | Synchronous parallel acquisition system for data of speed field and temperature field and application method | |
CN105738836B (en) | A kind of DC/DC converters automatization test system | |
CN115638952A (en) | Multipoint pressure data synchronous parallel acquisition system and application method | |
CN101520640A (en) | Time interval measuring instrument based on FPGA | |
CN101629995A (en) | Detection method of phasor measurement unit PMU | |
CN202710986U (en) | Backplane bus type structure of vibration monitoring protection device | |
CN108845281A (en) | A kind of online Determination system of equipment for monitoring power quality and its detection method | |
CN115077852A (en) | Dynamic measurement system and measurement method for continuous transonic wind tunnel | |
CN104459435A (en) | Wiring verification method and device for transformer substation | |
CN111934760A (en) | Signal processing transmission delay characteristic detection device and method and terminal equipment | |
CN115061082A (en) | Signal processing method and device for interferometer direction finding narrow-band receiver | |
CN209543097U (en) | Data collection system | |
CN110955179A (en) | Dual-channel shared clock trigger delay adjusting device based on PCI bus | |
CN201331680Y (en) | Time interval measuring instrument based on FPGA | |
CN104345284B (en) | A kind of radiation effect TT&C system of DC DC power modules | |
CN201766598U (en) | Multi-rate one-with-four bit code error analysis detector | |
CN115493801B (en) | Steady-state pressure and dynamic data phase synchronization parallel acquisition and preprocessing system | |
CN114460829B (en) | Accelerated aging reliability test device and method for chip atomic clock | |
CN207424123U (en) | Unmanned aerial vehicle onboard Aerial Electronic Equipment tests system | |
CN106019395B (en) | A kind of electrical method receiving processing system and its method based on Correlation Identification | |
CN211015495U (en) | Bus architecture and case type acquisition instrument with same | |
CN211699411U (en) | Circuit analysis experiment box with detection device | |
CN202329799U (en) | Large-scale distributed synchronous data acquisition and analysis system |
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 |