发明名称 一种采煤机摇臂振动传感器优化布置方法
摘要 本发明公开了一种监测采煤机摇臂振动信号的传感器优化布置方法。该方法首先建立采煤机摇臂壳体的有限元模型;再进行基于功率流方法的最优布置点的区域划分;初选出传感器的监测点,然后运用粒子群优化算法,基于模态置信准则,对传感器的数量和具体监测点进行优化计算,得到振动传感器的优化布置结果。该方法具有普适性,一方面保证了各传感器测量信号的有效性,另一方面也保证了各监测点的模态向量有较大的空间交角,可满足采煤机摇臂的在线监测和故障诊断的要求。
申请公布号 CN103413185B 申请公布日期 2016.04.06
申请号 CN201310342578.0 申请日期 2013.08.07
申请人 中国矿业大学 发明人 程刚;倪建明;张业林;刘后广;邱锦波;吴行标;周常飞
分类号 G06F17/00(2006.01)I 主分类号 G06F17/00(2006.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 一种采煤机摇臂振动传感器优化布置方法,其特征在于,包括以下步骤:(1)建立摇臂壳体的有限元模型;(2)基于功率流方法的可布置区域划分;21)运用有限元分析软件对摇臂壳体进行谐响应分析,得到每个单元的受力;22)提取每个单元的受力,根据板壳功率流理论计算出摇臂壳体中功率流分布;具体为:板壳单元横向振动的功率流表达式如下所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>P</mi><mrow><mi>x</mi><mn>1</mn></mrow></msub><mo>=</mo><mo>&lt;</mo><msub><mi>Q</mi><mi>x</mi></msub><mover><mi>w</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>M</mi><mi>x</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>y</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub><mo>-</mo><mo>&lt;</mo><msub><mi>M</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>x</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>y</mi><mn>1</mn></mrow></msub><mo>=</mo><mo>&lt;</mo><msub><mi>Q</mi><mi>y</mi></msub><mover><mi>w</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>M</mi><mi>y</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>x</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub><mo>-</mo><mo>&lt;</mo><msub><mi>M</mi><mrow><mi>y</mi><mi>x</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>y</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000900460300000011.GIF" wi="837" he="253" /></maths>式中,P<sub>x1</sub>为x方向截面的横向功率流,P<sub>y1</sub>为y方向截面的横向功率流,w为单元的z方向的位移,Q<sub>x</sub>、Q<sub>y</sub>分别为单元x方向和y方向的剪切力,M<sub>x</sub>、M<sub>y</sub>分别为单元x方向和y方向的弯矩,θ<sub>x</sub>、θ<sub>y</sub>分别为单元的x、y方向的角位移;板壳单元纵向振动的功率流表达式如下所示:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><msub><mi>P</mi><mrow><mi>x</mi><mn>2</mn></mrow></msub><mo>=</mo><mo>&lt;</mo><msub><mi>N</mi><mi>x</mi></msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>N</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>y</mi><mn>2</mn></mrow></msub><mo>=</mo><mo>&lt;</mo><msub><mi>N</mi><mi>y</mi></msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>N</mi><mrow><mi>y</mi><mi>x</mi></mrow></msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000900460300000012.GIF" wi="596" he="253" /></maths>式中,P<sub>x2</sub>为x方向截面的纵向功率流,P<sub>y2</sub>为y方向截面的纵向功率流,u、v分别为单元的x和y方向的位移,N<sub>x</sub>,N<sub>y</sub>和N<sub>xy</sub>=N<sub>yx</sub>为单元所受的薄膜力;板壳单元总振动功率流表达式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>P</mi><mi>x</mi></msub><mo>=</mo><mo>&lt;</mo><msub><mi>N</mi><mi>x</mi></msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>N</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>Q</mi><mi>x</mi></msub><mover><mi>w</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>M</mi><mi>x</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>y</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub><mo>-</mo><mo>&lt;</mo><msub><mi>M</mi><mrow><mi>x</mi><mi>y</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>x</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>P</mi><mi>y</mi></msub><mo>=</mo><mo>&lt;</mo><msub><mi>N</mi><mi>y</mi></msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>N</mi><mrow><mi>y</mi><mi>x</mi></mrow></msub><mover><mi>u</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>Q</mi><mi>y</mi></msub><mover><mi>w</mi><mo>&CenterDot;</mo></mover><msub><mo>&gt;</mo><mi>t</mi></msub><mo>-</mo><mo>&lt;</mo><msub><mi>M</mi><mi>y</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>x</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub><mo>+</mo><mo>&lt;</mo><msub><mi>M</mi><mrow><mi>y</mi><mi>x</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>y</mi></msub><msub><mo>&gt;</mo><mi>t</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000900460300000021.GIF" wi="1252" he="252" /></maths>式中,P<sub>x</sub>为x方向截面的总功率流,P<sub>y</sub>为y方向截面的总功率流;23)根据功率流计算结果,划分出振动传感器的可布置区域;具体为:首先在ANSYS后处理器中编制程序提取每个单元节点的应力,然后分别根据板壳单元横向功率流、纵向功率流和总功率流编制计算程序,计算壳体中的功率流,最后绘制功率流云图和矢量图;(3)基于粒子群算法的传感器数量及位置优化31)根据步骤21)分析出的传感器可布置区域选取m个监测点作为初始监测点;32)运用有限元分析软件对摇臂壳体进行模态分析,得到m个监测点的各阶模态振型,并将所得到的各阶振型组成模态向量;33)根据模态置信准则(MAC),从m个监测点中选取n个监测点作为摇臂壳体的监测点;具体为以MAC矩阵非对角元素最大值作为优化对象,以测点数为粒子,应用粒子群优化算法进行传感器位置寻优。
地址 221116 江苏省徐州市泉山区大学路1号