Abstract
The primary mode of premature failure for near-surface mounted strengthened beams is the concrete cover separation. Due to its complexity, most of the prediction methods for concrete cover separation tend to be empirical based, which can limit their usage to specific near-surface mounted strengthening configurations. In response to that, this article presents a mechanics-based design which uses the moment-rotation approach and the global energy balance approach which is less reliant on empirical formulations, as the mechanics of reinforced concrete beam such as tension stiffening and propagation of concrete cover separation debonding crack are directly simulated rather than empirically derived. The proposed design procedure was validated against published experimental results of reinforced concrete beams strengthened with near-surface mounted carbon fibre–reinforced polymer bars, near-surface mounted carbon fibre–reinforced polymer strips or side-near-surface mounted carbon fibre–reinforced polymer bars and show good accuracy. As it is less reliant on empirical formulations, the proposed design procedure should be applicable to various near-surface mounted reinforcement configurations and materials.
Keywords
Introduction
Structural reinforced concrete (RC) members can require strengthening to compensate for deficiencies in either flexural or shear strength. The reasons for structural strengthening are varied; most structures require strengthening to compensate for strength loss due to ageing, while some structures were damaged in some ways that result in loss of strength. The focus of this article is the near-surface mounted (NSM) method (Badawi and Soudki, 2009; Bilotta et al., 2011; Capozucca and Bossoletti, 2014; Galati and De Lorenzis, 2009; Sharaky et al., 2017; Wu et al., 2014), which is a type of flexural strengthening for RC beams or slabs. The NSM method involves drilling grooves within the concrete cover and placing an NSM reinforcement within it, after which the groove is filled with epoxy adhesive. Currently, there are also several new derivative methods based on the NSM method, such as the side-NSM (SNSM) method (Sharaky et al., 2017; Shukri et al., 2016a) where the NSM reinforcement is placed at the sides of the beam to allow NSM strengthening to be applied on beams with small width or with beam soffits which are inaccessible and the partially bonded NSM (Choi et al., 2011; Seo et al., 2016) where the high moment area of the NSM is left unbonded to increase ductility.
The main problem of the NSM method is its vulnerability to concrete cover separation (CCS) failure. The CCS involves a debonding crack which appears at the location of NSM reinforcement curtailment, which then propagates towards the higher moment area of the beam. It is a type of premature mode of failure, which means that NSM-strengthened beams that failed by CCS will have failed well below the design strength. There has not been much research done on CCS failure on NSM-strengthened beams, likely due to the presence of a large number of parameters involved. Several methods to predict or simulate CCS have been proposed using the finite element method (Al-Mahmoud et al., 2010; Zhang and Teng, 2014) or using the concrete tooth model (De Lorenzis and Nanni, 2003). Recently, Teng et al. (2016) proposed a strength model for NSM carbon fibre–reinforced polymer (CFRP) strips derived using finite element study while an analytical design approach was proposed by Rezazadeh et al. (2016), which was derived using concrete fracture mechanic. Most of these methods can be highly empirical, such as in terms of predicting crack spacing. Empirical methods that are formulated around a specific shape or material type of NSM reinforcement are only accurate within the regime of testing used to formulate them, which can limit their usage.
In response to all these problems, this article proposes a mechanics-based approach to design, which can prevent CCS failure while being less reliant on empirical means. An example of this can be seen in the work of Shukri and Jumaat (2016), where the moment-rotation (M/θ) approach (Knight et al., 2014; Mo et al., 2016; Oehlers et al., 2012; Shukri et al., 2015, 2016b; Visintin et al., 2013; Visintin and Oehlers, 2016, 2017) was used in conjunction with the global energy balance approach (GEBA) (Achintha and Burgoyne, 2013, 2011; Guan and Burgoyne, 2014) to simulate CCS failure. The M/θ approach is a mechanics-based method that applies the partial interaction theory (Gupta and Maestrini, 1990; Haskett et al., 2008; Muhamad et al., 2012; Visintin et al., 2013) to directly simulate concrete cracking, crack widening, and tension stiffening of RC beams. The GEBA, on the other hand, applies the fracture mechanics of concrete to predict whether the CCS debonding crack will propagate and cause failure. The proposed method by Shukri and Jumaat (2016) shows that the use of these two mechanics-based methods allows the CCS failure of NSM-strengthened beams to be predicted with good accuracy and applicable to various types of NSM reinforcement material and shape due to lower reliance on empirical means. However, it was not made for general design purpose and as such is difficult to use. Hence, in this article, closed-form solutions for crack spacing and load–slip relationships will be used to formulate a simpler design procedure based on the M/θ-GEBA method, which is then validated against published experimental results.
Fundamental theories
In this section, the fundamental theories used in the M/θ approach and GEBA will be presented. The aim of this article is to introduce a design procedure, hence only a brief description of the M/θ approach and GEBA that is relevant to this article are given while references to the original research are provided.
M/θ approach
The M/θ approach used is a segmental simulation focused on the behaviour of cracked RC beam, as shown in Figure 1(a), which feature an NSM-strengthened RC beam section with cross-section as shown in Figure 1(b). The spacing between these primary cracks is designated Scr, and due to symmetry of forces within the RC beam segment, it is possible to only consider half the crack spacing, Ldef, for analysis purpose.

RC beam: (a) NSM-strengthened beam segment and (b) cross-section of NSM-strengthened beam.
Slip occurs along Ldef between the steel and NSM reinforcements and adjacent concrete, causing a cumulative slip of δr and δf at the crack face, respectively; this occurrence can be idealized as prism with a reinforcement and its adjacent concrete as shown in Figure 2, where slip of steel reinforcement (δr) is shown as example. The size of the prisms is as shown in Figure 1(b), where cr refers to the distance from the centroid of the steel reinforcement to beam soffit and cf refers to the distance from the centroid of the NSM reinforcement to the beam soffit.

Partial interaction tension stiffening analysis of steel reinforcement.
The bond stress, τ, gradually transfers load Pr onto the adjacent concrete, causing both Pr and δr to reduce the further it gets from the crack face. Using the boundary condition of δr = 0 at length Ldef, the stresses and strains within the prism can be solved to determine the force Pr which causes the slip δr, which allows the load–slip (P-δ) relationship to be obtained as shown in Figure 3. This procedure is called the partial interaction tension stiffening analysis (Haskett et al., 2008; Muhamad et al., 2011, 2012); the benefit of this analysis is the resulting load–slip relationship has directly accounted for the effects of tension stiffening through the use of bond stress–slip relationship, such that other empirical indirect methods, such as the commonly used Branson’s equation (Branson, 1968), are not needed.

Load–slip relationship of steel or NSM reinforcement.
From Figure 3, it can be seen that there is an initial linear portion OA with stiffness K (Zhang et al., 2017), such that the Pr-δr and Pf-δr relationship can be written as Pr = Krδr and Pf = Kfδf, respectively. This assumption should be correct for steel reinforcements prior to steel yielding and for NSM FRP reinforcements, which do not yield. Closed-form solutions based on this theory have been proposed (Visintin and Oehlers, 2016; Zhang et al., 2017) and will be used extensively in this article.
GEBA
The GEBA is a fracture mechanics-based method where the primary assumption is that CCS debonding cracks will always appear on strengthened RC beams, causing a debonded length Ld as shown in Figure 4. It is then only a matter of determining whether there is enough energy for the debonding crack to propagate and cause failure. When the GEBA is applied on beams strengthened with FRP sheet (Achintha and Burgoyne, 2011) or NSM (Shukri et al., 2018; Shukri and Jumaat, 2016), it is usually assumed that the debonding crack starts to propagate at an angle of 45° to the beam axis until it reaches the shear link; with this assumption, Ld is considered equal to the concrete cover’s depth as shown in Figure 4(a). It was stated by Achintha and Burgoyne (2013) that the actual direction of the crack may be slightly varied from 45°, but it should not have a significant effect on the results. On the other hand, for SNSM-strengthened beams, the experimental studies on SNSM-strengthened beams (Hosen et al., 2015; Sharaky et al., 2017; Shukri et al., 2016a) had shown that the shear crack starts to propagate horizontally as it reaches the SNSM reinforcement, such that Ld = 0 as shown in Figure 4(b).

Initial debonded length of NSM- and SNSM-strengthened beams: (a) NSM strengthened beam; (b) SNSM strengthened beam.
As the GEBA is concerned with the start of the fracture process, mode II effects of concrete fracture such as aggregate interlock are not relevant; hence, the focus of the GEBA is mainly on mode I (Achintha and Burgoyne, 2008). For the fundamental procedure of the GEBA, consider the NSM-strengthened beam section as shown in Figure 4, with a CCS debonding crack already present at location La where moment Ma is acting. It is assumed that the debonding crack will propagate instantaneously, such that Ma remains the same as the change from strengthened to unstrengthened section due to CCS debonding occurs, as shown in Figure 5. The strain energy available at La that can cause the debonding crack to propagate, designated as Wa, can be determined from the difference between the moment–curvature (M/χ) relationship of the NSM-strengthened section (M/χ) s and the unstrengthened section (M/χ) u . The method to determine fracture energy available for debonding (Ga) is presented in section ‘CCS prediction’.

Moment–curvature relationship of strengthened and unstrengthened beam section.
Design procedure for CCS
The proposed mechanics-based design procedure will be presented in this section, while a flowchart is given in Figure 6 as an overview. A preliminary design for the NSM strengthening is first made. Assuming the maximum moment that the beam needs to withstand is known, the maximum moment at La, Ma is known. The rotation at the strengthened and unstrengthened section due to moment Ma is then determined and the GEBA is used to determine if the debonding crack will propagate. If the beam is predicted to fail by CCS, the design for the NSM strengthening is changed and the procedure is repeated. Once a suitable design for NSM strengthening is determined, the flexural strength of the NSM-strengthened beam is then determined.

Design procedure for NSM-strengthened beams.
The proposed design procedure is applicable for both virgin and cracked RC beams. It has been shown that the crack spacing near the support of the beams is usually similar regardless of whether the beam is virgin or cracked prior to NSM strengthening (Shukri et al., 2018), which allows the method for CCS prediction presented in Figure 6 to be used as it is for both types of beams. Furthermore, it was shown that the difference in ultimate load between virgin and cracked NSM-strengthened RC beam to be negligible when the failure is by flexure rather than CCS (Shukri et al., 2016a, 2018).
Primary crack spacing and load–slip stiffness
The primary crack spacing, Scr, and the length of deformation, Ldef, can be defined through the mechanically derived equation by Sturm et al. (2018)
where
where fct is the tensile strength of concrete, Lper is the total perimeter of a single tensile steel reinforcement of area Ar contained within the tension stiffening prism as shown in Figure 2. Acr is the area of adjacent concrete in the tension stiffening prism, which can be determined using Figure 1(b). The variables α, δ1 and τmax are the ascending branch of the non-linear bond stress–slip relationship by CEB-FIP (1993) where α = 0.4, δ1 = 1 mm and τmax = 1.25√fc.
Having determined the Scr, the load–slip stiffness parameter for the steel reinforcement, Kr can be determined as
where Er is the elastic modulus of steel and c2
The coefficient c2 allows for the effect of bond and was determined using semimechanical means by Zhang et al. (2017) using a numerical tension stiffening analysis shown in Figure 2 to perform parametric study, from which c2 was extracted. Since it was derived for steel reinforcement, it is unsuitable for NSM reinforcements. Hence, a parametric study was performed in this research to extract c2 for NSM CFRP bar and NSM CFRP strip.
Parametric study for c2
The objective of this brief parametric study is to obtain the relationship between Af/Acf and c2; Af is the area of NSM reinforcement and Acf is the area of adjacent concrete for the tension stiffening prism containing the NSM reinforcement. Three types of NSM reinforcement configuration will be used:
NSM CFRP bar with ratio of groove to bar diameter of 2;
NSM CFRP bar with ratio of groove to bar diameter of 1.5;
NSM CFRP strip with ratio of groove height to width of 2.75.
The NSM reinforcement configurations above were chosen as they are commonly used in the literature. The bond stress–slip relationship by De Lorenzis (2004) and Zhang et al. (2013) were used for NSM CFRP bar and NSM CFRP strip, respectively. The material properties used in the parametric study were fixed and are shown in Table 1. The Scr of a NSM-strengthened beam is assumed to be controlled by the steel reinforcement and not the NSM reinforcement (Shukri et al., 2015; Shukri and Jumaat, 2016); hence, for this parametric study of NSM reinforcements, the Scr is fixed to the value given in Table 1.
Fixed properties for parametric study.
NSM: near-surface mounted.
The full numerical procedure for the partial interaction tension stiffening analysis, as illustrated in Figure 2, has been shown in multiple published research papers (Haskett et al., 2008; Muhamad et al., 2011; Shukri et al., 2015, 2016b; Shukri and Jumaat, 2016) and so will not be repeated again here. The results of the parametric study are shown in Figures 7 to 9, where c2b2 is the coefficient of bond for NSM CFRP bar with ratio of groove to bar diameter of 2, c2b1.5 is the coefficient of bond for NSM CFRP bar with ratio of groove to bar diameter of 1.5 and c2s is the coefficient for NSM CFRP strip. From the results, it can be seen that the bond effect coefficient c2 decreases when Af/Acf increases.

Parametric study for NSM CFRP bars with ratio of groove to bar diameter of 2.

Parametric study for NSM CFRP bars with ratio of groove to bar diameter of 1.5.

Parametric study for NSM CFRP strips with ratio of groove height to width of 2.75.
Extracting the coefficients from the parametric study results using a linear relationship yield
The stiffness parameter for NSM CFRP bar can then be determined as
where Ef, Af and Acf are the elastic modulus of the NSM reinforcement, area of a single NSM reinforcement and area of adjacent concrete area within the tension stiffening prism of the NSM reinforcement, which can be determined using Figure 1(b). The coefficient c2b2 in equation (10) can be changed to either c2b1.5 or c2bs according to the NSM reinforcement configuration used. If some other NSM material or NSM configuration is used, to obtain the relevant coefficient c2 is only a matter of performing a numerical parametric study similar to what is shown here.
Depth of neutral axis for unstrengthened RC beam section
Having determined the Kr and Kf, the load acting on a single steel and NSM reinforcement for a given slip can be determined as
Let nr be the number of total steel reinforcement used in the RC beam. From equation (11), the sum of forces acting on all the steel reinforcements can be simplified as
where
Now consider the equation of equilibrium of the unstrengthened section
where Pc is the compressive force of concrete. Assuming a triangular-shaped concrete stress and inserting equation (13) into equation (15) leads to
where Ec is the elastic modulus of concrete, δT is the deformation of the topmost section of the beam and dna-u is the depth of neutral axis for the unstrengthened section. From Figure 10(a), the relationship between the rotation and the slips and deformations of a strengthened beam is

Moment-rotation of NSM-strengthened beam section: (a) beam section of length Ldef and deformation profile, (b) strain profile, (c) stress profile and (d) forces acting on the beam section.
where h is the height of the beam and dna is the depth of neutral axis. From equation (17), the relationship between the slip of reinforcement and concrete deformation at the topmost section of the beam for the unstrengthened section can be written as
Inserting equation (18) into equation (16) leads to
To simplify equation (19), let
Replacing equation (20) into equation (19) yields
Solving equation (21) for dna-u gives
Rotation of unstrengthened RC beam section
Since a triangular-shaped concrete stress was assumed, the lever arm for the steel reinforcement in the unstrengthened section, Zru, is
The equation of moment for the unstrengthened beam section is
Inserting equation (13) into equation (24) and rearranging leads to
where δru is the slip of steel reinforcement for the unstrengthened beam section. From equation (17), the rotation of the unstrengthened section can be determined as
Depth of neutral axis for strengthened beam section
From equation (17), the following relationships can be obtained
where dna-s is the depth of neutral axis for the strengthened section of the RC beam. The equilibrium equation for the strengthened beam section is
where Pf is the force acting on the NSM reinforcement and nf is the number of NSM reinforcement. Similar to equation (13), total NSM FRP force Pfnf can be written as
where
Assuming a triangular shape for the concrete stress and inserting equations (13), (27), (28) and (30) into equation (29) gives
Replacing equation (20) in equation (32) and simplifying it gives
Solving equation (33) for dna-s gives
Rotation of strengthened beam section
Since a triangular-shaped concrete stress was assumed, the lever arms for the steel and NSM reinforcements in the strengthened section can be written as
where Zrs and Zfs are the lever arm for the steel and NSM reinforcement in the strengthened beam section, respectively. The equation for moment in the strengthened beam section is
From equation (17), the following relationship can be obtained
where δf and δrs are the slips of NSM and steel reinforcement in the strengthened section, respectively. Replacing equations (13), (30) and (38) into equation (37) and solving for δrs leads to
From equation (17), the rotation of the strengthened section can be determined as
CCS prediction
The fracture strength of concrete can be determined using any appropriate model. Here, the CEB-FIP model (CEB-FIP, 1993) will be used
where Sa is maximum aggregate size used in the concrete in millimetres and Gfo is the base value for concrete fracture strength and can be taken as 0.037 for aggregate of size 20 mm. The fracture energy of the NSM-strengthened beam due to moment Ma can be determined as
where Δ L is the propagation of CCS debonding crack, bf is the total width of the tension stiffening prism for NSM reinforcements, which as shown in Figure 1(b) can be taken as the width of the beam (b) if two or more NSM reinforcement is used, or half the width of beam if only one NSM reinforcement is used (Shukri et al., 2015) and
If Ga > Gmax, the beam is predicted to fail by CCS. Assuming this occurs, to redesign the NSM-strengthened beam is only a matter of changing the properties and configuration of the NSM strengthening, such as diameter of the NSM reinforcement (which results in a smaller Af) and number of NSM reinforcement (nf). The design procedure is then repeated starting from section ‘Depth of neutral axis for strengthened beam section’.
With regard to the value of Δ L , as the GEBA is concerned with the start of the fracture process, the value Δ L used is usually a very small value. Previous researchers had used Δ L = 1 mm in their research (Achintha and Burgoyne, 2008; Shukri et al., 2018; Shukri and Jumaat, 2016). It was also noted that values of Δ L <1 mm can cause numerical convergence problem (Achintha and Burgoyne, 2008). In this article, Δ L = 1 mm will be used based on a sensitivity analysis of Δ L that will be presented in section ‘Sensitivity analysis for Δ L ’
Existing slip due to dead load
The dead load on RC beams prior to strengthening can affect the serviceability condition of the beam. Assuming that the moment due to existing permanent action (dead load), Mp, can be determined by the design engineer, equations derived for the unstrengthened beam section can be used to determine the slip of steel reinforcement due to permanent action
The lever arm should remain the same, as the equation used to obtain it and to obtain the depth of neutral axis is independent of applied load. Note that the equation above assumes a triangular shape for concrete stress. The existing crack width prior to strengthening would be equal to 2δp.
Design flexural strength of NSM-strengthened beam
The design flexural strength of NSM-strengthened beam will be limited by the concrete crushing. For simplicity, several assumptions are used:
A rectangular concrete stress block is assumed.
The concrete strain will be limited to 0.0035 to adhere to Eurocode 2 requirement.
The strain hardening of the steel reinforcements will be ignored.
Compression bars will be ignored.
The assumptions used are common in design practice. Since the concrete strain will be limited to 0.0035, the following relationship is obtained
The equilibrium equation at ultimate limit state is
where Py is the load at which the steel reinforcement yields. From equation (17)
where dna-d is the depth of neutral axis for the design strength of the NSM-strengthened beam and δfd is the slip of the NSM reinforcement at ultimate limit state. Inserting equation (47) into equation (49) yields
Replacing equations (30) and (50) into equation (48) and expanding yields
where concrete stress is assumed to be equal to fc at concrete strain of 0.0035, while αc and βc are ratio of equivalent concrete stress and equivalent concrete stress block height, respectively. From Eurocode 2, αc = 0.85 and βc = 0.8 for fc≤50 MPa. Equation (51) can be simplified into
where
Solving equation (52) for dna-d
With dna-d known, the slip of the NSM FRP bar, δf, can be determined using equation (49). The design flexural moment of the NSM-strengthened beam can be determined as
where Zrd and Zfd are the lever arm for steel and NSM reinforcements at ultimate limit state, respectively, and are determined as
It is also possible to determine the crack width of the beam at ultimate limit state. From equation (17)
Inserting equation (47) into equation (60) yields
where δrd is the slip of steel reinforcement due to applied load. The crack width, Δ r , of the beam is then determined as twice the slip due to permanent action and applied load
Validation
The proposed design procedure was validated using published experimental results (Al-Mahmoud et al., 2010; Barros et al., 2007; Ceroni, 2010; Sharaky et al., 2015; Shukri et al., 2015, 2016a; Teng et al., 2006) that failed by either CCS or flexure. In the validation process, the experimental and predicted failure mode will first be compared by determining the available fracture energy (Ga) and the fracture strength (Gmax). Comparisons between predicted and experimental moment of resistance at failure will also be given for beams that failed by flexure.
The geometric and material properties of the beam that are necessary to apply the design procedure are given in Tables 2 and 3. The beams were strengthened using NSM CFRP bars (Al-Mahmoud et al., 2010; Ceroni, 2010; Sharaky et al., 2015; Shukri et al., 2015), NSM CFRP strips (Barros et al., 2007; Teng et al., 2006) or with side-NSM CFRP bars (Shukri et al., 2016a). As the use of SNSM method only changes the lever arm of the NSM reinforcement, the proposed design procedure can be used without any changes. The maximum size of aggregates in concrete, Sa, was assumed to be 20 mm if not specified in the original research paper. Referring to Table 2, the moment Ma is the moment at the length La corresponding to the reported experimental failure load, Pmax, of each beam.
Geometric and loading properties of beams.
h: beam height; b: beam width; nr: number of tensile reinforcement bars; nf: number of FRP bars/strips; cr: distance from beam soffit to centre of tensile reinforcement bars; cf: distance from beam soffit to centre of FRP bars/strips; La: distance from end of beam to the location of curtailment of FRP bars/strips.
Material properties.
fc: concrete strength; σy: yield strength of steel; Er: elastic modulus of steel; Ef: elastic modulus of FRP bar/strip.
The result of the validation is given in Table 4, where the calculated Gmax, Ga, failure mode, experimental maximum moment (Me) and maximum moment obtained using the design procedure (Md) are presented. It should be noted that Md is only given when it is predicted that the beam fails by flexure. From Table 4, it can be seen that the proposed design procedure was able to correctly predict the CCS failure mode of the beams apart from beam A2, SNC10 and PSNC10. In the case of beam A2, the predicted failure mode is flexure, whereas the actual failure mode is CCS. The ratio of Md/Me, however, is very close, at 0.99, which shows that while the beam was reported to fail by CCS, the loss of strength due to the premature debonding is negligible. For the case of beams SNC10 and PSNC10, the predicted failure mode was CCS, while the actual failure mode was flexure. The ratio of Ga/Gmax for beam SNC10 shows that the fracture strength was only very slightly exceeded, hence the incorrect predicted failure mode can be due to difference in the calculated and actual fracture energy, Gmax, which is not uncommon due to the variable nature of concrete.
Predicted and experimental result comparison.
Ga: predicted fracture energy; Gmax: fracture strength; FMe: experimental failure mode; FMd: predicted failure mode; Me: experimental maximum moment at failure; Md: predicted maximum moment at failure; CCS: concrete cover separation failure mode; F: flexural failure mode α; n/a: not applicable.
Sensitivity analysis for Δ L
A sensitivity analysis was conducted using the properties of nearly all the beams in the used in the validation process that failed by CCS with the result as shown in Figure 11. As all the beams had experimentally failed by CCS, the Ga/Gmax should be more than 1. Where Ga/Gmax is less than 1, it shows that the value of Δ L used failed to provide a correct prediction of CCS failure for the beam.

Sensitivity analysis for Δ L .
It can be seen from Figure 11 that Δ L = 1 mm was able to give an accurate assessment of CCS failure for all beams. Values of Δ L <1 mm were also able to correctly predict the CCS failure; however, it should be noted that when Δ L <1 mm, the resulting Ga can be considerably higher. This can cause a very conservative design, as the design engineer may have to greatly reduce the amount of NSM reinforcement provided in order to ensure that Ga < Gmax and prevent CCS failure. Hence, the value of Δ L = 1 mm is suggested as it gives a reasonable balance between accuracy and conservativeness.
Conclusion
In this research article, a mechanics-based design procedure was proposed. The proposed design procedure uses the M/θ approach and the GEBA to predict the behaviour of NSM-strengthened RC beams and the CCS failure mode. Several conclusions can be made based on this study:
Published experimental results of beams strengthened with NSM CFRP bars, NSM CFRP strips or SNSM CFRP bars were used to validate the proposed design procedure and good correlation was found between the experimental and predicted results.
The proposed design approach should be more versatile compared to other existing design approach as it is less reliant on empirical formulations. Hence, it can easily be applicable to most types of NSM reinforcement material and configurations. The design approach can also accommodate any new innovations in terms of the NSM strengthening material or configurations. The coefficient of bond c2 for those new NSM material or configuration can be determined using a similar approach described in this article.
Footnotes
Appendix 1
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship and/or publication of this article: Financial support from the Fundamental Research Grant Scheme (FRGS), Ministry of Education Malaysia (project number FP004-2014B) is gratefully acknowledged.
