WO2022209303A1 - 作業機械 - Google Patents

作業機械 Download PDF

Info

Publication number
WO2022209303A1
WO2022209303A1 PCT/JP2022/004694 JP2022004694W WO2022209303A1 WO 2022209303 A1 WO2022209303 A1 WO 2022209303A1 JP 2022004694 W JP2022004694 W JP 2022004694W WO 2022209303 A1 WO2022209303 A1 WO 2022209303A1
Authority
WO
WIPO (PCT)
Prior art keywords
positioning
satellite
antenna
positioning accuracy
mask
Prior art date
Application number
PCT/JP2022/004694
Other languages
English (en)
French (fr)
Inventor
惠貞 李
伸一 小竹
枝穂 泉
新士 石原
Original Assignee
日立建機株式会社
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 日立建機株式会社 filed Critical 日立建機株式会社
Priority to US18/024,991 priority Critical patent/US20230374757A1/en
Priority to EP22779532.5A priority patent/EP4317617A1/en
Priority to KR1020237007658A priority patent/KR20230045068A/ko
Priority to CN202280005880.4A priority patent/CN116057417A/zh
Publication of WO2022209303A1 publication Critical patent/WO2022209303A1/ja

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02FDREDGING; SOIL-SHIFTING
    • E02F9/00Component parts of dredgers or soil-shifting machines, not restricted to one of the kinds covered by groups E02F3/00 - E02F7/00
    • E02F9/26Indicating devices
    • E02F9/264Sensors and their calibration for indicating the position of the work tool
    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02FDREDGING; SOIL-SHIFTING
    • E02F3/00Dredgers; Soil-shifting machines
    • E02F3/04Dredgers; Soil-shifting machines mechanically-driven
    • E02F3/28Dredgers; Soil-shifting machines mechanically-driven with digging tools mounted on a dipper- or bucket-arm, i.e. there is either one arm or a pair of arms, e.g. dippers, buckets
    • E02F3/30Dredgers; Soil-shifting machines mechanically-driven with digging tools mounted on a dipper- or bucket-arm, i.e. there is either one arm or a pair of arms, e.g. dippers, buckets with a dipper-arm pivoted on a cantilever beam, i.e. boom
    • E02F3/32Dredgers; Soil-shifting machines mechanically-driven with digging tools mounted on a dipper- or bucket-arm, i.e. there is either one arm or a pair of arms, e.g. dippers, buckets with a dipper-arm pivoted on a cantilever beam, i.e. boom working downwardly and towards the machine, e.g. with backhoes
    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02FDREDGING; SOIL-SHIFTING
    • E02F9/00Component parts of dredgers or soil-shifting machines, not restricted to one of the kinds covered by groups E02F3/00 - E02F7/00
    • E02F9/26Indicating devices
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C15/00Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/14Receivers specially adapted for specific applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/28Satellite selection
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60YINDEXING SCHEME RELATING TO ASPECTS CROSS-CUTTING VEHICLE TECHNOLOGY
    • B60Y2200/00Type of vehicle
    • B60Y2200/40Special vehicles
    • B60Y2200/41Construction vehicles, e.g. graders, excavators
    • B60Y2200/412Excavators
    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02FDREDGING; SOIL-SHIFTING
    • E02F3/00Dredgers; Soil-shifting machines
    • E02F3/04Dredgers; Soil-shifting machines mechanically-driven
    • E02F3/28Dredgers; Soil-shifting machines mechanically-driven with digging tools mounted on a dipper- or bucket-arm, i.e. there is either one arm or a pair of arms, e.g. dippers, buckets
    • E02F3/36Component parts
    • E02F3/42Drives for dippers, buckets, dipper-arms or bucket-arms
    • E02F3/43Control of dipper or bucket position; Control of sequence of drive operations
    • E02F3/435Control of dipper or bucket position; Control of sequence of drive operations for dipper-arms, backhoes or the like
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/53Determining attitude

Definitions

  • the present invention relates to a work machine whose position is detected using a positioning satellite system.
  • Information-aided construction refers to a series of construction processes consisting of research, design, construction, inspection, management, etc., which focuses on construction and utilizes electronic data and ICT (Information and Communication Technology) to facilitate construction. It is a system that achieves high efficiency.
  • Machines that support information-aided construction include a guidance function that displays the vehicle body position and the position and posture of the front working device (sometimes referred to as a working device) along with the position data of the work target plane on the monitor, and a bucket that displays the work target plane. 2. Description of the Related Art
  • Working machines typified by hydraulic excavators equipped with a machine control function for controlling a front working device to prevent over-digging are known.
  • Work machines that are compatible with such information-aided construction present information to the operator based on information-aided construction data that has three-dimensional coordinate data, and provide work support and operation support functions. For example, in the machine guidance of a hydraulic excavator, the position of the tip of the bucket is calculated from data on the position and attitude of the vehicle body and the attitude data of the front working device, and the position of the bucket relative to the target work surface is presented to the operator via a monitor.
  • This type of hydraulic excavator receives positioning signals from positioning satellites via a positioning antenna attached to the upper revolving structure in order to calculate the position of the upper revolving structure (body) in the global coordinate system (geographical coordinate system).
  • a satellite positioning system for example, GNSS (Global Navigation Satellite System)
  • GNSS Global Navigation Satellite System
  • the front parts such as the boom, arm, and bucket may exist above the positioning antenna of the satellite positioning system, and these may interfere with the reception of positioning signals on straight routes. be.
  • the positioning antenna receives the positioning signal as a diffracted wave or a reflected wave called multipath. If a diffracted wave or reflected wave is received and used for positioning calculation, there is a high possibility that an error will be included in the positioning result.
  • Patent Document 1 discloses a technique that attempts to reduce the influence of multipath.
  • a mask range (mask information) representing a selection rule of GPS satellites determined based on the arrangement of radio wave obstacles around a GPS receiver is stored in a database for each predetermined area, and GPS data is transmitted from the database.
  • GPS receiver control that acquires the mask range corresponding to the signal reception area and selects the GPS satellite to be used for positioning from among the GPS satellites that are located outside the mask range among the multiple GPS satellites that fly in the sky.
  • An apparatus is disclosed.
  • Patent Document 1 assumes only non-moving radio wave obstacles, and the mask range when there are movable radio wave obstacles around the positioning antenna, such as the front work equipment of a hydraulic excavator. not touched. In other words, in order to improve satellite positioning accuracy in the field of work machinery, it is important to set a mask range that takes into consideration the attitude of the front work equipment, which is a movable radio wave obstruction.
  • the mask range used for satellite selection is changed each time (or changed) according to changes in the attitude of the front working device, it would be possible to suppress the decrease in positioning accuracy, but the story is not so simple.
  • the posture of the front working device is frequently changed (eg, raising and lowering the boom). If the mask range is changed in accordance with the change in attitude, the satellites used for positioning calculations will switch frequently, and changing the mask range may result in a decrease in positioning accuracy.
  • the number of satellites that can be used for positioning calculation is extremely small, if the mask range is used according to the attitude of the front work device, there will be times when even the few satellites cannot be used for positioning. In some cases, not using it can suppress the decrease in positioning accuracy.
  • positioning accuracy degradation can also occur due to, for example, the number of satellites and the SN ratio (signal-to-noise ratio) of satellite signals.
  • the present invention has been made in view of the above circumstances, and its object is to provide a work machine that can suppress the deterioration of satellite positioning accuracy due to the operation of the work equipment.
  • the present application includes a plurality of means for solving the above-mentioned problems, and to give an example, an undercarriage, an upper revolving body rotatably mounted on the undercarriage, and the upper revolving body An attached articulated working device, a plurality of attitude sensors for detecting the attitudes of the working device and the upper slewing structure, and attached to the upper slewing structure for receiving satellite signals from a plurality of positioning satellites.
  • a receiver for calculating the position of the positioning antenna based on satellite signals received by the positioning antenna; and a receiver for calculating the position of the positioning antenna based on the satellite signals received by the positioning antenna; and a controller for calculating the attitude of the upper slewing structure, wherein the controller controls the working device and the positioning satellite for selecting the positioning satellite used for calculating the position of the positioning antenna by the receiver.
  • a plurality of mask ranges set with reference to the positioning antenna based on the attitude of the upper slewing body are stored, and the controller uses each of the plurality of mask ranges to select a positioning satellite.
  • the receiver is characterized by calculating the position of the positioning antenna based on the satellite signal of the positioning satellite selected when the positioning accuracy is the best among the plurality of positioning accuracies.
  • the mask range reflecting the posture of the front work device changes with the operation of the front work device
  • the mask range with the highest positioning accuracy is selected.
  • FIG. 1 is a side view of a hydraulic excavator 1 and a GNSS reference station 8 according to an embodiment of the present invention
  • FIG. FIG. 2 is a functional block diagram of an in-vehicle controller 40 mounted on the hydraulic excavator 1 of FIG. 1
  • FIG. 3 is a perspective view of a first antenna coordinate system based on the GNSS antenna 50A
  • FIG. 2 is a view (plan view) looking down on the hydraulic excavator 1 from the Z-axis in the first antenna coordinate system.
  • the figure (side view) which looked at the hydraulic excavator 1 from +Y-axis direction in a 1st antenna coordinate system.
  • FIG. 4 is an explanatory diagram of angles that define the posture of the front working device 6;
  • 5 is a flowchart of positioning processing of the GNSS antenna 50A by the controller 40 and the GNSS receiver 51 according to the present embodiment;
  • the present invention is applied to a crawler hydraulic excavator as a work machine. It also has a machine control function that limits the operation of the work equipment (that is, the operation of the actuator that drives the front member) so as not to exceed the construction target plane.
  • the same reference numerals are given to the same parts in each figure, and overlapping explanations are omitted as appropriate.
  • FIG. 1 is a side view of a hydraulic excavator 1 and a GNSS reference station 8 according to an embodiment of the invention.
  • the hydraulic excavator 1 shown in this figure includes a crawler-type running body (lower running body) 2, a revolving body (upper revolving body) 3 attached to the upper part of the running body 2 so as to be able to turn, and one end (base end) of which is A front working device (sometimes simply referred to as “working device”) 6 consisting of a multi-joint type link mechanism is attached to the front of the revolving body 3 .
  • Reference numeral 30 in the figure represents the ground.
  • the front working device 6 has a boom 6A whose one end is connected to the revolving body 3, an arm 6B whose one end is connected to the other end of the boom 6A, and a bucket 6C whose one end is connected to the other end of the arm 6B.
  • Each of these front members 6A, 6B, 6C is configured to rotate vertically.
  • a boom cylinder 11A, an arm cylinder 11B, and a bucket cylinder 11C are provided as actuators (hydraulic cylinders) for rotating the front members 6A, 6B, and 6C.
  • the revolving body 3 can be driven to turn left and right around a turning center axis O by a turning motor (not shown).
  • the boom 6A, arm 6B and bucket 6C operate on a common plane including the front working device 6, and hereinafter this plane may be referred to as the action plane.
  • the action plane is a plane perpendicular to the rotation axes of the boom 6A, arm 6B and bucket 6C.
  • the center of the boom 6A, arm 6B and bucket 6C in the width direction center of the rotation axis may be referred to as the action plane.
  • the hydraulic excavator 1 is provided with a plurality of attitude sensors 75A, 75B, 75C, and 23 for detecting the attitudes of the front working device 6 and the revolving body 3.
  • attitude sensors 75A, 75B, 75C, and 23 for detecting the attitudes of the front working device 6 and the revolving body 3.
  • IMU inertial measurement unit
  • a boom attitude sensor 75A is attached to the boom 6A
  • an arm attitude sensor 75B is attached to the arm 6B
  • a bucket attitude sensor 75C is attached to the bucket 6C (see FIG. 1).
  • a revolving structure attitude sensor 23 is attached to the revolving structure 3 (see FIG.
  • attitude sensors 75A, 75B, 75C, and 23 are input to the controller 40 via connection lines.
  • An angle sensor for example, a potentiometer or a rotary encoder
  • the bucket attitude sensor 75C may be attached to the bucket link instead of the bucket.
  • the revolving structure 3 has a plurality of control levers (not shown) operated by an operator, an operator's seat 4 having a monitor 60 displaying the positional relationship between the bucket 6C and the work target surface, and a plurality of positioning satellites.
  • a controller 40 which is a computer that calculates desired position coordinates on the front work device 6 based on the position and orientation calculated by the GNSS receiver 51 and detection signals from the plurality of attitude sensors 75A, 75B, 75C, and 23.
  • the positions of the two GNSS antennas 50A and 50B and the azimuth angle of the revolving body 3 are calculated by one GNSS receiver.
  • a configuration in which two GNSS receivers 51A and 51B are mounted may be adopted.
  • the GNSS reference station 8 that wirelessly transmits the GNSS correction data to the wireless device 7 of the hydraulic excavator 1 will be described.
  • the GNSS reference station 8 whose coordinate positions in the geographic coordinate system are known, has a GNSS antenna 80 for receiving satellite signals from a plurality of positioning satellites (GNSS satellites) and satellite signals received by the GNSS antenna 80 (satellite signals (includes satellite code, carrier, satellite orbit, satellite signal reception level, etc.);
  • a reference station controller 82 for generating GNSS correction data for wireless transmission to the radio 7 based on satellite signals, and a radio 87 for transmitting the GNSS correction data generated by the reference station controller 82 to the radio 7 are provided. .
  • a GNSS receiver 81 connected to a GNSS reference station antenna 80 wirelessly transmits GNSS correction data from a wireless device 87 via a reference station controller 82 . If the GNSS correction data received by the wireless device 7 is used for positioning by the GNSS receiver 51, highly accurate positioning on the centimeter level becomes possible.
  • ⁇ GNSS antenna 50> The two GNSS antennas 50A and 50B are fixed to the upper revolving body 3 via masts (antenna support members) 52a and 52b, respectively. are arranged at predetermined intervals.
  • the two masts 52a and 52b are pole-shaped support members for supporting the GNSS antennas 50A and 50B above the upper swing body 3, respectively.
  • the two masts 52a and 52b of this embodiment are arranged on the upper surface (first area) of the upper swing body 3, like the GNSS antennas 50A and 50B.
  • Base ends of the masts 52a and 52b are fixed to the upper surface of the upper rotating body 3, and the masts 52a and 52b extend substantially vertically from the base ends.
  • GNSS antennas 50A, 50B having a substantially disc-shaped external shape with the central part bulging in the axial direction are attached.
  • Each antenna 50A, 50B is supported so as to pass through the center axis of the GNSS antenna 50A, 50B.
  • the support members for the GNSS antennas 50A and 50B are not limited to the pole-shaped masts 52a and 52b, and can be supported by support members of various shapes.
  • the GNSS receiver 51 receives the two GNSS antennas 50A, 50B based on a plurality of satellite signals (including satellite codes, carriers, satellite orbits, satellite signal reception levels, etc.), the two GNSS antennas 50A, 50B. 50B of at least one GNSS antenna (for example, the GNSS antenna 50B) in the geographic coordinate system (global coordinate system) and the orientation between the two GNSS antennas 50A and 50B Azimuth (also called heading) is calculated.
  • a plurality of satellite signals including satellite codes, carriers, satellite orbits, satellite signal reception levels, etc.
  • the two GNSS antennas 50A, 50B. 50B of at least one GNSS antenna for example, the GNSS antenna 50B
  • Azimuth also called heading
  • Electromagnetic waves (satellite signals) containing transmission time information are transmitted from multiple positioning satellites.
  • the GNSS receiver 51 calculates the arrival time difference from the reception time of the electromagnetic wave from each GNSS satellite and the transmission time included in the electromagnetic wave, and based on the arrival time difference, the distance between each GNSS satellite and the GNSS antennas 50A and 50B. are estimated to calculate the positions of the GNSS antennas 50A and 50B.
  • the GNSS satellites are equipped with sophisticated clocks, and the distance between each GNSS satellite and the GNSS antenna is calculated by multiplying the arrival time difference obtained by demodulating the electromagnetic waves from each satellite by the speed of the electromagnetic waves.
  • An error may be included in the calculated distance between each GNSS satellite and each GNSS antenna. This error is caused by the fact that the speed change of the electromagnetic wave generated by the ionosphere and water vapor existing between the GNSS satellite and the GNSS antenna is different for each GNSS satellite with different azimuth and elevation angles, and the orbital information sent by electromagnetic waves from each GNSS satellite. is slightly different from the actual position, and there are some errors in the clock information between the GNSS satellites.
  • RTK-GNSS real-time kinematic GNSS
  • receives GNSS correction data transmitted from the GNSS reference station 8 and performs positioning For example, the positioning of the reference station GNSS antenna 80 with a known absolute position installed near the hydraulic excavator 1 (within several kilometers) and the calculation of the GNSS correction data are performed by the reference station GNSS receiver 81, and the correction data is transmitted by the wireless device 87. It is transmitted to the receiver 51 of the excavator 1 .
  • the relative position vector instead of the absolute position between the two GNSS antennas 50A (50B) and 80, errors can be reduced.
  • the correction data transmitted from the radio 87 of the GNSS reference station 8 is received by the radio 7 mounted on the hydraulic excavator 1 and transmitted to the GNSS receiver 51 .
  • the relative position between the reference station GNSS antenna 80 and the GNSS antenna 50A is calculated by comparing the satellite signal received by the GNSS antenna 50A (mobile station) and the signal of the reference station GNSS antenna 80 obtained from the correction data. Calculate (direction and distance).
  • the carrier wave phase information of the satellite signal from the satellite received by the base station antenna 80 is transmitted as correction information, and this is compared with the carrier wave phase information of the satellite signal received by the mobile station antenna 50A by the GNSS receiver 51.
  • the GNSS receiver 51 can output the positioning results of the GNSS antennas 50A and 50B in the NMEA format including the latitude, longitude and geoid height of each of the GNSS antennas 50A and 50B.
  • one GNSS antenna 50A can be regarded as a reference station and the other GNSS antenna 50B can be regarded as a mobile station.
  • Such a method is the moving base method. It is possible to measure the relative position (vector) between the two GNSS antennas 50A and 50B by using the correction data generated by the received signal of the GNSS antenna 50A to measure the relative position (vector) with the GNSS antenna 50B. becomes.
  • relative positions (vectors) can be calculated without using correction data transmitted from the wireless device 87 .
  • a direction calculation method there is also a method of calculating the positions of the GNSS antenna 50A and the GNSS antenna 50B from the reference station GNSS antenna 80 and obtaining the direction from the difference between the positions. Then, in the direction between the two GNSS antennas 50A and 50B calculated in this way, by considering the constant due to the mounting positions of the two GNSS antennas 50A and 50B on the excavator 1, and the azimuth (direction) of the front work device 6 can be calculated.
  • a system for calculating the directions of the upper rotating body 3 and the front work device 6 by wirelessly transmitting correction data from the reference station GNSS antenna 80 has been described, but VRS (virtual reference point system), quasi-zenith satellite, etc. have been described. You may use the service which distributes the correction
  • FIG. 2 is a functional block diagram of the controller 40 mounted on the hydraulic excavator 1 of FIG.
  • the controller 40 Based on the positions of the two GNSS antennas 50A and 50B and the azimuth (heading) of the revolving structure 3 calculated by the GNSS receiver 51, and the detection signals of the plurality of attitude sensors 75A, 75B, 75C, and 23, the controller 40: It is a computer that calculates the position coordinates of each of the front members 6A, 6B, and 6C that constitute the front working device 6.
  • FIG. 1 It is a computer that calculates the position coordinates of each of the front members 6A, 6B, and 6C that constitute the front working device 6.
  • the controller 40 includes an arithmetic processing unit (for example, CPU (not shown)), a storage device (for example, semiconductor memory such as ROM and RAM) 56, and an interface (input/output device (not shown)).
  • the program (software) stored in advance in 56 is executed by the arithmetic processing unit, and the arithmetic processing unit performs arithmetic processing based on the data specified in the program and the data input from the interface. output a signal (computation result) to Note that the GNSS receivers 51 and 81 can also have hardware similar to that of the controller 40 .
  • the storage device 56 may be a device independent of the controller 40 .
  • the controller 40 is connected to the GNSS receiver 51, the attitude sensors 75A, 75B, 75C, 23, the monitor 60, and the wireless device 7 via interfaces, and is capable of inputting and outputting data.
  • the storage device 56 of the controller 40 stores, for example, construction target surface data 55 defining the position of the construction target surface to be constructed by the hydraulic excavator 1, vehicle body shape and dimension data, and various programs executed by the arithmetic processing unit. is stored.
  • the controller 40 By executing a program stored in the storage device 56, the controller 40 performs a work device position/attitude calculation unit 41, a positioning result input unit 42, a satellite position extraction unit 43, an exclusion satellite determination unit 44, a mask range calculation unit 45 , an accuracy calculation unit 46 and a mask selection unit 47 .
  • the positioning result input unit 42 receives the position data of the two GNSS antennas 50A and 50B in the geographic coordinate system calculated by the GNSS receiver 51 and the azimuth data between the two GNSS antennas 50A and 50B (heading data of the upper rotating body 3). ).
  • the mask range calculation unit 45 calculates a plurality of mask ranges used for selecting positioning satellites used for calculation (positioning) of the positions of the two GNSS antennas 50A and 50B by the GNSS receiver 51 .
  • a plurality of mask ranges calculated by the mask range calculator 45 are stored in the storage device 56 within the controller 40 .
  • the mask range calculator 45 calculates (sets) a plurality of mask ranges based on each of the two GNSS antennas 50A and 50B.
  • the mask range calculated here includes a plurality of mask ranges (for example, a plurality of mask ranges with different sizes).
  • a mask range that changes based on the postures of the front work device 6 and the upper swing body 3 (at the time of detection) may be included.
  • the mask range of the GNSS antenna 50A includes a first mask range 21A set to a range in which the front working device 6 can become an obstacle when the GNSS antenna 50A receives satellite signals, and the first mask range 21A.
  • a second mask range 22A set in the range is included.
  • the first mask range 21A is based on the mounting position of the GNSS antenna 50A on the upper rotating body 3, the attitude of the upper rotating body 3 detected via the rotating body attitude sensor 23, and the maximum movable area of the front working device 6. Therefore, it may be set to the maximum range in which the front working device can become an obstacle when the GNSS antenna 50A receives satellite signals.
  • the second mask range 22A consists of the mounting position of the GNSS antenna 50A on the upper revolving body 3, the attitude of the upper revolving body 3 detected via the revolving body attitude sensor 23, and the attitude detected via the attitude sensors 75A-75C. Based on the attitude of the front work device 6, the front work device 6 becomes an obstacle when the attitude of the front work device 6 is calculated (when the attitude is detected) and when the GNSS antenna 50A receives the satellite signal. You can set the range.
  • the first mask range 21A and the second mask range 22A are set on, for example, a geographic coordinate system using a coordinate system based on the GNSS antenna 50A.
  • a coordinate system based on the GNSS antenna 50A may be referred to as a first antenna coordinate system
  • a coordinate system based on the GNSS antenna 50B may be referred to as a second antenna coordinate system.
  • the coordinate system for setting the first mask range 21A and the second mask range 22A is not limited to the geographic coordinate system, and may be set to, for example, a field coordinate system having an origin on the work site.
  • FIG. 3A and 3B are diagrams showing the first antenna coordinate system based on the GNSS antenna 50A.
  • FIG. 3A is a perspective view of the first antenna coordinate system
  • FIG. 3C is a view looking down (plan view)
  • FIG. 3C is a view (side view) of the hydraulic excavator 1 viewed from the +Y-axis direction in the first antenna coordinate system.
  • the first antenna coordinate system is a coordinate system whose origin is the center of the GNSS antenna 50A mounted on the revolving superstructure 3, and is fixed to the GNSS antenna 50A, that is, the revolving superstructure 3.
  • the X-axis of the first antenna coordinate system is a straight line extending along the longitudinal direction of the revolving body 3, and the front of the revolving body 3 is the positive direction.
  • the Y-axis is a straight line extending along the left-right direction of the revolving body 3, and the leftward direction of the revolving body 3 is the positive direction.
  • the Z-axis is orthogonal to the X-axis and the Y-axis, with upward being the positive direction.
  • the + direction of the rotation angle (roll angle, pitch angle, yaw angle (heading)) about each coordinate axis X, Y, Z is written in the figure.
  • the heading shown in FIG. 3B is the direction (azimuth) in which the revolving body 3 and the front work device 6 are facing, and is the line of the orthogonal projection of the X-axis of the first antenna coordinate system onto the horizontal plane and the true north. It is represented by the angle formed by In this embodiment, the true north direction is defined as 0 degrees, and the clockwise direction when looking down on the X-axis from the vertical direction is defined as positive. This defines heading between 0 and 360 degrees. That is, the heading matches the azimuth calculated by the GNSS receiver 51 .
  • the coordinate values in the geographic coordinate system consist of latitude, longitude and ellipsoidal height
  • the coordinate values in the planar rectangular coordinate system, geocentric rectangular coordinate system and field coordinate system are three-dimensional rectangular coordinate systems consisting of X, Y, Z coordinates, etc.
  • Geographical coordinate system coordinate values can be converted into a three-dimensional orthogonal coordinate system such as a planar rectangular coordinate system using the Gauss-Krugel conformal projection method or the like.
  • the planar rectangular coordinate system, the geocentric rectangular coordinate system, and the field coordinate system can be mutually transformed by using affine transformation or Helmert transformation.
  • the boom angle is based on a straight line extending horizontally from the center of the boom pin, which is the rotation axis of the boom 6A (the X-axis of the vehicle body (similar to the X-axis in Figure 3A)). It is the angle of rotation of a straight line 71 passing through the center of the arm pin, which is the driving axis.
  • the arm angle is the angle by which a straight line 72 passing through the center of the arm pin and the center of the bucket pin, which is the rotation axis of the bucket 6C, rotates with the straight line 71 as a reference.
  • the bucket angle is the angle by which a straight line 73 passing through the center of the bucket pin and the tip of the bucket rotates with respect to the straight line 72 .
  • the boom angle, arm angle, and bucket angle are detected by a boom orientation sensor 75A, an arm orientation sensor 75B, and a bucket orientation sensor 75C, respectively.
  • the mask range calculator 45 calculates the boom angle, arm angle, and bucket angle based on detection signals from the boom orientation sensor 75A, arm orientation sensor 75B, and bucket orientation sensor 75C. Next, using the calculated three angles, the three-dimensional model of the front work device 6 stored in the storage device 56 is appropriately rotated and translated, and the attitude of the three-dimensional model is changed to that of the actual front work device 6. position. In addition, the mask range calculation unit 45 calculates the roll angle and pitch angle, which are the inclination angles of the revolving structure 3, based on the detection signal of the revolving structure posture sensor 23, and creates a three-dimensional model of the front working device 6 in the same posture as the actual one.
  • the values of the calculated roll angle and pitch angle and the value of the heading (azimuth of the revolving body 3) calculated by the GNSS receiver 51 are added to rotate.
  • the coordinates of the center of the GNSS antenna 50A in the geographic coordinate system be (X0, Y0, Z0)
  • the coordinates of an arbitrary point Pn on the front working device 6 be (Xn, Yn, Zn).
  • n is a natural number
  • its maximum value is the number of vertices defining the three-dimensional model of the front working device 6 .
  • Equation 1 a vector V from the center of the GNSS antenna 50A to an arbitrary point Pn on the front working device 6 can be expressed by Equation 1 below.
  • X'n Xn-X0
  • Y'n Yn-Y0
  • Z'n Zn-Z0.
  • Coordinates (X'n, Y'n, Z'n) indicating the vector V are coordinates on the first antenna coordinate system.
  • a vector Vxy obtained by projecting this vector V onto the XY plane in the first antenna coordinate system can be expressed by Equation 2 below.
  • the angle from the north direction to the vector Vxy on the XY plane is represented as the azimuth angle ⁇ an (see FIG. 3B)
  • the angle is the coordinate values X′n, Y′n and the heading (heading angle) on the first antenna coordinate system. It can be represented by the following formula 3 using and.
  • the calculation formula for the azimuth angle ⁇ an is divided into five cases according to the combination of X'n and Y'n.
  • the angle formed by the vector Vxy and the vector V is represented as an elevation angle ⁇ en (see FIG. 3C)
  • the angle can be obtained by the following equation using the coordinate values X'n, Y'n, and Z'n on the first antenna coordinate system. 4 can be calculated. It should be noted that in Expression 4 below, the arithmetic expression for the elevation angle ⁇ en is divided into five cases according to the combination of X'n, Y'n, and Z'n.
  • the azimuth angle ⁇ an and the elevation angle ⁇ en can be calculated from the coordinate values of the point Pn on the front working device 6 in the first antenna coordinate system. Since the dimensions of the front work device 6 (boom 6A, arm 6B, bucket 6C) are already known (for example, they are stored in the storage device 56 in advance as vehicle body shape dimension data (details will be described later)), the attitude sensors 75A-75C If the posture (boom angle, arm angle, bucket angle) of the front work device 6 is specified by , the azimuth angle ⁇ an and elevation angle ⁇ en of an arbitrary point Pn on the front work device 6 in that posture can be calculated.
  • the mask range is defined by a combination of the range of the azimuth angle ⁇ an and the range of the elevation angle ⁇ en. That is, the mask range of this embodiment is defined by four parameters, two azimuth angles ⁇ an and two elevation angles ⁇ en.
  • FIG. 5 is a diagram showing an example of the first mask range 21A and the second mask range 22A set on the GNSS antenna 50A on a sky plot (satellite map).
  • the sky with the GNSS antenna 50A as a reference is represented by two-dimensional coordinates with two parameters, azimuth angle and elevation angle, and the mask range is represented in gray with dots.
  • the center of the circle indicates the center of the GNSS antenna 50A
  • the circumferential direction of the circle indicates the azimuth angle
  • the radial direction of the circle indicates the elevation angle.
  • a plurality of circles containing alphabets G, R and two-digit numbers in the figure respectively indicate the positions of positioning satellites captured by the GNSS receiver 51, and alphabets G, R and two-digit numbers represent the respective positioning satellites. number (satellite number).
  • FIG. 5A An example of the first mask range 21A is shown in FIG. 5A.
  • the azimuth angle is set to a range from ⁇ 1 to ⁇ 2 (where ⁇ 2> ⁇ 1) and the elevation angle is set to a range from ⁇ 3 to ⁇ 4 (where ⁇ 4> ⁇ 3).
  • ⁇ 1, ⁇ 2, ⁇ 3, and ⁇ 4 are determined based on the maximum movable range of the front work device 6.
  • ⁇ 1 is the minimum possible value of the azimuth angle ⁇ en when the revolving structure 3 operates the front working device 6 within the maximum movable range
  • ⁇ 2 is the maximum possible value of the azimuth angle ⁇ en under the same conditions. be.
  • ⁇ 3 is the minimum possible value of the elevation angle ⁇ en under the same conditions
  • ⁇ 4 is the maximum possible value of the elevation angle ⁇ en under the same conditions.
  • the first mask range 21A and the four parameters defining it are fixed as long as the excavator 1 (revolving body 3) does not perform a revolving motion or a traveling motion.
  • the second mask range 22A in this figure is set to a range of azimuth angles greater than or equal to ⁇ 1 and less than or equal to ⁇ 2 (where ⁇ 2> ⁇ 1) and an elevation angle greater than or equal to ⁇ 3 and less than or equal to ⁇ 4 (where ⁇ 4> ⁇ 3).
  • ⁇ 1, ⁇ 2, ⁇ 3, and ⁇ 4 are the postures of the front work device 6 at time t when the postures of the front work device 6 and the swing body 3 (boom angle, arm angle, bucket angle, pitch angle, roll angle, heading) are calculated. It is determined based on the range covering the sky of the GNSS antenna 50A.
  • the mask range calculation unit 45 utilizes the above equations (3) and (4) to calculate the attitudes of the front working device 6 and the revolving body 3 (boom angle, arm A combination of azimuth angle ⁇ an and elevation angle ⁇ en is obtained for all vertices Pn on the three-dimensional model of the front working device 6 at a certain time t at which the angle, bucket angle, pitch angle, roll angle, heading) have been calculated.
  • the combinations of azimuth angle ⁇ an and elevation angle ⁇ en at time t are ( ⁇ a1(t), ⁇ e1(t)), ( ⁇ a2(t), ⁇ e2(t)), ( ⁇ a3(t), ⁇ e3(t) ), . .
  • the mask range calculator 45 selects the largest azimuth angle ⁇ a_max(t), the largest elevation angle ⁇ e_max(t), the smallest azimuth angle ⁇ a_min(t), and the smallest elevation angle ⁇ e_min(t) from the obtained combinations.
  • the second mask range 22A and the four parameters defining it are calculated, for example, each time the posture of the front working device 6 is changed.
  • the second mask range 22A is inherently smaller than the first mask range 21A and included in the first mask range 21A (see FIG. 5B).
  • the heading and ⁇ 1 could be expressed with the same value, but in the case of FIG. 5B, the heading and ⁇ 1 are different. Therefore, the heading is not necessarily ⁇ a_min(t).
  • the GNSS antenna can be detected by running or turning the revolving superstructure 3. It moves with the GNSS antenna 50A even when the position of 50A on the geographic coordinate system is changed. That is, the mask ranges 21A and 22A are not fixed on the basis of the geographic coordinate system, but are fixed on the first antenna coordinate system shown in FIG. 3A.
  • the mask range is defined by four parameters.
  • the contour of the mask range 21 may be set.
  • the accuracy calculation unit 46 selects positioning satellites by using each of the plurality of mask ranges (for example, the first mask range 21A and the second mask range 22A) calculated by the mask range calculation unit 45. In the part that calculates (estimates) the positioning accuracy of the GNSS receiver 51 in a plurality of cases including when positioning satellites are selected without using any of the mask ranges (when satellite selection by mask ranges is interrupted) be.
  • the numerical value of N of the GNSS antenna 50B may differ from the numerical value of N of the GNSS antenna 50A.
  • N may also be referred to as the number of selected patterns in the mask range.
  • the positioning accuracy of the GNSS receiver 51 is the accuracy of the positioning results of the GNSS antennas 50A and 50B by the GNSS receiver 51, and may be simply referred to as "positioning accuracy” below.
  • This positioning accuracy includes data related to positioning satellites that can be selected after using each mask range (the number and arrangement of positioning satellites that can be used for positioning (for example, PDOP (Position Dilution of Precision)) and can be calculated from the SN ratio, level, etc. of satellite signals received from various positioning satellites (hereinafter sometimes referred to as "satellite-related data").
  • an index value indicating positioning accuracy is calculated from satellite-related data, and accuracy is evaluated based on the magnitude of the index value.
  • the satellite-related data may be received from the GNSS receiver 51, or may be calculated by the controller 40 based on satellite signals received by the GNSS antennas 50A and 50B.
  • the controller 40 may transmit a list of satellites to be excluded (excluded satellite list, which will be described later) when used to calculate satellite-related data.
  • the accuracy calculation unit 46 uses a certain mask range ("mask range” here includes cases where no mask is used) among a plurality of mask ranges set for each GNSS antenna 50.
  • the satellite-related data of satellites that can be selected later are obtained from the information of the remaining satellites, excluding the satellites included in the certain mask range, among the plurality of satellites input from the satellite position extraction unit 43 .
  • the accuracy calculator 46 calculates the positioning accuracy when the first mask range 21A is used, when the second mask range 22A is used, and when none of the mask ranges is used.
  • the positioning accuracy (first positioning accuracy) when selecting a positioning satellite using the first mask range 21A (first case) is the first This is calculated by the accuracy calculator 46 based on the satellite-related data of the satellites excluding the satellites (G17, G19, R12, R22) included in the mask range 21A.
  • the positioning accuracy (second positioning accuracy) when selecting a positioning satellite using the second mask range 22A (second case) is the second mask range from all satellites on the sky plot captured by the GNSS receiver 51.
  • 22A is calculated by the accuracy calculator 46 based on the satellite-related data of the satellites excluding the satellite (G17) included in 22A.
  • the positioning accuracy (third positioning accuracy) when selecting a positioning satellite without using any mask range (case 3) depends on the satellite-related data of all satellites on the sky plot supplemented by the GNSS receiver 51. Based on this, the accuracy calculation unit 46 calculates.
  • the positioning accuracy for N cases calculated by the accuracy calculation unit 46 is output to the mask selection unit 47 .
  • the mask selection unit 47 selects the case with the best positioning accuracy from among the N cases of positioning accuracy calculated by the accuracy calculation unit 46 .
  • the positioning accuracy is the best, not only is one selected from a plurality of mask ranges set for each of the GNSS antennas 50A and 50B, but there are cases where the mask range is not used.
  • Information on the mask range selected by the mask selection unit 47 is output to the exclusion satellite determination unit 44 .
  • the satellite position extraction unit 43 extracts the positions (elevation angle and azimuth angle in the geographic coordinate system) of multiple satellites whose satellite signals are captured by the GNSS receiver 51, and outputs them to the exclusion satellite determination unit 44 and the accuracy calculation unit 46. do.
  • the excluded satellite determination unit 44 is based on the mask range selected by the mask selection unit 47 and the position (elevation angle/azimuth angle) of the positioning satellite captured by the GNSS receiver 51 input from the satellite position extraction unit 43. , determines the excluded satellites that the GNSS receiver 51 does not use for the positioning calculation. Specifically, the excluded satellite determination unit 44 determines the satellites located in the mask range selected by the mask selection unit 47 from among the plurality of satellites extracted by the satellite position extraction unit 43 as excluded satellites, and determines the excluded satellites. Output the list of satellites to the GNSS receiver 51 .
  • the GNSS receiver 51 acquires the list of excluded satellites output from the excluded satellite determination unit 44, and determines the list of satellites excluded from the list from among the plurality of positioning satellites capable of acquiring satellite signals at that time. Based on the satellite signals, the position of at least one GNSS antenna 50 in the geographic coordinate system and the orientation between the two GNSS antennas 50A, 50B (orientation of the upper revolving body 3) are calculated. As a result, the GNSS receiver 51 selects at least one GNSS antenna in the geographical coordinate system based on the satellite signal of the positioning satellite selected when the positioning accuracy is the best among the N positioning accuracies calculated by the accuracy calculation unit 46. 50 and the azimuth between the two GNSS antennas 50A and 50B. The calculated position and orientation are input to the positioning result input unit 42 .
  • the position measurement results output from the GNSS receiver 51 are coordinates of at least one coordinate system of a planar rectangular coordinate system, a geocentric rectangular coordinate system, and a field coordinate system. It may be possible to output the value.
  • the work device position/orientation calculation unit 41 calculates the position of the GNSS antennas 50A and 50B and the orientation of the upper swing body 3 input from the positioning result input unit 42 and the outputs of the plurality of orientation sensors (75A, 75B, 75C, and 23). Based on the calculated angle values of the front members 6A, 6B, 6C, the inclination angle (pitch angle and roll angle) of the upper revolving body 3, and the vehicle body shape dimension data stored in the storage device 56, the work device 6 (for example, the tip position and attitude of the bucket 6C in the field coordinate system) are calculated.
  • the body shape data used to calculate the position and posture of the work device 6 include, for example, the length between two pins located at both ends of the boom 6A (boom pin length LB (see FIG. 4)), the arm 6B, the length between the two pins located at both ends of the bucket 6B (length between the arm pins), the length between the tip of the bucket 6C and the bucket pin (bucket tip length), and the two GNSS antennas 50A and 50B on the upper rotating body 3 GNSS mounting offset azimuth defined by the azimuth angle in the vehicle body coordinate system, the mounting position (antenna mounting position) of the two GNSS antennas 50A and 50B in the upper swing body 3 in the vehicle body coordinate system, and the two GNSS antennas There is a maximum movable area of the front work device 6 based on each of 50A and 50B.
  • the monitor 60 is based on the position and attitude data of the working device 6 in the field coordinate system calculated by the working device position/posture calculation unit 41 and the position data of the work target surface in the field coordinate system stored in the storage device 56. It is possible to display the positional relationship between the working device 6 and the work target surface calculated by This display allows the operator to easily grasp the position and attitude of the working device 6 with respect to the work target surface.
  • the position and orientation data of the working device 6 in the field coordinate system calculated by the working device position/posture calculation unit 41 and the position data of the work target surface in the field coordinate system stored in the storage device 56 are processed by the controller 40.
  • machine control is executed to restrict the operation of the work device 6 (that is, the operation of the actuators 11A, 11B, 11C that drive the front members 6A, 6B, 6C) so that the tip of the bucket does not exceed the target surface for construction.
  • the mask range calculator 45 defines the mask range using the azimuth angle and the elevation angle, the mask range may be defined using only the azimuth angle.
  • FIG. 5A represents a sky plot when the front work device 6 obstructs the sky view of the GNSS antenna 50A as an obstacle the most.
  • the first mask range 21A and the second mask range 22A are the same range, and can be collectively represented as the first mask range 21A.
  • the boom 6A is lower than in the case of FIG. 5A, and the range in which the front working device 6 blocks the sky view of the GNSS antenna 50A as an obstacle is narrower than in the case of FIG. 5A. .
  • the first mask range 21A and the second mask range 22A are different, the number of satellites selected when using the respective mask ranges 21A and 22A also changes. If the first mask range 21A is applied to positioning in the case of FIG. 5B, the first mask range 21A excludes a total of four satellites, R12, G19, R22, and G17, from satellites that can be used during positioning. However, if the second mask range 22A is applied to positioning in the case of FIG.
  • the accuracy calculation unit 46 extracts satellites that can be used for positioning when each mask range is applied (however, this may include cases where none of the mask ranges are used), and the satellite-related data of the extracted satellites. Based on this, the positioning accuracy when using each mask range is calculated.
  • the mask selection unit 47 compares a plurality of positioning accuracies calculated by the accuracy calculation unit 46, determines a mask range that produces the best positioning accuracy, and uses the mask range to perform positioning of the GNSS antenna 50 by the GNSS receiver 51. is done.
  • FIG. 6 is a flowchart of positioning processing of the GNSS antenna 50A by the controller 40 and the GNSS receiver 51 according to this embodiment. This processing flow is repeatedly calculated at regular intervals (for example, 100 ms).
  • regular intervals for example, 100 ms.
  • the same processing is also performed for the GNSS antenna 50B, and the azimuth between the two GNSS antennas 50A and 50B is determined based on the positioning results of both GNSS antennas 50A and 50B. shall be calculated.
  • step 51 the controller 40 (satellite position extraction unit 43 ) acquires from the GNSS receiver 51 data (including position data, etc.) of a plurality of satellites whose satellite signals are captured by the GNSS receiver 51 .
  • the controller 40 calculates the attitude (pitch angle, roll angle) of the revolving superstructure 3 based on the detection signal from the revolving superstructure position sensor 23 .
  • the controller 40 calculates the attitude of the front working device 6 based on the detection signals of the three attitude sensors 75A, 75B, and 75C.
  • step 54 the controller 40 (mask range calculation unit 45) calculates the posture of the swing body 3 calculated in step 52, the maximum movable range of the front working device 6 stored in the storage device 56, and the range calculated in the previous cycle.
  • the first mask range 21A is calculated based on the azimuth between the antennas 50A and 50B and stored in the storage device 56.
  • controller 40 calculates the attitude of the revolving body 3 calculated in step 52, the attitude of the front work device 6 calculated in step 53, and the position of the GNSS antenna 50A calculated in the previous cycle. , and the heading of the revolving body 3 calculated in the previous cycle, the second mask range 22A included in the first mask range 21A is calculated.
  • step 55 the controller 40 (accuracy calculation section 46) selects the remaining satellites by excluding the satellites located in the first mask range 21A.
  • a first positioning accuracy, which is accuracy, is calculated.
  • the satellite-related data of the satellites supplemented in the second case of excluding the satellites located in the second mask range 22A and selecting the remaining satellites is calculated based on the satellite-related data in the third case (which is assumed to be the same type as the The third positioning accuracy, which is the positioning accuracy in the third case, is calculated based on the same type as that used when calculating the first positioning accuracy.
  • the controller 40 compares the first positioning accuracy and the second positioning accuracy calculated at step 55 . If the first positioning accuracy is better than the second positioning accuracy, the process proceeds to step 57. If the first positioning accuracy is the same as or worse than the second positioning accuracy, the process proceeds to step 59.
  • step 57 the controller 40 (mask selection unit 47) determines the satellites located in the first mask range 21A from among the plurality of satellites whose positions are acquired in step 51 as excluded satellites, and lists the excluded satellites ( excluded satellite list) to the GNSS receiver 51.
  • FIG. 7 is a diagram showing an example of changes over time in the angle of the boom 6A (boom angle) and PDOP.
  • the angle of the boom 6A boost angle
  • PDOP PDOP
  • FIG. 7 shows a situation in which excavation and dumping (rotating motion) are repeatedly performed.
  • the satellites that are shielded will change greatly due to changes in the posture of the boom 6A and changes in the turning angle (orientation of the turning body 3).
  • the PDOP may increase before and after the boom is raised, that is, the positioning accuracy may decrease. Therefore, in such a case, using the wider first mask range 21A to suppress variations in the mask range may improve the positioning accuracy more than using the second mask range 22A. .
  • the GNSS receiver 51 excludes the satellites included in the excluded satellite list of step 57 from the satellites that transmitted data to the controller 40 in step 51, and based on the satellite signals of the remaining satellites, the GNSS antenna 50A. A position is calculated and it is output to the controller 40 as a positioning result of the GNSS antenna 50A.
  • the controller 40 compares the second positioning accuracy and the third positioning accuracy calculated at step 55 .
  • the process proceeds to step 60, and when the second positioning accuracy is the same as or worse than the third positioning accuracy, the process proceeds to step 62.
  • step 60 the controller 40 (mask selection unit 47) determines the satellites located in the second mask range 22A from among the plurality of satellites whose positions are acquired in step 51 as excluded satellites, and lists the excluded satellites ( excluded satellite list) to the GNSS receiver 51.
  • the GNSS receiver 51 excludes the satellites included in the excluded satellite list in step 60 from the satellites that transmitted data to the controller 40 in step 51, and based on the satellite signals of the remaining satellites, the GNSS antenna 50A. A position is calculated and it is output to the controller 40 as a positioning result of the GNSS antenna 50A.
  • step 62 the GNSS receiver 51 calculates the position of the GNSS antenna 50A based on the satellite signal of the satellite that transmitted the data to the controller 40 in step 51, and outputs it to the controller 40 as the positioning result of the GNSS antenna 50A.
  • FIG. 8 is a diagram showing the satellite arrangement at a certain time and the second mask range 22A.
  • the satellite G28 only one satellite (G28) exists in the azimuth of the front working device 6, and the number of satellites in the azimuth opposite to the azimuth of the front (180 degrees plus) is small.
  • the satellite G28 is excluded by the second mask range 22A, the number of satellites is insufficient, and the positioning accuracy tends to deteriorate.
  • the use of the second mask range 22A rather lowers the positioning accuracy, so the positioning accuracy is better when neither mask range is used.
  • the controller 40 when a plurality of mask ranges are set, the controller 40 performs positioning when each mask range is used based on satellite-related data when each mask range is used. Accuracy is calculated, and one mask range to be used for positioning of the GNSS antennas 50A and 50B is selected from among the plurality of mask ranges based on the calculated positioning accuracy. As a result, the positioning accuracy can be improved more than when the mask range is selected based on the attitude of the front working device 6 .
  • the controller 40 calculates the positioning accuracy in each case based on each satellite-related data when using a plurality of mask ranges of different sizes (including cases where mask ranges are not set), From among them, one mask range is selected for the case where the positioning accuracy is the best.
  • the controller 40 calculates the positioning accuracy in each case based on each satellite-related data when using a plurality of mask ranges of different sizes (including cases where mask ranges are not set), From among them, one mask range is selected for the case where the positioning accuracy is the best.
  • the satellite-related data in the selection of the mask range, it is possible to suppress the deterioration of the positioning accuracy more than when the mask range (in the above embodiment, the second mask range 22A) is changed along with the attitude of the front working device 6. Therefore, the working accuracy of the hydraulic excavator 1 can be improved as a result.
  • a first mask range 21A that covers the entire movable area of the front working device 6 and a second mask range that covers an area that is hidden by the posture of the front working device 6 at that time. 22A are included.
  • the attitude of the front working device 6 is taken into account in the selection of the mask range and the calculation of the positioning accuracy, and the influence of the shielding of the positioning signal by the front working device 6 on the positioning accuracy can be reduced.
  • step 56 if it is determined in step 56 that the first positioning accuracy is better than the second positioning accuracy, even if determination processing for comparing the first positioning accuracy and the third positioning accuracy is added. good.
  • this determination processing if it is determined that the first positioning accuracy is better, the process proceeds to step 57, and if it is determined that the third positioning accuracy is better, the same processing as in step 62 ( That is, the positioning process) may be performed without using any mask range.
  • the number of usable mask ranges is not limited as long as each mask range is different.
  • a range larger than the first mask range may be set.
  • the mask range is set in consideration of the front movement range in the same operation, and the operations performed by the hydraulic excavator are set. It is also possible to use the mask range that is estimated and set for the work. The mask range in the latter case is assumed to be smaller than the first mask range 21A and larger than the second mask range 22A.
  • the present invention is not limited to the above-described embodiments, and includes various modifications within a scope that does not deviate from the gist of the present invention.
  • the present invention is not limited to those having all the configurations described in the above embodiments, but also includes those with some of the configurations omitted. Also, it is possible to add or replace part of the configuration according to one embodiment with the configuration according to another embodiment.
  • each configuration related to the controller 40 and the functions and execution processing of each configuration are implemented partially or entirely by hardware (for example, logic for executing each function is designed by an integrated circuit).
  • the configuration related to the controller 40 may be a program (software) that implements each function related to the configuration of the controller 40 by being read and executed by an arithmetic processing unit (for example, a CPU).
  • Information related to the program can be stored, for example, in a semiconductor memory (flash memory, SSD, etc.), a magnetic storage device (hard disk drive, etc.), a recording medium (magnetic disk, optical disk, etc.), or the like.
  • control lines and information lines have been shown as necessary for the description of the embodiments, but not necessarily all the control lines and information lines related to the product does not necessarily indicate In reality, it can be considered that almost all configurations are interconnected.

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • Civil Engineering (AREA)
  • Mining & Mineral Resources (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mechanical Engineering (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Component Parts Of Construction Machinery (AREA)
  • Operation Control Of Excavators (AREA)

Abstract

作業機械に備えたコントローラには,受信機によるアンテナの位置の演算に利用する測位衛星を選択するための複数のマスク範囲であって,作業装置及び旋回体の姿勢に基づいてアンテナを基準に設定される複数のマスク範囲が記憶されている。コントローラは,複数のマスク範囲のそれぞれを利用して選択される測位衛星の衛星関連データと,複数のマスク範囲のいずれも利用することなく選択される測位衛星の衛星関連データとに基づいて複数の測位精度をそれぞれ演算する。受信機は,演算した前記複数の測位精度の中で最も測位精度が良い場合に選択される測位衛星の衛星信号に基づいてアンテナの位置を演算する。

Description

作業機械
 本発明は測位衛星システムを用いた位置検出が行われる作業機械に関する。
 近年,施工現場において情報化施工の導入が進められている。情報化施工とは,調査,設計,施工,検査,管理などからなる一連の建設工程のうち,施工に注目して,電子データとICT(Information and Communication Technology:情報通信技術)の活用により施工の高効率化を実現するシステムである。情報化施工に対応する機械としては,車体位置やフロント作業装置(作業装置と称すこともある)の位置及び姿勢を施工目標面の位置データとともにモニタに表示するガイダンス機能や,バケットが施工目標面を掘り過ぎないようにフロント作業装置を制御するマシンコントロール機能を搭載する油圧ショベルに代表される作業機械が知られている。このような情報化施工に対応した作業機械は,3次元座標データを持つ情報化施工データをもとに,オペレータに対して情報提示し,作業支援,運転支援する機能を提供する。例えば油圧ショベルのマシンガイダンスにおいては,車体の位置及び姿勢のデータやフロント作業装置の姿勢データからバケット先端位置が演算され,施工目標面に対するバケットの位置がモニタを介してオペレータに提示される。
 この種の油圧ショベルには,グローバル座標系(地理座標系)における上部旋回体(車体)の位置を演算する為に,上部旋回体に取り付けられた測位アンテナを介して測位衛星からの測位信号を受信することで上部旋回体の位置を演算する衛星測位システム(例えばGNSS(Global Navigation Satellite System:全球測位衛星システム))が搭載されることがある。しかし,油圧ショベルにおいては,ブーム,アーム,バケット等のフロント部材が,衛星測位システムの測位アンテナよりも上方に存在することがあり,これらが直線的なルートによる測位信号の受信を妨害することがある。このような場合,測位アンテナは,測位信号をマルチパスと呼ばれる回折波や反射波として受信する可能性が高い。回折波や反射波を受信して測位演算に使用すると測位結果に誤差を含む可能性が高くなる。
 そこで,マルチパスの影響の低減を試みた技術として,例えば特許文献1のものがある。特許文献1には,GPS受信機の周囲の電波障害物の配置などに基づいて定めたGPS衛星の選択規則を表すマスク範囲(マスク情報)を所定地域毎にデータベースに記憶し,当該データベースからGPS信号の受信地域に対応するマスク範囲を取得し,天空を航行する複数のGPS衛星のうち当該マスク範囲の外に位置するGPS衛星の中から測位に利用するGPS衛星を選択するGPS受信機の制御装置が開示されている。
特開2004-184121号公報
 特許文献1に記載の技術は,移動しない電波障害物のみを想定しており,油圧ショベルのフロント作業装置のように移動し得る電波障害物が測位アンテナの周囲に存在する場合のマスク範囲には触れられていない。つまり,作業機械の分野における衛星測位精度の向上には,移動し得る電波障害物であるフロント作業装置の姿勢を考慮したマスク範囲を設定することが重要である。
 そこでフロント作業装置の姿勢の変化に応じて衛星選択に利用するマスク範囲を都度変更すると(または変化させると)測位精度の低下を抑制できそうだが話はそう単純ではない。まず,作業機械ではフロント作業装置を動作させて作業を行うため,フロント作業装置の姿勢変更(例えばブームの上げ下げ動作)が頻繁に行われる。その姿勢変更に合わせてマスク範囲を変更すると,測位演算に利用される衛星が頻繁に切り替わってしまい,マスク範囲を変更した方が却って測位精度が低下する可能性がある。また,測位演算に利用可能な衛星が極端に少ない場合にフロント作業装置の姿勢に応じたマスク範囲を利用すると,その少ない衛星でさえも測位に利用できない時間帯が出てくるため,マスク範囲を利用しない方が却って測位精度の低下を抑えられる場合もある。
 このように電波障害物たり得るフロント作業装置の姿勢が作業中に頻繁に変化し得る作業機械では,フロント作業装置が存在する領域に位置する衛星を測位演算の際に単純に除外するだけでは作業中の測位精度の低下を抑制できない場合がある。
 なお,上記では主に衛星配置に起因した測位精度低下について触れたが,上記の他に例えば衛星数や衛星信号のSN比(信号とノイズの比)に起因した測位精度低下も起こりうる。
 本発明は上記の事情に鑑みてなされたものであり,その目的は,作業装置の動作に伴う衛星測位精度の低下を抑制し得る作業機械を提供することにある。
 本願は上記課題を解決する手段を複数含んでいるが,その一例を挙げるならば,下部走行体と,前記下部走行体の上に旋回可能に取り付けられた上部旋回体と,前記上部旋回体に取り付けられた多関節型の作業装置と,前記作業装置及び前記上部旋回体の姿勢を検出するための複数の姿勢センサと,前記上部旋回体に取り付けられ,複数の測位衛星からの衛星信号を受信するための測位アンテナと,前記測位アンテナで受信される衛星信号に基づいて前記測位アンテナの位置を演算する受信機と,前記複数の姿勢センサの検出信号に基づいて,前記作業装置の姿勢と前記上部旋回体の姿勢とを演算するコントローラとを備えた作業機械において,前記コントローラには,前記受信機による前記測位アンテナの位置の演算に利用する測位衛星を選択するために,前記作業装置及び前記上部旋回体の姿勢に基づいて前記測位アンテナを基準に設定される複数のマスク範囲が記憶されており,前記コントローラは,前記複数のマスク範囲のそれぞれを利用して測位衛星を選択した場合のそれぞれの測位精度と,前記複数のマスク範囲のいずれも利用することなく測位衛星を選択した場合の測位精度とを含む複数の測位精度を,選択した前記測位衛星の衛星関連データに基づいて演算し,前記受信機は,前記複数の測位精度のうち最も測位精度が良い場合に選択した前記測位衛星の衛星信号に基づいて前記測位アンテナの位置を演算することを特徴とする。
 本発明によれば,フロント作業装置の動作とともにフロント作業装置の姿勢を反映したマスク範囲が変化しても,最も測位精度が良いマスク範囲が選択される。これにより作業装置の動作に伴う測位精度の低下が抑制でき,作業機械による作業精度を向上させることができる。
本発明の実施形態に係る油圧ショベル1及びGNSS基準局8の側面図。 図1の油圧ショベル1に搭載された車載コントローラ40の機能ブロック図。 GNSSアンテナ50Aを基準とする第1アンテナ座標系の斜視図。 第1アンテナ座標系におけるZ軸上から油圧ショベル1を見下ろした図(平面図)。 第1アンテナ座標系における+Y軸方向から油圧ショベル1を見た図(側面図)。 フロント作業装置6の姿勢を規定する角度の説明図。 GNSSアンテナ50Aに設定される第1マスク範囲21Aの一例を衛星配置図上に示した図。 GNSSアンテナ50Aに設定される第2マスク範囲22Aの一例を衛星配置図上に示した図。 本実施形態に係るコントローラ40及びGNSS受信機51によるGNSSアンテナ50Aの測位処理のフローチャート。 ブーム6Aの角度(ブーム角度)とPDOPの時間変化の一例を示した図。 或る時刻の衛星配置と第2マスク範囲22Aを示す図。
 以下,本発明の実施の形態について図面を用いて説明する。以下の実施の形態は,作業機械としてクローラ式の油圧ショベルに本発明を適用したものであり,バケット先端と施工目標面の位置関係を運転室内のモニタに表示するマシンガイダンス機能と,バケット先端が施工目標面を超えないように作業装置の動作(すなわちフロント部材を駆動するアクチュエータの動作)に制限をかけるマシンコントロール機能とを有している。なお,各図において同じ部分には同じ符号を付し,重複した説明は適宜省略するものとする。
 <対象機械>
 図1は本発明の実施形態に係る油圧ショベル1及びGNSS基準局8の側面図である。この図に示す油圧ショベル1は,クローラ式の走行体(下部走行体)2と,走行体2の上部に旋回可能に取り付けられた旋回体(上部旋回体)3と,一端(基端)が旋回体3の前方に取り付けられた多関節型のリンク機構よりなるフロント作業装置(単に「作業装置」と称することもある)6とを備えている。図中の符号30は地面を表す。
 フロント作業装置6は,一端が旋回体3に連結されたブーム6Aと,一端がブーム6Aの他端に連結されたアーム6Bと,一端がアーム6Bの他端に連結されたバケット6Cとを有しており,これら各フロント部材6A,6B,6Cは,それぞれ上下方向に回動するように構成されている。
 また,各フロント部材6A,6B,6Cの回動を行うアクチュエータ(油圧シリンダ)として,ブームシリンダ11A,アームシリンダ11B,バケットシリンダ11Cが備えられている。旋回体3は図示しない旋回モータによって旋回中心軸Oを中心に左右に旋回駆動可能である。
 ブーム6A,アーム6B及びバケット6Cは,フロント作業装置6を含む共通の平面上で動作し,以下ではこの平面を動作平面と称することがある。つまり動作平面とは,ブーム6A,アーム6B及びバケット6Cの回動軸に直交する平面であり,例えばブーム6A,アーム6B及びバケット6Cの幅方向の中心(すなわち各フロント部材6A,6B,6Cの回動軸の中心)に設定できる。
 <姿勢センサ>
 油圧ショベル1には,フロント作業装置6と旋回体3の姿勢を検出するための複数の姿勢センサ75A,75B,75C,23が備えられている。本実施形態では各姿勢センサに,角度(または角速度)と加速度を検出可能な慣性計測装置(IMU:Inertial Measurement Unit)を用いている。これら姿勢センサのうち,ブーム6Aにはブーム姿勢センサ75Aが,アーム6Bにはアーム姿勢センサ75Bが,バケット6Cにはバケット姿勢センサ75Cが取り付けられている(図1参照)。また,旋回体3には旋回体姿勢センサ23が取り付けられており(図1参照),それにより旋回体3の傾斜角度(ピッチ角及びロール角),旋回速度及び旋回角度の計測が可能となっている。姿勢センサ75A,75B,75C,23の出力(検出信号)は,接続線を介してコントローラ40に入力されている。なお,フロント作業装置6の姿勢センサとしては,各フロント部材の回動角度を検出する角度センサ(例えば,ポテンショメータやロータリエンコーダ)を用いても良い。またバケット姿勢センサ75Cはバケットでなくバケットリンクに取り付けられていてもよい。
 旋回体3には,オペレータによって操作される複数の操作レバー(図示せず),バケット6Cと施工目標面の位置関係等が表示されるモニタ60が設置された運転席4と,複数の測位衛星(GNSS衛星)から衛星信号を受信するための2つのGNSSアンテナ(測位アンテナ)50A,50Bと,基準局8から送信されるGNSS補正データを受信するための無線機7と,2つのGNSSアンテナ50A,50Bのうち少なくとも1つのGNSSアンテナの地理座標系(グローバル座標系)における位置座標と,2つのGNSSアンテナ50A,50B間の方位(すなわち旋回体3の方位)とを演算するGNSS受信機51と,GNSS受信機51で演算された位置及び方位と,複数の姿勢センサ75A,75B,75C,23の検出信号とに基づいて,フロント作業装置6上の所望の位置座標を演算するコンピュータであるコントローラ40とが備えられている。なお,本実施形態では2つのGNSSアンテナ50A,50Bの位置及び旋回体3の方位角を1つのGNSS受信機で演算する構成を採っているが,2つのGNSSアンテナ50A,50Bのそれぞれに対応する2つのGNSS受信機51A,51Bを搭載する構成を採っても良い。
 <GNSS基準局>
 油圧ショベル1の無線機7に対してGNSS補正データを無線送信するGNSS基準局8について説明する。地理座標系における座標位置が既知であるGNSS基準局8には,複数の測位衛星(GNSS衛星)から衛星信号を受信するためのGNSSアンテナ80と,GNSSアンテナ80で受信された衛星信号(衛星信号には衛星のコード,キャリア,衛星軌道及び衛星信号受信レベル等が含まれる)に基づいてGNSSアンテナ80の地理座標系における位置座標を演算するGNSS受信機81と,GNSSアンテナ80で受信された複数の衛星信号に基づいて無線機7に無線送信するためのGNSS補正データを生成する基準局コントローラ82と,基準局コントローラ82で生成されたGNSS補正データを無線機7に送信する無線機87が備えられている。GNSS基準局アンテナ80に接続されたGNSS受信機81は,基準局コントローラ82を経由して無線機87よりGNSS補正データを無線送信する。無線機7で受信されたGNSS補正データをGNSS受信機51での測位に利用するとセンチメートル級の高精度な測位が可能となる。
 <GNSSアンテナ50>
 2つのGNSSアンテナ50A,50Bは,それぞれマスト(アンテナ支持部材)52a,52bを介して上部旋回体3に固定されており,上部旋回体3の上面にそれぞれ位置し,フロント作業装置6の前後方向に所定の間隔を介して配置されている。
 2本のマスト52a,52bはそれぞれ上部旋回体3の上方でGNSSアンテナ50A,50Bを支持するためのポール状の支持部材である。本実施形態の2本のマスト52a,52bは,GNSSアンテナ50A,50Bと同様に上部旋回体3の上面(第1領域)に配置されている。各マスト52a,52bの基端は上部旋回体3の上面に固定されており,各マスト52a,52bは当該基端から略垂直に伸びている。そして各マスト52a,52bの先端には,中心部が軸方向に膨らんだ略円盤状の外形を有するGNSSアンテナ50A,50Bが取り付けられており,各マスト52a,52bは自身の中心軸心が各GNSSアンテナ50A,50Bの中心軸心を通過するように各アンテナ50A,50Bを支持している。なお,GNSSアンテナ50A,50Bの支持部材は,ポール状のマスト52a,52bに限らず,種々の形状の支持部材による支持が可能である。
 <GNSS受信機51>
 GNSS受信機51は,2つのGNSSアンテナ50A,50Bで受信される複数の衛星信号(衛星のコード,キャリア,衛星軌道及び衛星信号受信レベル等が含まれる)に基づいて,2つのGNSSアンテナ50A,50Bのうち少なくとも1つのGNSSアンテナ(例えば,GNSSアンテナ50B)の地理座標系(グローバル座標系)における位置座標と,2つのGNSSアンテナ50A,50B間の方位(すなわち旋回体3やフロント作業装置6の方位(ヘディングとも称する))とを演算する。
 複数の測位衛星からは送信時刻情報を含んだ電磁波(衛星信号)が送信されている。GNSS受信機51は,各GNSS衛星からの電磁波の受信時刻とその電磁波に含まれた送信時刻とから到達時間差を演算し,その到達時間差を基に各GNSS衛星とGNSSアンテナ50A,50Bとの距離を推測してGNSSアンテナ50A,50Bの位置を算出する。GNSS衛星は精巧な時計を搭載しており,各衛星からの電磁波を復調して得られる到達時間差に電磁波の速度を乗算することにより各GNSS衛星とGNSSアンテナ間の距離が算出される。
 算出した各GNSS衛星と各GNSSアンテナとの距離に誤差が含まれ得る。この誤差は,GNSS衛星とGNSSアンテナ間に存在する電離層や水蒸気によって発生する電磁波の速度変化が方位や仰角が異なる各GNSS衛星の位置毎に異なることや,各GNSS衛星より電磁波で送られる軌道情報が実際の位置と若干異なることや,各GNSS衛星間の時計情報に若干の誤差があること等の要因により発生する。
 このような誤差は,GNSS基準局8から送信されるGNSS補正データを受信して測位するRTK-GNSS(リアルタイムキネマティックGNSS)を利用することで低減できる。例えば,油圧ショベル1の近くに(数km以内)設置した絶対位置が既知の基準局GNSSアンテナ80の測位とGNSS補正データの演算を基準局GNSS受信機81で行い,その補正データを無線機87にてショベル1の受信機51に送信する。そして2つのGNSSアンテナ50A(50B),80間の絶対位置ではなく相対位置(ベクトル)を測定することで誤差を低減することができる。
 GNSS基準局8の無線機87より送信された補正データは,油圧ショベル1に搭載された無線機7で受信されGNSS受信機51に送信される。GNSS受信機51ではGNSSアンテナ50A(移動局)で受信した衛星信号と補正データより得た基準局GNSSアンテナ80の信号を比較演算することにより,基準局GNSSアンテナ80とGNSSアンテナ50A間の相対的な位置(方向と距離)を算出する。このとき,補正情報として基地局アンテナ80が受信した衛星からの衛星信号の搬送波位相情報を送信し,これを移動局アンテナ50Aが受信した衛星信号の搬送波位相情報とGNSS受信機51で比較演算する。これにより数cmオーダーの移動局アンテナ50Aの測位が可能となり,ほぼ一点に収束した高精度の相対測位が可能となる。さらに,前述した補正データのなかに基準局GNSSアンテナ80の位置情報を含めることで,移動局であるGNSSアンテナ50Aの絶対位置を求めることが可能となる。また,基準局GNSSアンテナ80とGNSSアンテナ50Aの距離が近距離(一般的に数km以内)の場合は,前述した誤差要因(電磁波の速度変化,各GNSS衛星間の時計情報誤差)をよく相殺することが可能となる。2つのGNSSアンテナ50A,50B間の方位や,もう1つのGNSSアンテナ50Bの位置についても同様に演算できる。GNSS受信機51はそれぞれのGNSSアンテナ50A,50Bの緯度,経度,ジオイド高さを含むNMEAフォーマットなどでGNSSアンテナ50A,50Bの測位結果を出力可能である。
 ところで,本実施形態のGNSS受信機51の測位対象には2つのGNSSアンテナ50A,50Bが存在するため,一方のGNSSアンテナ50Aを基準局とし他方のGNSSアンテナ50Bを移動局とみなすことができる。このような手法がムービングベース方式である。GNSSアンテナ50Aの受信信号にて生成した補正データをGNSSアンテナ50Bとの相対位置(ベクトル)の測定に利用することで2つのGNSSアンテナ50A,50B間の相対位置(ベクトル)を測定することが可能となる。ムービングベース方式では無線機87から送信される補正データを利用することなく相対位置(ベクトル)を演算可能である。
 また,別の方向算出方法として,基準局GNSSアンテナ80からGNSSアンテナ50AとGNSSアンテナ50Bの位置をそれぞれ演算して,その位置の差分から方向を求める方法もある。そして,このようにして演算した2つのGNSSアンテナ50A,50B間の方向に,ショベル1における2つのGNSSアンテナ50A,50Bの取り付け位置に起因した定数を考慮することにより,上部旋回体2(車体)及びフロント作業装置6の方角(方向)が算出可能である。
 また,本実施形態では基準局GNSSアンテナ80から補正データを無線送信して上部旋回体3やフロント作業装置6の方向を演算するシステムについて説明したが,VRS(仮想基準点方式)や準天頂衛星等の補正データをネットワークで配信するサービスを用いても良い。なお,言うまでもないが,測位精度が許容される範囲であれば,基準局8等からの補正データを利用することなくGNSSアンテナ50A,50Bの測位をしても良い。
 <コントローラ40>
 図2は図1の油圧ショベル1に搭載されたコントローラ40の機能ブロック図である。
 コントローラ40は,GNSS受信機51で演算された2つのGNSSアンテナ50A,50Bの位置と旋回体3の方位(ヘディング)と,複数の姿勢センサ75A,75B,75C,23の検出信号に基づいて,フロント作業装置6を構成する各フロント部材6A,6B,6Cの位置座標を演算するコンピュータである。
 コントローラ40は,演算処理装置(例えばCPU(図示せず)),記憶装置(例えば,ROM,RAM等の半導体メモリ)56,インタフェース(入出力装置(図示せず))を備えており,記憶装置56内に予め保存されているプログラム(ソフトウェア)を演算処理装置で実行し,プログラム内で規定されているデータとインタフェースから入力されたデータに基づいて演算処理装置が演算処理を行い,インタフェースから外部に信号(演算結果)を出力する。なお,GNSS受信機51,81もコントローラ40と同種のハードウェアを備えることができる。また,記憶装置56はコントローラ40から独立した装置としても良い。
 コントローラ40は,インタフェースを介して,GNSS受信機51,姿勢センサ75A,75B,75C,23,モニタ60,及び無線機7と接続されており,データの入出力が可能になっている。
 コントローラ40の記憶装置56には,例えば,油圧ショベル1の施工対象である施工目標面の位置を定義した施工目標面データ55と,車体形状寸法データと,演算処理装置によって実行される各種プログラム等が記憶されている。
 コントローラ40は,記憶装置56内に格納されたプログラムを実行することで,作業装置位置・姿勢演算部41,測位結果入力部42,衛星位置抽出部43,除外衛星決定部44,マスク範囲演算部45,精度演算部46,及びマスク選択部47として機能する。
 (測位結果入力部42)
 測位結果入力部42は,GNSS受信機51で演算される地理座標系における2つのGNSSアンテナ50A,50Bの位置データと,2つのGNSSアンテナ50A,50B間の方位データ(上部旋回体3のヘディングデータ)とを入力する。
 (マスク範囲演算部45)
 マスク範囲演算部45は,GNSS受信機51による2つのGNSSアンテナ50A,50Bの位置の演算(測位)に利用する測位衛星を選択するために用いられる複数のマスク範囲を演算する。マスク範囲演算部45によって演算された複数のマスク範囲はコントローラ40内の記憶装置56に記憶される。マスク範囲演算部45は,2つのGNSSアンテナ50A,50Bのそれぞれを基準にして複数のマスク範囲を演算(設定)する。ここで演算されるマスク範囲には,複数のマスク範囲(例えば、大きさの異なる複数のマスク範囲)が含まれており,その中には演算時(言い換えると姿勢センサ75A-75C,23による姿勢検出時)のフロント作業装置6及び上部旋回体3の姿勢に基づいて変化するマスク範囲が含まれ得る。以下では,主にGNSSアンテナ50Aのマスク範囲の演算や設定について説明するが,GNSSアンテナ50Bについても同様にマスク範囲の演算や設定が行われる。
 (第1マスク範囲21A及び第2マスク範囲22A)
 GNSSアンテナ50Aのマスク範囲には,GNSSアンテナ50Aが衛星信号を受信する際にフロント作業装置6が障害物となり得る範囲に設定された第1マスク範囲21Aと,第1マスク範囲21Aに包含される範囲に設定された第2マスク範囲22Aが含まれる。
 第1マスク範囲21Aは,上部旋回体3におけるGNSSアンテナ50Aの取付位置と,旋回体姿勢センサ23を介して検出された上部旋回体3の姿勢と,フロント作業装置6の最大可動領域とに基づいて,GNSSアンテナ50Aが衛星信号を受信する際にフロント作業装置が障害物となり得る最大範囲に設定しても良い。
 第2マスク範囲22Aは,上部旋回体3におけるGNSSアンテナ50Aの取付位置と,旋回体姿勢センサ23を介して検出された上部旋回体3の姿勢と,姿勢センサ75A-75Cを介して検出されたフロント作業装置6の姿勢とに基づいて,フロント作業装置6の姿勢演算時(姿勢検出時)におけるフロント作業装置6が,GNSSアンテナ50Aが衛星信号を受信する際にフロント作業装置が障害物となる範囲に設定しても良い。
 本実施形態では,第1マスク範囲21A及び第2マスク範囲22Aを,GNSSアンテナ50Aを基準とする座標系を利用して,例えば地理座標系上に設定している。なお,GNSSアンテナ50Aを基準とする座標系を第1アンテナ座標系と称し,同様にGNSSアンテナ50Bを基準とする座標系を第2アンテナ座標系と称することがある。また,第1マスク範囲21A及び第2マスク範囲22Aを設定する座標系は地理座標系に限られず,例えば作業現場上に原点を有する現場座標系に設定しても良い。
 第1アンテナ座標系の定義について説明する。
 図3はGNSSアンテナ50Aを基準とする第1アンテナ座標系を表す図であり,図3Aは第1アンテナ座標系の斜視図,図3Bは第1アンテナ座標系におけるZ軸上から油圧ショベル1を見下ろした図(平面図),図3Cは第1アンテナ座標系における+Y軸方向から油圧ショベル1を見た図(側面図)である。
 図3Aに示すように,第1アンテナ座標系は,旋回体3に搭載されたGNSSアンテナ50Aの中心を原点とする座標系であり,GNSSアンテナ50A,すなわち旋回体3に固定されている。第1アンテナ座標系のX軸は旋回体3の前後方向に沿って延びる直線であり,旋回体3の前方が正の方向である。Y軸は旋回体3の左右方向に沿って延びる直線であり,旋回体3の左方向が正の方向である。Z軸はX軸及びY軸に直交しており,上方が正の方向である。各座標軸X,Y,Zを中心とした回転角(ロール角,ピッチ角,ヨー角(ヘディング))の+方向が図中に記入されている。
 図3B中に示したヘディング(heading)は,旋回体3やフロント作業装置6が向いている方向(方位)であり,第1アンテナ座標系のX軸を水平面に正射影した線と真北の成す角度で表される。本実施形態では,真北の方向を0度と定義し,X軸を鉛直上方から見下ろして時計回り方向が正とする。これによりヘディングは0度から360度の間で規定される。つまり,ヘディングはGNSS受信機51で演算される方位に一致する。
 地理座標系における座標値は緯度,経度及び楕円体高からなり,平面直角座標系,地心直交座標系及び現場座標系の座標値はX,Y,Z座標等からなる3次元直交座標系である。地理座標系座標値はガウス・クリューゲルの等角投影法などを用いて平面直角座標系などの3次元直交座標系に変換可能である。また,平面直角座標系,地心直交座標系及び現場座標系はアフィン変換またはヘルマート変換などを用いることで相互に変換可能である。
 図4を用いてフロント作業装置6の姿勢を規定する角度であって姿勢センサ75A-75Cによって検出される角度(ブーム角,アーム角,バケット角)について説明する。この図において,ブーム角は,ブーム6Aの回動軸であるブームピンの中心から水平に延びる直線(車体X軸(図3AのX軸と同様))を基準として,ブームピンの中心とアーム6Bの回動軸であるアームピンの中心を通る直線71が回転した角度である。アーム角は,直線71を基準として,アームピンの中心とバケット6Cの回動軸であるバケットピンの中心を通る直線72が回転した角度である。バケット角は,直線72を基準として,バケットピンの中心とバケット先端を通る直線73が回転した角度である。ブーム角,アーム角,バケット角は,それぞれ,ブーム姿勢センサ75A,アーム姿勢センサ75B,バケット姿勢センサ75Cによって検出される。
 マスク範囲演算部45は,ブーム姿勢センサ75A,アーム姿勢センサ75B,バケット姿勢センサ75Cの検出信号に基づいて,ブーム角,アーム角,バケット角を演算する。次に演算した3つの角度を利用して,記憶装置56に記憶されたフロント作業装置6の3次元モデルに回転移動と平行移動を適宜加え,その3次元モデルの姿勢を実際のフロント作業装置6の姿勢に合わせる。また,マスク範囲演算部45は旋回体姿勢センサ23の検出信号に基づいて旋回体3の傾斜角であるロール角及びピッチ角を演算し,実際と同じ姿勢のフロント作業装置6の3次元モデルに対して,演算したロール角及びピッチ角の値と,GNSS受信機51で演算されたヘディング(旋回体3の方位)の値とを加えて回転させる。ここでは地理座標系におけるGNSSアンテナ50Aの中心の座標を(X0,Y0,Z0)とし,同様にフロント作業装置6上の任意の点Pnの座標を(Xn,Yn,Zn)とおく。ただしnは自然数であり,その最大値はフロント作業装置6の3次元モデルを規定する頂点の数である。このときGNSSアンテナ50Aの中心からフロント作業装置6上の任意の点PnまでのベクトルVは下記の式1で表すことが出来る。ただし,X’n=Xn-X0, Y’n=Yn-Y0, Z’n=Zn-Z0とする。ベクトルVを示す座標(X’n,Y’n,Z’n)は第1アンテナ座標系上の座標である。
Figure JPOXMLDOC01-appb-M000001
 このベクトルVを第1アンテナ座標系におけるXY平面上に投影したベクトルVxyは下記の式2で表すことが出来る。
Figure JPOXMLDOC01-appb-M000002
 XY平面上で北方向からベクトルVxyまでの角を方位角θan(図3B参照)と表すと,その角度は第1アンテナ座標系上の座標値X’n,Y’nとヘディング(ヘディング角)とを利用して下記の式3で表すことができる。なお,下記式3では,X’n,Y’nの組み合わせに応じて方位角θanの演算式を5つの場合に分けている。
Figure JPOXMLDOC01-appb-M000003
 ベクトルVxyとベクトルVのなす角を仰角θen(図3C参照)と表すと,その角度は第1アンテナ座標系上の座標値X’n,Y’n,Z’nを利用して下記の式4の通り計算できる。なお,下記式4では,X’n,Y’n,Z’nの組み合わせに応じて仰角θenの演算式を5つの場合に分けている。
Figure JPOXMLDOC01-appb-M000004
 上記の式3,4より,方位角θan及び仰角θenはフロント作業装置6上の点Pnの第1アンテナ座標系における座標値から演算できる。フロント作業装置6(ブーム6A,アーム6B,バケット6C)の寸法は既知である(例えば,車体形状寸法データ(詳細は後述)として予め記憶装置56に格納されている)ため,姿勢センサ75A-75Cによってフロント作業装置6の姿勢(ブーム角,アーム角,バケット角)を特定すれば,その姿勢におけるフロント作業装置6上の任意の点Pnに係る方位角θan及び仰角θenを演算できる。
 本実施形態では,マスク範囲を方位角θanの範囲と仰角θenの範囲の組み合わせによって定義した。つまり,本実施形態のマスク範囲は2つの方位角θanと2つの仰角θenという4つのパラメータによって規定されている。
 図5はGNSSアンテナ50Aに設定される第1マスク範囲21A及び第2マスク範囲22Aの一例をスカイプロット(衛星配置図)上に示した図である。図5ではGNSSアンテナ50Aを基準とする天空を方位角及び仰角の2つをパラメータとする二次元座標で表しており,マスク範囲にドットを付してグレーで表している。円の中心がGNSSアンテナ50Aの中心を示し,当該円の周方向が方位角を,当該円の径方向が仰角を示している。図中のアルファベットG,Rと2桁の数字を内包する複数の円は,それぞれ,GNSS受信機51によって補足された測位衛星の位置を示し,アルファベットG,Rと2桁の数字は各測位衛星の番号(衛星番号)である。
 図5Aに第1マスク範囲21Aの一例を示す。この図の第1マスク範囲21Aは,方位角がα1以上かつα2以下の範囲(但しα2>α1)であって,仰角がα3以上かつα4以下の範囲(但しα4>α3)に設定されている。α1,α2,α3,α4はフロント作業装置6の最大可動領域に基づいて決定されている。α1は旋回体3がフロント作業装置6を最大可動領域内で動作させた場合に方位角θenが取り得る値の最小値であり,α2は同条件で方位角θenが取り得る値の最大値である。同様に,α3は同条件で仰角θenが取り得る値の最小値であり,α4は同条件で仰角θenが取り得る値の最大値である。例えば,図5Aの場合,α1=130,α2=210,α3=0,α4=75である。第1マスク範囲21A及びそれを規定する4つのパラメータは,油圧ショベル1(旋回体3)が旋回動作や走行動作をしない限り固定される。
 図5Bに第2マスク範囲22Aの一例を示す。この図の第2マスク範囲22Aは,方位角がβ1以上かつβ2以下の範囲(但しβ2>β1)であって,仰角がβ3以上かつβ4以下の範囲(但しβ4>β3)に設定されている。β1,β2,β3,β4は,フロント作業装置6と旋回体3の姿勢(ブーム角,アーム角,バケット角,ピッチ角,ロール角,ヘディング)を演算した時刻tにおけるフロント作業装置6の姿勢がGNSSアンテナ50Aの上空を覆う範囲に基づいて決定されている。
 具体的には,第2マスク範囲22Aの演算に際してマスク範囲演算部45は,上記の式(3),(4)を利用して,フロント作業装置6と旋回体3の姿勢(ブーム角,アーム角,バケット角,ピッチ角,ロール角,ヘディング)を演算した或る時刻tにおけるフロント作業装置6の3次元モデル上の全ての頂点Pnについて方位角θanと仰角θenの組み合わせを求める。ここでは,時刻tにおける方位角θanと仰角θenの組み合わせを,(θa1(t),θe1(t)),(θa2(t),θe2(t)),(θa3(t),θe3(t)),…(θan(t),θen(t))と表す。マスク範囲演算部45は,求めた組み合わせの中から最も大きい方位角θa_max(t)と最も大きい仰角θe_max(t)と最も小さい方位角θa_min(t)と最も小さい仰角θe_min(t)を選択する。第1マスク範囲21Aと同様,第2マスク範囲22Aも4つのパラメータ(β1,β2,β3,β4)によって規定されており,β1=θa_min(t),β2=θa_max(t),β3=θe_min(t),β4=θe_max(t)とする。例えば,図5Bの場合,β1=147,α2=185,α3=0,α4=38である。第2マスク範囲22A及びそれを規定する4つのパラメータは,例えば,フロント作業装置6の姿勢が変更される都度演算される。第2マスク範囲22Aは,その性質上,第1マスク範囲21Aよりも小さくかつ第1マスク範囲21Aに内包される(図5B参照)。
 なお,図5Aの場合にはヘディングとα1を同じ値で表現できたが,図5Bの場合には,ヘディングとβ1が異なる。よって,必ずしもヘディングがθa_min(t)になるとは限らない。
 上記のように生成されるマスク範囲21A,22Aは旋回体3に搭載されたGNSSアンテナ50Aの位置を基準に生成されるため,旋回体3を走行動作させたり旋回動作させたりすることでGNSSアンテナ50Aの地理座標系上における位置を変化させてもGNSSアンテナ50Aと共に動く。つまり,地理座標系を基準にしてマスク範囲21A,22Aが固定される訳ではなく,図3Aに示した第1アンテナ座標系に固定されることになる。
 なお,本実施形態では4つのパラメータによってマスク範囲を規定するが,例えば,フロント作業装置6Aを最大可動領域内で動作させた場合にGNSSアンテナ50Aの上空を覆うフロント作業装置6Aの輪郭を第1マスク範囲21の輪郭に設定しても良い。
 (精度演算部46)
 精度演算部46は,マスク範囲演算部45で演算された複数のマスク範囲(例えば,第1マスク範囲21A,第2マスク範囲22A)のそれぞれを利用して測位衛星を選択した場合と,当該複数のマスク範囲のいずれも利用することなく測位衛星を選択した場合(マスク範囲による衛星選択を中断した場合)とを含む複数の場合におけるGNSS受信機51の測位精度をそれぞれ演算(推定)する部分である。ここにおける「複数の場合」の数Nは,GNSSアンテナ50A,50Bごとに存在し,各GNSSアンテナ50A,50Bに設定されたマスク範囲の総数のそれぞれに1(つまり,マスク範囲を使用しない場合)を加えた数値が該当する。すなわち,GNSSアンテナ50Aに2つマスク範囲が設定されている場合にはN=3となる。なお,GNSSアンテナ50BのNの数値はGNSSアンテナ50AのNの数値と異なっていても良い。Nをマスク範囲の選択パターンの数と称することもある。
 GNSS受信機51の測位精度とは,GNSS受信機51による各GNSSアンテナ50A,50Bの測位結果の精度であり,以下では単に「測位精度」と称することもある。この測位精度は,各マスク範囲を利用した後に選択可能な測位衛星に関連するデータ(測位に利用可能な測位衛星の数や配置(例えば,PDOP(Position Dilution of Precision))と,測位に利用可能な測位衛星から受信した衛星信号のSN比やレベルなど(以下,これらを「衛星関連データ」と称することがある))から演算できる。測位精度は上記した衛星関連データのうち少なくとも1つに基づいて演算できるが,測位精度の演算に利用する衛星関連データはN個の場合全てで一致させることが好ましい。なお,上記のようにN=3の場合,精度演算部46によって演算される測位精度の個数も3となる。
 一般的に,衛星数が多いほど,PDOPの値が小さいほど,衛星信号のSN比やレベルが高いほど,測位精度は良くなる。本実施形態では,このことを勘案して衛星関連データから測位精度を示す指標値(測位精度指標値)を演算し,その指標値の高低により精度を評価するものとする。以下では測位精度指標値が高いほど精度が良いものとして説明する。なお,衛星関連データは,GNSS受信機51から受信しても良いし,GNSSアンテナ50A,50Bで受信した衛星信号に基づいてコントローラ40で演算しても良い。GNSS受信機51から受信する場合には,衛星関連データの演算に利用する際に除外する衛星のリスト(後述の除外衛星リスト)をコントローラ40から送信しても良い。
 精度演算部46には,GNSS受信機51が捕捉している衛星の情報が衛星位置抽出部43から入力される。精度演算部46は,各GNSSアンテナ50に設定された複数のマスク範囲のうちの或るマスク範囲(ここでいう「マスク範囲」にはマスクを使用しない場合も含まれるものとする)を利用した後に選択可能な衛星の衛星関連データを,衛星位置抽出部43から入力される複数の衛星の中から当該或るマスク範囲に含まれる衛星を除いた残りの衛星の情報から取得している。
 この点を具体例で説明する。図5の例では,精度演算部46は,第1マスク範囲21Aを利用した場合,第2マスク範囲22Aを利用した場合,いずれのマスク範囲も利用しない場合の測位精度を演算する。
 すなわち,第1マスク範囲21Aを利用して測位衛星を選択する場合(第1の場合)の測位精度(第1測位精度)は,GNSS受信機51が補足したスカイプロット上の全衛星から第1マスク範囲21Aに含まれる衛星(G17,G19,R12,R22)を除外した衛星の衛星関連データに基づいて精度演算部46によって演算される。第2マスク範囲22Aを利用して測位衛星を選択する場合(第2の場合)の測位精度(第2測位精度)は,GNSS受信機51が補足したスカイプロット上の全衛星から第2マスク範囲22Aに含まれる衛星(G17)を除外した衛星の衛星関連データに基づいて精度演算部46によって演算される。いずれのマスク範囲も利用せずに測位衛星を選択する場合(第3の場合)の測位精度(第3測位精度)は,GNSS受信機51が補足したスカイプロット上の全衛星の衛星関連データに基づいて精度演算部46によって演算される。
 精度演算部46で演算されたN個の場合の測位精度はマスク選択部47に出力される。
 (マスク選択部47)
 マスク選択部47は,精度演算部46で演算されたN個の場合の測位精度の中で,最も測位精度が良い場合を選択する。最も測位精度が良い場合には,各GNSSアンテナ50A,50Bに設定された複数のマスク範囲のうち1つが選択されるだけでなく,マスク範囲を利用しない場合が選択されることもある。マスク選択部47で選択されたマスク範囲の情報は除外衛星決定部44に出力される。
 (衛星位置抽出部43)
 衛星位置抽出部43は,GNSS受信機51が衛星信号を捕捉している複数の衛星の位置(地理座標系における仰角及び方位角)を抽出し,除外衛星決定部44と精度演算部46に出力する。
 (除外衛星決定部44)
 除外衛星決定部44は,マスク選択部47で選択されたマスク範囲と,衛星位置抽出部43から入力するGNSS受信機51が捕捉している測位衛星の位置(仰角・方位角)とに基づいて,GNSS受信機51が測位演算に利用しない除外衛星を決定する。具体的には,除外衛星決定部44は,衛星位置抽出部43で抽出された複数の衛星の中からマスク選択部47で選択されたマスク範囲に位置する衛星を除外衛星として決定し,その除外衛星のリストをGNSS受信機51に出力する。
 (GNSS受信機51)
 GNSS受信機51は,除外衛星決定部44から出力される除外衛星のリストを取得し,その時点で衛星信号を捕捉可能な複数の測位衛星の中から当該リストに含まれる衛星を除外した衛星の衛星信号に基づいて,地理座標系における少なくとも1つのGNSSアンテナ50の位置と,2つのGNSSアンテナ50A,50B間の方位(上部旋回体3の方位)とを演算する。これによりGNSS受信機51は,精度演算部46が演算したN個の測位精度の中で最も測位精度が良い場合に選択される測位衛星の衛星信号に基づいて地理座標系における少なくとも1つのGNSSアンテナ50の位置と,2つのGNSSアンテナ50A,50B間の方位を演算することとなる。演算された位置と方位は測位結果入力部42に入力される。
 なお,GNSS受信機51から出力される位置計測結果は,上記の地理座標系の他にも,平面直角座標系,地心直交座標系,及び現場座標系の少なくとも1つ以上の座標系の座標値を出力可能にしても良い。
 (作業装置位置・姿勢演算部41)
 作業装置位置・姿勢演算部41は,測位結果入力部42から入力するGNSSアンテナ50A,50Bの位置及び上部旋回体3の方位と,複数の姿勢センサ(75A,75B,75C,23)の出力から演算される各フロント部材6A,6B,6Cの角度値及び上部旋回体3の傾斜角(ピッチ角及びロール角)と,記憶装置56に記憶された車体形状寸法データとに基づいて,作業装置6の位置及び姿勢(例えば,現場座標系におけるバケット6Cの先端位置及び姿勢)を演算する。作業装置6の位置及び姿勢の演算に利用される車体形状寸法データとしては,例えば,ブーム6Aの両端に位置する2つのピン間の長さ(ブームピン間長LB(図4参照))と,アーム6Bの両端に位置する2つのピン間の長さ(アームピン間長)と,バケット6Cの先端とバケットピン間の長さ(バケット先端長)と,上部旋回体3における2つのGNSSアンテナ50A,50Bの位置関係を車体座標系における方位角で規定したGNSS取付オフセット方位角と,上部旋回体3における2つのGNSSアンテナ50A,50Bの車体座標系おける取付位置(アンテナ取付位置)と,2つのGNSSアンテナ50A,50Bのそれぞれを基準としたフロント作業装置6の最大可動領域がある。
 モニタ60は,作業装置位置・姿勢演算部41で演算された現場座標系における作業装置6の位置及び姿勢データと,記憶装置56に記憶された現場座標系における施工目標面の位置データとに基づいて演算される作業装置6と施工目標面の位置関係を表示できる。この表示によりオペレータは施工目標面に対する作業装置6の位置・姿勢を容易に把握することができる。
 なお,コントローラ40により,作業装置位置・姿勢演算部41で演算された現場座標系における作業装置6の位置及び姿勢データと,記憶装置56に記憶された現場座標系における施工目標面の位置データとに基づいて,バケット先端が施工目標面を超えないように作業装置6の動作(すなわちフロント部材6A,6B,6Cを駆動するアクチュエータ11A,11B,11Cの動作)に制限をかけるマシンコントロールを実行しても良い。
 また,上記のマスク範囲演算部45はマスク範囲を方位角と仰角で規定したが,方位角だけでマスク範囲を規定しても良い。
 次に,上述したマスク範囲について図5を用いて説明を追加する。
 図5Aが,フロント作業装置6が障害物としてGNSSアンテナ50Aの上空視界を最も遮った場合のスカイプロットを表しているとする。この場合,第1マスク範囲21Aと第2マスク範囲22Aが同じ範囲になり,第1マスク範囲21Aにまとめて表すことができる。
 しかし,図5Bの場合にはブーム6Aが図5Aの場合よりも下がっており,フロント作業装置6が障害物としてGNSSアンテナ50Aの上空視界を遮蔽する範囲が図5Aの場合よりも狭くなっている。この場合,第1マスク範囲21Aと第2マスク範囲22Aが異なるため,それぞれのマスク範囲21A,22Aを利用した場合に選択される衛星数も変わる。図5Bの場合に第1マスク範囲21Aを測位に適用すると,当該第1マスク範囲21Aによって衛星R12,G19,R22,G17の合計4つの衛星が測位時に使える衛星から除外される。しかし,図5Bの場合に第2マスク範囲22Aを測位に適用すると,当該第2マスク範囲22Aによって衛星G17のみが測位時に使える衛星から除外される。そのため,第1マスク範囲21Aを利用した場合よりも測位に使える衛星数が3つ増えて測位精度を上げられる。
 このように精度演算部46は各マスク範囲を適用した場合(但し,いずれのマスク範囲も利用しない場合を含めても良い)に測位に使える衛星を抽出し,抽出後の衛星の衛星関連データに基づいて各マスク範囲を利用した場合の測位精度を演算する。マスク選択部47は,精度演算部46が演算した複数の測位精度を比較し,最も良い測位精度を出すマスク範囲を決定し,そのマスク範囲を利用してGNSS受信機51によってGNSSアンテナ50の測位が行われる。以上の一連の流れを図6のフローチャートで説明する。
 <測位処理のフローチャート>
 図6は本実施形態に係るコントローラ40及びGNSS受信機51によるGNSSアンテナ50Aの測位処理のフローチャートである。この処理フローは,一定周期間隔(例えば100ms)で繰返し演算される。なお,ここではGNSSアンテナ50Aの測位についてのみ説明するが,GNSSアンテナ50Bについても同様の処理が行われ,両方のGNSSアンテナ50A,50Bの測位結果に基づいて2つのGNSSアンテナ50A,50B間の方位の演算が行われるものとする。
 ステップ51において,コントローラ40(衛星位置抽出部43)は,GNSS受信機51が衛星信号を捕捉している複数の衛星のデータ(位置データなどを含む)をGNSS受信機51から取得する。
 ステップ52において,コントローラ40(マスク範囲演算部45)は,旋回体姿勢センサ23の検出信号に基づいて旋回体3の姿勢(ピッチ角,ロール角)を演算する。
 ステップ53において,コントローラ40(マスク範囲演算部45)は,3つの姿勢センサ75A,75B,75Cの検出信号に基づいてフロント作業装置6の姿勢を演算する。
 ステップ54において,コントローラ40(マスク範囲演算部45)は,ステップ52で演算した旋回体3の姿勢と,記憶装置56に記憶されたフロント作業装置6の最大可動領域と,前回の周期で演算され記憶装置56に記憶されたGNSSアンテナ50Aの位置と,前回の周期で演算され記憶装置56に記憶された旋回体3の方位(この方位(ヘディング)は,前回の周期で演算された2つのGNSSアンテナ50A,50B間の方位に基づいて演算されて記憶装置56に記憶されている)とに基づいて,第1マスク範囲21Aを演算する。また,コントローラ40(マスク範囲演算部45)は,ステップ52で演算した旋回体3の姿勢と,ステップ53で演算したフロント作業装置6の姿勢と,前回の周期に演算したGNSSアンテナ50Aの位置と,前回の周期に演算した旋回体3の方位(ヘディング)とに基づいて,第1マスク範囲21Aに包含される第2マスク範囲22Aを演算する。
 ステップ55において,コントローラ40(精度演算部46)は,第1マスク範囲21Aに位置する衛星を除外して残りの衛星を選択する第1の場合に補足される衛星の衛星関連データ(衛星関連データは,例えば,衛星数,PDOP(衛星配置),衛星信号のSN比の少なくとも1つであり,衛星位置抽出部43を介してGNSS受信機51から取得できる)に基づいて第1の場合における測位精度である第1測位精度を演算する。同様に,第2マスク範囲22Aに位置する衛星を除外して残りの衛星を選択する第2の場合に補足される衛星の衛星関連データ(但し,第1測位精度を演算した際に利用したものと同じ種類とする)に基づいて第2の場合における測位精度である第2測位精度を演算し,いずれのマスク範囲も利用せずに衛星を選択する第3の場合における衛星関連データ(但し,第1測位精度を演算した際に利用したものと同じ種類とする)に基づいて第3の場合における測位精度である第3測位精度を演算する。
 ステップ56において,コントローラ40(マスク選択部47)は,ステップ55で演算した第1測位精度と第2測位精度を比較する。第1測位精度が第2測位精度よりも測位精度が良い場合にはステップ57に進み,第1測位精度が第2測位精度と測位精度が同じまたは悪い場合にはステップ59に進む。
 ステップ57において,コントローラ40(マスク選択部47)は,ステップ51で位置が取得された複数の衛星の中から第1マスク範囲21Aに位置する衛星を除外衛星として決定し,その除外衛星のリスト(除外衛星リスト)をGNSS受信機51に出力する。
 ここで,第1測位精度が第2測位精度よりも測位精度が良い場合,すなわち,第2マスク範囲22Aよりも大きい第1マスク範囲21Aを利用した方が測位精度が良くなる場合の具体例について図7を用いて説明する。図7はブーム6Aの角度(ブーム角度)とPDOPの時間変化の一例を示した図である。この図において短時間にブーム角度が上昇するタイミングは2回あるが,このタイミングではフロント作業装置6におけるブーム上げ操作と旋回体3の旋回動作が複合され,ダンプトラックなどへの放土作業が行われる。一方,ブーム角度が50度未満に保持される期間が3回あるが,このタイミングではフロント作業装置6による掘削作業が行われる。すなわち,図7は掘削と放土(旋回動作)が繰り返し行われている状況を示している。この図のように掘削と旋回動作を繰り返す場合に第2マスク範囲22Aの利用を継続すると,ブーム6Aの姿勢変化と旋回角度(旋回体3の方位)の変化によって遮蔽される衛星が大きく変化して図に示すようにブーム上げのタイミングの前後でPDOPが大きくなる,すなわち測位精度が低下することがある。そのため,このような場合には,より範囲の広い第1マスク範囲21Aを利用してマスク範囲の変動を抑制した方が第2マスク範囲22Aを利用するよりも測位精度が向上する可能性がある。
 ステップ58において,GNSS受信機51は,ステップ51でコントローラ40にデータを送信した衛星からステップ57の除外衛星リストに含まれる衛星を除外し,その残りの衛星の衛星信号に基づいてGNSSアンテナ50Aの位置を演算し,それをGNSSアンテナ50Aの測位結果としてコントローラ40に出力する。
 ステップ59において,コントローラ40(マスク選択部47)は,ステップ55で演算した第2測位精度と第3測位精度を比較する。第2測位精度が第3測位精度よりも測位精度が良い場合にはステップ60に進み,第2測位精度が第3測位精度と測位精度が同じ又は悪い場合にはステップ62に進む。
 ステップ60において,コントローラ40(マスク選択部47)は,ステップ51で位置が取得された複数の衛星の中から第2マスク範囲22Aに位置する衛星を除外衛星として決定し,その除外衛星のリスト(除外衛星リスト)をGNSS受信機51に出力する。
 ステップ61において,GNSS受信機51は,ステップ51でコントローラ40にデータを送信した衛星からステップ60の除外衛星リストに含まれる衛星を除外し,その残りの衛星の衛星信号に基づいてGNSSアンテナ50Aの位置を演算し,それをGNSSアンテナ50Aの測位結果としてコントローラ40に出力する。
 ステップ62において,GNSS受信機51は,ステップ51でコントローラ40にデータを送信した衛星の衛星信号に基づいてGNSSアンテナ50Aの位置を演算し,それをGNSSアンテナ50Aの測位結果としてコントローラ40に出力する。
 ここで,第3測位精度が第2測位精度よりも測位精度が良い場合,すなわち,いずれのマスク範囲も使用しない方が第2マスク範囲22Aを使用するよりも測位精度が良くなる場合(ステップ62に進む場合)の具体例について図8を用いて説明する。図8は或る時刻の衛星配置と第2マスク範囲22Aを示す図である。この図の場合,フロント作業装置6の方位に存在する衛星は1つ(G28)のみであり,フロントの方位の反対側の方位(180度加えた方位)に存在する衛星数も少ない。このような場合,第2マスク範囲22Aで衛星G28を除外してしまうと衛星数が充分でなく測位精度が低下し易い。つまり,第2マスク範囲22Aを利用すると却って測位精度が低下するため,いずれのマスク範囲も使用しない方が測位精度が良くなる。
 上記のように構成された油圧ショベル1では,複数のマスク範囲が設定されている場合に,各マスク範囲を利用した場合の衛星関連データに基づいてコントローラ40が各マスク範囲を利用した場合の測位精度を演算し,当該演算した測位精度に基づいて当該複数のマスク範囲の中からGNSSアンテナ50A,50Bの測位に利用するマスク範囲を1つ選択する。これによりフロント作業装置6の姿勢に基づいてマスク範囲を選択する場合よりも測位精度が向上し得る。
 図6のフローチャートでは,コントローラ40は,大きさの異なる複数のマスク範囲を利用した場合(マスク範囲を設定しない場合も含む)におけるそれぞれの衛星関連データに基づいて各場合の測位精度を演算し,その中から最も測位精度が良い場合に係るマスク範囲を1つ選択する。つまり,マスク範囲の選択に衛星関連データが考慮されることにより,フロント作業装置6の姿勢とともにマスク範囲(上記実施形態では第2マスク範囲22A)を変化させる場合よりも測位精度の低下を抑制できるので,結果的に油圧ショベル1による作業精度を向上させることができる。
 特に本実施形態では,マスク範囲の候補として,フロント作業装置6の可動領域の全てをカバーする第1マスク範囲21Aと,そのときのフロント作業装置6の姿勢で隠れる領域をカバーする第2マスク範囲22Aとが含まれる。これによりマスク範囲の選択及び測位精度の演算にフロント作業装置6の姿勢が考慮されることになり,フロント作業装置6が測位信号を遮蔽することが測位精度に与える影響も小さくできる。
 <その他>
 図6のフローにおいて,ステップ56で第1測位精度が第2測位精度より精度が良いと判定された場合には,さらに第1測位精度と第3測位精度を比較する判定処理を追加しても良い。この判定処理において,第1測位精度の方が精度が良いと判定された場合にはステップ57に進み,第3測位精度の方が精度が良いと判定された場合にはステップ62と同じ処理(すなわち,いずれのマスク範囲を使うことなく測位する処理)を行うようにしても良い。
 本実施形態では2つのマスク範囲21A,22Aを利用する場合について説明したが,各マスク範囲が異なっていれば,利用可能なマスク範囲の数に限定はない。例えば,第1マスク範囲より大きいものを設定しても良い。また,油圧ショベルで繰り返し行われる作業(例えば,掘削作業,運搬作業,放土作業)のうち,同一の作業におけるフロントの動作範囲を考慮してマスク範囲を設定し,油圧ショベルで行われる作業を推定して当該作業に設定されたマスク範囲を利用するようにしても良い。後者の場合のマスク範囲は,第1マスク範囲21Aより小さく,第2マスク範囲22Aより大きくなることが想定される。
 なお,本発明は,上記の実施の形態に限定されるものではなく,その要旨を逸脱しない範囲内の様々な変形例が含まれる。例えば,本発明は,上記の実施の形態で説明した全ての構成を備えるものに限定されず,その構成の一部を削除したものも含まれる。また,ある実施の形態に係る構成の一部を,他の実施の形態に係る構成に追加又は置換することが可能である。
 また,上記のコントローラ40に係る各構成や当該各構成の機能及び実行処理等は,それらの一部又は全部をハードウェア(例えば各機能を実行するロジックを集積回路で設計する等)で実現しても良い。また,上記のコントローラ40に係る構成は,演算処理装置(例えばCPU)によって読み出し・実行されることで当該コントローラ40の構成に係る各機能が実現されるプログラム(ソフトウェア)としてもよい。当該プログラムに係る情報は,例えば,半導体メモリ(フラッシュメモリ,SSD等),磁気記憶装置(ハードディスクドライブ等)及び記録媒体(磁気ディスク,光ディスク等)等に記憶することができる。
 また,上記の各実施の形態の説明では,制御線や情報線は,当該実施の形態の説明に必要であると解されるものを示したが,必ずしも製品に係る全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えて良い。
 1…油圧ショベル,2…走行体(下部走行体),3…旋回体(上部旋回体),4…運転席,6…フロント作業装置,6A…ブーム,6B…アーム,6C…バケット,7…無線機,8…GNSS基準局,8…基準局,11A…ブームシリンダ,11B…アームシリンダ,11C…バケットシリンダ,21A…第1マスク範囲,22A…第2マスク範囲,23…旋回体姿勢センサ,40…コントローラ,41…作業装置位置・姿勢演算部,42…測位結果入力部,43…衛星位置抽出部,44…除外衛星決定部,45…マスク範囲演算部,46…精度演算部,47…マスク選択部,50A…GNSSアンテナ(測位アンテナ),50B…GNSSアンテナ(測位アンテナ),51…GNSS受信機,56…記憶装置,60…モニタ,75A…ブーム姿勢センサ,75B…アーム姿勢センサ,75C…バケット姿勢センサ,80…基準局GNSSアンテナ,81…GNSS受信機,82…基準局コントローラ,87…無線機

Claims (9)

  1.  下部走行体と,
     前記下部走行体の上に旋回可能に取り付けられた上部旋回体と,
     前記上部旋回体に取り付けられた多関節型の作業装置と,
     前記作業装置及び前記上部旋回体の姿勢を検出するための複数の姿勢センサと,
     前記上部旋回体に取り付けられ,複数の測位衛星からの衛星信号を受信するための測位アンテナと,
     前記測位アンテナで受信される衛星信号に基づいて前記測位アンテナの位置を演算する受信機と,
     前記複数の姿勢センサの検出信号に基づいて,前記作業装置の姿勢と前記上部旋回体の姿勢とを演算するコントローラとを備えた作業機械において,
     前記コントローラには,前記受信機による前記測位アンテナの位置の演算に利用する測位衛星を選択するために,前記作業装置及び前記上部旋回体の姿勢に基づいて前記測位アンテナを基準に設定される複数のマスク範囲が記憶されており,
     前記コントローラは,前記複数のマスク範囲のそれぞれを利用して測位衛星を選択した場合のそれぞれの測位精度と,前記複数のマスク範囲のいずれも利用することなく測位衛星を選択した場合の測位精度とを含む複数の測位精度を,選択した前記測位衛星の衛星関連データに基づいて演算し,
     前記受信機は,前記複数の測位精度のうち最も測位精度が良い場合に選択した前記測位衛星の衛星信号に基づいて前記測位アンテナの位置を演算する
     ことを特徴とする作業機械。
  2.  請求項1の作業機械において,
     前記複数のマスク範囲には,
      前記測位アンテナが衛星信号を受信する際に前記作業装置が障害物となり得る範囲に設定された第1マスク範囲と,
      前記第1マスク範囲に包含される範囲に設定された第2マスク範囲と
     が含まれていることを特徴とする作業機械。
  3.  請求項2の作業機械において,
     前記第1マスク範囲は,演算した前記上部旋回体の姿勢と,前記作業装置の最大可動領域とに基づいて,前記測位アンテナが衛星信号を受信する際に前記作業装置が障害物となり得る最大範囲に設定されており,
     前記第2マスク範囲は,演算した前記上部旋回体の姿勢と,演算した前記作業装置の姿勢とに基づいて,前記作業装置の姿勢演算時における前記作業装置が,前記測位アンテナが前記複数の測位衛星から衛星信号を受信する際に障害物となる範囲に設定されている
     ことを特徴とする作業機械。
  4.  請求項2または3の作業機械において,
     前記コントローラは,
      前記第1マスク範囲を利用した第1の場合における前記衛星関連データに基づいて前記第1の場合における測位精度である第1測位精度を演算し,
      前記第2マスク範囲を利用した第2の場合における前記衛星関連データに基づいて前記第2の場合における測位精度である第2測位精度を演算し,
      いずれのマスク範囲も利用しない第3の場合における前記衛星関連データに基づいて前記第3の場合における測位精度である第3測位精度を演算し,
     前記受信機は,前記第1測位精度,前記第2測位精度および前記第3測位精度の中で最も測位精度が良い場合に選択される測位衛星の衛星信号に基づいて前記測位アンテナの位置を演算する
     ことを特徴とする作業機械。
  5.  請求項2または3の作業機械において,
     前記コントローラは,
      前記第1マスク範囲を利用した第1の場合における前記衛星関連データに基づいて前記第1の場合における測位精度である第1測位精度を演算し,
      前記第2マスク範囲を利用した第2の場合における前記衛星関連データに基づいて前記第2の場合における測位精度である第2測位精度を演算し,
      いずれのマスク範囲も利用しない第3の場合における前記衛星関連データに基づいて前記第3の場合における測位精度である第3測位精度を演算し,
     前記受信機は,前記第1測位精度が前記第2測位精度よりも測位精度が良い場合には,前記受信機は,前記第1マスク範囲を利用した場合に選択される測位衛星の衛星信号に基づいて前記測位アンテナの位置を演算する
     ことを特徴とする作業機械。
  6.  請求項5の作業機械において,
     前記第1測位精度が前記第2測位精度と測位精度が同じ又は悪い場合であって,前記第2測位精度が前記第3測位精度よりも測位精度が良い場合には,前記受信機は,前記第2マスク範囲を利用した場合に選択される測位衛星の衛星信号に基づいて前記測位アンテナの位置を演算する
     ことを特徴とする作業機械。
  7.  請求項6の作業機械において,
     前記第1測位精度が前記第2測位精度と測位精度が同じ又は悪い場合であって,前記第2測位精度が前記第3測位精度と測位精度が同じ又は悪い場合には,前記受信機は,いずれの衛星マスクも利用することなく選択される測位衛星の衛星信号に基づいて前記測位アンテナの位置を演算する
     ことを特徴とする作業機械。
  8.  請求項1の作業機械において,
     前記衛星関連データには,測位に利用可能な測位衛星の数と,測位に利用可能な測位衛星の配置と,測位に利用可能な測位衛星の衛星信号のSN比とのうち少なくとも1つが含まれることを特徴とする作業機械。
  9.  請求項1の作業機械において,
     前記測位アンテナは第1測位アンテナであり,
     前記上部旋回体に取り付けられ,前記複数の測位衛星からの衛星信号を受信するための第2測位アンテナをさらに備え,
     前記受信機は,前記第1測位アンテナ及び前記第2測位アンテナの2つの測位アンテナで受信される衛星信号に基づいて,前記2つの測位アンテナ間の方位を演算し,
     前記コントローラは,前記第1測位アンテナの位置と,前記2つの測位アンテナ間の方位とに基づいて,前記上部旋回体の方位を演算し,
     前記複数のマスク範囲は,前記作業装置及び前記上部旋回体の姿勢と,前記上部旋回体の方位とに基づいて設定されている
     ことを特徴とする作業機械。
PCT/JP2022/004694 2021-03-30 2022-02-07 作業機械 WO2022209303A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US18/024,991 US20230374757A1 (en) 2021-03-30 2022-02-07 Work machine
EP22779532.5A EP4317617A1 (en) 2021-03-30 2022-02-07 Work machine
KR1020237007658A KR20230045068A (ko) 2021-03-30 2022-02-07 작업 기계
CN202280005880.4A CN116057417A (zh) 2021-03-30 2022-02-07 作业机械

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2021058579A JP7039746B1 (ja) 2021-03-30 2021-03-30 作業機械
JP2021-058579 2021-03-30

Publications (1)

Publication Number Publication Date
WO2022209303A1 true WO2022209303A1 (ja) 2022-10-06

Family

ID=81214265

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2022/004694 WO2022209303A1 (ja) 2021-03-30 2022-02-07 作業機械

Country Status (6)

Country Link
US (1) US20230374757A1 (ja)
EP (1) EP4317617A1 (ja)
JP (1) JP7039746B1 (ja)
KR (1) KR20230045068A (ja)
CN (1) CN116057417A (ja)
WO (1) WO2022209303A1 (ja)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004184121A (ja) 2002-11-29 2004-07-02 Denso Corp Gps受信機の制御装置、及び、サーバ装置
JP2020008286A (ja) * 2018-07-02 2020-01-16 日立建機株式会社 作業機械
JP2020139933A (ja) * 2019-03-01 2020-09-03 日立建機株式会社 作業機械
JP2020144014A (ja) * 2019-03-06 2020-09-10 日立建機株式会社 作業機械
JP2020200597A (ja) * 2019-06-06 2020-12-17 日立建機株式会社 建設機械
WO2021060533A1 (ja) * 2019-09-26 2021-04-01 日立建機株式会社 作業機械
JP2021148467A (ja) * 2020-03-16 2021-09-27 日立建機株式会社 作業機械
JP2021195803A (ja) * 2020-06-16 2021-12-27 日立建機株式会社 建設機械

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004184121A (ja) 2002-11-29 2004-07-02 Denso Corp Gps受信機の制御装置、及び、サーバ装置
JP2020008286A (ja) * 2018-07-02 2020-01-16 日立建機株式会社 作業機械
JP2020139933A (ja) * 2019-03-01 2020-09-03 日立建機株式会社 作業機械
JP2020144014A (ja) * 2019-03-06 2020-09-10 日立建機株式会社 作業機械
JP2020200597A (ja) * 2019-06-06 2020-12-17 日立建機株式会社 建設機械
WO2021060533A1 (ja) * 2019-09-26 2021-04-01 日立建機株式会社 作業機械
JP2021148467A (ja) * 2020-03-16 2021-09-27 日立建機株式会社 作業機械
JP2021195803A (ja) * 2020-06-16 2021-12-27 日立建機株式会社 建設機械

Also Published As

Publication number Publication date
JP7039746B1 (ja) 2022-03-22
US20230374757A1 (en) 2023-11-23
CN116057417A (zh) 2023-05-02
EP4317617A1 (en) 2024-02-07
JP2022155192A (ja) 2022-10-13
KR20230045068A (ko) 2023-04-04

Similar Documents

Publication Publication Date Title
US9886038B2 (en) GNSS and optical guidance and machine control
US8594879B2 (en) GNSS guidance and machine control
US8639416B2 (en) GNSS guidance and machine control
US7831362B2 (en) Position measuring system for working machine
US20120059554A1 (en) Automatic Blade Control System during a Period of a Global Navigation Satellite System ...
JP2866289B2 (ja) 建設機械の位置及び姿勢表示方法
JP2007064853A (ja) 複合測位による移動体の測位制御装置、システムおよびプログラム
KR102606049B1 (ko) 건설 기계
US20220220695A1 (en) Work machine
JP7419119B2 (ja) 作業機械
WO2022209303A1 (ja) 作業機械
JP4202209B2 (ja) 作業機械の位置計測表示システム
JP7007313B2 (ja) 作業機械
JP7500420B2 (ja) 作業機械
JPH09196671A (ja) 双方向自動追尾システムによる位置決め方法及び装置
JP2022103603A (ja) 作業機械
JP7102071B2 (ja) 作業機械
JP2024050354A (ja) 作業機械の施工精度予測計算装置
JP2023012798A (ja) 作業機械
JP2023122748A (ja) 作業位置指示システム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22779532

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 20237007658

Country of ref document: KR

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 2022779532

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2022779532

Country of ref document: EP

Effective date: 20231030