Hole diameter effects on fiber-reinforced thermoplastic composite discs possessing distinct fiber arrays under thermal loading

Abstract

This study explores the thermomechanical behavior of steel fiber-reinforced high-density polyethylene matrix thermoplastic composite (TPC) discs possessing distinct fiber arrays and the influence of circular hole diameter, upon exposure to convective air cooling loading. The discs were manufactured and their thermomechanical properties were assessed. Models of the discs for different D/W ratios (where D is the circular hole diameter and W is the disc width) were constructed and analyzed. Two parameters, the circular-hole diameter and the fiber array configuration, were varied. Thermal stresses generated during the cooling of each disc were then determined. Residual thermal stress and equivalent plastic strain for each respective D/W ratio were presented in graphical form, enabling evaluation and comparison. A parametric study was conducted to identify the influence of both circular hole diameter and distinct fiber reinforcement upon residual thermal stress and plastic deformation for steel fiber/TPC discs subjected to thermal loading.

Keywords

Thermoplastic composite disc thermal residual stress plastic strain fiber arrangement finite element analysis

Introduction

Using composite technology, it is possible to design new, multifunctional materials with unique properties, unobtainable using any monolithic material. Thermal-induced internal stresses and their corresponding residual thermal stresses, arising during the production of thermoplastic composites (TPCs), have an important potential impact on the performance and properties of the composites. Fiber-reinforced TPCs have found widespread applications in a range of engineering fields, such as aircraft structure design¹ on account of their robustness, low sensitivity to moisture effects, and adaptability to mass production and repair. For the design and analytical modeling of such composites, accurate knowledge of residual stresses under loading is required since it can cause failures, including fiber misalignment,² transverse cracking,³ delamination,⁴ warpage,⁵ and fatigue behavior.⁶ Residual thermal stress is one form of residual stress. It arises due to greater shrinkage of the matrix compared with that of the fibers, following the processing of TPC laminates, as they cool from a high processing temperature to service temperature (22°C). Several researchers have studied residual thermal stresses and their causes under different conditions. Nairn and Zoller⁷ published an experimental study of thermal residual stress in TPC. According to Favre,⁸ if long-term and environmental parameters are excluded, the magnitude of thermal residual stresses in composite structures is dependent on four parameters: the temperature differential, the cooling thermal expansion/shrinkage coefficients of the composite constituents (or plies), the constituent elastic coefficients, and the fiber to volume ratio. These parameters are also dependent on other factors such as the thermoplastic matrix morphology (semicrystalline or amorphous), fiber types, fiber–matrix interface properties, fiber morphology (woven or unidirectional prepreg), and the processing conditions. Parlevliet et al.⁹ reviewed these parameters separately, starting at the micromechanical level and the fiber-dominated properties, followed by the matrix-dominated properties. Barnes and Byerly¹⁰ found that the further the temperature drops from the glass transition temperature, the more thermal residual stresses exist within all polymer composites. The study has shown that for thermoplastics matrix composites, the (radial) residual stress increases with an increasing cooling rate.¹¹ Sayman and Arman investigated the thermal residual stresses occurring in TPC discs reinforced by steel fibers with circular arrays, under steady state,¹² uniform,¹³ linear,¹⁴ and parabolic¹⁵ temperature distribution conditions by numerical solution and also finite element analyses. Hu and Weng¹⁶ analyzed the influence of thermal residual stresses on the deformation behavior of a composite using a new micromechanical method. Huang¹⁷ presented a formula for determining the strength of unidirectional composites with thermal residual stresses. Fiber arrangement also plays an important part in the formation of thermal residual stresses occurring in composite production processes operating under normal working conditions. Hobbiebrunken et al.¹⁸ examined the effects of nonuniform fiber arrangements on thermal residual stresses in carbon fiber-/epoxy-laminated composites. The effects of different fiber arrangements, including regular fiber arrays (square and hexagonal arrays) and a random fiber array, on residual thermal stresses in unidirectional composites of various fiber volume ratios were investigated by Jin et al.¹⁹ Tsai and Chi²⁰ studied the thermal residual stress effect on the constitutive behaviors of fiber composites with three different fiber arrangements, that is square edge packing, square diagonal packing, and hexagonal packing. Moreover, geometric stress concentration was found to arise due to holes, cutouts, or abrupt change in the section and to limit the load-carrying capabilities of structures under loading. As a result, stress distribution is not uniform throughout the whole cross section. Failure, such as cracking and plastic deformation, frequently occurs at points of geometric stress concentration, since stresses reach their maximum magnitude at those points. In the literature,^21
–23 various descriptions can be found of stress concentration in plates with a hole. Toubal et al.²⁴ have experimentally studied the stress concentration in composite plates with a circular hole. Ting et al.²⁵ presented a stress analysis theory, a method using alternating rhombic arrays for discs with multiple circular holes. Despite potentially wide-ranging applications, few studies have reported the influence of the hole diameter on thermal residual stresses within metal fiber TPC discs. The purpose of this research study is to investigate the influence of circular hole diameter on residual thermal stresses and plastic deformation, for steel fiber-reinforced TPC discs with a distinct fiber array exposed to convective air cooling condition.

Materials and methods

First, three TPC discs, each with different fiber array types (Figure 1) were manufactured. To produce each disc, raw granules of high-density polyethylene were initially put into molds and held for 5 min at 2.5 MPa. Electrical resistance was used to hold the temperature at 173°C. The steel wires were wound around the disc in patterns to form circular and radial fiber arrays. Having produced the matrix layer and its reinforcement, two matrix layers and one fiber array were bonded together under the same pressure and temperature conditions to obtain a TPC disc. Each disc had a diameter of 160 mm and a thickness of 2 mm, and the fiber volume ratio was measured at 8%.

Figure 1.

Woven, circular, and radial steel fiber arrays for TPC discs, respectively. TPC: thermoplastic composite.

An analysis was performed by the authors and the material and thermal properties of the discs were determined as displayed in Table 1. E, G, and $ν$ are the modulus of elasticity, modulus of rigidity, and Poisson’s ratio, respectively.

Table 1.

Average thermomechanical properties of each TPC disc.

Fiber array	E ₁ (MPa)	E ₂ (MPa)	$ν_{12}$	X (MPa)	Y (MPa)	$α_{r}$ 1 × 10⁻⁶	$α_{θ}$ 1 × 10⁻⁶	G (MPa)
Woven	16,058	16,058	0.378	20.45	20.45	12.185	12.185	245
	E _r (MPa)	E_θ (MPa)	$ν_{r θ}$	X (MPa)	Y (MPa)	$α_{r}$ 1 × 10⁻⁶	$α_{θ}$ 1 × 10⁻⁶	G _rθ (MPa)
Circular	2456	34,424	0.360	6	9	119	11.8	722
Radial	31,885	2763	0.296	3	2.5	34.2	4.77	847

TPC: thermoplastic composite.

Finite element method (FEM) and Tsai–Hill yield criterion were used to provide a numerical solution. It was important to ensure that the cooling process for each composite disc was consistent and that accurate data were recorded. ABAQUS engineering simulation software was used for this purpose. The numerical and experimental cooling curves of the discs were measured. Then the numerical and experimental curves were converged so that a precise value for the convective film coefficient could be calculated. Numerical cooling curves were obtained using data taken from the results of analysis solutions of each disc in ABAQUS simulation software. In order to obtain experimental cooling curves, each composite disc was warmed to operational temperature (94°C) in an oven and then cooled to service temperature (22°C). During this process, surface temperatures of the discs were measured at the same surface position at 20-s intervals by a digital thermometer as the discs cooled from 94°C to 22°C. Finally, the numerical and experimental cooling curves of each disc were plotted (Figure 2).

Figure 2.

Numerical and experimental cooling curves of TPC discs. TPC: thermoplastic composite.

The value for the convective film coefficient affecting the outer disc surfaces was determined as 6.67 W/m² °C h. Other factors influencing the convective film coefficient (h) of the composite disc on which analysis was carried out were its density (d), specific heat (C _p), and thermal conductivity (k) as shown in Table 2.

Table 2.

Physical properties of the TPC discs.

Fiber array	d (ton/mm³)	C _p (mJ/ton K)	k (mW/mm K)	h (mW/mm² K)
Woven	1317 × 10⁻¹²	2250 × 10⁶	0.5	6.67 × 10⁻³
Circular	1184 × 10⁻¹²	2250 × 10⁶	0.5	6.67 × 10⁻³
Radial	941 × 10⁻¹²	2250 × 10⁶	0.5	6.67 × 10⁻³

TPC: thermoplastic composite.

Finite element analysis

FEM was used as a basis for a numerical solution. Cartesian coordinates for woven fiber array and polar coordinates for circular and radial fiber array were employed in the analyses. Suitable element selection was the primary task in the finite element analysis. Multiple checks and convergence tests were made to guide the selection of suitable elements, including the determination of optimum element length. An elastoplastic material model, demonstrating anisotropic yield properties, was populated with data obtained in characterization tests. The models were established using quadrilateral plane solid elements, specified as CPS8T in the ABAQUS package. In this system, CPS8T is defined as “eight-node plane stress thermally coupled quadrilateral, biquadratic displacement, and bilinear temperature,” and it has an orthotropic plate element size of four. The symmetry of the disc geometry and loading conditions meant that assessing a quarter of each disc was sufficient. In order to construct graphical images and finite element modeling of the disc models for different D/W ratios, ABAQUS engineering simulation software was used. Figure 3 provides an example for graphical image and meshing the quarter disc model for a circular fiber array. This was one of the disc models used in this study, and in this case, it depicts a disc with a 0.3 D/W ratio. Mapped meshing is used for all models, so that more elements are employed near the surface–hole boundaries. After the mesh generation process was complete, all models were divided into elements and nodes.

Figure 3.

Graphical images (left) and finite element modeling (right) of the quarter TPC disc model possessing circular fiber array for D/W = 0.3 ratio. TPC: thermoplastic composite.

The discs were analyzed using simply supported boundary conditions, where both U_y and U_z are equal to zero. The Tsai–Hill theory was used to determine the yield, since the yield strengths in tension and compression were the same. According to this criterion, the equivalent stress in the fiber direction can be written as²⁶

σ_{e q} = \sqrt{σ_{θ}^{2} - σ_{r} σ_{θ} + \frac{σ_{r}^{2} X^{2}}{Y^{2}}}

where X, Y and $σ_{r}$ , $σ_{θ}$ are the yield strengths in the longitudinal and transverse directions for woven, circular, and radial fiber arrays, respectively.

For this study, the temperature at the beginning of the cooling process and the service temperature were 94°C and 22°C, respectively. Residual stresses and plastic strains that occurred were measured from the results of numerical solution of ABAQUS simulation software and then analyzed for the following ratios

\frac{R_{3}}{R_{2} - R_{1}} = \frac{D}{W} = \{[0], [0.1], [0.3], [0.5]\}

In case where D/W = 0 and $W = R_{2}$ , the disc does not possess a hole. The analysis in this case is extremely difficult; therefore, finite element model was adopted for whole analysis.

Results and discussion

The numerical results presented here are for fiber-reinforced thermoplastic orthotropic discs with woven, circular, and radial fiber arrays. Thermal stresses generated during convective air cooling were obtained, and residual stresses and plastic strains present in the discs at service temperature were evaluated for a number of D/W ratios.

According to Tsai–Hill yield criterion, equivalent stresses should be considered. Substituting the values of thermal stress components that occurred in the longitudinal and transverse directions during cooling into equation (1), equivalent stresses can be calculated for each disc employing software conducted by the authors. Equivalent thermal stresses revealed at small values compared to thermal stress components.

Figure 4 shows the change in equivalent thermal stresses present on the same point at the inner diameter in discs with woven, circular, and radial fiber arrays during cooling. The results were obtained subsequent to cooling the discs from the processing temperature (94°C) to the service temperature. Different coefficients of thermal expansion in radial and tangential directions for circular and radial discs caused higher thermal stresses than those experienced in woven discs for a uniform temperature distribution.

Figure 4.

Equivalent thermal stresses on the same point at inner diameter during convective air cooling for the TPC discs with distinct fiber array. TPC: thermoplastic composite.

On the other hand, thermal residual stresses attain their greatest values at the inner diameter edge in the discs, for all types of fiber array when the D/W = 0. As the D/W ratio is increased from 0 to 0.5, the numerical results indicate that thermal residual stresses occur at both the inner diameter edge and around each circular hole in the disc. For example, Figure 5 illustrates the thermal residual stress through the x direction distribution in the TPC disc with a circular fiber array at service temperature for each of the D/W ratios used in this study.

Figure 5.

Thermal residual stress distribution in the TPC disc with circular fiber array for each of D/W ratio. TPC: thermoplastic composite.

In order to correctly compare the thermomechanical behavior of the discs, it is necessary to introduce residual stress ratio (SR) for this study. The residual SR is defined as the maximum thermal residual stress in the circular or radial fiber array discs expressed in proportion to the same value for the disc with woven fiber array. To determine the highest thermal residual stress and equivalent plastic strain (not only strain in longitudinal and transverse directions) value occurring in each disc at service temperature, numerical results obtained from simulation software was employed. $σ_{x}_{m a x w o v e n}$ , $σ_{y}_{max w o v e n}$ , and $τ_{max w o v e n}$ are the highest values for $σ_{x}, σ_{y}$ , and $τ_{x y}$ residual stress values occurring in the woven disc. $σ_{r}_{max c i r c u l a r}$ , $σ_{θ}_{max c i r c u l a r}$ , and $τ_{max c i r c u l a r}$ are the highest values at service temperature for $σ_{r}$ , $σ_{θ}$ , and $τ_{r θ}$ residual stresses occurring in discs with circular arrays; and $σ_{r}_{max r a d i a l}$ , $σ_{θ}_{max r a d i a l}$ , and $τ_{max r a d i a l}$ indicate the maximum values for $σ_{r}$ , $σ_{θ}$ , and, $τ_{r θ}$ found in radial disc under these conditions.

$\frac{σ_{r}_{max C i r c u l a r}}{σ_{x}_{max W o v e n}}$ = SR1, $\frac{σ_{r}_{max R a d i a l}}{σ_{x}_{max W o v e n}}$ = SR2 $\frac{σ_{θ}_{max C i r c u l a r}}{σ_{y}_{max W o v e n}}$ =SR3, $\frac{σ_{θ}_{max R a d i a l}}{σ_{y}_{max W o v e n}}$ = SR4, and $\frac{τ_{}_{max C i r c u l a r}}{τ_{}_{max W o v e n}}$ = SR5, and $\frac{τ_{}_{max R a d i a l}}{τ_{max W o v e n}}$ = SR6 are the residual SRs used in this study, and they are SR1, SR2, SR3, SR4, SR5, and SR6, respectively.

Figures 6 to 11 show the effect of D/W on SR values and equivalent plastic strain occurring in the composite discs with woven, circular, and radial fiber array under uniformly distributed convective air cooling loading at service temperature. The following observations can be made from the results obtained. Figures 6 and 7 show the effect of D/W on SR1 and SR2 ratios and the $σ_{r}$ -induced plastic strain; Figure 6 for the circular fiber array and Figure 7 for the radial fiber array. Variation in SRs with respect to D/W has its greatest impact on the SR2 ratio and was significant for the SR1 ratio. The maximum value obtained for the SR2 ratio was 6, when D/W was equal to zero. Smaller stress occurred in the woven disc for all D/W ratios for both SR1 and SR2 ratios, when compared to the other types. In the case of D/W = 0, $σ_{r}_{max C i r c u l a r}$ was 1.6 times greater than $σ_{x}_{max W o v e n}$ and $σ_{r}_{m a x}_{R a d i a l}$ 6 times greater than $σ_{x}_{m a x}_{W o v e n}$ as seen in Figures 6 and 7, respectively.

Figure 11.

Influence of D/W ratio on plastic strain for disc with radial fiber array and SR6 ratio. SR: stress ratio.

Figure 7.

Influence of D/W ratio on plastic strain and SR2 ratio for discs with radial fiber array. SR: stress ratio.

Figure 6.

Influence of D/W ratio on plastic strain and SR1 ratio for discs with circular fiber array. SR: stress ratio.

The SR1 ratio continuously increased from 1.6 to 3.1 when D/W was increased from 0.1 to 0.5. That is, the value of $σ_{r max c i r c u l a r}$ became two times greater than $σ_{x}_{max w o v e n}$ as D/W was increased from 0.1 to 0.5. The SR2 ratio decreased from 6 to 2.9 when D/W was increased from 0 to 0.1 (surface holes were drilled). It then became constant at a value close to 3 when D/W was increased from 0.1 to 0.3. Then, the SR2 ratio increased to 5.5 when D/W increased from 0.3 to 0.5. No strain occurred in the woven disc, regardless of D/W ratio, but analysis found it present in circular and radial discs. Values of $σ_{r}$ −, $σ_{θ}$ −, and $τ_{r θ}$ −induced plastic strains were considered. The maximum strain value obtained was 7.4 × 10⁻² when D/W was equal to 0.5 for a circular array disc, as seen in Figure 6. In the case where D/W was equal to zero, $σ_{r}$ −induced strain value for a circular fiber array was 1 × 10⁻² as seen in Figure 6 and for a radial fiber array was 5 × 10⁻³ as seen in Figure 7. As D/W was increased from 0 to 0.1 through the drilling of surface holes, strain values remained constant for both the circular and radial discs. As D/W was increased from 0.1 to 0.5, the strain value continuously increased in the circular disc from 1 to 7. For the radial disc, it remained constant at 5 as D/W was increased from 0 to 0.1, then dropped to 3.6 as D/W reached 0.3, and then continuously increased from 3.6 to 11 as D/W was raised to 0.5.

Figures 8 and 9 show the effect of D/W on SR3, SR4 SRs, and $σ_{θ}$ -induced plastic strain; Figure 8 for a circular fiber array and Figure 9 for a radial fiber array. Variation in SR with respect to D/W was greatest for the SR4 ratio and was found to be significant for the SR3 ratio. The highest SR value obtained was 4.86 for SR4 when D/W was equal to 0.5. Minimum stress occurred in the woven disc for $σ_{y}$ stress for all D/W ratios. In the case where D/W was equal to zero, stress value for $σ_{y}$ in the woven disc was 1.5 times smaller than $σ_{θ}$ in the circular discs and 4.76 times smaller than $σ_{θ}$ in the radial discs. This can be seen in Figures 8 and 9.

Figure 9.

Influence of D/W ratio on plastic strain and SR4 ratio for discs with radial fiber array. SR: stress ratio.

Figure 8.

Influence of D/W ratio on plastic strain and SR3 ratio for discs with circular fiber array. SR: stress ratio.

The SR3 ratio continuously increased from 1.5 to 1.9 when D/W was increased from 0 to 0.1. It then remained constant at 2 when D/W was increased from 0.1 to 0.3, and then decreased to 1.25 when D/W was increased from 0.3 to 0.5. The SR4 ratio decreased from 4.76 to 4.72 when D/W was increased from 0 to 0.1 through the drilling of surface holes. It then sharply decreased from 4.72 to 4.55 when D/W was increased from 0.1 to 0.3 and continuously increased to 4.86 as D/W rose from 0.3 to 0.5. A maximum strain value for $σ_{y}$ was determined to be 7.4 × 10⁻² in the circular disc when the D/W ratio was equal to 0.5. Strain value for $σ_{y}$ with respect to D/W ratio was greatest in the circular disc. Maximum strain occurred for the circular disc when D/W was equal to 0.5. When D/W was equal to zero, strain value for $σ_{y}$ was 1 × 10⁻² in the circular disc and 4.25 × 10⁻³ in the radial disc. As D/W increased from 0 to 0.1 through the drilling of surface holes, strain value remained constant for both circular and radial discs. As D/W increased from 0.1 to 0.5, the strain value continuously increased in the circular disc from 1 × 10⁻² to 7.4 × 10⁻². As a comparison, in the radial disc, it remained constant at 4.25 × 10⁻³ as D/W was raised from 0 to 0.1, decreasing to 3.8 × 10⁻³ as D/W was increased to 0.3. It then continuously increased to 11 × 10⁻³ for as D/W reached a value of 0.5.

Figures 10 and 11 show the effect of D/W on SR5 and SR6 SRs for the highest $τ_{x y}, τ_{r θ}$ values and $τ_{r θ}$ -induced plastic strain; Figure 10 for circular fiber array and Figure 11 for radial fiber array. Variation in SR with respect to D/W was greatest for the SR6 ratio and significant for the SR5 ratio. The maximum SR value obtained was 7.85 for the SR6 ratio when D/W was equal to zero. Minimum stress for $τ_{x y}$ occurred in the disc with a woven fiber array for all D/W ratios. When D/W was equal to zero, the stress value for $τ_{x y}$ in the disc with woven fiber array was 3.4 times smaller than $τ_{r θ}$ in the disc with a circular fiber array and 7.85 times smaller than $τ_{r θ}$ in the disc with a radial fiber array. This can be seen in Figures 10 and 11.

Figure 10.

Influence of D/W ratio on plastic strain and SR5 ratio for discs with circular fiber array. SR: stress ratio.

The SR5 ratio continuously decreased from 3.4 to 0.1 when D/W was increased from 0 to 0.1 and then remained nearly constant at 0.1 as D/W was increased from 0.1 to 0.5. The value of $τ_{max c i r c u l a r}$ became approximately 0.1 times greater than $τ_{max w o v e n}$ as D/W was increased from 0 to 0.1. The SR6 ratio continuously decreased from 7.85 to 0.1 when D/W was increased from 0 to 0.1. It then remained nearly constant at 0.1 when D/W was increased from 0.1 to 0.5. A maximum plastic strain value of 7.4 × 10⁻² was obtained in the circular disk when D/W was equal to 0.5. This can be seen in Figure 10. In the case where D/W was equal to zero, $τ_{r}_{θ}$ -induced strain value for a disc with a circular fiber array was 1 × 10⁻² and for a disc with a radial fiber array was 4.5 × 10⁻³. As D/W was increased from 0 to 0.1 (surface hole drilled), the strain value remained nearly constant for both circular and radial discs. As D/W was increased from 0.1 to 0.5, the strain value continuously increased in the circular disc from 1 × 10⁻² to 7.4 × 10⁻². In the radial array disc, there was an initial decrease to a value of 3.6 × 10⁻³ as D/W was increased to 0.3, and from there a continuous increase to 11 × 10⁻³ as D/W rose from 0.3 to 0.5.

When D/W is equal to zero, the discs have no surface hole, the thermal stresses occurring in the disc with a woven fiber array were always smaller than the other discs, due to the same coefficient of thermal expansion in the radial and tangential directions under thermal loading.

In general, as D/W increased from 0 to 0.1 as surface holes were drilled, SR1 and SR3 ratios increased for $σ_{r}$ and $σ_{θ}$ stress values, as seen in Figures 6 and 8; whereas SR2 and SR4 ratios decreased, as seen in Figures 7 and 9. However, for $τ_{x y}$ and $τ_{r θ}$ stress values, both SR5 and SR6 ratios decreased, as seen in Figures 10 and 11. This means that when surface holes are drilled, $σ_{r}$ and $σ_{θ}$ stresses are diminishing in circular array discs and increasing in radial array discs in proportion to $σ_{x}$ stress occurring in the disc with woven fiber array. The lowest stress values appear in the disc with woven fiber array. On the other hand, $τ_{r θ}$ stresses increased in both the disc with the circular fiber array and the disc with the radial fiber array in proportion to $τ_{x}$ occurring in the disc with a woven fiber array. As D/W was increased from 0.1 to 0.3 due to surface holes being enlarged, the SR1 ratio increased (Figure 6), the SR2 ratio decreased (Figure 7), and both the SR5 and SR6 ratios remained constant (Figures 10 and 11). This means that $σ_{r}$ diminishes in discs with a circular fiber array as $σ_{θ}$ stress increases in the disc with the radial fiber array and that no change takes place for $τ_{x}$ and $τ_{r θ}$ stresses regardless of the D/W ratios. Based on this study executed on the developed composite, discs exposed to convectional cooling loading following salient main results can be revealed. Despite the results for the disc with woven fiber array, the thermal stresses exceeded the yield stress for discs with circular and radial fiber arrays. It can be explained by the properties of the disc that are equal in both tangential and longitudinal directions in woven fiber configuration. Thermal residual stress intensifies both inner diameter and around hole at the center of discs whatever fiber array has and their magnitudes are diminishing when D/W ratio is increased from 0 to 0.5. Plastic strains existed in the discs with both circular and radial fiber arrays except for woven fiber array. Magnitudes of the strains have smaller values in disc with radial fiber array rather than in disc with circular fiber array. Comparatively low thermal stress and plastic strain occurs in discs with radial fiber array compared to circular fiber array so it has more advantages in terms of thermal loading.

Conclusion

The values of the stresses and strains may be different depending on the distinct manufacturing technique and fiber volume concentration of TPC discs and diameter of the hole drilled on them. In order to compare their advantages and disadvantages of fiber array, same fiber volume ratio is considered in this study.

Woven fiber reinforcement is competent configuration in terms of reducing the formation of thermal stress and plastic deformation when compared to circular and radial fiber reinforcement. Comparison of TPC materials with radial and circular fiber array including same amount of fiber reinforcement confirmed that the radial is a better choice for thermal stress and strain point of view.

Existing of the circular holes on the discs have reducing effect on the thermal residual stresses and plastic deformations and as the diameter of hole increases, the values of the residual stresses and the plastic deformations decrease.

Finally, it can be understood that configuration of fiber array and hole dimension have an important influence on the thermal stresses and plastic deformations. It is possible to make new studies in which the effects of the different fiber arrays, composite geometries, materials, and load conditions may have been considered for TPC materials.

Footnotes

Funding

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

References

van Barschot

Airbus en Boeing vliegen met Ten Cate. In: NRC Handelsblad, Rotterdam, 5 January 2005, p. 17.

Kugler

Moon

. Identification of the most significant processing parameters on the development of the fiber wavines in thin laminates. J Compos Mater 2002; 36(12): 1451–1479.

Nairn

. Matrix microcracking in composites. In: Talreja

Manson

(eds) Polymer matrix composites. Oxford, UK: Elsevier Science Ltd., 2001, pp. 403–432.

Hsuch

Lee

Lin

. Analyses of mode I edge delamination by thermal stress in multilayer systems. Compos B: Eng 2006; 37(1): 1–9.

Wijskamp

Shape distortions in composite forming. PhD Thesis, Mech. Eng., University of Twente, Enschede, The Netherlands, 2005.

Unger

Hansen

HYS

. Method of reducing residual stress in thermoplastic laminates. European Patent, USA, 1993.

Nairn

Zoller

. Matrix solidification and the resulting residual thermal stress in composites. J Mater Sci 1985; 20(1): 355.

Favre

. Residual thermal stresses in fibre reinforced composite materials: a review. J Mech Behav Mater 1988; 1 (1–4): 37–53.

Parlevliet

Bersee

HEN

Beukers

. Residual stresses in thermoplastic composites—a study of the literature, part I: formation of residual stresses. Composite A 2007; 37(11): 1847–1857.

10.

Barnes

Byerly

. The formation of residual stress in laminated thermoplastic composites. Compos Sci Technol 1994; 51(4): 479–494.

11.

Kim

Hahn

Croman

. The effect of cooling rate on residual stress in a thermoplastic composites. J Compos Technol Res 1989; 11(2): 47–52.

12.

Sayman

Arman

. Thermal stresses in a thermoplastic composite disc under a steady state temperature distribution. J Reinforc Plast Compos 2006; 25: 1709.

13.

Sayman

. Stress analysis of a thermoplastic composite disc under uniform temperature distribution. J Thermoplast Compos Mater 2006; 19: 61.

14.

Sayman

Pekbey

. Elastic plastic stress analysis of a thermoplastic composite disc under linear temperature distribution. J Reinforc Plast Compos 2006; 25: 379.

15.

Şen

Pekbey

Sayman

. Elastic–plastic stress analyses of a thermoplastic composite disc under parabolic temperature distribution. Indıan J Eng Mater Sci 2007; 14: 282–288.

16.

Weng

. Influence of thermal residual stresses on the composite macroscopic behavior. Mech Mater 1998; 27: 229–240.

17.

Huang

. Strength formulae of unidirectional composites including thermal residual stresses. Mater Lett 2000; 43 (1–2): 36–42.

18.

Hobbiebrunken

Hojo

Jin

. Influence of non-uniform fiber arrangement on microscopic stress and failure initiation in thermally and transversely loaded CF/epoxy laminated composites. Compos Sci Technol, 2008; 68 (15–16): 3107–3113.

19.

Jin

. Effect of fiber arrangement on residual thermal stress distributions in a unidirectional composite. J Compos Mater 2007; 41: 2113–2130.

20.

Tsai

J-L

Chi

Y-K.

The thermal residual stress effect on the constitutive behaviors of fiber composites. Compos B: Eng 2008; 39(4): 714–721.

21.

On stress concentrations for isotropic/orthotropic plates and cylinders with a circular hole. Composites Part B 2003; 34(2): 127–134.

22.

Allam

MNM

Zenkour

. Stress concentration factor of a structurally anisotropic composite plate weakened by an oval opening. Compos Struct 2003; 61(3): 199–211.

23.

Mittelstedt

Becker

. A single-layer theory approach to stress concentration phenomena in layered plates. Compos Sci Technol 2004; 64(10-11): 1737–1748.

24.

Toubal

Karama

Lorrain

Stress concentration in a circular hole in composite plate. Compos Struct 2005; 68(1): 31–36.

25.

Ting

Chen

Yang

. Stress analysis of the multiple circular holes with the rhombic array using alternating method. Int J Pres Ves Pip 1999; 76(8): 503–514.

26.

Jones

RM.

Biaxial Strength Criteria for an Orthotropic Lamina, Mechanics of Composite Materials. Second Edition, 1999, 109–112, Taylor & Francis Inc., Philadelphia, PA, USA.