Performance Optimization, Prediction, and Adequacy by Response Surfaces Methodology with Allusion to DRF Technique

The RSM introduces statistically designed experiments for the purpose of making inferences from data. The second-order model is the most frequently used approximating polynomial model in RSM. The most common designs for the second-order model are the 3k factorial, Doehlert, Box-Behnken, and CCD. In this Box and Behnken design of three variables is selected as a representative of RSM and 70 : 30 polyester-wool DRF yarn knitted fabrics samples as a process representative. The survey reveals that secondorder model is the most frequently used approximating polynomial model in RSM.The Box-Behnken is the most suited design for optimization and prediction of data in textile manufacturing and this model is well-suited for DRF technique yarn knitted fabric. The trend was as higher wool fiber length shows higher fabric weight, abrasion, and bursting strength, correlation of TM was not visible; however, role of strands spacing is found dominant in comparison to other variables; at 14mm spacing it shows optimum behaviors. The optimum values were weight (gms/mt) 206 at length 75mm, TM 2.5 and 14mm spacing, abrasion (cycles) 1325 at length 70mm, TM 2.25 and 14mm spacing, bursting (kg/cm) 14.35 at length 70mm, and TM 2.00 and 18mm spacing. A selected variables, fiber length, TM, and strand spacing, have substantial influence. The adequacies of response surface equations are very high. The line trends of knitted fabric basic characteristics were almost the same for actual and predicted models. The difference (%) was in range of 1.21 to −1.45, 2.01 to −7.26, and 17.84 to −6.61, the accuracy (%) was in range of 101.45 to 98.79, 107.27 to 97.99, and 106.61 to 82.16, and the Discrepancy Factor (R-Factor) was noted to be 0.016, 0.002, and 0.229 for weight, abrasion, and bursting, respectively, between actual and predicted data. The L-estimation factors for actual and predicted data were that (i) the ratio were in range of 1.01 to 0.99, 1.02 to 0.93, and 1.22 to 0.94 for weight, abrasion, and bursting, respectively, (ii) the multiple-ratio was in range of 1.26 to 0.86, (iii) the ratio product was in range of 1.22 to 0.92, and (iv) the toting ratio was in range of 1.02 to 0.94.


Introduction
The response surface methodology (RSM) introduces statistically designed experiments for the purpose of making inferences from data.To achieve this goal, statistical considerations for preliminary planning of experiments, standard statistical designs for experiments, and underlying logic for using these designs are emphasized.It is a common but major error to view statistics as a tool to be used only after the experiments are completed.Even using their most sophisticated tools, researchers receiving data from improperly designed experiments can make only indistinct and approximate inferences.Therefore, it is unfortunate, because experimental data represent an expenditure of both time and money [1].
In general, the theoretical model that relates some controllable variables to a response either is not available or is very complex.Identifying and fitting from experimental data an appropriate response surface model requires some use of statistical experimental design fundamentals, regression modeling techniques, and optimization methods.As an important subject in the statistical design of experiments, RSM, introduced by Box and Wilson [2], comprises a group of mathematical and statistical techniques that is useful for empirical model building and analysis of problems in which a response of interest is influenced by several variables [3].

ISRN Textiles
It is important to fit a mathematical model equation in order to approximate a relationship between response and independent variables and determine the optimum settings of these variables that result in the maximum response.
The two important models that are commonly used in RSM, including the first-order model and second-order model [4], are as follows: where  is the response,  0 is the constant,   is the slope or linear effect of the factor   ,   is the quadratic effect of the factor   ,   is the interaction effect between the input factors   and   , and  is the residual term.The first-order models are inadequate to represent true functional relationships with independent variables.The second-order model is most suitable, highly structured, flexible, and diversified in order to locate the optimum point.
1.1.Designs for Fitting the Second-Order Model.The secondorder model is the most frequently used approximating polynomial model in RSM.The most common designs for the second-order model are the 3 factorial, Doehlert, Box-Behnken, and central composite designs (CCDs) [5,6].These symmetrical designs differ from one another with respect to their selection of experimental points, number of levels for variables, and number of runs and blocks.
1.2.3 Factorial Design.The 3 factorial design consists of all the combinations of the levels of the control variables with three levels each: low, medium or centre, and high [4].The number of experimental runs () required for this design is defined as  = 3, where  is the number of factors.The 3 factorial design needs a large number of experimental runs for large , which loses its efficiency in the modeling of quadratic functions.Therefore, a 3 factorial design is more appropriate having factors numbering less than five.Due to its requirement for more experimental runs it can usually be accommodated in practice; designs that present a smaller number of experimental points, such as Doehlert, Box-Behnken, and CCDs, are more often used.The application of 3 factorial design is not frequent, and the use of this design has been limited to the optimization of two variables, because its efficiency is very low for a higher number of variables.
1.3.Doehlert Design.The Doehlert (or uniform shell) design has been developed by Doehlert [7].The Doehlert design is for heterogeneous levels of variables.This property is important when some variables are subject to restrictions, such as cost and/or instrumental constraints, or when it is important to study a variable at major or minor levels.The intervals between each variable level must have a uniform distribution [6].The number of experiments required () for the Doehlert design is defined as  =  2 ++ 0 , where  is the number of factors and  0 is the number of centre points.For two variables, a central point surrounded by six points from a regular hexagon represents this design.For three variables, it is represented by a geometrical solid called a cub octahedron, and depending on how this solid is projected in the plane, it can generate some different experimental matrices.Although its matrices are neither orthogonal nor rotatable, it presents some advantages, such as requiring few experimental points for its application and high efficiency [8].
The number of experiments required () is given by  = 2( − 1) +  0 , where  is the number of variables and  0 is the number of central points.The design is represented as a cube and all points lie on a sphere of radius √ 2. In addition, this design does not contain any points at the vertices of the cubic region created by the upper and lower limits for each variable [10].The Box-Behnken design for three variables takes optimization with its 13 experimental points.This design is more economical and efficient in terms of the number of required runs than their corresponding 3 designs with 27 experiments.Therefore, this design is useful in avoiding experiments that would be performed under extreme conditions, for which unsatisfactory results might occur.However, it is ineffective for situations in which we would like to know the responses at extremes.The Box-Behnken design has been used for finding the optimum experimental conditions, leading to an optimal efficiency of different processes.[2] is the design most commonly used for fitting second-order models and it has been subjected to much attention in the theoretical development of its properties as in its practical use [10].This design combines a twolevel full or fractional factorial design with additional start points and at least one point at the centre of the experimental region.The CCD is widely used for the optimization of three variables.This design requires an experiment number according to  = 2+2+ 0 , where  is the number of factors and  0 is the number of central points.In CCD, all factors are studied in five levels.This  experiment is distributed as follows [4,10].

Central Composite Design. The CCD presented by Box and Wilson in 1951
(1) Full (or fractional) 2 factorial experiments, whose factors levels are coded as −1, +1: these experiments are the only points that contribute to the estimation of the two-factor interactions.
(3)  0 central points at (0, 0, . . ., 0): these experiments provide an estimation of pure error and contribute to the estimation of quadratic terms.The CCD is a rotatable and orthogonal design.A design is rotatable if the precision of the response estimation in all directions is equal and the orthogonality of the design means that different variable effects can be estimated independently.This design has been widely used for the optimization of several processes [11].
1.6.Optimization by Response Surface Methodology.In most production processes, the theoretical model that relates some controllable variables (factors) to a response either is not available or is very complex.In conventional methods used to determine this relationship, experiments are carried out varying systematically the studied parameter and keeping the others constant.This should be repeated for all the influencing parameters, resulting in an unreliable number of experiments.In addition, this exhaustive procedure is not able to find the combined effect of the effective parameters.In this way, the information about the relation between factors and response should be obtained in an empirical way [10,12].Using RSM, it is possible to estimate linear, interaction, and quadratic effects of the factors and to provide a prediction model for the response [13].
The textile industry is one of the largest and oldest industries worldwide and yarn manufacturing is the key process of it.The efficiency of yarn manufacturing depends on a number of factors, which are governed by the performance of fiber, yarn, and fabric initial characteristics and processing parameters of the experimental setup and also multiple pathways.Due to the complexity and variety of influencing factors, it is difficult to evaluate the relative significance of several affecting factors, especially in the presence of complex interactions [14].In the day by day innovations and introduction of latest technologies in yarn and fabric manufacturing, large numbers of textile scientists are developing so many advances.The development of double roving feed (DRF) techniques is one of them and widely accepted by the textile producers.The DRF yarn uses are increasing in the entire field including the knit-wears.
In the recent studies, only traditional one-factor-at-atime experiments were tested for evaluating the influence of operating factors on the DRF technique efficiency; however, very few researchers also used RSM.The DRF technique is not only time and work demanding but also completely lacks representation of the effect of interaction between different variables or factors.RSM allows an appropriate design of the experiments, which helps to decrease the number of runs.In addition, the modeling of the system facilitates the interpretation of multivariate phenomena and is valuable tool for scaling up [15].
The present endeavors reviewed the RSM techniques used for process optimization.The Box and Behnken design of three variables is selected as a representative of RSM.The DRF yarn knitted fabric production is chosen as a process for which the adequacy of the RSM is evaluated.Table 1.

Attachments to Produce DRF Yarn.
The following attachments are fitted in the conventional ring frame for the production of yarns by DRF technique: (1) rear roving guide; (2) double roving feeding attachments in drafting zone.

Experimental Design.
To study the individual and interactive effects of variables Box and Behnken factorial design was used for three variables.The following parameters are selected as prototype variables: (1) fiber length ( 1 ); (2) twist multiplier ( 2 ); (3) strand spacing ( 3 ).
Table 2 shows the coded and actual values of three parameters considered and fifteen sets of experimental combinations by DRF yarn and fabrics are knitted.

Measurement of Fabric Properties.
There are numbers of fabric properties that can be optimized by using this experimental approach; however, three fundamental properties are measured.The fabric weight per square meter (gms/mt 2 ) (weight) was evaluated by ASTM D-3776-79 method taking 10 × 10 cm sample from different places of knitted fabric.The abrasion cycle (abrasion) was measured by ASTM D-1966 method using martindale abrasion tester.The specimens were mounted on rectangular blocks of 1.5 × 2.5 inches with abrading material, which was itself fabric and then a number of rubs were counted by noting the number of cycles from counter in abrasion cycles.The bursting strength (bursting) is determined by ASTM D-3886-80 method using diaphragm type tester operated by hydrostatic pressure.The fabric samples are clamped by means of metal rings of internal diameter 30 ± 5 mm in the tester, by screwing the clamping ring too tight over the test piece.Thus, pressure was increased on the diaphragm until the test piece burst in between 7 and 20 seconds to increase pressure from zero to bursting point and then readings ware noted from the dial.

Development of Statistical Model.
To correlate the effects of variables and the response, the following secondorder standard polynomial was considered [16]: where  represents the responses and  0 ,  1 ,  2 , . . .,  23 are the coefficients of the model.The coefficients of main and interaction effects were determined by using the standard method.The response surface equations are calculated for prediction of responses.

Optimization of Fabric
Properties.The optimum fabric performance was predicted by using equations at all levels of variables drawing contours.
2.2.9.Adequacy of Models.The followings terms were studied for the adequacies of the models.
(a) Difference (%) is calculated by using the following equation: (b) Accuracy (%) is calculated by using the following equation: (c) Discrepancy Factor (-Factor) is calculated [17] by using the following equation: where  = Discrepancy Factor,   = actual values, and   = predicted values.

Results and Discussions
The actual observations, predicted values, and different calculated parameters for adequacy, response surface equations, and coefficient of correlation values are given in Table 3.The respective contours at different levels of variables were constructed and are given in Figures 1(a) to 1(c).The discussions are as follows.

Abrasion at Different Levels of Variables.
From Figure 1(b), as depicted in (1)(A), the trends were miscellaneous as at TM up to 2.25 and spacing near 14 mm the optimum is seen.(B) As TM and spacing increase the abrasion cycles are reduced to a certain limit.At TM and spacing levels 0 the optimum is found.(C) The trend was almost similar to (B).
As depicted in (2)(A), as length and spacing increase abrasion increases.The abrasion cycles are lowest at low fiber length.All parallel lines show similar trends in all spacings.(B) The parallel horizontal lines show that optimum could not be found in the range and the fiber length increases at all spacings as abrasion increases.(C) The trend is almost similar to (A) as length decreases and spacing increases from the decrease of the abrasion cycles.All lines are parallel and showing the same relationship at each spacing and length.
As depicted in (3)(A), as TM increases and length decreases the tendency of reduction in abrasion cycles is noted.(B) As length decreases and TM is below 2.25, the abrasion reduces.After TM is 2.25 as length increases, abrasion decreases.(C) As the length decreases, TM increases, and abrasion reduces, the trend is miscellaneous.(3) Weight at spacing (mm) (A) 10, (B) 14, and (C) 18     The fabric performance is proportional to the characteristics of fiber, yarn, and knit structures.In the study selected variables mainly have an impact on yarns and yarns are associated with knitted fabric.

Bursting at Different Levels of Variables. From
In DRF spinning as fiber length increases more length traps in drafting results, a more compact yarn.At optimum TM, the emerging fibers from front roller nip are trapped at higher binding force and gain better packing density because the packing density weight per unit length of yarn is higher.
In higher spacing, convergence angle is greater, which generates higher false twist as strand results in more compact trapping of surface fibers and ultimately more compact yarn structure.In DRF production after optimum TM, there is comparatively loosely packed yarn which shows less weight per unit length of fabric.Also after or before optimum spacing between strands, the trapping of fiber reduces causing loose structure of yarn and ultimately as fiber length increases the weight of fabric decreases, also the abrasion and bursting strength are reduced.With the increase of TM and strand spacing up to a certain limit (optimum condition of spinning), weight, abrasion, and bursting strength increase; however, after or before optimum condition adverse behaviors are seen.The probable reason may be that, up to a certain limit, it helps in better insertion of twist in single strand which causes better trapping of fibers in the yarn periphery and thereby improvement in various properties of yarn and respective fabrics.
In other studies, almost similar findings were reported by Ghasemi and other workers that wool/polyester blended worsted yarn is successfully produced by feeding two roving in spinning system and that yarns' specifications such as tensile strength, elongation, and abrasion resistance remain almost unchanged [18].
The literature reveals that as fiber length, TM, or strand spacing reaches the optimum value, the yarns produced are more compact due to better binding.The numbers of fibers in unit length are higher.The numbers of fibers present in the yarn structure are directly proportional to the ultimate product, that is, knitted fabric.The weight of unit area also increases or decreases, respectively, [19] due to fiber trappings.Similar results are admitted by other researchers that in same manufacturing conditions the breaking strength, tearing strength, abrasion resistance, and crease recovery properties of fabric are improved in case of DRF yarn than plies yarn [20].Figure 2 Figure 2(c) depicts that bursting strength trend of predicted values is different from the actual values.The actual bursting strength values were maximum at 14.35, minimum at 12.00, and average at 13.09; however, predicted bursting strength values were maximum at 12.90, minimum at 11.33, and average at 12.12, respectively.
The coefficients of correlation ( 2 ) between observed and predicted values were 0.96, 0.71, and 0.86 for weight, abrasion, and bursting, respectively, which shows significant influence.
The Discrepancy Factor (-Factor) was noted to be 0.016, 0.002, and 0.229 for weight, abrasion, and bursting, respectively.
The values under -estimation are as follows.
The values of ratio were maximum at 1.01, 1.02, and 1.22, minimum at 0.99, 0.93, and 0.94, and average at 1.00, 0.98, and 1.08 for weight, abrasion, and bursting, respectively.
The multiple-ratios were calculated maximum at 1.26, minimum at 0.86, and average at 1.05.
The values of ratio products were calculated maximum at 1.22, minimum at 0.92, and average at 1.06.
The values of toting ratio were calculated maximum at 1.02, minimum at 0.94, and average at 0.98.Box-Behnken is the most suited design for optimization and prediction of data in textile manufacturing and this model is well-suited for DRF technique yarn knitted fabric.

Conclusions
(2) The higher wool fiber length shows higher fabric weight, abrasion, and bursting strength.
(3) The correlation of TM is not visible.

3. 1 .
Weight at Different Levels of Variables.From Figure 1(a), as depicted in (1)(A), TM decreases from 2.25 and spacing increases as the weight reduces.The trend is miscellaneous.(B) The trend is miscellaneous; TM and spacing increase as the weight decreases consistently.(C) Up to 14 mm spacing, as TM decreases from 2.25 weight decreases; however, at slightly above 2.25 TM, spacing increases as the weight also increases.As depicted in (2)(A), two trends were found: first increasing weight as spacing increases in fiber length above coded level + 0.5 and second decreasing weight as length decreases from 70 mm.(B) As length decreases from 70 mm and spacing increases to 18 mm weight decreases.In length above 70 mm as spacing increases, weight increases and below reverse trend is visible.(C) In 70 to 75 mm in length as spacing increases, weight is reduced; however, when length is from 65 to 70 mm weight increases.As depicted in (3)(A): TM decreases and length increases up to 70 mm, while weight decreases; however, in further increases in length weight increases.(B) As TM decreases and length is up to 70 mm weight decreases; however, in further increases in length weight increases.(C) Optimums were found at TM 2.25 and fiber length 70 mm.
Figure 1(c), as depicted in (1)(A), as spacing and TM increase up to a certain TM bursting increases.(B) As TM increases and

3. 4 .
Adequacy of Models.The comparative analysis between actual and predicted performances of DRF yarn knitted fabric is shown by line diagram as given in Figures 2(a) to 2(c).The discussions are as follows.

( 1 )
The second-order model is the most frequently used approximating polynomial model in RSM.The

( 4 )
The role of strands spacing is dominant in comparison to other variables; at 14 mm spacing it shows optimum behaviors.(5) The optimum were weight (gms/mt 2 ) 206 at length 75 mm, TM 2.5 and 14 mm spacing, abrasion (cycles) 1325 at length 70 mm, TM 2.25 and 14 mm spacing, bursting (kg/cm 2 ) 14.35 at length 70 mm, and TM 2.00 and 18 mm spacing.

Table 2 :
Actual and coded values for independent variables and experimental design.

Table 3 :
Comparison of actual and predicted values.
(a) depicts that trend is almost similar in actual and predicted weight per square meter.The actual weight per square meter was maximum at 206, minimum at 184, and average at 192.27; however, predicted weight per square meter was maximum at 203.50, minimum at 184.54, and average at 192.95, respectively.Figure 2(b) depicts that abrasion cycles trend of predicted values is different from the actual values up to maximum values; therefore, it remains almost the same.The actual abrasion cycles were maximum at 1325, minimum at 1150, and average at 1236.67; however, predicted actual cycles were maximum at 1315.21, minimum at 1169.16, and average at 1263.78, respectively.