A lightning effect test waveform splicing method and system

Abstract

The present application relates to a lightning effect test waveform splicing method and system, the lightning effect test waveform splicing method is: obtaining a waveform discrete data combination of lightning impulse current component A wave, B wave, C wave, and D wave through a lightning impulse current component generating device and an acquisition device, performing data processing on the acquired waveform discrete data combination to obtain spliced data with uniform and continuous data distribution, and displaying the spliced data using logarithmic coordinates; the data processing includes the following steps: signal amplification, de-offset, filtering, downsampling, peak detection, edge pulse width calculation, data splicing, polynomial interpolation, etc. The present invention integrates the test of different current component combinations into one system, greatly improving the test efficiency; adopts a weighted average algorithm to solve the problem of uneven thickness caused by sampling differences of different current component data; adopts polynomial interpolation to solve the problem of discontinuous current component splicing; and adopts logarithmic coordinates for data display to solve the problem of uneven distribution of current components.

Description

Technical Field

[0001] The present application belongs to the technical field of aviation lightning protection experiments, and in particular relates to a lightning effect test waveform splicing method and splicing system. Background Art

[0002] The aircraft lightning direct effects test divides the simulated lightning impulse current into multiple different components, and selects different current components and their combinations for testing based on the lightning attachment area.

[0003] The current components include four types: A, B, C, and D. Current component A is the initial peak current, with a peak value of 200kA and a duration of less than or equal to 500μs; current component B is the intermediate current, with an average amplitude of 2kA and a duration of less than or equal to 5ms; current component C is the continuous current, with a duration of 0.25 to 1s; and current component D is the repetitive discharge current, with a peak value of 100kA and a duration of less than or equal to 500μs.

[0004] The lightning attachment area consists of five areas: 1A, 1B, 2A, 2B, and 3. Area 1A is a surface with a high probability of initial lightning attachment and a high probability of the attachment point position changing over time. It is tested using a combination of the A and B components. Area 1B is a surface with a high probability of initial lightning attachment and a low probability of the attachment point position changing over time. It is tested using a combination of the A, B, C, and D components. Area 2A is a surface with a high probability of lightning being blown through by airflow from the initial attachment point and a high probability of the attachment point position changing over time. It is tested using a combination of the B, C, and D components. Area 2B is a surface with a high probability of lightning being blown through by airflow from the initial attachment point and a low probability of the attachment point position changing over time. It is tested using a combination of the B, C, and D components. Area 3 covers the remaining surfaces and is tested using a combination of the A and C components.

[0005] The current existing technology uses four generators A, B, C, and D to output the four current components of the test respectively, and then uses the acquisition device to obtain the waveform to verify the reliability of the test. Generally, there are two acquisition methods:

[0006] Each generator uses a separate acquisition device. This method can capture the waveform characteristics of a single current component, but cannot obtain the characteristics of the component combination;

[0007] All generating devices use one acquisition device. This method can capture the characteristics of the current component combination, but the captured waveform is discontinuous and the waveform characteristics of each current component cannot be obtained. Summary of the Invention

[0008] In order to solve the above-mentioned deficiencies in the prior art, the present invention provides a lightning effect test waveform splicing method and system.

[0009] The technical solution adopted by the present invention to solve the above technical problems is:

[0010] A lightning effect test waveform splicing method is disclosed. The method obtains a waveform discrete data combination of lightning impulse current component A wave, B wave, C wave, and D wave through a lightning impulse current component generating device and a collecting device. The collected waveform discrete data combination is processed to obtain spliced data with uniform and continuous data distribution, and the spliced data is displayed using logarithmic coordinates.

[0011] Preferably, in the lightning effect test waveform splicing method of the present invention, the data processing comprises the following steps:

[0012] Signal amplification: converting the waveform discrete data combination into an actual waveform combination through A1=A0 / μ, where A1 represents the actual waveform data combination, A0 represents the collected waveform discrete data combination, and μ represents the current attenuation coefficient;

[0013] To remove the offset, use arithmetic averaging on the 0-10% data segment of A1 to obtain the average offset X, and subtract the average offset X from A1 to obtain the waveform data combination A2 after removing the offset;

[0014] Filtering: using a sliding average filtering algorithm on A2 to obtain a filtered waveform data combination A3;

[0015] Downsampling: Use the weighted average algorithm to reduce the number of samples of A3 to obtain the data combination A4;

[0016] Peak detection: Use the peak detection algorithm to find the peak point of A4. Using the peak point as the starting point, traverse A4 in reverse direction, taking the 90% peak point as tr90, the 50% peak point as tr50, and the 10% peak point as tr10. Using the peak point as the starting point, traverse A4 in forward direction, taking the 50% peak point as tt50.

[0017] Edge pulse width calculation: calculate the waveform rise time and pulse width based on the peak point in the peak detection step: rise time Tr = (tr90-tr10) / 0.8, pulse width Td = (tr90-tr10)*0.5 / 0.8+(tt50-tr50);

[0018] Data splicing, judging the current component type according to the current component peak value, rise time Tr and pulse width Td, and sorting them in the order of A wave, B wave, C wave and D wave;

[0019] Polynomial interpolation is performed by using peak detection to obtain the 0% peak point of the rising time period and the 0% peak point of the falling time period as the starting and ending points of the current components A and D waves. If the 0% peak point of the falling time period is not obtained, the end point of A4 is used as the 0% peak point of the falling time period, and the 95% peak point of the rising time period and the 95% peak point of the falling time period are obtained as the starting and ending points of the current components B and C waves. Polynomial interpolation is performed between the starting and ending points of the current components.

[0020] Preferably, in the lightning effect test waveform splicing method of the present invention, the polynomial interpolation is an order linear interpolation.

[0021] The present invention also provides a lightning effect test waveform splicing system, including a device driving module, a data processing module and a data display module. The device driving module is used to control the generating devices and collecting devices of four lightning impulse current components A wave, B wave, C wave and D wave to obtain a waveform discrete data set. The data processing module processes the collected waveform discrete data set to obtain spliced data with uniform and continuous data distribution. The data display module uses logarithmic coordinates to display the processed waveform.

[0022] Preferably, in the lightning effect test waveform splicing system of the present invention, the data processing module includes the following submodules:

[0023] The signal amplification submodule converts the waveform discrete data combination into an actual waveform combination through A1=A0 / μ, where A1 represents the actual waveform data combination, A0 represents the collected waveform discrete data combination, and μ represents the current attenuation coefficient;

[0024] The de-offset quantum module uses arithmetic averaging on the 0-10% data segment of A1 to obtain the average offset X, and subtracts the average offset X from A1 to obtain the waveform data combination A2 after de-offset;

[0025] The filtering submodule uses a sliding average filtering algorithm on A2 to obtain the filtered waveform data combination A3;

[0026] The downsampling submodule uses a weighted average algorithm to reduce the number of samples of A3 to obtain the data combination A4;

[0027] The peak detection submodule uses the peak detection algorithm to find the peak point of A4. It traverses A4 in reverse with the peak point as the starting point, taking the 90% peak point as tr90, the 50% peak point as tr50, and the 10% peak point as tr10. It traverses A4 in forward direction with the peak point as the starting point, taking the 50% peak point as tt50.

[0028] The edge pulse width calculation submodule calculates the waveform rise time and pulse width according to the peak point of the peak detection submodule. The rise time Tr = (tr90-tr10) / 0.8, and the pulse width Td = (tr90-tr10)*0.5 / 0.8+(tt50-tr50);

[0029] The data splicing submodule determines the current component type based on the current component peak value, rise time Tr and pulse width Td, and sorts them in the order of A wave, B wave, C wave and D wave;

[0030] The polynomial interpolation submodule uses peak detection to obtain the 0% peak point of the rising time period and the 0% peak point of the falling time period as the starting and ending points of the current components A wave and D wave. If the 0% peak point of the falling time period is not obtained, the end point of A4 is used as the 0% peak point of the falling time period, and the 95% peak point of the rising time period and the 95% peak point of the falling time period are obtained as the starting and ending points of the current components B wave and C wave, and polynomial interpolation is performed between the starting and ending points of the current components.

[0031] The beneficial effects of the present invention are:

[0032] The present invention integrates the testing of different current component combinations into one system, greatly improving the test efficiency. A weighted average algorithm is used to solve the problem of uneven coarseness and fineness caused by sampling differences in data of different current components. Polynomial interpolation is used to solve the problem of discontinuous splicing of current components. Logarithmic coordinates are used for data display to solve the problem of uneven distribution of current components. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] The technical solution of the present application is further described below with reference to the accompanying drawings and embodiments.

[0034] Figure 1 It is the current component waveform diagram in the background technology of this application;

[0035] Figure 2 This is a flow chart of the data processing steps of the lightning effect test waveform splicing method according to an embodiment of the present application;

[0036] Figure 3 This is a structural block diagram of a lightning effect test waveform splicing system according to an embodiment of the present application;

[0037] Figure 4 This is a structural block diagram of the data processing module of an embodiment of the present application. DETAILED DESCRIPTION

[0038] It should be noted that, unless there is any conflict, the embodiments and features in the embodiments of this application can be combined with each other.

[0039] The technical solution of the present application will be described in detail below with reference to the accompanying drawings and in combination with embodiments.

[0040] Example 1

[0041] This embodiment provides a lightning effect test waveform splicing method, which obtains a waveform discrete data combination of the lightning impulse current component A wave, B wave, C wave, and D wave through a lightning impulse current component generating device and a collection device, processes the collected waveform discrete data combination to obtain spliced data with uniform and continuous data distribution, and displays the spliced data using logarithmic coordinates.

[0042] The data processing steps in this embodiment are as follows: Figure 2 , as follows:

[0043] S1, signal amplification, converting the waveform discrete data combination into the actual waveform combination through A1=A0 / μ, where A1 represents the actual waveform data combination, A0 represents the collected waveform discrete data combination, and μ represents the current attenuation coefficient;

[0044] S2, offset removal, uses arithmetic averaging on the 0-10% data segment of A1 to obtain the average offset X, and subtracts the average offset X from A1 to obtain the waveform data combination A2 after offset removal;

[0045] S3, filtering, using a sliding average filtering algorithm on A2 to obtain the filtered waveform data combination A3;

[0046] S4, downsampling, uses the weighted average algorithm to reduce the number of samples of A3 to obtain the data combination A4;

[0047] S5, peak detection, use the peak detection algorithm to find the peak point of A4, traverse A4 in reverse with the peak point as the starting point, take the 90% peak point as tr90, the 50% peak point as tr50, and the 10% peak point as tr10, and traverse A4 in forward direction with the peak point as the starting point, take the 50% peak point as tt50;

[0048] S6, edge pulse width calculation, calculate the waveform rise time and pulse width according to the peak point in the peak detection step, rise time Tr = (tr90-tr10) / 0.8, pulse width Td = (tr90-tr10)*0.5 / 0.8+(tt50-tr50);

[0049] S7, data splicing, judging the current component type according to the current component peak value, rise time Tr and pulse width Td, and sorting them in the order of A wave, B wave, C wave and D wave;

[0050] S8, polynomial interpolation, using peak detection to obtain the 0% peak point of the rising time period and the 0% peak point of the falling time period as the starting and ending points of the current components A wave and D wave. If the 0% peak point of the falling time period is not obtained, the end point of A4 is used as the 0% peak point of the falling time period, and the 95% peak point of the rising time period and the 95% peak point of the falling time period are obtained as the starting and ending points of the current components B wave and C wave, and polynomial interpolation is performed between the starting and ending points of the current components.

[0051] The polynomial interpolation method described in this embodiment is a linear interpolation method. In addition, parabolic interpolation and cubic polynomial interpolation can also be used. The use of polynomial interpolation can make the current components spliced continuously and completely.

[0052] Example 2

[0053] This embodiment provides a lightning effect test waveform splicing system, such as Figure 3 shown.

[0054] The lightning effect test waveform splicing system includes a device driver module, a data processing module, and a data display module. The device driver module is used to control the generation and collection devices of four lightning impulse current components A wave, B wave, C wave, and D wave to obtain a discrete waveform data set of A wave, B wave, C wave, and D wave. The data processing module analyzes and processes the collected discrete waveform data set to obtain uniformly distributed and continuous spliced data. The data display module uses logarithmic coordinates to display the processed waveform.

[0055] The data processing module of this embodiment is as follows Figure 4 As shown, it specifically includes the following sub-modules:

[0056] The signal amplification submodule converts the waveform discrete data combination into an actual waveform combination through A1=A0 / μ, where A1 represents the actual waveform data combination, A0 represents the collected waveform discrete data combination, and μ represents the current attenuation coefficient;

[0057] The de-offset quantum module uses arithmetic averaging on the 0-10% data segment of A1 to obtain the average offset X, and subtracts the average offset X from A1 to obtain the waveform data combination A2 after de-offset;

[0058] The filtering submodule uses a sliding average filtering algorithm on A2 to obtain the filtered waveform data combination A3;

[0059] The downsampling submodule uses a weighted average algorithm to reduce the number of samples of A3 to obtain the data combination A4;

[0060] The peak detection submodule uses the peak detection algorithm to find the peak point of A4. It traverses A4 in reverse with the peak point as the starting point, taking the 90% peak point as tr90, the 50% peak point as tr50, and the 10% peak point as tr10. It traverses A4 in forward direction with the peak point as the starting point, taking the 50% peak point as tt50.

[0061] The edge pulse width calculation submodule calculates the waveform rise time and pulse width according to the peak point of the peak detection submodule. The rise time Tr = (tr90-tr10) / 0.8, and the pulse width Td = (tr90-tr10)*0.5 / 0.8+(tt50-tr50);

[0062] The data splicing submodule determines the current component type based on the current component peak value, rise time Tr and pulse width Td, and sorts them in the order of A wave, B wave, C wave and D wave;

[0063] The polynomial interpolation submodule uses peak detection to obtain the 0% peak point of the rising time period and the 0% peak point of the falling time period as the starting and ending points of the current components A wave and D wave. If the 0% peak point of the falling time period is not obtained, the end point of A4 is used as the 0% peak point of the falling time period, and the 95% peak point of the rising time period and the 95% peak point of the falling time period are obtained as the starting and ending points of the current components B wave and C wave, and polynomial interpolation is performed between the starting and ending points of the current components.

[0064] Since the pulse widths of the current components A, B, C, and D waves vary greatly, it is difficult to display the short pulse width current components (A and D current components) less than or equal to 500 μs using a general coordinate axis. Therefore, this embodiment uses a logarithmic coordinate axis to display the spliced waveform data to make the above current flow distribution more uniform.

[0065] Based on the above-mentioned ideal embodiments of this application, and in accordance with the above description, relevant personnel can make various changes and modifications without departing from the scope of the technical concept of this application. The technical scope of this application is not limited to the contents of the specification, but must be determined according to the scope of the claims.

Claims

1. A lightning effect test waveform splicing method, characterized in that: A waveform discrete data combination of the lightning impulse current component A wave, B wave, C wave, and D wave is obtained through a lightning impulse current component generating device and a collecting device. A weighted average algorithm is used on the waveform discrete data combination to obtain a current component with uniform thickness. After splicing the current components, polynomial interpolation is performed between the starting point and the end point of the current component. After data processing is performed on the collected waveform discrete data combination, spliced data with uniform and continuous data distribution is obtained, and the spliced data is displayed using logarithmic coordinates.

2. The lightning effect test waveform splicing method according to claim 1, characterized in that: The data processing includes the following steps: Signal amplification: converting the waveform discrete data combination into an actual waveform combination through A1=A0 / μ, where A1 represents the actual waveform data combination, A0 represents the collected waveform discrete data combination, and μ represents the current attenuation coefficient; To remove the offset, use arithmetic averaging on the 0-10% data segment of A1 to obtain the average offset X. Subtract the average offset X from A1 to obtain the waveform data combination A2 after removing the offset. Filtering: using a sliding average filtering algorithm on A2 to obtain a filtered waveform data combination A3; Downsampling: Use the weighted average algorithm to reduce the number of samples of A3 to obtain the data combination A4; Peak detection: Use the peak detection algorithm to find the peak point of A4. Using the peak point as the starting point, traverse A4 in reverse direction, taking the 90% peak point as tr90, the 50% peak point as tr50, and the 10% peak point as tr10. Using the peak point as the starting point, traverse A4 in forward direction, taking the 50% peak point as tt50. Edge pulse width calculation: Calculate the waveform rise time and pulse width based on the peak point in the peak detection step. Rise time Tr = (tr90-tr10) / 0.8, pulse width Td = (tr90-tr10)*0.5 / 0.8+(tt50-tr50); Data splicing, judging the current component type according to the current component peak value, rise time Tr and pulse width Td, and sorting them in the order of A wave, B wave, C wave and D wave; Polynomial interpolation is performed using peak detection to obtain the 0% peak point of the rising time period and the 0% peak point of the falling time period as the starting and ending points of the current components A and D waves. If the 0% peak point of the falling time period is not obtained, the end point of A4 is used as the 0% peak point of the falling time period. The 95% peak point of the rising time period and the 95% peak point of the falling time period are obtained as the starting and ending points of the current components B and C waves. Polynomial interpolation is performed between the starting and ending points of the current components.

3. The lightning effect test waveform splicing method according to claim 2, characterized in that: The polynomial interpolation is an order linear interpolation.

4. A lightning effect test waveform splicing system, characterized in that: The invention comprises a device driving module, a data processing module and a data display module. The device driving module is used to control the generating devices and collecting devices of four lightning impulse current components A wave, B wave, C wave and D wave to obtain a waveform discrete data set. The data processing module combines the waveform discrete data using a weighted average algorithm to obtain current components with uniform thickness. After splicing the current components, polynomial interpolation is performed between the starting point and the end point of the current components. After data processing the collected waveform discrete data set, the data distribution is uniform and continuous spliced data is obtained. The data display module uses logarithmic coordinates to display the processed waveform.

5. The lightning effect test waveform splicing system according to claim 4, characterized in that: The data processing module includes the following submodules: The signal amplification submodule converts the waveform discrete data combination into an actual waveform combination through A1=A0 / μ, where A1 represents the actual waveform data combination, A0 represents the collected waveform discrete data combination, and μ represents the current attenuation coefficient; The de-offset quantum module uses arithmetic averaging on the 0-10% data segment of A1 to obtain the average offset X. A1 minus the average offset X is used to obtain the waveform data combination A2 after de-offset. The filtering submodule uses a sliding average filtering algorithm on A2 to obtain the filtered waveform data combination A3; The downsampling submodule uses a weighted average algorithm to reduce the number of samples of A3 to obtain the data combination A4; The peak detection submodule uses the peak detection algorithm to find the peak point of A4. It traverses A4 in reverse direction starting from the peak point, taking the 90% peak point as tr90, the 50% peak point as tr50, and the 10% peak point as tr10. It then traverses A4 in forward direction starting from the peak point, taking the 50% peak point as tt50. The edge pulse width calculation submodule calculates the waveform rise time and pulse width based on the peak point of the peak detection submodule. The rise time Tr=(tr90-tr10) / 0.8, and the pulse width Td=(tr90-tr10)*0.5 / 0.8+(tt50-tr50); The data splicing submodule determines the current component type based on the current component peak value, rise time Tr and pulse width Td, and sorts them in the order of A wave, B wave, C wave and D wave; The polynomial interpolation submodule uses peak detection to obtain the 0% peak point of the rising time period and the 0% peak point of the falling time period as the starting and ending points of the current components A and D waves. If the 0% peak point of the falling time period is not obtained, the end point of A4 is used as the 0% peak point of the falling time period, and the 95% peak point of the rising time period and the 95% peak point of the falling time period are obtained as the starting and ending points of the current components B and C waves. Polynomial interpolation is performed between the starting and ending points of the current components.

Images