Improved numerical analysis of structures reinforced by composite FRP: Hygrothermal and prestressing loads with taper effect

Abstract

Fiber-reinforced polymer (FRP) composites are becoming suitable and substantial materials in repairing and replacing conventional metallic materials because of their high strength and stiffness. Steel beams can be strengthened in flexure using bonded FRP or using steel plates. In such plated beams, shear forces develop in the bonded beam and these will be transferred to the FRP plate via the adhesion technique. Thus, the interfacial shear stress and normal stress will develop consequently, and debonding may occur at the FRP plate ends due to high interfacial stress values in this area. This original research aims to study the debonding phenomenon using an analytical and a numerical finite element models, in order to identify the interfacial stresses of a steel beam strengthened by the FRP plate with taper model, taking into account a new coupled approach of prestressing force and hygrothermal effect. This article explores the effects of various parameters, such as geometrical and physical properties, on the stress behavior of FRP composites.

Keywords

FRP interfacial stress hygrothermal prestressing force taper effect

Introduction

Composite material made from two or more immiscible components together produce material with properties that exceed the individual components. Experts and researchers determined various solutions to deal with the crucial problem of structures aging caused by different loads (mechanical load, temperature, humidity, hygrothermal effect, corrosion, etc.), and they have been using new composite materials, fiber-reinforced polymers (FRPs), which have several practical applications. The FRP composites have been mainly used to reinforce concrete and steel beams; extensive research has shown that the bonding between an FRP plate to the face of concrete and the steel beam can effectively develop its capability force and ultimate strength; these types of reinforced models have been able to reduce the influence of the interfacial stresses effect, which are affected by the phenomenon of solicitations, and the success of this technique is due to the effect of the transferred interfacial stresses from the beam to the externally bonding FRP. Several works have been developed by Wu et al.,¹ Roberts,² Saadatmanesh and Malek,³ and Triantafillou and Antonopoulos⁴ to study the debonding phenomenon, while they have omitted several parameters in their studies. Smith and Teng⁵ presented an accurate and widely applicable solution, especially they have taken the flexural stiffness of the bonded beam and plate. Lau et al.,⁶ presented a simple theoretical model for estimating interfacial stresses, in which fiber orientations are considered in their analysis. However, this method does not take into account the effects of bending deformations in the plate and the thermal effect on the structures. In addition, an improved solution for calculating interfacial stresses in concrete beams reinforced with FRP plates is developed by Tounsi,⁷ in which shear deformations are considered, but the FRP fiber orientations are ignored. The solution developed by Kerboua et al.⁸ supports the shear deformation effects in their model under the thermal load using typical coefficients of thermal expansion. In some studies, the difference between the coefficients of thermal expansion is another important problem in steel beams reinforced with FRP plates. Denton⁹ found that the interfacial stresses can be developed in the adhesive layer at the ends of the FRP composite as a result of the difference between the coefficients of thermal expansion. Deng et al.¹⁰ presented a simple method for obtaining the interfacial shear and normal stress using a taper model in the composite plates under the thermal and mechanical loads, while Hart-Smith¹¹ proposed the use of lap joints to minimize stress concentrations. The use of the taper at the ends of the adherend has been shown to effectively limit delamination stresses that has been developed by Vinson and Sierakowski,¹² Hart-Smith,¹³ and Belabed et al.^14,15 The hygrothermal effect on the interfacial stresses for structures strengthened with FRP plates was presented by Philip and De Vries¹⁶ for a monolithic porous medium. In addition, there are several works to evaluate the couple (FRP-steel) and (FRP-concrete), which were experimentally investigated under the shear tensile and combined loading, such as Jiang et al.,¹⁷ Jiang et al.,¹⁸ Mendes et al.,¹⁹ Lin et al.,²⁰ and Jones et al.,²¹ the different experimental bonding studies give appreciable results. In recent applications, the use of prestressed laminates has been developed in several laboratories with tests on different materials and on various structures, such as Nordin,²² Bassetti,²³ Luke,²⁴ and De Lorenzis et al.²⁵ This technique of prestressing laminates can be exploited well, which gives a more effective strengthening, according to Benachour et al.²⁶ The problem of interfacial shear and normal stresses when prestressed laminates are used has treated by Benachour et al.,²⁶ Al-Emrani and Kliger,²⁷ and Brairi et al.²⁸ This original research aims to study the debonding phenomenon using an analytical model and a numerical finite element model, developed under ABAQUS software version 6.145 to identify the interfacial stresses of a steel beam strengthened by prestressed FRP plates, taking into account a new approach of prestressing force and hygrothermal effect with taper model. The numerical and analytical results are compared with the results developed by Deng et al.,¹⁰ This model gave a good agreement with those obtained by the experimental results done by Jones et al.²¹ After verifying the accuracy of the present work, a parametric study is used to identify the effects of various geometrical and material properties on the magnitude of stresses concentration. The present work gives an assessment of analytical and numerical models and a wide database including several experimental results for the evaluation of the bond strength between FRP plate composite and structures.

Basic assumptions

Figure 1 shows the geometry of the beam reinforced with FRP plate composite, and the configuration of the interface of the adhesive layer and the applied load are shown in Figure 2. The terms V, M, and N present, respectively, the shear force, the bending moment, and the longitudinal tension, while the beam and the FRP plate are symbolized, respectively, by 1 and 2, τ (x) is the shear stress, and (t) is the component thickness. This analytical and numerical model is madding by the following assumptions:

A linear elastic material is considered.

The shear stress in the adhesive layer does not vary through the thickness of the adhesive.

Shear deformations in the FRP plate and steel beam are neglected.

The deformations of the adhesive are neglected when the plate terminates further away from the support.

Hygrothermal and thermal load are taken into account.

The prestressed force model is considered.

Figure 1.

Geometry of structure.

Figure 2.

Distribution of forces.

Analytical equations

Shear stress formulation

From Figure 2, G is the shear modulus of the adhesive, u (x, y) and v (x, y), are, respectively, the longitudinal and the transverse displacement at any point in the adhesive layer. The corresponding shear stress is given by equation (1):

τ = - G (\frac{d u_{(x, y)}}{d_{y}} + \frac{d v_{(x, y)}}{d x})

Differentiating equation (1) gives equation (2):

\frac{d τ}{d x} = - G (\frac{d^{2} u_{(x, y)}}{d_{y} d_{x}} + \frac{d^{2} v_{(x, y)}}{d x^{2}})

The shear and the normal stresses in the adhesive layer do not vary through the thickness of the adhesive and are given by equations (3) and (4):

\frac{d u}{d y} = \frac{1}{t_{a}} (u_{1} - u_{2})

\frac{d^{2} u}{d y d x} = \frac{1}{t_{a}} (\frac{d u_{1}}{d x} - \frac{d u_{2}}{d x})

where t _a is the thickness of the adhesive layer and (u ₁, u ₂) present the longitudinal displacements at the bottom of the steel beam and at the top of the FRP plate, respectively. Integration of thermal and hygrothermal model is given by equations (5) and (6):

α_{1} Δ T + \frac{M_{1 t_{1} / 2}}{E_{1} I_{1}} + \frac{N_{1}}{E_{1} A_{1}} + β_{1} Δ c = ε_{1} (x)

α_{2} Δ T + \frac{M_{2 t_{2} / 2}}{E_{2} I_{2}} + \frac{N_{2}}{E_{2} A_{2}} + β_{2} Δ c - \frac{P_{0}}{E_{2} A_{2}} = ε_{2} (x)

where ε is the strain, β₁ is the humidity coefficient, α₁ is the thermal expansion coefficient, Δc is the humidity change, ΔT is the temperature change, P ₀ is the prestressing force, and A ₁ and A ₂ represent area. In calculating the interfacial shear stress, bending of the FRP plate is ignored.

Substituting equations (5) and (6) into equation (4) gives equation (7):

\frac{d^{2} u}{d y d x} = \frac{1}{t_{a}} (Δ T (α_{1} - α_{2}) + Δ C (β_{1} - β_{2}) + \frac{P_{0}}{E_{2} A_{2}} + \frac{M_{1 t_{1} / 2}}{E_{1} I_{1}} + \frac{N_{1}}{E_{1} A_{1}} - \frac{N_{2}}{E_{2} A_{2}})

Integrating equations (2) and (7) gives equation (8):

\frac{d^{2} u}{d y d x} = \frac{1}{t_{a}} (Δ T (α_{1} - α_{2}) + Δ C (β_{1} - β_{2}) + \frac{P_{0}}{E_{2} A_{2}} + \frac{(M (x) + N_{1} (t_{1} / 2 + t_{2} / 2) \frac{t_{1}}{2})}{E_{1} I_{1}} + N_{b} (\frac{1}{E_{1} A_{1}} + \frac{1}{E_{2} A_{2}}))

Substituting the longitudinal force from (N1 = −N2) gives equation (9):

\frac{d τ}{d x} = \frac{G}{t_{a}} (Δ T (α_{1} - α_{2}) + Δ C (β_{1} - β_{2}) + \frac{P_{0}}{E_{2} A_{2}} + \frac{(M (x) + N_{1} (t_{1} / 2 + t_{2} / 2) \frac{t_{1}}{2})}{E_{1} I_{1}} + N_{1} (\frac{1}{E_{1} A_{1}} + \frac{1}{E_{2} A_{2}}))

The equilibrium of the steel beam in the x direction can be expressed by equation (10):

\frac{1}{b} \frac{d N_{b}}{d x} = - τ

Differentiating equation (9) and substituting equation (10) gives the governing equation (11):

\frac{d^{2} τ}{d x^{2}} = λ^{2} τ - g λ^{2} M^{,} (x)

The solution of this equation has the general solution given by equation (12):

τ (x) = B_{1} e^{λ x} + B_{2} e^{- λ x} + g M^{,} (x)

and $λ = \sqrt{\frac{G 1}{t a} (\frac{(t_{1} / 2 + t_{2} / 2) t_{1} / 2}{E_{1} I_{1}} + \frac{1}{E_{1} A_{1}} + \frac{1}{E_{2} A_{2}})}$ , $g = \frac{G t_{1} / 2}{λ^{2} t_{a}} \frac{1}{E_{1} I_{1}}$ .

Considering the boundary condition due to symmetry, the shear stress at mid-span is zero:

τ (L_{p} / 2) = B_{1} e^{λ L_{2} / 2} + B_{2} e^{- λ L_{2} / 2} + g M^{,} (L_{2} / 2)

The longitudinal force is zero at the end of the plate, so equation (9) can be written by equation (14):

\frac{d τ (0)}{d x} = B_{1} λ - B_{2} λ + g M^{,,} (0) = - \frac{G}{t_{a}} (Δ T (α_{1} - α_{2}) + Δ C (β_{1} - β_{2}) + \frac{P_{0}}{E_{2} A_{2}} + \frac{M (0) t_{1} / 2}{E_{1} I_{1}})

where B ₁ and B ₂ are, respectively, given by equations (15) and (16):

B_{1} = \frac{\frac{- 1}{λ} e^{- λ L 2 / 2} (\frac{G}{t_{a}} (Δ T (α_{1} - α_{2}) + Δ C (β_{1} - β_{2}) + \frac{P_{0}}{E_{2} A_{2}} + \frac{M (0) t_{1} / 2}{E_{1} I_{1}}) + G M^{,,} (0)) - g M^{,} (\frac{L_{2}}{2})}{e^{\frac{λ L_{2}}{2} +} e^{\frac{- λ L_{2}}{2}}}

B_{2} = - B_{1}

Normal stress formulations

The interfacial normal stress is given by equation (17):

σ (x) = \frac{E_{a}}{t_{a}} (v_{2} (x) - v_{1} (x))

where v1 and v2 are the transverse displacements at the bottom of the steel beam and at the top of the FRP plate, respectively, which are given by equations (18) and (19):

\frac{d^{2} v_{1} (x)}{d x^{2}} = \frac{1}{E_{1} I_{1}} M_{1} (x)

\frac{d^{2} v_{2} (x)}{d x^{2}} = \frac{1}{E_{2} I_{2}} M_{2} (x)

Considering the moment equilibrium of an infinitesimal element of the steel beam:

(M_{1} + d M_{1}) - M_{1} - (V_{1} + d V_{1}) d x + \frac{1}{2} (σ b - q) d x^{2} + τ b (\frac{t_{1}}{2}) d x = 0

\frac{d V_{b}}{d x} = b σ - q

Substituting equation (20) into equation (21) gives equations (22) and (23):

\frac{d M_{1}}{d x} = V_{1} - \frac{τ b t_{1}}{2}

\frac{d M_{2}}{d x} = V_{2} - \frac{τ b t_{2}}{2}

Differentiating equation (17) three times gives equation (24):

\frac{d^{3} σ (x)}{d x^{3}} = \frac{E_{a}}{t_{a}} (\frac{1}{E_{2} I_{2}} (V_{2} - \frac{τ b t_{2}}{2}) - \frac{1}{E_{1} I_{1}} (V_{1} - \frac{τ b t_{1}}{2}))

Differentiating equation (24) gives the governing equation (25):

\frac{d^{4} σ (x)}{d x^{4}} + \frac{E_{a} b}{t_{a}} (\frac{1}{E_{2} I_{2}} + \frac{1}{E_{1} I_{1}}) σ (x) - \frac{E_{a} q}{t_{a} E_{1} I_{1}} - \frac{E_{a} b}{t_{a}} (\frac{\frac{t_{1}}{2}}{E_{1} I_{1}} - \frac{\frac{t_{2}}{2}}{E_{2} I_{2}}) \frac{d τ (x)}{d x} = 0

The general solution to equation (25) is given by equation (26):

σ (x) = e^{- γ x} (D_{1} cos (γ x) + D_{2} sin (γ x)) + m_{1} \frac{d τ}{d x} + m_{3} q

γ = \sqrt[4]{\frac{E_{a} b}{4 t_{a}} (\frac{1}{E_{2} I_{2}} + \frac{1}{E_{1} I_{1}})}

where $m_{1} = \frac{E_{2} I_{2} t_{1} / 2 - E_{1} I_{1} t_{2} / 2}{E_{1} I_{1} + E_{2} I_{2}}, m_{3} = \frac{E_{2} I_{2}}{b (E_{1} I_{1} + E_{2} I_{2})}$

At the end of the FRP plate, the bending moment and shear forces are zero. Hence, the boundary conditions can be written by equations (28) and (29):

\frac{d^{2} σ (x)}{d x^{2}} = \frac{E_{a}}{t_{a}} (\frac{1}{E_{2} I_{2}} M_{2} (0) - \frac{1}{E_{1} I_{1}} M_{1} (0)) = - \frac{E_{a}}{t_{a}} \frac{1}{E_{1} I_{1}} M (0)

\frac{d^{3} σ (0)}{d x^{3}} = m_{2} τ - \frac{E_{a}}{t_{a}} \frac{1}{E_{1} I_{1}} V (0)

and $m_{2} = \frac{E_{a} b}{t_{a}} (\frac{t_{1} / 2}{E_{1} I_{1}} - \frac{t_{2} / 2}{E_{2} I_{2}})$

When x = 0, equation (26) derives to equations (30) and (31):

\frac{d^{2} σ (0)}{d x^{2}} = - 2 D_{2} γ^{2} + m_{1} \frac{d^{3} τ}{d x^{3}}

\frac{d^{3} σ (0)}{d x^{3}} = 2 γ^{3} (D_{1} + D_{2}) + m_{2} \frac{d^{4} τ}{d x^{4}}

Combining equations (28) to (31) gives, respectively, equations (32) and (33):

D_{1} = \frac{1}{2 γ^{3}} (m_{2} τ (0) - m_{1} \frac{d^{4} τ (0)}{d x^{4}}) - \frac{m_{1}}{2 γ^{2}} \frac{d^{3} τ (0)}{d x^{3}} - \frac{E_{a}}{2 γ^{3} t_{a}} \frac{1}{E_{1} I_{1}} (V (0) + γ M (0))

D_{2} = \frac{m_{1}}{2 γ^{2}} \frac{d^{3} τ (0)}{d x^{3}} + \frac{1}{2 γ^{2}} \frac{E_{a}}{t_{a}} \frac{1}{E_{1} I_{1}} M (0)

Geometric and material properties

We consider in this research that the plate length is L _p = 5 m, it is symmetrically linked to its support. The steel beam has a UB122 profile of 533 × 210. A summary of material properties is given in Table 1.

Table 1.

A summary of the material.

Component	E11 (GPa)	E22 (GPa)	E33 (GPa)	υ
Steel beam	210	210	210	0.3
FRP plate	310	10	10	0.3
Adhesive layer	10	10	10	0.3

Comparison with existing analytical solutions and experimental results

The present prediction of the interfacial edge stresses is compared with that developed by Deng et al.¹⁰ for an applied thermal load (ΔT = 50°C) and with the experimental results done by Jones et al.²¹ for an uniformly distributed load (UDL) (q = 60 kN/m²). Figure 3, Figure 4, and Table 2 show an almost exact agreement between results.

Figure 3.

Interfacial shear stress concentration for different solutions (finite element model: ΔT = 50°C).

Figure 4.

Interfacial shear stress distribution: Jones et al.²¹ (q = 60 kN/m²).

Table 2.

Result of model (thermal load).

	Finite-element result (MPa)	Analytical result (MPa)
Deng et al.¹⁰	20.6	21.8
Present work	22.46	21.7

Finite element analysis (application of the taper and hygrothermal model)

The 3-D finite element model of the structural components has been employed to validate the results of the above procedure. The geometry of the reinforced beam is shown in Figure 5, in which only one half of the structure is presented, and the steel beam was modeled with shell elements (element S4 R) in the ABAQUS software (Figure 6), a half beam was considered because of symmetry, this half due to the loading conditions. Figure 7 shows the finite element results, where a prestressing force P ₀ = 100 kN and hygrothermal load (ΔT = 50°C, Relative humidity (RH) = 50%) are applied. Figure 8 plots the interfacial shear and normal stresses for the steel-plated beam, where a linear taper has been used in both analyses.

Figure 5.

Geometry of the plated beam.

Figure 6.

Finite element mesh of geometry.

Figure 7.

Stresses concentration (stage 4).

Figure 8.

Interfacial stress distribution (stage 4).

Different applied load

In this case, numerical results are presented to study the effect of the linear taper model on the distribution of interfacial stresses in beam strengthened with the composite FRP plate. Figure 9 shows the plots of the interfacial stresses for the following applied:

Stage 1: Prestressing force (P ₀ = 100 kN).

Stage 2: Thermal load (ΔT = 50°C).

Stage 3: Hygrothermal load (ΔT = 50°C, RH = 50%).

Stage 4: Prestressing force (P ₀ = 100 kN), hygrothermal load (ΔT = 50°C, RH = 50%).

Figure 9.

Interfacial shear stress distribution for different stages.

Parametric study

After verifying the accuracy of the present analytical and numerical model, a parameter study is carried out to better understand the effects of various parameters on the interfacial stresses under the coupling of different applied load (stage 4).

Effect of hygrothermal load

The effect of the hygrothermal load on the distribution of interfacial stresses at the end of the laminate is illustrated in Figures 10 and 11, and we notice that the higher values of shear and normal stresses correspond with the higher values of relative humidity and the temperature. Tables 3 presents the maximum values of interfacial stresses for the different values of RH and T for stage 4.

Figure 10.

Interfacial shear stress distribution: stage 4 (ΔT constant).

Figure 11.

Interfacial shear stress distribution: stage 4 (RH constant).

Table 3.

Hygrothermal effect on interfacial stresses (stage 4).

ΔT = 50°C	RH = 40%	RH = 50%	RH = 65%	RH = 80%
τ (MPa)	33.3	34.85	37.18	39.51
RH = 50%	ΔT = 25°C	ΔT = 40°C	ΔT = 50°C	ΔT = 60°C
τ (MPa)	24.01	30.51	34.85	39.19

Linear taper effect

The geometry of the taper is shown in Figure 12, where l is the length of the taper, t _end is the thickness at the end of the taper, and t _p is the plate thickness.

Figure 12.

Geometry of the linear taper.

The parameters of the taper (l) and (t _end) are very important. Figures 13 and 14 plot the interfacial shear stress distribution at the tapered end of the plate versus the length of the taper for different tapered end thicknesses. They indicate the beneficial effect of having a thin tapered end and a long taper, the benefit appears when the length of the taper is beyond 500 mm.

Figure 13.

Variation of the maximum shear stress (stage 4).

Figure 14.

Variation of the maximum normal stress (stage 4).

Effect of the prestressing force (P ₀)

In this case, numerical results of the present solution are presented to study the effect of the prestressing force P ₀ on the distribution of interfacial stress in a steel beam strengthened with bonded prestressed FRP plate. Four values of P ₀ are considered in this study (P ₀ =50, 100, 150, and 250 kN). Figure 15 plots the interfacial shear and normal stress for the steel beam strengthened with bonded prestressed FRP plate. From these results, it is seen that increasing the value of prestressing force P ₀ leads to high stress concentrations.

Figure 15.

Interfacial stresses distribution: prestressing force effect (stage 4).

Plate elasticity modulus effect

The effect of elasticity modulus E ₂ on the distribution of interfacial stresses near the end of the prestressed laminate is plotted in Figure 16. For four values of elasticity modulus (E ₂ = 110 GPa, E ₂ = 210 GPa, E ₂ = 310 GPa, and E ₂ = 410 GPa), the maximum values of shear and normal stresses are τ = 37.02 MPa and σ = 30.1 MPa for an elasticity modulus E ₂ = 410 GPa. So the increasing values of the plate elastics modulus will increase the values of maximum shear and normal stress.

Figure 16.

Interfacial stresses distribution: plate elasticity modulus effect (stage4).

Adhesive thickness effect

Figure 17 shows the effect of the thickness of the adhesive layer on interfacial stresses. It is seen that increasing the thickness of the adhesive layer leads to a significant reduction of interfacial stress. Thus, a thick adhesive layer is recommended.

Figure 17.

Interfacial stresses distribution: adhesive thickness effect (stage4).

Conclusion

A numerical and analytical solution for determining the shear and normal interfacial stresses between a composite FRP plate and beam submitted to prestressing force and hygrothermal load with a linear taper has been presented. In this research, a finite element model has been employed to validate the results comparatively done by Deng et al.¹⁰ and with the experimental results done by Jones et al.²¹ The results obtained from the different analysis is agreed, which demonstrates that the present model is simple and accurate. Some studies of the taper are carried out in this work to better understand the effects of various parameters on the edge interfacial stresses in a structure under different applied load. In this work, our model takes into account the parameters neglected in the previous studies. This work is of paramount importance for designers of materials and structures to optimize their operation and life in the context of sustainable development.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

References

Matsuzaki

Tanabe

. Interface crack propagation in FRP-strengthened concrete structures: non-metallic (FRP) reinforcement for concrete structures. Japan Concrete Inst 1997; 1(3): 19–26.

Roberts

. Approximate analysis of shear and normal stress concentrations in adhesive layer of plated RC beam. Struct Eng London 1989;67(12): 229–233.

Saadatmanesh

Malek

. Design guidelines for flexural strengthening of RC beams with FRP plates. J Compos Construct 1998; 2(4): 158–164.

Triantafillou

Antonopoulos

. Design of concrete flexural members trengthened with FRP. J Compos Construct 2000; 4(4): 198–205.

Smith

Teng

. Interfacial stresses in plate beams eng structure. 2001; 23(7): 857–871.

Lau

Zhou

Woo

. Strengthening and strain monitoring concrete structures using fiber glass composites and FBG sensor. J Mater Sci Res Int 1999; 5(3): 216–221.

Tounsi

. Improved theoretical solution for interfacial stresses in concrete beams strengthened with FRP plate. Inter J Solids Struct 2006; 43: 415–474.

Kerboua

Bensaid

Adda bedia

. New model of the laminate layers. J Adhes Sci Technol 2013; 27(17): 1853–1865.

Denton

SR.

Analysis of stresses developed in FRP plate beams due to thermal effects, Composites in civil engineering. Compos Civil Eng 2001; 1: 523–536.

10.

Deng

Marcus Lee

Stuart Moy

. Stress analysis of steel beams reinforced with a bonded CFRP plate. Compos Struct 2004; 65: 205–215.

11.

Hart-Smith

. Adhesive-bonded single-lap joints. Santa Monica: Douglas Aircraft Co., NASA Langley Contract Report CR-112–236, 1973.

12.

Vinson

Sierakowski

. The behaviour of structures composed of composite materials. Dordrecht: Kluwer, 1987.

13.

Hart-Smith

LJ.

Advances in the analysis and design of adhesive bonded joints in composite aerospace structures. In: SAMPE Process Engineering Series, Vol. 19. Azusa: SAMPE, Society for the Advancement of Material and Process Engineering, 1974, pp. 722–737.

14.

Belabed

Kerboua

Tarfaoui

. New design for reducing interfacial stresses of reinforced structures with FRP plates. Int J Build Pathol Adaptat 2019; 37(2): 196–207.

15.

Belabed

Kerboua

Tarfaoui

. New optimized numerical solution of interfacial stresses in steel strengthened structures with CFRP. Adv Civil Eng Mater 2019; 8(1): 117–133.

16.

Philip

De Vries

. Moisture movement in porous materials under temperature gradients. Trans Am Geophys Union 1957; 38: 222–232.

17.

Jiang

Kolstein

Bijlaard

FSK

. Experimental and numerical study, on mechanical behavior of an adhesively-bonded joint of FRP, steel composite bridge, under shear loading. Compos Struct 2014; 108: 387–399.

18.

Jiang

Kolstein

Bijlaard

FSK

. Experimental investigation on mechanical behaviour of FRP-to-steel adhesively-bonded joint under combined loading-part 2: after hygrothermal aging. Compos Struct 2015; 125: 687–697.

19.

Mendes

Phlippi

. A method for predicting heat and moisture transfer through multi layered walls based on temperature and moisture content gradients. Int J Heat Mass Tran 2005; 48: 37–51.

20.

Lin

Berman

Khoshbakht

. Modeling of moisture migration in an FRP reinforced masonry structure. Build Environ 2006; 41(5): 646–656.

21.

Jones

Swamy

Charif

. Plate separation and anchorage of reinforced concrete beams, strengthened by epoxy, bonded steel plates. Struct Eng 1988; 66(5): 85–94.

22.

Nordin

. Fibre reinforced polymers in civil engineering: flexural strengthening of concrete structures with prestressed near surface mounted CFRP rods. Licentiate Thesis, Lulea University of Technology, Sweden, 2003.

23.

Bassetti

Lamelles precontraintes en fibres de carbone pour le renforcement de ponts rivetes endommages par fatigue. Doctoral thesis No 2440, EPF-Lausanne, Switzerland, 2001.

24.

Luke

. The use of carbon fiber plates for the strengthening of two metallic bridges of an historic nature in the UK. In: Proceedings of the international conference on FR composites in civil engineering, Hong Kong, 12 December 2001.

25.

De Lorenzis

Micelli

Tegola

. Fiber reinforced polymer for structural strengthening: posttensioning of steel silos. Struct Eng Int 2003; 13(2): 124–127.

26.

Benachour

Benyoucef

Tounsi

. Interfacial stress analysis of steel beams reinforced with bonded prestressed FRP plate. Eng Struct 2008; 30(11): 3305–3315.

27.

Al-Emrani

Kliger

. Analysis of interfacial shear stresses in beams strengthened with bonded prestressed laminates. Compos B 2006; 37: 265–272.

28.

Brairi

Kerboua

Bensaid

. Improved stresses analysis of a functionally graded beam under prestressed CFRP plate. Adv Compos Lett 2018; 27(1): 10–22.

Improved numerical analysis of structures reinforced by composite FRP: Hygrothermal and prestressing loads with taper effect

Abstract

Keywords

Introduction

Basic assumptions

Analytical equations

Shear stress formulation

Normal stress formulations

Geometric and material properties

Comparison with existing analytical solutions and experimental results

Finite element analysis (application of the taper and hygrothermal model)

Different applied load

Parametric study

Effect of hygrothermal load

Linear taper effect

Effect of the prestressing force (P 0)

Plate elasticity modulus effect

Adhesive thickness effect

Conclusion

Footnotes

Declaration of conflicting interests

Funding

References

Effect of the prestressing force (P ₀)