Research on Dynamic Stability Optimization of Rock Spreader considering Prevention of Geological Subsidence and Local Collapse in Reclamation Area

The uneven stacking of waste rock materials caused by the dynamic instability of rock spreaders cannot be ignored for geological subsidence and local collapse in reclamation areas. Based on the dynamic stability of the existing rock spreader, combined with the actual application conditions, the load conditions of the rock spreader are analyzed. The static and dynamic analysis and optimization of the key structures such as the receiving arm, the discharging arm, and the main tower are carried out by using the topology optimization method. The optimized virtual prototype of the whole machine is established to verify its rationality. The research results show that the total weight of key parts such as receiving boom, discharging boom, and the main tower is reduced by about 4.5% (4650kg), the lower-order vibration frequencies below the ﬁfth order of those parts are not within the resonance range, and the stability of the whole machine is greatly improved, which fundamentally solves the problem of uneven rock stacking caused by impact stability. And, the maximum stress and deﬂection of the whole machine under the worst working conditions meet the actual requirements. The research ideas and conclusions can provide solutions to prevent the subsidence or local collapse of reclamation areas caused by uneven stacking materials.


Introduction
A self-moving rock spreader is a kind of equipment that can orderly discharge the rock and other materials except mineral resources to the fixed dump in the process of open-pit mining, which provides the basic conditions for the reclamation of the mining area after mining [1,2]. e development of open-pit mining toward high-speed, large capacity, and automation puts forward not only higher requirements for mining technology but also higher requirements for operating machinery and equipment [3,4]. In the early stage, the technology of loading and unloading waste materials by truck and long-distance transportation between mining area and dump has many problems, such as many potential safety hazards, low degree of automation, high emission of waste gas pollution, large consumption of gasoline, and chaotic dump. With the development of mining technology, most open-pit mines combine the rock spreader and belt conveyor to form the rock exhausting system of open-pit mines [5,6]. e combination of the two has the advantages of long conveying distance, high conveying efficiency, and good stability [7,8]. erefore, the self-moving rock spreader has become one of the indispensable types of equipment for continuous or semicontinuous high-efficiency mining technology of open-pit mines.
At present, under the action of continuous rainfall, weathering, and human disturbance, the phenomenon of uneven settlement or local collapse often occurs in the reclamation area of waste rock and soil [9][10][11]. However, the uneven stacking of rock materials by rock spreaders will make this phenomenon more serious and even lead to geological disasters and pollution [12][13][14].
ere are two main reasons for uneven rock materials piled up by rock spreaders. e first is the unreasonable position of lining and line technology. e second is the uneven rock stacking caused by the instability of the rock spreader. With the rapid development of intelligent control, accurate control has been achieved in the layout position and line technology [15]. However, the problem of dynamic stability of rock spreaders is easy to be ignored. erefore, it is urgent to prevent the subsidence or local collapse of reclamation area caused by uneven stacking materials due to the instability of the rock exhauster structure. Although the research on the rock spreader has achieved phased results, the design of the rock spreader mostly adopts the way of analogy design, and there is a technical bottleneck in the structure refinement design of the rock spreader [16,17]. Yang and Pu [18] optimized the limit device of the PL-2400 rock spreader support truck, which improved the reliability of the rock spreader. Dong and Wang [19] proposed the stacking and discharging process design of rock spreaders based on the scale parameters of waste dump and equipment parameters of rock spreader. Han [20] and Meng et al. [21] optimized the cab and control interface of the rock spreader according to the problems existing in the working environment and control room of the rock spreader. Hu et al. [22] and Zhang [23] have made a new design for the lubrication system of the rock spreader, which improves the working efficiency and service life of the rock spreader. Yao et al. [24] put forward a method to limit the vibration of the discharging boom of the rock spreader. erefore, on the basis of the existing research results of the rock exhauster, the formation of a more perfect calculation method for the structure optimization of the rock spreader can not only shorten the design cycle of the rock spreader and promote the standardization and serialization of the rock spreader but also improve the uniformity of the material of the rock spreader and prevent the serious subsidence or landslide in the reclamation area to a certain extent.
Structural topology optimization is an effective lightweight design method to solve complex working conditions [25][26][27]. In the field of mine engineering, relevant scholars have studied different types of structural parts by using topology optimization theory and achieved good results [28,29]. Wu et al. [30] applied the variable density method to optimize the topology of the main arm and the small arm of the hydraulic chisel loader. e result shows that the weight is reduced by 8%, but the stress is reduced and the safety factor is improved. Tang et al. [31] established the topology optimization model of the key structural parts of the feeder by using the variable density method and the substructure method and carried out the optimization. After the optimization, the natural frequency of the structure was increased by 4.17% and 6.30%, respectively, and the maximum dynamic stress was reduced by 46.4%, which improved the overall stiffness and structural reliability of the feeder. Pang et al. [32] used the topology optimization method to study the lightweight of the track frame of the WK-75 excavator based on multibody dynamics. After topology optimization, the thickness of the main support plate and rib plate of the track frame is reduced, and the material utilization efficiency is improved.
To sum up, relevant scholars have made some progress in structural stability and topology optimization of rock spreaders [33], but they have not formed a complete optimization design method and theory for the dynamic stability of key structures of rock spreaders. According to the field investigation and theoretical analysis, it is easy to produce large impact vibration due to the large self-weight of the key structure of the rock spreader and its low-order frequency is close to the natural frequency. is not only leads to the short service life of the rotary bearing of the receiving boom and the counterweight boom but also causes the failure of the hydraulic cylinder connecting the discharging boom from time to time [34]. More importantly, it will cause the uniform discharge of geotechnical materials by the rock spreader. In view of the above problems, firstly, the actual load conditions of rock spreaders are analyzed, based on this, a static and dynamic analysis of the existing structure was made, and the problems existing in its dynamic stability were discussed. en, the topology optimization analysis of the three main structural parts of the discharge arm, the receiving arm, and the main tower is carried out by using the variable density topology theory. And then, based on the topology results and combined with the field application conditions, the new design of each structure is carried out. Finally, the static and dynamic verification of the key structural parts after topology optimization is carried out, and the whole machine joint simulation is carried out to verify the correctness of the topology optimization results. Here, it is aimed to explore a set of topology optimization design theories and methods suitable for the key main structural parts of a large rock spreader, so as to improve the dynamic stability of its main structural parts, and make the stress on the structural parts uniform and give full play to the mechanical efficiency of materials. e research conclusion can provide an important guarantee for solving the discharge uniformity of rock spreaders and indirectly curb or prevent the serious settlement or local landslide in the reclamation area of waste rock and soil.

Analysis of Load Characteristics of
Rock Spreader e rock spreader is mainly composed of receiving system, discharging system, pitching system, counterweight system, main tower, and crawler walking device. e schematic diagram of the whole machine is shown in Figure 1. Taking a certain type of rock spreader as the research object, the main design parameters of this type of rock spreader are shown in Table 1.
In the process of discharging rock, the rock spreader will be subjected to the complex load, which is accompanied by the process of full load starting and braking, sudden unloading, and climbing with a load. In the calculation process, the load value is determined according to the comprehensive consideration of the material conveying capacity and the external environmental load [35,36]. According to the code for the design of steel structure of mobile continuous bulk material handling equipment, the load of rock spreader is quantified. en, it is designed and checked according to the following working conditions [37,38]: 2 Advances in Civil Engineering Condition I: the rock spreader is only subjected to selfweight load (no-load working condition) Condition II: normal operation of rock spreader (rated load with wind load) Condition III: abnormal operation (rated load with wind load and seismic load) To alleviate the impact of the rock spreader on the receiving boom in the working process, the left end of the receiving boom is connected with the receiving boom support trolley through the guide wheel group and guide rail. e length of the guide wheel group is 2500 mm, and the length of the guide rail at the left end of the receiving boom is 7500 mm. In the working process, the position of the guide wheel group determines the three typical working positions of the receiving boom. e three positions are as follows: the guide wheel of the support car of the unloading car is located in the left limit position, the middle position, and the right limit position of the guide rail. ere are three typical working positions for the discharging boom of rock discharging machine: the discharging boom is horizontal, the discharging boom is upward 13°, and the discharging boom is downward 4°.
Influence analysis of each part: ① e receiving boom of the rock spreader supports the whole receiving conveyor to complete the receiving operation. e left end of the receiving boom is connected with the receiving boom support trolley through the guide wheel guide rail, and the right end is suspended on the counterweight boom through the slewing bearing. e weight of the receiving boom directly affects the service life of the slewing bearing. ② e discharging boom supports the whole discharging conveyor to complete the discharging operation. e discharging boom mainly realizes the pitching movement of the discharging boom through the hydraulic cylinder. e quality of the discharging boom directly affects the service life of the hydraulic cylinder and other parts. ③ e main tower is used to fix the discharging boom and the counterweight boom. e above three key structural parts occupy a very important position in the upper structure of the rock spreader, and their structural rationality directly affects the working flexibility and energy consumption of the whole machine. e load and working condition combination of receiving boom are shown in Table 2. e load and working condition combination of the discharging boom are shown in Table 3.

Static and Dynamic Characteristics
Analysis of Key Structure of Rock Spreader  [39]. e dynamic equation of the structure is as follows:

Advances in Civil Engineering
For the static analysis of linear structure, the influence of time is ignored. e following relationship can be deduced by combining (1): where [F] is the static load without considering the time effect. e stress and strain values of the structure under load can be obtained by solving equation (2). According to the deformation value and stress value, we can check whether the structure meets the user requirements under different working conditions. e structural strain can be expressed as where ε { } is the strain matrix, f is the concentrated load matrix, [N] is the shape function matrix, δ { } e is the node displacement matrix, and [B] is the constant strain element matrix. e structural stress can be expressed as where σ { } is the stress matrix, [D] is the elastic matrix, ε { } is the strain matrix, ε { } e is the node displacement matrix, and [S] is the constant stress element matrix. Q345E is mainly used as the boom material of the rock spreader. In view of the complexity of the loads, the third strength theory or the fourth strength theory is needed to check the strength of the rock spreader structure. Because the fourth strength theory mainly considers the influence of the principal stress σ 2 , it is closer to reality than the third strength theory. erefore, the static strength of the rock spreader structure is checked according to the stress value σ e calculated by the fourth strength theory. e von Mises stress can be expressed as where σ 1 , σ 2 , and σ 3 are the three principal stresses of any particle, and σ 1 > σ 2 > σ 3 . e strength check standard is as follows: where σ e is the maximum equivalent stress and [σ] is the allowable stress. e structural stiffness can be expressed as where [K] e is the element stiffness matrix and V is the unit tetrahedral volume.  Advances in Civil Engineering e allowable tensile stress of Q345E used in the boom structure of the rock spreader is [σ] � 345 MPa, but the allowable safety factor is clearly different under different working conditions. According to the "code for the design of steel structure of mobile continuous bulk handling equipment," different safety factors should be selected for the structure of rock spreader under different working conditions [36], and the allowable stress table is shown in Table 4.
According to "the code for design of steel structures," the allowable deflection value of the boom of rock spreader is L/250 [39]. e length of the receiving boom is 50 m, so the allowable deflection is 250 mm. e length of discharging boom is 21.6 m, so the allowable deflection is 86.4 mm. e height of the main tower is 20 m, so the allowable deflection is 80 mm.

Establishment of the 3D Solid Model of the Key
Structure.
e assignment of structural material properties is as follows: elastic modulus E � 210 GPa, Poisson's ratio μ � 3, density ρ � 7850 kg/m 3 , and yield strength σ � 345 MPa. To improve the efficiency of simulation, on the basis of not affecting the accuracy of the key object, some secondary structural parts are simplified as follows [40]: (1) Some features such as threaded holes, fillets, and chamfers are simplified (2) e receiving conveyor, walking platform railing, driving device, and reversing device are added to the relevant boom in the form of load (3) e influence of the connecting plate on the structure is ignored (4) e weld and bolt connection are ignored (5) e wire rope is established by using the line body in the submodule "Geometry" of "ANSYS Mechanical" e models of receiving boom, discharging boom, and main tower are shown in Figure 2.

Finite Element Solution and Analysis of Key Parts
(1) Finite Element Analysis of Receiving Boom. e guide rail at the left end of the receiving boom is supported by the guide wheel of the supporting trolley, so the translational degrees of freedom along the direction of the guide rail are released, and the degrees of freedom in other directions are limited. e right end lug is hinged with the fixed end lug of the unloading boom, so the degree of freedom of rotation around the hanging shaft is released at the hinge, and the degree of freedom in other directions is limited.
As shown in Figure 3, under working condition III, the maximum deformation of the receiving boom is 66.74 mm, which occurs in the center of the receiving boom span; the maximum stress is 190.67 MPa, which occurs on the vertical steel plate at the left end of the receiving boom.
(2) Finite Element Analysis of Discharging Boom. According to the matching condition of the discharging boom and the rotary platform, its abutments can be set as a hinge joint; that is, it can limit the movement degree of freedom in XYZ direction and the rotation degree of freedom around X-and Z-axis and release the rotation degree of freedom around Yaxis. e connection between the strut and the hydraulic cylinder is hinged; that is, the movement degree of freedom in XYZ direction and the rotation degree of freedom around the X-and Z-axis are limited, and the rotation degree of freedom around the Y-axis is released.
When adding a load, in order to better simulate the load of walking platform railings and conveyors, the model is processed by the design modeler in ANSYS Workbench, and the imprint surface is added at the position where the load needs to be applied in the model, so as to apply the load more accurately [41].
It can be seen from Figure 4 that the total deformation displacement of the discharging boom under condition III is the largest, which occurs between the second and third hanging points, and the maximum deformation is 48.259 mm. Under condition III, the maximum stress of the discharging boom is 168.77 MPa, which occurs at the connection between the diagonal web member and the main chord between the two groups of suspension points, and the stress distribution in most other areas is about 90 MPa or less.
(3) Finite Element Analysis of Main Tower. It can be seen from Figure 5 that the maximum deformation of the tower under adverse conditions is 56.685 mm, and the maximum deformation occurs at the top of the main tower. e maximum stress value of the main tower is 218.06 MPa, and the maximum stress position occurs at the connection between the counterweight boom and the main tower. Because there is no stiffener at the junction of the original tower, the structure is relatively weak.
To sum up, by comparing the allowable values of stress in Table 4 and the allowable values of deflection and deformation mentioned above and analyzing Figures 3-5, it can be seen that the maximum stress value and maximum deflection and deformation value of key parts such as receiving boom, discharging boom, and main tower are within the safe range, but they are all abundant and have large optimization space. In addition, the distribution of stress and deflection in different parts is uneven, and the stress and deflection in many parts are very small, so it is difficult to play the maximum utility of materials.

Modal Analysis of Boom.
e receiving boom and discharging boom are equipped with belt conveyors, driving e periodic work of these parts will inevitably bring the excitation sources to the boom and affect the dynamic performance of the boom. e boom will produce complex vibration and deformation under the excitation of idler and motor, so the research on the dynamic characteristics of boom has a great impact on the safety and service life of boom. As the basis of dynamic analysis, modal analysis is a common method of structural dynamic characteristics analysis. Modal analysis can be used to solve the natural frequency and main mode of vibration of the boom of the rock spreader [42], and the frequency and resonance position of the boom of the rock spreader can be obtained.
According to the analysis results, the optimal design of the rock spreader can be carried out, which provides a theoretical basis for improving the dynamic characteristics of the boom.
When the speed of the driving motor is 1485 r/min, the excitation frequency of the driving motor is f � 1485/ 60 � 24.75 Hz. e roller rotation excitation should be calculated according to the belt speed. Under the rated working condition of the rock spreader, the belt speed v � 4 m/s, so the excitation frequency f � 1000v/πd � 4000/159π � 8 Hz. If the ±15% resonance band is selected [43], the resonance band caused by the idler is in the range of 6.8 Hz ∼ 9.2 Hz.   Figure 6. e frequency of the receiving boom is far from that of the motor, so it will not be affected by the vibration frequency of the motor. e fifth-order natural frequency of the receiving boom is 7.60 Hz, which is easy to cause resonance in the resonance band of the idler.
According to the vibration mode diagram of the first six modes in Figure 6 and the vibration mode description in Table 5, it can be seen that the stiffness of the head of the receiving boom is weak, and the multistage vibration modes show the maximum vibration mode at the head of the receiving boom. e first natural frequency of the receiving boom is low, so it is necessary to optimize the structure of the receiving boom, so as to improve the low natural frequency of the receiving boom, improve its dynamic performance, and reduce the risk of resonance.

Modal Analysis of Discharging Boom.
In the same way, constraints are imposed on the discharging boom and modal analysis is carried out to obtain the first six modes. e calculation results are shown in Figure 7.
It can be seen from the above Figure 7 and Table 6 that the natural frequencies of each order of the discharging boom are far different from the excitation frequency of the motor, so the motor excitation will not cause resonance. However, the fifth-and sixth-order natural frequencies of the discharging boom are 8.36 Hz and 8.93 Hz, respectively, which are all within the resonance band of the idler and are easy to cause resonance.

Topology Optimization Based on the Variable Density
Method [44,45].
e topology optimization module of ANSYS Workbench uses the variable density method to optimize the structure topology. Firstly, the continuum structure is discretized by the finite element method, and the relationship between relative density and elastic modulus is expressed in the form of the density interpolation function of continuous variables. en, taking the relative density of the elements as the optimization variable, each element corresponds to an optimization variable in the optimization process. By changing the value of each optimization variable, the elastic modulus of each element and the stiffness matrix of the structure are changed, so that the material is redistributed in the matrix model, and finally, the optimal structure meeting the stress, deformation, and other constraints is obtained. [44][45][46]. e Lagrange multiplier method is an important method to solve the extremum problem and conditional stationary value problem. Its central idea is to transform the extremum problem of original function under constraint condition into the extremum problem of the functional under unconstrained conditions, that is, to complete the transformation from a constrained optimization problem to an unconstrained optimization problem.  Advances in Civil Engineering introduction of the Lagrange multiplier method simplifies the mathematical calculation process, the Lagrange multiplier method is selected to solve the variable density topology optimization problems.

Mathematical Algorithm for Topology Optimization
It is assumed that the density of each virtual cell in the topology optimization space is ρ i , and the actual density of the optimized space is ρ 0 . e relationship between the virtual density ρ i and actual density ρ 0 can be expressed as follows: where x i is the relative density of elements and i is the total number of all elements in the optimization space. en, when the relative density x i � 1, it means that the entity is here; when the relative density x i � 0, it means that the hole is hollowed out. e relative density x i is defined as the threshold value in the software. Setting different threshold values, the topological results are different. Of course, it is necessary to ensure that 0 ≤ x i ≤ 1. To make the relative density x i converge to 0 or 1, a penalty factor p is introduced. K i is used to represent the relative stiffness of the i-th element, and K 0 represents the inherent stiffness of the topological space. e relationship between relative stiffness and inherent stiffness is as follows: According to Hooke's Law, where F is the load, K is the stiffness matrix, and U is the displacement. e flexibility and stiffness are reciprocal to each other. If the flexibility is C, that is, Substituting equation (11) into equation (10), we can obtain the following results: e objective of topology optimization is to minimize the flexibility of the structure, that is, to maximize the material stiffness. e volume ratio before and after optimization is taken as the constraint condition. e relative density of each unit in the topology optimization space is taken as the design variable. On the basis of the above, the mathematical model of topology optimization can be obtained as follows: where k is the set ratio of surplus material, which is the decisive parameter to be set before topology optimization, and v 1 is the total volume before topology optimization; x min and x max are the minimum value and the maximum value of the relative density value x i of the element, respectively. Limiting the value range of x i can ensure that the structure will not produce a singular total stiffness matrix.
Based on the mathematical model of topology optimization and Lagrange multiplier method, the constructed function can be expressed as follows: When the function takes the extremum, zL zx i � 0, i � 1, 2, . . . , n.
By deriving x i from both sides of equation (14) and simplifying it, the following equation can be obtained: By deriving x i from both sides of equation (10), the following equation can be obtained: Taking equation (17) into equation (16), the following equation can be obtained:  (9) and (18), when function (14) takes the extreme value, that is, zL/zx i � 0, the following can be obtained: Taking equation (9) into equation (19) and multiplying x m i at both ends of the equation, the following is obtained: Let ε i � (x i ) m− 1 U T i KU i , and substitute it into equation (20): Here, k represents the number of elements x i ≤ 1 and the iterative equation of x i can be expressed as follows: When x i ≤ x min , the relative density exceeds the lower limit, so x z+1 i � x min e above list is the mathematical calculation process of applying the Lagrange multiplier method to solve topology optimization problems based on the variable density method. As long as the volume ratio before and after topology optimization is limited, the virtual density of each unit in the design space can be calculated through a certain number of iterations. If the density of a place is 1, it means that it is reserved; if the density is zero, it means that it is set as a hole; if the density is between the two, it means that part of it is reserved, so the topology optimization problem can be solved. e topology optimization process in ANSYS Workbench is shown in Figure 8.

Topological Matrix Model.
Before topology optimization of structures, it is necessary to establish a reasonable matrix model for topology optimization. In general, the initial model should be larger than or equal to the outer dimension of the original structure. According to the design input of receiving boom, discharging boom, and main tower overall shape and actual working condition, the topology optimization matrix model of key parts as shown in Figure 9 is established by using SolidWorks software.

Topology Optimization Calculation.
Topology optimization is completed on the premise that the strength and stiffness meet the requirements. In ANSYS Workbench, the design variables, objective functions, and constraints are defined, and the value of residual material proportion k is set, which is consistent with the value of residual material proportion k set by the Lagrange multiplier method in the previous theory. e topology optimization module can iterate the relative density according to the input percentage of retained materials and finally remove the region of materials with a relative density of less than 1. After the previous calculation, the percentage of retained material is in the range of 70%∼90%, and the convergence effect is good. Five groups of topology optimization calculations with material removal percentages of 15%, 16%, 17%, 18%, 19%, and 20% are carried out for receiving boom and discharging boom structure. e topology structure diagram is shown in Figures 10 and 11.
Five groups of topology optimization calculation with material removal percentages of 17%, 18%, 19%, 20%, 21%, and 22% are carried out for the main tower structure. e topology structure is shown in Figure 12.
According to the uniformity and rationality of the structure distribution obtained by topology optimization of different k values, and combined with the particularity of field application conditions, k � 0.18 is adopted for the receiving arm, discharging arm, and the main tower as the basis for subsequent reconstruction.

Model Reconstruction after Optimization.
ere are many irregular parts in the model after topology optimization, so we need to follow the principle of manufacturability when reconstructing the model: (1) In the process of repairing the model, we need to consider manufacturability. We need to use the straight bar instead of the bending part in the topology optimization results. On the one hand, the bending web member is not easy to produce, which will increase the production cost; on the other hand, the bent web member may affect the installation of other parts. erefore, under the condition of keeping the whole structure unchanged, the bending web members after topology optimization are changed into straight web members (2) Because not all cell densities are 0 or 1 after topology optimization, there are some unreasonable local positions that need to be modified in the topology optimization results. When repairing the model, it should be noted that the color lighter part is the middle density unit, which is the removable unit (3) To ensure the lightweight of the boom, we try to use the bar to replace the plate, which can also maximize the lightweight of the boom (4) In the reconstruction process, we try to make the members form a triangle connection, because more triangles will make the boom structure more stable and the stiffness of the boom will be improved (5) e original structure of the connection part between the two ends of the boom and the rest of the rock spreader should not be changed as far as possible, so as to ensure the normal assembly of the parts

Receiving boom
Discharging boom Main tower e optimized models of receiving boom, discharging boom, and main tower are shown in Figure 13. e weight comparison of key parts before and after optimization is shown in Table 7.

Model Verification after Optimization
(1) Verification of Receiving Boom after Optimization. Referring to the previous method and the process of calculating the static and dynamic characteristics of the receiving boom, the static simulation and modal analysis of the optimized receiving boom are carried out. e deformation and stress distribution of the receiving boom under optimized condition III are shown in Figure 14; the comparison of deformation and stress before and after optimization is shown in Table 8. e low-order modal shapes of the optimized boom are shown in Figure 15, and the dynamic performance comparison before and after optimization is shown in Table 9.
It can be seen from Figure 14 and Table 8 that, under the worst condition III, the maximum deformation of the receiving boom after optimization is reduced by 10.8% compared with that before optimization. e maximum stress is increased by 6.2%, but it is still in the safe range, and the stress of the optimized boom tends to be uniform, so it is obvious that the new structure is more reasonable. According to Figure 15 and Table 9, after optimization, the multistage amplitude of the receiving boom decreases, and its stiffness is improved to a certain extent. After optimization, all low-order natural frequencies of the boom are improved to a certain extent, especially the 5th-order frequency is changed from 6.78 Hz to 9.32 Hz, which is 38.68% higher than that before optimization, avoiding the resonance band caused by the idler. e improvement of the natural frequency of the receiving boom means the improvement of the overall stiffness of the boom, which makes the operation of the rock spreader more stable.

(2) Verification of Discharging Boom after Optimization.
Referring to the previous method and the process of calculating the static and dynamic characteristics of the discharging boom, the static simulation and modal analysis of the optimized discharging boom are carried out. e deformation and stress distribution of the discharging boom under optimized condition III are shown in Figure 16. e comparison of deformation and stress before and after   Advances in Civil Engineering 13 optimization is shown in Table 9.
e low-order modal shapes of the optimized discharging boom are shown in Figure 17, and the dynamic performance comparison before and after optimization is shown in Table 10.
It can be seen from Figure 16 and Table 10 that under the worst condition III, the maximum deformation of the optimized discharging boom increases by 12.7%, and the maximum stress increases by 16.8%, but they are within the safe range. It is worth noting that the maximum stress and deformation of the new structure are more evenly distributed and the design is more reasonable.
It can be seen from Figure 17 and Table 11 that the fifthorder frequency of the discharging boom before the optimization is within the resonant frequency range of the idler, while the first five order modes after optimization avoid the resonant frequency of the idler, which improves the dynamic characteristics of the discharging boom to a certain extent.
(3) Verification of the Whole Machine Virtual Prototype. To verify whether the structure of the rock spreader meets the overall load requirements after optimization, the optimized receiving boom, discharging boom, and main tower are reconstructed and verified by finite element analysis. e action of the receiving boom is the most frequent during the operation of the rock spreader.
is section mainly focuses on the stress and deformation of the receiving boom of the rock spreader under the rated load, wind load, and seismic load at the upper, horizontal, and lower positions. Finite element analysis is carried out for the whole machine model assembled after topology optimization, in which the material properties, loads, and boundary conditions are set the same as those for single part analysis, and the finite element model is shown in Figure 18. e simulation analysis is carried out for the three main positions of the upper row, the lower row, and the horizontal row of the rock spreader under dangerous working          16 Advances in Civil Engineering are within the safe range and meet the engineering requirements. It can be seen from Table 7 that the weight of the receiving boom, discharging boom, and main tower is reduced by 5%, 4.1%, and 4.5%, respectively. It can be seen that the topology lightweight of each key structural part not only reduces the weight of the whole machine and reduces the production cost and power of the equipment but also reduces the grounding pressure of the rock spreader to a certain extent and improves the safety of the rock spreader.

Conclusion
(1) e actual load characteristic analysis, structural topology optimization and reconstruction, static and dynamic simulation of each optimized part, and static and dynamic performance analysis of the whole machine verify the correctness of the optimization of the key structure of the rock spreader. A set of topology optimization design theories and methods suitable for the key main structural parts of a large rock spreader is formed, which solves the problems of its large vibration and poor discharge uniformity (2) After the topology optimization of the key parts, the following conclusions are obtained in terms of stability. ① e multistage amplitude of the receiving boom decreases, which improves its stiffness to a certain extent. And its low-order natural frequencies are improved, avoiding the resonance band caused by the idler. ② e first five modes of the discharging boom avoid the resonance frequency of the idler, which improves the dynamic characteristics of the boom to a certain extent. ③ Compared with before optimization, the weight of receiving boom, discharging boom, and main tower can be reduced by 5.0%, 4.1%, and 4.5%, respectively, but their stress distributions are more uniform and their structures are more reasonable (3) e optimized key structures of the rock spreader are assembled and the virtual prototype is established. Under the worst working conditions, the maximum deflection deformation and the maximum equivalent stress of the receiving boom, discharging boom, and main tower after optimization are in the safe range, which can meet the working requirements of the rock spreader. e correctness of topology optimization lightweight design is verified.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.