<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">MME</journal-id><journal-title-group><journal-title>Modern Mechanical Engineering</journal-title></journal-title-group><issn pub-type="epub">2164-0165</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/mme.2017.74007</article-id><article-id pub-id-type="publisher-id">MME-79599</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Finite Element Model of a Filament-Wound Composite Tube Subjected to Uniaxial Tension
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Gyula</surname><given-names>Szabó</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Károly</surname><given-names>Váradi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Dávid</surname><given-names>Felhős</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Machine and Product Design, Budapest University of Technology and Economics, Budapest, Hungary</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>szabo.gyula@gt3.bme.hu(GS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>12</day><month>10</month><year>2017</year></pub-date><volume>07</volume><issue>04</issue><fpage>91</fpage><lpage>112</lpage><history><date date-type="received"><day>2,</day>	<month>August</month>	<year>2017</year></date><date date-type="rev-recd"><day>10,</day>	<month>October</month>	<year>2017</year>	</date><date date-type="accepted"><day>13,</day>	<month>October</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The aim of this paper is to demonstrate the mechanical behaviour of a filament-wound composite tube subjected to uniaxial tension by finite element analysis. Uniaxial tensile test experiments have been carried out on standard specimen and hose piece in order to verify finite element models and material properties and also to assess failure mode of composite plies. Composite reinforcement plies are modeled as linear orthotropic, while elastomer liners are described by hyperelastic material model. Results of finite element models and experiments show good agreement in the initial phase of uniaxial tension, which justifies utilized material models in the operating range. Results of finite element models show that transverse tension and shear load are dominant under tension. It is determined that principal failure mode of reinforcement plies is intra-ply yarn-matrix debonding caused by intensive shear of rubber matrix.
 
</p></abstract><kwd-group><kwd>Filament-Wound Composite Tube</kwd><kwd> Uniaxial Tension</kwd><kwd> Tensile Test</kwd><kwd> Debonding</kwd><kwd> Damage</kwd><kwd> Finite Element Analysis</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Composite tubes find extensive use in numerous different industrial fields; in transportation, oil industry, aerospace and aeronautical applications [<xref ref-type="bibr" rid="scirp.79599-ref1">1</xref>] . The most advantageous properties of polymeric-based tubes are high stiffness-to-weight ratio and corrosion resistance. Among manufacturing processes of composite tubes, filament winding is prominent because of high precision fibre positioning, high fibre content, low void content and good automation capability [<xref ref-type="bibr" rid="scirp.79599-ref2">2</xref>] . The most common operational loads are: internal pressure, axial and biaxial tension and buckling due to bending [<xref ref-type="bibr" rid="scirp.79599-ref3">3</xref>] . For balancing the internal forces, angles of &#177;Θ are generally adopted [<xref ref-type="bibr" rid="scirp.79599-ref4">4</xref>] . Optimal winding angles depend on the expected loading cases [<xref ref-type="bibr" rid="scirp.79599-ref5">5</xref>] . When solely internal pressure is applied the optimal winding angle is &#177;75˚, when loading is purely uniaxial winding angle needs to be as low as possible―taking into account manufacturing considerations, in case of combined internal pressure and axial loading―where ratio of circumferential to axial stress is 2:1, &#177;55˚ is adopted [<xref ref-type="bibr" rid="scirp.79599-ref6">6</xref>] , [<xref ref-type="bibr" rid="scirp.79599-ref7">7</xref>] . The lattermost is the most commonly encountered [<xref ref-type="bibr" rid="scirp.79599-ref8">8</xref>] . Rubber matrix gives the composite material favourable deformability, high flexibility and toughness.</p><p>One efficient method of obtaining material properties of a composite part is conducting uniaxial tensile test on a standard test specimen cut from a hose (in accordance with standard ASTM D638) or on a hose piece itself. Material properties in the operating range are to be determined from the initial section of the force-displacement curve characterized by constant slope. The following nonlinear section holds significance for investigating the failure procedure. [<xref ref-type="bibr" rid="scirp.79599-ref9">9</xref>] .</p><p>Bai et al. [<xref ref-type="bibr" rid="scirp.79599-ref10">10</xref>] manifested that in the initial phase of the uniaxial tension of a tensile specimen, microcracks appeared in the zone free of fibres, classified as mode I microcracks perpendicular to the loading direction. Microcracks are principally originated from porosity which can be either intra-ply or inter-ply in nature (intra-ply porosity is located inside the plies, inter-ply porosity is located in the interface layer between plies). Then, increasing tension brought about transverse cracking, which resulted in the debonding of fibre and matrix within plies.</p><p>Mechanical characterization of composite tubes is most frequently carried out by finite element models. [<xref ref-type="bibr" rid="scirp.79599-ref11">11</xref>]</p><p>The most widely used means of determining material properties of composite parts is usage of formulae of rules of mixture, which give estimation on elastic constants based on material properties of the reinforcing yarn and the matrix. In finite element models, plies are represented as laminates whose material properties are usually calculated with the use of the rules of mixtures. [<xref ref-type="bibr" rid="scirp.79599-ref12">12</xref>]</p><p>In this paper, material properties of a composite tube are elaborated with the above-mentioned approach (Chapter 2) followed by the tensile test experiment (Chapter 3) and the finite element simulation (Chapter 4) of the standard test specimen cut from the tube. Furthermore, tensile test experiment (Chapter 5) and finite element simulation (Chapter 6) of a hose piece are presented. Afterwards, comparison of results and conclusions end this article.</p><p>These analyses serve as a basis for verification of composite material properties and isotropic material properties as well, moreover, failure procedure of standard test specimen has been investigated. Stress states and strain states of finite element models corroborate assumed failure mode of composite plies. Experimental data and simulation results of standard test specimen and hose piece have been compared to one another, by which differences in mechanical behaviour of standard test specimen and hose piece have been demonstrated.</p></sec><sec id="s2"><title>2. Material Properties of Reinforcement Plies</title><p>Determination of material properties of reinforcement plies has been accomplished by the method illustrated in <xref ref-type="fig" rid="fig1">Figure 1</xref>. Firstly, photographs have been taken at the cross-section of the hose, in which yarns are respresented by their angular sections (end section). By using KLONK image measurement software [<xref ref-type="bibr" rid="scirp.79599-ref13">13</xref>] , average area of yarns has been calculated in their end sections. Average area of yarns in the end section has been utilized for the determination of fibre volume fraction. Then, average cross-sectional area of yarns has been calculated with the help of the orientation angle and average area of yarns in the end section. Afterwards, modulus of elasticity of fibre has been evaluated based on the load-strain curves of a reinforcement yarn and the average cross-sectional area of a reinforcement yarn. Two parameter Mooney-Rivlin material model has been fitted to stress-strain curves of rubber liners, therefore modulus of elasticity of matrix, being made of rubber, has been calculated. Finally, composite elastic constants have been derived based on modulus of elasticity of fibre, modulus of elasticity of matrix and fibre volume fraction.</p><p>Material properties of reinforcement plies can be determined in view of their geometric parameters and the tensile properties of the fibre and the matrix. As regards the geometric parameters, average area of yarns and fibre volume fraction are fundamental.</p><p>Computation of area has been carried out by KLONK image measurement software [<xref ref-type="bibr" rid="scirp.79599-ref13">13</xref>] .</p><p>In <xref ref-type="fig" rid="fig2">Figure 2</xref>, captured perpendicular to the axis of the hose (end section of yarns), yarns are surrounded by measurement splines. Average area of the yarns determined with the help of the software is 0.61 mm<sup>2</sup>. Outlines of the yarns might be regarded as ellipses with good approximation.</p><p>In <xref ref-type="fig" rid="fig2">Figure 2</xref>, yarns are represented by their angular cut surfaces because the included angle of the yarns and the axis of the hose is ω = 55˚. This angular cut is not representative when estimating elastic constants; therefore the cross-sec- tional area needs to be evaluated in this case (<xref ref-type="fig" rid="fig3">Figure 3</xref>). In the angular section of the yarn, the included angle of the ellipse and the normal to the cross-sectional plane is ω, hence</p><p>r = a ⋅ sin ω     while   r = b (1)</p><p>where r stands for the radius of the circular cross-section</p><p>a is the semi-major axis of the elliptical end section</p><p>b is the semi-minor axis of the elliptical end section</p><p>ω is the orientation angle</p><p>Therefore, area of the yarn in its cross-section equals:</p><p>A c = A e ⋅ sin ω (2)</p><p>where A<sub>c</sub> stands for the area of cross-section of yarn-perpendicular to its axis A<sub>e</sub> stands for the area of the angular end section of yarn-taken at the cross section of the hose.</p><p><xref ref-type="fig" rid="fig4">Figure 4</xref> shows orientation of the reinforcement plies related to the axis of the hose. The orientation angle is congruent with the included angle of the angular section of the yarn (end section) and the normal to the cross-sectional plane of the yarn.</p><sec id="s2_1"><title>2.1. Fibre Volume Fraction</title><p>Volume fraction has been calculated as the proportion of the end sectional area of yarns on the cross-section of the hose (A<sub>ae</sub>) and the cross-sectional area of the hose (A<sub>hc</sub>). The aggregate average end sectional area of yarns is (regarding 269 yarns):</p><p>A a e = 163.15   mm 2</p><p>Cross-sectional area of the hose has been calculated based on the outer (39.2 mm) and inner diameters (32.8 mm). Area of cross-section of the hose is A h c = 361.91   mm 2 . Therefore fibre volume fraction equals:</p><p>V f = A a e A h c = 45 % (3)</p></sec><sec id="s2_2"><title>2.2. Material Properties of Yarn</title><p>For obtaining material properties of the yarn, tensile test curves, provided by the manufacturer, are crucial (<xref ref-type="fig" rid="fig5">Figure 5</xref>). In this figure, two force-strain points were</p><p>selected (estimated approximate inflexion points, A and B), between which a line can be fitted for gaining characteristic stresses and therefore modulus of elasticity of yarn.</p><p>In actual practice of tensile testing, stress-strain curves can be used as a basis for acquiring modulus of elasticity of the fibre. By fitting a line to the load-strain curves and calculating stress values at these points, modulus of elasticity can be obtained by a simplified approach. With the help of the cross-sectional area of yarn ( A c = 0.496   mm 2 ) and force and strain values of point A and B ( F A = 78.57   N , F B = 221.4   N and ε A = 0.0654 , ε B = 0.163 ), the calculated stress values are: σ A = 158.1   MPa , σ B = 446   MPa ; the modulus of elasticity of the fibre is:</p><p>E f = σ B − σ A ε B − ε A = 2961   MPa (4)</p><p>σ B = 446   MPa is considered as the tensile strength of the yarn and ε B = 0.16 as the elongation at break.</p></sec><sec id="s2_3"><title>2.3. Material Properties of Rubber Liners</title><p>Material properties of rubber liners have been determined based on the curve fitted [<xref ref-type="bibr" rid="scirp.79599-ref14">14</xref>] to the uniaxial stress-strain curves provided by the manufacturer. The mechanical model of the rubber liners is 2 parameter Mooney-Rivlin hyperelastic model (<xref ref-type="fig" rid="fig6">Figure 6</xref>).</p><p>Mooney-Rivlin-parameters of rubber liners are:</p><p>C 10 = − 0.4982   MPa , C 01 = 1.523   MPa</p><p>D = 0 [1/MPa], so rubber liners are assumed to be incompressible in the FE model.</p><p>Equivalent shear modulus, Poisson’s ratio and modulus of elasticity of rubber liners are:</p><p>G r = 2 ⋅ ( C 10 + C 01 ) = 2.05   MPa (5)</p><p>ν r = 0.5</p><p>E r = 2 ⋅ ( 1 + ν r ) ⋅ G r = 6.15   MPa (6)</p></sec><sec id="s2_4"><title>2.4. Deriving Composite Elastic Constants</title><p>In the most general case, stress and strain states of composites can be described by an anisotropic material model. However, this requires 21 independent elastic constants the measurement of which is often too hard to achieve. Therefore, in most cases, the highest order of symmetry possible is taken into consideration. Firstly, composite parts are generally considered to be described as orthotropic having three mutually perpendicular planes of symmetry; filament-wound composite hoses, with their fibres aligned uniaxially, are usually regarded as transversely isotropic because the plane perpendicular to the fibre direction is a plane of isotropy ( E 2 = E 3 , G 12 = G 13 , υ 12 = υ 13 ). Transversely isotropic materials have five independent elastic constants (E<sub>1</sub>, E<sub>2</sub>, G<sub>12</sub>, G<sub>23</sub>, υ<sub>12</sub>) [<xref ref-type="bibr" rid="scirp.79599-ref15">15</xref>] .</p><p>Modulus of elasticity of the matrix is E m = E r = 6.15   MPa , its Poisson’s ratio is supposed to be υ m = υ r = 0.5 . Poisson’s ratio of the fibre is assumed to be υ f = 0.2 . Fibre volume fraction is V f = 45 % .</p><p>Matrix volume fraction and shear modulus of fibre are respectively as follows:</p><p>V m = 1 − V f = 1 − 0.45 = 0.55 (7)</p><p>G f = E f 2 ⋅ ( 1 + ν f ) = 1234   MPa (8)</p><p>Modulus of elasticity of one ply in direction 1 is:</p><p>E 1 = E f ⋅ V f + E m ⋅ V m = 1338   MPa (9)</p><p>Applying the formulae presented in [<xref ref-type="bibr" rid="scirp.79599-ref15">15</xref>] and [<xref ref-type="bibr" rid="scirp.79599-ref16">16</xref>] , the modulus of elasticity in directions 2 and 3 and shear moduli are:</p><p>E 2 = E 3 = E m 1 − V f ⋅ ( 1 − E m / E f ) = 19   MPa (10)</p><p>G 12 = G 13 = G 23 = G m 1 − V f ⋅ ( 1 − G m / G f ) = 6   MPa (11)</p><p>Since G f 12 = G f 23 , apparently G 12 = G 23 .</p><p>Poisson’s ratios of one ply are [<xref ref-type="bibr" rid="scirp.79599-ref12">12</xref>] :</p><p>υ 12 = υ 13 = V f ⋅ υ f + V m ⋅ υ m = 0.37 (12)</p><p>Contraction in plane 23 is principally governed by the mechanical behaviour of rubber, so</p><p>υ 23 = E 22 2 ⋅ G 23 − 1 = 0.498 (13)</p></sec></sec><sec id="s3"><title>3. Tensile Experiment on Standard Test Specimen Cut from the Hose</title><p>Test specimen used for uniaxial tension test has been cut from the hose in thickness direction in accordance with standard ASTM D638. The specimen, having a length of 75 mm, has been tested to complete failure. The experiment began with a pre-load of 1 N with a pre-load speed of 2 mm/min. The test was executed on a Zwick Z250 tensile testing machine with a tensile speed of 2 mm/min. The measuring limit of the load cell was 10 kN.</p><p>Firstly, fractures appeared in reinforcement plies, which led to yarn-matrix debonding. This may be attributed to yarns aligning themselves to the direction of load generating shear inside plies (<xref ref-type="fig" rid="fig7">Figure 7</xref>). As the tension progressed, gaps appeared between the two rubber liners.</p></sec><sec id="s4"><title>4. Finite Element Analysis of the Uniaxial Tension Test on Test Specimen</title><p>Finite element analysis of tensile test on standard test specimen consists of two separate sections, the first being the representation of fastening the specimen, while the second part simulates uniaxial tension.</p><sec id="s4_1"><title>4.1. Geometric Model</title><p>Inner diameter of the specimen is 28 mm, outer diameter is 44 mm, and thickness of the rubber liners equals 2.4 mm. There are 4 reinforcement plies consisting of natural yarns and rubber matrix. The material coordinate systems of the plies are cylindrical, the layup is [+55/−55/+55/−55] [˚] with respect to the axis of the hose based on technical specification provided by the manufacturer (<xref ref-type="fig" rid="fig8">Figure 8</xref>).</p><p>Geometry of the specimen can be observed in <xref ref-type="fig" rid="fig9">Figure 9</xref> along with the dimensions and the fibre orientation of ply 1 (its orientation angle is ω = +55˚).</p></sec><sec id="s4_2"><title>4.2. Connections</title><p>There are bonded contacts between the reinforcement plies and the outer rubber liner and the reinforcement plies and the inner rubber liner respectively because rubber, being the material of the matrix and also material of inner and outer liners, is vulcanized around yarns. The contacts defined between the rubber liners and the tensile jaws are frictional with a coefficient of friction &#181; = 0.8 based on [<xref ref-type="bibr" rid="scirp.79599-ref17">17</xref>] , and formulation “Augmented Lagrange”.</p></sec><sec id="s4_3"><title>4.3. Mesh</title><p>Finite element mesh (<xref ref-type="fig" rid="fig1">Figure 1</xref>0) consists of approximately 14,800 nodes and 13,100 elements.</p></sec><sec id="s4_4"><title>4.4. Loads, Boundary Conditions</title><p>The simulation consists of three time steps. In the first and the second one, test specimen is installed between the tensile jaws, in the third time step, tension is applied to the specimen.</p><p>In the first time step, the upper tensile jaw at the right end descends −2 mm in Y direction, while the lower tensile jaw opposite of the former stands still (<xref ref-type="fig" rid="fig1">Figure 1</xref>1).</p><p>In the second time step, left end of the specimen is being fixed. The upper tensile jaw sinks −2 mm in Y direction, while the lower tensile jaw is fixed.</p><p>Meanwhile, the upper right tensile jaw holds its position −2 mm from its initial horizontal position and the lower right tensile jaw is still fixed. In the third time step, tensile jaws situated at the right end are subjected to a prescribed displacement of 25 mm in positive X direction. In the meantime, tensile jaws situated at the left end hold their position.</p></sec><sec id="s4_5"><title>4.5. Material Properties</title><p>Material model of the reinforcement plies is linear orthotropic, the corresponding material properties can be found in Chapter 2.4. The rubber liners surrounding the reinforcement plies are described by a 2 parameter Mooney-Rivlin model, its parameters defined in Chapter 2.3.</p></sec><sec id="s4_6"><title>4.6. Results of the Simulation of Tensile Test on Tensile Specimen</title><p>Results of finite element analysis of uniaxial tension of tensile specimen at a displacement of 25 mm are disclosed below.</p><p>Elongation in the global coordinate system at the end of the third time step is significant in the whole cross-section, approximately ε<sub>x</sub> = 0.5 (<xref ref-type="fig" rid="fig1">Figure 1</xref>2).</p><p>Stress and strain results, being true stress and true strain results, regarding ply 1 are hereinafter presented, afterwards, mean values are disclosed for each ply in tables. Stress and strain results in ply 1 are to be considered based on the material orientation of <xref ref-type="fig" rid="fig9">Figure 9</xref>.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>3 shows normal strain distribution in material direction 1 (fibre direction) in fibre coordinate system for ply 1 at a prescribed displacement of 25 mm. The strain distribution tends to be uniform between the necks, while near the necks, it represents the material orientation of the ply. Strain values are low because of the high longitudinal modulus of elasticity, E<sub>1</sub> = 1338 MPa.</p><p>As regards other components of strain in the plane of ply1, normal strain in material direction 2 and shear strain in plane 12 have the same distribution as presented in <xref ref-type="fig" rid="fig1">Figure 1</xref>4.</p><p>The predominant strain component is shear strain (<xref ref-type="fig" rid="fig1">Figure 1</xref>4), whose absolute value is approximately the same for every ply, however sign of shear strains is opposite in adjacent plies, which can be attributed to the balanced layup (<xref ref-type="table" rid="table1">Table 1</xref>). The reason why the shear stress is so significant is that load is off-axis regarding the material orientation in each ply. The significant shear stress results in high shear strain considering fairly low shear moduli.</p><p>Average strains with the exception of longitudinal strains in the fibre coordinate system are considerable; <xref ref-type="table" rid="table1">Table 1</xref> indicates that global longitudinal strain is transformed primarily to transverse normal strain and shear strain in the fibre coordinate system.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>5 shows distribution of longitudinal stresses in the global coordinate system. Highest stress values can be recognized in reinforcement layers, while nearly one third of the peak value is present in rubber liners, so reinforcement plies are considered to be the load-bearing components.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Average strains in reinforcement plies in the vertical midplane</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >ply no.</th><th align="center" valign="middle" >ε<sub>x</sub> [−]</th><th align="center" valign="middle" >ε<sub>1</sub> [−]</th><th align="center" valign="middle" >ε<sub>2</sub> [−]</th><th align="center" valign="middle" >γ<sub>12</sub> [−]</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.4665</td><td align="center" valign="middle" >0.001</td><td align="center" valign="middle" >0.2493</td><td align="center" valign="middle" >−0.6349</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.4663</td><td align="center" valign="middle" >0.0015</td><td align="center" valign="middle" >0.2385</td><td align="center" valign="middle" >0.6504</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.4663</td><td align="center" valign="middle" >0.0016</td><td align="center" valign="middle" >0.2386</td><td align="center" valign="middle" >−0.6504</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.4664</td><td align="center" valign="middle" >0.0009</td><td align="center" valign="middle" >0.2488</td><td align="center" valign="middle" >0.6359</td></tr></tbody></table></table-wrap><p>Normal stress distribution in material direction 1 in fibre coordinate system is identical to normal strain distribution in the same direction; normal stress values (<xref ref-type="fig" rid="fig1">Figure 1</xref>6) are considerable. Shear stress (<xref ref-type="fig" rid="fig1">Figure 1</xref>7) tends to be the predominant stress component along with transverse normal stress (σ<sub>2</sub>) having the same distribution.</p><p>Average stress values can be observed in <xref ref-type="table" rid="table2">Table 2</xref>. The maximum values can be found in ply 2 and ply 3 for every stress component; however the stress variance in the direction normal to each ply tends to be slight. Sign of shear strains is opposite in adjacent plies, because the layup is balanced, similarly as in case of strains.</p></sec><sec id="s4_7"><title>4.7. Reviewing the Results</title><p>At the end of tension, longitudinal stress in the reinforcement plies at a displacement of 25 mm is approximately 8 MPa (<xref ref-type="fig" rid="fig1">Figure 1</xref>8) in the global coordinate system. Aggregate cross-sectional area of the reinforcement plies is 9 mm<sup>2</sup> (width: 3.1 mm, height: 2.9 mm). The respective average longitudinal stress value for the rubber liners is 1.9 MPa. Total cross-sectional area of rubber liners equals 11.8 mm<sup>2</sup> (width: 3.1 mm, height: 1.9 mm). Hence, the force exerted to the cross section is 94.3 N. Reaction force queried at this load step is 93.57 N.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>9 illustrates experimental data and simulation results of tensile test as force-displacement curves. There is good agreement between finite element re-</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Average stresses in reinforcement plies in the vertical midplane</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >ply no.</th><th align="center" valign="middle" >σ<sub>x</sub> [MPa]</th><th align="center" valign="middle" >σ<sub>1</sub> [MPa]</th><th align="center" valign="middle" >σ<sub>2</sub> [MPa]</th><th align="center" valign="middle" >τ<sub>12</sub> [MPa]</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >7.88</td><td align="center" valign="middle" >2.95</td><td align="center" valign="middle" >4.67</td><td align="center" valign="middle" >−4.09</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >8.05</td><td align="center" valign="middle" >3.68</td><td align="center" valign="middle" >4.46</td><td align="center" valign="middle" >4.1</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >8.06</td><td align="center" valign="middle" >3.69</td><td align="center" valign="middle" >4.46</td><td align="center" valign="middle" >−4.1</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >7.86</td><td align="center" valign="middle" >2.91</td><td align="center" valign="middle" >4.66</td><td align="center" valign="middle" >4</td></tr></tbody></table></table-wrap><p>sults and experimental data in the initial section of the force-displacement curves regarding the operating range. However, after a displacement of 10 mm, in the course of the failure procedure, substantial difference can be noticed between experimental and simulation results because failure mechanisms are not taken into account in the FE model.</p></sec></sec><sec id="s5"><title>5. Tensile Test Experiment of Hose Piece</title><p>Tensile test of a hose piece having length of 110 mm was performed (<xref ref-type="fig" rid="fig2">Figure 2</xref>0) to a displacement of 25 mm. Cross-section of the hose piece is in accordance with <xref ref-type="fig" rid="fig8">Figure 8</xref>. Uniaxial tension was realized with the help of two steel plugs fixed into the hose by pipe clamps and by bonding them to the inner liner of the hose.</p><p>Hose piece was pulled with a speed of 2 mm/min and a preload of 1 N. Tensile tests were performed three times.</p></sec><sec id="s6"><title>6. Finite Element Model of the Uniaxial Tension of Hose Piece</title><sec id="s6_1"><title>6.1. Connections</title><p>There are bonded contacts between the reinforcement plies and the outer rubber liner and the reinforcement plies and the inner rubber liner respectively. The plugs are connected to the inner lateral surface of the hose by bonded contacts.</p></sec><sec id="s6_2"><title>6.2. Mesh</title><p>Mesh of the finite element hose piece model, being a three-dimensional full model, consists of nearly 25,100 nodes and 21,000 elements (longitudinal section illustrated in <xref ref-type="fig" rid="fig2">Figure 2</xref>1). Reinforcement plies are represented by 4 adjacent layers of elements between 3 layers representing outer rubber liner and 3 layers representing inner rubber liner.</p></sec><sec id="s6_3"><title>6.3. Loads, Boundary Conditions</title><p>In the first time step, pressure of 1.5 MPa, which represents pipe clamps, is applied to sections of the outer side surface of the hose shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>2.</p><p>In the second time step, a prescribed displacement of 25 mm, being in line with the prescribed boundary condition of the test, is applied to the right end of the right plug. Meanwhile, the left end of the left plug is fixed.</p></sec><sec id="s6_4"><title>6.4. Material Properties</title><p>Material properties of the reinforcement plies, presented in detail in Chapter 2.4, and rubber liners, presented in detail in Chapter 2.3, are identical to material properties of the test specimen.</p><p>Steel plugs utilized for transmitting tensile load to the hose piece are modeled as rigid bodies in the present finite element model.</p></sec></sec><sec id="s7"><title>7. Results of FE Analysis and Experiment of Tensile Test of Hose Piece</title><p><xref ref-type="fig" rid="fig2">Figure 2</xref>3 shows experimental data of uniaxial tension of hose piece with the mean values along with simulation results.</p><p>Simulation results and experimental data are in good agreement, difference between simulation results and mean experimental data is less than 10% of mean experimental data, moreover, simulation results are located between minimum and maximum of experimental data at all data points (<xref ref-type="fig" rid="fig2">Figure 2</xref>3).</p><p>Deformation, strain and stress results at a displacement of 25 mm are displayed below.</p><p>Hose piece elongates as much as 0.33 (<xref ref-type="fig" rid="fig2">Figure 2</xref>4) at a displacement of 25 mm, which is considered significant. The deformed shape of <xref ref-type="fig" rid="fig2">Figure 2</xref>4 is noticeably similar to the deformed shape after tensile testing, shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>0.</p><p>Normal strain in material direction 1 in fibre coordinate system can be observed in <xref ref-type="fig" rid="fig2">Figure 2</xref>5. These strain values tend to be of low order of magnitude due to high stiffness in material direction 1.</p><p>Shear strain in plane 12 shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>6 has a predominant role in the structural behaviour of hose piece under tension because of the load being off-axis to the fibre direction and shear moduli being relatively low. Normal strain in material direction 2, having the same distribution as shear strain in plane 12, is also considerably significant (<xref ref-type="table" rid="table3">Table 3</xref>) being the other strain component into which global longitudinal strain is transformed. Strains are almost equally distributed along reinforcement plies.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Average strains in reinforcement plies in the vertical midplane</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >ply no.</th><th align="center" valign="middle" >ε<sub>x</sub> [−]</th><th align="center" valign="middle" >ε<sub>1</sub> [−]</th><th align="center" valign="middle" >ε<sub>2</sub> [−]</th><th align="center" valign="middle" >γ<sub>12</sub> [−]</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.3276</td><td align="center" valign="middle" >0.003</td><td align="center" valign="middle" >0.1694</td><td align="center" valign="middle" >−0.4531</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0.3279</td><td align="center" valign="middle" >0.0024</td><td align="center" valign="middle" >0.1675</td><td align="center" valign="middle" >0.4571</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0.3285</td><td align="center" valign="middle" >0.0006</td><td align="center" valign="middle" >0.1679</td><td align="center" valign="middle" >−0.4597</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0.3287</td><td align="center" valign="middle" >0.0018</td><td align="center" valign="middle" >0.1658</td><td align="center" valign="middle" >0.4571</td></tr></tbody></table></table-wrap><p><xref ref-type="fig" rid="fig2">Figure 2</xref>7 also proves that reinforcement plies are mainly responsible for the load bearing capability of the hose due to high stress values present in them. Local maximum stresses in the hose piece can be found near the contact to the plugs.</p><p>Stress distribution is very much alike strain distribution in reinforcement plies, so stresses are presented as average values in <xref ref-type="table" rid="table4">Table 4</xref>.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>8 illustrates experimental data and simulation results for the hose piece along with experimental data and simulation results for the standard test specimen proportioned to the hose piece. The scaling factor utilized for this is as follows.</p><p>f s c = P w = 28.27 (14)</p><p>where f<sub>sc</sub> is the scaling factor</p><p>P stands for the perimeter of the mid-diameter of hose piece</p><p>w stands for the width of standard test specimen</p><p>Experimental data and simulation results for hose piece and standard test specimen are all in good agreement with one another to a displacement of 10 mm. Afterwards, the difference between experimental data for the standard test specimen and experimental data for hose piece along with simulation results rises gradually. This increase of difference in reaction forces may be attributed to the difference in mechanical behaviour regarding standard test specimen and hose piece which is presented in Chapter 8 in detail.</p><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Average stresses in reinforcement plies in the vertical midplane</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >ply no.</th><th align="center" valign="middle" >σ<sub>x</sub> [MPa]</th><th align="center" valign="middle" >σ<sub>1</sub> [MPa]</th><th align="center" valign="middle" >σ<sub>2</sub> [MPa]</th><th align="center" valign="middle" >τ<sub>12</sub> [MPa]</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >6.63</td><td align="center" valign="middle" >5.44</td><td align="center" valign="middle" >2.94</td><td align="center" valign="middle" >−2.66</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >6.15</td><td align="center" valign="middle" >5.23</td><td align="center" valign="middle" >2.85</td><td align="center" valign="middle" >2.69</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >4.31</td><td align="center" valign="middle" >0.37</td><td align="center" valign="middle" >2.83</td><td align="center" valign="middle" >−2.7</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >3.63</td><td align="center" valign="middle" >2.58</td><td align="center" valign="middle" >2.84</td><td align="center" valign="middle" >2.73</td></tr></tbody></table></table-wrap></sec><sec id="s8"><title>8. Conclusions</title><p>Simulation results and experimental data of standard test specimen are in good agreement in the displacement interval of 0 - 10 mm. Relatively high shear strain, γ<sub>12</sub> = 0.26 and significant normal strain ε<sub>2</sub> = 0.11 belongs to the end of this displacement interval.</p><p>Tension exerted on the standard test specimen causes shear and transverse tension in reinforcement plies because load is aligned off-axis to the yarns. Load is distributed among yarns by the rubber matrix, having low modulus of elasticity, being susceptive to shear. In the course of uniaxial tensile loading, yarns and matrix debond gradually within the plies (<xref ref-type="fig" rid="fig7">Figure 7</xref>). This procedure, taking place between displacements of 15 and 25 mm is nonlinear (<xref ref-type="fig" rid="fig1">Figure 1</xref>9, <xref ref-type="fig" rid="fig2">Figure 2</xref>8). Yarns in the standard test specimen are regarded as short yarns because the thickness of the specimen is small, which makes failure procedure initiate at low strain values.</p><p>However, in the case of the hose piece, reinforcement plies being intact and continuous, in which yarns are much longer than in standard test specimen, are less susceptible to intra-ply yarn-matrix debonding. Continuous plies make the hose piece stiffer in radial direction (<xref ref-type="fig" rid="fig2">Figure 2</xref>0, <xref ref-type="fig" rid="fig2">Figure 2</xref>4), which raises resistance against tension, thus generating higher forces at the same displacement values (<xref ref-type="fig" rid="fig2">Figure 2</xref>8). Tension causes less transverse normal stress and shear stress in hose piece compared to the case of standard test specimen (<xref ref-type="table" rid="table2">Table 2</xref> and <xref ref-type="table" rid="table4">Table 4</xref> respectively).</p><p>With respect to the hose piece, the difference between simulation results and experimental data regarding force-displacement values is less than 10% of the experimental data (<xref ref-type="fig" rid="fig2">Figure 2</xref>3). Finite element results of standard test specimen and finite element results of hose piece are in good agreement in the initial phase of uniaxial tension. Experimental results and FE simulation results, regarding standard test specimen and hose piece, also show good agreement in this section of uniaxial tension.</p></sec><sec id="s9"><title>Acknowledgements</title><p>Authors would like to express their gratitude to Department of Polymer Engineering, Budapest University of Technology and Economics namely to G&#225;bor Szeb&#233;nyi and Bal&#225;zs Pinke for providing them with the testing environment, moreover to Gy&#246;rgy Szab&#243; at the Department of Machine and Product Design, Budapest University of Technology and Economics for test specimen preparation and to eCon Engineering Ltd., namely to D&#225;niel Farkas for valuable recommendations regarding the scope of this study.</p><p>The recent study and publication was realized within the Knorr-Bremse Scholarship Program supported by the Knorr-Bremse Rail Systems Budapest.</p></sec><sec id="s10"><title>Cite this paper</title><p>Szab&#243;, G., V&#225;radi, K. and Felhős, D. (2017) Finite Element Model of a Filament-Wound Composite Tube Subjected to Uniaxial Tension. Modern Mechanical Engineering, 7, 91-112. https://doi.org/10.4236/mme.2017.74007</p></sec></body><back><ref-list><title>References</title><ref id="scirp.79599-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Braiek, S., Zitoune, R., Ben Khalifa, A. and Zidi, M. (2017) Experimental and Numerical Study of Adhesively Bonded ±55° Filament Wound Tubular Specimens under Uniaxial Tensile Loading. Composite Structures, 172, 297-310.  
https://doi.org/10.1016/j.compstruct.2017.03.103</mixed-citation></ref><ref id="scirp.79599-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Almeida Jr., J.H.S., Ribeiro, M., Volnei, T. and Amico, S.C. (2017) Damage Modeling for Carbon Fiber/Epoxy Filament Wound Composite Tubes under Radial Compression. Composite Structures, 160, 204-210.  
https://doi.org/10.1016/j.compstruct.2016.10.036</mixed-citation></ref><ref id="scirp.79599-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Carroll, M., Ellyin, F., Kujawski, D. and Chiu, A.S. (1995) The Rate-Dependent Behaviour of ±55° Filament-Wound Glass-Fibre/Epoxy Tubes under Biaxial Loading. Composites Science and Technology, 55, 391-403.  
https://doi.org/10.1016/0266-3538(95)00119-0</mixed-citation></ref><ref id="scirp.79599-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Hinton, M.J., Soden, P.D. and Kaddour, A.S. (1995) Strength of Composite Laminates under Biaxial Loads, Proceedings of 10th Int. Conf. on Composite Materials (ICCM-10), IV, Whistler, British Columbia, Canada, 65-72.</mixed-citation></ref><ref id="scirp.79599-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Soden, P.D., Kitching, R., Tse, P.C. and Tsavalas, Y (1993) Influence of Winding angle on the Strength and Deformation of Filament-Wound Composite Tubes Subjected to Uniaxial and Biaxial Loads. Composites Science and Technology, 46, 363-378. https://doi.org/10.1016/0266-3538(93)90182-G</mixed-citation></ref><ref id="scirp.79599-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, L. and Li, H. (2011) Investigation on Mechanical Behavior of Filament-Wound CFRP Tubes. Proceedings of SPIE, 7978, 79781Y-6.  
https://doi.org/10.1117/12.880372</mixed-citation></ref><ref id="scirp.79599-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Soden, P.D., Kitching, R. and Tse, P.C. (1989.) Experimental Failure Stresses for ±55° Filament Wound Glass Fibre Reinforced Plastic Tubes under Biaxial Loads. Composites, 20, 125-135. https://doi.org/10.1016/0010-4361(89)90640-X</mixed-citation></ref><ref id="scirp.79599-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Hull, D., Legg, M.J. and Spencer, B (1978) Failure of Glass/Polyester Filament Wound Pipe. Composites, 9, 17-24. https://doi.org/10.1016/0010-4361(78)90513-X</mixed-citation></ref><ref id="scirp.79599-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Wang, Z., Liu, W. and Du, Y. (2012) Experimental Study of the Mechanical Properties of the Composite Materials Hose. Advanced Materials Research, 496, 3-6.  
https://doi.org/10.4028/www.scientific.net/AMR.496.3</mixed-citation></ref><ref id="scirp.79599-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Bai, J., Seeleuthner, P. and Bompard, P. (1997) Mechanical Behaviour of ± 55° Filament-Wound Glass-Fibre/Epoxy-Resin Tubes: I. Microstructural Analyses, Mechanical Behaviour and Damage Mechanisms of Composite Tubes under Pure Tensile Loading, Pure Internal Pressure, and Combined Loading. Composites Science and Technology, 57, 141-153. https://doi.org/10.1016/S0266-3538(96)00124-8</mixed-citation></ref><ref id="scirp.79599-ref11"><label>11</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Manach</surname><given-names> P.Y. </given-names></name>,<etal>et al</etal>. (<year>2008</year>)<article-title>Experimental Study of the Mechanical Properties of the Composite Materials Hose</article-title><source> Applied and Computational Mechanics</source><volume> 2</volume>,<fpage> 291</fpage>-<lpage>302</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.79599-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Kollár, L.P. and Springer, G.S. (2003) Mechanics of Composite Structures. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511547140</mixed-citation></ref><ref id="scirp.79599-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">KLONK User’s Manual. http://www.imagemeasurement.com/en/manual</mixed-citation></ref><ref id="scirp.79599-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Ansys Mechanical User’s Guide, Engineering Data User’s Guide, Material Data, Mechanical Material Curve Fitting.</mixed-citation></ref><ref id="scirp.79599-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Chawla, K.K. (2009) Composite Materials Science and Engineering. 3rd Edition, Springer, New York, London.</mixed-citation></ref><ref id="scirp.79599-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Chamis, C.C. (1983) NASA Tech. Memo. 83320. 38th Annual Conference of the Society of Plastics Industry, Houston.</mixed-citation></ref><ref id="scirp.79599-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Karger-Kocsis, J., Mousa, A., Major, Z. and Békési, N. (2008) Dry Friction and Sliding Wear of EPDM Rubbers against Steel as a Function of Carbon Black Content. Wear, 264, 359-367.</mixed-citation></ref></ref-list></back></article>