You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
<!-- metadata commented in wiki content
2
3
4
==Adaptive ensemble surrogate-based optimization and analysis of forklift pallet racks==
5
6
<big>Wei Zhang<sup>1</sup>, Yujun Lu<sup>1</sup></big>, Yaze Mei<sup>1</sup>, Liye Lv<sup>1</sup>
7
8
1 School of Mechanical Engineering, Zhejiang Sci-Tech University, Hangzhou, Zhejiang, China, 310018
9
2 Noblelift Intelligent Equipment Co.,Ltd., Huzhou, Zhejiang, China, 313000
10
3 Zhejiang Sci-Tech University Longgang Research Institute Co., Ltd., Wenzhou, Zhejiang, China, 325802
11
-->
12
==Abstract==
13
14
As an important part of the lifting platform of pallet forklift trucks, how to reduce the deformation of the pallet rack under working conditions while reducing the mass to ensure the safety performance of forklift trucks is the most concerning issue in the design of forklift truck structure. The pallet rack structure is complex, and optimizing simulation using traditional high-precision simulation models takes much time and effort. Therefore, this paper takes the lifting platform of an unmanned AVG forklift truck as the research object and establishes a parametric model of the pallet rack using the 3D modelling software SolidWorks and the finite element analysis software ANSYS to carry out static analysis of it. Optimization design variables are selected, a single surrogate model and ensemble surrogate model are chosen for various surrogate model techniques, a small number of sample points are used to construct a low-precision model, and adaptive infilling technology is used to improve the model accuracy, and the structure is optimized, and the results show that, while the pallet rack structure meets the requirements of light weight and strength, the mass is reduced by 1.2%, and the morphology is reduced by 17.2%. Moreover, a global sensitivity analysis of each design parameter was carried out under the guidance of the surrogate model, highlighting the most influential design variable as the height of the channel steel and establishing the design variables that should be taken into account in the structural design. This paper compares the performance of the mainstream single-surrogate model and ensemble-surrogate model as well as the adaptive infilling strategy in engineering design and proves that the surrogate model optimization method has a certain guiding significance for the structural optimization design of pallet racking.
15
16
'''Keywords''':Optimization design, surrogate models, pallet racks, DoEs, infilling
17
18
==1. Introduction==
19
20
AGV forklift is an important unmanned material handling robot under the trend of intelligent industry, and its mature technology makes it widely used in stations, ports, airports, factories and warehouses, AGV forklifts are highly efficient equipment for mechanised loading and unloading, stacking and short-distance transport. Among them, a pallet forklift can ensure the stability of the material itself in the case of uneven road surfaces and heavy loads of forklifts, which is an important part of the forklift branch. Therefore, AGV pallet forklift is the current research focus of unmanned material handling robots.As an important component of the pallet forklift working device, the pallet rack is the executor of the forklift to realize the handling function and the focus point of value embodiment. Its performance directly affects the efficiency of the forklift. Traditionally, the optimal design of pallet rack is mainly done through computer-aided design, and good results are obtained by simulation optimization, but there are many disadvantages in the simulation itself: the simulation process is time-consuming, and there is a lack of optimization algorithms for the time-consuming simulation, which results in the whole simulation process needing a lot of time and computational costs. The structural design of pallet racks involves multiple variables and highly nonlinear relationships, in this context, the surrogate modelling technique, which has a small amount of computation and is based on a small amount of sample data, is an excellent choice [1], and optimization based on the surrogate model technique can reduce the computation time and improve the optimization efficiency.
21
22
In recent years, the field of surrogate modelling has continued to develop, giving rise to numerous new methods and varieties, which have achieved excellent results in the fields of mechanical, materials, fluid and other engineering fields. Queipo et al. [2] provide a comprehensive discussion of the fundamental issues that arise from surrogate model-based analysis and optimization. In 2012, Song et al. [3] used multiple surrogate models to perform collisional foam-filled conical thin-walled structure optimization, and the results show that different surrogate model approaches produce different results for the final collision optimization design. Chen et al. [4] performed a multi-objective optimization design of an injection model based on the KRG model and the particle swarm optimization algorithm. In 2020, Jiang et al. [5,6] systematically elaborated on the application of surrogate models in the field of engineering and used a variable fidelity surrogate model to design the micro air vehicle fuselage as well as reinforced cylindrical shells. In 2020, Meng et al. [7] optimized the robust design of incomplete reinforcing bars based on an ensemble surrogate model. Yondo et al. [8] compared the surrogate model approach and the experimental design approach in the aerodynamic analysis of a flight vehicle try and introduced the optimization method based on the surrogate model. In 2021, Zhang et al. [9] proposed an optimization framework based on a multi-fidelity surrogate model by correlating the configuration parameters of an aircraft with its start-up performance.
23
24
Prediction is one of the most important functions of the surrogate model, and one of the important methods to improve the prediction accuracy of the surrogate model is to update the surrogate model by infilling the criterion. The more commonly used mainstream infilling criteria are the maximizing the probability of improvement (PoI) [10] criterion and maximizing expectation improvement (EI) [11] criterion, of which the EI method has better global search performance and a wide range of application scenarios and has become one of the most mainstream infilling criteria, which also has a lot of applications in engineering design. Zeng and Zhou [12] proposed a KRG model optimization method integrating LCB and trust domain and successfully applied it to the design optimization of I-beam. Kleijnen [13] proposed a bootstrap EI method, which can realize the multi-point infilling process based on EI and improve the modelling efficiency. Chaudhuri et al. [14] ensemble EI and PoI criteria to update the surrogate model every time by obtaining multiple sample points obtain multiple sample points to update the surrogate model and apply the method to the optimized design of vibrating wings. Liu et al. [16] converted the classical four criteria EI, PI, LCB, and MP into an infilling criterion with constraints and applied it to the wing and aerodynamic inverse design problems, and obtained better results.
25
26
In this paper, an unmanned AGV forklift truck is selected as the research object, focusing on the analysis and improvement of the pallet rack part of the forklift truck, to achieve lightweight under the premise of ensuring the structural strength and the optimization makes the forklift truck have a higher structural strength, a better loading capacity, and a lighter mass. At the same time, to compare the difference between the ensemble surrogate model and the single surrogate model in engineering applications, this paper chooses the single surrogate model represented by the Kriging(KRG) [10] model and the ensemble surrogate model represented by the extended adaptive hybrid functions(E-AHF) [15] to optimize the design and analysis of the forklift truck pallet racks, respectively. Among them, the KRG model, as one of the typical single-surrogate models, is essentially a linearly weighted combination of information about known points in order to predict the position information within a certain range. Unlike the rest of the single surrogate models, the KRG model has a better local prediction ability due to the correlation function, good continuity and derivability, and good approximation for nonlinear complex problems. The E-AHF ensemble surrogate model, one of the newly proposed ensemble surrogate models in recent years, has excellent predictive ability due to its extensive model library and the excellent weighting method used. It has better global performance for highly nonlinear problems. Therefore, the above two surrogate models are chosen to represent their respective fields. In order to further reduce the cost of the simulation model, this paper chooses to construct the low-precision surrogate model with less sample points. To further reduce the cost of the simulation model, this paper chooses to construct a low-precision surrogate model with fewer sample points and to improve the model accuracy by using the adaptive infilling criterion to reduce the number of times the simulation software is used. Finally, the most suitable model is selected from multiple error criteria.
27
28
The rest of the paper is arranged as follows: Section 2 briefly introduces the mathematical principles and related formulas involved in this paper, which are mainly related to the calculation of the surrogate model; Section 3 introduces the research object of this paper, the unmanned AGV stacker truck, and constructs the related finite element model; Section 4 constructs the relevant surrogates for the pallet rack with a single model and an ensemble model respectively; Section 5 compares the performance of the two surrogate models in this optimal design, selects the most suitable model, and based on which the design and optimization of the pallet rack are implemented; and the conclusions of this study are discussed in Section 6.
29
30
==2. Mathematical theories==
31
32
===2.1 KRG model===
33
34
The KRG model [10], originally proposed by the South African geologist Danie Krige in 1951 for the exploration of mineral reserves, is an uncertainty-free interpolation model for estimating standard deviation minimization, which uses a covariance-controlled Gaussian process to generate continuous functions.
35
36
A generic KRG model can be expressed as follows:
37
38
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
39
|-
40
| 
41
{| style="text-align: center; margin:auto;width: 100%;" 
42
|-
43
| style="text-align: center;" | <math>y\left( \boldsymbol{x}\right) =f\left( \boldsymbol{x}\right) +z\left( \boldsymbol{x}\right)</math> 
44
|}
45
|  style="text-align: center;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(1)
46
|}
47
48
where  <math display="inline">y\left( \boldsymbol{x}\right)</math>  is the response function for the region of interest,  <math display="inline">f\left( \boldsymbol{x}\right)</math>  is a polynomial function used as an approximation to the underlying functional model, and <math display="inline">z\left( \boldsymbol{x}\right)</math>  denotes a stochastic process with uncertainty in the mean of <math display="inline">f\left( \boldsymbol{x}\right)</math>. The KRG model is an interpolation method through each sample point. 
49
50
In this paper, we use the Gaussian correlation function:
51
52
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
53
|-
54
| 
55
{| style="text-align: center; margin:auto;width: 100%;" 
56
|-
57
| style="text-align: center;" | <math>R\left( {\boldsymbol{x}}^{i},{\boldsymbol{x}}^{j}\right) =\mathrm{exp}\,\left[ -\sum _{k=1}^{{n}_{dv}}{\theta }_{k}{\left| {x}_{k}^{i}-{x}_{k}^{j}\right| }^{{P}_{k}}\right]</math> 
58
|}
59
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(2)
60
|}
61
62
where <math display="inline">{n}_{dv}</math> is the number of design variables, <math display="inline">\theta</math>  and  <math display="inline">P</math> are parameters that vary flexibly with the variables, and  <math display="inline">\left| {x}_{k}^{i}-{x}_{k}^{j}\right|</math> is the distance between  <math display="inline">{\boldsymbol{x}}^{i}</math> and  <math display="inline">{\boldsymbol{x}}^{j}</math> concerning the k''-''th variable. 
63
64
The predicted value of  <math display="inline">\hat{y}\left( \boldsymbol{x}\right)</math>  at test point <math> x </math> is formulated as follows:
65
66
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
67
|-
68
| 
69
{| style="text-align: center; margin:auto;width: 100%;" 
70
|-
71
| style="text-align: center;" | <math>\hat{y}\left( x\right) =\hat{\beta }+\boldsymbol{r}\left( x\right) {\boldsymbol{R}}^{-1}\left( y-\boldsymbol{f}\hat{\beta }\right)</math> 
72
|}
73
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(3)
74
|}
75
76
where  <math display="inline">\boldsymbol{y}</math> is a column vector of sample response values, and when  <math display="inline">f\left( \boldsymbol{x}\right)</math>  is a constant vector,  <math display="inline">\boldsymbol{f}</math> is a column vector of length <math> n</math> with all elements of 1. <math display="inline">\boldsymbol{r}\left( x\right)</math>  is the correlation vector between the test point<math> x </math> and the sample points.
77
78
===2.2 Ensemble surrogate models===
79
80
Ensemble surrogate models are obtained by superimposing the weights of several individual surrogate models, and the key lies in solving for the appropriate weight coefficients. The current ensemble surrogate models can be divided into two types, namely, average weights and adaptive weights, according to the different forms of weight coefficients. This paper focuses on the adaptive weighted ensemble surrogate model used.
81
82
In 2018, Song et al. [15] proposed an ensemble surrogate model: Extended Adaptive Hybrid Functions (E-AHF). Three steps are proposed for model screening and adaptive weight calculation.
83
84
'''Step 1. Local measure estimation'''
85
86
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
87
|-
88
| 
89
{| style="text-align: center; margin:auto;width: 100%;" 
90
|-
91
| style="text-align: center;" | <math>{\hat{s}}^{2}(x)={\sigma }^{2}\left[ 1-{\boldsymbol{\psi }}^{T}{\boldsymbol{\Psi }}^{-1}\boldsymbol{\psi }+\right. </math><math>\left. \frac{\mathit{\boldsymbol{1}}-{\boldsymbol{1}}^{T}{\boldsymbol{\Psi }}^{-1}\boldsymbol{\psi }}{{\boldsymbol{1}}^{T}{\boldsymbol{\Psi }}^{-1}\boldsymbol{1}}\right]</math> 
92
|}
93
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(4)
94
|}
95
96
where  <math display="inline">\sigma</math> denotes the constant process variance of a Gaussian field, and <math display="inline">\boldsymbol{\Psi }</math> is the correlation matrix with the elements. 
97
98
'''Step 2. Probability estimation'''
99
100
The baseline model can represent the global trend of the hybrid surrogate model, due to its high accuracy across the entire design space. Therefore, this prediction by the baseline model can be deemed as an expected value of the hybrid model. With the local prediction estimation  <math display="inline">{s}^{2}</math> obtained from the previous local measure estimation, the probability coefficient of each individual surrogate model can be calculated by
101
102
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
103
|-
104
| 
105
{| style="text-align: center; margin:auto;width: 100%;" 
106
|-
107
| style="text-align: center;" | <math>{P}_{i}=exp\left\{ -\frac{{\left[ {y}_{ij}-{y}_{basej}\right] }^{2}}{2{s}_{j}^{2}}\right\}</math> 
108
|}
109
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(5)
110
|}
111
112
where  <math display="inline">{y}_{basej}</math> is the prediction value of the baseline model at the j-th sample point, Pi is the probability coefficient associated with the i-th individual surrogate model, and  <math display="inline">{s}_{j}^{2}</math> is the local prediction esti-mation at the j-th sample point.
113
114
'''Step 3. Local weight determination'''
115
116
With the probability coefficient assigned to each surrogate model, the weight factors can be computed by
117
118
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
119
|-
120
| 
121
{| style="text-align: center;margin:auto;width: 100%;"
122
|-
123
| <math>{\omega }_{i}=\frac{{P}_{i}}{\sum _{j=1}^{M}{P}_{j}}</math>
124
|}
125
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(6)
126
|}
127
128
As  <math display="inline">P</math> is a function of the input vector  <math display="inline">\boldsymbol{x}</math>, the weight factors obtained for each individual surrogate model is adaptive to the input vector as well.
129
130
===2.3 Maximizing EI===
131
132
Unlike the PoI strategy, which focuses on the probability of improvement, the EI strategy is related to the amount of improvement expected to be obtained when adding a filling point and given a filling point, the EI is expressed in terms of a Gaussian error function as follows:
133
134
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
135
|-
136
| 
137
{| style="text-align: center;margin:auto;width: 100%;"
138
|-
139
| <math>E\left[ I\left( x\right) \right] =\left( {y}_{min}-\hat{y}\left( x\right) \right) \left[ \frac{1}{2}+ \frac{1}{2}\mathrm{erf}\,\left( \frac{{y}_{min}-\hat{y}\left( x\right) }{\hat{s}\sqrt{2}}\right) \right] +</math><math>\hat{s}\frac{1}{\sqrt{2\pi }}exp\left[ \frac{-({{y}_{min}-\hat{y}(x))}^{2}}{2\hat{s}}\right]</math> 
140
|}
141
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(7)
142
|}
143
144
It is observed from equations <math display="inline">P\left[ I\left( x\right) \right]</math> and  <math display="inline">E\left[ I\left( x\right) \right]</math>  that the expressions for both PoI and EI strategies contain global terms (<math display="inline">{s}^{2}</math>) and local terms (<math display="inline">I=</math><math>{y}_{min}-Y\left( x\right)</math>), indicating that both approaches can balance global exploration and local exploitation to some extent.
145
146
===2.4 Performance criterion===
147
148
In this paper, we use two criteria to evaluate the accuracy of the model. That is, the coefficient of determination <math display="inline">\left( {R}^{2}\right)</math> is used as a global evaluation metric; the normalized max absolute error (NMAE) is used as a local evaluation metric to evaluate the global and local performance of different infilling strategies, respectively
149
150
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
151
|-
152
| 
153
{| style="text-align: center;margin:auto;width: 100%;"
154
|-
155
| <math>{R}^{2}=1-\displaystyle\frac{\displaystyle\sum\limits_{i=1}^{ntest}{\left( {y}_{i}-{\hat{y}}_{i}\right) }^{2}}{\displaystyle\sum \limits_{i=1}^{ntest}{\left( {y}_{i}-\overline{y}\right) }^{2}}</math>
156
|}
157
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(8)
158
|}
159
160
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
161
|-
162
| 
163
{| style="text-align: center;margin:auto;width: 100%;"
164
|-
165
| <math>NMAE=\max\left| y-\hat{y}\right| \sqrt{\displaystyle\frac{1}{ntest}\sum\limits _{i=1}^{ntest}{\left( {y}_{i}-{\hat{y}}_{i}\right) }^{2}}</math>
166
|}
167
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|(9)
168
|}
169
170
where  <math display="inline">{y}_{i}</math> denotes the true response at the i-th test point, and the predicted value of the model built by the surrogate model approach at the i-th test point is  <math display="inline">{\hat{y}}_{i}</math>,   and <math display="inline">\overline{y}</math> is the mean value of the true response at all test points, ntest denotes the number of test points. 
171
172
Usually, <math display="inline">{R}^{2}</math> ranges from 0 to 1, and a larger d or a smaller NMAE means more accurate surrogates. When  <math display="inline">{R}^{2}</math> is negative, it represents that the constructed surrogate model is trending opposite to the true model.
173
174
==3. Establishment of finite element model of pallet rack==
175
176
===3.1 Structure of pallet racks===
177
178
The AGV forklift involved in this paper is a battery-powered, unmanned material handling vehicle with automatic guidance, and its lifting platform consists of a pallet rack, two-wheel frame plates, a lifting plate, a guide rail, and a base plate, as shown in [[#img-1|Figure 1]]. Among them, the pallet rack, as the structural main body of the lifting platform, is used to bear all the weight of the pallet containing materials and the loaded materials and is all the force components of the lifting platform, as well as the key structure of the lifting platform. As an important part of the AGV forklift, the reasonable design of the lifting platform directly determines the working ability of the forklift. In daily operation, the pallet rack will be subjected to more complicated forces in repeated lifting operations. Therefore, the structural design and optimization of the pallet rack will directly lead to the degree of deformation and reliability of each component of the AGV forklift in the course of operation, affecting the safety performance and working efficiency of the lifting platform, which in turn seriously affects the overall working life of the forklift.
179
180
<div id='img-1'></div>
181
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
182
|-style="background:white;"
183
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image1.png|394px]]
184
|-
185
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 1'''. Structure of lifting platform
186
|}
187
188
===3.2 Simulation analysis of pallet racks===
189
190
A 3-D model of the pallet rack was built using SolidWorks and imported into ANSYS for static analysis. The material used for the pallet rack is 25MnV alloy steel, Set the material properties in Workbench: density 7850 <math display="inline">\frac{kg}{{m}^{3}}</math>, modulus of elasticity <math display="inline">2.06{_\ast}{10}^{5}</math>Mpa, Poisson's ratio of 0.3. The free mesh division of the pallet rack is carried out, with a grid cell size of 10mm, and local encryption is carried out on the parts that are prone to stress concentration, and a total of 191599 mesh cells are generated with 424581 nodes. The main loads and constraints on the pallet rack include: fixed constraints on the lifting plate of the fixed pallet rack, which is used to limit the vertical displacement of the pallet rack; the downward pressure acting on the side plates of the channel steel <math> F </math>; and the gravity of the pallet rack itself <math> G </math>. The overall load diagram of the pallet rack is shown in [[#img-2|Figure 2]].
191
192
<div id='img-2'></div>
193
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
194
|-style="background:white;"
195
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image2.png|382px]]
196
|-
197
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 2'''. Overall load diagram of pallet racks
198
|}
199
200
==4. Constructing surrogate model for pallet racking==
201
202
Surrogate model method is a mathematical analysis model method, which belongs to a branch of supervised machine learning, and has gained wide attention and research in many fields such as mechanical simulation, ecological modelling, geological statistics, artificial intelligence, etc., due to its ability to quickly and effectively achieve the approximate fitting of multidimensional variable complex systems based on limited data information. Surrogate modeling technology does not need to understand the physical meaning between input variables and output responses, and can only use the data itself to obtain an approximate mathematical model, so it can be applied to the optimization design and analysis where the physical meaning between the inputs and outputs is not clear, and the computational cost is high, and it has the advantages of good fitting accuracy, low cost, high efficiency, and so on. The application process of the surrogate model in engineering in general is shown in [[#img-3|Figure 3]], which mainly includes three modules: problem analysis, modelling analysis and result analysis.
203
204
<div id='img-3'></div>
205
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
206
|-style="background:white;"
207
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image3.jpeg|600px]]
208
|-
209
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 3'''. Flow of application of surrogate model in engineering
210
|}
211
212
213
===4.1 Identification of design variables===
214
215
As the main stress structure of the pallet rack, the channel steel and diagonal reinforcement is an important structure to support the pallet rack frame and transfer its stress, and its strength also determines whether the forklift truck can complete the engineering transport work more efficiently and safely. Therefore, the structural optimization of channel steel as well as diagonal tie bars can further reduce the stress and deformation while ensuring the strength, reduce the mass, and improve the economy and overall stability. After clarifying the overall structural weaknesses, this paper selects four design variables for the channel section and five design variables for the thickness of the diagonal bar as the optimization parameters, as shown in [[#img-4|Figure 4]], which are the height  <math display="inline">{x}_{1}</math>, width  <math display="inline">{x}_{2}</math>, web thickness  <math display="inline">{x}_{3}</math> and flange thickness of the channel  <math display="inline">{x}_{4}</math> and the thickness of the diagonal bar  <math display="inline">{x}_{5}</math>, respectively. The range and initial values of the corresponding design variables are shown in [[#tab-1|Table 1]].
216
217
<div id='img-4'></div>
218
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
219
|-style="background:white;"
220
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image4.jpeg|324px]] 
221
|-
222
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 4'''. Schematic diagram of the design variables of the pallet racks
223
|}
224
225
226
<div class="center" style="font-size: 75%;">'''Table 1'''. Design variable size parameters and range of values</div>
227
228
<div id='tab-1'></div>
229
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
230
|-style="text-align:center"
231
! Design variables !! Range of values (mm) !! Initial value (mm)
232
|-
233
|  style="text-align: center;"|Height of channel  <math display="inline">{x}_{1}</math>
234
|  style="text-align: center;"|90~120
235
|  style="text-align: center;"|105
236
|-
237
|  style="text-align: center;"|Width of channel  <math display="inline">{x}_{2}</math>
238
|  style="text-align: center;"|37~53
239
|  style="text-align: center;"|45
240
|-
241
|  style="text-align: center;"|Thickness of web  <math display="inline">{x}_{3}</math>
242
|  style="text-align: center;"|6.5~8.5
243
|  style="text-align: center;"|7.5
244
|-
245
|  style="text-align: center;"|Flange thickness  <math display="inline">{x}_{4}</math>
246
|  style="text-align: center;"|7.5~9
247
|  style="text-align: center;"|8.25
248
|-
249
|  style="text-align: center;"|Thickness of diagonal reinforcement  <math display="inline">{x}_{5}</math>
250
|  style="text-align: center;"|12~20
251
|  style="text-align: center;"|16
252
|}
253
254
255
===4.2 Design of experiment===
256
257
In the application of the surrogate model in engineering shown in [[#img-4|Figure 4]], the design of experiments (DoE) is the first step in constructing the surrogate model, and it is also the key step in constructing the surrogate model, the essence of which is to generate the sample points in the region of interest, and a good DoE method is able to extract more compliant sample points and to construct a more accurate model with a smaller amount of data. The DoE methods can be divided into static sampling and adaptive sampling. Static sampling, means that the training points needed to construct the surrogate model are extracted at one time. Currently, widely used static sampling methods include Latin hypercube sampling (LHS) [16], Fractional factorial design (FFD) [17] and central composite design (CCD) [18]. In static sampling, the training points are independent of the constructed surrogate model, i.e., static sampling cannot learn from known samples.
258
259
The adaptive infilling strategy is designed to overcome the limitations of static sampling independent models, using the results of the currently constructed surrogate model, and selecting new sample points that are more conducive to the construction of an accurate model through the "infilling strategy" to improve the accuracy of the surrogate model iteratively [19]. Typical infilling strategies include Statistical lower bound (SLB) [20] strategy, maximizing probability of improvement (PoI) strategy, maximizing expected improvement (EI) strategy, etc. The EI strategy has been used by most researchers, and the EI strategy has been used by most researchers. etc. Among them, the EI strategy has been confirmed by most studies to have better performance in local exploitation and global exploration, and it is one of the most widely used infilling strategies nowadays. In this paper, the EI strategy ischosen as the infilling strategy.
260
261
According to the process of constructing the surrogate model, after determining the range of design variables, LHS is used to select points for the experimental design, 60 points are selected, and the corresponding maximum stress and maximum deformations under the relevant loads are calculated using simulation, of which 10 groups of data are used to construct the low-precision surrogate model, and the remaining 50 groups of data are used as a database, which is used to provide response values for the new samples added at each iteration using the EI adaptive infilling strategy. The remaining 50 data sets are used as a database to provide response values for new sample points added at each iteration when using the EI adaptive infilling strategy, and the 50 data sets can also be used as a test set to examine the global and local accuracy (<math display="inline">{R}^{2}</math> and NMAE) of the constructed surrogate models. [[#img-5|Figure 5]] shows the distribution of the 60 groups of design variables in the design space, and it can be clearly seen that the training and test points are distributed in almost the whole design space with good space filling and projection, and the test and training points are interspersed with each other, which indicates that the sampling point analysis meets the requirements.
262
263
<div id='img-5'></div>
264
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
265
|-style="background:white;"
266
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image5.jpeg|600px]] 
267
|-
268
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 5'''. Distribution of training and test points
269
|}
270
271
===4.3 Constructing surrogate model for pallet racking===
272
273
In order to better reduce the cost brought by simulation experiments, this paper uses the EI infilling method to gradually guide the updating of the model and improve the accuracy of the surrogate model, therefore, when constructing the surrogate model, 10 sets of data will be used to construct the low-precision single surrogate model KRG and the low-precision ensemble surrogate model E-AHF, respectively, and according to the flow of the modelling in [[#img-3|Figure 3]], the model will be judged based on the accuracy judging criterion (<math display="inline">{R}^{2}</math> and NMAE) for the accuracy, and use the EI strategy to iteratively improve the accuracy until the surrogate model reaches a usable accuracy (usually <math display="inline">{R}^{2}\geq 0.8</math>), and in order to compare the performance of the single surrogate model and the ensemble surrogate model in this optimization, the construction of the surrogate model will be guided by the KRG model and the E-AHF model to construct the models of the design variables with respect to the deformation, mass, and stresses, respectively. The bar charts in [[#img-6|Figure 6]] show the performance of the two different models with 10 sample points with global metrics (<math display="inline">{R}^{2}</math>) and local metrics (NMAE) as references.
274
275
<div id='img-6'></div>
276
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
277
|-style="background:white;"
278
|align="center" | 
279
{|style="margin: 0em auto 0.1em auto;width:auto;" 
280
|+
281
|-
282
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image6-c.png|174px]]
283
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image7-c.png|174px]]
284
|style="text-align: center;padding:10px;"|[[Image:Review_120867196377-image8-c.png|174px]]
285
|-
286
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(a) Accuracy for deformation
287
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(b) Accuracy for stress
288
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(c) Accuracy for mass
289
|}
290
|-
291
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 6'''. Comparison of low precision surrogate models
292
|}
293
294
295
It is obvious through [[#img-6|Figure 6]] that 10 sample points are not capable of constructing a high-precision model, but in the case of limited data, the ensemble surrogate model E-AHF still exhibits higher accuracy, and the KRG model is superior in terms of scope of applicability and degree of accuracy in the single surrogate model, but the comprehensive comparison of the three response objectives of pallet racking's deformation, stress, and mass is not as good as the E-AHF model, whether it is in terms of global accuracy or local accuracy.
296
297
After the low-precision model is constructed, the EI strategy is used to guide the selection of new sample points based on the initial model, and one new sample point is selected to be added to the training set at each iteration until the model stops iterating when it reaches a usable accuracy. [[#img-7|Figure 7]] demonstrates a comparison of the accuracy of the KRG model and the E-AHF model after the convergence of the iterations. [[#img-8|Figure 8]] compares the number of sample points used to construct the model after reaching the usable accuracy for both models. It can be seen that the models constructed by both KRG and E-AHF reach high accuracy after the EI infilling strategy, but comparing the data in [[#img-7|Figure 7]], the E-AHF model still has an advantage in terms of global accuracy and local accuracy in deformation, stress and mass. Comparing the line graph in [[#img-8|Figure 8]], it can be seen that the same model that achieves <math display="inline">{R}^{2}\geq 0.8</math>, the E-AHF model can achieve the corresponding accuracy with fewer training points, which means that fewer times of the EI infilling strategy is used, and when put in the simulation experiments, it means that fewer times of the simulation are carried out, which saves the cost and time, and objectively proves the superiority of the ensemble surrogate model, which is represented by the E-AHF, compared with the single surrogate model. The exact number of additional sample points is shown in [[#tab-2|Table 2]].
298
299
<div id='img-7'></div>
300
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
301
|-style="background:white;"
302
|align="center" | 
303
{|style="margin: 0em auto 0.1em auto;width:auto;" 
304
|+
305
|-
306
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image9-c.png|174px]] 
307
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image10-c.png|174px]] 
308
|style="text-align: center;padding:10px;"|[[Image:Review_120867196377-image11-c.png|174px]]
309
|-
310
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(a) Accuracy for deformation
311
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(b) Accuracy for stress
312
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(c) Accuracy for mass
313
|}
314
|-
315
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 7'''. Comparison of high-precision surrogate models
316
|}
317
318
319
<div id='img-8'></div>
320
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
321
|-style="background:white;"
322
|align="center" | 
323
{|style="margin: 0em auto 0.1em auto;width:auto;" 
324
|+
325
|-
326
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image12-c.png|270px]]
327
|style="text-align: center;padding:10px;"| [[Image:Review_120867196377-image13-c.png|270px]]
328
|-
329
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(a) KRG
330
|style="text-align: center;font-size: 75%;padding-bottom:10px;"|(b) E-AHF
331
|}
332
|-
333
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 8'''. Comparison of the number of infilling points for the two surrogate models
334
|}
335
336
337
<div class="center" style="font-size: 75%;">'''Table 2'''. Comparison of the number of infilling points for the two surrogate models</div>
338
339
<div id='tab-1'></div>
340
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
341
|-style="text-align:center"
342
!  rowspan="2" | Response !!   colspan='2'  style="text-align: center;"|KRG !!   colspan='2'  style="text-align: center;"|E-AHF 
343
|-style="text-align:center"
344
! !! Initial points !! Additional points !!  Initial points !! Additional points
345
|-
346
|  style="text-align: center;"|Deformation
347
|  style="text-align: center;"|10
348
|  style="text-align: center;"|8
349
|  style="text-align: center;"|10
350
|  style="text-align: center;"|7
351
|-
352
|  style="text-align: center;"|Stress
353
|  style="text-align: center;"|10
354
|  style="text-align: center;"|11
355
|  style="text-align: center;"|10
356
|  style="text-align: center;"|11
357
|-
358
|  style="text-align: center;"|Mass
359
|  style="text-align: center;"|10
360
|  style="text-align: center;"|5
361
|  style="text-align: center;"|10
362
|  style="text-align: center;"|6
363
|}
364
365
==5. Optimization results and analysis==
366
367
===5.1 Definition of problem===
368
369
The KRG model and the E-AHF model presented above were used instead of the simulation model for optimization and design. In the optimization of the entire pallet rack, deformation and mass are the two objectives that need to be paid attention to during the optimization process, but in fact, these two optimization objectives are contradictory to each other, as a less deformed structure is more rigid and corresponds to a higher mass, so it is difficult to balance the two ends of the spectrum during the design process. Reviewing the specific situation in engineering applications, it is not difficult to find that compared with lightweight, higher structural strength is the focus of design considerations, higher strength to ensure the safety factor of machinery in the project, as well as a longer working life, placed in this optimization design process, smaller deformation is the priority of the experimental findings of the conflict between the object.
370
371
Therefore, in this paper, we will calculate the optimal solution of the two constraint problems for the optimization objective and compare the experimental results to conclude with the same number of total constraints, which are stress constraints; deformation constraints; and mass constraints, respectively. Where [T] and [M] are the permitted values of deformation and mass respectively. Considering the specific application and working conditions, the geometries of non-standard channels and diagonal ties are constrained. For the optimization problem described below, this paper uses a single-objective Genetic algorithm (GA), and the optimization problem for the whole pallet rack is described as:
372
373
Question 1: Deformation constraints, stress constraints and dimensional constraints.
374
375
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
376
|-
377
| 
378
{| style="text-align: center;margin:auto;width: 100%;"
379
|-
380
| <math>min\quad F({x}_{1},{x}_{2},{x}_{3},{x}_{4},{x}_{5})=M</math>
381
382
<math>Subject\, \, to\quad \, {\sigma }_{max}\leq \left[ \sigma \right] =648.1</math>
383
384
<math>{T}_{max}\leq 2.387</math>
385
386
<math>90\leq {x}_{1}\leq 120</math>
387
388
<math>37\leq {x}_{2}\leq 53</math>
389
390
<math>6.5\leq {x}_{3}\leq 8.5</math>
391
392
<math>7.5\leq {x}_{4}\leq 9</math>
393
394
<math>12\leq {x}_{5}\leq 20</math>
395
|}
396
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(10)</span>
397
|}
398
399
400
<span style="text-align: center; font-size: 75%;">Question 2: Mass constraints, stress constraints and dimensional constraints.</span>
401
402
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
403
|-
404
| 
405
{| style="text-align: center;margin:auto;width: 100%;"
406
|-
407
| <math>min\quad F({x}_{1},{x}_{2},{x}_{3},{x}_{4},{x}_{5})=T</math>
408
409
<math>Subject\, \, to\quad \, {\sigma }_{max}\leq \left[ \sigma \right] =648.1</math>
410
411
<math>{M}_{max}\leq 252.78</math>
412
413
<math>90\leq {x}_{1}\leq 120</math>
414
415
<math>37\leq {x}_{2}\leq 53</math>
416
417
<math>6.5\leq {x}_{3}\leq 8.5</math>
418
419
<math>7.5\leq {x}_{4}\leq 9</math>
420
421
<math>12\leq {x}_{5}\leq 20</math>
422
|}
423
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(11)</span>
424
|}
425
426
427
<span style="text-align: center; font-size: 75%;">where </span> <math display="inline">F(x)</math><span style="text-align: center; font-size: 75%;"> is the optimization objective, </span> <math display="inline">T</math><span style="text-align: center; font-size: 75%;"> is the total deformation and </span> <math display="inline">M</math><span style="text-align: center; font-size: 75%;"> is the total weight of the pallet rack structure.</span>
428
429
The specific results of the two optimization problems are shown in Tab. 3 and Tab. 4 respectively and compare the predicted and true responses of the objective function and constraints. From the information in the table, it can be seen that the structural value derived from minimizing mass as the optimization objective with deformation and stress as constraints, Although the mass reduction is larger, the optimized resultant deformation and stress reduction is small, and this degree of lightweight has limited cost reduction in practical engineering applications, so the results of the optimization problem I am not considered. On the other hand, the structural values derived by using mass and stress as constraints and minimizing deformation as the optimization objective not only better reduce the deformation of the pallet rack, but also significantly reduce the stress values, and even a small amount of reducing the mass of the pallet rack to meet the constraints is a better result. It is worth mentioning that in both optimization problems, the E-AHF model shows better performance in terms of all the indicators, and the constraint and objective functions written in its model give better guidance to the GA algorithm to get better results.
430
431
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
432
<span style="text-align: center; font-size: 75%;">Table 3 Comparison of optimization results of the KRG and E-AHF for question 1</span></div>
433
434
{| style="width: 100%;border-collapse: collapse;" 
435
|-
436
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Response</span>
437
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|
438
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{1}</math>
439
440
<span style="text-align: center; font-size: 75%;">(mm)</span>
441
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{2}</math>
442
443
<span style="text-align: center; font-size: 75%;">(mm)</span>
444
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{3}</math>
445
446
<span style="text-align: center; font-size: 75%;">(mm)</span>
447
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{4}</math>
448
449
<span style="text-align: center; font-size: 75%;">(mm)</span>
450
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{5}</math>
451
452
<span style="text-align: center; font-size: 75%;">(mm)</span>
453
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Deformation</span>
454
455
<span style="text-align: center; font-size: 75%;">(mm)</span>
456
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Stress</span>
457
458
<span style="text-align: center; font-size: 75%;">(Mpa)</span>
459
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Quality</span>
460
461
<span style="text-align: center; font-size: 75%;">(Kg)</span>
462
|-
463
|  colspan='2'  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Before optimization</span>
464
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">105.0</span>
465
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">45.0</span>
466
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.5</span>
467
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.25</span>
468
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">16.0</span>
469
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.39</span>
470
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">648.1</span>
471
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">252.7</span>
472
|-
473
|  rowspan='1' style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">After</span>
474
475
<span style="text-align: center; font-size: 75%;">optimization</span>
476
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">E-AHF</span>
477
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">107.6</span>
478
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">41.9</span>
479
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.5</span>
480
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.0</span>
481
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">12.3</span>
482
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.19</span>
483
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">446.0</span>
484
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">246.3</span>
485
|-
486
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">KRG</span>
487
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">105.8</span>
488
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">44.3</span>
489
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.9</span>
490
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.1</span>
491
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">12.2</span>
492
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.25</span>
493
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">428.2</span>
494
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">246.9</span>
495
|}
496
497
498
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
499
<span style="text-align: center; font-size: 75%;">Table 4 Comparison of optimization results of the KRG and E-AHF for question 2</span></div>
500
501
{| style="width: 100%;border-collapse: collapse;" 
502
|-
503
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Response</span>
504
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|
505
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{1}</math>
506
507
<span style="text-align: center; font-size: 75%;">(mm)</span>
508
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{2}</math>
509
510
<span style="text-align: center; font-size: 75%;">(mm)</span>
511
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{3}</math>
512
513
<span style="text-align: center; font-size: 75%;">(mm)</span>
514
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{4}</math>
515
516
<span style="text-align: center; font-size: 75%;">(mm)</span>
517
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math>{x}_{5}</math>
518
519
<span style="text-align: center; font-size: 75%;">(mm)</span>
520
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Deformation</span>
521
522
<span style="text-align: center; font-size: 75%;">(mm)</span>
523
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Stress</span>
524
525
<span style="text-align: center; font-size: 75%;">(Mpa)</span>
526
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Quality</span>
527
528
<span style="text-align: center; font-size: 75%;">(Kg)</span>
529
|-
530
|  colspan='2'  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Before optimization</span>
531
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">105.0</span>
532
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">45.0</span>
533
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.5</span>
534
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.25</span>
535
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">16.0</span>
536
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.39</span>
537
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">648.1</span>
538
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">252.7</span>
539
|-
540
|  rowspan='1' style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">After</span>
541
542
<span style="text-align: center; font-size: 75%;">optimization</span>
543
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">E-AHF</span>
544
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">109.8</span>
545
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">42.7</span>
546
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.9</span>
547
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">8.7</span>
548
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">13.2</span>
549
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.98</span>
550
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">407.7</span>
551
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">249.7</span>
552
|-
553
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">KRG</span>
554
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">110.6</span>
555
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">41.8</span>
556
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">6.7</span>
557
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">8.6</span>
558
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">13.9</span>
559
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.99</span>
560
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">443.6</span>
561
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">249.7</span>
562
|}
563
564
565
<span style="text-align: center; font-size: 75%;">Based on the structural values derived from the optimization objective of minimizing deformation, in order to check the accuracy of the optimization results (incorporating stress, deformation and mass) predicted by the two surrogate models for the values calculated by the optimization algorithms, the optimization results were brought into the simulation model to obtain the real results and the errors were obtained by comparing them with the predicted values of the models, which are shown in Tab. 5, Tab. 6 and Tab. 7, respectively.</span>
566
567
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
568
<span style="text-align: center; font-size: 75%;">Table 5 Deformation results of before and after optimization</span></div>
569
570
{| style="width: 100%;border-collapse: collapse;" 
571
|-
572
|  rowspan='1' style="border-top: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Optimal Model</span>
573
|  colspan='3'  style="border-top: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">After Optimization</span>
574
|  rowspan='1' style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Before Optimization</span>
575
|  rowspan='1' style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Decrease</span>
576
|-
577
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Prediction</span>
578
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">True</span>
579
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Error</span>
580
|-
581
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">E-AHF</span>
582
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.99</span>
583
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.98</span>
584
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.5%</span>
585
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.39</span>
586
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">17.2%</span>
587
|-
588
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">KRG</span>
589
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.98</span>
590
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.99</span>
591
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.5%</span>
592
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.39</span>
593
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">16.7%</span>
594
|}
595
596
597
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
598
<span style="text-align: center; font-size: 75%;">Table 6 Stress results of before and after optimization</span></div>
599
600
{| style="width: 100%;border-collapse: collapse;" 
601
|-
602
|  rowspan='1' style="border-top: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Optimal Model</span>
603
|  colspan='3'  style="border-top: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">After Optimization</span>
604
|  rowspan='1' style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Before Optimization</span>
605
|  rowspan='1' style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Decrease</span>
606
|-
607
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Prediction</span>
608
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">True</span>
609
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Error</span>
610
|-
611
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">E-AHF</span>
612
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">394.3</span>
613
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">407.7</span>
614
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.3%</span>
615
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">648.1</span>
616
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">37.1%</span>
617
|-
618
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">KRG</span>
619
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">430.5</span>
620
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">443.6</span>
621
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.0%</span>
622
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">648.1</span>
623
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">31.6%</span>
624
|}
625
626
627
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
628
<span style="text-align: center; font-size: 75%;">Table 7 Quality results of before and after optimization</span></div>
629
630
{| style="width: 100%;border-collapse: collapse;" 
631
|-
632
|  rowspan='1' style="border-top: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Optimal Model</span>
633
|  colspan='3'  style="border-top: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">After Optimization</span>
634
|  rowspan='1' style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Before Optimization</span>
635
|  rowspan='1' style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Decrease</span>
636
|-
637
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Prediction</span>
638
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">True</span>
639
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Error</span>
640
|-
641
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">E-AHF</span>
642
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">252.7</span>
643
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">249.7</span>
644
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.2%</span>
645
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">252.7</span>
646
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.2%</span>
647
|-
648
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">KRG</span>
649
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">252.7</span>
650
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">249.7</span>
651
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.2%</span>
652
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">252.7</span>
653
|  style="border-bottom: 2pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.2%</span>
654
|}
655
656
657
<span style="text-align: center; font-size: 75%;">From the prediction of the optimized results, the errors of the pallet rack shape variations of the E-AHF and KRG models are 0.5% and 0.5%, respectively; the errors of the stresses of the E-AHF and KRG models are 3.3% and 3.0%, respectively; and the errors of the pallet rack masses of the E-AHF and KRG models are 1.2% and 1.2%, respectively. Taken together, from the analysis of the accuracy of model prediction, the KRG model undoubtedly has better adaptability and stronger performance. From the comparison between before and after optimization, the reduction of pallet rack deformations by E-AHF and KRG models is 17.2%; 16.7%; the reduction of pallet rack stresses by E-AHF and KRG models is 37.1%; 31.6%; and the reduction of pallet rack mass by E-AHF and KRG models are 1.2%; 1.2%, respectively. Taken together, the optimization design based on the E-AHF model is more excellent and adaptable to engineering examples. Therefore, the optimal solution obtained from the E-AHF model is brought into the simulation model to obtain the optimized results, as shown in Fig. 9.</span>
658
659
{| style="width: 100%;border-collapse: collapse;" 
660
|-
661
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_120867196377-image14.jpeg|294px]] </span>
662
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_120867196377-image15.jpeg|294px]] </span>
663
|-
664
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">(a) Deformation after optimization</span>
665
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">(b) Stress after optimization</span>
666
|-
667
|  colspan='2'  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">Fig. 9 Optimized cloud maps</span>
668
|}
669
670
671
===5.2 Global sensitivity analysis===
672
673
Global sensitivity analysis (GSA) is the study of the effects of simultaneous changes in different design variables throughout the design space on a system or model, with the aim of implementing dimensionality reduction and simplifying the model in high-dimensional cases. The most representative of GSA, Sobol global sensitivity analysis, which is a standard deviation-based sensitivity analysis method, is used in this study [21]. This method can simultaneously calculate the first-order sensitivity index and global sensitivity index of design variables, in which the first-order sensitivity index indicates the impact of individual design variables, while the global sensitivity index takes into account the interactions between individual variables. Its specific calculation method is as follows.
674
675
Given any product function in n-dimensional space, it can be expressed as:
676
677
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
678
|-
679
| 
680
{| style="text-align: center;margin:auto;width: 100%;"
681
|-
682
| <math>f(x)={f}_{0}+\sum _{i}^{}{f}_{i}({x}_{i})+\sum _{i<j}^{}{f}_{ij}({x}_{i},{x}_{j})+</math><math>{f}_{1,2,\ldots ,n}({x}_{1},{x}_{2},\ldots ,{x}_{n})</math>
683
|}
684
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(12)</span>
685
|}
686
687
688
<span style="text-align: center; font-size: 75%;">where </span> <math display="inline">{x}_{i}</math><span style="text-align: center; font-size: 75%;"> denotes the i-th variable and </span> <math display="inline">{x}_{j}</math><span style="text-align: center; font-size: 75%;"> denotes the j-th variable.</span>
689
690
<span style="text-align: center; font-size: 75%;">The total standard deviation and partial standard deviation of </span> <math display="inline">f(x)</math><span style="text-align: center; font-size: 75%;"> are, respectively:</span>
691
692
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
693
|-
694
| 
695
{| style="text-align: center;margin:auto;width: 100%;"
696
|-
697
| <math>D=\int_{}^{}{f}^{2}(x)dx-{f}_{0}^{2}</math>
698
|}
699
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(13)</span>
700
|-
701
| 
702
{| style="text-align: center;margin:auto;width: 100%;"
703
|-
704
| <math>{D}_{{i}_{1}\ldots {i}_{s}}=\int_{}^{}{f}_{{i}_{1}\ldots {i}_{s}}^{2}(x)d{x}_{{i}_{1}}\ldots d{x}_{{i}_{s}}</math>
705
|}
706
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(14)</span>
707
|}
708
709
710
<span style="text-align: center; font-size: 75%;">The first-order sensitivity indices for each variable are:</span>
711
712
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
713
|-
714
| 
715
{| style="text-align: center;margin:auto;width: 100%;"
716
|-
717
| <math>{S}_{i}=\frac{{D}_{{i}_{1}\ldots {i}_{s}}}{D}</math>
718
|}
719
|  style="text-align: center;width: 5px;text-align: right;white-space: nowrap;"|<span style="text-align: center; font-size: 75%;">(15)</span>
720
|}
721
722
723
<span style="text-align: center; font-size: 75%;">where </span> <math display="inline">{S}_{i}</math><span style="text-align: center; font-size: 75%;"> denotes the first-order sensitivity index of the i-th variable. The global sensitivity index </span> <math display="inline">{S}_{ti}</math><span style="text-align: center; font-size: 75%;"> of the i-th variable is then the sum of the sensitivity indices of each order of the i-th variable.</span>
724
725
<span style="text-align: center; font-size: 75%;">Based on the E-AHF model, the GSA analysis of the pallet rack is carried out, and the results are shown in Fig. 10. It can be clearly seen that the channel height </span> <math display="inline">{x}_{1}</math><span style="text-align: center; font-size: 75%;"> has the greatest influence on the shape variables of the pallet rack, and both the first-order sensitivity index and the global sensitivity index are much higher than the remaining four design variables; whereas, the one that has a greater influence on the stress of the pallet rack is the channel flange thickness </span> <math display="inline">{x}_{4}</math><span style="text-align: center; font-size: 75%;">, and its sensitivity index is also much higher than that of the other four variables; whereas, for the pallet rack mass, it can be seen that, except for the web thickness </span> <math display="inline">{x}_{3}</math><span style="text-align: center; font-size: 75%;">, the rest of the design variables all have a greater influence. In summary, in the design of a pallet rack, the highest priority should be given to the height of channel steel </span> <math display="inline">{x}_{1}</math><span style="text-align: center; font-size: 75%;"> which has the greatest influence on its shape variables, to ensure the overall safety and pressure resistance of the lifting platform.</span>
726
727
{| style="width: 100%;border-collapse: collapse;" 
728
|-
729
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_120867196377-image16-c.png|180px]] </span>
730
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_120867196377-image17-c.png|174px]] </span>
731
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;"> [[Image:Review_120867196377-image18-c.png|180px]] </span>
732
|-
733
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">(1) GSA for deformation</span>
734
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">(2) GSA for stress</span>
735
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">(3) GSA for quality</span>
736
|-
737
|  colspan='3'  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Fig. 9 GSA based on E-AHF model</span>
738
|}
739
740
741
==6. Conclusion==
742
743
Taking the pallet rack of the lifting platform of the unmanned AGV stacker truck as the research object, the parametric design of the pallet rack is carried out by using SolidWorks and ANSYS, and then the design variables are selected, the design space is determined, and 60 sample points are selected by using the LHS method, of which 10 are used as a group for constructing an ensemble surrogate model represented by the E-AHF and a single surrogate model represented by KRG, the remaining 50 sample points are used as a group as a sample point library for the EI plus point strategy to improve the accuracy of the model by EI guidance, and the accuracy of two evaluation criteria can be tested. The remaining 50 sample points are used as the sample point bank for the EI infilling strategy, which is used to improve the accuracy of the model through the EI guidance, and at the same time, the accuracy of the constructed surrogate model can be examined by comparing the two evaluation criteria, which include the global and local indicators. Stress, mass and deformation are modelled using the two surrogate models, and two unused optimization problems with stress and deformation as constraints, and stress and mass as constraints are established, and Matlab Optimization Toolbox is used to optimize the models and compare the results of the two optimization problems, which show that the optimized mass has been reduced by 1.2%, and the deformation variable has been reduced by 17.2% compared with the initial deformation variable, which fully meets the objectives of the design, It fully meets the objectives of the optimized design. Finally, the E-AHF model is used as the basis for the global sensitivity analysis of the five design variables, and it is found that the most important shape variable is the length of the channel, which highlights the application of the surrogate model in the global sensitivity analysis, and at the same time provides a new research solution for the optimization of the forklift truck lifting platform.
744
745
==Acknowledgment==
746
747
The research was supported by the Longgang Institute of Zhejiang Sci-Tech University(LGYJY2021004).
748
749
==References==
750
751
[1] Chen, G.. Multi-objective optimization method based onmetamodeland its applications in vehicle body design [D].Hu’nan:Hu’nan University,2012 (6)
752
753
[2] Queipo, N.V., Haftka, R.T., Shyy, W. et al.. Surrogate-based analysis and optimization [J]. Progress in Aerospace Sciences,2005,41(1):1-28.
754
755
[3] Song, X., Sun, G. , Li, G. , Gao, W. , & Li, Q. . (2013). Crashworthiness optimization of foam-filled tapered thin-walled structure using multiple surrogate models. Structural and Multidisciplinary Optimization, 47(2), 221-231.
756
757
[4] Chen, W., Zhou, X., Wang, H. et al.. Multi-objective optimal approach for injection molding based on surrogate model and particle swarm optimization algorithm [J]. Journal of Shanghai Jiaotong University,2010,15(1): 88-93.
758
759
[5] Jiang, P., Zhou, Q., Shao, X.. Surrogate model-based engineering design and optimization [M]. Singapore : Springer,2020
760
761
[6] Ruan, X., Jiang, P., Zhou, Q. et al.. Variable-fidelity probability of improvement method for efficient global optimization of expensive black-box problems7 [J]. Structural and Multidisciplinary Optimization,2020, 62(6):3021-3052.
762
763
[7] Meng, Z., Zhang, Z., Zhou, H. et al.. Robust design optimization of imperfect stiffened shells using an active learning method and a hybrid surrogate model [J]. Engineering Optimization,2020,52(12):2044-2061.
764
765
[8] Yondo, R., Andres, E., Valero, E.. A review on design of experiments and surrogate models in aircraft real-time and many-query aerodynamic analyses [J]. Progress in Aerospace Sciences,2018,96:23-61.
766
767
[9] Zhang, X., Xie, F., Ji, T., Zhu, Z., Zheng, Y .. (2021). Multi-fidelity deep neural network surrogate model for aerodynamic shape optimization. Computer Methods in Applied Mechanics and Engineering, 373, 113485.
768
769
[10] Kushner, H. J.. A new method of locating the maximum point of an arbitrary multipeak curve in the  presence of noise [J]. Journal of Basic Engineering, 1964, 86(1): 97-106.
770
771
[11]Jones, D. R., Schonlau, M., Welch, W. J.. Efficient global optimization of expensive black-box functions [J]. Journal of Global Optimization, 1998, 13: 455-492.
772
773
[12] Zeng, F., Zhou, J.. Optimization strategy for dynamic metamodel integrating minimize lower confidence bound and trust region [J]. Journal of Mechanical Engineering,2017,53(13):170-178
774
775
[13] Kleijnen, J. P. C.. Expected improvement in efficient global optimization through bootstrapped kriging [J].Journal of Global Optimization,2012,54(1):59-73.
776
777
[14] Chaudhuria, A., Haftka, R. T. et al.. Experimental flapping wing optimization and uncertainty quantification using limited samples [J]. Structural & Multidisciplinary Optimization,2015,51(4):1-14.
778
779
[15] Song, X., Lv, L., LI, J., Sun, W., Zhang, J.. 2018. An Advanced and Robust Ensemble Surrogate Model: Extended Adaptive Hybrid Functions. Journal of Mechanical Design, 140.
780
781
[16] Mckay, M. D., Beckman, R. J., Conover, W. J.. A comparison of three methods for selecting values of input  variables in the analysis of output from a computer code [J]. Technometrics, 1979, 21(2): 239-245.
782
783
[17] GroempoIng, U.. 2014, “Fractional Factorial Design,” Wiley Interdiscipl. Rev.Comput. Stat., 1(2), pp. 234–244.
784
785
[18] Boddy, R., Smith, G.. 2010, Effective Experimentation: For Scientists and Technologists, John Wiley & Sons, New York.
786
787
[19] Turner, C. J., Crawford, R. H., Campbell, M. I.. Multidimensional sequential sampling for NURBs-based  metamodel development [J]. Engineering with Computers, 2007, 23: 155-174
788
789
[20] Cox, D. D., John, S.. 1992, “SDO: A Statistical Method for Global Optimization,” IEEE International Conference on Systems, Man and Cybernetics, Vol. 2, Chicago, IL, Oct. 18–21, IEEE Xplore, pp. 1241–1246.
790
791
[21] Sobol′, I. M.. 2001. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Mathematics and Computers in Simulation, 55, 271-280.
792

Return to Zhang et al 2023i.

Back to Top

Document information

Published on 27/10/23
Accepted on 17/10/23
Submitted on 14/09/23

Volume 39, Issue 4, 2023
DOI: 10.23967/j.rimni.2023.10.005
Licence: CC BY-NC-SA license

Document Score

0

Views 120
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?