CN107469702A - A kind of nonlinear chaotic vibration rod - Google Patents
A kind of nonlinear chaotic vibration rod Download PDFInfo
- Publication number
- CN107469702A CN107469702A CN201710657775.XA CN201710657775A CN107469702A CN 107469702 A CN107469702 A CN 107469702A CN 201710657775 A CN201710657775 A CN 201710657775A CN 107469702 A CN107469702 A CN 107469702A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- alpha
- rotating shaft
- chaos
- 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
Links
Classifications
-
- 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
- B01F31/44—Mixers with shaking, oscillating, or vibrating mechanisms with stirrers performing an oscillatory, vibratory or shaking movement
- B01F31/449—Stirrers constructions
Landscapes
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Transmission Devices (AREA)
- Coloring Foods And Improving Nutritive Qualities (AREA)
Abstract
The present invention relates to a kind of nonlinear chaotic vibration rod, and it belongs to the nonlinear chaotic vibration rod that the industries such as a kind of building, food, medicine use.The present invention mainly solve existing vibrating head existing for vibrational excitation be difficult to reach material region, leakiness of vibrating and the uneven technical problem of stirring comprehensively.The technical solution adopted by the present invention is:A kind of nonlinear chaotic vibration rod, including motor, driving snake and chaos caput, wherein:Described motor and the drive end of driving snake connect, rotating shaft is provided with the inner chamber of the chaos caput, one end of the rotating shaft is fixed on the opening of chaos caput inner chamber, the other end of the rotating shaft is provided with eccentric block, the end of the driving snake couples with the fixing end of rotating shaft, L section curve of the chaos caput inwall from rotating shaft fixing end to eccentric block working chamber along axis section is one section of parabolic segment, first set L sections curve as:G (x)=acos (bx)+d, then verify whether the L sections curve meets that eccentric block reaches the critical condition of chaotic motion.
Description
Technical field
The present invention relates to a kind of nonlinear chaotic vibration rod, it belong to the industries such as a kind of building, food, medicine use it is non-
Linear chaotic vibration rod.
Background technology
Vibrating head in present society using extensive, therefore, the selection of vibrating head whether rationally, stirring extent it is whether uniform
Directly affect the quality of product.Vibrating head is a kind of mechanization tamping tool, by vibrating the bubble excluded in material, is eliminated
Gap.But in actual use, the vibration rail of existing vibrating head is mostly that circle or regular figure, vibrational excitation are difficult to arrive comprehensively
Up to material region, bubble is caused to discharge completely, leakiness of vibrating, stirring is uneven, has had a strong impact on stirring of vibrating
Quality.
Macroscopical nonlinear system that chaos refers to determine shown under certain condition it is uncertain, similar with
The motion of machine.Using the method for machinery, theoretical by non-linear chaotic dynamics introduces vibrating head technology, is allowed and shaken by chaos
Dynamic rod is activated to each region up to material to be mixed so that stirring is especially uniform, and Vibratory Mixing quality greatly improves.
The content of the invention
It is difficult to reach material region comprehensively present invention aim to address vibrational excitation existing for existing vibrating head, vibrates not
A kind of technical problem closely knit and that stirring is uneven, there is provided nonlinear chaotic vibration rod.
In order to solve the above technical problems, the technical solution adopted by the present invention is:
A kind of nonlinear chaotic vibration rod, including motor, driving snake and chaos caput, wherein:Described motor
It is connected with the drive end of driving snake, horizontal direction is provided with rotating shaft in the inner chamber of the chaos caput, and one end of the rotating shaft leads to
The opening that bearing is fixed on chaos caput inner chamber is crossed, the other end of the rotating shaft is provided with eccentric block, the end of the driving snake
End couples with the fixing end of rotating shaft, L of the chaos caput inwall from rotating shaft fixing end to eccentric block working chamber along axis section
Section curve is the parabolic segment of one section of small curvature, first set the L sections curve as:G (x)=acos (bx)+d, in formula, g (x) is turns
The contact position of axle and chaos caput inwall is to the distance of chaos caput axis, the contact position of x-rotating shaft and chaos caput inwall
Put to the distance of fixing end;Verify whether the L sections curve meets that eccentric block reaches the critical condition of chaotic motion again, if L sections are bent
Line meets the critical condition of chaotic motion, then the curve is that the chaos caput inwall works from rotating shaft fixing end to eccentric block
Cross section curve of the chamber along axis, otherwise, it is necessary to separately set curve, repeated authentication step, until meeting that chaotic motion occurs for eccentric block
Critical condition.
Further, whether the checking L sections curve meets that the step of eccentric block reaches chaotic motion critical condition is as follows:
(1) rigidity of rotating shaft is calculated
In formula:The amplitude of y-rotating shaft free end;a1- Taylor series expansion constant term;b1- Taylor series expansion is once
Term coefficient;c1The secondary term coefficient of-Taylor series expansion;
(2) nonlinear dynamic equation of rotating shaft free end eccentric block is established
In formula:P=me ω2, m-eccentric block quality;E-eccentric throw;ω-angular speed,N-rotating speed, c-resistance
Buddhist nun's coefficient;
(1) formula is brought into (2) formula to obtain
In formula:ε is one a small amount of,
(3) α is judged by Melnikov functions1、α2、α3Whether meet that eccentric block reaches the critical condition of chaotic motion
In formula:η1The imaginary part of-the first limit,
η2The imaginary part of-the second limit,
(4) dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determination eccentric blocks is combined
Whether system occurs chaos
The eccentric distance e and damped coefficient c that are determined by critical condition formula (3) are brought into (2-1) formula, obtained power
System is
Phase plane trajectory figure, timeamplitude map and Poincare are drawn according to the dynamical system using matlab softwares to map
Whether figure, send out with reference to the dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determination eccentric blocks
Raw chaos;The L sections curve is verified if chaos occurs to produce the cross section curve of chaotic motion.
The present invention uses above-mentioned technical proposal, makes chaos caput inwall curve that the critical bar of chaotic motion occur according to rotating shaft
Part designs, and in high-speed rotation, eccentric block is reached chaotic motion state, under the exciting effect of eccentric block, from
And barred body is driven to produce nonlinear chaotic vibration.Because non-linear chaotic motion never repeats in finite region inner orbit, utilize
Motion state is extremely complex and the characteristics of ergodic, solves usual vibrating poker and is difficult to reach material region comprehensively, vibrates
The leakiness technical problem uneven with stirring, compared with background technology, the present invention, which has, significantly improves Vibratory Mixing effect
Fruit and the advantages of improving mixing quality.
Brief description of the drawings
Fig. 1 is the structural representation of the present invention;
Fig. 2 is the structural representation of the chaos caput of nonlinear chaotic vibration rod of the present invention;
Fig. 3 is the phase plane trajectory figure of eccentric block of the present invention;
Fig. 4 is the timeamplitude map of eccentric block of the present invention;
Fig. 5 is the Poincare mapping graphs of eccentric block of the present invention.
Embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.
As depicted in figs. 1 and 2, a kind of nonlinear chaotic vibration rod in the present embodiment, including motor 1, driving snake 2
With chaos caput 7, wherein:Described motor 1 is connected with the drive end of driving snake 2, the inner chamber reclaimed water of the chaos caput 7
The opening of the inner chamber of chaos caput 7, the rotating shaft 4 square are fixed on by bearing 3 to one end provided with rotating shaft 4, the rotating shaft 4
The other end be provided with eccentric block 6, the end of the driving snake 2 couples with the fixing end of rotating shaft 4, the chaos caput inwall 5
L section curve of 6 working chambers along axis section is the parabolic segment of one section of small curvature from rotating shaft fixing end to eccentric block, is first set described
L section curves are:G (x)=acos (bx)+d, in formula, g (x) is rotating shaft 4 and the contact position of chaos caput inwall 5 to chaos rod
The distance of first 7 axis, the contact position of x-rotating shaft 4 and chaos caput inwall 5 from fixing end with a distance from, a, b, d are constant and ﹥
0;Verify whether the L sections curve meets that eccentric block 6 reaches the critical condition of chaotic motion again, if L section curves meet chaotic motion
Critical condition, then the curve is that the chaos caput inwall 5 section along axis of 6 working chambers from fixing end to eccentric block is bent
Line, otherwise, it is necessary to separately set curve, repeated authentication step, until meeting that the critical condition of chaotic motion occurs for eccentric block 6.
Whether the checking L sections curve meets that the step of eccentric block 6 reaches chaotic motion critical condition is as follows:
(1) rigidity of rotating shaft 4 is calculated
In formula:The amplitude of y-rotating shaft free end;a1- Taylor series expansion constant term;b1- Taylor series expansion is once
Term coefficient;c1The secondary term coefficient of-Taylor series expansion;
(2) nonlinear dynamic equation of the free end eccentric block of rotating shaft 4 is established
In formula:P=me ω2, m-eccentric block quality;E-eccentric throw;ω-angular speed,N-rotating speed, c-
Damped coefficient;
(1) formula is brought into (2) formula to obtain
In formula:ε is one a small amount of,
(3) α is judged by Melnikov functions1、α2、α3Whether meet that eccentric block reaches the critical condition of chaotic motion
In formula:η1The imaginary part of-the first limit,
η2The imaginary part of-the second limit,
(4) dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determination eccentric blocks is combined
Whether system occurs chaos
The eccentric distance e and damped coefficient c that are determined by critical condition formula (3) are brought into (2-1) formula, obtained power
System is
Phase plane trajectory figure, timeamplitude map and Poincare are drawn according to the dynamical system using matlab softwares to map
Whether figure, send out with reference to the dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determination eccentric blocks
Raw chaos;The L sections curve is verified if chaos occurs to produce the cross section curve of chaotic motion.
The present invention implementation verification process be:
In the vibrating head course of work, the inwall of chaos caput plays a part of support, thus chaos caput inwall is from rotating shaft 4
Section L section curve g (x) of the free end to eccentric block working chamber along axis should be smaller than flexure curvature of a curve, now sets L section curvilinear equations
For g (x)=acos (bx)+d (a, b, d are constant), then verified according to calculating.
(1) rigidity of L section curve g (x) and rotating shaft 4 under the support of the inwall of chaos caput 7 is determined
Y=g (x)=acos (bx)+d is substituted into (1-1) formula, and obtained by Taylor series expansion:
Will(1-2) formula of substitution, is obtained:
Order
:
It is set to meet following condition:
In the vibrating head course of work, the inwall of caput plays a part of supporting rotating shaft, thus curve g (x) should be than without branch
Rotating shaft flexure curvature of a curve is small during support, therefore according to the deflection curve equation without rotating shaft when supporting
Condition in (1-4) takes a=0.001, b=0.5, d=-0.0016.That is g (x)=0.001cos (0.5x) -0.0016.
By a, b, d value substitute into a respectively1, b1And c1, obtain:
By a1, b1, and c1Value substitute into (1) formula, obtain:
(2) nonlinear dynamic equation of the free end eccentric block of rotating shaft 4 is established
In formula:C-damped coefficient;P-exciting force, P=me ω2;E-eccentric throw;ω-angular speed;N-rotating speed;
Take eccentric block quality m=0.1kg, exciting force P=50N, rotating speed n=2840r/min.
(1-5) formula substitution (2) formula is obtained into the nonlinear dynamic equation of eccentric block is:
Wherein
(3) α is judged by Melnikov functions1、α2、α3Whether meet that eccentric block 6 reaches the critical condition of chaotic motion
Defining Melnikov functions is:
Wherein
Then
It is computed:
In formula:
The real part of two limits of ζ-integrand,
η1The imaginary part of-the first limit,
η2The imaginary part of-the second limit,
When Melnikov functions have zero point, the perturbed system considered has the chaos under Smale horseshoe meanings, i.e. M
(τ)=0:
To make equation (3-4) have root, it is necessary to make sin ω (t+ τ)≤1, and due toTherefore, eccentric block is obtained to reach
Chaotic motion α1、α2、α3The critical condition must being fulfilled for:
By α1、α2And α3Value bring into (3) formula, by (3-5) formula understand exist e and c cause (3) formula set up, i.e. α1、α2With
α3Value can meet eccentric block occur chaotic motion critical condition.
(4) dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determinations eccentric block 6 is combined
Whether system occurs chaos.
E=0.005653m, c=0.1 are taken according to the critical condition that chaos occurs, are brought into (2-1) formula, what is obtained is dynamic
Force system is:
Matlab softwares are utilized according to the dynamical system, it is y=-0.0006 to take primary condition,Draw
Whether phase plane trajectory figure, timeamplitude map and Poincare mapping graphs, the dynamical system of synthetic determination eccentric block occur chaos.
The Poincare mapping graphs shown in timeamplitude map and Fig. 5 shown in phase plane trajectory figure as shown in Figure 3, Fig. 4
As can be seen that time-history curves have no rule, phase plane trajectory is even more disorderly and unsystematic just like inverted 8, Poincare mappings, with reference to
This group of parameter has fallen in the chaotic region that Melnikov functions give, and therefore, dynamical system corresponding to determining type (2-1) has been enter into
Chaos state, rotating shaft are issued to chaotic motion in the support of rod shell inwall.Therefore determine g (x)=0.001cos (0.5x) -0.0016
For the present embodiment chaos caput inwall (5) from rotating shaft fixing end to eccentric block L section curve of (6) working chamber along axis section.
Claims (2)
1. a kind of nonlinear chaotic vibration rod, including motor (1), driving snake (2) and chaos caput (7), it is characterised in that:
Described motor (1) is connected with the drive end of driving snake (2), and horizontal direction is provided with the inner chamber of the chaos caput (7)
Rotating shaft (4), one end of the rotating shaft (4) are fixed on the opening of chaos caput (7) inner chamber, the rotating shaft (4) by bearing (3)
The other end be provided with eccentric block (6), the end of the driving snake (2) couples with the fixing end of rotating shaft (4), the chaos caput
Inwall (5) the L sections curve of (6) working chamber along axis section from rotating shaft fixing end to eccentric block is the parabola of one section of small curvature
Section, first set the L sections curve as:G (x)=a cos (bx)+d, in formula, g (x) is rotating shaft (4) and chaos caput inwall (5)
Contact position is to the distance of chaos caput (7) axis, the contact position of x-rotating shaft (4) and chaos caput inwall (5) to fixing end
Distance;Verify whether the L sections curve meets that eccentric block (6) reaches the critical condition of chaotic motion again, if L sections curve meets to mix
The critical condition of ignorant motion, then the curve is chaos caput inwall (5) (6) working chamber from rotating shaft fixing end to eccentric block
Along the cross section curve of axis, otherwise, it is necessary to separately set curve, repeated authentication step, until meeting that chaotic motion occurs for eccentric block (6)
Critical condition.
A kind of 2. nonlinear chaotic vibration rod according to claim 1, it is characterised in that:Whether the checking L sections curve
Meet that the step of eccentric block (6) reaches chaotic motion critical condition is as follows:
(1) rigidity of rotating shaft (4) is calculated
<mrow>
<mi>k</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>y</mi>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mi>y</mi>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>></mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:The amplitude of y-rotating shaft free end;a1- Taylor series expansion constant term;b1Term system of-Taylor series expansion
Number;c1The secondary term coefficient of-Taylor series expansion;
(2) nonlinear dynamic equation of rotating shaft (4) free end eccentric block is established
<mrow>
<mi>m</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>+</mo>
<mi>c</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mi>k</mi>
<mi>y</mi>
<mo>=</mo>
<mi>P</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&omega;</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:P=me ω2, m-eccentric block quality, e-eccentric throw, ω-angular speed,N-rotating speed, c-damping system
Number;
(1) formula is brought into (2) formula to obtain
<mrow>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mi>y</mi>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
<msup>
<mi>y</mi>
<mn>3</mn>
</msup>
<mo>=</mo>
<msup>
<mi>e&omega;</mi>
<mn>2</mn>
</msup>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&omega;</mi>
<mi>t</mi>
<mo>-</mo>
<mi>&epsiv;</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:ε is one a small amount of,
(3) α is judged by Melnikov functions1、α2、α3Whether meet that eccentric block reaches the critical condition of chaotic motion
<mrow>
<mfrac>
<mi>e</mi>
<mi>c</mi>
</mfrac>
<mo>&GreaterEqual;</mo>
<mfrac>
<mrow>
<msqrt>
<mn>2</mn>
</msqrt>
<mrow>
<mo>(</mo>
<mrow>
<mi>c</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&omega;</mi>
</mrow>
<msqrt>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</msqrt>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>3</mn>
<msup>
<mi>&pi;&omega;</mi>
<mn>3</mn>
</msup>
<mi>m</mi>
<msqrt>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
</msqrt>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>ch&omega;&eta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>ch&omega;&eta;</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<msqrt>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</msqrt>
</mrow>
<mrow>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mn>9</mn>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<mn>9</mn>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
<msqrt>
<mrow>
<mn>2</mn>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mi>arcsin</mi>
<mfrac>
<mrow>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
<msqrt>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mn>9</mn>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
</mrow>
</msqrt>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula:
η1The imaginary part of-the first limit,
η2The imaginary part of-the second limit,
(4) it is with reference to the dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determination eccentric blocks
No generation chaos
The eccentric distance e and damped coefficient c that are determined by critical condition formula (3) are brought into (2-1) formula, obtained dynamical system
For
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>Z</mi>
<mo>=</mo>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>Z</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>&alpha;</mi>
<mn>3</mn>
</msub>
<msup>
<mi>y</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>e&omega;</mi>
<mn>2</mn>
</msup>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&omega;</mi>
<mi>t</mi>
<mo>-</mo>
<mi>&epsiv;</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Phase plane trajectory figure, timeamplitude map and Poincare mapping graphs are drawn using matlab softwares according to the dynamical system,
Whether mixed with reference to the dynamical system of phase plane trajectory figure, timeamplitude map and Poincare mapping graph synthetic determination eccentric blocks
It is ignorant;The L sections curve is verified if chaos occurs to produce the cross section curve of chaotic motion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710657775.XA CN107469702B (en) | 2017-08-03 | 2017-08-03 | Nonlinear chaotic vibrating spear |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710657775.XA CN107469702B (en) | 2017-08-03 | 2017-08-03 | Nonlinear chaotic vibrating spear |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107469702A true CN107469702A (en) | 2017-12-15 |
CN107469702B CN107469702B (en) | 2020-07-03 |
Family
ID=60598701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710657775.XA Active CN107469702B (en) | 2017-08-03 | 2017-08-03 | Nonlinear chaotic vibrating spear |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107469702B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113164950A (en) * | 2018-11-05 | 2021-07-23 | Hycor生物医学有限责任公司 | Apparatus and method for mixing fluids or media by vibrating a pipette using a non-concentric mass |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004126488A (en) * | 2002-10-02 | 2004-04-22 | Koichi Kameoka | Chaos vibration generating device |
CH694433A5 (en) * | 1999-07-28 | 2005-01-14 | Fabrizi Gmbh | Vibration assembly, to settle poured concrete in a mold, has a shrouding around the on/off switch and connection to the frequency converter housing to protect the electrical/electronic components against soiling |
CN101038304A (en) * | 2007-03-16 | 2007-09-19 | 哈尔滨工程大学 | Segmenting detecting method for pattern domain of Duffing vibrator phase diagram |
CN101793306A (en) * | 2010-03-12 | 2010-08-04 | 浣石 | Continuously variable-rigidity chaos vibration isolation device |
CN101972689A (en) * | 2010-10-08 | 2011-02-16 | 南京工程学院 | Secondary variable pendulum chaotic vibrating mill |
EP2347820A2 (en) * | 2008-10-21 | 2011-07-27 | Industry-Academic Cooperation Foundation Gyeongsang National University | An agitator vessel using a baffle and an agitator with improved stirring ability including the same |
CN201953080U (en) * | 2010-12-20 | 2011-08-31 | 王继龙 | Concrete vibrating rod |
JP2012058107A (en) * | 2010-09-09 | 2012-03-22 | Nittetsu Elex Co Ltd | Rotary machine abnormality diagnosis method |
CN202220446U (en) * | 2011-08-29 | 2012-05-16 | 孙曙光 | Pluggable concrete vibrating bar |
CN102622928A (en) * | 2012-03-01 | 2012-08-01 | 合肥工业大学 | Nonlinear buoyancy chaotic experiment device |
CN202778923U (en) * | 2012-08-04 | 2013-03-13 | 昆明理工大学 | End face type chaos vibration exciter |
CN103150483A (en) * | 2013-03-29 | 2013-06-12 | 王少夫 | Wiener noise-contained two-dimensional chaotic system |
CN103161866A (en) * | 2013-03-20 | 2013-06-19 | 王少夫 | Vibration damping chaotic system |
US20140169407A1 (en) * | 2012-03-26 | 2014-06-19 | John David Terry | Method and Apparatus for Multiple Signal Aggregation and Reception in Digital Chaos Network |
CN104516858A (en) * | 2014-12-31 | 2015-04-15 | 石家庄铁道大学 | Phase diagram matrix method for nonlinear dynamic behavior analysis |
CN105396496A (en) * | 2015-12-23 | 2016-03-16 | 哈尔滨杰曼机电设备有限公司 | Particulate material vibration stirring apparatus and vibration stirring method |
CN105822070A (en) * | 2016-04-07 | 2016-08-03 | 成都固特机械有限责任公司 | Eccentric concrete vibrating bar |
CN206245752U (en) * | 2016-09-30 | 2017-06-13 | 昆明理工大学 | A kind of plug-in type chaotic vibration vibrator |
-
2017
- 2017-08-03 CN CN201710657775.XA patent/CN107469702B/en active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CH694433A5 (en) * | 1999-07-28 | 2005-01-14 | Fabrizi Gmbh | Vibration assembly, to settle poured concrete in a mold, has a shrouding around the on/off switch and connection to the frequency converter housing to protect the electrical/electronic components against soiling |
JP2004126488A (en) * | 2002-10-02 | 2004-04-22 | Koichi Kameoka | Chaos vibration generating device |
CN101038304A (en) * | 2007-03-16 | 2007-09-19 | 哈尔滨工程大学 | Segmenting detecting method for pattern domain of Duffing vibrator phase diagram |
EP2347820A2 (en) * | 2008-10-21 | 2011-07-27 | Industry-Academic Cooperation Foundation Gyeongsang National University | An agitator vessel using a baffle and an agitator with improved stirring ability including the same |
CN101793306A (en) * | 2010-03-12 | 2010-08-04 | 浣石 | Continuously variable-rigidity chaos vibration isolation device |
JP2012058107A (en) * | 2010-09-09 | 2012-03-22 | Nittetsu Elex Co Ltd | Rotary machine abnormality diagnosis method |
CN101972689A (en) * | 2010-10-08 | 2011-02-16 | 南京工程学院 | Secondary variable pendulum chaotic vibrating mill |
CN201953080U (en) * | 2010-12-20 | 2011-08-31 | 王继龙 | Concrete vibrating rod |
CN202220446U (en) * | 2011-08-29 | 2012-05-16 | 孙曙光 | Pluggable concrete vibrating bar |
CN102622928A (en) * | 2012-03-01 | 2012-08-01 | 合肥工业大学 | Nonlinear buoyancy chaotic experiment device |
US20140169407A1 (en) * | 2012-03-26 | 2014-06-19 | John David Terry | Method and Apparatus for Multiple Signal Aggregation and Reception in Digital Chaos Network |
CN202778923U (en) * | 2012-08-04 | 2013-03-13 | 昆明理工大学 | End face type chaos vibration exciter |
CN103161866A (en) * | 2013-03-20 | 2013-06-19 | 王少夫 | Vibration damping chaotic system |
CN103150483A (en) * | 2013-03-29 | 2013-06-12 | 王少夫 | Wiener noise-contained two-dimensional chaotic system |
CN104516858A (en) * | 2014-12-31 | 2015-04-15 | 石家庄铁道大学 | Phase diagram matrix method for nonlinear dynamic behavior analysis |
CN105396496A (en) * | 2015-12-23 | 2016-03-16 | 哈尔滨杰曼机电设备有限公司 | Particulate material vibration stirring apparatus and vibration stirring method |
CN105822070A (en) * | 2016-04-07 | 2016-08-03 | 成都固特机械有限责任公司 | Eccentric concrete vibrating bar |
CN206245752U (en) * | 2016-09-30 | 2017-06-13 | 昆明理工大学 | A kind of plug-in type chaotic vibration vibrator |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113164950A (en) * | 2018-11-05 | 2021-07-23 | Hycor生物医学有限责任公司 | Apparatus and method for mixing fluids or media by vibrating a pipette using a non-concentric mass |
Also Published As
Publication number | Publication date |
---|---|
CN107469702B (en) | 2020-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107469702A (en) | A kind of nonlinear chaotic vibration rod | |
CN105986803B (en) | Underground vibrating monitoring method | |
CN106638558B (en) | Civil engineering high pressure jet grouting device | |
CN104285011A (en) | A variable moment resonance-free vibro hammer | |
CN104657611A (en) | Ship whole elastic-plastic movement response predication method and system under action of underwater explosion | |
CN103352414A (en) | Method for optimizing parameters of multi-frequency multi-amplitude vibratory roller | |
CN202900251U (en) | Drill hole angle-locating device for channeling hole | |
CN103511553B (en) | A kind of bat is shaken the control apparatus of noise, method, system and engineering machinery | |
JP5516531B2 (en) | Ground improvement excavator | |
CN104935229A (en) | Method of obtaining energy storage volute spiral spring real-time rotational inertia | |
CN106639123B (en) | Steel construction for mixed mud building block | |
CN107110105A (en) | Method and system for adjusting the quality in Wave-activated power generation equipment and the torque of spinning roller rotor | |
CN101550909A (en) | A wind power generator and its dynamic balance control system | |
JP2019078012A (en) | Pile placement management system | |
EP2067533A1 (en) | Vibrator for a vibratory pile driver | |
US3383531A (en) | Electric one-way unharmonical vibarator | |
CN106917459A (en) | A kind of half active mono-pendulum type tuned mass damper | |
CN103753691B (en) | A kind of even number parallel multi-shaft vertical orientation vibrational system | |
CN107497655A (en) | One kind reverses chaotic vibration vibrator | |
EP2392413A2 (en) | Vibration ram | |
CA1038649A (en) | Method of producing periodic shocks in one direction | |
CN103608518A (en) | Surface compactor and method of operation | |
Jebur et al. | Experimental investigation of under reamed pile subjected to dynamic loading in Sandy soil | |
Soldatenko et al. | Clarification of the methods used for calculating power of sieve separators | |
CN104657544B (en) | A kind of determination method of wind turbine gearbox torsion stiffness |
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 |