1. Introduction
Since flexible manufacturing systems rapidly developed in recent years, industrial robots play a more and more important role in various manufacturing scenarios [
1]. The modern application of industrial robots includes welding, assembly, painting, machining and so on. The continually expanding application of industrial robots brings new challenges, which include higher performance demands, such as faster responses, less energy consumption and greater robustness. The robot body design has a direct effect on its performance, and an excellent robot body design is generally regarded as a significant foundation for good performance [
2]. Therefore, the development of efficient industrial robot design methods has increasingly attracted wide attention. Either design optimization technology or simulation-based design technology can effectively support the design and development of industrial robots. In addition, the simulation technology and robot design optimization technology, as the main components of Industry 4.0 [
3,
4,
5], have a very important significance in the development of manufacturing industries.
The design optimization technology as an important measure has been broadly employed to improve the design and reduce the costs in many product designs, such as car design, airplane design, ship design, etc. Due to its outstanding advantages, more and more design optimization methods have been developed for robot body design. Hsiao et al. [
6] proposed a surrogate-based evolutionary optimization method, in which the response surface method is integrated with a multi-objective evolutionary algorithm. Their proposed method can not only improve the robot arm performance, but also alleviate the computational burden. Kouritem et al. [
7] proposed a multi-objective design mechanism for industrial robots. In this design mechanism, the initial and running costs, the stress analysis and the vibration analysis are considered, and therefore the optimum type of material and physical dimensions can be selected for industrial robot arms. Liang et al. [
8] proposed a structural optimization method for robot arms, which specially considers the influences of flexible joints on robot dynamics. Bugday and Karali [
9] analyzed different robot arms designs, in which both geometry and materials are changed, and their findings showed that the cost, efficiency and motor lifetime were all improved by adopting the optimization technology.
In addition to the optimization design, the simulation design is an advanced technology that is often used in modern product design. The simulation-based design is practically used to predict the behavior and results of products [
10]. Benefited from the development of computer technology, computer simulation design has been widely used in modern mechanical product design, greatly improving design efficiency and product quality. Similar to most mechanical product designs, industrial robot design also extensively adopts computer simulation methods. Ramasubramanian et al. [
11] used a commercial CAD platform to perform the initial design of a robot gripper finger. Brathikan et al. [
12] employed the dynamics simulation software ADAMS to analyze the kinematics of a five-axis industrial robotic arm. Sahu et al. [
13,
14] optimized the performance of a six-axis industrial robot on the basis of stress, modal, and vibration analyses by using FEM of the ANSYS workbench. Ghosh et al. [
15] modified the design of the robotic arm based on the analysis of von Mises stress and strain by using ANSYS simulation. Shanmugasundar et al. [
16] proposed an optimization framework to carry out a structural optimization design of a five-DOF welding robot considering the strength, stiffness, and robot weight by using ANSYS.
The above-mentioned research on the structural optimization and simulation design of robots was completed with a deterministic analysis. However, there inevitably exist various uncertainties in actual engineering [
17], for example, the variation of the coefficient related to a material’s properties derived from the discreteness of engineering material properties [
18], the variation of the dimensions related to manufacturing errors [
19], etc. The deterministic design ignores the influence of uncertainties, which generally yields unreliable design results, especially in the optimization design. In the design of an industrial robot, uncertainties are also an important factor to be considered.
Xu et al. [
20] proposed a dynamic and static structural topology optimization method with uncertain parameters for the lift arm of a parking robot. Wang et al. [
21] performed a trajectory optimization under non-probabilistic time-dependent reliability constraints to improve the soaring operating performance of robotic manipulators. Gao et al. [
22] realized the design optimization of clearance joints based on reliability sensitivity analysis. Lara-Molina et al. [
23] constructed the dynamic modeling of a one-link flexible manipulator using the stochastic finite element method in terms of the displacement of the manipulator’s tip and the frequency response function subjected to uncertainties. Lara-Molina et al. [
24] also proposed a novel robust optimal design for parallel manipulators to optimize the performance indices subject to the unavoidable effect of the uncertainties. Currently, the design optimization of a robot considering the uncertainties is mainly focused on the joints design and the trajectory optimization. As for the structural design of a robot, it is mainly performed by the deterministic design optimization.
Through the above analysis, the reliability design optimization for a robot is very important for improving the reliability of an industrial robot structural design. Moreover, improving the calculation efficiency of the method can effectively promote the application of the method in engineering. Therefore, in this paper, a simulation-based reliability design optimization method (SBRDOM) is prosed for an industrial robot structural design. In the proposed method, the simulation design is integrated with the reliability design optimization to obtain a robot structural design model, and therefore, the best balance between economy and safety could be achieved during the robot structural design. Furthermore, the SORA algorithm was employed to solve the problem of low computational efficiency of reliability design optimization in engineering applications.
The remainder of this paper is organized as follows.
Section 2 provides the proposed SBRDOM, including the SBRDOM procedure, the formulation of the reliability design optimization model and the SORA algorithm and its associated mathematical formulation. In
Section 3, the design of the big arm and small arm is used to demonstrate the validity of the proposed method. Finally,
Section 4 concludes this work and suggests future work.
2. Simulation-Based Reliability Design Optimization Method (SBRDOM)
2.1. SBRDOM Procedure
In the SBRDOM, the Latin hypercube sampling (LHS), simulation technology, response surface method (RSM) and SORA (Sequential Optimization and Reliability Assessment) algorithm were integrated to complete the structural design of the robot. The proposed method integrated the joint simulation of MTLAB and ANSYS into approximate technology so as to obtain the performance function of an industrial robot. In addition, the proposed method employed the SORA algorithm to solve the problem of low computational efficiency in engineering applications. The main procedure of the SBRDOM is shown in
Figure 1.
The main steps of the proposed method are summarized as follows.
Step 1: Determining the initial design parameters and variables of the industrial robot. In this step, the geometric configuration, load, material and working environment of the industrial robot were analyzed to define the uncertainty design parameters and variables.
Step 2: Sampling the design variables based on LHS. In this step, the LHS [
25] was used to obtain the samples of the design variables. The bounds of each variable were regarded as an interval. This interval was divided into
n small intervals, and then a sample point was randomly extracted from each small interval. The final samples were obtained by combining all the sample points.
Step 3: Constructing the performance function of the industrial robot based on computer simulation and RSM. In this step, the joint simulation of MTLAB and ANSYS was integrated with RSM to obtain the preformance function of the industrial robot. The detailed operation of this step is shown in
Figure 2 and explained as follows,
Step 3.1: Updating the simulation models. Firstly, the 3D parameter model of the industrial robot was constructed by using the ANSYS software. Then, the computer simulation model was updated with the sample point produced in step 2.
Step 3.2: Performing the simulations. The ANSYS simulation was performed with each sample, and then the response values corresponding to each sample were obtained.
Step 3.3: Approximating the performance function by using RSM. In this step, MATLAB was used to apply the approximation technology. Based on the samples and responses, the performance function could be obtained by using the approximation approach. The general used approximate model included a polynomial model, a kriging model and a support vector machine model.
Step 4: Constructing the reliability design optimization model of the industrial robot based on the probabilistic reliability theory. Based on the probability reliability theory [
26], the design parameters and variables were modeled as random parameters and variables, and then the reliability design optimization model of the industrial robot could be constructed.
Step 5: Solving the reliability design optimization problem by using the SORA algorithm. The reliability design optimization model of the industrial robot was carried out by employing the SORA algorithm in this step. The process of SORA is detailed in
Section 2.3.
Step 6: Obtaining the optimal design of the industrial robot. In this step, the optimization results obtained by SORA were evaluated. If the results were not reasonable, we turned back to Step 2 to perform the process again. The evaluation criterion was that the difference between two optimization results had to be within the allowable accuracy range.
2.2. Formulation of the Reliability Design Optimization Model
The general reliability design optimization model for the industrial robot is presented as follows:
where
f is the objective function,
is the vector of the random design variable,
is the vector of the mean values of the random design variable,
is the vector of the random design parameter,
stands for the probability,
is the
ith performance constraint function,
n is the number of the performance constraint function, and the superscripts
L and
U, respectively, stand for upper boundary and lower boundary.
The performance constraint function was obtained by using the approximate model, and therefore
was an implicit function. The expression of
depends on the selection of the approximate model. Since different approximate models have different accuracy and efficiency, the approximate model to be selected relies on the specific context. In this proposed method, the quadratic polynomial with cross terms was chosen for fitting the performance constraint function, which is given below.
where
stands for the approximation of the performance constraint function,
,
,
and
are coefficients,
xi and
xj is the
ith or the
jth random design variable,
m is the number of the random design variable.
2.3. SORA and Its Associated Mathematical Formulation
It has been practically testified that SORA [
27] is an efficient algorithm to solve reliability design optimization problems, which has been widely used due to its high efficiency. In SORA, the optimization solution process is decomposed into a deterministic optimization and an uncertainty analysis, which are executed alternately. The detailed process of SORA is shown in
Figure 3 and explained as follows.
Step 5.1: Initial value setting. The starting points including the design variables and the parameters were set. Let the shift vector , and set the iteration number k = 1. Set , , where the superscript (0) stands for the initial iteration, and and stand for the inverse most probable point (iMPP) of the random parameter and the iMPP of the random variable, respectively.
Step 5.2: Current shift vector
s calculation. The shift vector
s is calculated by
where
is the vector of the shift vector corresponding to the
ith performance constraint function,
k stands for the
kth iteration,
is the iMPP corresponding to the
ith performance constraint function.
Step 5.3: Deterministic design optimization. In this step, the deterministic design optimization was conducted, whose model is given as follows.
Step 5.4: Reliability assessment.
and
could be determined by using the reliability assessment, which was conducted through the following optimization model.
where
is the
ith performance function in standard normal space (u-space), and
is the vector of standard normal random variables and parameters, that is
.
and
are transformed from the original vector of random variables and random parameters, respectively. For the transformation formulas, please refer to Ref. [
28].
Step 5.5: Convergence judgment. Judging whether the objective function achieved convergence and whether the constraints were respected. If the convergence criterion os met, then go to step 6, otherwise let k = k + 1 and proceed to the next iteration.
4. Conclusions
Considering that the influence of uncertainties in the structural design of industrial robots is important for improving the performance of the robots, in this paper, a simulation-based reliability design optimization method is proposed to improve the reliability of robots. Four steps including LHS, approximating the performance function based on simulation, establishing the reliability design optimization model and optimization solution based on the SORA algorithm were involved in the proposed method. The effectiveness of the proposed method was verified by the design of the big arm and the small arm of an industrial robot with a load of 6 Kg. Compared with the deterministic design optimization, the proposed method appeared to be more efficient in ensuring performance reliability. The verification results showed that the proposed method can successfully realize the reliability design of industrial robots. The verification results also demonstrated that the proposed method has a better computational efficiency compared with the reliability design optimization of the double-loop method.
In this paper, the reliability design optimization of the big arm and small arm of an industrial robot was performed under conditions which did not consider the dynamic performance and the design optimization of other connecting parts, and therefore a whole reliability design optimization method can be developed on the basis of the proposed method in the future. The influence of the jump resonance on the reliability of the robotic arms was also not considered in this research. In addition, the study considered a single-disciplinary reliability design optimization. However, the structural design of industrial robots is a problem involving multiple disciplines such as mechanical and electrical system and control. Consequently, the development of a multidisciplinary reliability design optimization method simultaneously considering static performance, dynamic performance and jump resonance will be the focus of future work, so as to further improve the performance of the industrial robot.