Material dynamic test data processing method and system, electronic equipment and storage medium

Abstract

The invention discloses a material dynamic test data processing method and system, electronic equipment and a storage medium, and the method comprises the steps: after obtaining an incident rod signal curve and a transmission rod signal curve, firstly extracting an incident wave curve and a reflected wave curve from the incident rod signal curve, extracting a transmission wave curve from the transmission rod signal curve, and then extracting a reflected wave curve from the transmission rod signal curve; and then automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width, filtering and smoothing the incident wave curve, the reflected wave curve and the transmitted wave curve, and automatically calculating to obtain an engineering stress-strain curve. Due to the fact that the pulse widths of the incident wave, the reflected wave and the transmitted wave are automatically aligned, the consistency of the data time sequence is guaranteed, the accuracy of an engineering stress-strain curve is guaranteed, and the real dynamic mechanical property can be accurately reflected.

Description

Technical Field

[0001] The present invention relates to the technical field of material dynamic testing, and in particular, to a material dynamic testing data processing method and system, electronic equipment, and a computer-readable storage medium. Background Art

[0002] The mechanical properties of materials mainly refer to the properties of materials under external forces, including strength, hardness, plasticity, impact toughness, etc. These performance indicators are of great significance for engineering structure design, material selection and mechanical structure failure analysis. The real stress-strain curve of materials is mainly the relationship curve between stress and strain of materials during the stress process, which reflects the mechanical properties of materials. At present, there are three types of mechanical property tests of materials, namely low strain rate quasi-static test, medium strain rate dynamic test and high strain rate dynamic test. Among them, low strain rate quasi-static test (strain rate <0.01 / s) is carried out by electronic universal testing machine, dynamic mechanical properties of materials under medium strain rate (strain rate 0.1~500 / s) are carried out by hydraulic servo testing machine, and dynamic mechanical properties of materials under high strain rate (strain rate >1000 / s) are carried out by Hopkinson bar. Due to different test methods, the engineering stress-strain curve data of low strain rate quasi-static test and medium strain rate dynamic test can be directly measured, while the engineering stress-strain curve data of high strain rate dynamic test is obtained through existing automatic data processing program. However, the existing automated data processing program cannot automatically align the pulse widths of the incident wave, reflected wave, and transmitted wave, which results in poor accuracy of the stress-strain curve and makes it difficult to accurately reflect the true dynamic mechanical properties. In addition, the existing automated data processing program cannot automatically determine the validity of the test data and cannot guarantee the accuracy and reliability of the test results. In addition, the existing automated data processing program cannot perform plastic flow constitutive fitting, lacks fitting analysis of test data, and cannot allow the test data to be simulated and applied. Summary of the invention

[0003] The present invention provides a material dynamic test data processing method and system, electronic equipment, and computer-readable storage medium, which can automatically align the pulse widths of incident waves, reflected waves, and transmitted waves, ensuring the consistency of data timing, thereby ensuring the accuracy of engineering stress-strain curves, and can accurately reflect the true dynamic mechanical properties.

[0004] According to one aspect of the present invention, a material dynamic test data processing method is provided, comprising the following contents: Obtaining an incident rod signal curve and a transmission rod signal curve, extracting an incident wave curve and a reflected wave curve from the incident rod signal curve, extracting a transmitted wave curve from the transmission rod signal curve, and automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width; Filter and smooth the incident wave curve, reflected wave curve and transmitted wave curve; The engineering stress-strain curve is calculated based on the incident wave curve, reflected wave curve and transmitted wave curve.

[0005] Furthermore, the process of extracting the incident wave curve and the reflected wave curve from the incident rod signal curve, extracting the transmitted wave curve from the transmitted rod signal curve, and automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width includes the following: Find the extreme points of the incident wave and the reflected wave in the incident rod signal curve, and use them as the characteristic points of the incident wave curve and the reflected wave curve, respectively; find the extreme point of the transmitted wave in the transmitted rod signal curve, and use it as the characteristic point of the transmitted wave curve; Based on the characteristic point of the incident wave curve, traverse forward in the incident rod signal curve, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, and use it as the starting point of the incident wave waveform; then traverse backward in the incident rod signal curve based on the characteristic point, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all less than the measured value of this point, and use it as the ending point of the incident wave waveform, and determine the complete incident wave waveform curve based on the starting point and the ending point; Based on the characteristic point of the reflected wave curve, traverse forward in the incident rod signal curve to find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, and use it as the starting point of the reflected wave waveform; then traverse backward in the incident rod signal curve based on the characteristic point to find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all less than the measured value of this point, and use it as the ending point of the reflected wave waveform, and determine the complete reflected wave waveform curve based on the starting point and the ending point; Based on the characteristic points of the transmission wave curve, traverse forward in the transmission rod signal curve, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, use it as the starting point of the transmission wave waveform, and determine the end point based on the same pulse width as the incident wave waveform, so as to determine the complete transmission wave waveform curve; The starting points of the incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve are aligned, and the data with redundant time length in the reflected wave waveform curve are eliminated based on the principle of the same pulse width.

[0006] Furthermore, after obtaining the incident rod signal curve, the following contents are also included: The test type is determined according to the incident wave signal value. If the incident wave signal value is positive, it is determined to be a tensile test. If the incident wave signal value is negative, it is determined to be a compression test.

[0007] Furthermore, after obtaining the incident rod signal curve, the following contents are also included: Determine the validity of the incident rod signal and the transmitted rod signal.

[0008] Furthermore, the process of determining the validity of the incident rod signal includes the following: Obtain the theoretical incident rod signal curve, calculate the ratio of the difference between the maximum value of the measured incident rod signal and the maximum value of the theoretical incident rod signal to the maximum value of the theoretical incident rod signal, if the ratio is greater than a preset threshold, the incident rod signal is determined to be invalid, otherwise the incident rod signal is determined to be valid; calculate the correlation between the measured incident rod signal curve and the theoretical incident rod signal curve, if the correlation between the two is greater than a preset threshold, the incident rod signal is determined to be valid, otherwise the incident rod signal is determined to be invalid.

[0009] Furthermore, the following contents are also included: The calibration signal curve of the dynamic test calibration sensor is obtained, and the incident rod signal curve and the transmission rod signal curve are compared with the calibration signal curve respectively. The validity of the incident rod strain gauge signal and the transmission rod strain gauge signal are judged according to the comparison results.

[0010] Furthermore, the following contents are also included: The true stress-plastic strain curve is converted based on the engineering stress-strain curve, and then the constitutive model is used to perform fitting and extrapolation analysis based on the true stress-plastic strain curve to obtain the plastic flow curve.

[0011] In addition, the present invention also provides a material dynamic test data processing system, comprising: A waveform curve automatic alignment module is used to obtain the incident rod signal curve and the transmission rod signal curve, extract the incident wave curve and the reflected wave curve from the incident rod signal curve, extract the transmitted wave curve from the transmission rod signal curve, and automatically align the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width; A data filtering and smoothing processing module is used to filter and smooth the incident wave curve, the reflected wave curve and the transmitted wave curve; The engineering stress-strain curve generation module is used to calculate the engineering stress-strain curve based on the incident wave curve, the reflected wave curve and the transmitted wave curve.

[0012] In addition, the present invention also provides an electronic device, including a processor and a memory, wherein the memory stores a computer program, and the processor executes the steps of the above method by calling the computer program stored in the memory.

[0013] In addition, the present invention also provides a computer-readable storage medium for storing a computer program for processing material dynamic test data, wherein the computer program executes the steps of the method described above when running on a computer.

[0014] The present invention has the following beneficial effects: The material dynamic test data processing method of the present invention, after obtaining the incident rod signal curve and the transmission rod signal curve, first extracts the incident wave curve and the reflected wave curve from the incident rod signal curve, extracts the transmitted wave curve from the transmission rod signal curve, and then automatically aligns the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width, and then automatically calculates the engineering stress-strain curve after filtering and smoothing the incident wave curve, the reflected wave curve and the transmitted wave curve. Since the pulse widths of the incident wave, the reflected wave and the transmitted wave are automatically aligned, the consistency of the data time series is guaranteed, thereby ensuring the accuracy of the engineering stress-strain curve, and can accurately reflect the real dynamic mechanical properties.

[0015] In addition, the material dynamic test data processing system of the present invention also has the above advantages.

[0016] In addition to the above-described purposes, features and advantages, the present invention has other purposes, features and advantages. The present invention will be further described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The drawings constituting a part of this application are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings: Figure 1 It is a flow chart of a material dynamic test data processing method according to a preferred embodiment of the present application; Figure 2 It is a schematic diagram of a preferred embodiment of the present application in which a test piece is mounted on a Hopkinson bar to perform a high strain rate dynamic test; Figure 3 yes Figure 1 Schematic diagram of the sub-process of step S1; Figure 4 It is a schematic diagram of finding the extreme points of the incident wave curve, the reflected wave curve and the transmitted wave curve in the preferred embodiment of the present application; Figure 5 It is a schematic diagram of finding the starting point and the ending point of the incident wave curve, the reflected wave curve and the transmitted wave curve in the preferred embodiment of the present application; Figure 6 is a schematic diagram of automatically aligning an incident wave curve, a reflected wave curve, and a transmitted wave curve in a preferred embodiment of the present application; Figure 7is a schematic diagram of an engineering stress-strain curve generated in a preferred embodiment of the present application; Figure 8 is a schematic diagram of data signals for a tensile test in a preferred embodiment of the present application; Fig. 9 is a schematic diagram of data signals for a compression test in a preferred embodiment of the present application; Fig.10 is a schematic diagram of an incident rod signal curve and a theoretical incident rod signal curve measured when a tensile test is performed in a preferred embodiment of the present application; Fig.11 is a schematic diagram showing that the transmission rod signal curve in the preferred embodiment of the present application does not return to zero; Fig.12 It is a schematic diagram of the installation of the dynamic test calibration sensor in the preferred embodiment of the present application; Fig.13 It is a schematic diagram of force value curves of the incident rod strain gauge signal, the piezoelectric crystal force value signal, the calibration strain gauge signal and the transmission rod strain gauge signal obtained when the strain gauge is calibrated in the preferred embodiment of the present application; Fig.14 It is a schematic diagram of the alignment of the force curves of the incident rod strain gauge signal, the piezoelectric crystal force value signal, the calibration strain gauge signal and the transmission rod strain gauge signal obtained during calibration in the preferred embodiment of the present application; Fig.15 It is another flow chart of the material dynamic test data processing method of the preferred embodiment of the present application; Fig.16 It is a schematic diagram of plastic flow curves obtained by fitting simulation using different constitutive models in the preferred embodiment of the present application; Fig.17 It is a schematic diagram of the module structure of a material dynamic test data processing system according to another embodiment of the present application. DETAILED DESCRIPTION

[0018] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0019] Reference Figure 1 The preferred embodiment of the present application provides a material dynamic test data processing method, comprising the following contents: Step S1: obtaining an incident rod signal curve and a transmission rod signal curve, extracting an incident wave curve and a reflected wave curve from the incident rod signal curve, extracting a transmitted wave curve from the transmission rod signal curve, and automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width; Step S2: filtering and smoothing the incident wave curve, the reflected wave curve and the transmitted wave curve; Step S3: Calculate the engineering stress-strain curve based on the incident wave curve, the reflected wave curve and the transmitted wave curve.

[0020] It can be understood that the material dynamic test data processing method of this embodiment, after obtaining the incident rod signal curve and the transmission rod signal curve, first extracts the incident wave curve and the reflected wave curve from the incident rod signal curve, extracts the transmitted wave curve from the transmission rod signal curve, and then automatically aligns the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width, and then filters and smoothes the incident wave curve, the reflected wave curve and the transmitted wave curve, and then automatically calculates the engineering stress-strain curve. Since the pulse widths of the incident wave, the reflected wave and the transmitted wave are automatically aligned, the consistency of the data time series is guaranteed, thereby ensuring the accuracy of the engineering stress-strain curve, and can accurately reflect the real dynamic mechanical properties.

[0021] It is understandable that Figure 2 As shown, after the specimen is mounted on the Hopkinson bar, a high strain rate dynamic test can be performed. The test data that can be obtained through the ultra-dynamic strain gauge and the data acquisition system include time, incident rod signal and transmission rod signal, wherein the incident rod signal is the strain gauge signal on the incident rod, and the transmission rod signal is the strain gauge signal on the transmission rod. In the step S1, the collected test data is obtained from the data acquisition system to generate the incident rod signal curve and the transmission rod signal curve. When performing a high strain rate dynamic test, after the incident pulse signal is input into the incident rod, the strain gauge on the incident rod will collect the first incident wave signal. When the incident pulse signal propagates backward to the sample, the first reflected wave signal and the first transmitted wave signal will be generated. The first reflected wave signal will propagate forward along the incident rod, and the first transmitted wave signal will propagate backward along the transmission rod; then, the first reflected wave signal will be reflected to form the second incident wave signal, which will propagate backward again, and the wave will continue to attenuate during the reciprocating process. Therefore, it is necessary to accurately extract the first incident wave signal and the first reflected wave signal from the incident rod signal curve, that is, the incident wave curve and the reflected wave curve, and accurately extract the first transmitted wave signal from the transmission rod signal curve, that is, the transmitted wave curve. Then, the incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve are automatically aligned based on the pulse width, so as to accurately calculate the engineering stress-strain curve in the subsequent step.

[0022] Among them, Figure 3 As shown, the process of extracting the incident wave curve and the reflected wave curve from the incident rod signal curve, extracting the transmitted wave curve from the transmitted rod signal curve, and automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width includes the following: Step S11: finding the extreme points of the incident wave and the reflected wave in the incident rod signal curve respectively, and using them as the characteristic points of the incident wave curve and the reflected wave curve respectively; finding the extreme point of the transmitted wave in the transmitted rod signal curve, and using it as the characteristic point of the transmitted wave curve; Step S12: based on the characteristic point of the incident wave curve, traverse forward in the incident rod signal curve, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, and use it as the starting point of the incident wave waveform; then based on the characteristic point, traverse backward in the incident rod signal curve, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all less than the measured value of this point, and use it as the ending point of the incident wave waveform, and determine the complete incident wave waveform curve based on the starting point and the ending point; Step S13: based on the characteristic point of the reflected wave curve, traverse forward in the incident rod signal curve to find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, and use it as the starting point of the reflected wave waveform; then based on the characteristic point, traverse backward in the incident rod signal curve to find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all less than the measured value of this point, and use it as the ending point of the reflected wave waveform; and determine the complete reflected wave waveform curve based on the starting point and the ending point; Step S14: based on the characteristic points of the transmission wave curve, traverse forward in the transmission rod signal curve, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, use it as the starting point of the transmission wave waveform, and determine the end point based on the same pulse width as the incident wave waveform, so as to determine the complete transmission wave waveform curve; Step S15: aligning the starting points of the incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve, and removing the data of redundant duration in the reflected wave waveform curve based on the principle of the same pulse width.

[0023] Specifically, the present invention first finds the extreme points of the incident wave and the reflected wave in the incident rod signal curve. Since the wave is constantly attenuated, the extreme points can be used as the characteristic points of the first incident wave and the first reflected wave, that is, the characteristic points of the incident wave curve and the reflected wave curve, such as Figure 4 As shown, the minimum value in the incident rod signal curve is the characteristic point of the first incident wave waveform curve, and the maximum value is the characteristic point of the first reflected wave waveform curve.

[0024] Then, based on the characteristic points of the incident wave waveform curve, traverse forward in the incident rod signal curve to find the first data point whose measurement value is less than 0.05V and the measurement values ​​of subsequent data points are all greater than the measurement value of this point, and take the first data point that meets this condition as the starting point of the incident wave waveform. Based on the characteristic point, traverse backward in the incident rod signal curve to find the first data point whose measurement value is less than 0.05V and the measurement values ​​of subsequent data points are all less than the measurement value of this point, and take the first data point that meets this condition as the end point of the incident wave waveform, so that the complete incident wave waveform curve can be determined based on the starting point and the end point. Similarly, based on the characteristic point of the reflected wave curve, traverse forward in the incident rod signal curve, find the first data point whose measured value is less than 0.05V and the measured values ​​of subsequent data points are all greater than the measured value of this point, and use it as the starting point of the reflected wave waveform. Then, based on the characteristic point, traverse backward in the incident rod signal curve, find the first data point whose measured value is less than 0.05V and the measured values ​​of subsequent data points are all less than the measured value of this point, and use it as the end point of the reflected wave waveform, and determine the complete reflected wave waveform curve based on the starting point and the end point. It can be understood that due to the noise in the signal collected by the strain gauge, the measured value of the starting point of the wave will not be zero.

[0025] Next, based on the characteristic points of the transmission wave curve, traverse forward in the transmission rod signal curve to find the first data point whose measured value is less than 0.05V and the measured values ​​of subsequent data points are all greater than the measured value of this point, and use it as the starting point of the transmission wave waveform. Since the transmission wave passes the extreme point, it means that the specimen has broken or the waveform time has ended. In order to ensure the accuracy of the transmission wave data acquisition, the present invention determines the end point of the transmission wave waveform based on the same pulse width as the incident wave waveform, that is, the pulse width of the transmission wave waveform curve must be the same as the pulse width of the incident wave waveform curve, so as to determine the complete transmission wave waveform curve. Among them, the schematic diagram of the positions of the starting points and end points of the incident wave waveform curve, the reflected wave waveform curve and the transmission wave waveform curve is shown in Figure 5 shown.

[0026] Then, the starting points of the incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve are aligned, and the data of the redundant time length in the reflected wave waveform curve is removed based on the principle of the same pulse width. Since the end point of the transmitted wave waveform curve is determined based on the same principle as the pulse width of the incident wave waveform curve, it is no longer necessary to remove the data points of the redundant time length, and the incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve can be automatically aligned, such as Figure 6 Among them, the strain gauge measures the voltage value, which needs to be converted into the strain value. The specific conversion process belongs to the existing technology and will not be repeated here.

[0027] It can be understood that the present invention uses extreme points as characteristic points of the first incident wave, reflected wave and transmitted wave, and respectively searches for the starting point and ending point of the incident wave waveform curve, the reflected wave waveform curve, and the starting point of the transmitted wave waveform curve based on the characteristic points, and determines the ending point of the transmitted wave waveform curve based on the principle of the same pulse width. The incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve can be automatically aligned based on the pulse width, thereby ensuring the timing consistency of the three waveform curves, thereby ensuring the accuracy of the engineering stress-strain curve, and accurately reflecting the true dynamic mechanical properties.

[0028] It can be understood that in step S2, the signal collected by the HPKS (Hopkinson bar) data acquisition system is superimposed with a noise signal, and the noise signal contains periodic interference signals such as the 50 Hz power frequency and its octaves. Therefore, the three aligned waveform curves need to be filtered. Digital filtering is a processing method that selects a part of the signal of interest from the collected discrete signal through mathematical operations. Its main functions are to filter out noise or false components in the test signal, improve the signal-to-noise ratio, smooth analysis data, suppress interference signals, separate frequency components, etc. Optionally, the present invention uses an infinite length IIR impulse response digital filter (Infinite Impulse Response Digital Filter) for low-pass filtering. The filtering expression of the IIR filter can be defined as a differential equation: ,in, and represent the input and output time domain signal sequences respectively, and are all filter coefficients, N is the order of the IIR filter, or the number of poles of the filter system transfer function, and M is the number of zeros of the filter system transfer function. Its transfer function can be expressed as: , H() represents the transfer function. Commonly used prototype generation functions of analog low-pass filters include Butterworth filter prototype, Chebyshev type I and type II filter prototypes, elliptic filter prototype, Bessel filter prototype, etc. Since the Butterworth filter is characterized by the flattest amplitude characteristic in the passband, the low-pass filter of the present invention preferably adopts the Butterworth filter, and the MATLAB signal processing toolbox provides a function for directly designing the Butterworth filter, which greatly facilitates the design of general filters. The function butter is the design function of the Butterworth filter, and its input parameters mainly include the filter order n, the normalized inherent cutoff frequency W n , the output is the filter coefficient; the function buttord is the Butterworth filter order selection function, which can calculate the order and normalized inherent cutoff frequency required for the filter design function. The input parameters of this function are the passband and stopband cutoff frequencies W p and Ws , passband and stopband ripple coefficient R p and R s ; The cutoff frequency is the normalized frequency, and the frequency range is from 0 to 1, where 1 corresponds to the folding frequency, which is half the sampling frequency.

[0029] In addition, considering that the noise signal in the dynamic signal collected by the HPKS test machine also includes irregular random interference signals, the frequency band of this random interference signal is relatively wide, and sometimes the high-frequency component accounts for a relatively large proportion, which cannot be removed by filtering, resulting in many burrs on the curve drawn by the collected discrete data, which is very rough. In addition, sometimes the test system is subject to some unexpected interference, causing the sampling signals of individual measuring points to produce trend items that deviate greatly from the baseline and have irregular shapes. Therefore, the present invention also smoothes the waveform curve after filtering, wherein there are two ways of smoothing: a) removing abnormal areas and interpolating data from other areas, and the interpolation methods include piecewise linear interpolation, cubic interpolation, and cubic spline interpolation, etc.; b) using the sliding average method to perform multiple data smoothing processes on the measurement signal to remove irregular trend items, wherein the expression of the sliding average method is: , where x represents the sampled data, y is the smoothed data, and N is the sliding average order. The averaged data can be averaged again. The larger the number of sliding averages, the smoother the curve.

[0030] It can be understood that in step S3, after filtering and smoothing the data of the three waveform curves, the one-wave method or the three-wave method can be selected to automatically calculate the engineering stress-strain curve, for example Figure 7 As shown. Among them, the one-wave method and the three-wave method are both existing algorithms, and the specific engineering stress-strain curve calculation process is not repeated here. For example, the calculation formula of the one-wave method is: ; ;

[0031] in, represents the engineering strain, represents the engineering strain rate, represents the engineering stress, represents the reflected strain, represents the transmission strain, represents the incident strain, C0 represents the wave velocity, l s represents the length of the specimen, t represents the duration of the stress wave, A represents the cross-sectional area of ​​the transmission rod, and A s represents the cross-sectional area of ​​the specimen, and E represents the elastic modulus.

[0032] In addition, in step S3, according to the GB / T228.1 standard, the present invention also extracts the engineering stress data segment and engineering strain data segment of the engineering stress-strain curve in the 0.0005~0.0025 segment, obtains a series of numerical values, and performs linear fitting to obtain the slope of the linear segment, which is the dynamic elastic modulus of the specimen, and compares the dynamic elastic modulus with the static elastic modulus. If the difference between the two exceeds a preset value (e.g., 10%), the engineering stress-strain curve is determined to be invalid, and if the difference between the two does not exceed the preset value, the engineering stress-strain curve is determined to be valid. At the same time, the average engineering strain rate in the test process is also analyzed and calculated. If the deviation between the average engineering strain rate and the target engineering strain rate is not greater than a preset value (e.g., 30%), it is determined to be a true strain rate result. If the deviation between the two is greater than the preset value, the test is determined to be invalid. Therefore, the present invention can also judge the validity of the calculated engineering stress-strain curve, further ensuring the accuracy and reliability of the test results.

[0033] Optionally, in step S1, after obtaining the incident rod signal curve, the following contents are also included: The test type is determined according to the incident wave signal value. If the incident wave signal value is positive, it is determined to be a tensile test. If the incident wave signal value is negative, it is determined to be a compression test.

[0034] Specifically, after obtaining the incident rod signal curve, since the incident rod signal includes the incident wave signal and the reflected wave signal, and the reflected wave signal is after the incident wave signal, the test type can be judged according to the incident wave signal value. If the incident wave signal value is positive, it means that the incident rod has been stretched and deformed, such as Figure 8 As shown, it is determined to be a tensile test. If the incident wave signal is negative, it means that the incident rod is compressed and deformed. Fig. 9 Therefore, the present invention can quickly and intuitively determine the test type according to the incident wave signal value.

[0035] Optionally, the step S1 further includes the following contents: Determine the validity of the incident rod signal and the transmitted rod signal.

[0036] The process of judging the validity of the incident rod signal includes the following: Obtain the theoretical incident rod signal curve, calculate the ratio of the difference between the maximum value of the measured incident rod signal and the maximum value of the theoretical incident rod signal to the maximum value of the theoretical incident rod signal, if the ratio is greater than a preset threshold, the incident rod signal is determined to be invalid, otherwise the incident rod signal is determined to be valid; calculate the correlation between the measured incident rod signal curve and the theoretical incident rod signal curve, if the correlation between the two is greater than a preset threshold, the incident rod signal is determined to be valid, otherwise the incident rod signal is determined to be invalid.

[0037] Specifically, the theoretical incident rod signal curve is calculated based on the impact rod length and impact velocity of the Hopkinson rod, such as Fig.10 The blue line in the figure represents the incident rod signal curve obtained by experimental measurement. Then, the ratio between the difference between the maximum value of the measured incident rod signal and the maximum value of the theoretical incident rod signal and the maximum value of the theoretical incident rod signal is calculated. The calculation formula is: , It indicates the maximum value of the incident rod signal. represents the maximum theoretical incident rod signal, Represents a ratio, if If it exceeds 10%, the incident rod signal is considered invalid. If the incident rod signal does not exceed 10%, the incident rod signal is determined to be valid. In addition, the correlation between the measured incident rod signal curve and the theoretical incident rod signal curve is also calculated. Specifically, the correlation can be calculated by the existing cosine similarity algorithm or correlation coefficient algorithm. If the correlation is greater than 0.9, the incident rod signal is determined to be valid. If the correlation is not greater than 0.9, the incident rod signal is determined to be invalid. Therefore, the waveform of the incident rod signal curve needs to meet the above two conditions at the same time to determine that the incident rod signal is valid, otherwise the incident rod signal is determined to be invalid. If the incident rod signal is determined to be invalid, the test system needs to be re-tested, including the alignment of the impact rod, the incident rod and the transmission rod, the optimization of the shaping sheet, the strain gauge inspection, etc., and the static electricity of the measurement system needs to be eliminated.

[0038] In addition, the process of judging the validity of the transmission rod signal includes the following: In the initial stage of the pre-experiment, that is, when the test piece is not installed for testing, the transmission rod is used as the incident rod, and the impact rod is used to impact the transmission rod, and the impact waveform (i.e., the incident waveform) and the reflection waveform collected for the first time on the transmission rod are measured. If the difference between the first impact waveform and the first reflection waveform does not exceed 2%, for example, the correlation between the impact waveform curve and the reflection waveform curve exceeds 98%, then the transmission rod signal is determined to be valid, otherwise the transmission rod signal is determined to be invalid. In addition, it is determined whether the collected transmission rod signal curve eventually returns to zero, for example, Fig.11 As shown, if the transmission rod signal curve does not return to zero after the test, it is determined that there is a problem with the transmission rod signal and the strain gauge and other acquisition systems need to be rechecked.

[0039] It can be understood that the present invention can also judge the effectiveness of the incident rod signal and the transmitted rod signal, thereby ensuring the accuracy and reliability of the test data.

[0040] Optionally, after obtaining the incident rod signal curve and the transmitted rod signal curve, the step S1 further includes the following: The calibration signal curve of the dynamic test calibration sensor is obtained, and the incident rod signal curve and the transmission rod signal curve are compared with the calibration signal curve respectively. The validity of the incident rod strain gauge signal and the transmission rod strain gauge signal are judged according to the comparison results.

[0041] Specifically, before the test, in order to ensure the reliability of the strain gauge data collected, the incident rod strain gauge and the transmission rod strain gauge need to be calibrated. The present invention uses a dynamic test calibration sensor for calibration. The dynamic test calibration sensor consists of a piezoelectric crystal force sensor and a strain gauge connected in series, and its two ends are connected to the incident rod and the transmission rod by threads, such as Fig.12 As shown. Among them, the piezoelectric crystal force sensor is an element that uses the piezoelectric effect to realize the sensing function. It has a forward piezoelectric effect and a reverse piezoelectric effect. The forward piezoelectric effect is: when mechanical pressure or stress is applied to the piezoelectric crystal, the crystal will generate an electric charge. This is because the distribution of positive and negative charges inside the crystal changes under the action of pressure, thereby forming a charge on the surface of the crystal and generating a voltage signal; the reverse piezoelectric effect is: when an electric field is applied to the piezoelectric crystal, the crystal will deform. This deformation is due to the electric field causing the distribution of positive and negative charges inside the crystal to change, resulting in a slight change in the size of the crystal. Typical piezoelectric crystal force sensors are made of piezoelectric materials, such as quartz, strontium tungstate and other crystals. These materials can generate measurable charge or voltage signals when external force or electric field is applied. Piezoelectric crystal force sensors can generate large electrical signal responses to small pressures or strains, and have a fast response speed, which can meet the dynamic performance test of high strain rates. When the bullet hits the incident rod, four signals can be collected, namely: the incident rod strain gauge signal, the piezoelectric crystal force value signal, the calibration strain gauge signal, and the transmission rod strain gauge signal. The four signal curves obtained after converting the strain gauge signal into the force value signal are as follows: Fig.13 Then align the time starting points of the four signal curves, as shown in Fig.14 As shown, a comparative analysis is performed. If the following conditions are met, the incident rod strain gauge signal and the transmitted rod strain gauge signal are determined to be valid: 1) The correlation between the incident rod strain gauge signal and the transmitted rod strain gauge signal and the piezoelectric crystal force value signal and the calibration strain gauge signal is not less than a preset value, such as 0.9; 2) The platform value deviation does not exceed 5%; 3) The total pulse width deviation does not exceed 5%.

[0042] If the incident rod strain gauge signal and the transmitted rod strain gauge signal do not meet the above three conditions at the same time, they are judged invalid and the pasting of the strain gauge needs to be checked to ensure the accuracy and reliability of the strain gauge acquisition signal.

[0043] In addition, the piezoelectric crystal force signal is basically consistent with the force signal of the calibration strain gauge, so the calibration strain gauge can be omitted, that is, the dynamic test calibration sensor only uses the piezoelectric crystal force sensor.

[0044] Optionally, the present invention can also perform dynamic sensitivity measurement on the incident rod strain gauge and the transmission rod strain gauge according to the piezoelectric crystal force value signal, specifically based on the following formula: , where a represents the dynamic sensitivity, F load represents the force value of the piezoelectric crystal force sensor, A represents the cross-sectional area of ​​the incident rod or the transmission rod, E represents the elastic modulus of the incident rod or the transmission rod, and f represents the gain coefficient of the incident rod strain gauge or the transmission rod strain gauge. It can be understood that, usually, the sensitivity coefficient of the strain gauge is obtained by testing under low-speed static conditions, and the sensitivity coefficient under high-speed dynamic conditions will be different. The present invention can accurately obtain the sensitivity coefficient of the strain gauge under high-speed dynamic conditions.

[0045] Alternatively, if Fig.15 As shown, the material dynamic test data processing method also includes the following contents: Step S4: The true stress-plastic strain curve is converted based on the engineering stress-strain curve, and then the constitutive model is used to perform fitting and extrapolation analysis based on the true stress-plastic strain curve to obtain the plastic flow curve.

[0046] Specifically, the true stress-plastic strain curve is first converted based on the engineering stress-strain curve, and the conversion formula is: ,in, represents the true strain, represents the engineering strain, represents the true stress, represents the engineering stress, Represents plastic strain, and E represents elastic modulus. After obtaining the true stress-plastic strain curve, the data after the necking point is removed, and the data before necking are used by programming to use different simulation models, such as JC, swfit, Ludwick, Voce, etc., to uniformly fit and extrapolate multiple strain rate curves of the same material to a strain of 1. Among them, the same strain rate is generally taken from the middle curve of the three test curves, and the number of curve points is 100, such as Fig.16 After fitting the extrapolated plastic flow curve, the extrapolated curve can be directly applied to simulation software such as LS-DYNA and AQAQUS to simulate the dynamic mechanical properties of materials, thus improving the usability of HPKS test data.

[0047] In addition, if Fig.17 As shown, another embodiment of the present invention further provides a material dynamic test data processing system, preferably using the material dynamic test data processing method as described above, including: A waveform curve automatic alignment module is used to obtain the incident rod signal curve and the transmission rod signal curve, extract the incident wave curve and the reflected wave curve from the incident rod signal curve, extract the transmitted wave curve from the transmission rod signal curve, and automatically align the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width; A data filtering and smoothing processing module is used to filter and smooth the incident wave curve, the reflected wave curve and the transmitted wave curve; The engineering stress-strain curve generation module is used to calculate the engineering stress-strain curve based on the incident wave curve, the reflected wave curve and the transmitted wave curve.

[0048] It can be understood that the material dynamic test data processing system of this embodiment, after obtaining the incident rod signal curve and the transmission rod signal curve, first extracts the incident wave curve and the reflected wave curve from the incident rod signal curve, and extracts the transmitted wave curve from the transmission rod signal curve, and then automatically aligns the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width, and then automatically calculates the engineering stress-strain curve after filtering and smoothing the incident wave curve, the reflected wave curve and the transmitted wave curve. Since the pulse widths of the incident wave, the reflected wave and the transmitted wave are automatically aligned, the consistency of the data time series is guaranteed, thereby ensuring the accuracy of the engineering stress-strain curve, and can accurately reflect the real dynamic mechanical properties.

[0049] In addition, the material dynamic test data processing system also includes: The fitting and extrapolation analysis module is used to convert the true stress-plastic strain curve based on the engineering stress-strain curve, and then use the constitutive model to perform fitting and extrapolation analysis based on the true stress-plastic strain curve to obtain the plastic flow curve.

[0050] In addition, another embodiment of the present invention further provides an electronic device, including a processor and a memory, wherein the memory stores a computer program, and the processor executes the steps of the above method by calling the computer program stored in the memory.

[0051] In addition, another embodiment of the present invention further provides a computer-readable storage medium for storing a computer program for processing material dynamic test data, wherein the computer program executes the steps of the method described above when running on a computer.

[0052] Common computer-readable storage media include: floppy disks, flexible disks, hard disks, magnetic tapes, any other magnetic media, CD-ROMs, any other optical media, punch cards, paper tapes, any other physical media with patterns of holes, random access memory (RAM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), flash-erasable programmable read-only memory (FLASH-EPROM), any other memory chip or cartridge, or any other medium that can be read by a computer. Instructions can further be transmitted or received by a transmission medium. The term transmission medium may include any tangible or intangible medium that can be used to store, encode or carry instructions for execution by a machine, and includes digital or analog communication signals or intangible media that facilitate communication of the above instructions. Transmission media include coaxial cables, copper wire, and optical fiber, including the wires of a bus used to transmit a computer data signal.

[0053] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The schemes in the embodiments of the present application may be implemented in various computer languages, for example, object-oriented programming language Java and literal scripting language JavaScript, etc.

[0054] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0055] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0056] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0057] Although the preferred embodiments of the present application have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present application.

[0058] Obviously, those skilled in the art can make various changes and modifications to the present application without departing from the spirit and scope of the present application. Thus, if these modifications and variations of the present application fall within the scope of the claims of the present application and their equivalents, the present application is also intended to include these modifications and variations.

[0059] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A material dynamic test data processing method, characterized in that: Includes the following: Obtaining an incident rod signal curve and a transmission rod signal curve, extracting an incident wave curve and a reflected wave curve from the incident rod signal curve, extracting a transmitted wave curve from the transmission rod signal curve, and automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width; Filter and smooth the incident wave curve, reflected wave curve and transmitted wave curve; The engineering stress-strain curve is calculated based on the incident wave curve, reflected wave curve and transmitted wave curve; The process of extracting the incident wave curve and the reflected wave curve from the incident rod signal curve, extracting the transmitted wave curve from the transmitted rod signal curve, and automatically aligning the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width includes the following: Find the extreme points of the incident wave and the reflected wave in the incident rod signal curve, and use them as the characteristic points of the incident wave curve and the reflected wave curve, respectively; find the extreme point of the transmitted wave in the transmitted rod signal curve, and use it as the characteristic point of the transmitted wave curve; Based on the characteristic point of the incident wave curve or the reflected wave curve, traverse forward in the incident rod signal curve to find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of the point, and use it as the starting point of the incident wave waveform or the reflected wave waveform; then traverse backward based on the characteristic point to find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all less than the measured value of the point, and use it as the ending point of the incident wave waveform or the reflected wave waveform, and determine the complete incident wave waveform curve or the reflected wave waveform based on the starting point and the ending point; Based on the characteristic points of the transmission wave curve, traverse forward in the transmission rod signal curve, find the first data point whose measured value is less than the preset value and the measured values ​​of subsequent data points are all greater than the measured value of this point, use it as the starting point of the transmission wave waveform, and determine the end point based on the same pulse width as the incident wave waveform, so as to determine the complete transmission wave waveform curve; The starting points of the incident wave waveform curve, the reflected wave waveform curve and the transmitted wave waveform curve are aligned, and the data with redundant time length in the reflected wave waveform curve are eliminated based on the principle of the same pulse width.

2. The material dynamic test data processing method according to claim 1, characterized in that: After obtaining the incident rod signal curve, the following contents are also included: The test type is determined according to the incident wave signal value. If the incident wave signal value is positive, it is determined to be a tensile test. If the incident wave signal value is negative, it is determined to be a compression test.

3. The material dynamic test data processing method according to claim 1, characterized in that: After obtaining the incident rod signal curve, the following contents are also included: Determine the validity of the incident rod signal and the transmitted rod signal.

4. The material dynamic test data processing method according to claim 3, characterized in that: The process of judging the validity of the incident rod signal includes the following: Obtain the theoretical incident rod signal curve, calculate the ratio of the difference between the maximum value of the measured incident rod signal and the maximum value of the theoretical incident rod signal to the maximum value of the theoretical incident rod signal, if the ratio is greater than a preset threshold, the incident rod signal is determined to be invalid, otherwise the incident rod signal is determined to be valid; calculate the correlation between the measured incident rod signal curve and the theoretical incident rod signal curve, if the correlation between the two is greater than a preset threshold, the incident rod signal is determined to be valid, otherwise the incident rod signal is determined to be invalid.

5. The material dynamic test data processing method according to claim 1, characterized in that: Also included are the following: The calibration signal curve of the dynamic test calibration sensor is obtained, and the incident rod signal curve and the transmission rod signal curve are compared with the calibration signal curve respectively. The validity of the incident rod strain gauge signal and the transmission rod strain gauge signal are judged according to the comparison results.

6. The material dynamic test data processing method according to claim 1, characterized in that: Also included are the following: The true stress-plastic strain curve is converted based on the engineering stress-strain curve, and then the constitutive model is used to perform fitting and extrapolation analysis based on the true stress-plastic strain curve to obtain the plastic flow curve.

7. A material dynamic test data processing system, using the material dynamic test data processing method according to any one of claims 1 to 6, characterized in that: include: A waveform curve automatic alignment module is used to obtain the incident rod signal curve and the transmission rod signal curve, extract the incident wave curve and the reflected wave curve from the incident rod signal curve, extract the transmitted wave curve from the transmission rod signal curve, and automatically align the incident wave curve, the reflected wave curve and the transmitted wave curve based on the pulse width; A data filtering and smoothing processing module is used to filter and smooth the incident wave curve, the reflected wave curve and the transmitted wave curve; The engineering stress-strain curve generation module is used to calculate the engineering stress-strain curve based on the incident wave curve, the reflected wave curve and the transmitted wave curve.

8. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the memory stores a computer program, and the processor executes the steps of the method according to any one of claims 1 to 6 by calling the computer program stored in the memory.

9. A computer-readable storage medium for storing a computer program for processing material dynamic test data, characterized in that: When the computer program is run on a computer, the computer program executes the steps of the method according to any one of claims 1 to 6.

Images