The present invention relates to a construction machine, and more particularly to a construction machine that controls, on the basis of position information of the own machine, a work device or the like that changes in position and posture.
In a field of construction machines such as hydraulic excavators, the introduction of intelligent construction, which efficiently utilizes various information by applying an information communication technology to construction work, and thereby streamlines the work, has recently been under way. For example, there are construction machines having functions of assisting in operator’s operation such as machine guidance that displays, to the operator, the position and posture of an articulated work device formed by coupling together a plurality of members such as a boom, an arm, and a bucket, and machine control that performs control such that the work device moves along a target construction surface.
Intelligent construction using the coordinates of an own machine on a construction site when performing work support such as the machine guidance or the machine control is referred to as three-dimensional intelligent construction (hereinafter referred to as 3D intelligent construction). A construction machine supporting the 3D intelligent construction includes a satellite positioning system (GNSS: Global Navigation Satellite System) in order to obtain the position of the own machine. The GNSS receives positioning signals from a plurality of satellites, and measures the three-dimensional position (latitude, longitude, and altitude) of the own machine. A construction machine with a work device, such as a hydraulic excavator, necessitates, for the work support, not only the position of the own machine but also a direction in which the work device is facing (orientation of the work device). Therefore, a construction machine is publicly known which is equipped with two GNSS antennas that receive the positioning signals, and identifies the orientation of the work device on the basis of the positioning signals received by the GNSS antennas.
Because of the use of the GNSS, positioning is not normally performed when positioning signals from a sufficient number of satellites cannot be captured. In this case, it is not possible to obtain posture information in a three-dimensional space of the construction machine, the posture information indicating the posture (the three dimensional position coordinates and orientation or the like) of a machine body and the posture of the work device, so that the work support such as the machine guidance or the machine control needs to be stopped. A high frequency of stopping the work support leads to a significant decrease in work efficiency.
In addition, it is known that a small number of satellites available for the positioning can cause variation in positioning result even when positioning using the GNSS can be performed. In this case, there is a fear that highly accurate position information cannot be obtained, and that the performance of the machine guidance or the machine control is consequently decreased.
For such a problem, a technology described in Patent Document 1 has been proposed. A control system of a work machine described in Patent Document 1 is intended to reduce an effect of variation in the positioning result of a positioning device (corresponding to the GNSS) by smoothing a first position (positioning result of the GNSS) measured by the positioning device, using working information of the work machine detected by a state device. In addition, the control system is intended to make it possible to continue the machine guidance or the machine control, even if an abnormality occurs in positioning by the positioning device, by continuing to use position information obtained by the smoothing processing before the abnormality in the positioning as long as a traveling or a swinging (action) of the machine body is not performed.
Patent Document 1: WO 2015/186845
The work device of a hydraulic excavator as one of construction machines is formed by metallic members. Thus, the positioning signals from the satellites may be reflected or interrupted by the work device. When the hydraulic excavator performs excavating, up and down motion of the work device are repeated, so that the reflection or interruption of the positioning signals may occur repeatedly. Under such an environment, though the number of positioning signals that can be received by the GNSS antennas changes, a sufficient number of satellites available for positioning computation are usually secured. Therefore, the GNSS can continue highly accurate positioning computation (calculation of an RTK-FIX solution).
However, change of the number and arrangement of the satellites available for the positioning (satellite positioning environment) may cause different positioning results of the GNSS. For example, when the reflection or interruption of the positioning signals due to motion of the work device decreases the number of available satellites and causes an imbalance in the satellite arrangement, the positioning results may change even though the positioning results are highly accurate. In order to implement the machine guidance and the machine control of the hydraulic excavator, a strict required accuracy specified in implementation guidelines of intelligent construction (for example ±5 cm in ordinary civil engineering work) needs to be satisfied. However, a change in the positioning result, which is caused by a change in the satellite positioning environment, may adversely affect the implementation of intelligent construction. Therefore, there is a desire for a method of obtaining highly accurate position information even when a change in the posture of the hydraulic excavator itself causes a change in the computation result of satellite positioning.
In the control system described in Patent Document 1, a positioning condition of the positioning device (corresponding to the GNSS) is determined from the reception of the positioning signals on the positioning device, or the like. It is determined that the positioning condition is normal when positioning accuracy is excellent (Fix). On the other hand, it is determined that the positioning condition is abnormal when positioning is not possible or when positioning is possible but the positioning accuracy is poor. When a determination result represents normality, smoothing processing is performed while the work machine is in a non-travelling state and a non-swinging state. When the determination result represents abnormality, on the other hand, the smoothing processing is not performed.
Also in the control system described in Patent Document 1, when motion of the work device causes reflection or interruption of a part of the positioning signals, the accuracy of positioning of the positioning device is assumed to be maintained in an excellent condition (Fix). In this case, it is determined that the positioning condition is normal, and therefore the smoothing processing is performed under conditions where the work machine is not travelling or swinging. Depending on a method of the smoothing processing, variation in the positioning result of the positioning device, which is caused by a change in the satellite positioning environment, may not be able to be suppressed sufficiently. For example, when the variation in the positioning result is large, smoothing processing of a low-pass filter with a small time constant, a moving average, or the like may not be able to smooth the positioning results sufficiently, so that the position information to be used for the machine guidance or the like may vary.
Even if the smoothing processing used by the control system described in Patent Document 1 can suppress variation in the positioning result caused by a change in the satellite positioning environment, it appears to be difficult to obtain accurate position information of the own machine when movement of the work machine occurs without operator’s operation. For example, under an environment where the work machine tends to slip, such as an environment with a ground after rainy weather or the like, an underbody may slip only by work of the work device. In this case, the position of the own machine (GNSS antennas) is shifted. The control system determines a non-travelling state and a non-swinging state of the work machine according to the presence or absence of operation of control levers, and determines whether to perform or discontinue the smoothing processing according to a result of the determination. When the work machine slips, the control levers are not operated, and therefore the smoothing processing is performed on the computation results (positioning results) of the positioning device. Therefore, even though the displacement of the GNSS antennas (own machine) occurs actually, the smoothing processing causes a response delay, and it is difficult to obtain accurate position information of the own machine.
The present invention has been made based on the above-described circumstances. It is an object of the present invention to provide a construction machine that can obtain posture information of the construction machine based on accurate position information of the own machine even when variation occurs in computation results of satellite positioning or even when movement of the own machine occurs without operator’s operation.
The present application includes a plurality of means for solving the above-described problems. To cite an example thereof, there is provided a construction machine including: a machine body; a work device attached to the machine body in a raiseable and lowerable manner; an antenna that is attached to the machine body and receives positioning signals from a plurality of satellites; a first sensor that senses information on a posture and motion of the machine body; a second sensor that senses information on a posture of the work device; and a computing device that computes posture information indicating the postures of the machine body and the work device. The computing device is configured to perform positioning computation that is to compute a position of the machine body and a variance value of the position of the machine body on a basis of the positioning signals from the plurality of satellites received by the antenna, subject the position of the machine body obtained by the positioning computation to first smoothing processing that is to increase a degree of smoothing as magnitude of the variance value of the position of the machine body obtained by the positioning computation becomes larger, and compute the posture information based on a result of the first smoothing processing, the information sensed by the first sensor, and the information sensed by the second sensor.
According to the present invention, even when variation in the position as a result of the positioning computation becomes large due to a degradation in the satellite positioning environment, the variance value as a result of the positioning computation is increased, and therefore, by correspondingly strengthening smoothing on the position as a result of the positioning computation, it is possible to obtain position information in which the variation in the position as a result of the positioning computation is suppressed. In addition, when a movement of the machine body occurs without operator’s operation, the variance value as a result of the positioning computation hardly changes unless the satellite positioning environment is degraded, and therefore position information following the movement of the machine body can be obtained by smoothing the position as a result of the positioning computation with a strength corresponding to the magnitude of the variance value as a result of the positioning computation. That is, even when variation occurs in computation result of satellite positioning or even when the movement of the own machine occurs without operator’s operation, it is possible to obtain the posture information of the construction machine based on accurate position information of the own machine.
Problems, configurations, and effects other than those described above will be made apparent by the following description of embodiments.
Construction machines according to embodiments of the present invention will hereinafter be described with reference to the drawings. In the present embodiments, description will be made by taking a hydraulic excavator as an example of the construction machine.
A general configuration of a hydraulic excavator as a construction machine according to a first embodiment of the present invention will first be described with reference to
In
The front work device 1 is an articulated work device formed by coupling together a plurality of driven members rotatably in a vertical direction. The plurality of driven members include, for example, a boom 6, an arm 7, and a bucket 8 as a work tool. A proximal end portion of the boom 6 is rotatably supported by a front portion of the upper swing structure 3. A proximal end portion of the arm 7 is rotatably supported by a distal end portion of the boom 6. The bucket 8 is rotatably supported by a distal end portion of the arm 7. The boom 6, the arm 7, and the bucket 8 are respectively driven by a boom cylinder 10, an arm cylinder 11, and a bucket cylinder 12 as hydraulic actuators.
The lower track structure 2, for example, has a crawler type track device 14 (only a track device on one side is shown). The track device 14 is driven by a travelling hydraulic motor 14a as a hydraulic actuator.
The upper swing structure 3 is, for example, configured to be swing-driven with respect to the lower track structure 2 by a swing hydraulic motor 4 as a hydraulic actuator. The upper swing structure 3 includes a cab 16 where an operator boards and a machine room 17 housing various kinds of apparatuses.
The cab 16 is provided with operation devices 18a and 18b for operating the hydraulic actuators 4, 10, 11, 12, and 14a. Each of the operation devices 18a and 18b is, for example, an electric control lever device with a control lever that can be inclined forward and rearward and leftward and rightward. The electric control lever devices 18a and 18b have a sensor (not shown) that electrically senses an inclination direction and an inclination amount, that is, an operation direction and an operation amount of the control lever. The electric control lever devices 18a and 18b output operation signals corresponding to the sensed operation directions and the sensed operation amounts to a controller 40 (see
A prime mover 21 such as an engine or a motor, a hydraulic pump device 22 driven by the prime mover 21, and the like are arranged in the machine room 17. A hydraulic fluid delivered from the hydraulic pump device 22 is supplied to each of the hydraulic actuators 4, 10, 11, 12, and 14a, thereby driving each of the hydraulic actuators 4, 10, 11, 12, and 14a. The driving of each of the hydraulic actuators 4, 10, 11, 12, and 14a is controlled by a control valve unit 23 as an aggregate of control valves corresponding to the respective hydraulic actuators 4, 10, 11, 12, and 14a. The control valves constituting the control valve unit 23 control directions and flow rates of the hydraulic fluid supplied from the hydraulic pump device 22 to the corresponding hydraulic actuators 4, 10, 11, 12, and 14a. The driving of each control valve is, for example, controlled by an operation pilot pressure output from a pilot pump (not shown) via a solenoid proportional valve (not shown). Each solenoid proportional valve is controlled by the controller 40 based on the operation signals from the operation devices 18a and 18b. Drives of the respective hydraulic actuators 4, 10, 11, 12, and 14a are thereby controlled via the respective control valves of the control valve unit 23.
As shown in
In the boom 6, the arm 7, and the bucket 8 as constituent members of the front work device 1, inertial measurement units 26, 27, and 28 that measure the angular velocities and accelerations of the respective constituent members are respectively installed. A measuring device constituted by these three inertial measurement units 26, 27, and 28 functions as a sensor that senses information on the posture of the front work device 1. In order to distinguish these three inertial measurement units 26, 27, and 28 from one another, the inertial measurement unit 26 for the boom 6 will be referred to as a boom IMU, the inertial measurement unit 27 for the arm 7 will be referred to as an arm IMU, and the inertial measurement unit 28 for the bucket 8 will be referred to as a bucket IMU. The boom IMU 26, the arm IMU 27, and the bucket IMU 28 each output a sensing result to the controller 40. The measuring device including the three IMUs, that is, the boom IMU 26, the arm IMU 27, and the bucket IMU 28 functions as a sensor that senses information on the posture of the front work device 1. In addition, a measuring device having a configuration obtained by adding the machine body IMU 25 to the three IMUs, that is, the boom IMU 26, the arm IMU 27, and the bucket IMU 28 functions as a sensor that senses information on the posture in a three-dimensional space of the construction machine.
In addition, two GNSS antennas 31 and 32 capable of receiving positioning signals from satellites are attached to the upper swing structure 3. The positioning signals received by the respective GNSS antennas 31 and 32 are input to a GNSS receiver 33 shown in
The GNSS receiver 33 (positioning system 30) can perform RTK (Real Time Kinematic) positioning by connecting to a fixed station of GNSS installed within a site via wireless communication. In a case of a site without a fixed station of GNSS, positioning using network type RTK that obtains the information of an electronic reference station via the Internet can be performed. In the following, the GNSS receiver 33 is assumed to be able to perform RTK positioning irrespective of the presence or absence of the fixed station within the site.
As shown in
The controller 40 has a positioning computing section 51, a construction target surface computing section 52, a monitor display control section 53, and a hydraulic system control section 54 as a part of the functions performed by the processor 42.
The positioning computing section 51 computes posture information in the three-dimensional space of the hydraulic excavator, which indicates the posture (three-dimensional position coordinates and an orientation or the like) of the machine body and the posture of the front work device 1 within the work site, on the basis of the computation results of the GNSS receiver 33 and the sensing results of the machine body IMU 25, the boom IMU 26, the arm IMU 27, and the bucket IMU 28. The positioning computing section 51 outputs a computation result to the construction target surface computing section 52, the monitor display control section 53, and the hydraulic system control section 54. Accuracy of the computation result of the positioning computing section 51 with respect to the actual position and posture of the hydraulic excavator affects the computation or control of these sections. Therefore, the computation result of the positioning computing section 51 is desired to be maintained at high accuracy. Details of the configuration of the positioning computing section 51 will be described later.
The construction target surface computing section 52 computes a construction target surface that defines the target shape of a construction object, based on construction information of a three-dimensional working drawing or the like, which is stored in the storage device 41 in advance, and the posture information of the hydraulic excavator as the computation result of the positioning computing section 51. The construction information is, for example, input by a construction manager via the monitor 19 or the like as an input device. The construction target surface computing section 52 further calculates a distance between the construction target surface and a reference point (for example a claw tip of the bucket 8). Incidentally, the construction target surface can also be set by the operator on the spot by operating the monitor 19. The construction target surface computing section 52 outputs a computation result to the monitor display control section 53 and the hydraulic system control section 54.
The monitor display control section 53 controls the display of the monitor 19 within the cab 16. The monitor display control section 53 computes instruction contents of operation assistance for the operator based on the construction target surface as the computation result of the construction target surface computing section 52 and the posture information of the hydraulic excavator as the computation result of the positioning computing section 51. The monitor display control section 53 displays a computation result on the monitor 19. The monitor display control section 53, for example, performs a part of functions as a machine guidance system that assists in operator’s operation by displaying the posture of the front work device 1 and the tip position and angle of the bucket 8 of the front work device 1 on the monitor 19. When an abnormality occurs in the positioning result of the GNSS receiver 33, an interruption of the guidance functions can be notified to the operator by displaying a warning on the monitor 19. The monitor 19 can be configured to give a work instruction or call attention by not only making screen display but also outputting sound. In addition, the monitor 19 is desirably usable as an input device by including a touch panel, rather than being a mere display device. The monitor 19 can, for example, be configured to be installed within the cab 16, or configured as a portable tablet terminal detachably installed within the cab 16.
The hydraulic system control section 54 controls a hydraulic system including the hydraulic pump device 22, the plurality of hydraulic actuators 4, 10, 11, 12, and 14a, and the control valve unit 23. The hydraulic system control section 54, for example, computes an action of the hydraulic excavator based on the construction target surface as the computation result of the construction target surface computing section 52 and the posture information of the hydraulic excavator as the computation result of the positioning computing section 51, and controls the hydraulic system so as to realize the action as a computation result. Specifically, the hydraulic system control section 54 performs a part of functions as a machine control system that restricts the action such that a tip of the bucket 8 does not approach the construction target surface over a certain distance, or performs control such that the bucket 8 moves along the construction target surface. The hydraulic system control section 54 preferably stops the machine control functions when an abnormality occurs in the positioning result of the GNSS receiver 33 and the monitor display control section 53 displays a warning on the monitor 19.
In the following, a system combining the functions of both of the construction target surface computing section 52 and the monitor display control section 53 will be referred to as a machine guidance system. In addition, a system combining the functions of both of the construction target surface computing section 52 and the hydraulic system control section 54 will be referred to as a machine control system.
In the present embodiment, the controller 40 including the positioning computing section 51 and the GNSS receiver 33 constitute a computing device that computes the posture information of the hydraulic excavator. In addition, while in the present embodiment, one controller 40 is configured to execute all of the functional sections shown in
Next, referring to
The positioning computing section 51 of the controller 40 uses the computation results of the GNSS receiver 33 (positioning system 30) using the positioning signals of the satellites, in order to compute the posture information of the hydraulic excavator. The posture information as the computation result of the positioning computing section 51 is used in the monitor display control section 53 and the hydraulic system control section 54. When positions as computation results of the positioning system 30 vary despite no actual displacement of the GNSS antennas 31 and 32, the variation affects the computation result of the positioning computing section 51. Consequently, there is a fear that the claw tip position of the bucket 8 of the hydraulic excavator displayed on the monitor 19 by the monitor display control section 53 may vary from an actual position, and that it may therefore become difficult to perform appropriate machine guidance. In addition, there is a fear that the claw tip position of the bucket 8 controlled by the hydraulic system control section 54 may become different from the actual position, and that a finished surface may therefore have an undulating shape. The following are cited as factors that vary the computation results of the positioning system 30 despite no actual displacement.
The driven members of the boom 6, the arm 7, and the bucket 8 constituting the front work device 1 are formed by metallic members. Thus, the positioning signals from the satellites may be reflected or interrupted due to a raising motion of the front work device 1. In excavation motion of the hydraulic excavator, up and down motion (elevating motion) of the front work device 1 are repeated. There is thus a fear that the reflection or interruption of the positioning signals may occur repeatedly. Under such conditions, the number of positioning signals that can be received by the GNSS antennas 31 and 32 is increased or decreased. Even under such conditions, a sufficient number of satellites available for the positioning computation of the positioning system 30 can be usually secured, and therefore the positioning system 30 can continue highly accurate positioning computation (calculation of an RTK-FIX solution). However, the change of the number or arrangement of satellites, that is, the change of the satellite positioning environment may result in positioning results with different values (variation in positioning result) despite no actual displacement of the GNSS antennas 31 and 32. This phenomenon will be described in the following.
On the other hand,
When the imbalance in the satellite arrangement is small as in the celestial diagram shown in
However, in order to implement the machine guidance and the machine control of the hydraulic excavator, a required accuracy of implementation guidelines for intelligent construction (±5 cm in ordinary civil engineering work) needs to be satisfied. Thus, even a slight difference in results of the positioning computation cannot be ignored. When the machine guidance or the like is performed without consideration being given to differences between the positioning computation result (the average value and the error ellipse) obtained under the satellite positioning environment shown in
In addition, in order to implement the machine guidance or the like, when the machine body (GNSS antennas) is actually moved (displaced) according to operator’s operation, various kinds of information to be used in the machine guidance or the like are desired to follow the actual displacement. For example, when the hydraulic excavator travels or swings, the various kinds of information are desired to follow a displacement corresponding to the travelling or the swinging. In addition, also when the machine body is slipped due to only a motion of the front work device 1 without the hydraulic excavator travelling or swinging, and the GNSS antennas are thereby displaced without operator’s operation, the various kinds of information to be used in the machine guidance or the like are desired to follow an actual displacement.
In addition, even when the hydraulic excavator is not travelling or swinging, the actual position of the GNSS antennas may be slightly displaced by vibration caused by motion of the front work device 1 or the like. At this time, the positioning result of the positioning system 30 may vary due to the actual displacement caused by the vibration of the hydraulic excavator. The variation in the positioning result of the positioning system 30 caused by the vibration of the hydraulic excavator needs to be suppressed in order to implement the machine guidance or the like.
The various kinds of information to be used for the machine guidance and the machine control are thus desired to be maintained at high accuracy with respect to the actual information. Accordingly, the controller 40 according to the present embodiment is configured to correct the positioning results of the positioning system 30 according to a change in the satellite positioning environment and the presence or absence of motion of the machine body (upper swing structure 3).
Next, referring to
In
The position correction computing section 61 corrects the position as a computation result of the positioning system 30 according to the satellite positioning environment (the number and arrangement of satellites or the like from which reception can be performed by the GNSS antennas 31 and 32) and the presence or absence of motion of the machine body (travelling, swinging, or the like) to obtain the corrected position of a predetermined part of the machine body (GNSS antennas 31 and 32). The position correction computing section 61 is characterized by subjecting the position as a computation result of the positioning system 30 to two stages of smoothing processing, that is, a smoothing processing according to the satellite positioning environment and a smoothing processing according to the presence or absence of machine body motion. Specifically, the position correction computing section 61, for example, includes: a motion determining section 64 that determines the presence or absence of motion (displacement) of the machine body; a first smoothing processing section 65 that subjects the position as a computation result of the positioning system 30 to the smoothing processing according to the satellite positioning environment (first smoothing processing); and a second smoothing processing section 66 that performs smoothing processing according to a determination result of the motion determining section 64 (second smoothing processing). The position correction computing section 61 outputs the corrected position as a computation result to the three-dimensional posture computing section 62.
The motion determining section 64 determines the presence or absence of motion of the machine body (upper swing structure 3) based on the angular velocity and acceleration of the machine body (upper swing structure 3) as sensing results of the machine body IMU 25. The motion determining section 64 outputs a determination result to the second smoothing processing section 66. When the machine body is in motion, for example when the machine body is travelling or swinging, the actual position of the GNSS antennas 31 and 32 attached to the machine body changes with the motion of the machine body. On the other hand, it is assumed that the actual position of the GNSS antennas 31 and 32 does not change when the machine body is motionless (stationary).
The motion determining section 64 for example determines that the machine body is in motion when the angular velocity and the acceleration as sensing results of the machine body IMU 25 each exceed corresponding threshold values, whereas otherwise (when equal to or lower than the threshold values), the motion determining section 64 determines that the machine body is motionless (stationary). More specifically, it is determined that the machine body is in motion when a composite value Acc of three-axis accelerations sensed by the machine body IMU 25 does not coincide with a gravitational acceleration g, or when an angular velocity of any axis of three axes is not zero. However, because the sensed values of the machine body IMU 25 have a bias (steady error) and include sensor noise, threshold values are desirably set for the determination. For example, it is determined that the machine body is in motion when any of the following inequations of (Formula 1) to (Formula 5) holds, where an acceleration threshold value is set as Ath, and an angular velocity threshold value is set as ωth.
Incidentally, ωx, ωy, and ωz are of angular velocities of an x-axis, a y-axis, and a z-axis as sensed values of the machine body IMU 25, respectively.
The first smoothing processing section 65 subjects a position p (position of the GNSS antennas 31 and 32) as one computation result of the positioning system 30 to the first smoothing processing that is to increase a degree of smoothing as a position variance value Rp as another result of the positioning computation of the positioning system 30 becomes larger. As described earlier, when the index (DOP) indicating the quality of the satellite positioning environment for the positioning system 30 is degraded, variation in the position p as a result of the positioning computation of the positioning system 30 is increased, and the position variance value Rp as another result of the positioning computation is also increased. The first smoothing processing is to utilizes this fact, and to change the strength of smoothing according to the magnitude of the variance value Rp as a computation result of the positioning system 30 to suppress the variation in the position p as a computation result of the positioning system 30 caused according to the satellite positioning environment. The first smoothing processing section 65 outputs a first corrected position pc1 obtained by the first smoothing processing to the second smoothing processing section 66.
Specifically, as shown in
In addition, a second example of the first smoothing processing section 65 can, for example, be constituted by a Kalman filter (hereinafter referred to as a KF) that uses the position as a computation result of the positioning system 30 and the variance thereof. Using of the KF allows the position p as a computation result of the positioning system 30 to be subjected to the first smoothing processing that is to increase the degree of smoothing as the variance value Rp as a computation result of the positioning system 30 becomes larger.
The KF generally handles a state equation of (Formula 6) and an observation equation of (Formula 7) in the following.
Here, xk is a state vector, yk is an observation vector, wk-1 is process noise, vk is observation noise, Fk-1 is a transition matrix, and Hk is an observation matrix. Incidentally, suffixes k and k-1 denote time. Suppose that the variance of process noise wk is given by Qk, and the variance of the observation noise vk is given by Rk.
After the above preparation, an update equation of the KF can be given by (Formula 8), and a gain used for an update can be given by (Formula 9).
Incidentally, a prediction value ^xk|k-1 (“^” is a circumflex) can be given by (Formula 10) as the following equation with respect to an update value “xk-1|k-1 preceding by one time interval.
In addition, supposing that a position vector at a certain time k is pk, and supposing that the position pk-1 at time k-1 preceding by one time interval is driven by an effect of certain noise, the state equation can be expressed by the dynamics of (Formula 11).
Here, wpk-1 is process noise related to arrangement. The variance of this process noise wpk-1 corresponds to Qk-1 in the above description.
In addition, the observation equation can be given by (Formula 12) in the following, supposing that position information pk at time k is affected by observation noise vpk.
Here, the variance of the observation noise vpk corresponds to Rk in the above description.
(Formula 11) is obtained supposing that the state vector is xk = pk, the transition matrix is Fk-1 = I, and the process noise is wk-1 = wpk-1 in (Formula 6). In addition, (Formula 12) is obtained supposing that the observation matrix is Hk = I, and the observation noise is vk = vpk in (Formula 7). Incidentally, I is a unit matrix.
As a result of the above preparation, the KF can be applied as the first smoothing processing section 65 that performs the first smoothing processing, which is to increase the degree of smoothing as the variance value as a computation result of the positioning system 30 becomes larger.
Here, the variance Rk is a certain value set by a designer when the KF is designed. In general, the variance Rk is often set at a fixed value. On the other hand, in the present embodiment, the variance as a computation result of the positioning system 30 is used as the variance Rk of the KF. Therefore, in the KF applied as the first smoothing processing section 65 according to the present embodiment, the variance Rk is sequentially changed according to changes in the satellite positioning environment.
According to (Formula 8), at each time, a value obtained by multiplying the position information yk received from the positioning system 30 by the update gain Kk is reflected in the estimated value ^xk|k. In addition, according to (Formula 9), the update gain Kk is calculated on the basis of an inverse matrix (division) of the variance Rk of the observation noise. This means that when the variance Rk of the observation noise (variance as a result of the positioning computation of the positioning system 30) becomes larger, the update gain Kk is correspondingly decreased. From the above, it is understood that the position information yk from the positioning system 30 is less reflected in the estimated value ^xk|k when the variance Rk of the observation noise becomes larger. That is, this KF can realize the first smoothing processing that is to increase the degree of smoothing as the variance value as a computation result of the positioning system 30 becomes larger, and can provide a smoothing effect similar to that of the configuration including the time constant changing LPF 652 in the first example shown in
In addition, a first modification of the second example of the first smoothing processing section 65 can be constituted by a KF that uses also speed information in addition to the position as input information from the positioning system 30.
Supposing that a position vector is pk, that a speed vector is vk, and that a computation cycle of the controller 40 is Δt, dynamics can be expressed by (Formula 13) in the following.
Here, wpk-1 is process noise related to the position, and wvk-1 is process noise related to the speed.
The variance Qk-1 of these process noises is given by (Formula 14) in the following.
Incidentally, Qpk-1 is the variance of wpk-1, and Qvk-1 is the variance of wvk-1.
In addition, the observation equation can be given by (Formula 15) in the following.
Here, vvk is observation noise related to the speed.
The variance Rk of the observation noises related to the position and the speed is given by (Formula 16) in the following.
Incidentally, Rpk is the variance of vpk, and Rvk is the variance of vvk. The variances of computation results of the positioning system 30 are used as Rpk and Rvk.
(Formula 13) is obtained when the state vector xk, the transition matrix Fk, and the process noise wk-1 in (Formula 6) as a state equation are set as follows.
In addition, (Formula 15) is obtained when the observation matrix Hk and the observation noise vk in (Formula 7) as an observation equation are set as follows.
From the above, the KF adopting the above-described equations can be applied as the first smoothing processing section 65 with a configuration that uses both of the position and the speed as computation results of the positioning system 30. A characteristic of not reflecting the position as a computation result of the positioning system 30 when the satellite positioning environment is not excellent (the variance of the positioning result is large) is maintained by using, as the variance Rk of the observation noises of the KF, the variances as computation results of the positioning system 30. That is, this KF smooths the position of the positioning result to which the speed as a computation result of the positioning system 30 is related after changing the strength of smoothing according to the variances of the position and the speed as another result of the positioning computation.
Incidentally, in (Formula 13), unlike (Formula 11), speed information is considered in computing the position. Therefore, even when the positions as results of the positioning computation momentarily vary, grounds that the actual position is not changed are strengthened when the speeds as results of the positioning computation do not vary. The first modification of the second example (KF) of the first smoothing processing section 65 is therefore expected to increase the effect of suppressing variation in the position as results of the positioning computation.
In addition, as a second modification of the second example including the KF of the first smoothing processing section 65, a KF is applied which uses information of both of results (information of position and speed) of the positioning computation of the positioning system 30 and a sensing result (acceleration) of the machine body IMU 25. This considers the dynamics of (Formula 17) in which an acceleration vector ak sensed by the IMU is added to the position vector pk and the speed vector vk.
However, the positioning system 30 cannot compute the acceleration of the machine body (GNSS antennas 31 and 32). Thus, the observation equation is set as (Formula 18).
In this case, it is to be noted that the variance Rk of the observation noises is set to be the same as in (Formula 16).
The first smoothing processing that is to increase the degree of smoothing as the variance value as a computation result of the positioning system 30 becomes larger can be realized by using the dynamics of (Formula 17) to configure the KF according to the above-described method. This KF takes into consideration information that the machine body IMU 25 does not sense acceleration when the machine body (upper swing structure 3) is motionless (stationary). Thus, the KF can be expected to further improve the effect of suppressing variation in the position as computation results of the positioning system 30, which occurs according to a change in the satellite positioning environment.
It is to be noted that the second example of the first smoothing processing section 65 is not limited to the KF as long as the second example of the first smoothing processing section 65 uses certain dynamics as in (Formula 6) and (Formula 7). For example, when nonlinearity is included in the dynamics, a particle filter or a GSF (Gaussian Sum Filter) may be used. These filters all perform statistical processing, and have a function of similar smoothing processing in that the position information yk is less reflected in the estimated value ^xk|k as the variance value as a computation result of the positioning system 30 becomes larger. In the present description, filters that perform smoothing processing including these pieces of statistical processing described above will be referred to as statistical filters.
The second smoothing processing section 66 subjects the first corrected position pc1 as a processing result of the first smoothing processing section 65 to the second smoothing processing that is to change the strength of smoothing according to the presence or absence of motion of the machine body (the lower track structure 2 and the upper swing structure 3) as the determination result of the motion determining section 64. The second smoothing processing section 66 outputs a second corrected position pc2 as a processing result to the three-dimensional posture computing section 62. The second smoothing processing section 66 makes smoothing stronger when the motion determining section 64 determines that the machine body is motionless (stationary) than when the motion determining section 64 determines that the machine body is in motion. The actual position of the GNSS antennas 31 and 32 is assumed to be unchanged when the machine body is stationary. In this case, relatively strengthened smoothing of the input position information (first corrected position pc1) suppresses variation in the position information. When the machine body is in motion, on the other hand, the actual position of the GNSS antennas 31 and 32 changes with the motion of the machine body. In this case, the second corrected position pc2 output from the second smoothing processing section 66 to the three-dimensional posture computing section 62 is desired to follow the actual displacement of the GNSS antennas 31 and 32. Accordingly, relatively weakened smoothing of the input position information (first corrected position pc1) suppresses a response delay in the second corrected position pc2 as the processing result of the second smoothing processing section 66.
Specifically, as shown in
In addition, as shown in
As shown in
Next, referring to
In
Next, the first smoothing processing section 65 of the positioning computing section 51 subjects the position as a computation result of the positioning system 30 to the first smoothing processing that is to increase the degree of smoothing as the variance value of a positioning computation of the positioning system 30 becomes larger (step S20).
Next, the motion determining section 64 of the positioning computing section 51 determines whether or not the machine body is in motion (step S30). When the motion determining section 64 determines that the machine body is in motion (YES), the processing proceeds to step S40. When the motion determining section 64 determines that the machine body is motionless (NO), on the other hand, the processing proceeds to step S50.
In step S40, the second smoothing processing section 66 of the positioning computing section 51 subjects the first corrected position as a processing result of the first smoothing processing section 65 to relatively weakened smoothing processing. In the case of the first example of the second smoothing processing section 66 shown in
In addition, in step S50, the second smoothing processing section 66 of the positioning computing section 51 subjects the first corrected position as the processing result of the first smoothing processing section 65 to relatively strengthened smoothing processing. In the case of the first example of the second smoothing processing section 66 shown in
In step S60, the three-dimensional posture computing section 62 of the positioning computing section 51 computes the posture information (three-dimensional coordinates of various kinds of parts, which are necessary for the machine guidance and the machine control, such as the central position of the claw tip of the bucket 8) of the hydraulic excavator by using the second corrected position pc2 as a processing result of the second smoothing processing section 66.
A return is made to the start after the execution of the processing of step S60. The positioning computing section 51 can obtain highly accurate position information necessary for the machine guidance and the machine control by repeatedly performing the computation cycle of steps S10 to S60 many times.
Next, actions and effects of the construction machine according to the first embodiment of the present invention will be described while compared with a technology according to a comparative example. First, referring to
On the other hand,
It is to be noted that while the position information as a processing result of the first smoothing processing section 65 according to the present embodiment and the position information as a processing result of the technology according to the comparative example represent similar results, computation contents differ. Next, referring to
In order to resolve the response delay caused by the smoothing processing of the technology according to the comparative example and make the corrected position follow the actual displacement, the smoothing processing needs to be weakened. However, the smoothing processing of the technology according to the comparative example does not distinguish the change in the positioning result of the positioning system 30, which is caused by the degradation in the satellite positioning environment as shown in
On the other hand,
Next, a processing result of the first smoothing processing section according to the present embodiment in a case where the condition shown in
On the other hand,
As described above, the hydraulic excavator (construction machine) according to the first embodiment of the present invention includes: the machine body (the lower track structure 2 and the upper swing structure 3); the front work device 1 (work device) attached to the machine body 2 and 3 in a raiseable and lowerable manner; the GNSS antennas 31 and 32 (antenna) that are attached to the machine body 2 and 3, and receive positioning signals from the plurality of satellites; the machine body IMU 25 (first sensor) that senses information on the posture and motion of the machine body 2 and 3; the IMUs 26, 27, and 28 (second sensor) that sense information on the posture of the front work device 1 (work device); and the computing device (the GNSS receiver 33 and the controller 40) that computes posture information indicating the postures of the machine body 2 and 3 and the front work device 1 (work device). The computing device (the GNSS receiver 33 and the controller 40) is configured to perform positioning computation that is to compute the position of the machine body 2 and 3 and the variance value of the position of the machine body 2 and 3 based on the positioning signals from the plurality of satellites received by the GNSS antennas 31 and 32 (antenna), subject the position of the machine body 2 and 3 obtained by the positioning computation to the first smoothing processing that is to increase the degree of smoothing as the magnitude of the variance value of the position of the machine body 2 and 3 obtained by the positioning computation becomes larger, and compute the posture information based on a result of the first smoothing processing, the information sensed by the machine body IMU 25 (first sensor), and the information sensed by the IMUs 26, 27, and 28 (second sensor).
According to this configuration, even when variation in the position as a result of the positioning computation becomes large due to a degradation in the satellite positioning environment, the variance value as a result of the positioning computation is increased, and therefore, by correspondingly strengthening smoothing on the position as a result of the positioning computation, it is possible to obtain position information in which the variation in the position as a result of the positioning computation is suppressed. In addition, when a movement of the machine body 2 and 3 occurs without operator’s operation, the variance value as a result of the positioning computation hardly changes unless the satellite positioning environment is degraded, and therefore position information following the movement of the machine body 2 and 3 can be obtained by smoothing the position as a result of the positioning computation with a strength corresponding to the magnitude of the variance value as a result of the positioning computation. That is, even when variation occurs in computation result of the satellite positioning or even when a movement of the own machine occurs without operator’s operation, it is possible to obtain the posture information of the hydraulic excavator (construction machine) based on accurate position information of the own machine.
In addition, the positioning computation of the GNSS receiver 33 (part of the computing device) according to the present embodiment is further to compute the speed of the machine body 2 and 3 and the variance value of the speed of the machine body 2 and 3 on the basis of the positioning signals from the plurality of satellites received by the GNSS antennas 31 and 32. In addition, the first smoothing processing of the controller 40 (part of the computing device) is to set, as input information, the position of the machine body 2 and 3 and the speed of the machine body 2 and 3 obtained by the positioning computation, and to perform smoothing processing with the strength of smoothing changed according to the magnitude of the variance value of the position of the machine body 2 and 3 obtained by the positioning computation, and the magnitude of the variance value of the speed of the machine body 2 and 3 obtained by the positioning computation.
According to this configuration, information on the speed and the variance value thereof as results of the positioning computation is added in the first smoothing processing. Thus, in addition to the presence or absence of reliability of the position as a result of the positioning computation, the presence or absence of reliability of the speed as a result of the positioning computation is reflected in the result of the first smoothing processing. Hence, when a position as a result of the positioning computation changes momentarily, and there is no variation in the speed as a result of the positioning computation, probability that there is no actual displacement is increased, so that an effect of suppressing variation in the position information can be enhanced.
In addition, the controller 40 (part of the computing device) according to the present embodiment further makes a motion determination that is to determine whether or not the machine body 2 and 3 is in motion based on a sensing result of the machine body IMU 25 (first sensor), and performs the second smoothing processing that is to further smooth the first smoothing processing. The posture information of the controller 40 is computed based on the result of the second smoothing processing, the information sensed by the machine body IMU 25 (first sensor), and the information sensed by the IMUs 26, 27, and 28 (second sensor). In addition, the strength of smoothing of the second smoothing processing is set to be stronger in a case where it is determined in the motion determination that the machine body 2 and 3 is motionless than in a case where it is determined in the motion determination that the machine body 2 and 3 is in motion.
According to this configuration, by changing the strength of smoothing according to the presence or absence of motion of the machine body 2 and 3, it is possible to suppress variation in result of the positioning computation when the machine body 2 and 3 is motionless, and obtain position information following an actual displacement during motion of the machine body 2 and 3 (suppress a response delay caused by smoothing processing).
In addition, the first smoothing processing of the controller 40 (part of the computing device) according to the present embodiment uses a statistical filter (for example a KF, a GSF, or a particle filter) that performs smoothing processing including statistical processing, and the variance value as a result of the positioning computation is used as an input value of the statistical filter.
According to this configuration, this statistical filter inherently has a function of changing the strength of smoothing according to the magnitude of the variance value. Thus, the configuration including the statistical filter is made simpler than that including a filter (for example, an LPF) that does not perform statistical processing and does not inherently have the function.
A construction machine according to a modification of the first embodiment of the present invention will next be described with reference to
The construction machine according to the modification of the first embodiment of the present invention is different from the first embodiment in that the positioning computation performed by the GNSS receiver 33 (see
A correcting computation (smoothing processing) of the position correction computing section 61 on the computation results of the positioning computing section 56 is similar to that in the case of the first embodiment. In addition, a procedure for computation processing of the positioning computing section 51 of the controller 40A is similar to the flowchart showing the procedure for the computation processing of the positioning computing section 51 according to the first embodiment (see
Incidentally, the GNSS antennas 31 and 32 and the positioning computing section 56 of the controller 40A constitute a positioning system 30A. In addition, the controller 40A including the positioning computing section 56 constitutes a computing device that computes the posture information of the hydraulic excavator.
According to the above-described construction machine according to the modification of the first embodiment of the present invention, as in the foregoing first embodiment, it is possible to obtain the posture information of the hydraulic excavator (construction machine) based on accurate position information of the own machine even when variation occurs in computation result of the satellite positioning or even when a displacement of the own machine occurs without operator’s operation.
Construction machines according to a first example and a second example of a second embodiment of the present invention will next be described with reference to
The construction machine according to the second embodiment of the present invention shown in
When the first smoothing processing section 65B of the construction machine according to the first example of the second embodiment of the present invention shown in
Specifically, as shown in
Incidentally, any configuration is possible as long as the configuration can implement a correction that is to increase the variance values as computation results of the positioning system 30 only when the front work device 1 is assumed to be positioned in the specific region. For example, a configuration is possible in which the selecting section 654 outputs zero at a time of boom lowering, and outputs the positive number Kup at a time of boom raising, and an adding section is used in place of the multiplying section 655. In addition, the KF 656 according to the present embodiment can be replaced with the time constant changing LPF 652 shown in
Incidentally, when the hydraulic excavator travels while maintaining the boom-raised posture, the first smoothing processing section 65B according to the first example makes a correction that is to increase the variance values as computation results of the positioning system 30 according to the boom-raised posture. In this case, a result is obtained such that smoothing on the position as a computation result of the positioning system 30 is strengthened though the actual position is changing due to travelling. In this case, the position corrected by the first smoothing processing section 65B may be shifted with respect to the actual displacement of the GNSS antennas 31 and 32.
Accordingly, the first smoothing processing section 65B according to the second example of the second embodiment shown in
Specifically, the first smoothing processing section 65B according to the second example, for example, includes: a gain table 657 that computes a gain on the basis of the sensing result of the boom IMU 26 and the determination result of the motion determining section 64; a multiplying section 655 that multiplies the variance values ∑p and ∑v as computation results of the positioning system 30 by an output value (gain) from the gain table 657; and a KF 656 that sets, as input values thereof, information of the position p and speed v as computation results of the positioning system 30 and an output value from the multiplying section 655. The KF 656 is the same as the KF in the first example of the present embodiment. The gain table 657 is, for example, illustrated in the following Table 1. The positive number Kup larger than one is output when the posture of the front work device 1 is a boom-raised posture and the determination result of the motion determining section 64 represents that the machine body is motionless. One is output otherwise, that is, when the posture of the front work device 1 is a boom-lowered posture, or when the posture of the front work device 1 is a boom-raised posture and the determination result of the motion determining section 64 represents that the machine body is in motion.
Thus, only when the posture of the front work device 1 is a boom-raised posture (the satellite positioning environment is degraded) and the machine body is motionless (there is no displacement of the GNSS antennas 31 and 32), the variance values as computation results of the positioning system 30 are corrected to be larger. On the other hand, otherwise, the variances as computation results of the positioning system 30 are not corrected. Hence, when the hydraulic excavator travels while maintaining a boom-raised posture, the variances as computation results of the positioning system 30 are input to the filter performing smoothing processing as they are not corrected to be larger, so that the position information corrected by the first smoothing processing section 65 follows the actual position displacement.
According to the above-described construction machine according to the first example and second example of the second embodiment of the present invention, as in the foregoing first embodiment, it is possible to obtain the posture information of the hydraulic excavator (construction machine) based on accurate position information of the own machine even when variation occurs in computation result of the satellite positioning or even when a displacement of the own machine occurs without operator’s operation.
In addition, the controller 40B (part of the computing device) according to the first example of the present embodiment further determines on the basis of a sensing result of the IMUs 26, 27, and 28 (second sensor) whether or not the front work device 1 (work device) is positioned in the specific region in which the front work device 1 degrades the reception of the positioning signals on the GNSS antennas 31 and 32. In addition, when it is determined that the front work device 1 (work device) is positioned in the specific region, the first smoothing processing of the controller 40B corrects the variance value of the position of the machine body 2 and 3 obtained by the positioning computation such that the variance value of the position becomes larger than in a case where the front work device 1 (work device) is not positioned in the specific region, and then performs smoothing processing with a strength corresponding to the magnitude of the corrected variance value.
According to this configuration, variation in the positioning result of the positioning system 30 caused by a degradation in the satellite positioning environment can be further suppressed by determining the condition of the satellite positioning environment based on the posture of the front work device 1 (work device).
In addition, the controller 40B (part of the computing device) according to the second example of the present embodiment further makes a motion determination that is to determine whether or not the machine body 2 and 3 is in motion based on a sensing result of the machine body IMU 25 (first sensor). In addition, the first smoothing processing of the controller 40B is to perform smoothing with a strength corresponding to the magnitude of the variance value of the position of the machine body obtained by the positioning computation, without correcting the variance value of the position of the machine body obtained by the positioning computation, when it is determined that the front work device 1 (work device) is positioned in the specific region and it is determined that the machine body 2 and 3 is in motion. The first smoothing processing of the controller 40B is to correct, when it is determined that the front work device 1 (work device) is positioned in the specific region and it is determined that the machine body 2 and 3 is motionless, the variance value of the position of the machine body 2 and 3 obtained by the positioning computation such that the variance value of the position becomes larger than in the case where the front work device 1 (work device) is not positioned in the specific region, and then perform smoothing with a strength corresponding to the magnitude of the corrected variance value.
According to this configuration, in other than a case where the satellite positioning environment is degraded and the position of the GNSS antennas 31 and 32 is not changed, the variance value as a result of the positioning computation is not corrected to be larger. Thus, a response delay in the processing result (corrected position) of the first smoothing processing section 65 can be prevented when the hydraulic excavator travels while maintaining a boom-raised posture.
A construction machine according to the third embodiment of the present invention will next be described with reference to
The construction machine according to the third embodiment of the present invention shown in
In addition, when swing motion is performed at an extremely slow speed, changes in the sensed values of the machine body IMU 25 are very small, and therefore it may be difficult to determine the swing motion even when any of the conditional expressions of (Formula 1) to (Formula 5) is used. On the other hand, the motion determining section 64C determines the presence or absence of motion of the machine body based on operation of the operation devices 18a and 18b, and therefore swing motion at an extremely slow speed can be determined correctly. Consequently, excessive smoothing of the second smoothing processing section 66 can be prevented. It is therefore possible to suppress a response delay caused by the smoothing of the second smoothing processing section 66.
According to the above-described construction machine according to the third embodiment of the present invention, as in the foregoing first embodiment, it is possible to obtain the posture information of the hydraulic excavator (construction machine) based on accurate position information of the own machine even when variation occurs in computation result of the satellite positioning or even when a movement of the own machine occurs without operator’s operation.
In addition, the hydraulic excavator (construction machine) according to the present embodiment further includes: the operation devices 18a and 18b for operating the machine body (the lower track structure 2 and the upper swing structure 3); and a third sensor (operation devices 18a and 18b) that senses operation of the operation devices 18a and 18b. Further, the motion determination of the controller 40C (part of the computing device) is to determine that the machine body 2 and 3 is in motion when the machine body IMU 25 (first sensor) senses motion of the machine body 2 and 3 or when the operation devices 18a and 18b (third sensor) senses operation of the operation devices 18a and 18b.
According to this configuration, motion of the machine body 2 and 3 “a few time intervals ahead” is takes into consideration by sensing operation of the operation devices 18a and 18b, which reflects an intention of the operator who operates the machine body 2 and 3. This allows the second smoothing processing section 66, which changes the strength of smoothing according to the presence or absence of motion of the machine body 2 and 3, to weaken smoothing before the machine body 2 and 3 is actually displaced. Consequently, a response delay caused by the smoothing of the second smoothing processing section 66 can be suppressed. It is therefore possible to obtain position information that follows an actual displacement of the machine body 2 and 3.
It is to be noted that while an example in which the present invention is applied to a hydraulic excavator has been illustrated in the foregoing embodiments, the present invention is widely applicable to various kinds of construction machines that control, on the basis of the position information of the own machines, a work device or the like that changes in position and posture.
In addition, the present invention is not limited to the present embodiments, but includes various modifications. The foregoing embodiments have been described in detail in order to describe the present invention in an easily understandable manner, and are not necessarily limited to embodiments including all of the described configurations. A part of a configuration of a certain embodiment can be replaced with a configuration of another embodiment, and a configuration of another embodiment can be added to a configuration of a certain embodiment. In addition, for a part of a configuration of each embodiment, another configuration can be added, deleted, or substituted.
For example, in the foregoing embodiments, an example has been illustrated in which the operation devices 18a and 18b are constituted by operation devices of an electric type. However, the operation devices can also be constituted by operation devices of a hydraulic type. In a case of this configuration, the hydraulic actuators 4, 10, 11, 12, and 14a are driven by supplying an operation pilot pressure corresponding to the operation direction and the operation amount of each operation device as a driving signal to the control valves of the control valve unit 23. The operation direction and the operation amount of the operation device can be sensed by a pressure sensor that senses the operation pilot pressure generated by the operation device. That is, the pressure sensor functions as a sensor that senses operation (the operation direction and the operation amount) of the operation device.
In addition, in the foregoing second embodiment, a configuration has been illustrated in which it is determined, on the basis of a sensing result of the boom IMU 26, whether or not the front work device 1 is positioned in the specific region where the front work device 1 degrades the reception of the positioning signals on the GNSS antennas 31 and 32. However, a configuration is also possible in which it is determined, on the basis of sensing results of the three IMUs 26, 27, and 28 that sense the posture of the front work device 1, whether or not the front work device 1 is positioned in the specific region.
1:
2:
3:
25:
26:
27:
28:
31, 32:
33:
40:
18
a, 18b:
Number | Date | Country | Kind |
---|---|---|---|
2020-103572 | Jun 2020 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2021/021974 | 6/9/2021 | WO |