Amplitude acquisition method, device and equipment of crystallizer and storage medium

Abstract

The invention relates to an amplitude acquisition method, device and equipment for a crystallizer and a storage medium, and the method comprises the steps: carrying out the band-pass filtering processing of collected acceleration data of the crystallizer, and obtaining the filtering acceleration data; performing double-inversion integration processing on the filtered acceleration data to obtain displacement data; based on the vibration periodic boundary condition of the crystallizer, boundary correction is conducted on the displacement data, and corrected displacement data are obtained; calculating amplitude data of the crystallizer based on the correction displacement data; compared with the prior art, according to the technical scheme, high-precision amplitude calculation is achieved through band-pass filtering denoising, double-inversion integral trend term drift suppression and periodic boundary correction deviation correction.

Description

Technical Field

[0001] This application relates to the technical field of the iron and steel metallurgical industry, and in particular to a method, apparatus, equipment and storage medium for obtaining the amplitude of a crystallizer. Background Technology

[0002] As the core equipment in continuous casting production, the vibration state of the crystallizer directly affects the surface quality of the cast billet and the stability of the continuous casting process. Therefore, high-precision monitoring of crystallizer amplitude data is a key link to ensure the efficient operation of continuous casting production. In the existing technology, the core idea of ​​obtaining crystallizer amplitude data is to filter and integrate the acceleration data collected during the crystallizer vibration process, and obtain the amplitude data through the conversion relationship of acceleration-velocity-displacement. However, this basic processing method has significant defects: the amplitude data after direct processing has large fluctuations and poor waveform reproduction, which cannot meet the actual needs of high-quality monitoring of the crystallizer vibration state. The data processing process needs to be optimized and improved.

[0003] Furthermore, to address the aforementioned issues, current methods generally involve processing the acceleration data of the crystallizer, such as using multiple rounds of filtering to optimize the amplitude results, or proposing a segmented calculation of the average trend value to reduce the influence of the trend term.

[0004] However, existing technical solutions for processing crystallizer acceleration data to obtain high-precision amplitude generally suffer from problems such as incomplete trend term processing, complex data processing procedures, and insufficient amplitude calculation accuracy. These solutions cannot fully meet the needs of the iron and steel metallurgical industry for high-quality monitoring of crystallizer vibration status. There is an urgent need for a technical solution that can effectively suppress trend term drift, simplify the processing procedure, and ensure high-precision amplitude results. Summary of the Invention

[0005] This application provides a method, apparatus, device, and storage medium for obtaining the amplitude of a crystallizer. High-precision amplitude calculation is achieved through bandpass filtering for noise reduction, double inverted integral to suppress trend term drift, and periodic boundary correction to correct deviation.

[0006] In a first aspect, this application provides a method for obtaining the amplitude of a crystallizer, comprising: performing bandpass filtering on the collected acceleration data of the crystallizer to obtain filtered acceleration data; performing double inverted integral processing on the filtered acceleration data to obtain displacement data; performing boundary correction on the displacement data based on the vibration periodic boundary conditions of the crystallizer to obtain corrected displacement data; and calculating the amplitude data of the crystallizer based on the corrected displacement data.

[0007] In one possible implementation, the step of performing double inverted integration on the filtered acceleration data to obtain displacement data specifically includes: performing a first integration on the filtered acceleration data to obtain velocity data; performing a first inversion on the velocity data to obtain inverted velocity data; performing a second integration on the inverted velocity data to obtain preliminary displacement data; and performing a second inversion on the preliminary displacement data to obtain displacement data.

[0008] In one possible implementation, the step of performing boundary correction on the displacement data based on the periodic boundary conditions of the crystallizer to obtain corrected displacement data specifically includes: obtaining the initial and final displacement values ​​corresponding to the displacement data based on the periodic boundary conditions of the crystallizer; constructing a correction linear function based on the initial and final displacement values; and subtracting the correction linear function from the displacement data to obtain the corrected displacement data.

[0009] In one possible implementation, the step of obtaining the initial and final displacement values ​​corresponding to the displacement data based on the periodic boundary conditions of the crystallizer's vibration; constructing a correction linear function based on the initial and final displacement values; and subtracting the correction linear function from the displacement data to obtain corrected displacement data specifically includes: dividing the displacement data into displacement signal quantities and displacement deviation quantities; obtaining the initial and final displacement signal values ​​corresponding to the displacement signal quantities based on the periodic boundary conditions of the crystallizer's vibration; and constructing a correction linear function based on the initial and final displacement values; and subtracting the correction linear function from the displacement data to obtain corrected displacement data. A linear correction function for displacement signal quantity is constructed using the displacement signal value at the final moment; the linear correction function is subtracted from the displacement signal quantity to obtain the corrected displacement signal quantity; based on the periodic boundary conditions of the crystallizer's vibration, the initial moment displacement deviation value and the final moment displacement deviation value corresponding to the displacement deviation amount are obtained; a linear correction function for displacement deviation amount is constructed based on the initial moment displacement deviation value and the final moment displacement deviation value; the linear correction function is subtracted from the displacement deviation amount to obtain the corrected displacement deviation amount; the corrected displacement signal quantity and the corrected displacement deviation amount are integrated to obtain the corrected displacement data.

[0010] In one possible implementation, after obtaining the corrected displacement data, the method further includes: performing standard first-order integration on the filtered acceleration data to obtain standard velocity data, and performing standard second-order integration on the standard velocity data to obtain standard displacement data; performing boundary correction on the standard displacement data based on the periodic boundary conditions of the crystallizer to obtain standard corrected displacement data; and quantifying the error accuracy of the corrected displacement data based on the standard corrected displacement data and the corrected displacement data.

[0011] In one possible implementation, before performing bandpass filtering on the collected acceleration data of the crystallizer, the method further includes: collecting the acceleration data of the crystallizer based on an accelerometer, wherein the accelerometer is fixed on the vibration table of the crystallizer and vibrates synchronously with the vibration table; and uploading the collected acceleration data to a host computer based on a wireless transmission module.

[0012] In one possible implementation, calculating the amplitude data of the crystallizer based on the corrected displacement data specifically includes: obtaining the peak and valley values ​​of the corrected displacement data within a single vibration cycle, and calculating the amplitude data of the crystallizer based on the peak and valley values.

[0013] Secondly, this application provides an amplitude acquisition device for a crystallizer, comprising: a bandpass filtering module, a double inverted integral module, a correction module, and an amplitude data calculation module; wherein, the bandpass filtering module is used to perform bandpass filtering on the collected acceleration data of the crystallizer to obtain filtered acceleration data; the double inverted integral module is used to perform double inverted integral processing on the filtered acceleration data to obtain displacement data; the correction module is used to perform boundary correction on the displacement data based on the periodic boundary conditions of the crystallizer's vibration to obtain corrected displacement data; and the amplitude data calculation module is used to calculate the amplitude data of the crystallizer based on the corrected displacement data.

[0014] Thirdly, embodiments of this application also provide a computer device, which includes a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the above-described method.

[0015] Fourthly, embodiments of this application also provide a computer-readable storage medium storing a computer program that, when executed by a processor, can implement the above-described method.

[0016] This application provides a method, apparatus, device, and storage medium for obtaining the amplitude of a crystallizer, which has the following advantages compared with the prior art:

[0017] This method involves bandpass filtering the collected acceleration data from the crystallizer to obtain filtered acceleration data; performing double inverted integration on the filtered acceleration data to obtain displacement data; performing boundary correction on the displacement data based on the periodic boundary conditions of the crystallizer's vibration to obtain corrected displacement data; and calculating the crystallizer's amplitude data based on the corrected displacement data. Compared with existing technologies, the technical solution of this application can accurately select the effective vibration frequency band in the acceleration data through bandpass filtering, filtering out irrelevant noise interference and laying a high-quality foundation for subsequent data processing; the double inverted integration can effectively offset the trend generated during the integration process. The term drift avoids the error accumulation caused by traditional single integration, significantly improving the accuracy of displacement data. Combined with boundary correction of vibration periodic boundary conditions, the deviation of displacement data at periodic nodes can be further corrected, ensuring a high degree of fidelity between the displacement waveform and the actual vibration state of the crystallizer. Finally, the amplitude is calculated based on the corrected displacement data. The simplified operation steps are achieved through a simple process of filtering-integration-correction. Furthermore, the multi-stage error control solves the problems of incomplete trend term processing, complex process, and insufficient accuracy in the existing technology. This achieves efficient acquisition of high-precision amplitude results, which can fully meet the needs of the iron and steel metallurgical industry for high-quality monitoring of the crystallizer vibration state. Attached Figure Description

[0018] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, for those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0020] One or more embodiments are illustrated by way of example with reference numerals in the accompanying drawings. These illustrations do not constitute a limitation on the embodiments. Elements with the same reference numerals in the drawings are denoted as similar elements. Unless otherwise stated, the figures in the drawings are not to be limited by scale.

[0021] Figure 1 This is a flowchart illustrating one embodiment of a method for obtaining the amplitude of a crystallizer provided in this application;

[0022] Figure 2 This is a schematic diagram of one embodiment of the amplitude acquisition device for a crystallizer provided in this application;

[0023] Figure 3 This is a schematic diagram of the drift of the integral displacement curve under standard double integration;

[0024] Figure 4 This is a schematic diagram showing the relationship between the standard integral and the double inverted integral correction signal.

[0025] Figure 5 This is a schematic diagram of the structure of a computer device provided in this application. Detailed Implementation

[0026] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0027] The following disclosure provides numerous different embodiments or examples for implementing various structures of this application. To simplify the disclosure, specific examples of components and arrangements are described below. These are merely examples and are not intended to limit the scope of this application. Furthermore, reference numerals and / or letters may be repeated in different examples. Such repetition is for simplification and clarity and does not in itself indicate a relationship between the various embodiments and / or arrangements discussed.

[0028] It should be understood that, when used in this specification and the appended claims, the terms "comprising" and "including" indicate the presence of the described features, integrals, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components and / or collections thereof.

[0029] It should also be understood that the terminology used in this specification is for the purpose of describing particular embodiments only and is not intended to limit the scope of the application. As used in this specification and the appended claims, the singular forms “a,” “an,” and “the” are intended to include the plural forms unless the context clearly indicates otherwise.

[0030] It should also be further understood that the term “and / or” as used in this application specification and the appended claims means any combination of one or more of the associated listed items and all possible combinations, and includes such combinations.

[0031] As used in this specification and the appended claims, the term "if" may be interpreted, depending on the context, as "when," "once," "in response to determination," or "in response to detection." Similarly, the phrases "if determined" or "if [described condition or event] is detected" may be interpreted, depending on the context, as "once determined," "in response to determination," "once [described condition or event] is detected," or "in response to detection of [described condition or event]."

[0032] Example 1, see Figure 1 , Figure 1 This is a schematic flowchart of an embodiment of a method for obtaining the amplitude of a crystallizer provided in this application, as shown below. Figure 1 As shown, the method includes steps 101-104, as detailed below:

[0033] Step 101: Perform bandpass filtering on the collected acceleration data of the crystallizer to obtain filtered acceleration data.

[0034] In one embodiment, acceleration data of the crystallizer is collected based on an accelerometer, wherein the accelerometer is fixed on the vibration table of the crystallizer and vibrates synchronously with the vibration table; the collected acceleration data is uploaded to a host computer based on a wireless transmission module.

[0035] Specifically, the acceleration sensor includes, but is not limited to, a piezoelectric acceleration sensor; wherein, the piezoelectric acceleration sensor has high sensitivity to low-frequency vibrations and can accurately capture minute acceleration changes of the vibration table, providing high-quality raw data for subsequent amplitude calculations.

[0036] Specifically, the accelerometer sensor is fixed using a neodymium magnet. By mounting a neodymium magnet at the bottom of the accelerometer sensor, the sensor is fixed to the crystallizer vibration table through strong magnetic attraction. On the one hand, the attraction strength of the neodymium magnet can be adapted to the vibration intensity of the crystallizer vibration table, ensuring that the sensor maintains absolute synchronization with the vibration table during vibration, avoiding data distortion caused by relative displacement due to loose installation. On the other hand, the magnetic attraction method is detachable, which facilitates the installation, debugging, periodic calibration, maintenance and replacement of the sensor, and is suitable for the needs of rotating monitoring of multiple crystallizers in industrial scenarios.

[0037] Specifically, the acceleration sensor belongs to the data acquisition system, which further includes a wireless transmission module; preferably, the wireless transmission module includes, but is not limited to, a WIFI module.

[0038] Preferably, the data acquisition system further includes a FLASH storage chip, an MCU, a WIFI module, a battery, and a charge / discharge management circuit.

[0039] In one embodiment, bandpass filtering is performed on the collected acceleration data of the crystallizer to extract effective vibration signals from the original acceleration data and suppress irrelevant interference signals, thus laying a high-quality data foundation for subsequent integral calculation and amplitude acquisition.

[0040] In one embodiment, the frequency range of the bandpass filter is set to 1-50Hz. This is because the vibration frequency of the crystallizer vibration table is generally below 6Hz. Setting the frequency range to 1-50Hz can prevent low-frequency drift signals below 1Hz from entering subsequent processing. The upper limit of the frequency range is much higher than the highest vibration frequency of the crystallizer, which can completely cover the fundamental frequency and harmonic components related to the crystallizer vibration, ensuring that the effective vibration signal is not filtered out and avoiding signal truncation distortion caused by the filter range being too narrow.

[0041] In one embodiment, the collected acceleration data from the crystallizer is bandpass filtered to obtain filtered acceleration data, as shown below:

[0042] a(t) = A sin(ωt+φ) + σ(t);

[0043] In the formula, a(t) is the filtered acceleration data at time t, φ is the initial phase, ω=2πf is the angular velocity, ε(t) is the deviation and noise, and A is the acceleration amplitude.

[0044] Since high-frequency noise is suppressed after bandpass filtering from 1 to 50 Hz, ε(t) mainly manifests as low-frequency deviation caused by sensor manufacturing process and filter characteristics. Considering the limiting case, i.e., limε(t)=ε0, the drift caused by the low-frequency deviation on the integral is the largest. At this time, the obtained filtered acceleration data can be expressed as:

[0045] a(t) = A sin(ωt+φ) + ε0;

[0046] In the formula, ε0 is the initial low-frequency deviation constant remaining after the acceleration data is bandpass filtered, and A sin(ωt+φ) is the effective acceleration data of the crystallizer vibration.

[0047] Step 102: Perform double inverted integration on the filtered acceleration data to obtain displacement data.

[0048] In one embodiment, the filtered acceleration data is integrated once to obtain velocity data.

[0049] Specifically, the filtered acceleration data a(t) is integrated from the initial time 0 to any time t to obtain the velocity data v(t), as shown in the following process:

[0050]

[0051] In the formula, ε0(t) is the deviation term at time t. This is an effective speed signal.

[0052] Specifically, by integrating the filtered acceleration data once, the physical quantity conversion from acceleration to velocity is achieved.

[0053] In one embodiment, the speed data is first inverted to obtain inverted speed data.

[0054] Specifically, the inversion process refers to the time axis reversal, which replaces the variable t in the velocity data that changes with time with the acquisition duration (Tt). That is, through the mapping from v(t) to v(Tt), the inverted velocity data u(t) is obtained.

[0055] Specifically, during the first inversion of the velocity data, t in v(t) is replaced with Tt, and substituted into the calculation formula for the velocity data to obtain the inverted velocity data u(t), as shown below:

[0056]

[0057] Specifically, by reversing the time in the first step, the time dependency of the deviation term is adjusted, so that after the subsequent second integration, the form of the deviation term changes from a quadratic function to a form that can be further suppressed by boundary correction.

[0058] In one embodiment, the inverted velocity data is integrated twice to obtain preliminary displacement data.

[0059] Specifically, by performing a second integration on the inverted velocity data, the final physical quantity conversion from velocity to displacement is achieved, resulting in preliminary displacement data without undergoing a second inversion and boundary correction.

[0060] Specifically, the inverted velocity data u(t) is integrated from the initial time 0 to any time t to obtain the preliminary displacement data w(t), as shown below:

[0061]

[0062] Specifically, after obtaining the preliminary displacement data w(t), the calculation formula for the preliminary displacement data w(t) is further integrated and decomposed according to the effective signal term and the deviation term, resulting in:

[0063] w(t) = w1(t) + w2(t);

[0064]

[0065] In the formula, w1(t) is the effective signal term of the preliminary displacement data, and w2(t) is the deviation term of the preliminary displacement data.

[0066] Specifically, for the effective signal term w1(t) of the initial displacement data: Let θ = ωt + φ, then cos(ω(T-τ) + φ) = cos(θ-ωτ), we can obtain:

[0067]

[0068] Specifically, for the deviation term w2(t) of the initial displacement data, the following can be calculated using the constant and polynomial integration rules:

[0069]

[0070] Specifically, by integrating the effective signal term w1(t) and the deviation term w2(t) of the preliminary displacement data, the final expression for the preliminary displacement data w(t) is obtained, as shown below:

[0071]

[0072] Specifically, the effective signal term w1(t) of the preliminary displacement data w(t) has initially reflected the displacement change of the crystallizer vibration, but there is still a phase shift due to time inversion; the deviation term w2(t) of the preliminary displacement data it contains is a quadratic function drift term generated by the second integral of the low-frequency deviation ε0. The deviation needs to be further adjusted through a subsequent second inversion.

[0073] In one embodiment, the preliminary displacement data is inverted a second time to obtain the displacement data.

[0074] Specifically, the second inversion is logically consistent with the first inversion, and it is still a time axis reversal: the time variable t in the initial displacement data is replaced with Tt, that is, the final displacement data s(t) is obtained through the mapping from w(t) to w(Tt); the two inversions form a closed loop of double inversion integral of acceleration integral - velocity inversion - velocity integral - displacement inversion.

[0075] Specifically, when performing the second inversion process on the initial displacement data w(t), t in w(t) is replaced with Tt, and this Tt is substituted into the final expression of the initial displacement data w(t) to obtain the displacement data s(t), as shown below:

[0076]

[0077] Specifically, because the acceleration data collected by the accelerometer contains several cycles, therefore Given sin(ωT+φ)=sinφ; cos(ωT+φ)=cosφ; therefore, simplifying the above displacement data s(t), we obtain the final expression for the displacement data s(t), as shown below:

[0078]

[0079] Step 103: Based on the vibration periodicity boundary conditions of the crystallizer, perform boundary correction on the displacement data to obtain corrected displacement data.

[0080] In one embodiment, based on the periodic boundary conditions of the crystallizer's vibration, the initial displacement value and the final displacement value corresponding to the displacement data are obtained; based on the initial displacement value and the final displacement value, a correction linear function is constructed; the correction linear function is subtracted from the displacement data to obtain the correction displacement data.

[0081] Specifically, the displacement data is divided into displacement signal quantity and displacement deviation quantity; based on the periodic boundary conditions of the crystallizer's vibration, the initial and final displacement signal values ​​corresponding to the displacement signal quantity are obtained; based on the initial and final displacement signal values, a displacement signal quantity correction linear function is constructed; the displacement signal quantity is subtracted from the displacement signal quantity correction linear function to obtain the corrected displacement signal quantity; based on the periodic boundary conditions of the crystallizer's vibration, the initial and final displacement deviation values ​​corresponding to the displacement deviation quantity are obtained; based on the initial and final displacement deviation values, a displacement deviation quantity correction linear function is constructed; the displacement deviation quantity is subtracted from the displacement deviation quantity correction linear function to obtain the corrected displacement deviation quantity; the corrected displacement signal quantity and the corrected displacement deviation quantity are integrated to obtain the corrected displacement data.

[0082] Specifically, the displacement data is divided into displacement signal quantity and displacement deviation quantity, as shown below:

[0083] s(t)=s signal (t)+s ε (t);

[0084]

[0085] In the formula, s signal (t) represents the displacement signal quantity, s ε (t) represents the displacement deviation.

[0086] Specifically, since the vibration of the crystallizer's vibration table is periodic, meaning the initial and final displacements are zero (i.e., s(0) = 0 and s(T) = 0), the s in s(t) can be corrected using forced boundary conditions. signal(t) and s ε Subtract the corresponding linear function from each (t) to satisfy the boundary conditions.

[0087] Specifically, let the corrected displacement signal be S. signal_corr (t):

[0088]

[0089] S signal (T) = 0;

[0090] Based on the initial displacement signal value and the final displacement signal value, the linear function for displacement signal correction can be obtained as follows:

[0091]

[0092] The corrected displacement signal is:

[0093]

[0094] In the above formula, s signal (0) represents the initial displacement signal value, s signal (T) represents the displacement signal value at the final moment.

[0095] Specifically, let the corrected displacement deviation be s. ε_corr (t):

[0096]

[0097] s ε (T) = 0;

[0098] Based on the initial displacement deviation value and the final displacement deviation value, a linear function for displacement deviation correction is constructed as follows:

[0099]

[0100] The corrected displacement deviation is:

[0101]

[0102] In the above formula, s ε (0) represents the displacement deviation value at the first initial moment, s ε (T) represents the displacement deviation value at the final moment.

[0103] Specifically, the corrected displacement signal and the corrected displacement deviation are integrated to obtain the corrected displacement data, as shown below:

[0104]

[0105] Step 104: Calculate the amplitude data of the crystallizer based on the corrected displacement data.

[0106] In one embodiment, the peak and valley values ​​of the corrected displacement data within a single vibration cycle are obtained, and the amplitude data of the crystallizer is calculated based on the peak and valley values.

[0107] Specifically, the average value between the peak value and the valley value is calculated, and the average value is used as the amplitude data of the crystallizer.

[0108] In one embodiment, the filtered acceleration data is further subjected to standard first-order integration to obtain standard velocity data, and the standard velocity data is subjected to standard second-order integration to obtain standard displacement data; based on the vibration periodic boundary conditions of the crystallizer, the standard displacement data is subjected to boundary correction to obtain standard corrected displacement data; based on the standard corrected displacement data and the corrected displacement data, the error accuracy of the corrected displacement data is quantified.

[0109] Specifically, after performing two integrations on the filtered acceleration data, the standard displacement data x(t) is obtained as follows:

[0110]

[0111] Specifically, based on the periodic boundary conditions of the crystallizer's vibration, the standard displacement data is boundary-corrected to obtain standard corrected displacement data. Similarly, the standard displacement data is divided into standard displacement signal quantity and standard displacement deviation quantity. Based on the periodic boundary conditions of the crystallizer's vibration, the standard initial moment displacement signal value and the standard final moment displacement signal value corresponding to the standard displacement signal quantity are obtained. Based on the standard initial moment displacement signal value and the standard final moment displacement signal value, a standard displacement signal quantity correction linear function is constructed. The standard displacement signal quantity is subtracted from the standard displacement signal quantity correction linear function to obtain the standard corrected displacement signal quantity. Based on the periodic boundary conditions of the crystallizer's vibration, the standard initial moment displacement deviation value and the standard final moment displacement deviation value corresponding to the standard displacement deviation quantity are obtained. Based on the standard initial moment displacement deviation value and the standard final moment displacement deviation value, a standard displacement deviation quantity correction linear function is constructed. The standard displacement deviation quantity is subtracted from the standard displacement deviation quantity correction linear function to obtain the standard corrected displacement deviation quantity. The standard corrected displacement signal quantity and the standard corrected displacement deviation quantity are integrated to obtain the standard corrected displacement data.

[0112] Specifically, the standard displacement data is divided into standard displacement signal quantities x. signal (t) and standard displacement deviation x ε(t), as shown below:

[0113] x(t)=x signal (t)+x ε (t);

[0114]

[0115] Specifically, for the standard displacement signal quantity x in the standard displacement data x(t) signal (t) Correction

[0116] x signal (0) = 0;

[0117]

[0118] Based on the standard initial displacement signal value and the standard final displacement signal value, a standard displacement signal correction linear function is constructed as follows:

[0119]

[0120] Obtain the standard correction displacement signal quantity x signal_corr (t):

[0121]

[0122] Based on the above corrected displacement signal quantity, it can be known that: x signal_corr (t)=-s signal_corr (t); This indicates that after boundary condition correction of the signals x(t) and s(t), their amplitudes are equal. The waveform of x(t) is in phase with the standard displacement integral signal when there is no deviation, while s(t) is out of phase with it. The main difference is caused by the deviation.

[0123] s ε_corr (t) and x ε (t) compare the maximum values ​​corresponding to each, where the standard displacement deviation x under the standard integral is... ε (t): When t = T, the deviation term is at its maximum. For the corrected displacement deviation s under the double inverted integral method ε_corr (t): When The maximum deviation term is Based on the above comparison, it can be seen that the double inverted integral method reduces the maximum error from Reduce to It has been reduced by 75%.

[0124] To verify the above theory, assume that the data acquisition system performs continuous data acquisition for 5 seconds at a sampling rate of Fs = 4000Hz. The expression for the acceleration data a(t) of the crystallizer after bandpass filtering is as follows:

[0125]

[0126] That is, the maximum amplitude A = 0.003m, the vibration frequency f = 3Hz, and the initial phase ε0 is a constant deviation.

[0127] Simultaneously, the acceleration data a(t) is subjected to standard quadratic integration and correction to obtain standard corrected displacement data x(t), where t∈[0,5]:

[0128]

[0129] When ε0 = 0.0000003, it will cause the integral displacement curve to drift, and the amount of drift will increase as time continues to increase. Figure 3 As shown, Figure 3 This is a schematic diagram of the drift of the integral displacement curve under standard double integration.

[0130] The acceleration data a(t) is then subjected to double inverted integration and correction to obtain the corrected displacement data s(t), where t∈[0,5]:

[0131]

[0132] For the s signal corresponding to s(t) respectively signal_corr (t) and the semaphore x corresponding to x(t) signal_corr (t) Perform graph comparison, such as Figure 4 As shown, Figure 4 This is a schematic diagram showing the relationship between the standard integral and the double inverted integral corrected signal quantities; based on the diagram, it can be verified that the amplitudes of the two signal quantities are equal and their phases are opposite.

[0133] The maximum value of the deviation between the standard second integral and the double inverted integral is then verified. Let the unbiased acceleration be a′(t), as shown below:

[0134]

[0135] After performing a standard quadratic integral and correcting it, we get x. ′ (t), the expression is as follows:

[0136]

[0137] Based on the theoretical analysis of the above maximum deviation, the maximum value of the standard displacement integral deviation can be determined. and the maximum value of the integral deviation of the double inverted displacement The expression is as follows:

[0138]

[0139] Easy to obtain:

[0140]

[0141] This ratio directly proves that the maximum deviation of the double inverted integral method is only 25% of that of the standard double integral method, a reduction of 75%. Therefore, it is believed that the double inverted integral can reduce the influence of the deviation term.

[0142] Example 2, see Figure 2 , Figure 2 This is a schematic diagram of an embodiment of the crystallizer amplitude acquisition device provided in this application. Corresponding to the above-described crystallizer amplitude acquisition method, this application also provides a crystallizer amplitude acquisition device. The crystallizer amplitude acquisition device includes modules for executing the above-described crystallizer amplitude acquisition method, and can be configured in a desktop computer, tablet computer, laptop computer, or other terminal. Specifically, the crystallizer amplitude acquisition device includes a bandpass filter module 201, a double inverted integral module 202, a correction module 203, and an amplitude data calculation module 204.

[0143] The bandpass filter module 201 is used to perform bandpass filtering on the collected acceleration data of the crystallizer to obtain filtered acceleration data.

[0144] The dual inverted integral module 202 is used to perform dual inverted integral processing on the filtered acceleration data to obtain displacement data.

[0145] The correction module 203 is used to perform boundary correction on the displacement data based on the vibration periodic boundary conditions of the crystallizer to obtain corrected displacement data.

[0146] The amplitude data calculation module 204 is used to calculate the amplitude data of the crystallizer based on the corrected displacement data.

[0147] In one embodiment, the dual inverted integration module 202 is used to perform dual inverted integration processing on the filtered acceleration data to obtain displacement data. Specifically, it includes: performing a first integration on the filtered acceleration data to obtain velocity data; performing a first inversion processing on the velocity data to obtain inverted velocity data; performing a second integration on the inverted velocity data to obtain preliminary displacement data; and performing a second inversion processing on the preliminary displacement data to obtain displacement data.

[0148] In one embodiment, the correction module 203 is used to perform boundary correction on the displacement data based on the periodic boundary conditions of the crystallizer to obtain corrected displacement data. Specifically, it includes: obtaining the initial displacement value and the final displacement value corresponding to the displacement data based on the periodic boundary conditions of the crystallizer; constructing a correction linear function based on the initial displacement value and the final displacement value; and subtracting the correction linear function from the displacement data to obtain the corrected displacement data.

[0149] In one embodiment, the correction module 203 is used to obtain the initial displacement value and the final displacement value corresponding to the displacement data based on the periodic boundary conditions of the crystallizer's vibration; construct a correction linear function based on the initial displacement value and the final displacement value; subtract the correction linear function from the displacement data to obtain corrected displacement data. Specifically, obtaining the corrected displacement data includes: dividing the displacement data into displacement signal quantity and displacement deviation quantity; obtaining the initial displacement signal value and the final displacement signal value corresponding to the displacement signal quantity based on the periodic boundary conditions of the crystallizer's vibration; and subtracting the correction linear function from the initial displacement data to obtain corrected displacement data. Based on the initial displacement signal value and the final displacement signal value, a linear correction function for the displacement signal is constructed. The linear correction function is subtracted from the displacement signal to obtain the corrected displacement signal. Based on the periodic boundary conditions of the crystallizer's vibration, the initial displacement deviation value and the final displacement deviation value corresponding to the displacement deviation are obtained. A linear correction function for the displacement deviation is constructed based on the initial displacement deviation value and the final displacement deviation value. The linear correction function is subtracted from the displacement deviation to obtain the corrected displacement deviation. The corrected displacement signal and the corrected displacement deviation are integrated to obtain the corrected displacement data.

[0150] In one embodiment, the crystallizer amplitude acquisition device provided in this application further includes: an error accuracy quantization module.

[0151] In one embodiment, the error accuracy quantization module is used to perform standard first-order integration processing on the filtered acceleration data to obtain standard velocity data, and perform standard second-order integration processing on the standard velocity data to obtain standard displacement data; based on the vibration periodic boundary conditions of the crystallizer, perform boundary correction on the standard displacement data to obtain standard corrected displacement data; and quantify the error accuracy of the corrected displacement data based on the standard corrected displacement data and the corrected displacement data.

[0152] In one embodiment, the crystallizer amplitude acquisition device provided in this application further includes a data acquisition module.

[0153] In one embodiment, the data acquisition module is used to collect acceleration data of the crystallizer based on an accelerometer, wherein the accelerometer is fixed on the vibration table of the crystallizer and vibrates synchronously with the vibration table; and the collected acceleration data is uploaded to a host computer based on a wireless transmission module.

[0154] In one embodiment, the amplitude data calculation module 204 is used to calculate the amplitude data of the crystallizer based on the corrected displacement data, specifically including: obtaining the peak value and valley value of the corrected displacement data in a single vibration cycle, and calculating the amplitude data of the crystallizer based on the peak value and the valley value.

[0155] The above-described crystallizer amplitude acquisition device can implement the crystallizer amplitude acquisition method of the above method embodiment. The options in the above method embodiment are also applicable to this embodiment, and will not be detailed here.

[0156] like Figure 5 As shown, Figure 5 This is a schematic diagram of the structure of a computer device provided in this application; it includes a processor 111, a communication interface 112, a memory 113 and a communication bus 114, wherein the processor 111, the communication interface 112 and the memory 113 communicate with each other through the communication bus 114, and the memory 113 is used to store computer programs.

[0157] In one embodiment of this application, the processor 111, when executing the program stored in the memory 113, implements the crystallizer amplitude acquisition method provided in any of the foregoing method embodiments.

[0158] It will be understood by those skilled in the art that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program may be stored in a storage medium, which is a computer-readable storage medium. The computer program is executed by at least one processor in the computer system to implement the process steps of the embodiments of the above methods.

[0159] Therefore, embodiments of this application also provide a computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the steps of the crystallizer amplitude acquisition method provided in any of the foregoing method embodiments.

[0160] The storage medium is a physical, non-transient storage medium, such as a USB flash drive, external hard drive, read-only memory (ROM), magnetic disk, or optical disk, or any other physical storage medium capable of storing program code. The computer-readable storage medium can be non-volatile or volatile.

[0161] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this application.

[0162] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For example, the division of each unit is merely a logical functional division, and there may be other division methods in actual implementation. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed.

[0163] The steps in the methods of this application embodiment can be adjusted, merged, or deleted according to actual needs. The units in the apparatus of this application embodiment can be merged, divided, or deleted according to actual needs. Furthermore, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0164] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, a terminal, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application.

[0165] In the above embodiments, the descriptions of each embodiment have different focuses. For parts that are not described in detail in a certain embodiment, please refer to the relevant descriptions in other embodiments.

[0166] Obviously, those skilled in the art can make various modifications and variations to this application without departing from the spirit and scope of this application. Since these modifications and variations fall within the scope of the claims and their equivalents, this application also intends to include these modifications and variations.

[0167] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any person skilled in the art can easily conceive of various equivalent modifications or substitutions within the technical scope disclosed in this application, and these modifications or substitutions should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for obtaining the amplitude of a crystallizer, characterized in that, include: The collected acceleration data from the crystallizer is bandpass filtered to obtain filtered acceleration data. The filtered acceleration data is subjected to double inverted integration to obtain displacement data; Based on the vibrational periodicity boundary conditions of the crystallizer, the displacement data is boundary-corrected to obtain corrected displacement data; The amplitude data of the crystallizer is calculated based on the corrected displacement data.

2. The method as described in claim 1 above, characterized in that, The process of performing double inverted integral processing on the filtered acceleration data to obtain displacement data specifically includes: Integrate the filtered acceleration data once to obtain velocity data; The speed data is first inverted to obtain inverted speed data; By performing a second integration on the inverted velocity data, preliminary displacement data is obtained; The preliminary displacement data is then inverted a second time to obtain the final displacement data.

3. The method as described in claim 1 above, characterized in that, The displacement data is corrected based on the vibrational periodicity boundary conditions of the crystallizer to obtain corrected displacement data, specifically including: Based on the vibration periodic boundary conditions of the crystallizer, the initial displacement value and the final displacement value corresponding to the displacement data are obtained; Based on the initial displacement value and the final displacement value, a correction linear function is constructed; Subtract the correction linear function from the displacement data to obtain the corrected displacement data.

4. The method as described in claim 3 above, characterized in that, Based on the periodic boundary conditions of the crystallizer's vibration, the initial and final displacement values ​​corresponding to the displacement data are obtained; a correction linear function is constructed based on the initial and final displacement values. Subtracting the correction linear function from the displacement data yields the corrected displacement data, which specifically includes: The displacement data is divided into displacement signal quantity and displacement deviation quantity; Based on the periodic boundary conditions of the crystallizer, the initial displacement signal value and the final displacement signal value corresponding to the displacement signal quantity are obtained; Based on the initial displacement signal value and the final displacement signal value, a linear function for displacement signal correction is constructed. Subtract the displacement signal correction linear function from the displacement signal quantity to obtain the corrected displacement signal quantity; Based on the vibration periodic boundary conditions of the crystallizer, the initial displacement deviation value and the final displacement deviation value corresponding to the displacement deviation amount are obtained. Based on the initial displacement deviation value and the final displacement deviation value, a linear function for displacement deviation correction is constructed. Subtract the displacement deviation correction linear function from the displacement deviation to obtain the corrected displacement deviation. By integrating the corrected displacement signal and the corrected displacement deviation, corrected displacement data is obtained.

5. The method as described in claim 1, characterized in that, After obtaining the corrected displacement data, the process further includes: The filtered acceleration data is subjected to standard first-order integration to obtain standard velocity data, and the standard velocity data is subjected to standard second-order integration to obtain standard displacement data. Based on the vibration periodic boundary conditions of the crystallizer, the standard displacement data is boundary corrected to obtain standard corrected displacement data; Based on the standard corrected displacement data and the corrected displacement data, the error accuracy of the corrected displacement data is quantified.

6. The method as described in claim 1, characterized in that, Before performing bandpass filtering on the collected acceleration data of the crystallizer, the process further includes: The acceleration data of the crystallizer is collected based on an accelerometer sensor, wherein the accelerometer sensor is fixed on the vibration table of the crystallizer and vibrates synchronously with the vibration table. The collected acceleration data is uploaded to the host computer using a wireless transmission module.

7. The method as described in claim 1, characterized in that, The calculation of the crystallizer amplitude data based on the corrected displacement data specifically includes: The peak and valley values ​​of the corrected displacement data within a single vibration cycle are obtained, and the amplitude data of the crystallizer is calculated based on the peak and valley values.

8. An amplitude acquisition device for a crystallizer, characterized in that, include: Bandpass filter module, double inverted integral module, correction module and amplitude data calculation module; The bandpass filter module is used to perform bandpass filtering on the collected acceleration data of the crystallizer to obtain filtered acceleration data. The double inverted integral module is used to perform double inverted integration on the filtered acceleration data to obtain displacement data; The correction module is used to perform boundary correction on the displacement data based on the vibration periodic boundary conditions of the crystallizer to obtain corrected displacement data. The amplitude data calculation module is used to calculate the amplitude data of the crystallizer based on the corrected displacement data.

9. A computer device, characterized in that, The computer device includes a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the method as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The storage medium stores a computer program that, when executed by a processor, can implement the method as described in any one of claims 1-7.

Images