Improved Wave Equation Analysis for Piles in Soil-Based Intermediate Geomaterials with LRFD Recommendations and Economic Impact Assessment

: The Wave Equation Analysis of Pile Driving (WEAP) has been widely used to determine drivability, predict static resistance, and assure the integrity of piles in soils. Assigning static and dynamic properties of Soil-based Intermediate Geomaterials (S-IGMs) remains a challenge in WEAP, partly attributed to IGMs that act as transition geomaterials between soil and hard rock. Furthermore, reliable static analysis methods for unit resistance predictions are rarely available for driven piles in S-IGMs in the default WEAP method. To alleviate these challenges, this study presents improved WEAP methods for steel piles driven in S-IGMs, including proposed damping parameters and Load and Resistance Factor Design (LRFD) recommendations based on newly developed static analysis methods and the classification of S-IGMs. A back calculation approach is used to generate the appropriate damping parameters for S-IGMs for three distinct subsurface conditions utilizing a database of 34 steel H-and pipe piles. Newly developed WEAP and LRFD procedures are also recommended. Additional independent 22 test pile data are used to compare and evaluate the accuracy and efficiency of the proposed WEAP methods with the default WEAP method. Compared with the default WEAP, bearing graph analysis results revealed that the selected proposed WEAP method, on average, reduces the underprediction of pile resistances by 6% and improves the reliability with a 43% reduction in the coefficient of variation (COV). Calibrated resistance factors for the proposed WEAP method increase to as high as 0.75 compared to the current AASHTO recommendation of 0.50. An economic impact assessment reveals that the proposed WEAP method is more efficient than the default WEAP method as the average difference in steel weight for 32 test piles is 0.06 kg/kN, almost close to zero, reducing the construction challenges in the current engineering practice.


Introduction
In rocky mountain areas, pile foundations are frequently driven in intermediate geomaterials (IGMs) to achieve higher resistances to meet the need for loads, satisfy the requirements of the Load and Resistance Factor Design (LRFD), increase design effectiveness, and reduce construction difficulties.Transitional geomaterials, known as IGMs, have a high level of natural variability.However, the definition of these geomaterials varies depending on the different applications [1].According to Clarke and Smith [2], Marinos [3], De Freitas [4], and Martin and Stacey [5], most definitions are based on unconfined compressive strength.The design and construction of bridge foundations are difficult and challenging due to significant geological variability and a need for standardized definitions of IGM [6].The lack of dependable design methodologies and measured IGM characteristics hindered reliable prediction of geotechnical resistance of driven piles in IGMs [7][8][9][10][11].
During the construction stage, pile performance in IGM is frequently assessed utilizing dynamic testing, such as the Pile Driving Analyzer (PDA) with a subsequent signalmatching analysis using the Case Pile Wave Analysis Program (CAPWAP) [12,13].However, once a test or production pile is equipped with PDA sensors and driven by a pile hammer during construction, the PDA-CAPWAP analysis can be carried out.Before construction, pile-driving analysis using the 1D Wave Equation Analysis Program (WEAP) was carried out to determine pile driving and acceptance criteria to identify the appropriate penetration depth that would accomplish the desired pile resistance using GRLWEAP14 software.Furthermore, WEAP has been frequently employed as an affordable construction approach due to the increased expense of a pile load test.WEAP is commonly used to assess pile drivability, predict static resistance, and ensure the integrity of piles in geomaterials.This program simulates pile penetration for each hammer strike.It estimates total pile resistance, compressive and tensile stresses in a pile, pile drivability, and pile integrity by modeling the driving system, pile, and geomaterials.One of WEAP's key outputs, the bearing graph, calculates the ultimate pile resistance as a function of blow count.The pile type, drive mechanism, percentage, and distribution of shaft resistance, and, most crucially, the dynamic soil factors, quake, and damping values all affect the bearing graph analysis.The proper selection of damping and quake values in the bearing graph analysis is crucial for the creation of pile driving criteria.Despite the quality of the modeling of the pile and drive systems, the damping and quake values are neither inherent geomaterial qualities nor quantifiable from laboratory measurements.The damping and quake values have been established from historical pile load test data for soils but not for the soil-based IGM (S-IGM).Previous studies examined pile drivability on soil to predict geotechnical resistances [14][15][16][17].In those studies, the observed pile resistance from static load tests [18,19] or dynamic load tests [20][21][22][23] were matched with the ultimate pile resistance from WEAP using a back-analysis approach.Those approaches are only limited to driven piles in soil, not S-IGM, which is different from soils as the classification system has been created for distinguishing S-IGMs from soil to improve the prediction of pile resistances [7,24].The study conducted by Masud et al. [7] shows that if the undrained shear strength (s u ) of any fine-grained geomaterials is equal to or greater than 129 kPa, then it would be classified as a fine-grained soil-based IGM (FG-IGM).Furthermore, when an SPT N-value is equal to or greater than 58 blows/0.3m, coarse-grained soil-based IGM (CG-IGM) is classified for pile resistance predictions [24].
The current WEAP procedure possesses challenges with assigning the static unit resistances and dynamic parameters for S-IGMs.This challenge is particularly true for CG-IGMs that frequently have Standard Penetration Test (SPT) N-values higher than the WEAP maximum permissible 60 blows/0.3m.Additionally, there are no reliable static analysis methods in default GRLWEAP14 software for estimating unit resistance for driven piles in S-IGMs, and the suggested Smith parameters were first created using pile load test data in soils.In this study, the newly developed static analysis methods for driven piles in S-IGMs are adopted from the study conducted by Masud et al. [7,24].New static analysis methods for predicting unit shaft resistance (q s ) and unit end bearing (q b ) in S-IGM were developed for different subgroups of CG-IGM and FG-IGM to reduce the uncertainties associated with the prediction of pile resistances in S-IGMs [7,24].For instance, CG-IGM is further differentiated based on pile types (steel H piles and steel pipe piles), and FG-IGM is further divided into clay-based IGM and silt-based IGM based on grain size.This paper presents a comprehensive wave equation analysis of piles driven in S-IGM using GRLWEAP14 software version 2010-2014 [25] by incorporating classification criteria and proposed static analysis methods for S-IGMs developed by Masud et al. [7,24].
The back-calculation of the dynamic parameters for the S-IGMs is conducted using 34 dynamic pile load test results in S-IGMs collected from four different states: Wyoming, Montana, Idaho, and North Dakota.This study introduces newly developed dynamic parameters for the S-IGMs in the improved WEAP analysis to increase the reliability of pile resistance prediction based on the proposed classification and new static analysis methods.Four different WEAP methods, namely WEAP SAD (default WEAP method), WEAP UWD (newly developed static analysis methods are used as input parameters for S-IGMs), WEAP UWR (newly developed static analysis methods and damping parameters for S-IGMs are used), and WEAP SAR (default SA method but newly developed damping parameters are used), are proposed in this study to evaluate the pile resistance predictions in S-IGMs.An additional 22 test piles are collected from the literature as an independent dataset for validating the proposed WEAP methods.The calibrated LRFD resistance factors (φ) are recommended for the proposed WEAP methods using the three probability-based methods, Factor of Safety Method (FOSM), First Order Reliability Method (FORM), and Monte Carlo Simulation (MCS), to achieve target reliabilities.
An economic impact study was conducted to assess the pile performances in shale using the GRLWEAP14 software by Islam et al. [10].However, similar studies have yet to be conducted to evaluate the economic impact analysis for driven piles in S-IGMs.This study evaluates the economic impact of WEAP SAD and WEAP UWR methods for driven piles in S-IGM to reduce construction challenges in future engineering practice.

Existing Studies on Wave Equation Analysis of Pile Driving
The numerical solution of the one-dimensional wave equation was presented by Smith [26,27] and subsequently used by Goble and Rausche [28] to develop the wave equation analysis of pile driving.They used a combination of lumped masses, springs, and dashpots to simulate the pile, hammer, driving system, and soil conditions.The Smith model is the most common soil model used in WEAP.This model accounts for the static and dynamic soil behaviors, and the total soil resistance during pile driving is the summation of static and dynamic resistances.Since its introduction, WEAP has been widely used to determine drivability, predict static and ultimate pile resistances, evaluate hammer performance, and assure the integrity of piles in soils.The compatibility of driving equipment with the pile and geomaterial conditions can be easily analyzed using WEAP [29].WEAP has proven to be a reliable method for dynamic pile analysis and an economical construction control method.
The accuracy of pile resistance prediction determined from the wave equation analysis is highly dependent on the dynamic soil parameters: quake values (Q s and Q t ) and damping factors (J s and J t ).The quake and damping parameters can be predicted from the results of Static Load Tests (SLT) and dynamic load testing utilizing the PDA with a subsequent signal-matching analysis using the CAPWAP.Several combinations of quake and damping factors can be obtained to match results from load tests and WEAP [30].Several past studies have been performed to determine the dynamic parameters summarized in Table 1.Liang and Sheng [31] provided theoretical expressions for predicting the dynamic parameters of soil in terms of pile penetration velocities, accelerations, pile sizes, and soil properties.McVay and Kuo [32] expressed the dynamic soil parameters regarding the SPT N-value and hammer energy (Er).It should be noted that these recommended dynamic parameters were developed for piles driven in soils.However, the dynamic parameters for piles driven in S-IGMs have yet to be investigated. .v p = Pile penetration acceleration; v p = Pile penetration velocity; A, B = Regression constants [27]; N s = Weighted average skin SPT N-value; N t = Weighted average toe SPT N-value; and E r = Rated energy of SPT hammer.
Using 24 steel H-piles in Wyoming, WEAP analysis was performed to predict pile resistances based on Q t values of D/120 for S-IGMs [35].Since J s for IGMs were unavailable, the authors assumed J s values of 0.16 s/m for non-cohesive geomaterials, 0.65 s/m for cohesive geomaterials, and 0.33 s/m for silt-like geomaterials.Q s = 2.5 mm and J t = 0.5 s/m were used independent of the geomaterial and pile types.The analysis yielded a mean resistance bias (ratio of CAPWAP measured pile resistance to the WEAP predicted pile resistance) of 1.0 but a relatively high Coefficient of Variation (COV) of resistance biases of 0.30.However, this study was limited to only one state and used no new static analysis methods.Similarly, Islam et al. [10] used 46 test pile data from Kansas, Montana, Iowa, and Wyoming to back-calculate the dynamic parameters for piles driven in shale.However, no similar studies have been conducted to develop a WEAP procedure for driven piles in S-IGMs.
The existing SPT N-value-based (SA) method is widely used as a geomaterial input method in WEAP.However, the SA method poses a challenge as there is no option for assigning S-IGM as a geomaterial.Although the "Other" option can be used, the program does not assign q s or q b to the S-IGMs.The q s and q b need to be manually input into WEAP; however, this creates another challenge as reliable static analysis methods for predicting these unit pile resistances in S-IGMs are unavailable in the default WEAP method.S-IGMs often exhibit SPT N-values greater than the maximum allowable input of 60 in the SA method.In those cases, other geomaterial properties, such as the friction angle (ϕ) or unconfined compressive strength (q u ), can be used to describe the geomaterial, but these properties are usually unavailable.In addition, WEAP limits the maximum q s and q b assigned by the SA method to the geomaterials (Table 2), and for a past study performed for the Wyoming Department of Transportation (WYDOT), the catalog of geomaterial properties and CAPWAP-measured resistances prepared by Adhikari [35] suggested that the CAPWAP-measured q s and q b of the cohesive soil-based IGMs ranged from 52 kPa to 172 kPa and 3687 kPa to 10,677 kPa, respectively.These average q s and q b were 78 kPa and 6608 kPa, respectively, which are higher than the maximum q s and q b recommended for clay limited by the SA method (Table 2).Applying maximum q s and q b values could lead to the underprediction of pile resistances by WEAP for piles driven in S-IGMs.Therefore, a user-defined input of the unit resistances for soil-based IGMs in WEAP is recommended.Table 3 summarizes the WEAP-recommended dynamic parameters for soil and hard rocks [35].These recommended dynamic parameters were initially developed based on pile load test data in soils [36].However, recommendations have yet to be provided for determining the dynamic parameters and unit pile resistance for S-IGMs.Thus, improved WEAP analysis methods are proposed to incorporate unit pile resistances using new static analysis methods and back-calculated dynamic parameters to predict ultimate pile resistances in S-IGMs.J s = Shaft damping; J t = Toe damping; Q s = Shaft quake; Q s = Toe quake; D = Pile diameter or width.

Pile Load Test Data
This study is conducted using 34 usable historical test pile data collected from Wyoming, Montana, Idaho, and North Dakota (Table 4).The test pile data are considered usable because they contain relevant subsurface information, PDA and CAPWAP results, and pile, hammer, and driving information.Table 4 summarizes the state, project, pile location, pile type, total pile penetration, stroke height, blow count, and hammer type of 34 test piles.The relevant geomaterial properties like SPT N-value, unit weight (γ), s u , ϕ, and q u are obtained from the respective boring logs and geotechnical reports.Among 34 test piles, 25 are steel H-piles and 9 are steel pipe piles.The total pile penetration ranges from 5.9 m to 41.6 m.The stroke height ranges from 1.68 to 3.05 m.These 34 test piles are used to conduct the bearing graph analyses, improve WEAP input methods based on the developed new static analysis methods for S-IGMs by Masud et al. [7,24], and recommend new damping parameters for both CG-IGM and FG-IGM for three different Subsurface conditions.

Four WEAP Procedures for Bearing Graph Analyses
Bearing graph analyses are performed on all useable test piles using the default WEAP SAD method and three proposed WEAP UWD, WEAP SAR, and WEAP UWR methods.The term "SA" refers to the default SPT N-value-based method, and "UW" refers to the new static analysis methods developed for the S-IGMs by Masud et al. [7,24].The "D" refers to the GRLWEAP14 software recommended dynamic parameters, and the "R" refers to newly recommended dynamic parameters determined from the back-calculation procedure for S-IGMs.These four WEAP methods are explicitly described in the following subsections, and a summary of different WEAP methods is given in Table 5.

Method
Geomaterial Input for q s and q b Quake Damping Toe quake; J s = Shaft damping; J t = Toe damping; D = Pile Diameter (mm).

WEAP SAD Method
The SA method is programmed in WEAP to define the geomaterial profile and geomaterial properties using the SPT N value (N), unit weight (γ), friction angle (ϕ), and unconfined compressive strength (q u ).Since the WEAP SA method only allows input of N ≤ 60 blows/0.3m, it is difficult to define the S-IGMs that generally have N > 60 blows/0.3m.An alternative procedure is thus proposed to define the S-IGMs in the default SA method for WEAP analysis as follows: • N ≤ 60: N and γ are used to define the geomaterial.

•
N > 60 for fine-grained soil-based IGM (FG-IGM): q u and γ are used to define the geomaterial.
WEAP assigns the default "D" dynamic parameters based on geomaterial and pile size (Table 3).The percentage shaft resistance obtained from the drivability analysis is used in the bearing graph analysis.

WEAP UWD Method
The default SA input procedure is used to input soil properties to predict the shaft resistance of overburdened soils along each test pile.For the S-IGMs, the unit resistances determined using the newly developed equations (Table 6) are manually input via the "SA input form" for individual layers.The WEAP-recommended dynamic parameters (Table 3) are assigned to the soil and S-IGMs based on the description of S-IGMs and pile type.The new percentage shaft resistance determined from a drivability analysis is used as an input in a bearing graph analysis.

WEAP SAR Method
The default SA input procedure is used to input soil properties to predict the shaft resistance of overburdened soils along each test pile.For the S-IGMs, the geomaterial input procedures explained in the WEAP SAD method are used.The significant difference between the WEAP SAR and WEAP SAD methods is that in the WEAP SAR method, newly recommended damping parameters are used, but in WEAP SAD, default WEAP damping parameters are used.The newly proposed damping parameters are manually input for each geomaterial layer using the "pile segment" and "soil segment damping/quake" options in the WEAP.

WEAP UWR Method
Since the dynamic parameters and static resistance equations used by WEAP were initially developed for soils, the existing static analysis methods might not accurately determine the q s and q b for S-IGMs.Hence, the new static analysis methods (Table 6) developed by Masud et al. [7,24], authors from the University of Wyoming "UW", are adopted to predict the q s and q b for the S-IGMs that are manually input into WEAP.The new percentage shaft resistance obtained from the drivability analysis is used in the bearing graph analysis.The recommended "R" dynamic parameters determined from the backcalculation procedure discussed in the next section are manually input into the bearing graph analysis through the "pile segment" and "soil segment damping/quake" options.For all WEAP methods, bearing graph analysis is conducted to predict the ultimate pile resistance using the field-reported hammer blow count and stoke height for each test pile shown in Table 5.

Back-Calculation Procedure for Dynamic Parameters
Pile resistance prediction from WEAP depends on selecting the quake values and damping factors.However, matching the pile resistance from CAPWAP with the predicted ultimate resistance from WEAP's bearing graph analysis can yield a combination of quake and damping parameters.Thus, back-calculation becomes an indeterminate problem that requires a reduction in unknown quake and damping parameters.A parametric study has been performed to determine the effect of the quake and damping parameters on the bearing graph [10].The study found that the J s have the most effect on the bearing graph.
On the other hand, the Q s value has little effect on the bearing graph, and thus, the default Q s value of 2.5 mm recommended in WEAP is adopted in this study for S-IGMs.Since S-IGMs are more rigid and stiffer than soil, it is reasonable to assume the WEAPrecommended Q t value of D/120 (hard and stiff soils) for the S-IGMs.For the damping factors, a J s -to-J t ratio of one is assumed to create a determinate problem as pile resistance from WEAP depends on the relative effect of both damping factors.The bearing graph is generated by changing J s and J t for the S-IGM layers while considering the WEAPrecommended values for the overburdened soil layers.The predicted resistance from WEAP is determined using the stoke height and blow count at the End of Driving (EOD) given in pile driving reports.The best back-calculated damping factors are determined by matching the predicted ultimate resistance from WEAP with the resistance determined from CAPWAP at EOD until the difference is less than 0.1%.A sample calculation is illustrated using the HP 310×79 test pile (pile ID-26 from Table 4) at the Abutment 1 location of the Pine Bluff-Beech Street bridge project in Wyoming.The CAPWAP resistance of the test pile at EOD is 1272 kN at this particular location.The pile length, penetration, hammer, stroke height, and blow counts at the EOD are 15.3 m, 13.6 m, Delmeg D 16-32, 2.07 m, and 62 b/0.3 m, respectively.Pile resistance from WEAP methods depends on quake and damping values, creating an indeterminate problem.To eliminate this indeterminacy, the shaft quake Q s value of 2.5 mm, toe quake Q t value of 2.58 mm, and J s -to-J t ratio of 1 are chosen for the bearing graph analysis.After several trials and errors, the damping factor J s = J t = 0.5 s/m was finally selected for the bearing graph analysis, and the WEAP predicted a resistance of 1270.9 kips.The resistance difference between WEAP and CAPWAP at EOD is 0.09% in this test pile.Based on this methodology, all the back-calculated damping parameters for 34 test piles are determined.

Back-Calculation Results
The back-calculations show that the back-calculated damping parameters depend upon the bearing geomaterial and the overburden subsurface condition.Three subsurface conditions, I, II, and III, are proposed, along with the recommended damping factors for driven piles in S-IGMs.A flowchart to determine the different subsurface conditions based on geotechnical and boring log information is provided in Figure 1.To eliminate this indeterminacy, the shaft quake Qs value of 2.5 mm, toe quake Qt value of 2.58 mm, and Js-to-Jt ratio of 1 are chosen for the bearing graph analysis.After several trials and errors, the damping factor Js = Jt = 0.5 s/m was finally selected for the bearing graph analysis, and the WEAP predicted a resistance of 1270.9 kips.The resistance difference between WEAP and CAPWAP at EOD is 0.09% in this test pile.Based on this methodology, all the back-calculated damping parameters for 34 test piles are determined.

Back-Calculation Results
The back-calculations show that the back-calculated damping parameters depend upon the bearing geomaterial and the overburden subsurface condition.Three subsurface conditions, I, II, and III, are proposed, along with the recommended damping factors for driven piles in S-IGMs.A flowchart to determine the different subsurface conditions based on geotechnical and boring log information is provided in Figure 1.

Subsurface Condition I (Fine-Grained Soil-Based IGM as Bearing Layer)
Out of the 34 test piles, 15 test piles are driven onto fine-grained soil-based IGMs (FG-IGMs) as the bearing layer.Subsurface condition I can be identified as having FG-IGM as the bearing IGM layer and overlying with one or more of either soil layers or FG-IGMs.No coarse-grained soil-based IGM (CG-IGM) is encountered in subsurface condition I.The back-calculated damping factors varied from 0.10 s/m to 1.23 s/m.The average back-calculated damping factor of 0.59 s/m for the FG-IGM is higher than the J t of 0.50 s/m for cohesive soils, but lower than the J s of 0.65 s/m for cohesive soils.The higher average back-calculated J t and lower average back-calculated J s indicate that the current WEAP-recommended J t and J s for cohesive soils will underpredict the shaft resistance and overpredict the end bearing of FG-IGMs.A relationship shown in Figure 2a is developed for the back-calculated damping factor based on the s u (kPa) of FG-IGM at the pile tip and slenderness ratio (i.e., the ratio of pile penetration L (m) to pile dimension D p (mm)).The fitted linear regression function for predicting Ĵt (s/m) and Ĵs (s/m) is given by Equation ( 1) with a relatively high coefficient of determination (R 2 ) of 0.74.Equation ( 1) infers that those piles with a deeper bearing FG-IGM layer will demand higher damping factors, but the damping factors will decrease with increasing pile size.Equation ( 1

Subsurface Condition II (Coarse-Grained Soil-Based IGM as Bearing Layer)
Nine test piles are driven onto coarse-grained soil-based IGMs (CG-IGM) as the bearing layer.Subsurface condition II can be identified as having CG-IGM as the bearing IGM layer and overlying with one or more of either soil layers or CG-IGMs.No FG-IGM is encountered in this subsurface condition II.The back-calculated damping factors ranged from 0.075 s/m to 0.82 s/m with the average back-calculated damping factor of 0.33 s/m.The average back-calculated damping is lower than the Jt, but higher than the Js recommended in WEAP for non-cohesive soils.This comparison indicates that the current WEAP recommendations will underpredict the end bearing and overpredict the shaft resistance in CG-IGMs. Figure 2b shows a non-linear relationship between the back-calculated damping and the combined term of pile penetration length, L (m), and corrected N, (N1)60 (blows/0.3m) below the pile tip.The nonlinear trend can be best described by the exponential relationship given by Equation ( 2) .

Subsurface Condition II (Coarse-Grained Soil-Based IGM as Bearing Layer)
Nine test piles are driven onto coarse-grained soil-based IGMs (CG-IGM) as the bearing layer.Subsurface condition II can be identified as having CG-IGM as the bearing IGM layer and overlying with one or more of either soil layers or CG-IGMs.No FG-IGM is encountered in this subsurface condition II.The back-calculated damping factors ranged from 0.075 s/m to 0.82 s/m with the average back-calculated damping factor of 0.33 s/m.The average back-calculated damping is lower than the J t , but higher than the J s recommended in WEAP for non-cohesive soils.This comparison indicates that the current WEAP recommendations will underpredict the end bearing and overpredict the shaft resistance in CG-IGMs. Figure 2b shows a non-linear relationship between the backcalculated damping and the combined term of pile penetration length, L (m), and corrected N, (N 1 ) 60 (blows/0.3m) below the pile tip.The nonlinear trend can be best described by the exponential relationship given by Equation ( 2) The predicted Ĵt (s/m) and Ĵs (s/m) from Equation (2) show that the damping factor decreases with increasing (N 1 ) 60 of CG-IGMs below the pile tip.Because of limited data points with 8.5 m ≥ L ≤ 14.8 m, 55 blows/0.3m ≥ (N 1 ) 60 ≤ 151 blows/0.3 m.The pile types used in the development of the prediction equation in Equation ( 2) include HP 310×79, HP 360×108, and 406-mm open-ended pipe piles (OEP).Future studies are recommended to improve the prediction equation for damping factors by including additional data with longer pile lengths and variable pile types and sizes.

Subsurface Condition III (Either Soil or Soil-Based IGM as Bearing Layer)
Out of the 34 test piles, 10 fall under subsurface condition III, with the bearing layer identified as either a fine-grained soil-based IGM (FG-IGM), coarse-grained soil-based IGM (CG-IGM), or soil.For the test piles with the FG-IGM or CG-IGM as the bearing layer, this subsurface condition differs from the subsurface conditions I and II due to the presence of both FG-IGM and CG-IGMs along the pile shaft.For piles with soil as the bearing layer, subsurface condition III is identified when FG-IGM, CG-IGM, or both are present as overburden layers.Figure 1 clearly illustrates the identification of subsurface condition III from subsurface conditions I and II.Since damping parameters are found to depend on the subsurface conditions, the application of Equations ( 1) and ( 2) is found to be impractical for this condition.The average back-calculated damping factors of 0.33 s/m and 0.59 s/m are recommended for the CG-IGM and FG-IGM layers, respectively, to facilitate the proposed WEAP methods.In addition, the recommended dynamic parameters from the back-calculation study for S-IGMs, along with the default dynamic parameters for soils, are summarized in Table 7 to facilitate the implementation of improved WEAP methods summarized in Table 5.
Table 7. Quake and damping parameters for subsurface conditions I, II, and III.

Determination of Pile Resistances from Bearing Graph Analysis
Bearing graph analyses are conducted to determine the ultimate pile resistances considering all four WEAP methods and three subsurface conditions.The bearing graph analysis procedure for different WEAP methods is illustrated in a flowchart in Figure 3.Each predicted ultimate pile resistance from each WEAP method is compared with the measured pile resistance from CAPWAP at EOD in terms of a resistance bias (λ) where λ > 1 indicates that WEAP underpredicts the pile resistance and λ < 1 indicates that WEAP overpredicts the pile resistance.The comparison of the predicted ultimate resistances from each WEAP method and resistances determined from CAPWAP for 34 test piles is shown in Figure 4. Figure 4 shows that WEAP SAD and WEAP UWD, with a mean resistance bias of 1.14, underpredict the pile resistances, on average, by approximately 12%, while the WEAP UWR and WEAP SAR, with a mean bias of 1.02, provide better pile resistance prediction as the mean bias is close to 1.Among the four WEAP methods, WEAP UWR provides better pile resistance prediction as the mean bias is close to unity and the coefficient of variation (COV) of bias is 0.16, the lowest among the four methods.

Validation of Improved WEAP Methods
Twenty-two independent pile data obtained from the pile load test program in Illinois by Long and Anderson [37], as summarized in Table 8, are utilized to validate the improved WEAP methods for predicting pile resistances in subsurface conditions for S-IGMs.Of the 22 test piles, 9 piles are steel H-piles and the remaining 13 test piles are 356-mm diameter close-ended pipe piles.For the bearing IGM layers, 8 test piles are driven in FG-IGM, and 14 test piles are driven in CG-IGM.The classification for CG-IGMs is conducted based on N ≥ 58 blows/0.3m, and the classification of FG-IGMs is based on the reported s u values ≥ 129 kPa following the criteria conducted by Masud et al. [6,22].Bearing graph analyses are performed using the reported stoke heights and blow counts at the EOD for 22 independent test pile data.The predicted ultimate pile resistances from four different WEAP methods are compared with those from CAPWAP in Figure 5.The mean and COV of resistance biases for the four WEAP methods are included in Figure 5 for comparative study.The highest COV of 0.25 is obtained from the WEAP SAD and WEAP UWD methods.Having the mean resistance bias of 1.04, the pile resistances are underpredicted by the WEAP SAD and WEAP UWD, on average, by only approximately 4%.However, the WEAP UWR and WEAP SAR improve the pile resistance predictions by reducing the underprediction of pile resistances by around 5% (from 1.04 to 0.99) and the COV by 20% (from 0.25 to 0.20).
Considering the different environmental conditions outside the study's geographical scope, WEAP UWR is the best method of four for pile resistance prediction with a mean bias closer to 1 and the lowest COV of 0.20.
This WEAP UWR method incorporates the newly proposed damping parameters and static analysis methods and provides more consistent resistance estimation compared to other methods (Figure 5).The same observations of the mean close to unity and lower COV are obtained for WEAP UWR methods based on 34 training test pile data (Figure 4).This WEAP UWR method incorporates the newly proposed damping parameters and static analysis methods and provides more consistent resistance estimation compared to other methods (Figure 5).The same observations of the mean close to unity and lower COV are obtained for WEAP UWR methods based on 34 training test pile data (Figure 4).

LRFD Resistance Factors
The resistance factors (φ) for WEAP, as recommended by the AASHTO [38], are developed initially for piles driven in soils.To ascertain the reliability of WEAP on pile resistance predictions, LRFD resistance factors are calibrated to reflect the uncertainty associated with the resistance prediction of driven piles in S-IGMs.Improved pile design reli-

LRFD Resistance Factors
The resistance factors (φ) for WEAP, as recommended by the AASHTO [38], are developed initially for piles driven in soils.To ascertain the reliability of WEAP on pile resistance predictions, LRFD resistance factors are calibrated to reflect the uncertainty associated with the resistance prediction of driven piles in S-IGMs.Improved pile design reliability is obtained by limiting the database to a particular geologic region and corresponding construction practices [38].LRFD resistance factor values are calibrated using FOSM, FORM, and MCS.
The efficiency of the different WEAP methods Is assessed using efficiency factors ( φ (x) ), which are established at reliability indexes (β T ) of 2.33 and 3.0 for a redundant pile group and a non-redundant pile group, respectively [39].The calibration of resistance factors using FOSM, FORM, and MCS required the mean, standard deviation, and COV of resistance biases.The details of the calibration procedures are discussed in the literature [39,40].
The φ values are calibrated using the 34 training test piles and 22 independent test piles used in the validation.The Shapiro-Wilk (SW) [41] test and Anderson-Darling (AD) [42] test are performed to assess whether the resistance biases meet the assumptions of Normality or Log-Normality.A significance level of p value = 0.05 is used for these normality tests.The results of the normality tests are summarized in Table 9.For the WEAP SAD and WEAP UWD, no violation of the Normal and Log-Normal distribution assumptions is detected because the P-values are greater than 0.05.However, for WEAP UWR and WEAP SAR, both the Normal and the Log-Normal distribution assumptions are rejected by SW and AD tests as their P-values are less than 0.05.Test piles are identified as outliers if the difference between λ and exceeds twice the standard deviation.These outliers are then excluded to improve the quality of the calibrated φ [39].The normality test results after removing the outliers (Table 9) showed that both the Normal and Log-Normal distribution assumptions were not violated.Thus, the Log-Normal distribution is selected in the LRFD resistance factor calibration based on the higher Log-likelihood value.Table 10 summarizes the resistance and efficiency factors for the four WEAP methods.The φ values of WEAP SAD and WEAP UWD vary from 0.70 to 0.82, while higher φ values ranging from 0.72 to 0.94 are calibrated for WEAP UWR and WEAP SAR at β T = 2.33.Compared with φ = 0.50 recommended in the AASHTO [38] for soils, the calibrated φ values using FOSM at β T = 2.33 for piles in S-IGMs are 40%, 40%, 44%, and 50% higher for WEAP SAD, WEAP UWD, WEAP SAR, and WEAP UWR, respectively.The obtained resistance factors are also higher than the φ = 0.40 recommended by Paikowsky et al. [40].These higher calibrated φ values are attributed to mean bias closer to one and lower COV This overestimation of the number of piles could create cost overrun during the construction phase due to excess steel piles used.Alternatively, if the WEAP UWR is used, there is a slight overestimation of the pile resistance, leading to a slight underestimation of several WEAP SAD methods; the WEAP UWR method seems to be more efficient as the average difference in steel weight per unit load is closer to zero.Following these steps, the economic impact assessment is conducted only for the first 32 test piles shown in Table 5 because pile numbers 33 and 34 do not have the reported load demand.The validation dataset for the piles driven in S-IGMs is not used as the structural load demand is not reported.Figure 7 shows the difference in steel weight for 32 test piles for both the WEAP SAD and WEAP UWR methods.The average differences in steel weight per unit load for the WEAP SAD and WEAP UWR are −0.54kg/kN and 0.06 kg/kN, respectively.If the WEAP SAD method had been used as the construction control method instead of CAPWAP, it would have highly underestimated pile resistance and resulted in overestimating the number of piles.This overestimation of the number of piles could create cost overrun during the construction phase due to excess steel piles used.Alternatively, if the WEAP UWR is used, there is a slight overestimation of the pile resistance, leading to a slight underestimation of several WEAP SAD methods; the WEAP UWR method seems to be more efficient as the average difference in steel weight per unit load is closer to zero.

Summary and Conclusions
Thirty-four test piles driven into S-IGMs from four states are used in this study to develop a new procedure for WEAP analysis to improve the prediction of pile resistances in S-IGMs.The back-calculation of the dynamic parameters for the S-IGMs is performed by matching the predicted ultimate pile resistance from WEAP with the resistance determined from CAPWAP at the EOD.This study focuses on improved WEAP analysis for steel piles driven into soil-based IGMs only.However, the same methodologies can be adopted for other pile materials and types to develop new damping parameters and WEAP methods.The economic impact assessment conducted in this study is limited to performance analysis in terms of the number of differences in steel weight.However, different pile types, maintenance, and long-term performance can be considered in future studies.Subsurface variability and uncertainty are also important factors that are not considered in this study but are suggested for future studies.The study yields the following conclusions:

•
Quake values of 2.5 mm and D/120 for Qs and Q t , respectively, are adequate for the S-IGMs.Smith damping factors are found to depend on the pile driving and subsurface conditions.Single values of Smith damping factors are inadequate for different pile and driving conditions.Hence, new Smith damping factors are proposed for three different subsurface conditions.

•
Using 34 training pile test data and 22 independent test pile data, it is found that WEAP UWR is the most efficient as it provides a mean resistance bias of 1.02 closer to 1 and the lowest COV of 0.18.• A φ value of 0.75 for WEAP UWR calibrated based on FOSM at β T = 2.33 for piles driven into S-IGMs is higher than the value of φ of 0.5 recommended in AASHTO [38] for piles in soils.Compared with the default WEAP SAD, the selected proposed WEAP UWR method, on average, reduces the underprediction of pile resistances by 6% and improves the reliability with a 43% reduction in the coefficient of variation (COV) for β T = 2.33.• The economic impact assessment reveals that the average differences in steel weight per unit load for the WEAP SAD and WEAP UWR are −0.54kg/kN and 0.06 kg/kN, respectively.Compared with the WEAP SAD method, the WEAP UWR method seems to be more efficient as the average difference in steel weight per unit load is closer to zero, which will reduce construction challenges.

Figure 1 .
Figure 1.Flowchart to classify the subsurface conditions I, II, and III.

5. 2 .
Subsurface Condition I (Fine-Grained Soil-Based IGM as Bearing Layer) Out of the 34 test piles, 15 test piles are driven onto fine-grained soil-based IGMs (FG-IGMs) as the bearing layer.Subsurface condition I can be identified as having FG-IGM as the bearing IGM layer and overlying with one or more of either soil layers or FG-IGMs.No coarse-grained soil-based IGM (CG-IGM) is encountered in subsurface condition I.The back-calculated damping factors varied from 0.10 s/m to 1.23 s/m.The average backcalculated damping factor of 0.59 s/m for the FG-IGM is higher than the Jt of 0.50 s/m for

Figure 1 .
Figure 1.Flowchart to classify the subsurface conditions I, II, and III.

1 )Figure 2 .
Figure 2. (a) Relationship between back-calculated damping factor and the product of undrained shear strength and slenderness ratio for FG-IGMs and (b) relationship between damping factor and the ratio of pile penetration to corrected SPT N-value for CG-IGMs.

Figure 2 .
Figure 2. (a) Relationship between back-calculated damping factor and the product of undrained shear strength and slenderness ratio for FG-IGMs and (b) relationship between damping factor and the ratio of pile penetration to corrected SPT N-value for CG-IGMs.

Figure 3 .
Figure 3. Flowchart comparing the four different WEAP methods.Figure 3. Flowchart comparing the four different WEAP methods.

Figure 3 .
Figure 3. Flowchart comparing the four different WEAP methods.

Figure 4 .
Figure 4. Comparison of measured pile resistances from CAPWAP and predicted pile resistance from WEAP SAD, WEAP UWD, WEAP SAR, and WEAP UWR for the 34-training pile dataset.

Figure 4 .
Figure 4. Comparison of measured pile resistances from CAPWAP and predicted pile resistance from WEAP SAD, WEAP UWD, WEAP SAR, and WEAP UWR for the 34-training pile dataset.

Figure 5 .
Figure 5.Comparison of measured pile resistances from CAPWAP and predicted pile resistances from WEAP SAD, WEAP UWD, WEAP SAR, and WEAP UWR for the 22 independent pile data.

Figure 5 .
Figure 5.Comparison of measured pile resistances from CAPWAP and predicted pile resistances from WEAP SAD, WEAP UWD, WEAP SAR, and WEAP UWR for the 22 independent pile data.

Figure 6 .
Figure 6.A flow chart showing the methodology for economic impact assessment.

Figure 7 .
Figure 7. Plot of the differences in steel weight based on the WEAP SAD and WEAP UWR methods.

Figure 6 .
Figure 6.A flow chart showing the methodology for economic impact assessment.

Figure 7 .
Figure 7. Plot of the differences in steel weight based on the WEAP SAD and WEAP UWR methods.

Figure 7 .
Figure 7. Plot of the differences in steel weight based on the WEAP SAD and WEAP UWR methods.

Table 1 .
Summary of dynamic parameters for soils from previous studies.

Table 2 .
[35]mum unit shaft resistance and unit end bearing used in the SA method[35].

Table 4 .
Summary of historical test pile data and relevant End of Driving (EOD) information.

Table 5 .
Summary of WEAP methods.

Table 8 .
Summary of independent test pile data and relevant End of Driving (EOD) information for validation.

Table 9 .
Results from the Shapiro-Wilk and Anderson-Darling tests.