Experimental and Numerical Study of Behaviour of Reinforced Masonry Walls with NSM CFRP Strips Subjected to Combined Loads

: Near surface mounted (NSM) carbon fibers reinforced polymer (CFRP) reinforcement is one of the techniques for reinforcing masonry structures and is considered to provide significant advantages. This paper is composed of two parts. The first part presents the experimental study of brick masonry walls reinforced with NSM CFRP strips under combined shear-compression loads. Masonry walls have been tested under vertical compression, with different bed joint orientations 90° and 45° relative to the loading direction. Different reinforcement orientations were used including vertical, horizontal, and a combination of both sides of the wall. The second part of this paper comprises a numerical analysis of unreinforced brick masonry (URM) walls using the detailed micro-modelling approach (DMM) by means of ABAQUS software. In this analysis, the non-linearity behavior of brick and mortar was simulated using the concrete damaged plasticity (CDP) constitutive laws. The results proved that the application of the NSM-CFRP strips on the masonry wall influences significantly strength, ductility, and post-peak behavior, as well as changing the failure modes. The adopted DMM model provides a good interface to predict the post peak behavior and failure mode of unreinforced brick masonry walls.


Introduction
Masonry is an anisotropic material, considering only plane homogenous stress states, it is characterized by many different failure modes and performances [1,2]. This property strongly affects the response of masonry wall subjected to in-plane seismic forces. From the uniaxial compressive load tests on masonry panels in which bed joints are oriented at different angles to the direction of the applied load, the combined shear-compression behavior of masonry have also been reported by some researchers [3,4]. The failure of masonry under uniaxial compression combined shear and compression has been extensively studied in the past by many researches [5]. These failures represent particular points on the general failure surface. The influence of the orientation of the applied stresses to the joints has also been discussed in the literature concerning shear wall behavior [6].
Around the world, there are many old unreinforced masonry (URM) structures such as historical cultural monuments and bridges. Many of these structures are exposed either to excessive loads or to unpredicted loading scenarios such as earthquakes [7,8]. Different techniques have been developed to prolong the life of these structures. Several of these reinforcement systems include the application of fiber-reinforced polymer (FRP) composites to strengthen masonry structures. Using FRPs offers several advantages such as increasing strength and ductility of walls subjected to either out-of-or inplane lateral loading. El-Gawady et al. [9] and Korany et al. [10] have reviewed the techniques of retrofit/repair masonry walls with FRP composites. Such composites offer the possibility of either adhering them on the outside surface (externally bonded-EB) or inserting them inside a groove of a masonry element by the NSM technique. Seracino and Wylie [11] recommended that the placement of NSM FRP in the vertical and/or horizontal direction depends on the loading and boundary conditions of the wall. Petersen et al. [12] carried out tests to study the bond behavior of clay brick masonry prisms reinforced with vertical strips and horizontal strips (CFRP perpendicular and parallel to bed joints). They applied CFRP strips with a rectangular form to maximize the confinement from the surrounding concrete. Tests noted that the principal failure mode was the debonding of CFRP from the masonry wall for both orientations of CFRP strips. For solid bricks with vertical NSM-CFRP strips inserted into brick only, they found that the bond strength decreased by a rate of 8%. If the vertical NSM-CFRP strips passed through mortar head joint, a reduction in the bond strength of 11% was observed. Further decrease (31%) in bond strength was recorded when the FRP was aligned horizontally. A number of researchers investigated the effect of single-sided (or non-symmetric) strengthening. Valluzzi et al. [13] examined a series of unreinforced brick masonry walls strengthened by different FRP composites such as CFRP, GFRP, and AFRP fibers of different forms and subjected to a diagonal compression test. The results showed that the application of FRP at one side of the wall produces a significant failure mode of out-of-plane deformation. The magnitude of this failure intensified by the limited restraint at the bottom and top corners of the wall. Petersen et al. [14] investigated the in-plane shear behavior of masonry walls reinforced with NSM CFRP using different reinforcement orientations such as vertical strips, horizontal strips, and a mixture of both. It was observed that the application of vertical NSM-FRP strips applied to both front and backsides of the masonry walls led to 46% improvement of the load capacity. In addition, when both horizontal and vertical NSM strips were used, the horizontal strips prevented the opening of diagonal cracks while the vertical strips prevented the sliding failure. Dizhur et al. [15] performed several tests with various NSM-CRFP retrofitting and repair schemes on URM walls loaded in diagonal compression. They reported that retrofitted walls showed an increase in the maximum shear strength ranging from 1.3 to 2.6 times. For repaired walls, the increase ranging from 1.3 to 3.7 times compared to the masonry wall. Furthermore, they observed a considerable increase in ductility of 2.6 times for walls retrofitted on one side and 25.5 times for the walls retrofitted on both sides. The effect of different parameters including the groove dimensions, the dimensions of CFRP, shape of CFRP, and the adhesive type in NSM-strengthened brick masonry has been investigated by Maljaee et al. [16]. Existing experimental and analytical research has led to two major approaches for masonry modeling, namely heterogeneous and homogeneous modeling. In the former, the brick units and joint mortar are considered separately. While in the latter, the brick units, joint mortar and interfaces between them are assumed by an isotropic or anisotropic composite material. According to Lourenco [17], there are two main modelling approaches macro-modelling and micro-modelling. The micro-modeling can be divided into two techniques: detailed micro-modeling and simplified micro modeling [18,19]. In the detailed micro-modeling, each component (units, mortar) is modeled separately with unit-mortar interface. The masonry units and mortar are modeled with true thickness, whereas zero thickness is attributed to unit-mortar interfaces. The units and mortar are modelled as continuous elements and unit-mortar interfaces are modelled as discontinuous elements.
To understand the behavior of masonry, the literature has focused on the importance of including all rupture mechanisms of the masonry in the modeling in terms of ultimate load and ductility on all the damage, which is normally concentrated at the mortar interface [20,21]. Modeling the complex behavior of masonry is a difficult task in the finite element analysis of civil engineering structures. Only tensile brick rupture is often considered in microscopic models, relying mainly on linear elastic behavior. The authors took that model to introduce the mortar's post-peak softening behavior.
There is no published literature on reinforced masonry walls that discusses the shear test when the loading is applied by means of compression force only. In such a scenario, the bed joints are at an angle with the loading direction. However, the effect of combined stresses to the strength and deformation of this reinforced masonry system is not widely available in the literature. With a view to bridging this gap, experimental studies on the behavior including ductility and failure mode of NSM CFRP reinforced brick masonry wall under combined shear-compression were carried out.
This paper presents masonry panels with two bed joint orientations that were constructed and tested under vertical compression to generate a wide range of compression to shear ratios from one to infinity at each of the bed joints of the panels. The effect of the type of joint mortar, the efficiency of NSM-FRP technique and the position of the CFRP strips for improving the shear strength and ductility of the reinforced wall are discussed in detail. An analytical approach is developed to simulate the behavior of unreinforced masonry walls under combined load, using the detailed micromodelling approach (DMM) implemented in ABAQUS. The material nonlinearities of units, mortar and unit-mortar interfaces in both tensile and compression regimes have been assigned. In this analysis, the non-linearities behavior of brick and mortar was simulated using the Concrete Damaged Plasticity (CDP) constitutive laws.

Brick Units
To evaluate the compressive strength and elastic modulus of masonry units, nine specimens of perforated clay bricks with dimensions 220 mm long, 105 mm wide and 55 mm high were tested according to EN 771-1 [22] (Figure 1). The compressive strength and elastic modulus values are presented in Table 1.  Two types of mortar consisting of cement, lime and sand were prepared. The proportions of cement, lime and sand were 1:1:3 and 1:1:5 respectively. Both types of mortar had a water to cement ratio of 0.5. Flexural and compressive resistance were determined in accordance with EN1015-11 [23] and the results are presented in Table 1 The CFRP strip has an elastic modulus equal to approximately 165,000 MPa, and rupture strain equal to 1.7%. The properties of composite materials of CFRP applied in the reinforcement of masonry wall panels are presented in Table 2.

Methods
The masonry panels were categorized into two series. In the first series, the compressive behavior is defined on a panel subjected to uniaxial compression perpendicular to bed joint θ = 90° with two different types of mortar (A, B). In the second series, a combined shear-compression behavior of un-strengthened and strengthened masonry (θ = 45°) has been studied using the same types of mortar.
The experimental program consisted of a total of 20 tests, 6 un-strengthened masonry panels as control specimens (MCA, MCB) with θ = 90° and (MTA, MTB) with θ = 45°, and 12 strengthened panels. Four of the strengthened panels have only horizontal NSM CFRP strips (parallel to bed joints) with both types of mortar (MRHA, MRHB), and four having only vertical NSM CFRP strips (perpendicular to bed joints) with both types of mortar (MRVA, MRVB). The remaining four panels were reinforced with two vertical strips on one side of the panels and two horizontal strips on the other side MR2A and MR2B (Table 3). Full view of walls specimens is shown in Figure 2.  The test specimens were made according to RILEM technical recommendation for testing small walls under compression [24]. The specimen was a masonry panel made of perforated bricks with dimensions 220 × 105 × 55 mm. The specimens were tested under uniaxial loading with different orientations of the bed joints. The varied inclinations were 90° and 45°.
The specimens were aligned between the platens of the testing machine carefully to avoid eccentricity in loading. Furthermore, to avoid any contact between the steel plates and the surfaces of the masonry specimen, support plates were placed to ensure a better distribution of forces without any eccentricity. Full contact between upper and lower surfaces of specimens and surfaces of steel plates of the testing machine was assured. A computer-controlled servo hydraulic frame with a capacity of 160 kN was used to carry out the tests. The displacement measurements were carried out using comparators that have been installed on both sides of the panel. All test data were collected by a data acquisition system.
The tests were performed under displacement control in order to obtain the complete stressstrain curve of the panels. All tested wall panels were of similar dimensions in order to allow direct comparison of their failure loads. The FRP reinforcement was glued, using two-part epoxy adhesive, into rectangular grooves cut in the surface of the masonry with a circular saw. Full view of the fabrication of specimens and installation of CFRP is shown in Figure 3.

Unreinforced Masonry Panels Subjected to Uniaxial Compression θ = 90° (MCA, MCB)
Panels subjected to uniaxial compression perpendicular to bed joint (θ = 90°) collapsed due to vertical cracking at the level of face shells of the masonry wall ( Figure 4). The face-shell cracking occurred at a load near to 95% of the ultimate load. The vertical cracks on the face shell perpendicular to the joints can be due to different levels of lateral expansion of the units and mortar under compressive stresses, which causes tensile splitting of brick and perpendicular joints. The splitting cracks result in two face shells that deform individually and become more fragile. This failure pattern has been reported by other researchers [25,26]. A reduction of mortar strength with an increase in the proportion of sand in mortar is observed, whereas a slight increase was observed in the value of compressive strength for masonry panels MCB compared to masonry panels MCA (2.23 MPa and 2.16 MPa respectively). This indicates the insignificant effect of mortar strength on the compressive strength of the masonry wall (3%). The elastic modulus of the masonry wall is calculated as the second modulus of the stress-strain curve obtained during an experimental prism test within stress levels of 0.05 and 0.33 times σm.  The values of Young's modulus of masonry panels obtained for each wall MC are summarized in Table 4. The results show that the value of the modulus of linear deformation is low compared with the correlations between E and fc provided by the codes: E/fc = 1000 as in CSA 2004 [27] and E/fc = 850 in EN 1052-1 [28]. Most of the formulae that calculate the elastic modulus of the masonry give a value greater than the experimental value. Augenti et al. [29] found that the elastic modulus varies between 250-1100 times the compressive strength of masonry. They proposed an average value of Young's modulus equal to 550 times the compressive strength. The results obtained show that the measured values are close to those found by FEMA306 [30], which propose E, mac ≈ 550σc, mac, where, σc, mac is the compressive strength of masonry units.  [29,30] 1245.2 1224.3 Venkatarama Reddy and Uday Vyas [31] investigated the influence of bond strength on stressstrain characteristics of masonry employing a cement-lime mortar and soil-cement blocks. These studies show that when the masonry unit is stiffer than that of mortar (Eblock/Emortar ratio greater than one), the masonry compressive strength is not sensitive to bond strength variation and the modulus decreases with an increase in bond strength. In addition, the modulus of masonry is less than that of the block and the mortar when Eunit/Emortar ratio is less than one. However, the results of the present study indicated that the modulus of masonry is less than that of the block and mortar, although the Eunit/Emortar ratio is greater than one (Table 5). The combined shear and normal stresses play a significant role in the failure modes and deformation on the levels of the bed and head joints. For a uniaxial state of stress that is inclined relative to the x-axis at an angle θ, if angle θ varies, the normal stress decreases and shear stress τn increases on an inclined plane. The maximum shear stress of magnitude occurs on the planes oriented at 45° to the x-axis. Using equilibrium Equations (1-3) along the bed and perpendicular joints respectively, the applied vertical compressive stress can be converted to compressive and shear stresses (σn, τnp) and (σp, τnp). The three linear strains (εx, εy, ε45) can be used to determine the normal and parallel strains (εn, εp) and the shear strain (γnp) on the plane of the specimen using the strain transformation equations [32], in which θ is the angle of bed joint to the x axis ( Figure 5).
In case θ equals 45°, the shear stress (τ) equals the normal stress ( ) on the sliding surface. The shear modulus (G) is calculated using: = The shear angle is calculated following: In wall panels MTA and MTB, cracking occurred primarily in a diagonal through the mortar joints followed by a rapid reduction in load capacity. With load increasing, the wall exhibited a gradual increase in the width of a predominately diagonally oriented crack, followed by sliding along the formed cracks (Figure 6a). With further increase in load, multiple cracks were observed in the panel MTB before failure as shown in Figure 6b. The ultimate loads of MTA and MTB panels were measured as 69.2 kN and 91.0 kN, respectively ( Table 6). The response of the tested unreinforced wall panels (MTA and MTB) to shear stress-strain is summarized in Figure 6c. The shear strength of the masonry panel constructed with mortar A was lower than the shear strength of the corresponding masonry panel constructed with mortar B. It is observed that both MTA and MTB wall panels presented an approximately linear shear stress-strain relationship before cracking, followed by rapid degradation of shear strength once cracking had propagated, followed by a slight increase in shear strength and deformation before rupture (Figure 6c).  The wall panels tested with an orientation of bed joint by 45° failed at lower load compared to the wall panel tested when its bed joint makes 90° with loading axis. The same observation was found for the Young's modulus. For these orientations, the load bearing capacity was affected by the brick strength and the shear bond characteristics of the joints. Since the masonry joints act as planes of weakness in the traditional masonry, they tend to degrade the orthotropic behavior under different orientations of bed joints to the direction of loading. The ratio of the compressive strength at 45° (f45°) to the compressive strength perpendicular to the bed joint (f90°), was found at a range of 0.38 to 0.48, which is highly orthotropic. This is what led to a reduction of strength by up to 62%. Consequently, the failure load of the unreinforced wall panels was highly dependent on the bed joint orientation. As the orientation changes from 90° to 45°, the average strength value reduced from 2.22 MPa to 1.08 MPa for masonry panels MTB and from 2.16 MPa to 0.82 MPa for masonry panels MTA, which represented about 40% of the strength reduction. The shear modules of wall panel A and wall panel B were determined as 125 MPa and 137.3 MPa, respectively.

Reinforced Masonry Walls Subjected to Uniaxial Compression with θ = 45°
The reinforcement schemes used for the strengthening of panels are summarized in Figure 1b: • Panels MRVA and MRVB were reinforced with two vertical strips on one side; • Panels MRHA and MRHB were reinforced with two horizontal strips on one side; • Panels MR2A and MR2B were reinforced with two vertical strips on one side and two horizontal strips on the other side. The distance between the staggered vertical reinforcement was 135 mm. The criteria established by Park [32] were used to calculate the displacement ductility factor of all tested panels. The displacement ductility factor ( ) is defined as the ratio between the ultimate displacement to the yield displacement: δ = displacement at ultimate load δ = displacement at the load causing yield condition In unreinforced walls, the tensile stress causes an appearance of cracks leading to a complete destruction. However, in the case of walls reinforced with CFRP strips, the tensile stresses are transferred to these strips resulting in the masonry wall being greatly reduced in stress. As the load increases, the maximum tensile stress occurs in the corner of the wall. Figure 5 illustrates the observed failure modes for all strengthened masonry panels. Wall panels reinforced by NSM CFRP strips presented a ductile behavior compared to the control wall panels. The results in Table 4 showed that, when the ratio Eunit/Emortar was doubled, the compressive strength of masonry panels increased by 102% and the shear strength increased by 132%. Furthermore, the Eunit/Emortar ratio was found to have a significant impact on the in-plane shear capacity of the shear walls. Likewise, it was found that the improvement in shear strength for strengthened wall panels type A ranged from 123% to 185%, whereas for wall panels type B it ranged between 142% to 196% (Table 6). Concerning the behavior of the reinforced panels, the ultimate load in the (MRVA, MRVB) panels reinforced with vertical CFRP strips on one side only, had an increase of up to 149%. Furthermore, the displacement increased to 200%, which led to an increase of ductility by 384% (Table 7). These results show that all the reinforced walls exhibited a significant improvement of ductility when compared with the corresponding control wall panel (μu/μo). They showed a substantial increase in deformation capacity, which remained between 2.26 and 4.27 times (between 222% and 392% as a percentage). Similar findings were also reported in literature by Dizhur et al. [15]. The vertical reinforcement contributed more to strength enhancement as compared to horizontal reinforcement. Wall panel reinforced using horizontal CFRP strips as shown in Figure 7f,g failed by sliding along the bed joint (sliding shear failure mode), which resulted in a substantially lower increase in shear strength when compared with wall panels having vertical-oriented reinforcement. The reinforced walls with two horizontal CFRP strips showed a brittle failure with a sudden loss of strength. On the other hand, the reinforced walls using vertical CFRP strips revealed less brittle failure mode with larger deformation (Figure 7h,i). Similarly, the studies performed by Parvin and Syed Shah [33] and Seracino and Wylie [34] investigated the efficiency of using vertical-oriented discrete FRP strips. These studies have shown that the panels reinforced with vertical CFRP strips showed an increase in vertical moment capacity and deformation capacity of walls under out-of-plane loading or in-plane loading. The reinforced walls (MR2A and MR2B) exhibited a vertical splitting of the interior webs followed by a gradual increase in the load up to the peak load. After reaching the peak load, the webs of the blocks completely fractured and delaminated from the face shells as shown in (Figure 7e). Petersen et al. [14] suggested that the vertical reinforcement was very effective in restraining sliding and diagonal cracking and hence preventing the URM failure mode. The non-symmetrical reinforcement regimes resulted in out-of-plane deformations and, therefore, it could not be avoided by reinforcing the panels at both sides (Figure 7a-d). Higher strains in vertical and horizontal directions were recorded for the MR2B panel while lower values were obtained for MRHA panel as shown in Figure 8. The strain values recorded in type B masonry panels were much higher than that in type A panels. The higher shear stress-strain values for masonry panels MR2 were primarily due to confinement of masonry on both sides of the wall. No rupture of the CFRP strips was observed during testing or when the CFRP strips were exposed during demolition of the wall panels. With increasing wall panel deformation, the debonding failure mode and pull out of the middle CFRP strips was not observed.

Material Parameters
In this study, the FE software ABAQUS was used to evaluate the validity of Detailed Micro-Modeling (DMM) as a numerical model for brick masonry wall analysis. The units and mortar joints are modelled using eight-nodded 3D continuum elements with hourglass control and reduced integration (C3D8R), and the unit-mortar interface with zero thickness. The penalty friction formulation was introduced to model tangential behavior with a friction coefficient value of 0.6. The normal behavior of the interface was modeled using "hard contact". Furthermore, the nonlinear behavior of brick and mortar was simulated by using the CDP model while the interface between unit and mortar was modeled through surface-to-surface contact. The adopted modeling strategy used in this study is illustrated in Figure 9.

Presentation of the Numerical Model
The numerical study will be limited to analysing the same unreinforced brick masonry wall (MCB and MTB) that was studied in the experimental part. Figure 10 shows the numerical simulation, the geometry and loading condition for the FE model that has been implemented using ABAQUS.

Constitutive Behavior of Units and Mortar a. Concrete Damage Plasticity (CDP) Model
The CDP model provides a general capability for modeling concrete and other quasi-brittle materials, which was developed by Lubliner et al. [35], and is available in the ABAQUS [36] software library. The CDP model uses the concepts of isotropic damaged elasticity in combination with isotropic tensile and compressive plasticity [36]. This approach has been developed to predict failure modes as well as tensile fracture and the compressive crushing; it is selected in this study to simulate the nonlinear behavior of the masonry units. The level of damage is evaluated by damage parameters dc and dt. The CDP model was based on the uniaxial stress-strain relationship: The damaged parameter ( ) is calculated by equation (6), this parameter ranging from 0 to 1: σ′ The compressive strength of masonry Damaged parameter ( ) can be calculated by equation (7), this parameter ranging from 0 to 1.
′ Masonry tensile strength The uniaxial compressive and tensile responses of concrete with respect to the concrete damage plasticity model subjected to compression and tension load were given by: The effective uniaxial stress and were derived as follows: b. Constitutive Behavior of Units-Mortar Interface In this model, brick-mortar interface is modeled by means of Mohr-Coulomb frictional models in order to adequately reproduce the shear response of joint. Cohesion and friction angle values for Coulomb type friction model are given by the shear strength envelope for different normal stress values, which was obtained by the experimental triplet test ( Figure 11). The interaction module of Abaqus/Explicit analysis was used to make the contact between units and mortar through the option surface-to surface contact. In this step, it is necessary to define two contact properties: Normal contact and tangential behavior. Under normal behavior: The hard contact behavior normal to the surfaces is chosen. The goal is to prevent interpenetration of surfaces and to permit a separation between them once a contact has been established.
Tangential behavior: The analysis needs to take frictional forces because when the surfaces are in contact, they habitually transmit shear and normal forces athwart their interface, which resist the relative sliding of the surfaces.

Model Input Parameter
The input parameters used to reproduce the mechanical behavior of masonry wall using the detailed micro-modeling approach are presented. The mechanical parameters for brick and mortar employed in this analysis were obtained from tests carried out in this research (see Tables 8-10).  Figure 12 shows the numerical and experimental curves of the stress-strain relationship of unreinforced walls (MCB, MTB). The numerical results show a good agreement with the experimental results concerning not only at the initial rigidity of the elastic phase but also from the non-linear phase to the post peak corresponding response. However, the values of numerical modeling are somewhat higher than that of the experimental results except the result of maximum strain for the wall sample MTB (θ = 45°). The numerical predicted deformations and peak loads for each panel compared to the average experimental results are shown in Table 11.  For the wall subjected to uniaxial compression perpendicular to bed joint (MCB), the initial cracking occurs along the vertical mortar joint as an initial response. As the vertical loads increase, more cracks occur in the vertical mortar joints of wall from top to down. After that, cracks appear and propagate in the brick units directly as a final response, which causes failure due to vertical cracking of the face shells of the masonry wall ( Figure 13). The same mode of rupture was found experimentally. In MTB wall panels, cracking took place primarily through the bed and head mortar joints, and then extended diagonally. With load increasing, the wall exhibited a gradual increase in the width of the main diagonally-oriented crack, followed by sliding along the diagonal parallel to bed joints. On the other hand, many cracks appeared at the level of brick (Figure 14).  By comparing the numerically predicted damage with the observed failure pattern in the experimental test, tensile damage and thus cracking of the brick are clearly visible in the center of the wall. Otherwise, the tensile stress in the head mortar joint identified in the numerical model has not been observed by experimental test (show Figure 15). Before crushing occurred, the crack pattern appeared in the head joint (initial step), cracks through the mortar joints continued to develop vertically but with a low damage factor value. This model crack pattern was unlike those observed in the experiment test. This can be explained by the existence of micro-crack at the mortar joint but cannot be seen in the experiment test with the overall scale.  Figure 16 displays a comparison of the numerical crack pattern developed in the masonry wall MTB with that which shown by experimental tests. The damage pattern obtained by the developed numerical model is represented in term of compressive damage (DAMAGEC) and tensile damage (DAMAGET) contour plot. This comparison showed us that the crack patterns predicted in the mortar joint by the FE model and crack patterns observed during the experiment resemble each other to a good extent. A good confrontation was found not only for the crack pattern of mortar but also for the brick units. However, sometimes the position of numerical cracks at the brick differs from the experimental results, among them, crushing in the two extremities of the diagonal. This can be explained by the numerical simplification, which comprises considering that all the mortar joints have the same thickness, the same mechanical characteristics, which is not assured experimentally.

Conclusions
Strengthening of structures with NSM CFRP reinforcement is a technique that has attracted considerable attention as a feasible and economical method compared to externally bonded FRP reinforcement. The following conclusions were drawn from the experimental and numerical study: • The increase in the proportion of sand in the mortar (from 3 to 5) led to an increase in the shear strength and the ductility of the masonry panels, especially in the case of reinforcement.

•
The shear strength of unreinforced and reinforced masonry wall panel oriented by 45° was mainly related to the compressive strength of the masonry wall panel.

•
The compressive strength and shear strength of masonry panels were affected by Eunit/Emortar ratio even when this ratio was greater than one.

•
As the horizontal reinforcement restrained the opening of diagonal cracks, the sliding failures along single mortar bed joints caused by the horizontal reinforcement were prevented by the vertical reinforcement.

•
The application of CFRP strips improves and increases the ductility and the bond strength of wall masonry especially in the case of reinforcing of both sides of the panel.

•
The improvement in shear strength of strengthened wall panels with NSM CFRP strips ranged from 1.3 to 2 times (123% to 196%). • All reinforced wall panels showed a substantial increase in the ductility by 2.68 to 4.7 times compared with unreinforced wall panels.

•
The most significant increase in ductility was achieved by vertical NSM-FRP reinforcement in both sides of the wall panel.

•
The Finite Element model proposed in this paper shows a considerable accuracy for the prediction of maximum strength and failure mode.

•
The behavior of the masonry is strongly affected by the behavior of the interface; Coulomb friction criterion are important to simulate correctly the load transmission brick-mortar. • The CDP model can well reflect the damage and the non-linear behavior of brick and mortar. • The developed model has been proved to capture the crack patterns and the stress distribution patterns in both brick and mortar. However, the behaviour of the NSM strengthened masonry wall has not been validated since more experimental data and a particular model would be needed to substantiate this. Therefore, future work should focus on developing numerical models to analyze the wall behavior reinforced with NSM-CFRP under biaxial loading. In addition, several factors such as type of loading and the characteristics of the single component brick and mortar should be considered in future researches in order to develop a suitable and optimal design approach for strengthening of in-plane loaded masonry walls.