CN112620101B - 单质体四机倍频自同步驱动振动机及其参数确定方法 - Google Patents
单质体四机倍频自同步驱动振动机及其参数确定方法 Download PDFInfo
- Publication number
- CN112620101B CN112620101B CN202011353808.XA CN202011353808A CN112620101B CN 112620101 B CN112620101 B CN 112620101B CN 202011353808 A CN202011353808 A CN 202011353808A CN 112620101 B CN112620101 B CN 112620101B
- Authority
- CN
- China
- Prior art keywords
- mass
- vibration
- frequency
- formula
- vibration exciters
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000001360 synchronised effect Effects 0.000 claims abstract description 23
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 230000006698 induction Effects 0.000 claims abstract description 4
- 210000002706 plastid Anatomy 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 abstract description 4
- 239000010802 sludge Substances 0.000 abstract description 4
- 239000002002 slurry Substances 0.000 abstract description 3
- 239000012535 impurity Substances 0.000 abstract 1
- 238000012795 verification Methods 0.000 abstract 1
- 238000006073 displacement reaction Methods 0.000 description 23
- 238000004088 simulation Methods 0.000 description 15
- 238000012360 testing method Methods 0.000 description 9
- 230000004044 response Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000012216 screening Methods 0.000 description 4
- 230000035939 shock Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005056 compaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011178 precast concrete Substances 0.000 description 3
- 239000004567 concrete Substances 0.000 description 2
- 230000018044 dehydration Effects 0.000 description 2
- 238000006297 dehydration reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000007493 shaping process Methods 0.000 description 2
- 238000007873 sieving Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005495 investment casting Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011819 refractory material Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B07—SEPARATING SOLIDS FROM SOLIDS; SORTING
- B07B—SEPARATING SOLIDS FROM SOLIDS BY SIEVING, SCREENING, SIFTING OR BY USING GAS CURRENTS; SEPARATING BY OTHER DRY METHODS APPLICABLE TO BULK MATERIAL, e.g. LOOSE ARTICLES FIT TO BE HANDLED LIKE BULK MATERIAL
- B07B1/00—Sieving, screening, sifting, or sorting solid materials using networks, gratings, grids, or the like
- B07B1/42—Drive mechanisms, regulating or controlling devices, or balancing devices, specially adapted for screens
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B07—SEPARATING SOLIDS FROM SOLIDS; SORTING
- B07B—SEPARATING SOLIDS FROM SOLIDS BY SIEVING, SCREENING, SIFTING OR BY USING GAS CURRENTS; SEPARATING BY OTHER DRY METHODS APPLICABLE TO BULK MATERIAL, e.g. LOOSE ARTICLES FIT TO BE HANDLED LIKE BULK MATERIAL
- B07B1/00—Sieving, screening, sifting, or sorting solid materials using networks, gratings, grids, or the like
- B07B1/28—Moving screens not otherwise provided for, e.g. swinging, reciprocating, rocking, tilting or wobbling screens
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Apparatuses For Generation Of Mechanical Vibrations (AREA)
Abstract
本发明公开了单质体四机倍频自同步驱动振动机及其参数确定方法,该振动系统中质体通过弹簧A和弹簧B与地基相连,弹簧对称分布于质体上;四个激振器两两分别安装在质体的上侧和下侧,每个激振器中各有一偏心转子,偏心转子由各自的感应电动机驱动,分别绕着旋转轴线中心旋转,同侧的两个频率相同的激振器关于y轴对称分布,但旋转方向相反,自倍频同步振动驱动设备工作;利用倍频振动自同步原理,通过建立动力学模型和运动微分方程、进行倍频同步理论解析、推导二倍频以及三倍频同步性条件和稳定性条件,并进行实验验证。能够有效提高振动系统工作效率。特别适用于黏湿性大等难筛分物料的分级,也适合粘性大杂质多的工程泥浆及污泥的脱水等。
Description
发明领域
本发明属于振动装置技术领域,涉及一种单质体四机倍频自同步驱动双直线运动轨迹振动机及参数确定方法。
背景技术
振动筛分/脱水/密实/成型设备,是一种利用振动对物料进行分级,实现泥浆或污泥的固液分离,以及实现预制混凝土及精密铸造构件的振动密实与成型,适用于砂石骨料、矿山、钢厂、食品、化工、石油、建筑、隧道工程及轨道交通等行业。该类设备主要实现各类干式物料的分级,固液的分离以及泥浆或污泥的脱水,也能实现建筑及轨道交通行业用预制混凝土的振动密实与成型等。本发明属于振动倍频驱动实现系统双频运动的振动机械。普通的单频振动机械会产生许多如下问题:
1.单频驱动振动机械,如振动筛很容易产生筛堵或筛糊现象,导致筛分效率较低;振动脱水筛容易造成脱水效率低下;混凝土密实成型效果差及效率低下。
2.采取提高抛掷指数来提高设备振动强度,某种意义上可以提高设备功能,但会使整个设备寿命降低,不利于工业生产中的使用性能要求。
3.为了提高工作效率,现在工程上普遍采用其他一些方式来改善机器功能及性能,但这些方法会导致设备结构复杂、体积大、加工成本高。
随着倍频同步理论的不断完善,有必要应用先进的倍频振动同步技术,设计一款既能提高设备功能,又能保证设备性能要求,且结构紧凑的振动机械,使其既能提高效率又能长期稳定运行。
发明内容
本发明为了克服现有技术中存在的问题,通过以下技术方案实现:
单质体四机倍频自同步驱动振动机的动力学模型包括:四个激振器、一个质体、弹簧A与弹簧B;其中质体通过弹簧A与弹簧B与地基相连,弹簧对称分布于质体上;四个激振器两两分别安装在质体的上侧和下侧,每个激振器中各有一偏心转子,偏心转子由各自的感应电动机驱动,分别绕着旋转轴线中心旋转,同侧的两个频率相同的激振器关于y轴对称分布,但旋转方向相反,倍频自同步驱动,实现设备双频双直线轨迹运动功能;以此提高振动筛分/脱水/密实/成型设备的工作效率,改善该类设备的工作质量。
所述振动机的四个激振器的参数确定方法,包括如下步骤:
步骤1,建立动力学模型和系统运动微分方程
如图1所示,建立如图所示的坐标系。四个激振器分别绕着旋转中心轴o1,o2,o3和o4旋转。分别是四个转子(URs)的旋转角。而且,四个激振器与x轴的夹角分别用β1,β2,β3和β4表示。根据这个模型,整个系统展现出三个自由度:直线运动x,y和摆动角ψ。
根据Lagrange方程,得系统的运动微分方程如下:
其中
式中,
m——质体质量;
mj——激振器j偏心质量,m1=m2,m3=m4;
Ji——激振器i转动惯量,J1=J2,J3=J4;
Jm——质体自身转动惯量;
g——重力加速度;
lj——激振器j回转轴心oj到质体中心O的距离,l1=l2,l3=l4;
rj——激振器j偏心距,r1=r2,r3=r4;
fi——电机i轴阻尼系数,i=1,2,3,4;
le——系统绕质心当量回转半径;
Tei——电机i电磁输出转矩,i=1,2,3,4;
βj——激振器j回转轴心oj到质体中心O的连线与x轴夹角,β1+β2=-π,β3+β4=π;
kx,ky,kψ——x,y和ψ方向弹簧刚度;
fx,fy,fψ——x,y和ψ方向阻尼系数;
步骤2,四激振器倍频同步理论解析
假设四个激振器实现同步运转,激振器1和2的转速相同,激振器3和4的转速是他们的整数倍(分别用n3、n4表示),以顺时针旋转方向为正方向,有:
考虑到激振器是以远远高于系统固有频率的速度旋转,将式(1)等号左边的第二项和第三项省略掉。此外,考虑到x,y和ψ都是很小的并且四个激振器的运转状态接近稳定,代入式(1)的最后一个表达式,去掉式(1)的x,y和ψ并且取其二阶导,得到的近似表达式为:
其中
σ1=σ3=-1,σ2=σ4=1
式中,小参数ε是激振器偏心块质量与系统总质量的比值。σi(i=1,2,3,4)值的正负表示偏心块的旋转方向,正值代表顺时针方向旋转,负值代表逆时针方向旋转。
把旋转相位角作如下表示:
其中
τ=ωt,n1=n2=1
式中,Δi是由于系统的运动而随着激振器偏心转子的产生阶段而缓慢变化的函数。
将式(4)代入式(3)中得:
其中,
ψij +=(σini+σjnj)τ+σiΔi+σjΔj+βi+βj,ψij -=(σini-σjnj)τ+σiΔi-σjΔj+βi-βj式(5)为激振器实现同步的基本表达式。
把式(5)写成标准形式:
式(5)和(6)关于未知参数Δi和νi建立的一阶微分方程表达式如下:
其中
σini+σjnj≠0时pij=1/(σini+σjnj),σini+σjnj=0时pij=0
σini-σjnj≠0时qij=1/(σini-σjnj),σini-σjnj=0时qij=0
同样改进第二近似解:
Δi=Δi,i=1,2,3,4
将式(9)代入式(7)等号的左边,Ωi和Δi作为固定值并取关于τ=0~2π上的平均值。考虑到相同转速的激振器会互相反向旋转,得到如下关系式:
其中
σini+σjnj=0时us=1,ψij (1)=σiΔi+σjΔj+βi+βj否则us=0
σini+σrnr=0时uh=1,ψij (2)=σiΔi+σrΔr+βi+βr否则uh=0
σini-2σjnj=0时ul=1,γij (1)=σiΔi-2σjΔj+βi-βj否则ul=0
σini+2σrnr=0时um=1,γij (2)=σiΔi+2σjΔj+βi+βj否则um=0
σini-2σjnj+σrnr=0时ud1=1,ηijr (1)=σiΔi-2σjΔj+σrΔr+βi-2βj+βr否则ud1=0
σini-2σjnj-σrnr=0时ud2=1,ηijr (2)=σiΔi-2σjΔj-σrΔr+βi-2βj-βr否则ud2=0
σini+2σjnj+σrnr=0时ud3=1,ηijr (3)=σiΔi+2σjΔj+σrΔr+βi+2βj+βr否则ud3=0
σini+2σjnj-σrnr=0时ud4=1,ηijr (4)=σiΔi+2σjΔj-σrΔr+βi+2βj-βr否则ud4=0
这里,通过求出稳定解,由式(10)可知在第项和第ε项中可以确定转速相等的激振器间相位关系,可知次的第二项是在n3=n4=2的情况下有转速比为1:2的同步相位关系式,第三项是在n3=n4=3的情况下有1:3的同步相位关系式。因此,首先可以通过取到式(10)的ε次为止来确定转速相等激振器间的相位关系,将该相位关系代入式(10),然后从次项求出转速比为1:2或1:3的激振器之间的同步相位关系。另外,考虑到系统的结构是对称的,所以有:
a12=a21=1,a34=a43,a13=a14=a23=a24,α1=α2,α3=α4,
k1=k2,k3=k4,l1=l2,l3=l4,
A11=A12=A21=A22=A1,A33=A34=A43=A44=A2,A13=A14=A24=A23
步骤三,推导四激振器同步性及稳定性条件
(1)当转速相等时,取到式(10)的ε项为止,得:
当系统处于稳定状态时,式(11)中参数的表达式为:
因此,激振器实现相同频率(转速比为1:1)的同步性条件为:
在式(14)中,考虑到稳定状态,二倍频同步条件写为下式:
假定初始相位Δi0和Ωi0都具有小偏差,需要做如下设定:
Δi=Δi0+δi,Ωi=Ωi0+ξi,i=1,2,3,4 (16)
将式(16)代入式(11)中得到系统微分方程表达式为:
其中
整理式(17)得到关于δi(i=1,2,3,4)的表达式为:
取特征值为λ,得到的特征方程如下:
根据Routh-Hurwitz准则分析并整理得到以下稳定性判据:
式中,ε,α1 (1),α3 (1),a31,a34,k1,A1,A2均大于0,且ε值很小,ε2K21和ε2K23无限接近于0,所以在计算中可认为εA1+ε2K21,4εA2+ε2K23,1.5ε2a31k1均为正值,经过分析可确定:
cos(Δ20-Δ10)>0 (21)
(3)当n3=n4=3时,激振器3和4的稳定转速是激振器1和2的三倍,系统实现三倍频同步。同步条件表达式(13)改为:
考虑到稳定状态,转速比为1:3的同步性条件可写为下式:
为了寻求稳定相位角,与二倍频同步分析方法一样,求出系统在稳定状态下的微分方程表达式,引出它的特征方程为:
其中
根据Routh-Hurwitz准则分析得到稳定性判据如下:
和二倍频一样方法整理求解式(25)得:
cos(Δ20-Δ10)>0 (26)
本发明的有益效果:
1)本专利采用四机倍频自同步驱动,无论是2倍频还是3倍频,都能实现机体的双频双直线运动轨迹,能够有效提高设备的处理效率和工作质量,如提高筛分及脱水效率,有效提高振动密实成型效果和质量等(例如,增强预制混凝土或精密铸造构件的密实度,改善生产量,同时也能提高设备振动密实成型作业的工作效率等)。
2)筛机双频双直线轨迹运动,有效提高筛机的处理量和效率,特别适合黏湿性物料的分级,工程泥浆或污泥的脱水,以及密实成型效果和质量。
附图说明
图1为反向回转四机驱动单质体机械系统动力学模型图。
图中:1.激振器3;2.质体;3.激振器4;4.弹簧A;5.激振器2;6.激振器1;7.弹簧B。
图中各参数含义:
oxy--绝对坐标系
O--整个系统的中心;
O1--激振器1旋转中心;
O2--激振器2旋转中心;
O3--激振器3旋转中心;
O4--激振器4旋转中心;
m1--激振器1质量;
m2--激振器2质量;
m3--激振器3质量;
m4--激振器4质量;
r1--激振器1偏心距;
r2--激振器2偏心距;
r3--激振器3偏心距;
r4--激振器4偏心距;
m--质体质量;
kx--弹簧A在x方向的刚度系数;
ky--弹簧B在y方向的刚度系数;
β1--激振器1与x轴的夹角;
β2--激振器2与x轴的夹角;
β3--激振器3与x轴的夹角;
β4--激振器4与x轴的夹角;
l1--激振器1回转轴心o1到质体中心O的距离;
l2--激振器2回转轴心o2到质体中心O的距离;
l3--激振器3回转轴心o3到质体中心O的距离
l4--激振器4回转轴心o4到质体中心O的距离
lx--弹簧A与质体连接点到系统中心O的距离;
ly--弹簧B与质体连接点到系统中心O的距离;
ψ--质体绕中心轴摆动的角度。
图2为在不同rl1和rl3下四个激振器的二倍频同步稳定性区域(η1=η2=η3=η4=1)。
图3为在不同rl1和rl3下四个激振器的三倍频同步稳定性区域(η1=η2=η3=η4=1)。
图4为η1=η2=1,η3=η4=0.5时二倍频条件下超远共振状态下的仿真结果:
(a)四电机转速;
(b)电机输出转矩;
(c)激振器1和2相位差;
(d)激振器2和3相位差;
(e)激振器3和4相位差;
(f)x和y方向位移;
(g)摆动角。
图5为η1=η2=η3=η4=1时二倍频条件下超远共振状态下的仿真结果:
(a)激振器1和2相位差;
(b)激振器2和3相位差;
(c)激振器3和4相位差;
(d)x和y方向位移;
(e)摆动角。
图6为η1=η2=1,η3=η4=0.5时三倍频条件下超远共振状态下的仿真结果:
(a)四电机转速;
(b)电机输出转矩;
(c)激振器1和2相位差;
(d)激振器2和3相位差;
(e)激振器3和4相位差;
(f)x和y方向位移;
(g)摆动角。
图7为η1=η2=η3=η4=1时三倍频条件下超远共振状态下的仿真结果:
(a)激振器1和2相位差;
(b)激振器2和3相位差;
(c)激振器3和4相位差;
(d)x和y方向位移;
(e)摆动角。
图8为四机驱动振动同步试验台。
图9为η1=η2=η3=η4=1时四激振器二倍频同步试验结果。
具体实施方案
实施例1:
假定振动系统的参数:m=1430kg,m0=20kg,J=1050kg·m2,kx=ky=180kN/m,kψ=120kN/rad,fx=fy=3.83kN·s/m,fψ=3.2kN·s/rad,r1=r2=r3=r4=0.15m,l1=l2=l3=l4=1.1m,β1=-2π/3,β2=-π/3,β3=2π/3,β4=π/3,设激振器偏心块质量为mi=ηim0(i=1,2,3,4,m0为激振器偏心块的标准质量),根据振动系统的参数,容易求出主要固有频率:ωn=ωx=ωy=10.92rad/s,ωψ=10.69rad/s。电动机的类型:三相鼠笼式,50Hz,380V,6-pole,0.75kW,额定转速:980r/min。设置了电机参数:转子电阻Rr=3.40Ω,定子电阻Rs=3.35Ω互感系数Lm=164mH,转子电感Lr=170mH,定子电感Ls=170mH。
(a)四个激振器在二倍频同步条件下的稳定区域
为了便于讨论激振器在倍频同步条件下的动力学特性,引入了无量纲参数rli=li/le,i=1,2,3,4,因为系统的结构是对称的,有rl1=rl2,rl3=rl4。根据理论分析可知相同转速的激振器1和2,激振器3和4的稳定相位差均为0度左右。由二倍频同步稳定性条件表达式(20)分析可以得到二倍频同步稳定性区域,如图2所示。分别改变无量纲参数rl1和rl3,可以得到两组特性曲线。在曲线下方,二倍频相位差及的值稳定在0°附近;曲线上方,二倍频相位差值稳定在180°附近。此外,还可以看出小参数ε随rl1及rl3的增大呈增长趋势,并且激振器回转中心与质心间距离越大,增长趋势越明显。
(b)四个激振器在三倍频同步条件下的稳定区域
同样方法分析三倍频同步,根据稳定性判据表达式(25)可以得到三倍频同步稳定区域,如图3所示。改变无量纲参数rl1和rl3的值,得到两组小参数ε随rl1和rl3的变化曲线。根据分析可知在曲线下方,三倍频相位关系及稳定在0°附近;而曲线上方,三倍频相位差值稳定在180°附近。其余分析和二倍频一样,这里就不再赘述。
实施例2
(a)四个激振器在二倍频同步条件下的仿真
设定参数m0=20kg,rl1=rl3=1.3,其余系统参数和电机参数在上述部分已给出。仿真过程中,需要通过控制感应电机的工作频率来得到不同的转速,因此模型中电机1和2的工作频率均为25Hz,电机3和4的工作频率均为50Hz。改变激振器间的质量关系获得系统在不同偏心质量距条件下的动力学特性,应用Runge-Kutta程序给出了二倍频条件下的电脑仿真结果,两组仿真均在20s时给电机2一个π/4的干扰。
如图4所示,仿真中激振器间的质量关系为η1=η2=1,η3=η4=0.5,电机1和2的稳定转速约为490r/min,电机3和4的稳定转速约为980r/min,恰好为电机1和2的二倍,激振器间实现了二倍频同步。此时ω≈98.4rad/s,ω>ωn,因此系统在超远共振状态下运转。
从图4(c)(d)(e)可以看出,大约7s后相位差达到稳定。激振器间稳定相位关系为:干扰后,激振器间的稳定相位差出现短时间的波动,然后迅速恢复到原状态,表明该状态下系统是稳定运转的,其运动状态不受外部干扰的影响。
图4(f)表示质体在x和y方向的位移曲线图,可看出x方向上稳态时质体的振幅约为0mm,y方向上质体的振幅约为3.0mm。根据图4(g)可看出摆动角的大小几乎为0度,可认为质体没有摆动,这表明系统的主要运动形式为y方向的直线运动,其具体运动形式可以根据图4(f)的放大图清晰地看出。此外,在20s时给电机2加一个π/4相位的干扰,x和y方向的位移有一个迅速增大,然后又回到原来稳定状态,其最大位移值保持不变。
改变激振器间的质量关系η1=η2=η3=η4=1,即四个激振器的偏心块质量相等,得到的仿真结果如图5所示。电机的转速不变,其稳定相位关系也基本不变,仍然是说明激振器的质量几乎不影响相位差的稳定值。不同的是系统在y方向的振幅有所改变,此时振幅约为4.0mm,这表明系统的位移与激振器的偏心质量距有关系,质量距越大,系统的最大位移越大。观察图5(d)的放大图,可以看到质体在y方向的运动过程中会有位移的突然增大,这对工程中是有用的,可以实现更大的振幅,更好地为振动摇筛的设计提供理论指导。
(b)四个激振器在三倍频同步条件下的仿真
改变电机的工作频率,使模型中质体下方电机1和2的频率均为16Hz,质体上方电机3和4的频率为48Hz,其余仿真参数不变,得到的仿真结果展示在图6和图7中。如图6(a)和7(a)所示,电机1和2的稳定转速约为320r/min,电机3和4的转速约为电机1和2的三倍,即:960r/min。在20s时同样给电机2一个π/4相位的干扰,如图6(c)(d)(e)所示,大约7s后系统达到稳定,此时激振器间的稳定相位关系为:和二倍频同步相比变化不大。这表明不管是二倍频同步还是三倍频同步,激振器间的稳定相位关系基本不受影响,主要变化的是质体的位移。
质体在x方向的位移几乎为零,可以忽略,因此系统实现的是y方向上的直线运动,其运动的最大位移约为3.3mm。改变激振器的偏心质量,使四个激振器的质量相等,在35s给电机2加一个干扰,得到的仿真结果如图7所示。容易发现质体在y方向的振幅变大了,约为4.7mm,表明激振器的偏心质量距越大,质体在y方向的位移越大。同样地,当四个激振器偏心转子质量相等时,质体在y方向运动过程中会有位移突然的增加,类似于一个冲击波。这种响应冲击波对于工程上是有用的,可以实现振幅的迅速增大,为新型振动成型密实装备的设计提供参考,实现对预制混泥土的构建或耐火材料的密实成型。
实施例3
为了进一步验证理论及数值分析的正确性,根据该模型搭建了试验台,进行试验研究。四机振动同步试验台如图8所示,选取的四个电机型号上述部分已给出,其额定转速均为980r/min。双频振动同步试验系统的参数为:m=345kg,m0=4kg,J=44.5kg·m2,kx=ky=110.87kN/m,kψ=12.65kN/rad,fx=fy=0.37kN·s/m,fψ=0.22kN·s/rad,r1=r2=r3=r4=0.05m,l1=l2=l3=l4=0.45m,β1=-3π/4,β2=-π/4,β3=3π/4,β4=π/4,rl1=rl3=1.3。根据振动系统的参数,容易求出主要固有频率:ωx′=ωy′=17.52rad/s,ωψ′=16.86rad/s。电机1和2对称安装在质体的下方,其旋转方向相反;电机3和4对称安装在质体的上方,旋转方向同样相反。试验中通过变频器调节电机的供电频率来得到不同的电机转速,可通过调整其偏心块夹角的大小来调整激振器的激振力,一般有夹角越大,其偏心力越大。利用霍尔传感器的脉冲触发点来测电机的转速及相位,质体的位移可由加速度传感器间接测得,把得到的加速度曲线进行二次积分可得位移曲线。将智能信号分析仪采集的数据导入Matlab软件进行编程处理,最后通过OriginPro 8成像,可以得到转速、相位差、位移响应图等。
试验过程中,调节电机1和2的供电频率为25Hz,电机3和4的供电频率为50Hz,采样时间为100s,得到二倍频同步试验结果如图9所示。电机转速展示在9(a)中,可以清晰地看出电机1和2的转速相同,约为500r/min,电机3和4的转速均为1000r/min左右,是电机1和2转速的二倍。
四个电机刚开始启动时,各电机的转动惯量相等,其转速也接近相同。当电机转速到达共振点时会激起共振响应,此时振幅达到最大。一段时间后,通过调节电机间的相位差来使各电机的负载扭矩达到同步,从而使电机转速快速稳定下来,激振器间实现同步并稳定运转。
根据图9(b)(c)(d)可知,二倍频相位差达到稳定状态所需的过渡时间较长,大约30s后才能稳定;而相同频率间相位差只需10s左右就能稳定下来。稳定运转时,低频电机1和2间的相位差大约稳定在0°,高频和低频电机间的二倍频相位差大约稳定在30°左右,两个高频电机3和4的相位差稳定在0°附近,与仿真结果相比,略有偏差,但定性上是一致的。出现偏差的原因可能是即使选取四个电机型号完全相同,其输出转矩也不可能完全相同。这对于转速不同的电机间相位差影响比较大,使其不能恰好稳定在0度左右,还有可能是霍尔传感器的布置不精确导致测得的相位出现偏差。
由加速度传感器采集并经过二次积分可以得到试验台机体在x,y和ψ方向的位移,其随时间的变化曲线分别如图9(e)(f)(g)所示。根据图9(f)的放大图可明显看出稳态时质体在y方向的运动状态,其运动形式大致与仿真结果相同。而x方向的最大位移约为0.12mm,摆动最大角约为0.5°,均可以被忽略,因此可认为稳态时质体的运动形式为y方向的直线运动。在稳定运转过程中,y方向会有类似脉冲响应冲击波的产生,导致位移迅速增大,这正是工程中所需要的,可以为振动成型装备的密实过程提供借鉴。
Claims (2)
1.单质体四机倍频自同步驱动振动机,其特征在于,该振动机的动力学模型包括:四个激振器、一个质体、弹簧A与弹簧B;其中质体通过弹簧A与弹簧B与地基相连,弹簧对称分布于质体上;四个激振器两两分别安装在质体的上侧和下侧,每个激振器中各有一偏心转子,偏心转子由各自的感应电动机驱动,分别绕着旋转轴线中心旋转,同侧的两个频率相同的激振器关于y轴对称分布,但旋转方向相反,倍频自同步驱动,实现设备双频双直线轨迹运动功能。
2.权利要求1所述的单质体四机倍频自同步驱动振动机的参数确定方法,其特征在于,所述的四个激振器的参数确定方法,包括如下步骤:
步骤1,建立动力学模型和系统运动微分方程;
建立坐标系:四个激振器分别绕着旋转中心轴o1,o2,o3和o4旋转;分别是四个转子的旋转角;四个激振器与x轴的夹角分别用β1,β2,β3和β4表示;整个系统展现出三个自由度:直线运动x,y和摆动角ψ;
根据Lagrange方程,得系统的运动微分方程如下:
其中
式中,
M——系统总质量;
m——质体质量;
mj——激振器j偏心质量,m1=m2,m3=m4;
Ji——激振器i转动惯量,J1=J2,J3=J4;
Jm——质体自身转动惯量;
g——重力加速度;
lj——激振器j回转轴心oj到质体中心O的距离,l1=l2,l3=l4;
rj——激振器j偏心距,r1=r2,r3=r4;
fi——电机i轴阻尼系数,i=1,2,3,4;
le——系统绕质心当量回转半径;
Tei——电机i电磁输出转矩,i=1,2,3,4;
βj——激振器j回转轴心oj到质体中心O的连线与x轴夹角,β1+β2=-π,β3+β4=π;
kx,ky,kψ——x,y和ψ方向弹簧刚度;
fx,fy,fψ——x,y和ψ方向阻尼系数;
步骤2,四激振器倍频同步理论分析
四个激振器实现同步运转,激振器1和2的转速相同,激振器3和4的转速是激振器1和2转速的整数倍,分别用n3、n4表示,以顺时针旋转方向为正方向,有:
其中
σ1=σ3=-1,σ2=σ4=1
式中,小参数ε是激振器偏心块质量与系统总质量的比值;σi(i=1,2,3,4)值的正负表示偏心块的旋转方向,正值代表顺时针方向旋转,负值代表逆时针方向旋转;
把旋转相位角作如下表示:
其中
τ=ωt,n1=n2=1
式中,Δi是由于系统的运动而随着激振器偏心转子的产生阶段而缓慢变化的函数;
将式(4)代入式(3)中得:
其中,
ψij +=(σini+σjnj)τ+σiΔi+σjΔj+βi+βj,ψij -=(σini-σjnj)τ+σiΔi-σjΔj+βi-βj
式(5)为激振器实现同步的基本表达式;
把式(5)写成标准形式:
式(5)和(6)关于未知参数Δi和νi建立的一阶微分方程表达式如下:
Δi=Δi,i=1,2,3,4
其中
σini+σjnj≠0时pij=1/(σini+σjnj),σini+σjnj=0时pij=0
σini-σjnj≠0时qij=1/(σini-σjnj),σini-σjnj=0时qij=0
同样改进第二近似解:
Δi=Δi,i=1,2,3,4
将式(9)代入式(7)等号的左边,Ωi和Δi作为固定值并取关于τ=0~2π上的平均值;考虑到相同转速的激振器会互相反向旋转,得到如下关系式:
其中
σini+σjnj=0时us=1,ψij (1)=σiΔi+σjΔj+βi+βj否则us=0
σini+σrnr=0时uh=1,ψij (2)=σiΔi+σrΔr+βi+βr否则uh=0
σini-2σjnj=0时ul=1,γij (1)=σiΔi-2σjΔj+βi-βj否则ul=0
σini+2σrnr=0时um=1,γij (2)=σiΔi+2σjΔj+βi+βj否则um=0
σini-2σjnj+σrnr=0时ud1=1,ηijr (1)=σiΔi-2σjΔj+σrΔr+βi-2βj+βr否则ud1=0
σini-2σjnj-σrnr=0时ud2=1,ηijr (2)=σiΔi-2σjΔj-σrΔr+βi-2βj-βr否则ud2=0
σini+2σjnj+σrnr=0时ud3=1,ηijr (3)=σiΔi+2σjΔj+σrΔr+βi+2βj+βr否则ud3=0
σini+2σjnj-σrnr=0时ud4=1,ηijr (4)=σiΔi+2σjΔj-σrΔr+βi+2βj-βr否则ud4=0
a12=a21=1,a34=a43,a13=a14=a23=a24,α1=α2,α3=α4,
k1=k2,k3=k4,l1=l2,l3=l4,
A11=A12=A21=A22=A1,A33=A34=A43=A44=A2,A13=A14=A24=A23
步骤三,推导四激振器同步性及稳定性条件
(1)当转速相等时,取到式(10)的ε项为止,得:
当系统处于稳定状态时,式(11)中参数的表达式为:
激振器实现相同频率,转速比为1:1的同步性条件为:
Ωi0=0,i=1,2,3,4
在式(14)中,考虑到稳定状态,二倍频同步条件为:
假定初始相位Δi0和Ωi0都具有小偏差,做如下设定:
Δi=Δi0+δi,Ωi=Ωi0+ξi,i=1,2,3,4 (16)
将式(16)代入式(11)中得到系统微分方程表达式为:
其中
整理式(17)得到关于δi(i=1,2,3,4)的表达式为:
取特征值为λ,得到的特征方程如下:
根据Routh-Hurwitz准则分析并整理得到以下稳定性判据:
式中,ε,α1 (1),α3 (1),a31,a34,k1,A1,A2均大于0,且ε值很小,ε2K21和ε2K23无限接近于0,所以认为εA1+ε2K21,4εA2+ε2K23,1.5ε2a31k1均为正值,经过分析确定:
cos(Δ20-Δ10)>0(21)
(3)当n3=n4=3时,激振器3和4的稳定转速是激振器1和2的三倍,系统实现三倍频同步;同步条件表达式(13)改为:
考虑到稳定状态,转速比为1:3的同步性条件为:
为了寻求稳定相位角,求出系统在稳定状态下的微分方程表达式,引出它的特征方程为:
其中
根据Routh-Hurwitz准则分析得到稳定性判据如下:
整理求解式(25)得:
cos(Δ20-Δ10)>0(26)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011353808.XA CN112620101B (zh) | 2020-11-27 | 2020-11-27 | 单质体四机倍频自同步驱动振动机及其参数确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011353808.XA CN112620101B (zh) | 2020-11-27 | 2020-11-27 | 单质体四机倍频自同步驱动振动机及其参数确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112620101A CN112620101A (zh) | 2021-04-09 |
CN112620101B true CN112620101B (zh) | 2021-11-02 |
Family
ID=75306419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011353808.XA Active CN112620101B (zh) | 2020-11-27 | 2020-11-27 | 单质体四机倍频自同步驱动振动机及其参数确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112620101B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113158365B (zh) * | 2021-04-12 | 2023-10-03 | 东北大学 | 一种单机驱动摆锤类直线运动振动机及其参数确定方法 |
CN113111457B (zh) * | 2021-04-12 | 2023-10-27 | 东北大学 | 一种双机驱动单摆类平面运动振动破碎机及参数确定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2465821Y (zh) * | 2001-01-03 | 2001-12-19 | 西南石油学院 | 变椭圆变直线振型振动筛 |
CN101662248A (zh) * | 2009-09-28 | 2010-03-03 | 东北大学 | 空间三向自同步振动筛及其结构参数的确定方法 |
CN109499696A (zh) * | 2018-12-17 | 2019-03-22 | 东北大学 | 一种多机驱动高频振动磨机的参数确定方法 |
CN109794329A (zh) * | 2018-12-17 | 2019-05-24 | 东北大学 | 一种双质体四机驱动圆周运动高频振动磨机的参数确定方法 |
WO2019200426A1 (en) * | 2018-04-20 | 2019-10-24 | S and V Brown Holdings Pty Ltd | Mountings for vibrating machines and methods of isolating vibrations |
CN210788090U (zh) * | 2019-08-31 | 2020-06-19 | 河南威猛振动设备股份有限公司 | 一种十字形振动筛 |
-
2020
- 2020-11-27 CN CN202011353808.XA patent/CN112620101B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2465821Y (zh) * | 2001-01-03 | 2001-12-19 | 西南石油学院 | 变椭圆变直线振型振动筛 |
CN101662248A (zh) * | 2009-09-28 | 2010-03-03 | 东北大学 | 空间三向自同步振动筛及其结构参数的确定方法 |
WO2019200426A1 (en) * | 2018-04-20 | 2019-10-24 | S and V Brown Holdings Pty Ltd | Mountings for vibrating machines and methods of isolating vibrations |
CN109499696A (zh) * | 2018-12-17 | 2019-03-22 | 东北大学 | 一种多机驱动高频振动磨机的参数确定方法 |
CN109794329A (zh) * | 2018-12-17 | 2019-05-24 | 东北大学 | 一种双质体四机驱动圆周运动高频振动磨机的参数确定方法 |
CN210788090U (zh) * | 2019-08-31 | 2020-06-19 | 河南威猛振动设备股份有限公司 | 一种十字形振动筛 |
Also Published As
Publication number | Publication date |
---|---|
CN112620101A (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112620101B (zh) | 单质体四机倍频自同步驱动振动机及其参数确定方法 | |
CN101997357B (zh) | 单电机单向非周期振动产生方法及单向非周期振动装置 | |
CN109794329B (zh) | 一种双质体四机驱动圆周运动高频振动磨机的参数确定方法 | |
CN107252780A (zh) | 一种三机驱动双质体亚共振自同步概率等厚筛及参数确定方法 | |
CN109499697B (zh) | 一种三磨筒高频振动磨机的参数确定方法 | |
CN106503307B (zh) | 双激振电机驱动的振动系统的同步特性分析方法及装置 | |
CN101690929A (zh) | 一种四机驱动自同步振动筛及结构参数确定方法 | |
CN109499696B (zh) | 一种多机驱动高频振动磨机的参数确定方法 | |
CN109499695B (zh) | 一种多机驱动大型高频振动磨机的参数确定方法 | |
CN112620102B (zh) | 双机双频自同步驱动双圆周运动轨迹振动筛及其参数确定方法 | |
CN109635453B (zh) | 一种双机驱动高频振动球磨机的设计参数确定方法 | |
CN112604955B (zh) | 三机倍频自同步驱动可变轨迹振动筛及参数确定方法 | |
CN103595177A (zh) | 单向非周期振动装置及非周期振动控制方法 | |
CN107187810A (zh) | 一种四机驱动自同步振动给料机及其结构参数确定方法 | |
CN109701697B (zh) | 一种四机驱动双质体振动冲击破碎机及其参数确定方法 | |
CN109614725B (zh) | 一种紧凑型振动干燥/冷却流化床的参数确定方法 | |
CN113032924B (zh) | 一种单质体复合同步驱动振动机及其参数确定方法 | |
CN109649965B (zh) | 一种亚共振双质体四机驱动振动筛分输送机的参数确定方法 | |
CN106441359A (zh) | 一种高精度的微机械陀螺现场模态匹配方法 | |
CN112604954B (zh) | 双质体四机倍频自同步驱动振动机及其参数确定方法 | |
Li et al. | Multi-frequency vibration synchronization and stability of the nonlinear screening system | |
CN113158365B (zh) | 一种单机驱动摆锤类直线运动振动机及其参数确定方法 | |
CN117181581A (zh) | 一种内驱式自同步振动机及参数确定方法 | |
Zhang et al. | Comments on the stability of the synchronous states of three vibrators in a vibrating system with two rigid frames | |
CN113111460B (zh) | 一种摆驱动自具清筛孔功能的振动筛及其参数确定方法 |
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 |