Boundary element analysis of FRP-concrete delamination

Similar documents
Modeling of Interfacial Debonding Induced by IC Crack for Concrete Beam-bonded with CFRP

Finite Element Analysis of FRP Debonding Failure at the Tip of Flexural/Shear Crack in Concrete Beam

FRACTURE MECHANICS APPROACHES STRENGTHENING USING FRP MATERIALS

NUMERICAL SIMULATION OF THE FATIGUE BEHAVIOUR OF FRP STRIPS GLUED TO A BRITTLE SUBSTRATE

181/A, I Parma, Italy DRAFT

Heterogeneous structures studied by interphase elasto-damaging model.

Abstract. 1 Introduction

Fluid driven cohesive crack propagation in quasi-brittle materials

Influence of bond stress-slip relationship on bond strength prediction

HIGHLY ADAPTABLE RUBBER ISOLATION SYSTEMS

THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS.

3D Finite Element analysis of stud anchors with large head and embedment depth

LOCAL BOND STRESS SLIP RELATIONS FOR FRP SHEETS-CONCRETE INTERFACES

On characterising fracture resistance in mode-i delamination

Discrete Element Modelling of a Reinforced Concrete Structure

WP6 - Thought for Eurocodes Upgrade

Numerical analysis of the mechanical response of wood glulam beams reinforced through the thickness by FRP rods.

Bond-Slip Characteristics between Cold-Formed Metal and Concrete

Mechanics PhD Preliminary Spring 2017

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element

Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas

Modeling issues of the FRP detachment phenomenon

POST-PEAK BEHAVIOR OF FRP-JACKETED REINFORCED CONCRETE COLUMNS

ULTIMATE SHEAR OF BEAMS STRENGTHENED WITH CFRP SHEETS

Modeling the bond of GFRP and concrete based on a damage evolution approach

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS

Fracture Mechanics of Non-Shear Reinforced R/C Beams

Department of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil

Basic Energy Principles in Stiffness Analysis

Numerical simulation of delamination onset and growth in laminated composites

CFRP. FRP FRP. Abaqus

Homework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2

SIMULATION OF BOND FAILURE IN RC BEAMS STRENGTHENED WITH FRP SHEETS SUBJECTED TO DYNAMIC/IMPACT LOADINGS

An orthotropic damage model for crash simulation of composites

The Effects of Transverse Shear on the Delamination of Edge-Notch Flexure and 3-Point Bend Geometries

Optimal Slope of Dramix Type Fibers in Reinforced Concrete

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

Non-uniqueness of FRP bond stress-slip relationships in the presence of steel. Mehdi Taher Khorramabadi and Chris J. Burgoyne

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering

Influence of bond-slip on the behaviour of reinforced concrete beam to column joints

Piezoelectric Bimorph Response with Imperfect Bonding Conditions

MODELING OF THE WEDGE SPLITTING TEST USING AN EXTENDED CRACKED HINGE MODEL

Aim of the study Experimental determination of mechanical parameters Local buckling (wrinkling) Failure maps Optimization of sandwich panels

MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN LAYERED STRUCTURES USING GEOMETRICALLY NONLINEAR BEAM FINITE ELEMENTS

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS

KINK BAND FORMATION OF FIBER REINFORCED POLYMER (FRP)

Chapter Two: Mechanical Properties of materials

Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach

Hypersingular formulation for boundary stress evaluation revisited. Part 1: Smooth boundaries.

Mechanical Properties of Materials

International Journal of Pure and Applied Mathematics Volume 58 No ,

Modelling of non linear soil-structure interface behaviour using M. Boulon & P. G arnica

ASSESSMENT OF NONLINEAR BOND LAWS FOR NEAR- SURFACE-MOUNTED SYSTEMS IN CONCRETE ELEMENTS

Powerful Modelling Techniques in Abaqus to Simulate

Size effect in the strength of concrete structures

An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK

INTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011

Calibration and Experimental Validation of LS-DYNA Composite Material Models by Multi Objective Optimization Techniques

ISSN: ISO 9001:2008 Certified International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 2, Issue 4, July 2013

MECHANICS OF MATERIALS. Prepared by Engr. John Paul Timola

Multi-scale representation of plastic deformation in fiber-reinforced materials: application to reinforced concrete

Structural behaviour of traditional mortise-and-tenon timber joints

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS

L e c t u r e. D r. S a s s a n M o h a s s e b

Chapter 3 Variational Formulation & the Galerkin Method

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS

QUANTITATIVE ANALYSIS OF ACOUSTIC EMISSION WAVEFORMS AND. M.Ohtsu Department of Civil and Environmental Eng. Kumamoto University Kumamoto 860, Japan

Finite Element Method in Geotechnical Engineering

A Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK

Cracked concrete structures under cyclic load

Sensitivity and Reliability Analysis of Nonlinear Frame Structures

NUMERICAL INVESTIGATION OF DELAMINATION IN L-SHAPED CROSS-PLY COMPOSITE BRACKET

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

PERIDYNAMICS WITH ADAPTIVE GRID REFINEMENT

Adhesive Joints Theory (and use of innovative joints) ERIK SERRANO STRUCTURAL MECHANICS, LUND UNIVERSITY

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

University of Sheffield The development of finite elements for 3D structural analysis in fire

Dynamic analysis of a reinforced concrete shear wall with strain rate effect. Synopsis. Introduction

Module 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression

3D dynamic crack propagation analysis with PDS-FEM

ICM11. Simulation of debonding in Al/epoxy T-peel joints using a potential-based cohesive zone model

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Isotropic Elastic Models: Invariant vs Principal Formulations

FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS

arxiv:cond-mat/ v1 25 Feb 1994

MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS

Stress-strain response and fracture behaviour of plain weave ceramic matrix composites under uni-axial tension, compression or shear

Computational Analysis for Composites

Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework

STRAIN LOCALIZATION AS BIFURCATION ELASTO-PLASTIC SOFTENING MATERIALS

Experimental and numerical study on GFRP-glass adhesively bonded joints

Analytical formulation of Modified Upper Bound theorem

MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK

Test Results and Nonlinear Analysis of RC T-beams Strengthened by Bonded Steel Plates

Interface properties between a steel pre-stressing strand and an epoxy matrix.

Prediction of Delamination Growth Behavior in a Carbon Fiber Composite Laminate Subjected to Constant Amplitude Compression-Compression Fatigue Loads

5 ADVANCED FRACTURE MODELS

An alternative multi-region BEM technique for layered soil problems

Fracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design

Continuation methods for non-linear analysis

Transcription:

Boundary element analysis of FRP-concrete delamination F. Freddi 1, A. Salvadori 2 &M.Savoia 3 1 Department of Civil Engineering., University of Parma, Italy 2 Department of Civil Engineering., University of Brescia, Italy 3 DISTART-Structural Engineering,University of Bologna, Italy Abstract The FRP-concrete delamination problem is studied making use of BEM. Cohesive crack model has been adopted for the interface, whereas linear elasticity is used for the two materials outside the process zone. Numerical analyses are carried out by means of symmetric Galerkin boundary element method for cohesive interfaces, adopting the arc-length technique to follow the equilibrium path beyond its critical point. Two different test setups have been numerically simulated and results are compared with experimental tests. For bond lengths longer than minimum anchorage length, tests may exhibit a snap-back branch after the attainment of the maximum force. Keywords: SG-BEM, cohesive interface, FRP-concrete delamination, arc-length. 1 Introduction The effectiveness of strengthening concrete structures with composite materials can be compromised if debonding of Fiber-Reinforced Polymer (FRP) plate occurs. It is then very important to accurately predict the debonding failure load as well as the post-failure behavior. Recently, boundary element method (BEM) has received great attention, especially in fracture mechanics. This method is particularly attractive in the present case, since all non linearities are localized on a boundary, the interface, whereas the two materials, FRP plates and concrete, can be considered as elastic domains. The role of the interface is essential to the stress transfer between FRP and concrete. A comprehensive way of modelling FRP-concrete debonding phenomena is not yet available. At present, experimental data can be used to define mode II 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

336 Boundary Elements XXVI shear stress - tangential slip law [1], [2]. However these laws must be modified to be used for 2D or 3D problems, since compliance of concrete cover contributing to interface slip is implicitly included in these laws. The incremental symmetric boundary integral formulation for cohesive interface problems is briefly described in section 3 (see [3] for details). Riks arc-length technique with local control is used for the solution of non linear problems [4]. Several numerical simulations have been performed, concerning two different setups of delamination tests of FRP-plates bonded on concrete substrate. Due to the presence of a softening branch in shear stress-slip law, the behaviour of a specimen undergoing delamination may be strongly non linear, causing very brittle failure mechanisms, and possibility of snap-back behaviour after the attainment of maximum load. Results are found in good agreement with experimental results reported in the literature. Moreover, crack nucleation and growth can be numerically followed, until complete delamination. 2 Interface law In the so-called cohesive crack model, the fracture process zone is modelled as a fictitious crack, representing a region where stresses can be transmitted, whose values depend on crack amplitude. In this study, a mode II cohesive crack law is adopted for the FRP-concrete interface. This law must correctly reproduce elastic stiffness for low slip levels, maximum shear stress and softening branch for high slips. An important starting point is the correct definition of kinematical variables involved, that is interface slips and dual static variables (interface stresses). Hence, different models, where slips represent different kinematical variables, must adopt different interface laws. For instance, for 1D bond - slip models based on the assumption of plane strain profile for both concrete and plate (see for instance [1]), the non linear interface law must involve compliances of both adhesive layer and external concrete cover. A mode II power fractional interface law has been proposed in [2] by post-processing experimental data of delamination tests, giving the local bond stress τ as a function of slip s (see Fig. 1): τ τ = s s n (n 1) + (s/ s) n (1) where ( τ, s) denote peak shear stress and corresponding slip, and n is a free parameter mainly governing the softening branch. According to (1), linear compliances of both adhesive and concrete cover for low slip values as well as failure and subsequent post-failure softening behaviors are modelled at the interface level. A non linear constitutive interface law for FRP-concrete suitable for multidimensional domain decomposition method, such as FEM of BEM, is not yet available in literature. In the present study, the interface law has been defined starting from that proposed in [2]. That law has been modified for 2D analysis by subtracting the elastic contribution of concrete cover compliance, since this contribution is 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

Boundary Elements XXVI 337 8 [ N/mm 2 6 4 1D Model 2D Model Bilinear 2.1.2.3.4.5 s (mm) Model 1D 2D τ (N/mm 2 ) 6.8753 6.8753 s (mm).48.42 n 2.52 2.444 Gf (N mm) 1.575 1.575 K (N/mm 3 ) 237.47 276.35 Figure 1: Mode II non-linear FRP - concrete interface laws: the proposed (2D) interface law, a 1D bond - slip law (from [2]) and a bilinear law. explicitly taken into account in the BEM. Hence, the initial stiffness of 2D interface law, obtained from (1) as: τ (s =)=K2D = τ s n s n 1 (2) is given by: 1/K 2D =1/K 1D 1/K c,wherek 1D is the initial stiffness of 1D bond - slip law, and K c = G c /h c is the stiffness of external concrete cover, with G c being shear modulus and h c =26mm the equivalent thickness of concrete contributing to interface compliance [1]. Moreover, the value of fracture energy of 1D law has been preserved, since that value was calibrated from experimental tests in order to give the value of maximum transmissible force by an anchorage of infinite length [5]. Fracture energy of interface law (1) can be written as: G f = τ sπ [1/(n 1)] (1 2 n) csc (2π/n) (3) Hence, for prescribed values of K 2D and G f, the parameters s and n of 2D interface law can be easily obtained from eqns (2), (3). The interface law for 2D analyses is reported in Fig. 1. In the same figure, a bilinear law is also depicted, obtained by prescribing the same peak values ( τ, s) and fracture energy G f of 2D power fractional law. It is worth noting that, at present, no experimental information are available to calibrate a non linear law for normal stress-crack opening displacement. Therefore, in the present analysis, a linearly elastic law is adopted, σ = k p w, with compliance due to adhesive deformation, i.e., k p = E g /h g, E g and h g being the elastic modulus and the thickness of the adhesive layer respectively. Denoting with p = def {σ, τ} the cohesive tractions vector and with w = def {w, s} the relative opening displacements vector, the interface constitutive equation can 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

338 Boundary Elements XXVI be written in an incremental form making use of the tangent stiffness matrix of cohesive law, denoted with D T : ṗ = D T (w) ẇ. 3 Boundary integral formulation Single-zone boundary integral equations - Consider a homogeneous solid with domain Ω R 2 and boundary Γ=Γ u Γ p. Assuming small strains and displacements, consider the following quasi-static external actions: tractions p(x) on Γ p, displacements ū(x) on Γ u ; the constitutive law is assumed to be isotropic and linear elastic. Boundary integral equations (BIEs) for the linear elastic problem can be derived from Somigliana s and hypersingular identities by performing the boundary limit Ω x x o Γ. Somigliana s identity, which is the boundary integral representation of displacements in the interior of the domain x Ω, is based on Green s functions which represent components u i of the displacement vector u in a point x due to: i) a unit force concentrated in space (point y) and acting on the unbounded elastic space Ω in direction j (such functions are gathered in matrix G uu ); ii) a unit relative displacement concentrated in space (at a point y), crossing a surface with normal l(y) and acting on the unbounded elastic space Ω (in direction j) (gathered in matrix G up )). The hypersingular identity, which is the boundary integral representation of tractions on a surface of normal n(x) in the interior of the domain, involves Green s functions (collected in matrices G pu and G pp ) which describe components (p i ) of the traction vector p on a surface of normal n(x) due to: i) a unit force concentrated in space (point y) and acting on the unbounded elastic space Ω in direction j; ii) a unit relative displacement concentrated in space (at a point y), crossing a surface with normal l(y) and acting on the unbounded elastic space Ω (in direction j). The two integral equations are usually referred to as displacement and traction equations and read as follows [6] for x Γ, havingsetr = x y: C(x)u(x)+ G up (r; l(y))u(y)dy = G uu (r)p(y)dy (4) Γ D(x)p(x)+ = G pp (r; n(x); l(y))u(y)dy = G pu (r; n(x))p(y)dy (5) Γ Γ As seen from equations (4)-(5), singularities of Green s functions are triggered off in the limit process, thus generating free terms C(x) and D(x) and giving to the integrals the distributional connotation of finite part of Hadamard and of Cauchy s principal value. Incremental multi-zone formulation - Consider N domains Ω n, n =1, 2,..., N connected to each other by M interfaces. Quasi-static external tractions t n are imposed on the Neumann boundary Γ n p of each domain Ωn, whereas displacements ū n are imposed on the Dirichlet boundary Γ n u of each domain Ω n. The incremental BIE formulation of the problem has been proposed in [3] and for paucity of space will be here merely briefly summarized starting from its compact form: Γ M u =ṗ (6) 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

Boundary Elements XXVI 339 Equation (6) represents a set of BIEs obtained from the incremental form of equations (4), (5) on the Dirichlet and Neumann boundaries of each domain, as well as suitable linear combinations, by means of the transposed tangent stiffness matrix D T T, of equations (4), (5) written on the two faces of each interface (see [3] for further details). It has been proved [3] that if D T = D T T then the integral operator M is symmetric with respect to a suitable bilinear form A, thatis:a(m u, v) = A( u, M v). It is then possible to give a variational formulation to problem (6), whose solution is a critical point of the functional Ψ[ u] = 1 2A( u, M u) A( u, ṗ) which does not possess a sign definiteness property unless the interface law is stable in a second order sense (D T negative definite). Symmetric Galerkin approximation - The Galerkin approximation scheme [7] ensures stability and convergence of the numerical solution in suitable functional spaces. Let h>be a parameter and let U h denote a family of finite dimensional subspaces of U such that u U,inf uh U h u u h as h. The symmetric Galerkin approximation of (6) consists in finding the critical point u h U h of the functional Ψ[ u h ]= 1 2 A( u h, M u h ) A( u h, ṗ) (7) From the algebraic point of view, let ψ u (y) be a matrix of shape functions and u h (y) =ψ u (y) u be discrete approximations of u(y). From the stationarity condition of (7), the linear system M u = ṗ (8) can be obtained, where M = A(ψ u, Mψ u ) and ṗ = A(ψ u, ṗ). Integration in time - Problem (8) is a set of first order differential equations. Supplied with initial conditions, a unique solution u exists for any given term ṗ provided that discrete operator M(u) is invertible. Such a condition, in the analysis of non-linear responses of structures, is not a priori fulfilled if the interface law is not stable in a second order sense (i.e. D T is not negative definite). Solution u of (8) in the range of D T negative definite may be approximated [4] by standard integration schemes: explicit Runge-Kutta methods of order two and four have been adopted in the present paper. Arc-length method - Consider the non-linear incremental problem (6) where the structural response u is obtained for p varying from a null initial value and evolving quasi statically in time by means of a load factor λ(t): ṗ = λ p, supposing that λ at least initially increases. The problem under consideration presents a limit points and conventional incremental algorithms, based on a fixed value of load, may fail to overcome such points. Hence, path-following techniques must be used, whereby the load factor λ is added to the set of unknowns. Equation (8) becomes therefore: M u = λ p (9) The number of unknowns increases by one and the set of governing equations is given by (9) plus an additional constraint equation which accounts for the increment of the load factor. The arc-length type method is based on the introduction 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

34 Boundary Elements XXVI of a control function giving a measure of the evolution of the loading process. The pioneering form of the constraint equation, proposed originally in [8], reads as follows: u 2 + λ 2 = c(λ, u) 2 (1) Huge efforts have been made in the last two decades - see [9] for an exhaustive review - in order to improve computational efficiency of constraint equation (1). In the present paper, a local control function is used, analogous that proposed in [1], with arc-length suitably adapted by the curvature along the equilibrium path. Differently from all previously cited papers, the present formulation is differential in time whereby the solution is reached by (explicit) time integration strategies. Details of the adopted approach, coupling boundary integral equations and arclength technique, are reported in [4]. 4 Simulation of delamination problems via BEM The proposed BE model for cohesive interfaces has been used to simulate some FRP - concrete delamination tests. Two different test setups are considered, depicted in Figs. 2a and 7a, together with mechanical and geometrical properties of specimens. For FRP-concrete interface, the non linear interface law reported in Fig. 1 has been adopted. 4.1 Simulation of Setup N. 1 delamination test Fig. 2a shows a typical configuration for pull-pull delamination tests. Left and bottom sides have been constrained in order to have no displacements in the direction normal to the surface and free displacement tangent to it. Two different numerical tests have been performed. 4.1.1 Comparison with experimental tests Experimental test performed by Chajes et al. [11] have been considered. Experimental results were post-processed in Savoia et al. [1]: FRP - concrete shear stresses and slips along the bond length have been derived from experimental strains in FRP plate. In these tests, FRP plate width (25.4 mm) was much smaller than concrete width. Nevertheless, in the present simulation a plane stress problem has been considered. Comparison between experimental data and numerical results obtained through the present BE model are reported in Figs. 3a-b, for different bonding lengths. In each figure, strains in FRP plate along the bonding length are reported, for different values of applied force. The highest load corresponds to experimental failure load. All figures show good agreement between numerical and experimental results both at low and high loading levels. Agreement at low loadings (the interface behavior being almost linear) means that the initial interface stiffness is correctly defined. The agreement is very good also for high loadings, where plate delamination is occurring. 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

Boundary Elements XXVI 341 t (mm) 1.16 h (mm) 152.4 E p (MPa) 23 E c (MPa) 3364 f c (MPa) 36.4 b p (mm) 25.4 b c (mm) 152.4 d (mm) 5 Figure 2: Geometrical and mechanical properties of specimens tested in [11]. Strain ( ) 5 4 3 2 L 5.8 mm 1934 N 459 N 69 N 895 N Strain ( ) 5 4 3 2 L 23.2mm 219 N 4624 N 71 N 9198 N 1157 N 1 1 2 4 6 x (mm) 5 1 15 2 25 x(mm) Figure 3: Strains along FRP plates: experimental results [11] and numerical simulations. Moreover, delamination failure load as a function of bonding length is reported in Fig. 6a. Experimental results are compared with numerical predictions obtained from the proposed BE solution and those given by a recently proposed 1D bondslip model [2]. The two different models give practically the same results, for both small and long bonding lengths. This circumstance confirms that the corresponding interface laws are actually equivalent. In the same figure, dashed line indicates failure loads obtained by adopting the bilinear interface law (also reported in Fig. 1) together with 2D BE solution. It is clearly shown that failure loads obtained using the bilinear law are significantly overestimated. The asymptotic value (for bonding length approaching infinity) is practically the same, since the same value of fracture energy of interface law has been adopted. 4.1.2 Other numerical tests With reference to the same test setup and mechanical properties of specimens depicted in Fig. 2a, different bond lengths have been considered (L = 5 and 2 mm), whereas the plate width is 5 mm. In this case load-displacement curves obtained from numerical simulations are reported in Figs. 4a-b, compared with results given by 1D model described in [1]. Abscissa refers to relative axial displacement u between the initial and the final section of FRP plate (axial elongation of the plate). It is worth noting that, for long bonding lengths (see Fig. 4b), snap-back branch occurs after the attainment of maximum load, due to elastic shortening of FRP 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

342 Boundary Elements XXVI P(kN) 2 15 1 5 SG-Bem 1D Model.2.4.6.8 u(mm) P(kN) 4 35 3 25 2 15 1 5 SG-Bem 1D Model.2.4.6.8 u (mm) Figure 4: Pull-pull delamination test for L =5mm and L = 2 mm..4 8.3 6 s (mm).2.1 A B C s=s (N/mm 2 ) 4 2 A B C 4 8 12 16 2 x (mm) 4 8 12 16 2 x (mm) Figure 5: FRP-concrete slip s and shear stress τ distributions along the interface. plate when axial loading is decreasing during delamination. Of course, this simulation cannot be performed with a standard force or displacement control procedure, and more refined techniques are required, as the arc-length technique described in Section 3. Snap-back branch is not present in the case of short bonding lengths (see Fig. 4a). For the 2 mm bonding length case, slip and shear stress profiles along the anchorage are reported in Figs. 5a-b. Curves refer to three different equilibrium points (A, B, C) underlined in Fig. 4b. Straight line in Fig. 5a shows the slip value corresponding to peak shear stress of interface law (s =.421 mm). Hence, when slip is higher, softening branch in shear stress profile occurs (see Fig. 5b). Finally, for the 2 mm bonding length, load-plate elongation curve is reported in Fig. 6a, adopting the proposed power fractional law and the bilinear law for the interface. The figure confirms that, adopting the bilinear law, higher values of delamination load are predicted; moreover, the snap-back branch is given by a straight line until complete delamination. 4.2 Simulation of Setup N. 2 test In Setup N. 2, depicted in Fig. 7a, the specimen is symmetric and the load is applied at two external FRP plates. Moreover, both concrete and FRP plates are restrained at the left end section. Bond length is L = 6 mm. 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

Boundary Elements XXVI 343 P max (kn) 2 15 1 5 Chajes et al. 1D Model SG-Bem Bilinear 1 2 3 4 Anchorage Length (mm) P(kN) 5 45 4 35 3 25 2 15 1 5 Bilinear Proposed law.2.4.6.8 u (mm) Figure 6: a) Maximum transmissible force vs. bond length: experimental and numerical results obtained with different interface laws (b p =25.4 mm). b) Pull-pull delamination test for L = 2 mm: solution obtained with different interface laws (b p =5mm). The corresponding load-displacement curve is reported in Fig. 7b. The figure clearly shows the transition between State 1 condition for low level loads (both stiffness contributions of concrete and FRP plates are present) and State 2 condition after complete FRP delamination (FRP plates only contribute to specimen stiffness). Hence, in this case the post-delamination branch is stable and mono- 6 5 4 State 1 B C P(kN) 3 2 A State 2 1 1 2 3 L (mm) Figure 7: Geometrical properties of the specimen and structural response. 2.5 8 s(mm) 2 1.5 1 C B A (N mm 2 ) 6 4 C B A.5 2 2 4 6 x (mm) 2 4 6 x (mm) Figure 8: Setup N.2: slip s and shear stress τ distributions along the interface. 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6

344 Boundary Elements XXVI tonically increasing. Figs. 8a-b show the profiles along the bonded length of FRPconcrete slip and shear stresses, with reference to three different equilibrium points A, B, C underlined in Fig. 7b. Point A indicates the beginning of delamination process, point C a condition where almost the whole plate is delaminated, and point B an intermediate condition. Fig. 8b clearly shows that, during delamination, the portion of plate contributing to the anchorage translates along the specimen from the (right) loaded end to the (left) restrained section. Acknowledgements The financial supports of (Italian) MIUR (PRIN 23 Grant, FIRB Grant) and C.N.R., PAAS Grant 21 are gratefully acknowledged. References [1] M. Savoia, B. Ferracuti, C. Mazzotti, Delamination of frp plate/ sheets used for strengthening of r/c elements, insecond International Structural Engineering and Construction (ISEC-2) Conference Proceeding, Rome (23). [2] M. Savoia, B. Ferracuti, C. Mazzotti, Non linear bond-slip law for frpconcrete delamination, ink.tan,ed.,frprcs-6 Conference Proceedings, Singapore (23), pp. 1 1. [3] A. Salvadori, A symmetric boundary integral formulation for cohesive interface problems, Computational Mechanics, 32(4-6), (23), 381 391. [4] F. Freddi, Cohesive interface analysis via boundary integral equations,ph.d. thesis, University of Bologna, Italy (24). [5] C. Mazzotti, B. Ferracuti, M. Savoia, An experimental study on frp-concrete delamination, inframcos-5, Fracture Mechanics of Concrete and Concrete Structures, Colorado, U.S.A (24). [6] K. Hong, J. Chen, Derivations of integral equations of elasticity, ASCEJ. Engrg. Mechanics, 114, (1988), 128 144. [7] W. McLean, Strongly elliptic systems and boundary integral equations,cambridge University Press, New York (2). [8] E. Riks, An incremental approach to the solution of snapping and buckling problems, Int. J. Solid Structures, 15, (1979), 529 551. [9] M. Geers, Enhanced solution control for physically and geometrically nonlinear problems. part 1. the subplane control approach, Int. J. Numer. Methods Eng., 46, (1999), 177 24. [1] I. May, Y. Duan, A local arc-length procedure for strain softening., Computers and Structures, 64, (1997), 297 33. [11] M.Chajes,W.J.Finch,T.Januska,T.T.jr.,Bond and force transfer of composite material plates bonded to concrete, ACI Struct. J., 93, (1996), 28 217. 24 WIT Press, www.witpress.com, ISBN 1-85312-78-6