Using improved simulated anneal arithmetic to

Special Issue Article

Using improved simulated anneal arithmetic to improve semiconductor packaging ball placement process capability

Advances in Mechanical Engineering 2018, Vol. 10(2) 1–12 Ó The Author(s) 2018 DOI: 10.1177/1687814017746514 journals.sagepub.com/home/ade

Shen-Tsu Wang1 and Meng-Hua Li2

Abstract This study used improved simulated annealing arithmetic to improve the semiconductor packaging ball placement process capability. The simulated annealing computing process was combined with an artificial neural network, Metropolis algorithm, and sequential Gaussian simulation. Grey relation analysis was used as the target value, which is intended to obtain the optimum parameter design of process capability. The results showed that the substrate design-Pad Open is 0.30 mm ball size, 0.275 mm Pad open, Profile type B, maximum temperature 236 (°C), solder melting time 53 (s), preheating temperature 176 (°C), and a solder ball component, as SAC405 can result in a parameter design for optimal solder ball shear strength of Cpk = 2.38 and G0i = 0:87. This optimum parameter design process is better than other methods, and it reduces the cost waste of extra defectives. Keywords Semiconductor packaging, ball placement, improved simulated annealing, grey relation analysis, optimum parameter design process

Date received: 12 July 2017; accepted: 10 November 2017 Handling Editor: Peter Nielsen

Introduction The purpose of packaging is to provide the finished product (package) with an interface, in order to connect the internal electrical signal via the packaging material to the system and to prevent the silicon chip from being damaged and corroded by external force, water, moisture, and chemicals. The main process of semiconductor packaging is as shown in Figure 1. This research scope is Ball Placement, and the purpose of Ball Placement is to place the Solder Ball on the substrate. The Solder Ball is the circuit connection and heat sinking conductor between substrate and mainboard. The ball placement process consists of a Ball Mount, IR-Reflow, and Cleaner. In terms of ball placement, the pad of the substrate is coated with flux, and the Solder Ball is placed on the Pad of the substrate at the high temperature

(150°C–245°C) of IR-Reflow furnace, the Solder Ball reacts with flux and half of it melts on the substrate pad. The Cleaner removes the excess flux and the oxide after reaction from the substrate to avoid contaminating the substrate. In terms of the research setting range of ball placement, this study uses the reflow parameter as the research subject and uses solder ball strength 1

Department of Commerce Automation and Management, National Pingtung University, Pingtung, Taiwan, R.O.C. 2 Department of Industrial Management, National Formosa University, Huwei, Taiwan, R.O.C. Corresponding author: Shen-Tsu Wang, Department of Commerce Automation and Management, National Pingtung University, 51 Min Sheng E. Road, Pingtung 900, Taiwan, R.O.C. Email: [email protected]

Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License (http://www.creativecommons.org/licenses/by/4.0/) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).

2

Advances in Mechanical Engineering

Figure 1. Assembly process and research process.

process stability as the target response value. The thin and fine ball grid array (TFBGA) (Figure 2) product is validated, the solder ball shear strength (Figure 3), the fatigue life of the TFBGA can be obtained (larger-thebetter, LTB), warpage (smaller-the-better), residual stresses (smaller-the-better) are the quality characteristics, and a higher solder ball shear is better, thus, it is a LTB characteristic. However, higher strength may be accompanied by IMC (intermetallic compound) layer embrittlement, which may result in another reliability problem; thus, it is not discussed in this study for the moment.1,5–8 The purpose of ball placement is to lay the solder ball on the substrate. The solder ball is the conductor of electrical connection and heat dissipation between the substrate and the mainboard. The process of ball placement can be divided into three procedures: Ball

Mount, IR-Reflow, and Cleaner. As for Ball Mount, it coats a layer of flux on the pad base of the substrate first, and then it mounts the solder ball on the substrate. Through the high temperature of the reflow oven (about 150°C–245°C), half of the solder ball will be fused to the pad of the substrate after it reacts with the flux. Cleaner is for removing the excess flux and the oxide after the reaction on the substrate to avoid contamination of the substrate. This study used the reflow parameter as the study subject and used the processing stability of the solder ball strength as the target response value. Taking the TFBGA product as the subject for verification, the thrust strength of the solder ball (Figure 3) was used as the quality characteristics, in which case the higher the thrust is of the solder ball, the better it is, and so it exhibits the LTB characteristic.

Wang and Li

3 The optimum architecture of the relation of the standardized data was realized using the back-propagation neural network method. The approximately global optimum was achieved by employing a genetic algorithm and the improved simulated anneal arithmetic. Finally, the optimal solutions obtained by the aforesaid methods were validated by a confirmation experiment to find out the optimum parameter combination.9–24 The higher the thrust of the solder ball, the better it is, and therefore, it exhibits an LTB characteristic25 SNLTB =  10log10 ðMSDÞ " # n 1X 1 =  10log10 n i = 1 y2i

Figure 2. Solder bump structure. References: Lau and Lee,1 Lee et al.,2 Tee et al.,3 and Akbari et al.4

SNLTB is SN (signal-to-noise) ratio of the LTB characteristic, MSD is mean square error of deviation from the target value, and yi is observation data.

IMC IMC is a chemical composition, and so the energy must be given to the formation of IMC, which is the reason why the solder paste needs to be heated during the welding process. In addition, in the composition of the solder paste, only the pure tin (Sn) has a diffusion reaction with the copper base (e.g. organic solderability preservative, I-Ag, I-Sn surface-treated board) or nickel base (Electroless nickel immersion gold surface-treated board) in the strong heat, and thus a solid interface IMC is generated.1–4

ANOVA

Figure 3. Solder ball shear test. References: Lau and Lee,1 Lee et al.,2 Tee et al.,3 and Akbari et al.4

However, when the strength rises, it may be associated with the embrittlement of the IMC layer, and it may cause another reliability problem. Thus, it is not included in the research scope of this study.1–4 This study used the Taguchi method and multiple quality characteristics to calculate the signal-to-noise (SN) ratio and factor level combination from the original data measured after the soldering experiment and took multiple quality characteristics to work out analysis of variance (ANOVA) and contribution rate data, which are integrated into a cross-over analysis table in order to obtain four visualized optimal solutions. The SN ratio was used as an input value, the priority value of various schemes was found by the conversion mode of grey relation analysis, and a rough optimal solution was obtained.

It is a common statistical model in data analysis, and it mainly explores the relationship between the dependent variable of the continuous data type and the independent variable of the categorical data type.25 Literature review of this study is shown in Table 1.

Research method Current condition analysis and parameter setting The current condition data were the calculated process capability index Cpk = 1.34 ( 1:33 acceptable) as shown in Table 2, meaning the process capability is unstable and insufficient, and there is a large room for improvement. The experiment is designed according to the ball placement process parameters of the case company, and the six most important controlled factors in the experiment are selected from the process parameters according to the field engineer’s advice, which are substrate design-Pad Open, Profile type, maximum temperature (°C), solder melting time (s), preheating temperature (°C), and solder ball composition, and the

4

Advances in Mechanical Engineering

Table 1. Literature review of this study. Title

Main contents

Source

Process and decisionmaking method

A design optimization method based on reliability for an improved multitarget system is developed to explore the design concept of the conceptual car body under uncertainty. Dimension optimization is combined with a reliability-based design optimization method to improve the performance of the accelerometer. The Six Sigma Method is used in parametric design optimization applications in the aluminum extrusion plant. A parameter design from the integrated production process for the model of the forecasted product is proposed, which takes into account the design, assembly, and process. A genetic algorithm based on the genetic algorithm with two different local search strategies is proposed to maximize the module density, and it uses a more general objective function, where tunable parameters can solve the resolution limit. A heuristic method based on simulated annealing (SA) for solving the open location-routing problem is proposed. The goal is to minimize the total cost, including equipment operating costs, vehicle fixed costs, and travel costs. A hybrid real coded genetic algorithm with particle swarm optimization and mixed artificial immune algorithm for solving constrained global optimization problems is proposed. The ‘‘Analytic Hierarchy Process Method’’ and the grey relational analysis based on degree of similarity are used for empirical investigation. An enhanced resonance search algorithm is developed, which makes it possible to quickly escape the local optimal solution. A powerful feature extraction method is designed for classifying multiple types of signals, identifying valuable features from the original display graph data, and finding the effective classifier of the feature. A flexible process scheduling combinatorial optimization study of nonhomogeneous parallel machines minimizes the maximum completion time or manufacturing time. The Grey Decision Theory is used for management application basis. This book describes the decision-making methods that apply to different situations. It observes a large number of pictures, deals with data quality problems, uses customer background and knowledge to develop appropriate statistical models, and analyzes different management implications. A multi-objective (minimization of availability and maintenance cost) preventive maintenance scheduling model of the continuous operating system is proposed. A multi-objective genetic algorithm is used to optimize the objective function. This reference presents shear stress and shear creep strain hysteresis loops at the fillet welds, shear stress history, shear creep strain history, and creep strain density range, so as to better understand the thermodynamic behavior of the wafer level chip scale packages of the lead-free solder bumps. The fatigue life of the fine-pitch BGA package under thermal cycling is mainly analyzed. A detailed drop test and a simulation of the thin-profile fine-pitch BGA and the very-thin-profile fine-pitch BGA package are performed at the board level with the use of an internally developed test program. A double cantilever beam is configured with the test BGA and the printed circuit boards assembled with the microelectronics package. As for the scanning probe lithography, it proposes that the anodic oxidation of the atomic force microscope should be used to define the line thickness directly on the surface of the silicon chip. It improves the optical diffraction of conventional optical lithography. In consideration of the uncertain system parameters, the optimal design of the base isolation system used to control seismic vibration is usually the unconditional expected value through minimizing the mean square response of the structure, regardless of the effect due to system parameters’ uncertainty.

Duan et al.5

Discussion on process parameters

Teves et al.6 Ketan and Nassir7 Boorla and Howard8 Mu et al.9

Yu and Lin11

Wu12 Yang et al.13 Maheri and Narimani17 Siuly and Li18 Dai et al.23 Deng26 Chien27 Hamada and Hamada28 Adhikary et al.29

Lau and Lee1

Lee et al.2 Tee et al.3 Akbari et al.4 Kuo et al.10

Roy and Chakraborty14

(continued)

Wang and Li

5

Table 1. Continued Title

Main contents

Source

Due to difficulties in the conceptual design of the compact heat exchanger, it adopts a multi-level, multi-factor, and non-structural fuzzy optimal decision model to select the compact heat exchanger. An application of the harmony search–based algorithm in the optimal detail design of the special seismic moment–reinforced concrete structure under seismic load based on the American standard is introduced. This reference used the global optimization framework based on real coding genetic algorithm (GA) as the common template for designing the optimum proportional-derivative controller to control four different multiwing chaotic systems. The emulational auxiliary power model and the scalable component model are used to improve the accuracy of the design process by means of tradeoff analysis and evaluation methods for multi-objective design problems. In order to improve the solidification speed of soft processing technology, the silicone rubber composite mold material is designed based on multiobjective optimization (MOO), which is contradictory to the target. The reference performed statistical analysis of the realization of microelectronic electro-discharge machining process parameter measurement and acquisition system and noted its effect on the process performance. The linear and non-linear regression methods are used to obtain the most important predictive equation for the micro-electroplating electro-discharge machining process. The choice of processing parameters is important for obtaining the required accuracy for good surface quality in cylindrical grinding. This literature focuses on the optimization of AISI 4140 steel’s continuous and interrupted cylindrical grinding. The emulational annealing and the genetic algorithmic computing technology are integrated to find a set of optimal cutting condition values. The main objective is to minimize the annual cost of the hybrid system, the environmental costs, and the loss of power supply probability due to the use of hybrid power generation systems to avoid pollutant emissions.

Zhou et al.15 Akin and Saka16 Das et al.20

Sarioglu et al.21 Nandi et al.22 D’Urso et al.24

Ko¨klu¨30

Zain et al.31 Shayeghi and Hashemi32

BGA: ball grid array.

Table 2. Current data (solder ball strength). 634.2 646.3 634.3 .. . 559.8

496.7 492.0 703.5 .. . 665.4

431.1 464.5 426.6 .. . 539.0

levels of the control factors are set as shown in Table 3, and the improved simulated annealing (SA) is used for experimental design.11,13 As for the temperature characteristics of the alloy only, the difference is not big, but the slurry interval of SAC105 is larger. This means that when SAC305/405 are completely in the liquid phase, SAC105 is still partially not fully liquified. Therefore, if SAC105 is used in the solder paste, then when the solder ball uses SAC305/405, it will have a higher probability of generating a void. Kim et al.19 used a sandwich joint test piece to carry out the relevant tensile test on SAC305. The experimental results show that the tensile strength

559.8 478.9 582.5 .. . 506.6

623.3 524.1 539.7 .. . 588.2

of the test piece is larger than that of the block. It shows that the thickness of the medium metal phase has a considerable correlation with the anti-shear strength of the joint surface, and the thickened medium metal phase will make the breaking strength of the joint surface decrease significantly.1–4 The solder ball shear strength is a quality characteristic, in which a higher solder ball shear is better, and so it presents a LTB characteristic. However, the high strength may be accompanied by IMC layer embrittlement. The viscoplastic strain energy density is substituted in the model proposed by Darveaux, thus obtaining the fatigue life of TFBGA (LTB).2 The

6

Advances in Mechanical Engineering

Table 3. Control factors and levels. Standard

Level 1

Level 2

Level 3

A. Substrate design—Pad open

0.30 mm ball size, 0.275 mm Pad open

0.20 mm ball size, 0.275 mm Pad open

–

230 40 140 SAC105

235 60 160 SAC305

245 120 180 SAC405

B. Profile type

C. Maximum temperature (°C) D. Solder melting time (s) E. Preheating temperature (°C) F. Solder ball composition

References: Lau and Lee,1 Lee et al.,2 Tee et al.,3 and Akbari et al.4

Table 4. Cross-over analysis of quality characteristics.

Solder ball shear strength SN ratio significance Contribution rate after weighting (%) The fatigue life SN ratio significance Contribution rate after weighting (%) Warpage SN ratio significance Contribution rate after weighting (%) Residual stresses SN ratio significance Contribution rate after weighting (%) Optimum parameter level combination

A

B

C

D

E

F

A1 (F = 22.36)* 23.82 A1 * 12.53 A1 * 6.29 A2

B2 (F = 9.87)* 10.29 B3

C2 (F = 10.21)* 9.31 C2

D3

F1

0.93 D1

E2 (F = 5.19)* 4.09 E3

0.91 B2 * 6.26 B3

6.12 C1 * 4.32 C3

1.08 D1

3.91 E1

1.32 F2 (F = 3.29)* 1.68 F1

6.28 D2

0.87 E1

1.06 F2

2.62 A1

5.19 B2

1.23 C1 C2

3.26 D1 D2

1.99 E2

0.87 F2

*F-value is relative efficiency of factors effect in the Taguchi method. F-value \ 1: the factor effect is small; F-value . 2: the factor effect is medium; F-value . 4: the factor effect is large.25,30

literature has paid increasing attention to the warpage (smaller-the-better) and residual stresses (smaller-thebetter). Previous studies indicated that the first cause for the warpage of various components of an IC package is the non-uniform warpage resulting from different material coefficient of thermal expansion, or the curing shrinkage of epoxy molding compound, or the nonuniform temperature change.3,4 According to the specific quality characteristic SN ratio’s significance of the Taguchi method and the specific quality characteristic contribution rate obtained by ANOVA, the cross-over analysis comparison table can be used for judgment, so as to obtain the optimum parameter combination of this part, as shown in Table 4. In terms of significance, according to the contribution rate after weighting, among the ABCDEF factors, the levels reaching the optimum contribution degree are A1B2C2D1E2F2. However, the maximum weighted ratio of solder ball shear strength is 40% in this experiment, and the Factor C solder ball shear

strength test has a significant impact. Therefore, aside from considering C2, C1 can also be considered. The quality parameters of Factor D are insignificant, but the contribution rate in warpage is the highest, and the residual stresses rank second. Thus, D1 and D2 factors shall be considered as well. The four optimum combinations are hence A1B2C1D1E2F2, A1B2C1D2E2F2, A1B2C2D1E2F2, and A1B2C2D2E2F2.

Grey relation analysis optimum parameter level combination The optimum parameter level combination considering the G0i of grey relation analysis of cost and quality aspects.11 The degree of relationship among sub-systems or elements could be evaluated through grey relational analysis,26–28 and important influential factors to the development trend are then determined to learn the major features of the system through the following steps:

Wang and Li

7

Table 5. Grey relation analysis. Divisor ½R0i ðkÞ

A

B

C

D

E

F

Level 1 Level 2 Level 3 Best level

0.89 0.63 – A1

0.88 0.66 0.59 B1

0.76 0.61 0.58 C1

0.52 0.69 0.59 D2

0.66 0.72 0.51 E2

0.83 0.76 0.53 F1

Step 1. Normalize original data. Normalize by dividing the original data zi (k) with the mean value of the sequence shown in equation (1)

zi ðk Þ ; i = a, :::, d; k = A, :::, N R i ðk Þ = N P zi ðk Þ k =1 N

ð1Þ

Step 2. Designate the standard sequence and calculate the difference sequence. Take the mean value as a standard sequence, that is, sequence 0; the difference sequence D0i (k) indicates the absolute difference of elements k between the other sequence i and the standard sequence 0, as expressed in equation (2) D0i ðk Þ = jR0 ðk Þ  Ri ðk Þj; i = 1, 2, 3:::; k = A, :::, N

Dmax = Max D0i (k)

ð3Þ

Dmin = Max D0i (k)

ð4Þ

i, k

Step 4. Calculate grey relational coefficient: R0i (k). The relational coefficient: R0i (k) is defined below, of which w is the adjustment factor, as shown in equation (5) R0i ðk Þ =

Dmin + wDmax D0i ðk Þ + wDmax

ð5Þ

Step 5. Calculate the grey relationship G0i between each sequence and the standard sequence. The grey relationship G0i is defined as in equation (6) G0i =

N X R0i ðk Þ N k =A

Network model

Training

Testing

6-1-1 6-2-1 6-3-1 6-4-1 6-5-1 6-6-1 6-7-1 6-8-1 6-9-1

0.0856 0.0923 0.0026 0.000938 0.000834 8.76E–09 9.28E–09 1.263E–08 2.234E–08

0.0613 0.0823 0.0734 0.0265 0.0346 0.0123 0.0862 0.0923 0.0723

MSE: mean squared error.

Table 5 lists the optimum parameter level combination considering G0i of grey relation analysis for the cost and quality aspects.

ð2Þ

Step 3. Calculate maximal difference Dmax and minimal difference Dmin as expressed in equations (3) and (4) i, k

Table 6. Artificial neural network MSE.

Artificial neural network setting and execution The mode is selected as TanhAxon, the lower bound of the learning rate is set as 0.03, the upper bound is 0.2; the lower bound of the initial value of momentum is set as 0.6, the upper limit is 0.9; the number of training cycles is set as 10,000; the number of layers of the artificial neural network is set as 1–9; the performance index is mean squared error (MSE) for determining network quality, as expressed as equation (7) n  P

MSE =

i=1

0

Si  Si

2

n

ð7Þ

where n is the number of samples, Si is the actual value, 0 and Si is the prediction value. The standard criterion is that the MSE value of training and testing is lower the better in the test values; thus, the 6-6-1 artificial neural network is used as the executed network architecture according to Table 6.

ð6Þ

Step 6. Conduct sequencing according to the grey relationship.

Optimum parameter analysis of genetic algorithm According to the trained neural network model, the model has six input layers, one output layer, six nodes of a hidden layer, the learning rate is 0.15, and the

8

Advances in Mechanical Engineering

Table 7. Genetic algorithm implementation results. Network model 6-6-1

G0i

Generation number

Maximum

Minimum

Mean

0.71 0.78 0.82 0.86 0.83 0.85

0.49 0.39 0.43 0.46 0.56 0.59

0.53 0.59 0.69 0.68 0.65 0.63

500 1000 1500 2000 3000 4000

Table 8. Summary analysis of SA computed performance. Parameter

Numerical value

Initial temperature (T0 ), S0 Final temperature (Ti ), Si Cooling factor a Total number of cooling Optimal strategy temperature layer Optimal strategy target value

160, 0.42 1, 0.69 0.96 38 26 0.86

SA: simulated annealing.

Table 9. Improved SA implementation results. Network model

G0i Maximum

Minimum

Mean

6-6-1

0.62 0.66 0.71 0.87 0.68 0.59

0.39 0.38 0.41 0.42 0.32 0.31

0.46 0.52 0.56 0.69 0.55 0.46

SA: simulated annealing.

Figure 4. Improved SA parameter convergence map.

Markov chain 300 400 500 600 700 800

momentum is 0.95. The G0i value is taken as the target, and the input layer is the input parameter value for gene screening. According to the genetic algorithm principle, the data have approached to convergence in Generation 2000, the minimum value of G0i of 3000 generations and 4000 generations increases slightly, but it is unnecessary to waste this cost according to practice. This study lays emphasis on G0i value LTB; thus, the executed standard is defined as 2000 generations, and the results are shown in Table 7.15–17,28,29

Optimum parameter analysis of improved SA The neural network training model is called, the G0i value is set as the target value, the initial temperature is 160, the stop temperature is 1, and cooling factor a is 0.96. It is obvious that the data have approached to convergence when the Markov Chain is 600; thus, the test stops, and the results are as shown in Tables 8 and 9 and Figure 4. The SA structure of this study consists of the Metropolis algorithm, sequential Gaussian simulation, and an annealing process, as described below.11–16,31 The Metropolis algorithm is as follows: Step 1. When iteration t = 0, an arbitrary feasible solution S(0) = S is generated. The sequential Gaussian simulation generates S(t). Step 2. A neighborhood solution is generated according to the state S of current solution S(t) using an effective disturbance mechanism 0

S = N ½S ðtÞ  S, N ½S ðtÞ 6¼ F (it can contain S or not contain S)

Wang and Li

9 0

0

Step 3. Calculate increment DE = C(S )  C(S). 0 C(S ) as the current solution and C(S) is the value of the iteration. 0 Step 4. If DE \0, the solution of the new state is 0 better, S is accepted as the new current solution.

function of the actual value and the distribution of the simulation values in the space. The simulation procedure is as follows: 1.

0

If DE  0, the probability EXP( DE =kT ) is used to 0 judge whether or not to accept S as new current solution. A random number u subject to U (0, 1) is generated, 0 0 if EXP( DE =kT).u, and S is temporarily accepted as new current solution; otherwise it is not, where k is the Boltzmann constant and T is the current temperature. Step 5. t = t + 1 and judge whether the algorithm stops according to the preset convergence criteria (stop search rule), if yes, the computational process is stopped; otherwise proceed to Step 2. The sequential Gaussian simulation is as follows. The sequential Gaussian simulation is based on the stochastic simulation theory, where the known data points, and all the simulation values in simulated positions, are combined with the conditional probability to simulate the new location point. This simulated new value will be put in the condition data of the next simulation, and the remaining can be deduced accordingly till all unknown point values in the simulation region are simulated. A random variable Z(x) of arbitrary statistical distribution is converted by normal score into a multivariate Gaussian function. The execution steps are described as follows:17,18,20,21,32 1. 2.

3. 4. 5. 6. 7.

Define a random path. Decide whether to use simple kriging or ordinary kriging and normal variation functions. While the simple kriging is generally used, if the data are large and the local average varies significantly with the region, then ordinary kriging is used. The variation function is selected by cross-validation. The selected kriging method and variation function are used to determine the complementary cumulative distribution function (CCDF) parameter of Y (x) on x coordinate. The CCDF simulates the Y (x) value of the mth point in the simulation region. The simulated value obtained in Step 3 is added to the data. Execute the next point simulation till all the simulated points are simulated. Convert the simulated normal y(x) value into z(x) value. Execute simulation once and discuss whether the variation functions of various simulation values conform to the theoretical variation

2.

3. 4.

5.

An initial numerical model is constructed, this model simulates the initial melting of an actual annealing phenomenon, and a random value is allocated to each grid. If the problem to be handled is of co-simulation, that is, the secondary variable can be used, then the initial values can be allocated according to the appropriate condition distribution, as obtained from validation. The objective function is defined as the measurement of the gap between the wanted spatial characteristics and the implementation result characteristics, for example, the gap between the variation function of the implementation result and the theoretical variation function. The image is corrected by giving a random location a new value. The objective function value decreases, and the correction is always accepted. If the objective function value increases, the correction is accepted with a specific probability. When the probability is accepted with an inappropriate exchange, and the value decreases, the disturbance process continues until a low objective function state is far.

The annealing process is as follows: Step 1. Select an initial solution S0 as the initial current solution, and let S(0) = S0 , set the initial temperature T0 and iteration i = 0. Step 2. Let T = Ti , and use T and Si to call the Metropolis algorithm, in order to judge whether current solution S is changed to Si . Step 3. Reduce temperature T according to a specific model, that is, Ti + 1\T , i = i + 1. Step 4. Check whether the annealing process is finished. If yes, the current solution Si is used as the final solution; otherwise proceed to Step 2.

Validation experiment According to Table 8, the grey relation analysis, artificial neural network, and the improved SA make up the optimum combination of the validation experiment. Its G0i is 0.87, better than other combinations, and the G0i value is 0.87 in the confidence interval, proving that the confirmation experiment results match the experimental data. A total of 10 confirmatory experiments were conducted in this study. The 95% confidence interval was

SN: signal to noise.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

Experiment number

1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2

A

1 1 1 2 2 2 1 1 1 2 2 2 1 1 1 2 2 2

B

Impact factor

1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3

C

1 2 3 1 2 3 2 3 1 3 1 2 2 3 1 3 1 2

D 1 2 3 2 3 1 1 2 3 3 1 2 3 1 2 2 3 1

E 1 2 3 2 3 1 3 1 2 2 3 1 1 2 3 3 1 2

F 496.4 600.2 586.6 536.9 650.8 585.5 551.3 580.2 490.6 495.4 608.4 570.8 533.2 512.3 468.9 532.6 549.6 501.6

608.2 555.6 459.5 469.3 660.2 458.4 512.3 511.3 530.8 536.6 480.1 560.3 599.6 464 529.2 518.4 483.2 532.6

526.1 533.2 622.6 561.4 660.2 621.5 505.4 610.1 620.3 528 580.1 537.8 580.3 559.8 510.6 646.3 588.2 519.3

546.9 613.9 512.1 610.1 653.2 511.2 611.6 500.6 540.2 542.8 501.2 616.8 620.2 534.5 500.8 529.1 493.5 503.9

599.8 570.2 499.1 540.8 651.9 460.2 499.9 589.2 508.1 510.7 530.9 577.6 530.1 513.9 466.1 537.8 533.7 499.9

Experimental response value (solder ball strength)

Table 10. Experimental response value and SN ratio of the L18 orthogonal array.

512.9 560.3 560.3 498.6 631.6 561.1 619.8 509.4 555.8 499.9 496.3 584.3 518.3 508.6 480.2 616.9 612 511.6

612.8 531.2 524.6 557.3 660.6 522.4 510.9 614.8 512.9 501.9 522.8 538.4 538.1 556.8 510.8 639.5 540.6 616.9

536.2 610.2 614.6 617.2 655.8 618.4 493.8 498.2 530.1 538 501.3 499.8 488.2 609.2 600.9 530.8 500.1 502.8

510.9 570.6 498.9 497.3 650.9 512.3 549.2 512.3 508.1 610.7 498.6 573.5 506.2 514.9 478.3 577.8 496.3 499.8

520.8 555.6 583.2 575.6 652.4 500.9 512.6 598.3 520.1 526.3 581.8 560.4 577.2 515.6 482.1 526.3 580.9 501.3

54.69 55.09 54.63 54.66 56.29 54.43 54.52 54.75 54.46 54.43 54.41 54.96 54.72 54.41 53.96 54.96 54.53 54.26

SN

0.775 0.832 0.713 0.768 0.873 0.732 0.742 0.779 0.722 0.721 0.732 0.799 0.746 0.582 0.591 0.812 0.749 0.659

G0i

10 Advances in Mechanical Engineering

Wang and Li

11

Table 11. Validation experiment judgment value. Optimum experimental combination

A

B

C

D

E

F

G0i

Cpk

Taguchi Cross-over analysis table for quality characteristics (1) Cross-over analysis table for quality characteristics (2) Cross-over analysis table for quality characteristics (3) Cross-over analysis table for quality characteristics (4) Grey relation analysis Integration of neural network and genetic algorithm Integration of neural network and SA Integration of neural network and improved SA

A1 A1 A1 A1 A1 A1 A1 A2 A1

B2 B2 B2 B2 B2 B1 B2 B2 B2

C2 C1 C1 C2 C2 C1 238(°C) 241 (°C) 236 (°C)

D2 D1 D2 D1 D2 D2 56 (s) 59 (s) 53 (s)

E2 E2 E2 E2 E2 E2 178 (°C) 179 (°C) 176 (°C)

F3 F2 F2 F2 F2 F1 F3 F3 F3

0.41 0.68 0.62 0.69 0.61 0.51 0.86 0.69 0.87

1.39 2 1.98 2.03 1.97 1.48 2.32 2.03 2.38

SA: simulated annealing.

verified before the G0i confirmatory experiment. The 95% confidence interval value for this experiment is [0.568, 0.893]. Through the L18 orthogonal array, a higher SN ratio represents better quality (smaller loss), and so the highest SN ratio is the best combination of the parameter level. The variation of the product produced at this parameter level is also at a minimal. The factor effect graph shows the best significant level. From the SN ratio data of the experiment, we can find that the best level of factor combination should be A1B2C2D2E3F3. The confirmatory experiment can verify whether or not the average of the predicted optimal conditions is valid, as shown in Table 10. Therefore, the confidence interval of the expected SN ratio of the confirmatory experiment is [56.28 6 3.23] = [53.05, 59.51]. Through the above confirmatory experiments, the SN ratio falls within the 95% confidence interval, indicating that the experimental results are successful.

Results and discussion According to Table 11, this study is the optimum combination of a validation experiment, the G0i value is 0.87, Cpk = 2:38, which is better than other combinations, and the confidence interval proves that the confirmation experiment results match the experimental data.

Conclusion and future research According to the case analysis in this study, substrate design-Pad Open as 0.30 mm ball diameter, 0.275 mm Pad open, Profile type B, maximum temperature 236 (°C), solder melting time 53 (s), preheating temperature 176 (°C), and solder ball component SAC405 can result in the parameter design of the optimal solder ball shear strength of Cpk = 2.38 and G0i = 0:87, thus reducing the waste costs of extra defectives and improving the process. In the future, as for multiple processes before and after Ball Placement, with an overall consideration of more of its parameters, the decision-making method

in this study can be used to determine the multiple parameters of its multiple processes. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding The author(s) received no financial support for the research, authorship, and/or publication of this article.

References 1. Lau JH and Lee SWR. Reliability of wafer level chip scale package (WLCSP) with 96.5Sn-3.5Ag lead-free solder joints on build-up microvia printed circuit board. In: Proceedings of the international symposium on electronic materials and packaging, 2000 (EMAP 2000), Hong Kong, China, 30 November–2 December 2000, vol. 23, pp.19–27. New York: IEEE. 2. Lee KL, Liu HH and Pan WF. Fatigue life estimation of thin fin-pitch ball grid array subjected to thermal cyclic loading. J Chung Cheng Inst Technol 2013; 42: 223–232. 3. Tee TY, Ng HS, Lim CT, et al. Impact life prediction modeling of TFBGA packages under board level drop test. Microelectron Reliab 2004; 44: 1131–1142. 4. Akbari S, Nourani A and Spelt JK. Effect of adhesive fillet geometry on bond strength between microelectronic components and composite circuit boards. Compos Part A: Appl S 2016; 87: 228–236. 5. Duan LB, Li GY, Cheng AG, et al. Multi-objective system reliability-based optimization method for design of a fully parametric concept car body. Eng Optimiz 2017; 49: 1247–1263. 6. Teves AC, De Lima CR, Passaro A, et al. Structural design of high-performance capacitive accelerometers using parametric optimization with uncertainties. Eng Optimiz 2017; 49: 365–380. 7. Ketan H and Nassir M. Aluminium hot extrusion process capability improvement using Six Sigma. Adv Prod Eng Manag 2016; 11: 59–69.

12 8. Boorla SM and Howard TJ. Production monitoring system for understanding product robustness. Adv Prod Eng Manag 2016; 11: 159–172. 9. Mu CH, Xie J, Liu Y, et al. Memetic algorithm with simulated annealing strategy and tightness greedy optimization for community detection in networks. Appl Soft Comput 2015; 34: 485–501. 10. Kuo CFJ, Su TL and Hung SJ. Process parameter optimization of scanning probe lithography for anodic oxidation of nanowire. Mater Manuf Process 2015; 30: 669–676. 11. Yu VF and Lin SY. A simulated annealing heuristic for the open location-routing problem. Comput Oper Res 2015; 62: 184–196. 12. Wu JY. Solving constrained global optimization problems by using hybrid evolutionary computing and artificial life approaches. Math Probl Eng 2012; 2012: 841410-1–841410-36. 13. Yang SH, Chung CC and Lee HS. Containership flag selection: the opening of direct shipping between Taiwan and China. Math Probl Eng 2014; 2014: 425808-1–425808-11. 14. Roy BK and Chakraborty S. Robust optimum design of base isolation system in seismic vibration control of structures under random system parameters. Struct Saf 2015; 55: 49–59. 15. Zhou GY, Wu E and Tu ST. Optimum selection of compact heat exchangers using non-structural fuzzy decision method. Appl Energ 2014; 113: 1801–1809. 16. Akin A and Saka MP. Harmony search algorithm based optimum detailed design of reinforced concrete plane frames subject to ACI 318-05 provisions. Comput Struct 2015; 147: 79–95. 17. Maheri MR and Narimani MM. An enhanced harmony search algorithm for optimum design of side sway steel frames. Comput Struct 2014; 136: 78–89. 18. Siuly S and Li Y. Designing a robust feature extraction method based on optimum allocation and principal component analysis for epileptic EEG signal classification. Comput Meth Prog Bio 2015; 119: 29–42. 19. Kim KS, Huha SH and Suganumab K. Effects of fourth alloying additive on microstructures and tensile properties of Sn–Ag–Cu alloy and joints with Cu. Microelectron Reliab 2003; 43: 259–267.

Advances in Mechanical Engineering 20. Das S, Acharya A and Pan I. Simulation studies on the design of optimum PID controllers to suppress chaotic oscillations in a family of Lorenz-like multi-wing attractors. Math Comput Simulat 2014; 100: 72–87. 21. Sarioglu IL, Czapnik B, Bostanci E, et al. Optimum design of a fuel-cell powertrain based on multiple design criteria. J Power Sources 2014; 266: 7–21. 22. Nandi AK, Deb K and Datta S. Genetic algorithm–based design and development of particle-reinforced silicone rubber for soft tooling process. Mater Manuf Process 2013; 28: 753–760. 23. Dai M, Tang D, Zheng K, et al. An improved geneticsimulated annealing algorithm based on a hormone modulation mechanism for a flexible flow-shop scheduling problem. Adv Mech Eng. Epub ahead of print 26 January 2015. DOI: 10.1155/2013/124903. 24. D’Urso G, Maccarini G, Quarto M, et al. Micro-electro discharge machining drilling of stainless steel with copper electrode: the influence of process parameters and electrode size. Adv Mech Eng 2016; 8: 1–16. 25. Su CT. Quality engineering. Taipei, Taiwan: Chinese Society for Quality, 2009. 26. Deng JL. Control problems of grey systems. Syst Control Lett 1982; 1: 288–294. 27. Chien CF. Decision analysis and management—a unison framework for total decision quality enhancement. Taipei, Taiwan: Yeh Yeh Book Gallery, 2015. 28. Hamada CA and Hamada MS. Illustrating the practice of statistics as a stepping stone to statistical engineering. Qual Eng 2012; 24: 184–192. 29. Adhikary DD, Bose GK, Jana DK, et al. Availability and cost-centered preventive maintenance scheduling of continuous operating series systems using multi-objective genetic algorithm: a case study. Qual Eng 2016; 28: 352–357. 30. Ko¨klu¨ U. Optimisation of machining parameters in interrupted cylindrical grinding using the Grey-based Taguchi method. Int J Comput Integ M 2013; 26: 696–702. 31. Zain AM, Haron H and Sharif S. Integration of simulated annealing and genetic algorithm to estimate optimal solutions for minimising surface roughness in end milling Ti-6AL-4V. Int J Comput Integ M 2011; 24: 574–592. 32. Shayeghi H and Hashemi Y. Application of fuzzy decision-making based on INSGA-II to designing PV– wind hybrid system. Eng Appl Artif Intel 2015; 45: 1–17.