Development of Rock Embedded Drilled Shaft Resistance Factors in Korea based on Field Tests

: Load and resistance factor design (LRFD) is a limit state design method that has been applied worldwide. Because the data for determining LRFD factors in Korea has been insu ﬃ cient, the resistance factors suggested by American Association of State Highway and Transportation O ﬃ cials (AASHTO) in the US have been used for design in Korea; however, these resistance factors were deﬁned based on the characteristics of the predominant bedrock types in the U.S. As such, it remains necessary to determine resistance factors that reﬂect the bedrock conditions in Korea. Accordingly, in this study, LRFD resistance factors were determined using 13 sets of drilled shaft load test data. To obtain accurate resistance factors, calibration of the elastic modulus of the drilled shaft and the equivalent load–displacement curve considering the axial load and elastic settlement was conducted. After determining accurate resistance values, a reliability analysis was performed. The resistance factors were determined to be within 0.13–0.32 of the AASHTO factors for the shaft resistance, 0.19–0.29 for the base resistance, and 0.28–0.42 for the total resistance. This is equivalent to being 30–60% of the AASHTO-recommended values for the shaft resistance and 40–60% of the AASHTO-recommended values for the base resistance. These di ﬀ erences in resistance factors were entirely the result of discrepancies in the conditions of the rock in the US and Korea in which the shafts were founded.


Introduction
The load resistance factor design (LRFD) method can evaluate the resistance of a structure considering uncertainty based on the results of a reliability analysis. Using a reliability analysis, the probability of failure can be quantitatively evaluated, and the design of consistent and reliable geotechnical structures can be made possible. In recent decades, LRFD has replaced the allowable stress design (ASD) and working stress design (WSD) methods worldwide for the design of foundations [1,2]. The LRFD method was introduced in South Korea about a decade ago, and has been reflected in the design code for foundations in the country since 2012. The new Korean Highway Bridge Design Code (KHBDC) (limit state design) was published in 2012 and took effect for field applications starting in 2015. The resistance factors for foundations suggested by KHBDC [3], however, are based on American Association of State Highway and Transportation Officials (AASHTO) [4] in the U.S.; indeed, almost all the resistance factors were simply borrowed from AASHTO. Resistance factors for drilled shaft in rock from AASHTO were suggested as 0.50~0.55 for shaft resistance factors, and 0.50 for base resistances.
The resistance factors in AASHTO, however, reflect the characteristics of the bedrock found in the U.S., which may differ from the characteristics of the bedrock found in Korea. Notably, the resistance factors have already been calibrated for a number of states in U.S. based on area-specific bedrock types [5][6][7]. Also, the representative bedrock types in Korea are granite and gneiss classified as moderate or hard rock by strength criterion, they exhibit significant discontinuity and weathering. As such, these rocks are classified as weathered or soft rock by the weathering degree criterion [8]. This classification constitutes a discrepancy between the bedrock conditions in the U.S. and Korea, and as such, the direct application of the AASHTO resistance factors to Korean domestic design is inappropriate.
Considering these discrepancies, a great deal of research into resistance factors has been conducted in Korea. In 2008, the Korea Institute of Construction Technology [9] suggested total resistance factors for driven steel pipe piles (0.29~0.38 for target reliability index of 2.33) and drilled shafts (0.20~0.70 for target reliability index of 3.0). Jung [10] proposed total, shaft, and base resistance factors for drilled shafts using a numerical analysis. From Jung's research, resistance factors for target reliability index of 3.0 were determined as 0.55~0.80 for total resistance, 0.20 for shaft resistance, and 0.20~0.70 for base resistance. Kim et al. [11] suggested the base and shaft resistance factors for axially loaded driven piles in sands based on the Imperial College pile (ICP) design method. The calculated resistance factors (target reliability index of 3.0) were 0.61~0.86 for base resistance and 0.55~0.96 for shaft resistance. Park [12,13] updated the total resistance factors for driven steel pipe piles based on Bayesian theory to incorporate the results of proof tests. The updated resistance factors were determined as 0.19~0.71 for target reliability index of 3.0. Recently, there were some researches to suggest the resistance factors for deep foundations using load test data in domestic region [14,15] and researches for suggestion of resistance factors and comparison according to design methods for shallow foundations [16,17].
Most of this previous research, however, has been limited to suggesting the total resistance factors, and the validity of the shaft and base resistance factors suggested by Jung are somewhat doubtful as they are derived from the results of a numerical analysis. Therefore, it remains necessary to determine accurate resistance factors for drilled shafts considering the discrepancy between the rock conditions in the U.S. and Korea using reliable load test results.
In this study, 13 sets of drilled shaft load test data obtained using strain gauges were collected and a load transfer analysis was conducted to determine accurate values for the shaft and base resistances. When the load tests were conducted as bi-directional load tests, the equivalent load-displacement curve was drawn to determine the total resistance. To obtain sufficiently accurate resistance values, calibrations of the elastic modulus of the concrete drilled shaft and of the equivalent load-displacement curve considering the axial load and elastic settlement were conducted. After determining accurate resistance values, a reliability analysis was performed to determine the target reliability index and resistance factors using the advanced first-order second-moment (AFOSM) reliability method.

Collected Load Test Data
It is critical that accurate resistance values be determined. In the procedure for reliability-based resistance factor calibration, the measured resistance value can significantly influence the results of the reliability analysis and thus the subsequent the resistance factor calibration. Therefore, load test data obtained by strain gauges located in test drilled shafts were collected to provide accurate measurements of the shaft and base resistance. Based on these accurate shaft, base, and total resistance values obtained from the load transfer analysis, sufficiently accurate resistance factors can be determined.
In this study, a total of 13 sets of in-situ drilled shaft load test results were collected from three sites: shaft tests TP1-TP4 and TP5-TP8 were conducted on two different bridge construction sites and shaft  tests TP9-TP13 were conducted on a building construction site (see Figures 1-3 for site characteristics).  Sites 1 and 2, corresponding to TP1-TP4 and TP5-TP8, respectively, contained large-diameter shafts,  while Site 3 contained medium-diameter shafts. The drilled shaft profiles and rock properties, including  the unconfined strength, RQD, and TCR were obtained from the site investigations, and data from the  drilled shaft rock sockets were utilized to capture the conditions and effects of weathered, soft rock.  Detailed information describing shaft construction and rock conditions can be found in Table 1.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 18 In this study, a total of 13 sets of in-situ drilled shaft load test results were collected from three sites: shaft tests TP1-TP4 and TP5-TP8 were conducted on two different bridge construction sites and shaft tests TP9-TP13 were conducted on a building construction site (see Figure 1-3 for site characteristics). Sites 1 and 2, corresponding to TP1-TP4 and TP5-TP8, respectively, contained largediameter shafts, while Site 3 contained medium-diameter shafts. The drilled shaft profiles and rock  properties, including the unconfined strength, RQD, and TCR were obtained from the site  investigations, and data from the drilled shaft rock sockets were utilized to capture the conditions  and effects of weathered, soft rock. Detailed information describing shaft construction and rock  conditions can be found in Table 1.   In this study, a total of 13 sets of in-situ drilled shaft load test results were collected from three sites: shaft tests TP1-TP4 and TP5-TP8 were conducted on two different bridge construction sites and shaft tests TP9-TP13 were conducted on a building construction site (see Figure 1-3 for site characteristics). Sites 1 and 2, corresponding to TP1-TP4 and TP5-TP8, respectively, contained largediameter shafts, while Site 3 contained medium-diameter shafts. The drilled shaft profiles and rock properties, including the unconfined strength, RQD, and TCR were obtained from the site investigations, and data from the drilled shaft rock sockets were utilized to capture the conditions and effects of weathered, soft rock. Detailed information describing shaft construction and rock conditions can be found in Table 1.       Over the last few decades, the bi-directional load test has been widely used to test drilled shafts due to its advantages compared to conventional top-down pile load tests. Among the load tests considered in this study, bi-directional load tests were used in Sites 1 and 2 in order to apply sufficient loading for the large-diameter drilled shafts. However, the shaft diameters in Site 3 were smaller, so bi-directional tests were not necessary to provide sufficient loading. Therefore, conventional pile (top-down) static loading tests were applied. The test programs are summarized in Table 2.

Procedures for Determination of Measured Resistance
Measured resistances were determined using the ultimate bearing capacities obtained from the load tests. In general cases, the measured resistance indicates the total measured resistance from the load-displacement curve of a static load test or the equivalent load-displacement curve of a bi-directional load test. However, load-displacement curves are generally drawn using the applied load and measured displacement of the shaft head, so they cannot include detailed ground properties or axial load transfer characteristics of a drilled shaft in the loading stage. To overcome these limitations and check the shaft and base resistances according to depth, a load transfer analysis using the measured strain data was performed.
The procedure for the calculation of measured resistance can be summarized as follows: From the applied load and displacement measured during the load test, draw the loaddisplacement curve and determine the total resistance using this curve. If the load test was a bi-directional load test, draw the equivalent load-displacement curve If it is necessary to calculate the shaft and base resistances, a load transfer analysis is conducted using data from strain gauges installed along the pile shaft that report the strain at each evaluated depth.
Calculate the axial load for each depth using a load transfer analysis as follows: where P is the calculated axial load at each depth, ε is measured strain at each depth, E pile is the elastic modulus of pile, and A pile is the area of the test shaft. From the calculated axial load and displacement for each depth, draw the shaft resistancedisplacement curve (f-w curve) and base resistance-displacement curve (q-w curve) and determine the shaft and base resistances using these curves.

Calibration for Load Test Transfer Analysis
The axial loads and the shaft and base resistance values can be determined by a load transfer analysis using strain data in conjunction with applied load levels. In a load transfer analysis, the elastic modulus of concrete is one of the most important parameters. Typically, a constant elastic modulus (E 50 ) has been applied in the past. However, the elastic modulus of concrete exhibits nonlinear stress-strain characteristics, and as such, nonlinearity should be considered in any load transfer analysis. According Appl. Sci. 2019, 9,2201 6 of 17 to Fellenius [18], the stress-strain relationship of a concrete pile shows nonlinear behavior when loaded, as can be seen in Figure 4. Therefore, the tangent modulus of the stress-strain relationship of a concrete pile is not constant but rather decreases with the downward slope of the curve with increasing strain.
The axial loads and the shaft and base resistance values can be determined by a load transfer analysis using strain data in conjunction with applied load levels. In a load transfer analysis, the elastic modulus of concrete is one of the most important parameters. Typically, a constant elastic modulus (E ) has been applied in the past. However, the elastic modulus of concrete exhibits nonlinear stress-strain characteristics, and as such, nonlinearity should be considered in any load transfer analysis. According to Fellenius [18], the stress-strain relationship of a concrete pile shows nonlinear behavior when loaded, as can be seen in Figure 4. Therefore, the tangent modulus of the stress-strain relationship of a concrete pile is not constant but rather decreases with the downward slope of the curve with increasing strain.
(a) stress-strain relationship (b) tangent modulus of the stress-strain relationship In this study, the elastic modulus was calculated using the method of Fellenius, assuming the stress-strain relation of concrete to be a quadratic function. Then the calculated elastic modulus was applied to the load transfer analysis. For the calibration of the elastic modulus, the strain data near the load cell was used due to the large strain range present in the data. An example of the calibrated elastic modulus (TP 1) is shown in Figure 5. The solid line denotes the typically applied constant elastic modulus, while the dotted line describes the calibrated elastic modulus of the test drilled shaft. In this study, the elastic modulus was calculated using the method of Fellenius, assuming the stress-strain relation of concrete to be a quadratic function. Then the calculated elastic modulus was applied to the load transfer analysis. For the calibration of the elastic modulus, the strain data near the load cell was used due to the large strain range present in the data. An example of the calibrated elastic modulus (TP 1) is shown in Figure 5. The solid line denotes the typically applied constant elastic modulus, while the dotted line describes the calibrated elastic modulus of the test drilled shaft.

Calibration for Load Test Transfer Analysis
The axial loads and the shaft and base resistance values can be determined by a load transfer analysis using strain data in conjunction with applied load levels. In a load transfer analysis, the elastic modulus of concrete is one of the most important parameters. Typically, a constant elastic modulus (E ) has been applied in the past. However, the elastic modulus of concrete exhibits nonlinear stress-strain characteristics, and as such, nonlinearity should be considered in any load transfer analysis. According to Fellenius [18], the stress-strain relationship of a concrete pile shows nonlinear behavior when loaded, as can be seen in Figure 4. Therefore, the tangent modulus of the stress-strain relationship of a concrete pile is not constant but rather decreases with the downward slope of the curve with increasing strain.
(a) stress-strain relationship (b) tangent modulus of the stress-strain relationship In this study, the elastic modulus was calculated using the method of Fellenius, assuming the stress-strain relation of concrete to be a quadratic function. Then the calculated elastic modulus was applied to the load transfer analysis. For the calibration of the elastic modulus, the strain data near the load cell was used due to the large strain range present in the data. An example of the calibrated elastic modulus (TP 1) is shown in Figure 5. The solid line denotes the typically applied constant elastic modulus, while the dotted line describes the calibrated elastic modulus of the test drilled shaft. After calibration, the elastic moduli of the shafts are applied to load transfer analyses, which were conducted using the measured strain data. First, the axial loads at each depth were determined using the data from the strain gauges as follows: where P is the calculated axial load at each depth, ε is measured strain at each depth, E Fellenius is the calibrated elastic modulus, and A pile is the area of the test shaft. Using the calibrated load from Equation (2) for each depth, the load transfer analysis was performed to determine the measured shaft and base resistances.

Calibration for the Equivalent Load-Displacement Curve
As described in the previous section, eight load tests were bi-directional and five were conventional. However, as the bi-directional drilled shaft load test relies upon a loading mechanism entirely different from that of the conventional pile load test, many investigators and practicing engineers have expressed concern that the bi-directional load test could produce inaccurate results, especially describing drilled shaft head settlement behavior. According to Kwon et al. [19], the top-down equivalent curve constructed from the bi-directional load test results predicted the drilled shaft head displacement to be about half that predicted by the conventional top-down load test. This is due to the difference in the load transferred to the soil layer. In the case of the bi-directional load test, the applied load is transferred to the rock layer, whose elastic modulus is larger than that of the soil layer, while in the case of the top-down load test, the load is transferred to the soil layer, which has a small elastic modulus and is much softer than the rock. The neglected elastic displacement is another reason for the potential inaccuracy of bi-directional load test results. In the procedure for obtaining the equivalent load-displacement curve, it is assumed that the drilled shaft is incompressible, so the elastic compression of the shaft is neglected during the construction of the equivalent load-displacement curve. As a result, many researchers including Yu et al. [20] have insisted that the additional elastic compression of the shaft be considered in a complete analysis. Accordingly, Kwon et al. suggested the use of the newly developed plotting method for the equivalent load-displacement curve considering the load transferred to the soil layer and the elastic compression.
To plot the equivalent load-displacement curves in this study, the axial load distributions were determined for each test drilled shaft and the load ratios quantifying the difference in load distribution between the top-down test and the bi-directional load test were calculated according to Kwon et al. The shaft displacement can be calculated from the end bearing load as follows: where Q pile is the end bearing load for the downward movement of the pressure cell base w; A pile and L pile are the cross-sectional area and the length of the shaft, respectively; and E pile is the Young's modulus of the shaft. Taking as many random points as required, the top-down equivalent curve can be constructed by connecting the points described by the calculated load-displacement coordinates. For example, in the case of shaft TP 2, the load ratio (= ABCEA/ ABCDA) was approximately 3.82, indicating that the shaft head settlement should be w + 3.82∆w s to induce the downward movement of w at the shaft toe when the shaft is loaded from the top ( Figure 6). The calculated calibrated elastic moduli and load ratios for all test cases are summarized in Table 3.    Using the resulting axial load distribution and appropriately reduced elastic modulus, the equivalent load-displacement curves for each test shaft were constructed and compared with the results of the conventional load testing method. In this study, the equivalent load-displacement curves used by Kwon et al. and the calibrated elastic modulus used by Fellenius were selected to measure the resistance of the test drilled shafts. An example of the equivalent load-displacement curve is shown in Figure 7.  Using the resulting axial load distribution and appropriately reduced elastic modulus, the equivalent load-displacement curves for each test shaft were constructed and compared with the results of the conventional load testing method. In this study, the equivalent load-displacement curves used by Kwon et al. and the calibrated elastic modulus used by Fellenius were selected to measure the resistance of the test drilled shafts. An example of the equivalent load-displacement curve is shown in Figure 7.

Measured Resistance from the Load Test Results
A load transfer analysis was conducted to measure the shaft and base resistance based on the drilled shaft load test results. The strain data were utilized to calculate the axial loads at each depth along the test shaft. In addition, the f-w curves, describing the unit shaft displacement for each depth, and the q-w curves, describing the unit displacement at the base of the drilled shafts, were established using the calculated axial loads. Then, according to O'Neill and Reese [21], the ultimate shaft and base resistance values were determined when the displacement of the shaft was equal to 1% and 5% of the shaft diameter, respectively. If the load tests were terminated before the applied loads reached these ultimate resistance values, the f-w and q-w curves were extrapolated using the hyperbolic method [22], as shown in Figure 8. along the test shaft. In addition, the f-w curves, describing the unit shaft displacement for each depth, and the q-w curves, describing the unit displacement at the base of the drilled shafts, were established using the calculated axial loads. Then, according to O'Neill and Reese [21], the ultimate shaft and base resistance values were determined when the displacement of the shaft was equal to 1% and 5% of the shaft diameter, respectively. If the load tests were terminated before the applied loads reached these ultimate resistance values, the f-w and q-w curves were extrapolated using the hyperbolic method [22], as shown in Figure 8.  To determine the ultimate total resistance values, equivalent load-displacement curves were applied using the ASCE 20-96 method [23]. In this method, the ultimate total resistance is defined as the load at which the displacement of the shaft is equal to the elastic displacement + 3.81 + D/100 (mm). If a load test was terminated before the applied load reached this ultimate resistance value, the equivalent load-displacement curve was extrapolated using the hyperbolic method shown in Figure  9.

Predicted Resistance Values Obtained using Bearing Capacity Equations
Once the measured resistance values were obtained from the load test data, the predicted resistance values were used to evaluate the pile resistance uncertainty. While the measured resistance values were determined by analyzing the load test data, the predicted resistance values were calculated using bearing capacity equations. In this study, because the test piles were embedded in rock, the predicted resistance was calculated using a bearing capacity equation and the unconfined To determine the ultimate total resistance values, equivalent load-displacement curves were applied using the ASCE 20-96 method [23]. In this method, the ultimate total resistance is defined as the load at which the displacement of the shaft is equal to the elastic displacement + 3.81 + D/100 (mm). If a load test was terminated before the applied load reached this ultimate resistance value, the equivalent load-displacement curve was extrapolated using the hyperbolic method shown in Figure 9.
along the test shaft. In addition, the f-w curves, describing the unit shaft displacement for each depth, and the q-w curves, describing the unit displacement at the base of the drilled shafts, were established using the calculated axial loads. Then, according to O'Neill and Reese [21], the ultimate shaft and base resistance values were determined when the displacement of the shaft was equal to 1% and 5% of the shaft diameter, respectively. If the load tests were terminated before the applied loads reached these ultimate resistance values, the f-w and q-w curves were extrapolated using the hyperbolic method [22], as shown in Figure 8.  To determine the ultimate total resistance values, equivalent load-displacement curves were applied using the ASCE 20-96 method [23]. In this method, the ultimate total resistance is defined as the load at which the displacement of the shaft is equal to the elastic displacement + 3.81 + D/100 (mm). If a load test was terminated before the applied load reached this ultimate resistance value, the equivalent load-displacement curve was extrapolated using the hyperbolic method shown in Figure  9.

Predicted Resistance Values Obtained using Bearing Capacity Equations
Once the measured resistance values were obtained from the load test data, the predicted resistance values were used to evaluate the pile resistance uncertainty. While the measured resistance values were determined by analyzing the load test data, the predicted resistance values were calculated using bearing capacity equations. In this study, because the test piles were embedded in rock, the predicted resistance was calculated using a bearing capacity equation and the unconfined Figure 9. Determination of the ultimate total resistance value at TP 1.

Predicted Resistance Values Obtained using Bearing Capacity Equations
Once the measured resistance values were obtained from the load test data, the predicted resistance values were used to evaluate the pile resistance uncertainty. While the measured resistance values were determined by analyzing the load test data, the predicted resistance values were calculated using bearing capacity equations. In this study, because the test piles were embedded in rock, the predicted resistance was calculated using a bearing capacity equation and the unconfined compressive strengths of the rock cores, as these values could be easily obtained through field and laboratory tests. Note that the selected bearing capacity equations are currently widely applied in actual design and specifications. The selected bearing capacity equations, required parameters, and remarks for the predicted shaft, base, and total resistance values are summarized in Table 4 [24][25][26][27].

Load Statistics
Contrary to the resistance statistics, it is known that regional variations in the load statistics are negligible. As a result, in this study, it was assumed that the load statistics in Korea are not significantly different from those in the U.S., meaning that the statistical load characteristics and load factors from the current AASHTO LRFD Specifications could be used. The dead and live load combination for the Strength I case was selected because this combination is considered the most conservative for the calibration of the resistance factors of axially loaded piles and drilled shafts, and is often the governing design load case in bridge design. The load factors and load statistical parameters for the Strength I case are summarized in Table 5. Table 5. Statistical characteristics of the loads.

Resistance Statistics
The uncertainty of the ground conditions, design parameters, bearing capacity equations, and the drilled shaft resistance is established using a bias factor. The statistical characteristics of resistance are represented in terms of this bias factor, which is defined as the ratio of the measured resistance to the predicted resistance. The measured and predicted resistances were determined from the load test and bearing capacity equations as mentioned in previous sections. In this study, 13 load test cases were collected and used. These load tests were planned to install the strain gauge by depth and perform the load transfer analysis. From load transfer analysis, shaft and base resistances were separately determined and summarized as measured resistances in Tables 6-8. However, due to some reasons like site condition, pile construction quality and proficiency of engineer, load transfer analysis was not able to perform in some cases especially Site 3(TP 9~TP 13). The strain gauge condition of test piles in Site 3(TP 9~TP 13) were not good to measure the strain and perform the load transfer analysis, it was impossible to determine the shaft resistances. So, only base resistances for TP 9~TP13 were used for this study. Consequently, 21 bias factors for shaft resistance, 10 for base resistance, and 6 for total resistance were calculated and are summarized in Tables 6-8.  These calculated bias factors were then used to determine the resistance statistics, such as the mean, standard deviation, and coefficient of variation. These resistance statistics are summarized in Table 9. In this procedure, extremely outlying data points may not be representative of the resistance due to some large error; therefore, the bias factor values outside the boundary defined by the mean plus two times the standard deviation were omitted according to Paikowsky [28]. In case of Carter & Kulhawy, statistics based on all data were calculated as follows; mean as 1.55, Standard deviation as 1.11, and COV as 0.72. However, boundary of bias factors was calculated as 3.77 (mean plus times of standard deviation), so outlying data (bias factor of TP 2 as 3.80 in Table 7) was omitted for statistical analysis. Consequently, statistics were calculated based on nine data (excluding TP 2), and results were summarized in Table 9.

Determination of Target Reliability Index
To consider the transition from ASD to LRFD, a reliability analysis was performed and the reliability indices of shafts designed using the current design practice (ASD) were estimated. The Korean design standard for foundation structures requires a safety factor of 3.0 for pile bearing capacity design; therefore, the reliability analysis was performed for safety factors of 2.0-5.0 using the advanced first-order second moment (AFOSM) method.
The AFOSM method, named so because it is based on a first-order Taylor series approximation of the limit state function and because it uses only the second-moment statistics (mean and COV) of the variables, was proposed by Hasofer and Lind [29] and is known as one of the most accurate methods of reliability analysis. When the AFOSM method is used, the limit state function is linearized at the design point on the nonlinear failure surface, so the full distribution of information can be used rather than simply the mean value of the random variables. The results of the AFOSM reliability analysis are summarized in Table 10. Recently, AASHTO LRFD [1] selected a target reliability index of 3.0 for a typical drilled shaft. Moreover, KICT [5] suggested a target reliability index between 3.0-3.5 for drilled shafts. Based on the reliability analysis results and previous research [1,5], target reliability indices (β T ) of 2.5, 3.0, and 3.5 were selected for the calibration of the resistance factors in this study.

Determination of Resistance Factors
Considering three different target reliability indices (β T = 2.5, 3.0, and 3.5), the resistance factors were calibrated using the AFOSM method. The calibration results were expressed as the relationship between the reliability index and the resistance factor. The results are shown in Figure 10, and the resistance factors for each bearing capacity equation and target reliability index are summarized in Table 11.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 15 of 18 Considering three different target reliability indices (β = 2.5, 3.0, and 3.5), the resistance factors were calibrated using the AFOSM method. The calibration results were expressed as the relationship between the reliability index and the resistance factor. The results are shown in Figure 10, and the resistance factors for each bearing capacity equation and target reliability index are summarized in Table 11.      These differences in resistance factors are the result of the discrepancy in the conditions of the rocks in which the drilled shafts are founded. The resistance factors suggested by AASHTO are determined based on sedimentary and intact rock conditions, while according to geologic maps from Korea Structure and Foundation Structure Design Specifications [30], the representative types of bedrock in Korea are granite and gneiss, classified as weathered or soft rock by the weathering degree criterion. Clearly, this difference in rock condition has a clear and detrimental effect on drilled shaft performance that must be appropriately accounted for in the LRFD design resistance factors used in Korea. Also, suggested resistance factors were based on only 13 load test cases considering regional characteristics in Korea. But, if additional load test data will be collected, LRFD methods can be updated and calibrated the resistance factors based on reliability analysis. So, results in this study are expected to be helpful for future reference.

Conclusions
The load resistance factor design (LRFD) method can be used evaluate the resistance of a structure considering uncertainty based on a reliability analysis. In recent decades, LRFD has replaced allowable stress design (ASD) and working stress design (WSD) for the design of foundations worldwide, and is reflected in the design code for foundations in South Korea. However, the resistance factors used in Korea are based on those recommended by AASHTO in the U.S., which reflect characteristics of the bedrock in the U.S. that may differ from those of the bedrock in Korea. Based on such discrepancies in rock conditions, the direct application of the AASHTO-recommended resistance factors for design in Korea is inappropriate. Therefore, it was necessary to determine accurate resistance factors for drilled shafts considering the rock conditions in Korea based on reliable load test results.
In this study, 13 sets of drilled shaft load test data were obtained using strain gauges, and a load transfer analysis was performed to determine the accurate shaft and base resistance values. For bi-directional load tests, the equivalent load-displacement curve was drawn to determine the total resistance. Calibrations of the elastic modulus of the drilled shaft concrete and of the equivalent load-displacement curve considering the axial load and elastic settlement were conducted to obtain more accurate resistance values. After determining accurate resistance values, a reliability analysis was performed to determine the target reliability index and the resistance factors using the advanced first-order second-moment (AFOSM) reliability method.
For the AASHTO-recommended target reliability index of 3.0, the shaft resistance factors were found to be within 0.13-0.32 of the AASHTO-based values, the base resistance factors were within 0.19-0.29 of the AASHTO-based values, and the total resistance factors were within 0.28-0.42 of the AASHTO-based values for each bearing capacity equation evaluated. The resistance factors were accordingly determined to be 30-60% of the AASHTO-recommended values for the shaft resistance and 40-60% of the AASHTO-recommended values for the base resistance. These differences in resistance factors are the result of discrepancies in the condition of the rocks in which the drilled shafts were founded. The resistance factors suggested by AASHTO are determined using sedimentary and intact rock conditions prevalent in the US, while the bedrock in Korea tends to be weathered or soft granite and gneiss. This study has thus determined adjusted resistance factors to account for this difference in rock type and condition, providing improved design safety and load resistance accuracy when using drilled shafts in Korea. Also, suggested resistance factors were based on only 13 load test cases considering the regional characteristics in Korea. But, if additional load test data is collected, LRFD