ESD.77 Lecture 18, Robust Design
Total Page:16
File Type:pdf, Size:1020Kb
ESD.77 – Multidisciplinary System Design Optimization Robust Design main two-factor interactions III Run a resolution III Again, run a resolution on effects n-k on noise factors noise factors. If there is an improvement, in transmitted Change variance, retain the change one factor k b b a c a c If the response gets worse, k n k b go back to the previous state 2 a c B k 1 Stop after you’ve changed C b 2 a c every factor once A Dan Frey Associate Professor of Mechanical Engineering and Engineering Systems Research Overview Concept Outreach Design to K-12 Adaptive Experimentation and Robust Design 2 2 n x1 x2 1 2 1 2 2 INT 2 2 2 1 x 2 ME (n 2) INT erf 1 e 2 1 n 2 Pr x x 0 INT dx dx 12 1 2 12 ij 2 1 2 2 2 1 2 0 x2 (n 2) INT ME INT 2 Complex main effects Methodology Systems A B C D Validation two-factor interactions AB AC AD BC BD CD three-factor interactions ABC ABD ACD BCD ABCD four-factor interactions Outline • Introduction – History – Motivation • Recent research – Adaptive experimentation – Robust design “An experiment is simply a question put to nature … The chief requirement is simplicity: only one question should be asked at a time.” Russell, E. J., 1926, “Field experiments: How they are made and what they are,” Journal of the Ministry of Agriculture 32:989-1001. “To call in the statistician after the experiment is done may be no more than asking him to perform a post- mortem examination: he may be able to say what the experiment died of.” - Fisher, R. A., Indian Statistical Congress, Sankhya, 1938. Estimation of Factor Effects Say the independent (bc) (abc) experilimental error o f (b) (ab) observations + (a)(), (ab), et cetera is σε. B (c) (ac) We define the main effect + - estimate Α to be - C (1) (a) - A + 1 (abc )()()()()()()(1)≡A [ ab + +ac+ − a− b− c− ] bc 4 The standard deviation of the estimate is 1 1 A factor of two improvement in σ=8 σ = 2 σ A 4 ε 2 ε efficiency as compared to “single question methods” Fractional Factorial Experiments “It will sometimes be advantageous deliberately to sacrifice all possibility of obtaining information on some points, these being confidently believed to be unimportant … These comparisons to be sacrificed will be deliberately confounded with certain elements of the soil heterogeneity… Some additional care should, however, be taken…” Fisher, R. A., 1926, “The Arrangement of Field Experiments,” Journal of the Ministry of Agriculture of Great Britain, 33: 503-513. Fractional Factorial Experiments + B + - - C - A + 3 1 2III Fractional Factorial Experiments Trial A B C D E F G FG=-A 1 -1 -1 -1 -1 -1 -1 -1 +1 2 -1 -1 -1 +1 +1 +1 +1 +1 3 -1 +1 +1 -1 -1 +1 +1 +1 4 -1 +1 +1 +1 +1 -1 -1 +1 5 +1 -1 +1 -1 +1 -1 +1 -1 6 +1 -1 +1 +1 -1 +1 -1 -1 7 +1 +1 -1 -1 +1 +1 -1 -1 8 +1 +1 -1 +1 -1 -1 +1 -1 27-4 Design (aka “orthogonal array”) Every factor is at each level an equal number of times (balance). High replication numbers provide precision in effect estimation. Resolution III. Robust Parameter Design Robust Parameter Design … is a statistical / engineering methodology that aims at reducing the performance variation of a system (i.e. a product or process) by choosing the setting of its control factors to make it less sensitive to noise variation. Wu, C. F. J. and M. Hamada, 2000, Experiments: Planning, Analysis, and Parameter Design Optimization, John Wiley & Sons, NY. Cross (or Product) Arrays Noise Factors Control Factors a -1 -1 +1 +1 b -1 +1 -1 +1 A B C D E F G c -1 +1 +1 -1 1 -1 -1 -1 -1 -1 -1 -1 2 -1 -1 -1 +1 +1 +1 +1 3 -1 +1 +1 -1 -1 +1 +1 27 4 4 -1 +1 +1 +1 +1 -1 -1 III 5 +1 -1 +1 -1 +1 -1 +1 6 +1 -1 +1 +1 -1 +1 -1 7 +1 +1 -1 -1 +1 +1 -1 8 +1 +1 -1 +1 -1 -1 +1 3 1 724 3 1 2III III2III Taguchi, G., 1976, System of Experimental Design. Step 1 Step 1 Summary:Step 5 Step 5 Summary: • Form cross function team of experts.• Determine noise factors and levels. Identify Project • Clearly define project objective. • Determine noise strategy • Assign Noise and Team Define roles and responsibilities to team - Surrogatemembers. Noise Strategy • TranslateFactors customer to Outer intent non-technical - termsCompound into technical Noise terms. • Identify productArray quality issues. - Treat Noise Individually • Isolate the boundary conditions and• describe Establish the outer system noise in arrayterms matrix.of its inputs and outputs. Step 2 Step 2 Summary:Step 6 Step 6 Summary: • Formulate Select a response function(s). • Cross functional team will need to develop a step-by-step plan to carry out the • Select a signal parameter(s). logistical activities necessary for successful completion of the data collection Engineered • Conduct Step 4 Summary: Determine if problem is static or dynamicphase parameter. of the optimization Note: Static experiment. has a single input System: Ideal and oneExperiment or more responses. and Dynamic• hasIdentify multiple Facility input Constraints. signals and multiple output Function / Quality signals. Collect Data • Determine logistical/run order. • Determine the S/N function. See section Step2 e.g. Nominal-the-best, etc. • Determine control factor levels Characteristic(s) • Calculate the DOF • Determine if there are any interactions Step 3 Step 3 Summary:Step 7 Step 7 Summary: • Select control factor(s). • Calculate the Following Values: Formulate • Rank control factors. • Mean • Select the appropriate orthogonal array Engineered • SelectAnalyze noise factors(s). Data and • Variance • Signal to Noise Ratio System: Select Optimal • ANOVA Parameters Design • Factor Effects • Interpret the Results • Select Factor Levels providing the largest improvement in S/N should be • Use the mean effect values to predict the new mean of the improved system • If needed, select a value for the scaling factor Step 4 Step 4 Summary: • DetermineStep control factor8 levels. Step 8 Summary: Assign Control • Calculate the DOF. • TBD - Canice • Determine if there are any interactions between control factors. Factors to Inner • Select thePredict appropriate and Orthogonal Array. Array Confirm Majority View on “One at a Time” One way of thinking of the great advances of the science of experimentation in this century is as the final demise of the “one factor at a time” method, although it should be said that there are still organizations which have never heard of factorial experimentation and use up many man hours wandering a crooked path. Logothetis, N., and Wynn, H.P., 1994, Quality Through Design: Experimental Design, Off-line Quality Control and Taguchi’s Contributions, Clarendon Press, Oxford. Minority Views on “One at a Time” “…the factorial design has certain deficiencies … It devotes observations to exploring regions that may be of no interest…These deficiencies … suggest that an efficient design for the present purpose ought to be sequential; that is, ought to adjust the experimental program at each stage in light of the results of prior stages.” Friedman, Milton, and L. J. Savage, 1947, “Planning Experiments Seeking Maxima”, in Techniques of Statistical Analysis, pp. 365-372. “Some scientists do their experimental work in single steps. They hope to learn something from each run … they see and react to data more rapidly …If he has in fact found out a good deal by his methods, it must be true that the effects are at least three or four times his average random error per trial.” Cuthbert Daniel, 1973, “One-at-a-Time Plans”, Journal of the American Statistical Association, vol. 68, no. 342, pp. 353-360. My Observations of Industry • Farming equipent company has reliability problems • Large blocks of robustness experiments had been planned at outset of the design work • More than 50% were not finished • Reasons given – Unforseen changes “Well, in the third experiment, we – Resource pressure found a solution that met all our needs, so we cancelled the rest – Satisficing of the experiments and moved on to other tasks…” More Observations of Industry • Time for design (concept to market) is going down • Fewer physical experiments are being conducted • Greater reliance on computation / CAE • Poor answers in computer modeling are common – Right model → Inaccurate answer – Right model → No answer whatsoever – Not-so right model → Inaccurate answer • Unmodeled effects • Bugs in coding the model Human Subjects Experiment • Hypothesis: Engineers using a flawed simulation are more likely to detect the flaw while using OFAT than while using a more complex design. • Method: Between-subjects experiment with human subjects (engineers) performing parameter design with OFAT vs. designed experiment. Treatment: Design Space Sampling Method Trial A B C D E F G 1 - - - - - - - • Adaptive OFAT 2 + - - - - - - 3 + - - - - - – One factor changes in 4 + - - - - each trial 5 + - - - 6 filled in as + - - 7 + - 8 required to adapt + Trial A B C D E F G • Plackett-Burman L8 1 - - - - - - - 2 + - + - - + + – Four factors change 3 + + + - + - - 4 - + - - + + + between any two trials 5 + - - + + - + 6 - - + + + + - 7 - + + + - - + 8 + + - + - + - Using a 27 system avoids possible confounding factor of number of trials. Increasing number of factors likely means increasing discrepancy in detection. Larger effect sizes require fewer test subjects for given Type I and II errors. Parameter Design of Catapult to Hit Target • Modeled on XPultTM • Commonly used in DOE demonstrations • Extended to 27 system by introducing – Arm material – Air relative humidity – Air temperature Control Factors and Settings Control Factor Nominal Setting Alternate Setting Relative Humidity 25% 75% • Control factor tied directly to Pullback 30 degrees 40 degrees simulation mistake Type of Ball Orange (Large-ball TT) White (regulation TT) Arm Material Magnesium Aluminum Launch Angle 60 degrees 45 degrees • Arm material selected for its Rubber Bands 3 2 moderate effect size Ambient 72 F 32 F Temperature 50 • Computer simulation equations 100 45 are “correct”, but intentional 40 80 mistake is that arm material 35 properties are reversed 30 60 Effect Size 25 Cumulative Percent 20 40 15 • Control factor ordering is not 10 20 random, to prevent variance due 5 to learning effect 0 0 • “Bad” control factor placed in 4th Pullback Type of Ball Arm MaterialLaunch Angle Relative Humidity column in both designs No.