Experimental and numerical studies on the flexural behavior of fibre reinforced concrete beams with innovative hybrid FRP-wrapped steel bars

ABSTRACT The durability and lifespan of concrete structures are significantly affected by the corrosion of steel reinforcement. Though epoxy-coated steel reinforcement is being used for practical applications, slight damage in the coating during the construction process is enough to initiate the corrosion process. The shortcomings in the current practices to prevent corrosion of steel reinforcement have led the researchers to look for alternatives. Innovative hybrid reinforcing elements with solid steel core and fibre-reinforced polymer (FRP) wrapping (called steel–FRP composite bars (SFCB)) can be considered as a suitable alternative to conventional steel reinforcement. In this research paper, experimental and numerical studies that investigated the flexural behaviour of concrete beam specimens reinforced with SFCB produced using glass fibres or carbon fibres are presented and discussed. Experimental study consisted of flexural testing of six concrete beam specimens with conventional steel reinforcement or SFCB with glass fibres or SFCB with carbon fibres. Two different spacings of shear stirrups were considered. Failure mode, first crack load, ultimate flexural load and load-deflection behaviour were recorded during the experiment. The test results showed that beam specimens with SFCB reinforced achieved less ultimate flexural load as compared to beam specimens with conventional steel reinforcement. The failure mode was found to be influenced by the spacing of the shear stirrups. The results of the numerical analysis were close to the experimental results. .


Introduction
Reinforced concrete (RC) is inevitable in the construction of civil infrastructure facilities such as bridges and buildings. The durability of RC structures is significantly affected due to corrosion of steel reinforcement, and this reduces their life span. It has been reported that the United States spends around 276 USD billion (3.1% of GDP) annually as a direct cost to corrosion (URL-1). In order to curb the cost incurred due to corrosion in the future, new innovative reinforcing elements such as fibrereinforced polymer (FRP), plasma rebars, corrosionresistant rebars (steel rebars with higher content of chromium) and epoxy-coated rebars are being adopted for practical applications.
research studies (Toutanji and Saafi (2000), Toutanji and Young (2003), Masmoudi, Theriault, and Benmokrane (2001)) have been conducted that investigated the influence of using FRP bars as reinforcement on the flexural behaviour of concrete beams. However, the use of FRP rebars for practical applications is limited primarily due to lower elastic modulus, brittle nature and high cost. In order to reduce the cost and achieve durable and ductile RC structures, researchers have proposed using combined steel and FRP bars as reinforcement (called hybrid reinforcement) in concrete members (Xu et al. (2019), Sun et al. (2019)). This type of reinforcement would reduce the cost of construction as compared to that incurred by using only FRP bars. While using hybrid reinforcement, FRP bars are placed near the outer surface of the concrete members (with lesser cover thickness) and steel bars are placed at inner areas (with higher cover thickness) thereby prolonging the time of initiation of corrosion and enhancing the lifespan of the structures. It is to be noted that, even though adopting hybrid reinforcement would increase the concrete cover for tension steel, steel stirrups that would have lesser cover would be used. Steel stirrups would therefore provide an electrical contact of the inner steel bars, and there are possibilities that the envisaged lifespan of the RC members might not be achieved. Towards improving the possibilities of combined use of steel bars and FRP bars for reinforcing concrete members, attempts have been made as early as the 1980s for developing rebars consisting of steel core and FRP wrapping around the core (herein called as steel-FRP composite bars (SFCB)). A typical photograph of SFCB is shown in Figure 1. The FRP shielding would prevent or significantly reduce corrosion of steel core, thereby improving the lifespan of RC structures reinforced with SFCB and subjected to corrosive environment, especially in coastal regions. Zhou et al. (2019) have conducted experiments for studying the corrosion behaviour of SFCB using X-ray microcomputed tomography. Test results have showed that the type of fibre, micropore structure of the fibre coating layer and manufacturing process arethe factors that influenced the corrosion resistance of the rebars. The corrosion rate of SFCB was found to be significantly less than that of bare steel rebars. The authors have also found through tensile tests that the elastic modulus of the SFCB was more than that of GFRP and the tensile strength was more than that of steel rebars. Zhao et al. (2020) have conducted experimental studies to determine the tensile properties of SFCB and their bonding performance in concrete. Test results have showed that SFCB with round steel core and FRP wrapping showed better tensile behaviour as compared to those with ribbed steel core and FRP wrapping. The authors have observed that the rule of mixtures can be used to describe the tensile behaviour of SFCB. The stress-strain curve of SFCB showed a bilinear trend. Bond strength of SFCB was influenced by the diameter and the surface treatment.
Fibre-reinforced concrete consisting of steel or synthetic fibres is reported to enhance the ductility of concrete, and it has also been shown that the addition of fibres improves the resistance of concrete to formation and propagation of cracks (Zhao et al. (2016), Yehia et al. (2016), Aldabagh, Abed, and Yehia (2016), Farrag and Yehia (2015), Yehia et al. (2015), Yehia (2012)).
In the literature, experimental research studies that determined the behaviour of concrete structural elements with SFCB are reported. However, experimental and numerical studies that investigated the flexural behaviour of fibre-reinforced concrete beams consisting of SFCB are not found to be reported. Numerical studies reported in the literature consisted of simulating the flexural behaviour of concrete beams with FRP bars as reinforcing elements or hybrid steel and FRP bars as reinforcing elements. Abed et al. (2020) have carried out numerical simulations on the behaviour of concrete elements reinforced with FRP bars under different loading conditions. The authors have observed that the predicted behaviour of the elements compared well with the experimental results. Hawileh (2014) conducted finite element studies to predict the load vs deflection behaviour of concrete beams reinforced with hybrid steel and FRP bars produced with aramid fibers. The authors observed that the load vs deflection behaviour could be predicted well using the model developed. The author also remarked that only limited numerical studies that investigated the behaviour of FRPreinforced concrete elements are available in the literature. Bencardino, Condello, and Ombres (2016) have carried out numerical studies on the behaviour of concrete beams with steel, FRP and hybrid FRP-steel reinforcements. The authors have observed that the tension stiffening model they have considered predicted the behaviour well for beams with low and normal reinforcement ratio; however, for beams with higher reinforcement ratio, the numerical predictions were not comparable to the test results. Similar numerical studies can be found elsewhere in the literature (Al-Rahmani and Abed (2013), D D, Pisano, andFuschi (2014), Ferreira et al. (2001), Rafi, Nadjai, and Ali (2008), Nour et al. (2007)). To the authors' knowledge, numerical simulation studies that predicted the flexural behaviour of concrete beams consisting of SFCB as reinforcing elements are not available in the literature.
The combined use of SFCB and steel fibres as reinforcement can be considered to be a novel technique for achieving higher durability and improved cracking resistance in concrete structural members. For developing design guidelines for practical applications, flexural behavior of concrete beams reinforced with SFCB and fibres has to be investigated experimentally and numerically, and that motivated the research presented in this paper.
In order to fill these gaps in the literature, in this paper, experimental studies are conducted to investigate the flexural behaviour of fibre-reinforced concrete beam specimens consisting of SFCB as tensile reinforcement. Two types of SFCB produced with glass fibers and carbon fibres are considered. The flexural behaviour of concrete beam specimens with SFCB is compared with control concrete beam specimens reinforced with conventional steel bars. The beam specimens are tested under 4-point bending with simply supported boundary conditions. Numerical analyses consisted of developing models of the tested beam specimens in general-purpose finite-element software ABAQUS and validating with test results.

Accelerated corrosion test
A pilot experimental study (Accelerated corrosion test) was conducted to determine the corrosion behaviour (gravimetric weight loss) of the proposed SFCB and compare it with that of conventional steel bars. It is known that corrosion of steel is an electrochemical process, and it involves anode, cathode and an electrolytic solution. The process of corrosion is shown in Figure 2. (URL-2). The corroding portion of steel is the anode and the un-corroded portion is the cathode. Electrolyte solution provides electrical contact between the anode and the cathode. Under favourable conditions, ferrous ions from anode (steel) enter into the electrolyte. The electrons that are released from the corroding part of the steel (anode) move towards the cathode and form hydroxyl ions by combining them with water and oxygen. These hydroxyl ions react with the ferrous ions that are released from the anodic portion of steel and form hydrated ferrous oxide. The ferrous oxide further gets oxidized to form ferric oxide which is known as the "red dust".
Three bare specimens such as (i) steel, (ii) SFCB with glass fibers and (iii) SFCB with carbon fibers were considered for accelerated corrosion test. The SFCB were manufactured manually. Details of the manufacturing procedure of SFCB are available in the literature (Sivaraman and Chinnaraju 2019).
The acceleration of the corrosion process was achieved by impressing an anodic potential between the rebar anode and steel plate (Saraswathy and Song (2008), Güneyisi, Özturan, and Gesolu (2006), Song et al. (2009)). The rebars were immersed in a 3.5% sodium chloride (NaCl) solution during the test. An impressed voltage of 10 V was applied to accelerate the corrosion process. The accelerated corrosion process was carried out for 24 h. The schematic diagram of the accelerated corrosion test setup is shown in Figure 3.
The SFCB specimens for the corrosion tests are shown in Figure 4.
The weight of steel and SFCB specimens was determined before the corrosion test and after the corrosion test and listed in Table 1. Figure 5 shows the steel and SFCB specimens after the accelerated corrosion test.
It was observed from the test results that the conventional bare steel rebar lost most of its material due to corrosion. However, SFCB rebars with glass fibers or carbon fibers were found to be less affected. It was observed that the percentage loss of conventional steel rebar was very high as compared to the percentage loss of steel in SFCB with glass or carbon fibers. The corrosion test results showed that the durability of concrete structures could be improved enormously by using SFCB as reinforcing elements. It was noted that the accelerated corrosion tests were conducted on bare steel bars and SFCB. However, in practical applications, the steel reinforcement would be embedded in concrete. Further research studies on the corrosion behaviour of SFCB embedded in concrete are required to determine their actual corrosion behaviour.

Details of beam specimens and flexure test setup
The experimental program consisted of testing six concrete beam specimens of dimension 130x230x1400mm (bxdxl) and reinforced with conventional steel rebars or    SFCB with glass fibers (GFRP wrapping) or SFCB with carbon fibers (CFRP wrapping) in the tensile zone. In SFCB, GFRP or CFRP wrapping was provided over the entire length of the steel core. The diameter of the conventional steel rebars was 10 mm, and the diameter of SFCB produced by wrapping GFRP or CFRP for a thickness of 1 mm over the steel core of 10 mm diameter was 12 mm. Fe500 grade steel rebars were used. The hanger bars consisted of 2#10 mm diameter steel rebars, and stirrups were made of 8 mm steel rebars spaced at 150 mm c/c or 300 mm c/c. The effective cover for the tensile reinforcement was kept as 30 mm. Fibrereinforced concrete with an average 28-day compressive strength of 26.3 was used for casting the specimens. The mix proportion of the concrete mix used was 1:1.51:2.77:0.45 in the order of cement (Portland Pozzolana Cement), fine aggregates (river sand) and coarse aggregates (20 mm gravel). Water-cement ratio of 0.5 was considered. Hooked type steel fibres with an aspect ratio of 60 were used in the concrete mix. The amount of fibres considered was 2% of cement content. The details of the beam specimens tested are given in Table 2.
In Table 2, the effective reinforcement ratios are determined using the equation given below, where ρ s = A s /bh o , ρ f = A f /bh o , A s is area of steel bars, A f is the area of FRP bars, b is section width, h o is section effective depth, f y is yield strength of steel, σ cg or σ cc is tensile strength of FRP using glass or carbon fiber. The corresponding balanced reinforcement ratio can be calculated as the ratio where concrete crushing and steel yield occur simultaneously given by ACI 440-1 R-06.
where f fu is tensile strength of FRP, ε cu is ultimate strain in concrete, ε fu is rupture strain of FRP bars, E f is modulus of elasticity of FRP bars, β1 is the strength reduction factor taken as 0.85 for concrete strength up to 27.6 Mpa. This factor can be reduced if the strength is in excess of 27.6 Mpa, but it should not be taken less than 0.65.
The beams were subjected to 4-point bending. The beams were simply supported at both edges with a clear span of 1200 mm. The loading was applied using a 300 kN hydraulic jack and the load was transferred using a rigid transfer beam to form two point loads on the beam specimens. The distance between one of the loading points and the nearest support was 400 mm, and the shear span-to-depth ratio was kept as 1.95. A schematic diagram of the beam specimen with reinforcement details is shown in Figure 6(a) and a typical photograph of the test setup used is shown in Figure 6(b).

Details of numerical model
The FE model was validated first with beam specimens with conventional steel reinforcement ( Figure  12a) and only then used for modelling beam specimens with SFCB. Numerical models were developed for all the tested beams consisting of shear stirrups at a spacing of 150 mm in the finite element software ABAQUS. Concrete and reinforcement (steel and SFCB) were modelled as 3D parts using C3D8Relements. SFCB was modelled with a steel core of diameter 10 mm and GFRP/CFRP wrapping over the steel core for a thickness of 1 mm. The bond between the steel core and the GFRP/CFRP wrapping was assumed to be perfect, and it was achieved by using Tie constraint. Test results indicated that the tested beam specimens with SFCB as tensile reinforcement achieved lesser ultimate flexural load as compared to the beams specimens with conventional steel reinforcement. The reason for this observation was attributed due to the lesser bonding action between the SFCB and the surrounding concrete due to the smoother surface of the SFCB reinforcement (see Section 4.3). This behaviour was simulated in the model by specifying interaction property at the SFCB-concrete interface. Guo et al. (2020) have observed that the friction coefficient at the steel-concrete interface could vary in the range 0.377-0.458. In this numerical model developed, the steel-concrete friction coefficient at the interface was kept as 0.3 for conventional steel reinforcement. For the beam specimens with SFCB, the friction coefficient of the SFCB-concrete interface was adjusted so as to get load- deflection curves comparable to that of test results. Iterations were made to get numerical load-deflection curves close to experimental load-deflection curves. Based on iteration, a friction coefficient of 0.1 was finalized as it resulted in numerical loaddeflection curves close to experimental load-deflection curves. Unidirectional tensile test specimen was prepared according to ASTM D3916. UTM was used for the tensile test. The loading rate was 5 mm/min. After the rods are cured, additional epoxy resin impregnated glass fabrics were wrapped around the ends to increase their end diameter at least 1 mm to serve as tabs for gripping in a tensile testing machine. The rods are approximately 10 mm in diameter and 600 mm in length. The hybrid rebars using glass fibers consist of 30 numbers of glass fiber strands and hybrid rebars using carbon fibers consist of 14 numbers of carbon fiber strands. Figure 7 shows the tensile specimen and unidirectional test set up. Figure 8 shows the stressstrain graph of the hybrid rebars.
The material property of steel was assumed to be bilinear. An earlier study by Zhao et al. (2020) indicated that the tensile strength of FRP based on the rule of mixtures was compared well with the experimental results.
In the present study, the yield strength of steel is 510 Mpa and the ultimate tensile strength of steel is 640 Mpa; tensile strength of GFRP and CFRP wrapping was obtained based on the rule of mixtures, and it was 1157MPa and 882MPa, respectively (Sivaraman and Chinnaraju 2019). The elastic modulus of GFRP and CFRP wrapping was also obtained based on the rule of mixtures and were found to be 37,920 MPa and 41,400 MPa, respectively (Sivaraman and Chinnaraju 2019).
The constitutive relationship for concrete was modelled using Concrete Damaged Plasticity option available in ABAQUS. The uni-axial stress-strain curve of  concrete in compression was modelled using the equation proposed by Saenz (1964) and is shown in Figure 9.
The behaviour of concrete under tension was modelled using the fracture energy criterion. The fracture energy of the concrete was calculated using Fib Model Code (2010) equation as given below.
Fracture lab energy The tensile strength of concrete was determined using the formula suggested by Graham and Scanlon (1986) as given below.
Tensile strength, where f ck is the characteristic compressive strength and f cm is the mean compressive strength of concrete.

Failure mode
The photograph of the tested beam specimens is shown in Figure 10. It is observed from Figure 8 that the cracking pattern of the tested beam specimens was influenced by the type of reinforcement used and the spacing of shear stirrups. More number of flexural cracks occurred in the beam specimens BS150 and BS300 (as compared to other beam specimens tested), and the final failure was due to the opening of one flexural crack formed near the mid-span. The number of cracks in the beam specimen BS150 was more than that formed in the beam specimen BS300. It was also observed that, at failure, in beam specimen BS150, crushing of concrete in the compression zone occurred in the mid-span. The beam specimens BS150 and BS300 failed in flexural mode.
In the beam specimens HBG150 and HBG300, the number of cracks formed at failure was significantly influenced by the spacing of the shear stirrups. The number of cracks in the beam specimen HBG150 was more than that in the beam specimen HBG300. The beam specimen HBG150 failed by the widening of the crack formed in the shear span region. The crack was initially oriented vertically such that it occurred  primarily due to flexural stress; however, it became inclined with an increased magnitude of applied flexural loading. The crack tended to join one of the loading points and the bottom edge of a cross-section that was located slightly away from the nearest support.
The beam specimen HBG300 failed due to an inclined shear crack that tended to join one of the loading points and the nearest support. At failure, the concrete in the compression zone of the beam specimen HBG300 also failed due to crushing. It is worth to note the difference in the observed failure modes of the beam specimens HBG150 and HBG300. The difference was that the inclined shear crack tended to join one of the loading points and bottom edge of the cross-section that was located slightly away from the support in the beam specimen HBG150 (called flexure- shear interaction failure mode); however, the inclined shear crack tended to join one of the loading points and the support in the beam specimen HBG300 (called as shear failure mode). The failure mode of the beam specimen HBG300 was exactly similar to the failure of deep beams. These observations clearly indicated that the spacing of shear stirrups influenced the failure mode of the beams reinforced using SFCB produced using glass fibre wrapping.
Except for the crushing failure of concrete in the compression zone in HBG300, the failure modes of the beam specimens HBC150 and HBC300 were similar to that of the beam specimens HBG150 and HBG300, respectively.

First crack load
The first crack load of the tested beam specimens is given in Table 3.
It is observed from Table 3 that the first crack load of all the tested beam specimens was influenced by the spacing of the shear connectors. For the beam specimens reinforced with conventional steel reinforcement, the first crack load increased with a decrease in the spacing of the shear connectors. However, for the beam specimens reinforced with SFCB, test results showed that the first crack load increased with an increased spacing of the shear connectors.
The first crack load of the beam specimen BS150 was 18% more than that of the beam specimen BS300. The first crack load of the beam specimens HBG300 and HBC300 were 100% and 35%, respectively, more than that of the beam specimens HBG150 and HBC150. This observation indicated that the first crack load was significantly influenced by the spacing of the shear stirrups in the beam specimens reinforced with SFCB produced using glass fibres as compared to those produced using carbon fibres.

Ultimate flexural load
The ultimate flexural load of the tested beam specimens are given in Table 3. It is observed from Table 3 that the ultimate flexural load of the tested beam specimens (BS150, HBG150 and HBC150) with 150 mm shear stirrups spacing was 5%, 3% and 7%, respectively, more than that of the corresponding beam specimens (BS300, HBG300 and HBC300) with 300 mm shear stirrups spacing. This observation indicated that the ultimate flexural load of the tested beams specimens was not significantly influenced by the spacing of the shear stirrups considered in this study.
The beam specimens with conventional steel reinforcement achieved higher ultimate flexural load as compared to the beam specimens with SFCB are tensile reinforcement. The ultimate flexural load of the beam specimens HBG150 and HBC150 was nearly 50% of the ultimate flexural load of the beam specimen BS150. Similarly, the ultimate flexural load of the beam specimens HBG300 and  Figure 11. Load-deflection curves of tested beam specimens. HBC300 was nearly 50% of the ultimate flexural load of the beam specimen BS300. The reason for less ultimate flexural load of the tested beams specimens with SFCB as tensile reinforcement could be attributed to the following. Visual observations indicated that the surface of SFCB reinforcement was relatively smooth as compared to conventional steel reinforcement. Due to the smoother surface of the SFCB reinforcement, its bond with concrete could have been less as compared to that of conventional steel reinforcement, and this could have resulted in lesser ultimate flexural load.

Load-deflection behavior
The load-deflection curves of the tested beam specimens are shown in Figure 11. It is observed from Figure 11 that the load-deflection curves of all the tested beam specimens showed tri-linear behaviour. The stiffness of the tested beam specimens was found to be influenced by the type of reinforcement and spacing of the shear stirrups.
It is observed from Figure 11 that the initial stiffness of the tested beams BS150&, BS300, and HBC150 & HBC300 can be assumed to be nearly the same irrespective of the spacing of the shear stirrups. However, Figure 11 shows that the initial stiffness of the tested beams HBG150 and HBG300 was found to be influenced by the spacing of the shear stirrups. The initial stiffness of the beam specimen HBG300 was more than that of the beam specimen HBG150.
With the increase in the applied flexural loading beyond the first crack load of the tested beam specimens, the flexural stiffness of the tested beams decreased linearly due to the formation and widening of cracks. Critical observation of Figure 11 indicated that even with increased loading magnitude (beyond first crack load), the stiffness of the beam specimen HBG150 was more than that of the beam specimen HBG300 till ultimate load. However, test results ( Figure 11) showed that, in general, the stiffness of the beam specimens HBG300 and HBC300 was more than the beam specimens HBG150 and HBC150, respectively. The reason for this observation could be attributed due to the difference in the failure modes of the beam specimens HBG300 and HBC300 as compared to that of the beam specimens HBG150 and HBC150. Since the failure mode of the beam specimens HBG300 and HBC300 was similar to that of deep beams, their behaviour was slightly stiffer (due to relatively less ductile behaviour) than the beam specimens HBG150 and HBC150, respectively.

Flexural toughness
The flexural toughness of the tested beam specimens was obtained as the area under the load-deflection curves and is given in Table 3. The variation in the flexural toughness of the tested beams is also shown in Figure 12.
It is observed from Table 3 and Figure 12 that, irrespective of spacing of shear stirrups, the flexural toughness of the beam specimens BS150 and BS300 with conventional steel reinforcement was more than that of other beam specimens consisting of SFCB as reinforcement. Among the beam specimens with SFCB as reinforcement, HBC150 achieved higher flexural toughness.

Numerical analysis results
Comparison of failure modes of beam specimens BS150, HBG150 and HBC150 obtained from the numerical simulation and experiment is shown in Figure 13.
It is observed from Figure 13 that the numerical models developed are able to predict the flexural failure of beam specimen BS150 and failure due to combined shear-flexural stresses of the beam specimens HBG150 and HBC150. Comparison of load-deflection curves of the beam specimens BS150, HBG150 and HBC150 obtained from numerical simulation and experiment is shown in Figure 14.
It is observed from Figure 14(a) that the loaddeflection curves of beam specimen BS150 obtained from numerical analysis and experiment are comparable. However, it is observed from Figure 14(b and c) that the initial stiffness of the beam specimens HBG150 and HBC150 as obtained from the numerical analyses was more than that obtained experimentally. The deviation was more for beam specimen HBG150 than that for the beam specimen HBC150. Comparison of experimental ultimate flexural load and as obtained from the numerical analyses of the tested beams BS150, HBG150 and HBC150 is shown in Table 4.
It was observed from Table 4 the predicted ultimate flexural load of beam specimen BS150 was 7% less than the experimental results. The predicted ultimate flexural load of the beam specimens HBG150 and HBC150 was 10% and 4%, respectively, more than the experimental results. In general, the results indicated that the predicted ultimate flexural load of the beam specimens was within ±10%.

Summary and conclusions
In this paper, experimental and numerical studies conducted on the flexural behaviour of fibre-reinforced concrete beams consisting of conventional steel bars or innovative hybrid bars in the tension zone were presented. Experimental study consisted of testing six concrete beam specimens with conventional steel bars, hybrid bars produced with glass fibres and hybrid bars produced with carbon fibres. Two different spacings of shear stirrups were considered. The conclusions based on the experimental and numerical studies presented in this paper are summarized below.
• Hybrid reinforcement consisting of FRP wrapping over steel core area can be considered as a suitable alternative to conventional steel reinforcement towards improving the lifespan of reinforced concrete structures located in a corrosive environment. • The flexural behaviour of concrete beams specimens was influenced by the type of reinforcement used in the tension zone. • The beam specimens with hybrid reinforcement failed either due to shear-flexural stresses interaction or due to shear depending upon the spacing of the shear stirrups. • The numerical models developed were able to predict the ultimate flexural load comparable to the concrete beam specimens BS150, HBG150 and HBC150 with reasonable accuracy. • Further research is required in this area towards developing design guidelines.

Disclosure statement
No potential conflict of interest was reported by the authors.

Notes on contributors
Sivaramakrishnan Subbaram, corresponding author, is working as a assistant professor in Sri Sairam Institute of Technology Chennai.
Chinnaraju Komarasamy, is working as Professor in College of Engineering Guindy. The first author is doing Ph.D under the guidance of second author.