# Multiobjective Optimization for Fixture Locating Layout of Sheet Metal Part Using SVR and NSGA-II

Authors: Yuan Yang; Zhongqi Wang; Bo Yang; Zewang Jing; Yonggang Kang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Fixture plays a significant role in determining the sheet metal part (SMP) spatial position and restraining its excessive deformation in many manufacturing operations. However, it is still a difficult task to design and optimize SMP fixture locating layout at present because there exist multiple conflicting objectives and excessive computational cost of finite element analysis (FEA) during the optimization process. To this end, a new multiobjective optimization method for SMP fixture locating layout is proposed in this paper based on the support vector regression (SVR) surrogate model and the elitist nondominated sorting genetic algorithm (NSGA-II). By using ABAQUS™ Python script interface, a parametric FEA model is established. And the fixture locating layout is treated as design variables, while the overall deformation and maximum deformation of SMP under external forces are as the multiple objective functions. First, a limited number of training and testing samples are generated by combining Latin hypercube design (LHD) with FEA. Second, two SVR prediction models corresponding to the multiple objectives are established by learning from the limited training samples and are integrated as the multiobjective optimization surrogate model. Third, NSGA-II is applied to determine the Pareto optimal solutions of SMP fixture locating layout. Finally, a multiobjective optimization for fixture locating layout of an aircraft fuselage skin case is conducted to illustrate and verify the proposed method.


---

## Body

## 1. Introduction

Sheet metal part (SMP), having the advantages of high strength, light weight, and fluent geometric modeling, widely exists in various fields of aerospace, vehicle, and so on [1, 2]. However, it is easy to deform during its assembly, machining, and inspection processes because of the properties of thin wall, large size, and low rigidity. Recently, flexible fixture is widely used during the whole manufacturing process to locate and constrain SMP accurately to reduce its dimensional and form errors. Moreover, the “N-2-1” (N>3) locating principle, originally proposed by Cai et al. [3], is extensively recognized in the fixture design to not only ensure the locating deterministic but also restrain the excessive deformation of SMP. However, the key to the application of “N-2-1” principle in flexible fixture design is how to find the optimal fixture locating layout to minimize the SMP deformation.To determine the optimal SMP fixture locating layout effectively and efficiently, many scholars and technicians have conducted a lot of research. In the beginning, finite element analysis (FEA) and evolutionary algorithm were integrated to model the SMP fixture system and search for the optimal SMP fixture layout. Kaya [4] proposed a combined GA and FEA approach to optimizing the fixture layout to decrease the workpiece deformation during its machining process. Prabhaharan et al. [5] used FEA to compute the workpiece deformation and applied ant colony algorithm (ACA) for the workpiece fixture layout optimization to minimize its dimensional and form errors. Dou et al. [6] established an FEA model implemented by ANSYS to calculate the elastic deformation of the part and employed four evolutionary approaches including GA, PSO, improved GA, and improved PSO for fixture layout optimization. Cheng et al. [7] presented a hierarchical fixture layout model to optimize the locating points of aeronautical thin-walled structure to decrease the assembly variation by genetic algorithm and ants algorithm. Xiong et al. [8] put forward an “N-2-1-1” locating principle and then solved the fixture layout optimization of a compliant aerospace part by integrating FEA and GA.We can see that the optimization method integrating FEA with evolutionary algorithm has become the main popular method to optimize the fixture locating layout of SMP. However, since there is no explicit functional relation between fixture layout scheme and the corresponding SMP deformation, the above optimization methods for SMP fixture layout often involve a lot of time-consuming FEA which results in unbearable computational costs and burdens. Therefore, the surrogate or prediction models, such as artificial neural network (ANN) and response surface methodology (RSM), were developed and applied to approximate the mapping relation between the design variables and objective values so as to improve the optimization efficiency and to decrease the computational cost. Li et al. [9] developed three different design models of an ANN-enhanced RSM, a nonlinear programming model, and a polynomial RSM to assist optimization and design of fixture layout during the SMP spot welding assembly. Rex and Ravindran [10] also established an ANN-based model to approximate the workpiece fixture system and determine the optimal fixture layout to reduce the maximum elastic deformation by incorporating full factorial design of experiments (DOE) and statistical analysis. Lu and Zhao [11] built a back propagation neural network (BPNN) so as to approximate the deformation of SMP under a given fixture layout and then employed GA to find the optimal fixture layout based on the BPNN prediction model. Wang et al. [12, 13] proposed a radial basis function neural network (RBFNN) prediction model to predict the deformation of SMP and then carried out the follow-up work to search for the optimal fixture layout by integrating RBFNN and bat algorithm to improve the location quality and optimization efficiency. Furthermore, Yang et al. [14] presented an integrated method to determine the optimum fixture locating layout to minimize the overall deformation of SMP by combining kriging with cuckoo search algorithm.All researches above only involve single objective of fixture layout optimization. However, in most cases, the fixture layout design should satisfy the requirements of multiple objectives. So as to decrease the workpiece deformation and improve the deformation distributing uniformity, Chen et al. [15] put forward a multiobjective method for fixture design optimization, where FEA was used to calculate the deformation and GA was employed to solve the optimization model. By regarding fixture locating layout as a multiobjective optimization problem, Wang and Nee [16] applied nondominated GA and ACA to solve the multiobjective problem with the objective of maximizing the distance of the locators and minimizing the locating error and the deviation of the locating error. Liu et al. [17] developed an accelerated FEA model to calculate the nodal displacements of the workpiece under external forces and proposed an optimization method to find the optimal fixture layout based on the multiobjective GA. Considering multiple requirements for the locating layout of a checking fixture such as detachability, stability, and accuracy, Jiang et al. [18] presented a multiobjective optimization method to continuous search for a detachable, stable, and accurate locating layout in checking fixture design. By combining nondominated sorting algorithm and the social radiation algorithm, a nondominated sorting social radiation algorithm was proposed by Xing et al. [19] to conduct the multiobjective optimization of compliant parts fixture scheme to simultaneously meet the multiple requirements of gravity deformation and assembly tolerance.In this paper, a new multiobjective optimization method by integrating support vector regression (SVR) with the elitist nondominated sorting genetic algorithm (NSGA-II) is proposed for the SMP fixture locating layout design. The objective here is twofold. One is to minimize the overall deformation of SMP, and the other is to minimize the maximum deformation. A parametric FEA model implemented by Python scripting interface of ABAQUS is developed to calculate the deformation for a given fixture layout under given forces. Then NSGA-II is integrated, and the NSGA-II toolbox based on MATLAB™ is developed to solve the optimization problem. The remainder of this paper is organized as follows. The mathematical model for SMP fixture layout multiobjective optimization is represented in Section2. The fundamental theory related to SVR is introduced in Section 3. Section 4 describes the construction flowchart for the multiobjective optimization of SMP fixture layout using SVR and NSGA-II. In Section 5, an experiment is conducted to verify the proposed optimization method by comparing the simulated results with the experimental data. Finally, the major conclusions are provided in Section 6.


## 2. Problem Formulation

Since it is easy to deform in many manufacturing operations, SMP is often located under an overconstraint condition in the flexible fixture based on the “N-2-1” (N>3) locating principle. Figure 1 depicts the typical “N-2-1” (N=4) fixture locating scheme for a curved SMP. It can be seen thatN locating points (LPs), generally more than three, are required on the primary datum to prevent excessive deformation and supply more reinforcements for SMP. And two LPs on the secondary datum and one on the tertiary datum are expected to locate SMP uniquely in space. Obviously, the locating effectiveness based on “N-2-1” principle is strongly related to the layout of theN LPs where the numberN is determined by the SMP dimensional specifications. In this paper, our interests focus on the flexible fixture layout optimization of SMP on the primary datum to minimize its overall deformation and maximum deformation synchronously. Therefore, the fixture layout design is regarded as a multiobjective optimization problem in this work.Figure 1
The “N-2-1” locating principle for a curved SMP.In SMP fixture layout optimization problem, two goals (the overall deformation accuracy and the maximum deformation accuracy) should be synchronously satisfied to obtain an optimal solution. In this paper, FEA method is employed to model the SMP fixture system and to compute the objective values. Here, in order to minimize the human assistance and realize the automation of the whole iterative optimization and FEA, this paper performs the secondary development of ABAQUS and establishes a parametric FEA model for SMP fixture layout optimization. Thus, the multiobjective optimization problem can be formulated as(1)Find:X=x1,x2,…,xi,…,xj,…,xNMinimize:F=F1,F2F1X=∑i=1Lwi2XLF2X=maxwiXSubject  to:xi∈S,xj∈S,xi≠xj,where X represents the design variables on behalf of various fixture layout schemes; xi and xj represent the coordinates of any two different LPs, respectively, where i,j=1,…,N; F1 and F2 denote the two goals, the overall deformation, and maximum deformation, respectively; wi is the displacement deviation of theith finite element node after deformation;L depicts the total number of the SMP finite element nodes;S represents the set of all finite element nodes. Here, the design variableX has to be in the predetermined domainS, and in each fixture layout scheme, any two LPs cannot overlap.


## 3. SVR

As a powerful machine learning method, support vector machine (SVM) was initially presented by Vapnik [20, 21] and applied to solve regression and classification problems having small samples and high dimensions on the base of structural risk minimization principle and statistical learning theory. By introducing Vapnik’s ε-insensitive loss function, the regression version of SVM, named SVR, has been used as a powerful and effective method to solve the nonlinear regression problems.Given a training setT=Xj,Fj∣j=1,2,…,l of size l, where Xj are the input vectors (representing fixture locating layout schemes) and Fj=Fj1,Fj2 are the associated output values (representing the SMP overall deformation Fj1 and maximum deformation Fj2) of Xj. The aim of SVR is to provide a nonlinear mapping function to map the training data to a higher dimensional feature space. Now, the nonlinear mapping relation between the input data and the output values can be expressed as follows:(2)F^X=WTφX+b,W,X∈RN, where φX is the feature that is nonlinearly mapped from the input vectorX;W andb are the coefficients to be adjusted. Moreover, Vapnik’sε-insensitive loss function [20] is employed in this paper to panelize the training error between the real training value Fj and the regression value F^j.As displayed in Figure2, thisε-insensitive loss function formulizes a tube with a radius ofε around the assumptive regression function so that if a training sample point lies within the tube, the loss function equals zero, while if a training sample point falls on or outside the tube, the loss is in proportion to the magnitude of the Euclidean difference between the training sample point and the radiusε of the tube. Thus, the SVR regression problem can be solved by calculating the following constrained optimization function:(3)minimize:12WTW+C∑j=1lξj+ξj∗subject  to:Fj-WTφXj+b≤ε+ξjWTφXj+b-Fj≤ε+ξj∗ε,ξi,ξj∗≥0,j=1,2,…,l.Figure 2
SVR and Vapnik’sε-insensitive loss function.Here, the penalty factorC>0 defines the trade-off between the model complexity of F^X=WTφX+b and the amount up to which deviations larger than the radiusεcan be tolerated. ξi and ξi∗ denote the nonnegative slack variables measuring the error of the up and down sides, separately. Then, by solving the optimization problem above and introducing Lagrange multipliers, the coefficientW of (1) can be obtained by(4)W=∑j=1lαj-αj∗φXj, where αi and αi∗0≤αi,αi∗≤C are the Lagrange multipliers. Then the SVR regression function can be reformulated to a dual equation:(5)F^X=∑j=1lαj-αj∗φTXjφX+b=∑j=1lαj-αj∗KXj,X+b.Here,KXj,X stands for the so-called kernel function. In this paper, the most commonly used kernel function, namely, the Gaussian radial basis function (RBF), is adopted in the SVR modeling. The RBF kernel function can be expressed as(6)KXj,X=exp⁡-γXj-X2, where γ denotes the width of the Gaussian RBF function; Xj represents thejth input of the training sample set, and X represents the input of the prediction data. Here, the coefficient of aj, aj∗, andb can be obtained by several iterations of the training sample set. Now, by tuning the penalty factorC, width of the RBF kernel function γ, and tube size ofε-insensitive loss functionε, the SVR model can be constructed to deal with the optimization problem of the SMP fixture locating layout.


## 4. Multiobjective Optimization Based on SVR and NSGA-II

Multiobjective optimization always typically involves multiple conflicting objectives and does not have a feasible solution that minimizes or maximizes all objectives synchronously. Hence, the final outcome of such a multiobjective optimization problem is a set of trade-off solutions between different objectives. These trade-off points are termed as Pareto optimal solutions which are not dominated by any other solution and cannot be improved in the case of at least one other objective without worsening. The set of all these feasible nondominated solutions is termed as Pareto optimal solution set, and the corresponding objective values are named the Pareto front [22]. Figure 3 demonstrates the Pareto front for the optimization of two objectives of Obj-1 and Obj-2.Figure 3
Pareto front for a two-objective optimization.Recently, it has been found that NSGA-II can converge near to the real Pareto front and keep a good diversity of the population in the final Pareto optimal solution set [23, 24]. In this work, the Pareto optimal solutions and the corresponding Pareto front are obtained by NSGA-II which was originally proposed by Deb et al. [24] based on the nondominated sorting genetic algorithm (NSGA) [25]. The three main features of NSGA-II are elitist strategy, crowding distance, and fast ranking of nondominated solutions. Therefore, so as to verify the feasibility and effectiveness of the proposed multiobjective optimization method integrating SVR and NSGA-II, based on the ABAQUS [26] secondary development, this paper uses Python language to compile FEA model of the parametric SMP fixture system to generate the limited training and testing samples for SVR surrogate models and finally uses MATLAB language to compile the NSGA-II program to do the multiobjective optimization of SMP fixture layout.ABAQUS can execute the finite element calculation by running a script file with the  .PY extension and store the postprocessing results in a database file with the  .ODB extension. MATLAB can read and modify the PY and ODB files to achieve the connection and data transmission between MATLAB language and ABAQUS software. Thus, the information from optimization algorithm can be sent to finite element model by modifying the PY files, and the displacement can be received from ABAQUS by reading the ODB files. The main steps of the multiobjective optimization method are described as follows.Step 1.
Generate the limited training and testing samples by Latin hypercube design (LHD) [27] and then calculate the responding overall deformation and maximum deformation of the SMP under external forces with the help of the parametric FEA model.Step 2.
Scale the training data set and testing data set into the range of0,1 using the following normalization formula:(7)xi′=xi-xminxmax-xmin,where xi is theith original data value, xi′ is theith scaled data value, and xmax and xmin are the upper and lower bound of the original data, respectively.Step 3.
Tune the SVR hyperparameters (C, γ and ε) and establish the SVR prediction models to depict the nonlinear mapping relation between different fixture layout and the corresponding SMP deformation.Step 4.
Calculate and analyze the relative root mean squared error (RRMSE) of the SVR models on the testing sample set. Then increase training samples and repeat Steps1–4 until the prediction accuracy is satisfied.Step 5.
Initialize the parent populationP0 of NSGA-II randomly and evaluate the population by calculating the objectives of each individual using the established SVR prediction models.Step 6.
Perform the nondominated sorting onP0 and rank each individual based on the dominance criteria.Step 7.
Compute the crowding distances and sort the individuals.Step 8.
Select individuals using a tournament selection and store them in a temporary mating pool. Then apply the crossover and mutation operators to produce the offspring populationOt where subscriptt denotes the count of current generation.Step 9.
Merge the parent populationPt and offspring population Ot to generate a larger population and assign a fitness value to each individual by a fast nondominated sorting.Step 10.
Perform elitism to find the individuals with better fitness from the merged population as the new parent populationPt+1.Step 11.
Repeat Steps8–10 until the maximum number of iterations is reached (t=tmax).When the optimization process terminates, the nondominated solutions of the final population are the approximate Pareto optimal solutions with a good diversity. Figure4 explains the implementation of the proposed multiobjective optimization of SMP fixture layout.Figure 4
Flowchart for multiobjective optimization of SMP fixture layout using SVR and NSGA-II.


## 5. Case Study

In this section, a simplified aircraft fuselage skin is chosen as an example to illustrate and evaluate the multiobjective optimization method for SMP fixture locating layout by combining SVR and NSGA-II. The basic dimension, force condition, and initial fixture layout of the sheet metal fuselage skin are demonstrated in Figure5, and the physical properties of the material are listed in Table 1. “N” (N=4) LPs on the primary datum plane of the fuselage skin are LP-1, LP-2, LP-3, and LP-4, and the “2” LPs on the secondary datum are LP-5 and LP-6, while the “1” LP on the tertiary datum is LP-7. The coordinates of the fixed locating points (FLPs) of LP-1, LP-2, LP-3, LP-5, LP-6, and LP-7 are preset as (900,1500),(300,1500),(300,500),(0,1500),(0,500), and (600,2000), respectively. The LP to be optimized (OLP) is LP-4 and its coordinate is expressed by (x,y). The sheet metal fuselage skin is subjected to its own gravity (G) and a concentrated external force (F).Table 1
The material physical properties.
Material properties
Values
Mass density
2.8
×
10
3 kg/m3
Poisson ratio
0.33
Young’s modulus
7.12
×
10
4 MPaFigure 5
The initial “N-2-1” fixture layout for the simplified aircraft fuselage skin.The Pareto optimal position of LP-4 is to be determined by the method proposed in this paper by combining SVR and NSGA-II for a trade-off between the two goals, the minimum overall deformation and the minimum maximum deformation of the aluminum alloy sheet metal fuselage skin. According to the workflow in Figure4, the training and testing sample points are produced by LHD and the two-objective function values of the fuselage skin under its dead weight and a concentrated force is calculated by ABAQUS. In this work, after several trials and a normalization of the sample data sets, two SVR models are built with 800 sets of known input and output values for training and 40 sets of known values for testing. Figure 6 depicts the comparison between the SVR output and the expected FEA output for the same testing samples. As a result, the RRMSE of the SVR prediction models for the overall deformation and maximum deformation are 3.68% and 3.55%, respectively; both are less than the need of general engineering of 10%. The parameter setting and the RRMSE of the two SVR models are listed in Table 2. Therefore, the response surface models depicting the mapping relation between fixture layout and the corresponding SMP deformation are built, as demonstrated in Figure 7. Now, for a given fixture locating layout, the SMP deformation (the overall deformation and the maximum deformation) can be obtained.Table 2
Experimental results of SVR prediction models.
Penalty factor (C)
Width of RBF kernel function (γ)
Tube size ofε-insensitive loss function (ε)
RRMSE
SVR model forF1
6.7
213
0.0171
3.68%
SVR model forF2
17
164.5
0.02
3.55%Figure 6
Output curves of SVR prediction models for the overall and maximum deformation.
(a)
(b)Figure 7
Response surfaces of SVR prediction models for the overall and maximum deformation.
(a)
(b)The crossover and mutation probabilities of NSGA-II used in this example are 0.9 and 0.1, respectively. The Pareto front by NSGA-II is illustrated in Figure8, with 100 population individuals and the maximum number of iterations is 1000. (The coordinate values and the corresponding overall and maximum deformation of the trade-off points on the Pareto front are given in Table 4.)Figure 8
Pareto optimal solutions by NSGA-II.In order to further verify the feasibility and effectiveness of the proposed multiobjective optimization method, the multipoint flexible fixture system, as shown in Figure9, was designed, and the validation experiment for the Pareto optimal solutions of the sheet metal fuselage skin was conducted. The experimental data were measured to assess the validity of the simulated results. In this experiment, the high precision coordinate measurement system of API Laser Tracker 3 was used to determine the locators’ spatial position of the flexible fixture system, and the HandySCAN 3D™ handheld scanner was employed to measure the deformation condition of the sheet metal fuselage skin after being located, as shown in Figure 10.Figure 9
Multipoint flexible locating fixture system.Figure 10
SMP fixture system and the measurement operation.In this paper, the two ends (numbers 1 and 100 in Table4) of the Pareto front and an arbitrary point (number 57) in the middle were taken for experimental verification. The experimental comparisons are given in Table 3.Table 3
Comparison of simulated results with experimental data.
Number
Coordination
Simulated results
Experimental results
Relative errors
F
1/mm
F
2/mm
F
1/mm
F
2/mm
F
1
F
2
1
(857.1263, 640.515)
1.727699
6.430427
1.6582
6.125
4.1895%
4.9866%
57
(871.6877, 641.0927)
1.75448
6.208433
1.9243
6.562
8.825%
5.3881%
100
(883.6917, 651.7052)
1.819898
6.161721
1.978
6.443
7.9917%
4.3657%Table 4
Pareto optimal solution set.
Number
Coordination
F/mm
F
1
F
2
1
(857.1263, 640.515)
1.727699
6.430427
2
(857.4084, 640.3446)
1.72771
6.423829
3
(857.558, 640.5037)
1.727724
6.420588
4
(857.9121, 640.012)
1.727788
6.41228
5
(858.0873, 640.1298)
1.727824
6.408521
6
(858.2955, 640.4132)
1.727882
6.404211
7
(858.4601, 640.1936)
1.727936
6.400485
8
(858.6489, 639.8772)
1.728013
6.39622
9
(858.8504, 640.4498)
1.728097
6.39238
10
(858.9442, 639.7643)
1.728146
6.389922
11
(859.1689, 639.9301)
1.728254
6.385363
12
(859.2113, 640.0448)
1.728275
6.384564
13
(859.4965, 639.8107)
1.728445
6.378577
14
(859.6698, 639.8657)
1.728555
6.375112
15
(859.8157, 640.0276)
1.728655
6.372288
16
(859.9493, 639.4548)
1.728765
6.369304
17
(860.1582, 640.032)
1.728912
6.365522
18
(860.4934, 639.3426)
1.729205
6.358663
19
(860.6067, 638.5751)
1.729375
6.356181
20
(860.8294, 638.8946)
1.729543
6.352104
21
(861.2121, 638.8173)
1.729931
6.344983
22
(861.3658, 638.5834)
1.730117
6.342108
23
(861.8739, 639.0724)
1.73066
6.333174
24
(862.0308, 639.0256)
1.730856
6.330407
25
(862.1956, 639.5911)
1.731055
6.327726
26
(862.8157, 639.1991)
1.731914
6.317078
27
(863.0007, 639.3575)
1.732186
6.314057
28
(863.2695, 639.3298)
1.732599
6.309667
29
(863.5583, 640.197)
1.733112
6.305263
30
(863.66, 640.1977)
1.733282
6.303648
31
(863.9331, 639.6238)
1.733705
6.299226
32
(864.1062, 639.6468)
1.734011
6.296564
33
(864.293, 639.0971)
1.734331
6.293644
34
(864.828, 639.2486)
1.735339
6.285737
35
(864.9587, 639.1982)
1.735594
6.283842
36
(865.2114, 640.2559)
1.736193
6.280348
37
(865.4428, 639.7651)
1.736616
6.277022
38
(865.7169, 639.772)
1.737202
6.273236
39
(865.7257, 639.9955)
1.737247
6.273134
40
(865.9755, 639.4373)
1.73774
6.269726
41
(866.0171, 639.4455)
1.737833
6.269171
42
(866.2835, 639.4555)
1.738436
6.265652
43
(866.4907, 638.627)
1.738879
6.263041
44
(866.5304, 638.8366)
1.738974
6.262505
45
(867.1737, 639.0451)
1.740534
6.254456
46
(867.3659, 639.0463)
1.741016
6.252137
47
(868.4293, 639.2489)
1.74385
6.239908
48
(868.4349, 639.2294)
1.743863
6.239851
49
(868.6918, 638.9492)
1.744555
6.237148
50
(868.9649, 638.973)
1.745335
6.234259
51
(869.5347, 639.8578)
1.747138
6.228198
52
(869.9662, 639.4221)
1.748383
6.224126
53
(870.3028, 639.4324)
1.74944
6.220968
54
(870.9122, 639.925)
1.751508
6.215329
55
(871.0192, 639.9366)
1.751865
6.214398
56
(871.1483, 640.2595)
1.752376
6.213178
57
(871.6877, 641.0927)
1.75448
6.208433
58
(872.2787, 640.9801)
1.75653
6.203838
59
(872.5152, 641.3372)
1.757522
6.201948
60
(872.6415, 641.3698)
1.757997
6.20101
61
(873.039, 641.3271)
1.759454
6.198197
62
(873.4614, 641.3901)
1.761079
6.195316
63
(873.9644, 641.7911)
1.763208
6.191943
64
(874.2811, 641.796)
1.764462
6.190021
65
(874.3631, 641.9308)
1.764855
6.189481
66
(874.7061, 640.6866)
1.765692
6.188154
67
(875.1304, 642.3306)
1.768186
6.185035
68
(875.3949, 642.0009)
1.769109
6.183821
69
(876.0423, 641.8134)
1.771753
6.18083
70
(876.1436, 641.8305)
1.772197
6.180367
71
(876.6645, 641.9238)
1.77451
6.178109
72
(876.9324, 641.9135)
1.775686
6.177062
73
(877.2183, 642.3166)
1.777186
6.175749
74
(877.2129, 644.6902)
1.77883
6.17469
75
(877.5359, 645.1398)
1.780681
6.173351
76
(877.7426, 644.8819)
1.781407
6.172697
77
(878.1743, 645.8126)
1.784253
6.170986
78
(878.4108, 646.2754)
1.78582
6.170133
79
(878.666, 646.2975)
1.787051
6.169366
80
(878.8959, 646.2787)
1.788129
6.168725
81
(879.0174, 646.8048)
1.789262
6.168275
82
(879.4558, 646.3675)
1.790925
6.167285
83
(879.6992, 646.4447)
1.792194
6.166714
84
(880.1559, 646.2866)
1.794284
6.165852
85
(880.6159, 646.3753)
1.796674
6.165043
86
(880.8007, 646.5716)
1.797813
6.164695
87
(880.6268, 648.522)
1.799213
6.164473
88
(881.1824, 646.3025)
1.799467
6.164314
89
(881.6068, 646.5678)
1.801923
6.163765
90
(881.4987, 648.1459)
1.803189
6.163329
91
(881.941, 648.5292)
1.805964
6.162812
92
(882.3562, 648.1199)
1.80762
6.162642
93
(882.6308, 648.0689)
1.809005
6.162536
94
(882.5089, 649.601)
1.81038
6.162209
95
(882.6627, 649.5111)
1.811069
6.162146
96
(882.9777, 649.8114)
1.813167
6.161979
97
(883.1503, 650.466)
1.815044
6.161847
98
(883.3552, 650.5182)
1.816223
6.161807
99
(883.6135, 650.6113)
1.817757
6.161787
100
(883.6917, 651.7052)
1.819898
6.161721It is known from Table3 that the maximum relative error of the simulated results is not more than the need of general engineering of 10% from the experimental data. The results show that the established FEA model can be efficiently applied to analyze the SMP overall deformation and maximum deformations, and the proposed multiobjective optimization method as well as the given Pareto front (see Figure 8) can be used to guide the designers to evaluate and determine the final optimal acceptable fixture layout design with a proper trade-off based on multiple objectives.


## 6. Conclusions

In order to realize the multiobjective optimization of SMP fixture locating layout and reduce the cost of FEA calculation during the iterative optimization, a new multiobjective optimization design method is proposed, which combines SVR surrogate model and NSGA-II in fixture layout design. The two major conclusions are as follows:(1)
On the basis of a few FEA calculations, the SVR prediction models of the SMP overall deformation and maximum deformation are constructed, respectively, which greatly reduce the cost of finite element calculation of fixture layout optimization and improve the fixture design effectiveness. At the same time, the RRMSE of the prediction data satisfies the general engineering requirement (10%), and the prediction speed is fast and the prediction result is stable.(2)
A new multiobjective optimization model for SMP fixture locating layout is established by combining SVR with NSGA-II. The Pareto distribution is obtained, which can provide a theoretical guidance for multiobjective optimization design of SMP fixture locating layout.



---

*Source: 7076143-2017-08-17.xml*