CN109974832B - 一种高速摇振系统振幅的算法 - Google Patents
一种高速摇振系统振幅的算法 Download PDFInfo
- Publication number
- CN109974832B CN109974832B CN201910266563.8A CN201910266563A CN109974832B CN 109974832 B CN109974832 B CN 109974832B CN 201910266563 A CN201910266563 A CN 201910266563A CN 109974832 B CN109974832 B CN 109974832B
- Authority
- CN
- China
- Prior art keywords
- displacement
- shaking
- amplitude
- sampling points
- maximum
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/12—Measuring characteristics of vibrations in solids by using direct conduction to the detector of longitudinal or not specified vibrations
- G01H1/16—Amplitude
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种高速摇振系统振幅的算法,包括以下步骤:第一步:在高速摇振系统内设置位移传感器;第二步:所述第一步位移传感器测定一个振动周期内至少4个取样点位移;所述相邻两取样点对应的时间差相同;第三步:根据第二步所得取样点,摇振位移取样设置的最大硬件扫描周期计算公式为:Tscanmax=Tshake/4第四步:在一个振动周期内测定数据读取受不确定因素影响大,采用队列法测定多个摇振周期的数值,在这多个周期的数据中选出最大值和最小值,计算出振幅;第五步:根据第四步所得的最大位移Smax和最小位移Smin,通过Smax‑Smin计算得出振幅。可以设置合理的振幅编码器取样点的采样周期;在确保CPU正常使用负荷,不影响其他软件使用的同时,测定数据更精准。
Description
技术领域
本发明涉及高速摇振系统领域,具体为一种高速摇振系统振幅的算法。
背景技术
目前,在车速250m/min以下的造纸机上,网案摇振几乎是必须配备的装置,在这样的低车速下,网案摇振对成形匀度的影响是显著的。随着自动化技术发展,在造纸行业,低速纸机已经逐渐被取缔和淘汰,取而代之的中高速纸机同样有改善纸页匀度的需求。而对于中高速纸机,早期的网案摇振系统对于改变纸张特性已经无能为力了。
对于中高速纸机,高速摇振装置能具有更高的频率和更宽的振幅,并且只对胸辊进行摇振,对高速纸机的网案设计和脱水元件的排列不产生任何影响,能有效改善纸张的物理特性。
对于摇振装置的振幅检测,一般都是安装一个振幅编码器,实时读取摇振装置的振动位移。最大位移减去最小位移,就是振幅。
但是,在振幅的计算过程中,往往会遇到以下两点问题:1、应该如何设置合理的振幅编码器取样点的采样周期。设置过大,无法采样到振动位移的最大值和最小值,导致振幅计算不准确。设置过小,会增加CPU的负荷,影响其他程序的正常运行,严重的情况下会导致CPU停机。2、由于程序运行、信号传输等问题,无法保证取样点的最大值和最小值分别取值在振动位移的波峰和波谷,所以,该用什么方法去处理取样点,使取样点更加接近波峰和波谷,使振幅的计算值更加精确。这两个难点决定了摇振系统的运行稳定性。目前,市场上还没有一套行之有效的标准算法来解决这两个问题。
发明内容
本发明的目的就在于为了解决上述问题而提供一种高速摇振系统振幅的算法,包括以下步骤:
第一步:在高速摇振系统内设置位移传感器;
第二步:所述第一步位移传感器测定一个振动周期内至少4个取样点位移;所述相邻两取样点对应的时间差相同;摇振位移的取样点,即为硬件的扫描时间;
第三步:根据第二步所得取样点,摇振位移取样设置的最大硬件扫描周期计算公式为:
Tscanmax=Tshake/4=(60/shakemax*1000)/4;
所述Tshake为一个摇振周期时间;
所述shakemax为摇振系统的最大摇振频率;
所述Tscanmax为最大硬件扫描时间;
第四步:在一个振动周期内测定数据读取受不确定因素以及误差影响大,因此采用队列法测定多个摇振周期的数值,并且在这多个周期的数据中选出最大值和最小值,计算出振幅;
N=C*Tshake/Tscan
=C*(60/shakemax*1000)/[(60/shakemax*1000*4n)]
=C*4n
其中Tscan为硬件扫描周期;
C是一个常数,为整数个摇振周期;
n为大于等于1的整数;
N为取样点数;
根据队列法运算,从中选出最大位移Smax和最小位移Smin;
第五步:根据第四步所得的最大位移Smax和最小位移Smin,通过公式PV_STROKE=Smax-Smin计算得出振幅。
作为上述技术方案的进一步改进:
所述第四步不确定因素包括传感器读取数据所需要时间、程序扫描所需要时间和阻力。
所述摇振位移为PLC通过模拟量通道直接从传感器读取,或者为PLC通过PROFIBUS网络从传感器读取。
所述第二步取样点以时间为轴相对应成倍增加,取样点越接近波峰值和波谷值。
有益效果:
1、可以设置合理的振幅编码器取样点的采样周期。
2、在确保CPU正常使用负荷,不影响其他软件使用的同时,测定数据更精准。
3、更加快捷方便准确地测定振幅。
4、由于在摇振系统的控制中,频率在一定的车速范围内都是保持不变的,所以,摇振系统的控制关键就在于振幅的调节。振幅实际值的计算精确程度直接决定了摇振系统的换相电机是否能准确的进行换相工作。所以,振幅的精确计算,对于摇振系统的稳定运行有着重要的意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是摇振最少取样点理想位移示意图。
图2是摇振最少取样点实际位移示意图。
图3是摇振合理取样点理想位移示意图。
图4是摇振合理取样点实际位移示意图。
图5是队列法计算示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行详细的描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式,都属于本发明所保护的范围。
如图1-5所示,本实施例的高速摇振系统振幅的算法,包括以下步骤:
第一步:在高速摇振系统内设置位移传感器;
第二步:第一步位移传感器测定一个振动周期内至少4个取样点位移;相邻两取样点对应的时间差相同;摇振位移的取样点,即为硬件的扫描时间;
第三步:根据第二步所得取样点,摇振位移取样设置的最大硬件扫描周期计算公式为:
Tscanmax=Tshake/4=(60/shakemax*1000)/4;
Tshake为一个摇振周期时间;
shakemax为摇振系统的最大摇振频率;
Tscanmax为最大硬件扫描时间;
第四步:在一个振动周期内测定数据读取受不确定因素以及误差影响大,因此采用队列法测定多个摇振周期的数值,并且在这多个周期的数据中选出最大值和最小值,计算出振幅;
N=C*Tshake/Tscan
=C*(60/shakemax*1000)/[(60/shakemax*1000*4n)]
=C*4n
其中Tscan为硬件扫描周期;
C是一个常数,为整数个摇振周期;
n为大于等于1的整数;
N为取样点数;
根据队列法运算,从中选出最大位移Smax和最小位移Smin;
第五步:根据第四步所得的最大位移Smax和最小位移Smin,通过公式PV_STROKE=Smax-Smin计算得出振幅。
第四步不确定因素包括传感器读取数据所需要时间、程序扫描所需要时间和阻力。
摇振位移为PLC通过模拟量通道直接从传感器读取,或者为PLC通过PROFIBUS网络从传感器读取。
第二步取样点以时间为轴相对应成倍增加,取样点越接近波峰值和波谷值。
实施例一、
在高速摇振系统上,一般是由位移传感器去检测左右摆动的位移,再由最大位移减去最小位移,得出振幅。在一个振动周期内,如果要取样到最大值和最小值,那么最少要取样4个点,如图1所示,t1-t4分别为4个振动位移取样点,S2和S4分别为一个摇振周期内位移的最大值和最小值,则实际振幅PV_STROKE=Smax-Smin=S2-S4。
图1是摇振最少取样点理想位移图,但是由于传感器读取数据需要时间,程序扫描需要时间,以及阻力、生产等一些其他不确定因素,所以,根本没法保证在一个摇振周期内,程序读取到的位移最大值,就是波峰值S2,同理,也无法保证程序读取到的位移最小值,就是波谷值S4。
实际上,真正读取到的数据应该是如图2所示。
从图2可以看出,如果一个振动周期内只取样4次,得到的位移最大值和最小值,有可能是远离波峰值和波谷值的。
所以,为了使计算更加精确,在波峰两边各增加一个取样点,波谷两边各增加一个取样点,那么一个振动周期内就是取样8次,其理想位移图如图3所示。
再结合实际情况的取样点位移图,如图4所示。
从图4可以看出,一个振动周期内至少取样8次,才可以保证取样到的位移最大值和最小值,是比较接近波峰值和波谷值的。
在实际检测过程中,摇振位移是由PLC通过模拟量通道直接从传感器读取,或者是PLC通过PROFIBUS等网络从传感器读取。所以,不管通过何种方式,振动位移的取样点,就是硬件的扫描时间。
所以,对于高速摇振系统,振动位移取样可以设置的最大硬件扫描周期计算公式为:
Tscanmax=Tshake/8=(60/shakemax*1000)/8,单位:ms。
式中,Tshake指的是一个摇振周期时间,单位是ms。
shakemax为摇振系统的最大摇振频率,单位是shk/min。
Tscanmax为最大硬件扫描时间,单位是ms。
进一步考虑,如果想要取样点的最大值和最小值,更加接近波峰值和波谷值,可以在波峰两边各增加两个取样点,波谷两边各增加两个取样点,即一共取样12个点,以此类推,可以得到最终的硬件扫描周期计算公式:
Tscan=(60/shakemax*1000)/(8+4n')
式中,Tscan为硬件扫描周期,单位是ms。
n'为自然数。
n'越大,则取样点越接近波峰和波谷,但是同时也要兼顾CPU的负荷能力。
实际振幅的进一步精确计算
因为传感器读取数据需要时间,程序扫描需要时间,以及阻力、生产等其他不确定因素,所以每个周期读取的数据都是不同的。甚至,每个周期的波峰和波谷值,都不一定是一样的。
那么,为了尽可能使读取的最大值接近波峰,最小值接近波谷,需要怎么做呢?
可以利用队列法,如图5。需要多取几个摇振周期的数值,并且在这几个周期的数据中选出最大值和最小值,计算出振幅。
根据图5可以看出,
Smax=max(Sn,Sn+1…Sn+39)
Smin=min(Sn,Sn+1…Sn+39)
然后在下个取样时间,将取样点Sn删除,读取新的取样点Sn+40。则会得出一个新的摇振位移最大值和最小值,分别为
Smax=max(Sn+1,Sn+2…Sn+40)
Smin=min(Sn+1,Sn+2…Sn+40)
以此类推,根据先进先出的原则,加入每个新的取样点,都可以计算出一个新的Smax和Smin。再根据公式PV_STROKE=Smax-Smin计算出振幅。
根据多次试验,取4-5个振动周期的取样点作为计算对象。进行队列法运算的取样点个数的计算公式:
N=C*Tshake/Tscan=C*(60/shakemax*1000)/[(60/shakemax*1000)/(8+4n')]=C*(8+4n')
式中,shakemax为摇振系统的最大摇振频率,单位是shk/min。
Tshake为一个振动周期,单位是ms。
Tscan为硬件扫描周期,单位是ms。
C是一个常数,表示取多少个摇振周期。该值取值太小无法达到精确计算的效果,取值太大会增加数据的滞后性。一般C=4或者C=5。
n'为自然数。
根据这个公式,就可以计算出,选取N个取样点,进行队列法运算,从中选出最大位移Smax和最小位移Smin。
由于在摇振系统的控制中,频率在一定的车速范围内都是保持不变的,所以,摇振系统的控制关键就在于振幅的调节。振幅实际值的计算精确程度直接决定了摇振系统的换相电机是否能准确的进行换相工作。所以,振幅的精确计算,对于摇振系统的稳定运行有着重要的意义。
本发明中各实施例的技术方案可进行组合,实施例中的技术特征亦可进行组合形成新的技术方案。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (4)
1.一种高速摇振系统振幅的算法,其特征在于,包括以下步骤:
第一步:在高速摇振系统内设置位移传感器;
第二步:所述第一步位移传感器测定一个振动周期内至少4个取样点位移;所述相邻两取样点对应的时间差相同;摇振位移的取样点,即为硬件的扫描时间;
第三步:根据第二步所得取样点,摇振位移取样设置的最大硬件扫描周期计算公式为:
Tscanmax=Tshake/4=(60/shakemax*1000)/4;
所述Tshake为一个摇振周期时间;
所述shakemax为摇振系统的最大摇振频率;
所述Tscanmax为最大硬件扫描时间;
第四步:在一个振动周期内测定数据读取受不确定因素以及误差影响大,因此采用队列法测定多个摇振周期的数值,并且在这多个周期的数据中选出最大值和最小值,计算出振幅;
N=C*Tshake/Tscan
=C*(60/shakemax*1000)/[(60/shakemax*1000*4n)]
=C*4n
其中Tscan为硬件扫描周期;
C是一个常数,为整数个摇振周期;
n为大于等于1的整数;
N为取样点数;
根据队列法运算,从中选出最大位移Smax和最小位移Smin;
第五步:根据第四步所得的最大位移Smax和最小位移Smin,通过公式PV_STROKE=Smax-Smin计算得出振幅。
2.根据权利要求1所述的高速摇振系统振幅的算法,其特征在于:所述第四步不确定因素包括传感器读取数据所需要时间、程序扫描所需要时间和阻力。
3.根据权利要求1所述的高速摇振系统振幅的算法,其特征在于:所述摇振位移为PLC通过模拟量通道直接从传感器读取,或者为PLC通过PROFIBUS网络从传感器读取。
4.根据权利要求1所述的高速摇振系统振幅的算法,其特征在于:所述第二步取样点以时间为轴相对应成倍增加,取样点越接近波峰值和波谷值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910266563.8A CN109974832B (zh) | 2019-04-03 | 2019-04-03 | 一种高速摇振系统振幅的算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910266563.8A CN109974832B (zh) | 2019-04-03 | 2019-04-03 | 一种高速摇振系统振幅的算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109974832A CN109974832A (zh) | 2019-07-05 |
CN109974832B true CN109974832B (zh) | 2021-04-02 |
Family
ID=67082742
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910266563.8A Active CN109974832B (zh) | 2019-04-03 | 2019-04-03 | 一种高速摇振系统振幅的算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109974832B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111649978B (zh) * | 2020-06-05 | 2022-04-29 | 吉林大学 | 基于形状记忆合金驱动的土壤采样器 |
CN114280921B (zh) * | 2021-12-15 | 2023-11-14 | 浙江华章科技有限公司 | 一种摇振装置参数优化方法、装置及存储介质 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4318617A (en) * | 1979-12-14 | 1982-03-09 | Keuffel & Esser Company | DC Shift error correction for electro-optical measuring system |
JPS6123559A (ja) * | 1984-07-12 | 1986-02-01 | Nippon Kokan Kk <Nkk> | 鋼の連続鋳造用鋳型の振動方法 |
JP3183237B2 (ja) * | 1997-11-21 | 2001-07-09 | 株式会社ケイ | 可変型振動テーブル用振動装置 |
DE10124760A1 (de) * | 2001-05-21 | 2003-02-20 | Siemens Ag | Verfahren zur kontaktlosen, linearen Positionsmessung |
JP4209709B2 (ja) * | 2003-03-20 | 2009-01-14 | 株式会社キーエンス | 変位計 |
CN1995936A (zh) * | 2005-12-31 | 2007-07-11 | 中国科学院半导体研究所 | 波长可调谐激光器的自动化测试控制系统及方法 |
JP5421669B2 (ja) * | 2009-06-23 | 2014-02-19 | アズビル株式会社 | 振幅算出装置および振幅算出方法 |
CN203768726U (zh) * | 2013-11-22 | 2014-08-13 | 广东理文造纸有限公司 | 一种造纸网部摇振系统 |
CN108411678B (zh) * | 2018-04-24 | 2023-06-16 | 华南理工大学 | 一种振幅可调的高速摇振箱装置及振幅调节方法 |
-
2019
- 2019-04-03 CN CN201910266563.8A patent/CN109974832B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109974832A (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109974832B (zh) | 一种高速摇振系统振幅的算法 | |
CN108776460B (zh) | 基于pc的高速图像检测分拣生产线控制系统及其控制方法 | |
CN101750046A (zh) | 一种角度测量装置和方法及工程机械 | |
CN110879298B (zh) | 基于通讯式编码器的速度获取方法 | |
CN103675324A (zh) | 一种具有防抖功能的转速测量方法 | |
CN101893457B (zh) | 连续运动高精度全参数检测方法 | |
CN114229385A (zh) | 基于双转速传感器的火电厂给煤机转速输出方法及装置 | |
CN110273260B (zh) | 一种基于无传感滚筒洗衣机偏心和称重的采集方法 | |
CN109764897A (zh) | 一种正余弦编码器高速信号采集及细分方法和系统 | |
CN101478285B (zh) | 电机速度检测方法 | |
CN111630772B (zh) | 伺服系统机械参数辨识方法、伺服控制系统和存储装置 | |
CN113091797A (zh) | 一种监测编码器故障状态的方法及装置 | |
CN105763170B (zh) | 一种电力信号数字滤波方法 | |
CN103439660A (zh) | 一种用于测试汽车仪表步进电机转动的方法 | |
CN115453140A (zh) | 一种增量式编码器的检测方法及检测装置 | |
CN201561749U (zh) | 一种角度测量装置及工程机械 | |
CN215727577U (zh) | 一种颗粒计数系统 | |
CN104444671B (zh) | 电梯测速方法、系统以及电梯低速启动控制方法 | |
CN104678301B (zh) | 一种基于脉冲反馈检测电机运动学参数的方法 | |
CN112560762A (zh) | 车体回转角度数据处理方法、装置、控制器以及介质 | |
CN112301641A (zh) | 一种带mems传感器的洗衣机转速测量方法及位移数据检测方法 | |
US6310458B1 (en) | Blended velocity estimation | |
CN113722376B (zh) | 一种用于船用发电机组瞬态转速采集的数据处理方法 | |
CN117990170B (zh) | 用于质量流量计的振动传感器驱动控制方法、系统及设备 | |
CN1458500A (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 |