CN105322891B - A kind of multiaxis sinusoidal vibration controls drive signal generation method - Google Patents
A kind of multiaxis sinusoidal vibration controls drive signal generation method Download PDFInfo
- Publication number
- CN105322891B CN105322891B CN201510280799.9A CN201510280799A CN105322891B CN 105322891 B CN105322891 B CN 105322891B CN 201510280799 A CN201510280799 A CN 201510280799A CN 105322891 B CN105322891 B CN 105322891B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- frequency
- drive signal
- phase
- 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.)
- Expired - Fee Related
Links
Landscapes
- Apparatuses For Generation Of Mechanical Vibrations (AREA)
- Particle Accelerators (AREA)
Abstract
The invention discloses a kind of multiaxis sinusoidal vibration to control drive signal generation method.The drive signal generation method carries out data retrieval according to calculation procedure according to the frequency of sinusoidal signal to be generated, amplitude, phase and four parameters of biasing in sinusoidal data list, and the data generated after retrieval are sent with fixed frequency by DA data collecting plate cards.The sinusoidal drive signals of drive signal generation method generation are highly stable, and signal can reach continuous, smooth, smoothly transition when frequency, amplitude or phase arbitrary parameter change.
Description
Technical field
The invention belongs to spacecraft dynamic environmental test field, and in particular to a kind of multiaxis sinusoidal vibration controls drive signal
Generation method.
Background technology
In multiaxis sinusoidal vibration frequency sweep test, control system requirement is in real time to the amplitude, frequency and phase of drive signal
Position is modified, so that response meets the requirement of reference spectrum.Therefore, continuous sinusoidal drive signals how are produced, and ensure letter
It is multiaxis sinusoidal vibration control software exploitation asking of must studying and solve number to reach smooth smoothly transition when parameter updates
One of topic.Multiaxis sinusoidal vibration control drive signal generation method proposed by the present invention based on data list retrieval type, according to
Frequency, amplitude, phase and four parameters of biasing of sinusoidal signal to be generated enter line number according to calculation procedure in sinusoidal data list
According to retrieval, the data generated after retrieval are sent into board by DA data with fixed frequency and sent.This method generation
Sinusoidal drive signals are highly stable, and signal can reach when the parameters such as frequency, amplitude or phase change it is continuous, smooth,
Smoothly transition.
The content of the invention
The technical problem to be solved in the present invention is that providing a kind of multiaxis sinusoidal vibration controls drive signal generation method, with
Make the sinusoidal drive signals of generation highly stable, and signal can reach company when the parameters such as frequency, amplitude or phase change
Continuous, smooth, smoothly transition.
In order to solve the above-mentioned technical problem, the invention provides a kind of multiaxis sinusoidal vibration control drive signal generation side
Method, comprise the following steps:
Four key elements of the generation of drive signal based on sinusoidal signal:Frequency, amplitude, phase and biasing,
If current, wherein the frequency of drive signal is f ' all the wayi, amplitude is A 'i, phase is α 'i, it is biased to ζ 'i, D/A adopts
Sample frequency is Fs, for the feedback by a closed loop, it is necessary to generate new drive signal, its parameter is as follows:Frequency fi, amplitude Ai, phase
Position is αi, bias ζi(i=1,2 ...), it is in particular:
(a) it is N first by the points of computer generation one, amplitude is 1 sinusoidal fundamental wave list complete cycle:
Wherein, N is points sine wave complete cycle, value 0x80000, pt=0,1,2......, N-1.
(b) the frequency f ' of certain current drive signal all the way is obtainedi, amplitude A 'i, phaseWith biasing ζ 'iFour basic parameters
The frequency f of value and new drive signal to be generatedi, amplitude Ai, phaseWith biasing ζiFour base values, wherein phaseWith points
To represent:
Wherein, αiFor the true phase of signal, unit is radian;
(c) buffer time of driving signal frequency, amplitude and phase takes the 1/4 of the new signal cycle to be generated, biasing
Buffer time takes a steady state value, i.e.,
(d) several frequency F are sent according to DAsDetermine signal frequency, amplitude, phase and the buffer point of biasing respectively with buffer time
Number Nf(i), NA(i),Nζ(i), its value is calculated with following equation respectively:
Nf(i)=Tf(i)·Fs, NA(i)=TA(i)·Fs…………………………(5)
(e) incremental change that drive signal parameters are buffered to signal parameter value to be generated by currency is calculated, uses symbol respectively
Δf(i), ΔA(i),Δζ(i) represent, calculation formula is as follows:
(f) each frequency f ' for waiting to send several points of a new frame data is determinedi(j), amplitude A 'i(j), phaseBias ζ 'i
(j) (j=0,1,2 ...), calculation formula is as follows:
f′i(j)=f 'i+(j+1)·Δf(i), A 'i(j)=A 'i+(j+1)·ΔA(i)……··…··…(9)
Wherein i is that the way of drive signal is numbered;J is the discrete digital signal point numbering of the i roads drive signal of generation.
(g) each index in sinusoidal fundamental wave for sending several points is calculated, formula is as follows:
(h) each actual value for waiting to send several points is finally calculated:
xi(j)=A 'i(j)Sine[Sj]+ζ′i(j)…………………………··(12)
The discrete series for the drive signal that formula (12) is calculated is converted into the digital quantity under DA reference voltages respectively, so
Cyclically sequence forfeited in successively afterwards it is several in a manner of be stored in DA cachings (FIFO), DA boards are to send several frequency FsBy each point according to
The smoothing processing of device can produce continuous sinusoidal drive signals after secondary submitting and after filtering.
Multiaxis sinusoidal vibration provided by the invention controls drive signal generation method, according to the frequency of sinusoidal signal to be generated
Rate, amplitude, phase and four parameters of biasing carry out data retrieval according to calculation procedure in sinusoidal data list, by the number after retrieval
According to being handled and sent with fixed frequency by DA data collecting plate cards.The sinusoidal drive signals of this method generation are non-
It is often stable, and signal can reach continuous, smooth, smoothly transition when the parameters such as frequency, amplitude or phase change.
Brief description of the drawings
The drive signal of the tunnel difference amplitudes of Tu1Wei tri-.
The drive signal of the tunnel outs of phase of Tu2Wei tri-.
The drive signal that the tunnel frequencies of Tu3Wei tri-, amplitude and phase change.
Embodiment
A kind of multiaxis sinusoidal vibration control drive signal generation method to the present invention is further described below.
(1) drive signal generation method
Multiaxis sinusoidal vibration control needs to generate multichannel independent drive signals.The generation of drive signal is based on sinusoidal signal
Four key elements:Frequency, amplitude, phase and biasing.
If current, wherein the frequency of drive signal is f ' all the wayi, amplitude is A 'i, phase is α 'i, it is biased to ζ 'i, D/A adopts
Sample frequency is Fs, for the feedback by a closed loop, it is necessary to generate new drive signal, its parameter is as follows:Frequency fi, amplitude Ai, phase
Position is αi, bias ζi(i=1,2 ...).It is transformed into new drive signal by current drive signal, there is one section of smooth transition centre
Journey.Wait after smooth transition process terminates, DA (digital-to-analogue conversion) ends are just with new drive signal parameters stably output signal.Its is detailed
Thin formation mechanism is as follows:
(a) it is N first by the points of computer generation one, amplitude is 1 sinusoidal fundamental wave list complete cycle:
Wherein N is points sine wave complete cycle, value 0x80000, pt=0,1,2......, N-1.
(b) the frequency f ' of certain current drive signal all the way is obtainedi, amplitude A 'i, phaseWith biasing ζ 'iFour basic parameters
The frequency f of value and two-way drive signal to be generatedi, amplitude Ai, phaseWith biasing ζiFour base values, wherein phaseWith point
Count to represent (αiFor the true phase of signal, unit is radian):
(c) change of any driving signal frequency all the way, four amplitude, phase and biasing basic parameters is not to be mutated
Interior generation, but one section of buffering course is needed, to reach the smooth transition of parameter.This section of buffering course is with parameter buffer time
To represent.The buffer time of driving signal frequency, amplitude and phase takes the 1/4 of new signal cycle to be generated, the buffering of biasing
Time takes a steady state value.I.e.
(d) several frequency F are sent according to DAsDetermine signal frequency, amplitude, phase and the buffer point of biasing respectively with buffer time
Number Nf(i), NA(i),Nζ(i).Its value is calculated with following equation respectively:
Nf(i)=Tf(i)·Fs, NA(i)=TA(i)·Fs…………·………………(5)
(e) incremental change that drive signal parameters are buffered to signal parameter value to be generated by currency is calculated, uses symbol respectively
Δf(i), ΔA(i),Δζ(i) represent.Calculation formula is as follows:
(f) each frequency f ' for waiting to send several points of a new frame data is determinedi(j), amplitude A 'i(j), phase(j) ζ ', is biasedi
(j) (j=0,1,2 ...).Calculation formula is as follows:
f′i(j)=f 'i+(j+1)·Δf(i), A 'i(j)=A 'i+(j+1)·ΔA(i)…··……··…(9)
Wherein i is that the way of drive signal is numbered;J is the discrete digital signal point numbering of the i roads drive signal of generation.
(g) each index in sinusoidal fundamental wave for sending several points is calculated.Formula is as follows:
It can be seen from formula (7)~(11) when two parameters of frequency and phase of signal change, buffer time
Interior sampling rate is also what is changed over time.Only as each frequency f ' for waiting to send several pointsiAnd phase (j)Reach fiWithWhen, sampling rate just reaches a steady state value.Now the index value of each sample point and time are linear.
(h) each actual value for waiting to send several points is finally calculated:
xi(j)=A 'i(j)Sine[Sj]+ζ′i(j)…………………………(12)
The digital quantity discrete series of the above-mentioned drive signal calculated being converted into respectively under DA reference voltages, Ran Houyi
It is secondary cyclically sequence is forfeited in it is several in a manner of be stored in DA cachings (FIFO), DA plates are to send several frequency FsEach point is sent out successively
Afterwards and the smoothing processing of device can produce continuous sinusoidal drive signals after filtering.
(2) drive signal generation emulation
According to sinusoidal drive signals formation mechanism, drive signal generation program function module has been write, and in Matlab rings
The emulation of three road drive signals has been carried out under border.Different drivings can be generated when driving signal frequency, amplitude, phase parameter regulation
Signal, general main point of following three kinds of situations:
(a) drive signal of three tunnel difference amplitudes
The tunnel difference amplitude drive signal parameters table of table 1 three
Frequency (Hz) | Amplitude (V) | Phase (°) | |
Drive1 | 5.0 | 2.0 | 0 |
Drive2 | 5.0 | 4.0 | 0 |
Drive3 | 5.0 | 6.0 | 0 |
(b) drive signal of three tunnel outs of phase
The tunnel out of phase drive signal parameters table of table 2 three
Frequency (Hz) | Amplitude (V) | Phase (°) | |
Drive1 | 5.0 | 4.0 | 0 |
Drive2 | 5.0 | 4.0 | 90 |
Drive3 | 5.0 | 4.0 | 180 |
(c) drive signal that three tunnel frequencies, amplitude and phase change
The drive signal parameters table that the tunnel frequency of table 3 three, amplitude and phase change
Frequency changes (Hz) | Amplitude changes (V) | Phase place change (°) | |
Drive1 | 5.0→5.5 | 2.0→2.5 | 0→0 |
Drive2 | 5.0→5.5 | 4.0→4.5 | 0→45 |
Drive3 | 5.0→5.5 | 6.0→5.5 | 0→90 |
1~Fig. 3 of accompanying drawing sets forth simulation result of the sinusoidal signal under different frequency, amplitude and phase condition.From
Simulation result can be seen that generation driving curve continuous and derivable at Parameters variation, reach smooth transition.
Although giving detailed description and explanation to the embodiment of the present invention above, it should be noted that
We can carry out various equivalent changes and modification to above-mentioned embodiment according to the conception of the present invention, and its caused function is made
, all should be within protection scope of the present invention during with the spirit still covered without departing from specification and accompanying drawing.
Claims (2)
1. a kind of multiaxis sinusoidal vibration controls drive signal generation method, it is characterised in that the multiaxis sinusoidal vibration control is driven
Dynamic signal creating method includes:
Four key elements of the generation of drive signal based on sinusoidal signal:Frequency, amplitude, phase and biasing,
If current, wherein the frequency of drive signal is f ' all the wayi, amplitude is A 'i, phase is α 'i, it is biased to ζ 'i, D/A sample frequencys
For Fs, for the feedback by a closed loop, it is necessary to generate new drive signal, its parameter is as follows:Frequency fi, amplitude Ai, phase be
αi, bias ζi(i=1,2 ...), it is in particular:
(a) it is N first by the points of computer generation one, amplitude is 1 sinusoidal fundamental wave list complete cycle:
<mrow>
<mi>S</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mo>&lsqb;</mo>
<mi>p</mi>
<mi>t</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mi>p</mi>
<mi>t</mi>
</mrow>
<mi>N</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>...</mn>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, N is points sine wave complete cycle, value 0x80000, pt=0,1,2......, N-1;
(b) the frequency f ' of certain current drive signal all the way is obtainedi, amplitude A 'i, phaseWith biasing ζ 'iFour basic parameter values and
The frequency f of two-way drive signal to be generatedi, amplitude Ai, phaseWith biasing ζiFour base values, wherein phaseWith points come table
Show:
Wherein, αiFor the true phase of signal, unit is radian;
(c) buffer time of driving signal frequency, amplitude and phase takes the 1/4 of the new signal cycle to be generated, the buffering of biasing
Time takes a steady state value, i.e.,
<mrow>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>T</mi>
<mi>A</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>4</mn>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mn>...</mn>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
(d) several frequency F are sent according to DAsDetermine signal frequency, amplitude, phase and the buffering of biasing points N respectively with buffer timef
(i), NA(i),Nζ(i), its value is calculated with following equation respectively:
Nf(i)=Tf(i)·Fs, NA(i)=TA(i)·Fs…………·………………(5)
(e) incremental change that drive signal parameters are buffered to signal parameter value to be generated by currency is calculated, respectively with symbol Δf
(i), ΔA(i),Δζ(i) represent, calculation formula is as follows:
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>f</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>&Delta;</mi>
<mi>A</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>A</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>A</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mn>...</mn>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
(f) each frequency f ' for waiting to send several points of a new frame data is determinedi(j), amplitude A 'i(j), phaseBias ζ 'i(j)(j
=0,1,2 ...), calculation formula is as follows:
f′i(j)=f 'i+(j+1)·Δf(i), A 'i(j)=A 'i+(j+1)·ΔA(i)…………·…(9)
ζ′i(j)=ζ 'i+(j+1)·Δζ(i)……………(10)
Wherein i is that the way of drive signal is numbered;J is the discrete digital signal point numbering of the i roads drive signal of generation;
(g) each index in sinusoidal fundamental wave for sending several points is calculated, formula is as follows:
(h) each actual value for waiting to send several points is finally calculated:
xi(j)=A 'i(j)Sine[Sj]+ζ′i(j)……………………………(12)
The discrete series for the drive signal that formula (12) is calculated is converted into the digital quantity under DA reference voltages, Ran Houyi respectively
It is secondary cyclically sequence is forfeited in it is several in a manner of be stored in DA cachings, DA plates are to send several frequency FsAfter each point is sent out successively and pass through
The smoothing processing of wave filter can produce continuous sinusoidal drive signals.
2. multiaxis sinusoidal vibration as claimed in claim 1 controls drive signal generation method, it is characterised in that in the frequency of signal
When two parameters of rate and phase change, sampling rate is also what is changed over time in buffer time, waits to send several points when each
Frequency f 'iAnd phase (j)Reach fiWithWhen, sampling rate just reaches a steady state value, now the index value of each sample point
It is linear with the time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510280799.9A CN105322891B (en) | 2015-05-28 | 2015-05-28 | A kind of multiaxis sinusoidal vibration controls drive signal generation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510280799.9A CN105322891B (en) | 2015-05-28 | 2015-05-28 | A kind of multiaxis sinusoidal vibration controls drive signal generation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105322891A CN105322891A (en) | 2016-02-10 |
CN105322891B true CN105322891B (en) | 2018-04-06 |
Family
ID=55249629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510280799.9A Expired - Fee Related CN105322891B (en) | 2015-05-28 | 2015-05-28 | A kind of multiaxis sinusoidal vibration controls drive signal generation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105322891B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111256927B (en) * | 2019-12-31 | 2022-03-25 | 杭州亿恒科技有限公司 | Self-adaptive double-closed-loop time adjustment sinusoidal vibration control method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303610A (en) * | 2008-07-04 | 2008-11-12 | 西北工业大学 | Signal discrete method for signal generating device sampling step length complicated change |
CN101741318A (en) * | 2009-12-31 | 2010-06-16 | 中国航空工业集团公司北京长城计量测试技术研究所 | Method for controlling degree of distortion of sine vibrating table |
CN102645585A (en) * | 2012-05-14 | 2012-08-22 | 中国航空工业集团公司北京长城计量测试技术研究所 | Synchronous rapid measuring method and device for ultralow frequency sinusoidal signal phase difference |
CN105092197A (en) * | 2015-06-19 | 2015-11-25 | 北京航天斯达科技有限公司 | Multi-degree-of-freedom sine vibration control method and controller |
CN106130324A (en) * | 2016-06-29 | 2016-11-16 | 韩伟 | A kind of SPWM waveform controlling method based on FPGA |
CN106502308A (en) * | 2016-09-20 | 2017-03-15 | 江苏大学 | A kind of wave generator system of ultrasonic pulse pumping signal and production method |
-
2015
- 2015-05-28 CN CN201510280799.9A patent/CN105322891B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101303610A (en) * | 2008-07-04 | 2008-11-12 | 西北工业大学 | Signal discrete method for signal generating device sampling step length complicated change |
CN101741318A (en) * | 2009-12-31 | 2010-06-16 | 中国航空工业集团公司北京长城计量测试技术研究所 | Method for controlling degree of distortion of sine vibrating table |
CN102645585A (en) * | 2012-05-14 | 2012-08-22 | 中国航空工业集团公司北京长城计量测试技术研究所 | Synchronous rapid measuring method and device for ultralow frequency sinusoidal signal phase difference |
CN105092197A (en) * | 2015-06-19 | 2015-11-25 | 北京航天斯达科技有限公司 | Multi-degree-of-freedom sine vibration control method and controller |
CN106130324A (en) * | 2016-06-29 | 2016-11-16 | 韩伟 | A kind of SPWM waveform controlling method based on FPGA |
CN106502308A (en) * | 2016-09-20 | 2017-03-15 | 江苏大学 | A kind of wave generator system of ultrasonic pulse pumping signal and production method |
Also Published As
Publication number | Publication date |
---|---|
CN105322891A (en) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104537260B (en) | The dynamic soft-measuring method and system returned based on slow feature | |
CN103900610B (en) | MEMS gyro random error Forecasting Methodology based on Lycoperdon polymorphum Vitt wavelet neural network | |
CN102968995B (en) | A kind of sound mixing method and device of audio signal | |
Lessen | On the stability of the free laminar boundary layer between parallel streams | |
CN106546400B (en) | A kind of multiple-input and multiple-output non-gaussian random vibration test system and test method | |
Zhao et al. | Discussions on observer design of nonlinear positive systems via T–S fuzzy modeling | |
CN102915036A (en) | Method for suppressing limit cycle of inclination angle control system of aircraft with parameter uncertainty | |
CN107832870A (en) | Project of transmitting and converting electricity cost forecasting method based on grey correlation analysis and SVMs | |
CN102110010A (en) | Hardware-in-the-loop (HIL) real-time simulation platform of permanent magnet linear synchronous motor | |
CN104750982A (en) | Reliability interval estimation method based on proportion grouping re-sampling | |
CN105322891B (en) | A kind of multiaxis sinusoidal vibration controls drive signal generation method | |
Rizwan et al. | Implementation of ANN-based embedded hybrid power filter using HIL-topology with real-time data visualization through node-RED | |
CN109782588A (en) | The tight format non-model control method of the different factor of MIMO | |
CN104240263B (en) | A kind of motion subtree system for Delta parallel manipulators | |
Resch et al. | Stability analysis of power hardware-in-the-loop simulations for grid applications | |
CN103279030A (en) | Bayesian framework-based dynamic soft measurement modeling method and device | |
Hirase et al. | Dynamics analysis using koopman mode decomposition of a microgrid including virtual synchronous generator-based inverters | |
Kokovin et al. | Intelligent power electronic converter for wired and wireless distributed applications | |
Zhang et al. | Vibration control of disturbed all-clamped plate with an inertial actuator based on cascade active disturbance rejection control | |
Kengne et al. | Optimal phase control in a Remoissenet–Peyrard substrate potential: numerical and analogical investigations | |
Andrs | Using Industry 4.0 technologies for teaching and learning in education process | |
Ribeiro et al. | Brownian regime of finite-N corrections to particle motion in the XY Hamiltonian mean field model | |
CN201887775U (en) | Simulation system for radio altitude echoed signals | |
CN103647635B (en) | Step wave based multiple scroll wave circuit | |
CN103150475B (en) | The flexible measurement method of melt index |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180406 Termination date: 20210528 |