The present invention relates to a control system for a construction machine, the construction machine, and a control method for the construction machine.
In a technical field related to construction machines, there is known a control system for a construction machine that moves a bucket of working equipment along a design surface that indicates a target shape of a construction target, as disclosed in Patent Literature 1.
The construction machine has a swing body. The working equipment is supported by the swing body. In some cases, the construction machine has a tilt bucket that is capable of tilt rotation. The construction machine is configured to form a slope on the construction target only by operating the working equipment so that the tilt bucket approaches the swing body. In some cases, when forming the slope on the construction target by the tilt bucket, the working equipment is operated while swinging the swing body so as to prevent the tilt bucket from digging below the design surface. When the swing body swings on the basis of operation of an operation device by a driver, if the swing body has a high swing speed, the tilt bucket may dig below the design surface, preventing excavation of the construction target into a desired shape.
An object of an aspect of the present invention is to excavate a construction target into a desired shape.
According to an aspect of the present invention, a control system for a construction machine including working equipment that includes an arm and a tilt bucket, and a swing body supporting the working equipment, the control system comprises: a design surface acquisition unit configured to acquire a design surface indicating a target shape of a construction target; an operation data acquisition unit configured to acquire operation data of an operation device operating the working equipment and the swing body; a target value generation unit configured to generate a target value of an amount of control relating to tilt rotation of the tilt bucket and a target value of an amount of control of the swing body, based on the operation data obtained when the working equipment and the swing body are operated and the design surface; a prediction unit configured to calculate a predicted value of the amount of control relating to the tilt rotation of the tilt bucket and a predicted value of an amount of control relating to swing of the swing body, based on the target values and a prediction model for the construction machine, and configured to calculate amounts of drive to control the working equipment and the swing body, based on the predicted values and the design surface so as to maintain a distance from a predetermined portion of the tilt bucket to the design surface and an attitude of the tilt bucket; and a command unit configured to output control commands to control the working equipment and the swing body based on the amounts of drive.
According to the aspect of the present invention, the construction target can be excavated into the desired shape.
Embodiments according to the present invention will be described below with reference to the drawings, but the present invention is not limited to this description. Component elements according to the embodiments described below can be appropriately combined with each other. Furthermore, in some cases, some of the component elements are not used.
In the following description, a three-dimensional vehicle body coordinate system (X, Y, Z) is defined, and a positional relationship between respective component elements will be described. The vehicle body coordinate system represents a coordinate system that has the origin fixed to a construction machine. The vehicle body coordinate system is defined by an X-axis that extends in a defined direction based on the origin set to the construction machine, a Y-axis that is orthogonal to the X-axis, and a Z-axis that is orthogonal to each of the X-axis and the Y-axis. A direction parallel to the X-axis is an X-axis direction. A direction parallel to the Y-axis is a Y-axis direction. A direction parallel to the Z-axis is a Z-axis direction. A rotation or tilting direction about the X-axis is a OX direction. A rotation or tilting direction about the Y-axis is a OY direction. A rotation or tilting direction about the Z-axis is a OZ direction.
[Construction Machine]
As illustrated in
The traveling body 3 has a pair of tracks 3C. The excavator 100 travels by the rotation of the tracks 3C. Note that the traveling body 3 may have tires.
The working equipment 1 is supported by the swing body 2. The working equipment 1 has a boom 6 that is connected to the swing body 2, an arm 7 that is connected to a distal end of the boom 6, and a bucket 8 that is connected to a distal end of the arm 7. The bucket 8 has a blade edge 9. In the present embodiment, the blade edge 9 of the bucket 8 is the edge of a straight blade. The blade edge 9 of the bucket 8 may be the tip of a protruded tooth provided at the bucket 8.
The boom 6 is rotatable about a boom axis AX1 relative to the swing body 2. The arm 7 is rotatable about an arm axis AX2 relative to the boom 6. In the present embodiment, the bucket 8 is a tilt bucket. The bucket 8 is rotatable relative to the arm 7 about each of a bucket axis AX3, a tilt axis AX4, and a rotating axis AX5. The boom axis AX1, arm axis AX2, and bucket axis AX3 are parallel to the Y-axis. The tilt axis AX4 is orthogonal to the bucket axis AX3. The rotating axis AX5 is orthogonal to each of the bucket axis AX3 and the tilt axis AX4. The swing axis RX is parallel to the Z-axis. The X-axis direction is a front-rear direction of the swing body 2. The Y-axis direction is a vehicle width direction of the swing body 2. The Z-axis direction is a vertical direction of the swing body 2. A direction of the working equipment 1 relative to the driver seated on the seat 4S is a front direction.
As illustrated in
The working equipment 1 is operated by power generated by the hydraulic cylinders 10. Each of the hydraulic cylinders 10 is driven on the basis of the hydraulic oil supplied from the hydraulic pump 17. The hydraulic cylinder 10 includes a boom cylinder 11 that is configured to operate the boom 6, an arm cylinder 12 that is configured to operate the arm 7, and a bucket cylinder 13, a tilt cylinder 14, and a rotating cylinder 15 that are configured to operate the bucket 8. The boom cylinder 11 generates power to rotate the boom 6 about the boom axis AX1. The arm cylinder 12 generates power to rotate the arm 7 about the arm axis AX2. The bucket cylinder 13 generates power to rotate the bucket 8 about the bucket axis AX3. The tilt cylinder 14 generates power to rotate the bucket 8 about the tilt axis AX4. The rotating cylinder 15 generates power to rotate the bucket 8 about the rotating axis AX5.
In the following description, the rotation of the bucket 8 about the bucket axis AX3 is appropriately referred to as bucket rotation, and the rotation of the bucket 8 about the tilt axis AX4 is appropriately referred to as tilt rotation, and the rotation of the bucket 8 about the rotating axis AX5 is appropriately referred to as rotation.
The swing body 2 swings by power generated by the swing motor 16. The swing motor 16 is a hydraulic motor and is driven on the basis of the hydraulic oil supplied from the hydraulic pump 17. The swing motor 16 generates power to cause the swing body 2 to swing about the swing axis RX.
The engine 5 is mounted on the swing body 2. The engine 5 generates power for driving the hydraulic pump 17.
The hydraulic pump 17 discharges the hydraulic oil for driving the hydraulic cylinder 10 and the swing motor 16.
The valve device 18 has a plurality of valves configured to distribute the hydraulic oil supplied from the hydraulic pump 17 to the plurality of hydraulic cylinders 10 and the swing motor 16. The valve device 18 adjusts the flow rate of the hydraulic oil supplied to each of the plurality of hydraulic cylinders 10. The adjustment of the flow rate of the hydraulic oil supplied to the hydraulic cylinder 10 adjusts the operating speed of the hydraulic cylinder 10. The valve device 18 regulates the flow rate of the hydraulic oil supplied to the swing motor 16. The adjustment of the flow rate of the hydraulic oil supplied to the swing motor 16 adjusts the rotation speed of the swing motor 16.
The position calculation device 20 calculates the positional data of the swing body 2. The positional data of the swing body 2 includes the position of the swing body 2, the attitude of the swing body 2, and the orientation of the swing body 2. The position calculation device 20 has a position calculator 21 configured to calculate the position of the swing body 2, an attitude calculator 22 configured to calculate the attitude of the swing body 2, and an orientation calculator 23 configured to calculate the orientation of the swing body 2.
The position calculator 21 calculates the position of the swing body 2 in a global coordinate system, as the position of the swing body 2. The position calculator 21 is arranged in the swing body 2. The global coordinate system represents a coordinate system that has the origin fixed to the earth. The global coordinate system is a coordinate system defined by a Global Navigation Satellite System (GNSS). The GNSS is a global navigation satellite system. A global positioning system (GPS) is an example of the global navigation satellite system. The GNSS has a plurality of positioning satellites. The GNSS detects a position defined by coordinate data of a latitude, a longitude, and an altitude. The swing body 2 is provided with a GPS antenna. The GPS antennas receives a radio wave from a GPS satellite and outputs a signal generated on the basis of the received radio wave to the position calculator 21. The position calculator 21 calculates the position of the swing body 2 in the global coordinate system, on the basis of the signal supplied from the GPS antenna. The position calculator 21 calculates the position of a representative point O of the swing body 2, for example, as illustrated in
The attitude calculator 22 calculates the inclination angle of the swing body 2 relative to a horizontal plane in the global coordinate system, as the attitude of the swing body 2. The attitude calculator 22 is arranged in the swing body 2. The attitude calculator 22 includes an inertial measurement unit (IMU). The inclination angle of the swing body 2 relative to the horizontal plane includes a roll angle α indicating an inclination angle of the swing body 2 in the vehicle width direction and a pitch angle β indicating an inclination angle of the swing body 2 in the front-rear direction.
The orientation calculator 23 calculates the orientation of the swing body 2 relative to a reference orientation in the global coordinate system, as the orientation of the swing body 2. The reference orientation is, for example, north. The orientation calculator 23 is arranged in the swing body 2. The orientation calculator 23 includes a gyroscope sensor. The orientation calculator 23 may calculate the orientation on the basis of the signal supplied from the GPS antenna. The orientation of the swing body 2 relative to the reference orientation includes a yaw angle γ indicating an angle formed by the reference orientation and the orientation of the swing body 2.
The angle detection device 30 detects the angles θ of the working equipment 1. The angle detection device 30 is arranged in the working equipment 1. As illustrated in
The angle detection device 30 has a boom angle detector 31 that detects the boom angle θ1, an arm angle detector 32 that detects the arm angle θ2, a bucket angle detector 33 that detects the bucket angle θ3, a tilt angle detector 34 that detects the tilt angle θ4, and a rotation angle detector 35 that detects the rotation angle θ5. The angle detection device 30 may include a stroke sensor configured to detect a stroke of the hydraulic cylinder 10, or may include an angle sensor, such as a rotary encoder, configured to detect the angles θ of the working equipment 1. When the angle detection device 30 includes the stroke sensor, the angle detection device 30 calculates the angles θ of the working equipment 1 on the basis of detection data of the stroke sensor.
The operation device 40 is operated by the driver to drive the hydraulic cylinder 10 and the swing motor 16. The operation device 40 is arranged in the cab 4. Operation of the operation device 40 by the driver operates the working equipment 1. The operation device 40 includes a lever that is operated by the driver of the excavator 100. The lever of the operation device 40 includes a right operating lever 41, a left operating lever 42, and a tilt operating lever 43.
When the right operating lever 41 in a neutral position is operated forward, the boom 6 is operated to lower, and when the right operating lever 41 is operated backward, the boom 6 is operated to rise. When the right operating lever 41 in the neutral position is operated rightward, the bucket 8 dumps, and when the right operating lever 41 is operated leftward, the bucket 8 digs.
When the left operating lever 42 in a neutral position is operated forward, the arm 7 dumps, and when the left operating lever 42 is operated backward, the arm 7 digs. When the left operating lever 42 in the neutral position is operated rightward, the swing body 2 swings rightward, and when the left operating lever 42 is operated to leftward, the swing body 2 swings leftward.
Operation of the tilt operating lever 43 causes a tilt rotation or rotation of the bucket 8.
[Control Device]
The positional data acquisition unit 51 acquires the positional data of the swing body 2 from the position calculation device 20. The positional data of the swing body 2 includes the position of the swing body 2, the attitude of the swing body 2, and the orientation of the swing body 2.
The angle data acquisition unit 52 acquires angle data indicating the angles θ of the working equipment 1 from the angle detection device 30. The angle data of the working equipment 1 includes the boom angle θ1, the arm angle θ2, the bucket angle θ3, the tilt angle θ4, and the rotation angle θ5.
The operation data acquisition unit 53 acquires operation data of the operation device 40 operating the working equipment 1 and the swing body 2. The operation data of the operation device 40 includes an amount of operation of the operation device 40. The operation device 40 is provided with an operation amount sensor that detects an amount of operation of each lever. The operation data acquisition unit 53 acquires the operation data of the operation device 40 from the operation amount sensor of the operation device 40.
The design surface acquisition unit 54 acquires a design surface that indicates a target shape of a construction target. The design surface indicates a three-dimensional target shape after construction by the excavator 100. In the present embodiment, a design surface data supply device 70 generates design surface data that indicates the design surface. The design surface acquisition unit 54 acquires the design surface data from the design surface data supply device 70. The design surface data supply device 70 may be provided at a distant place from the excavator 100. The design surface data generated by the design surface data supply device 70 may be transmitted to the control device 50 via a communication system. Note that the design surface data generated by the design surface data supply device 70 may be stored in the storage unit 60. The design surface acquisition unit 54 may acquire the design surface data from the storage unit 60.
The target value generation unit 55 generates a target value of an amount of control of the working equipment 1. In the present embodiment, the amount of control of the working equipment 1 includes one or both of a moving speed of the bucket 8 and the position of a predetermined portion of the bucket 8. The predetermined portion of the bucket 8 includes the blade edge 9 of the bucket 8. The moving speed of the bucket 8 includes a moving speed of the blade edge 9. The position of the predetermined portion of the bucket 8 includes the position of the blade edge 9. The target value generation unit 55 generates the target value of the amount of control of the working equipment 1, on the basis of the operation data acquired by the operation data acquisition unit 53. Furthermore, the target value generation unit 55 generates a target value of an amount of control of the swing body 2. In the present embodiment, the amount of control of the swing body 2 includes a swing speed of the swing body 2. The target value generation unit 55 generates the target value of the amount of control of the swing body 2, on the basis of the operation data acquired by the operation data acquisition unit 53. In the present embodiment, the target value generation unit 55 generates a target value of an amount of control relating to the tilt rotation of the bucket 8 and a target value of an amount of control relating to the swing of the swing body 2, on the basis of at least the operation data obtained upon operation of the working equipment 1 and the swing body 2, and the design surface.
The amount of control relating to the tilt rotation of the bucket 8 includes at least one of the tilt angle θ4 indicating the position of the bucket 8 in a tilt direction, a tilt speed indicating a rotation speed of the bucket 8 in the tilt direction, and a tilt acceleration indicating a rotational acceleration of the bucket 8 in the tilt direction. The amount of control relating to the swing of the swing body 2 includes the swing speed of the swing body 2.
In the following description, it is assumed that the predetermined portion of the bucket 8 is the blade edge 9. Note that the predetermined portion of the bucket 8 may not necessarily be the blade edge 9. The predetermined portion of the bucket 8 may be a floor surface (bottom surface) of the bucket 8.
The moving speed of the bucket 8 includes a translational speed and a rotation speed of the bucket 8. The translational speed of the bucket 8 represents a moving speed in each of the X-axis direction, the Y-axis direction, and the Z-axis direction. The rotation speed of the bucket 8 represents a rotational angular speed in each of the OX direction, the OY direction, and the θZ direction. In the present embodiment, the target value generation unit 55 includes a target translational speed calculation unit 551 configured to calculate a target translational speed vtarget that is a target value of the translational speed, and a target rotation speed calculation unit 552 configured to calculate a target rotation speed ωtarget that is a target value of the rotation speed. The target value generation unit 55 calculates each of the target translational speed vtarget and the target rotation speed ωtarget, on the basis of the angle data acquired by the angle data acquisition unit 52, the operation data acquired by the operation data acquisition unit 53, and the design surface acquired by the design surface acquisition unit 54.
The target translational speed calculation unit 551 calculates the target translational speed vtarget of the bucket 8 so as not to dig below the design surface. The target translational speed vtarget of the bucket 8 is calculated on the basis of formulas (1) to (6).
n ∈ R3 is a unit normal vector of the design surface closest to the blade edge 9, wR1 ∈ R3×3 is a rotation matrix for transformation from the vehicle body coordinate system to the global coordinate system, vsagyo ∈ R3 is a translational speed component of the bucket operated by the boom 6 and the arm 7 on a working equipment plane (XZ plane in the vehicle body coordinate system) of the translational speed upon operation of the working equipment 1 based on the operation of the operation device 40, and VMAX is a maximum speed of the bucket 8 in a direction normal to the design surface, for prevention of digging below the design surface. Jv ∈ R3×5 and Jω ∈ R3×5 represent a translational speed component and a rotation speed component of the Jacobian matrix, respectively.
The target translational speed calculation unit 551 is configured to calculate the distance between the blade edge 9 and the design surface, on the basis of the positional data of the swing body 2 acquired by the positional data acquisition unit 51, the angle data of the working equipment 1 acquired by the angle data acquisition unit 52, and working equipment data stored in the storage unit 60. As illustrated in
The target translational speed calculation unit 551 calculates positional data of each of the outline points RP. The target translational speed calculation unit 551 calculates a relative position between the representative point O of the swing body 2 and each of the plurality of outline points RP, in the vehicle body coordinate system. The target translational speed calculation unit 551 is configured to calculate the relative position between the representative point O of the swing body 2 and each of the plurality of outline points RP of the bucket 8, in the vehicle body coordinate system, on the basis of the working equipment data including the boom length L1, the arm length L2, the bucket length L3, the tilt length L4, the bucket width L5, and the bucket profile data, and the angle data of the working equipment 1 including the boom angle θ1, the arm angle θ2, the bucket angle θ3, the tilt angle θ4, and the rotation angle θ5. Setting of the outline points RP to the blade edge 9 makes it possible for the target translational speed calculation unit 551 to calculate a relative position between the representative point O and the blade edge 9. The design surface is defined in the vehicle body coordinate system. Therefore, the target translational speed calculation unit 551 is allowed to calculate the distance between the blade edge 9 and the design surface in the vehicle body coordinate system. In addition, the target translational speed calculation unit 551 calculates a position of each of the plurality of outline points RP in the global coordinate system. The target translational speed calculation unit 551 is allowed to calculate the position of each outline point RP of the bucket 8 in the global coordinate system, on the basis of an absolute position of the representative point O of the swing body 2 and the relative position between the representative point O of the swing body 2 and the position of the outline point RP of the bucket 8.
The speed limit calculation unit 551B determines a speed limit of the boom 6 in the direction normal to the design surface by using a speed limit table indicating a relationship of the distance between the bucket 8 and the design surface to a speed limit of the working equipment 1.
The rotation speed ω′target is calculated on the basis of formulas (7) to (10).
ΔTtarget is a parameter corresponding to a time required to correct the attitude of the bucket 8. The P control unit 552D calculates the target rotation speed ωtarget by performing P control on the basis of the rotation speed ωtarget calculated by the rotation speed calculation unit 552C.
The model predictive control unit 56 calculates a predicted value of the amount of control of the working equipment 1 and a predicted value of the amount of control of the swing body 2, on the basis of the target value of the amount of control of the working equipment 1 and the target value of the amount of control of the swing body 2 which are generated by the target value generation unit 55 and a prediction model for the excavator 100 including the working equipment 1 and the swing body 2. The model predictive control unit 56 calculates an amount of drive to control the working equipment 1 and the swing body 2, on the basis of the predicted value of the amount of control of the working equipment 1, the predicted value of the amount of control of the swing body 2, and the design surface so as to maintain the distance between the predetermined portion of the bucket 8 and the design surface and the attitude of the bucket 8, that is, so that the bucket 8 follows a target design surface in a predetermined attitude. The predicted value of the amount of control of the working equipment 1 also includes a predicted value relating to the tilt rotation of the bucket 8.
A prediction model storage unit 561 stores the prediction model for the excavator 100 including the working equipment 1 and the swing body 2. The prediction model includes a dynamic model for the excavator 100. The prediction model includes a model for the swing body 2 that swings about the swing axis RX, a model for the boom 6 that rotates about the boom axis AX1, a model for the arm 7 that rotates about the arm axis AX2, and a model for the bucket 8 that rotates about the bucket axis AX3, tilt axis AX4, and rotating axis AX5.
The prediction model is represented by a discrete state equation and an output equation. The state equation of the prediction model for control of the excavator 100, discretized with a sampling time ΔT is shown in formula (11). Matrices of the state equation is shown in formulas (12) and (13). The output equation of the prediction model is shown in formula (14).
Each of M ∈ R5×5 and Co ∈ R5 is an inertial matrix of an equation of motion and a Coriolis force/gravity vector. Ctay ∈ R2nP is a constant term when a tailor expansion of n·p around an angle θ is performed at the predetermined time t. np is the number of design surfaces to be considered. Outputs from the output equation of the prediction model are an angle θ, an angular speed, the target translational speed vtarget, the target rotation speed ωtarget, the distance d between the blade edge 9 and the design surface, and the flow rate Q of the hydraulic oil.
A prediction unit 562 performs an optimization operation on the basis of the prediction model and calculates the predicted value of the amount of control of each of the working equipment 1 and the swing body 2. As described above, in the present embodiment, the amount of control of the working equipment 1 includes one or both of the moving speed of the bucket 8 and the position of the predetermined portion of the bucket 8. The predetermined portion of the bucket 8 includes the blade edge 9. Furthermore, the amount of control of the working equipment 1 includes an angular speed of the boom 6, an angular speed of the arm 7, and an angular speed of the bucket 8. The angular speed of the bucket 8 includes an angular speed about the bucket axis AX3, an angular speed about the tilt axis AX4, and an angular speed about the rotating axis AX5.
In the present embodiment, the prediction unit 562 calculates a predicted value of the amount of control relating to the tilt rotation of the bucket 8 on the basis of the target value of the amount of control relating to the tilt rotation of the bucket 8 and the prediction model, and controls the working equipment 1 on the basis of the predicted value. As described above, the amount of control relating to the tilt rotation of the bucket 8 includes at least one of the tilt angle θ4, the tilt speed, and the tilt acceleration of the bucket 8. Furthermore, the prediction unit 562 calculates the predicted value of the amount of control of the swing body 2 on the basis of the target value of the amount of control of the swing body 2 and the prediction model, and controls the swing body 2 on the basis of the predicted value so that the bucket 8 follows the target design surface in the predetermined attitude.
The prediction unit 562 predicts a value on the left side of formula (14) several steps ahead from the present time.
The prediction unit 562 calculates the amount of drive to control the working equipment 1, on the basis of at least one of a predicted value of the moving speed of the bucket 8, a predicted value of the angular speed of each axis, a predicted value of the position of the blade edge 9 of the bucket 8, and a predicted value of the flow rate of the hydraulic oil. The prediction unit 562 calculates the amount of drive so that the predicted value of the amount of control follows the target value thereof. The prediction unit 562 calculates the amount of drive to control the swing body 2 on the basis of a predicted value of the swing speed of the swing body 2.
In the present embodiment, the prediction unit 562 calculates the amount of drive so that the bucket 8 in a predetermined attitude follows a target design surface, on the basis of the predicted value of the moving speed of the bucket 8, the predicted value of the angular speed of each axis, the predicted value of the position of the blade edge 9 of the bucket 8, the predicted value of the flow rate of the hydraulic oil, a predicted value of a swing speed of the swing body 2 and the design surface. In other words, the prediction unit 562 calculates the amount of drive so that the bucket 8 does not dig below the design surface and the position of the blade edge 9 and the position of the design surface coincide with each other.
The prediction unit 562 calculates the amount of drive to control the working equipment 1 and the swing body 2 so that an evaluation function has a minimum value and each constraint condition is satisfied.
In model predictive control, the evaluation function as shown in formula (15) is generally used.
E(t)=Ey(t)±Eu(t)+EΔu(t)+Ec(t) (15)
Ey(t) is a difference between a target value and a predicted value in output, Eu(t) is a difference between a target value and a predicted value in input, EΔu(t) is a magnitude of the change in input, and Ec(t) is a penalty function that is imposed when the constraint conditions which are described later are not satisfied. In the present embodiment, Eu(t)=0 and EΔu(t)=0, and a tracking error in output with respect to the target value in output is used as the evaluation function. The evaluation function is shown in formulas (16) and (17).
r(t+i|t) is a target value of time t+i at the time t, y(t+i|t) is a plant output at the time t+i predicted at the time t, Hp is a prediction horizon that determines how many steps ahead a prediction is made, and W is a diagonal matrix that weights to make a prediction variables.
The constraint condition calculation unit 57 calculates the constraint conditions. The constraints include a first constraint condition relating to the performance of the excavator 100 and a second constraint condition relating to the position of the bucket 8. The prediction unit 562 calculates the amount of drive so as to satisfy the constraint conditions calculated by the constraint condition calculation unit 57.
The excavator 100 as a control target has a limit on the angles θ, angular speed, angular acceleration, and the flow rate of the hydraulic oil, of the working equipment 1. For example, each angle θ through which the working equipment 1 is movable has a limit. Likewise, each of the angular speed and angular acceleration of the working equipment 1 has a limit. Furthermore, the flow rate of the hydraulic oil discharged from the hydraulic pump 17 has a limit. As described above, the excavator 100 has a limit on hardware. Therefore, it is necessary to consider the first constraint condition indicating limitation on the hardware of the excavator 100, also in the model predictive control. The constraint condition calculation unit 57 calculates the first constraint condition including the angles θ, angular speed, angular acceleration, and flow rate of the hydraulic oil, of the working equipment 1. The prediction unit 562 calculates the amount of drive so as to satisfy the first constraint condition.
The constraint conditions of the angles θ, the angular speed, and the flow rate of the hydraulic oil are shown in formulas (18) to (21).
The constraint condition of the angular acceleration is shown in formula (22).
{umlaut over (θ)}min≤{umlaut over (θ)}(t)≤{umlaut over (θ)}max (22)
In the present embodiment, the constraint condition calculation unit 57 converts the constraint condition of the angular acceleration into the constraint condition of torque. The constraint condition of the angular acceleration after conversion is shown in formula (23).
{umlaut over (θ)}min≤M−1(t)τ(t)−M−1(t)Co(t)≤{umlaut over (θ)}max (23)
In control of the working equipment 1, it is necessary to prevent the bucket 8 from digging below the design surface. In other words, the bucket 8 has a limit on position so as not to dig below the design surface. Therefore, it is necessary to consider the second constraint condition indicating limitation on the position of the bucket 8, also in the model predictive control. The constraint condition calculation unit 57 calculates the second constraint condition including the position of the bucket 8 relative to the design surface. The prediction unit 562 calculates the amount of drive so as to satisfy the second constraint condition.
An output d(t) indicates the distance between the blade edge 9 and the design surface. An equation of the i-th design surface is represented by a unit normal vector ni as n·p+di=0. Formulas (24) and (25) show the constraint conditions for preventing the right and left ends of the blade edge 9 from digging below the design surface.
ni·pL(t)≥−di (24)
ni·pR(t)≥−di (25)
The coordinates of the blade edge 9 are non-linear with respect to an angle θ in a state variable. Therefore, a linear approximation is applied as shown in formulas (26) and (27).
The prediction unit 562 uses the evaluation functions shown in formulas (16) and (17) to perform the optimization operation in the model predictive control so as to satisfy the constraint conditions shown in formulas (18) to (27). An optimization problem in the present embodiment is shown in formula (28). For example, quadratic programming (QP) is used for the optimization operation, but another calculation method may be used.
τ(t) is a control input torque for a control plant and is a solution of the optimization operation. Hu is a control horizon that determines how many steps ahead the inputs are to be handled in the optimization problem.
The prediction unit 562 outputs a control command to control the swing speed of the swing body 2, on the basis of the predicted value of the amount of control relating to the tilt rotation of the bucket 8 and a predicted value of the amount of control relating to the swing of the swing body 2 that have been predicted, so that the distance between the blade edge 9 of the bucket 8 and the design surface and the attitude of the bucket 8 are maintained. The swing body 2 swings on the basis of the operation of the operation device 40 by the driver. When the swing speed of the swing body 2 swung on the basis of the operation of the operation device 40 is excessively high, the prediction unit 562 outputs a restriction command to control the valve device 18 on the basis of the predicted value of the swing speed of the swing body 2 so as to reduce the swing speed of the swing body 2. As described above, the valve device 18 is configured to adjust the flow rate of the hydraulic oil supplied to the swing motor 16. The adjustment of the flow rate of the hydraulic oil supplied to the swing motor 16 adjusts the rotation speed of the swing motor 16 and the swing speed of the swing body 2. The output of the restriction command from the prediction unit 562 to the valve device 18 reduces the flow rate of the hydraulic oil supplied to the swing motor 16. The reduction of the flow rate of the hydraulic oil supplied to the swing motor 16 reduces the swing speed of the swing body 2.
Each of
As illustrated in
The swing body 2 swings on the basis of the operation of the operation device 40 by the driver. When the swing body 2 swings by operation of the operation device 40, if the swing body 2 has a high swing speed, the bucket 8 may dig below the design surface, preventing excavation of the construction target into a desired shape. In other words, when the swing speed of the swing body 2 is high, the tilt rotation of the bucket 8 at a maximum speed may cause digging below the design surface by the bucket 8 because of delay in raising one side of the blade edge 9 as illustrated in
In the present embodiment, the prediction unit 562 calculates a predicted value of the tilt speed of the bucket 8 on the basis of the operation data of the operation device 40 operating the working equipment 1 and the swing body 2. Furthermore, the prediction unit 562 calculates the amount of drive to control the operation of the working equipment 1 including the tilt rotation of the bucket 8, on the basis of the operation data of the operation device 40 operating the working equipment 1 and the swing body 2 so that the blade edge 9 follows the design surface IS.
The command unit 58 outputs the control command to control the working equipment 1 and the swing body 2 (valve device 18) on the basis of the amount of drive calculated by the prediction unit 562.
The command unit 58 outputs the control command to control the swing of the swing body 2, on the basis of the amount of drive of the swing body 2 calculated by the prediction unit 562 when a predicted value of the tilt speed of the bucket 8 is equal to or larger than a predetermined defined value. The predicted value of the tilt speed of the bucket 8 is calculated on the basis of the operation data of the operation device 40 operating the swing body 2 so that the distance between the blade edge 9 and the design surface IS and the attitude of the bucket 8 are maintained.
The defined value relating to the tilt speed may be the maximum speed (maximum value) of the tilt speed, or may be a speed lower than the maximum speed. The maximum speed of the tilt speed is a maximum speed of the tilt rotation of the bucket 8 and is known data determined on the basis of, for example, the performance of the tilt cylinder 14 or the performance of the hydraulic pump 17.
In the present embodiment, the command unit 58 outputs the control command to control the swing speed of the swing body 2, on the basis of the amount of drive of the swing body 2 calculated by the prediction unit 562 so as to limit the swing speed of the swing body 2 when the predicted value of the tilt speed of the bucket 8 exceeds the maximum speed. The predicted value of the tilt speed of the bucket 8 is calculated based on the operation data of the operation device 40 operating the swing body 2 so that the blade edge 9 follows the design surface IS.
The bucket 8 is not able to perform tilt rotation beyond the maximum speed of the tilt speed. When the swing body 2 continues to swing at high speed on the basis of the operation of the operation device 40 even though the predicted value of the tilt speed of the bucket 8 exceeds the maximum speed, the bucket 8 may dig below the design surface IS. The prediction unit 562 is configured to calculate the predicted value of the swing speed of the swing body 2, on the basis of the operation data of the operation device 40 swinging the swing body 2. When the tilt speed of the bucket 8 exceeds the maximum speed, the command unit 58 outputs the control command so as to limit the swing speed of the swing body 2. The tilt speed of the bucket 8 is calculated on the basis of the predicted value of the swing speed of the swing body 2 so that the blade edge 9 of the bucket 8 follows the design surface IS. In other words, the command unit 58 outputs the control command to swing the swing body 2 at a swing speed that is calculated so as to prevent the bucket 8 from digging below the design surface IS when the tilt speed of the bucket 8 exceeds the maximum speed. The tilt speed of the bucket 8 is calculated so that the blade edge 9 of the bucket 8 follows the design surface IS. This inhibits the bucket 8 from digging below the design surface IS.
[Control Method]
The design surface acquisition unit 54 acquires the design surface data (step S1).
The positional data acquisition unit 51 acquires the positional data of the swing body 2 from the position calculation device 20, as a current value. Furthermore, the angle data acquisition unit 52 acquires the angle data and angular speed data of the working equipment 1 from the angle detection device 30, as the current values (step S2).
The driver operates the operation device 40 so that the bucket 8 moves from the first position P1 to the second position P2. Furthermore, the driver operates the operation device 40 so that the swing body 2 turns while the bucket 8 is moving from the first position P1 to the second position P2. The operation data acquisition unit 53 acquires the operation data from the operation device 40. The target value generation unit 55 generates the target value of the amount of control of the working equipment 1 and the target value of the amount of control of the swing body 2 on the basis of at least the operation data of the operation device 40 (step S3).
The target value of the amount of control of the working equipment 1 includes a target value of the moving speed of the bucket 8. The target value of the amount of control of the swing body 2 includes a target value of the swing speed of the swing body 2. The target value of the moving speed of the bucket 8 includes the target translational speed vtarget of the bucket 8, described with reference to
The constraint condition calculation unit 57 calculates the first constraint condition relating to the performance of the excavator 100 and the second constraint condition relating to the position of the bucket 8, on the basis of the operation data of the operation device 40, the angle data indicating each angle θ of the working equipment 1 that changes by operating the operation device 40 and the angular speed data indicating the amount of change in each angle θ per unit time, and the design surface data (step S4).
The prediction unit 562 calculates the amount of drive to control the working equipment 1 and the amount of drive to control the swing body 2 so as to satisfy the constraint condition calculated in step S4, on the basis of the target value of the amount of control of the working equipment 1, the target value of the amount of control of the swing body 2, and the prediction model stored in the prediction model storage unit 561 (step S5).
For the amount of drive, an amount of drive of each valve of the valve device 18 is exemplified. For the amount of drive, a torque of the working equipment 1 may be calculated.
The prediction unit 562 calculates the amount of drive of the working equipment 1 and the amount of drive of the swing body 2, for example, to ten steps ahead from the present time.
The prediction unit 562 calculates the predicted value of the amount of control of the working equipment 1 and the predicted value of the amount of control of the swing body 2, on the basis of the amount of drive calculated in step S5 and the current value acquired in step S2 (step S6).
The prediction unit 562 calculates the predicted value of the moving speed of the working equipment 1, the predicted value of the position of the blade edge 9, and the predicted value of the swing speed of the swing body 2, for example, to ten steps ahead from the present time.
The prediction unit 562 determines whether a predicted value of a bucket speed that is calculated so that the blade edge 9 of the bucket 8 follows the design surface IS exceeds a maximum speed, on the basis of the operation data of the operation device 40 operating the working equipment 1 and the swing body 2 (step S7).
In step S7, if it is determined that the predicted value of the bucket speed does not exceed the maximum speed (step S7: No), the prediction unit 562 recalculates the amount of drive so that the predicted value of the amount of control follows the target value thereof (step S5).
The prediction unit 562 recalculates the amount of drive so as to minimize the value of the evaluation function defined by the target value and the current value of the amount of control. The prediction unit 562 recalculates the amount of drive so as to satisfy the first constraint condition and the second constraint condition.
In step S7, if it is determined that the predicted value of the bucket speed exceeds the maximum speed (step S7: Yes), the prediction unit 562 determines whether the evaluation function has the minimum value (step S8).
The speed of the bucket 8 may be the angular speed or angular acceleration of each axis of the working equipment 1 or swing body 2. The maximum speed may be an upper limit. In other words, in step S7, the prediction unit 562 may determine whether a predicted value of the angular acceleration of each axis exceeds the upper limit angular acceleration.
In step S8, if it is determined that the evaluation function does not have the minimum value (step S8: No), the prediction unit 562 recalculates the amount of drive so that the predicted value of the amount of control follows the target value (step S5).
The prediction unit 562 repeats the processing of step S5, step S6, step S7, and step S8, until the evaluation function has the minimum value.
In step S8, if it is determined that the evaluation function has the minimum value (step S8: Yes), the command unit 58 outputs the control command to control the swing body 2 so as to limit the swing speed of the swing body 2, on the basis of the amount of drive to control the swing body 2 which is calculated so that the evaluation function has the minimum value. The command unit 58 outputs the control command to the valve device 18 so as to drive the valve on the basis of the amount of drive calculated by the prediction unit 562. Furthermore, the command unit 58 outputs the control command to control the working equipment 1, on the basis of the amount of drive to control the working equipment 1 which is calculated in step S5 (step S9).
As described above, the amount of drive is calculated from the present time, for example, to ten steps ahead. The command unit 58 outputs, as the control command, the amount of drive having been calculated in the first step, from the amounts of drive having been calculated ten steps ahead.
The swing body 2 swings with the swing speed limited, on the basis of the control command output from the command unit 58. This makes it possible for the bucket 8 moving from the first position P1 to the second position P2 to perform tilt rotation so as to follow the design surface IS, as described with reference to
[Effects]
As described above, according to the present embodiment, moving the working equipment 1 and the swing body 2 by the model predictive control makes it possible for the control device 50 to output the control command to limit the swing speed of the swing body 2 so that the bucket 8 can follow the design surface IS, when the swing speed of the swing body 2 defined by the operation data of the operation device 40 is excessively high. Therefore, the bucket 8 is inhibited from digging below the design surface IS, and the construction target is excavated into the desired shape.
[Computer System]
According to the embodiments described above, the computer system 1000 is configured to acquire the operation data of the operation device 40 operating the working equipment 1 and the swing body 2, to generate the target value of the amount of control relating to the tilt rotation of the bucket 8 and the target value of the amount of control of the swing body on the basis of the operation data obtained upon operation of the working equipment 1 and the swing body 2, and the design surface IS indicating the target shape of the construction target, to calculate the predicted value of the amount of control relating to the tilt rotation of the bucket 8 and the predicted value of the amount of control relating to the swing of the swing body 2 on the basis of the target values and the prediction model for the excavator 100, to calculate the amounts of drive to control the working equipment 1 and the swing body 2 on the basis of the predicted values and the design surface IS so as to maintain the distance from the predetermined portion of the bucket 8 to the design surface IS and the attitude of the bucket 8, and to output the control commands to control the working equipment 1 and the swing body 2 on the basis of the amounts of drive.
In the embodiments described above, the target value generation unit 55 is configured to employ the speed (translational speed and rotation speed) of the bucket 8 to generate the target value for the model predictive control unit 56. The target value generation unit 55 may employ the position and attitude of the bucket 8 to generate the target value for the model predictive control unit 56.
Note that in the embodiments described above, some or all of the functions of the control device 50 may be provided in an external computer system for the excavator 100. For example, the target value generation unit 55 and the model predictive control unit 56 may be provided in the external computer system so that the amount of drive calculated by the external computer system is transmitted to the excavator 100 via a wireless communication system.
Note that, in the embodiments described above, the construction machine 100 is the excavator. The component elements described in the above embodiments are applicable to a construction machine with working equipment, which is different from the excavator.
Note that in the embodiments described above, the swing motor 16 configured to swing the swing body 2 does not need to be the hydraulic motor. The swing motor 16 may be an electric motor configured to be driven by power supply. Furthermore, the working equipment 1 may be operated not by the hydraulic cylinder 10 but by power generated by an electric actuator such as an electric motor.
Number | Date | Country | Kind |
---|---|---|---|
2019-017393 | Feb 2019 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2019/049874 | 12/19/2019 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2020/158235 | 8/6/2020 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
9464406 | Matsuyama et al. | Oct 2016 | B2 |
11015319 | Izumikawa | May 2021 | B2 |
20090083206 | Shigemori | Mar 2009 | A1 |
20140099178 | Nomura | Apr 2014 | A1 |
20160097184 | Matsuyama | Apr 2016 | A1 |
20160244950 | Kami | Aug 2016 | A1 |
20160251835 | Kitajima | Sep 2016 | A1 |
20160321763 | Shike | Nov 2016 | A1 |
20170130428 | Matsuzaki | May 2017 | A1 |
20170342679 | Iwamura | Nov 2017 | A1 |
20180002900 | Morimoto | Jan 2018 | A1 |
20180016768 | Izumikawa | Jan 2018 | A1 |
20180016771 | Izumikawa | Jan 2018 | A1 |
20190017248 | Okada | Jan 2019 | A1 |
20190127946 | Sakai | May 2019 | A1 |
20190292747 | Takehara | Sep 2019 | A1 |
20200011029 | Sano | Jan 2020 | A1 |
Number | Date | Country |
---|---|---|
20170131484 | Nov 2017 | KR |
20180135939 | Dec 2018 | KR |
2014167718 | Oct 2014 | WO |
2016158779 | Oct 2016 | WO |
2018030220 | Feb 2018 | WO |
Number | Date | Country | |
---|---|---|---|
20220106773 A1 | Apr 2022 | US |