Material constants and correlation coefficients -60 °C ≤ T ≤ 23 °C
Fiber reinforced composite materials have been increasingly used as structural members in many structures such as airplane. The advantages of these materials are derived from their high strength, stiffness and damping together with low specific weight. Low temperature mechanical properties of glass fiber-reinforced epoxy have to be assessed, because of composite materials are subjected to low temperatures in service. Experimental or analytical investigation on the tensile failure behavior of glass/epoxy laminated composite with/or without stress concentration subjected to thermo-mechanical static loadings at low temperatures has not been done yet. In the present work, a model was developed to perform the progressive failure analysis of quasi isotropic composite plates at low temperatures. The initial failure load is calculated by means of an elastic stress analysis. The load is increased step by step. For each given load, the stresses are evaluated and the appropriate failure criterion is applied to inspect for possible failure. For the failed element, material properties are modified according to the failure mode using a non-zero stiffness degradation factor. Then, the modified Newton–Raphson iteration is carried out until convergence is reached. This analysis is repeated for each load increment until the final failure occurs and the ultimate strength is determined. The present method yields results in a reasonable agreement with the experimental data at room temperature and -60
2. Materials and mechanical testing
2.1. Material properties
Unidirectional glass fibers have been used in this investigation as reinforcing material, while epoxy resin has been utilized as the matrix material. Hand lay-up method was used to fabricate thin laminates with epoxy resin ML-506 with hardener HA11. Test specimens were cut from laminates according to relevant standard codes. All specimens had a constant cross section with tabs bonded to the ends. The fiber volume fraction of the composites was 55%.
2.2. Mechanical test equipments
The specimens are tested from beginning to the failure mode both at room and low temperatures. All tests were conducted under displacement control condition using Instron 5582 machine adaptable for cryogenic service by an environmental chamber. The displacement rate was 2 mm/min. Environmental chamber has the ability to cool down the temperature to -196 °C by evaporating liquid cryogenic Nitrogen. During the tests, a pressurizing device is used to control the cooling time from room temperature to -20 °C and -60 °C and maintain an evaporating pressure of 152 kPa. In order to reach thermal equilibrium in mechanical parts of machine one hour prior to first time of testing, temperature of environmental chamber was kept at desired low temperatures. Since composite materials require more time than metal materials to reach thermal equilibrium states in low temperatures, specimens stayed in a constant low temperature for at least an hour prior to start loading. Environmental chamber was equipped with a fan to perform a uniform temperature. A digital thermometer was mounted inside of environmental chamber to monitor active temperature during the tests. Experimental set up for mechanical tests at different temperature is shown in Figure 1.
3. Material characterization
In order to apply progressive damage modeling to a laminated composite, mechanical properties of unidirectional laminates should be assessed. For this purpose test specimens was fabricated for tensile, compressive and shear loading in longitudinal and transverse direction. In this section, for each case, specimen dimensions and experimental results are reported. Also some master curves are illustrated to predict mechanical properties in the assumed temperature range. Finally all of mechanical properties are assessed by a regression function.
3.1. Tensile test
The unidirectional glass fiber-reinforced epoxy which composed of ten plies was used in this test; the total thickness of the laminate is 2
The purpose of the static tensile tests was to determine the low temperature effects on the tensile strength, Young modulus and ultimate strain of unidirectional laminates in longitudinal and transverse directions. Figure 3 show the typical stress-strain behavior of unidirectional laminate under tensile loads in longitudinal direction at room temperature, -20
Figure 4 show the load-deflection curve of unidirectional laminates under tensile loads in transverse direction. The behavior is the same as described in the previous tests. The only difference is that the stress-strain curve in this case showed insignificant nonlinearity before reaching the maximum stress, which is due to the nature of the plasticity of epoxy. It was observed because the composite becomes more brittle in low temperatures, this nonlinearity behavior of matrix decreased by decreasing temperature in Figure 4 and Ref. . Typical results was also reported that both transverse tensile strength and stiffness increased by decreasing temperature from RT for glass-epoxy composites .
Failure regions of the composite specimens at different temperatures in longitudinal and transverse directions under tensile loading are shown in Figure 5 and Figure 6 respectively. It is found that in the case of longitudinal tension, the failure mechanism changes with temperature while no significant change was observed on the transverse direction.
For the longitudinal direction, from a visual inspection of the damage surfaces, a significant change could be observed in the fracture surface at room temperature to -60
3.2. Compression test
The specimens used for compression test had the same material as those of the tensile test except that these specimens were made of fifteen plies so that the total thickness is 3
In order to determine the effect of low temperature service on the strength and strain to failure of composite specimens in longitudinal and transverse directions, a set of experimental tests were performed. For this reason, a special fixture recommended in ASTM D3410 was used to prevent possible buckling. Figure 8 illustrate typical stress-strain behavior of composite laminate under compressive loading in longitudinal direction at various temperatures. As shown, the stress-strain relation is slightly nonlinear and final failure occurs catastrophically in all cases. The main specific property of longitudinal compressive stress-strain curves is that against other nonlinear stress-strain curves of UD composite in longitudinal or transverse directions, these curves exhibit positive curvature. This may be due to interface bonding between matrix and fibers. Imperfect fiber-matrix bonding is a manufacturing defect that is created during curing process. This defect influences the role of the matrix that supports the fibers when the composite is under compression , . From stress-strain curves, it seems that by increasing compression in longitudinal direction, this bonding becomes grater and the role of matrix against microbuckling of fibers is more highlighted. As shown in the figure, significant non-linear deformation was often observed before the maximum load , , and this behavior was associated to the plastic deformation of the polymeric matrix. Compressive strength in longitudinal direction increases significantly by decreasing temperature (from 570.37
Figure 9 depict stress-strain curves of unidirectional composite under transverse compressive loading at different temperatures. Again significant increase in strength of composite was observed by decreasing temperature (from 122.12 MPa at room temperature, to 186.22 MPa at -60 °C). Typical behavior was reported by  for glass-epoxy composites. But in this case against other directions of loading, by decreasing temperature, transverse compressive stiffness was decreased and strain to failure was increased significantly. The reason of the observed phenomenon in this type of loading is unknown and it is interested to investigate more about this behavior in future works.
The main failure mode of specimens in this case is micro buckling of fibers in shear or extensional mode. Also, in this case, there is a small amount of tab deboning near the gauge area for room temperature specimens. The change in the mechanical response with decreasing temperature is associated with a change in failure modes. As shown in Figure 11, failure occurred at the middle of specimen under micro-buckling failure mode. The difference is that, in low temperature cases, more fiber pull-out was occurred. Also, extensive deboning between fibers and the matrix was also observed because of weaker fiber and matrix bonding at low temperatures. In this case, only at low temperatures failed specimens fractured in two separate parts. In the case of compression failure in the transverse direction at low temperature, the failure area is much more than the room temperature specimens. Also at room temperature, specimen is not fractured in two separate parts. This is because of major effects of low temperature on mechanical properties of epoxy matrix of composites. In this case, there is no significant change in failure mechanism for specimens at -20
3.3. In-plane shear test
In the following, method of static experiments for characterizing material properties under in-plane shear loading is illustrated. Among numerous testing methods, the three-rail shear test method, described by the ASTM A4255 , is a fairly reliable technique. This method is modified by Shokrieh . He developed new specimen by inserting notches into the edges of specimen. The reason behind the insertion of notches at the locations of stress concentrations is to replace a very sharp crack with a blunt crack with much lower stress concentration factor. To increase the stability of the specimen during the test, a [0,90]s configuration is selected in this test. The specimen and the dimension are shown in Figure 12.
Electric strain gauges were used to measure the strains of the specimens during the shear tests. The strain gauges are attached to the specimens in the direction of 45°. The strain gauges and lead wires were manufactured by TML Corporation for low temperature ranges. However, since the changing in temperature would also affect the strain in strain gauges, Wheatstone half bridge circuit has been used to reduce this unfavorable effect.
Figure 13 show the stress-strain behavior of modified specimen under in-plane shear loadings at room temperature, -20
For metallic materials, the linear region of the stress-strain curve corresponds to elastic deformation and the nonlinear region to plastic deformation. However, when working with non-crystalline substances (such as epoxy), the glass transition temperature (
3.4. Master curves for mechanical properties
Figure 15 show the experimental values of tensile and compressive strengths for unidirectional laminates at different temperatures in longitudinal and transverse directions with their calculated standard deviations.
In each case in the above figures, by a polynomial curve fitting to experimental data, the increasing trend of strengths by decreasing the temperature is shown. However, the increasing rate of compressive strength by decreasing temperature is more for that of the tensile strength in both longitudinal and transverse directions. Figure 16 show experimental values of tensile modulus for unidirectional laminates at different temperatures in longitudinal and transverse directions. Both longitudinal and transverse tensile modulus increased by decreasing temperature to -60
Decreasing temperature, in this case, has lower effect on the final failure mode of the compression tests of the unidirectional plies in longitudinal and transverse directions at various temperatures in comparison with the previous one. Experimental magnitudes of shear strength and modulus at different temperatures are compared in Figure 17. In the case of shear loading, low temperature has major effect on mechanical properties of composites.
|Longitudinal compressive strength||0.0093||1.5706||604.05||0.9575|
|Transverse compressive strength||0.0121||0.3784||124.04||0.9678|
|Longitudinal tensile strength||0.0098||0.3808||718.52||0.9642|
|Transverse tensile strength||0.0010||0.0789||70.390||0.9622|
|In-plane shear strength||0.0058||0.0341||65.339||0.9682|
|In-plane shear modulus||0.0002||0.0184||2.477||0.9769|
4. Specimen geometry
Quasi-isotropic lay-up ([0/±45/90]2s) was used in this study for tensile tests at room temperature and -60 °C. For this reason, thin laminate composed of ten plies of reinforcement with epoxy resin were fabricated with considered configuration, giving a laminate approximately 2 mm in thickness. For laminate with stress concentration, a central hole was made by a machine. Woven glass/epoxy tabs with tapered ends were locally bonded on each side of the specimens. These tabs allow a smooth load transfer from the grip to the specimen especially for low temperature test. All specimens had a constant cross section with tabs bonded to the ends. The geometry of the specimen with stress concentration for tensile tests is shown in Figure 18.
5. Progressive Damage Modeling (PDM)
5.1. Stress analysis
The first component of the PDM is finite element stress analysis. Consider a composite plate with/without stress concentration (Figure 18). The plate has width w = 25 mm, length l =170 mm, thickness t = 2 mm and central hole D = 6 mm (if applicable). The plate is a laminated composite with quasi-isotropic ply orientation. A two dimensional macro code by APDL of ANSYS  is developed to perform finite element analysis. In this paper, the 8-node layered element SHELL 99 is adopted to model the laminates, which allows up to 250 different material layers in the thickness direction in each element without much increase of counting time. Mechanical properties of a unidirectional laminate at room temperature and -60°C, which is used as initial values in FEM and listed in Table 2, were tested by the present authors .
|Longitudinal tensile strength
|Longitudinal compression strength
|Transverse tensile strength
|Transverse compression strength
In this table, the script x refers to the fiber direction, and y refers to direction perpendicular to the fiber direction. In this study, the stress resultants are defined as follow:
where and are the mechanical and thermal stresses respectively. Thermal stresses are due to decreasing temperature from room temperature to -60 °C. These stress resultants will be calculated by the following relation:
In the above equation, is transformed reduced stiffness matrix for a laminated composite and is thermal strain vector which is defined as follow:
where and are temperature difference from room temperature and coefficients of thermal expansion for an angle ply laminate, respectively. can be given in terms of the coefficients of thermal expansion for a unidirectional laminate as:
and is a transformation matrix for an angle ply laminate .
5.2. Failure criterion
The second part of PDM is failure analysis. By using finite element results, at layer level in each element, stiffness reduction is carried out considering five types of damages: fiber and matrix in tension and compression and fiber-matrix shearing modes. To detect them, a set of two dimensional stress based failure criterion is selected. The following Hashin criteria  are used to detect five different failure modes. The first two of the failure modes are catastrophic and the others are not.
|Fiber tensile failure|
|Fiber compressive failure|
|Matrix tensile failure|
|Matrix compressive failure|
|Fiber-Matrix shearing failure|
5.3. Material properties degradation rules
The last component of PDM is material properties degradation. As failure occurs in a unidirectional ply of a laminate, material properties of that failed ply are changed by a set of sudden material property degradation rules. In the present method after failure occurrence in a ply of the laminate, instead of inducing real crack, the failed region of the unidirectional ply is replaced by an intact ply of lower material properties. A complete set of sudden material property degradation rules for all various failure modes of a unidirectional ply under a uni-axial static stress is explained in the following. The rules must be carefully applied to avoid numerical instabilities during computation by the computer program.
Fiber tension failure
Fiber tension failure mode of a ply is a catastrophic mode of failure and when it occurs, the failed material cannot sustain any type or combination of stresses. Thus, all material properties of the failed ply are reduced, as follows:
Fiber compression failure
Fiber compression failure mode of a unidirectional ply is a catastrophic mode of failure and when it occurs, the failed material cannot sustain any type or combination of stresses. Thus, all material properties of the failed ply are reduced. Equations 7 and 8 show this degradation rule.
As mentioned, these two modes of failure are catastrophic, therefore if it occurs, the other modes of failure do not need to also be verified.
Matrix tension failure
Matrix compression failure
Matrix compression failure mode results in the same type of damage to the composite ply as the matrix tension failure mode. This mode of failure is not catastrophic; therefore, other material properties are left unchanged:
Fiber-Matrix shear failure
In fiber-matrix shearing failure modes of a ply, the material can still carry load in the fiber and matrix directions, but in-plane shear stress can no longer be carried. This is modeled by reducing the in-plane shear material properties of the failed ply, as follows:
The PDM is an integration of the three important components: stress analysis, failure analysis and material property degradation. The model is capable of simulating the first and final failure load of composite laminates with arbitrary geometry and stacking sequence under tensile static loading at room temperature and -60 °C.
A computer program, the algorithm of which is shown in Figure 19, is established to analyze the failure mechanism of composite plates at low temperatures using APDL of ANSYS. All material properties are set to initial values which are experimentally evaluated by present authors . The initial failure load is calculated by means of an elastic stress analysis. The load is increased step by step. For each given load, the stresses at each integration point are evaluated and the appropriate failure criterion is applied to inspect for possible failure. At the point with failure, the material properties are modified according to the failure mode using a non-zero stiffness degradation factor. Then, the modified Newton–Raphson iteration is carried out until convergence is reached. The convergence tolerance is assumed to be 0.001. This analysis is repeated for each load increment until the final failure occurs and the ultimate strength is determined.
Theoretically, the smaller load increment between successive steps, the more accurate analysis result can be achieved. However, a reasonable load increment should be prescribed to avoid too much analysis time and to ensure accuracy. After sensitivity analysis on load increment, 1 KN is applied in the current model.
6. Results and discussion
The specimens are tested under static tensile loading at room and low temperatures. In each case (with or without stress concentration) at room temperature four specimens and at -60 ºC five coupons were tested to show statistic scatter of experiments. By statistical evaluation (mean values and standard deviation) reliability of results were appraised. The experimental setup for low temperature tests using an environmental chamber is shown in Figure 20.
The average tensile properties such as first ply failure (FPF) load, final failure (FF) load and ultimate strain to failure (USF) for quasi-isotropic laminates with/without stress concentration determined based on results of stress–strain curves from experimental tests and numerical analysis are summarized in Table 4.
Figure 21 and Figure 22 show the failure process predicted by the model at room temperature and -60 ºC, respectively. At the FPF load, a mainly obvious damage around the hole of plate is matrix cracking (see Figure 21a and Figure 22a). By increasing the load, other failure modes are also occurred (see Figure 21b and Figure 22b). At the final failure load, the plate breaks along the transverse direction through the central hole edge, the same as noticed in the experimental tests. In this load, the mainly failure mode is fiber breakage (see Figure 21c and Figure 22c). As shown in the following figures, the failure regions of specimens at -60 ºC at each step is much more than room temperature.
In all cases, major failure mode was fiber breakage, matrix cracking and fiber-matrix shearing which are shown in the figures by red, pink and blue paints respectively. Other failure modes are also occurred in the final failure but can not be shown in the figures.
Figure 23 illustrates mean values of tensile strength for quasi-isotropic laminate at room temperature and -60 °C with and without stress concentration. This figure also compares experimental results with those obtained from the present finite element model. Results show that strength of laminate increased by decreasing temperature. This is because of change of micromechanical properties of composites at low temperature.
Figure 24 shows typical stress-strain curve for the laminate with stress concentration based on experimental results at room temperature and -60 °C.
Failure mechanism of tested specimens with central hole at room temperature and -60 °C are different. Figure 25 shows failed specimens at two different temperatures. From a visual inspection, there is a small amount of tab de-bonding near the gage area for both two cases with more fiber pull-out for low temperature specimen. At low temperature, because of the interface between fiber and matrix are much weaker and the fiber debond from the matrix, synchronous with fiber breakage, matrix cracking and a few fiber-matrix shearing were occurred. However, the mainly failure mode for all cases is fiber breakage and matrix cracking.
Tensile failure behavior of glass/epoxy laminated composite subjected to thermo-mechanical loadings at low temperatures with/without stress concentration was investigated experimentally and numerically. A finite element code was utilized to model the progressive failure analysis of quasi-isotropic composite plates at low temperatures under static loading. For each given load step, the stresses at each integration point are evaluated and the appropriate failure criterion is applied to inspect for possible failure by using Hashin failure criteria. At the point with failure, the material properties are modified according to the failure mode using a non-zero stiffness degradation factor. In case of failure detection, because of nonlinear behavior, the modified Newton–Raphson method was carried out until convergence is reached. This analysis is repeated for each load increment until the final failure occurs and the ultimate strength is determined. Based on the results of the present study, the following conclusions can be drawn:
The stress-strain behavior of UD laminate under tensile loads in longitudinal direction was linear elastic until breakage and the slope of the stress-strain curve and the strength increased about 12% as the temperature decreased to -60 °C. On the other hand, by decreasing temperature, strain to failure decreased slightly about 10%. However, the laminate under transverse tensile loading exhibited insignificant nonlinearity (especially at room temperature) before reaching the maximum stress, which is due to the nature of the plasticity of epoxy.
The UD laminate showed a nonlinear elastic relation between stress and strain under compressive loads in both directions until final failure occurred catastrophically. The compressive strength increased by decreasing temperature (for longitudinal direction 28% and transverse direction 50%) while the strain to failure decreased. But in transverse direction, strain to failure increased by decreasing temperature and against other nonlinear curves, positive curvature was observed.
The UD laminate under in-plane shear loading behaves highly nonlinear for all temperature tests until final failure. Again, in this case, both shear strength and stiffness increased by decreasing temperature in about 32% and 70% respectively. Also nonlinear region decreased by decreasing temperature due to increasing brittleness of epoxy matrix.
Failure type of UD laminates under various loadings was affected by low temperature. It was found that, by decreasing temperature a small amount of tab deboning occurred near the gauge area of specimens in longitudinal tensile loading. Also, because of the interface between fiber and matrix was much weaker at low temperature, fibers deboned the matrix in all test cases. Therefore, it may be concluded that the low temperature affects the micro mechanisms of damage in composite specimen.
From general master curves and illustrated regression function, mechanical properties of unidirectional glass fiber polymeric composites at temperature range of -60 °C to 23 °C can be evaluated.
The stress-strain behavior of laminate under tensile loads was linear elastic until first ply failure (FPF). After this, the behavior of laminate was nonlinear until final failure occurred. This trend was observed for laminated composite with/without stress concentration at both temperatures.
The slope of the stress-strain curve and the strength of laminate increased as the temperature decreased to -60 °C. On the other hand, by decreasing temperature, strain to failure decreased slightly. So, in spite of improvement in strength and stiffness of composites under static loading at low temperatures in comparison with room temperature, their strain to failure under these environmental conditions becomes weaker.
The failure mode of laminated composite at low temperature changes from matrix cracking at FPF to mixed mode failure (fiber breakage, fiber matrix shearing and matrix cracking) at final failure load.
Failure type of laminates under various loadings was affected by low temperature. It was found that, by decreasing temperature a small amount of tab deboning occurred near the gage area with more fiber pull-out. Also, due to weakness of the interface between fiber and matrix at low temperature, fiber debones the matrix. Therefore, it may be concluded that the lower temperature affects the micro mechanisms of damage.
Good agreement was achieved between results from experimental and analytical calculation at room temperature and -60 °C. This agreement also showed the validity of model.
D 3039M-95a, Standard test method for tensile properties of polymer matrix composite materials, Astm D 1997
D 3410M-95, Standard test method for compressive properties of polymer matrix composite materials with unsupported gage section by shear loading, Astm D 1997
D 4255M-83, Standard guide for testing in-plane shear properties of composite laminates, Astm D 1994
Properties of composite materials for cryogenic application, Cryogenics. 38, 3-12 ( Schutz J. S 1998
and A. Wanki Jun, Compressive strength of unidirectional fiber composites with matrix non-linearity, Composite Science and Technology, C. T Sun 52 1994 577 587
Damage trends in cryogenically cycled carbon/polymer composites, Composite Science and Technology. 64, 1773-1784 ( Bechel V. T Kim R. Y 2004
Experimental and analytical studies on the damage initiation in composite laminates at cryogenic temperature, Composite Structure. 76, 62-66 ( Kim R. Y Donaldson S. L 2006
Residual stress and thermal expansion of graphite epoxy laminates subjected to cryogenic temperatures, Composite Science and Technology. 66, 2449-2455 ( Ifju P Myers D Schulz W 2006
Mechanical response of a unidirectional graphite fiber/polyimide composite as a function of temperature, Composite Science and Technology. 66, 1045-1055 ( Rupnowski P Gentz M Kumosa M 2006
Tensile response of graphite/epoxy composite at low temperatures, Composite Structures. 79(1), 84-89 ( Kim M. G Kang S. G Kim C. G Kong C. W 2007
Three-dimensional thermoelastic analysis of cracked plain weave glass/epoxy composites at cryogenic temperatures, Composite Science and Technology. 64, 2353-2362 ( Takeda T Shindo Y Narita F 2004
Double cantilever beam measurement and finite element analysis of cryogenic Mode I interlaminar fracture toughness of glass-cloth/epoxy laminates, Journal of Engineering Materials and Technology. 123, 191-197 ( Shindo Y Horiguchi K Wang R Kudo H 2001
Mode I fracture toughness of an adhesively bonded composite-composite joint in a cryogenic environment, Composite Science and Technology. 67(3-4), 501-6 ( Melcher R. J Johnson W. S 2007
Characterization of Mode I fatigue crack growth in GFRP woven laminates at low temperatures, Acta Materialia. 53, 1389-1396 ( Shindo Y Inamoto A Narita F 2005
Mode I fatigue delamination growth in CFRP woven laminates at low temperatures, Engineering Fracture Mechanics. 73, 2080-2090 ( Shindo Y Inamoto A Narita F Horiguchi K 2006
Tension-tension fatigue behavior of GFRP woven laminates at low temperatures, Cryogenics. 45, 123-128 ( Kumagai S Shindo Y Inamoto A 2005
Cryogenic fatigue behavior of plain weave glass/epoxy composite laminates under tension-tension cycling, Cryogenics. 46, 794-798 ( Shindo Y Takano S Horiguchi K Sato T 2006
Interaction of damage failure and post-buckling behavior of composite plates with cut-outs by progressive damage modeling, Composites: Part B. 39(2), 304-15 ( Labeas G Belesis S Stamatelos D 2008
Progressive failure analysis of bonded composite repairs, Composite Structures. 81, 331-340 ( Liu X Wang G 2007
Progressive failure of triaxial woven fabric (TWF) composites with open holes, Composite Structures. 65, 419-431 ( Zhao Q Hoa S. V Ouellette P 2004
Deformation and progressive failure behavior of woven-fabric-reinforced glass/epoxy composite laminates under tensile loading at cryogenic temperatures, Composite Science and Technology. 65, 1691-1702 ( Takeda T Takano S Shindo Y Nurita F 2005
Tensile and damage behavior of plain weave glass/epoxy composites at cryogenic temperatures, Fusion Engineering and Design. 81, 2479-2483 ( Shindo Y Takano S Narita F Horiguchi K 2006
Progressive failure analysis of thick composite plates using spline finite strip method, Composite Structures. 79, 34-43 ( Akhras G Li W. C 2007
Progressive failure analysis of composite plates, Proceedings of 8th Iranian aerospace society conference. Esfahan, Iran ( Shokrieh M. M Torabizadeh M. A Fereidoon A 2009
ANSYSVer. 10, Canonsburg (PA): SAS IP; ( 2005
An investigation on damage of quasi-isotropic laminated composite, Proceedings of 18th annual international conference on mechanical engineering (In Persian), Tehran, Iran ( Shokrieh M. M Torabizadeh M. A Fereidoon A 2010
Mechanics of composite materials, Taylor and Francis Group, LLC ( Kaw A. W 2006
Progressive fatigue damage modeling of composite materials, Ph.D. Thesis, McGill University ( Shokrieh M. M 1996
Longitudinal compressive strength of a unidirectional fibrous composite, Composite Science and Technology, A. P Wilczynski 45 1992 37 41
Inelastic behavior of an AS4/PEEK composite under combined transverse compression and shear. Part I: Experiments, International Journal of Plastics, T. J Vogler S Kyriakides 15 1999 783 806
Strain rate behavior of composite materials, Composites: Part B, H. M Hsiao I. M Daniel 29 1998 521 533
Mechanical behavior of unidirectional fiber-reinforced polymers under transverse compression: Microscopic mechanisms and modeling, Composite Science and Technology, C Gonzalez J Llorca 67 2007 2795 2806
ANSYSVer. 10, Canonsburg (PA): SAS IP; 2005
An investigation on damage of quasi-isotropic laminated composite,” Proceedings of 18th annual international conference on mechanical engineering (In Persian), Tehran, Iran, ISME M. M Shokrieh M. A Torabizadeh A Fereidoon 2010
Mechanics of composite materials, Taylor and Francis Group, LLC, A. K Kaw 978-0-84931-343-1 2006
W.D., Materials Science and Engineering: An Introduction, Third edition, John Wiley & Sons, New York, J. R Callister 480 1994
Low temperature mechanical testing of carbon-fiber/epoxy-resin composite materials, NASA Technical Paper, 3663, and A. T Nettles E. J Biss 1996
Experimental study on mechanical behavior of laminates at low temperatures, Cryogenics, M Gong X. F Wang J. H Zhao 47 2007 1 7
Low temperature properties of a unidirectionally reinforced epoxy fiberglass composite, Cryogenics, R. P Walsh J. D Colskey R. P Reed 35 1995 723 725
Physically-based failure models and criteria for laminated fiber-reinforced composites with emphasis on fiber-kinking. Part I: development, Composites: Part A, S. T Pinho L Iannucci P Robinson 37 2006 63 73
Failure analysis of FRP laminates by means of physically based phenomenological models, Composite Science and Technology, A Puck H Schurmann 62 2002 1633 1662
Fracture micromechanisms in C/epoxy composites under transverse compression, Master thesis, Universidad Politecnica de Madrid, D Aragones 2007