m (Onate moved page Zarate et al 2023a to Zarate et al 2023b) |
|||
(49 intermediate revisions by 2 users not shown) | |||
Line 26: | Line 26: | ||
david.cendon.franco@upm.es''' | david.cendon.franco@upm.es''' | ||
--> | --> | ||
− | |||
==Abstract== | ==Abstract== | ||
Line 33: | Line 32: | ||
'''Keywords:''' 3D solid finite elements, rotation-free beam elements, polymer fiber rebars | '''Keywords:''' 3D solid finite elements, rotation-free beam elements, polymer fiber rebars | ||
− | ==1 Introduction== | + | ==1. Introduction== |
There is an increasing interest in the use of non-metallic materials, such as polymers, in the building and construction (B&C) sector as an alternative to standard materials, such as steel and cement <span id='citeF-1'></span><span id='citeF-2'></span>[[#cite-1|[1,2]]]. These new materials have many advantages, compared to traditional ones, in terms of strength-weight and stiffness-weigh ratios, corrosion resistance to chemical agents, and fatigue performance, to name a few. All these advantages can extend the life-frame of structures and minimize the maintenance required, thus contributing to improve the sustainability of the sector. Current research efforts aim to develop resistant, sustainable and cost-effective materials incorporating polymers for a number of applications in building and construction. | There is an increasing interest in the use of non-metallic materials, such as polymers, in the building and construction (B&C) sector as an alternative to standard materials, such as steel and cement <span id='citeF-1'></span><span id='citeF-2'></span>[[#cite-1|[1,2]]]. These new materials have many advantages, compared to traditional ones, in terms of strength-weight and stiffness-weigh ratios, corrosion resistance to chemical agents, and fatigue performance, to name a few. All these advantages can extend the life-frame of structures and minimize the maintenance required, thus contributing to improve the sustainability of the sector. Current research efforts aim to develop resistant, sustainable and cost-effective materials incorporating polymers for a number of applications in building and construction. | ||
− | Polymers and composite materials can be used in B&C in many forms: as external reinforcements of concrete beams and slabs <span id='citeF-3'></span><span id='citeF-4'></span>[[#cite-3|[3,4]]], as short fibers embedded in concrete <span id='citeF-5'></span>[[#cite-5|[5]]], or as pultruded elements which can be used as rebars for concrete reinforcement <span id='citeF-6'></span>[[#cite-6|[6]]] or as beams elements to construct different structures such as bridges <span id='citeF-7'></span>[[#cite-7|[7]]], or lighthouses <span id='citeF-8'></span>[[#cite-8|[8]]] ( | + | Polymers and composite materials can be used in B&C in many forms: as external reinforcements of concrete beams and slabs <span id='citeF-3'></span><span id='citeF-4'></span>[[#cite-3|[3,4]]], as short fibers embedded in concrete <span id='citeF-5'></span>[[#cite-5|[5]]], or as pultruded elements which can be used as rebars for concrete reinforcement <span id='citeF-6'></span>[[#cite-6|[6]]] or as beams elements to construct different structures such as bridges <span id='citeF-7'></span>[[#cite-7|[7]]], or lighthouses <span id='citeF-8'></span>[[#cite-8|[8]]] ([[#img-1|Figures 1]]a and [[#img-1|1]]b), to name some examples. |
<div id='img-1a'></div> | <div id='img-1a'></div> | ||
<div id='img-1b'></div> | <div id='img-1b'></div> | ||
<div id='img-1'></div> | <div id='img-1'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |align="center" | | ||
+ | {|style="margin: 1em auto 0.1em auto;width:auto;" | ||
+ | |+ | ||
|- | |- | ||
− | |[[Image:Draft_Onate_870101189-Figure1-intro-a.png|393px|]] | + | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure1-intro-a.png|393px|]] |
− | |[[Image:Draft_Onate_870101189-Figure1-intro-b.png|121px|]] | + | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure1-intro-b.png|121px|]] |
− | |- style="text-align: center; font-size: 75%;" | + | |- |
− | | (a) | + | |style="text-align: center;font-size: 75%;padding-bottom:10px;"|(a) |
− | | | + | |style="text-align: center;font-size: 75%;padding-bottom:10px;"|(b) |
− | + | ||
− | + | ||
|} | |} | ||
+ | |- | ||
+ | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 1'''. Applications of pultruded beams in construction <span id='citeF-7'></span><span id='citeF-8'></span>[[#cite-7|[7,8]]] | ||
+ | |} | ||
+ | |||
The use of pultruded FRP as reinforcement in concrete has several advantages compared to steel reinforcement such as corrosion resistance, high longitudinal (unidirectional) strength (2 or 3 times higher than steel for equivalent bar diameter), high fatigue endurance (which depends on the type of reinforcing fiber and bar), magnetic transparency, lightweight (25% of steel weight) and low thermal and electric conductivity (for glass and aramid fibers) <span id='citeF-9'></span>[[#cite-9|[9]]]. Glass fibers are the common choice due to their good performance and reduced cost, compared to other synthetic fibers that have improved properties but are more expensive, such as carbon. Research efforts are currently carried out for developing new families of rebars with synthetic fibers, as reinforcement for concrete structures <span id='citeF-10'></span><span id='citeF-11'></span>[[#cite-10|[10,11]]]. This type of synthetic fiber rebars could be an alternative to current FRP rebars using glass fibers, and even to steel rebars for specific applications. | The use of pultruded FRP as reinforcement in concrete has several advantages compared to steel reinforcement such as corrosion resistance, high longitudinal (unidirectional) strength (2 or 3 times higher than steel for equivalent bar diameter), high fatigue endurance (which depends on the type of reinforcing fiber and bar), magnetic transparency, lightweight (25% of steel weight) and low thermal and electric conductivity (for glass and aramid fibers) <span id='citeF-9'></span>[[#cite-9|[9]]]. Glass fibers are the common choice due to their good performance and reduced cost, compared to other synthetic fibers that have improved properties but are more expensive, such as carbon. Research efforts are currently carried out for developing new families of rebars with synthetic fibers, as reinforcement for concrete structures <span id='citeF-10'></span><span id='citeF-11'></span>[[#cite-10|[10,11]]]. This type of synthetic fiber rebars could be an alternative to current FRP rebars using glass fibers, and even to steel rebars for specific applications. | ||
Line 63: | Line 68: | ||
The content of the paper is as follows. In the next section we present the overall finite element approach for non-linear analysis of FRP rebars combining solid elements and rotation-free beam elements. Then we describe the 3D solid finite element formulation for modelling the behavior of the matrix material. This is followed by the finite element formulation for the fibers using rotation-free beam elements based on a cell-vertex approach. The coupling of 3D-solid elements and rotation-free beam elements is detailed. The finite element formulation is validated in the non-linear analysis up to failure of FRP bar specimens formed by an epoxy matrix and carbon fibers tested under axial, bending and shear loads. Numerical results are compared with experimental results for the same tests. Excellent agreement is obtained in the comparison. | The content of the paper is as follows. In the next section we present the overall finite element approach for non-linear analysis of FRP rebars combining solid elements and rotation-free beam elements. Then we describe the 3D solid finite element formulation for modelling the behavior of the matrix material. This is followed by the finite element formulation for the fibers using rotation-free beam elements based on a cell-vertex approach. The coupling of 3D-solid elements and rotation-free beam elements is detailed. The finite element formulation is validated in the non-linear analysis up to failure of FRP bar specimens formed by an epoxy matrix and carbon fibers tested under axial, bending and shear loads. Numerical results are compared with experimental results for the same tests. Excellent agreement is obtained in the comparison. | ||
− | ==2 Combination of 3D solid elements and rotation-free elements for analysis of polymer fiber rebars== | + | ==2. Combination of 3D solid elements and rotation-free elements for analysis of polymer fiber rebars== |
− | + | [[#img-2|Figure 2]] depicts the modelling strategy for analysis of FRP rebars with the FEM <span id='citeF-14'></span><span id='citeF-15'></span>[[#cite-14|[14,15]]]. The procedure proposed dissociates the longitudinal contribution of fibers to the composite strength and stiffness, from their contribution in the transversal and shear directions. This is done by modeling a transverse reinforced resin, that accounts for the resin and the transverse and shear performance of the fibers, with standard 3D solid elements; and by modeling the longitudinal contribution of fibers to the composite with a collection of beams embedded in the 3D solid model. | |
<div id='img-2'></div> | <div id='img-2'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure1.png|600px|Finite element modelling of FRP rebar]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 2'''. Finite element modelling of FRP rebar |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
Both the 3D solid elements and the beam elements are formulated using an updated Lagrangian description in order to account for geometrically non-linear effects in the deformation of the FRP bars. Following standard finite element procedures, the discretized equilibrium equation for the bar structure, expressed in the current configuration at time <math display="inline">t</math>, can be written in the following quasi-static form <span id='citeF-15'></span>[[#cite-15|[15]]] | Both the 3D solid elements and the beam elements are formulated using an updated Lagrangian description in order to account for geometrically non-linear effects in the deformation of the FRP bars. Following standard finite element procedures, the discretized equilibrium equation for the bar structure, expressed in the current configuration at time <math display="inline">t</math>, can be written in the following quasi-static form <span id='citeF-15'></span>[[#cite-15|[15]]] | ||
Line 247: | Line 253: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>{}^t {\boldsymbol K}_{T_r} = \int _{V_r} {}^t{\boldsymbol B}_r^T {}^t {\boldsymbol D}_{T_r} {}^t{\boldsymbol B}_r dV\, , \qquad \mbox{with } r=s,b </math> | + | | style="text-align: center;" | <math>{}^t {\boldsymbol K}_{T_r} = \int _{V_r} {}^t{\boldsymbol B}_r^T {}^t {\boldsymbol D}_{T_r} {}^t{\boldsymbol B}_r dV\, , \qquad \mbox{with }\,\,\, r=s,b </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (13) | | style="width: 5px;text-align: right;white-space: nowrap;" | (13) | ||
Line 256: | Line 262: | ||
The computation of these matrices is detailed in the next sections. | The computation of these matrices is detailed in the next sections. | ||
− | ==3 Formulation of 3D solid elements for the transverse reinforced resin material== | + | ==3. Formulation of 3D solid elements for the transverse reinforced resin material== |
===3.1 Definition of displacements, strains and stresses=== | ===3.1 Definition of displacements, strains and stresses=== | ||
Line 416: | Line 422: | ||
{| style="text-align: left; margin:auto;width: 100%;" | {| style="text-align: left; margin:auto;width: 100%;" | ||
|- | |- | ||
− | | style="text-align: center;" | <math>\hat{\boldsymbol D}^{-1} = V_m [{\boldsymbol D}_m(1-d_m)]^{-1} + V_f {\boldsymbol D}_f(1-d_f)]^{-1} </math> | + | | style="text-align: center;" | <math>\hat{\boldsymbol D}^{-1} = V_m \left[{\boldsymbol D}_m(1-d_m)\right]^{-1} + V_f \left[{\boldsymbol D}_f(1-d_f)\right]^{-1} </math> |
|} | |} | ||
| style="width: 5px;text-align: right;white-space: nowrap;" | (24) | | style="width: 5px;text-align: right;white-space: nowrap;" | (24) | ||
Line 449: | Line 455: | ||
|} | |} | ||
− | ==4 Formulation of rotation-free beam elements for the fibers== | + | ==4. Formulation of rotation-free beam elements for the fibers== |
The elimination of the rotations as nodal degrees of freedom in bending-type elements (such as beams, plates and shells) leads to the so-called ''rotation-free'' elements <span id='citeF-20'></span><span id='citeF-21'></span><span id='citeF-22'></span>[[#cite-20|[20,21,22]]]. Rotations can be simply eliminated by expressing the curvature field in a node, or a point within an element, in terms of the displacements of nodes contained in an appropriate ''patch of elements'' surrounding the node or the point. This procedure emanates from the efforts of finite difference practitioners for analysis of thin plates based on Kirchhoff plate theory using displacement variables only <span id='citeF-20'></span>[[#cite-20|[20]]]. The idea evolved in the last decades of the 20th century and lead to the formulation of different families of rotation-free elements for analysis of plates and shells. Oñate and Zárate <span id='citeF-20'></span>[[#cite-20|[20]]] developed a unified formulation for this type of rotation-free plate and shell elements combining ideas from the finite volume method with the FEM. Some of the rotation-free shell elements have enjoyed much popularity for non-linear analysis of shells structures under dynamic and impact loads and sheet stamping processes using explicit and implicit dynamic formulations <span id='citeF-23'></span>[[#cite-23|[23]]]. | The elimination of the rotations as nodal degrees of freedom in bending-type elements (such as beams, plates and shells) leads to the so-called ''rotation-free'' elements <span id='citeF-20'></span><span id='citeF-21'></span><span id='citeF-22'></span>[[#cite-20|[20,21,22]]]. Rotations can be simply eliminated by expressing the curvature field in a node, or a point within an element, in terms of the displacements of nodes contained in an appropriate ''patch of elements'' surrounding the node or the point. This procedure emanates from the efforts of finite difference practitioners for analysis of thin plates based on Kirchhoff plate theory using displacement variables only <span id='citeF-20'></span>[[#cite-20|[20]]]. The idea evolved in the last decades of the 20th century and lead to the formulation of different families of rotation-free elements for analysis of plates and shells. Oñate and Zárate <span id='citeF-20'></span>[[#cite-20|[20]]] developed a unified formulation for this type of rotation-free plate and shell elements combining ideas from the finite volume method with the FEM. Some of the rotation-free shell elements have enjoyed much popularity for non-linear analysis of shells structures under dynamic and impact loads and sheet stamping processes using explicit and implicit dynamic formulations <span id='citeF-23'></span>[[#cite-23|[23]]]. | ||
Line 461: | Line 467: | ||
A particular class of rotation-free beam element can be derived by computing the curvature ''at each node'' using a finite difference scheme. The resulting element is termed CVB (for Cell-Vertex Beam element) <span id='citeF-14'></span><span id='citeF-20'></span>[[#cite-14|[14,20]]]. | A particular class of rotation-free beam element can be derived by computing the curvature ''at each node'' using a finite difference scheme. The resulting element is termed CVB (for Cell-Vertex Beam element) <span id='citeF-14'></span><span id='citeF-20'></span>[[#cite-14|[14,20]]]. | ||
− | For the sake of simplicity we will show the derivation of the stiffness equation of a rotation-free straight beam element that bends in the <math display="inline">x-z</math> plane ( | + | For the sake of simplicity we will show the derivation of the stiffness equation of a rotation-free straight beam element that bends in the <math display="inline">x-z</math> plane ([[#img-3|Figure 3]]). A similar process can be followed for deriving the stiffness equations for bending in the <math display="inline">x-y</math> plane. Due to the particular one-dimensional (1D) features of the rebars considered in this work, we will neglect torsional effects in the fibers. |
<div id='img-3'></div> | <div id='img-3'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure2.png|600px|Modelling strategy for analysis of FRP with the FEM]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 3'''. Modelling strategy for analysis of FRP with the FEM |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
Let us formulate the deformation of a straight beam element for a current configuration at time <math display="inline">t</math>. The bending behavior of the beam is modelled by considering the control domain formed by half of the lengths of the elements adjacent to a node. The curvature at node <math display="inline">i</math> is computed as | Let us formulate the deformation of a straight beam element for a current configuration at time <math display="inline">t</math>. The bending behavior of the beam is modelled by considering the control domain formed by half of the lengths of the elements adjacent to a node. The curvature at node <math display="inline">i</math> is computed as | ||
Line 484: | Line 491: | ||
|} | |} | ||
− | where <math display="inline">l_i={1\over 2} (l^e+ l^{e-1})</math> and subscripts <math display="inline">R_i</math> and <math display="inline">L_i</math> denote the midpoints of the elements located at the right and left of node <math display="inline">i</math> (Fig. [[#img-4|4]]). The curvature <math display="inline">\kappa _i</math> is assumed to be constant in the control domain <math display="inline">l_i</math> assigned to node <math display="inline">i</math> ( | + | where <math display="inline">l_i={1\over 2} (l^e+ l^{e-1})</math> and subscripts <math display="inline">R_i</math> and <math display="inline">L_i</math> denote the midpoints of the elements located at the right and left of node <math display="inline">i</math> (Fig. [[#img-4|4]]). The curvature <math display="inline">\kappa _i</math> is assumed to be constant in the control domain <math display="inline">l_i</math> assigned to node <math display="inline">i</math> ([[#img-4|Figure 4]]). |
<div id='img-4'></div> | <div id='img-4'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig1_15.png|440px|CVB element. Control domains around two nodes i and i+1]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 4'''. CVB element. Control domains around two nodes <math>i</math> and <math>i+1</math> |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
The rotations <math display="inline">\left(\displaystyle {\partial w\over \partial x}\right)_{\!\scriptscriptstyle R_i}</math> and <math display="inline">\left(\displaystyle {\partial w\over \partial x}\right)_{\!\scriptscriptstyle L_i}</math> are expressed in terms of the nodal deflections as | The rotations <math display="inline">\left(\displaystyle {\partial w\over \partial x}\right)_{\!\scriptscriptstyle R_i}</math> and <math display="inline">\left(\displaystyle {\partial w\over \partial x}\right)_{\!\scriptscriptstyle L_i}</math> are expressed in terms of the nodal deflections as | ||
Line 634: | Line 642: | ||
The prescribed values for the deflection are imposed when solving the global system of equations, as usual. We briefly explain next how to compute the curvature at a free end and at simply supported and clamped nodes. For details see <span id='citeF-14'></span>[[#cite-14|[14]]]. | The prescribed values for the deflection are imposed when solving the global system of equations, as usual. We briefly explain next how to compute the curvature at a free end and at simply supported and clamped nodes. For details see <span id='citeF-14'></span>[[#cite-14|[14]]]. | ||
+ | The curvature at a ''free end node and at a simply supported node'' is ''zero''. This condition is implemented by neglecting the contribution of the boundary node to the element stiffness matrix. For a boundary element with a free or simply supported node the stiffness matrix is ([[#img-5|Figure 5]]a) | ||
<span id="eq-36"></span> | <span id="eq-36"></span> | ||
− | |||
− | |||
− | |||
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
|- | |- | ||
Line 659: | Line 665: | ||
<div id='img-5'></div> | <div id='img-5'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:65%;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig1_16.png|480px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 5'''. CVB element. (a) Boundary condition of zero curvature at a free or simply supported node <math>B</math>. (b) Control domain for a clamped or symmetry node |
− | + | ||
− | + | ||
|} | |} | ||
− | The curvature at a ''clamped or symmetry left-end node'' <math display="inline">i</math> is computed as ( | + | |
+ | The curvature at a ''clamped or symmetry left-end node'' <math display="inline">i</math> is computed as ([[#img-5|Figure 5]]b) | ||
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 680: | Line 687: | ||
where <math display="inline">w_{\scriptscriptstyle i-1}</math> is an auxiliary fictitious deflection. | where <math display="inline">w_{\scriptscriptstyle i-1}</math> is an auxiliary fictitious deflection. | ||
− | Similarly, for a right-end clamped or symmetry node ( | + | Similarly, for a right-end clamped or symmetry node ([[#img-5|Figure 5]]b) |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 730: | Line 737: | ||
===4.4 Extension to curved rotation-free CVB elements=== | ===4.4 Extension to curved rotation-free CVB elements=== | ||
− | A curved rebar will be modelled by a collection of straight rotation-free CVB elements | + | A curved rebar will be modelled by a collection of straight rotation-free CVB elements. [[#img-6|Figure 6]] shows a curved rebars under bending loads in a plane <math display="inline">x-z</math>, discretized with three CVB elements. |
<div id='img-6'></div> | <div id='img-6'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:50%;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig9_26.png|390px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 6'''. Cell-vertex rotation-free beam element. Control domains for computing the average curvature. A, B and C are element mid-points |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
The stiffness matrix is obtained as the sum of the axial stiffness matrix (<math display="inline">{\boldsymbol K}_A</math>) and the bending stiffness matrix (<math display="inline">{\boldsymbol K}_B</math>). | The stiffness matrix is obtained as the sum of the axial stiffness matrix (<math display="inline">{\boldsymbol K}_A</math>) and the bending stiffness matrix (<math display="inline">{\boldsymbol K}_B</math>). | ||
Line 776: | Line 784: | ||
The procedure for computing the bending stiffness matrix is an extension of that followed for straight CVB elements in Section 4.2. <span id="eq-41"></span> | The procedure for computing the bending stiffness matrix is an extension of that followed for straight CVB elements in Section 4.2. <span id="eq-41"></span> | ||
− | The curvature at node <math display="inline">i</math> is computed in terms of the gradients of the normal deflection at points <math display="inline">A</math> and <math display="inline">B</math> ( | + | The curvature at node <math display="inline">i</math> is computed in terms of the gradients of the normal deflection at points <math display="inline">A</math> and <math display="inline">B</math> ([[#img-6|Figure 6]]) |
<span id="eq-41.a"></span> | <span id="eq-41.a"></span> | ||
Line 791: | Line 799: | ||
|} | |} | ||
− | In Eqs.([[#eq-41.a|41.a]]) <math display="inline">{w'}_j</math> are the nodal deflections in the local axes of the element defined in | + | In Eqs.([[#eq-41.a|41.a]]) <math display="inline">{w'}_j</math> are the nodal deflections in the local axes of the element defined in [[#img-7|Figure 7]] and |
<span id="eq-41.b"></span> | <span id="eq-41.b"></span> | ||
Line 818: | Line 826: | ||
<div id='img-7'></div> | <div id='img-7'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure3.png|600px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 7'''. Local coordinate system of curved rotation-free CVB element bending in the <math>xz</math> plane |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
Similarly, for node <math display="inline">i+1</math> | Similarly, for node <math display="inline">i+1</math> | ||
Line 866: | Line 875: | ||
|} | |} | ||
− | The bending stiffness matrix for the CVB element is expressed as the sum of the stiffness contributions from the two control subdomains 1 and 2 that form the element ( | + | The bending stiffness matrix for the CVB element is expressed as the sum of the stiffness contributions from the two control subdomains 1 and 2 that form the element ([[#img-6|Figure 6]]), i.e. |
{| class="formulaSCP" style="width: 100%; text-align: left;" | {| class="formulaSCP" style="width: 100%; text-align: left;" | ||
Line 897: | Line 906: | ||
The assembly process follows the general rule of the FEM <span id='citeF-14'></span>[[#cite-14|[14]]]. | The assembly process follows the general rule of the FEM <span id='citeF-14'></span>[[#cite-14|[14]]]. | ||
− | ==5 Coupling of 3D solid elements and rotation-free beam elements== | + | ==5. Coupling of 3D solid elements and rotation-free beam elements== |
− | + | [[#img-8|Figure 8]] shows a conceptual representation of a FRP bar of rectangular section containing of a number of fibers embedded in the (bulk) matrix material. The bar geometry is discretized into 3D solid elements (8-noded hexahedral are shown). The fibers contained within each solid element are assumed to be concentrated in one-dimensional (1D) packs of fibers distributed within the 3D solid elements. [[#img-8|Figure 8]] shows the concentration of fibers in four 1D packs placed at the <math display="inline">2\times 2</math> Gauss points locations in each cross section <span id='citeF-14'></span>[[#cite-14|[14]]]. Each fiber pack is modelled as a rotation-free beam element with a cross-section area equivalent to the amount of fibers considered in the pack. | |
The displacements of the rotation-free beam element nodes are related to the nodal displacement of the 3D solid element by a simple interpolation. | The displacements of the rotation-free beam element nodes are related to the nodal displacement of the 3D solid element by a simple interpolation. | ||
<div id='img-8'></div> | <div id='img-8'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure4.png|500px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 8'''. Patch of 3D solid and rotation-free beam elements |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
The discretized equilibrium equations for the coupled solid-beam system can be written as explained in Section [[#2 Combination of 3D solid elements and rotation-free elements for analysis of polymer fiber rebars|2]]. The numerical response of the FRP bars for different incremental load up to failure is computed using the standard Newton-Raphson iterative procedure described in Eqs.[[#eq-4|(4)]]-[[#eq-13|(13)]]. | The discretized equilibrium equations for the coupled solid-beam system can be written as explained in Section [[#2 Combination of 3D solid elements and rotation-free elements for analysis of polymer fiber rebars|2]]. The numerical response of the FRP bars for different incremental load up to failure is computed using the standard Newton-Raphson iterative procedure described in Eqs.[[#eq-4|(4)]]-[[#eq-13|(13)]]. | ||
− | ==6 Examples of application== | + | ==6. Examples of application== |
===6.1 Composite material properties=== | ===6.1 Composite material properties=== | ||
Line 927: | Line 937: | ||
===6.2 Axial tension test=== | ===6.2 Axial tension test=== | ||
− | + | [[#img-9|Figure 9]] shows the bar modelled with 960 8-noded hexahedral elements and 1197 nodes. The fibers have been modelled using 1140 RFB elements. The length of the sample is 100 mm long with 9.58 mm diameter. Symmetry conditions have been considered by prescribing a zero axial displacement in the normal direction to the surface at a bar end, while the other end has a prescribed axial displacement condition to induce tensile stresses in the bar. The bar axis has the displacements in the normal plane yz constrained to zero. | |
<div id='img-9'></div> | <div id='img-9'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure_8.png|400px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 9'''. Tensile test geometry: 3D solid and RFB elements mesh |
− | + | ||
− | + | ||
|} | |} | ||
− | + | ||
+ | [[#img-10|Figure 10]] shows a typical bar tested at UPM after failure. | ||
<div id='img-10'></div> | <div id='img-10'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_2.png|220px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 10'''. Tensile test of the PFR performed at UPM |
− | + | ||
− | + | ||
|} | |} | ||
− | The numerical results obtained are shown in | + | |
+ | The numerical results obtained are shown in [[#img-11|Figures 11]] and [[#img-12|12]]. The displacement field obtained presents a linear distribution, until the rebar reaches brittle failure. | ||
<div id='img-11'></div> | <div id='img-11'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_3.png|500px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 11'''. Displacement field [m] in deformed geometry |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
<div id='img-12'></div> | <div id='img-12'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_4.png|484px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 12'''. Axial tension tests stress-strain curves. Numerical and experimental values |
− | + | ||
− | + | ||
|} | |} | ||
− | |||
− | + | The numerical load-displacement curve obtained is converted into a stress-strain curve to be compared with the experimental results. A good correlation can be seen between the numerical and experimental values ([[#img-12|Figure 12]]). | |
+ | [[#table-1|Table 1]] shows the good correlation between the numerical and experimental stiffness and strength magnitudes. | ||
− | + | <div class="center" style="font-size: 75%;">'''Table 1'''. Tensile test. Comparison of numerical and experimental results</div> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
+ | <span id='table-1'></span> | ||
+ | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" | ||
+ | |-style="text-align:center" | ||
+ | ! Property !! Lab tests !! Numerical | ||
+ | |- style="text-align:center" | ||
+ | | Tensile strength (MPa) | ||
+ | | 1150 <math>\boldsymbol{\pm }</math> 34 (MPa) | ||
+ | | 1173.33 (MPa) | ||
+ | |- style="text-align:center" | ||
+ | | Strain to max strength | ||
+ | | 2.3 <math>\pm </math> 0.1 (%) | ||
+ | | 2.22 (%) | ||
+ | |- style="text-align:center" | ||
+ | | Young Modulus (GPa) | ||
+ | | 55 <math>\pm </math> 4 (GPa) | ||
+ | | 53.03 (GPa) | ||
|} | |} | ||
===6.3 Bending test=== | ===6.3 Bending test=== | ||
− | The experiment performed for measuring the bending properties, corresponds to the standard three-point bending test showed in | + | The experiment performed for measuring the bending properties, corresponds to the standard three-point bending test showed in [[#img-14|Figure 14]]. The tests were carried out on bars of full circular section, instead of using cut specimens as indicated in the ASTMD-4476 norm. [[#img-13|Figure 13]] shows the bar modelled with 5400 8-noded hexahedral elements (6327 nodes) and 6156 rotation-free beam elements. The length of the bar sample is 100 mm long with 9.58 mm diameter. Symmetry and boundary conditions have been considered, as shown in [[#img-13|Figure 13]]. To avoid problems with the stress concentration at the loading points, a small part of the beam (the blue zone in [[#img-13|Figure 13]]) has been defined with an unbreakable material. The elastic properties of this material are the same as for the rest of the bar. |
<div id='img-13'></div> | <div id='img-13'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:50%;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_5.png|450px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 13'''. Three point bending test for PFR rebar. Dimensions, FEM mesh and boundary conditions |
− | + | ||
− | + | ||
|} | |} | ||
− | + | ||
+ | [[#img-14|Figure 14]] shows a snapshot of the laboratory test for the FRP rebar in a three-point bending configuration. Note the rollers used as a pin-supports of the beam. | ||
<div id='img-14'></div> | <div id='img-14'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_6.png|434px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 14'''. Three-point bending test of FRP rebars performed at UPM |
− | + | ||
− | + | ||
|} | |} | ||
− | + | ||
+ | [[#img-15|Figure 15]] shows the strain distribution over the deformed bar ([[#img-15|Figure 15]]a) and the mesh of hexahedra ([[#img-15|Figure 15]]b) in a load state close to the onset of damage. The damage distribution for each element is displayed in [[#img-16|Figure 16]]. | ||
<div id='img-15'></div> | <div id='img-15'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:50%;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_7.png|450px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 15'''. Three-point bending tests. Numerical results. Deformed bar. (a) Strain distribution for RFB elements. (b) Strain distribution in 3D solid elements |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
<div id='img-16'></div> | <div id='img-16'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:50%;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure_15.png|450px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 16'''. Contour field of damage variable in longitudinal and transverse sections on the deformed geometry (1x) |
− | + | ||
− | + | ||
|} | |} | ||
− | |||
− | + | In [[#img-17|Figure 17]], the numerical load-deflection curve is shown together with that obtained in the laboratory for comparison purposes. The maximum loading force and the flexural stiffness captured by the numerical simulation are in good agreement with the experimental results. | |
<div id='img-17'></div> | <div id='img-17'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:60%;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_9.png|480px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 17'''. Three-point bending tests. Force-deflection curves for shearing testing. Numerical results compared with laboratory tests |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
+ | |||
+ | At the laboratory, samples of 240mm of total length were cut for the 10mm diameter full cylinder tests. The span of the tests was set to 200mm. Tests were carried out under displacements control, with a crosshead displacement velocity of 3mm/min. [[#img-18|Figure 18]] shows the test set-up and one specimen after failure. | ||
<div id='img-18'></div> | <div id='img-18'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_10.png|450px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 18'''. Three-point bending tests. Failed specimen of PFR rebar |
− | + | ||
− | + | ||
|} | |} | ||
===6.4 Shear test=== | ===6.4 Shear test=== | ||
− | The experiment carried out for evaluating the shear properties, corresponds to the ASTM-7617 standard test showed in | + | The experiment carried out for evaluating the shear properties, corresponds to the ASTM-7617 standard test showed in [[#img-19|Figure 19]] <span id='citeF-25'></span>[[#cite-25|[25]]]. For the numerical analysis, the symmetry of the geometry and the imposed displacement have been taken into account as shown in [[#img-20|Figure 20]]a in order to speed up the calculations. [[#img-20|Figure 20]]b displays the bar modelled with 5225 8-noded hexahedral elements (4608 nodes) and 5016 rotation-free beam elements. The length of the sample is 100 mm long with 9.58 mm diameter. |
<div id='img-19'></div> | <div id='img-19'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |align="center" | | ||
+ | {| | ||
+ | |+ | ||
|- | |- | ||
− | |[[Image:Draft_Onate_870101189-Fig6_11a.png| | + | |style="text-align: center;padding:10px;"|[[Image:Draft_Onate_870101189-Fig6_11a.png|180px|]] |
− | |[[Image:Draft_Onate_870101189-Fig6_11b.png| | + | |style="text-align: center;padding:10px;"|[[Image:Draft_Onate_870101189-Fig6_11b.png|162px|]] |
|- | |- | ||
− | | colspan="2" style=" padding:10px;" |[[Image:Draft_Onate_870101189-Fig6_11c.png| | + | |colspan="2" style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_11c.png|330px]] |
− | + | ||
− | + | ||
|} | |} | ||
+ | |- | ||
+ | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 19'''. Shear test. Double shear fixture according to ASTM D-7617 standard <span id='citeF-25'></span>[[#cite-25|[25]]] | ||
+ | |} | ||
+ | |||
<div id='img-20'></div> | <div id='img-20'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |align="center" | | ||
+ | {| | ||
+ | |+ | ||
|- | |- | ||
− | |[[Image:Draft_Onate_870101189-Fig6_12a.png| | + | |style="text-align: center;padding:10px;"|[[Image:Draft_Onate_870101189-Fig6_12a.png|500px|]] |
− | |[[Image:Draft_Onate_870101189-Fig6_12b.png| | + | |- |
− | |- style="text-align: center; font-size: 75%;" | + | |style="text-align: center;font-size: 75%;"|(a) |
− | | | + | |- |
+ | |style="text-align: center;padding-bottom:10px;"| [[Image:Draft_Onate_870101189-Fig6_12b.png|500px]] | ||
+ | |- | ||
+ | |style="text-align: center;font-size: 75%;padding-bottom:10px;"|(b) | ||
+ | |} | ||
+ | |- | ||
+ | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 20'''. Shear test geometry. (a) Boundary conditions. (b) Mesh and dimensions | ||
|} | |} | ||
− | |||
− | + | In [[#img-20|Figure 20]] a gap of 0.001 m can be observed between the loaded and supported zones of the specimen. This zone is necessary to capture accurately the shear forces and to avoid stress concentration in the numerical model. | |
− | + | [[#img-21|Figure 21]] shows a representative configuration of the numerical results. [[#img-21|Figure 21]]a displays the deformed bar and the localized damage in the resin material at the shear zone. [[#img-21|Figure 21]]b is a close-up of the shear zone were the RFB are presented. The curvature of these elements due to the shear deformation is depicted. | |
+ | |||
+ | [[#img-22|Figure 22]] shows the numerical load-deflection curve superimposed to a set of curves obtained in the laboratory tests for comparison purposes. In the curves the non-linear load phenomena of the fiber reorientation can be clearly observed. For an applied load slightly larger than 5kN there is an initial plateau in the experimental tests, also reproduced by the numerical model, that corresponds to the cracking of the matrix. At this point there is a fiber reorientation that provides the stiffness increase and the final composite failure, due to fiber failure, for a load close to 35kN. The good agreement obtained between the numerical and the experimental results, especially in this case, highly non-linear, proves the validity of the proposed procedure to simulate pultruded fiber reinforced polymer rebars. | ||
<div id='img-21'></div> | <div id='img-21'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Fig6_13.png|600px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 21'''. Shear test. Contour field of damage variable in the matrix. Deformed geometry (1x) |
− | + | ||
− | + | ||
|} | |} | ||
+ | |||
<div id='img-22'></div> | <div id='img-22'></div> | ||
− | {| class=" | + | {| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;width:auto;" |
+ | |-style="background:white;" | ||
+ | |style="text-align: center;padding:10px;"| [[Image:Draft_Onate_870101189-Figure_21.png|600px]] | ||
|- | |- | ||
− | | | + | | style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 22'''. Shear test. Force-deflection curves. Numerical results compared with laboratory tests |
− | + | ||
− | + | ||
|} | |} | ||
− | ==7 Concluding remarks== | + | ==7. Concluding remarks== |
The combination of 3D solid finite elements with rotation-free beam elements has shown to be adequate for modelling axial, bending and shear deformation states in fiber reinforced polymer rebars. The fact that both types of elements are formulated in terms of displacement DOFs only simplifies the combination of the best features of both elements for modelling the 3D behavior of the matrix material and the 1D behavior of the fibers, as well as the coupling interactions during complex deformation modes. | The combination of 3D solid finite elements with rotation-free beam elements has shown to be adequate for modelling axial, bending and shear deformation states in fiber reinforced polymer rebars. The fact that both types of elements are formulated in terms of displacement DOFs only simplifies the combination of the best features of both elements for modelling the 3D behavior of the matrix material and the 1D behavior of the fibers, as well as the coupling interactions during complex deformation modes. | ||
Line 1,105: | Line 1,140: | ||
The approach presented in this work for combining 3D solid elements with rotation-free beam elements can be implemented in more sophisticated formulations for modelling fiber reinforced materials, such as that proposed in <span id='citeF-26'></span>[[#cite-26|[26]]]. | The approach presented in this work for combining 3D solid elements with rotation-free beam elements can be implemented in more sophisticated formulations for modelling fiber reinforced materials, such as that proposed in <span id='citeF-26'></span>[[#cite-26|[26]]]. | ||
− | == | + | ==Data Availability Statement== |
Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request. | Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request. | ||
− | == | + | ==Acknowledgements== |
This research was carried out as a part of the project “Development of polymer fiber reinforced rebars” funded by Saudi Aramco. The authors also acknowledge the financial support from the Spanish Ministry of Economy and Competitiveness, through the “Severo Ochoa Programme for Centres of Excellence in R&D” (CEX2018-000797-S). | This research was carried out as a part of the project “Development of polymer fiber reinforced rebars” funded by Saudi Aramco. The authors also acknowledge the financial support from the Spanish Ministry of Economy and Competitiveness, through the “Severo Ochoa Programme for Centres of Excellence in R&D” (CEX2018-000797-S). | ||
− | === | + | ==References== |
+ | |||
+ | <div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;"> | ||
<div id="cite-1"></div> | <div id="cite-1"></div> | ||
− | + | [[#citeF-1|[1]]] Pacheco-Torgal F., Jalali S., Labrincha J., John V.M. Eco-efficient concrete. Woodhead Publishing Series in Civil and Structural Engineering, pp. 624, 2013. | |
<div id="cite-2"></div> | <div id="cite-2"></div> | ||
− | + | [[#citeF-2|[2]]] Bank L.C. Composites for construction: Structural design with FRP materials. John Wiley & Sons, Inc., pp. 560, 2006. | |
<div id="cite-3"></div> | <div id="cite-3"></div> | ||
− | + | [[#citeF-3|[3]]] Marí A.R., Oller E., Bairán J.M. Predicting the response of FRP-strengthened reinforced-concrete flexural members with nonlinear evolutive analysis models. J. Compos. Constr., 15(5):799–809, 2011. | |
<div id="cite-4"></div> | <div id="cite-4"></div> | ||
− | + | [[#citeF-4|[4]]] Martínez X., Oller S., Rastellini F., Barbat A.H. A numerical procedure for simulating RC structures reinforced with FRP using the serial/parallel mixing theory. Computers & Structures, 86(15-16):1604-1618, 2008. doi:[https://doi.org/10.1016/j.compstruc.2008.01.007 10.1016/j.compstruc.2008.01.007] | |
<div id="cite-5"></div> | <div id="cite-5"></div> | ||
− | + | [[#citeF-5|[5]]] Nogales A., Tosić N., de la Fuente A. Rotation and moment redistribution capacity of fiber-reinforced concrete beams: Parametric analysis and code compliance. Struct. Concr., 23(1):220-239, 2021. | |
<div id="cite-6"></div> | <div id="cite-6"></div> | ||
− | + | [[#citeF-6|[6]]] Pultron Composites [Online]. Available: https://www.pultron.com/products/features/. [Accessed: 17-Apr-2020]. | |
<div id="cite-7"></div> | <div id="cite-7"></div> | ||
− | + | [[#citeF-7|[7]]] Lleida GFRP Pedestrian Bridge [Online]. Available: http://www.pedelta.com/lleida-gfrp-pedestrian-bridge-p-52-en. [Accessed: 17-Apr-2020]. | |
<div id="cite-8"></div> | <div id="cite-8"></div> | ||
− | + | [[#citeF-8|[8]]] Hernández S., et al. Design and development of an innovative structural system based on composites applied on the constrction of maritme lighthouses with lower maintenance requeriments (FAROCOMP). Mater. Compuestos, 3(2):54–58, 2019. | |
<div id="cite-9"></div> | <div id="cite-9"></div> | ||
− | + | [[#citeF-9|[9]]] GangaRao H.V.S., Taly N., Vijay P.V. Reinforced concrete design with FRP composites. 1st Edition, CRC Press, pp. 400, 2019. | |
<div id="cite-10"></div> | <div id="cite-10"></div> | ||
− | + | [[#citeF-10|[10]]] Sharifianjazi F., Zeydi P., Bazli M., Esmaeilkhanian A., Rahmani R., Bazli L., Khaksar S. Fibre-reinforced polymer reinforced concrete members under elevated temperatures: A review on structural performance. Polymers, 14(3):472, 2022. [https://doi.org/10.3390/polym14030472 10.3390/polym14030472] | |
<div id="cite-11"></div> | <div id="cite-11"></div> | ||
− | + | [[#citeF-11|[11]]] Sagar B., Sivakumar M.V.N. Performance evaluation of basalt fibre-reinforced polymer rebars in structural concrete members–a review. Innovative Infrastructure Solutions, 6(2):1-18, 2021. [https://doi.org/10.1007/s41062-020-00452-2 10.1007/s41062-020-00452-2] | |
<div id="cite-12"></div> | <div id="cite-12"></div> | ||
− | + | [[#citeF-12|[12]]] Strau S., Senz A., Ellinger J. Comparison of the processing of epoxy resins in pultrusion with open bath impregnation and closed-injection pultrusion. J. Compos. Sci., 3(3):87, 2019. | |
<div id="cite-13"></div> | <div id="cite-13"></div> | ||
− | + | [[#citeF-13|[13]]] Starr T.F. Pultrusion for engineers. Woodhead Publishing Series in Composites Science and Engineering, pp. 336, 2000. | |
<div id="cite-14"></div> | <div id="cite-14"></div> | ||
− | + | [[#citeF-14|[14]]] Oñate E. Structural analysis with the finite element method. Linear statics. Vol 1: Basis and solids, 2009. Vol. 2: Beams, plates and shells, 2013, Springer-CIMNE. | |
<div id="cite-15"></div> | <div id="cite-15"></div> | ||
− | + | [[#citeF-15|[15]]] Zienkiewicz O.C., Taylor R.L. The finite element method for solids and structural mechanics. 6th edition, Elsevier, 2005. | |
<div id="cite-16"></div> | <div id="cite-16"></div> | ||
− | + | [[#citeF-16|[16]]] Martínez X., Rastellini F., Oller S., Flores F., Oñate E. Computationally optimized formulation for the simulation of composite materials and delamination failures. Composites Part B: Engineering, 42(2):134-144, 2011. [https://doi.org/10.1016/j.compositesb.2010.09.013 10.1016/j.compositesb.2010.09.013] | |
<div id="cite-17"></div> | <div id="cite-17"></div> | ||
− | + | [[#citeF-17|[17]]] Rastellini F., Oller O., Salomón O., Oñate E. Composite materials non-linear modelling for long fibre-reinforced laminates: Continuum basis, computational aspects and validations. Computers & Structures, 86(9):879-896, 2008. [https://doi.org/10.1016/j.compstruc.2007.04.009 10.1016/j.compstruc.2007.04.009] | |
<div id="cite-18"></div> | <div id="cite-18"></div> | ||
− | + | [[#citeF-18|[18]]] Oller S., Oñate E., Oliver J., Lubliner J. Finite element non-linear analysis of concrete structures using a plastic-damage model. Int. Journal of Engineering Fracture Mechanics, 35(1/2/3):219-231, 1990. [https://doi.org/10.1016/0013-7944(90)90200-Z 10.1016/0013-7944(90)90200-Z] | |
<div id="cite-19"></div> | <div id="cite-19"></div> | ||
− | + | [[#citeF-19|[19]]] Lubliner J., Oller S., Oliver J., Oñate E. A plastic damage model for concrete. Int. Journal of Solids and Structures, 25(3):299-326, 1989. [https://doi.org/10.1016/0020-7683(89)90050-4 10.1016/0020-7683(89)90050-4] | |
<div id="cite-20"></div> | <div id="cite-20"></div> | ||
− | + | [[#citeF-20|[20]]] Oñate E., Zárate F. Rotation-free triangular plate and shell elements. Int. J. Numer. Meth. Engng., 47(1-3):557-603, 2000. [https://doi.org/10.1002/(SICI)1097-0207(20000110/30)47:1/3<557::AID-NME784>3.0.CO;2-9 10.1002/(SICI)1097-0207(20000110/30)47:1/3<557::AID-NME784>3.0.CO;2-9] | |
<div id="cite-21"></div> | <div id="cite-21"></div> | ||
− | + | [[#citeF-21|[21]]] Oñate E., Zárate F. Extended rotation-free plate and beam elements with shear deformation effects. Int. J. Numer. Meth. Engng., 83(2):196-227, 2010. [https://doi.org/10.1002/nme.2836 10.1002/nme.2836] | |
<div id="cite-22"></div> | <div id="cite-22"></div> | ||
− | + | [[#citeF-22|[22]]] Jovicevic J., Oñate E. Analysis of beams and shells using a rotation - free finite element - finite volume formulation. Monograph M43, CIMNE, ISBN: 84-89925-36-4, Barcelona, 1999. | |
<div id="cite-23"></div> | <div id="cite-23"></div> | ||
− | + | [[#citeF-23|[23]]] Oñate E., Flores F.G. Advances in the formulation of the rotation-free basic shell triangle. Computer Methods in Applied Mechanics and Engineering, 194(21-24):2406-2443, 2005. [https://doi.org/10.1016/j.cma.2004.07.039 10.1016/j.cma.2004.07.039] | |
<div id="cite-24"></div> | <div id="cite-24"></div> | ||
− | + | [[#citeF-24|[24]]] Flores F., Oñate E. Rotation-free finite element for the non-linear analysis of beams and axisymmetric shells. Computer Methods in Applied Mechanics and Engineering, 195(41-43):5297-5315, 2006. [https://doi.org/10.1016/j.cma.2005.08.021 10.1016/j.cma.2005.08.021] | |
<div id="cite-25"></div> | <div id="cite-25"></div> | ||
− | + | [[#citeF-25|[25]]] ASTM-International, ASTM D7617/D7617M - 11(2017) Standard Test Method for Transverse Shear Strength of Fiber-reinforced Polymer Matrix Composite Bars, West Conshohocken, PA, 2017. | |
<div id="cite-26"></div> | <div id="cite-26"></div> | ||
− | + | [[#citeF-26|[26]]] Khristenko U., Schuß S., Krüger M., Schmidt F., Wohlmuth B., Hesch C. Multidimensional coupling: A variationally consistent approach to fiber-reinforced materials. Computer Methods in Applied Mechanics and Engineering, 382:113869, 2021. [https://doi.org/10.1016/j.cma.2021.113869 10.1016/j.cma.2021.113869] |
We present a combination of three dimensional (3D) solid elements and rotation-free beam elements for non-linear analysis of fiber-reinforced polymer (FRP) of rebars. The matrix material is modelled by 3D solid elements while the fibers are modelled with rotation-free beam elements. The absence of rotation variables in the beam elements allows the straightforward coupling of 3D solid and beam elements using a formulation with displacement nodal degrees of freedom only. Both solid and beam elements are formulated in an updated Lagrangian description. The behavior of the matrix and the fiber material are modelled with an elastic-damage model. The efficiency and accuracy of the combined 3D-beam element formulation are verified in examples of application to the analysis of FRP rebars up to fracture in axial, bending and shear tests for which experimental results are available.
Keywords: 3D solid finite elements, rotation-free beam elements, polymer fiber rebars
There is an increasing interest in the use of non-metallic materials, such as polymers, in the building and construction (B&C) sector as an alternative to standard materials, such as steel and cement [1,2]. These new materials have many advantages, compared to traditional ones, in terms of strength-weight and stiffness-weigh ratios, corrosion resistance to chemical agents, and fatigue performance, to name a few. All these advantages can extend the life-frame of structures and minimize the maintenance required, thus contributing to improve the sustainability of the sector. Current research efforts aim to develop resistant, sustainable and cost-effective materials incorporating polymers for a number of applications in building and construction.
Polymers and composite materials can be used in B&C in many forms: as external reinforcements of concrete beams and slabs [3,4], as short fibers embedded in concrete [5], or as pultruded elements which can be used as rebars for concrete reinforcement [6] or as beams elements to construct different structures such as bridges [7], or lighthouses [8] (Figures 1a and 1b), to name some examples.
| ||||
Figure 1. Applications of pultruded beams in construction [7,8] |
The use of pultruded FRP as reinforcement in concrete has several advantages compared to steel reinforcement such as corrosion resistance, high longitudinal (unidirectional) strength (2 or 3 times higher than steel for equivalent bar diameter), high fatigue endurance (which depends on the type of reinforcing fiber and bar), magnetic transparency, lightweight (25% of steel weight) and low thermal and electric conductivity (for glass and aramid fibers) [9]. Glass fibers are the common choice due to their good performance and reduced cost, compared to other synthetic fibers that have improved properties but are more expensive, such as carbon. Research efforts are currently carried out for developing new families of rebars with synthetic fibers, as reinforcement for concrete structures [10,11]. This type of synthetic fiber rebars could be an alternative to current FRP rebars using glass fibers, and even to steel rebars for specific applications.
Pultrusion beams are obtained by pulling the reinforcement fibers, which have previously impregnated with a matrix resin, through a heated die where the composite is cured [12]. The shape of the die defines the shape of the final beam produced. Due to the fabrication process, the final composite contains mainly fibers aligned with the beam axis, and the efficiency of the method allows having fiber participations that can reach 80% of the total composite weight content [13].
This paper presents a new computational approach for assessing the non-linear behavior of pultruded composites beams using a combination of 3D solid finite elements with rotation-free beam elements. These last ones are included in the model to characterize the longitudinal performance of the composite.
The formulation presented is of special interest for pultruded beams, as in these elements fibers are mainly oriented in the longitudinal direction, which simplifies the distribution of the rotation-free beams in the numerical model. As it will be shown, the adaptability of the beam elements in the model geometry allows characterizing complex failure mechanisms such as shear failure, which depends highly on the tensile forces taken by the fibers once the matrix has failed. To do so, the geometrical non-linear effects in the solid and beam elements are taken into account by using an updated Lagrangian description, allowing for large displacements effects that can occur in specific loading cases, such as in shear tests and in post-failure situations. The fact that both the 3D solid elements and the rotation-free beam element are formulated in terms of displacement degrees of freedom (DOFs) only facilitates the integration of both types of elements in a unified formulation. An elastic-damage model is used for characterizing the behavior of both the matrix and the fibers.
The content of the paper is as follows. In the next section we present the overall finite element approach for non-linear analysis of FRP rebars combining solid elements and rotation-free beam elements. Then we describe the 3D solid finite element formulation for modelling the behavior of the matrix material. This is followed by the finite element formulation for the fibers using rotation-free beam elements based on a cell-vertex approach. The coupling of 3D-solid elements and rotation-free beam elements is detailed. The finite element formulation is validated in the non-linear analysis up to failure of FRP bar specimens formed by an epoxy matrix and carbon fibers tested under axial, bending and shear loads. Numerical results are compared with experimental results for the same tests. Excellent agreement is obtained in the comparison.
Figure 2 depicts the modelling strategy for analysis of FRP rebars with the FEM [14,15]. The procedure proposed dissociates the longitudinal contribution of fibers to the composite strength and stiffness, from their contribution in the transversal and shear directions. This is done by modeling a transverse reinforced resin, that accounts for the resin and the transverse and shear performance of the fibers, with standard 3D solid elements; and by modeling the longitudinal contribution of fibers to the composite with a collection of beams embedded in the 3D solid model.
![]() |
Figure 2. Finite element modelling of FRP rebar |
Both the 3D solid elements and the beam elements are formulated using an updated Lagrangian description in order to account for geometrically non-linear effects in the deformation of the FRP bars. Following standard finite element procedures, the discretized equilibrium equation for the bar structure, expressed in the current configuration at time , can be written in the following quasi-static form [15]
|
(1) |
with
|
(2) |
In Eqs.(1) and (2), is the vector of residual forces, is the vector of nodal displacement, is the vector of internal nodal forces, is the strain matrix relating the strains with the nodal displacements, is the Cauchy stress vector, is the vector of equivalent nodal forces, is the volume of the FRP bar and the upper left index denotes the time .
The vector of internal forces is computed by sum of the contributions from the solid elements and the beam elements, i.e.
|
(3) |
where indices and denote the contribution of the 3D solid elements and the beam elements to vector , respectively.
The equilibrium solution at time is found by solving Eq.(1) using a standard Newton-Raphson procedure briefly outlined below.
For each iteration of the Newton-Raphson scheme [15] we perform the following computations:
Computation of displacement increments
|
(4) |
with
|
(5) |
Displacement update
|
(6) |
Stress update
|
(7) |
|
(8) |
|
(9) |
Residual force vector computation
|
(10) |
Convergence check
The iteration stop when
|
(11) |
where is a prescribed error value. In our computations we have taken .
In Eq.(4) is the tangent stiffness matrix. We have computed this matrix as
|
(12) |
where indices and denote the contribution of the solid elements and the beam elements to the tangent stiffness matrix, respectively. These contributions are computed as
|
(13) |
where and with denote the deformation matrix and the tangent constitutive matrix for the solid elements and the beam elements, respectively at time .
The computation of these matrices is detailed in the next sections.
Let us consider a FRP rebar discretized using 3D-solid elements of nodes. The displacement field within the element can be expressed in terms of the nodal displacements in the standard finite element form as [14,15]
|
(14) |
where is the displacement vector, and are the displacements along the cartesian axes , respectively, and are the shape functions matrix and the nodal displacement vector for the element given by
|
(15) |
where
|
(16) |
are the nodal components of matrix and vector , and is the shape function matrix of node .
The strains within the element can be readily expressed in terms of the nodal displacements as
|
(17) |
where
|
(18) |
where and are the strain matrices for the element and the th node, respectively.
The stress-strain relationship is written in incremental form (neglecting initial strains and initial stressed) as
|
(19) |
where
|
(20) |
is the Cauchy stress vector and is the tangent constitutive matrix. In our computations we have modelled the matrix material with an elastoplastic model based on a Drucker-Prager yield functions and an elastic-plastic-damage model. The damage model is described in the next section.
The expressions of the strain matrix and the Cauchy stress vector are used for computing the nodal forces (Eq.(2)) and the tangent constitutive matrix.
The validation examples presented in the paper have been solved using three-linear 8-noded hexahedral elements [14].
The material model defined for the solid elements has to reproduce the performance of what has been called transverse reinforced resin. This is, a resin with improved properties in its transversal and shear directions, due to the presence of fibers. The transversal directions are defined as those directions perpendicular to the fiber longitudinal orientation.
As stated by Rastellini et al. when defining the serial parallel mixing theory [17] and used by Martinez et al. in the simplified model presented in [16], the contribution of the fibers to the composite in the transversal and shear directions can be captured by applying the inverse mixing theory. This states that the constitutive matrix of the composite () can be obtained by the inverse combination of the stiffness matrix of its constituents ( and ), proportionally to their volumetric participation ( and ), being and matrix and fiber, respectively. This is
|
(21) |
The constitutive matrices and are shown in Eq.(22)
|
(22) |
with
|
(23) |
We note that the constitutive matrix accounts for the effect of the fiber matrix where has a null value in the fiber direction (in practice a small value is chosen to avoid numeric instability). The stiffness in the fiber direction is provided by the RFB elements.
The distinct form of matrices and is needed in order to accurately reproduce the shear stiffness in the resulting matrix . In this manner, the combination of solid elements and RFB elements can reproduce a parallel mixed material behavior along the longitudinal direction of the RFB elements (deformations in the fiber direction are compatible between all materials) and a serial mixed material behavior in the orthogonal direction to the fibers (stresses in the normal fiber direction are compatible between all materials), according to the inverse mixing theory [17]. An advantage of this approach versus the standard serial-parallel material mixing model is the compatibility between the stresses in the concrete and the fibers in the orthogonal direction to the fibers, which facilitates the convergence of the nonlinear iterative process.
Non-linear effects in the response of the material are modelled by a simple isotropic damage model [18,19]. According to this, Eq.(21) is modified as
|
(24) |
where is the damaged constitutive matrix for the solid element and the damage parameters are an scalar numbers between 0 and 1 defined as
|
(25) |
In Eq.(25) is the elastic limit, is the maximum principal stress, is the fracture energy, is the Young's modulus and index () refers to the matrix or fiber material, respectively. In addition, is the characteristic length of the element.
The above constitutive model can accurately reproduce axial and pure bending stress states in which the RFB elements are mainly subjected to tensile forces. The model has proved to work well in shear-dominated regions in which two different stress situations appear in the shear zone and in the region where the loads are applied. In the shear zone the matrix material is soon fully damaged (i.e. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat{\boldsymbol D}\to 0} ) and the stiffness is mainly contributed by the fibers as they rotate and change their orientation. In zones where the loads are applied in a direction orthogonal to the fibers, the fibers are agglomerated in packages that allow the transmission of compression stresses to the resin material, inducing damage. The contribution of matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \hat{\boldsymbol D} }
is important in these zones. This distinct phenomenon has been modelled by limiting the maximum value of the damage parameters in Eq.(25) in these regions as
|
(26) |
The elimination of the rotations as nodal degrees of freedom in bending-type elements (such as beams, plates and shells) leads to the so-called rotation-free elements [20,21,22]. Rotations can be simply eliminated by expressing the curvature field in a node, or a point within an element, in terms of the displacements of nodes contained in an appropriate patch of elements surrounding the node or the point. This procedure emanates from the efforts of finite difference practitioners for analysis of thin plates based on Kirchhoff plate theory using displacement variables only [20]. The idea evolved in the last decades of the 20th century and lead to the formulation of different families of rotation-free elements for analysis of plates and shells. Oñate and Zárate [20] developed a unified formulation for this type of rotation-free plate and shell elements combining ideas from the finite volume method with the FEM. Some of the rotation-free shell elements have enjoyed much popularity for non-linear analysis of shells structures under dynamic and impact loads and sheet stamping processes using explicit and implicit dynamic formulations [23].
The use of rotation-free beam elements has been less popular in practice. The derivation of this type of elements based on Euler-Bernoulli beam theory can be found in [14,22,24]. Extensions of these elements to account for shear deformation effect was reported by Zarate and Oñate [21].
A distinct feature of all rotation-free elements is that the discretized equilibrium equations can be expressed in terms of displacement degrees of freedom only. This makes these elements advantageous for their coupling with solid elements, that are also formulated in terms of nodal displacement DOFs, in a straightforward manner. In this work we will model the mechanical behavior of the fibers with rotation-free Euler-Bernoulli beam elements formulated using an updated Lagrangian description and a cell-vertex approach [14,20]. We present below details of the element formulation.
A particular class of rotation-free beam element can be derived by computing the curvature at each node using a finite difference scheme. The resulting element is termed CVB (for Cell-Vertex Beam element) [14,20].
For the sake of simplicity we will show the derivation of the stiffness equation of a rotation-free straight beam element that bends in the plane (Figure 3). A similar process can be followed for deriving the stiffness equations for bending in the plane. Due to the particular one-dimensional (1D) features of the rebars considered in this work, we will neglect torsional effects in the fibers.
![]() |
Figure 3. Modelling strategy for analysis of FRP with the FEM |
Let us formulate the deformation of a straight beam element for a current configuration at time . The bending behavior of the beam is modelled by considering the control domain formed by half of the lengths of the elements adjacent to a node. The curvature at node is computed as
|
(27) |
where and subscripts and denote the midpoints of the elements located at the right and left of node (Fig. 4). The curvature is assumed to be constant in the control domain assigned to node (Figure 4).
![]() |
Figure 4. CVB element. Control domains around two nodes and |
The rotations and are expressed in terms of the nodal deflections as
|
(28) |
Substituting Eqs.(28) into (27) gives
|
(29) |
where
|
(30) |
Similarly, the curvature at node is found as
|
(31) |
with
|
(32) |
The internal virtual work over the element is obtained by adding up the contributions from the two control domains and as
|
(33) |
Substituting Eqs.(27) and (31) into (33) and following the standard assembly process of the FEM yields the element stiffness matrix as
|
(34.a) |
where
|
(34.b) |
In Eqs.(33) and (34) and are the Young modulus of the fiber material and the moment of inertia of the fiber section with respect the bending axis, respectively. In the computation of the above integrals, we have assumed that and are constant over the control domains.
Note that, as it is usual in rotation-free beam, plate and shell elements, the stiffness matrix of an element involves the nodal deflections of the adjacent elements.
The global equivalent nodal force is computed
|
(35.a) |
with
|
(35.b) |
where is the external nodal force acting at node and and and are uniformly distributed loads acting in the vertical direction over elements and , respectively.
The prescribed values for the deflection are imposed when solving the global system of equations, as usual. We briefly explain next how to compute the curvature at a free end and at simply supported and clamped nodes. For details see [14].
The curvature at a free end node and at a simply supported node is zero. This condition is implemented by neglecting the contribution of the boundary node to the element stiffness matrix. For a boundary element with a free or simply supported node the stiffness matrix is (Figure 5a)
|
(36.a) |
|
(36.b) |
![]() |
Figure 5. CVB element. (a) Boundary condition of zero curvature at a free or simply supported node . (b) Control domain for a clamped or symmetry node |
The curvature at a clamped or symmetry left-end node is computed as (Figure 5b)
|
(37) |
where is an auxiliary fictitious deflection.
Similarly, for a right-end clamped or symmetry node (Figure 5b)
|
(38) |
where is an auxiliary fictitious deflection.
The element stiffness matrix is computed in both cases by Eq.(36) with or equal to zero, as appropriate. The fictitious nodal deflections or (and indeed ) are prescribed to a zero value when solving the global system of equations.
The above rotation-free formulation can be extended to account for the bending behavior of the beam in the plane in a straightforward manner. The nodal degrees of freedom are the vertical displacement along the and axes [14].
The extension of the bending formulation of rotation-free beam element to account for transverse shear deformation effect is possible. The details can be found in [21].
The material behavior of the rotation-free CVB element is simply characterized by the Young modulus. Accounting for non-linear effects in the fiber material under axial stresses has been modelled with a standard one-dimensional elastic-damage constitutive law, i.e.
|
where and are the Young modulus of the fiber in the undamaged and damaged states, respectively.
The evolution of the damage parameter for the beam element is governed the maximum deformation of the fiber at failure, i.e.
|
where is the axial deformation of the fiber, is the deformation at the elastic limit, , and is the deformation at the ultimate (fracture) strength value.
A curved rebar will be modelled by a collection of straight rotation-free CVB elements. Figure 6 shows a curved rebars under bending loads in a plane , discretized with three CVB elements.
![]() |
Figure 6. Cell-vertex rotation-free beam element. Control domains for computing the average curvature. A, B and C are element mid-points |
The stiffness matrix is obtained as the sum of the axial stiffness matrix () and the bending stiffness matrix ().
The axial strain matrix is obtained as
|
(39) |
with and denotes variables in the local coordinates system.
The axial stiffness matrix is computed as
|
(40) |
The procedure for computing the bending stiffness matrix is an extension of that followed for straight CVB elements in Section 4.2.
The curvature at node is computed in terms of the gradients of the normal deflection at points and (Figure 6)
|
(41.a) |
In Eqs.(41.a) are the nodal deflections in the local axes of the element defined in Figure 7 and
|
(41.b) |
and
|
(41.c) |
![]() |
Figure 7. Local coordinate system of curved rotation-free CVB element bending in the plane |
Similarly, for node
|
(42.a) |
with
|
(42.b) |
and
|
(42.c) |
The bending stiffness matrix for the CVB element is expressed as the sum of the stiffness contributions from the two control subdomains 1 and 2 that form the element (Figure 6), i.e.
|
(43) |
The total stiffness matrix for the element is finally obtained by
|
(44) |
where is given by Eq.(40).
We note that the size of for the CVB element bending in the plane is , as it involves the two DOFs of nodes and . A similar procedure will be followed for the case of bending in the plane [14].
The assembly process follows the general rule of the FEM [14].
Figure 8 shows a conceptual representation of a FRP bar of rectangular section containing of a number of fibers embedded in the (bulk) matrix material. The bar geometry is discretized into 3D solid elements (8-noded hexahedral are shown). The fibers contained within each solid element are assumed to be concentrated in one-dimensional (1D) packs of fibers distributed within the 3D solid elements. Figure 8 shows the concentration of fibers in four 1D packs placed at the Gauss points locations in each cross section [14]. Each fiber pack is modelled as a rotation-free beam element with a cross-section area equivalent to the amount of fibers considered in the pack.
The displacements of the rotation-free beam element nodes are related to the nodal displacement of the 3D solid element by a simple interpolation.
![]() |
Figure 8. Patch of 3D solid and rotation-free beam elements |
The discretized equilibrium equations for the coupled solid-beam system can be written as explained in Section 2. The numerical response of the FRP bars for different incremental load up to failure is computed using the standard Newton-Raphson iterative procedure described in Eqs.(4)-(13).
The three examples shown hereafter demonstrate the good behavior of the formulation described for analysis of FRP rebars. The examples include the analysis of a FRP rebar in three standard laboratory tests: an axial test, a three-point bending test and a shear test. The numerical results are compared with results from laboratory tests carried out at the experimental lab for structural materials of the Centro de Investigación en Materiales Estructurales (CIME) of the Technical University of Madrid (UPM), using rebars made of ECR-Glass and a matrix of DERAKANE8084EpoxyVinylester.
The rebar nominal diameter is 10mm, including helicoidal ribs. As these ribs are manufactured by machining the originally pultruded rebars, the effective diameter of the bars is reduced to 9.58 mm, on average. This diameter is the one adopted for creating the model of the bar geometries.
The mechanical parameters of the fibers are: Young modulus: 74.61 GPa, Poisson ratio: 0.22, Elastic limit (tension): 1,653 MPa, Elastic limit (compression): 1,322 MPa, Fracture energy: 126,778 J/m2, Deformation at failure: 0.0277, Volume fraction: 70%
The mechanical parameters of the matrix material are: Young modulus: 3.17 GPa, Poisson coefficient: 0.38, Elastic limit: 72 MPa, Fracture energy: 45,657 J/m2, Deformation at failure: 0.11, Volume fraction 30%
Figure 9 shows the bar modelled with 960 8-noded hexahedral elements and 1197 nodes. The fibers have been modelled using 1140 RFB elements. The length of the sample is 100 mm long with 9.58 mm diameter. Symmetry conditions have been considered by prescribing a zero axial displacement in the normal direction to the surface at a bar end, while the other end has a prescribed axial displacement condition to induce tensile stresses in the bar. The bar axis has the displacements in the normal plane yz constrained to zero.
![]() |
Figure 9. Tensile test geometry: 3D solid and RFB elements mesh |
Figure 10 shows a typical bar tested at UPM after failure.
![]() |
Figure 10. Tensile test of the PFR performed at UPM |
The numerical results obtained are shown in Figures 11 and 12. The displacement field obtained presents a linear distribution, until the rebar reaches brittle failure.
![]() |
Figure 11. Displacement field [m] in deformed geometry |
![]() |
Figure 12. Axial tension tests stress-strain curves. Numerical and experimental values |
The numerical load-displacement curve obtained is converted into a stress-strain curve to be compared with the experimental results. A good correlation can be seen between the numerical and experimental values (Figure 12).
Table 1 shows the good correlation between the numerical and experimental stiffness and strength magnitudes.
Property | Lab tests | Numerical |
---|---|---|
Tensile strength (MPa) | 1150 34 (MPa) | 1173.33 (MPa) |
Strain to max strength | 2.3 0.1 (%) | 2.22 (%) |
Young Modulus (GPa) | 55 4 (GPa) | 53.03 (GPa) |
The experiment performed for measuring the bending properties, corresponds to the standard three-point bending test showed in Figure 14. The tests were carried out on bars of full circular section, instead of using cut specimens as indicated in the ASTMD-4476 norm. Figure 13 shows the bar modelled with 5400 8-noded hexahedral elements (6327 nodes) and 6156 rotation-free beam elements. The length of the bar sample is 100 mm long with 9.58 mm diameter. Symmetry and boundary conditions have been considered, as shown in Figure 13. To avoid problems with the stress concentration at the loading points, a small part of the beam (the blue zone in Figure 13) has been defined with an unbreakable material. The elastic properties of this material are the same as for the rest of the bar.
![]() |
Figure 13. Three point bending test for PFR rebar. Dimensions, FEM mesh and boundary conditions |
Figure 14 shows a snapshot of the laboratory test for the FRP rebar in a three-point bending configuration. Note the rollers used as a pin-supports of the beam.
![]() |
Figure 14. Three-point bending test of FRP rebars performed at UPM |
Figure 15 shows the strain distribution over the deformed bar (Figure 15a) and the mesh of hexahedra (Figure 15b) in a load state close to the onset of damage. The damage distribution for each element is displayed in Figure 16.
![]() |
Figure 15. Three-point bending tests. Numerical results. Deformed bar. (a) Strain distribution for RFB elements. (b) Strain distribution in 3D solid elements |
![]() |
Figure 16. Contour field of damage variable in longitudinal and transverse sections on the deformed geometry (1x) |
In Figure 17, the numerical load-deflection curve is shown together with that obtained in the laboratory for comparison purposes. The maximum loading force and the flexural stiffness captured by the numerical simulation are in good agreement with the experimental results.
![]() |
Figure 17. Three-point bending tests. Force-deflection curves for shearing testing. Numerical results compared with laboratory tests |
At the laboratory, samples of 240mm of total length were cut for the 10mm diameter full cylinder tests. The span of the tests was set to 200mm. Tests were carried out under displacements control, with a crosshead displacement velocity of 3mm/min. Figure 18 shows the test set-up and one specimen after failure.
![]() |
Figure 18. Three-point bending tests. Failed specimen of PFR rebar |
The experiment carried out for evaluating the shear properties, corresponds to the ASTM-7617 standard test showed in Figure 19 [25]. For the numerical analysis, the symmetry of the geometry and the imposed displacement have been taken into account as shown in Figure 20a in order to speed up the calculations. Figure 20b displays the bar modelled with 5225 8-noded hexahedral elements (4608 nodes) and 5016 rotation-free beam elements. The length of the sample is 100 mm long with 9.58 mm diameter.
| ||||
Figure 19. Shear test. Double shear fixture according to ASTM D-7617 standard [25] |
| ||||
Figure 20. Shear test geometry. (a) Boundary conditions. (b) Mesh and dimensions |
In Figure 20 a gap of 0.001 m can be observed between the loaded and supported zones of the specimen. This zone is necessary to capture accurately the shear forces and to avoid stress concentration in the numerical model.
Figure 21 shows a representative configuration of the numerical results. Figure 21a displays the deformed bar and the localized damage in the resin material at the shear zone. Figure 21b is a close-up of the shear zone were the RFB are presented. The curvature of these elements due to the shear deformation is depicted.
Figure 22 shows the numerical load-deflection curve superimposed to a set of curves obtained in the laboratory tests for comparison purposes. In the curves the non-linear load phenomena of the fiber reorientation can be clearly observed. For an applied load slightly larger than 5kN there is an initial plateau in the experimental tests, also reproduced by the numerical model, that corresponds to the cracking of the matrix. At this point there is a fiber reorientation that provides the stiffness increase and the final composite failure, due to fiber failure, for a load close to 35kN. The good agreement obtained between the numerical and the experimental results, especially in this case, highly non-linear, proves the validity of the proposed procedure to simulate pultruded fiber reinforced polymer rebars.
![]() |
Figure 21. Shear test. Contour field of damage variable in the matrix. Deformed geometry (1x) |
![]() |
Figure 22. Shear test. Force-deflection curves. Numerical results compared with laboratory tests |
The combination of 3D solid finite elements with rotation-free beam elements has shown to be adequate for modelling axial, bending and shear deformation states in fiber reinforced polymer rebars. The fact that both types of elements are formulated in terms of displacement DOFs only simplifies the combination of the best features of both elements for modelling the 3D behavior of the matrix material and the 1D behavior of the fibers, as well as the coupling interactions during complex deformation modes.
The good agreement of the numerical results for several benchmark laboratory tests with experimental data for the same tests have, validated the accuracy of the combined 3D-1D model presented. The combined model can be therefore be used confidently in practical applications of the fiber polymer rebars as reinforcement of concrete structures.
The approach presented in this work for combining 3D solid elements with rotation-free beam elements can be implemented in more sophisticated formulations for modelling fiber reinforced materials, such as that proposed in [26].
Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request.
This research was carried out as a part of the project “Development of polymer fiber reinforced rebars” funded by Saudi Aramco. The authors also acknowledge the financial support from the Spanish Ministry of Economy and Competitiveness, through the “Severo Ochoa Programme for Centres of Excellence in R&D” (CEX2018-000797-S).
[1] Pacheco-Torgal F., Jalali S., Labrincha J., John V.M. Eco-efficient concrete. Woodhead Publishing Series in Civil and Structural Engineering, pp. 624, 2013.
[2] Bank L.C. Composites for construction: Structural design with FRP materials. John Wiley & Sons, Inc., pp. 560, 2006.
[3] Marí A.R., Oller E., Bairán J.M. Predicting the response of FRP-strengthened reinforced-concrete flexural members with nonlinear evolutive analysis models. J. Compos. Constr., 15(5):799–809, 2011.
[4] Martínez X., Oller S., Rastellini F., Barbat A.H. A numerical procedure for simulating RC structures reinforced with FRP using the serial/parallel mixing theory. Computers & Structures, 86(15-16):1604-1618, 2008. doi:10.1016/j.compstruc.2008.01.007
[5] Nogales A., Tosić N., de la Fuente A. Rotation and moment redistribution capacity of fiber-reinforced concrete beams: Parametric analysis and code compliance. Struct. Concr., 23(1):220-239, 2021.
[6] Pultron Composites [Online]. Available: https://www.pultron.com/products/features/. [Accessed: 17-Apr-2020].
[7] Lleida GFRP Pedestrian Bridge [Online]. Available: http://www.pedelta.com/lleida-gfrp-pedestrian-bridge-p-52-en. [Accessed: 17-Apr-2020].
[8] Hernández S., et al. Design and development of an innovative structural system based on composites applied on the constrction of maritme lighthouses with lower maintenance requeriments (FAROCOMP). Mater. Compuestos, 3(2):54–58, 2019.
[9] GangaRao H.V.S., Taly N., Vijay P.V. Reinforced concrete design with FRP composites. 1st Edition, CRC Press, pp. 400, 2019.
[10] Sharifianjazi F., Zeydi P., Bazli M., Esmaeilkhanian A., Rahmani R., Bazli L., Khaksar S. Fibre-reinforced polymer reinforced concrete members under elevated temperatures: A review on structural performance. Polymers, 14(3):472, 2022. 10.3390/polym14030472
[11] Sagar B., Sivakumar M.V.N. Performance evaluation of basalt fibre-reinforced polymer rebars in structural concrete members–a review. Innovative Infrastructure Solutions, 6(2):1-18, 2021. 10.1007/s41062-020-00452-2
[12] Strau S., Senz A., Ellinger J. Comparison of the processing of epoxy resins in pultrusion with open bath impregnation and closed-injection pultrusion. J. Compos. Sci., 3(3):87, 2019.
[13] Starr T.F. Pultrusion for engineers. Woodhead Publishing Series in Composites Science and Engineering, pp. 336, 2000.
[14] Oñate E. Structural analysis with the finite element method. Linear statics. Vol 1: Basis and solids, 2009. Vol. 2: Beams, plates and shells, 2013, Springer-CIMNE.
[15] Zienkiewicz O.C., Taylor R.L. The finite element method for solids and structural mechanics. 6th edition, Elsevier, 2005.
[16] Martínez X., Rastellini F., Oller S., Flores F., Oñate E. Computationally optimized formulation for the simulation of composite materials and delamination failures. Composites Part B: Engineering, 42(2):134-144, 2011. 10.1016/j.compositesb.2010.09.013
[17] Rastellini F., Oller O., Salomón O., Oñate E. Composite materials non-linear modelling for long fibre-reinforced laminates: Continuum basis, computational aspects and validations. Computers & Structures, 86(9):879-896, 2008. 10.1016/j.compstruc.2007.04.009
[18] Oller S., Oñate E., Oliver J., Lubliner J. Finite element non-linear analysis of concrete structures using a plastic-damage model. Int. Journal of Engineering Fracture Mechanics, 35(1/2/3):219-231, 1990. 10.1016/0013-7944(90)90200-Z
[19] Lubliner J., Oller S., Oliver J., Oñate E. A plastic damage model for concrete. Int. Journal of Solids and Structures, 25(3):299-326, 1989. 10.1016/0020-7683(89)90050-4
[20] Oñate E., Zárate F. Rotation-free triangular plate and shell elements. Int. J. Numer. Meth. Engng., 47(1-3):557-603, 2000. <557::AID-NME784>3.0.CO;2-9 10.1002/(SICI)1097-0207(20000110/30)47:1/3<557::AID-NME784>3.0.CO;2-9
[21] Oñate E., Zárate F. Extended rotation-free plate and beam elements with shear deformation effects. Int. J. Numer. Meth. Engng., 83(2):196-227, 2010. 10.1002/nme.2836
[22] Jovicevic J., Oñate E. Analysis of beams and shells using a rotation - free finite element - finite volume formulation. Monograph M43, CIMNE, ISBN: 84-89925-36-4, Barcelona, 1999.
[23] Oñate E., Flores F.G. Advances in the formulation of the rotation-free basic shell triangle. Computer Methods in Applied Mechanics and Engineering, 194(21-24):2406-2443, 2005. 10.1016/j.cma.2004.07.039
[24] Flores F., Oñate E. Rotation-free finite element for the non-linear analysis of beams and axisymmetric shells. Computer Methods in Applied Mechanics and Engineering, 195(41-43):5297-5315, 2006. 10.1016/j.cma.2005.08.021
[25] ASTM-International, ASTM D7617/D7617M - 11(2017) Standard Test Method for Transverse Shear Strength of Fiber-reinforced Polymer Matrix Composite Bars, West Conshohocken, PA, 2017.
[26] Khristenko U., Schuß S., Krüger M., Schmidt F., Wohlmuth B., Hesch C. Multidimensional coupling: A variationally consistent approach to fiber-reinforced materials. Computer Methods in Applied Mechanics and Engineering, 382:113869, 2021. 10.1016/j.cma.2021.113869
Published on 04/05/23
Accepted on 04/05/23
Submitted on 16/03/23
Volume 39, Issue 2, 2023
DOI: 10.23967/j.rimni.2023.05.001
Licence: CC BY-NC-SA license
Are you one of the authors of this document?