Abstract
A novel rescue accessory with spreading and supporting functions is proposed to promote the performance of usage continuity and improve the rescue efficiency of traditional rescue devices. To optimally design this rescue device, we apply related working principles to analyze the motion of the whole device and its spreading mechanism. In optimization, three important parameters of the spreading width, the jaw rotation angle, and the spreading force are simultaneously chosen as the objective functions. Meanwhile, specific mathematical models are established for estimating the performance of the novel rescue accessory. Then, an improved algorithm based on the nondominated sorting genetic algorithm (INSGA-II) is proposed to solve the above multiobjective optimization problem of the rescue accessory. The INSGA-II optimizes the population diversity, the crossover operator, and the mutation operation in calculation procedures. Furthermore, it can dynamically adjust the crowded distance and enhance the global as well as the local searching capabilities of INSGA-II. Results from computational tests show that INSGA-II has the best performance among other popular algorithms. More importantly, comparison results show that the compromise solution of the Pareto front for the above three objective functions obtained by INSGA-II is superior to those by other algorithms.
1. Introduction
Earthquake is a sudden and destructive natural disaster. Once happened, it usually brings about huge damages and casualties to humankind due to its unpredictability [1, 2]. Therefore, earthquake defenses and rescue operations have always been the focus of governments around the world [3,4]. Researches on the high-performance earthquake rescue equipment become very important.
In recent years, the rescue robot played an important role in the postdisaster rescue status [5, 6]. This is because it can enter the narrow space of a building to carry out rescue actions. Bai et al. [7] designed a transformable wheel-legged mobility rescue robot that can be switched between the wheel and leg functions to adapt to different terrains. Deng et al. [8] proposed a control strategy based on rolling gait and trajectory planning, which might enable the hexapod robot to walk through dynamic environments and complete rescue missions. Zhang and Gao [9] designed a mine rescue robot with a head mechanism of the groundhog-like and studied its mobility and inverse kinematics. In addition, the structure and motion parameters of the hybrid head mechanism were optimized, and the optimal dexterity was obtained. Li et al. [10] designed a tracked rescue robot with an explosion-proof multidrive crawler and a lightweight explosion-proof body support system. Zang et al. [11] designed a cutting robot and a Jack robot to carry out works of cutting obstacles and jacking up heavy debris. However, in general, the rescue robot is too small to cross large obstacles and provide enough applying force. That is why this kind of robot is usually used to carry out detection rather than rescue missions.
On the other hand, the machinery of large-scale rescue equipment and general large-scale construction equipment played an important role in the postdisaster rescue [12, 13]. However, after the earthquake, roads and bridges are often interrupted, and it is difficult to transport the large-scale rescue equipment to the quake-stricken area over a long distance and perform in-time (within the golden 72 hours) rescue. To overcome this problem, lots of general large-scale construction equipment were designed and proposed with features of economic production and simple construction. And, they may reach the disaster area as quickly as possible. However, due to the lack of specialized rescue accessories such as an end-effector, their functions are still limited, and they cannot be effectively used to perform deeper rescue missions [14, 15].
For rescue equipment, some specific rescue accessories with proper functions are essential in performing different specific rescue missions. For example, in the postearthquake rescue scene, when there is a gap in the ruin, we may overcome it by using a spreader which may help to expand the whole rescue equipment to a proper height [16, 17]. Then, the supporting accessory can be used to support or further expand, so as to complete the rescue operation [18, 19]. In short, due to limited space at the rescue scene, it is difficult to accommodate large-scale construction machinery to perform rescue operations. Therefore, spreading and supporting rescue accessories are essential and commonly used in rescue operations.
As regards, in this paper, we intend to design a novel rescue accessory with spreading and supporting functions, which can be quickly docked with general large-scale construction equipment. Then, to better solve the optimal design problem for this accessory, an improved algorithm is proposed based on the nondominated sorting genetic algorithm (INSGA-II) [20] which has been widely used in the field of multiobjective optimization [21–23].
The remainder of the paper is structured as follows. Section 2 describes the integration of spreading and supporting functions, and the working principle of the proposed novel rescue accessory. The kinematics analysis of the spreading mechanism is made, and the multiobjective optimization problem for the accessory is formulated in Section 3. An improved strategy of NSGA-II is proposed, and a comparison of algorithm performance is made in Section 4. Results and discussions are drawn in Section 5. The conclusions of this work are given in Section 6.
2. Novel Design of Rescue Accessory with Spreading and Supporting Functions
2.1. Mechanism
As mentioned above, the spreading and supporting functions are often required in rescue operations. In order to perform these two operations coherently so as to improve the rescue efficiency, we innovatively design a rescue accessory with spreading and supporting functions, as shown in Figure 1. The designed rescue accessory consists of a spreading hydraulic actuator, two supporting hydraulic actuators, two translational hydraulic actuators, a spreading mechanism, and a supporting mechanism. The spreading and supporting forces are provided by a translational hydraulic actuator and two supporting hydraulic actuators, respectively. This rescue accessory may be quickly docked with general large-scale construction machinery, and its hydraulic power is supplied from the construction machinery.

2.2. Working Principle
During rescuing, if the gap in the ruin is small, we may use the spreading mechanism first to expand this gap, as shown in Figure 2(a). Then, we use the supporting mechanism to support or further expand the gap, as shown in Figure 2(b). During the spreading operation, the supporting and translational hydraulic actuators must be retracted to the preset limit positions and then locked. In addition, the supporting actuator must be aligned in parallel with the translational hydraulic actuator so as not to interfere with the spreading operation, as shown in Figure 2(a). When the spreading mechanism expands the gap size to a sufficient extent, we may then adjust and lock the supporting hydraulic actuator at an appropriate angle to perform the supporting operation, as shown in Figure 2(b).

(a)

(b)
3. Kinematics Analysis and Problem Formulation
The spreading mechanism plays an important role in the rescue accessory. Its maximum spreading width and bearing load directly determine the whole rescue efficiency. Therefore, in this section, the kinematics analysis of the spreading mechanism will be made. The spreading width, jaw rotation angle, and spreading force are chosen as the objective functions.
3.1. Kinematics Analysis of Spreading Mechanism
Due to geometric symmetry, only half of the spreading mechanism is shown here in Figure 3, which majorly consists of a four-bar linkage. For clear visualization, the supporting mechanism is hidden. The linkage FD and the axis of symmetry are used as the base to construct a rectangular coordinate system for describing motions. The ith linkage and its length are represented by the vector notation and , respectively. The initial and final positions of the spreading mechanism are shown in black and blue vector lines, respectively. The vector loop closure equations for the initial and final positions are expressed as follows:

Related angles between each link and AD are shown in Figure 4. Then, we may deduce thatwhere j is the imaginary unit.

In addition, (2) can be rewritten based on the Eulerian formulation, as follows:
By solving (3), we can obtainwhere , , , and .
3.2. Mathematical Model of Spreading Width and Jaw Rotational Angle
As known, the efficiency of spreading can be improved when the spreading width becomes large. To enlarge the spreading width, the projection length in the initial and final positions of link 5 on the X-axis should be as long as possible. Here, we investigate three cases of the spreading mechanism with different linkage lengths in the initial and final positions, as shown in Figure 5. For these cases, we can getwhere

For all three cases of different spreading widths, it is observed that the angle between the linkage 5 and the linkage 3 always remains at a fixed value. This means that this angle of the spreading mechanism in the initial position is equal to that in the final position. Therefore, for the first case, we may obtain that
For the second case, we obtain that
For the third case, we obtain that
In total, we have
According to the geometric relationship of the spreading mechanism shown in Figure 5, we can easily derive thatwhere is the spreading width and is the rotational angle of the jaw.
3.3. Mathematical Model of Spreading Force
The force distribution for a spreading operation of the accessory is shown in Figure 6. We assume that the force on the tip of the spreading mechanism is perpendicular to the axis of symmetry during the spreading operation. And, the thrust of the spreading hydraulic actuator is denoted as and the pressure at the point E′ is . The mechanical stability is a necessary condition in spreading operations of rescue missions. We now assume the system is under an equilibrium state during the spreading operation. Thus,

From (12), we obtain
Similarly, we can find that the spreading force is also a function of the length of each linkage. Here, the thrust of the spreading hydraulic actuator is set as .
3.4. Multiobjective Optimization Formulation
As stated above, the spreading mechanism must not only be able to move between the initial and final positions but also have the largest possible spreading width and force. However, a large spreading width is always accompanied by a large rotational angle of the jaw. With the increase of jaw rotational angle, the contact area between the jaw tips and the object gradually decreases, which makes the spreading operation more unstable and is prone to secondary collapse. That is, the spreading width and jaw rotational angle have an opposite relationship. Therefore, our target is to maximize the spreading width and the spreading force and minimize the jaw rotational angle. Thus, three objectives are defined aswhere is the design variable which is defined as
As previously stated, the motion of the four-bar linkage is driven by the spreading hydraulic actuator. Although the length of each linkage is determined at the initial design stage, it should be carefully chosen, not just simply being greater than zero or arbitrarily. In order to enhance the calculation speed and avoid excluding possible optimal solutions, we set the upper and lower bounds of the design variables based on the actual working conditions as
In addition, a four-bar linkage mechanism can be well constructed by properly selecting different design variables. The suggested restrictions for the lengths of linkages are
As shown in Figure 4, when the spreading mechanism moves to the final position, the spreading hydraulic actuator must extend upward, and the linkage BE must expand to the right side. Then, the , , and are restricted as
In many encountering rescue cases, the gap is usually very small, so we expect that the initial spreading width of the spreading mechanism should be as small as possible. Therefore, the is restricted as
In addition, in (4), since the angles and cannot be imaginary numbers, these two angles should satisfy the following requirements:
4. Methodology
On the basis of the so-called “genetic algorithm (GA)” [24, 25], the traditional NSGA-II is a popular method for solving multiobjective optimization problems [20]. However, in practical applications, the NSGA-II has the defects of being prone to premature convergence and having a nonuniform distribution of the Pareto solution set. To improve its defects so as to successfully apply to our study, we propose an improved strategy of NSGA-II which modifies the crowded distance, the crossover operator, and the mutation operation that appeared in traditional NSGA-II. The followings are the details of INSGA-II.
4.1. Improving Crowded Distance
The sorting strategy of NSGA-II is based on the individual crowded distance. The crowded distance is used to estimate the density of solutions around a particular individual in the population. The crowded distance of the ith individual is the sum of the distances between two neighboring individuals in its layer for each object. In the same layer, the individual with a large crowded distance will be preferred. Although this strategy retains individuals within a small crowded distance, it has some limitations. Since the crowded distance is only calculated once, if multiple individuals with a small crowded distance are concentrated in a certain region, these individuals will be deleted at the same time.
To solve the above-mentioned problem, a dynamically adjusting method of the crowded distance is proposed as follows.
Step 1. Calculate the Euclidean distance between two adjacent individuals. Assuming there are n individuals, the crowded distance between the kth and (k + 1)th individuals, , can be expressed aswhere m is the number of objective functions, is the value of the kth individual on the ith object, is the minimum value of the ith objective function, and is the maximum value of the ith objective function.
Step 2. Find the minimum crowded distance and compare with . If is small, then the kth individual is eliminated; otherwise, the (k + 1)th individual is eliminated.
Step 3. Update the crowded distance. Determine whether to meet the number of individuals required. If it is not satisfied, return to Step 2, otherwise, terminate.
4.2. Improving Crossover Operator
Traditional NSGA-II adopts the simulated binary crossover in generating an individual’s offspring. This operation has limited ability in dealing with the optimization problem with complex Pareto fronts, such as slow convergence and low accuracy.
In order to improve the global searching capability of NSGA-II, the normal-distribution crossover is introduced. The normal distribution crossover operator is defined aswhere is the ith crossover variable for the offspring chromosome, and are the ith crossover variables for two paternal chromosomes, and is a normally distributed random variable.
To further enhance the spatial-searching capability of NSGA-II, we introduce a discrete recombination operation to the previous normal distribution crossover operator, as follows:where is a uniformly distributed random number on the interval (0,1).
4.3. Improving Mutation Operation
To improve the local-searching capability of NSGA-II, a cross-mutation algorithm based on the close relative coefficient is proposed. Firstly, we define two variables: the close relative coefficient and the close relative individuals.
Definition 1. The close relative coefficient of two individuals is defined aswhere l is the length of the individual gene, is the kth gene of the individual , and is the weight of the kth allele.
Definition 2. If the close relative coefficient of two individuals is greater than the set threshold, the two individuals are considered to be close relative individuals.
Secondly, the implementation process of this proposed cross-mutation algorithm is as follows.
4.4. Test of INSGA-II
The relevant parameters of these algorithms are shown in Table 1, and the calculation results via these algorithms are shown in Figures 7–14. For both two test functions of ZDT and DTLZ, it is seen from Figures 7–14 that INSGA-II may well converge to the true Pareto front, much better than the other three algorithms. In addition, the distribution of the Pareto solution set obtained by INSGA-II is more uniform than that obtained by the other three algorithms. This test identifies the effectiveness of the proposed INSGA-II.








In order to further verify the conditions of convergence and uniformity for INSGA-II, we now further select the statistical parameters of generational distance (GD) [26] and spread (SP) [27] as the measurement indicators. They can be respectively expressed as follows:where n is the number of individuals in the Pareto front, is the Euclidean distance between the ith objective function vector in the Pareto solution set and the nearest individual in the true Pareto solution set, and is the average of . All four algorithms are now running 30 times independently, and the obtained data are then averaged. The obtained simulation results are given in Table 2.
It is found in Table 2 that the mean values of GD and SP obtained by INSGA-II are smaller than those of the other three algorithms. In other words, INSGA-II has better convergence and uniformity. In addition, when using the test functions of ZDT2 and ZDT3, the standard deviation of GD and SP obtained by INSGA-II is smaller than those of the other three algorithms; that is, this method has better robustness. Specifically, the simulation results were obtained via any of the aforementioned algorithms, adopting ZDT3 as the test function, which usually has a noncontinuous Pareto front distribution. It is seen that the distribution performance of the Pareto front for INSGA-II is better than those for the other three algorithms when using ZDT3 as the test function. In short, the test results indicate that the proposed INSGA-II has significant improvement in calculation convergence and results from consistency when compared with other traditional algorithms.
5. Results and Discussions
5.1. Optimization Analysis with Different Methods
We now perform the optimization analysis for our designed spreading mechanism using INSGA-II, traditional NSGA-II, PESAII, and SPEA2 methods with the related fundamental parameters as shown in Table 1. This multiobjective optimization problem of the spreading mechanism has been already defined in Section 3.4. Through calculations, we obtain the distributions of the Pareto front using the above four methods, as shown in Figure 15.

(a)

(b)

(c)

(d)
Since it is difficult to observe the distribution of the solution set in the three-dimensional space, we now extract and compare the data of extreme values of obtained objective functions, the maximum , the maximum , and the minimum , for better understanding, as shown in Table 3. The comparison shows that the Pareto front distribution obtained by INSGA-II is more uniform than those obtained by the other algorithms. And, larger and and smaller occur in the case by using INSGA-II than those by using the other algorithms. This proves that INSGA-II has better population diversity and global and local search abilities and can provide more suitable solutions for designers to choose.
5.2. Comparison of the Compromise Solutions
In solving our multiobjective optimization problem of the designed rescue device, there always exists some contradiction between objects. These usually lead to different solution qualities of the Pareto set for each objective function. However, it is difficult for designers to observe and select suitable solutions from Pareto set in the multidimensional space. A method for converting the Pareto set from the multidimensional space to the two-dimensional space was proposed by Walker [28]. In the two-dimensional space, these solutions are mapped into corresponding points restricted to a specific polygon. Also, the vertices of this polygon are determined from the following formula:where is the number of objectives. These points are further mapped bywhere is the normalized value of solution j′s ith objective, is the number of solutions in the Pareto front.
In order to observe the characteristics of these points more clearly, we mark them with different colors by the following formula:
Then, the Pareto fronts obtained by different algorithms may be mapped into two-dimensional spaces, as shown in Figure 16. Blue and yellow points represent the minimum and maximum of c, respectively. If a point is closer to the vertex, it indicates that the value of the corresponding objective function is better, and the values of other objective functions are worse. The point closest to the center of the polygon has the same performance for each objective function, which is defined as the compromise solution. Eventually, the compromise solutions by different algorithms are obtained as shown in Figure 16 and Table 4.

(a)

(b)

(c)

(d)
6. Conclusion
For improving the efficiency of traditional rescue equipment, a novel rescue accessory with spreading and supporting functions is proposed in this paper. The working principle for the main components of the rescue accessory is discussed, and the kinematics analysis of the spreading mechanism is made. In analysis, the spreading width, jaw rotational angle, and spreading force are simultaneously chosen as the objective functions, and some mathematical models are established for estimating the performance of the accessory. The proposed INSGA-II modifies the crowded distance, crossover operator, and mutation operation to improve the population diversity and search capabilities of traditional NSGA-II. For comparison, a multiobjective optimization test of the rescue accessory is performed with the methods of INSGA-II, NSGA-II, PESAII, and SPEA2. Simulation results show that the Pareto front obtained by the proposed INSGA-II distributes more uniformly than those by the other algorithms, proving that INSGA-II has the highest effectiveness among all algorithms. These matters provide valuable design guides for designers. Especially, the compromise solutions of different algorithms are obtained and compared, and the comparison results show that the compromise solution of three objective functions obtained by INSGA-II is superior to those obtained by the other algorithms. This again proves that our designed novel rescue accessory has good performance compared with traditional ones.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Special Project of Central Government Guiding Local Science and Technology Development (no. 2021L3029), the Leading Project of Fujian Province (no. 2020H0048), the Leading Technology Project of Sanming City (no. 2021-G-7), the National Natural Foundation Cultivation Project of Sanming University (no. PYT2006), and the Department of Education Science and Technology Program of Fujian Province (no. JAT200641).