Civil Engineering Applications of Polymer CompositesView this Special Issue
Research Article | Open Access
Minkwan Ju, Hongseob Oh, "Experimental Assessment on the Flexural Bonding Performance of Concrete Beam with GFRP Reinforcing Bar under Repeated Loading", International Journal of Polymer Science, vol. 2015, Article ID 367528, 11 pages, 2015. https://doi.org/10.1155/2015/367528
Experimental Assessment on the Flexural Bonding Performance of Concrete Beam with GFRP Reinforcing Bar under Repeated Loading
This study intends to investigate the flexural bond performance of glass fiber-reinforced polymer (GFRP) reinforcing bar under repeated loading. The flexural bond tests reinforced with GFRP reinforcing bars were carried out according to the BS EN 12269-1 (2000) specification. The bond test consisted of three loading schemes: static, monotonic, and variable-amplitude loading to simulate ambient loading conditions. The empirical bond length based on the static test was 225 mm, whereas it was 317 mm according to ACI 440 1R-03. Each bond stress on the rib is released and bonding force is enhanced as the bond length is increased. Appropriate level of bond length may be recommended with this energy-based analysis. For the monotonic loading test, the bond strengths at pullout failure after 2,000,000 cycles were 10.4 MPa and 6.5 MPa, respectively: 63–70% of the values from the static loading test. The variable loading test indicated that the linear cumulative damage theory on GFRP bonding may not be appropriate for estimating the fatigue limit when subjected to variable-amplitude loading.
Fiber-reinforced polymer (FRP) bars have been used widely in reinforced concrete structures due to their many advantages; FRP bars have superior material properties, such as a high tensile strength and corrosion resistance. From a maintenance point of view, using FRP bars inside concrete structures is clearly cost-effective with respect to life cycle costs. Although the bond performance of FRP bars is weaker than that of steel bars, FRP bars have generated interest as an advanced substitute material for reinforced concrete structures.
The bond in reinforced concrete is important in transferring stress from the reinforcing bar to the concrete. To have a composite action in reinforced concrete, perfect bond capacity is required, and conventional steel bars are usually considered to satisfy this bond performance. Unlike conventional steel reinforcing bars, however, determining the performance bond capacity of FRP reinforcing bars is difficult. The bond between concrete and FRP reinforcing bars is complex and several factors influence the bond characteristics. According to previous research, the key factors are the concrete compressive strength, bar diameter, embedment length, geometry and surface treatment scheme of the FRP reinforcing bar, temperature changes, and environmental conditions .
Many experimental studies have been conducted to investigate the bond performance of FRP reinforcing bars using pullout tests [2–7]. Most of these have been monotonic and uniaxial tests based on the ASTM standard. Beam tests are believed to more realistically simulate the stress transfer of reinforced concrete subjected to flexural bending, but few studies have been reported determining the fatigue bond characteristics of FRP reinforcing bars, which is essential to understand reinforced concrete structures in a service state under repeated monotonic and amplitude loading. Investigating bond characteristics of the bond surface is important, wherein longitudinal friction of the FRP reinforcing bar is usually concentrated when subjected to repeated external loading of the reinforced concrete beam.
This study investigated the bond performance of GFRP reinforcing bars subjected to fatigue loading. The beam tests for bonding were carried out according to the BS EN 12269-1  specification. The test variables were the three types of bond length. For static and monotonic loading tests, maximum bond strength and slip were analyzed, and a reduction factor in bond strength was suggested experimentally to evaluate safety in bond design subjected to fatigue loading. For variable-amplitude loading tests, the fatigue limit was evaluated with an S-N curve based on Miner’s theory, and its application to the nonlinear cumulative damage for the bonding behavior of glass fiber-reinforced polymer (GFRP) reinforcing bars subjected to variable-amplitude loading is discussed.
2. Experimental Program
2.1. Description of Properties of GFRP Reinforcing Bars, Steel Reinforcing Bars, and Concrete
The GFRP reinforcing bar, as depicted in Figure 1, used consisted of continuous longitudinal glass fibers of 67% volume fraction in a thermosetting epoxy . To mold the surface pattern and to enhance the bond performance of the GFRP reinforcing bar, its external layer was manufactured by mixing milled glass fiber and epoxy at a ratio of 1 : 1, and it was cured for 15 min at a temperature above 160°C. To enhance bond performance, the surface pattern of the developed GFRP reinforcing bar was treated to have ribs similar to those of a steel reinforcing bar using epoxy resin containing milled glass fiber by mechanical pressing after pultrusion of the FRP core section. The rib section was a mixture of epoxy resin and milled glass fibers. The external deformed-rib arrangement angle was set at 80°. The rib height and spacing were chosen based on a literature review to obtain optimal mechanical performance in terms of bond and tension characteristics . The GFRP reinforcing bar used for the tensile area had a nominal diameter of 9.53 mm and a design tensile strength which is a guaranteed tensile strength multiplied by the environmental reduction factor in compliance with ACI 440 1R-06 of fu (GFRP) = 616 MPa; its modulus of elasticity was 42.9 GPa .
(a) Conventional steel reinforcing bar
(b) Newly developed GFRP reinforcing bar
(c) Design detail of GFRP reinforcing bar
2.2. Flexural Bond Strength from the Beam Test
The test method and details for evaluating the flexural bond strength of the GFRP reinforcing bar adopted the scheme of the British Standard . The beam test in bonding is considered to more realistically mimic the stress conditions of the reinforced concrete (RC) structure subjected to bending. Tighiouart et al.  conducted beam tests in bonding according to the RILEM specification. The specimen consisted of two rectangular blocks joined at the top by a steel ball joint. The design concrete compressive strength is 27 MPa. This specimen detail may not be appropriate for fatigue but for a monotonic loading test due to the possibility of the stress concentration in the concrete at the ball joint. Thus, for the fatigue bond test in this study, the test scheme of BS EN 12269-1  was adopted, whereby a concrete compressive block was used to create the stress distribution in the fatigue test.
The size of the test specimens was 180 × 200 × 1300 mm, and the major variables according to bond length () are shown in Figure 2. The lower section of the mid-span in the beam was formed as a semicircular type to prevent stress concentration caused by flexural cracks tips. Bond lengths were set at (45 mm), (90 mm), and (135 mm), where is the diameter of the FRP reinforcing bar, for simulating pullout failure of the GFRP reinforcing bar. As shown in Figure 2, the unbounded region of the FRP reinforcing bar in the concrete beam was secured by PVC pipes that had a slightly larger inner diameter than the diameter of the GFRP reinforcing bar. The test beam was loaded using a 250 kN hydraulic actuator (MTS) and was executed with a four-point loading method. The unloaded slip was measured at both ends of the specimens (Figure 2).
Figure 3 illustrates the force equilibrium condition, and its relationship is shown in (1). Once the tensile load () is applied, the tensile stress () at the load level is calculated by dividing by the bond surface area according to BS EN 12269-1 . The optimal bond length is defined as the minimum length to transfer the ultimate load from the GFRP reinforcing bar to the concrete through the bond surface. Considerwhere = tensile load of the GFRP reinforcing bar (kN), = ultimate applied load (kN), = measured strain of the steel or GFRP reinforcing bar, = modulus of elasticity of the steel or GFRP reinforcing bar (MPa), = shear span (mm), and = distance between the resulting tensile and compressive loads (mm).
As mentioned above, the fatigue loading applied maximum stress level with the ratio to the static test results as follows: 60%, 70%, 75%, 80%, 85%, and 90% of the peak load of each variable. The minimum stress level was considered to be 10% of the applied maximum load. All specimens were tested at 3 Hz, and sinusoidal loading was controlled up to fatigue failure or 2,000,000 loading cycles. Data acquisition for calculating bond stress was performed after finishing the cycle controls.
3. Test Results of the Static Loading Test
The relationship between flexural bond stress and slip from the static test, performed according to Oh et al. , is summarized in Figure 4 and Table 1. The bond stress of each specimen was expressed under force equilibrium conditions and reinforcing bar strain. According to the static loading test, the bond strength determined using the measured strain of the GFRP reinforcing bar was found experimentally to be lower than that using equilibrium theory. This difference might have been caused by incomplete bond stress transfer to the concrete. The reason was analyzed in terms of the bond behavior of the GFRP reinforcing bar and the change in the neutral axis of the beam section during loading.
As reported by other researchers, the bond strength calculated from force equilibrium at the beam section of the mid-span of the specimen is slightly higher than the strength calculated from FRP strain. The work energy of the pullout action of FRP reinforcing bar indicates the capacity for energy dissipation at the interface between concrete and the surface of reinforcing bar. While the bond strength decreases, according to the increase in bond length, total work energy, defined by the area of loads and slip curves of the specimen, increases (Figure 5).
The required development length to ensure reinforcing bar strength should be calculated under force equilibrium conditions for the beam specimen and also demonstrated by the equation , as reported in ACI 440 1R-03, where is the bond length, is the nominal diameter of the GFRP reinforcing bar, and is the designed tensile strength of the GFRP reinforcing bar. The experimental length obtained from the regression data based on the test results and development length according to ACI 440 is shown in Figure 6. Both the regression curve from test results and the analytical bond length from (1) curve tended to decrease the bond strength as the bond length increased. The large discrepancy at 50 mm bond length was due to the short bond length and bond perimeter, which could be conservative as the bond length is increased further. The minimum bond length that the GFRP reinforcing bar can afford to contribute to the tensile strength of the concrete was 225 mm from the crossing point between the two curves by regression and (1). The bond length, however, was calculated as 317 mm according to ACI 440 1R-03.
4. Test Results of the Monotonic Fatigue Loading Test
4.1. Bond Strength When Subjected to Fatigue Loading
Each specimen was tested for the designated level of stress as summarized in Table 2, and 10 specimens in total were tested. Two types of failure pattern were investigated: typical pullout failure and concrete splitting failure at the end of bonded region of the reinforcing bar () (Figures 7 and 8). Pullout failure occurs when the bond stress between the concrete and GFRP reinforcing bar exceeds the stress resulting from the tensile load on the bar. Concrete splitting failure must be caused by insufficient concrete strength, even though the embedded length of the reinforcing bar is sufficient to perform the composite action. As shown in Figure 6, no shear-off of the ribs on the GFRP reinforcing bar occurred with the pullout failure, indicating that the interlocking mechanism with the concrete is valid.
|Pullout failure and concrete splitting failure.|
The bond strength was determined at the fatigue limit of 2,000,000 cycles, at which point a failure test was finally applied. The and specimens showed the fatigue limit at a stress level of 70%. The number of cycles increased as the development length was longer. For the pullout failure specimens of and , the bond strengths at failure testing after 2,000,000 cycles were 10.5 MPa and 6.4 MPa, respectively, which is 63–70% of those from the static bond test. This difference was caused by repeated stress on the bond surface so that fatigue stress weakened the adhesive capacity between the concrete and surface of the GFRP bar, including ribs, so that it degraded the bonding performance. For the specimen, the mode of failure was concrete splitting. The bond length over can be concluded to be close to having sufficient bonding performance for the GFRP reinforcing bar. This is consistent with the results of the static loading test, whereby the bond strength difference between the calculated and measured strain on the GFRP reinforcing bar and that using equilibrium theory keeps decreasing as the bond length increases.
4.2. Fatigue Bond Strength and Slip Relationship
Bond stress and slip relationships under fatigue loading are shown in Figures 9 and 10, including the results of the static bond tests. The specimen was not included for slip behavior after 2,000,000 cycles due to the concrete splitting mode of failure. The slip amount at 2,000,000 cycles for the specimen was not measured due to a data recording error. According to Oh et al. , the bond performance of the steel reinforcing bar was better than that of the GFRP reinforcing bar, considering bars of the same diameter, due to enhanced shear action of the GFRP rib. For the fatigue test, as the number of repeated cycles increased, the residual slip amount also increased, although its scale was very small. The specimen showed less than 50% slip versus the specimen at the ultimate bond strength. As the bond length increased, the bond strength decreased [13, 14]. The bond strength can be defined as a mechanical performance of the rib. As the bond length is shorter, the bond strength at the rib is increased. If the bond length is longer, however, the bond strength at the rib is decreased due to the stress distribution. Consequently, the ribs with longer bond lengths can stably contribute stress transference to the concrete with much higher tensile stress of the reinforcing bar than with a short bond length.
4.3. Fatigue Limit of GFRP Reinforcing Bar Bonding
Comparing the specimen with the specimen at 2,000,000 cycles in Figures 8 and 9, the maximum slip of the specimen was found to be less than that of the specimen. This indicated that as the bond length was longer, the maximum bond strength resulted in the shorter slip. This was confirmed with experimental results of maximum slips of 0.175 mm for and 0.083 mm for . The slip limit for longer bond length specimen of the GFRP reinforcing bar at 2,000,000 cycles can be estimated; it will diminish with a decreasing rate close to that for the and specimens.
For the fatigue test, the fatigue limit state can be estimated from the S-N curve based on Miner’s theory, which is based mostly on linear cumulative damage concepts, as proposed by Palmgren and Miner . Figure 11 shows the log-scaled S-N curve and regression analysis from the results of the fatigue tests for the and specimens. For S-N curve analysis, the fatigue limit of the bonding of the GFRP reinforcing bar can be investigated for pullout failure. The fatigue limit state of bonding is useful in designing and analyzing flexural reinforced concrete members, especially those reinforced with GFRP reinforcing bars. Table 3 shows the estimated fatigue limit stress by regression analysis.
5. Test Results of the Variable-Amplitude Loading Test
Due to the uncertain characteristic of the bond strength of the FRP reinforcing bar, investigating only monotonic loaded fatigue test is not sufficient. Authors conducted the flexural fatigue bond test under variable-amplitude loading condition to evaluate the real bond fatigue performance. One can simulate the accelerated effect of accumulated damage on bond strength of the GFRP reinforcing bar by subjecting it to variable-fatigue loading conditions. Table 4 shows the applied load levels and cases for the variable-amplitude loading test. Each specimen was tested for the designated variable loading case, and in total six specimens were tested for and . Loading was applied based on the designated proportion of the maximum bond strength, as in the monotonic loading test. The loading sequence, however, was set as shown in Table 4. The minimum load level was considered to be 15% of the failure number of cycles () and the variable-amplitude load was initiated at this level. Variable load cases 1 and 2 are the upward (low 75% to high 80%) and downward (high 80% to low 75%) load cases, respectively. These loading cases were performed to investigate the accumulated damage effect according to the loading sequence magnitude. Variable load case 3 is an upward case including two steps of the minimum load level of 0.15 . The experimental results from these load cases were analyzed using the linear damage theory of fatigue.
Figure 11 shows the relationship of applied bond stress and unbounded end slip for test specimens with bond lengths of 45 mm and 135 mm. The two specimens have different performance limits, such as pullout and concrete splitting failure. For this reason, the unloaded end slips before the concrete splitting failure of the specimen were compared with those of the specimen. The specimen showed a sudden increase in slip in the downward load case (80 to 75) as it approached failure. The specimen, however, showed no such sudden increase and the slip at failure was stable at 0.2 mm. As the applied load level was increased, the specimen with the shorter bond length was more sensitive to causing a larger slip. Figure 12 shows the relationship between cumulative end slip and normalized relative fatigue life ().
Figure 13 illustrates the cumulative end-slip increase according to the relative fatigue life for the monotonic and the variable-amplitude loading test. The cumulative end slip under variable-amplitude loading condition, compared with that of monotonic loading condition, exhibited a sudden increase of slip at the moment of the specific fatigue life. In the case of the FRP reinforcing bar, verifying the bond characteristics based on bond length is more important than the conventional steel reinforcing bar subjected to the variable-amplitude load regarding the real loading frequency of the traffic volume. Table 5 summarizes the result of estimation of fatigue life by the Palmgren-Miner rule, which is commonly used to calculate the cumulative fatigue damage on the following equation: where = number of loading cycles at a given stress level sigma , = number of cycles to failure at sigma , and = total damage on fatigue (usually considered to be 1 at fatigue failure).
6. Results and Discussions
The bond length of the GFRP reinforcing bar in this study is “safer” than that mandated by ACI 440 1R-03, and the outer ribs can sufficiently resist the shear friction against the tensile stress of the GFRP reinforcing bar. For purposes of structural designs, the bond characteristics should also be confirmed when subjected to repeated loading conditions.
The work energy of the pullout action of FRP reinforcing bar indicates the capacity for energy dissipation at the interface between concrete and the surface of reinforcing bar. While the bond strength decreases, according to the increase in bond length, total work energy, defined by the area of loads and slip curves of the specimen, increases (Figure 5). This means that each bond stress on the rib is released and bonding force is enhanced as the bond length is increased. Appropriate level of bond length may be recommended with this energy-based analysis.
Compared with the static loading test, the bond strength decreased markedly for the and specimens, up to a maximum 63–70%, as noted above. This degraded bond strength was caused by the residual slip resulting from the repeated loading. Repeated loading of a flexural member reinforced with GFRP reinforcing bar was concluded to affect the loss of bond strength, which is a factor usually considered in the design of the bond.
At the assumed fatigue cycles of 10,000,000, the fatigue limit stresses for and specimens were 68.0% and 65.0%, respectively. These results showed that the GFRP reinforcing bar had sufficient bond performance in a repeated loading state, such as vehicular traffic, to be used in the design of concrete flexural members.
Initial damage was calculated by considering the applied load level of 0.15. For the total damage on fatigue, two apparent differences were detected. The low-high cases (cases 1 and 3) showed total damage over 1.0 according to test results of the fatigue failure. The high-low case (case 2), however, showed total damage still under 1.0 even though the specimen had already failed. This result indicates that the fatigue performance was more vulnerable when a high load is applied prior to a low load on the structure. The linear cumulative damage theory based on the Palmgren-Miner rule may not be appropriate for estimating the fatigue limit when subjected to variable-amplitude loading.
6.1. Reduction Factor of the Fatigue Bond Strength
From the fatigue bond test, the relationship between the reduction in bond strength and the residual slip at maximum bond strength in fatigue at 2,000,000 cycles was investigated. As the slip increased, the bond strength was found to decrease. Thus, the reduction in bond strength was proportional to the amount of residual slip, which means that bond stiffness, according to repeated loading, decreases. This relationship can be redefined using an energy concept to calculate half the area of the graph.
In this study, a reduction factor for bond stiffness of the FRP reinforcing bar, , can be defined by an energy concept using the area of the shaded section in Figure 14. The conceptual formula is as follows: is suggested for evaluating the reduced capacity in bonding of an FRP reinforcing bar reinforced in a concrete flexural member when it is serviced under an ambient repeated loading state such as vehicle load. The performance standard for fatigue bonding of the FRP reinforcing bar should be used in designing the member. The reduced capacity under fatigue loading must be evaluated quantitatively because of the various types of rib shapes for the FRP reinforcing bar so that the rib shapes directly affect determining the bond capacity, especially in fatigue loading. By applying (3) to the results of the and specimens, for the and specimens was found to be 0.41 and 0.16, respectively, which means that as the bond length is longer, the reduction in bond stiffness decreases. Thus, a quantitative evaluation indicated that the bond stress can be transferred to the concrete well under fatigue loading.
Using this factor for the reduced capacity of bonding of FRP reinforcing bars, one can check the reduction rate of bonding of FRP reinforcing bars under fatigue loading. For a more reliable evaluation, should be formulated as a function of various experimental data so that the nonlinearity of the relationship between the reduced bond strength and the residual slip can be confirmed.
In this study, we investigated the bond performance of GFRP reinforcing bars under fatigue loading. The fatigue test was conducted until 2,000,000 cycles, and the bond strength and slip relationship were examined. The conclusions are as follows.
(1) For the static loading test, the bond length was evaluated as 225 mm, whereas it was 317 mm according to ACI 440 1R-03. The bond length of the GFRP reinforcing bar in this study well satisfied the ACI 440 1R-03 limit, and the outer ribs sufficiently resisted the shear friction against the tensile stress of the GFRP reinforcing bar. Each bond stress on the rib is released and bonding force is enhanced as the bond length is increased. Appropriate level of bond length may be recommended with this energy-based test result. For the purposes of structural design, these bond characteristics should be further confirmed when subjected to repeated loading conditions.
(2) For the pullout failure specimens, the bond strengths at failure testing after 2,000,000 cycles were 10.4 MPa and 6.5 MPa, respectively, which is 63–70% when compared with those from the static bond test. This was caused by repeated stress on the bond surface so that fatigue stress weakened the adhesive capacity between the concrete and the surface of the GFRP reinforcing bar. Thus, FRP bonding in designs should be checked with respect to the fatigue behavior of the flexural member.
(3) In this study, a reduction factor for bond stiffness for the FRP reinforcing bar, , was suggested experimentally. Using this factor for the reduced capacity of bonding of FRP reinforcing bars, one can quantitatively check the expected reduced bonding of the FRP reinforcing bar under fatigue. For a more reliable evaluation, should be formulated as a function of various experimental data so that the nonlinearity of the relationship between the reduced bond strength and residual slip can be confirmed.
(4) For the variable loading test, the linear cumulative damage theory on GFRP bonding was found to perhaps not be appropriate to estimate the fatigue limit when subjected to variable-amplitude loading. The GFRP reinforcing bar was concluded to have sufficient bond performance in a repeated loading state, such as vehicular traffic; that is, it can be used in designing concrete flexural members. In future studies, fatigue limit estimation should be researched for bonding of GFRP reinforcing bars under ambient loading conditions.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by research grants of Korea Institute of Marine Science & Technology Promotion (PJT200493) and Korea Institute of Energy Technology Evaluation and Planning (0000000015513).
- E. Cosenza, G. Manfredi, and R. Realfonzo, “Behavior and modeling of bond of FRP rebars to concrete,” Journal of Composites for Construction, vol. 1, no. 2, pp. 40–51, 1997.
- E. Cosenza, G. Manfredi, and R. Realfonzo, “Development length of FRP straight rebars,” Composites, Part B: Engineering, vol. 33, no. 7, pp. 493–504, 2002.
- Z. Achillides and K. Pilakoutas, “Bond behavior of fiber reinforced polymer bars under direct pullout conditions,” Journal of Composites for Construction, vol. 8, no. 2, pp. 173–181, 2004.
- M. Baena, L. Torres, A. Turon, and C. Barris, “Experimental study of bond behaviour between concrete and FRP bars using a pull-out test,” Composites Part B: Engineering, vol. 40, no. 8, pp. 784–797, 2009.
- Z. He and G.-W. Tian, “Probabilistic evaluation of the design development length of a GFRP rod pull-out from concrete,” Journal of Engineering Structures, vol. 33, no. 10, pp. 2943–2952, 2011.
- C. G. Papakonstantinou, P. N. Balaguru, and Y. Auyeung, “Influence of FRP confinement on bond behavior of corroded steel reinforcement,” Cement and Concrete Composites, vol. 33, no. 5, pp. 611–621, 2011.
- Y. Ding, X. Ning, Y. Zhang, F. Pacheco-Torgal, and J. B. Aguiar, “Fibres for enhancing of the bond capacity between GFRP rebar and concrete,” Construction and Building Materials, vol. 51, pp. 303–312, 2014.
- British Standard, “Determination of the bond behavior between reinforcing steel and autoclaved aerated concrete by the beam test,” BS EN EN 12269-1, 2000.
- J. Sim, H. Oh, D. Moon, and M. Ju, “A hybrid fiber reinforced plastic rebar having a optic sensor for concrete,” Patent No. 100709292, 2006 (Korean).
- D. Moon, Bond behavior of newly developed deformed GFRP bars [Ph.D. thesis], Hanyang University, Seoul, Republic of Korea, 2004.
- H. Oh, J. Sim, T. Kang, and H. Kwon, “An experimental study on the flexural bonding characteristic of a concrete beam reinforced with a GFRP rebar,” KSCE Journal of Civil Engineering, vol. 15, no. 7, pp. 1245–1251, 2011.
- B. Tighiouart, B. Benmokrane, and D. Gao, “Investigation of bond in concrete member with fibre reinforced polymer (FRP) bars,” Construction and Building Materials, vol. 12, no. 8, pp. 453–462, 1998.
- E. Makitani, I. Irisawa, and N. Nishiura, “Investigation of bond inconcrete member with fibre reinforced polymer bars,” in Proceedings of the International Symposium Fibre-Reinforced-Plastic Reinforcement for Concrete Structures, ACI SP-138, pp. 315–331, 1993.
- A. Nanni, M. M. Al-Zahrani, S. U. Al-Dulaijan, C. E. Bakis, and T. E. Boothby, “Bond of FRP reinforcement to concrete,” in Proceedings of the 2nd International RILEM Symposium (FRPRCS-2 '95), pp. 135–145, London, UK, 1995.
- J. Schijve, Fatigue of Structures and Materials, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2001.
Copyright © 2015 Minkwan Ju and Hongseob Oh. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.