CN104268414A - 一种故障自修复的鲁棒的煤矿巷道剪应力采集方法 - Google Patents
一种故障自修复的鲁棒的煤矿巷道剪应力采集方法 Download PDFInfo
- Publication number
- CN104268414A CN104268414A CN201410516131.5A CN201410516131A CN104268414A CN 104268414 A CN104268414 A CN 104268414A CN 201410516131 A CN201410516131 A CN 201410516131A CN 104268414 A CN104268414 A CN 104268414A
- Authority
- CN
- China
- Prior art keywords
- foil gauge
- strain
- rod
- bending
- prime
- 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
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种故障自修复的鲁棒的煤矿巷道剪应力采集方法。煤矿现场通常用杆体来检测巷道的应变程度,主要测量杆体的剪应力分布,一般采用在杆体上安装应变片的方法,但是应变片数量多且容易损坏,而应变片数据采集都是自动采集,无人员实时监控,所以应变片发生故障损坏时不能被及时发现,得到错误的弯矩拟合函数和剪应力分布,这对于煤矿巷道观察被测杆体应变特性是非常不利的;其次,在煤矿现场安装杆体需要用特种水泥将其固定在巷道里,应变片发生故障损坏时将无法更换新的应变片,如不进行故障修复,少数应变片故障将造成杆体的浪费。
Description
技术领域
本发明专利涉及一种故障自修复的鲁棒的煤矿巷道剪应力采集方法,属于测控和自动化技术领域。
背景技术
剪应力是指物体由于外因(受力、湿度变化等)而变形时,在物体内各部分之间产生相互作用的内力,来抵抗外因的作用,并力图使物体从形变后的位置回复到形变前的位置。对于杆体,准确掌握整个长度范围内的剪应力分布,是研究整个杆体结构特性的重点。
煤矿巷道中,对于一些大型结构(桥梁等),通常测量内置杆体的整体形变来反映外部结构的形变,准确掌握大型结构的形变,对于及时做好预防措施、减少安全隐患具有重要意义。因此,准确检测杆体的形变,特别是准确掌握杆体剪应力分布至关重要。
测量杆体的剪应力分布,通常采用在杆体上安装应变片的方法,由于杆体一般用于测量大型结构,所以杆体比较长,所需要的应变片数量很大。由于应变片体积小且薄,容易在煤矿巷道的环境中损坏,而应变片数据采集都是自动采集,无人员实时监控,所以当应变片发生故障损坏时不能被及时发现,故障应变片的错误采集值仍参与计算,得到错误的弯矩函数拟合和剪应力分布,因此,工作人员检查结构时看到错误的杆体性质,会相应的采取错误的措施,影响工作效率,并且杆体一旦浇筑进大型结构中无法拿出,如果出现故障不能解决,那么杆体检测失去了意义。本发明提出一种故障自修复的鲁棒的煤矿巷道剪应力采集方法,通过实时监测所有应变片是否有故障,当应变片发生故障时对其进行相应的采集数据剔除或者进行误差模型补偿,使得被测杆体在无人监控的环境下实现故障自修复,保证弯矩拟合函数的准确性,得到准确的剪应力的函数,从而获得整个被测杆体上的剪应力分布。
发明内容
本发明解决的问题是:为了克服工业中杆体上应变片发生故障时检测出错误的杆体剪应力分布,影响杆体结构性质判断的缺点,结合以上背景和需求,本发明提供一种故障自修复的鲁棒的煤矿巷道剪应力采集方法,通过检测出有故障的应变片,对其采集数值进行误差补偿,用补偿后的采样数值进行拉伸应变、弯曲应变、弯矩的计算,然后对每个测量点的弯矩进行函数拟合获得整个被测杆体上的弯矩函数,由弯矩拟合函数得到沿测点长度范围剪应力的函数,根据这个函数可以很容易得到杆体上每个测点处的剪应力。
本发明的技术解决方案是:
1、一种故障自修复的鲁棒的煤矿巷道剪应力采集方法,其特征在于,包括
以下步骤:
步骤1:设煤矿巷道剪应力测试杆体长度L,在该杆体上贴N对应变片,每对应变片之间的间隔d=L/(N-1),每对应变片正对贴在杆体上,设上面一层应变片的编号为Xi1,正对面一层应变片编号为Xi2(i=0、1……N-1,i为第i对应变片),将该杆体安装在对其进行弯曲、拉伸的相关设备上;
步骤2:误差模型的输出值的评估:
步骤2.1:利用相关设备对装有2N个应变片的煤矿巷道剪应力测试杆体以时间间隔t1分别进行C次拉伸应变和C次弯曲应变,采集并记录弯曲应变、拉伸应变数值,记为εi1 拉伸,εi2 拉伸,εi1 弯曲,εi2 弯曲(i=0、1……N-1),此处时间间隔t1时间间隔较短,模拟的是煤矿巷道应变突变的情况下的数据采集;
步骤2.2:分别计算煤矿巷道剪应力测试杆体模拟煤矿巷道应变突变的情况下每一对应变片的拉伸应变和弯曲应变数值误差δij 拉伸=εi1 拉伸-εi2 拉伸,δij 弯曲=εi1 弯曲+εi2 弯曲(i=0、1……N-1,j=0、1……C-1);
步骤2.3:利用相关设备对装有2N个应变片的煤矿巷道剪应力测试杆体以时间间隔t2分别进行C′次拉伸应变和C′弯曲应变,采集并记录弯曲应变、拉伸应变数值,记为(i=0、1……N-1),此处时间间隔t2时间间隔较长,模拟的是煤矿巷道应变变化缓慢的情况下的数据采集;
步骤2.4:分别计算设煤矿巷道剪应力测试杆体模拟煤矿巷道应变变化缓慢的情况下每一对应变片的拉伸应变和弯曲应变数值误差(i=0、1……N-1,j′=0、1……C′-1);
步骤2.5:分别计算设煤矿巷道剪应力测试杆体每一对应变片拉伸应变和弯曲应变应的数值综合工况下误差平均值,(i=0、1……N-1),该平均值即为煤矿巷道综合工况下的应变误差模型输出值的估算值,该估算值即为应变偏差补偿的数据源;
步骤3:将剪应力测试杆体安装在煤矿巷道里,进行煤矿巷道综合工况下的数据采集。
步骤4:根据采集到的数据判断应变片是否故障:
步骤4.1:每隔时间t对杆体上的应变片采集数值进行读取;
步骤4.2:读取杆体上2N个应变片的采集数值εi1,εi2(i=0、1……N-1);
步骤4.3:如果连续三次εi1或者εi2值为0,则认为应变片发生故障,进入步骤5;
步骤4.4:如果εi1,εi2(i=0、1……N-1)都不为0,即应变片未出现故障,则此时有效应变片对数N′=N,进入步骤7;
步骤5:应变片故障类型判断:
步骤5.1:如果应变片Xi1出现故障,应变片Xi2正常工作或者应变片Xi2出现故障,应变片Xi1正常工作,进入步骤6.1;
步骤5.2:如果应变片Xi1和Xi2都出现故障,则记录应变片Xi1和Xi2都出现故障的应变片对数为P,进入步骤6.2;
步骤6:根据故障类型进入相应的应变片故障处理流程:
步骤6.1:对于故障应变片的采集数值进行误差补偿,如果Xi1出现故障,则计算拉伸应变时对Xi1应变片进行补偿,令则计算弯曲应变和弯矩时令则可以得到第i个应变片的弯矩:Mi 弯矩=π·E·R3·(ε′i1-εi2)/8,其中π是圆周率,E为弹性系数,R为杆体的截面半径;如果Xi2出现故障,则计算拉伸应变时对Xi2进行补偿,令则计算弯曲应变和弯矩时令则可以得到第i个应变片的弯矩:Mi 弯矩=π·E·R3·(εi1-ε′i2)/8,进入步骤8;
步骤6.2:如果第i对两个应变片都是有故障的,则剔除这一对应变片的采样数值,将剩下的正常工作的应变片作为有效测点,则修正后有效应变片对数N′=N-P进入步骤7;
步骤7:被测杆体应变参数计算,E为弹性系数,R为杆体的截面半径,A为杆体横截面积,i′为故障处理之后有效应变片的序号,εi′1,εi′2为步骤6中进行误差补偿之后的应变值:
步骤7.1:计算被测杆体上每个应变片测点拉伸应变参数,拉伸应变为拉伸应力为(i′=0、1……N′-1);
步骤7.2:;计算被测杆体上每个应变片测点弯曲应变参数,弯曲应变为弯曲应力为(i′=0、1……N′-1);
步骤7.3:;计算被测杆体每个有效应变片点弯矩:
Mi′=π·E·R3·(εi′1-εi′2)/8,(i′=0、1……N′-1) (1)
步骤8:对被测杆体的N′对有效应变片的弯矩Mi′进行函数拟合,在被测杆体由应变片分割的N′-1个区间段上使用三次样条插值拟合获取被测杆体上的弯矩函数,每个区间段上的拟合函数的一阶导数和二阶导数均为连续函数,则获得的拟合弯矩函数曲线是光滑曲线,具体拟合步骤如下:
步骤8.1:被测杆体的N′对有效应变片,相邻两个应变片之间是一个区间,所以有(N′-1)个区间,设每个区间段[zk,zk+1](k=0,1,...,N′-2),k表示(N′-1)个区间的序号。设第k个区间[zk,zk+1]的三次样条插值拟合函数为Sk(x),Sk(x)表示距离被测杆体端点为x处的弯矩,即:
Sk(x)=ak+bk(x-zk)+ck(x-zk)2+dk(x-zk)3 (2)
其中x是距离被测杆端点的距离,zk=k×d(k=0,1,...,N′-2)是当应变片i=k序号时,该应变片距离被测杆端点的距离,ak,bk,ck,dk为三次样条插值拟合函数需要求解的参数。
根据三次样条插值拟合函数的要求,其插值条件为:
当x=zk时,Sk(zk)=Mk (3)
其中,Mk是当i′=k时Mi′的值,Mi′的值从公式(1)得到。
三次样条插值拟合函数的连续性条件为:
三次样条插值拟合函数的一阶导数连续条件为:
其中,Sk′(x)是Sk(x)的一次导数。
三次样条插值拟合函数的二阶导数连续条件为:
其中,Sk″(x)是Sk(x)的二次导数。
以上插值条件和连续性条件总共能够得到(4N′-2)条等式。而三次样条插值在每一个区间上有4个未知参数,因此总共有4N′个待求解未知参数。被测杆体两端的端点位置没有受到让它们弯曲的力,因此端点条件为自由边界,所以边界端点约束提供了两个等式,为:
通过以上4N′条等式,能够求解出Sk(x)中4N′个未知参数ak,bk,ck,dk(k=0、1……N′-2)。
步骤8.2:求解三次样条插值拟合函数为Sk(x)。
步骤8.2.1:将拟合函数Sk(x)的二阶导数作为待定参数,令Tk=Sk″(x);
步骤8.2.2:因为Sk(x)是三次多项式,所以其二阶导数是一次多式,
令
其中,hk=zk+1-zk。
对二阶导数积分,得到一阶导数:
其中,pk为常量参数。
对一阶导数积分,得到函数Sk(x):
其中,qk为常量参数。
将x=zk带入公式(9)中,可以得到Sk(zk)=qk,qk即当x=zk时Sk(x)的值,根据公式(3)可知Sk(zk)等于序号为k的应变片的弯矩值,即Mk,因此qk是已知量。
因为拟合曲线平滑时端点处函数值相等,所以
Sk(zk+1)=Sk+1(zk+1)=qk+1(11)将公式(11)代入公式(10)中,可以得到:
由公式(12)得到
其中,
将公式(13)带入公式(10)中可以得到三次样条插值拟合函数为:
公式(15)中,Tk,Tk+1是未知量,需要继续求解。
步骤8.2.3:确定二阶导数待定参数Tk。
因为拟合曲线一次导数平滑所以端点处相邻两个函数值相等:
S′k-1(zk)=S′k(zk) (16)
将公式(16)代入公式(9)中,即:
化简公式(17),可以得到:
令:
则公式(18)可以化简为:
μkTk-1+2Tk+λkTk+1=6dk (19)
又因为边界条件公式(7)可以得到
可以得到求解Tk(k=1,2......N′-2)的矩阵方程,(N′-2)个方程可以求解出(N′-2)个Tk的值:
步骤8.3:将求得的Tk(k=0,1......N′-1)代入三次样条插值拟合函数Sk(x)公式(15)中,获得每一个区间段[zk,zk+1](k=0,1,...,N′-2)对应的弯矩函数:
步骤9:沿测点长度范围剪应力分布函数:
步骤9.1:对公式(22)的函数Sk(x)进行求导运算得到S′k(x),根据拟合的弯矩函数的导数S′k(x)得到每个区间段剪应力分布函数:
τk(x)=4·S′k(x)/(3π·R2) (23)
其中x是距离被测杆端点的距离,R是杆体横截面的半径;
步骤9.2:将被测杆体上任意一点距离杆体端点的距离z代入公式(23)中,可以得到该点的剪应力。
本发明的有益效果包括:
(1)本发明及时对应变片进行故障检测以及相应的故障处理,该发明将故障分为一对应变片中有一片发生故障及一对应变片中两片均发生故障,根据不同的故障类型采用不同的故障处理流程,保证了整个杆体剪应力分布的准确性,避免了因为应变片故障导致整个杆体应变检测失效的问题,延长了杆体检测的使用寿命;
(2)本发明采用自适应的故障处理策略,当判断故障类型为一对应变片中只有一个应变片发生故障的情况时将采用本发明中提出的煤矿巷道应变误差模型偏差补偿法对故障进行修复,若判断为一对应变片都发生故障的情况采用剔除该应变片数据,将剩下的应变片数据作为煤矿巷道剪应力计算的数据源,保证在分析整个杆体剪应力分布时不会受到错误数据干扰,并对部分情况下错误数据进行合理修复,提高煤矿巷道剪应力采集的鲁棒性,同时提高杆体的使用寿命;
(3)本发明对于一对应变片中只有一个应变片发生故障的情况采用了一种误差模型补偿的故障修复策略,该误差模型是基于模拟煤矿巷道应力突变及模拟煤矿巷道应力平缓综合数据采集的情况下建立的,保证了该误差模型更加接近煤矿巷道剪应力检测杆体在煤矿综合工况下的误差真实模型,保证了故障修复的鲁棒性及准确性;
(4)本发明对被测杆体的弯矩进行三次样条插值拟合,插值函数节点处的波动只对节点两侧的分段有影响,离该点较远的分段影响逐渐变小,有较好的稳定性。同时,由于个别端点的测量值异常不会引起被测杆体整体数据的异常,提高精度。
附图说明
图1为本发明的工作原理流程图;
具体实施方法
一种故障自修复的鲁棒的煤矿巷道剪应力采集方法,其特征在于,包括以下
步骤:
步骤1:设煤矿巷道剪应力测试杆体长度L,在该杆体上贴N对应变片,每对应变片之间的间隔d=L/(N-1),每对应变片正对贴在杆体上,设上面一层应变片的编号为Xi1,正对面一层应变片编号为Xi2(i=0、1……N-1,i为第i对应变片),将该杆体安装在对其进行弯曲、拉伸的相关设备上;
步骤2:误差模型的输出值的评估:
步骤2.1:利用相关设备对装有2N个应变片的煤矿巷道剪应力测试杆体以时间间隔t1分别进行C次拉伸应变和C次弯曲应变,采集并记录弯曲应变、拉伸应变数值,记为εi1 拉伸,εi2 拉伸,εi1 弯曲,εi2 弯曲(i=0、1……N-1),此处时间间隔t1时间间隔较短,模拟的是煤矿巷道应变突变的情况下的数据采集;
步骤2.2:分别计算煤矿巷道剪应力测试杆体模拟煤矿巷道应变突变的情况下每一对应变片的拉伸应变和弯曲应变数值误差δij 拉伸=εi1 拉伸-εi2 拉伸,δij 弯曲=εi1 弯曲+εi2 弯曲(i=0、1……N-1,j=0、1……C-1);
步骤2.3:利用相关设备对装有2N个应变片的煤矿巷道剪应力测试杆体以时间间隔t2分别进行C′次拉伸应变和C′弯曲应变,采集并记录弯曲应变、拉伸应变数值,记为(i=0、1……N-1),此处时间间隔t2时间间隔较长,模拟的是煤矿巷道应变变化缓慢的情况下的数据采集;
步骤2.4:分别计算设煤矿巷道剪应力测试杆体模拟煤矿巷道应变变化缓慢的情况下每一对应变片的拉伸应变和弯曲应变数值误差(i=0、1……N-1,j′=0、1……C′-1);
步骤2.5:分别计算设煤矿巷道剪应力测试杆体每一对应变片拉伸应变和弯曲应变应的数值综合工况下误差平均值,(i=0、1……N-1),该平均值即为煤矿巷道综合工况下的应变误差模型输出值的估算值,该估算值即为应变偏差补偿的数据源;
步骤3:将剪应力测试杆体安装在煤矿巷道里,进行煤矿巷道综合工况下的数据采集。
步骤4:根据采集到的数据判断应变片是否故障:
步骤4.1:每隔时间t对杆体上的应变片采集数值进行读取;
步骤4.2:读取杆体上2N个应变片的采集数值εi1,εi2(i=0、1……N-1);
步骤4.3:如果连续三次εi1或者εi2值为0,则认为应变片发生故障,进入步骤5;
步骤4.4:如果εi1,εi2(i=0、1……N-1)都不为0,即应变片未出现故障,则此时有效应变片对数N′=N,进入步骤7;
步骤5:应变片故障类型判断:
步骤5.1:如果应变片Xi1出现故障,应变片Xi2正常工作或者应变片Xi2出现故障,应变片Xi1正常工作,进入步骤6.1;
步骤5.2:如果应变片Xi1和Xi2都出现故障,则记录应变片Xi1和Xi2都出现故障的应变片对数为P,进入步骤6.2;
步骤6:根据故障类型进入相应的应变片故障处理流程:
步骤6.1:对于故障应变片的采集数值进行误差补偿,如果Xi1出现故障,则计算拉伸应变时对Xi1应变片进行补偿,令则计算弯曲应变和弯矩时令则可以得到第i个应变片的弯矩:Mi 弯矩=π·E·R3·(ε′i1-εi2)/8,其中π是圆周率,E为弹性系数,R为杆体的截面半径;如果Xi2出现故障,则计算拉伸应变时对Xi2进行补偿,令则计算弯曲应变和弯矩时令则可以得到第i个应变片的弯矩:Mi 弯矩=π·E·R3·(εi1-ε′i2)/8,进入步骤8;
步骤6.2:如果第i对两个应变片都是有故障的,则剔除这一对应变片的采样数值,将剩下的正常工作的应变片作为有效测点,则修正后有效应变片对数N′=N-P进入步骤7;
步骤7:被测杆体应变参数计算,E为弹性系数,R为杆体的截面半径,A为杆体横截面积,i′为故障处理之后有效应变片的序号,εi′1,εi′2为步骤6中进行误差补偿之后的应变值:
步骤7.1:计算被测杆体上每个应变片测点拉伸应变参数,拉伸应变为拉伸应力为(i′=0、1……N′-1);
步骤7.2:;计算被测杆体上每个应变片测点弯曲应变参数,弯曲应变为弯曲应力为(i′=0、1……N′-1);
步骤7.3:;计算被测杆体每个有效应变片点弯矩:
Mi′=π·E·R3·(εi′1-εi′2)/8,(i′=0、1……N′-1) (1)
步骤8:对被测杆体的N′对有效应变片的弯矩Mi′进行函数拟合,在被测杆体由应变片分割的N′-1个区间段上使用三次样条插值拟合获取被测杆体上的弯矩函数,每个区间段上的拟合函数的一阶导数和二阶导数均为连续函数,则获得的拟合弯矩函数曲线是光滑曲线,具体拟合步骤如下:
步骤8.1:被测杆体的N′对有效应变片,相邻两个应变片之间是一个区间,所以有(N′-1)个区间,设每个区间段[zk,zk+1](k=0,1,...,N′-2),k表示(N′-1)个区间的序号。设第k个区间[zk,zk+1]的三次样条插值拟合函数为Sk(x),Sk(x)表示距离被测杆体端点为x处的弯矩,即:
Sk(x)=ak+bk(x-zk)+ck(x-zk)2+dk(x-zk)3 (2)
其中x是距离被测杆端点的距离,zk=k×d(k=0,1,...,N′-2)是当应变片i=k序号时,该应变片距离被测杆端点的距离,ak,bk,ck,dk为三次样条插值拟合函数需要求解的参数。
根据三次样条插值拟合函数的要求,其插值条件为:
当x=zk时,Sk(zk)=Mk (3)
其中,Mk是当i′=k时Mi′的值,Mi′的值从公式(1)得到。
三次样条插值拟合函数的连续性条件为:
三次样条插值拟合函数的一阶导数连续条件为:
其中,Sk′(x)是Sk(x)的一次导数。
三次样条插值拟合函数的二阶导数连续条件为:
其中,Sk″(x)是Sk(x)的二次导数。
以上插值条件和连续性条件总共能够得到(4N′-2)条等式。而三次样条插值在每一个区间上有4个未知参数,因此总共有4N′个待求解未知参数。被测杆体两端的端点位置没有受到让它们弯曲的力,因此端点条件为自由边界,所以边界端点约束提供了两个等式,为:
通过以上4N′条等式,能够求解出Sk(x)中4N′个未知参数ak,bk,ck,dk(k=0、1……N′-2)。
步骤8.2:求解三次样条插值拟合函数为Sk(x)。
步骤8.2.1:将拟合函数Sk(x)的二阶导数作为待定参数,令
Tk=Sk″(x);
步骤8.2.2:因为Sk(x)是三次多项式,所以其二阶导数是一次多式,
令
其中,hk=zk+1-zk。
对二阶导数积分,得到一阶导数:
其中,pk为常量参数。
对一阶导数积分,得到函数Sk(x):
其中,qk为常量参数。
将x=zk带入公式(9)中,可以得到Sk(zk)=qk,qk即当x=zk时Sk(x)的值,根据公式(3)可知Sk(zk)等于序号为k的应变片的弯矩值,即Mk,因此qk是已知量。
因为拟合曲线平滑时端点处函数值相等,所以
Sk(zk+1)=Sk+1(zk+1)=qk+1 (11)
将公式(11)代入公式(10)中,可以得到:
由公式(12)得到
其中,
将公式(13)带入公式(10)中可以得到三次样条插值拟合函数为:
公式(15)中,Tk,Tk+1是未知量,需要继续求解。
步骤8.2.3:确定二阶导数待定参数Tk。
因为拟合曲线一次导数平滑所以端点处相邻两个函数值相等:
S′k-1(zk)=S′k(zk) (16)
将公式(16)代入公式(9)中,即:
化简公式(17),可以得到:
令:
则公式(18)可以化简为:
μkTk-1+2Tk+λkTk+1=6dk (19)
又因为边界条件公式(7)可以得到
可以得到求解Tk(k=1,2......N′-2)的矩阵方程,(N′-2)个方程可以求解出(N′-2)个Tk的值:
步骤8.3:将求得的Tk(k=0,1......N′-1)代入三次样条插值拟合函数Sk(x)公式(15)中,获得每一个区间段[zk,zk+1](k=0,1,...,N′-2)对应的弯矩函数:
步骤9:沿测点长度范围剪应力分布函数:
步骤9.1:对公式(22)的函数Sk(x)进行求导运算得到S′k(x),根据拟合的弯矩函数的导数S′k(x)得到每个区间段剪应力分布函数:
τk(x)=4·S′k(x)/(3π·R2) (23)
其中x是距离被测杆端点的距离,R是杆体横截面的半径;
步骤9.2:将被测杆体上任意一点距离杆体端点的距离z代入公式(23)中,可以得到该点的剪应力。
以上的实施方法只是已实现的有效的具体实施方式之一,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。
Claims (1)
1.一种故障自修复的鲁棒的煤矿巷道剪应力采集方法,其特征在于,包括以下步骤:
步骤1:设煤矿巷道剪应力测试杆体长度L,在该杆体上贴N对应变片,每对应变片之间的间隔d=L/(N-1),每对应变片正对贴在杆体上,设上面一层应变片的编号为Xi1,正对面一层应变片编号为Xi2(i=0、1……N-1,i为第i对应变片),将该杆体安装在对其进行弯曲、拉伸的相关设备上;
步骤2:误差模型的输出值的评估:
步骤2.1:利用相关设备对装有2N个应变片的煤矿巷道剪应力测试杆体以时间间隔t1分别进行C次拉伸应变和C次弯曲应变,采集并记录弯曲应变、拉伸应变数值,记为εi1 拉伸,εi2 拉伸,εi1 弯曲,εi2 弯曲(i=0、1……N-1),此处时间间隔t1时间间隔较短,模拟的是煤矿巷道应变突变的情况下的数据采集;
步骤2.2:分别计算煤矿巷道剪应力测试杆体模拟煤矿巷道应变突变的情况下每一对应变片的拉伸应变和弯曲应变数值误差δij 拉伸=εi1 拉伸-εi2 拉伸,δij 弯曲=εi1 弯曲+εi2 弯曲(i=0、1……N-1,j=0、1……C-1);
步骤2.3:利用相关设备对装有2N个应变片的煤矿巷道剪应力测试杆体以时间间隔t2分别进行C′次拉伸应变和C′弯曲应变,采集并记录弯曲应变、拉伸应变数值,记为εi1 ′拉伸,εi2 ′拉伸,εi1 ′弯曲,εi2 ′弯曲(i=0、1……N-1),此处时间间隔t2时间间隔较长,模拟的是煤矿巷道应变变化缓慢的情况下的数据采集;
步骤2.4:分别计算设煤矿巷道剪应力测试杆体模拟煤矿巷道应变变化缓慢的情况下每一对应变片的拉伸应变和弯曲应变数值误差δij ′拉伸=εi1 ′拉伸-εi2 ′拉伸,δij ′弯曲=εi1 ′弯曲+εi2 ′弯曲(i=0、1……N-1,j′=0、1……C′-1);
步骤2.5:分别计算设煤矿巷道剪应力测试杆体每一对应变片拉伸应变和弯曲应变应的数值综合工况下误差平均值, 该平均值即为煤矿巷道综合工况下的应变误差模型输出值的估算值,该估算值即为应变偏差补偿的数据源;
步骤3:将剪应力测试杆体安装在煤矿巷道里,进行煤矿巷道综合工况下的数据采集。
步骤4:根据采集到的数据判断应变片是否故障:
步骤4.1:每隔时间t对杆体上的应变片采集数值进行读取;
步骤4.2:读取杆体上2N个应变片的采集数值εi1,εi2(i=0、1……N-1);
步骤4.3:如果连续三次εi1或者εi2值为0,则认为应变片发生故障,进入步骤5;
步骤4.4:如果εi1,εi2(i=0、1……N-1)都不为0,即应变片未出现故障,则此时有效应变片对数N′=N,进入步骤7;
步骤5:应变片故障类型判断:
步骤5.1:如果应变片Xi1出现故障,应变片Xi2正常工作或者应变片Xi2出现故障,应变片Xi1正常工作,进入步骤6.1;
步骤5.2:如果应变片Xi1和Xi2都出现故障,则记录应变片Xi1和Xi2都出现故障的应变片对数为P,进入步骤6.2;
步骤6:根据故障类型进入相应的应变片故障处理流程:
步骤6.1:对于故障应变片的采集数值进行误差补偿,如果Xi1出现故障,则计算拉伸应变时对Xi1应变片进行补偿,令则计算弯曲应变和弯矩时令则可以得到第i个应变片的弯矩:Mi 弯矩=π·E·R3·(ε′i1-εi2)/8,其中π是圆周率,E为弹性系数,R为杆体的截面半径;如果Xi2出现故障,则计算拉伸应变时对Xi2进行补偿,令则计算弯曲应变和弯矩时令则可以得到第i个应变片的弯矩:Mi 弯矩=π·E·R3·(εi1-ε′i2)/8,进入步骤8;
步骤6.2:如果第i对两个应变片都是有故障的,则剔除这一对应变片的采样数值,将剩下的正常工作的应变片作为有效测点,则修正后有效应变片对数N′=N-P进入步骤7;
步骤7:被测杆体应变参数计算,i′为故障处理之后有效应变片的序号,εi′1,εi′2为步骤6中进行误差补偿之后的应变值:
步骤7.1:计算被测杆体上每个应变片测点拉伸应变参数,拉伸应变为拉伸应力为
步骤7.2:;计算被测杆体上每个应变片测点弯曲应变参数,弯曲应变为弯曲应力为
步骤7.3:;计算被测杆体每个有效应变片点弯矩:
Mi′=π·E·R3·(εi′1-εi′2)/8,(i′=0、1……N′-1) (1)
步骤8:对被测杆体的N′对有效应变片的弯矩Mi′进行函数拟合,在被测杆体由应变片分割的N′-1个区间段上使用三次样条插值拟合获取被测杆体上的弯矩函数,每个区间段上的拟合函数的一阶导数和二阶导数均为连续函数,则获得的拟合弯矩函数曲线是光滑曲线,具体拟合步骤如下:
步骤8.1:被测杆体的N′对有效应变片,相邻两个应变片之间是一个区间,所以有(N′-1)个区间,设每个区间段[zk,zk+1](k=0,1,...,N′-2),k表示(N′-1)个区间的序号。设第k个区间[zk,zk+1]的三次样条插值拟合函数为Sk(x),Sk(x)表示距离被测杆体端点为x处的弯矩,即:
Sk(x)=ak+bk(x-zk)+ck(x-zk)2+dk(x-zk)3 (2)
其中x是距离被测杆端点的距离,zk=k×d(k=0,1,...,N′-2)是当应变片i=k序号时,该应变片距离被测杆端点的距离,zk+1=(k+1)×d(k=0,1,...,N′-2)是当应变片i=k+1序号时,该应变片距离被测杆端点的距离,ak,bk,ck,dk为三次样条插值拟合函数需要求解的参数。
根据三次样条插值拟合函数的要求,其插值条件为:
当x=zk时,Sk(zk)=Mk (3)
其中,Mk是当i′=k时Mi′的值,Mi′的值从公式(1)得到。
三次样条插值拟合函数的连续性条件为:
三次样条插值拟合函数的一阶导数连续条件为:
其中,Sk'(x)是Sk(x)的一次导数。
三次样条插值拟合函数的二阶导数连续条件为:
其中,Sk″(x)是Sk(x)的二次导数。
以上插值条件和连续性条件总共能够得到(4N′-2)条等式。而三次样条插值在每一个区间上有4个未知参数,因此总共有4N′个待求解未知参数。被测杆体两端的端点位置没有受到让它们弯曲的力,因此端点条件为自由边界,所以边界端点约束提供了两个等式,为:
其中,z0=0×d=0表示当应变片i=0序号时,该应变片距离被测杆端点的距离为0,zN′-1=(N′-1)×d是当应变片i=N′-1序号时,该应变片距离被测杆端点的距离,通过以上4N′条等式,能够求解出Sk(x)中4N′个未知参数ak,bk,ck,dk(k=0、1……N′-2)。
步骤8.2:求解三次样条插值拟合函数为Sk(x)。
步骤8.2.1:将拟合函数Sk(x)的二阶导数作为待定参数,令Tk=Sk″(x);
步骤8.2.2:因为Sk(x)是三次多项式,所以其二阶导数是一次多式,
令
其中,hk=zk+1-zk。
对二阶导数积分,得到一阶导数:
其中,pk为常量参数。
对一阶导数积分,得到函数Sk(x):
其中,qk为常量参数。
将x=zk带入公式(9)中,可以得到Sk(zk)=qk,qk即当x=zk时Sk(x)的值,根据公式(3)可知Sk(zk)等于序号为k的应变片的弯矩值,即Mk,因此qk是已知量。
因为拟合曲线平滑时端点处函数值相等,所以
Sk(zk+1)=Sk+1(zk+1)=qk+1 (11)
将公式(11)代入公式(10)中,可以得到:
由公式(12)得到
其中,
将公式(13)带入公式(10)中可以得到三次样条插值拟合函数为:
公式(15)中,Tk,Tk+1是未知量,需要继续求解。
步骤8.2.3:确定二阶导数待定参数Tk。
因为拟合曲线一次导数平滑所以端点处相邻两个函数值相等:
S′k-1(zk)=S′k(zk) (16)
将公式(16)代入公式(9)中,即:
化简公式(17),可以得到:
令:
则公式(18)可以化简为:
μkTk-1+2Tk+λkTk+1=6dk(19)
又因为边界条件公式(7)可以得到
可以得到求解Tk(k=1,2......N′-2)的矩阵方程,(N′-2)个方程可以求解出(N′-2)个Tk的值:
步骤8.3:将求得的Tk(k=0,1......N′-1)代入三次样条插值拟合函数Sk(x)公式(15)中,获得每一个区间段[zk,zk+1](k=0,1,...,N′-2)对应的弯矩函数:
步骤9:沿测点长度范围剪应力分布函数:
步骤9.1:对公式(22)的函数Sk(x)进行求导运算得到S′k(k),根据拟合的弯矩函数的导数S′k(x)得到每个区间段剪应力分布函数:
τk(x)=4·S′k(x)/(3π·R2) (23)
其中x是距离被测杆端点的距离,R是杆体横截面的半径;
步骤9.2:将被测杆体上任意一点距离杆体端点的距离z代入公式(23)中,可以得到该点的剪应力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410516131.5A CN104268414B (zh) | 2014-09-29 | 2014-09-29 | 一种故障自修复的鲁棒的煤矿巷道剪应力采集方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410516131.5A CN104268414B (zh) | 2014-09-29 | 2014-09-29 | 一种故障自修复的鲁棒的煤矿巷道剪应力采集方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104268414A true CN104268414A (zh) | 2015-01-07 |
CN104268414B CN104268414B (zh) | 2017-04-19 |
Family
ID=52159935
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410516131.5A Active CN104268414B (zh) | 2014-09-29 | 2014-09-29 | 一种故障自修复的鲁棒的煤矿巷道剪应力采集方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104268414B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760938A (zh) * | 2016-03-18 | 2016-07-13 | 江苏省特种设备安全监督检验研究院 | 采集周期自适应的鲁棒的电梯曳引钢丝绳应变检测方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4579007A (en) * | 1983-11-04 | 1986-04-01 | Sfernice Societe Francaise De L'electro-Resistance | Dynamometer comprising an elastic bar provided with strain gages |
CN201497630U (zh) * | 2009-09-23 | 2010-06-02 | 河海大学 | 基于无线转矩测量的有载分接开关传动部分故障诊断系统 |
CN102072926A (zh) * | 2010-11-30 | 2011-05-25 | 浙江大学 | 一种发动机机体疲劳裂纹诊断的方法 |
CN102466531A (zh) * | 2010-11-10 | 2012-05-23 | 三一电气有限责任公司 | 一种螺栓故障监测系统和监测方法 |
CN102730571A (zh) * | 2012-06-01 | 2012-10-17 | 华中科技大学 | 一种起重机在线监测与故障诊断系统 |
CN103411659A (zh) * | 2013-08-12 | 2013-11-27 | 国电联合动力技术有限公司 | 一种风力发电机叶片与塔筒状态监测方法及系统 |
CN203704983U (zh) * | 2014-03-08 | 2014-07-09 | 湖南科技大学 | 一种基于无线传输的风力机叶根应力和叶片振动检测装置 |
-
2014
- 2014-09-29 CN CN201410516131.5A patent/CN104268414B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4579007A (en) * | 1983-11-04 | 1986-04-01 | Sfernice Societe Francaise De L'electro-Resistance | Dynamometer comprising an elastic bar provided with strain gages |
CN201497630U (zh) * | 2009-09-23 | 2010-06-02 | 河海大学 | 基于无线转矩测量的有载分接开关传动部分故障诊断系统 |
CN102466531A (zh) * | 2010-11-10 | 2012-05-23 | 三一电气有限责任公司 | 一种螺栓故障监测系统和监测方法 |
CN102072926A (zh) * | 2010-11-30 | 2011-05-25 | 浙江大学 | 一种发动机机体疲劳裂纹诊断的方法 |
CN102730571A (zh) * | 2012-06-01 | 2012-10-17 | 华中科技大学 | 一种起重机在线监测与故障诊断系统 |
CN103411659A (zh) * | 2013-08-12 | 2013-11-27 | 国电联合动力技术有限公司 | 一种风力发电机叶片与塔筒状态监测方法及系统 |
CN203704983U (zh) * | 2014-03-08 | 2014-07-09 | 湖南科技大学 | 一种基于无线传输的风力机叶根应力和叶片振动检测装置 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760938A (zh) * | 2016-03-18 | 2016-07-13 | 江苏省特种设备安全监督检验研究院 | 采集周期自适应的鲁棒的电梯曳引钢丝绳应变检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104268414B (zh) | 2017-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106556498B (zh) | 桥梁结构损伤识别方法及系统 | |
CN101221104B (zh) | 基于分布式应变动态测试的结构健康监测方法 | |
TWI449883B (zh) | 結構體安全性之分析方法 | |
CN104457681B (zh) | 一种基于应变模态的梁结构动挠度监测方法 | |
CN101782372B (zh) | 基于梁端纵向位移的桥梁伸缩缝损伤诊断智能方法 | |
CN105716814B (zh) | 一种评估桁架结构损伤的实时监测系统及其方法 | |
CN106932135B (zh) | 一种基于加权窄带搜峰法识别振动频率的柔性拉索索力测试方法 | |
WO2019001016A1 (zh) | 一种基于温度与位移关系模型的桥梁伸缩缝性能预警方法 | |
CN101806668B (zh) | 一种基于索力监测的索结构健康监测方法 | |
Ding et al. | Assessment of bridge expansion joints using long-term displacement measurement under changing environmental conditions | |
CN110728089B (zh) | 基于botda技术的大跨桥梁斜拉索结构损伤诊断方法 | |
CN109858112B (zh) | 基于结构应力监测结果的数值反演分析方法 | |
RU2645903C1 (ru) | Способ контроля напряженно-деформированного состояния конструктивных элементов массивных бетонных сооружений при длительной эксплуатации | |
CN107356417B (zh) | 一种融合时序分析与信息熵的栓接结合部损伤识别方法 | |
CN202382726U (zh) | 复合梁应变测试系统 | |
CN110672154A (zh) | 土木工程建筑监测系统 | |
CN104880217A (zh) | 一种基于测量值关联度的故障传感器信息重构方法 | |
CN105158300A (zh) | 一种桥梁线状钢构件检测方法 | |
Iliopoulos et al. | Continuous fatigue assessment of an offshore wind turbine using a limited number of vibration sensors | |
CN107421672B (zh) | 一种基于振动频率全域搜峰的加权索力计算方法 | |
CN113483723B (zh) | 基于被动激励的桥梁结构应变监测系统在线校准方法 | |
CN105091923B (zh) | 用于操作测量点的方法 | |
Nyman et al. | Smart condition monitoring of a steel bascule railway bridge | |
CN117852122A (zh) | 一种pc桥梁结构内钢绞线现存应力检测的方法 | |
CN104268414A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |