Art[,' c ""j) 01/1 /

Similar documents
5.1 Introduction. Learning Objectives

Chapter 5: Methods and Philosophy of Statistical Process Control

Example 1: One Way ANOVA in MINITAB

Equation 1: F spring = kx. Where F is the force of the spring, k is the spring constant and x is the displacement of the spring. Equation 2: F = mg

Robert Jones Bandage Report

Announcements. Lecture 19: Inference for SLR & Transformations. Online quiz 7 - commonly missed questions

Dwell Time and Immersion Depth Optimization For Immediate Defect Reduction

Running head: DATA ANALYSIS AND INTERPRETATION 1

LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/12

Analysis of Variance. Copyright 2014 Pearson Education, Inc.

Statistical Analysis of PGA Tour Skill Rankings USGA Research and Test Center June 1, 2007

LOW PRESSURE EFFUSION OF GASES adapted by Luke Hanley and Mike Trenary

Citation for published version (APA): Canudas Romo, V. (2003). Decomposition Methods in Demography Groningen: s.n.

A Method Quantitatively Evaluating on Technical Progress of Students in Ship Handling Simulator Training ABSTRACT

The Discussion of this exercise covers the following points:

Exercise 3. Power Versus Wind Speed EXERCISE OBJECTIVE DISCUSSION OUTLINE. Air density DISCUSSION

Legendre et al Appendices and Supplements, p. 1

Thermal Profiling the Reflow Process

y ) s x x )(y i (x i r = 1 n 1 s y Statistics Lecture 7 Exploring Data , y 2 ,y n (x 1 ),,(x n ),(x 2 ,y 1 How two variables vary together

Importance of Wave Height Measurement in Wave Solder Process Control

Novel empirical correlations for estimation of bubble point pressure, saturated viscosity and gas solubility of crude oils

Design of Experiments Example: A Two-Way Split-Plot Experiment

Fail Operational Controls for an Independent Metering Valve

Minimal influence of wind and tidal height on underwater noise in Haro Strait

Characterizers for control loops

ROSE-HULMAN INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering. Mini-project 3 Tennis ball launcher

SPIRIT III Radiometer Saturation Effect

Gerald D. Anderson. Education Technical Specialist

Statistical Process Control Basics. LeanSix LLC

A Depletion Compensated Wet Bath Simulator For Calibrating Evidential Breath Alcohol Analyzers

WATER OIL RELATIVE PERMEABILITY COMPARATIVE STUDY: STEADY VERSUS UNSTEADY STATE

M. Hussain, J. Gu, R. Engel, and D. Shortt *General Motors Corporation, GMPT Headquarters, Pontiac, Michigan USA

Section I: Multiple Choice Select the best answer for each problem.

Unit 4: Inference for numerical variables Lecture 3: ANOVA

INSTITUTE AND FACULTY OF ACTUARIES. Curriculum 2019 AUDIT TRAIL

Linear Compressor Discharge Manifold Design For High Thermal Efficiency

Naval Postgraduate School, Operational Oceanography and Meteorology. Since inputs from UDAS are continuously used in projects at the Naval

Impact of Lead Free on Automated X-Ray Inspection

WMO LABORATORY INTERCOMPARISON OF RAIN INTENSITY GAUGES

Calculation of Trail Usage from Counter Data

Investigating the effects of interchanging components used to perform ripple assessments on calibrated vector network analysers

Wing-Body Combinations

29 Pressure, Temperature relationship of a gas

Improving the Serving Motion in a Volleyball Game: A Design of Experiment Approach

Competitive Performance of Elite Olympic-Distance Triathletes: Reliability and Smallest Worthwhile Enhancement

International Journal of Technical Research and Applications e-issn: , Volume 4, Issue 3 (May-June, 2016), PP.

Hardware Triaxial and Consolidation Testing Systems Pressure Measurement and Control

The Simple Linear Regression Model ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD

Lab 4: Root Locus Based Control Design

CONTROL VALVE WHAT YOU NEED TO LEARN?

Journal of Applied Fluid Transients, Vol 1-1, April 2014 (3-1)

S-CO 2 Brayton Recompression Loop Design and Control

CORESTA RECOMMENDED METHOD N 6

Essential ski characteristics for cross-country skis performance

Time Pressure Dispensing

Exercise 2-2. Second-Order Interacting Processes EXERCISE OBJECTIVE DISCUSSION OUTLINE. The actual setup DISCUSSION

Midterm Exam 1, section 2. Thursday, September hour, 15 minutes

INTERACTION BETWEEN WIND-DRIVEN AND BUOYANCY-DRIVEN NATURAL VENTILATION Bo Wang, Foster and Partners, London, UK

AC : MEASUREMENT OF HYDROGEN IN HELIUM FLOW

You should know how to find the gradient of a straight line from a diagram or graph. This next section is just for revision.

Week 7 One-way ANOVA

Robust specification testing in regression: the FRESET test and autocorrelated disturbances

Chapter 12 Practice Test

A COMPARATIVE STUDY OF MIX FLOW PUMP IMPELLER CFD ANALYSIS AND EXPERIMENTAL DATA OF SUBMERSIBLE PUMP

Tip Seal Behavior in Scroll Compressors

SHOT ON GOAL. Name: Football scoring a goal and trigonometry Ian Edwards Luther College Teachers Teaching with Technology

Hydraulic and Economic Analysis of Real Time Control

Policy Management: How data and information impacts the ability to make policy decisions:

An Impeller Blade Analysis of Centrifugal Gas Compressor Using CFD

Design of a Microcontroller-Based Pitch Angle Controller for a Wind Powered Generator

Valve Losses in Reciprocating Compressors

Hydronic Systems Balance

COMPUTER-AIDED DESIGN AND PERFORMANCE ANALYSIS OF HAWT BLADES

Queue analysis for the toll station of the Öresund fixed link. Pontus Matstoms *

Tightening Evaluation of New 400A Size Metal Gasket

ARA PROJECT Prepared by: Peter T. Dzwilewski Applied Research Associates, Inc Shaffer Parkway Littleton, Colorado 80127

E. Agu, M. Kasperski Ruhr-University Bochum Department of Civil and Environmental Engineering Sciences

E.I. Kugushev, 7.6. Jaroshevskij Institute of Applied Mathematics, the USSR Academy of Sciences, Moscow, A-47, Miusskaya Sq», 4

Habit Formation in Voting: Evidence from Rainy Elections Thomas Fujiwara, Kyle Meng, and Tom Vogl ONLINE APPENDIX

Power Prediction Model for Milling 618 Stainless Steel Using Response Surface Methodology

IDENTIFYING SUBJECTIVE VALUE IN WOMEN S COLLEGE GOLF RECRUITING REGARDLESS OF SOCIO-ECONOMIC CLASS. Victoria Allred

Automating Injection Molding Simulation using Autonomous Optimization

Aalborg Universitet. Published in: Proceedings of Offshore Wind 2007 Conference & Exhibition. Publication date: 2007

ESTIMATION OF THE DESIGN WIND SPEED BASED ON

Grade: 8. Author(s): Hope Phillips

Impact Fatigue on Suction Valve Reed: New Experimental Approach

Measurement of Representative Landfill Gas Migration Samples at Landfill Perimeters: A Case Study

Denise L Seman City of Youngstown

Procedure 1: Volume vs. Pressure 1.) Using the lap tops, go to the Physics Education Technology from the University of Colorado at:

Lecture 22: Multiple Regression (Ordinary Least Squares -- OLS)

INCLINOMETER DEVICE FOR SHIP STABILITY EVALUATION

Accurate Measurement of Steam Flow Properties

PRESENTS. Solder & Oven Profiles Critical Process Variables

Paper 2.2. Operation of Ultrasonic Flow Meters at Conditions Different Than Their Calibration

Influence of Fly Ash Content on Compaction Characteristics of Fly Ash Clay Mixture

Volume-to-Capacity Estimation of Signalized Road Networks for Metropolitan Transportation Planning

INSTRUMENT INSTRUMENTAL ERROR (of full scale) INSTRUMENTAL RESOLUTION. Tutorial simulation. Tutorial simulation

EEF. Fatigue Testing Unit PROCESS DIAGRAM AND UNIT ELEMENTS ALLOCATION. Engineering and Technical Teaching Equipment

Experiment. THE RELATIONSHIP BETWEEN VOLUME AND TEMPERATURE, i.e.,charles Law. By Dale A. Hammond, PhD, Brigham Young University Hawaii

Available online at ScienceDirect. Energy Procedia 53 (2014 )

Transcription:

J-fhf.sfJrJ /~ - Ir pu--<vzy>c..f, IIMT-4T ~. JZ~l- ~7'cn-td C7' Art[,' c ""j) 01/1 / APPLICATION OF DESIGN OF EXPERIMENT IN PRINTED CUIRCUIT BOARD ASSEMBLY INDUSTRY V.Mathaeswarit, U.Safian 2, & A.Zalila 3 ICelestica Malaysia Sdn. BM, Kulim, Kedah,Malaysia 2,3School ofmathematical Sciences, Universiti Sains Malaysia, 11800 Minden Malaysia Imlingam@celestica.com;2safian@cs,usm,my; 3zalila@cs.usm.my Abstract Printed circuit board assembly (PCBA) is a process of mounting and inserting components to printed circuit board (PCB) by means of several principle processes such as Manual Insertion, Auto Insertion (Axial & Radial), Wave Soldering, Surface Mount Technology and Box Build. Manufacturers of PCBA are forced to produce better quality with less cost in order to stay in the business as globalization taking shape. To achieve such a goal, manufacturers are forced to strive to attain and maintain high first pass-yield. Hence, Design of Experiments (DOE) is one of the statistical tools that can be used by the manufacturers of PCBA to purposefully evaluate and optimize designs and processes; so that the manufacturers can achieve greater first pass yield and lower wastage which attributes to decrease the product cost. A project was implemented in a PCBA manufacturing company using the Full Factorial Design with two levels for the wave soldering process with the aid of MINITAB Software. And the objectives of the experiment are to determine the key factors of the process and optimize the process for higher yield. The results of the experiments were excellent, the yield at the wave soldering increased from 9 to 97.4. At the same time the key factors were identified and monitored closely to minimize variations. Keywords: printed cuircuit board, full factorial design, wave soldering; 1. Introduction Design of Experiments (DOE) is the statistical methodology of planning investigations and research, collecting data and data analysis. DOE identifies the key inputs or factors of a process and these factors are consciously changed to observe the consequences or outputs. A process is series of activities that convert inputs such as man, material, machine, method, etc. into the desired stage or outputs. The main reasons for implementing DOE are to: understand and reduce variability, determine optimum operating conditions or setting, perform a comparative study, ascertain the most influential inputs on response or process, acquire knowledge of the system, reduce wastage and save costs, reduce product development cycle (see Tony et ai., 2004).

Mathaeswari,Safian,Zaliia A good yardstick to identify the process stability is the Statistical Process Control (SPC) technique. Since, before an experiment can be carried out the process variations got to be addressed. If the process is in control, the performance is predictable over time, then the consequence can be determined by varying the principle factors of the process.the common reasons for poor experimental results: too much noise, incapable or unstable metrology, undefined scope, perplexed effects, interference with uncontrolled effects, inadequate sample size, inappropriate factor setting (Harriet et ai., 2003, Stephen et al., 1999). For full factorial design, the number of runs (n) for each combinations of the k factors is given by n = 2 k. The effect of a factor or interactions of factors is the average of all the data when the term is at the high level minus the average of all the data when the term is at the low level i.e, Effect = Min Value Response high - Min Value Response low 2.0 DOE Application The project presented in this dissertation was undertaken in Celestica (M) Sdn. Bhd. which is a Multinational company situated in the Kulim High Tech Park. The company is a PCBA contract manufacturer. Design of Experiment (DOE) was carried for the wave soldering process of a PCBA line and the objectives of the experiment were to: determine the key factors of the process, optimize the process for higher yield, reduce the conversion time. A team of people who have technical knowledge of the wave soldering machine and the processes carried out the project. 4 3 16 11 7 6 8 5 9 12 10 15 13 3 4 5 6 7 8 9 10 11 12 13 14 15 16 97.3 96.4 93.3 89.4 97.7 92.3 89.9 91.7 90.2 93.9 97.5 95.2 90.5 91.4 Table 2.1 : Results of the Experiments The team identified 3 key factors that may effect process i.e. machine conveyor speed, main wave pump and pre-heater. The experiment was carried out based on the 2

Proceeding of The IRCMSA 2005 combinations given in table 2.1. For each run, 100 units of Printed Circuit boards were evaluated and the corresponding yield (results) were recorded In respect to the results of analysis in figure 2.1, the P values for Main Effects in the Analysis of Variance (ANOVA) table for Yield is 0.000 which is less than ex = 0.05, where ex (alpha) is the significance level used for this experiment. From the ANOVA table the interaction effects are not significant. However, from the Effects and Coefficients for Yield table suggest that Main Wave, Pre-heat and Main Wave * Pre- heat (interaction) are significant using ex = 0.05. Factorial Fit: Yield versus Speed, Main Wave, Pre-heat Estimated Effects Term Constant Speed Main Wav Pre-heat Speed*Main Wav Speed*Pre-heat Main Wav*Pre-heat Speed*Main Wav* Pre-heat and Coefficients Effect Coef 93.150-0.125-0.062-0.125-0.400 Analysis of Variance for Yield Source DF Seq SS Main Effects 3 100.065 2-Way Interactions 3 11.265 3-Way Interactions 1 0.160 Residual Error 8 10.930 Pure Error 8 10.930 Total 15 122.420 for Yield SE Coef (coded units) Adj SS Adj MS 100.065 33.3550 11.265 3.7550 0.160 0.1600 10.930 1.3662 10.930 1.3662 (coded T 318.77-0.21 3.25-7.91-0.21-8 -2.78-0.34 F 24.41 2.75 0.12 Estimated Coefficients for Yield using data in uncoded units Term Coef Constant 44.156 Speed -52.281 Main Wav 0.164062 Pre-heat 0.140937 Speed*Main Wav 0.154688 Speed*Pre-heat 0.212500 Main Wav*Pre-heat -0.00051562 Speed*Main Wav*Pre-heat -0.00062500 Figure 2.1 : Minitab window outputfor Data Analysis units) P 0.000 0.836 0.836 0.513 111111 0.741.. P 0.112 0.741 3

Mathaeswari,Safian,Zalila Hence, the combination of those results, the following notion being interpreted: Main Wave has an effect of 1.900 which means that high factor setting i.e. rpm results in a better yield as compared to the low factor setting i.e. rpm. The intensity of this effect is lower than Pre-heat but higher than of the interaction of Main Wave * Pre- heat. Pre-heat has the greatest effect i.e. - 4.625 as compared to the Main Wave and interaction of Main Wave * Pre- heat. Its effect is significant at low factor setting i.e. C. Main Wave * Pre- heat, the value of this interaction is -1.625, which means it has impact during the low factor settings. The Pareto in figure 2.2 (a) was used to compare the relative magnitude and statistical significance for both main and interaction effects. Based on the Pareto chart, the 3 effects identified earlier (i.e. Main Wave, Pre-heat and Main Wave * Pre- heat) passes the reference line. Besides the Pareto, the normal effects plot in figure 2.2 (b) also reinforces the notion that Main Wave, Pre-heat and Main Wave * Pre- heat are significant. 2.2 (a) PareloChart alhe Standll'dlzed Effects (~I8Yie1d,~",.05) 012345678 k Spood B: J.WnWaI/ C:~"'1I1 1.'.B ~~wav c: I f>nt.hlial 1.0 0.' j.. 0.0 ~....Be -1.0.e.1.5 B., 0 SWldardized Effecl Figure 2.2: Pareto Chart & Normal Probability Plot of the Standardized Effects 2.1 Fitting a Reduced Model Next, we fitted a reduced model; for the terms which were identified significant in earlier analysis. The results for the reduced model are shown in figure 2.3. In the ANOVA table for yield, the P values for both Main Effects and Two-Way interactions are both lower than the ex value of 0.05. Therefore, the ANOVA suggests that both Main Effects and Two-Way interactions are significant. Again in the "Estimated Effects and Coefficients for Yield (coded units)" table Main Wave, Pre-heat and Main Wave * Pre-heat read P value which is inferior to the ex value of 0.05. This means that these factors are significant. The reduced model is further checked by using Residual Histogram for Yield - Figure 2.4 (a), Normal plot of Residuals for Yield - Figure 2.4 (b), Residuals vs. Fits for Yield - Figure 2.4 (c), and Residuals vs. Order for Yield - Figure 2.4 (d) 4

Proceeding of The IRCMSA 2005 Factorial Fit: Yield versus Main Wave, Pre-heat Estimated Effects and Coefficients for Yield (coded units) Term Effect Coef SE Coef T P Constant 93.150 0.2485 374.87 0.000 MainWav 1.900 0.950 0.2485 3.82 Pre-heat -4.625-2.312 0.2485 Main Wav*Pre-heat -1.625-0.812 0.2485 Analysis of Variance for Yield (coded units) Source DF Seq SS Adj SS AdjMS F P Main Effects 2 100.002 100.002 50.0012 51 2-Way Interactions 1 10.562 10.562 10.5625 19 Residual Error 12 11.855 11.855 0.9879 Pure Error 12 11.855 11.855 0.9879 Total 15 122.420 Estimated Coefficients for Yield using data in uncoded units Term Constant Main Wav Pre-heat Main Wav*Pre-heat Coef 2.3312 0.287813 0.310938-0.00101563 Figure 2.3 : Minitab window outputfor Reduced Model The Histogram in figure 2.4 (a) is approximately normally distributed, therefore it is basically can be interpreted as that the data do not contain any major unusual observations. The normal plot in figure 2.4 (b) was used to determine the distribution of the residuals. The points are very close to in forming a straight line, which means that the residuals are normally distributed. Residuals versus fits in figure 2.4 (c) illustrates there is no fixed pattern and scatted randomly in the plot, therefore it can be assumed that there is no unusual factor and it has non-constant variance. The residual versus order plot in figure 2.4 (d) is randomly scatted closed to the zero line and no specific pattern detected. This pattern indicates that there are no systematic effects of time or data collection order. 5

Mathaeswari,Safian,Zalila 2.4 (a) Histogram of the Residuals (IBSpOIVl&is Y"'ld) 2.4 (b) Normal Probability Plot of the Residuals {'~eisyillid) j, " ~ : z 2.4 (c) Residuals Versus the Fitted Values (IHpoIISeillYIIIll) R... 2.4 (d) ResIduals Versus the Order of the Data {lftpomuiliv-...,. " J. ----------'--------~ ~, J FIUedVaIue " ObMrvalIonOrder Figure 2.4 3.0 Drawing Conclusion Factorial plots for main effects (figure 3.I(a)) and interaction effects (figure 3.1(b)) were generated to compare the relative strength of the effects. From figure 3.1 it can be interpreted that Pre-heat has the highest impact since the line connecting the mean responses for Pre-heat C to c has a steeper slope. Followed by the Main Wave since there is a slope between Main Wave and. As for the Speed, the line connecting the mean responses from to 1.0 relatively parallel to X-axis, which interpreted as no main effect present... II 3.1 (a): Main effects plot 3.1 (b): Interaction plot (data means) for yield (data means) for yield........ I.' " #'.. i'., Spe<>d...-r-- 90.0.1.---- eo.8 Main Wave "'- 95.0 ~ 92.' ~.... -. ---+----..._._..- '. 95.0 Wl -- -- """'"...,38Q., ~ 91 Pte-heat 92.' 90.0 Figures 3.1(a) & 3.1 (b) 6

Proceeding of The IRCMSA 2005 On the other hand, from figure 3.1 (b) the interaction plot, only the interaction between Main Wave and Pre-heat has some degree of divergence, which suggests there is some effects due to the interactions of these factors. Other plots are relatively parallel which means there is no interaction present. Besides that, cube plot in figure 4.8 was drawn to illustrate the relationship between factors and the response. The interpretation derived from the cube plot is that Main Wave set at rpm and Pre-heat set at C produce the best yield. Nevertheless, varying speed between and 1.0 shows a very small improvement; however, in practice Speed may not exhibit any difference. Main Wave = Pre-Heat = Speed = Cube Plol (data means) for Yield Main Wave = Pre Heat = Speed = 1.0 MalnWave,! llo.80j-,...,... "'''",,'-' -V"-. 1.0 9O,~ao Figure 3.2: Cube Plot Based upon the above analysis and confirmation carried out (through out the experiment), the following mathematical model derived to identify the critical factors by using the coefficients: Yield =93.150 + 0.950 (main wave) - 2.312 (Pre-heat) - 0.812 (Wave * Pre-heat) Further to the above mathematical expression of the Yield, the epsilon squared (10 2 ) values (illustrated in figure 3.3) were calculated to assess the significance in terms of practical. Based on the epsilon squared values, pre-heat has the most impact, which controls 69.89 %; followed by main wave (11.80 %) and the interaction of pre-heat and main wave 8.63%. However, there is an error of 9.68 % which is contributed by other factors that were not covered in the experiment; but for practical application this value is within the acceptable level. Analysis ofvariance for Yield Source DF SS Pre-heat 1 85.563 Main Wav 1 14.440 Pre-heat*Main Wav 1 10.563 Error 12 11.855 Total 15 122.420 MS 85.563 14.440 10.563 0.988 F 86.61 14.62 19 P 0.000 0.002 0.007 Figure 3.3: Epsilon Squared values for the Reduced Model 7

Mathaeswari,Safian,Zaliia 4.0 Project Conclusion The Project described in section 2 was carried out in a Canadian company operating in Kulim High Tech Park. The aim of the project was to carry out Design of Experiment (DOE) to the wave soldering process of a PCBA line. With the objectives to determine the key factors of the process and optimizes the process for higher yield. Wave soldering process was selected because it was performing below the established target of98 percent for the Yield. The average yield from work week (ww) 30 to 41 is 9 percent; which is very far below the established target. The high rejection rate translated into high rework & repair cost and sometimes delays in delivery to customer. Yield (Actual vs Target) - Wave Soldering (SOLTEe 2) 10,000 8,000 6,000 4,000 2.000 100% 96% " 92% ~ 88% ~ ~ 84% 80% _Sal 62 8.750 8.970 8.Kl] 9,055 7.934 8,805 8.926 9.D} 8.021 5,034 5.201 6Ai06 _Sale>1 839 673 714 873 646 801 904 827 799 471 529 817 -...-.Ssles3 98,0% 98.0% 98.0% 98.0% 98.0% 98.0% 98.0% 98.0% 98.0% 98,0% 98.0% 98.0% """"""-Sala;,4 90.41% 90,27% 91.19% 90.36 91.87% 90,90% 89.87% 90.91% 90.04% 94 89.83% 90.40% Series 2: Total OutDut Series I: Total Reiect Series 3: Tamet Yield Series 4: Actual Yield Figure 4.1: Wave Solder Performance before DOE Figure 4.1 illustrates the performance of the wave soldering machine after design of experiments being carried out (data collected from work week 47 to 51). The average yield had improved from 9 percent to 97.4 percent; which is an improvement of 6.8 percent. The most crucial improvement was the reduction in number of PCBAs for rework; previously on the average 749 PCBAs due for rework every month and with improvement the average monthly quantity dropped to 119 PCBAs (which is approximately 84 percent reduction). This proves the potential of Design of Experiments in the PCBA industry. Although, the target of 98 percent yet to be achieved, but it has come to close to the target line. Besides that, the DOE enabled to identify the key factors that influence the wave soldering process. This knowledge is the crucial to the next path of continuous improvements. 8

Proceeding of The IRCMSA 2005 Yield (Actual vs Target) - Wave Soldering (SOLTEe 2) 10,000 T".... _..,,_...._.._........,,_ _ _.._ _ _ _ ~.. _..T 8,000 6,000 4,000 2,000 100% 98% CD 96% Cl ~ CD 94% I:! If 92% 90% IIIIIIIIIIISeries 1 _S erles 3 _8erles4 Serie~ 2: Totnl Ontnnt Serie~ 1: Totnl Reiect Serie~ 1: Tnmet Yiele! Serie~ 4: Actnnl Yiele! Figure 4.2: Wave Solder Performance after implementation ofdoe 5.0 Scope for Further Work While DOE proved to be effective in the PCBA industry but there are also other advanced tools like Response Surface Method. Which can used to examine the relationship between a response and a setof quantitative experimental variables or factors. Response surface method can be engaged to: find factor settings (operating conditions) that produce the "best" response, find factor settings that satisfy operating or process specifications, identify new operating conditions that produce demonstrated improvement in product quality over the quality achieved by current conditions, and model a relationship between the quantitative factors and the response Acknowledgements This paper is part of the first author's MSc dissertation. The second author would like to thank Universiti Sains Malaysia for the support and financial assistance. Such contributions are gratefully acknowledged. 9

Mathaeswari,Safian,Zalila References 1. Harriet, B. N. & Rene, V. V., "Integrating Experimental Design and Statistical Control for Quality Improvement", Journal of Quality Technology, 35 (2003), 403-406. 2. MINITAB Statistical Software, Release 13.2 by MINITAB Inc. 3. Stephen, R. S. & Robert G. L., " Understanding Industrial Designed Experiments", (Colorado: Air Academy Press, 1997). 4. Tony, L. & Belli, c., "Modifying 2 2 Factorial Designs to Accommodate a Restricted Design Space", Quality Engineering, 16 (2004), 99-104. 10