Abstract
To address the various instability problems in cold-formed steel members, many researchers have mainly focused on developing innovative sectional profiles wherein geometry of the section plays a vital role in enhancing the inherent resistance of such sections against premature buckling. However, the process of forming such innovative shapes is not only complex and time-consuming but sometimes such sections fail to mobilize their complete reserve strength. Hence, a stiffening arrangement of weaker zones for mobilizing the untapped reserve strength is suggested. The contribution of this simple, effective and partly stiffening arrangements, aimed at eliminating/delaying the premature local buckling, is studied both experimentally and numerically and also compared with existing codes. Experimental study was carried out on different simply supported cold-formed steel beams with judiciously proposed stiffening arrangements under four-point loading. An equivalent hot-rolled steel beam was also tested to compare the efficiency of the cold-formed steel beams. The cold-formed steel beams investigated had different width-to-thickness ratio, different geometries and different stiffening arrangements. The test strengths, failure modes, deformed shapes, load versus mid-span displacements and geometric imperfections were measured and reported. The test strengths of the beam models are also compared with the design strength predicted by North American Standards and Eurocode for cold-formed steel structures. To validate the test results further, a numerical study was carried out on such stiffened cold-formed steel beams using finite element software ABAQUS. All these results show that the proposed strengthening system is efficient and economical and allow cold-formed steel beams to reach greater load carrying capacity.
Introduction
In most of the developing countries, structural designers face the global competition of designing fast-track and cost-effective structural systems to meet the huge deficit of various infrastructural systems. Such a development is crucial for bringing such countries at par with global standards of living. For this purpose, cold-formed steel (CFS) sections provide the best choice not only for cost-effective structural systems but also an ideal choice for the desired fast-track construction (Anbarasu and Sukumar, 2013, 2014). The main members in any structure need to have an adequate safety to avoid catastrophic or total collapse, thus use of hot-rolled steel sections for such vital members are unavoidable, owing to their superior performance against premature buckling, thus can be justified. However, use of hot-rolled steel sections for moderate to lightly loaded members such as floor beams and purlins generally remains under-utilized (Keerthan et al., 2014; Wang et al., 2014; Wang and Young, 2016). Moreover, optimization techniques in design are not fruitful because of the limited availability of rolled section in the market place (Subramanian and Venugopal (1977). This contributes to the highly conservative use of such a precious construction material with limited reserves. Keeping in view the importance of steel as an ideal construction material in challenging situations, it is of paramount importance to make the most optimum use of such vital construction material (Valsa Ipe et al., 2013). CFS provides an ideal choice to avoid such a wasteful use of steel. Unlike hot rolling, the cold forming process permits an almost infinite variety of shapes to be produced which can serve desired needs efficiently. Generally, the width-to-thickness ratio of individual components of CFS members is high, hence are prone to premature buckling at moderate compressive stress levels (much below the yield stress). This problem can be tackled effectively either by developing innovative sectional profiles with intrinsic resistance against premature buckling and/or by appropriate stiffening arrangement at vulnerable locations. Past research has stressed mainly on the development of innovative sections.
In the past decade, due to advances in the manufacturing technology, many attempts were made to make changes in the cross section of the member in order to attain economical and efficient sectional profiles (Hancock, 2016; Schafer, 2011). SudhirSastry et al. (2015) carried out a numerical study to investigate the effect of different flange configurations on the buckling of CFS channel beams. It was observed that the beams with extended open flanges and rounded flanges had enhanced critical buckling moments compared to beams with dropped flanges. Obst et al. (2016) carried out tests to study the behaviour of non-standard channel beams with single- and double-box flanges. Beams with double-box flanges had higher capacities than the ones with single-box flanges. It was also observed that reinforced beams performed better than the unreinforced ones. Ye et al. (2018) conducted a study to optimize the CFS channel beam sections for higher energy dissipation and ductility. It was found that by incorporating intermediate web and flange stiffeners to slender channel beams, there was a significant improvement in their energy dissipating capacity. However, in stockier channel beams, such introduction of stiffeners did not help. Instead, increasing the cross-sectional depth and decreasing the flange width helped in attaining higher energy dissipating capacity. Paczos and Wasilewicz (2009) tested anti-symmetrical CFS I-beam sections fabricated out of a single steel sheet. It was observed that there was quick loss in the stability of beams that were loaded with a concentrated force at the mid-span. Also, the critical loads for lipped beams were higher than that of the non-lipped ones. Trahair and Papangelis (2018) studied the lateral distortional buckling in hollow flanged beams with corrugated web plates. It was observed that the distortion in the hollow flanges resulted in the reduction of their torsional rigidities as well as their lateral buckling resistance. However, corrugation of web plate helped in preventing web distortion and significantly reduced the lateral buckling resistance, which is generally observed in beams with flat webs. Laím et al. (2015) conducted a series of tests to study the flexural behaviour of CFS beams with sigma profiles at ambient and elevated temperatures. It was seen that the web stiffeners behave differently under elevated temperatures and that behaviour depends upon the sectional profile. Axially unrestrained beams perform better than the restrained ones. Under ambient temperature, chances of excessive non-uniform compressive stress distribution are very likely compared to that in tension regions. Ye et al. (2016) carried out a numerical study to develop more efficient CFS channel sections under flexure. Just by optimizing the relative dimensions of flat plates and inclination of lips, the bending resistance can increase by nearly 30%. Double folded lips can substantially improve the flexural resistance, whereas intermediate web stiffeners may not necessarily improve flexural performance. Adequately designed CFS channel beam sections have the potential to reach their plastic moment capacity (Kumar and Sahoo, 2016). Siahaan et al. (2016a) studied innovative rectangular hollow flanged channel beams to develop optimum sections which can delay the buckling failure in them. Dar et al. (2015a) investigated various innovative CFS beam sectional profiles under flexure so as to find an ideal replacement for conventional hot-rolled steel sections. Wang and Young (2015) carried out an experimental and numerical investigation to study the local buckling and/or distortional buckling behaviour of the built-up open and closed sections under flexure. All these studies indicate that innovative profiles were successful in postponing the buckling failures in CFS beams; however, the process of forming such innovative sections of complex profiles requires a lot of time and effort, thus making the process very difficult and time-consuming. Appropriate stiffening of vulnerable zones is an effective alternate solution that can eliminate/delay premature buckling (Laı’m et al., 2013; Paczos, 2014; Moen et al., 2013). Hence, there is an urgent need to develop judicious stiffening arrangements that can be effectively used in simple CFS sections to overcome the complex problem of premature buckling failure, thus making CFS construction fast, simple and efficient.
The primary objective of this research is to conduct experimental investigation on different simply supported CFS beams with judiciously proposed stiffening arrangements under four-point loading. The tests would compensate the lack of experimental data on this form of construction and act as a data base for numerical models to be developed. An equivalent hot-rolled steel beam was also tested to compare the efficiency and the structural performance of these CFS beams. The CFS beams investigated had different width-to-thickness ratio, different geometries and different stiffening arrangements. The test strengths, failure modes, deformed shapes, load versus mid-span displacements and geometric imperfections were measured and reported in this article. The test strengths of the beam models are compared with the design strength predicted by North American Specifications and Eurocode for CFS structures. In order to validate the test results, a numerical study was also carried out on CFS beams using finite element software ABAQUS.
Experimental investigation
Test models
To achieve the well-defined objective of this study, four CFS beam models have been fabricated with and without appropriate stiffening arrangements. The tested CFS beams comprised two channel sections connected back to back by black bolts of size 5 mm and class 4.6. Two rows of bolts were provided with 100 mm centre to centre spacing along the depth as shown in Figure 1 and 200 mm centre to centre spacing along the length of the beam. To compare the efficiency and effectiveness of the proposed stiffening arrangements, one hot-rolled steel section was tested. The dimensional details of all the models are given in Table 1. A digital vernier calliper was used to measure the dimensions of various components of the models. The details of various models fabricated and tested are described below.

Arrangement of bolts in the cross section.
Model nominal and measured dimensions.
Model I: unstiffened model
The sectional geometry of this model consisting of I shape using two-lipped channel sections bolted back to back are shown in Figure 2(I). A cold-formed sheet having a thickness of 2 mm only has been used for the fabrication of this model. The nominal and measured dimensions of the various elements of the model are shown Table 1. All the beams had a span of 2.1 m.

Cross sectional details of the models.
Model II: angle-stiffened model
During the testing of Model I, localized lip buckling failure on compression side was observed at a low magnitude of loading. Therefore, it was expected that effective stiffening of compression area falling under high bending moment zone would have arrested such localized buckling, thereby considerably improving its structural performance. Accordingly, Model II, which is a modified version of Model I was stiffened by attaching two hot-rolled angle stiffeners (25 × 25 × 5) as shown in Figure 2(II). This hot-rolled angle stiffener was placed in the central 1.5 m length of the beam falling under high bending moment zone. The steel angle stiffener was welded to the inside of the compression flange lip by using Shielded Metal Arc Welding process. Rest of the details of Model II were kept strictly identical to Model I so as to investigate the contribution of the proposed variation in Model II towards improved structural performance. The dimensions of hot-rolled angle stiffener used in Model II are given in Table 2.
Dimensions of hot-rolled angle stiffener used in Model II.
Model III: unstiffened lightest model
Keeping in view the importance of achieving steel economy, it was deemed appropriate to fabricate lighter models using thinnest possible steel sheets. Accordingly, Model III with symmetrical I shape (using two-lipped channel sections back to back similar to Model I) involving simple fabrication was fabricated as shown in Figure 2(III). A cold-formed sheet having a thickness of 1 mm only was used for fabrication of this model.
Model IV: stiffened lightest model
During the testing of Model III, localized lip buckling failure (similar to Model I) on the compression side was observed in high bending moment zone. It was again expected that strengthening of flange lips in the compression zone could have prevented/delayed such a failure. Accordingly, Model IV which is a modified version of Model III was fabricated by stiffening the compression flange lips using a small-channel section as shown in Figure 2(IV). The small-channel section was bolted to the compression flange of the beam using 5 mm black bolts of class 4.6. For each small-channel stiffener, a single row of bolts were provided with centre to centre spacing of 150 mm along the length.
By developing structurally efficient CFS beam sections, steel economy was the main consideration. However, it was equally important to evaluate these proposed CFS sections from structural performance consideration too. Accordingly, Model V, that is, ISMB-150 (hot-rolled steel section), was chosen as the reference model for meaningful comparison with the various CFS beam models. The dimensions of ISMB-150 are given in Table 3.
Dimensions of hot-rolled beam ISMB 150.
Material properties
The test specimens were fabricated from locally available structural steel. Tensile coupon tests were used to determine the mechanical properties of the same. Since two categories of steel sheet thicknesses were used to fabricate the models. Three coupons were prepared from the centre of the flange in the longitudinal direction from each sheet. Various standards exist which specify the requirements for testing of tensile specimens. However, the dimensions of the coupons, as conforming to the Indian Standards (IS1608:2005), were used for material testing. A computerized universal testing machine was used to conduct the tensile tests of the coupons. The relevant material properties of the steel obtained from the material testing are given in Table 4. A typical stress–strain curve of CFS used in this study is shown in Figure 3. Since hot-rolled steel angle was used in Model II, three tensile test coupons were prepared from the flat portion of the angle along the longitudinal portion. The relevant properties of this hot-rolled steel angle are given in Table 5.
Material properties of steel used.

Stress–strain curve of CFS used in this study.
Material properties of hot-rolled steel angle used.
Geometric imperfections
Prior to testing, the initial overall geometric imperfections were measured. The imperfections were measured at the bottom flange web junctions near the centre both along the longitudinal as well as transverse directions. An optical theodolite and a calibrated vernier calliper were used to obtain the readings at the mid-length and near both ends of the models. The imperfections measured at the mid-lengths along the model in two orthogonal directions (Figure 4) are given in Table 6. The maximum geometric imperfection measured at the mid-length in δ1 and δ2 directions was 1/2167 mm and 1/2131 mm, respectively, and found in Model III. The minimum geometric imperfection measured at the mid-length in δ1 and δ2 directions was 1/4112 and 1/4346, respectively. As a comparison, the magnitude of the maximum and minimum imperfections measured by Yuan et al. (2015) was 1/100 and 1/4856, respectively.

Directions of measured geometric imperfections.
Measured geometric imperfections at mid-length.
Test setup
The testing of models was carried out on a loading frame of 500 kN capacity as shown in Figure 5. A hydraulic loading jack of 500 kN capacity was used to transmit load on a rigid spreader beam to ensure four-point loading as shown in Figure 6. Dial gauges of least count of 0.01 and 75 mm travel were used to record the vertical displacements. Identical bearing plates of size 150 mm ×150 mm × 15 mm were placed under concentrated loading points to prevent punching failure. To prevent web buckling under concentrated loading points, bearing stiffeners comprising two angles ISA 50 × 50 × 6 (SP 6-1 :2003) back to back were bolted to the web of the models on both sides. Simply supported end conditions with one end hinged and the other end pinned was adopted as shown in Figure 6. The moment span was laterally unrestrained. To assess the contribution of stiffening arrangements towards favourable structural performance, it was important to ensure strict uniformity of various parameters which include span of the beam, support conditions, location of applied point loads and bearing/stiffening arrangement under concentrated applied load/reaction points (Dar et al., 2015b, 2017). The supports were restrained laterally at the supports.

Model mounted on the testing rig.

Loading arrangement.
Test results and discussions
Sometimes, semi-log curves are plotted for better interpretation of results (Dar et al., 2018; Manikandan et al., 2014). Figure 7 shows the log of load versus displacement (at mid-span) curve of various models. The Model I resisted the load until it experienced premature local buckling failure (as seen in Figure 8) in the lip within the central one-third span with the maximum bending moment. The said failure was noticed corresponding to a maximum load of 44.1 kN and mid-span displacement of 11.67 mm.

Combined load–displacement curves of various models.

Lip buckling of Model I.
Except a localized lip failure confined to a small-length segment, rest of the model was in sound condition and still possessed enough reserve strength. To exploit this reserved strength, appropriate stiffening arrangement (as mentioned earlier) was introduced in Model II. The Model II resisted a higher load (68.3 kN) compared to Model I. The extra load carried by Model II was partly due to reserve capacity of the section and partly due to strain hardening. The mode of failure was initiated with local buckling of a compression flange shown in Figure 9 falling within the central one-third length and was under highest compressive stress corresponding to a failure load of 68.3 kN and a maximum deflection of 20.9 mm at mid-span. It is worth highlighting here that the stiffening arrangement adopted over a full length of the vulnerable zone has greatly contributed to a much-improved load carrying capacity from 44.1 to 68.3 kN (i.e. an increase of 39%). This encouraging experimental result confirms the important role of judiciously provided stiffening arrangements in CFS construction.

Compression flange buckling in Model II.
The Model III experienced premature local lip buckling failure within the middle third zone of the maximum bending moment as shown in Figure 10. The said failure was observed at a lower load of 12.7 kN and mid-span displacement of 5.0 mm. At this stage, the model suddenly stopped resisting any further load; hence, the model was said to have reached its limit state. Except localized lip failure confined to a small length, the rest of the model seemed in good condition and expected to possess appreciable reserve strength. To exploit this reserved strength, appropriate stiffening (as mentioned earlier) was introduced in Model IV.

Lip buckling failure in Model III.
The stiffened model (Model IV) failed at a load of 15 kN against 12.7 kN resisted by the unstiffened model (Model III; i.e. only a meagre increase of 2.3 kN). The mode of failure was again local lip buckling within the middle third zone of high bending moment and is prominently noticeable as seen in Figure 11. It is therefore concluded that using thin sheets with a thickness of 1 mm for the fabrication of CFS beams is highly vulnerable to premature buckling at very low loads even after adopting some stiffening arrangements; hence, such thin cold-formed sheets are not recommended. During testing of specimens, no distortional and lateral torsional buckling was visually observed in any of the selected cross sections.

Lip buckling failure in Model IV.
As mentioned earlier, it would be appropriate to have experimental results of a comparable hot-rolled section for meaningful comparison. Accordingly, a hot-rolled beam of ISMB-150 at 15 kg/m, named Model V, was chosen as the only reference beam model, with all other conditions remaining identical. The curve for ISMB-150 shows a linear response till it experienced lateral buckling as seen in Figure 12 corresponding to failure load of 83.65 kN and mid-span displacement of 16 mm. It needs to be highlighted here that the promising results of partly stiffened CFS Model II are comparable with reference hot-rolled section (i.e. Model V), thus confirming the vital role of proper stiffening arrangements in enhancing the structural performance of CFS construction. Hence, judiciously selected stiffening arrangements can be adopted with confidence in new CFS construction as well as in the strengthening of existing CFS structures (which demand upgradation).

Lateral buckling in ISMB-150.
The flat-width-thickness ratio of the flange had an effect on both strengths as well as stiffness as shown in Figure 13. The stiffness of beams was calculated by the method adopted by Deepak and Shanthi (2018). The effect of stiffening on the strength and initial stiffness of CFS beams is shown in Figure 14. As the flat-width-thickness ratio (b/t) of the flange reduced from 62.5 to 31.25 in unstiffened beams and from 31.25 to 18.75 in stiffened beams, the strength increased by 228% and 354%, respectively, as shown in Figure 13(a). For the same reduction in flat-width-thickness ratio (b/t) of the flange, the initial stiffness increased by 35% and 106% for unstiffened and stiffened beams, respectively, as shown in Figure 13(b).

Effect of flange’s flat-to-width ratio on ultimate load and stiffness: (a) ultimate load variation and (b) stiffness variation.

Effect of stiffening on ultimate load and stiffness: (a) ultimate load variation and (b) stiffness variation.
Design rules
Since one of the objectives of this research was to assess the strengths of the CFS beams against the current design standards. Accordingly, the un-factored design strengths predicted by European code and North American Standards were calculated and are summarized in Table 7.
Summary of all model results.
Design rules specified in EC3
The unfactored design strength of flexural members depends on the minimum effective section modulus depending upon the position of the neutral axis. According to EC3 (BS-EN1993-1-3), the unfactored design strength (MEC3) is calculated as follows
where fyb is the basic yield strength and Weff,y is the effective section modulus given by
where Weff, y,c and Weff, y,t are the section moduli with regard to compression and tension flanges, respectively; Z is the position of neutral axis from respective flange and Ieff,y is the second moment of area of the effective section.
The effective section properties are calculated by using a reduction factor for elements in compression
where ρ is the width reduction factor, λp,h is the relative slenderness, Kσ is the buckling factor, ψ is the stress ratio, hp is the nominal dimension of cross-sectional element and hc is the distance of the point of maximum compressive stress in element from neutral axis.
Design rules specified in AISI-S100
The unfactored design strength (Mn) of flexural members using the AISI specification (2016) is calculated as follows
where Fy is the nominal yield strength, Se is the elastic section modulus relative to top fibre, ycg is the depth of neutral axis with respect to the compression flange and Ix is the second moment of area of the effective section, determined by using a reduction factor, given by
where beff is the effective design width, w is the width of compression element, ρ is the reduction factor, k is the plate buckling co-efficient, t is the thickness of compression element, E is the modulus of elasticity and f is the maximum compressive edge stress in the element.
Comparison with design rules
Figure 15 shows the comparison between test strengths and design strength predictions of North American Standards and Eurocode. From Table 4 and Figure 15, it can be concluded that North American Standards are conservative for beams with a wall thickness of 2 mm, but un-conservative for the stiffened CFS beam with a wall thickness of 1 mm. A similar behaviour was indicated by Eurocode, except for Model III (where a slight degree of unconservativeness was observed). In all the four models, non-dimensional slenderness for lateral torsional buckling is less than limiting slenderness 0.4; therefore, there is no possibility of occurrence of lateral torsional buckling.

Comparison between test and design strengths for CFS beams.
Numerical calibration
The finite element analysis using ABAQUS (2004) version 6.14 was conducted to simulate an experimental behaviour of CFS lipped channel beams under four-point loading as shown in Figure 16. Nominal cross-sectional dimensions, material imperfections and initial geometric imperfections of the test specimens were incorporated. S4R5 shell element was selected to develop the finite element model. This element is thin, shear flexible, isometric quadrilateral shell with four nodes and five degrees of freedom per node, utilizing reduced integration and bilinear interpolation scheme (Ammash, 2017; Anbarasu, 2016; Keerthan and Mahendran, 2013). Mesh convergence study was carried out to find the optimum mesh size. The model with 5 mm mesh size provided reasonable accuracy and was hence used in all the finite element models (FEMs). The size of the element adopted was 5 mm × 5 mm (25) mm2. Elastic perfectly plastic model with a modulus of elasticity of 210 GPa and yield stress of 450 MPa was used in this study. To ensure proper distribution of concentrated forces on to the beams, the load was applied at the centre of the rigid plate attached to the beams as shown in Figure 16. Idealized simply supported end condition was modelled by restraining the displacements in x, y and z directions and rotations in z directions at the pinned support. The displacements were restrained in y and z directions and rotations in z directions at the roller support. To avoid contact problems in-between the layers, general hard surface contact was adopted. In the assembly, various instances of master–slave surfaces were created between the surfaces. Frictionless hard contact was adopted. One major problem while meshing sections with contact faces is that there are penetrations of layers during analysis. From various trials, it was found that square mesh of size 5 mm × 5 mm can be adopted for modelling all parts of the assembled sections to avoid any penetrations (Deepak and Shanthi, 2018). Residual stresses have a negligible effect on the strength (Schafer and Pekoz, 1998) and hence were ignored. Local, distortional and global geometric imperfections were incorporated in the model. The maximum magnitude of local, distortional and global imperfections adopted was 0.34 × t, 0.94 × t and L/1000, respectively (Schafer and Pekoz, 1998). In the nonlinear analysis, RIKS method was used.

Numerical model used in the finite element analysis.
The comparison of load versus deflection curves for test and FEM studies is shown in Figure 17. The mean value of the (PTest/PFEM) ratio is 0.987 with the corresponding standard deviation of 0.037, as shown in Table 8 and is also presented in Figure 18. Figure 19 shows the comparison of deformed shapes (test vs FEM) in Model II and Model III. Table 5 and Figures 17 to 19 indicate that the FEM-predicted results are in good agreement with test results, thus confirming the accuracy of the experimental results.

Comparison of test and FEM results.
Summary of all test and FEM results.
FEM: finite element model.
Average = 0.987.
Standard deviation = 0.037.

Verification of test results.

Comparison of deformed shapes (test vs FEM): (a) Model II and (b) Model III.
Conclusion
Based on experimental as well as numerical investigations carried out to study the effectiveness of stiffening arrangements in mobilizing the untapped reserve strength in CFS beams (comprising two channel sections connected back to back by bolts), the following important conclusions are drawn:
The judiciously provided steel angle stiffener which was welded to the inside of the compression flange lip only over the vulnerable zone considerably improved the load carrying capacity from 43 to 67.8 kN (i.e. increase in strength by 40%) and the initial stiffness from 3.58 to 5.71 kN/mm (i.e. increase by 60%). Hence, instead of using a heavy hot-rolled section or larger CFS section, a suitable smaller CFS section can be partially stiffened in the compression zone to avoid premature local buckling. This type of stiffening can be adopted both in existing structures (which demand stiffening) as well as in the structures to be built to obtain economy.
The stiffening of the compression flange lips in sections with flat width-to-thickness ratio of compression flange greater than 62.5 using small-channel section did not perform satisfactorily. It shows that the proposed stiffening arrangement may be suitable only when the flat width-to-thickness ratio of compression flange is not greater than 32.5.
The finite element models developed were in good agreement with the experimental results. The ratio of FEM-predicted strengths and test results had an average PTest/PFEM of 0.987 and standard deviation of 0.037. It shows that ABAQUS software can be used to predict the behaviour of such partly stiffened CFS beams.
Design strengths computed as per North American Specifications and European Code for CFS structures were found to be conservative, except for the stiffened CFS beam with a wall thickness of 1 mm. Hence, these provisions can be safely used for the design of such CFS beams.
Lip buckling and flange buckling in the compression zone were the primary types of observed failures, which can be effectively controlled by the proposed stiffening.
The minimum and maximum geometric imperfections observed at mid-lengths in the various models were 1/4346 mm and 1/2131 mm, respectively.
It has to be noted that in this study, angles were used to stiffen the compression flange of CFS beams. However, we can optimize the dimensions of the stiffener in order to obtain economy. On this aspect, a parametric study being carried out by the authors is under progress.
Footnotes
Appendix 1
Acknowledgements
The authors would like to thank the Civil Engineering Department of National Institute of Technology Srinagar for their support by permitting the testing of the models in their Structural Engineering Laboratory. Prior to joining IIT Delhi, M. Adil Dar was working as an MTech scholar in Structural Engineering at Kurukshetra University and wishes to thank the University for their support.
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: The experimental work described in this paper has been supported by a grant from Consulting Engineers, PVT. LTD (Project No. CES2015/8360).
