CN102680082A - 一种结晶器振幅测试方法和装置 - Google Patents

一种结晶器振幅测试方法和装置 Download PDF

Info

Publication number
CN102680082A
CN102680082A CN2012101645090A CN201210164509A CN102680082A CN 102680082 A CN102680082 A CN 102680082A CN 2012101645090 A CN2012101645090 A CN 2012101645090A CN 201210164509 A CN201210164509 A CN 201210164509A CN 102680082 A CN102680082 A CN 102680082A
Authority
CN
China
Prior art keywords
crystallizer
vibration
value
integrated
numerical
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.)
Pending
Application number
CN2012101645090A
Other languages
English (en)
Inventor
李凯
翟翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongye Continuous Casting Technology Engineering Co Ltd
CCTec Engineering Co Ltd
Original Assignee
Zhongye Continuous Casting Technology Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhongye Continuous Casting Technology Engineering Co Ltd filed Critical Zhongye Continuous Casting Technology Engineering Co Ltd
Priority to CN2012101645090A priority Critical patent/CN102680082A/zh
Publication of CN102680082A publication Critical patent/CN102680082A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明实施例提供一种结晶器振幅测试方法和装置,该方法包括:对结晶器振动过程中的加速度传感器得到的振动加速度值进行采样;将采集到的结晶器的振动加速度利用数值积分方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分方法进行积分得到结晶器振动的振幅;在对结晶器的振动加速度和振动速度进行数值积分之前,将被积分的数值按照时间顺序划分成多段;在每段内计算该段内的被积分数值的趋势值的平均值,并将该段内的每个被积分数值都减去获得的趋势值的平均值。

Description

一种结晶器振幅测试方法和装置
技术领域
本发明属于连铸技术领域,尤其涉及一种结晶器振幅测试方法和装置。
背景技术
结晶器时连铸机的“心脏”,需要不断地振动,这样铸坯才能强制脱模和消除粘连。目前结晶器的振动主要有正弦波式振动和非正弦波式振动两种方式。
结晶器的振动参数中的振幅可以衡量结晶器页面的稳定性,判断浇注是否易于控制,考量铸坯表面是否平滑。
结晶器振动测试通常利用加速度传感器获得原始的加速度采样值,再经过软件积分一次得到速度值,最后再积分一次得到位移值。由于现有的软件积分算法存在初始值不为零、加速度均值不为零、算法误差等问题,会使得经过上面两次积分得到的位移会产生趋势项,无法正确反映结晶器的振幅参数。
发明内容
有鉴于此,本发明的目的在于提供一种结晶器振幅测试方法和装置,能够得到更为准确的振幅值。
为实现上述目的,本发明的一个实施例提供一种结晶器振幅测试方法,包括:
对结晶器振动过程中的加速度传感器得到的振动加速度值进行采样;
将采集到的结晶器的振动加速度利用数值积分方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分方法进行积分得到结晶器振动的振幅;
在对结晶器的振动加速度和振动速度进行数值积分之前,将被积分的数值按照时间顺序划分成多段;
在每段内计算该段内的被积分数值的趋势值的平均值,并将该段内的每个被积分数值都减去获得的趋势值的平均值。
另一方面,本发明还提供一种结晶器振幅测试装置,包括:
加速度采样单元,用于对结晶器振动过程中的加速度传感器得到的振动加速度值进行采样;
积分单元,用于将采集到的结晶器的振动加速度利用数值积分的方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分的方法进行积分得到结晶器振动的振幅;
去趋势单元,用于在积分单元对结晶器的振动加速度和振动速度进行数值积分之前,将被积分的数值按照时间顺序划分成多段,并在每段内计算该段内的被积分数值的趋势值的平均值,并将该段内的每个被积分数值都减去获得的趋势值的平均值。
根据本发明实施例,在得到结晶器振幅的过程中,能够去除积分过程产生的趋势项,使得最终获得的结晶器振幅更正确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是积分含义的示意图;
图2是利用现有的数值积分方法计算得到的结晶器振幅的结果示意图;
图3是本发明提供的一种结晶器振幅测试方法的示意图;
图4是利用本发明实施例的方法获得的结晶器振幅的结果的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下简单说明利用数值积分即软件积分的方法计算结晶器振动的振幅的原理。
根据数学上的定律:
设函数F(x)在x轴上存在两点x1、x2,如图1围绕起来的封闭图形的面积S即是函数f(x)在[x1,x2]区间里的积分
Figure BDA00001679765900031
设x2-x1≥0,那么则有
Figure BDA00001679765900032
(可以称为梯形公式)。
但是在利用软件进行数值积分的过程中会存在如下问题:
1、积分初始值问题:无法获得初速度的值,在实际数值计算中只能假设初速度为0;
2、加速度总和不为零问题:类似sin曲线的加速度围绕起来的封闭图形的面积S不为0,加速度可能偏向于正象限,也有可能偏向于负象限,由此造成加速度总和不为零;
3、基于上述梯形公式计算积分本身也会具有误差,误差大小为: ∫ x 1 x 2 f ( x ) dx - x 2 - x 1 N [ f ( x 1 ) + f ( x 2 ) 2 + Σ k = 1 N f ( x 1 + k x 2 - x 1 N ) ] , 这里的N是在x1至x2之间所划分的梯形的个数,上述误差是在x1至x2之间因为梯形累加造成的误差总和。
由于上述问题的存在,会使得数值积分的结果中引入趋势项,例如,在计算结晶器振动振幅的时候,会得到如图2中示出的检测结果(蓝色曲线是加速度,绿色曲线是速度,红色曲线是振幅),由该检测结果可知第条线所反映的振幅带有趋势项的影响。
考虑到上述问题的存在,本发明提供一种结晶器振幅测试方法,如图3所示,该方法包括如下步骤:
步骤S301:对结晶器振动过程中加速度传感器得到的振动加速度值进行采样。
根据实际的需要上述采样频率或者在规定的时间内采样的数据点数可以灵活设定,根据实验发现,在采样频率为100个/秒时,最终得到的结晶器振动振幅更准确。
步骤S302:将采集到的结晶器的振动加速度利用数值积分方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分方法进行积分得到结晶器振动的振幅。
在步骤S302中如果直接对采集到的原始振动加速度采样值进行数值积分,就会在振动速度中引入趋势项,进一步地以该振动速度直接进行数值积分又会再次引入新的趋势项,两次数值积分的过程会造成与图2中相类似的问题。
因此本发明在步骤S302中在对结晶器的振动加速度和振动速度进行数值积分之前都要对将被积分的数值减去相应的趋势值,可以称为去趋势的过程,由此可见,本发明中在对原始的振动加速度进行数值积分之前,要对振动加速度去趋势,而将去趋势后的振动加速度值作为计算振动速度的数值积分过程的被积分数值;同理,在之后利用获得的振动速度进行数值积分之前,也要将直接计算得到的振动速度去趋势,将去趋势后的振动速度作为计算振动振幅的数值积分过程的被积分数值。
以下详细说明上述去趋势的过程,由于在两次积分过程之前分别需要对振动加速度和振动速度进行去趋势,而实际上两次去趋势的过程是相类似的,以下不妨以每次被去趋势的值(也就是数值积分过程的被积分的数值)来说明。
将被积分的数值按照时间顺序划分成多段,例如划分成N段,N>1。每段内独立进行去趋势,可以将每段称为一个去趋势窗体,每段内的数据点数量(可以称为窗体大小)可以记为WindowSize。
在实际中,可以通过一维先进先出(FIFO)链表TrendList将被积分的数值按照时间顺序划分成多段,该FIFO链表引入了趋势项的被积分数值。
针对每段的被积分数值,计算该段内的被积分数值的趋势值的平均值,该平均值的计算过程可以通过如下公式表示:
Σ 0 WindowSize Trend ( i ) / WindowSize
其中Trend(i)表示该段内的每个被积分数值的趋势值。
每个被积分数值的趋势值可以根据如下算法确定:
将多个Trend(i)求和,除以窗体大小。
将该段内的每个被积分数值都减去前面获得的趋势值的平均值。
至此完成被积分数值的去趋势过程。
在实际中为了防止一维FIFO链表TrendList耗尽内存,还在采集到1024个采样点时清空所述FIFO链表中的数据量为该链表大小减每段的数据量的数据,即清空数量为TrendList.size()-WindowSize的数据。
由此可见看出,本发明实施例提供的测试方法在得到结晶器振幅的过程中,能够去除积分过程产生的趋势项,使得最终获得的结晶器振幅更正确,例如图4示出了利用本发明实施例的方法获得的结晶器振幅的结果,由图4可以看出,本发明得到的结晶器振幅与客观情况比较相符,更为准确。
此外,本发明实施例还相应提供一种结晶器振幅测试装置,包括:
加速度采样单元,用于对结晶器振动过程中的加速度传感器得到的振动加速度值进行采样;
积分单元,用于将采集到的结晶器的振动加速度利用数值积分的方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分的方法进行积分得到结晶器振动的振幅;
去趋势单元,用于在积分单元对结晶器的振动加速度和振动速度进行数值积分之前,将被积分的数值按照时间顺序划分成多段,并在每段内计算该段内的被积分数值的趋势值的平均值,并将该段内的每个被积分数值都减去获得的趋势值的平均值。
在一个优选的实施例中,将被积分的数值按照时间顺序划分成多段具体实现为:通过一维先进先出FIFO链表将被积分的数值按照时间顺序划分成多段。
考虑到为了防止一维FIFO链表TrendList耗尽内存,上述装置还包括清空单元,用于在。。。。时清空所述FIFO链表中的数据量为该链表大小减每段的数据量的数据。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (6)

1.一种结晶器振幅测试方法,其特征在于,包括:
对结晶器振动过程中的加速度传感器得到的振动加速度值进行采样;
将采集到的结晶器的振动加速度利用数值积分方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分方法进行积分得到结晶器振动的振幅;
在对结晶器的振动加速度和振动速度进行数值积分之前,将被积分的数值按照时间顺序划分成多段;
在每段内计算该段内的被积分数值的趋势值的平均值,并将该段内的每个被积分数值都减去获得的趋势值的平均值。
2.根据权利要求1所述的方法,其特征在于,所述将被积分的数值按照时间顺序划分成多段具体实现为:通过一维先进先出FIFO链表将被积分的数值按照时间顺序划分成多段。
3.根据权利要求2所述的方法,其特征在于,所述方法还包括:在采集到1024个采样点时清空所述FIFO链表中的数据量为该链表大小减每段的数据量的数据。
4.一种结晶器振幅测试装置,其特征在于,包括:
加速度采样单元,用于对结晶器振动过程中的加速度传感器得到的振动加速度值进行采样;
积分单元,用于将采集到的结晶器的振动加速度利用数值积分的方法积分得到结晶器的振动速度,并将获得的振动速度利用数值积分的方法进行积分得到结晶器振动的振幅;
去趋势单元,用于在积分单元对结晶器的振动加速度和振动速度进行数值积分之前,将被积分的数值按照时间顺序划分成多段,并在每段内计算该段内的被积分数值的趋势值的平均值,并将该段内的每个被积分数值都减去获得的趋势值的平均值。
5.根据权利要求4所述的装置,其特征在于,所述将被积分的数值按照时间顺序划分成多段具体实现为:通过一维先进先出FIFO链表将被积分的数值按照时间顺序划分成多段。
6.根据权利要求5所述的装置,其特征在于,还包括:
清空单元,用于在采集到1024个采样点时清空所述FIFO链表中的数据量为该链表大小减每段的数据量的数据。
CN2012101645090A 2012-05-24 2012-05-24 一种结晶器振幅测试方法和装置 Pending CN102680082A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012101645090A CN102680082A (zh) 2012-05-24 2012-05-24 一种结晶器振幅测试方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012101645090A CN102680082A (zh) 2012-05-24 2012-05-24 一种结晶器振幅测试方法和装置

Publications (1)

Publication Number Publication Date
CN102680082A true CN102680082A (zh) 2012-09-19

Family

ID=46812354

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012101645090A Pending CN102680082A (zh) 2012-05-24 2012-05-24 一种结晶器振幅测试方法和装置

Country Status (1)

Country Link
CN (1) CN102680082A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102970003A (zh) * 2012-11-16 2013-03-13 中冶连铸技术工程股份有限公司 一种结晶器振幅信号的滤波方法及装置
CN106323451A (zh) * 2015-06-26 2017-01-11 陕西重型汽车有限公司 通过加速信号获得位移信号的方法及其装置
CN109084885A (zh) * 2018-08-14 2018-12-25 中国科学院上海高等研究院 一种机械设备振动检测系统、方法、可读存储介质及终端
CN110775758A (zh) * 2019-10-09 2020-02-11 浙江大学 一种基于轿厢加速度信号解析的电梯运行健康度评价方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007155335A (ja) * 2005-11-30 2007-06-21 Institute Of National Colleges Of Technology Japan プレス金型の動的変形測定法及び動的変形測定システム
CN101113936A (zh) * 2007-07-20 2008-01-30 广州市计量检测技术研究院 一种虚拟振动台检测信号处理方法及其设备
CN101221066A (zh) * 2008-01-25 2008-07-16 太原理工大学 工程非线性振动检测方法
CN101707512A (zh) * 2009-11-25 2010-05-12 华为技术有限公司 自适应调整数据速率的方法及装置
CN101722284A (zh) * 2009-09-17 2010-06-09 杭州谱诚泰迪实业有限公司 一种连铸结晶器振动检测装置及检测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007155335A (ja) * 2005-11-30 2007-06-21 Institute Of National Colleges Of Technology Japan プレス金型の動的変形測定法及び動的変形測定システム
CN101113936A (zh) * 2007-07-20 2008-01-30 广州市计量检测技术研究院 一种虚拟振动台检测信号处理方法及其设备
CN101221066A (zh) * 2008-01-25 2008-07-16 太原理工大学 工程非线性振动检测方法
CN101722284A (zh) * 2009-09-17 2010-06-09 杭州谱诚泰迪实业有限公司 一种连铸结晶器振动检测装置及检测方法
CN101707512A (zh) * 2009-11-25 2010-05-12 华为技术有限公司 自适应调整数据速率的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘业颜: "结晶器振动装置超低频振动分析系统研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》, no. 5, 15 May 2011 (2011-05-15) *
李东文等: "振动加速度信号处理探讨", 《机电工程技术》, no. 9, 15 December 2008 (2008-12-15) *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102970003A (zh) * 2012-11-16 2013-03-13 中冶连铸技术工程股份有限公司 一种结晶器振幅信号的滤波方法及装置
CN102970003B (zh) * 2012-11-16 2015-05-27 中冶连铸技术工程股份有限公司 一种结晶器振幅信号的滤波方法及装置
CN106323451A (zh) * 2015-06-26 2017-01-11 陕西重型汽车有限公司 通过加速信号获得位移信号的方法及其装置
CN106323451B (zh) * 2015-06-26 2019-10-11 陕西重型汽车有限公司 通过加速信号获得位移信号的方法及其装置
CN109084885A (zh) * 2018-08-14 2018-12-25 中国科学院上海高等研究院 一种机械设备振动检测系统、方法、可读存储介质及终端
CN110775758A (zh) * 2019-10-09 2020-02-11 浙江大学 一种基于轿厢加速度信号解析的电梯运行健康度评价方法
CN110775758B (zh) * 2019-10-09 2020-12-18 浙江大学 一种基于轿厢加速度信号解析的电梯运行健康度评价方法

Similar Documents

Publication Publication Date Title
CN103674062B (zh) 基于Allan方差与ARMA模型分析提高陀螺仪测量精度的方法
CN104375912B (zh) 移动终端卡顿的测量方法与装置
CN102680082A (zh) 一种结晶器振幅测试方法和装置
CN106523928B (zh) 基于声波实时数据二级筛选的管道泄漏检测方法
CN102735971A (zh) 一种基于fpga的电力系统同步向量测量计算的装置
EP1785702A3 (en) Method and device for measuring gas meter and volume corrector accuracy
CN102043091B (zh) 数字化高精度相位检测器
CN101930351B (zh) 基于cordic的信号处理方法和装置
CN110019386A (zh) 一种流数据处理方法及设备
CN103823177A (zh) 基于窗函数设计的滤波器性能检测方法和系统
CN102043090A (zh) 轨道电路移频信号参数测量系统
CN104048165B (zh) 管道泄漏诊断的方法
CN104880713B (zh) Gbas完好性评估处理方法和系统
KR101041990B1 (ko) 레이더 모의표적장치의 도플러 주파수 구현방법
EP3370074B1 (en) Method for detecting background noise of sensor, and device thereof
CN105628187B (zh) 一种振动信号的零偏校正方法、装置和振动监测仪
KR102062792B1 (ko) 회귀 분석 기반의 추정 알고리즘에 의하여 변형된 칼만 필터 알고리즘을 적용한 파고 데이터 연산을 사용하는 간이형 파고 측정 장치
CN109870404B (zh) 一种雨棚结构损伤识别方法、装置及终端设备
CN102970259A (zh) 一种测量td-lte终端矢量幅度误差的方法及装置
CN111177923A (zh) 一种用于蒸发器中结垢维护的预测方法及装置
CN106352974B (zh) 一种数字式的声级计脉冲计权方法及装置
CN111130681B (zh) 一种噪声传递特性确定方法和装置
CN102591616B (zh) 浮点计算性能确定装置和方法
CN208669907U (zh) 一种提高起重机制动下滑量精度的测量装置
CN106970357A (zh) 高精度的实时超声波定位追踪方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20120919