Research Article  Open Access
Validating the Classical Failure Criteria for Applicability to the Notched WovenRoving Composite Materials
Abstract
The classical failure criteria are phenomenological theories as they ignore the actual failure mechanism and do not concentrate on the microscopic events of failure. The main objective of the current investigation is to modify the classical failure theories to comprise the essential failure mechanism (interfacial shear failure) in the thinlayered wovenroving composite materials. An interfacial shear correction factor (MH_{6}) is introduced into the nondimensional shear terms in the studied classical failure criteria. Thus the validity of applying these theories to the investigated material will be augmented. The experimental part of the current study is conducted on thinlayered circular specimens. The specimens are fabricated from two plies of fiber Eglass wovenroving fabric reinforced with polyester. The fabrics are laid to have [45°] or [0°, 90°] fiber orientation. The specimens used are plain, where no macroscopic sources of stress concentration exist or having circular notches of five, seven, or nine mm radii. The specimens are subjected to low cycle completely reversed fatigue bending loading where the SN and the R.D.N curves are plotted for each group of specimens.
1. Introduction
The failure modes in wovenroving composite materials are more complex than those of other manmade materials or metals [1]. To improve the performance of these materials, the source of damage must be insightfully understood to prevent its subsequent evolution and the deterioration of fundamental mechanical properties. The adhesion integrity between fiber and matrix is the justification for superior composite materials properties over their constituent components as discussed by a lot of researchers [2–5]. Their investigations were focusing on the fiber/matrix interaction at micro level scale [2–5]. Moreover, different techniques have been proposed to measure the fiber/matrix adhesion levels [6–9]. Swan et al. [5] and Jayaraman et al. [10, 11] have proposed different analytical methodologies and experimental techniques to improve the interpretation of the interface region. Dasgupta and Sirkis [12] studied the effect of coating modulus and thickness on the stress distribution around an optical fiber subjected to axial loading. Carman and Case [13] demonstrated that an optimal coating exists for the longitudinal shear problem. Carman et al. [14] analytically investigated the influence of finite size interphase coatings applied to the structural fiber on the stress distribution in continuous unidirectional composite subjected to transverse loading. They also proposed an optimization of fiber coating procedure to reduce the stress concentration at the fiber matrix interaction locations. At the macrolevel scale, Madhukar and Drazal [15–18] published a series of papers describing the considerable effect of the interphase/interface on macrolevel strength and stiffness properties.
In 1960s, most of the classical failure criteria were proposed prior to the above mentioned studies. In the current investigation, a bridge between the classical failure theories which is phenomenological in nature and the interfacial deficiencies as the main cause of failure on the insightful interpretation is built. An interfacial shear correction factor denoted MH_{6} is used to augment the shear stress in the nondimensional shear term in the classical criteria. The modified classical criteria can fairly estimate the relative damage of the tested material under the assigned loading type.
The rest of the paper is organized as follows. Section 2 presents the analytical analysis. Section 3 provides the experimental work. The experimental data used in our previous investigation [1] are also used for the current study. The validation procedure of the studied classical criteria is explained in Section 4. Finally, the paper is concluded in Section 5.
2. Analytical Analysis
2.1. Stress Representation
As the specimens were subjected to completely reversed bending stress, the only existing stress component, , is the stress in the global axis (see (1)): where , for completely reversed loading and .
Based on the global stress calculated above, the local stress components could be calculated based on fiber orientation as follows. For stress element of [0°, 90°]_{2s } specimens, , . For stress element of [±45°]_{2s } specimens . Figure 1(a) shows these two stress states. Notice that the stress element is located at the midpoint of the middle length of the plain specimens as shown in Figure 1(a). For the notched specimen, the stress element was selected to be on the midline of the specimen length, above or below the notch where the moment loading value is uniform. These points are denoted with stars as shown in Figure 1(b). Notice that the distance between these points and the neutral axis is almost equal to the mean radius of the specimen cross section, such that the global stress at these points (denoted with stars in Figure 1(b)) is equal to (see (1)).
(a)
(b)
2.2. Theories of Failure
The diversity of the wovenroving composite material application made the need for the simple failure criteria to account for the microscopic failure modes in the composite structures a demanding task. Most of the modern failure criteria of composite material are proposed to account for the complexity of failure mechanisms in the composite materials which does not heavily exist in the thinlayered specimens used in the current study. Besides, the fatigue process in thinlayered wovenroving fiber reinforced with ductile matrices is similar to that in metals in the sense that it consists of two stages: crack initiation followed by crack propagation [1]. For the previous reasons, the classical failure criteria listed in Table 1 are focused on and the modern criteria are out of the current scope.

The classical failure criteria are not based on the physical mechanisms that are blamed for failure on the microscopic level at the time these theories were introduced. These theories are also known as phenomenological theories. Some of the recent forms of these classical criteria include normal stress interaction terms “” and stability condition factors which are used to account for the material anisotropy; see Table 1. Nevertheless, they did not reflect the dominant microscopic failure mechanism (interfacial shear failure) in the notched Wovenroving Fiber Reinforced Composite (WRFRC) thinlayered materials. Hence, an augmentation of the shear stress component is required to validate these criteria for applicability to the studied material.
As seen in Table 1, these criteria are based on the principal static stresses together with their corresponding static strength. In the present study, the principal fatigue stress components (, , and ) and the corresponding fatigue strengths (, , and ) are used in the failure criteria [19]. The values of the fatigue strengths , , and are calculated from tests conducted on fivelayer plates fabricated from the same constituent materials of the studied materials. The common observation in these theories is that the nondimensional shear stress component is not coupled with any of the other nondimensional local principal stresses.
These criteria could be further simplified based on the fiber orientation as shown in Table 2. The simplified criteria for [0°, 90°]_{2s } specimens do not contain the nondimensional shear stress term as the local shear stress component is zero. When applying those theories to [0°, 90°]_{2s } and [±45°]_{2s } specimens, we will consider the following.(i) = = the fatigue strength of plain [0°, 90°]_{2s } specimens subjected to completely reversed bending test. Consequently, Norris distortional energy [20], Hoffman [21], TsaiWu [22], TsaiHahn [23], Cowin [24], and Puppo [25] criteria will yield to the TsaiHill yield criterion [26] as shown in Table 2.(ii) = the shear endurance strength of plain [0°, 90°]_{2s } specimens subjected to completely reversed torsion test.
Notice that the simplified failure criteria of the [0°, 90°]_{2s } fiber orientation yield to the same form as shown in Table 2. The reason for this is the zero values of the local stress components and as discussed before.
3. Experimental Work
3.1. Testing Machine
A strain controlled testing machine designed by Abouelwafa et al., [27] Sharara [28], and Elshabasy [29] is used to conduct the current experiments. Figure 2 shows the two views of the testing machine. The testing machine was designed to perform pure bending, pure twisting moment loadings, or combined loading, whether in phase or out of phase, with a constant frequency of 8.75 Hz. The loading system of the machine is independent and is designed to produce the load in form of sine wave with variable stress ratio () between the minimum and maximum applied stresses. The bending system gives a fourpoint configuration. Therefore, the specimen is acted upon by a constant bending moment over its entire length. The applied bending moment is measured via a load cell. This cell is fixed on one of the two grippers and it has four active strain gauges, forming a full Wheatstone bridge.
3.2. Specimens
The thinlayered cylinder specimens were made of two layers of wovenroving E glass fabric, with equal glass intensity in both fill and warp directions. The reinforcing resin is SIROPOL 8330 Polyester. The Resin was treated with the methyl ethyl ketone peroxide as a catalyst in a percentage of 2.0 to 2.5% by volume, and the cobalt naphthenate (6% solution) is used as an accelerator in a percentage of 0.2% by volume [29–31]. Crosslinking and curing took place at the ambient conditions. The nominal dimensions of the tested specimen are shown in Figure 3. Chamis [32] showed that optimum load transfer between fiber and matrix occurs when the volume fraction varies between 0.55 and 0.65. Therefore, the manufactured specimens with fiber volume fraction () in the range of 0.5 to 0.64 by weight are accepted [1, 29]. Figure 3 shows the layout of the specimen with the two wooden plugs at the specimen end and rubber sleeve to relief and distributes the contact stress between the specimen ends and the copper gripper. The only counted failed specimens are those having failure signs at the middle third of the specimen length to avoid the failure emanating from the gripped ends.
3.3. Test Results
Before conducting the fatigue tests, five flat specimens made from five plies of the same constituent materials of the tubular specimens were subjected to static bending loading. The bending strengths for both fill and warp directions are very close to the ASTM standards of glass composite [33]. The dimensions of the specimen and the test results are listed in Table 3 [1].

In order to avoid any misleading results, only the specimens that had their failure features within the gauge section (the middle third of the specimen length) were considered. More than 110 specimens for each orientation were used in running the fatigue tests. Each group of the same orientation is divided into two groups, plain and notched. The specimens are subjected to completely reversed bending fatigue loading with different amplitudes. Each experimental plotted point is obtained by considering the average of three to five specimens tested under the same conditions [33]. The test results are used in plotting the corresponding  curve where they are fitted to the exponential form shown in
The test results of [±45°]_{2s } specimens are used in plotting their corresponding  curves, shown in Figure 4. Figure 5 shows the  curves of [0°, 90°]_{2s } specimens. Figures 6, 7, 8, and 9 are plotted for the sake of comparison between the [±45°]_{2s } specimens and the corresponding ones having (0°, 90°) fiber orientations.
3.4. Failure Modes
The failure modes for the different studied specimens are clearly shown in the zoomedin photos; see Figures 10 and 11.
For [0°, 90°]_{2s } plain specimen, the observed failure mode was a matrixcracking mode normal to the loading axis and along the transverse fiber direction as shown in Figure 10. This may be attributed to brittleness of the matrix compared to the filling fiber fabrics. The matrix cracks, propagates, and stops at the interface as the fiber has a higher strength than the polyester matrix. In case of [0°, 90°]_{2s } specimens with holes, the observed failure mode is also matrix cracking and slight interfacial shear failure along the fiber tangent to the notch and parallel to the loading axis, as shown in Figure 12.
For [±45°]_{2s } plain specimens, the observed failure mode is a combined mode of slight interfacial shear associated with slight matrix cracking; see Figure 11. For the [±45°]_{2s } specimens with notch, a clear white zone extended along the uncut fiber tangent to the notch edges is shown in Figure 13. Although the visual inspection on the microscopic level was missing, the white damaged zone around the notch was visually inspected carefully. The debond between the fiber filaments and matrix was clear. The justification for this is that the load is mainly transferred by the interfacial region between the fiber and matrix and the existence of shear stress component [1, 15–18]. Another attestation for our interpretation is the Xradiograph taking a graphiteepoxy [0_{2}/90_{2}]_{s} laminate having a center notch by Binienda et al. [34]. They noticed that the driving force was the interfacial shear due to load transfer from the fiber bundle cut by the hole to the fiber bundle which is uncut; see Figure 14.
4. Classical Theories Validation
All discussed failure criteria have their right hand side to be unity, and the left hand side is composed of normalized local stress components added together. The left hand side of each criterion will be called the calculated relative damage (R.D.). The validity of the manipulated criteria depends on the propinquity of its calculated R.D. value from unity. Using the reduced forms of the failure criteria, the calculated R.D. is plotted against the cycles to failure for the two fiber orientations.
For [0°, 90°]_{2s } specimen, the calculated R.D. is unity regardless of the specimen life and the used failure criteria; see Figure 15. The justification is that, in [0°, 90°] fiber orientation, the all theories of failure yielded to the same form shown in Table 2. Besides, the shear component which is blamed for the specimens’ failure does not exist in the reduced classical failure criteria for [0°, 90°]_{2s } fiber orientation.
In the [±45°] fiber orientation, the calculated relative damage values calculated by the above criteria listed in Table 2 are different. Figure 16 shows that the calculated R.D. using the Hill criterion [35] is very small compared to unity and almost constant regardless of the notch diameter or the number of cycles. The calculated R.Ds. are exactly the same for the criteria [20–26] listed in second row of Table 2; see Figure 17. The reason is that these theories yield to the same reduced form as shown in Table 2. Figure 18 shows that the R.D. calculated from the Norris interaction theory is the largest between the calculated R.D. values discussed before. It is also noticed that this value is almost constant regardless of the number of cycles or the notch diameter. From Figure 19, it is noticed that Ashkenazi criterion is overestimating the failure where the relative damage values are greater than one regardless of the number of cycles. Besides, the relative damage varies with the notch radii. Accordingly and for simplicity, this criterion is not considered for further investigation in the following section.
4.1. Interfacial Shear Correction Factor for [±45°]_{2s } Specimens
The simple classical phenomenological theories discussed above will be modified to account for the interfacial shear failure modes visibly shown in the failed specimens; recall Figures 10–14.
In the Hill criterion [35], it is noticed that the only surviving term is the , which is common in the all reduced forms of the [±45°]_{2s } specimens. Accordingly, the calculated R.D. value from this theory is expected and observed to be the smallest value as the local shear stress component does not include the augmentation of its value due to the overlooked interfacial shear failure mechanism. The R.D. value calculated from Norris interaction reduced criterion [36] is the greatest between the all calculated R.D. values except the refused Ashkenazi [37] as mentioned before. The justification for this is that the double value of the normalized local normal stress contributes to the R.D. value beside the shear term.
Also recall that the visual inspection of Figures 11 and 13 shows a white zone extended along the fiber tangent to the notch that fortifies the existence of interfacial shear failure, although the microscopic inspection is not implemented. Accordingly, the normalized local shear stress component will be augmented by a proposed correction factor. This factor is called interfacial shear correction factor and is denoted MH_{6}. The MH_{6} will be simply multiplied by the local shear stress component in the above discussed theories. Thus, the corrected failure criteria could fairly be used with the thinlayered [±45°]_{2s } WRFRC materials subjected to flexural fatigue stress. The only surviving term in the reduced criteria applied to the [0°, 90°]_{2s } specimens is ; hence the MH_{6} factor is no longer needed for validating these theories to the thinlayered [0°, 90°]_{2s } WRFRC materials subjected to flexural fatigue stress.
As the R.D. value differs from theory to another, the values of the MH_{6} factors are expected to be different as shown in Figure 20. From Figure 20, it is noticed that the smaller the R.D. value, the greater the interfacial shear correction factor. Besides, as the shear correction factor is used to account for the overlooked interfacial shear correction failure mechanism which is microscopic in nature, thus it should be independent of the dimension of the macroscopic source of stress concentration (circular notch). This is obvious in Figure 20 as the MH_{6} is independent of the normalized notch radius.
5. Conclusion
One of the microscopic failure mechanisms that the classical failure theories did not consider is the interfacial shear failure. The effect of overlooking this microscopic failure mechanism leads to poor applicability of these criteria to the wovenroving composite materials with [±45°] fiber orientation. As the local shear stress component, , is the driving force of the propagating interfacial debonds, the interfacial shear correction factor (MH_{6}) is proposed. This factor is independent of the notches’ dimensions. The classical failure criteria applied to the [0°, 90°]_{2s } specimens are not modified to account for the interfacial shear failure. The justification is that the [0°, 90°]_{2s } specimen behaves like isotropic materials under the loading conditions in the current investigation and the classical failure criteria of composite materials were based on isotropiclike material assumptions. Finally, under the completely reversed bending loading, the [0°, 90°]_{2s } fiber oriented wovenroving composite materials are more reliable than those of [±45°] fiber orientations.
Nomenclature
:  Maximum stress 
Number of cycles to failure  
:  Local normal stress in fiber longitudinal direction 
:  Local normal stress in fiber transverse direction 
:  Local shear stress 
:  Global normal stress in the loading axis 
:  Global normal stress normal to the loading axis 
:  Global torsion stress 
:  Fatigue strength in fiber longitudinal orientation 
:  Fatigue strength in diber transverse direction 
:  Compressive fatigue strength in fiber orientation, where or 2 
:  Tensile fatigue strength in fiber orientation, where or 2 
Applied bending moment  
:  Mean bending moment 
:  Variable bending moment 
Amplitude of applied bending moment  
Angular displacement  
Fiber orientation direction  
Mean radius of specimen  
Second polar moment of area of specimen  
Fatigue strength prediction factor for FRP  
Stress ratio between minimum and maximum value  
Radius of the circular hole  
Mean diameter of the thinlayered specimens  
:  Fiber volume fraction. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 M. M. Y. B. Elshabasy, M. N. Abou ElWafa, A. H. Hamdy, and A. A. Elmidany, “Comparison between inclusions and holes as stress kaisers in composite materials used in wind turbine components,” in Proceedings of the 45th AIAA Aerospace Sciences Meeting and Exhibit, AIAA20070819, Reno, Nev, USA, 2007. View at: Google Scholar
 P. S. Chua and M. R. Piggott, “The glass fiberpolymer interface: II work of fracture and shear stresses,” Composites Science and Technology, vol. 22, no. 2, pp. 107–119, 1985. View at: Google Scholar
 P. S. Chua and M. R. Piggott, “The glass fibrepolymer interface: III pressure and coefficient of friction,” Composites Science and Technology, vol. 22, no. 3, pp. 185–196, 1985. View at: Google Scholar
 P. S. Chua and M. R. Riggott, “The glass fibrepolymer interface: IV controlled shrinkage polymers,” Composites Science and Technology, vol. 22, no. 4, pp. 245–258, 1985. View at: Google Scholar
 R. Swan, K. L. Reifnider, K. Jayaraman, and M. ElZein, “Interface/interphase concepts in composite materials systems,” Journal of Thermoplastic Composite Materials, vol. 3, no. 1, pp. 13–23, 1990. View at: Publisher Site  Google Scholar
 E. J. Chen and J. C. Young, “The microdebonding test systems a method of quantifying adhesion in real composites,” Composites Science and Technology, vol. 42, no. 1–3, pp. 189–206, 1991. View at: Google Scholar
 L. J. Broutman, “Measurement of the fiberpolymer matrix interfacial strength,” in Interfaces in Composites, vol. 452 of ASTM Specical Technical Publication, pp. 27–41, ASTM International, 1969. View at: Google Scholar
 A. Kelly and W. R. Tyson, “Tensile properties of fibrereinforced metals: copper/tungsten and copper/molybdenum,” Journal of the Mechanics and Physics of Solids, vol. 13, no. 6, pp. 329–350, 1965. View at: Google Scholar
 D. H. Grande, Microdebonding Test for Measuring Shear Strength of the Fiber/Matrix Interface in Composite Materials [M.S. thesis], Department of Materials Science and Engineering, Massachusetts Institute of Technology, Boston, Mass, USA, 1983.
 K. Jayaraman and K. L. Reifsnider, “The interphase in unidirectional fiberreinforced epoxies: effect on residual thermal stresses,” Composites Science and Technology, vol. 47, no. 2, pp. 119–129, 1993. View at: Google Scholar
 J. K. Jim and Y. W. Mai, Engineered Interfaces in Fiber Reinforced Composites, Elsevier, London, UK, 1st edition, 1998.
 A. Dasgupta and J. S. Sirkis, “The importance of coatings to structurally embedded optical fiber sensors in smart structures,” AIAA journal, vol. 30, no. 5, pp. 1337–1343, 1992. View at: Google Scholar
 G. P. Carman and S. W. Case, “Minimizing stress concentrations in material systems with appropriate fiber coatings,” in Proceedings of the 7th Technical Conference, American Society for Composites, pp. 889–898, Technomic, 1992. View at: Google Scholar
 G. P. Carman, R. C. Averill, K. L. Reifsnider, and J. N. Reddy, “Optimization of fiber coatings to minimize stress concentrations in composite materials,” Journal of Composite Materials, vol. 27, no. 6, pp. 589–612, 1993. View at: Google Scholar
 M. S. Madhukar and L. T. Drazal, “Fibermatrix adhesion and its effect on composite mechanical properties. I. Inplane and interlaminar shear behavior of graphite/epoxy composites,” Journal of Composite Materials, vol. 25, no. 8, pp. 932–957, 1991. View at: Google Scholar
 M. S. Madhukar and L. T. Drazal, “Fibermatrix adhesion and its effect on composite mechanical properties. II. Longitudinal (0°) and transverse (90°) tensile and flexure behavior of graphite/epoxy composites,” Journal of Composite Materials, vol. 25, no. 8, pp. 958–991, 1991. View at: Google Scholar
 M. S. Madhukar and L. T. Drazal, “Fibermatrix adhesion and its effect on composite mechanical properties. III. Longitudinal (0°) compressive properties of graphite/epoxy composites,” Journal of Composite Materials, vol. 26, no. 3, pp. 310–333, 1992. View at: Google Scholar
 M. S. Madhukar and L. T. Drazal, “Fibermatrix adhesion and its effect on composite mechanical properties: IV. Mode I and mode II fracture toughness of graphite/epoxy composites,” Journal of Composite Materials, vol. 26, no. 7, pp. 936–968, 1992. View at: Google Scholar
 K. E. Remple and T. Niu, “Graphite/epoxy [±45]s tubes, their static axial & shear properties and their fatigue behavior under completely reversed loading,” Journal of Composite Materials, vol. 16, no. 3, pp. 172–187, 1982. View at: Publisher Site  Google Scholar
 C. B. Norris, “Strength of orthotropic materials subjected to combined stress US forest products,” Laboratory Report 1816, 1962. View at: Google Scholar
 N. O. Hoffman, “The brittle strength of orthotropic materials,” Journal of Composite Materials, vol. 1, no. 2, pp. 200–206, 1967. View at: Publisher Site  Google Scholar
 S. W. Tsai and E. M. Wu, “A general theory of strength for anisotropic materials,” Journal of Composite Materials, vol. 5, pp. 58–80, 1971. View at: Google Scholar
 S. W. Tsai and H. T. Hahn, Introduction to Composite Materials, Technomic, Westport, Conn, USA, 1980.
 S. C. Cowin, “On the strength anisotropy of bone and wood,” Journal of Applied Mechanics, vol. 46, no. 4, pp. 832–838, 1979. View at: Google Scholar
 A. H. Puppo and H. A. Evensen, “Strength of anisotropic materials under combined stresses,” AIAA Journal, vol. 10, no. 4, pp. 468–474, 1972. View at: Google Scholar
 V. D. Azzi and S. W. Tsai, “Anisotropic strength of composites,” Experimental Mechanics, vol. 5, no. 9, pp. 283–288, 1965. View at: Publisher Site  Google Scholar
 M. N. Abouelwafa, A. H. Hamdy, and E. A. Showaib, “A new testing machine for fatigue under combined bending and torsion acting outofphase,” Alexandria Engineering Journal, vol. 28, no. 4, pp. 113–130, 1989. View at: Google Scholar
 A. I. Sharara, Effect of Stress Ratio on Fatigue Characteristics of WovenRoving Glass Reinforced Polyester [M.S. thesis], Faculty of Engineering, Alexandria University, Alexandria, Egypt, 1997.
 M. M. Y. B. Elshabasy, The Inclusion Effect on the Fatigue Strength of Woven Roving GRP Composite Materials [M.S. thesis], Faculty of Engineering, Alexandria University, Alexandria, Egypt, 2001.
 A. J. Elkhateeb, Delamination Growth of Glass Fiber Reinforced Composite under Cyclic Torsional Moments [Ph.D. dissertation], Faculty of Engineering, Alexandria University, Alexandria, Egypt, 2001.
 A. A. ElMidany, Fatigue of WovenRoving Glass Fiber Reinforced Polyester Under Combined Bending and Torsion [Ph.D. dissertation], Faculty of Engineering, Alexandria University, Alexandria, Egypt, 1995.
 C. C. Chamis, Mechanics of Load Transfer at the Fiber/ Matrix Interface NASA TN D6588, 1972.
 “Tentative specification for woven roving glass fabric for polyester glass laminates ASTM designation,” Tech. Rep. D2150 63T, 1963. View at: Google Scholar
 W. Binienda, A. S. D. Wang, Y. Zhong, and E. S. Reddy, “A criterion for mixedmode matrix cracking in graphiteepoxy,” in Composites Composite Materials: Testing and Design, S. P. Garbo, Ed., vol. 9 of Special Technical Publication, pp. 287–300, American Society for Testing and Materials, Philadelphia, Pa, USA, 1990. View at: Google Scholar
 R. A. Hill, “Theory of the yielding and plastic flow of anisotropic metals,” Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, vol. 193, pp. 281–297, 1948. View at: Google Scholar
 C. B. Norris and P. F. McKinnon, “Supplement to compression, tension, and shear tests on yellow poplar plywood panels of sizes that do not buckle with tests made at various angles to the face grain (shear tests) US forest products,” Laboratory Report 1328, 1956. View at: Google Scholar
 E. K. Ashkenazi, “On the problem of strength for anisotropic materials,” Technical Physics, vol. 4, no. 3, pp. 333–338, 1959. View at: Google Scholar
Copyright
Copyright © 2014 Mohamed Mostafa Yousef Bassyouny Elshabasy. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.