本发明属于大型设备筛分输送技术领域,主要应用振动同步理论,利用振动系统不同共振类型的振动状态对物料进行筛分输送的大型设备的开发研制,并最终实现其工程应用价值。
背景技术:
在振动筛分输送领域,很多设备已经应用到工程实际中,而本专利提出一种新的振动筛分输送设备模型。以双质体四机驱动动力学模型为研究对象,应用平均法和哈密顿最小作用量的原理,分别得到四激振器最终实现同步的同步性判据,分析了系统实现同步的耦合机理,定义了同步性和稳定性能力系数,数值方面,给出了两质体相对运动的幅-频曲线关系、系统同步能力系数曲线、无量纲耦合力矩最大值及稳定性系数,界定出系统不同共振区间下的三类相位关系:激振器间、质体间以及质体与激振器间的相位关系。而三类相位关系就是机械设备最终功能的体现。仿真方面,验证了数值结论的正确性。可根据双质体四机驱动振动同步理论为新型、长距离、大输送量的振动筛分输送设备的研制开发提供理论指导。
技术实现要素:
发明目的:针对目前振动筛分输送设备的短距离、小型化、耗能等弊端,本发明提出了亚共振双质体四机驱动振动筛的设计方法,理论上论述了该动力学模型的同步性条件及同步状态下的稳定性判据,并通过仿真分析验证数值分析的正确性,最终确定了亚共振区为系统的合理工作点,以实现两激振器正向叠加,进而,为新型、长距离、大输送量的振动输送机的研制开发提供理论指导。
本发明是通过以下技术方案实现的,将该振动筛分输送机转化成动力学模型为:包括两个质体、四个激振器及两组弹簧,上下两部分分别称为主振系统和隔振系统;
主振系统上分别刚性连接四个激振器,分别记为左右两组,每组激振器均沿着弹簧对称分布于质体m1上下两侧,相互对称的一组激振器转向相反;质体m1通过一组弹簧k1将其与质体m2连接,并在两个质体之间垂直于弹簧伸缩的方向上设置导向板,使系统只有沿x方向一个自由度;
隔振系统质体m2上方连接于质体m1,下方通过弹簧k2和导向板固连于基础上;所述的激振器和振动系统的参数确定方法,包括如下步骤:
步骤1,动力模型的建立
振动筛分输送机动力学模型见图1,分别以质体1和质体2的质心为原点建立两个坐标系o1-x1、o2-x2。
在广义坐标系o-x,根据系统的动能t,势能v以及能量函数,得出两质体在x方向上的运动微分方程如下
其中
m1=m1+4m0,m2=m2,joi=m0iri2
式中,m1——质体1质量;
m2——质体2质量;
m0i——激振器i质量(i=1~4);
r——激振器偏心距;
k1,k2——x方向上弹簧刚度;
f1,f2——x方向上阻尼系数。
joi——激振器i的转动惯量(i=1~4);
令系统中的四个激振器的所有参数保持一致,取m01=m02=m03=m04=m0。假定四个激振器的平均相位为
设定四个偏心转子在稳态时同步角速度为ωm0,当系统稳态运行时,可得到位移与加速度的关系,如下,
当系统处于稳态运行时,其角加速度几乎为0,即
对于式(6)和(7),由于k2<<k1,有
做如下处理,
其中,
式中,m为诱导质量。
由式(8),可得两质体在x方向相对运动固有频率ω0(也称其为主振系统固有频)为
式(8)的响应为
其中
工程上,对于小阻尼振动机,f1=2ξ′xmω0,其中ξ′x为等效相对阻尼系数,ξ′x≤0.07。由式(10)可知,当z0=1(即ωm0=ω0)时,a12在共振点处取最大值,说明ω0对应的是两质体x方向反相位相对运动固有频率。通过对式(10)进一步化简,可以得到反相位相对运动响应幅值λ12,表示为λ12=|a12·s|(11)
其中
s=(η12+η22+η32+η42+2η1η2cos(2α1)+2η2η3cos(2α2)+2η3η4cos(2α3)+2η1η3cos(2α1+2α2)+2η1η4cos(2α1+2α2+2α3)+2η2η4cos(2α2+2α3))1/2
反相位相对运动响应幅值λ12在工程中具有很大应用价值。
步骤2,系统同步性分析
使用传递函数法,对式(1)与(2)进行求解,令x1=x1(s),
其中
式中,γ1——质体1在x方向上的滞后角;
γ2——质体2在x方向上的滞后角。
根据公式(1),可得两质体在x方向上质量矩阵、刚度矩阵及特征方程如下
式中,m′为质量耦合矩阵,k′为刚度耦合矩阵,δ(ω2)为特征方程。
令δ(ω2)=0,可以得到
一般情况下,隔振系统的刚度k2远远小于主振系统的刚度k1,即k2<<k1,因此将ωinv中的k2忽略,则有ωinv=ω0,可以看出,ωinv对应的是两质体在x方向反相位相对运动固有频率。显然,ωsa是两质体在x方向同相位相对运动固有频率。
当四个激振器能够同步运转时,有
其中
在上述积分过程中,需注意的是固有频率ω1为质体2的固有频率,在后续的数值与仿真分析中会重点说明。对比
各激振器间输出力矩之差为
对式(16)~(18)进行整理,可得,
其中
在上述推导中,
综上,结合(19),(20)与(21),可得到四激振器的同步判据为
通过式(25),(26)与(27),可以看出相邻两个激振器的无量纲残余力矩之差的绝对值小于或等于无量纲耦合力矩的最大值。进一步可得到τc41max的表达式如下:
τc41max=-(τc12max+τc23max+τc34max)(28)
对
四激振器平均无量纲负载力矩的约束函数如下
定义同步能力系数为ζij(i,j=1,2,3,4),可得,
同步能力系数越大,系统的同步能力越强,越容易达到同步。
步骤3,推导稳定性条件
对于双质体四激振器系统,系统的动能(t)与动能(v)如下,
在一个周期内平均动能(et)与平均势能(ev)分别为
其中
在一个周期内系统的hamilton平均作用量(i)为
两种同步状态中的稳定相位差解
其中
令,
为了使i的hesse矩阵正定,即h矩阵正定,应满足
h1>0,h2>0,h3>0(39)
将h1,h2与h3定义为系统的同步条件下的稳定能力系数,式(39)即为系统的稳定性能力的表达式,当满足式(39)时,系统稳定。
本发明的有益效果:
(1)本发明在模型上进行创新,选用两个质体,其中一质体上安有四个电机,采用四机驱动,另一质体作隔振体,并且两质体之间以及隔振体与地基之间也同样通过弹簧相互连接,在模型上进行创新,更接近工程实践。
(2)本发明应用振动同步理论,采用四机驱动实现系统的同步工作。区别以往工作点选择在超远工作区的振动筛分输送类设备,本专利提出的模型将工作区域选择在亚共振区域,在该区域内,系统在相同振幅的条件下,亚共振区域内激起的同样的振幅所需的激振力是其超远共振条件下的1/5~1/3。因而,在亚共振状态工作的振动系统所需驱动电机功率会相应降低,进而可以实现能源的节约。
(3)本发明的研究内容对于工程上大型筛分输送的机械设备,即新型、长距离、大产量大型振动筛分输送设备,对其结构参数设计以及工作区域的选择具有重大指导作用。
附图说明
图1双质体四机驱动振动系统动力学模型
图中:1.激振器2、2.质体1、3.弹簧1、4.质体2、5.导板、6.弹簧2、7.激振器4、8.激振器1、9.激振器3
图中各参数含义:
o1-x1质体1坐标系;
o2-x2质体2坐标系;
m01——激振器1质量;
m02——激振器2质量;
m03——激振器3质量;
m04——激振器4质量;
m1——质体1质量;
m2——质体2质量;
k1——弹簧1在x1方向上弹簧刚度;
k2——弹簧2在x2方向上弹簧刚度。
图2反相位相对运动幅-频响应曲线。
图3系统同步性能力曲线。
图4系统稳定性能力系数;
(a)稳定性能力系数h1;
(b)稳定性能力系数h2;
(c)稳定性能力系数h3。
图5四激振器之间的相位关系。
图6系统响应及其与激振器之间的相位关系。
图7相对ω1的亚共振状态(区域i)的仿真结果;
(a)激振器1,2稳定相位差;
(b)激振器2,3稳定相位差;
(c)激振器3,4稳定相位差;
(d)质体1,2的相对位移(x1-x2);
(e)位移x1与x2。
图8相对ω0的亚共振状态(区域ii)的仿真结果;
(a)激振器1,2稳定相位差;
(b)激振器2,3稳定相位差;
(c)激振器3,4稳定相位差;
(d)质体1,2的相对位移(x1-x2);
(e)位移x1与x2。
图9相对ω0的超共振状态(区域iii)的仿真结果;
(a)激振器1,2稳定相位差;
(b)激振器2,3稳定相位差;
(c)激振器3,4稳定相位差;
(d)质体1,2的相对位移(x1-x2);
(e)位移x1与x2。
具体实施方式
实施例1:一种亚共振双质体四机驱动振动筛分输送机。其动力学模型见图1,包括1.激振器1;2.质体1;3.弹簧1;4.质体2;5.导板;6.弹簧2;7.激振器4;8.激振器2;9.激振器3。该模型由两个质体、四个激振器及两组弹簧组成,四组激振器两两反向安装在质体1上,通过弹簧1将其与质体2连接,并在两个质体之间的弹簧垂直方向上设置导板,使系统只有沿x方向一个自由度,同时,通过弹簧2将质体2连接在地基上。如图1,四个激振器回转中心分别用o1,o2,o3及o4表示,回转半径均为r。激振器1与激振器4反向回转,激振器2与激振器3反向回转。整个系统在x方向上产生位移,并且每个激振器绕自身回转轴旋转,以
双质体四机驱动振动系统的数值分析
基于四机驱动(两两反向)双质体动力学模型,对系统参数给定一些数值,进行分析,已验证理论推导的正确性。具体的系统参数如下:m1=600kg,m2=1500kg,m0=10kg,r=0.15m,k1=8000kn/m,k2=100kn/m。选取四个电机型号一致,三相鼠笼式(型号vb-1082-w,380v,50hz,6-极,δ-连接,0.75kw,转速980r/min,39kg)。
结合前面内容的分析,可以看出,对于四机驱动(两两反向)双质体系统,共存在三组重要的固有频率,质体2的固有频率ω1,反相位相对运动固有频率ωinv以及相对运动固有频率ω0(也称其为主振系统固有频)。当忽略隔振系统的刚度k2时,ωinv=ω0。将固有频率ω1与ω0作为激振频率的划分值,将整个系统的频率区间分成三部分,当ωm0<ω1时,该区域为相对ω1的亚共振区域,记作区域i;当ω1≤ωm0≤ω0时,该区域为相对ω1的超共振区域和相对ω0的亚共振区域,记作区域ii;当ωm0>ω0时,该区域为相对ω0的超共振区域,记作区域iii。
下面,对系统的稳态幅-频特性、同步性能力、系统最大耦合力矩、同步状态的稳定性以及稳态时系统相位关系在数值上加以分析。
(a)稳态幅-频特性
图2为系统相对运动幅-频响应曲线,根据系统的固有特性,将质体相对运动幅-频响应曲线分为三部分。
可以看出,在区域i,激振频率ωm0由0rad/s开始增加,两质体的相对位移为0mm。当激振频率逐渐增加,逐渐接近ω1,即a点,频率为73.03rad/s。系统进入区域ii,响应幅值不断增大,在b点处达到最大的相对幅值。此处,对应系统的另一固有频率ω0两质体相对运动的固有频率,大小为133.48rad/s。随着ωm0继续增加,系统进入区域iii,响应幅值开始下降,并在c点处出现了断点,随后,相位位移又回到0mm。由此,分析出了在不同共振区域,两质体相对运动振幅的变化趋势,可与仿真形成对比,同时还能够清晰地看出双质体的相对运动状态。
经过上述分析可以看出,区域i与区域iii的相对振幅基本为零,而响应曲线中的区域ii,即相对ω1的超共振区域和相对ω0的亚共振区域,为有用区域,在该段区域内系统具有较大并且稳定的幅值,可应用于工程实际中。
(b)同步性能力
同步性能力系数用ζ表示,是衡量每个电机达到同步的指标,同步性能力系数越大,系统的同步性能力越好。由于四个电机的型号一致,因此,四电机的同步性能力系数大小一致,如图3,四条曲线完全重合。同样,系统的两个固有频率ω1和ω0将同步性能力曲线分为三部分。在区域i中,初始阶段,当频率增高,同步性能力系数逐渐增大,随后,同步性能力系数随着激振频率ωm0的增加而减小,在点a处达到最小,此处为共振点ω1。随着激振频率继续增加,同步性能力曲线进入区域ii,在该区域内,同步性能力系数先增加,后随着系统接近另一共振点b,同步性能力系数又降低至最小,b点为共振点ω0。通过对激振器间的相位差的数值分析可知,此段中稳定相位为0。之后,同步性能力曲线进入区域iii,并且同步性能力系数随着激振频率的增加而不断增大。
其中,同步性能力曲线中的a点和b点,a点为ω1共振点,频率为73.03rad/s,b点为ω0共振点,频率为133.48rad/s。当激振频率处于a,b点时,同步性能力系数最小。可见,当系统处于共振点时,系统的同步性能力最小。
(c)系统最大耦合力矩
根据式(19)~(21)四激振器之间的无量纲耦合力矩,将系统参数带入,得到四激振器间的最大耦合力矩,可知激振器1,2与激振器2,3与激振器3,4以及激振器4,1间的最大耦合力矩均相等。
随着激振频率ωm0的增加,无论是在区域i,ii或者iii,稳定后,四激振器之间的无量纲最大耦合力矩均在在5.2附近,除了a与b点。同样,如同同步性能力曲线,a点为ω1共振点,频率为73.03rad/s,b点为ω0共振点,频率为133.48rad/s。在a,b两点,即共振点处,四激振器之间的无量纲最大耦合力矩均产生了下降。
(d)同步状态的稳定性
将系统参数带入式(38)的h1,h2与h3中,得到系统的稳定性能力系数,如图4所示,在区域i中,图(a),(b),(c)的稳定能力曲线具有相同的变化趋势,随着激振频率的增加,稳定性能力系数均为0。当激振频率逐渐经过固有频率ω1,稳定性能力曲线进入区域ii,稳定性能力系数均产生明显的增大,同时,在该区域具有稳定的相位差。随着频率的继续增加,进入区域iii,稳定性能力系数h1一直大于0,而h2则出现了多个稳定的值,h3则又回到0。
该结果与图5的激振器之间的相位差关系相互对应,可以说明系统出现了多样性的现象,该现象将在图5中详细描述。通过数值分析结果可以得出存在非线性系统多样性情况的条件有2个:其一是系统在同步条件下的稳定性系数为0;其二是存在多个相位差的稳定解。
由于非线性系统的多样性,在区域i与iii中,稳定性系数为0或者出现多个稳定值。而在区域ii,即相对固有频率ω0的亚共振区域,系统具有较强的稳定性。因此,该区域具有将强的应用价值。
(e)稳态时系统相位关系
系统的稳态相位关系主要包括三类相位关系,分别是激振器之间的相位差,系统响应与激振器间的相位关系以及系统响应之间的相位关系。
如图5所示,表示四激振器之间的相位关系,2α12表示激振器1与2的相位差,同理,2α23表示激振器2与3的相位差,2α34表示激振器3与4的相位差。从图中可以看出,根据两个共振点ω1和ω0将相位差分为三个区域,当激振频率处于区域i和区域iii时,四激振器间的相位差在[-180,180]之间存在多值,即存在多个稳定解,说明了系统的多样性。而当激振频率处于区域ii时,四激振器之间具有稳定的相位差,均为0,即2α12=2α23=2α34=0rad/s。该结论与下一章节的仿真结果相对应。由此可见,四激振器间的相位差在区域ii是稳定的,也就是说在相对ω0的亚共振区是稳定的。
如图6所示,表示系统响应与激振器间的相位关系以及系统响应之间的相位关系,γ1表示质体1响应滞后于激振器的角度,γ2表示质体2响应滞后于激振器的角度,γ12表示质体1响应滞后于质体2的角度。由下图可以看出,在区域i,质体1与质体2均滞后于激振器180°,即与激振器呈反相位。在区域ii,质体1滞后于激振器的角度逐渐减小;而质体2滞后于激振器的角度先减小后又逐渐增大。在区域iii中,质体1与激振器的滞后角度接近0度,即质体1与激振器呈同相位运动;质体2与激振器的滞后角度仍为180°,即与激振器呈反相位运动。
通过γ12可以看出系统响应之间的关系,在区域i,γ12处于0附件,说明质体1与质体2之间不存在角度滞后,呈同相位运动的状态。在区域ii,二者之间的之后角度逐渐增加,逐渐接近180°。在区域iii,二者的之后角度滞后180°,呈反相位运动的状态。
上述分析说明了ω1为质体1与质体2的同反相位运动的临界同步转速。当激振频率ωm0小于ω1时,两质体呈同相位运动,滞后角为0°;当ω1≤ωm0≤ω0时,两质体的响应开始滞后;当ωm0>ω0时,两质体呈反相位运动,滞后角为180°。此结果与下一节的仿真结果保持一致。
实施例2:双质体四机驱动振动系统的仿真分析
双质体四机驱动振动系统动力学模型,对其系统的微分方程,即式(1)~(3)利用四阶rouge-kutta程序进行仿真。在特性分析中,根据系统的固有特性,将其运行状态分为三个区域,如图2~6的区域i,ii及iii所示。同理,在仿真分析中,取三组参数,分别使系统处于三个区域内。通过仿真分析出在对应区域内系统在稳态时电机的转速、相位差以及质体1与2的位移曲线。
实际工程应用中,一般取相同的激振器,四电机的参数相同,即
η1=η2=η3=η4=1.0。系统整体参数选用如下:转子电阻rr=3.40ω,定子电阻rs=3.35ω,转子电感lr=170mh,定子电感ls=170mh,互感lm=164mh,f1=f2=0.05。振动系统的其他参数:r=0.15m,m1=600kg,m2=1500kg,m0=10kg,fx=7.6kn·s/m,ξnx=0.07。调整参数,使系统分别处于亚共振状态和超共振状态。
(a)相对ω1的亚共振状态(区域i)的仿真结果
取k1=80000kn/m,k2=100kn/m,zx=0.45,系统处于区域i,即相对于ω1的亚共振状态。下面,对系统在区域i的稳态时电机的转速、激振器间的相位差、质体1与2的位移进行分析。随着时间增加,在前20s,四个激振器达到了同步和稳定。但是在受到干扰后,系统出现了非线性多样性情况。
如图7中的曲线均是在稳态运行时四电机的同步转速稳定在980r/min~985r/min范围内绘制而成。当系统处于区域i稳态运行时,各电机的转速接近激振转速983r/min。
如图7(a),(b),(c)表示区域i中系统的稳定相位差曲线,可以看出,在未受干扰前,即前20s,激振器1与激振器2的稳定相位差2α12=-180°,激振器2与激振器3的稳定相位差为2α23=0°,激振器3与激振器4的稳定相位差均为0°,即2α34=0°。在20s与35s处分别给电机2一个干扰,可见,系统相位出现了变化。该结果与数值分析的结果相对应,体现出了非线性系统的多样性。
图7(d),(e),表示在区域i中质体1,2的位移及二者相对位移曲线。前20s,质体1的位移x1稳定在-0.4mm~0.4mm,质体2的位移x2稳定在-0.7mm~0.7mm,则二者的相对位移稳定在-0.3mm~0.3mm。分别于20s和35s处给系统干扰后,质体1与2的位移出现了变化,均变为0mm,即在给定干扰后,系统由小振幅运动变为了完全静止。
如图7(e),清晰地表示出了质体1的位移x1与质体2的位移x2之间的关系,当质体1的位移处于波峰时,质体2的位移也处于波峰,可见,质体1与质体2的位移在x方向上呈同向运动,相对位移相互抵消。此结果与数值分析结果保持一致。
(b)相对ω0的亚共振状态(区域ii)的仿真结果
取k1=8000kn/m,k2=100kn/m,zx=0.77,系统处于区域ii,相对于固有频率ω0的亚共振状态。下面,对系统在区域ii的稳态时电机的转速、激振器间的相位差、质体1与2的位移进行分析。随着时间的增加,在前20s,系统很快达到了同步和稳定,并在给干扰后,系统仍旧能够保持稳定。
如图8中的曲线均是在稳态时短时间内各电机同步转速基本稳定在750r/min~800r/min范围内时的数据绘制而成。当系统处于亚共振时,虽然激振转速为980r/min,但是系统稳态运行时,各电机的转速均小于激振转速。
如图8(a),(b),(c)表示区域ii中系统的稳定相位差曲线,可以看出在干扰前后稳定相位差不变,此时系统处于稳定状态。
如图8(d),(e)表示质体1,2的位移以及二者相对位移的曲线,前20s,质体1的位移x1稳定在-12.5mm~12.5mm,质体2的位移x2稳定在-6.5mm~6.5mm,则二者的相对位移稳定在-19mm~19mm。于20s和35s处给系统干扰后,质体1与质体2的位移很快稳定到受干扰前的状态。可见,系统具有较强的稳定性。
如图9(e),清晰地表示出了质体1的位移x1与质体2的位移x2之间的关系,当质体1的位移处于波峰时,质体2的位移正处于波谷,可见,质体1与质体2的位移在x方向上呈反向运动,相对位移相互叠加。
(c)相对ω0的超共振状态(区域iii)的仿真结果
取k1=3000kn/m,k2=100kn/m,zxo=1.27,振动系统处于区域iii,即相对于固有频率ω0的超共振状态。下面,对系统在区域iii的稳态时电机的转速、激振器间的相位差、质体1与2的位移进行分析。随着时间增加,在前20s,经过系统共振区后,四个激振器达到了同步和稳定。但是在受到干扰后,系统出现了多样性的情况。
如图9中的曲线均是在稳态时电机转速,可以得到,在运行初始阶段,各电机的转速很快达到稳定,电机1的转速稳定在975r/min-985r/min,其余电机基本稳定在983r/min。
如图9(a),(b),(c)表示在超共振状态下系统稳定的相位差,可以看出,在未受干扰前,激振器1与激振器2的稳定相位差为2α12=-180°,激振器2与激振器3的稳定相位差为2α23=0°,激振器3与激振器4的稳定相位差均为2α34=180°。即在系统达到稳定同步运行时,激振器2与3,1与4同相位,激振器1与2,3与4反相位。在20s处和35s处分别给激振器2一个干扰后,相位差出现了渐变,与区域i有同样的情况,体现出了非线性系统的多样性。
如图9(d),(e)表示质体1,2的位移以及二者相对位移的曲线,在开始阶段系统经过共振区,位移x1与x2具有较大幅值,大约3s后,质体1与质体2的位移达到稳定。质体1的位移x1稳定在-2.7mm~2.7mm,质体2的位移x2稳定在-0.25mm~0.25mm,则二者的相对位移稳定在-0.3mm~0.3mm。在20s处和35s处分别给激振器2干扰后,质体1与2的位移均变为0mm,即在给定干扰后,系统由小振幅运动变为了完全静止。
如图9(e)所示,清晰地表示了受干扰前质体1的位移x1与质体2的位移x2之间的关系,当质体1的位移处于波峰时,质体2的位移则处于波谷,可见,质体1与质体2的位移在x方向上呈反向运动,位移相互叠加。
实施例3:一款振动筛分输送机的示例数据参数。本发明并不仅限于此设计参数。
质体1质量m1=600kg,质体2质量m2=1500kg,激振器偏心块质量m0=10kg,激振器回转半径r=0.15m,质体1,2间的弹簧刚度k1=8000kn/m,质体2与地基间的弹簧刚度k2=100kn/m,电机转速ω=980r/min=102.57rad/s,此时