The Effectiveness of CFRP- and Auxetic Fabric-Strengthened Brick Masonry under Axial Compression: A Numerical Investigation

Bonded brickwork used for loadbearing walls is widely found in heritage structures worldwide. The evaluation of bonded masonry structures and their strengthening strategies against dynamic actions require appropriate understanding under cyclic loading. Subsequently, a simplified 3D microscale numerical model is developed in this paper to analyse bonded brickwork under cyclic compression. A plasticity-based damage constitutive model to represent damage in masonry bricks under cyclic compression loading was employed, and zero-thickness interfaces were considered with non-linear damage properties to simulate the mechanical behaviour of masonry. A threshold strain level was used to enact the element deletion technique for initiating brittle crack opening in the masonry units. The developed model was validated against the experimental results published by the authors in the past. The models were able to accurately predict the experimental results with an error limit of 10% maximum. Mainly, two types of strengthening materials, possessing (1) high energy absorption characteristics (auxetic fabric) and (2) high strength properties (carbon fibre reinforced polymer composites/CFRP) were employed for damage mitigation under cyclic compression. Results show that the CFRP-strengthened masonry failure was mainly attributed to de-bonding of the CFRP and crushing under compression. However, the auxetic strengthening is shown to significantly minimise the de-bonding phenomenon. Enhanced energy dissipation characteristics with relatively higher ductility (up to ~50%) and reduced damages on the bonded brickwork were observed as compared to the CFRP-strengthened brickwork under cyclic compression loading. Additionally, the auxetic fabric application also increased the compressive resistance of brickwork by 38–60% under monotonic loading, which is comparably higher than with the CFRP strengthening technique.


Introduction
Unreinforced masonry (URM) has been widely used in heritage structures across the world as a construction material that needs continuous monitoring and intervention to withstand any unforeseen actions. The URM walls of these buildings, on many occasions, were constructed of bonded brickwork walls, consisting of two or more bricks layers in thickness. These bonded URM buildings usually suffer severe damage under seismic actions due to brittleness, poor connectivity/tensile strength and low energy dissipation capacity. Therefore, it is necessary to evaluate the optimum retrofitting techniques to minimise damage during seismic actions. Extensive research studies have been performed on the strengthening of URM walls to mitigate damages against seismic and other extreme actions. More than 50,000 material types have been used in retrofitting, strengthening, design and as an energy-absorbing protective layer on the surface of structures [1]. Carbon modelling when the focus of the study is to evaluate only the global behaviour [38]. The micro-modelling approach, although computationally expensive, can be utilised for smallscale studies to understand local failure behaviours such as debonding or delamination [57]. As the main aim of this paper is to predict the effectiveness of CFRP fabric and auxetic fabric application on masonry wallettes under compression loading, the macro-modelling technique provided a simpler model with reasonable failure characterisation to represent the delamination of strengthening fabrics from the masonry wallettes.
Thus, in this study, a simplified finite element (FE)-based micro-modelling technique was developed to analyse the bonded brickwork with and without strengthening materials subjected to monotonic and cyclic compression loads. The formulated FE modelling technique of the bonded brickwork was verified with the experimental results published earlier by the authors [16]. Later, the validated numerical models of the masonry wallettes were expanded with a single layer of CFRP laminates on the surfaces of the wallettes. The model of CFRP-strengthened masonry developed was validated with the experimental results obtained by the authors [17]. Finally, the validated CFRP strengthening layer on the brickwork masonry wallette was replaced with auxetic fabric layers, and then the compression resistance and mechanical characteristics of the wallettes under both monotonic and cyclic compression loadings were examined. Thus, this study contributes to the knowledge on the effectiveness of these two externally bonded strengthening techniques (i.e., CFRP layup and auxetic fabric layup) on masonry under axial compression.

Numerical Modelling of URM Axial Compression
A simplified 3D micro-mechanical numerical modelling technique was adopted in this paper to analyse double brick-bonded masonry wallettes under both monotonic and cyclic compression. An implicit finite-element modelling method was developed to simulate the response of masonry wallettes under axial monotonic and cyclic displacement. The FE model incorporated the material and contact non-linearities, and the input parameters were calibrated using the datasets of the wallette properties obtained through experimental study [16]. Figure 1 shows the FE model developed for the double-brick wallettes constructed of two types of brick units (namely B1 and B2, described in the following sub-section).
stage is where large-scale, discrete element models of staggered composites with statistical variation are applied to capture the failure of composites. Considering all these modelling techniques, the homogenous macro-modelling approach is usually adapted for large-scale structure modelling when the focus of the study is to evaluate only the global behaviour [38]. The micro-modelling approach, although computationally expensive, can be utilised for small-scale studies to understand local failure behaviours such as debonding or delamination [57]. As the main aim of this paper is to predict the effectiveness of CFRP fabric and auxetic fabric application on masonry wallettes under compression loading, the macro-modelling technique provided a simpler model with reasonable failure characterisation to represent the delamination of strengthening fabrics from the masonry wallettes.
Thus, in this study, a simplified finite element (FE)-based micro-modelling technique was developed to analyse the bonded brickwork with and without strengthening materials subjected to monotonic and cyclic compression loads. The formulated FE modelling technique of the bonded brickwork was verified with the experimental results published earlier by the authors [16]. Later, the validated numerical models of the masonry wallettes were expanded with a single layer of CFRP laminates on the surfaces of the wallettes. The model of CFRP-strengthened masonry developed was validated with the experimental results obtained by the authors [17]. Finally, the validated CFRP strengthening layer on the brickwork masonry wallette was replaced with auxetic fabric layers, and then the compression resistance and mechanical characteristics of the wallettes under both monotonic and cyclic compression loadings were examined. Thus, this study contributes to the knowledge on the effectiveness of these two externally bonded strengthening techniques (i.e., CFRP layup and auxetic fabric layup) on masonry under axial compression.

Numerical Modelling of URM Axial Compression
A simplified 3D micro-mechanical numerical modelling technique was adopted in this paper to analyse double brick-bonded masonry wallettes under both monotonic and cyclic compression. An implicit finite-element modelling method was developed to simulate the response of masonry wallettes under axial monotonic and cyclic displacement. The FE model incorporated the material and contact non-linearities, and the input parameters were calibrated using the datasets of the wallette properties obtained through experimental study [16]. Figure 1 shows the FE model developed for the double-brick wallettes constructed of two types of brick units (namely B1 and B2, described in the following subsection). The mortar joints were modelled as zero-thickness interfaces with non-linear damaging properties, while the bricks were enlarged to incorporate half of the thickness of mortar joints all around the bricks, as shown in Figure 1. This type of modelling has been adopted earlier to model bonded brickwork under monotonic compression [64]. The loading plates as illustrated in Figure 1 could move in the direction of compression loading and unloading at the rate of 0.25 mm/min. The experimental programme, model geometry, constitutive material models and validation results are discussed in detail in the following sub-sections.

Experimental Programme
Experimental results for masonry under monotonic and cyclic compression established by the authors [16] were used to validate the applicability of the established numerical modelling technique. The masonry wallettes in this experimental programme were constructed using two types of clay brick units, i.e., B1 and B2 bricks. The difference in the mechanical properties of the B1 and B2 bricks lay mainly in the manufacturing process, i.e., controlled and uncontrolled firing in the kiln. However, the B1 bricks with low strength and modulus of elasticity are very common in many developing countries. In total, 14 masonry wallettes were constructed using B1 and B2 clay bricks with double brickwork thickness to assess the monotonic and cyclic behaviour of bonded masonry. The wallettes were built and tested according to BS EN 1052-1 [65] and loaded under the monotonic and cyclic displacement-controlled protocols. The thickness of the mortar was strictly maintained at 10 mm with the help of an experienced mason. The dimensions and the average compressive strengths of the wallettes are presented in Table 1 along with their coefficients of variation (COV). The experimental representation with the numerical model of the double layer-bonded brickworks constructed of B1 and B2 brick types are shown in Figure 2.
The wallettes were tested using a 1000 kN-capacity servo-controlled Universal Testing Machine (UTM) to assess the monotonic and cyclic compressive strengths and deformation characteristics. Five-millimetre-thick plywood was used between the top and bottom platens of the UTM to uniformly distribute the compression loading over the wallettes. In total, eight displacement transducers were attached to record the axial and lateral deformation in each wallette. Four transducers were attached vertically on both faces (two on each face) of the wallette to record the axial deformation, as shown in Figure 2a,b. Similarly, four remaining transducers (two on each face) were placed horizontally to record lateral deformation on each face of the wallette. A uniform monotonic displacement of 10 mm was applied on the top of the loading plate, with boundary constraints marked with a reference node to simulate the compression tests, as shown in Figure 2. For cyclic testing, a typical cyclic loading protocol in terms of displacement was applied using the same reference node as shown in Figure 3. The cyclic displacement procedure was taken from the load-displacement responses of similar wallettes tested under monotonic loads, as presented in previously published articles [16,17,64]. The elastic, hardening and peak limits were identified from the load-displacement response. The elastic limit point was observed at one-third of the peak load, and the hardening limit was considered as 0.8 times the peak load [16,17,64]. The considered number of cycles within each limit are mentioned in Table 2.  The wallettes were tested using a 1000 kN-capacity servo-controlled Universal Testing Machine (UTM) to assess the monotonic and cyclic compressive strengths and deformation characteristics. Five-millimetre-thick plywood was used between the top and bottom platens of the UTM to uniformly distribute the compression loading over the wallettes. In total, eight displacement transducers were attached to record the axial and lateral deformation in each wallette. Four transducers were attached vertically on both faces (two on each face) of the wallette to record the axial deformation, as shown in Figure  2a,b. Similarly, four remaining transducers (two on each face) were placed horizontally to record lateral deformation on each face of the wallette. A uniform monotonic displacement of 10 mm was applied on the top of the loading plate, with boundary constraints marked with a reference node to simulate the compression tests, as shown in Figure 2. For cyclic testing, a typical cyclic loading protocol in terms of displacement was applied using the same reference node as shown in Figure 3. The cyclic displacement procedure was taken from the load-displacement responses of similar wallettes tested under monotonic loads, as presented in previously published articles [16,17,64]. The elastic, hardening and peak limits were identified from the load-displacement response. The elastic limit point was observed at one-third of the peak load, and the hardening limit was considered as 0.8 times the peak load [16,17,64]. The considered number of cycles within each limit are mentioned in Table 2.     [16,17,64]).

Load-Displacement Response Number of Steps Number of Cycles
Elastic Limit (one-third of peak load) 4 steps 2 cycles at each step Hardening Limit (0.8 times the peak load) 4 Steps 2 cycles at each step Peak Limit 3 Steps 2 cycles at each step

Model Geometry Details
For the simplified 3D micro-mechanical numerical modelling method, the homogenised brick units B1 and B2 were modelled using 3D solid elements (C3D8R) with 8 nodes and 24 degrees of freedom. The unit size was represented by an enlarged brick consisting of full-scale brick enveloped by a half-thickness of the mortar bedding layer, as shown in Figure 2. These masonry bricks were arranged in multiple layers using zero-thickness cohesive interface elements to simulate the bond behaviour under shear, tension, and compression actions. A mesh convergence study was performed for brick units (B1 and B2) to obtain accurate and computationally efficient results. Subsequently, a mesh size of 20 mm × 20 mm × 20 mm was employed for the brick elements after the mesh convergence studies had been carried out.

Material Properties and Constitutive Model
A plasticity-based damage constitutive model to represent the damage in the masonry under cyclic compression loading was employed. The concrete damage plasticity (CDP) model in-built in ABAQUS FE software [66] was used to simulate the mechanical behaviour of smeared brick and mortar together. This model can predict two main failure modes: tensile cracking and compression crushing. Beyond the failure limit, the formation of minor and major cracks was represented with a softening stress-strain response, as shown in Figure 4. After calibrating with the experimental results, a threshold strain level of 0.05 was used to enact the element deletion technique for initiating the brittle crack opening in the masonry units. The element deletion limit value was defined from the studies performed in the past by the authors [36,[66][67][68]. The constitutive model used for the B1 and B2 brick wallettes is shown in Figure 4. The elastic modulus of the B1 series wallette was input as 200 MPa, while for the B2 series wallette, an average value of 995 MPa was adopted from the authors' previous experimental study [16]. The Poisson's ratio was kept as 0.18, as was observed from the experiments. In addition, the CDP model depends on five main parameters: dilation angle (ψ) to account for the volume change caused by the inelastic behaviour of bricks; flow potential eccentricity (e), which defines the rate at which the hyperbolic flow potential approaches its asymptote; the ratio of the biaxial compressive strength of the material to its uniaxial compressive strength ratio (f b0 /f c0 ); shape factor (K c ), which is the ratio of the tensile to the compressive meridian for describing the shape of the yield surface, and viscosity (µ), which represents the relaxation time of the visco-plastic system and is used for the visco-plastic regularisation. More details of these parameters can be found in [69,70]. The assigned CDP parameters are listed in Table 3. the rate at which the hyperbolic flow potential approaches its asymptote; the ratio of the biaxial compressive strength of the material to its uniaxial compressive strength ratio (fb0/fc0); shape factor (Kc), which is the ratio of the tensile to the compressive meridian for describing the shape of the yield surface, and viscosity (µ), which represents the relaxation time of the visco-plastic system and is used for the visco-plastic regularisation. More details of these parameters can be found in [69,70]. The assigned CDP parameters are listed in Table 3.

Zero Thickness Interface: Cohesive Zone Model
The zero-thickness interface between the homogenised units was simulated using a cohesive zone model. The interaction was modelled using the simple bi-linear separation constitutive law as shown in Figure 5. The traction law is graphically represented as the traction force (τ) resulting in the opening displacement (δ). The interface allowed damage initiation followed by damage evolution until the total degradation of surface bonding. The elastic behaviour of the interface relates the normal and the shear (tractions) to the normal and shear displacements (separation) across the interface.

Zero Thickness Interface: Cohesive Zone Model
The zero-thickness interface between the homogenised units was simulated using a cohesive zone model. The interaction was modelled using the simple bi-linear separation constitutive law as shown in Figure 5. The traction law is graphically represented as the traction force (τ) resulting in the opening displacement (δ). The interface allowed damage initiation followed by damage evolution until the total degradation of surface bonding. The elastic behaviour of the interface relates the normal and the shear (tractions) to the normal and shear displacements (separation) across the interface.  The normal tensile and shear , stiffness of the interface are related to the properties of mortar used for connecting the bricks. The lime mortar adopted in this research has low tensile strength and slightly higher shear stiffness. The input values for the interface modelling are provided in Table 4. The damage initiation in the interface was assumed to occur when quadratic traction factors involving the combined normal stress and shear stress , ratios reached the value 1 as shown in Equation (1) The normal tensile (K nn ) and shear (K ss , K ss ) stiffness of the interface are related to the properties of mortar used for connecting the bricks. The lime mortar adopted in this research has low tensile strength and slightly higher shear stiffness. The input values for the interface modelling are provided in Table 4. The damage initiation in the interface was Polymers 2022, 14, 1800 9 of 26 assumed to occur when quadratic traction factors involving the combined normal stress (σ n ) and shear stress (τ s , τ t ) ratios reached the value 1 as shown in Equation (1) where σ 0 n is the limiting (maximum) normal strength and τ 0 s , τ 0 t are the limiting (maximum) shear strength of the interface mortar. The interface damage evolution is expressed in terms of fracture energy release as shown in Equation (2) where G n , G s and G t are work done in normal and shear directions while G c t and G c s are the corresponding maximum fracture energies, which cause failures in normal and shear directions, respectively. The dependency of the fracture energy η is a material parameter proposed by Benzeggagah & Kenane [71] and was taken as 1.5. The Benzeggagah-Kenane [71] fracture criterion defines the critical fracture energies when the deformation and the first and second shear directions are the same. The properties of the zero-thickness interface to account for cohesion, traction separation law and friction are listed in Table 4. Since the masonry assemblages considered to validate the numerical models were made of lime-mortared clay bricks with low tensile strength, the adopted interfacial properties are relatively lower than for the conventional masonry made with cement mortars [70][71][72].

Validation Results
This section presents the validation of the FE model developed for the double-brick masonry wallettes under monotonic and cyclic compression loading using the experimental datasets presented in [16]. The validation includes the failure modes of bonded brickwork wallettes, stress-strain curve validation with the experimental results and the strain variation in the region of significant damage.

Failure Modes
The validation of the FE results under monotonic loading is discussed in this section. However, the detailed study of the monotonic testing results, including the effect of geometry and thickness of wallettes, have been presented in the authors' previous research [64]. The failure modes of the double-thick brickwork models of the B1 and B2 series are presented in terms of maximum principal strain in Figure 6. The numerically predicted failure patterns of the double-thick brick masonry wallettes are shown to be in close agreement with the experimental failure modes (as shown in Figure 6a,b). The FE result also predicted similar vertical cracking on the front face and tensile splitting on the side faces. The tensile cracking phenomena of masonry under axial compression are well understood to be due to incompatible deformation characteristics between unit and mortar and subsequent lateral tensile strains caused by a Poisson's ratio effect which induces vertical tensile splitting cracks [73,74]. The maximum principal strain contour plot of magnitude 0.009 shows that the masonry wallettes of brick type B2 ( Figure 6b) were more brittle than the B1 brick wallettes. These failure strain limit values can also be correlated with the axial stress-strain plot of the wallettes under the monotonic compression load shown in below Section 2.5.2. The failure modes of double-thick bonded brickwork models of the B1 and B2 series under cyclic compression are shown in Figure 7 also in terms of principal strain. The failure patterns observed in the numerical model for both wallettes are shown to be in close agreement with the experimental failure modes. The wallettes under cyclic loading observed similar single splitting cracks through the thickness of the wallette, and parallel tensile cracks were observed along the width of the wallettes. The vertical cracks were mainly observed in the central area and near the edges of the faces. After the peak load, these cracks propagated diagonally and opened, leading to spalling from the edge of the wallettes. The post-peak stage under cyclic loading resulted in wider crack propagation within a short period, associated with the spalling of brick-and-mortar pieces. The failure modes of double-thick bonded brickwork models of the B1 and B2 series under cyclic compression are shown in Figure 7 also in terms of principal strain. The failure patterns observed in the numerical model for both wallettes are shown to be in close agreement with the experimental failure modes. The wallettes under cyclic loading observed similar single splitting cracks through the thickness of the wallette, and parallel tensile cracks were observed along the width of the wallettes. The vertical cracks were mainly observed in the central area and near the edges of the faces. After the peak load, these cracks propagated diagonally and opened, leading to spalling from the edge of the wallettes. The post-peak stage under cyclic loading resulted in wider crack propagation within a short period, associated with the spalling of brick-and-mortar pieces.

Monotonic and Cyclic Stress-Strain Relationship
The axial compressive stress vs strain curves obtained from the numerical models and experimental data are presented in Figure 8 for B1 and B2 brick wallettes. To plot the numerical stress-strain curves, the axial (vertical) deflection from the nodes equivalent to the gauge length used in the experiments [16] was used to calculate the axial strain, while the reaction force was used to find the axial stress by dividing it with the cross-section area of the respective specimens. The lateral strains represent the horizontal strain measurement on the front face of the wallette samples, as indicated in Figure 2. The predicted stress-strain curves observed from the numerical models match reasonably well with the experimental curves for both types of brick, B1 and B2. Figure 9 shows the cyclic stress-strain curves of the bonded brickwork of the B1 and B2 series. The axial stress-strain curves derived from the numerical models were also compared with the experimental curves for both brickwork specimens. From the cyclic envelope stress-strain curves, the elastic moduli were derived, which varied from 135 MPa for B1 series wallettes to 506 MPa for B2 series wallettes. It can be clearly seen that the non-linear behaviour under cyclically loaded wallettes started nearly around 40% to 50% of peak compressive stress, which was associated with the initiation of cracks in the wallettes. The accumulation of non-reversible axial strains in the masonry wallettes resulted in the degradation of stiffness in each cycle, especially after about 40% of peak stress, for both B1 and B2 series wallettes. This degradation indicates the gradual damage that has occurred in each cycle in the masonry wallettes. Therefore, the pattern of cyclic stress-strain curves of wallettes shown in Figure 9 justifies the overall decrease in the compressive strength of masonry wallettes under cyclic compression when compared to monotonic loading [16].

Monotonic and Cyclic Stress-Strain Relationship
The axial compressive stress vs strain curves obtained from the numerical mod and experimental data are presented in Figure 8 for B1 and B2 brick wallettes. To plot t numerical stress-strain curves, the axial (vertical) deflection from the nodes equivalent the gauge length used in the experiments [16] was used to calculate the axial strain, wh the reaction force was used to find the axial stress by dividing it with the cross-secti area of the respective specimens. The lateral strains represent the horizontal strain me urement on the front face of the wallette samples, as indicated in Figure 2. The predict stress-strain curves observed from the numerical models match reasonably well with t experimental curves for both types of brick, B1 and B2.

Monotonic and Cyclic Stress-Strain Relationship
The axial compressive stress vs strain curves obtained from the numerical models and experimental data are presented in Figure 8 for B1 and B2 brick wallettes. To plot the numerical stress-strain curves, the axial (vertical) deflection from the nodes equivalent to the gauge length used in the experiments [16] was used to calculate the axial strain, while the reaction force was used to find the axial stress by dividing it with the cross-section area of the respective specimens. The lateral strains represent the horizontal strain measurement on the front face of the wallette samples, as indicated in Figure 2. The predicted stress-strain curves observed from the numerical models match reasonably well with the experimental curves for both types of brick, B1 and B2.  lettes. The accumulation of non-reversible axial strains in the masonry wallettes resulted in the degradation of stiffness in each cycle, especially after about 40% of peak stress, for both B1 and B2 series wallettes. This degradation indicates the gradual damage that has occurred in each cycle in the masonry wallettes. Therefore, the pattern of cyclic stressstrain curves of wallettes shown in Figure 9 justifies the overall decrease in the compressive strength of masonry wallettes under cyclic compression when compared to monotonic loading [16]. Additionally, Figure 9 shows that the B2 series wallettes have significantly lower axial and lateral deformability than the B1 series wallettes. This phenomenon is due to high strength and lower deformability properties of B2 bricks, which influenced the overall stress-strain behaviour of the masonry. Hence, the accumulation of strain and the time history of the strains is significant for the compressive strength of the masonry wallettes, as well as for defining mitigating strategies under cyclic loading.

Numerical Modelling of Strengthened Brick Wallettes under Axial Compression
This section entails the mitigating methods to enhance the performance of brick masonry under axial compression loading. The technique involves a layup of CFRP laminates on the faces of the bonded brickwork wallettes, as shown in Figure 10, and its effectiveness was studied using numerical modelling. The fabric material properties and the adhesion properties of the epoxy resin are discussed. This section also describes the adhesive interface modelling between the surface of masonry and CFRP-composite rendering under monotonic and cyclic compression testing. Moreover, the validation of the modelling technique for CFRP layup with the experimental results under monotonic loading, and then predicting the behaviour of CFRP layup for the same kind of masonry wallettes under cyclic compression testing, are presented in this section. Additionally, Figure 9 shows that the B2 series wallettes have significantly lower axial and lateral deformability than the B1 series wallettes. This phenomenon is due to high strength and lower deformability properties of B2 bricks, which influenced the overall stress-strain behaviour of the masonry. Hence, the accumulation of strain and the time history of the strains is significant for the compressive strength of the masonry wallettes, as well as for defining mitigating strategies under cyclic loading.

Numerical Modelling of Strengthened Brick Wallettes under Axial Compression
This section entails the mitigating methods to enhance the performance of brick masonry under axial compression loading. The technique involves a layup of CFRP laminates on the faces of the bonded brickwork wallettes, as shown in Figure 10, and its effectiveness was studied using numerical modelling. The fabric material properties and the adhesion properties of the epoxy resin are discussed. This section also describes the adhesive interface modelling between the surface of masonry and CFRP-composite rendering under monotonic and cyclic compression testing. Moreover, the validation of the modelling technique for CFRP layup with the experimental results under monotonic loading, and then predicting the behaviour of CFRP layup for the same kind of masonry wallettes under cyclic compression testing, are presented in this section.

CFRP Fabric: Geometric Detailing
The fabric geometry was modelled using the conventional planar 3D shell element (S4R) with six DOF at each node [38,54,75,76], as shown in Figure 10. The fabric was assigned a homogenous laminar property at each layer using five integration points through

CFRP Fabric: Geometric Detailing
The fabric geometry was modelled using the conventional planar 3D shell element (S4R) with six DOF at each node [38,54,75,76], as shown in Figure 10. The fabric was assigned a homogenous laminar property at each layer using five integration points through the thickness. A single layer of fabric layup on the front and rear face of the masonry wallettes was employed. The dimension of the CFRP fabric used for bonded bricks wallettes of the B1 and B2 series were 410 mm (L) × 730 mm (H) × 1 mm (T) and 430 mm (L) × 750 mm (H) × 1 mm (T), respectively, as shown in Figure 10.

CFRP Fabric: Material Properties
The Hashin's damage model available in ABAQUS [66] was used to simulate the characteristics of the CFRP fabric [77]. Hashin's damage initiation and evolution laws were used to define the progressive damage of CFRP fabrics. Equations (3)-(6) were used to define the damage initiation and the damage evolution, which were adapted from the study by Asad et al. [38]. The maximum degradation was considered for the fabric when the damage variable D reaches a magnitude of 1. The element deletion technique was only activated when the damage variable achieved this value. A failure is initiated in the fabric laminate when one of the following indices is greater than or equal to 1: fibre tension, fibre compression, matrix tension and matrix compression, as represented by Equations (3)-(6), respectively.
Matrix Tension (σ yy ≥ 0) Matrix Compression (σ yy < 0) where X T denotes the tensile strength in the fibre direction, X c denotes the compressive strength in the fibre direction (which is assumed to be negligible in the analyses), Y T and Y C are the tensile and compressive strengths in the direction perpendicular to the fibres (Y C also assumed negligible in the analyses), S T denotes the transverse shear strength, S L denotes the longitudinal shear strength and α h is the coefficient that determines the contribution of shear stress to the fibre tensile failure initiation criteria. Furthermore, σ xx , σ yy & τ xy are the components of the effective stress tensor. The input properties of the CFRP laminate were taken from past experimental studies [73,74,78,79]. The input properties of the CFRP fabric are listed in Tables 5-7. It has to be mentioned that the compressive strength of CFRP fabric is not explicitly considered for design or analyses. Although FRP fabrics are compressed in the strengthened masonry, the FRP failures were mainly dominated by tensile rupture or delamination due to the dilation of the strengthened assembly under axial compression.

CFRP Fabric Brick Wallettes: Interface Modelling
From the literature [17,37,38], delamination/de-bonding was observed between CFRP surfaces and the masonry wallettes. Such delamination or de-bonding at the interfaces was mainly due to the adhesive failure between the epoxy-coated CFRP and masonry substrates. In numerical models, the de-bonding was simulated using the cohesive zone interface model. The interaction features in the cohesive model were evaluated using the same bi-linear traction separation constitutive law used for zero-thickness mortar bedding, as shown in Figure 5b. The cohesive properties between the CFRP fabric and the masonry wallettes are summarised in Table 8.

Validation of the CFRP Fabric-Strengthened Brick Wallettes under Monotonic Compression Loading
The developed numerical model for CFRP-strengthened brickwork wallettes was validated with the experimental results presented in Thamboo et al. [17]. The experimental programme involved testing of six CFRP-strengthened masonry wallettes under monotonic axial static loading. The masonry wallettes considered for the testing were made from the same types of bricks (B1 and B2 clay brick units) as discussed in Section 2. Hence, the constitutive material model was duplicated, as explained in detail in Section 2.3 for the masonry wallettes under monotonic and cyclic loading. However, the application of the CFRP and the testing procedures were explained in detail in Thamboo et al. [17]. The unidirectional CFRP sheet layups were applied using saturated epoxy resin as an adhesive layer on the front and the back surface of the masonry wallettes, as shown in Figure 11. It must be mentioned that the unidirectional CFRP sheets were applied parallel to the bed joints, a technique which is commonly employed to improve the shear/tensile resistance of masonry walls. Additionally, a confining of walls with fabrics as is usually carried out for column elements is not possible due to their considerably larger length than columns and connections with crossing walls. The panels were then tested under axial compression to verify the influence of the conventional strengthening method (with unidirectional CFRP parallel to bed joints) on the compression strength characteristics of masonry, which is a prime input parameter for the macro-modelling of masonry structures. The testing was carried out at a constant displacement rate of 0.25 mm/min using the same UTM of 1000 kN capacity to capture the complete post-crack behaviour of masonry wallettes under axial compression. The numerical model developed corresponding to the experimental test setup is shown in Figure 10. Displacement transducers were attached to the wallettes to capture the vertical deformation at mid-height of the wallettes, as shown in Figure 11c. The typical failure modes of monotonically loaded wallettes of B1 and B2 brick masonry strengthened with CFRP sheets are shown in Figure 12. Due to the CFRP sheets attached to the front and back face of the masonry wallettes, no visible cracks were noted in the experimental test, as shown in Figure 12. Previous similar research studies on multiple thin CFRP strips layups on masonry walls [80,81] were verified to confirm the failure modes developed on the side and front/rear face of the wallettes. The failure modes on the sides with visible CFRP deformation buckling are shown in Figure 13a. A premature interaction failure at the side edge between the CFRP strip and the masonry surface is shown in Figure 13b.  The testing was carried out at a constant displacement rate of 0.25 mm/min using the same UTM of 1000 kN capacity to capture the complete post-crack behaviour of masonry wallettes under axial compression. The numerical model developed corresponding to the experimental test setup is shown in Figure 10. Displacement transducers were attached to the wallettes to capture the vertical deformation at mid-height of the wallettes, as shown in Figure 11c. The typical failure modes of monotonically loaded wallettes of B1 and B2 brick masonry strengthened with CFRP sheets are shown in Figure 12. Due to the CFRP sheets attached to the front and back face of the masonry wallettes, no visible cracks were noted in the experimental test, as shown in Figure 12. Previous similar research studies on multiple thin CFRP strips layups on masonry walls [80,81] were verified to confirm the failure modes developed on the side and front/rear face of the wallettes. The failure modes on the sides with visible CFRP deformation buckling are shown in Figure 13a. A premature interaction failure at the side edge between the CFRP strip and the masonry surface is shown in Figure 13b. The testing was carried out at a constant displacement rate of 0.25 mm/min using the same UTM of 1000 kN capacity to capture the complete post-crack behaviour of masonry wallettes under axial compression. The numerical model developed corresponding to the experimental test setup is shown in Figure 10. Displacement transducers were attached to the wallettes to capture the vertical deformation at mid-height of the wallettes, as shown in Figure 11c. The typical failure modes of monotonically loaded wallettes of B1 and B2 brick masonry strengthened with CFRP sheets are shown in Figure 12. Due to the CFRP sheets attached to the front and back face of the masonry wallettes, no visible cracks were noted in the experimental test, as shown in Figure 12. Previous similar research studies on multiple thin CFRP strips layups on masonry walls [80,81] were verified to confirm the failure modes developed on the side and front/rear face of the wallettes. The failure modes on the sides with visible CFRP deformation buckling are shown in Figure 13a. A premature interaction failure at the side edge between the CFRP strip and the masonry surface is shown in Figure 13b.  The testing was carried out at a constant displacement rate of 0.25 mm/min using the same UTM of 1000 kN capacity to capture the complete post-crack behaviour of masonry wallettes under axial compression. The numerical model developed corresponding to the experimental test setup is shown in Figure 10. Displacement transducers were attached to the wallettes to capture the vertical deformation at mid-height of the wallettes, as shown in Figure 11c. The typical failure modes of monotonically loaded wallettes of B1 and B2 brick masonry strengthened with CFRP sheets are shown in Figure 12. Due to the CFRP sheets attached to the front and back face of the masonry wallettes, no visible cracks were noted in the experimental test, as shown in Figure 12. Previous similar research studies on multiple thin CFRP strips layups on masonry walls [80,81] were verified to confirm the failure modes developed on the side and front/rear face of the wallettes. The failure modes on the sides with visible CFRP deformation buckling are shown in Figure 13a. A premature interaction failure at the side edge between the CFRP strip and the masonry surface is shown in Figure 13b.  The failure modes obtained from the FE model of both types of brick (B1 and B2) masonry wallettes strengthened with CFRP fabric are shown in Figure 14. The strain contour plots show that the maximum stress concentration occurred at the vertical edge of the wall, which resulted in the vertical cracks due to splitting. The magnitude of the strain observed in the B2 brick-type masonry wallettes was smaller than in the B1-brick type masonry wallettes, due to premature failure observed at the interface. This was mainly due to the brittle nature of the B2 bricks, which failed at a lower strain, compared to the B1 brick-type masonry wallettes. The failure modes obtained from the FE model of both types of brick (B1 and B2) masonry wallettes strengthened with CFRP fabric are shown in Figure 14. The strain contour plots show that the maximum stress concentration occurred at the vertical edge of the wall, which resulted in the vertical cracks due to splitting. The magnitude of the strain observed in the B2 brick-type masonry wallettes was smaller than in the B1-brick type masonry wallettes, due to premature failure observed at the interface. This was mainly due to the brittle nature of the B2 bricks, which failed at a lower strain, compared to the B1 brick-type masonry wallettes. The numerically predicted stress-strain curves were also compared with the experimental curves for the CFRP-strengthened masonry wallettes, as shown in Figure 15. For plotting the numerical stress-strain curves, the axial (vertical) deflection from the nodes for the same gauge length that was used in the experiments [17] was used to calculate the axial strain, while the reaction force on the loading plate was used to find the axial stress by dividing it with the cross-section area of the respective specimens. The lateral strains shown in Figure 15 represent horizontal strain measured on the front face of the CFRP sheet samples. The same gauge length as in the experiments was used to calculate the horizontal lateral strains from the models. The numerical model results matched the experimental curves reasonably well, despite the variation in brick strengths. The stressstrain curve of the CFRP-strengthened wallette of B2 bricks, as shown in Figure 15b, confirmed the sudden or premature failure in these wallettes due to the brittleness of the B2 brick types. The numerically predicted stress-strain curves were also compared with the experimental curves for the CFRP-strengthened masonry wallettes, as shown in Figure 15. For plotting the numerical stress-strain curves, the axial (vertical) deflection from the nodes for the same gauge length that was used in the experiments [17] was used to calculate the axial strain, while the reaction force on the loading plate was used to find the axial stress by dividing it with the cross-section area of the respective specimens. The lateral strains shown in Figure 15 represent horizontal strain measured on the front face of the CFRP sheet samples. The same gauge length as in the experiments was used to calculate the horizontal lateral strains from the models. The numerical model results matched the experimental curves reasonably well, despite the variation in brick strengths. The stress-strain curve of the CFRP-strengthened wallette of B2 bricks, as shown in Figure 15b, confirmed the sudden or premature failure in these wallettes due to the brittleness of the B2 brick types.

Prediction of Cyclic Behaviour of CFRP Fabric-Strengthened Brickwork Wallettes
The numerical models were modified with the cyclic loading protocol (as shown in Figure 3) to study the effectiveness of strengthening under cyclic compression. This subsection presents the performance of the CFRP fabric-strengthened masonry wallettes under cyclic compression loading. The failure modes of the CFRP fabric-strengthened masonry wallettes constructed of B1 and B2 brick types are shown in Figure 16. The failure modes are presented in terms of maximum principal strain, and the images shown in Figure 16 are at the ultimate failure stage of the CFRP-strengthened masonry wallettes. Debonding in the CFRP-strengthened masonry wallettes of B2 type brickwork was observed at a magnitude 83% lower than the B1 type brickwork-constructed masonry wallettes. This reduction in strain in B2 bricks wallettes has already been discussed in the previous section under monotonic static loading. The CFRP-strengthened masonry wallettes of B2 bricks showed a 67% decrease in strain compared to B1 brick wallettes strengthened with CFRP fabric. Based on these results, it can be said that the formulated modelling technique appropriately predicted the failure pattern of the CFRP-strengthened masonry wallettes using different types of brick units.

Prediction of Cyclic Behaviour of CFRP Fabric-Strengthened Brickwork Wallettes
The numerical models were modified with the cyclic loading protocol (as shown in Figure 3) to study the effectiveness of strengthening under cyclic compression. This subsection presents the performance of the CFRP fabric-strengthened masonry wallettes under cyclic compression loading. The failure modes of the CFRP fabric-strengthened masonry wallettes constructed of B1 and B2 brick types are shown in Figure 16. The failure modes are presented in terms of maximum principal strain, and the images shown in Figure 16 are at the ultimate failure stage of the CFRP-strengthened masonry wallettes. Debonding in the CFRP-strengthened masonry wallettes of B2 type brickwork was observed at a magnitude 83% lower than the B1 type brickwork-constructed masonry wallettes. This reduction in strain in B2 bricks wallettes has already been discussed in the previous section under monotonic static loading. The CFRP-strengthened masonry wallettes of B2 bricks showed a 67% decrease in strain compared to B1 brick wallettes strengthened with CFRP fabric. Based on these results, it can be said that the formulated modelling technique appropriately predicted the failure pattern of the CFRP-strengthened masonry wallettes using different types of brick units. The stress-strain curves of the CFRP-strengthened cyclically loaded wallettes are presented in Figure 17. The cyclic hysteresis stress-strain curves predicted that the CFRPstrengthened masonry wallettes behaved much more stiffly than the control (un-strengthened) masonry under cyclic loads. However, the CFRP-strengthened masonry wallettes of B2 brick showed a premature failure due to the brittle and stiff properties of the B2 bricks, as shown in Figure 17b. The reason could be the early de-bonding of the fabric The stress-strain curves of the CFRP-strengthened cyclically loaded wallettes are presented in Figure 17. The cyclic hysteresis stress-strain curves predicted that the CFRPstrengthened masonry wallettes behaved much more stiffly than the control (un-strengthened) masonry under cyclic loads. However, the CFRP-strengthened masonry wallettes of B2 brick showed a premature failure due to the brittle and stiff properties of the B2 bricks, as shown in Figure 17b. The reason could be the early de-bonding of the fabric from the surface of the masonry wallettes. This phenomenon was further investigated for both types of brick wallettes by studying the strain levels at the interface of fabric and masonry wallette. The stress-strain curves of the CFRP-strengthened cyclically loaded wallettes are presented in Figure 17. The cyclic hysteresis stress-strain curves predicted that the CFRPstrengthened masonry wallettes behaved much more stiffly than the control (un-strengthened) masonry under cyclic loads. However, the CFRP-strengthened masonry wallettes of B2 brick showed a premature failure due to the brittle and stiff properties of the B2 bricks, as shown in Figure 17b. The reason could be the early de-bonding of the fabric from the surface of the masonry wallettes. This phenomenon was further investigated for both types of brick wallettes by studying the strain levels at the interface of fabric and masonry wallette. A simplified illustration of the axial strain distribution with respect to time on the surfaces of the masonry wallettes and the CFRP sheet under cyclic loading is plotted in Figure 18. The axial strain obtained from the numerical model was plotted at mid-height of the wallettes. The difference in strain levels of masonry and CFRP sheet in Figure 18a,b is clear evidence of debonding at the interface between the CFRP fabric and the masonry wallettes. It is also proven from Figure 18 that the damage in B2 brick wallettes occurred at a much lower strain compared to the B1 bricks wallettes, due to the brittle nature of the B2 bricks. A simplified illustration of the axial strain distribution with respect to time on the surfaces of the masonry wallettes and the CFRP sheet under cyclic loading is plotted in Figure 18. The axial strain obtained from the numerical model was plotted at mid-height of the wallettes. The difference in strain levels of masonry and CFRP sheet in Figure 18a,b is clear evidence of debonding at the interface between the CFRP fabric and the masonry wallettes. It is also proven from Figure 18 that the damage in B2 brick wallettes occurred at a much lower strain compared to the B1 bricks wallettes, due to the brittle nature of the B2 bricks.

Alternative Strengthening of Masonry Using Auxetic Fabrics
From the results of Section 3, it can be easily concluded that the CFRP strengthening of the masonry wallettes was not fully effective due to the early debonding of the CFRP fabric from the surface of the masonry wallettes under monotonic and cyclic compression loading. Alternative to these failure modes, it was observed from the literature [37,54] that an innovative auxetic fabric performed reasonably well when compared with CFRPs due to better bonding properties of the auxetic fabric. The same auxetic fabric was employed here to evaluate the behaviour of the auxetic fabric layup on the surface of the brick masonry wallettes under monotonic and cyclic loading.

Alternative Strengthening of Masonry Using Auxetic Fabrics
From the results of Section 3, it can be easily concluded that the CFRP strengthening of the masonry wallettes was not fully effective due to the early debonding of the CFRP fabric from the surface of the masonry wallettes under monotonic and cyclic compression loading. Alternative to these failure modes, it was observed from the literature [37,54] that an innovative auxetic fabric performed reasonably well when compared with CFRPs due to better bonding properties of the auxetic fabric. The same auxetic fabric was employed here to evaluate the behaviour of the auxetic fabric layup on the surface of the brick masonry wallettes under monotonic and cyclic loading.

Auxetic Fabric-Strengthened Masonry Wallettes under Compression Loading
The numerical model of auxetic fabric-strengthened masonry wallettes was developed by replacing the CFRP fabric with the auxetic fabric, which was also modelled using Hashin's damage criteria [77] inbuilt in ABAQUS [66]. The properties of auxetic fabric are presented in Tables 9-11, which were taken from Asad et al. [37,38,76], and were validated earlier in [38] with experimental results [52]. It can be mentioned that for the masonry applications under axial compression, the fabrics mostly fail by tensile rupture or delamination rather than in compression, due to dilation of the strengthened masonry. Furthermore, the interface modelling and cohesive properties were assumed to be the same as used to validate the CFRP-strengthened masonry wallettes for both types of brick units, due to the same epoxy used for the bonding.   The failure modes of the auxetic fabric-strengthened masonry wallettes under monotonic loading for both types of bricks are shown in Figure 19 in terms of maximum principal strain. The predicted numerical failure patterns on the surface of the auxetic fabric, as shown in Figure 19a, imply effectiveness over the CFRP-strengthened masonry wallettes ( Figure 14) under monotonic loading. The uniform distribution of the principal strains on the surface of the auxetic fabric compared to the CFRP-strengthened masonry wallettes increased the ductility of the wallettes. The maximum increase in the ductility (corresponding to an increase in ultimate failure strain) was, for the B1 brickwork wallettes, due to the non-debonding of the fabric. However, compared to the B1 brickwork wallettes, the B2 brick wallettes, as shown in Figure 19b, exhibited lower ductility due to premature debonding of the auxetic fabric. However, the results show that B2 brick assemblies were still more efficient when compared to the CFRP-strengthened B2 brick wallettes.
The numerically derived axial stress-strain curves of auxetic-strengthened wallettes are shown in Figure 20, and are compared to the CFRP-strengthened wallette results. Figure 20a shows an enhancement in strength with a significant increase in strain for the auxeticstrengthened masonry wallettes constructed of B1 type bricks. The enhancement is mainly due to the non-debonding behaviour of auxetic fabric on the surface of masonry wallettes. In comparison, the B2 type auxetic-strengthened masonry wallettes exhibited a slight increase in strength with a 67% increment in strain compared to the CFRP-strengthened masonry wallettes. Again, the reason for the slight increment in strength could lie in the stiffness properties of the B2 bricks. However, the bonding between the auxetic fabric and B2 brick masonry wallette surfaces increased due to the negative Poisson's ratio of the fabric. cipal strain. The predicted numerical failure patterns on the surface of the auxetic fabric, as shown in Figure 19a, imply effectiveness over the CFRP-strengthened masonry wallettes ( Figure 14) under monotonic loading. The uniform distribution of the principal strains on the surface of the auxetic fabric compared to the CFRP-strengthened masonry wallettes increased the ductility of the wallettes. The maximum increase in the ductility (corresponding to an increase in ultimate failure strain) was, for the B1 brickwork wallettes, due to the non-debonding of the fabric. However, compared to the B1 brickwork wallettes, the B2 brick wallettes, as shown in Figure 19b, exhibited lower ductility due to premature debonding of the auxetic fabric. However, the results show that B2 brick assemblies were still more efficient when compared to the CFRP-strengthened B2 brick wallettes.

Auxetic Fabric-Strengthened Masonry Wallettes under Cyclic Compression
The auxetic fabric-strengthened masonry wallette models were also analysed under the cyclic loading protocol (as shown in Figure 3) to study the effectiveness of auxetic strengthening under cyclic compression. The failure modes of the auxetic fabric-strengthened masonry wallettes constructed of B1 and B2 brick types are shown in Figure 21. The failure modes are presented in terms of maximum principal strain, and the images shown in Figure 21 are at the ultimate failure stage. Figure 21a shows the initiation of the debonding of the auxetic-strengthened masonry wallettes at the maximum principal strain of 0.055, which is 67% larger than with the CFRP sheet ( Figure 16a). Additionally, a reduction of 60% in the principal strain was observed in the B2-type wallettes when compared with the B1-type brickwork masonry wallettes (as shown in Figure 21a,b).

Auxetic Fabric-Strengthened Masonry Wallettes under Cyclic Compression
The auxetic fabric-strengthened masonry wallette models were also analysed under the cyclic loading protocol (as shown in Figure 3) to study the effectiveness of auxetic strengthening under cyclic compression. The failure modes of the auxetic fabric-strengthened masonry wallettes constructed of B1 and B2 brick types are shown in Figure 21. The failure modes are presented in terms of maximum principal strain, and the images shown in Figure 21 are at the ultimate failure stage. Figure 21a shows the initiation of the debonding of the auxetic-strengthened masonry wallettes at the maximum principal strain of 0.055, which is 67% larger than with the CFRP sheet ( Figure 16a). Additionally, a reduction of 60% in the principal strain was observed in the B2-type wallettes when compared with the B1-type brickwork masonry wallettes (as shown in Figure 21a,b).
The stress-strain curves of the cyclically loaded auxetic fabric-strengthened wallettes are presented in Figure 22. The cyclic hysteresis stress-strain curves predicted that the auxetic-strengthened masonry wallettes behaved in a softer and more ductile way than the control and CFRP-strengthened masonry wallettes (presented in Figure 17). Additionally, as expected, the auxetic-strengthened masonry wallettes of B2 brick showed an early failure compared to the B1 brick wallettes, due to the brittle and stiff nature of the brick, as shown in Figure 22b. However, compared to CFRP-strengthened brickwork, the auxeticstrengthened masonry wallettes showed a significant improvement in terms of ductility. It can be noted that there is an increase of 50% in compression strains from un-strengthened to auxetic-strengthened wallettes of the B1 brick type, as shown in Figure 22a. In contrast, when compared with the CFRP-strengthened wallettes (presented in Figure 17), the same compression strain showed a significant increase of about 80%. A similar observation of enhancement in compression strain was observed for the auxetic-strengthened masonry wallettes of B2 brick types when compared to the corresponding un-strengthened masonry wallettes. Hence, to confirm the non-debonding characteristics of the auxetic fabric, an investigation was performed. The stress-strain curves of the cyclically loaded auxetic fabric-strengthened wallettes are presented in Figure 22. The cyclic hysteresis stress-strain curves predicted that the auxetic-strengthened masonry wallettes behaved in a softer and more ductile way than the control and CFRP-strengthened masonry wallettes (presented in Figure 17). Additionally, as expected, the auxetic-strengthened masonry wallettes of B2 brick showed an early failure compared to the B1 brick wallettes, due to the brittle and stiff nature of the brick, as shown in Figure 22b. However, compared to CFRP-strengthened brickwork, the auxetic-strengthened masonry wallettes showed a significant improvement in terms of ductility. It can be noted that there is an increase of 50% in compression strains from unstrengthened to auxetic-strengthened wallettes of the B1 brick type, as shown in Figure  22a. In contrast, when compared with the CFRP-strengthened wallettes (presented in Figure 17), the same compression strain showed a significant increase of about 80%. A similar observation of enhancement in compression strain was observed for the auxetic-strengthened masonry wallettes of B2 brick types when compared to the corresponding unstrengthened masonry wallettes. Hence, to confirm the non-debonding characteristics of the auxetic fabric, an investigation was performed. A simplified illustration of the axial strain distribution with respect to time on the surfaces of the masonry wallettes and the auxetic fabric sheet is plotted in Figure 23. The axial strain was plotted at mid-height of the wallettes. Figure 23a,b show that the strain distribution on the surface of the auxetic fabric and at the masonry surface completely overlapped with each other, with variation in time. The plotted strain distribution time history provides clear evidence of non-debonding at the interface between the auxetic fabric and the masonry wallettes. However, comparing Figure 23a,b, the damage in the auxetic-strengthened masonry wallettes constructed of B2 brick occurred at lower strain (0.005) compared to the B1 bricks wallettes (strain = 0.02) due to the stiff and brittle nature of the B2 bricks. A simplified illustration of the axial strain distribution with respect to time on the surfaces of the masonry wallettes and the auxetic fabric sheet is plotted in Figure 23. The axial strain was plotted at mid-height of the wallettes. Figure 23a,b show that the strain distribution on the surface of the auxetic fabric and at the masonry surface completely overlapped with each other, with variation in time. The plotted strain distribution time history provides clear evidence of non-debonding at the interface between the auxetic fabric and the masonry wallettes. However, comparing Figure 23a,b, the damage in the auxetic-strengthened masonry wallettes constructed of B2 brick occurred at lower strain (0.005) compared to the B1 bricks wallettes (strain = 0.02) due to the stiff and brittle nature of the B2 bricks.

Summary and Conclusions
A simplified micro-modelling method of un-strengthened and CFRP-/auxetic fabricstrengthened masonry was developed for predicting the behaviour of bonded brickwork wallettes under monotonic and cyclic compression loading. The element deletion technique was applied for the 3D solid and shell elements to allow for damage in the brick units and the fabric sheets, respectively. The material and contact/interface nonlinearities were incorporated into the simplified micro-model developed. The modelling methodologies were validated using different experimental results of control/un-strengthened and CFRP-strengthened masonry wallettes under cyclic and monotonic loading. The validated models were used to systematically predict the monotonic and cyclic compression  A simplified illustration of the axial strain distribution with respect to time on the surfaces of the masonry wallettes and the auxetic fabric sheet is plotted in Figure 23. The axial strain was plotted at mid-height of the wallettes. Figure 23a,b show that the strain distribution on the surface of the auxetic fabric and at the masonry surface completely overlapped with each other, with variation in time. The plotted strain distribution time history provides clear evidence of non-debonding at the interface between the auxetic fabric and the masonry wallettes. However, comparing Figure 23a,b, the damage in the auxetic-strengthened masonry wallettes constructed of B2 brick occurred at lower strain (0.005) compared to the B1 bricks wallettes (strain = 0.02) due to the stiff and brittle nature of the B2 bricks.

Summary and Conclusions
A simplified micro-modelling method of un-strengthened and CFRP-/auxetic fabricstrengthened masonry was developed for predicting the behaviour of bonded brickwork wallettes under monotonic and cyclic compression loading. The element deletion technique was applied for the 3D solid and shell elements to allow for damage in the brick units and the fabric sheets, respectively. The material and contact/interface nonlinearities were incorporated into the simplified micro-model developed. The modelling methodologies were validated using different experimental results of control/un-strengthened and CFRP-strengthened masonry wallettes under cyclic and monotonic loading. The validated models were used to systematically predict the monotonic and cyclic compression

Summary and Conclusions
A simplified micro-modelling method of un-strengthened and CFRP-/auxetic fabricstrengthened masonry was developed for predicting the behaviour of bonded brickwork wallettes under monotonic and cyclic compression loading. The element deletion technique was applied for the 3D solid and shell elements to allow for damage in the brick units and the fabric sheets, respectively. The material and contact/interface nonlinearities were incorporated into the simplified micro-model developed. The modelling methodologies were validated using different experimental results of control/un-strengthened and CFRPstrengthened masonry wallettes under cyclic and monotonic loading. The validated models were used to systematically predict the monotonic and cyclic compression behaviour of auxetic fabric-strengthened masonry under monotonic and cyclic compression. Following are the main conclusions drawn from this research: The simplified micro-modelling technique for the masonry wallettes constructed with different brick types was found to be computationally efficient without compromising the accuracy in predicting compression behaviour under the monotonic and cyclic compression loading protocols. Due to the de-bonding nature of failure in CFRP-strengthened masonry wallettes, the increases in the compressive resistance and ductility of masonry wallettes were limited. This phenomenon was more pronounced in the cyclic loading conditions, where with progressive damage, premature compression failures were observed due to early de-bonding. The CFRP strengthening increased the compressive resistance of the masonry wallettes compared to the un-strengthened masonry wallettes by about 10-20%. In contrast, the auxetic fabric application increased the compressive resistance of brickwork by 38-60% under monotonic loading. In cyclic compression loading, the auxetic fabric application to the surface of the masonry wallettes increased the ultimate strains without compromising the strength, compared to the un-strengthened masonry wallettes, which significantly enhanced the ductility (up to~50%) of the masonry wallettes.
From the evaluated results of this study, it can be concluded that the application of the auxetic fabric on the surface of the masonry wallettes is more efficient than conventional CFRP fabric in enhancing compressive resistance and ductility under monotonic and cyclic compression loading. This enhancement was achieved due to the negative Poisson's ratio effect of the auxetic fabric, which ultimately improves the bonding nature. However, more experimental and numerical studies are needed to confirm the parametric analyses and to further experimentally verify the effectiveness of auxetic fabric application in strengthening masonry under different stress states in order to comprehend the complete characteristics. Data Availability Statement: The raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.