[go: up one dir, main page]

CN111323011A - Cooperative positioning device and positioning method of shearer body and rocker arm - Google Patents

Cooperative positioning device and positioning method of shearer body and rocker arm Download PDF

Info

Publication number
CN111323011A
CN111323011A CN202010327977.XA CN202010327977A CN111323011A CN 111323011 A CN111323011 A CN 111323011A CN 202010327977 A CN202010327977 A CN 202010327977A CN 111323011 A CN111323011 A CN 111323011A
Authority
CN
China
Prior art keywords
rocker arm
fuselage
inertial navigation
strapdown inertial
shearer
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
CN202010327977.XA
Other languages
Chinese (zh)
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.)
Xijing University
Original Assignee
Xijing University
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 Xijing University filed Critical Xijing University
Priority to CN202010327977.XA priority Critical patent/CN111323011A/en
Publication of CN111323011A publication Critical patent/CN111323011A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

A co-location device and a co-location method for a coal mining machine body and a rocker arm comprise a body strapdown inertial navigation device and a wireless sensing anchor node which are arranged on the coal mining machine body; the rocker arm strapdown inertial navigation device and the wireless sensing mobile node are arranged on the rocker arm on the coal cutting side of the coal cutter; the wireless sensing anchor node receives the wireless signal; during work, the machine body moves along the track, the rocker arm cuts a coal wall to move, the rocker arm strapdown inertial navigation device measures angular velocity and acceleration information of the rocker arm of the coal mining machine, the machine body strapdown inertial navigation device measures the angular velocity and acceleration information of the machine body of the coal mining machine in real time, and the strapdown navigation algorithm is used for positioning and attitude determination; the invention has the advantages of high positioning and attitude determination precision, high fully mechanized mining efficiency and remote automatic control.

Description

一种采煤机机身与摇臂协同定位装置及定位方法Cooperative positioning device and positioning method of shearer body and rocker arm

技术领域technical field

本发明属于采煤机定位技术领域,具体涉及一种采煤机机身与摇臂协同定位装置及定位方法。The invention belongs to the technical field of shearer positioning, and in particular relates to a cooperative positioning device and a positioning method of a shearer body and a rocker arm.

背景技术Background technique

采煤机、液压支架和刮板输送机是井下综采工作面最重要的三种设备,相互配合完成割煤、运煤和支护工作。其中采煤机为主导设备,是综采工作面割煤和装煤的主要设备,是一种高集成度综采设备。采煤机在工作时沿刮板输送机的轨道进行往复式割煤,液压支架支护顶板并推进工作面。为了实现综采工作面自动化、远程自动控制,需要对采煤机精确动态定位,而且,针对自动化采煤需求,需要同时对采煤机的机身和摇臂进行定位,所以采煤机机身和摇臂的定位技术是煤矿生产装备自动化的关键技术。Shearer, hydraulic support and scraper conveyor are the three most important equipments in underground fully mechanized mining face, which cooperate with each other to complete coal cutting, coal transportation and support work. Among them, the shearer is the main equipment, which is the main equipment for coal cutting and coal loading in the fully mechanized mining face, and is a highly integrated fully mechanized mining equipment. During operation, the shearer performs reciprocating coal cutting along the track of the scraper conveyor, and the hydraulic support supports the roof and pushes the working face. In order to realize the automation and remote automatic control of the fully mechanized mining face, it is necessary to precisely and dynamically position the shearer. Moreover, for the needs of automatic coal mining, the fuselage and rocker arm of the shearer need to be positioned at the same time. And the positioning technology of rocker arm is the key technology of coal mine production equipment automation.

煤矿综采工作面工况复杂,同时空间封闭,因此采煤机机身与摇臂定位是一个典型的复杂封闭环境下室内定位问题。在采煤机的采煤作业中,采煤机机身沿轨道移动;而摇臂进行割煤工作,其相对于采煤机机身存在相对运动。为了实现自动化采煤,必须同时测量采煤机机身和摇臂的位置和姿态。The working conditions of the fully mechanized coal mining face are complex and the space is closed, so the positioning of the shearer body and rocker arm is a typical indoor positioning problem in a complex closed environment. In the coal mining operation of the shearer, the body of the shearer moves along the track; while the rocker arm performs the coal cutting work, there is a relative movement relative to the body of the shearer. In order to realize automatic coal mining, the position and attitude of the shearer body and rocker arm must be measured at the same time.

目前采煤机所用的定位方法主要包括捷联惯导定位法、红外定位法、超声波定位法、齿轮计数定位法、无线传感器网络定位法等。At present, the positioning methods used by shearers mainly include strapdown inertial navigation positioning method, infrared positioning method, ultrasonic positioning method, gear counting positioning method, wireless sensor network positioning method and so on.

红外定位法由安装在采煤机上红外发射装置发射信号,液压支架上安装的接收装置接收信号,利用红外测距定位采煤机的位置,但其红外信号易受粉尘影响,且存在固有定位盲区,定位精度不高,因此使用有局限性。In the infrared positioning method, the infrared transmitting device installed on the shearer transmits the signal, the receiving device installed on the hydraulic support receives the signal, and uses the infrared ranging to locate the position of the shearer, but the infrared signal is easily affected by dust, and there is an inherent positioning blind spot. , the positioning accuracy is not high, so the use is limited.

超声波定位法将超声波发射装置安装在工作面巷道中,当采煤机经过时,机身发射超声波,根据各位置超声波接收装置接收信号,利用超声波测距定位采煤机的位置,超声波的优点是可以穿透粉尘,但由于工作面较长,信号失信严重,其定位精度不高,因此使用有局限性。The ultrasonic positioning method installs the ultrasonic transmitting device in the roadway of the working face. When the shearer passes by, the fuselage emits ultrasonic waves, and the ultrasonic receiving device receives the signal according to each position, and uses the ultrasonic ranging to locate the position of the shearer. The advantages of ultrasonic waves are: It can penetrate dust, but due to the long working surface, serious signal loss, and low positioning accuracy, there are limitations in its use.

齿轮计数定位法是对采煤机行走齿轮的转动圈数进行计数,根据转动圈数和齿轮圆周长计算采煤机沿输送机轨道方向上的位移;但该方法只能用于定位采煤机沿轨道方向的一维位置,且受到齿轮计数误差影响,不能满足三维定位需要。The gear counting and positioning method is to count the number of rotations of the walking gear of the shearer, and calculate the displacement of the shearer along the conveyor track according to the number of rotations and the circumference of the gear; but this method can only be used to locate the shearer The one-dimensional position along the track direction, and affected by the gear count error, cannot meet the needs of three-dimensional positioning.

无线传感器网络定位法是在液压支架上布置多个位置已知的无线传感器(称为锚节点),在采煤机上布置待定位节点(称为移动节点),移动节点发射无线信号,锚节点接收无线信号监测采煤机与液压支架间的位置关系,解算出采煤机的位置;但由于工作面环境复杂,无线定位数据不稳定,且锚节点随液压支架移动后会导致自身位置变化,需要更新锚节点位置信息,同时不能进行采煤机定姿,不能满足实时定位定姿需求。The wireless sensor network positioning method is to arrange multiple wireless sensors (called anchor nodes) with known positions on the hydraulic support, and arrange the nodes to be located (called mobile nodes) on the shearer, the mobile nodes transmit wireless signals, and the anchor nodes receive The wireless signal monitors the positional relationship between the shearer and the hydraulic support, and calculates the position of the shearer; however, due to the complex working face environment, the wireless positioning data is unstable, and the anchor node will change its position after moving with the hydraulic support. The anchor node position information is updated, and at the same time, the shearer attitude cannot be determined, which cannot meet the needs of real-time positioning and attitude determination.

捷联惯导定位方法是一种全自主导航定位方法,无需借助外部信息,利用捷联惯导装置的三轴陀螺和三轴加速度计实时测量采煤机的角速度和线加速度,结合初始装订信息,通过姿态更新算法先解算出采煤机的运动姿态,然后将加速度根据姿态信息投影到导航坐标系,通过积分和二次积分获得采煤机的速度和位置等信息;高精度惯导虽然满足精度要求,但由于体积、重量等原因,一般只安装在采煤机的机身,无法用于采煤机摇臂的定位和定姿,或采煤机摇臂只能安装小型低精度惯导,其定位、定姿精度不能满足要求。The strapdown inertial navigation positioning method is a fully autonomous navigation and positioning method, without external information, using the three-axis gyroscope and three-axis accelerometer of the strapdown inertial navigation device to measure the angular velocity and linear acceleration of the shearer in real time, combined with the initial binding information , through the attitude update algorithm, the motion attitude of the shearer is first solved, and then the acceleration is projected to the navigation coordinate system according to the attitude information, and the speed and position of the shearer are obtained through integration and quadratic integration. Accuracy requirements, but due to volume, weight and other reasons, it is generally only installed on the fuselage of the shearer, and cannot be used for the positioning and attitude determination of the shearer rocker arm, or the shearer rocker arm can only be installed with small and low-precision inertial navigation. , its positioning and attitude accuracy cannot meet the requirements.

发明内容SUMMARY OF THE INVENTION

为了克服上述现有技术存在的问题,本发明提供了一种采煤机机身与摇臂协同定位装置及定位方法,能够利用采煤机机身和摇臂上的捷联惯导装置和无线测距传感器,同时实现采煤机机身与摇臂的定位定姿,具有定位定姿精度高、综采工作效率高及远程自动控制的优点。In order to overcome the above-mentioned problems in the prior art, the present invention provides a cooperative positioning device and a positioning method of the shearer body and the rocker arm, which can utilize the strapdown inertial navigation device and wireless The distance measuring sensor can realize the positioning and attitude of the shearer body and the rocker arm at the same time, and has the advantages of high positioning and attitude accuracy, high working efficiency of fully mechanized mining and remote automatic control.

为了实现上述目的,本发明采用如下的技术方案:In order to achieve the above object, the present invention adopts the following technical scheme:

一种采煤机机身与摇臂协同定位装置,包括在采煤机机身5上安装的机身捷联惯导装置1和无线传感锚节点2;以及在采煤机割煤一侧的摇臂6上安装的摇臂捷联惯导装置3和无线传感移动节点4;其中无线传感移动节点4发射无线信号,无线传感锚节点2接收无线信号。A cooperative positioning device for a shearer body and a rocker arm, comprising a body strapdown inertial navigation device 1 and a wireless sensor anchor node 2 installed on the shearer body 5; and a coal cutting side of the shearer The rocker arm strapdown inertial navigation device 3 and the wireless sensor mobile node 4 are installed on the rocker arm 6; wherein the wireless sensor mobile node 4 transmits wireless signals, and the wireless sensor anchor node 2 receives wireless signals.

所述机身捷联惯导装置1包括三只陀螺和三只加速度计。The fuselage strapdown inertial navigation device 1 includes three gyroscopes and three accelerometers.

基于上述的采煤机机身与摇臂协同定位装置的定位方法,具体包括以下步骤:Based on the above-mentioned positioning method of the shearer fuselage and the rocker arm cooperative positioning device, it specifically includes the following steps:

步骤一、定义坐标系Step 1. Define the coordinate system

1)定义采煤机机身坐标系b1:Ob1Xb1Yb1Zb1:坐标系原点Ob1固连在机身捷联惯导装置1中心,Xb1轴正向由采煤机指向煤壁,Yb1轴垂直于Xb1轴向上,Zb1轴与Xb1轴、Yb1轴构成右手坐标系,构成前上右坐标系,下标1表示该坐标系为机身坐标系;1) Define the shearer fuselage coordinate system b1: O b1 X b1 Y b1 Z b1 : the origin of the coordinate system O b1 is fixedly connected to the center of the fuselage strapdown inertial navigation device 1, and the positive X b1 axis points from the shearer to the coal wall, the Y b1 axis is perpendicular to the X b1 axis upward, the Z b1 axis, the X b1 axis, and the Y b1 axis form a right-hand coordinate system, which constitutes the front upper right coordinate system, and the subscript 1 indicates that the coordinate system is the fuselage coordinate system;

2)定义采煤机摇臂坐标系b2:Ob2Xb2Yb2Zb2;坐标系原点Ob2固连在摇臂捷联惯导装置3中心,Xb2轴正向由摇臂指向煤壁,Yb2轴垂直于Xb2轴向上,Zb2轴与Xb2轴、Yb2轴构成右手坐标系,构成前上右坐标系,下标2表示该坐标系为摇臂坐标系;2) Define the shearer rocker arm coordinate system b2: O b2 X b2 Y b2 Z b2 ; the origin of the coordinate system O b2 is fixedly connected to the center of the rocker arm strapdown inertial navigation device 3, and the X b2 axis is directed to the coal wall by the rocker arm. , the Y b2 axis is perpendicular to the X b2 axis upward, the Z b2 axis, the X b2 axis, and the Y b2 axis form a right-hand coordinate system, forming a front upper right coordinate system, and the subscript 2 indicates that the coordinate system is a rocker arm coordinate system;

3)定义导航坐标系n:OnXnYnZn:北天东地理坐标系,Xn轴指向地理北向,Yn轴指向天向,Zn轴指向地理东向;3) Define the navigation coordinate system n: O n X n Y n Z n : the north celestial east geographic coordinate system, the X n axis points to the geographic north direction, the Y n axis points to the sky direction, and the Z n axis points to the geographic east direction;

4)导航坐标系n经过三次旋转后与采煤机机身坐标系b1重合,三次旋转的角度即为机身航向角ψ1、俯仰角θ1和滚动角γ1;类似地,导航坐标系n经过三次旋转后与采煤机摇臂坐标系b2重合,三次旋转的角度即为摇臂航向角ψ2、俯仰角θ2和滚动角γ24) The navigation coordinate system n coincides with the shearer fuselage coordinate system b1 after three rotations, and the angles of the three rotations are the fuselage heading angle ψ 1 , pitch angle θ 1 and roll angle γ 1 ; similarly, the navigation coordinate system After three rotations, n coincides with the shearer rocker arm coordinate system b2, and the angles of the three rotations are rocker arm heading angle ψ 2 , pitch angle θ 2 and roll angle γ 2 ;

步骤二、计算采煤机机身信息Step 2. Calculate the information of the shearer fuselage

采煤机机身捷联惯导装置1利用陀螺和加速度计测量信息,解算采煤机机身姿态俯仰角θ1、航向角ψ1和滚动角γ1,计算机身的速度

Figure BDA0002463896770000041
位置
Figure BDA0002463896770000042
和经度λ1、纬度
Figure BDA0002463896770000043
高度h1;The shearer fuselage strapdown inertial navigation device 1 uses the gyro and accelerometer measurement information to calculate the shearer fuselage attitude pitch angle θ 1 , heading angle ψ 1 and roll angle γ 1 , and calculate the speed of the body
Figure BDA0002463896770000041
Location
Figure BDA0002463896770000042
and longitude λ 1 , latitude
Figure BDA0002463896770000043
height h 1 ;

1)机身捷联惯导装置1接收初始装订信息:初始位置

Figure BDA0002463896770000044
和初始姿态
Figure BDA0002463896770000045
其中上标1表示机身捷联惯导,N,U,E分别表示导航坐标系下北向、天向和东向位置;同时由于采煤机初始处于停止状态,初始速度V1=[0 0 0]T;1) Body SINS 1 receives initial binding information: initial position
Figure BDA0002463896770000044
and initial posture
Figure BDA0002463896770000045
The superscript 1 represents the fuselage strapdown inertial navigation, and N, U, and E represent the north, sky and east positions in the navigation coordinate system, respectively. 0] T ;

2)机身捷联惯导装置1的三只陀螺测量采煤机机身三个轴的角速度向量

Figure BDA0002463896770000046
三只加速度计测量采煤机机身三个轴向的加速度向量
Figure BDA0002463896770000051
2) The three gyroscopes of the fuselage strapdown inertial navigation device 1 measure the angular velocity vectors of the three axes of the shearer fuselage
Figure BDA0002463896770000046
Three accelerometers measure the acceleration vectors in three axial directions of the shearer fuselage
Figure BDA0002463896770000051

3)机身捷联惯导装置1首先更新采煤机机身b1坐标系相对于导航坐标系n的角速度向量

Figure BDA0002463896770000052
3) The fuselage strapdown inertial navigation device 1 first updates the angular velocity vector of the shearer fuselage b1 coordinate system relative to the navigation coordinate system n
Figure BDA0002463896770000052

Figure BDA0002463896770000053
Figure BDA0002463896770000053

其中,

Figure BDA0002463896770000054
表示机身捷联惯导装置1纬度,
Figure BDA0002463896770000055
为机身捷联惯导装置1姿态矩阵,ωie是地球自转角速度,R为地球半径,
Figure BDA0002463896770000056
Figure BDA0002463896770000057
表示机身捷联惯导装置1东向速度和北向速度;in,
Figure BDA0002463896770000054
Indicates the 1 latitude of the fuselage strapdown inertial navigation device,
Figure BDA0002463896770000055
is the attitude matrix of the fuselage SINS1, ω ie is the angular velocity of the earth's rotation, R is the earth's radius,
Figure BDA0002463896770000056
and
Figure BDA0002463896770000057
Indicates the east speed and north speed of the fuselage strapdown inertial navigation device 1;

4)更新四元数q1=[q0 q1 q2 q3]T 4) Update quaternion q 1 =[q 0 q 1 q 2 q 3 ] T

Figure BDA0002463896770000058
Figure BDA0002463896770000058

其中,

Figure BDA0002463896770000059
T表示姿态计算周期;in,
Figure BDA0002463896770000059
T represents the attitude calculation period;

5)更新机身捷联惯导装置1相对于导航坐标系的姿态矩阵

Figure BDA00024638967700000510
5) Update the attitude matrix of the fuselage strapdown inertial navigation device 1 relative to the navigation coordinate system
Figure BDA00024638967700000510

Figure BDA00024638967700000511
Figure BDA00024638967700000511

其中q1=[q0 q1 q2 q3]T,令

Figure BDA00024638967700000512
其中i,j=1,2,3where q 1 =[q 0 q 1 q 2 q 3 ] T , let
Figure BDA00024638967700000512
where i,j=1,2,3

6)根据

Figure BDA00024638967700000513
计算机身姿态角:6) According to
Figure BDA00024638967700000513
Computer body attitude angle:

俯仰角θ1=sin-1(C12)Pitch angle θ 1 =sin −1 (C 12 )

航向角

Figure BDA00024638967700000514
Heading
Figure BDA00024638967700000514

滚动角

Figure BDA0002463896770000061
roll angle
Figure BDA0002463896770000061

7)利用

Figure BDA0002463896770000062
将加速度计输出的加速度向量
Figure BDA0002463896770000063
投影到导航坐标系
Figure BDA00024638967700000614
7) Utilize
Figure BDA0002463896770000062
Convert the acceleration vector output by the accelerometer
Figure BDA0002463896770000063
Project to Navigation Coordinate System
Figure BDA00024638967700000614

Figure BDA0002463896770000064
Figure BDA0002463896770000064

8)更新机身捷联惯导装置1的速度

Figure BDA0002463896770000065
位置
Figure BDA0002463896770000066
和经度λ1、纬度
Figure BDA0002463896770000067
高度h1,8) Update the speed of the fuselage SINS 1
Figure BDA0002463896770000065
Location
Figure BDA0002463896770000066
and longitude λ 1 , latitude
Figure BDA0002463896770000067
height h 1 ,

Figure BDA0002463896770000068
Figure BDA0002463896770000068

Figure BDA0002463896770000069
Figure BDA0002463896770000069

Figure BDA00024638967700000610
Figure BDA00024638967700000610

步骤三、计算采煤机摇臂信息Step 3. Calculate the rocker arm information of the shearer

采煤机摇臂捷联惯导装置3利用陀螺和加速度计测量信息,解算采煤机摇臂姿态俯仰角θ2、航向角ψ2和滚动角γ2,计算摇臂的速度

Figure BDA00024638967700000611
位置
Figure BDA00024638967700000612
和经度λ2、纬度
Figure BDA00024638967700000613
高度h2;The shearer rocker arm strapdown inertial navigation device 3 uses the gyro and accelerometer measurement information to solve the shearer rocker arm attitude pitch angle θ 2 , heading angle ψ 2 and roll angle γ 2 , and calculate the rocker arm speed
Figure BDA00024638967700000611
Location
Figure BDA00024638967700000612
and longitude λ 2 , latitude
Figure BDA00024638967700000613
height h 2 ;

摇臂捷联惯导装置3接收初始装订信息:初始位置

Figure BDA0002463896770000071
和初始姿态
Figure BDA0002463896770000072
其中上标2表示摇臂捷联惯导,N,U,E分别表示导航坐标系下北向、天向和东向位置;同时由于采煤机初始处于停止状态,初始速度V2=[0 0 0]T;Rocker arm SINS 3 receives initial binding information: initial position
Figure BDA0002463896770000071
and initial posture
Figure BDA0002463896770000072
The superscript 2 represents the rocker arm strapdown inertial navigation, and N, U, and E represent the north, sky and east positions in the navigation coordinate system, respectively. 0] T ;

摇臂捷联惯导装置3与机身捷联惯导装置1采用相同的位置和姿态解算方法,将步骤二中1)至8)中变量上标1更换为2,即为摇臂捷联惯导装置3的位置和姿态解算方法;The rocker arm strapdown inertial navigation device 3 and the fuselage strapdown inertial navigation device 1 adopt the same position and attitude calculation method, and the variable superscript 1 in step 2 in 1) to 8) is replaced by 2, that is, the rocker arm The position and attitude calculation method of the coupled inertial navigation device 3;

步骤四、计算采煤机机身与摇臂的相对距离Step 4. Calculate the relative distance between the shearer body and the rocker arm

摇臂无线传感移动节点4相对于机身无线传感锚节点2测距:The rocker arm wireless sensor mobile node 4 measures the distance relative to the fuselage wireless sensor anchor node 2:

摇臂无线传感移动节点4发射无线信号,机身无线传感锚节点2接收移动节点的无线信号,利用RSSI算法进行测距,基于RSSI的测距方法将接收到的信号强度转化为无线传感移动节点4与无线传感锚节点2之间的距离,使用公知的对数-常数无线信号传播模型,The rocker arm wireless sensor mobile node 4 transmits wireless signals, the fuselage wireless sensor anchor node 2 receives the wireless signal of the mobile node, and uses the RSSI algorithm for ranging, and the RSSI-based ranging method converts the received signal strength into wireless transmission. sense the distance between the mobile node 4 and the wireless sensing anchor node 2, using the well-known log-constant wireless signal propagation model,

S=A-10m lg(d) (8)S=A-10m lg(d) (8)

其中,S为信号强度,d为距离,A为接收端距离发射端为1m时接收到的无线信号强度,m为路径损耗,A和m为已知参数,根据S计算距离d;Among them, S is the signal strength, d is the distance, A is the received wireless signal strength when the receiving end is 1m away from the transmitting end, m is the path loss, A and m are known parameters, and the distance d is calculated according to S;

根据公式(8),根据无线传感锚节点2接收到的信号强度S,可得摇臂相对于机身的距离d;According to formula (8), according to the signal strength S received by the wireless sensor anchor node 2, the distance d of the rocker arm relative to the fuselage can be obtained;

在无线传感移动节点4发射无线信号时,摇臂捷联惯导装置3记录该时刻惯导解算的位置信息,并发送给机身捷联惯导装置1;机身无线传感锚节点2接收无线信号时,接收摇臂捷联惯导装置3发送的测距时刻位置信息,同时记录机身捷联惯导装置1解算的机身位置姿态信息及摇臂位置姿态信息。When the wireless sensor mobile node 4 transmits a wireless signal, the rocker arm SINS device 3 records the position information calculated by the inertial navigation at that moment, and sends it to the fuselage SINS device 1; the fuselage wireless sensor anchor node 2. When receiving wireless signals, receive the ranging moment position information sent by the rocker arm strapdown inertial navigation device 3, and simultaneously record the body position and attitude information and the rocker arm position and attitude information calculated by the fuselage strapdown inertial navigation device 1.

步骤五、利用步骤四中在无线测距信号发射时刻记录的机身位置姿态信息、摇臂位置姿态信息、机身与摇臂的相对距离,进行基于卡尔曼滤波器的协同导航滤波解算,估计摇臂捷联惯导装置3的位置、速度和姿态误差;Step 5: Use the fuselage position and attitude information, the rocker arm position and attitude information, and the relative distance between the fuselage and the rocker arm recorded at the time of transmitting the wireless ranging signal in step 4, to carry out the Kalman filter-based collaborative navigation filtering solution, Estimate the position, velocity and attitude errors of the rocker arm SINS 3;

1)协同导航滤波器的状态变量为

Figure BDA0002463896770000081
其中
Figure BDA0002463896770000082
表示摇臂捷联惯导装置3位置误差;
Figure BDA0002463896770000083
表示摇臂捷联惯导装置3速度误差;φ=[φN φU φE]T表示摇臂捷联惯导装置3姿态误差;1) The state variable of the collaborative navigation filter is
Figure BDA0002463896770000081
in
Figure BDA0002463896770000082
Indicates the position error of the rocker arm strapdown inertial navigation device 3;
Figure BDA0002463896770000083
Represents the speed error of the rocker arm SINS device 3; φ=[φ N φ U φ E ] T represents the attitude error of the rocker arm SINS device 3;

2)协同导航滤波状态方程由摇臂捷联惯导装置3的速度误差方程、位置误差方程和姿态误差方程组成:2) The cooperative navigation filtering state equation is composed of the velocity error equation, position error equation and attitude error equation of the rocker arm SINS device 3:

速度误差方程:Velocity Error Equation:

Figure BDA0002463896770000084
Figure BDA0002463896770000084

其中,in,

Figure BDA0002463896770000085
Figure BDA0002463896770000086
Figure BDA0002463896770000085
Figure BDA0002463896770000086

位置误差方程Position Error Equation

Figure BDA0002463896770000087
Figure BDA0002463896770000087

姿态误差方程Attitude Error Equation

Figure BDA0002463896770000091
Figure BDA0002463896770000091

将方程(9)、(10)、(11)进行离散化并改写为协同导航滤波状态方程形式,其中tk表示滤波时间点;The equations (9), (10) and (11) are discretized and rewritten into the form of cooperative navigation filtering state equation, where t k represents the filtering time point;

X(tk)=F(tk-1)X(tk-1) (12)X(t k )=F(t k-1 )X(t k-1 ) (12)

3)协同导航滤波量测方程如下:3) The collaborative navigation filtering measurement equation is as follows:

Figure BDA0002463896770000092
Figure BDA0002463896770000092

4)修正机身捷联惯装置1导同步误差4) Correct the 1-lead synchronization error of the fuselage strapdown inertial device

机身捷联惯导装置1定位数据更新率高于无线测距的定位数据更新率,可采用线性插值方法,在摇臂无线传感移动节点4发射无线信号时,机身捷联惯导装置1记录时间间隔Δt,Δt<ΔT,采用插值方法,计算用于协同导航的机身捷联惯导装置1的位置P1The update rate of the positioning data of the fuselage strapdown inertial navigation device 1 is higher than that of the wireless ranging, and the linear interpolation method can be used. When the rocker arm wireless sensor mobile node 4 transmits wireless signals, the fuselage strapdown inertial navigation 1. Record the time interval Δt, Δt<ΔT, and use the interpolation method to calculate the position P 1 of the fuselage strapdown inertial navigation device 1 used for cooperative navigation:

Figure BDA0002463896770000093
Figure BDA0002463896770000093

5)计算量测量Z5) Calculate the amount of measurement Z

在无线测距信号发射时刻,机身无线传感锚节点2位置(即机身捷联惯导装置位置)为

Figure BDA0002463896770000094
在摇臂无线传感移动节点4位置即摇臂捷联惯导装置3位置为
Figure BDA0002463896770000095
根据几何关系有:At the moment when the wireless ranging signal is transmitted, the position of the fuselage wireless sensor anchor node 2 (that is, the position of the fuselage strapdown inertial navigation device) is
Figure BDA0002463896770000094
At the position of the rocker arm wireless sensor mobile node 4, that is, the position of the rocker arm strapdown inertial navigation device 3 is:
Figure BDA0002463896770000095
According to the geometric relationship there are:

Figure BDA0002463896770000096
Figure BDA0002463896770000096

其中,Y表示利用摇臂捷联惯导装置3位置与机身捷联惯导装置1位置计算的相对距离。那么协同导航滤波的量测量Z计算如下:Among them, Y represents the relative distance calculated by using the position of the rocker arm SINS device 3 and the position of the fuselage SINS device 1 . Then the quantity measurement Z of the collaborative navigation filtering is calculated as follows:

Z=Y-dZ=Y-d

Figure BDA0002463896770000097
Figure BDA0002463896770000097

6)机身捷联惯导装置1利用公知的卡尔曼滤波算法,可得协同导航滤波状态量估计值

Figure BDA0002463896770000101
6) The fuselage strapdown inertial navigation device 1 uses the well-known Kalman filter algorithm to obtain the estimated value of the cooperative navigation filter state quantity
Figure BDA0002463896770000101

步骤六、利用协同导航滤波状态估计值

Figure BDA0002463896770000102
修正摇臂捷联惯导装置3的位置、速度和姿态误差,得到误差修正后的摇臂的位置、速度和姿态:Step 6. Use collaborative navigation to filter state estimates
Figure BDA0002463896770000102
Correct the position, speed and attitude errors of the rocker arm strapdown inertial navigation device 3, and obtain the position, speed and attitude of the rocker arm after error correction:

机身捷联惯导装置1将

Figure BDA0002463896770000103
发送给摇臂捷联惯导装置3,对其误差进行闭环点修正:The fuselage strapdown inertial navigation device 1 will be
Figure BDA0002463896770000103
Send it to the rocker arm strapdown inertial navigation device 3, and perform closed-loop point correction for its error:

P2(tk)=P2(tk)-ΔP(tk) (17)P 2 (t k )=P 2 (t k )-ΔP(t k ) (17)

V2(tk)=V2(tk)-ΔV(tk) (18)V 2 (t k )=V 2 (t k )-ΔV(t k ) (18)

Figure BDA0002463896770000104
Figure BDA0002463896770000104

得到误差修正后的摇臂捷联惯导装置3位置、速度和姿态矩阵,对修正误差后的

Figure BDA0002463896770000105
按照步骤二的6),将变量上标1更换为2,计算误差修正后的摇臂姿态角:俯仰角θ2、航向角ψ2和滚动角γ2。The position, velocity and attitude matrix of rocker arm SINS device 3 after error correction is obtained.
Figure BDA0002463896770000105
According to 6) of step 2, replace the variable superscript 1 with 2, and calculate the attitude angle of the rocker arm after error correction: pitch angle θ 2 , heading angle ψ 2 and roll angle γ 2 .

本发明的有益效果是:提供了一种采煤机机身5与摇臂6协同定位方法,高精度的机身捷联惯导装置1解算机身位置和姿态,摇臂捷联惯导装置3解算摇臂位置和姿态,通过机身无线传感锚节点2与摇臂无线传感移动节点4之间的无线测距,利用协同导航滤波方法,估计并修正摇臂位置和姿态误差,实现采煤机机身5和摇臂6的高精度定位定姿,使其满足综采工作面自动化、远程自动控制的需求,具有定位、定姿精度高、综采工作效率高及远程自动控制的优点。The beneficial effects of the present invention are as follows: a method for co-locating the fuselage 5 and the rocker arm 6 of the shearer is provided; The device 3 calculates the position and attitude of the rocker arm, and estimates and corrects the position and attitude errors of the rocker arm through the wireless ranging between the fuselage wireless sensor anchor node 2 and the rocker arm wireless sensor mobile node 4, using the collaborative navigation filtering method , to achieve high-precision positioning and attitude determination of the shearer fuselage 5 and rocker arm 6, so that it can meet the needs of automatic and remote automatic control of fully mechanized mining face, with high positioning and attitude accuracy, high fully mechanized mining work efficiency and remote automatic control. Advantages of control.

附图说明Description of drawings

图1为采煤机机身与摇臂协同定位示意图。Figure 1 is a schematic diagram of the cooperative positioning of the shearer body and the rocker arm.

图2为机身(摇臂)坐标系、导航坐标系、姿态示意图。Figure 2 is a schematic diagram of the fuselage (rocker) coordinate system, the navigation coordinate system and the attitude.

图3为采煤机机身与摇臂协同定位工作流程图。Figure 3 is a flow chart of the co-location of the shearer body and the rocker arm.

图4为采煤机惯性/无线传感器网络组合定位信息时序图。Figure 4 is a time sequence diagram of shearer inertial/wireless sensor network combined positioning information.

图中:1、机身捷联惯导装置;2、机身无线传感锚节点;3、摇臂捷联惯导装置;4、摇臂无线传感移动节点;5、机身;6、摇臂。In the figure: 1. Airframe strapdown inertial navigation device; 2. Airframe wireless sensor anchor node; 3. Rocker arm strapdown inertial navigation device; 4. Rocker arm wireless sensor mobile node; 5. Airframe; 6. rocker arm.

具体实施方式Detailed ways

下面结合附图对本发明做进一步详细说明。The present invention will be further described in detail below in conjunction with the accompanying drawings.

参见图1,一种采煤机机身与摇臂协同定位装置,包括在采煤机机身上安装的机身捷联惯导装置1和无线传感锚节点2;以及在采煤机割煤一侧的摇臂6上安装的摇臂捷联惯导装置3和无线传感移动节点4;其中无线传感移动节点4发射无线信号,无线传感锚节点2接收无线信号。Referring to Fig. 1, a shearer fuselage and rocker arm co-positioning device, comprising a fuselage strapdown inertial navigation device 1 and a wireless sensor anchor node 2 installed on the shearer fuselage; The rocker arm strapdown inertial navigation device 3 and the wireless sensor mobile node 4 are installed on the rocker arm 6 on the coal side; the wireless sensor mobile node 4 transmits wireless signals, and the wireless sensor anchor node 2 receives wireless signals.

所述机身捷联惯导装置1包括三只陀螺和三只加速度计。The fuselage strapdown inertial navigation device 1 includes three gyroscopes and three accelerometers.

基于上述的采煤机机身与摇臂协同定位装置的定位方法,具体实施方式如下:Based on the above-mentioned positioning method of the shearer fuselage and the rocker arm cooperative positioning device, the specific implementation is as follows:

步骤一、定义坐标系Step 1. Define the coordinate system

1)定义采煤机机身坐标系b1:Ob1Xb1Yb1Zb1:坐标系原点Ob1固连在机身捷联惯导装置1中心,Xb1轴正向由采煤机指向煤壁,Yb1轴垂直于Xb1轴向上,Zb1轴与Xb1轴、Yb1轴构成右手坐标系,构成前上右坐标系,下标1表示该坐标系为机身坐标系;1) Define the shearer fuselage coordinate system b1: O b1 X b1 Y b1 Z b1 : the origin of the coordinate system O b1 is fixedly connected to the center of the fuselage strapdown inertial navigation device 1, and the positive X b1 axis points from the shearer to the coal wall, the Y b1 axis is perpendicular to the X b1 axis upward, the Z b1 axis, the X b1 axis, and the Y b1 axis form a right-hand coordinate system, which constitutes the front upper right coordinate system, and the subscript 1 indicates that the coordinate system is the fuselage coordinate system;

2)定义采煤机摇臂坐标系b2:Ob2Xb2Yb2Zb2;坐标系原点Ob2固连在摇臂捷联惯导装置中心,Xb2轴正向由摇臂指向煤壁,Yb2轴垂直于Xb2轴向上,Zb2轴与Xb2轴、Yb2轴构成右手坐标系,构成前上右坐标系,下标2表示该坐标系为摇臂坐标系;2) Define the shearer rocker arm coordinate system b2: O b2 X b2 Y b2 Z b2 ; the origin of the coordinate system O b2 is fixedly connected to the center of the rocker arm strapdown inertial navigation device, and the positive X b2 axis is directed by the rocker arm to the coal wall, The Y b2 axis is perpendicular to the X b2 axis upward, the Z b2 axis, the X b2 axis, and the Y b2 axis form a right-hand coordinate system, which constitutes the front upper right coordinate system, and the subscript 2 indicates that the coordinate system is a rocker arm coordinate system;

3)定义导航坐标系n:OnXnYnZn:北天东地理坐标系,Xn轴指向地理北向,Yn轴指向天向,Zn轴指向地理东向;3) Define the navigation coordinate system n: O n X n Y n Z n : the north celestial east geographic coordinate system, the X n axis points to the geographic north direction, the Y n axis points to the sky direction, and the Z n axis points to the geographic east direction;

4)导航坐标系n经过三次旋转后与采煤机机身坐标系b1重合,三次旋转的角度即为机身航向角ψ1、俯仰角θ1和滚动角γ1,如图2所示;类似地,导航坐标系n经过三次旋转后与采煤机摇臂坐标系b2重合,三次旋转的角度即为摇臂航向角ψ2、俯仰角θ2和滚动角γ24) The navigation coordinate system n coincides with the shearer fuselage coordinate system b1 after three rotations, and the angles of the three rotations are the fuselage heading angle ψ 1 , pitch angle θ 1 and roll angle γ 1 , as shown in Figure 2; Similarly, the navigation coordinate system n coincides with the shearer rocker arm coordinate system b2 after three rotations, and the angles of the three rotations are the rocker arm heading angle ψ 2 , pitch angle θ 2 and roll angle γ 2 ;

步骤二、计算采煤机机身信息Step 2. Calculate the information of the shearer fuselage

采煤机机身捷联惯导装置1利用陀螺和加速度计测量信息,解算采煤机机身姿态俯仰角θ1、航向角ψ1和滚动角γ1,计算机身的速度

Figure BDA0002463896770000121
位置
Figure BDA0002463896770000122
和经度λ1、纬度
Figure BDA0002463896770000123
高度h1;The shearer fuselage strapdown inertial navigation device 1 uses the gyro and accelerometer measurement information to calculate the shearer fuselage attitude pitch angle θ 1 , heading angle ψ 1 and roll angle γ 1 , and calculate the speed of the body
Figure BDA0002463896770000121
Location
Figure BDA0002463896770000122
and longitude λ 1 , latitude
Figure BDA0002463896770000123
height h 1 ;

1)机身捷联惯导装置1接收初始装订信息:初始位置

Figure BDA0002463896770000124
和初始姿态
Figure BDA0002463896770000125
其中上标1表示机身捷联惯导,N,U,E分别表示导航坐标系下北向、天向和东向位置;同时由于采煤机初始处于停止状态,初始速度V1=[0 0 0]T;1) Body SINS 1 receives initial binding information: initial position
Figure BDA0002463896770000124
and initial posture
Figure BDA0002463896770000125
The superscript 1 represents the fuselage strapdown inertial navigation, and N, U, and E represent the north, sky and east positions in the navigation coordinate system, respectively. 0] T ;

2)机身捷联惯导装置1的三只陀螺测量采煤机机身三个轴的角速度向量

Figure BDA0002463896770000126
三只加速度计测量采煤机机身三个轴向的加速度向量
Figure BDA0002463896770000127
2) The three gyroscopes of the fuselage strapdown inertial navigation device 1 measure the angular velocity vectors of the three axes of the shearer fuselage
Figure BDA0002463896770000126
Three accelerometers measure the acceleration vectors in three axial directions of the shearer fuselage
Figure BDA0002463896770000127

3)机身捷联惯导装置1首先更新采煤机机身b1坐标系相对于导航坐标系n的角速度向量

Figure BDA0002463896770000131
3) The fuselage strapdown inertial navigation device 1 first updates the angular velocity vector of the shearer fuselage b1 coordinate system relative to the navigation coordinate system n
Figure BDA0002463896770000131

Figure BDA0002463896770000132
Figure BDA0002463896770000132

其中,

Figure BDA0002463896770000133
表示机身捷联惯导装置1纬度,
Figure BDA0002463896770000134
为机身捷联惯导装置1姿态矩阵,ωie是地球自转角速度,R为地球半径,
Figure BDA0002463896770000135
Figure BDA0002463896770000136
表示机身捷联惯导装置1东向速度和北向速度;in,
Figure BDA0002463896770000133
Indicates the 1 latitude of the fuselage strapdown inertial navigation device,
Figure BDA0002463896770000134
is the attitude matrix of the fuselage SINS1, ω ie is the angular velocity of the earth's rotation, R is the earth's radius,
Figure BDA0002463896770000135
and
Figure BDA0002463896770000136
Indicates the east speed and north speed of the fuselage strapdown inertial navigation device 1;

4)更新四元数q1=[q0 q1 q2 q3]T 4) Update quaternion q 1 =[q 0 q 1 q 2 q 3 ] T

Figure BDA0002463896770000137
Figure BDA0002463896770000137

其中,

Figure BDA0002463896770000138
T表示姿态计算周期;in,
Figure BDA0002463896770000138
T represents the attitude calculation period;

5)更新机身捷联惯导装置1相对于导航坐标系的姿态矩阵

Figure BDA0002463896770000139
5) Update the attitude matrix of the fuselage strapdown inertial navigation device 1 relative to the navigation coordinate system
Figure BDA0002463896770000139

Figure BDA00024638967700001310
Figure BDA00024638967700001310

其中q1=[q0 q1 q2 q3]T,令

Figure BDA00024638967700001311
其中i,j=1,2,3where q 1 =[q 0 q 1 q 2 q 3 ] T , let
Figure BDA00024638967700001311
where i,j=1,2,3

6)根据

Figure BDA00024638967700001312
计算机身姿态角:6) According to
Figure BDA00024638967700001312
Computer body attitude angle:

俯仰角θ1=sin-1(C12)Pitch angle θ 1 =sin −1 (C 12 )

航向角

Figure BDA00024638967700001313
Heading
Figure BDA00024638967700001313

滚动角

Figure BDA00024638967700001314
roll angle
Figure BDA00024638967700001314

7)利用

Figure BDA00024638967700001315
将加速度计输出的加速度向量
Figure BDA0002463896770000141
投影到导航坐标系
Figure BDA0002463896770000142
7) Utilize
Figure BDA00024638967700001315
Convert the acceleration vector output by the accelerometer
Figure BDA0002463896770000141
Project to Navigation Coordinate System
Figure BDA0002463896770000142

Figure BDA0002463896770000143
Figure BDA0002463896770000143

8)更新机身捷联惯导装置1的速度

Figure BDA0002463896770000144
位置
Figure BDA0002463896770000145
和经度λ1、纬度
Figure BDA0002463896770000146
高度h1,8) Update the speed of the fuselage SINS 1
Figure BDA0002463896770000144
Location
Figure BDA0002463896770000145
and longitude λ 1 , latitude
Figure BDA0002463896770000146
height h 1 ,

Figure BDA0002463896770000147
Figure BDA0002463896770000147

Figure BDA0002463896770000148
Figure BDA0002463896770000148

Figure BDA0002463896770000149
Figure BDA0002463896770000149

步骤三、计算采煤机摇臂信息Step 3. Calculate the rocker arm information of the shearer

采煤机摇臂捷联惯导装置3利用陀螺和加速度计测量信息,解算采煤机摇臂姿态俯仰角θ2、航向角ψ2和滚动角γ2,计算摇臂的速度

Figure BDA00024638967700001410
位置
Figure BDA00024638967700001411
和经度λ2、纬度
Figure BDA00024638967700001412
高度h2;The shearer rocker arm strapdown inertial navigation device 3 uses the gyro and accelerometer measurement information to solve the shearer rocker arm attitude pitch angle θ 2 , heading angle ψ 2 and roll angle γ 2 , and calculate the rocker arm speed
Figure BDA00024638967700001410
Location
Figure BDA00024638967700001411
and longitude λ 2 , latitude
Figure BDA00024638967700001412
height h 2 ;

摇臂捷联惯导装置3接收初始装订信息:初始位置

Figure BDA00024638967700001413
和初始姿态
Figure BDA00024638967700001414
其中上标2表示摇臂捷联惯导,N,U,E分别表示导航坐标系下北向、天向和东向位置;同时由于采煤机初始处于停止状态,初始速度V2=[0 0 0]TRocker arm SINS 3 receives initial binding information: initial position
Figure BDA00024638967700001413
and initial posture
Figure BDA00024638967700001414
The superscript 2 represents the rocker arm strapdown inertial navigation, and N, U, and E represent the north, sky and east positions in the navigation coordinate system, respectively. 0] T ;

摇臂捷联惯导装置3与机身捷联惯导装置1采用相同的位置和姿态解算方法,将步骤二中1)至8)中变量上标1更换为2,即为摇臂捷联惯导装置3的位置和姿态解算方法;The rocker arm strapdown inertial navigation device 3 and the fuselage strapdown inertial navigation device 1 adopt the same position and attitude calculation method, and the variable superscript 1 in step 2 in 1) to 8) is replaced by 2, that is, the rocker arm The position and attitude calculation method of the coupled inertial navigation device 3;

步骤四、计算采煤机机身与摇臂的相对距离Step 4. Calculate the relative distance between the shearer body and the rocker arm

摇臂无线传感移动节点4相对于机身无线传感锚节点2测距:The rocker arm wireless sensor mobile node 4 measures the distance relative to the fuselage wireless sensor anchor node 2:

摇臂无线传感移动节点4发射无线信号,机身无线传感锚节点2接收移动节点的无线信号,利用RSSI算法进行测距,基于RSSI的测距方法将接收到的信号强度转化为无线传感移动节点4与无线传感锚节点2之间的距离,使用公知的对数-常数无线信号传播模型,The rocker arm wireless sensor mobile node 4 transmits wireless signals, the fuselage wireless sensor anchor node 2 receives the wireless signal of the mobile node, and uses the RSSI algorithm for ranging, and the RSSI-based ranging method converts the received signal strength into wireless transmission. sense the distance between the mobile node 4 and the wireless sensing anchor node 2, using the well-known log-constant wireless signal propagation model,

S=A-10m lg(d) (8)S=A-10m lg(d) (8)

其中,S为信号强度,d为距离,A为接收端距离发射端为1m时接收到的无线信号强度,m为路径损耗,A和m为已知参数,根据S计算距离d;Among them, S is the signal strength, d is the distance, A is the received wireless signal strength when the receiving end is 1m away from the transmitting end, m is the path loss, A and m are known parameters, and the distance d is calculated according to S;

根据公式(8),根据无线传感锚节点2接收到的信号强度S,可得摇臂相对于机身的距离d;According to formula (8), according to the signal strength S received by the wireless sensor anchor node 2, the distance d of the rocker arm relative to the fuselage can be obtained;

在无线传感移动节点4发射无线信号时,摇臂捷联惯导装置3记录该时刻惯导解算的位置信息,并发送给机身捷联惯导装置1;机身无线传感锚节点2接收无线信号时,接收摇臂捷联惯导装置3发送的测距时刻位置信息,同时记录机身捷联惯导装置1解算的机身位置姿态信息及摇臂位置姿态信息;When the wireless sensor mobile node 4 transmits a wireless signal, the rocker arm SINS device 3 records the position information calculated by the inertial navigation at that moment, and sends it to the fuselage SINS device 1; the fuselage wireless sensor anchor node 2. When receiving the wireless signal, receive the position information at the time of ranging sent by the rocker arm strapdown inertial navigation device 3, and simultaneously record the body position and attitude information and the rocker arm position and attitude information calculated by the fuselage strapdown inertial navigation device 1;

步骤五、利用步骤四中在无线测距信号发射时刻记录的机身位置姿态信息、摇臂位置姿态信息、机身与摇臂的相对距离,进行基于卡尔曼滤波器的协同导航滤波解算,估计摇臂捷联惯导装置3的位置、速度和姿态误差;Step 5: Use the fuselage position and attitude information, the rocker arm position and attitude information, and the relative distance between the fuselage and the rocker arm recorded at the time of transmitting the wireless ranging signal in step 4, to carry out the Kalman filter-based collaborative navigation filtering solution, Estimate the position, velocity and attitude errors of the rocker arm SINS 3;

1)协同导航滤波器的状态变量为

Figure BDA0002463896770000161
其中ΔPI=[ΔPN ΔPUΔPE]T表示摇臂捷联惯导装置3位置误差;ΔV=[ΔVN ΔVU ΔVE]T表示摇臂捷联惯导装置3速度误差;φ=[φN φU φE]T表示摇臂捷联惯导装置3姿态误差;1) The state variable of the collaborative navigation filter is
Figure BDA0002463896770000161
Wherein ΔP I =[ΔP N ΔP U ΔP E ] T represents the position error of the rocker arm SINS 3; ΔV=[ΔV N ΔV U ΔV E ] T represents the speed error of the rocker arm SINS 3; φ= [φ N φ U φ E ] T represents the attitude error of the rocker arm strapdown inertial navigation device 3;

2)协同导航滤波状态方程由摇臂捷联惯导装置3的速度误差方程、位置误差方程和姿态误差方程组成:2) The cooperative navigation filtering state equation is composed of the velocity error equation, position error equation and attitude error equation of the rocker arm SINS device 3:

速度误差方程:Velocity Error Equation:

Figure BDA0002463896770000162
Figure BDA0002463896770000162

其中,in,

Figure BDA0002463896770000163
Figure BDA0002463896770000164
Figure BDA0002463896770000163
Figure BDA0002463896770000164

位置误差方程Position Error Equation

Figure BDA0002463896770000165
Figure BDA0002463896770000165

姿态误差方程Attitude Error Equation

Figure BDA0002463896770000166
Figure BDA0002463896770000166

将方程(9)、(10)、(11)进行离散化并改写为协同导航滤波状态方程形式,其中tk表示滤波时间点;The equations (9), (10) and (11) are discretized and rewritten into the form of cooperative navigation filtering state equation, where t k represents the filtering time point;

X(tk)=F(tk-1)X(tk-1) (12)X(t k )=F(t k-1 )X(t k-1 ) (12)

3)协同导航滤波滤波量测方程如下:3) The collaborative navigation filtering filtering measurement equation is as follows:

Figure BDA0002463896770000171
Figure BDA0002463896770000171

4)修正机身捷联惯装置1导同步误差4) Correct the 1-lead synchronization error of the fuselage strapdown inertial device

机身捷联惯导装置1定位数据更新率高于无线测距的定位数据更新率,采用如图4所示的插值方法,在摇臂无线传感移动节点4发射无线信号时,机身捷联惯导装置1记录时间间隔Δt,Δt<ΔT,采用插值方法,计算用于协同导航的机身捷联惯导装置1的位置P1The update rate of the positioning data of the fuselage strapdown inertial navigation device 1 is higher than the update rate of the positioning data of the wireless ranging, and the interpolation method shown in Figure 4 is used. The strap-down inertial navigation device 1 records the time interval Δt, Δt<ΔT, and uses the interpolation method to calculate the position P 1 of the fuselage strap-down inertial navigation device 1 for cooperative navigation:

Figure BDA0002463896770000172
Figure BDA0002463896770000172

5)计算量测量Z5) Calculate the amount of measurement Z

在无线测距信号发射时刻,机身无线传感锚节点2位置(即机身捷联惯导装置位置)为

Figure BDA0002463896770000173
在摇臂无线传感移动节点4位置即摇臂捷联惯导装置3位置为
Figure BDA0002463896770000174
根据几何关系有:At the moment when the wireless ranging signal is transmitted, the position of the fuselage wireless sensor anchor node 2 (that is, the position of the fuselage strapdown inertial navigation device) is
Figure BDA0002463896770000173
At the position of the rocker arm wireless sensor mobile node 4, that is, the position of the rocker arm strapdown inertial navigation device 3 is:
Figure BDA0002463896770000174
According to the geometric relationship there are:

Figure BDA0002463896770000175
Figure BDA0002463896770000175

其中,Y表示利用摇臂捷联惯导装置3位置与机身捷联惯导装置1位置计算的相对距离。那么协同导航滤波的量测量Z计算如下:Among them, Y represents the relative distance calculated by using the position of the rocker arm SINS device 3 and the position of the fuselage SINS device 1 . Then the quantity measurement Z of the collaborative navigation filtering is calculated as follows:

Z=Y-dZ=Y-d

Figure BDA0002463896770000176
Figure BDA0002463896770000176

6)机身捷联惯导装置1利用公知的卡尔曼滤波算法,可得协同导航滤波状态量估计值

Figure BDA0002463896770000177
6) The fuselage strapdown inertial navigation device 1 uses the well-known Kalman filter algorithm to obtain the estimated value of the cooperative navigation filter state quantity
Figure BDA0002463896770000177

步骤六、利用协同导航滤波状态估计值

Figure BDA0002463896770000178
修正摇臂捷联惯导装置3的位置、速度和姿态误差,得到误差修正后的摇臂的位置、速度和姿态:Step 6. Use collaborative navigation to filter state estimates
Figure BDA0002463896770000178
Correct the position, speed and attitude errors of the rocker arm strapdown inertial navigation device 3, and obtain the position, speed and attitude of the rocker arm after error correction:

机身捷联惯导装置1将

Figure BDA0002463896770000181
发送给摇臂捷联惯导装置3,对其误差进行闭环点修正:The fuselage strapdown inertial navigation device 1 will be
Figure BDA0002463896770000181
Send it to the rocker arm strapdown inertial navigation device 3, and perform closed-loop point correction for its error:

P2(tk)=P2(tk)-ΔP(tk) (17)P 2 (t k )=P 2 (t k )-ΔP(t k ) (17)

V2(tk)=V2(tk)-ΔV(tk) (18)V 2 (t k )=V 2 (t k )-ΔV(t k ) (18)

Figure BDA0002463896770000182
Figure BDA0002463896770000182

得到误差修正后的摇臂捷联惯导装置3位置、速度和姿态矩阵,对修正误差后的

Figure BDA0002463896770000183
按照步骤二的6),将变量上标1更换为2,计算误差修正后的摇臂姿态角:俯仰角θ2、航向角ψ2和滚动角γ2。The position, velocity and attitude matrix of rocker arm SINS device 3 after error correction is obtained.
Figure BDA0002463896770000183
According to 6) of step 2, replace the variable superscript 1 with 2, and calculate the attitude angle of the rocker arm after error correction: pitch angle θ 2 , heading angle ψ 2 and roll angle γ 2 .

参见图3,本发明的工作原理为:Referring to Figure 3, the working principle of the present invention is:

在采煤机的采煤作业中,采煤机机身5沿轨道移动,而摇臂6进行截割煤壁运动;摇臂捷联惯导装置3测量采煤机摇臂的角速度和加速度信息,利用捷联导航算法进行定位定姿;摇臂无线传感移动节点4发射无线信号,同时摇臂捷联惯导装置3将该时刻的惯性信息发送给机身捷联惯导装置1;机身捷联惯导装置1实时测量采煤机机身的角速度和加速度信息,利用捷联导航算法进行定位定姿;机身无线传感锚节点2接收摇臂无线传感移动节点4的无线信号,利用RSSI算法进行测距解算;机身捷联惯导装置1解算无线信号到达时刻的机身惯性信息,接收摇臂捷联惯导装置3的惯性信息,与机身无线传感锚节点2测距定位结果进行协同导航滤波,估计摇臂位置、姿态误差;根据误差估计结果,将修正后的摇臂位置、姿态信息发送给摇臂捷联惯导装置,修正摇臂的定位和定姿误差,从而同时获得采煤机机身和摇臂的精确位置和姿态信息。In the coal mining operation of the shearer, the shearer body 5 moves along the track, while the rocker arm 6 moves to cut the coal wall; the rocker arm strapdown inertial navigation device 3 measures the angular velocity and acceleration information of the shearer rocker arm , using the strapdown navigation algorithm for positioning and attitude determination; the rocker arm wireless sensor mobile node 4 transmits wireless signals, and the rocker arm strapdown inertial navigation device 3 sends the inertial information at this moment to the fuselage strapdown inertial navigation device 1; The strapdown inertial navigation device 1 measures the angular velocity and acceleration information of the shearer fuselage in real time, and uses the strapdown navigation algorithm for positioning and attitude determination; the fuselage wireless sensor anchor node 2 receives the wireless signal of the rocker arm wireless sensor mobile node 4 , using the RSSI algorithm to calculate the distance; the fuselage SINS 1 calculates the inertial information of the fuselage at the time of arrival of the wireless signal, receives the inertial information of the rocker SINS 3, and communicates with the fuselage wireless sensor anchor The node 2 ranging and positioning results are subjected to collaborative navigation filtering to estimate the position and attitude error of the rocker arm; according to the error estimation results, the corrected rocker arm position and attitude information are sent to the rocker arm strapdown inertial navigation device, and the positioning and attitude of the rocker arm are corrected. The attitude error is determined, so as to obtain the precise position and attitude information of the shearer fuselage and rocker arm at the same time.

需要指出的是,本发明虽然以捷联惯导装置作为具体实施实例,但本发明只需对步骤一至步骤六的导航解算方式稍作更改,同样适用于平台惯导及其它形式的惯导装置。It should be pointed out that although the present invention takes the strapdown inertial navigation device as a specific implementation example, the present invention only needs to slightly change the navigation solution method in steps 1 to 6, and is also applicable to platform inertial navigation and other forms of inertial navigation. device.

显然,以上具体实施方式中实施例仅用于说明本发明的技术方案而非对其限制,尽管参照上述具体实施方式对本发明进行了详细说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本权利要求范围当中。Obviously, the embodiments in the above specific embodiments are only used to illustrate the technical solutions of the present invention and not to limit them. Although the present invention has been described in detail with reference to the above specific embodiments, those of ordinary skill in the art should understand that: the present invention can still be The specific embodiments of the present invention are modified or equivalently replaced, and any modification or equivalent replacement that does not depart from the spirit and scope of the present invention shall be included in the scope of the present claims.

Claims (3)

1.一种采煤机机身与摇臂协同定位装置,其特征在于:包括在采煤机机身(5)上安装的机身捷联惯导装置(1)和无线传感锚节点(2);以及在采煤机割煤一侧的摇臂(6)上安装的摇臂捷联惯导装置(3)和无线传感移动节点(4);其中无线传感移动节点(4)发射无线信号,无线传感锚节点(2)接收无线信号。1. A shearer fuselage and a rocker arm cooperative positioning device, characterized in that: comprising the fuselage strapdown inertial navigation device (1) and the wireless sensor anchor node (1) installed on the shearer fuselage (5). 2); and a rocker arm strapdown inertial navigation device (3) and a wireless sensor mobile node (4) installed on the rocker arm (6) on the coal cutting side of the shearer; wherein the wireless sensor mobile node (4) The wireless signal is transmitted, and the wireless sensor anchor node (2) receives the wireless signal. 2.根据权利要求1所述的一种采煤机机身与摇臂协同定位装置,其特征在于:所述机身捷联惯导装置(1)包括三只陀螺和三只加速度计。2. A shearer body and rocker arm cooperative positioning device according to claim 1, characterized in that: the body strapdown inertial navigation device (1) comprises three gyroscopes and three accelerometers. 3.基于权利要求1或2所述的采煤机机身与摇臂协同定位装置的定位方法,具体包括以下步骤:3. Based on the positioning method of the shearer fuselage and the rocker arm cooperative positioning device according to claim 1 or 2, the method specifically comprises the following steps: 步骤一、定义坐标系Step 1. Define the coordinate system 1)定义采煤机机身坐标系b1:Ob1Xb1Yb1Zb1:坐标系原点Ob1固连在机身捷联惯导装置1中心,Xb1轴正向由采煤机指向煤壁,Yb1轴垂直于Xb1轴向上,Zb1轴与Xb1轴、Yb1轴构成右手坐标系,构成前上右坐标系,下标1表示该坐标系为机身坐标系;1) Define the shearer fuselage coordinate system b1: O b1 X b1 Y b1 Z b1 : the origin of the coordinate system O b1 is fixedly connected to the center of the fuselage strapdown inertial navigation device 1, and the positive X b1 axis points from the shearer to the coal wall, the Y b1 axis is perpendicular to the X b1 axis upward, the Z b1 axis, the X b1 axis, and the Y b1 axis form a right-hand coordinate system, which constitutes the front upper right coordinate system, and the subscript 1 indicates that the coordinate system is the fuselage coordinate system; 2)定义采煤机摇臂坐标系b2:Ob2Xb2Yb2Zb2;坐标系原点Ob2固连在摇臂捷联惯导装置中心,Xb2轴正向由摇臂指向煤壁,Yb2轴垂直于Xb2轴向上,Zb2轴与Xb2轴、Yb2轴构成右手坐标系,构成前上右坐标系,下标2表示该坐标系为摇臂坐标系;2) Define the shearer rocker arm coordinate system b2: O b2 X b2 Y b2 Z b2 ; the origin of the coordinate system O b2 is fixedly connected to the center of the rocker arm strapdown inertial navigation device, and the positive X b2 axis is directed by the rocker arm to the coal wall, The Y b2 axis is perpendicular to the X b2 axis upward, and the Z b2 axis, the X b2 axis, and the Y b2 axis form a right-hand coordinate system, forming the front upper right coordinate system, and the subscript 2 indicates that the coordinate system is a rocker arm coordinate system; 3)定义导航坐标系n:OnXnYnZn:北天东地理坐标系,Xn轴指向地理北向,Yn轴指向天向,Zn轴指向地理东向;3) Define the navigation coordinate system n: O n X n Y n Z n : the north celestial east geographic coordinate system, the X n axis points to the geographic north direction, the Y n axis points to the sky direction, and the Z n axis points to the geographic east direction; 4)导航坐标系n经过三次旋转后与采煤机机身坐标系b1重合,三次旋转的角度即为机身航向角ψ1、俯仰角θ1和滚动角γ1;类似地,导航坐标系n经过三次旋转后与采煤机摇臂坐标系b2重合,三次旋转的角度即为摇臂航向角ψ2、俯仰角θ2和滚动角γ24) The navigation coordinate system n coincides with the shearer fuselage coordinate system b1 after three rotations, and the angles of the three rotations are the fuselage heading angle ψ 1 , pitch angle θ 1 and roll angle γ 1 ; similarly, the navigation coordinate system After three rotations, n coincides with the shearer rocker arm coordinate system b2, and the angles of the three rotations are rocker arm heading angle ψ 2 , pitch angle θ 2 and roll angle γ 2 ; 步骤二、计算采煤机机身信息Step 2. Calculate the information of the shearer fuselage 采煤机机身捷联惯导装置1利用陀螺和加速度计测量信息,解算采煤机机身姿态俯仰角θ1、航向角ψ1和滚动角γ1,计算机身的速度
Figure FDA0002463896760000021
位置
Figure FDA0002463896760000022
和经度λ1、纬度
Figure FDA0002463896760000023
高度h1
The shearer fuselage strapdown inertial navigation device 1 uses the gyro and accelerometer measurement information to calculate the shearer fuselage attitude pitch angle θ 1 , heading angle ψ 1 and roll angle γ 1 , and calculate the speed of the body
Figure FDA0002463896760000021
Location
Figure FDA0002463896760000022
and longitude λ 1 , latitude
Figure FDA0002463896760000023
height h 1 ;
1)机身捷联惯导装置1接收初始装订信息:初始位置
Figure FDA0002463896760000024
和初始姿态
Figure FDA0002463896760000025
其中上标1表示机身捷联惯导,N,U,E分别表示导航坐标系下北向、天向和东向位置;同时由于采煤机初始处于停止状态,初始速度V1=[0 0 0]T
1) Body SINS 1 receives initial binding information: initial position
Figure FDA0002463896760000024
and initial posture
Figure FDA0002463896760000025
The superscript 1 represents the fuselage strapdown inertial navigation, and N, U, and E represent the north, sky and east positions in the navigation coordinate system respectively. 0] T ;
2)机身捷联惯导装置1包括三只陀螺和三只加速度计,三只陀螺测量采煤机机身三个轴的角速度向量
Figure FDA0002463896760000026
三只加速度计测量采煤机机身三个轴向的加速度向量
Figure FDA0002463896760000027
2) The fuselage strapdown inertial navigation device 1 includes three gyroscopes and three accelerometers, and the three gyroscopes measure the angular velocity vectors of the three axes of the shearer fuselage
Figure FDA0002463896760000026
Three accelerometers measure the acceleration vectors in three axial directions of the shearer fuselage
Figure FDA0002463896760000027
3)机身捷联惯导装置1首先更新采煤机机身b1坐标系相对于导航坐标系n的角速度向量
Figure FDA0002463896760000028
3) The fuselage strapdown inertial navigation device 1 first updates the angular velocity vector of the shearer fuselage b1 coordinate system relative to the navigation coordinate system n
Figure FDA0002463896760000028
Figure FDA0002463896760000029
Figure FDA0002463896760000029
其中,
Figure FDA00024638967600000210
表示机身捷联惯导装置1纬度,
Figure FDA00024638967600000211
为机身捷联惯导装置1姿态矩阵,ωie是地球自转角速度,R为地球半径,
Figure FDA00024638967600000212
Figure FDA00024638967600000213
表示机身捷联惯导装置1东向速度和北向速度;
in,
Figure FDA00024638967600000210
Indicates the 1 latitude of the fuselage strapdown inertial navigation device,
Figure FDA00024638967600000211
is the attitude matrix of the fuselage SINS1, ω ie is the angular velocity of the earth's rotation, R is the earth's radius,
Figure FDA00024638967600000212
and
Figure FDA00024638967600000213
Indicates the east speed and north speed of the fuselage strapdown inertial navigation device 1;
4)更新四元数q1=[q0 q1 q2 q3]T 4) Update quaternion q 1 =[q 0 q 1 q 2 q 3 ] T
Figure FDA0002463896760000031
Figure FDA0002463896760000031
其中,
Figure FDA0002463896760000032
T表示姿态计算周期;
in,
Figure FDA0002463896760000032
T represents the attitude calculation period;
5)更新机身捷联惯导装置1相对于导航坐标系的姿态矩阵
Figure FDA0002463896760000033
5) Update the attitude matrix of the fuselage strapdown inertial navigation device 1 relative to the navigation coordinate system
Figure FDA0002463896760000033
Figure FDA0002463896760000034
Figure FDA0002463896760000034
其中q1=[q0 q1 q2 q3]T,令
Figure FDA0002463896760000035
其中i,j=1,2,3
where q 1 =[q 0 q 1 q 2 q 3 ] T , let
Figure FDA0002463896760000035
where i,j=1,2,3
6)根据
Figure FDA0002463896760000036
计算机身姿态角:
6) According to
Figure FDA0002463896760000036
Computer body attitude angle:
俯仰角θ1=sin-1(C12)Pitch angle θ 1 =sin −1 (C 12 ) 航向角
Figure FDA0002463896760000037
Heading
Figure FDA0002463896760000037
滚动角
Figure FDA0002463896760000038
roll angle
Figure FDA0002463896760000038
7)利用
Figure FDA0002463896760000039
将加速度计输出的加速度向量
Figure FDA00024638967600000310
投影到导航坐标系
Figure FDA00024638967600000311
7) Utilize
Figure FDA0002463896760000039
Convert the acceleration vector output by the accelerometer
Figure FDA00024638967600000310
Project to Navigation Coordinate System
Figure FDA00024638967600000311
Figure FDA00024638967600000312
Figure FDA00024638967600000312
8)更新机身捷联惯导装置1的速度
Figure FDA00024638967600000313
位置
Figure FDA00024638967600000314
和经度λ1、纬度
Figure FDA00024638967600000315
高度h1
8) Update the speed of the fuselage SINS 1
Figure FDA00024638967600000313
Location
Figure FDA00024638967600000314
and longitude λ 1 , latitude
Figure FDA00024638967600000315
height h 1 ,
Figure FDA0002463896760000041
Figure FDA0002463896760000041
Figure FDA0002463896760000042
Figure FDA0002463896760000042
Figure FDA0002463896760000043
Figure FDA0002463896760000043
Figure FDA0002463896760000044
Figure FDA0002463896760000044
Figure FDA0002463896760000045
Figure FDA0002463896760000045
步骤三、计算采煤机摇臂信息Step 3. Calculate the rocker arm information of the shearer 采煤机摇臂捷联惯导装置3利用陀螺和加速度计测量信息,解算采煤机摇臂姿态俯仰角θ2、航向角ψ2和滚动角γ2,计算摇臂的速度
Figure FDA0002463896760000046
位置
Figure FDA0002463896760000047
和经度λ2、纬度
Figure FDA0002463896760000048
高度h2
The shearer rocker arm strapdown inertial navigation device 3 uses the gyro and accelerometer measurement information to solve the shearer rocker arm attitude pitch angle θ 2 , heading angle ψ 2 and roll angle γ 2 , and calculate the rocker arm speed
Figure FDA0002463896760000046
Location
Figure FDA0002463896760000047
and longitude λ 2 , latitude
Figure FDA0002463896760000048
height h 2 ;
摇臂捷联惯导装置3接收初始装订信息:初始位置
Figure FDA0002463896760000049
和初始姿态
Figure FDA00024638967600000410
其中上标2表示摇臂捷联惯导,N,U,E分别表示导航坐标系下北向、天向和东向位置;同时由于采煤机初始处于停止状态,初始速度V2=[0 0 0]T
Rocker arm SINS 3 receives initial binding information: initial position
Figure FDA0002463896760000049
and initial posture
Figure FDA00024638967600000410
The superscript 2 represents the rocker arm strapdown inertial navigation, and N, U, and E represent the north, sky and east positions in the navigation coordinate system, respectively. 0] T ;
摇臂捷联惯导装置3与机身捷联惯导装置1采用相同的位置和姿态解算方法,将步骤二中1)至8)中变量上标1更换为2,即为摇臂捷联惯导装置3的位置和姿态解算方法;The rocker arm strapdown inertial navigation device 3 and the fuselage strapdown inertial navigation device 1 adopt the same position and attitude calculation method, and the variable superscript 1 in step 2 in 1) to 8) is replaced by 2, that is, the rocker arm The position and attitude calculation method of the coupled inertial navigation device 3; 步骤四、计算采煤机机身与摇臂的相对距离Step 4. Calculate the relative distance between the shearer body and the rocker arm 摇臂无线传感移动节点4相对于机身无线传感锚节点2测距:The rocker arm wireless sensor mobile node 4 measures the distance relative to the fuselage wireless sensor anchor node 2: 摇臂无线传感移动节点4发射无线信号,机身无线传感锚节点2接收移动节点的无线信号,利用RSSI算法进行测距,基于RSSI的测距方法将接收到的信号强度转化为无线传感移动节点4与无线传感锚节点2之间的距离,使用公知的对数-常数无线信号传播模型,The rocker arm wireless sensor mobile node 4 transmits wireless signals, the fuselage wireless sensor anchor node 2 receives the wireless signal of the mobile node, and uses the RSSI algorithm for ranging, and the RSSI-based ranging method converts the received signal strength into wireless transmission. sense the distance between the mobile node 4 and the wireless sensing anchor node 2, using the well-known log-constant wireless signal propagation model, S=A-10mlg(d) (8)S=A-10mlg(d) (8) 其中,S为信号强度,d为距离,A为接收端距离发射端为1m时接收到的无线信号强度,m为路径损耗,A和m为已知参数,根据S计算距离d;Among them, S is the signal strength, d is the distance, A is the received wireless signal strength when the receiving end is 1m away from the transmitting end, m is the path loss, A and m are known parameters, and the distance d is calculated according to S; 根据公式(8),根据无线传感锚节点2接收到的信号强度S,可得摇臂相对于机身的距离d;According to formula (8), according to the signal strength S received by the wireless sensor anchor node 2, the distance d of the rocker arm relative to the fuselage can be obtained; 在无线传感移动节点4发射无线信号时,摇臂捷联惯导装置3记录该时刻惯导解算的位置信息,并发送给机身捷联惯导装置1;机身无线传感锚节点2接收无线信号时,接收摇臂捷联惯导装置3发送的测距时刻位置信息,同时记录机身捷联惯导装置1解算的机身位置姿态信息及摇臂位置姿态信息;When the wireless sensor mobile node 4 transmits a wireless signal, the rocker arm SINS device 3 records the position information calculated by the inertial navigation at that moment, and sends it to the fuselage SINS device 1; the fuselage wireless sensor anchor node 2. When receiving the wireless signal, receive the position information at the time of ranging sent by the rocker arm strapdown inertial navigation device 3, and simultaneously record the body position and attitude information and the rocker arm position and attitude information calculated by the fuselage strapdown inertial navigation device 1; 步骤五、利用步骤四中在无线测距信号发射时刻记录的机身位置姿态信息、摇臂位置姿态信息、机身与摇臂的相对距离,进行基于卡尔曼滤波器的协同导航滤波解算,估计摇臂捷联惯导装置3的位置、速度和姿态误差;Step 5: Use the fuselage position and attitude information, the rocker arm position and attitude information, and the relative distance between the fuselage and the rocker arm recorded at the time of transmitting the wireless ranging signal in step 4, and perform collaborative navigation filtering based on Kalman filter. Estimate the position, velocity and attitude errors of the rocker arm SINS 3; 1)协同导航滤波器的状态变量为
Figure FDA0002463896760000051
其中ΔPI=[ΔPN ΔPU ΔPE]T表示摇臂捷联惯导装置3位置误差;ΔV=[ΔVN ΔVU ΔVE]T表示摇臂捷联惯导装置3速度误差;φ=[φN φU φE]T表示摇臂捷联惯导装置3姿态误差;
1) The state variable of the collaborative navigation filter is
Figure FDA0002463896760000051
Wherein ΔP I =[ΔP N ΔP U ΔP E ] T represents the position error of the rocker arm SINS 3; ΔV=[ΔV N ΔV U ΔV E ] T represents the speed error of the rocker arm SINS 3; φ= [φ N φ U φ E ] T represents the attitude error of the rocker arm strapdown inertial navigation device 3;
2)协同导航滤波状态方程由摇臂捷联惯导装置3的速度误差方程、位置误差方程和姿态误差方程组成:2) The cooperative navigation filtering state equation is composed of the velocity error equation, position error equation and attitude error equation of the rocker arm SINS device 3: 速度误差方程:Velocity Error Equation:
Figure FDA0002463896760000061
Figure FDA0002463896760000061
其中,in,
Figure FDA0002463896760000062
Figure FDA0002463896760000062
Figure FDA0002463896760000063
Figure FDA0002463896760000063
位置误差方程Position Error Equation
Figure FDA0002463896760000064
Figure FDA0002463896760000064
姿态误差方程Attitude Error Equation
Figure FDA0002463896760000065
Figure FDA0002463896760000065
将方程(9)、(10)、(11)进行离散化并改写为协同导航滤波状态方程形式,其中tk表示滤波时间点;The equations (9), (10) and (11) are discretized and rewritten into the form of cooperative navigation filtering state equation, where t k represents the filtering time point; X(tk)=F(tk-1)X(tk-1) (12)X(t k )=F(t k-1 )X(t k-1 ) (12) 3)协同导航滤波滤波量测方程如下:3) The collaborative navigation filtering filtering measurement equation is as follows:
Figure FDA0002463896760000066
Figure FDA0002463896760000066
4)修正机身捷联惯装置1导同步误差4) Correct the 1-lead synchronization error of the fuselage strapdown inertial device 机身捷联惯导装置1定位数据更新率高于无线测距的定位数据更新率,可采用线性插值方法,在摇臂无线传感移动节点4发射无线信号时,机身捷联惯导装置1记录时间间隔Δt,Δt<ΔT,采用插值方法,计算用于协同导航的机身捷联惯导装置1的位置P1The update rate of the positioning data of the fuselage strapdown inertial navigation device 1 is higher than that of the wireless ranging, and the linear interpolation method can be used. When the rocker arm wireless sensing mobile node 4 transmits wireless signals, the fuselage strapdown inertial navigation 1. Record the time interval Δt, Δt<ΔT, and use the interpolation method to calculate the position P 1 of the fuselage strapdown inertial navigation device 1 used for cooperative navigation:
Figure FDA0002463896760000071
Figure FDA0002463896760000071
5)计算量测量Z5) Calculate the amount of measurement Z 在无线测距信号发射时刻,机身无线传感锚节点2位置(即机身捷联惯导装置位置)为
Figure FDA0002463896760000072
在摇臂无线传感移动4节点位置即摇臂捷联惯导装置3位置为
Figure FDA0002463896760000073
根据几何关系有:
At the moment when the wireless ranging signal is transmitted, the position of the fuselage wireless sensor anchor node 2 (that is, the position of the fuselage strapdown inertial navigation device) is
Figure FDA0002463896760000072
When the rocker arm wireless sensor moves 4 nodes, that is, the rocker arm strapdown inertial navigation device 3 position is:
Figure FDA0002463896760000073
According to the geometric relationship there are:
Figure FDA0002463896760000074
Figure FDA0002463896760000074
其中,Y表示利用摇臂捷联惯导装置3位置与机身捷联惯导装置1位置计算的相对距离。那么协同导航滤波的量测量Z计算如下:Among them, Y represents the relative distance calculated by using the position of the rocker arm SINS device 3 and the position of the fuselage SINS device 1 . Then the quantity measurement Z of the collaborative navigation filtering is calculated as follows: Z=Y-dZ=Y-d
Figure FDA0002463896760000075
Figure FDA0002463896760000075
6)机身捷联惯导装置1利用公知的卡尔曼滤波算法,可得协同导航滤波状态量估计值
Figure FDA0002463896760000076
6) The fuselage strapdown inertial navigation device 1 uses the well-known Kalman filter algorithm to obtain the estimated value of the cooperative navigation filter state quantity
Figure FDA0002463896760000076
步骤六、利用协同导航滤波状态估计值
Figure FDA0002463896760000077
修正摇臂捷联惯导装置3的位置、速度和姿态误差,得到误差修正后的摇臂的位置、速度和姿态:
Step 6. Use collaborative navigation to filter state estimates
Figure FDA0002463896760000077
Correct the position, speed and attitude errors of the rocker arm strapdown inertial navigation device 3, and obtain the position, speed and attitude of the rocker arm after error correction:
机身捷联惯导装置1将
Figure FDA0002463896760000078
发送给摇臂捷联惯导装置3,对其误差进行闭环点修正:
The fuselage strapdown inertial navigation device 1 will be
Figure FDA0002463896760000078
Send it to the rocker arm strapdown inertial navigation device 3, and perform closed-loop point correction for its error:
P2(tk)=P2(tk)-ΔP(tk) (17)P 2 (t k )=P 2 (t k )-ΔP(t k ) (17) V2(tk)=V2(tk)-ΔV(tk) (18)V 2 (t k )=V 2 (t k )-ΔV(t k ) (18)
Figure FDA0002463896760000081
Figure FDA0002463896760000081
得到误差修正后的装置3位置、速度和姿态矩阵,对修正误差后的
Figure FDA0002463896760000082
按照步骤二的6),将变量上标1更换为2,计算误差修正后的摇臂姿态角:俯仰角θ2、航向角ψ2和滚动角γ2
The position, velocity and attitude matrix of device 3 after error correction is obtained, and the error correction
Figure FDA0002463896760000082
According to 6) of step 2, replace the variable superscript 1 with 2, and calculate the attitude angle of the rocker arm after error correction: pitch angle θ 2 , heading angle ψ 2 and roll angle γ 2 .
CN202010327977.XA 2020-04-23 2020-04-23 Cooperative positioning device and positioning method of shearer body and rocker arm Pending CN111323011A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010327977.XA CN111323011A (en) 2020-04-23 2020-04-23 Cooperative positioning device and positioning method of shearer body and rocker arm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010327977.XA CN111323011A (en) 2020-04-23 2020-04-23 Cooperative positioning device and positioning method of shearer body and rocker arm

Publications (1)

Publication Number Publication Date
CN111323011A true CN111323011A (en) 2020-06-23

Family

ID=71172000

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010327977.XA Pending CN111323011A (en) 2020-04-23 2020-04-23 Cooperative positioning device and positioning method of shearer body and rocker arm

Country Status (1)

Country Link
CN (1) CN111323011A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111894582A (en) * 2020-08-04 2020-11-06 中国矿业大学 A kind of control method of coal shearer
CN112902952A (en) * 2021-01-21 2021-06-04 中国煤炭科工集团太原研究院有限公司 Development machine pose detection method based on millimeter wave radar and inertial navigation
CN114662259A (en) * 2020-12-07 2022-06-24 北斗天地股份有限公司 Coal mining machine three-dimensional positioning method and system based on high-precision inertial navigation
CN115326064A (en) * 2022-07-19 2022-11-11 西京学院 Coal mining machine body and rocker arm cooperative high-precision positioning device and positioning method
CN116698077A (en) * 2023-04-25 2023-09-05 天地上海采掘装备科技有限公司 Coal mining machine posture detection method and device and computer readable storage medium
CN118746297A (en) * 2024-07-26 2024-10-08 北京中航天佑科技有限公司 A positioning method and system for anchor digging machine
CN120063289A (en) * 2025-04-25 2025-05-30 北京星际导控科技有限责任公司 Dynamic navigation positioning method and system for intelligent fully-mechanized mining of coal mine

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010151459A (en) * 2008-12-24 2010-07-08 Seiko Epson Corp Method and device for calculating position
CN103413185A (en) * 2013-08-07 2013-11-27 中国矿业大学 Optimal arrangement method of coal mining machine rocker arm vibration sensors
CN108387227A (en) * 2018-02-22 2018-08-10 北京航空航天大学 The multinode information fusion method and system of airborne distribution POS
CN109974715A (en) * 2019-04-30 2019-07-05 中国矿业大学(北京) Autonomous navigation system and method for roadheader combining strapdown inertial navigation and light spot recognition
CN110118550A (en) * 2019-05-14 2019-08-13 中国矿业大学 A kind of generation system and method for fully-mechanized mining working high-precision three-dimensional navigation map
CN110162036A (en) * 2019-04-09 2019-08-23 中国矿业大学 A kind of development machine Camera calibration system and method
CN110702109A (en) * 2019-06-05 2020-01-17 西京学院 Coal mining machine inertial navigation/wireless sensor network combined positioning method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010151459A (en) * 2008-12-24 2010-07-08 Seiko Epson Corp Method and device for calculating position
CN103413185A (en) * 2013-08-07 2013-11-27 中国矿业大学 Optimal arrangement method of coal mining machine rocker arm vibration sensors
CN108387227A (en) * 2018-02-22 2018-08-10 北京航空航天大学 The multinode information fusion method and system of airborne distribution POS
CN110162036A (en) * 2019-04-09 2019-08-23 中国矿业大学 A kind of development machine Camera calibration system and method
CN109974715A (en) * 2019-04-30 2019-07-05 中国矿业大学(北京) Autonomous navigation system and method for roadheader combining strapdown inertial navigation and light spot recognition
CN110118550A (en) * 2019-05-14 2019-08-13 中国矿业大学 A kind of generation system and method for fully-mechanized mining working high-precision three-dimensional navigation map
CN110702109A (en) * 2019-06-05 2020-01-17 西京学院 Coal mining machine inertial navigation/wireless sensor network combined positioning method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SCHIFFBAUER W H 等: ""Navigation and control of continuous mining systems for coal mining"", 《CONFERENCE RECORD OF THE 1996 IEEE INDUSTRY APPLICATIONS CONFERENCE THIRTY-FIRST IASANNUAL MEETING》 *
巩剑波等: "采煤机姿态矫正系统研究综述", 《煤炭技术》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111894582A (en) * 2020-08-04 2020-11-06 中国矿业大学 A kind of control method of coal shearer
CN111894582B (en) * 2020-08-04 2021-09-24 中国矿业大学 A kind of control method of coal shearer
CN114662259A (en) * 2020-12-07 2022-06-24 北斗天地股份有限公司 Coal mining machine three-dimensional positioning method and system based on high-precision inertial navigation
CN112902952A (en) * 2021-01-21 2021-06-04 中国煤炭科工集团太原研究院有限公司 Development machine pose detection method based on millimeter wave radar and inertial navigation
CN112902952B (en) * 2021-01-21 2023-11-10 中国煤炭科工集团太原研究院有限公司 Heading machine pose detection method based on millimeter wave radar and inertial navigation
CN115326064A (en) * 2022-07-19 2022-11-11 西京学院 Coal mining machine body and rocker arm cooperative high-precision positioning device and positioning method
CN116698077A (en) * 2023-04-25 2023-09-05 天地上海采掘装备科技有限公司 Coal mining machine posture detection method and device and computer readable storage medium
CN118746297A (en) * 2024-07-26 2024-10-08 北京中航天佑科技有限公司 A positioning method and system for anchor digging machine
CN120063289A (en) * 2025-04-25 2025-05-30 北京星际导控科技有限责任公司 Dynamic navigation positioning method and system for intelligent fully-mechanized mining of coal mine

Similar Documents

Publication Publication Date Title
CN110702109B (en) Coal mining machine inertial navigation/wireless sensor network combined positioning method
CN111323011A (en) Cooperative positioning device and positioning method of shearer body and rocker arm
CN109813311B (en) A UAV formation collaborative navigation method
CN110792430B (en) A method and device for measuring inclination while drilling based on multi-sensor data fusion
CN103410512B (en) Merge coal-winning machine absolute positioning apparatus and the method for geological environment information
EP1478903B1 (en) Device for use with a portable inertial navigation system (pins) and method for processing pins signals
US11274788B2 (en) Gimbal pose correction method and device
CN103335647B (en) A kind of attitude of shield machine measuring system and measuring method thereof
CN109540130A (en) A kind of continuous milling machine inertial navigation positioning and orienting method
CN114753841B (en) Coal mining machine posture correction method and device
CN112378399A (en) Coal mine tunnel tunneling robot precise positioning and orientation method based on strapdown inertial navigation and digital total station
CN114739425A (en) Coal mining machine positioning calibration system based on RTK-GNSS and total station and application method
WO2017092181A1 (en) Combined initial alignment system and alignment method for strapdown inertial navigation system of underground coal mining machine
CN110672099A (en) Course correction method and system for indoor robot navigation
CN113236363B (en) Mining equipment navigation positioning method, system, equipment and readable storage medium
CN104515527A (en) Anti-rough error integrated navigation method under non-GPS signal environment
CN113075650A (en) Underground roadway tunneling equipment real-time positioning method based on UWB and inertial unit
CN110847905A (en) Autonomous navigation system and method for coal mining machine
CN108868772B (en) A fast collimation control method for continuous miner
CN110095135A (en) A kind of method and device for development machine positioning and directing
CN111207743A (en) Method for realizing centimeter-level accurate positioning based on close coupling of encoder and inertial equipment
CN110274591A (en) Deep diving manned underwater vehicle ADCP assists SINS air navigation aid
CN115655267A (en) Underground positioning and orientation system based on millimeter wave radar ranging of inertial navigation odometer
CN115655268A (en) Underground mine heading machine navigation positioning method based on graph optimization algorithm
CN115326064A (en) Coal mining machine body and rocker arm cooperative high-precision positioning device and positioning method

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