Research on the Prediction Method of Braking Rotation Angle for Remote-Controlled Excavator

To calculate, analyze, and predict the rotation angle during the deceleration and braking process of large remote-controlled excavators, this article established a spatial coordinate system based on a simplified model of a hydraulic excavator’s upper structure. Using the D-H parameter method, a mathematical model of the working device’s center of gravity and its rotational inertia was established. Based on the characteristics of the excavator’s hydraulic system and the relationship between brake torque variations, a prediction model was developed to forecast the stopping position (brake rotary angle) of the excavator’s bucket after braking. Subsequently, the predicted results were validated using simulation and compared with existing experimental data to assess the accuracy of the model. The findings demonstrate that the predictive model exhibited high precision with minimal error. The utilization of this model enabled effective forecasting of the excavator’s braking position changes, providing a theoretical foundation for the intelligent remote control of excavators.


Introduction
Excavators, as one of the most prominent and widely used engineering machinery, play a crucial role in energy exploration, infrastructure construction, and emergency relief [1]. They have also gradually become indispensable equipment in fields such as transportation, energy, and national defense. However, with the increasingly complex working environment and specific working conditions, modern excavators have put forward stricter requirements in terms of personnel safety, construction accuracy, and work efficiency so that traditional excavators can no longer meet the development needs of future society. Therefore, research on intelligent technology [2], energy efficiency improvement technology [3,4], comfort and safety technology [5][6][7], and environmental protection technology for excavators is of great practical significance. Among them, the development of intelligent technology for excavators is of paramount importance, mainly including autonomous control [8,9], automatic navigation [10,11], and automatic scheduling of excavators. For example, through machine learning and other technologies, excavators can recognize and adapt to complex terrain and obstacles [12][13][14], achieving more efficient operations.
Under normal circumstances, excavators work in harsh environments and even in dangerous places [15]. Operators face the problems of high labor intensity and high safety risks [16]. Remote-controlled and unmanned operation modes have become feasible solutions to these problems and have become the main direction of intelligent excavator development [17]. Among them, the remote-controlled operation mode has been widely used [18]. For example, China Railway Engineering Group used remote-controlled excavators to improve construction quality and efficiency in railway tunnel projects, while the ECR235CL excavator of Swedish company Volvo can be equipped with a remote-controlled system to allow operators to work in dangerous locations. In this mode of operation, a corresponding spatial coordinate system was established to determine the initial center of gravity position of each component of the excavator in the static state. The starting state was then iteratively analyzed and calculated to determine the dynamic changes of the center of gravity position and rotational inertia of the excavator's upper structure during the relative rotation of the working device and the rotary platform around the swing axis. The braking rotation angle of the vehicle itself was predicted in combination with the torque changes of the engine and the reducer during the rotation process, providing a theoretical basis for the positioning of the bucket for the remote operators and exploring a new method for analyzing the dynamic operation of construction machinery. It also provides a new direction for the intelligent development of construction machinery. Figure 1 shows a schematic diagram of the excavator's braking rotation angle, where ω represents the angular velocity of the vehicle during the rotation process, and θ represents the braking rotation angle, which is the angle between the current position and the predicted stopping position.

Decomposition and Simplification of Excavator Modules
To analyze the dynamic changes in the center of gravity and rotational inertia of the hydraulic excavator during the rotation process, the first step was to determine the center

Decomposition and Simplification of Excavator Modules
To analyze the dynamic changes in the center of gravity and rotational inertia of the hydraulic excavator during the rotation process, the first step was to determine the center of gravity coordinates of each component involved in the rotation. The commonly used method is to install corresponding sensors on the specific excavator, and then use related algorithms and software for measurement and reading. However, the noise, temperature, and humidity on the construction site can greatly affect the accuracy of the sensors. Therefore, this paper attempted to use a method based on the spatial positions of each component to determine their center of gravity coordinates. The key was to determine a parameter set that can characterize the size of each component of the excavator, to effectively screen and express the data points of the 3D structure. Based on this, mathematical methods can be used to reconstruct its 3D structure and obtain the initial center of gravity position.

Decomposition Module Hinge Point Coordinate Collection
The Volvo company's EC650 medium-sized hydraulic excavator was selected as the research object, and the components involved in the rotary motion were identified, including the bucket, dipper, boom, cockpit, base, bucket cylinder, dipper cylinder, and boom cylinder, as shown in Figure 2. The maximum lifting state of the excavator (where each working device reaches the highest and farthest position) was chosen as the initial state. The three-dimensional structure of the model shown in Figure 1 was converted into effective data points. Firstly, the excavator's upper structure was disassembled, and each module was a homogeneous component made of steel material. Secondly, we defined feature points as points that can accurately or approximately form the external contour of the model by connecting them one by one. The edge points of the model can be selected as the initial feature points, if it is a plane, they can be directly connected, and if it is a curved surface, the corresponding curvature needs to be determined. We collected the initial feature points of the disassembled modules and organized their spatial coordinates, as shown in Table 1.  The three-dimensional structure of the model shown in Figure 1 was converted into effective data points. Firstly, the excavator's upper structure was disassembled, and each module was a homogeneous component made of steel material. Secondly, we defined feature points as points that can accurately or approximately form the external contour of the model by connecting them one by one. The edge points of the model can be selected as the initial feature points, if it is a plane, they can be directly connected, and if it is a curved surface, the corresponding curvature needs to be determined. We collected the initial feature points of the disassembled modules and organized their spatial coordinates, as shown in Table 1. The center of gravity coordinates of a three-dimensional model can be directly read using Solidworks. However, the calculation accuracy of the algorithm used by the software was affected by the modeling method of the model. When the geometric shape, mass distribution, or material density of the model changes, manual adjustments and corrections must be made. Obviously, different brands and models of excavators correspond to different models. In practical use, only the center of gravity coordinates of simple objects are calculated through software, and other methods are needed to assist in calculating the center of gravity coordinates of complex and asymmetric objects. In order to calculate the initial center of gravity of the excavator, interpolation based on the initial feature points can be used, that is to predict the values of unknown data points using known data points, thereby expanding the number of effective data points to ensure the accuracy of the calculation result. In theory, the more initial points collected, the better the corresponding interpolation effect. Considering that most of the components in the excavator structure can be represented by straight lines to describe their surfaces or contours, linear interpolation can also achieve good results for data conversion.

Simplification of the Upper Structure and Composite Position Model of Oil Cylinder
For the convenience of subsequent calculations, the excavator's upper structure was projected onto the symmetrical center plane of the working device to be simplified, as shown in Figure 3. For the convenience of subsequent calculations, the excavator's upper structure was projected onto the symmetrical center plane of the working device to be simplified, as shown in Figure 3.  In the figure above, A~F and H~J are the hinge points between various components of the excavator; K represents the tip of the bucket tooth; G a , G b , G c , G BC , G DE , and G H I represent the center of gravity positions of the boom, dipper, bucket, and their respective oil cylinders.
Since the center of gravity position of the cylinder barrel and piston rod of the working cylinder changed with the relative motion of the working device, we needed to simplify its model: the center of gravity position of the cylinder barrel and piston rod was equivalent to a composite center of gravity for calculation. The equivalent model is shown in Figure 4, where 1 is the cylinder barrel and 2 is the piston rod. Since the center of gravity position of the cylinder barrel and piston rod of the working cylinder changed with the relative motion of the working device, we needed to simplify its model: the center of gravity position of the cylinder barrel and piston rod was equivalent to a composite center of gravity for calculation. The equivalent model is shown in Figure 4, where is the cylinder barrel and is the piston rod. The functional relationship of the composite center of gravity position is as follows: (1) • represents the gravity of the cylinder barrel; • represents the gravity of the piston rod; • represents the distance between cylinder and the adjacent hinge point; • represents the distance between piston rod and the adjacent hinge point; • represents the distance between two hinge points of the cylinder; • represents the distance between composite center of gravity and hinge point of the cylinder.
The mass parameters of the cylinder barrel and the piston rod of the bucket, dipper, and boom cylinders are obtained by using Solidworks, as shown in Table 2.  The functional relationship of the composite center of gravity position is as follows: • G 1 represents the gravity of the cylinder barrel; • G 2 represents the gravity of the piston rod; • L 1 represents the distance between cylinder and the adjacent hinge point; • L 2 represents the distance between piston rod and the adjacent hinge point; • L represents the distance between two hinge points of the cylinder; • L c represents the distance between composite center of gravity and hinge point of the cylinder.
The mass parameters of the cylinder barrel and the piston rod of the bucket, dipper, and boom cylinders are obtained by using Solidworks, as shown in Table 2. The distance from the center of gravity of each working device oil cylinder to its lower connecting hinge point in the initial state is: In which, L BC , L DE , L H I represent the composite center of gravity of the oil cylinder of the bucket, dipper, and boom, respectively. For each cylinder, the center of gravity coordinates in the initial state can be easily obtained based on their geometric relationship, which is represented as: in which, x BC , y BC , z BC represent the center of gravity coordinates of the bucket oil cylinders, and the center of gravity coordinates of the dipper oil cylinders and the boom oil cylinders can be calculated similarly. All the initial center-of-gravity coordinates are shown in Table 3. Ignoring the mass of oil entering and leaving the oil cylinder during the actual operation and only considering the dynamic changes of its center of gravity coordinates, a corresponding D-H transformation coordinate system can be established for calculation.

Calculation of the Center of Gravity Coordinates of Excavator
Excavators' upper structure can be abstractly viewed as a large engineering manipulator. The most commonly used method for describing their joint motion is the D-H parameter method, which is often used to describe articulated coordinate systems, but for non-joint motion, conventional transformation matrices are necessary. Therefore, the overall transformation matrix for the center of gravity coordinates can be obtained by multiplying the homogeneous transformation matrices determined by the D-H parameter method. In addition, if the excavator itself is not on level ground, the effect of the gravity component on the entire vehicle must be taken into consideration. In this study, we assumed that the initial static state of the excavator was perpendicular to the ground with no obvious tilt or sway [25] and that work scenarios involving uphill or downhill operations are not within the scope of this paper.

Excavator's Link Coordinate System
Suppose that two adjacent links in a robot are as shown in Figure 5. We denoted the joint closer to the base as the driving joint and the joint closer to the end effector as the driven joint. no obvious tilt or sway [25] and that work scenarios involving uphill or downhill operations are not within the scope of this paper.

Excavator's Link Coordinate System
Suppose that two adjacent links in a robot are as shown in Figure 5. We denoted the joint closer to the base as the driving joint and the joint closer to the end effector as the driven joint.  Generally, the D-H coordinate system is established at the driven joint, which is the side closer to the end effector. That is, the coordinate system {O i−1 } is fixed to Link i−1 and Axis i−1 corresponds to the driving axis of the Link i−1 , Axis i corresponds to the driven axis of the Link i−1 and driving axis of Link i , Axis i+1 corresponds to the driven axis of the Link i .
In Figure 5, the two pairs of red lines are parallel to each other: • θ i -The angle of the x-axis between the coordinate system {O i−1 } and the coordinate system {O i }, which is also the rotational angle of the Axis i ; • d i -The offset of the coordinate system {O i } relative to the coordinate system {O i−1 } along the z i−1 -axis, which can be considered as a translation transformation; • α i -The angle between the driving link and the driven link; • a i -The straight-line length between the driving link and the driven link.
In this context, the rotary coordinate system was fixed to the excavator chassis, and the rotation of the working device relative to their respective pivot points can be considered as the motion of a robot arm with joints. The coordinate system at the boom pivot point A was chosen as the base coordinate system for D-H transformations. Accordingly, the corresponding D-H coordinate systems for the excavator links [26] can be established as shown in Figure 6. Generally, the D-H coordinate system is established at the driven joint, which is the side closer to the end effector. That is, the coordinate system is fixed to and corresponds to the driving axis of the , corresponds to the driven axis of the and driving axis of , corresponds to the driven axis of the . In Figure 5, the two pairs of red lines are parallel to each other: • -The angle of the -axis between the coordinate system and the coordinate system , which is also the rotational angle of the ; • -The offset of the coordinate system relative to the coordinate system along the -axis, which can be considered as a translation transformation; • -The angle between the driving link and the driven link; • -The straight-line length between the driving link and the driven link.
In this context, the rotary coordinate system was fixed to the excavator chassis, and the rotation of the working device relative to their respective pivot points can be considered as the motion of a robot arm with joints. The coordinate system at the boom pivot point A was chosen as the base coordinate system for D-H transformations. Accordingly, the corresponding D-H coordinate systems for the excavator links [26] can be established as shown in Figure 6.  The parameters l AF , l FJ , and l JK are the distances between the pivot points of the boom, dipper, and bucket, respectively (corresponding to the sizes of the D-H parameters).

The Homogeneous Transformation Matrix of the Excavator
The relevant D-H parameters were determined according to the specific dimensions of the excavator's upper structure components (links) shown in Figure 6, as shown in Table 4. Table 4. Parameters of the linkage of the excavator's upper structure.
According to the parameters in Table 4, the D-H base coordinate system A − x 1 y 1 z 1 can be obtained by rotating and translating the rotary coordinate system O − x 0 y 0 z 0 around the z-axis and x-axis, the corresponding rotation matrix is 0 Here, the rotation matrix is determined by Equation (3), and ϕ is the angle parameter corresponding to this rotation method.
The translation matrix is 0 P 1 = 0 0 h 0 T . Therefore, the homogeneous transformation matrix 0 T 1 is: Similarly, the reference coordinate system F − x 2 y 2 z 2 can be obtained by rotating and translating the D-H base coordinate system A − x 1 y 1 z 1 , with the corresponding rotation matrix as 1 R 2 and the translation matrix as 1 P 2 . The homogeneous transformation matrix 1 T 2 is: By analogy, the transformation matrices between different coordinate systems can be obtained as follows: Sensors 2023, 23, 6780 9 of 27 Continuing from the previous statement, the total homogeneous transformation matrix of the center of gravity of each component relative to the revolving coordinate system can be obtained by multiplying the homogeneous transformation matrices of each coordinate system. Therefore, the total transformation matrices of the joint are: The 3 × 1 sub-matrix in the upper-right corner of the total transformation matrix represents the spatial coordinates of the point ( x i y i z i T ). Combining Figures 3 and 6, it can be seen that the dynamic coordinates of each feature point were difficult to obtain relatively accurately through analytical expressions [27]. Therefore, the dynamic center of gravity of the excavator working device during the rotation process can be obtained by homogeneously transforming its initial center of gravity coordinates.

Calculation of Center of Gravity Coordinates for Each Excavator Module
Let x c , y c , and z c be the composite center of gravity coordinate of the excavator working device (excluding the corresponding oil cylinder). It was assumed that the mass of each component is uniformly distributed, so the formula for calculating the composite center of gravity is shown as Equation (11): • m sep -The mass of each decomposed module of the excavator (kg); • tol-Total number of decomposed modules; • x s , y s , z s -Coordinates of each valid data point.
According to Section 2.1.1, 117 characteristic points were collected for the dipper, 83 characteristic points for the boom, 176 characteristic points for the cockpit, and 168 characteristic points for the base in the initial state. As the bucket is composed of irregular surfaces, the use of linear interpolation would harm the modeling process [28]. Therefore, filling the enclosed area surrounding the bucket surface with additional points can be used as new valid data points [29]. The excavator's decomposed modules were created by connecting these initial characteristic points, as shown in Figure 7.
Once the initial characteristic points of the 3D model were determined, they were interpolated using the Griddata Method (a method that fits the scattered data points in vector (x, y, z) to a surface v = f (x, y)), and the interpolation results were checked for how well they matched the original model to allow the center of gravity position to be calculated using the expanded valid data points. The interpolation effect is shown in Figure 8, and it is apparent that the expanded valid data points were sufficient to fit the shape of each component. Once the initial characteristic points of the 3D model were determined, they were interpolated using the Griddata Method (a method that fits the scattered data points in vector to a surface ), and the interpolation results were checked for how well they matched the original model to allow the center of gravity position to be calculated using the expanded valid data points. The interpolation effect is shown in Figure 8, and it is apparent that the expanded valid data points were sufficient to fit the shape of each component. Combining Equation (11), the center of gravity coordinates of each module of the excavator were calculated based on the set of valid data points and were compared with the measured results from Solidworks software, as shown in Table 5 (where blue represents calculated values, orange represents measured values, yellow represents outliers, and Dv represents the difference).  Once the initial characteristic points of the 3D model were determined, they were interpolated using the Griddata Method (a method that fits the scattered data points in vector to a surface ), and the interpolation results were checked for how well they matched the original model to allow the center of gravity position to be calculated using the expanded valid data points. The interpolation effect is shown in Figure 8, and it is apparent that the expanded valid data points were sufficient to fit the shape of each component. Combining Equation (11), the center of gravity coordinates of each module of the excavator were calculated based on the set of valid data points and were compared with the measured results from Solidworks software, as shown in Table 5 (where blue represents calculated values, orange represents measured values, yellow represents outliers, and Dv represents the difference). Combining Equation (11), the center of gravity coordinates of each module of the excavator were calculated based on the set of valid data points and were compared with the measured results from Solidworks software, as shown in Table 5 (where blue represents calculated values, orange represents measured values, yellow represents outliers, and Dv represents the difference).
It can be seen that the calculated values of the center of gravity coordinates of each component of the excavator had some outliers (absolute position error > 30 mm, marked in yellow), which may have been due to the large overall mass of the corresponding component, and the asymmetry in the front, back, left, and right also led to the calculation deviation of the center of gravity coordinates. The algorithm used by Solidworks will make the results more biased toward the symmetrical reference plane of the model, so in actual applications, the characteristics of both methods should be considered to obtain more accurate results. The maximum absolute error between the remaining results and the corresponding software reading values did not exceed 30 mm (excluding the Z-axis, as it is not the sensitive direction of the rotational inertia), which was acceptable relative to the distance of the working device from the vehicle body for several meters. It also had little effect on the final calculation result of the rotational inertia, proving that this method can not only be used to calculate the center of gravity coordinates of engineering vehicles but also overcome various defects in software reading methods and have stronger adaptability. The maximum bucket capacity of the Volvo EC650 excavator was approximately 3.3 m 3 , and its main targets during excavation operations were soil, sand, and finely crushed rocks. The mixed density of these materials fell within the range of 1.35 to 1.45 g/cm 3 . For the purpose of analysis, we assumed a density of 1.4 g/cm 3 for the loaded materials and investigated the variation of the excavator's load center of gravity.
In Solidworks, we loaded the excavator bucket with sand and gravel (density set to the aforementioned 1.4 g/cm 3 ), and measured the variation of the load center of gravity by increasing and decreasing the volume of sand and gravel (ensuring that the total volume remained within the bucket's capacity range), as shown in Figure 9. When the volume of the material reached 3.3 m 3 , the coordinates of the load center of gravity underwent the following changes: • -The center of gravity of the bucket module when loaded with material; When the volume of the material reached 3.3 m 3 , the coordinates of the load center of gravity underwent the following changes: The center of gravity of the bucket module when loaded with material; • X Unload -The center of gravity of the bucket module when unloaded.
The representation method for the Y-axis and Z-axis coordinates followed similarly.
It is important to note the aforementioned center-of-gravity measurement results assumed that the material was evenly distributed inside the bucket (which is the typical working condition in most cases) and not tilted to one side, as shown in Figure 10. When the volume of the material reached 3.3 m 3 , the coordinates of the load center of gravity underwent the following changes: • -The center of gravity of the bucket module when loaded with material; • -The center of gravity of the bucket module when unloaded.
The representation method for the -axis and -axis coordinates followed similarly.
It is important to note the aforementioned center-of-gravity measurement results assumed that the material was evenly distributed inside the bucket (which is the typical working condition in most cases) and not tilted to one side, as shown in Figure 10. When the bucket capacity gradually increased, the center-of-gravity variation along the -axis can be ignored because the even distribution of material did not cause any changes in the center of gravity for a symmetric structure. However, for the variation of When the bucket capacity gradually increased, the center-of-gravity variation along the X-axis can be ignored because the even distribution of material did not cause any changes in the center of gravity for a symmetric structure. However, for the variation of the load center of gravity along the Y-axis and Z-axis, it can be calculated using the following formula: • 4620-Total weight at full load.
Indeed, the numerical changes along the Z-axis will not have an impact on the moment of inertia. In the subsequent calculation process, it is sufficient to consider only the variation along the Y-axis.

The Measurement Method of the Center of Gravity of Whole Excavator
Existing excavator's center-of-gravity measurement techniques mostly focus on measuring the whole vehicle's center of gravity rather than studying the center of gravity position of individual components separately. This may require data collection during the manufacturing of vehicle components. In the measurement method for the whole vehicle center of gravity, the measurement site is limited by the location of the truck scale, and there are issues such as the excavator having a large tilt angle, which can lead to dangerous situations, significant errors in traditional line-marking methods for center of gravity measurement, and difficulties in obtaining repeatable measurement data. Given the limitations mentioned above, this paper proposed a method for measuring the center of gravity of the entire excavator as follows: (1) The whole vehicle weight of the excavator M is measured using a truck scale; (2) The excavator is placed with one side on the knife edge while the other side is suspended with a sling to achieve a horizontal position. A tension detection device is installed on the sling, measure the horizontal distance from the sling to the knife edge, and the value from the tension detection device is recorded. Based on the principle of torque balance, the horizontal distance from the blade edge to the excavator's center of gravity D is calculated and converted into horizontal coordinates.
(3) When suspending the excavator to an inclined position, measure its tilt angle a and record the value from the tension detection device F. Measure the horizontal distance from the sling to the knife edge d. Then, according to the formula: Certainly, the vertical distance from the excavator's center of gravity to the knife edge can be calculated when the excavator is in a horizontal position. This value can then be converted into vertical coordinates.
(4) In step (2), when 'one side of excavator' refers to the rear side, the measurement is the horizontal distance D x , the tension detection device gives the value F x , and the calculated center of gravity coordinates is X = F x · D x /(M · g); when 'one side of excavator' refers to the left side, In the end, corresponding values of D y , F y , and Y = F y · D y /(M · g) will be obtained.
A graphical illustration of the above measurement steps is shown in Figure 11.

Calculation of the Rotational Inertia for Excavator
The attitude of the working device of the excavator constantly changed during the actual working process, such as excavation, lifting, rotation, unloading, and other movements or adjustments, which changed the overall vehicle's rotational inertia [30]. In addition, due to the large distance between the working devices and the rotation center during the working process, coupled with the mass of the rotation platform itself, the excavator as a whole had large rotational inertia. In order to study the influence of the rotational inertia on the braking rotation angle of the excavator, a functional relationship between the rotational inertia relative to the swing axis of the excavator's upper structure and the angle of the boom , the angle of the dipper , and the angle of the bucket was established based on the simplified structure and D-H coordinate system of the upper structure.
(14) Figure 11. Method of measuring center of gravity of excavator.

Calculation of the Rotational Inertia for Excavator
The attitude of the working device of the excavator constantly changed during the actual working process, such as excavation, lifting, rotation, unloading, and other movements or adjustments, which changed the overall vehicle's rotational inertia [30]. In addition, due to the large distance between the working devices and the rotation center during the working process, coupled with the mass of the rotation platform itself, the excavator as a whole had large rotational inertia. In order to study the influence of the rotational inertia on the braking rotation angle of the excavator, a functional relationship between the rotational inertia relative to the swing axis O Z 0 of the excavator's upper structure and the angle of the boom α, the angle of the dipper β, and the angle of the bucket γ was established based on the simplified structure and D-H coordinate system of the upper structure.
The above equation decomposes the rotational inertia of the entire excavator into two parts: the rotational inertia relative to the axis of each component of the working device, which changes with the position of the center of gravity, and the rotational inertia relative to its axis of the rotation platform during operation, which changes less with rotation.

Calculation of Center-of-Gravity Coordinates for Each Excavator Module
Before calculating the rotational inertia of the boom, dipper, and bucket, their working space should be first determined. With the help of Figure 3, a simplified model of the excavator's upper structure was established using Matlab 2019b software, as shown in Figure 12. According to the stroke of each cylinder of the excavator, the range of motion of the boom was determined to be , the range of motion of the dipper was , and the range of motion of the bucket was . The actual rotation range of the entire vehicle was , but the rotational inertia was not affected by the rotation of the vehicle around the slewing axis. To facilitate demonstration, some discrete values were taken as an example. Assuming that the working device of the excavator rotates around the axis each time, we recorded the spatial coordinates of the points that the bucket teeth could reach. There were a total of 270 points, as shown in Figure 13 (The red arrows indicate the direction of rotation of each axis and the black arrow indicates the working radius). According to the stroke of each cylinder of the excavator, the range of motion of the boom was determined to be 60~80 • , the range of motion of the dipper was 20~75 • , and the range of motion of the bucket was 20~150 • . The actual rotation range of the entire vehicle was 360 • , but the rotational inertia was not affected by the rotation of the vehicle around the slewing axis. To facilitate demonstration, some discrete values were taken as an example. Assuming that the working device of the excavator rotates around the axis 5 • each time, we recorded the spatial coordinates of the points that the bucket teeth could reach. There were a total of 270 points, as shown in Figure 13 (The red arrows indicate the direction of rotation of each axis and the black arrow indicates the working radius).
When the excavator performed rotation alone, the distance between each working device and the slewing center remained unchanged, so the rotational inertia of each component did not change. According to the parallel axis theorem for rotational inertia, the rotational inertia of the working device relative to the slewing axis O Z 0 can be calculated as follows [31]: • m i -The mass of each component of the working device (kg); • tol-The distance from the center of gravity of each component to the slewing axis (mm); • n-The number of components and hydraulic cylinders in the working device.
vehicle was , but the rotational inertia was not affected by the rotation of the vehicle around the slewing axis. To facilitate demonstration, some discrete values were taken as an example. Assuming that the working device of the excavator rotates around the axis each time, we recorded the spatial coordinates of the points that the bucket teeth could reach. There were a total of 270 points, as shown in Figure 13 (The red arrows indicate the direction of rotation of each axis and the black arrow indicates the working radius). When the excavator performed rotation alone, the distance between each working device and the slewing center remained unchanged, so the rotational inertia of each component did not change. According to the parallel axis theorem for rotational inertia, the rotational inertia of the working device relative to the slewing axis can be calculated as follows [31]: The variable d i in Equation (15) can be calculated using the following equation. The calculation method for the related distances in the subsequent equations is the same and will not be repeated here.
Here, i T j (1, 4) is the value in the first row and fourth column of the corresponding homogeneous transformation matrix of the component, which represents the value of the center of gravity coordinate x. Similarly, i T j (2, 4) represents the value of the center of gravity coordinate y. In addition, J gi is the rotational inertia of each component relative to its center of gravity. Since this part of the inertia does not change much with the working posture of the component, it can be assumed that the rotational inertia changes linearly from the extreme lifting posture to the extreme lowering posture, and a function for J gi (α, β, γ) can be established as follows: • ∆J i -The total change in rotational inertia of the component from the extreme lifting posture to the extreme lowering posture kg · m 2 ; • ϕ -The rotational angle of each component relative to the hinge point during the working process, where α, β, and γ represent the rotational angles of the boom, dipper, and bucket, respectively ( • ); • ϕ-The range of motion for each component ( • ). It was assumed that the composite center of gravity of each working device cylinder is iteratively obtained from the center-of-gravity coordinates in the initial state. Therefore, the rotational inertia of the bucket, dipper, boom, and their respective cylinders relative to the slewing axis can be calculated according to Tables 2 and 5, and Equations (15) and (16), as shown in Figure 14.
According to Equation (17), the rotational inertia of the bucket, dipper, boom, and their respective oil cylinder relative to their center of gravity was calculated, as shown in Figure 15. According to Equation (17), the rotational inertia of the bucket, dipper, boom, and their respective oil cylinder relative to their center of gravity was calculated, as shown in Figure 15.

Calculation of the Rotational Inertia of the Rotation Platform
Since the composite center of gravity of the rotation platform almost does not change in the rotating space, the rotational inertia of the rotation platform is taken as a fixed value: According to Equation (17), the rotational inertia of the bucket, dipper, boom, and their respective oil cylinder relative to their center of gravity was calculated, as shown in Figure 15.

Calculation of the Rotational Inertia of the Rotation Platform
Since the composite center of gravity of the rotation platform almost does not change in the rotating space, the rotational inertia of the rotation platform is taken as a fixed value:  Figure 15. Rotational inertia of the excavator working device relative to its center of gravity.

Calculation of the Rotational Inertia of the Rotation Platform
Since the composite center of gravity of the rotation platform almost does not change in the rotating space, the rotational inertia of the rotation platform is taken as a fixed value: The term J gj is similar to J gi , which refers to the rotational inertia of each component of the rotation platform relative to its center of gravity (the rotation platform's center of gravity does not coincide with the slewing axis). This value can be taken as a fixed proportion of the excavator's overall rotation inertia. Generally, this coefficient is assumed to be k p = 0.1~0.2.
According to Table 5 and Equation (18), by substituting the mass of the cockpit and the rotation platform as well as their respective distances from the slewing axis, let k p = 0.15, the rotational inertia of the rotation platform can be calculated J pm = 14239 kg · m 2 .

Calculation of the Rotational Inertia of the Excavator
During the actual operation, the mass and mass distribution of each component of the excavator did not change (excluding fuel). The change in the rotational inertia of the whole excavator was caused by changes in the slewing angle of the boom, dipper, and bucket. By combining the previous analyses of the rotational inertia of each component and summing them, the following functional relationship can be obtained: In Equation (19):

R Boom
Cylinder · = (1.53c α + 7.43s α ) 2 + 6.64 m Bucket , m Dipper , m Boom , m Load , m Bucket Cylinder , m Dipper Cylinder , and m Boom Cylinder represent the mass of the bucket, dipper, boom, load, bucket cylinder, dipper cylinder, and boom cylinder, respectively; s α , c α , s β , c β , s γ , and c γ represent the sin α, cos α, sin β, cos β, sin γ, and cos γ respectively. According to Equation (19), the relationship between the rotational inertia of excavator J E and the variations in α, β, and γ can be graphically represented as shown in Figure 16.

Calculation of the Braking Rotation Angle
The braking rotation angle of an excavator refers to the angle when the vehicle runs at the maximum rotation speed and starts braking by cutting off the motor oil supply until the speed reaches zero. Currently, there is no effective method for accurately calculating the exact value of this angle, because it is difficult to accurately measure the excavator's rotational inertia and output torque [32].

Calculation of the Braking Rotation Angle
The braking rotation angle of an excavator refers to the angle when the vehicle runs at the maximum rotation speed and starts braking by cutting off the motor oil supply until the speed reaches zero. Currently, there is no effective method for accurately calculating the exact value of this angle, because it is difficult to accurately measure the excavator's rotational inertia and output torque [32].

Analysis of the Excavator Rotation Motion
As mentioned earlier, the initial braking position of the excavator is at the horizontal plane, and its rotation system is a large inertial system, so its rotary motion characteristics determine the braking angle of the entire vehicle [33]. When the rotation motor and the working oil cylinder perform a combined action, the rotation process can be divided into three stages: the starting acceleration stage, the constant speed slewing stage, and the deceleration braking stage. Among them, the acceleration and deceleration stages have a greater impact on the braking angle. Therefore, considering the torque output of the hydraulic system during these two stages can more accurately predict the excavator's braking rotation angle. The rotary motion characteristic curve of the Volvo EC650 excavator was plotted by collecting angular velocity data [34], as shown in Figure 17. From Figure 17, it can be seen that the segment is the start-up acceleration phase of the excavator's rotary platform, which can be further divided into a variable acceleration phase and a uniform acceleration phase, and they can be simplified as a

Analysis of the Excavator Rotation Motion
As mentioned earlier, the initial braking position of the excavator is at the horizontal plane, and its rotation system is a large inertial system, so its rotary motion characteristics determine the braking angle of the entire vehicle [33]. When the rotation motor and the working oil cylinder perform a combined action, the rotation process can be divided into three stages: the starting acceleration stage, the constant speed slewing stage, and the deceleration braking stage. Among them, the acceleration and deceleration stages have a greater impact on the braking angle. Therefore, considering the torque output of the hydraulic system during these two stages can more accurately predict the excavator's braking rotation angle. The rotary motion characteristic curve of the Volvo EC650 excavator was plotted by collecting angular velocity data [34], as shown in Figure 17.

Calculation of the Braking Rotation Angle
The braking rotation angle of an excavator refers to the angle when the vehicle runs at the maximum rotation speed and starts braking by cutting off the motor oil supply until the speed reaches zero. Currently, there is no effective method for accurately calculating the exact value of this angle, because it is difficult to accurately measure the excavator's rotational inertia and output torque [32].

Analysis of the Excavator Rotation Motion
As mentioned earlier, the initial braking position of the excavator is at the horizontal plane, and its rotation system is a large inertial system, so its rotary motion characteristics determine the braking angle of the entire vehicle [33]. When the rotation motor and the working oil cylinder perform a combined action, the rotation process can be divided into three stages: the starting acceleration stage, the constant speed slewing stage, and the deceleration braking stage. Among them, the acceleration and deceleration stages have a greater impact on the braking angle. Therefore, considering the torque output of the hydraulic system during these two stages can more accurately predict the excavator's braking rotation angle. The rotary motion characteristic curve of the Volvo EC650 excavator was plotted by collecting angular velocity data [34], as shown in Figure 17. From Figure 17, it can be seen that the segment is the start-up acceleration phase of the excavator's rotary platform, which can be further divided into a variable acceleration phase and a uniform acceleration phase, and they can be simplified as a From Figure 17, it can be seen that the 0~t 1 segment is the start-up acceleration phase of the excavator's rotary platform, which can be further divided into a variable acceleration phase and a uniform acceleration phase, and they can be simplified as a uniform acceleration process for simplicity. t 1~t2 segment is the uniform rotation phase, during which the excavator enters this phase when the speed of the rotation platform reaches its maximum, and its speed-time curve approximates a trapezoid. t 2~t3 segment is the deceleration braking phase, during which the speed of the excavator gradually decreases to zero and generates a braking rotation angle. Typically, the state in which the excavator begins to brake is referred to as the "stalled state", meaning that the excavator's kinetic energy is gradually consumed until it comes to a complete stop, which occurs when the brake pedal is pressed or the throttle pedal is released.

Typical characteristics Analysis of Excavator Hydraulic System
The hydraulic excavator performs hydraulic braking through the motor when braking rotation. The hydraulic principle of the rotary motor is illustrated in Figure 18. Assuming the reducer rotates in a certain direction, high-pressure oil is supplied to port A while low-pressure oil returns to port B. When braking, the hydraulic circuit from the hydraulic pump is cut off, i.e., ports A and B are blocked. As a result, the rotary motor operates in a pump mode under the influence of the inertia torque of the turntable. The oil pressure at port A decreases while the pressure at port B increases until it reaches the set pressure of the relief valve B. Consequently, the relief valve B opens, allowing the overflow to flow back to port A. reaches its maximum, and its speed-time curve approximates a trapezoid. segment is the deceleration braking phase, during which the speed of the excavator gradually decreases to zero and generates a braking rotation angle. Typically, the state in which the excavator begins to brake is referred to as the "stalled state", meaning that the excavator's kinetic energy is gradually consumed until it comes to a complete stop, which occurs when the brake pedal is pressed or the throttle pedal is released.

Typical Characteristics Analysis of Excavator Hydraulic System
The hydraulic excavator performs hydraulic braking through the motor when braking rotation. The hydraulic principle of the rotary motor is illustrated in Figure 18. Assuming the reducer rotates in a certain direction, high-pressure oil is supplied to port while low-pressure oil returns to port . When braking, the hydraulic circuit from the hydraulic pump is cut off, i.e., ports and are blocked. As a result, the rotary motor operates in a pump mode under the influence of the inertia torque of the turntable. The oil pressure at port decreases while the pressure at port increases until it reaches the set pressure of the relief valve . Consequently, the relief valve opens, allowing the overflow to flow back to port . Calibrating an accurate model for hydraulic braking equipment involves gathering a substantial amount of data, which can be expensive and time-consuming. Additionally, it is challenging to obtain a pure physical model that suits different behaviors of various braking equipment by precisely modeling each individual component. This is primarily due to the following issues [35]: (1) For fluid systems, hydraulic equipment consists of components that exhibit highly nonlinear behavior, and these components often generate strong discontinuities in the solutions. Moreover, the nonlinear behavior is not only influenced by the values of the considered system states (such as pressure and flow rate) but also by their derivatives. Calibrating an accurate model for hydraulic braking equipment involves gathering a substantial amount of data, which can be expensive and time-consuming. Additionally, it is challenging to obtain a pure physical model that suits different behaviors of various braking equipment by precisely modeling each individual component. This is primarily due to the following issues [35]: (1) For fluid systems, hydraulic equipment consists of components that exhibit highly nonlinear behavior, and these components often generate strong discontinuities in the solutions. Moreover, the nonlinear behavior is not only influenced by the values of the considered system states (such as pressure and flow rate) but also by their derivatives. This type of behavior is particularly common in electromechanical and electro-hydraulic actuators [35].
(2) Internal friction affects the response of the brake in terms of applied torque/force. In practical operation, both friction and more general hysteresis effects are strongly influenced by the travel speed [36].
To reduce computational load, certain parameters involved in the calculations can be simplified. During the braking process, the pressure at port B remains constant at the set pressure of the relief valve B (also known as constant-pressure relief). Referring to Figure 17, during the acceleration phase 0~t 1 , the pressure at port A rapidly increased and then remained almost constant at high pressure (i.e., constant-pressure relief), while the pressure at port B almost reached zero. When the speed is constant in a certain time period t 1~t2 , the pressures at ports A and B are nearly zero. In the braking phase of the time period t 2~t3 , the pressure at port B rapidly rises and then remains at constant high pressure, while the pressure at port A remains at low pressure. By considering the acceleration and deceleration phases 0~t 1 and t 2~t3 as uniform acceleration and uniform deceleration motion, with equal magnitude and opposite direction of acceleration, according to the torque formula M = J × ω (where M is the torque (N · m), J is the rotational inertia kg · m 2 , and ω is the angular acceleration rad/s 2 ), the braking torque is equal in magnitude to the rotational torque, enabling further calculations.

Prediction of the Braking Rotation Angle
When the excavator needs to change its posture during the rotation process, the rotational inertia of the whole vehicle is not fixed. Therefore, the continuous function for calculating the braking angle needs to be discretized, and the time axis is divided into finer intervals for integration [34]. When the excavator does not have any additional actions during the rotation process, the rotational inertia of the whole vehicle is constant. It was previously known that the rotational torque and the braking torque are the same magnitude, then: • Predicting the braking rotation angle by combining Equation (20) and acceleration formula: • θ-The braking rotation angle ( • ); • ω B -The angular acceleration during the braking process rad/s 2 ; • ω max -The maximum rotation speed of excavator (rad/s). In Equation (21), ω max = ω R i s , ω R is the operating rotation speed of the reducer (rad/s). By combining Equations (20) and (21), we can obtain: In general, the transmission efficiency of excavator slewing support ρ s is above 90%, while the transmission ratio from the output gear of the reducer to the slewing support shaft i s may vary depending on the model and manufacturer. Usually, this ratio is 3~6. Therefore, in a specific excavator, the selection of the transmission ratio needs to be determined based on its working load, speed, torque, and other parameters. For the same model (EC650) of the excavator, if ρ s = 90% and i s = 4.5 are given, then λ can be considered a constant. The predicted braking rotation angle θ is positively correlated with the excavator's rotational inertia J E and the square of the working speed of the reducer ω R , and negatively correlated with the working output torque of the reducer.
The output torque of the reducer M R is generally determined according to Equation (23).
• 9550-The torque calculation formula coefficient; • P M -The real-time engine power (kw); • ω M -The engine output speed (r/min); • i R -Gear ratio, which is the ratio of engine output speed to the output speed of the gearbox; • ξ-The gearbox mechanical efficiency. Substituting Equation (23) into Equation (22), we obtain: Before entering the deceleration and braking process, the excavator's engine and gearbox operate stably, allowing the entire vehicle to enter the "idle rotation" state, during which the speed of the gearbox gradually decreases from its initial value to its minimum value. Consult the relevant manual and select the parameters corresponding to the excavator under full load as follows: P M = 393, ω M = 1800, i R = 9600, ξ = 0.83, and ω M = 600~800. However, these parameters will not remain fixed during the rotation process but will fluctuate within a certain range, which will result in a disturbance factor ∆ in Equation (24). The mean value of the ω R was taken as 700 and ∆ conformed to the normal distribution. The excavator's braking rotation angle change relationship under the stationary state was calculated by substituting random numbers, as shown in Figure 19. From Figure 19, it can be seen that although the braking rotation angle fluctuated continuously as the reducer speed decreased, it will gradually decrease to its minimum value. To more scientifically evaluate the braking performance of the excavator, further consideration was given to how to offset most of the effects caused by the difference in rotation speed and torque. According to Equation (25), the theoretical values of the excavator's rotation speed and rotation torque are represented by and , respectively, which are calculated as: Generally speaking, the maximum flow rate of the swing motor for medium-sized excavators is around , while for large excavators, the maximum flow rate may reach or even higher. In addition to the maximum flow rate, the From Figure 19, it can be seen that although the braking rotation angle fluctuated continuously as the reducer speed decreased, it will gradually decrease to its minimum value. To more scientifically evaluate the braking performance of the excavator, further consideration was given to how to offset most of the effects caused by the difference in rotation speed and torque. According to Equation (25), the theoretical values of the excavator's rotation speed ω E and rotation torque M E are represented by ω Et and M Et , respectively, which are calculated as: Generally speaking, the maximum flow rate of the swing motor for medium-sized excavators is around 30~60 (L/min), while for large excavators, the maximum flow rate may reach 100 (L/min) or even higher. In addition to the maximum flow rate, the volumetric efficiency of the swing motor is also a very important indicator, which refers to the ratio of the flow rate to the speed of the swing motor during operation. It affects the rotation speed and output power of the excavator and usually ranges from 0.8~2.0 (mL/r). Furthermore, the mechanical efficiency of the swing motor refers to the ratio of the mechanical power output to the hydraulic power input, which is directly related to the energy consumption of the excavator, and is generally around 80~90%. The excavator's displacement refers to the amount of oil output by the excavator per unit of time, usually expressed by (mL/r), and the size of the displacement determines the torque of the entire vehicle, generally ranging from 10 to 100 (mL/r). Finally, the working pressure of the swing motor refers to the pressure of the hydraulic oil supplied by the hydraulic system to the swing motor, which is measured in MPa, and an appropriate value is around 10~35 (MPa).
Substituting (26), compare θ and θ before and after correction, as shown in Figure 20.

Results
Therefore, selecting , , , , and , the predicted braking rotation angle was corrected according to Equation • -Actual rotation speed of excavator ; • -Actual rotation torque of excavator .
Substituting and into Equation (26), compare and before and after correction, as shown in Figure 20. Based on Figure 20, the comparison between the braking rotary angles before and after the correction showed a numerical difference of (i.e., the corrected angle was different from the angle before correction). However, there is currently no clear reference range for the braking rotation angle of excavators, and the precise values of the vehicle's rotational inertia and output torque cannot be obtained [30]. Therefore, it is necessary to complete the vehicle coordination in actual application scenarios, collect key parameters that reflect its ability status, and use reinforcement learning to make the predicted braking angle more reasonable.
The instrument used to measure the braking rotation angle is shown in Figure 21. Based on Figure 20, the comparison between the braking rotary angles before and after the correction showed a numerical difference of 13.67% (i.e., the corrected angle was 13.67% different from the angle before correction). However, there is currently no clear reference range for the braking rotation angle of excavators, and the precise values of the vehicle's rotational inertia and output torque cannot be obtained [30]. Therefore, it is necessary to complete the vehicle coordination in actual application scenarios, collect key parameters that reflect its ability status, and use reinforcement learning to make the predicted braking angle more reasonable.
The instrument used to measure the braking rotation angle is shown in Figure 21.
was different from the angle before correction). However, there is currently no clear reference range for the braking rotation angle of excavators, and the precise values of the vehicle's rotational inertia and output torque cannot be obtained [30]. Therefore, it is necessary to complete the vehicle coordination in actual application scenarios, collect key parameters that reflect its ability status, and use reinforcement learning to make the predicted braking angle more reasonable.
The instrument used to measure the braking rotation angle is shown in Figure 21. Figure 21. Experimental instrument. Figure 21. Experimental instrument.
The angular velocity sensor (CRS03 gyroscope) was fixed at a certain position on the upper platform of the excavator, and its output signal was connected to the Ncode high-low level layer. The measurement range was ±100( • )/s (turning right means positive, and turning left means negative), and resolution was 0.05( • )/s. When conducting the test for the braking rotary angle, it was necessary to quickly operate the pilot control lever to its maximum stroke, causing the excavator to rotate at its highest speed. When the excavator completed one full rotation and returned to the starting point (this process requires assistance from personnel observing and signaling to the operator), the pilot control lever was swiftly released to its neutral position, allowing the excavator to naturally brake and come to a stop. Then, the data collected by Ncode can be processed through integration to obtain the braking rotary angle. Table 6 shows the data obtained from measuring the braking rotary angle of a batch of the same model (50-ton class hydraulic excavator) from a certain company. In theory, the braking rotary angle for the same model should be the same. However, the measured angle difference rate was as high as 19%. After analysis, the potential reasons for this discrepancy may include variations in the rotary speed, differences in rotary torque, the positioning of the excavator during the experiment, operator's actions, mechanical efficiency, and reading accuracy, all of which can influence the predicted results [37].
It was obvious that the predicted whole vehicle braking rotary angle after correction was closer to the measured value (the average corrected braking rotary angle for 8 excavators in the experiment was 61.401( • ), while the predicted braking rotary angle when the excavator starts braking and the reducer reached maximum speed was 64.288( • )). This demonstrates that the correction method can compensate for some of the differences in rotary speed and torque, validating the scientific and effective nature of the method proposed in this paper. The method is also applicable for comparing prototypes with benchmarks. When manufacturing companies upgrade their products based on market demands, the excavator's rotary speed and torque may undergo changes. In such cases, this method can be utilized to predict the whole vehicle braking rotary angle after the upgrade to evaluate its braking performance. For instance, if it is determined that the braking is too fast, adjustments can be made by increasing the maximum rotary speed or reducing the whole vehicle's rotary torque. Figure 22 shows an another application scenario of the excavator braking rotation angle, where people can determine the distance D between the final position of the bucket and the current position based on the calculated value and the distance between the bucket and the cockpit, providing visual assistance for the remote-controlled operation of construction machinery.

Discussion
The simulation result showed that the excavator had a simplified and clear structure, and the calculation model of the center of gravity and rotational inertia was accurate and reliable. The predicted value of braking rotation angle was in good agreement with the experimental value. In actual application, the prediction frequency can be determined by the gradually reducing speed of the gearbox during the rotation braking process, and the use of multiple dynamic corrections can ensure that the predicted position is consistent with the final stopping position. This result also proved that the prediction of the excavator's braking rotation angle based on the dynamic change of the rotational inertia is feasible. In engineering applications, calibration or learning modes can be introduced to optimize certain parameters selected in the entire method and improve the final prediction accuracy.
The methods and problem-solving strategies employed in this study can also be applied to the prediction of other actions of remote-controlled excavators and related research on the structures of similar large-scale construction machinery. With the acquisition of relevant parameters , , and , the dynamic changes in the rotational inertia during the rotation process of the vehicle at any posture can be analyzed. Combining this with the performance parameters of the gearbox, the final braking rotation angle after the rotation process can be predicted, providing a reliable solution for remotecontrolled operations of construction machinery.
We hope to conduct the following experiments or solve the following problems in the future: (1) In subsequent research, when conditions permit, it would be advantageous to choose specific components of the excavator and measure their actual center of gravity positions. These measured data can then be compared with the results obtained from the simulation method proposed in this paper. It is important to note that the center-ofgravity position only needs to be measured once, and in subsequent calculations, it can be compared with the values obtained from specific software reading values, allowing for the selection of data that best suits the research needs. Once we obtain the relative position

Discussion
The simulation result showed that the excavator had a simplified and clear structure, and the calculation model of the center of gravity and rotational inertia was accurate and reliable. The predicted value of braking rotation angle was in good agreement with the experimental value. In actual application, the prediction frequency can be determined by the gradually reducing speed of the gearbox during the rotation braking process, and the use of multiple dynamic corrections can ensure that the predicted position is consistent with the final stopping position. This result also proved that the prediction of the excavator's braking rotation angle based on the dynamic change of the rotational inertia is feasible. In engineering applications, calibration or learning modes can be introduced to optimize certain parameters selected in the entire method and improve the final prediction accuracy.
The methods and problem-solving strategies employed in this study can also be applied to the prediction of other actions of remote-controlled excavators and related research on the structures of similar large-scale construction machinery. With the acquisition of relevant parameters α, β, and γ, the dynamic changes in the rotational inertia during the rotation process of the vehicle at any posture can be analyzed. Combining this with the performance parameters of the gearbox, the final braking rotation angle after the rotation process can be predicted, providing a reliable solution for remote-controlled operations of construction machinery.
We hope to conduct the following experiments or solve the following problems in the future: (1) In subsequent research, when conditions permit, it would be advantageous to choose specific components of the excavator and measure their actual center of gravity positions. These measured data can then be compared with the results obtained from the simulation method proposed in this paper. It is important to note that the center-of-gravity position only needs to be measured once, and in subsequent calculations, it can be compared with the values obtained from specific software reading values, allowing for the selection of data that best suits the research needs. Once we obtain the relative position of the center of gravity with respect to the corresponding coordinate system for a particular component, we can conduct homogeneous transformation studies to investigate the dynamic changes in the center of gravity.
(2) The excavator's load varies with changes in the bucket capacity, and the increase or decrease in oil in each working device's cylinder can also be considered as part of the overall load of the vehicle. The latter was neglected in this paper. While the former was analyzed to some extent, the calculations were based on full load conditions. In practical applications, it is recommended to employ visual methods to reconstruct a 3D model of the material pile and estimate its volume and mass to improve the accuracy of predictions.
(3) For the braking rotation angle, its influencing factors are concentrated in various stages of the braking rotation process. This involves complex velocity and energy transformations, making accurate modeling and utilization challenging. Combining the model with the hydraulic system could potentially enhance the prediction accuracy and applicability. Many assumptions made in this paper were based on ideal and simplified conditions, and the selection of data requires further discussion and limitations. Once precise control of the vehicle's rotation speed during the rotation process is achieved, all parameters related to speed can be separated from Equation (24). The entire prediction model can then be divided into constant (determined by the mechanical characteristics of the excavator) and variable (rotational inertia, hydraulic system, rotation speed, etc.) parts for separate calculations.
(4) Currently, there are very limited experimental conditions for testing the braking rotation angle of engineering vehicles. The experimental method presented in this paper effectively approximated the theoretical requirements, and the obtained prediction results were also satisfactory. Subsequent research should focus on certain factors such as errors due to human factors, the applicability of this method to different vehicles, and improvements to the correction formula. It should also promote the construction, maintenance, and development of experimental equipment and facilities.