Multiobjective Optimization of Injection Molding Process Parameters for the Precision Manufacturing of Plastic Optical Lens

Injectionmolding process parameters (IMPP) have a significant effect on the optical performance and surface waviness of precision plastic optical lens.This paper presents a set of procedures for the optimization of IMPP, with haze ratio (HR) reflecting the optical performance and peak-to-valley 20 (PV20) reflecting the surface waviness as the optimization objectives. First, the orthogonal experiment was carried out with the Taguchi method, and the results were analyzed by ANOVA to screen out the IMPP having a significant effect on the objectives. Then, the 3 full-factor experiment was conducted on the key IMPP, and the experimental results were used as the training and testing samples. The BPNN algorithm and the M-SVR algorithm were applied to establish the mapping relationships between the IMPP andobjectives. Finally, themultiple-objective optimizationwas performedby applying the nondominated sorting genetic algorithm (NSGA-II), with the builtM-SVRmodels as the fitness function of the objectives, to obtain a Pareto-optimal set, which improved the quality of plastic optical lens comprehensively. Through the experimental verification on the optimization results, the mean prediction error (MPE) of HR and PV20 is 7.16% and 9.78%, respectively, indicating that the optimization method has high accuracy.


Introduction
In recent years, with the rapid development of the photovoltaic industry and information technology, the demand for high-quality optical components has shown a sharp upward trend.Optical lens is an indispensable component of these systems, so it is particularly important to carry out an indepth study on this aspect.The conventional optical lens processing materials are generally made of glass.However, with the development of new polymer materials with features such as low density, impact resistance, and low cost, as well as the production technology, more and more optical lenses are being made of special resins.The manufacturing technology of plastic optical lens generally involves injection molding process and compression methods.The injection molding process can produce complex products, with the advantages of high efficiency, low cost, and automation, so it is mainly used for producing plastic optical lens.Although the production of plastic optical lens has these advantages, compared with the glass optical lens, it also has some shortcomings due to constraints from material characteristics and the injection molding process, such as low light transmittance, birefringence, moisture absorption, poor heat resistance, and low geometry dimensional accuracy.In order to further enhance the quality of plastic optical lenses, more suitable new materials need to be developed and a new injection molding process needs to be studied.For a specific type of plastic, the quality of its injection molding products is affected by the injection machine, mold, and IMPP; the IMPP may be adjusted easily and therefore it has been widely investigated by researchers.
In the field of plastic optical lens processing, many researchers mainly adopt the geometric accuracy and defects of the optical lens as the research objectives, putting forward the corresponding optimization method and successfully improving the size, shape, and surface accuracy of the optical lens.Pazos et al. [1] estimated the ultimate thickness of PC biconcave and biconvex lenses using the CAE tool and evaluated the quality of the optical lens based on the number and location of weld lines and air-traps.J. Wen and P. Wen [2], Bensingh et al. [3], and Shieh et al. [4] optimized PMMA aspheric lens, bi-aspheric lens, and large diameter aspheric plastic lens, respectively, through design of experiments (DOE) and numerical simulation, reducing the volumetric shrinkage of the lens effectively.With the linear shrinkage of a specific point of the optical lens as an optimization objective, Hu and Cui [5] studied the effect of holding pressure and gate size on the objective.Lu and Khim [6] and Spina et al. [7] investigated some effects of the molding conditions on the surface contours and birefringence patterns of the lens, using a statistical experimental analysis method and grey relational component analysis coupled to the Taguchi design.Taking light transmission, surface waviness, and surface finish as the optimization objectives, Tsai et al. [8] established linear, exponential, and nonlinear regression models for the objectives and predicted the optimal IMPP; they concluded that the nonlinear regression model had higher prediction accuracy.They established an inverse model of injection molding for form accuracy of optical lens by combining the artificial neural network (ANN) with a genetic algorithm (GA) [9].They also investigated the correlation between melt pressure and cavity pressure in different runner positions [10].Other researchers, taking the optical and mechanical performance of the optical lens as the optimization objectives, proposed different methods to adjust the IMPP, reducing the birefringence and stress of the optical lens effectively.They mainly focused on the optimization of birefringence, used the numerical simulation method to study a PMMA CD-ROM lens [11], and verified the simulation results of a cyclic-olefincopolymers aspheric lens by combining with the Taguchi method [12].Through the adjustment of IMPP, the birefringence of the optical lens was effectively reduced.Using the shortening of the processing cycle time and reduction of optical retardation as the optimization objectives, Turng et al. [13] achieved the high-quality PMMA lens via CAE simulation and experimental verification and evaluated the generated stresses and birefringence.
As mentioned above, although some researchers studied multiple quality objectives in the existing optimization of the plastic optical lens processing, they separately studied each individual objective while determining the optical IMPP, without considering the mutual effect of multiple objectives.However, the overall quality of the optical lens is obviously determined by its geometric accuracy, optical performance, and mechanical performance.Only optimizing the IMPP of an individual objective may lead to large defects in other optical quality objectives.From the practical point of view, therefore, it is necessary to optimize multiple quality objectives of the optical lens at the same time.The multiobjective optimization problem (MOP) has been widely used in the optimization of injection molding process and machine.For injection molding machine, Kuo and Chang [14] proposed a turbo injection mode (TIM) for an axial flux motor, and the particle swarm optimization (PSO) was used to find out the multiobjective optimal design solution.The product quality, energy-saving efficiency, and molding efficiency were optimized, and the corresponding multiple-objective Pareto-optimal set was obtained.Zhao et al. [15] took warp, volumetric shrinkage, and sink marks as the optimization objectives and achieved the multiobjective optimization by establishing the Kriging surrogate model and mapping the nonlinear relationships between the IMPP and objectives, using the NSGA-II algorithm.Yin et al. [16] presented an investigation on birefringence and warpage of thermoplastic lenses using microinjection molding, and the simultaneous control of both properties was accomplished through multiobjective optimization using the desirability function approach.Wang et al. [17] proposed a two-stage optimization approach to improve heating efficiency and temperature distribution uniformity for a new water-assisted rapid heat cycle molding, using RSM and the NSGA-II algorithm.The multiobjective optimization method mentioned above can be called model-based optimization (MBO) methods, and the quality of models plays an important role in the optimization results.In addition, there is another important multiobjective optimization method called model-free optimization (MFO), which integrates the simplex method.The MFO has been successfully applied to control and optimize part weight, part dimensions, and the focal length of molded products.Kong et al. [18] and Zhu et al. [19] proposed a novel MFO method for batch processes with short cycle times and low operational costs.In order to eliminate defects such as shrinkage and flash, Yang et al. [20] proposed a systematic method that combined digital image processing and MFO to solve such problems, and the MFO used online measurement as feedback to determine the optimal settings.
In order to obtain high-quality plastic optical lens, this paper takes the surface waviness and optical performance into account and presents a set of technical procedures for optimization of IMPP, providing reference for manufacturing high-quality optical lens.The existing literature contained the studies on the optimization of IMPP in terms of birefringence and light transmittance of plastic optical lens.In this paper, another important optical quality indicator, haze ratio (HR), which has not been studied, is one of the optimization objectives.It refers to the percentage of the intensity of transmitted light deviating from the incident light by 2.5 ∘ in the overall transmitted light intensity.The greater the HR is, the more significantly lens' gloss and transparency will decrease.Meanwhile, another indicator reflecting the surface waviness, namely, peak-to-valley 20 (PV 20 ), was also taken as the optimization objective.The technical procedures adopted in this study are shown in Figure 1.First, the DOE was conducted, and then the key IMPP having significant effect on the two quality objectives were screened out via the Taguchi method combined with ANOVA, and the optimal combination of IMPP was initially determined.After that, a full factorial experiment was performed by the key IMPP.In order to obtain the nonlinear mapping relationships between the key IMPP and objectives, this paper used the data of the full factorial experiment as the sample, adopted multioutput support vector machine for regression (M-SVR) algorithm and the back-propagation neural network (BPNN) algorithm, and compared their prediction accuracies.Finally, the NSGA-II algorithm was used to optimize the two objectives, and a Pareto-optimal set was obtained and verified through experiment.

Optical Lens and Experimental Setup
Plastic optical lens needs to meet higher manufacturing requirements compared to that of common injection molding products, especially the surface waviness, and it must reach certain optical performance.Therefore, the equipment and process used for the manufacture of optical lens shall also meet more stringent requirements.Any large deviation in the product will seriously affect the quality of the lens.The lens studied in this paper, as shown in Figure 2, has a diameter of 20 mm, and its material is PMMA-HT50Y made by SUMITOMO.The precise injection machine adopts ARBURG ALL-ROUNDER 307 S 600-100.Its maximum injection weight is 230 grams, and the maximum injection pressure is 230 MPa.In this paper, the experiment on the lens involves HR and PV 20 , in which the haze ratio measurement device is the WGT-S transparent/haze tester, with the measurement accuracy of up to 0.01%; and the PV 20 measurement device is the Bruker-Contour GT-X3 three-dimensional surface profiler, which has the measurement accuracy of up to 0.1 nm and can accurately measure PV 20 , which is the average of the ten top points minus the ten lowest points on the wave surface in the measurement area.Evaluating the surface waviness with PV 20 can eliminate the errors caused by some anomalous points or defects in the test data to a certain extent.

Optimization Methodologies
The optimization of precision injection molding process needs to be carried out generally by DOE.Due to the high-standard requirement of product quality, it is difficult to meet the requirement through only conducting trialand-error experiment by experienced technicians for long time.This has the characteristics of lack of planning and time consuming.The optimization of optical lens mainly relies on DOE, which has overcome the shortcomings of the trial-and-error method.If there are multiobjective optimization requirements to the product, especially when the IMPP are conflicting with each other, it is hard to obtain the global optimal solution only by DOE.Therefore, the intelligent algorithm was introduced on the basis of DOE.A set of optimization procedures for manufacturing plastic optical lens has been presented, as shown in Figure 1.The study conducted the double-objective optimization on the HR and PV 20 to obtain the optimum IMPP for precise injection molding.The optimization process was divided into four stages.In the first stage, the ANOVA was carried out by using the statistical tool Minitab 15 with the Taguchi method as the DOE, in order to screen out the IMPP having a great influence on the two objectives, and then significant IMPP were obtained for full factorial experiments.In the second stage, the M-SVR and BPNN models were built and then trained and tested according to the data from full factorial experiments, to obtain the mapping relationships between key IMPP and objectives.In the third stage, a Pareto-optimal set was obtained through selection, crossover, and mutation, with the mapping relationships constructed by M-SVR as the fitness function of the NSGA-II algorithm.In the last stage, the obtained optimum IMPP were tested and verified to determine the effectiveness of the optimization method.

Taguchi Method.
The HR and PV 20 of plastic optical lens are affected by many IMPP.In order to study these indicators effectively and easily, the IMPP having little influence on the objectives need to be excluded using DOE.Taguchi method is an excellent DOE method.By arranging fewer experiments, we can obtain the influence degrees of all IMPP on the optimization objective, with combination of ANOVA.Taguchi method also introduces the concept of Signal to Noise ratio (S/N ratio), depending on which, we can ensure the product has stable quality and small fluctuation and is not sensitive to noise during production.The smaller the two optimization objectives are, the higher the quality of the lens will be, showing a smaller-the-better (STB) characteristic.The S/N ratio of the STB is defined as [21] where   is the response and  is the number of replications.
According to the calculated average value of the S/N ratio for IMPP at each level, we can determine their optimum level and finally obtain the optimum level combination.
It should be noted that only a few discrete levels were selected for the experiment, so the optimization scheme obtained by the Taguchi method is not necessarily the global optimum.However, the initial optimized results can be used to adjust the range of process parameters, which can provide reference for further intelligent algorithm optimization.In order to verify the validity of the Taguchi method, the optimal combination was obtained according to S/N ratio and then used to conduct the experiment.Compare the results after optimization with the previous ones.If the two objectives are smaller, it indicates that the Taguchi method is effective.Otherwise, it is necessary to reselect the IMPP to repeat the experiment.

BPNN Algorithm and M-SVR Algorithm.
The ANN is an information processing method developed from inspiration of the biological nervous system.The network model built through the study of samples and analysis of the data relations has the prediction ability.The BPNN algorithm was first proposed by Rumelhart et al. [22], and it is a representative neural network, involving back propagation of errors and forward propagation of information.With the BPNN algorithm, the nonlinear relationships between the objectives and IMPP can be established.In this paper, a three-layer network structure was built, as shown in Figure 3; the output of each neuron is equal to the sum of the weights of all neurons from the previous layer plus a bias, and then the activation function generated the output where   is the output of the th neuron of the hidden layer or the output layer,   is the network weight from the th neuron to the th neuron,   is the output of the th neuron, and   is the bias of the th neuron.(⋅) is an activation function of the hidden layer and the output layer, and a sigmoid function is selected as follows: where  is the gain factor of neurons and(⋅) is within the range of (−1, 1).The model gradually adjusts the network weight by minimizing the actual output and network output where () is the difference between the actual output and the network output and  is the learning rate.As the BPNN algorithm is based on the gradient descent method, there are some shortcomings such as large computation burden, long training time, demand for a large number of samples, and slow convergence.In this paper, the small number of experiments and training samples may result in that the prediction ability of BPNN is not so strong.Therefore, the M-SVR algorithm was also adopted to establish a nonlinear model in terms of objectives and IMPP, so that we could compare the prediction ability of the two algorithms.Pérez-Cruz et al. [23] improved the SVR algorithm for single objective and obtained the M-SVR algorithm which can handle multiple output objectives.They used the iterative reweight least square (IRWLS) method to solve the dual problem by redefining the -insensitive loss function in the hyperspherical space and finally obtained the prediction model of each output terminal.It is a type of effective method for Multi-Input Multi-Output (MIMO) modelling of small samples.
Suppose that the given sample data is  , = {(  ,   ) |  = 1, 2, . . ., },   ∈ R  , map the input data   into the highdimensional Hilbert space, and the mapped data (  ) ∈ R  may be obtained.The problem to be solved by SVR is how to select vectors W and b, so that the regression results obtained from ỹ = 0(  )W + b have the minimum error, where The -insensitive loss function can be extended into a multidimensional form.Supersede the original singledimensional error with  2 norm of the output variable error, and (  ) was defined as where   = ‖  ‖ = √     ,   =   − ỹ ,  is the allowable deviation, and  is the penalty factor.Through the first-order Taylor expansion, the approximate quadratic programming form of Formula ( 5) is obtained: where   and  are the constants independent of W, b, and then, we can minimize    (W, b) with the IRWLS method to calculate the optimal solution W s , b s .Finally, the relation between input and output is obtained: In order to evaluate the ability of BPNN and M-SVR models of fitting and predicting, the coefficient of determination ( 2 ) was used for analysis in this paper.The closer  2 is to 1, the more data the model can interpret, and the better the model performance will be. 2 is defined as follows: where ESS is the error sum of squares, TSS is the total sum of squares, and RSS is the regression sum of squares, and they satisfy the relation TSS = ESS + RSS.

NSGA-II Algorithm.
Since MOP often involves optimizing several objectives with trade-off relationships, it has no solution that can achieve the optimization of all objectives at the same time.The purpose of the multiobjective optimization algorithm is to obtain a representative Pareto-optimal set in the decision space, so that their distribution in the objective space Pareto front is broad and uniform as much as possible.There are different algorithms to solve MOP; the most widely used one is the NSGA-II algorithm developed by Deb et al. [24].The algorithm enables individuals in a Paretooptimal set to extend into the entire objective space uniformly through the introduction of the fast nondominated sorting algorithm and the elite retention strategy, which ensures the diversity of the population.The basic steps of the NSGA-II algorithm are as follows.
Step 1. Initialize parameters, such as population size, control variable limits, crossover and mutation probability, and termination criteria.
Step 2. Nondominated sorting: sort the population according to nondomination criteria.
Step 3. Calculate the crowding distance for each individual.
Step 4. Execute the genetic operator.Generate the intermediate population by selection, crossover, and mutation operators.
Step 5. Selection and combination: consolidate descendant and the contemporary populations; the descendant individuals are reset through rank and crowding distance.
Step 6.Return to Step 2 and repeat until termination criteria are reached.

Screening Significant IMPP.
There may be a number of IMPP that can affect the quality of the injection products.Taking all of them into consideration to carry out experiments means a very large workload, and it is not necessary.We can conduct the screening experiment on some IMPP that may have significant effect on the product quality, according to our production practice experience.In this paper, a total of eight IMPP were selected for experiments, including melt temperature, injection velocity, injection pressure, VP switchover, packing time, packing pressure, molding temperature, and cooling time.The Taguchi method was taken as the DOE method for the  18 (2 1 × 3 7 ) orthogonal experiment.The diameter of the optical lens is 20 mm, which belongs to a small size product.In order to improve production efficiency, six cavities of the same size have been designed.According to the suggestion of PMMA-HT50Y supplier, the number of cavities, and the size of the product, the levels of IMPP should be initially set out as shown in Table 1.
In order to reduce the experimental and measurement errors, each experimental treatment was carried out twice, totaling 36 experiments.At the same time, the experiments were arranged in random order to eliminate the interference of some hidden factors.The mold has six identical cavities.To ensure the tested data can reflect the whole batch of optical lens, this paper adopted a sequential sampling method.After each injection, one of optical lenses was selected to a box marked with serial numbers.Then, the next optical lens was taken in clockwise direction.The mean of the experimental results and the S/N ratio of the two objectives, HR and PV 20 , are shown in Table 2. ANOVA of HR and PV 20 is shown in Tables 3 and 4, respectively.The results show that the corresponding  value of some IMPP is not more than 0.05, indicating it has a significant effect on the objectives, which must be considered during the optimization.The IMPP which have a significant effect on HR include injection pressure, packing pressure, and molding temperature.The IMPP which have a significant effect on PV 20 include packing pressure, molding temperature, and cooling time.
HR and PV 20 have the smaller-the-better characteristic; the main effects plots for S/N ratios are shown in Figures 4  and 5, respectively, which further show the degrees effect of each process parameter on the objectives.At the same time, the initial optimization results of the Taguchi method were obtained, and the minimum process parameter combinations of HR and PV 20 are  = 220,  = 60,  = 80,  = 5.2,  = 4,  = 110,  = 60,  = 35 and  = 220,  = 70,  = 90,  = 5.6,  = 6,  = 100,  = 75,  = 25, respectively.In order to verify the initial optimization effect of the Taguchi method on the objectives, the experiments were conducted based on the initial optimization results, and the mean values of HR and PV 20 were 4.36% and 0.303 m, respectively.In comparison with the results in Table 2, the optimized objective is minimal, indicating that the Taguchi method has played a role in the initial optimization.

Building BPNN and M-SVR Model.
The optimization on objectives using the Taguchi method cannot guarantee  the global optimal value.In order to further improve the quality of optical lens, this paper established the mapping relationships between the key IMPP and objectives and conducted a global optimization using the multiobjective intelligent optimization algorithm, NSGA-II.The methods used for establishing the mapping relationships were M-SVR and BPNN algorithms, and their prediction abilities were compared.In this paper, the 3 4 full factors experiment was conducted with injection pressure, packing pressure, molding temperature, and cooling time as key IMPP.The experimental results were used as training and test samples for M-SVR and BPNN.The IMPP that have no significant effect on the quality objectives were initialized to the mean value of best levels according to their S/N ratios.For example, the best values of injection velocity for HR and PV 20 were 60 and 70, respectively, so its value was set to the 65.In addition, it was necessary to adjust the level value for one of the key IMPP.From Figure 4, it can be seen that the injection pressure only has a great effect on the HR and that the S/N ratios have a negative correlation; namely, the quality of optical lens decreases more quickly as a result of the increase in HR, with the increase in the injection pressure.Therefore, the level of injection pressure should be adjusted to smaller values.The IMPP and levels for 3 4 full factorial experiments are shown in Table 5.Similarly, the experiments were arranged in random order, and the experimental results are shown in Table 6.
The M-SVR and BPNN models were built with the data of 81 samples as in Table 6, of which 60 samples were randomly selected as training data, and those of the remaining 21 samples were used as test data.This paper used MATLAB R2014b as a modelling tool, the M-SVR algorithm was Libsvm toolbox developed by Chang and Lin [25], the kernel function was the radial basis function (RBF), and the performance goal was 1−4.The Levenberg-Marquardt algorithm [26,27] with a higher network convergence rate was selected in the BPNN algorithm, and the training variables were set as epochs 1000, were obviously smaller than those of the M-SVR model, indicating that the effect was not stable in the training of small samples.Especially in the testing stage, the prediction effort at some points was very large, completely unsuitable for precise prediction.Therefore, the M-SVR algorithm was selected to establish the relationships between the key IMPP and the two objectives.

Multiobjective Optimization.
In this paper, the optical performance and surface waviness of the optical lens needed to be optimized at the same time; this belongs to the typical multiobjective optimization.The two objectives have the smaller-the-better characteristic.Therefore, the issue of multiobjective optimization is defined as follows: Find: x Minimize:  (x) = ( 1 (x) ,  2 (x))    where  1 (x) and  2 (x) denote HR and PV 20 , respectively.x is the key IMPP vector (   )  .According to Table 5, it is known that b lower = (60 90 60 15)  and b upper = (80 110 90 35)  .There are many effective methods to solve multiobjective optimization problems, in which the traditional linear weighted method (LWM) and the popular NSGA-II algorithms are representative.In order to compare the optimization effects of these two methods, they have been tested separately in MATLAB R2014b.LWM is the method that can transform multiobjective optimization into single objective optimization.Each objective function was weighted and then summed up.Since the two objective functions have different units, they need to be dimensionless.The value of weights can reflect the importance of (x), and each set of weights,  1 and  2 , can obtain an optimized key IMPP by minimizing (x).In this paper, one objective reflects the optical performance, while the other reflects the surface waviness.Theoretically speaking, if the two objectives are equally important, the weights should be set as  1 =  2 = 0.5.In order to explore the overall optimization effect of LWM, 9 sets of weights have been set for testing.Using constrained nonlinear minimization Solver as the optimization tool, the operating parameters were set as follows: max iterations = 1500, function tolerance = 1 − 5, and the other stopping criteria were set as default value.The optimal key IMPP and two objectives have been achieved, as shown in Table 8. Figure 10 clearly shows that there is a dominating relationship between some optimization points.For example, under the condition of  = ( 1  2 ) = (0.9 0.1), the optimal solution is inferior to the optimal solution under the conditions of  = (0.7 0.3), (0.6 0.4).Therefore, if the weight of HR is greater, the optimization result can be unreliable.This is  mainly because the key IMPP and the two objectives have strong nonlinear relations, and (x) is a nonconvex function.
In fact, the optimal solution is highly likely to be the local optimum.
The NSGA-II algorithm has excellent multiobjective optimization ability.The functions  1 (x),  2 (x) established using the M-SVR algorithm were taken as the two fitness functions of the NSGA-II algorithm.The running parameters were set as follows: population size 100, maximum number of generations 200, stall generation limits 200, tolerance of fitness function 1−100, mutation probability 0.25, crossover probability 0.8, and Pareto fraction 0.1.The running results of the NSGA-II algorithm are shown in Figure 10.The Pareto-optimal set has 100 × 0.1 = 10 individuals.After the approximate treatment was performed on the optimized IMPP and the corresponding objectives, the results are as shown in Table 9.The range of HR is 3.826 to 5.766 and that of PV 20 is 0.347 to 0.791.As can be seen from Table 6, the known minimum values of HR and PV 20 are 3.82 and 0.334, respectively.The Pareto-optimal set does not contain the case in which the two objectives are minimum, indicating there is a conflict between the two objectives and the optimal status cannot be achieved at the same time.Figure 10 further shows that PV 20 gradually increases with the decrease of HR.In contrast to LWM and NSGA-II, all the Pareto-optimal  combination of Experiment number 5 was selected, and the injection plastic optical lens is shown in Figure 11.

Conclusion
In order to consider the optical performance and surface waviness of the plastic optical lens at the same time, this paper presents a set of optimization procedures with four stages, taking HR and PV 20 as the quality objectives.The NSGA-II algorithm was applied to solve this multioptimization problem.In particular, the introduction of the M-SVR algorithm as the fitness function played an important role, through which we improved the prediction accuracy of the quality objectives.The Pareto-optimal set obtained in this study improved the quality of plastic optical lens comprehensively.The conclusions are as follows: (1) The IMPP having a significant effect on both objectives, HR and PV 20 , include packing pressure and molding temperature.The IMPP only having a significant effect on HR is injection pressure.With the increase in injection pressure, the optical performance of the lens decreases continuously.The IMPP only having a significant effect on PV 20 is cooling time.With the increase in cooling time, the lens surface waviness was increased substantially at the beginning and decreased slightly subsequently.
(2) When establishing the relationships between the objectives and the key IMPP, the 3 4 full factorial experimental results of the IMPP were used as samples.At the same time, the BPNN algorithm and the M-SVR algorithm were used to establish the relationships.The results show that the determinant coefficients  2 of the M-SVR model, which are close to 1, are bigger than that of the BPNN model in the training and testing stages, indicating that the M-SVR algorithm has high prediction accuracy when dealing with low-dimensional small samples.
(3) Using the NSGA-II algorithm can effectively find the Pareto-optimal set of IMPP, so that the two objectives HR and PV 20 are as small as possible, to achieve the optimal balance.Through the verification on the optimization results, the MPE of HR and PV 20 reached a high accuracy of 7.16% and 9.78%, respectively.Finally, the optimization results may be selected from the Pareto-optimal set according to the bias of the two quality objectives.
row vector, b  ( = 1, . . ., ) is a constant, and 0(  ) is a nonlinear mapping function.Thus, the minimized objective function   is established to solve W and b: min W,b

Figure 6 :
Figure 6: Comparison of HR prediction results in training set.

Figure 7 :
Figure 7: Comparison of PV 20 prediction results in training set.

Figure 9 :
Figure 9: Comparison of PV 20 prediction results in testing set.

Figure 11 :
Figure 11: The injection plastic optical lens.

Table 1 :
The IMPP and their levels for the Taguchi experiment.

Table 2 :
Experimental results and S/N ratio of HR and PV 20 .

Table 7 :
The  2 of BPNN model and M-SVR model at training and testing stage.

Table 10 :
Confirmation experimental results for Pareto-optimal set of HR and PV 20 .