Statistical Process Control For Monitoring Nonlinear Profiles A Six Sigma Project On Curing Process
Statistical Process Control For Monitoring Nonlinear Profiles A Six Sigma Project On Curing Process
Shing I. Chang , Tzong-Ru Tsai , Dennis K. J. Lin , Shih-Hsiung Chou & Yu-Siang
Lin
To cite this article: Shing I. Chang , Tzong-Ru Tsai , Dennis K. J. Lin , Shih-Hsiung Chou & Yu-
Siang Lin (2012) Statistical Process Control for Monitoring Nonlinear Profiles: A Six Sigma Project
on Curing Process, Quality Engineering, 24:2, 251-263, DOI: 10.1080/08982112.2012.641149
FIGURE 4 Patterns of the air temperature profiles during the cool-down stage: (a) the duration of vent open (X3), (b) the duration of
spray 1 (X4), (c) the duration of rest 1 (X5), and (d) the duration of spray 2 (X6).
making. Currently, the front-end databases are fed by they are not exactly random variables in the usual
the PLC. Our tasks involve the development of the sense because a PLC is used to make sure that each
middleware, SPC, and decision-making components cycle achieves the targets for QC2 and QC3; that is,
that will be integrated into the existing quality assur- 500 F for 480 time units. When the air temperature
ance database. reading gets close to 500 F, the PLC starts the timer
PH Corporation, like many other manufacturers, is for 480 time units to achieve the required fixed dur-
experiencing a modern-day factory phenomenon; ation. Based on the data recording tag shown in the
that is, a data tsunami. Lots of data are generated by database, the PLC would start the curing timer when
various operations, but their use in decision making a temperature reading first reaches 492 F. Each air
is ignored. The PLC that controls the curing process temperature profile exhibits a different pattern initi-
is capable of generating a large volume of data, most ally as shown in Figure 3. The temperature would
of which are simply stored, archived, and eventually eventually reach the steady state of 500 F. Therefore,
deleted. None or very little of it is transformed into the true quality characteristic for the curing time
useful information or statistics on which insightful requirement is how much time each profile stays at
decision makings can be made. There is really a miss- or above 500 F.
ing link between the statistical thinking=methods that We define X2 as the cumulative time that tem-
suggested collecting these data and engineering perature readings exceeding the target temperature
practice. In this study, we propose to integrate IT (500-d), where d is the standard deviation of the
and the proposed SPC procedure. steady-state temperature readings. Note that d is small
because of the PLC control actions. Based on the
steady-state segment of the temperature profiles
ESTABLISH MEASURE OF
shown in Figure 3, d is estimated to be 0.05 F. This
CRITICAL-TO-QUALITY quality characteristic can be obtained by an algorithm
CHARACTERISTICS written in MATLAB to process any profile during the
curing stage. It captures the key characteristics impor-
This measure step is the second step of the DMAIC
tant for the curing recipe; that is, QC2 and QC3.
process. Based on the CTQs QC1 to QC4, we would
like to establish statistics that can directly measure
the CTQs and be fed into control charts. It is straight- A SOLUTION TO PROFILE ANALYSIS
forward to quantify QC1 as X1, which is the duration
between the start of a curing stage and the start of the
DURING THE COOL-DOWN STAGE
countdown timer of QC3. This statistic is relatively We propose the use of profile monitoring techni-
easy to obtain via an algorithm written in MATLAB ques for tackling the SPC implementation issues for
(Manhattan, KS). QC4. Profile monitoring techniques have gained
It is a challenge to derive statistics based on QC2 much attention since Jin and Shi (1999) introduced
and QC3. A closer look at QC2 and QC3 reveals that their application to the force of a stamping process.
cycles. Specifically, cycles 1 and 98 are below the control limit of the Hotelling T2 control chart changes
mean curing duration but above the minimal 360 from 11.594 to 11.565 (based on Tracy et al. 1992).
time units. Another cycle 54 also belongs to this
group. Although this point falls within the IX chart STUDY DURING THE
control limits and is above 360 time units, it is far
from the mean point of the ellipsoid, causing it to
COOL-DOWN STAGE
be plotted outside the control limit of the T2 chart. The Hotelling’s T2 control chart for the cool-down
We keep the cycles of these two groups when estab- stage has four variables X3 to X6, which are the aver-
lishing the control limits of the T2 chart. age deviation statistics from the B-spline fitted nom-
The group of out-of-control points that cause con- inal profile. This is displayed in Figure 8. Cycles 41,
cerns consists of cycles 60, 73, 131, and 137. 136, and 141 are above the control limit of 15.764.
Although they reach the target temperature like the Figure 9 shows detailed plots of each segment with
majority of the other cycles, they fail to keep the these three profiles highlighted. In various segments,
minimal of 360 time units of curing time above the shapes of these cycles are very different from the
500 F. This result is forwarded to process engineers rest of the profiles. In addition, we can also use a
for more detailed study. After discarding the cycles matrix plot to show the correlation structure among
of this group from the Phase I data set, the new X3 to X6. As shown in Figure 10, the outliers corre-
spond to cycles 41, 136, and 141 in various places.
After removing these points, the control limit of the
Hotelling’s T2 control chart changes from 15.764 to
15.732.
prevent major process catastrophes from taking We now summarize how all of the statistical
place and provide useful information and statistics methods and tools are used in this study. The main
for informative decision making. process monitoring tools proposed are a pair of
FIGURE 9 Cool-down process profiles by (a) segment X3, (b) segment X4, (c) segment X5, and (d) segment X6 (out-of-control profiles 41,
136, 141).
Hotelling T2 control charts. The studies described operators would use the control charts to judge
above show how the control limits for both charts whether the process is in control or not. When a
are established based on the Phase I data set. During point plots outside the control limits of either one
the process monitoring phase, an algorithm written of the T2 control charts, process engineers are then
in MATLAB is used to obtain the quality characteris- notified to diagnose the process.
tics X1 to X6 from a sample temperature profile. Then
variables X1 and X2 form an input vector to be used
in the first T2 control chart as shown in Figure 6 and
IMPROVE THE CURING PROCESS
the rest of the variables are put into the second T2 During the fourth step of the DMAIC process, we
control chart as shown in Figure 8. If the points identify and generate potential solutions for process
plotted on both T2 control charts fall below their improvement. Currently air temperature profiles are
respective control limits, the process is deemed in not used for any process-related analysis or diag-
control. Otherwise, a biplot shown in Figure 7 is nosis. The PH Corporation will benefit from decision
used on variables X1 and X2 for fault diagnoses when making based on objective data analysis as demon-
the first T2 control chart triggered the out-of-control strated in this case study. Furthermore, during this
signal. Furthermore, a matrix plot of cool-down vari- Phase I study, we have identified a process adjust-
ables X3 to X6 should be used when the second T2 ment opportunity that may potentially save produc-
control chart triggered the out-of-control signal. tion cost. The distribution of the curing duration
Operators are responsible for monitoring both T2 above the target temperature reveals that it is poss-
control charts at the end of each curing cycle, and ible to trim the target duration to a smaller value as
quality engineers should use either a biplot or a shown in Figure 7. The curing recipe specifies 480
matrix plot for fault diagnoses when an out-of- time units of curing to accommodate all types of pro-
control signal is triggered by either or both of the ducts so that the actual curing time above 500 F may
T2 control charts. During a routine curing operation, be at the minimum of 360 time units. Most X2 values
FIGURE 11 Diagram of current process flow and the proposed embedded SPC flow. (Color figure available online.)