Next Article in Journal
Synthesis and Properties of Moisture-Cured Reactive Polyurethane Containing Castor Oil and Oxime Compounds
Previous Article in Journal
Production of Biodegradable Palm Oil-Based Polyurethane as Potential Biomaterial for Biomedical Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence Mechanism and Optimization Analysis of Technological Parameters for the Composite Prepreg Tape Winding Process

1
Laboratory of Aero-engine High Performance Manufacturing, School of Mechanical Engineering, Northwestern Polytechnical University, Xi’an 710072, China
2
Xi’an Electronic Engineering Research Institute, Xi’an 710100, China
*
Author to whom correspondence should be addressed.
Submission received: 20 July 2020 / Revised: 9 August 2020 / Accepted: 13 August 2020 / Published: 17 August 2020

Abstract

:
Composite prepreg tape winding technology has proven to be an effective method for manufacturing revolving body composite structures in aerospace field. Process parameters including heating temperature, tape tension and roller pressure have an important impact on the winding products’ mechanical property such as tensile strength. The aim of this study is to investigate the influence mechanism and optimization analysis of parameters for the composite prepreg tape winding process. Firstly, the sensitivity analysis for single parameter had be employed to reveal the influence mechanism of each winding parameter change on tensile strength. Secondly, iso-surfaces analysis for parameter range had be applied to describe the distribution law of parameter with continuous distribution characteristics. Then the coupling analysis for process parameters was carried out employing response surface methodology. The analysis results showed that tape tension has the most significant effect on the winding products’ tensile strength. And the outstanding parameter combination with the heating temperature of 72 °C, tape tension of 307 N and roller pressure of 1263 N was provided by response surface design software via desirability function method. The validation experiments showed that the optimal parameter combination has a positive guiding significance for improving the quality of winding products.

Graphical Abstract

1. Introduction

Composite materials have been widely used in engineering applications given their high specific intensity and specific rigidity [1,2,3,4,5,6]. Composite prepreg tape winding—a process utilized to fabricate revolving body composite structures—is recognized as an effective method for producing aerospace products, such as solid rocket motor nozzles, satellite fairings and missile noses [7,8,9]. The winding process can be described as the prepreg tapes’ hoop winding process designed to fabricate the rotary body composite [10,11]. Under the pressure and heating of compaction roller, the tensioned prepreg tapes are steadily wound around the mandrel rotating at a constant speed. The wound product will be cured and demolded to become the final product. Studies have shown that the process parameters during the winding process, such as heating temperature, tape tension and roller pressure and winding speed have a significant effect on the mechanical property of the composite wound products [12,13,14,15,16,17]. As an important index to characterize the mechanical properties of wound products, hoop tensile strength is employed to reflect the ability of composite cylindrical parts to resist hoop tensile failure. To improve the hoop tensile strength of the wound products, the study of the influence mechanism and optimization analysis of technological parameters in the composite prepreg tape winding process is necessary.
In recent years, more and more scholars have made great contributions to the research of composite winding technology. Some new filament winding technologies and equipment have been appeared as the technology develops. For example, Sorrentino et al. [18] used the robotic filament winding to manufacture a complex shape structural part belonging to aeronautic field. Libonati et al. [19] provided a modified classic filament winding named squeeze-winding to fabricate the bone-like structure characterized by axial cylindrical features. With the advantages of scouring resistance and high temperature resistance, the fiber-reinforcement composite tape winding products are usually used for aerospace vehicle. Therefore, the physical or mechanical properties of the winding products need to be carefully studied. Such as Wang et al. [20] researched the void growth and fiber volume fraction in the filament winding process. And Eggers et al. [21] evaluated the mechanical response of filament winding composite rings under tension and compression. Moreover, the process parameters play an indispensable important role in improving the winding products’ various performances. In terms of process controlling, the control strategies for winding tension [22] and compaction pressure [23] are the main object of study. Some scholars have studied the influence of process parameter on the properties of an end winding product. According to Mazumdar et al. [12] and Mack et al. [13], several parameters such as heat intensity, winding velocity, roving tension, consolidation pressure, of tape and tool material and cooling rate can affect the properties of an end product. In Nath’s article [14], four control factors machine speed, roller pressure, tape tension and tape temperature were investigated for the tape winding process. Meanwhile, processing temperature, winding speed and compaction pressure were selected as the key parameters to study for tape winding process according to Dai et al. [15] and Zaami et al. [16]. However, three parameters including tape tension, roller pressure and temperature were analyzed for developing a NC tape winding machine from a research of Shi et al. [17]. In order to save experimental resources unfortunately, three main process parameters tape tension, roller pressure and heating temperature are considered in this manuscript. At last, a reasonable parameters combination will contribute to improving the performances of winding product. For the parameter optimization method, Nath [14] used Taguchi’s robust design methodology to optimize the winding parameters to enhance the performance of solid rocket nozzle. For another example, Colombo et al. [24] presented an analytical study for the optimization of the filament winding parameters of a plain pipe in glass reinforced polymers composite.
According to the multitudinous researches, we can clearly understand that parameter influencing, process controlling and parameter optimizing on the winding process have been widely studied. These scientific research results have actively promoted the process development and industrial application of composite winding technology. However, the influence mechanism of single factor, multi-factor coupling and parameter range of the process parameters in the composite prepreg tape winding process on the winding products lack a systematic summary.
The aim of this study is to investigate the influence mechanism and optimization analysis of parameters for the composite prepreg tape winding process. First of all, heating temperature, tape tension and roller pressure are selected as the key important process parameters. And then the hoop tensile strength will be selected as the mechanical property of winding products.
In our engineering practice, we devoted to developing a quality control system for the prepreg tape winding products. For a particular prepreg tape, we need to give a suggested parameter combination via theoretical model, empirical model and experimental verification. And in the actual manufacturing process, the types of prepreg tapes are usually limited. Our work in this article will provide an analysis process for how we studied and obtained the parameter combination.
At first, we intended to research importance degree and distribution characteristics of the single parameter in the winding process model. According to the studies of Nikishova et al. [25] and Saltelli et al. [26], sensitivity analysis can be used to identify the key parameters which affect model performance. It plays an important role in model parameterization, calibration, optimization and uncertainty quantification [27]. The iso-surface method is a method of computer graphics. In three-dimensional space, a set of curved surface graphs are fitted by points with the same magnitude in a certain physical quantity. Then the curved surface can be applied to describe the distribution law of those physical quantities with continuous distribution characteristics [28,29]. In the paper, sensitivity analysis for single parameter will be applied to research the influence mechanism of each winding parameter change on the tensile strength. Then the iso-surfaces analysis for parameter range will be applied to describe the distribution law of parameters with continuous distribution characteristics.
Then response surface methodology will be employed to carry out the coupling analysis for process parameters. The response surface methodology, as a collection of mathematical and statistical techniques, has been widely used to develop meta-models due to its low computational effort [30]. The method is effective for developing, optimizing and improving productive processes when many factors affect one or more particular response values [31,32].
In addition, the desirability function method has been widely applied to the multiple response optimization process for the advantage of simplicity [33,34,35]. The desirability function method can convert each response variable to an individual desirability function [33]. Here, we want to obtain the optimal parameter combination by the aid of the response surface design software based on desirability function method. At last, the authors wish these researches will be useful to enhance the winding products’ quality in the composite tape winding manufacturing process.

2. Composite Prepreg Tape Winding Process

Composite prepreg tape winding technology can be defined as a composite fabrication process by which resin-impregnated tapes are consolidated around a rotating mandrel shape with a geometric path and tension to create a structural solid of revolution [7,36]. Figure 1 gives a three-dimensional schematic of the composites prepreg tape winding process. In the tape winding process, the prepreg tapes are pulled and placed on a series of specific rollers from the tray as shown in Figure 1. Generally, the material tray will provide the prepreg tape with a certain pretension to ensure the smooth conveyance. At the same time, the magnetic powder brake provides the required tensile tension for the tape during the winding process. The force sensor will monitor the tension value in real time and feed back to the control system. Meanwhile, we can see that the high precision partial adjust mechanism has been adopted to maintain the correct winding trajectory instead of deviation. At first, the unwinding tapes need to be heated to the suitable temperature via the electric heating device from the inside of compaction roller. At the same time, the cylinder provided compaction roller pressure on the surface of tape. Subsequently, the unwinding prepreg tapes can be continuously winded to the outer layers of the mandrel with an invariable speed. During the composite tape winding process, the heating temperature will contribute to reducing the resin viscosity and improving the degree of interlaminar contact and positive force provided by a compaction roller is conducive to squeezing out the bubbles among the interlaminar contact interfaces and enhancing the contact degree of interlayers. In addition, a suitable tape tension can greatly strengthen the hoop tensile strength of tape winding products [36,37].

3. Experimental Design and Procedure

Tensile strength is a representative mechanical property index for the composite tape winding products. As the composite tape winding is a complicated process, an ideal theory model describing the objective law between winding process and technological parameters cannot be easily established. In this work, we were devoted to building a mathematical simulation between parameters and responses via response surface methodology. However, the simple linear functions are difficult to explain the behavior of the complex system. To reduce the cost of computing, in practice, we choose to select the second order polynomial regression to simulate and analyze the composite tape winding process.
In this paper, a second order polynomial regression model was employed to describe the changing rule of tensile strength with the process parameter. The tape winding process model of the influence mechanism of technology parameters on the tensile strength can be presented:
T S = f ( x 1 , x 2 , , x n ) = α 0 + i = 0 n α i x i + i = 0 , j = 0 n α i j x i x j + i = 0 n α i i x i 2 + ε
where TS is the tensile strength of winding product, n is the number of variables, xi is the ith winding process parameter; α0 is the constant term, αi represents the coefficients of the linear parameters, αij represents the coefficients of the interaction parameters and ε is the residual associated with the experiments.

3.1. Experimental Design

The response surface methodology is generally applied to experimental situations where several independent variables influence a response variable [34]. As a type of response surface methodology, the Box–Behnken design is rotatable second-order designs based on three-level incomplete factorial designs, allowing the number of design points to increase at the same rate as the number of polynomial coefficients [34,38]. In this study, the Box–Behnken design was chosen to find out the relationship between the response functions and variables [39]. Therefore, the three-factor and three-level Box–Behnken design based on response surface methodology theory was used to design this multi-factor and single-target experimental process.
In the experiments, the range of heating temperature T is between 40 and 100 °C, the tape tension F is between 100 and 500 N and the roller pressure P is between 500 and 2000 N. Meanwhile, the interval of process parameters in this paper came from an approximate range of experience and manufacture practice. To facilitate the latter part of comparative analysis process, we used a formula to converted each parameter into a coding system in the first place. Assuming that x1, x2 and x3 denotes the temperature, tension and pressure, respectively, and that βi1, βi0 and βi−1 represent the variable levels 1, 0 and −1 for each process parameter, respectively. Then the three variables can be described:
x i = β i j β i 0 Δ i   ( i = 1 , 2 , 3 ; j = 1 , 0 , 1 )
where xi is the code of the ith variable, βij denotes the jth variable level of the ith processing parameter, βi0 represents the 0 level of each processing parameter. Δi represents the interval range of the ith processing parameter, which can be calculated as Δi = (βi1βi0).
The test levels and coding results of processing parameters were combined and displayed in Table 1. The experiments were performed 17 times based on different parameters’ combinations.

3.2. Sample Preparation and Measurement

In the experiments, all the test specimens were fabricated on the Automate Tape Winding KUKA Robot (XGD-1200) (KUKA, Augsburg, Germany) independently designed and manufactured by Northwestern Polytechnical University, China. The compaction roller was made of 45 steel with the size of 160 mm in external diameter and 150 mm in width. The winding mandrel was a 45 steel cylinder with a size of 150 mm in external diameter and 1200 mm in length. In addition, a high precision partial adjust mechanism has been adopted to maintain the correct winding trajectory instead of deviation.
The materials used in the tests were carbon fiber-reinforced epoxy resin unidirectional prepreg tapes. The carbon fiber and epoxy resin were T-300 coming from Toray Industries, Inc. (Tokyo, Japan) and YH-69 provided by Wuxi Resin Factory of Bluestar New Chemical Materials Co., Ltd. (Wuxi, China), respectively. The width of the prepreg tape was 50 mm and the thickness was 0.2 mm.
All the winding specimens were manufactured by way of hoop winding. During the whole winding process, the ambient temperature was 20 ± 2 °C. And the mandrel rotated steadily at a speed of 10 revolutions per minute from beginning to end. Finally, the winding products were cured in the autoclave manufactured by TEDA Industrial Equipment Co., Ltd. (Tianjin, China). Meanwhile, the heating rate in the autoclave was strictly controlled at 2.5 °C/min. In addition, the curing temperature need to be consistently retained for 120 min after the temperature rose from room temperature to 150 °C. At the same time, the curing pressure must be kept around 0.15 MPa.
Since we had obtained the cured winding products, a split disc test method was developed for the determination of the apparent hoop tensile strength of composite winding ring-shaped specimens. As shown in Figure 2a, the standard testing rings for split disc tensile tests were turned along the tape winding direction on the universal lathe. Figure 2b displays the standard testing rings which have a same dimension of 150 ± 0.2 mm in inner diameter, 6 ± 0.2 mm in width and 3 ± 0.1 mm in thickness. The measurement for the hoop tensile strength of tape winding specimens were carried out following the standard GB/T 1458-2008 [40]. All the tensile strength tests were accomplished on the electronic universal testing machine DDL100 (Figure 2c) produced by Changchun Research Institute for Mechanical Science Co., Ltd (Changchun, China). According to the standard, the computational formula for the hoop tensile strength of composite tape winding ring can be given as
σ t = F b 2 b h
where σt denotes the tensile strength TS; Fb means the maximum load and b and h are the width and thickness of specimens, respectively.

3.3. Experimental Results and Analysis

Based on Table 1, the turning test was carried out 17 times. Figure 2d–f shows the typical images specimens after tensile failure. From overall view of standard testing ring, some fractures occurred locally just as the red box selection in Figure 2d-1. This means that a considerable weak part exists in the wound ring. By contrast, some specimens had a global decomposition-type failure mode as shown in Figure 2d-2. This is probably because the adhesive force between the fiber and the substrate in the corresponding specimen was relatively unsubstantial. As a result, more dispersion than fracture occurred among the fibers of specimen. In terms of the split form of the fiber bundle, a part of testing rings as shown in Figure 2e-1 had cracks along the ring direction and breakages in different locations. However, another part as shown in Figure 2e-2 had cracks in the radial direction and breakages in different positions. From the basic shape of the crack, some fractures were irregular and destroyed in an interlaced manner in Figure 2f-1. Still, some specimens’ fractures look comparatively flatter than the irregular ones as shown in Figure 2f-2.
All the test results of tensile strength for standard testing winding rings were recorded in Table 2. In the table, T is the heating temperature, F is the tape tension and P is the roller pressure. From the table, the maximum and minimum values of tensile strength are 1124.31 and 721.39 MPa, respectively.
The experimental factors and tensile strength were converted into matrix form by the test parameter transformation. The coefficients of Equation (1) were obtained from the Design Expert 10 (Stat-Ease, Inc., Minneapolis, MN, USA). Then the polynomial regression model of the influence mechanism of technology parameters on the tensile strength is given as follows:
Y T S = f ( x 1 , x 2 , x 3 ) = 1121.008 + 10.31 x 1 + 12.488 x 2 + 8.758 x 3 + 21.555 x 1 x 2 26.980 x 1 x 3 96.525 x 2 x 3 71.604 x 1 2 182.739 x 2 2 101.6889 x 3 2
To find out the significant influencing factors and reveal the interactions between the factors, the analysis of variance for samples was carried out based on the response surface design software Design Expert 10 (Stat-Ease, Inc., Minneapolis, MN, USA). The analysis results are shown in the Table 3. According to the analysis of variance, the F value is a main technical parameter to describe the impact of each input variable on the response value. In short, the larger F value, the more significant the corresponding variable is. From the Table 3, the F value of polynomial regression model is 2567.59, far exceeding the F-test critical value F0.05 (9, 7) = 3.677, which means the developed model is significant. Meanwhile, the value of ‘Prob > F’ is less than 0.05 meaning the corresponding model term is significant. So, in this model, T, F, P, TF, TP, FP, T2, F2 and P2 are all the significant model terms. Here, T, F and P represents the heating temperature, tape tension and roller pressure, respectively. At last, the ‘Prob > F’ value of ‘Lack of Fit’, 0.2749, means the lack of fit terms were not significant relative to the pure error. Therefore, the corresponding model can be considered a complete model including all the necessary items.
In addition, the residual analysis had been employed to validate the reliability of corresponding polynomial regression model. Figure 3a shows the diagram of internally studentized residuals. In the graph, the residuals scatters are an approximate linear distribution along a straight line. The result meets the requirement of normal distribution meaning that the fitted model for tensile strength is trustworthy. According to the residual operation diagram in Figure 3b, the residuals scatters is closer to the random distribution pattern. Meanwhile, Figure 3c gives the comparison between the actual and predicted value of the tensile strength. The max error between the actual and predicted value is less than 5%. Therefore, the developed polynomial regression model can be considered as reasonable for the composite prepreg tape winding process. In brief, the established second order polynomial regression model could be effectively utilized to predict the tensile strength of winding products.

4. Influence Mechanism and Optimization of Parameters

4.1. Sensitivity Analysis for Single Parameter

Parameters sensitivity denotes the sensitive degree of target value to each design variable change [25,26]. In this study, sensitivity analysis for single parameter were applied to make clear the influence mechanism of each winding parameter changes on tensile strength. In a mathematical sense, the sensitivity coefficient is the first-order derivative of model output with reference to the model parameter [41]. Therefore, let S(xi) denote the sensitivity coefficients, we can obtain the sensitivity of single parameter as follows:
S ( x i ) = l i m Δ x i 0 f ( x i + Δ x i ) f ( x i ) Δ x i = f ( x ) x i
where S(xi) denotes the sensitivity coefficients, f(x) is the model output and xi (i = 1, 2, 3, …, n) is the ith input parameter.
In the experiment, the tensile strength model was fitted by a battery of discrete experimental data. To obtain a given single factor model, the other two factors were simultaneously set as the zero level. The sub models of heating temperature, tape tension and roller pressure are as follows,
y T = f ( x 1 , x ¯ 2 , x ¯ 3 ) = 1121.008 + 10.31 x 1 71.604 x 1 2 y F = f ( x ¯ 1 , x 2 , x ¯ 3 ) = 1121.008 + 12.488 x 2 182.739 x 2 2 y P = f ( x ¯ 1 , x ¯ 2 , x 3 ) = 1121.008 + 8.758 x 3 101.6889 x 3 2
where the yT, yF and yP denote the sub models of heating temperature, tape tension and roller pressure, respectively; and x ¯ 1 , x ¯ 2 and x ¯ 3 represent the zero level of each variable.
According to Equation (5), the sensitivity models of three winding parameter including heating temperature, tape tension and roller pressure can be described as follows:
y T = f ( x 1 , x ¯ 2 , x ¯ 3 ) / x 1 = 10.31 143.208 x 1 y F = f ( x ¯ 1 , x 2 , x ¯ 3 ) / x 2 = 12.488 365.478 x 2 y P = f ( x ¯ 1 , x ¯ 2 , x 3 ) / x 3 = 8.758 203.3778 x 3
Figure 4 shows the effects of processing parameters on the tensile strength model. As shown in Figure 4a, the regression curves of sub-model are parabola going downwards. The results show that the tensile strength increases first and then decreases with the increase of each process parameter. Meanwhile, x2 has significant effect on the tensile strength of tape winding products. The value of tensile strength will change greatly as the value of x2 changes. The tensile strength, relatively speaking, are not so sensitive to the x1 and x3. According to Figure 4b, the images of sensitivity model of each factor are three linear functions with different slopes. The results show that, for the composite tape winding process, the tensile strength is most sensitive to the variation of tape tension, next sensitive to roller pressure and least sensitive to heating temperature. All the process parameters have a positive impact first and then a negative one on the tensile strength. In addition, the impact degree is gradually weakened with the increasing of the variable first and then gradually strengthened.

4.2. Iso-Surfaces Analysis for Parameter Range

The previous sensitivity analysis shows that the tensile strength of composite winding products is deeply influenced by the processing parameters. In three-dimensional space, a set of curved surface graphs are fitted by points with the same magnitude in a certain physical quantity. Then the curved surface can be applied to describe the distribution law of those physical quantities with continuous distribution characteristics [28].
According to the experiment’s results, the range of tensile strength value is 721.39~1124.31 MPa. For a more intuitive analysis of process parameters, four representative tensile strength values 800, 900, 1000 and 1100 MPa were selected as the certain objective value to complete drawing of three-dimensional iso-surfaces. The iso-surfaces analysis of different tensile strength values in Equation (4) were carried out on MATLAB (R2017a, The MathWorks Company, Natick, MA, USA) software platform. The results of iso-surfaces analysis are displayed in the Figure 5a–d. For the iso-surfaces of Ts = 800, 900 and 1000 MPa, x1 has the same range of [−1, 1] as shown in Figure 5a–c. Only in the case of Ts = 1100 MPa in Figure 5d, x1 falls into the [−0.49, 0.64]. When the tensile strength value is 800 MPa, x2 has the compound interval [−1, −0.54] and [0.72, 1] in Figure 5a. Similarly, x2 has a compound interval of [−1, −0.02] and [0.33,1] when Ts = 900 MPa in Figure 5b. As shown in Figure 5c,d, x2 falls into the [−0.85, 0.92] and [−0.34, 0.41], respectively. For the iso-surfaces of Ts = 900 and 1000 MPa, x3 has the same interval [−1, 1] in Figure 5b,c. However, in Figure 5c, x3 falls into the compound interval [−1, −0.37] and [0.45, 1] when Ts = 800 MPa. In addition, in Figure 5d, x3 falls into the range of [−0.49, 0.52]. The iso-surfaces analysis indicate that in the composite tape winding process, for a particular value of Ts, the parameter interval of x2 is the narrowest, followed by x3 and x1. Meanwhile, the results indicate that x2 has the largest impact on the value of Ts. Therefore, to get a higher tensile strength value, the parameter interval and control strategy of x2 need to be seriously considered.

4.3. Coupling Analysis for Process Parameters

In the actual winding process of composite materials prepreg tape, individual process parameters have a crucial impact on the performance of the winding product. In addition, the coupling between various process parameters will have a non-negligible interactive effect on the quality of the product such as tensile strength. The influence mechanism of each process parameter on tensile strength is affected by other parameters value. The process parameters including heating temperature, tape tension and winding pressure are coupled with each other to determine the final tensile strength of winding product. According to the response surface analysis results, the polynomial regression model of the tensile strength in composite tape winding process can be obtained as:
T S = f ( T , F , P ) = 290.644 + 11.903 T + 3.356 F + 0.741 P + 3.592 × 10 3 T F 1.199 × 10 3 T P 6.435 × 10 4 F P 7.956 × 10 2 T 2 4.568 × 10 3 F 2 1.808 × 10 4 P 2
where TS represents the tensile strength of winding product, MPa; T denotes the heating temperature, °C; F is the tape tension, N; and P is the roller pressure, N.
The coupling effect of process parameters on the tensile strength of winding products were analyzed using the response surface design software Design Expert 10 (Stat-Ease, Inc., Minneapolis, MN, USA). And Figure 6 shows the response surfaces and contour lines for the effect of parameters coupling on tensile strength. Figure 6a demonstrates the coupling effect of heating temperature and tape tension on the tensile strength of winding product. During the winding molding process, prepreg tapes can be heated to enhance the viscous fluidity of the resin. The heated resin is conducive to improving the degree of interface bonding between the resin and the fiber and between the tape and the laminate. When the temperature is low, the resin has poor fluidity. At this time, a small tape tension cannot provide a suitable deformation for the prepreg tape. However, excessive tension can easily cause the breaking of fibers in the tape, resulting in uneven stress distribution. When the temperature is too high, the resin will undergo a local curing reaction, which will reduce the overall mechanical properties of the winding product. The highest value of tensile strength locates in the area of 55~85 °C for heating temperature and 250~350 N for tape tension.
Figure 6b displays the coupling effect of heating temperature and roller pressure on the tensile strength of winding product. The positive pressure produced by the compaction roller can be employed to squeeze out the bubbles between the tape and laminated layers. The resin on the contact surface of the prepreg tape and the laminate will fuse into one under the suitable temperature and positive pressure. When the temperature is low, the polymer matrix does not reach the molten state. At this moment, a large pressure is difficult to improve the surface bonding strength for the poor degree of molecular chain diffusion. With the increase of temperature, the tensile strengths of the winding products had improved significantly. However, an excessive high temperature will cause the curing reaction of the polymer matrix, which will seriously reduce the products’ mechanical properties. When the temperature is suitable, too small pressure cannot make the tape and the laminate come into intimate contact. Conversely, extreme high positive pressure will excessively extrude the resin in a molten state from the laminates in large amounts. Looking at the whole picture, the tensile strength has a generally high value. As shown in Figure 6b, the ideal tensile strength appears in the area of 55~85 °C for heating temperature and 1000~1500 N for roller pressure.
Figure 6c shows the coupling effect of tape tension and roller pressure on the tensile strength of winding product. As a whole, the span of the tensile strength range is relatively large. The lager tensile strength values appear in the area of 250~350 N for tape tension and 900~1600 N for roller pressure. When both the tape tension and roller pressure are very low or high, the value of tensile strength is the smallest. Meanwhile, the value of tensile strength will first increase rapidly and then decrease greatly with the increase of the tape tension. When the tape tension is small, a large number of gaps appear between the tape and the base laminate. If the roller pressure is also small, at this time, these bubbles will be difficult to squeeze out. The degree of bonding between the winding layers becomes lower, resulting in a significant decrease for tensile strength. When the tape is wound around the mandrel through the hot pressing roller, the tension will generate a part of the normal pressure. If the tension and pressure are both large, the tape will be severely deformed, which will severely affect the mechanical properties of the winding product.

4.4. Optimization and Validation

To obtain a greater tensile strength value, the excellent parameter combination needs to be provided by optimizing process based on the desirability function. In this paper, the response was optimized using the general linear scale desirability function. The desirability function was beneficial for searching for the most favorable point in the design space [34,42]. Here, the response of tensile strength was the larger the better. Let ŷ represents the target experimental response. When the ŷ exceeds a particular criteria value, the desirability value equals to 1. If the ŷ is less than a particular criteria value, the desirability value equals to 0. Then the desirability function for the response to be maximized can be described as follow [35]:
d i = { 0 , y ^ y m i n ( y ^ y m i n y m a x y m i n ) r , y m i n y ^ y m a x 1 , y ^ y m i n
where di is the individual desirability of the response to be maximized; ymin and ymax denote the lowest and highest limits for each response, respectively; and r represents the weight, which can be specified by the user, r ≥ 0.
A global desirability function, D, also be called as overall desirability, is defined as the geometric mean of individual desirability values in Equation (10) [30]:
D = ( i = 1 n d i ) 1 / n
where the n is the number of the experimental responses.
In this article, the response surface design software Design Expert 10 (Stat-Ease, Inc., Minneapolis, MN, USA) was applied to determine the maximum desirability. The results were sorted by the desirability value from largest to smallest. Then the combination of process condition with a highest desirability value D will be selected as the optimum condition for the response.
According to the analysis results of the response surface software, the optimal parameter combination was selected as heating temperature with 72 °C, tape tension with 307 N and roller pressure with 1263 N. The maximum tensile strength value predicted by the response surface model was 1121.68 MPa. Figure 7 shows the desirability values for the response value. As shown in Figure 7, the highlighted dot from each ramp represents the recommended value for corresponding process parameter. Meanwhile, the desirability value was 0.993, which means that the target value and the response value have a high degree of closeness. To verify the effectiveness and reliability of the optimization results, three groups of repeated winding experiments were carried out using the same developed parameter combination. The materials, environment, testing process and so forth, in the verification experiments were the same as the previous experiments. Table 4 shows the results of validation experiments using the optimal parameter combination. The average tensile strength value of the testing specimens, 1131.26 MPa, was 9.58 MPa higher than the predicted value. Meanwhile, the average relative error of the verification tests was 0.85%. The comparison results showed that the quadratic regression model provided by response surface methodology was useful and reliable for the prepreg tape winding process. In addition, the optimal parameter combination has a positive guiding significance for improving the quality of winding products.

5. Conclusions

This article was dedicated to researching the influence mechanism and optimization analysis of parameters for the composite prepreg tape winding process. Process parameters including heating temperature, tape tension and roller pressure have an important impact on the winding products’ mechanical property such as tensile strength.
According to sensitivity analysis for single parameter, the tensile strength is most sensitive to the variation of tape tension, next sensitive to roller pressure and least sensitive to heating temperature. The tape tension has the most significant effect on the tensile strength of winding products. Meanwhile, the parameter interval of the processing parameters can be effectively determined via iso-surfaces for different tensile strength values in three-dimensional factor space. For a particular value of tensile strength, the parameter interval of tape tension is the narrowest according to the iso-surface diagrams of tensile strength. When only considering the enhancement of the tensile strength, the parameter interval and control strategy of tape tension needs to be confirmed first. In addition, three key process parameters, including heating temperature, tape tension and roller pressure, are coupled with each other to determine the tensile strength of the winding products. Then the outstanding parameter combination with the heating temperature of 72 °C, tape tension of 307 N and roller pressure of 1263 N was provided by response surface design software via desirability function method. The validation experiments showed that the optimal parameter combination has a positive guiding significance for improving the quality of the winding product.

Author Contributions

Methodology and writing—original draft preparation, B.D.; writing—review and editing, Y.S.; project administration, T.Y.; funding acquisition, P.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China with Grant number 51475377.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lukaszewicz, D.H.J.A.; Ward, C.; Potter, K.D. The engineering aspects of automated prepreg layup: History, present and future. Compos. Part B Eng. 2012, 43, 997–1009. [Google Scholar] [CrossRef]
  2. Oromiehie, E.; Prusty, B.G.; Compston, P.; Rajan, G. Automated fibre placement based composite structures: Review on the defects, impacts and inspections techniques. Compos. Struct. 2019, 224, 110987. [Google Scholar] [CrossRef]
  3. Rafiee, R. On the mechanical performance of glass-fibre-reinforced thermosetting-resin pipes: A review. Compos. Struct. 2016, 143, 151–164. [Google Scholar] [CrossRef]
  4. Zu, L.; Xu, H.; Zhang, B.; Li, D.; Zi, B.; Zhang, B. Design and production of filament-wound composite square tubes. Compos. Struct. 2018, 191, 202–208. [Google Scholar] [CrossRef]
  5. Saenz-Castillo, D.; Martin, M.I.; Calvo, S.; Rodriguez-Lence, F.; Guemes, A. Effect of processing parameters and void content on mechanical properties and ndi of thermoplastic composites. Compos. Part A Appl. Sci. Manuf. 2019, 121, 308–320. [Google Scholar] [CrossRef]
  6. Leon, A.; Perez, M.; Barasinski, A.; Abisset-Chavanne, E.; Defoort, B.; Chinesta, F. Multi-scale modeling and simulation of thermoplastic automated tape placement: Effects of metallic particles reinforcement on part consolidation. Nanomaterials 2019, 9, 695. [Google Scholar] [CrossRef] [Green Version]
  7. Henriquez, R.G.; Mertiny, P. 3.21 Filament winding applications. In Comprehensive Composite Materials II; Elsevier: Amsterdam, Holland, 2018; Volume 3, pp. 556–577. [Google Scholar]
  8. Han, Z.Y.; Sun, S.Z.; Li, W.Q.; Zhao, Y.X.; Shao, Z.X. Experimental study of the effect of internal defects on stress waves during automated fiber placement. Polymers 2018, 10, 413. [Google Scholar] [CrossRef] [Green Version]
  9. Sun, S.; Han, Z.; Fu, H.; Jin, H.; Dhupia, J.S.; Wang, Y. Defect characteristics and online detection techniques during manufacturing of frps using automated fiber placement: A review. Polymers 2020, 12, 1337. [Google Scholar] [CrossRef]
  10. Dalibor, I.H.; Lisboa, T.V.; Marczak, R.J.; Amico, S.C. Optimum slippage dependent, non-geodesic fiber path determination for a filament wound composite nozzle. Eur. J. Mech. A Solids 2020, 82, 103994. [Google Scholar] [CrossRef]
  11. Stokes-Griffin, C.M.; Compston, P.; Matuszyk, T.I.; Cardew-Hall, M.J. Thermal modelling of the laser-assisted thermoplastic tape placement process. J. Thermoplast. Compos. Mater. 2015, 28, 1445–1462. [Google Scholar] [CrossRef]
  12. Mazumdar, S.K.; Hoa, S. Determination of manufacturing conditions for hot-gas-aided thermoplastic tape winding. J. Thermoplast. Compos. Mater. 1996, 9, 35–53. [Google Scholar] [CrossRef]
  13. Mack, J.; Schledjewski, R. Filament winding process in thermoplastics. In Manufacturing Techniques for Polymer Matrix Composites (PMCS); Elsevier: Amsterdam, Holland, 2012; pp. 182–208. [Google Scholar]
  14. Nath, N.K. Optimization of tape winding process parameters to enhance the performance of solid rocket nozzle throat back up liners using taguchi’s robust design methodology. J. Inst. Eng. Ser. C 2017, 98, 479–484. [Google Scholar] [CrossRef]
  15. Dai, S.C.; Ye, L. Characteristics of cf/pei tape winding process with on-line consolidation. Compos. Part A Appl. Sci. Manuf. 2002, 33, 1227–1238. [Google Scholar] [CrossRef]
  16. Zaami, A.; Schäkel, M.; Baran, I.; Bor, T.C.; Janssen, H.; Akkerman, R. Temperature variation during continuous laser-assisted adjacent hoop winding of type-iv pressure vessels: An experimental analysis. J. Compos. Mater. 2020, 54, 1717–1739. [Google Scholar] [CrossRef] [Green Version]
  17. Shi, Y.-Y.; Tang, H.; Yu, Q. Development of a NC tape winding machine. In Global Design to Gain a Competitive Edge; Springer: London, UK, 2008; pp. 753–762. [Google Scholar]
  18. Sorrentino, L.; Anamateros, E.; Bellini, C.; Carrino, L.; Corcione, G.; Leone, A.; Paris, G. Robotic filament winding: An innovative technology to manufacture complex shape structural parts. Compos. Struct. 2019, 220, 699–707. [Google Scholar] [CrossRef]
  19. Libonati, F.; Vellwock, A.E.; El Louizi, F.; Hoffmann, R.; Colombo, C.; Ziegmann, G.; Vergani, L. Squeeze-winding: A new manufacturing route for biomimetic fiber-reinforced structures. Compos. Part A Appl. Sci. Manuf. 2020, 132, 105839. [Google Scholar] [CrossRef]
  20. Wang, Q.; Li, T.; Wang, B.; Liu, C.; Huang, Q.; Ren, M. Prediction of void growth and fiber volume fraction based on filament winding process mechanics. Compos. Struct. 2020, 246, 112432. [Google Scholar] [CrossRef]
  21. Eggers, F.; Almeida, J.H.S.; Azevedo, C.B.; Amico, S.C. Mechanical response of filament wound composite rings under tension and compression. Polym. Test. 2019, 78, 105951. [Google Scholar] [CrossRef]
  22. Xu, X.-M.; Zhang, W.-X.; Ding, X.-L.; Zhang, M.; Wei, S.-H. Design and analysis of a novel tension control method for winding machine. Chin. J. Mech. Eng. 2018, 31, 101. [Google Scholar] [CrossRef] [Green Version]
  23. He, X.D.; Shi, Y.Y.; Kang, C.; Yu, T. Analysis and control of the compaction force in the composite prepreg tape winding process for rocket motor nozzles. Chin. J. Aeronaut. 2017, 30, 836–845. [Google Scholar] [CrossRef] [Green Version]
  24. Colombo, C.; Vergani, L. Optimization of filament winding parameters for the design of a composite pipe. Compos. Part B Eng. 2018, 148, 207–216. [Google Scholar] [CrossRef]
  25. Nikishova, A.; Comi, G.E.; Hoekstra, A.G. Sensitivity analysis based dimension reduction of multiscale models. Math. Comput. Simul. 2020, 170, 205–220. [Google Scholar] [CrossRef] [Green Version]
  26. Saltelli, A.; Aleksankina, K.; Becker, W.; Fennell, P.; Ferretti, F.; Holst, N.; Li, S.S.; Wu, Q.L. Why so many published sensitivity analyses are false: A systematic review of sensitivity analysis practices. Environ. Model. Softw. 2019, 114, 29–39. [Google Scholar] [CrossRef]
  27. Song, X.; Zhang, J.; Zhan, C.; Xuan, Y.; Ye, M.; Xu, C. Global sensitivity analysis in hydrological modeling: Review of concepts, methods, theoretical framework, and applications. J. Hydrol. 2015, 523, 739–757. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, A.M.; Han, Y.; Pan, Y.H.; Xing, H.W.; Li, J.Z. Optimum surface roughness prediction for titanium alloy by adopting response surface methodology. Results Phys. 2017, 7, 1046–1050. [Google Scholar] [CrossRef]
  29. Nagaraj, S.; Natarajan, V. Relation-aware isosurface extraction in multifield data. IEEE Trans. Vis. Comput. Graph. 2010, 17, 182–191. [Google Scholar] [CrossRef]
  30. Hernandez-Vazquez, J.-M.; Garitaonandia, I.; Fernandes, M.H.; Muñoa, J.; Lacalle, L.N.L.d. A consistent procedure using response surface methodology to identify stiffness properties of connections in machine tools. Materials 2018, 11, 1220. [Google Scholar] [CrossRef] [Green Version]
  31. Dharma, S.; Masjuki, H.H.; Ong, H.C.; Sebayang, A.H.; Silitonga, A.S.; Kusumo, F.; Mahlia, T.M.I. Optimization of biodiesel production process for mixed jatropha curcas-ceiba pentandra biodiesel using response surface methodology. Energy Convers. Manag. 2016, 115, 178–190. [Google Scholar] [CrossRef]
  32. Danmaliki, G.I.; Saleh, T.A.; Shamsuddeen, A.A. Response surface methodology optimization of adsorptive desulfurization on nickel/activated carbon. Chem. Eng. J. 2017, 313, 993–1003. [Google Scholar] [CrossRef]
  33. Lee, D.H.; Jeong, I.J.; Kim, K.J. A desirability function method for optimizing mean and variability of multiple responses using a posterior preference articulation approach. Qual. Reliab. Eng. Int. 2018, 34, 360–376. [Google Scholar] [CrossRef]
  34. Cho, H.J.; Lee, D.W.; Marasini, N.; Poudel, B.K.; Kim, J.H.; Ramasamy, T.; Yoo, B.K.; Choi, H.G.; Yong, C.S.; Kim, J.O. Optimization of self-microemulsifying drug delivery system for telmisartan using box—Behnken design and desirability function. J. Pharm. Pharmacol. 2013, 65, 1440–1450. [Google Scholar] [CrossRef] [PubMed]
  35. Jenarthanan, M.; Jeyapaul, R. Optimisation of machining parameters on milling of GFRP composites by desirability function analysis using taguchi method. Int. J. Eng. Sci. Technol. 2013, 5, 22–36. [Google Scholar] [CrossRef] [Green Version]
  36. Deng, B.; Shi, Y. Modeling and optimizing the composite prepreg tape winding process based on grey relational analysis coupled with bp neural network and bat algorithm. Nanoscale Res. Lett. 2019, 14, 296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Deng, B.; Shi, Y. Modeling and simulation of voids in composite tape winding process based on domain superposition technique. Appl. Compos. Mater. 2017, 25, 1219–1236. [Google Scholar] [CrossRef]
  38. Aslan, N.; Cebeci, Y. Application of box—Behnken design and response surface methodology for modeling of some turkish coals. Fuel 2007, 86, 90–97. [Google Scholar] [CrossRef]
  39. Ye, G.; Ma, L.; Li, L.; Liu, J.; Yuan, S.; Huang, G. Application of box—Behnken design and response surface methodology for modeling and optimization of batch flotation of coal. Int. J. Coal Prep. Util. 2020, 40, 131–145. [Google Scholar] [CrossRef]
  40. Standardization Administration of the People’s Republic of China. Test Method for Mechanical Properties of Ring of Filament-Winding Reinforced Plastics; Standardization Administration of the People’s Republic of China: Beijing, China, 2008; GT/B 1458-2008.
  41. Dezan, D.J.; Rocha, A.D.; Ferreira, W.G. Parametric sensitivity analysis and optimisation of a solar air heater with multiple rows of longitudinal vortex generators. Appl. Energy 2020, 263, 114556. [Google Scholar] [CrossRef]
  42. Derringer, G.; Suich, R. Simultaneous optimization of several response variables. J. Qual. Technol. 1980, 12, 214–219. [Google Scholar] [CrossRef]
Figure 1. Schematic of composites prepreg tape winding process.
Figure 1. Schematic of composites prepreg tape winding process.
Polymers 12 01843 g001
Figure 2. Experimental test samples machining and tensile strength testing process. (a) Turning machining process; (b) Standard testing rings; (c) Split disc tensile tests on electronic universal testing machine DDL100; (df) Typical failure mode of the ring specimen.
Figure 2. Experimental test samples machining and tensile strength testing process. (a) Turning machining process; (b) Standard testing rings; (c) Split disc tensile tests on electronic universal testing machine DDL100; (df) Typical failure mode of the ring specimen.
Polymers 12 01843 g002
Figure 3. Residual analysis results. (a) Residuals normal plot of tensile strength model; (b) Residual versus run plot of tensile strength model; (c) Predicted versus actual plot of tensile strength model.
Figure 3. Residual analysis results. (a) Residuals normal plot of tensile strength model; (b) Residual versus run plot of tensile strength model; (c) Predicted versus actual plot of tensile strength model.
Polymers 12 01843 g003
Figure 4. Effects of processing parameters. (a) The regression curve of sub-model; (b) The sensitivity model of each factor.
Figure 4. Effects of processing parameters. (a) The regression curve of sub-model; (b) The sensitivity model of each factor.
Polymers 12 01843 g004
Figure 5. Iso-surfaces for different tensile strength values in three-dimensional factor space. (a) Iso-surface diagram of TS = 800 MPa; (b) Iso-surface diagram of TS = 900 MPa; (c) Iso-surface diagram of TS = 1000 MPa; (d) Iso-surface diagram of TS = 1100 MPa.
Figure 5. Iso-surfaces for different tensile strength values in three-dimensional factor space. (a) Iso-surface diagram of TS = 800 MPa; (b) Iso-surface diagram of TS = 900 MPa; (c) Iso-surface diagram of TS = 1000 MPa; (d) Iso-surface diagram of TS = 1100 MPa.
Polymers 12 01843 g005
Figure 6. The response surface and contour lines map of process parameters on the tensile strength. (a) Coupling effect of temperature and tension; (b) Coupling effect of temperature and pressure; (c) Coupling effect of tension and pressure.
Figure 6. The response surface and contour lines map of process parameters on the tensile strength. (a) Coupling effect of temperature and tension; (b) Coupling effect of temperature and pressure; (c) Coupling effect of tension and pressure.
Polymers 12 01843 g006
Figure 7. Ramp function graph of desirability and the optimized result for tensile strength.
Figure 7. Ramp function graph of desirability and the optimized result for tensile strength.
Polymers 12 01843 g007
Table 1. Level of process parameters.
Table 1. Level of process parameters.
Levels
Experimental ParametersSymbolCodeUnitsLevel −1Level 0Level 1
Heating temperatureTx1°C4070100
Tape tensionFx2N100300500
Roller pressurePx3N50012502000
Table 2. The experiment results.
Table 2. The experiment results.
No.T/(°C)F/(N)P/(N)x1x2x3TS/(MPa)
11005001250110911.33
27030012500001124.31
310030050010−1976.92
47030012500001119.68
540300500−10−1899.4
6401001250−1−10865.11
7705002000011758.72
87030012500001123.53
9403002000−101972.47
10701005000−1−1721.39
117030012500001120.32
12405001250−110850.54
137010020000−11930.36
1410010012501−10839.68
151003002000101942.07
167030012500001117.2
177050050001−1935.85
Table 3. ANOVA results for response surface quadratic model of tensile strength.
Table 3. ANOVA results for response surface quadratic model of tensile strength.
SourceSSDFMSF ValueProb > F
Model2.7 × 105930,002.312567.59<0.0001
T850.371850.3772.77<0.0001
F1247.511247.5106.76<0.0001
P613.551613.5552.510.0002
TF1858.4711858.47159.05<0.0001
TP2911.6812911.68249.18<0.0001
FP37,268.3137,268.33189.41<0.0001
T221,587.93121,587.931847.49<0.0001
F21.406 × 10511.41 × 10512,032.87<0.0001
P243,539.59143,539.593726.1<0.0001
Residual81.8711.69
Lack of Fit47.79315.931.870.2749
Pure Error3448.5
Cor Total2.701 × 10516
DF: Degrees of Freedom; SS: Sum of Squares; MS: Mean Square
Table 4. Results of validation experiments using the optimal parameter combination.
Table 4. Results of validation experiments using the optimal parameter combination.
No.T/(°C)F/(N)P/(N)Tensile Strength/(MPa)Relative Error
PredictedExperiment
17230712631121.681129.160.67%
21133.831.08%
31130.790.81%

Share and Cite

MDPI and ACS Style

Deng, B.; Shi, Y.; Yu, T.; Zhao, P. Influence Mechanism and Optimization Analysis of Technological Parameters for the Composite Prepreg Tape Winding Process. Polymers 2020, 12, 1843. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12081843

AMA Style

Deng B, Shi Y, Yu T, Zhao P. Influence Mechanism and Optimization Analysis of Technological Parameters for the Composite Prepreg Tape Winding Process. Polymers. 2020; 12(8):1843. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12081843

Chicago/Turabian Style

Deng, Bo, Yaoyao Shi, Tao Yu, and Pan Zhao. 2020. "Influence Mechanism and Optimization Analysis of Technological Parameters for the Composite Prepreg Tape Winding Process" Polymers 12, no. 8: 1843. https://0-doi-org.brum.beds.ac.uk/10.3390/polym12081843

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop