CN113144997A - Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing - Google Patents

Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing Download PDF

Info

Publication number
CN113144997A
CN113144997A CN202110308316.7A CN202110308316A CN113144997A CN 113144997 A CN113144997 A CN 113144997A CN 202110308316 A CN202110308316 A CN 202110308316A CN 113144997 A CN113144997 A CN 113144997A
Authority
CN
China
Prior art keywords
acceleration
fluctuation rate
mixing
sample
mixed
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110308316.7A
Other languages
Chinese (zh)
Other versions
CN113144997B (en
Inventor
马宁
陈松
孙晓朋
张哲�
谢中元
朱士富
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Modern Chemistry Research Institute
Original Assignee
Xian Modern Chemistry Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Modern Chemistry Research Institute filed Critical Xian Modern Chemistry Research Institute
Priority to CN202110308316.7A priority Critical patent/CN113144997B/en
Publication of CN113144997A publication Critical patent/CN113144997A/en
Application granted granted Critical
Publication of CN113144997B publication Critical patent/CN113144997B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01F—MIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F35/00—Accessories for mixers; Auxiliary operations or auxiliary devices; Parts or details of general application
    • B01F35/20—Measuring; Control or regulation
    • B01F35/21—Measuring
    • B01F35/211—Measuring of the operational parameters
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01F—MIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F31/00—Mixers with shaking, oscillating, or vibrating mechanisms
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01F—MIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F2215/00—Auxiliary or complementary information in relation with mixing
    • B01F2215/04—Technical information in relation with mixing
    • B01F2215/0413—Numerical information
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01F—MIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F2215/00—Auxiliary or complementary information in relation with mixing
    • B01F2215/04—Technical information in relation with mixing
    • B01F2215/0413—Numerical information
    • B01F2215/0436—Operational information
    • B01F2215/0477—Numerical time values
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01F—MIXING, e.g. DISSOLVING, EMULSIFYING OR DISPERSING
    • B01F2215/00—Auxiliary or complementary information in relation with mixing
    • B01F2215/04—Technical information in relation with mixing
    • B01F2215/0413—Numerical information
    • B01F2215/0436—Operational information
    • B01F2215/0481—Numerical speed values

Landscapes

  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

The invention discloses a method and a system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing, wherein the method comprises the following steps: acquiring the resonant sound mixed acceleration of the mixed sample according to the time sequence to obtain an acceleration sample set; dividing samples in the acceleration sample set into a plurality of sample subsets according to a time sequence, wherein each sample subset comprises N continuous acceleration values; calculating the acceleration fluctuation rate of each sample subset, determining the corresponding acceleration fluctuation rate when the mixed samples are uniformly mixed as the critical acceleration fluctuation rate through sampling physical and chemical analysis, acquiring the acceleration fluctuation rate of each sample subset in real time for the samples to be mixed, and stopping mixing when k continuous acceleration fluctuation rates are less than or equal to the critical fluctuation rate, wherein k is more than or equal to 2, and the samples to be mixed reach a uniform state.

Description

Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing
Technical Field
The invention belongs to the technical field of material mixing degree detection, relates to resonance sound mixing uniformity detection, and particularly relates to a method and a system for determining critical fluctuation rate and judging mixing uniformity in resonance sound mixing.
Background
The resonant sound mixing technology is a technology based on vibration macro mixing and sound field micro mixing coupling action which is started in recent years, and relies on low frequency (about 60 Hz) and large acceleration (100g, 9.8m/s for g) of resonant sound mixing equipment2) And vibrating in a reciprocating manner. Under the condition of low-frequency large acceleration vibration, the materials to be mixed are fluidized to generate a macroscopic vibration mixing vortex; meanwhile, the large acceleration vibration excites a sound field (pressure wave) at the bottom of the mixing container, and the sound field generates a couple effect on the materials when propagating inside the materials to form micro-scale mixing (the scale is less than 50 um). The resonant sound mixing has the advantages of whole field mixing, no paddle intervention, high mixing efficiency and the like, is widely applied to aspects of energetic materials, food, pharmacy, chemical eutectic and the like, and can improve the mixing efficiency by 2-10 times compared with the traditional mixing mode.
The same as the problems faced by all mixing modes, the resonance sound mixing also faces the problem of evaluating the mixing uniformity so as to accurately judge the mixing end point and prevent resource waste caused by insufficient mixing or over mixing. The characterization modes of the mixing uniformity comprise direct sampling component test, sampling electron microscope observation, indirect sampling density test, sampling hardness test, sampling chromaticity test and the like. However, there are three main problems brought by sampling test, the first is long period of off-line sampling test, the second is large error of sampling test, and the third is that the local sampling point is used for judging the whole mixing uniformity, which has a certain limitation.
In order to improve the efficiency and accuracy of the mixing uniformity judgment, an online mixing uniformity judgment method is always an object of the efforts of scientific researchers. At the present stage, a better online mixing uniformity judgment method is developed and is a near infrared spectroscopy method, the uniformity is judged by utilizing different degrees of reflection and absorption of near infrared spectra due to different material states in the mixing process, and when the reflection and absorption parameters of the spectra are consistent, the mixing is considered to be uniform, so that the method is more applied in the pharmaceutical industry. However, the application of near infrared spectroscopy to some industrial applications has certain condition limitation, which is mainly expressed in three aspects. Firstly, the near-infrared device is sensitive to environmental influence, the noise, temperature and humidity, vibration and the like of an installation occasion of the near-infrared device are stable, and otherwise, the test accuracy is influenced; secondly, the near-infrared emitting and receiving device must be installed within 2cm of the surface of the measured material, and a sapphire glass window is needed, so that the installation process is complex and the cost is high; third, the uniformity of the near infrared spectroscopy detection is only in response to the surface of the material, since light cannot penetrate the interior of the material, the detection may be affected by the material being stuck to the surface of the sapphire glass window.
The resonance sound mixing is characterized by large acceleration mixing, and the process mainly has three problems in the aspect of uniformity detection: firstly, the mixing speed is high, generally 10-20min, and for materials such as liquid phase, powder and the like, the time for achieving uniform mixing is only minutes or tens of seconds, so that the offline detection is extremely inconvenient for judging the uniformity; secondly, the mixing process container is vibrated with large amplitude and large acceleration, the installation of the near-infrared device is inconvenient, and the test accuracy is difficult to ensure.
Disclosure of Invention
In order to solve the defects in the prior art, the invention provides a method and a system for determining the critical fluctuation rate and judging the mixing uniformity in resonant sound mixing, and solves the defects of long period and large error in the traditional sampling uniformity judgment, harsh installation conditions, obvious environmental influence on the test accuracy and the like in the online method.
In order to solve the technical problems, the invention adopts the following technical scheme:
the method for determining the critical fluctuation rate in the resonant sound mixing comprises the following steps:
step 1, acquiring resonance sound mixed acceleration of a mixed sample according to a time sequence to obtain an acceleration sample set;
step 2, collecting acceleration samplesThe samples in (1) are divided into a plurality of sample subsets according to the time sequence, each sample subset comprises N continuous acceleration values, and the ith sample subset is { G }j,Gj+1,...,Gj+N-1J is 1,2, N +1-N, i is j, N is more than or equal to 2, and N is the sampling number;
step 3, calculating the acceleration fluctuation rate delta of each sample subsetiAcceleration fluctuation rate deltaiCalculated by any one of the formulas (1) to (3),
Figure BDA0002988562550000031
Figure BDA0002988562550000032
Figure BDA0002988562550000033
Gi1for the first acceleration value in the ith sample subset,
Figure BDA0002988562550000035
is the average of the accelerations in the ith subset of samples; gimaxIs the maximum value of acceleration, G, in the ith subset of samplesiminIs the minimum of the accelerations in the ith subset of samples;
Figure BDA0002988562550000034
is the jth acceleration value;
generating an acceleration fluctuation rate curve by the acceleration fluctuation rate values of all the sample subsets;
step 4, determining a corresponding time point when the mixed sample is uniformly mixed through sampling physical and chemical analysis, and taking the time point as a critical time point; according to the acceleration fluctuation rate curve, taking the acceleration fluctuation rate corresponding to the critical time point as the critical acceleration fluctuation rate deltac。
Preferably, the number n of the acceleration samples is more than or equal to 10, and the sampling step length of the acceleration is 1-10 s.
Preferably, in the step 2, N is more than or equal to 3 and less than or equal to 6.
Preferably, in the step 5, k is more than or equal to 5 and less than or equal to 15.
The invention also discloses a method for judging the mixing uniformity of the resonance sound, which comprises the following steps: acquiring acceleration in real time in the mixing process of the materials to be mixed to obtain an acceleration data set; then calculating the acceleration fluctuation rate according to the method from the step 2 to the step 3, and when k continuous acceleration fluctuation rates are less than or equal to the critical fluctuation rate deltacAnd when k is more than or equal to 2, stopping mixing when the mixed sample reaches a uniform state.
The invention also discloses a system for judging the mixing uniformity of the resonance sound, which comprises: an acceleration acquisition module, a data processing module, an acceleration fluctuation rate curve generation module and a mixed state judgment module, wherein,
the acceleration acquisition module is used for acquiring the resonance sound mixed acceleration of the mixed sample in real time to obtain an acceleration data set;
the data processing module comprises a data dividing module and an acceleration fluctuation rate calculating module;
the data dividing module is used for dividing data in the acceleration data set into a plurality of sample subsets according to a time sequence, each sample subset comprises N continuous acceleration values, and the ith sample subset is { G }j,Gj+1,...,Gj+N-1J is 1,2, N +1-N, i is j, N is more than or equal to 2, and N is the sampling number;
the acceleration fluctuation rate calculation module is used for calculating the acceleration fluctuation rate delta of each sample subsetiCalculated by any one of the formulas (1) to (3),
Figure BDA0002988562550000041
Figure BDA0002988562550000042
Figure BDA0002988562550000043
Gi1for the first acceleration value in the ith sample subset,
Figure BDA0002988562550000044
is the average of the accelerations in the ith subset of samples; gimaxIs the maximum value of acceleration, G, in the ith subset of samplesiminIs the minimum of the accelerations in the ith subset of samples;
Figure BDA0002988562550000045
is the jth acceleration value;
the acceleration fluctuation rate curve generation module is used for generating acceleration fluctuation rate curves from the acceleration fluctuation rate values of all the sample subsets;
the mixing state judgment module is used for judging whether the resonant sound mixed material is uniform or not according to the acceleration fluctuation rate and the critical acceleration fluctuation rate, and when the k continuous acceleration fluctuation rates are deltaiLess than or equal to the critical acceleration fluctuation rate deltacWhen k is more than or equal to 2, the resonant sound mixed material reaches a uniform state, and the mixing is stopped.
Specifically, the acceleration acquisition module adopts an acceleration sensor installed on the mixing container to acquire.
Preferably, the number n of the acceleration data of the acceleration acquisition module is more than or equal to 10, and the sampling step length of the acceleration is 1-10 s.
Preferably, in the data dividing module, N is more than or equal to 3 and less than or equal to 6.
Preferably, in the mixed state judgment module, k is greater than or equal to 5 and less than or equal to 15.
Compared with the prior art, the invention has the beneficial effects that:
the method of the invention combines the parameter acceleration of the resonance sound mixing process to judge the uniformity of the mixed material, and has the advantages of online, real-time and accuracy; moreover, no third-party test equipment is introduced, the system is simple and reliable, and is not influenced by environments such as vibration, temperature and humidity; the method of the invention judges the uniformity of the whole mixed flow field, and avoids the influence of wall surfaces, sampling areas and the like on the uniformity judgment of local sampling, near infrared and other methods.
Drawings
Fig. 1 is a graph of the acceleration and its fluctuation rate acquired in example 1.
Fig. 2 is a graph of the acceleration and its fluctuation rate acquired in example 2.
Fig. 3 is a graph of the acceleration and its fluctuation rate acquired in example 3.
Detailed Description
The resonant sound mixing mode has a remarkable characteristic that the vibration acceleration along with the whole mixing process is obviously influenced by the material state besides the material quality, the material type, the input energy and the like. When the mass, the type and the input energy of the materials are fixed, the change of the material state is the only factor influencing the mixing acceleration. In the initial stage of mixing, the state change of the materials is obvious along with the mixing due to the unevenness, and the fluctuation of the mixing acceleration is severe; along with the mixing, the materials are gradually uniform, the state tends to be consistent and stable, and the fluctuation of the mixing acceleration is reduced, so that the vibration acceleration is used as the index for judging the mixing uniformity. On the basis, the resonance sound mixing uniformity judgment method based on the acceleration fluctuation rate is provided, the method firstly needs to determine the critical fluctuation acceleration, and specifically comprises the following steps of 1-4:
step 1, collecting resonance sound mixed acceleration of a mixed sample according to a time sequence, specifically collecting the resonance sound mixed acceleration by adopting an acceleration sensor arranged on a mixing container to obtain an acceleration sample set, wherein the number n of the acceleration samples is preferably more than or equal to 10, and the acceleration sampling step length delta t is preferably 1-10 s.
Step 2, dividing samples in the acceleration sample set into a plurality of sample subsets according to a time sequence, wherein each sample subset comprises N continuous acceleration values, and expressing the ith sample subset as { G }j,Gj+1,...,Gj+N-1J 1,2, N +1-N, i j, N ≧ 2, preferably 3 ≦ N ≦ 6；
Step 3, calculating the acceleration fluctuation rate delta of each sample subsetiThe invention can adopt any one of the following three formulas to calculate,
Figure BDA0002988562550000061
Figure BDA0002988562550000062
Figure BDA0002988562550000063
Gi1the acceleration value is the first acceleration value in the ith sample subset, so that each acceleration sample can only participate in the calculation of the primary fluctuation rate;
Figure BDA0002988562550000064
the average value of the acceleration in the ith sample subset is the sum of all acceleration values in the ith sample subset divided by N; gimaxIs the maximum value of the acceleration in the ith sample subset; giminIs the minimum of the accelerations in the ith subset of samples;
Figure BDA0002988562550000071
is the jth acceleration value and the value is located in the ith sample subset, the acceleration fluctuation rate delta in the formula (3)iIs the standard deviation of the samples of the acceleration in the ith subset of samples.
Generating an acceleration fluctuation rate curve by using the acceleration fluctuation rate values of all the sample subsets;
step 4, determining the critical fluctuation rate delta of the acceleration through sampling analysiscThe method specifically comprises the following steps:
sequentially sampling at different sampling times to perform physical and chemical analysis until a time point of uniform mixing is determined, and taking the time point as a critical time point; according to the curve of acceleration fluctuation rate, the critical point is determinedThe acceleration fluctuation rate corresponding to the time point is taken as the critical acceleration fluctuation rate deltac。
After the critical fluctuation rate of the acceleration is obtained, when the mixing uniformity of other samples is judged, the following method is adopted for judging, specifically:
for a certain sample to be mixed, acquiring acceleration in real time in the mixing process to obtain an acceleration data set; then calculating the acceleration fluctuation rate according to the method from the step 2 to the step 3, and when k continuous acceleration fluctuation rates are less than or equal to the critical fluctuation rate deltacWhen k is more than or equal to 2, preferably, k is more than or equal to 5 and less than or equal to 15, which indicates that the sample to be mixed reaches a uniform state, and the mixing can be stopped.
The sample material to be mixed in the present invention is the same as the sample material to be mixed.
An example of a specific cluster analysis method employed by the present invention is given below. It should be noted that the present invention is not limited to the following embodiments, and all equivalent changes based on the technical solutions of the present application fall into the protection scope of the present invention.
Example 1
The embodiment discloses a process for determining critical fluctuation rate in resonant sound mixing based on acceleration fluctuation rate, which specifically comprises the following steps:
the mixed material in this embodiment is a poured PBX explosive with a solid content of 84%, an acceleration sensor mounted on the mixing equipment is used to collect acceleration, a sampling step length Δ t of the mixing acceleration is 2s, and a total sample number n is 182;
dividing the sample set into a plurality of sample subsets, wherein each sample subset is composed of 3 acceleration samples, namely N is 3, and the sample subset is { G1,G2,G3},{G2,G3,G4},…,{G180,G181,G182}。
Acceleration fluctuation rate delta of each sample subsetiCalculated according to the following formula:
Figure BDA0002988562550000081
at acceleration fluctuation rate deltaiDuring calculation, only the first time fluctuation rate calculation is involved in each acceleration sample in this embodiment, for example, for the first sample subset, the acceleration sample G is included1、G2、G3In the calculation of the acceleration fluctuation rate, only G is used1Making a difference value with the average value; for the second acceleration cell, the sample containing acceleration is G2、G3、G4In the calculation of the acceleration fluctuation rate, only G is used2The difference is made from the average.
The solid line in fig. 1 is an acceleration fluctuation rate curve with N equal to 3. The time point of obtaining uniform mixing is about 280s through sampling test; determining the fluctuation rate corresponding to the point on the fluctuation rate curve to be about 0.02, and setting the critical fluctuation rate deltac＝0.02。
In the same manner as described above, the acceleration fluctuation rate is calculated once more with N-6, and the dashed-dotted line in fig. 1 is an acceleration fluctuation rate curve with N-6. Likewise, the corresponding critical fluctuation ratio is 0.02.
The test material of this example is also a cast PBX explosive with 84% solids; in the mixing process, acquiring acceleration, dividing samples and calculating the fluctuation rate of the acceleration in real time according to the above thought; when the 5 consecutive acceleration fluctuation rates (i.e., k-5) are 0.02 or less, the mixing is stopped, and the corresponding time in fig. 1 is about 286 s. It can also be seen from fig. 1 that when the mixing time reaches around 280s, the acceleration fluctuation becomes significantly smaller, tending to a steady value.
Example 2
The mixed material in this embodiment is a poured PBX explosive with a solid content of 84%, a sampling step length Δ t of a mixing acceleration in a resonant acoustic mixing process is 2s, and a total number of samples n is 182; each acceleration unit consists of 3 acceleration samples, namely N is 3; acceleration fluctuation rate delta of each sample subsetiCalculated according to the following formula:
Figure BDA0002988562550000091
in the same manner, the acceleration fluctuation rate was calculated again in accordance with N ═ 6.
Fig. 2 shows the acceleration and the acceleration fluctuation rate thereof in the present embodiment, in which the solid line shows the acceleration fluctuation rate with N-3, and the dashed dotted line shows the acceleration fluctuation rate with N-6.
Through sampling test, the time point of uniform mixing is about 280s, the fluctuation rate corresponding to the point is determined to be about 0.05 on the fluctuation rate curve, and the critical fluctuation rate delta is setc＝0.05。
In this embodiment, if the acceleration fluctuation rate is continuously 5 times (i.e., k is 5) and is 0.05 or less, the mixture is uniformly mixed, and the time period in 2 is about 290 s.
Example 3
The mixed material in this embodiment is a poured PBX explosive with a solid content of 84%, a sampling step length Δ t of a mixing acceleration in a resonant acoustic mixing process is 2s, and a total number of samples n is 182; each acceleration unit consists of 3 acceleration samples, namely N is 3; fluctuation rate δ of acceleration of each sample subsetiIs the standard deviation of the samples for that subset of samples,
Figure BDA0002988562550000092
at acceleration fluctuation rate deltaiDuring calculation, only the first time fluctuation rate calculation is involved in each acceleration sample in this embodiment, for example, for the first sample subset, the acceleration sample G is included1、G2、G3In the calculation of the acceleration fluctuation rate, only G is used1Making a difference value with the average value; for the second acceleration cell, the sample containing acceleration is G2、G3、G4In the calculation of the acceleration fluctuation rate, only G is used2The difference is made from the average.
In the same manner, the acceleration fluctuation rate was calculated again with N ═ 6.
Fig. 3 shows the acceleration and the acceleration fluctuation rate thereof in the present embodiment, in which the solid line shows the acceleration fluctuation rate with N-3, and the dashed dotted line shows the acceleration fluctuation rate with N-6.
Through sampling test, the time point of uniform mixing is about 280s, the fluctuation rate corresponding to the time point is determined to be about 1.4 on the fluctuation rate curve, and the critical fluctuation rate delta is setc＝1.4。
In this embodiment, if the continuous k is 5 waves with a fluctuation ratio of 1.4 or less, the mixture is uniform, and the corresponding time in fig. 3 is about 290 s.
Example 4
The embodiment discloses a resonance sound mixing uniformity judgment system based on acceleration fluctuation rate, which comprises: an acceleration acquisition module, a data processing module, an acceleration fluctuation rate curve generation module and a mixed state judgment module, wherein,
the acceleration acquisition module is used for acquiring the resonant sound mixed acceleration of the mixed sample according to the time sequence to obtain an acceleration data set; specifically, the sampling is performed using an acceleration sensor mounted on the mixing container. The acceleration sampling step length of the embodiment is 1-10 s.
The data processing module comprises a data dividing module and an acceleration fluctuation rate calculating module; specifically, the data dividing module is configured to divide the samples in the acceleration sample set into a plurality of sample subsets according to a time sequence, each sample subset includes N consecutive acceleration values, and an ith sample subset is { G }j,Gj+1,…,Gj+N-1J is 1,2, …, N +1-N, i is j, in this embodiment, N is equal to or greater than 3 and equal to or less than 6, and N is the number of samples;
the acceleration fluctuation rate calculation module is used for calculating the acceleration fluctuation rate delta of each sample subsetiThe calculation formula adopts any one of formulas (1) to (3);
the acceleration fluctuation rate curve generation module is used for generating an acceleration fluctuation rate curve from the acceleration fluctuation rate values of all the sample subsets and displaying the curve;
the mixing state judgment module is used for judging whether the resonant sound mixed material is uniform or not according to the acceleration fluctuation rate and the critical acceleration fluctuation rate, and particularly when k acceleration fluctuation rates delta are continuousiLess than or equal to the critical acceleration fluctuation rate deltacWhen the material is mixed, the resonance sound is uniform; k is more than or equal to 2, preferably, k is more than or equal to 5 and less than or equal to 15. And then according to the judgment instruction output by the module, ending the resonance mixing process of the mixing container.

Claims (10)

1. The method for determining the critical fluctuation rate in the resonant sound mixing is characterized by comprising the following steps of:
step 1, acquiring resonance sound mixed acceleration of a mixed sample according to a time sequence to obtain an acceleration sample set;
step 2, dividing samples in the acceleration sample set into a plurality of sample subsets according to a time sequence, wherein each sample subset comprises N continuous acceleration values, and the ith sample subset is { G }j,Gj+1,...,Gj+N-1J is 1,2, N +1-N, i is j, N is more than or equal to 2, and N is the sampling number;
step 3, calculating the acceleration fluctuation rate delta of each sample subsetiAcceleration fluctuation rate deltaiCalculated by any one of the formulas (1) to (3),
Figure FDA0002988562540000011
Figure FDA0002988562540000012
Figure FDA0002988562540000013
Gi1for the first acceleration value in the ith sample subset,
Figure FDA0002988562540000014
is the average of the accelerations in the ith subset of samples; gimaxIs the maximum value of acceleration, G, in the ith subset of samplesiminIs the minimum of the accelerations in the ith subset of samples;
Figure FDA0002988562540000015
is the jth acceleration value;
generating an acceleration fluctuation rate curve by the acceleration fluctuation rate values of all the sample subsets;
step 4, determining a corresponding time point when the mixed sample is uniformly mixed through sampling physical and chemical analysis, and taking the time point as a critical time point; according to the acceleration fluctuation rate curve, taking the acceleration fluctuation rate corresponding to the critical time point as the critical acceleration fluctuation rate deltac。
2. The method for determining the critical fluctuation rate in resonant acoustic mixing according to claim 1, wherein in the step 1, the number n of acceleration samples is greater than or equal to 10, and the sampling step length of the acceleration is 1-10 s.
3. The method for determining critical fluctuation ratio in resonant acoustic mixing according to claim 1, wherein in step 2, N is 3. ltoreq. N.ltoreq.6.
4. The method for determining critical fluctuation ratio in resonant acoustic mixing according to claim 1, wherein in step 5, k is 5. ltoreq. k.ltoreq.15.
5. The method for determining the mixing uniformity of the resonant sound is characterized by comprising the following steps: acquiring acceleration in real time in the mixing process of the materials to be mixed to obtain an acceleration data set; then calculating the acceleration fluctuation rate according to the method of step 2 to step 3 in claim 1, when k successive acceleration fluctuation rates less than or equal to the critical fluctuation rate δ occurcAnd when k is more than or equal to 2, stopping mixing when the mixed sample reaches a uniform state.
6. A system for determining the uniformity of a mixture of resonant sounds, comprising: an acceleration acquisition module, a data processing module, an acceleration fluctuation rate curve generation module and a mixed state judgment module, wherein,
the acceleration acquisition module is used for acquiring the resonance sound mixed acceleration of the mixed sample in real time to obtain an acceleration data set;
the data processing module comprises a data dividing module and an acceleration fluctuation rate calculating module;
the data dividing module is used for dividing data in the acceleration data set into a plurality of sample subsets according to a time sequence, each sample subset comprises N continuous acceleration values, and the ith sample subset is { G }j,Gj+1,...,Gj+N-1J is 1,2, N +1-N, i is j, N is more than or equal to 2, and N is the sampling number;
the acceleration fluctuation rate calculation module is used for calculating the acceleration fluctuation rate delta of each sample subsetiCalculated by any one of the formulas (1) to (3),
Figure FDA0002988562540000021
Figure FDA0002988562540000022
Figure FDA0002988562540000023
Gi1is the first acceleration value, G, in the ith sample subsetiIs the average of the accelerations in the ith subset of samples; gimaxIs the maximum value of acceleration, G, in the ith subset of samplesiminIs the minimum of the accelerations in the ith subset of samples;
Figure FDA0002988562540000031
is the jth acceleration value;
the acceleration fluctuation rate curve generation module is used for generating acceleration fluctuation rate curves from the acceleration fluctuation rate values of all the sample subsets;
the mixed state judging module is used for judging the critical acceleration according to the acceleration fluctuation rateJudging whether the resonant sound mixed material is uniform or not according to the degree fluctuation rate, and judging whether the resonant sound mixed material is uniform or not when k acceleration fluctuation rates delta are continuousiLess than or equal to the critical acceleration fluctuation rate deltacWhen k is more than or equal to 2, the resonant sound mixed material reaches a uniform state, and the mixing is stopped.
7. The system for determining the uniformity of a resonant acoustic mixing according to claim 6, wherein the acceleration acquisition module acquires the signals using an acceleration sensor mounted on the mixing container.
8. The system for determining the mixing uniformity of the resonant sound according to claim 6, wherein the number n of the acceleration data of the acceleration acquisition module is greater than or equal to 10, and the sampling step length of the acceleration is 1-10 s.
9. The system for determining the uniformity of a resonant acoustic mixing according to claim 6, wherein N is 3. ltoreq. N.ltoreq.6 in the data partitioning module.
10. The system for determining the uniformity of mixing of a resonating sound according to claim 6, wherein k is 5 ≦ k ≦ 15 in the mixing status determining module.
CN202110308316.7A 2021-03-23 2021-03-23 Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing Active CN113144997B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110308316.7A CN113144997B (en) 2021-03-23 2021-03-23 Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110308316.7A CN113144997B (en) 2021-03-23 2021-03-23 Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing

Publications (2)

Publication Number Publication Date
CN113144997A true CN113144997A (en) 2021-07-23
CN113144997B CN113144997B (en) 2022-10-25

Family

ID=76888193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110308316.7A Active CN113144997B (en) 2021-03-23 2021-03-23 Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing

Country Status (1)

Country Link
CN (1) CN113144997B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113731279A (en) * 2021-08-31 2021-12-03 华中科技大学 Online evaluation method and equipment for mixing state in acoustic resonance mixing process

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100294113A1 (en) * 2007-10-30 2010-11-25 Mcpherson Michael D Propellant and Explosives Production Method by Use of Resonant Acoustic Mix Process
US20150059746A1 (en) * 2012-03-30 2015-03-05 Vectura Limited Method and apparatus
CN107051296A (en) * 2017-01-23 2017-08-18 西安近代化学研究所 A kind of electromagnetic excitation resonance mixing arrangement and its control method
CN108287465A (en) * 2018-01-03 2018-07-17 湖北航鹏化学动力科技有限责任公司 Control system and method for acoustic resonance material mixer
CN109694292A (en) * 2018-12-25 2019-04-30 西安近代化学研究所 A kind of cast PBX explosive resonance hybrid technique
CN109939600A (en) * 2019-03-12 2019-06-28 湖北航鹏化学动力科技有限责任公司 Resonance mixing method and resonance mixing equipment for fire transfer agent
CN111024214A (en) * 2019-12-25 2020-04-17 华中科技大学 Method for acquiring natural frequency of acoustic resonance mixer in real time in operation process

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100294113A1 (en) * 2007-10-30 2010-11-25 Mcpherson Michael D Propellant and Explosives Production Method by Use of Resonant Acoustic Mix Process
US20150059746A1 (en) * 2012-03-30 2015-03-05 Vectura Limited Method and apparatus
CN107051296A (en) * 2017-01-23 2017-08-18 西安近代化学研究所 A kind of electromagnetic excitation resonance mixing arrangement and its control method
CN108287465A (en) * 2018-01-03 2018-07-17 湖北航鹏化学动力科技有限责任公司 Control system and method for acoustic resonance material mixer
CN109694292A (en) * 2018-12-25 2019-04-30 西安近代化学研究所 A kind of cast PBX explosive resonance hybrid technique
CN109939600A (en) * 2019-03-12 2019-06-28 湖北航鹏化学动力科技有限责任公司 Resonance mixing method and resonance mixing equipment for fire transfer agent
CN111024214A (en) * 2019-12-25 2020-04-17 华中科技大学 Method for acquiring natural frequency of acoustic resonance mixer in real time in operation process

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张毅铭等: "固液两相声共振混合数值模拟", 《化工进展》 *
陈松等: "声共振混合设备用于百克量级火炸药制备的能量转换特性", 《火炸药学报》 *
马宁等: "塑料粘结炸药模拟物共振声混合工艺过程变化规律研究", 《兵工学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113731279A (en) * 2021-08-31 2021-12-03 华中科技大学 Online evaluation method and equipment for mixing state in acoustic resonance mixing process

Also Published As

Publication number Publication date
CN113144997B (en) 2022-10-25

Similar Documents

Publication Publication Date Title
JP3151629B2 (en) Method of measuring distribution of component concentration in centrifugal separation system by transmitting / receiving mechanical signal
Comte et al. Microprobe-type measurement of Young's modulus and Poisson coefficient by means of depth sensing indentation and acoustic microscopy
Schork et al. On‐line measurement of surface tension and density with applications to emulsion polymerization
Williams et al. Multi-phase flow measurements in powder processing
CN113144997B (en) Method and system for determining critical fluctuation rate and judging mixing uniformity in resonant sound mixing
Best et al. Modelling harmonic generation measurements in solids
Toledo et al. Piezoelectric MEMS resonators for monitoring grape must fermentation
CN109443406B (en) Acoustic fluid sensor
Fathi-Haftshejani et al. Nondestructive evaluation of clad rods by inversion of acoustic scattering data
Álvarez-Arenas et al. Characterization of suspensions of particles in water by an ultrasonic resonant cell
JP4121460B2 (en) Acoustic cell for material analysis
US8307711B2 (en) Apparatus for inspection of a fluid and method
US7331232B2 (en) Measurement method and biosensor apparatus using resonator
Campos-Pozuelo et al. Comparative study of the nonlinear behavior of fatigued and intact samples of metallic alloys
Jang et al. Techniques for cure monitoring of thermoset resins and composites-a review
JP2006052996A (en) Determination method of extreme value frequency
Addali Monitoring gas void fraction in two-phase flow with acoustic emission
US5376276A (en) In situ primary froth quality measurements using microwave monitor
Parker et al. Ultrasonic Sensor for the characterization of colloidal slurries
RU2217742C2 (en) Method of prediction of residual life of metal of mining equipment
RU2485482C2 (en) Method to determine porosity of ceramic and silicate materials using magnetic fluid
Villa-Arango et al. Acoustic spectrometer: Resonant sensing platform for measuring volumetric properties of liquid samples
Tsahalis Modelling of fluid damping in thermal insulation
RU2431135C1 (en) Method of ultrasound tube surface roughness measurement (versions)
Liang et al. Coupling of Induction with Damping Behavior for Viscosity Sensing via Design of Magnetized Oscillator

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant