TAMU-POST: An analysis tool for vehicle impact on in-line pile group

Abstract The tragic consequences of vehicles running into infrastructures have raised the need for perimeter protection. One common perimeter barrier is a set of piles or posts in an in-line geometry as an efficient way to contain or redirect errant vehicles. To date, the design of such barriers relies mostly on performing full-scale crash tests. These crash tests are expensive, and it is often practical to run such tests. In this paper, a general yet simple analysis-design model called TAMU-POST was developed to predict the response of a group of in-line piles impacted by a vehicle. TAMU-POST is based on the finite difference solution to the governing differential equation for a beam supported by piles. The piles are represented by single degrees of freedom consisting of a dashpot, a lumped mass, a spring, and a slider. A large number of computer simulations using a non-linear finite element program LS-DYNA as well as the data obtained from two full-scale crash tests were used to calibrate the proposed model. The design quantities are the barrier deflection, the vehicle dynamic penetration defined as the maximum vehicle intrusion into the barrier, as well as other parameters including the bending moment in the piles and in the beams. A Monte Carlo Simulation analysis was conducted using TAMU-POST to evaluate the probability of failure of a group of in-line piles under a given vehicular impact when considering the inherent uncertainties associated with the input parameters and the model coefficients.


PUBLIC INTEREST STATEMENT
Infrastructures are constantly exposed to extreme loading conditions that would cause tremendous damage, injuries, and fatalities. Advancements in computational modeling and design technologies have provided tools to protect our infrastructure against hazards such as vehicle crashes, explosions, security threats, etc. Cost of real-scale testing for certifying infrastructure designs is drastically high. Predictive models and computational techniques are needed to understand the behavior of materials/ components under extreme loads for a resilient and sustainable design. This research focuses on one of the most common extreme hazards on infrastructures, vehicle impacts, and design of various crashworthy geo-structural systems through real-scale, component tests, and advanced numerical simulations. First, the response of soil and soil-structure interaction under extreme loads is characterized through laboratory and in-situ tests, pendulum tests, bogie tests, and real-scale crash tests. Computational simulations are developed to replicate these highly nonlinear problems and then validated against the tests. Numerical results coupled with experiments were used to develop a simplified macro-element model to predict pile deflection when subjected to impacts.

Asadollahi Pajouh Mojdeh
Briaud Jean-Louis deflection, the vehicle dynamic penetration defined as the maximum vehicle intrusion into the barrier, as well as other parameters including the bending moment in the piles and in the beams. A Monte Carlo Simulation analysis was conducted using TAMU-POST to evaluate the probability of failure of a group of in-line piles under a given vehicular impact when considering the inherent uncertainties associated with the input parameters and the model coefficients.

Introduction
Proper design of anti-ram perimeter barriers becomes crucial when their failure most likely involves tremendous loss of lives and destruction. Groups of in-line piles have been used for such barriers yet little is known about their behavior under lateral impact. Numerous studies have been focused on lateral response of piles under cyclic loading, earthquake loading, and other types of vibration but much less have been done on impact loading. Early work on dynamic soil-pile interaction was first reported by Wolf (1980). Then, a number of analytical approaches were suggested to evaluate soil-pile interaction under cyclic loading (Brown, Morrison, & Resse, 1988;El Naggar & Novak, 1996;Nogami, Konagai, & Otani, 1991;Wu & Fin, 1997). These studies made a number of simplifying assumptions. Other studies followed with more advanced input parameters and analytical models (Dobry & Gazetas, 1988;Wolf & Somaini, 1986).
Most of the impact loading studies have been conducted in the field of roadside safety to develop guardrail barrier systems. These barriers consist of widely spaced piles, referred to as posts in the transportation field, tied together with a rail aimed at safely redirecting errant vehicles. In this work, various design configurations have been studied including post embedment depth, post spacing, and post materials (either wood or steel) under different impact conditions. However, most of the work has consisted of crash testing at full scale without much attention paid to the soil behavior including the development of a simple analysis program as is proposed in this article. This work is described in references such as Abu-Odeh et al. (2014), Abu-Odeh, McCaskey, Bligh, Menges, andKuhn (2015), Alberson, Menges, and Schoeneman (2001), Bligh, Abu-Odeh, and Menges (2011), Bligh and Menges (1997), Faller et al. (2004. These studies indicated that the guardrail performance in redirecting impacting vehicles largely depends mostly on the soil strength and the pile stiffness. The work described in this article was motivated by the lack and thus the need for a simple analytical approach to design soil-embedded piles subjected to impact loading. The current approach is to perform a full-scale crash test with the right vehicle hitting the right barrier. However, the approach is costly. Numerical simulations are useful to augment the value of the fullscale crash tests but are computationally expensive and require to be validated against full-scale crash tests before they can be used with improved confidence. In previous studies by the authors, two full-scale crash tests, test nos. PU60 and M50 were performed on in-line pile groups (Asadollahi Pajouh, Briaud, Lim, & Mirdamadi, 2017). In test no. PU60, a 2300-kg pickup truck impacted head on an in-line group of four HSS6x6x1/2 piles embedded 2 m in stiff clay, spaced 5.2 m apart, tied together with two steel HSS8x8x1/2 beams, and impacted at 96.55 km/h. In test no. M50, a 6800-kg medium-duty truck impacted head on an in-line group of eight steel W14x90 piles embedded 3 min loose sand, spaced 5.2 m and tied together with two steel HSS8x8x1/2 beams, and impacted at 80 km/h. In both cases, the in-line group of piles successfully stopped the test vehicle (Figure 1). A number of detailed LS-DYNA simulations were performed to replicate the crash tests (Hallquist, 2007). These experiments and the simulation results were then used to develop TAMU-POST. TAMU-POST is a simple code to predict the behavior of a group of in-line piles embedded in the soil, connected by a beam, and subjected to a vehicle impact.

Model theory
The development of the analytical model had three goals: to be as fundamentally correct as possible, to be simple enough to use without advanced computing, and to predict reasonably well the actual deflection of the barrier. Indeed, the primary parameter to be predicted was the dynamic penetration of the vehicle into the barrier. The dynamic penetration is defined in ASTM F2656-07, "Standard Test Method for Vehicle Crash Testing of Perimeter Barriers" (ASTM Standard F2656, 2007, as the distance between the leading edge of the flatbed of the vehicle at maximum displacement during the impact and the undeflected position of the barrier (Figure 2). According to the standard ASTM F2656-07, the test designation P1 requires that the dynamic penetration is less than 1 m for the barriers protecting sensitive buildings.
In early studies of piles subjected to static lateral loads, the beam-on-springs foundation approach was used and assumed that the pile was an Euler-Bernoulli beam on a series of springs (Broms, 1964a(Broms, , 1964bMatlock & Reese, 1960). The governing differential equation for the beam deflection was developed and the finite difference method was used to solve this equation. The beam-on-springs foundation approach was then modified and extended to the p-y curve method (Matlock, 1970;O'Neill, Reese, & Cox, 1990;Reese & Welch, 1975). For the dynamic behavior, a number of simplified models were used to represent the soil. One of the most common models was proposed by Smith (1960), as shown in Figure 3(a,b). This model consisted of a lumped mass, a spring, a dash pot, and a slider. In the current study, this lumped-parameter model (Smith, 1960) was adopted to simulate the soil-pile interaction (Figure 3(c)). In this Single Degree of Freedom (SDOF) model, the mass represents  the mass of the beam tributary to one pile plus the mas of the engaged surrounding soil. The damper and the spring represent the damping coefficient and the static stiffness, respectively. The soil surrounding the pile fails when the soil deflection reaches the soil-failure deflection. The failure load (P ult ) and the associated failure deflection (y f ) are discussed later in the paper. When the soil fails, the spring stiffness becomes zero.
In order to derive the governing differential equation for the barrier including the piles and the beam connecting the piles together, the following simplification was chosen. The beam was modeled as an elastic beam supported at discrete intervals by soil-embedded piles; these soilembedded piles were each represented by an SDOF element as shown in Figure 3(c). The complete beam model was thus as shown in Figure 4. With this assumption, the Euler-Bernoulli beam theory equation relating the beam deflection (y) to the applied lateral load q(x, t) on the beam is expressed as Equation (1): where E b is the elastic modulus of the beam, I b is the second moment of area (moment of inertia) of the beam with respect to the axis perpendicular to the applied load, y, and x are the beam deflection and the distance along the beam, respectively, q(x, t) is the transverse load applied at a distance x and at a time t. The load itself was not constant and varied with time; therefore, the equation controlling the behavior of the SDOF mass-spring-dashpot model with an applied load varying with the time t was:

M
@ 2 y @t 2 þ C @y @t þ Ky ¼ Àqðx; tÞ where M is the mass, K the stiffness of the system and C the dashpot constant. Assuming that the bending stiffness E b I b of the beam is constant, Equations (1) and (2) can be combined to give the governing differential equation for the system.
Equation (3) is solved for the deflection of the beam (y) using the central finite difference method together with the corresponding initial values and boundary conditions. The beam was discretized into a number of equal length elements h. For each node i at time j, the deflection of the beam y is computed explicitly using Equation (4).
In order to facilitate the numerical solution, the following parameters m, e, and c were introduced.
Equation (4) is written in matrix form as Equation (6).
: : : : Àe 0 0 : : : : : : : : : Solving these equations for all the node deflections requires solving a simple diagonal matrix equation without the need to inverse the matrices. As such, the beam deflection y at node i at any time step was computed based on the known beam deflections of that node and the adjacent nodes from the previous time steps as follows: ðm þ cÞy i; jþ1 ¼ ðc À mÞy i; jÀ1 À ey iþ2; j þ 4ey iþ1; j À ð6e À 2m þ kÞy i; j þ 4ey iÀ1; j À ey iÀ2; j (7)

Initial conditions
The node under impact displaces y i during the first time step Δt and the displacement is given by Equation (8).
where V v is the velocity of the approaching vehicle at the time of impact. This results in an extremely large acceleration at that node. To avoid convergence issues associated with this problem, a very small-time step of 10 −5 sec was used in the explicit solution.

Boundary conditions
Boundary conditions express that the shear force proportional to @ 3 y=@x 3 À Á and the bending moment proportional to @ 2 y=@x 2 À Á at both ends of the beam are zero. Therefore, the following conditions should be met: 3. TAMU-POST model formulation The theory described above was applied to a beam supported by piles subjected to a vehicle impact. To optimize the combination of reasonable accuracy and required computational effort, the beam was discretized into 15 nodes and therefore 14 elements ( Figure 5). This combination was convenient for various in-line group configurations (three to eight piles) and allowed the user to choose the impacted node either at a pile location or at mid-span between two piles. For node i at time j the beam deflection y was computed explicitly for a time increment equal to Δt, using the finite difference expression of the differential equation (Equation (4)). A summary of the model input parameters is listed in Table 1. Note that although the beam was assumed to behave elastically, the model parameters were identified using data obtained from fully nonlinear dynamic numerical analyses (using LS-DYNA simulations) and full-scale crash tests.

Numerical simulations
Full-scale crash tests provide the best verification for a given system. Numerical simulations are less expensive, typically less time-consuming, while reasonably reliable once they are calibrated  & Conover, 1979) was used to select the best cases for the numerical work. The stratification and random sampling within the strata incorporated in SLHS guarantee a better convergence and ensure that the selected random samples best represent the real variability.
To determine the strata, the most influential input variables were identified based on experience and engineering judgment. They were the vehicle mass, the vehicle velocity, the soil strength, the pile embedment depth, and the pile spacing each having a range of values (Table 2). With respect to soil strength, a statistical sampling was performed for both clay and sand soil types. With respect to the barrier configuration, two barriers of four in-line piles and eight in-line piles were considered as they had been previously crash tested. Tables 3 and 4 give a list of the soil strength categories and parameters for clay and sand, respectively. These tables lead to a total number of cases equal to 90 simulations. Note that for modeling the sand behavior in LS-DYNA, a small nonzero value for cohesion was used to avoid numerical instabilities.

Lateral Static Stiffness
The lateral soil-pile stiffness (K) in the TAMU-POST model (Equation (3)) is the ratio of the load applied at the top of the pile (groundline) and the associated deflection. The behavior of piles subjected to static lateral loading has received a lot of attention since the work of Briaud and Ballouz (1996), Broms (1964aBroms ( , 1964b, Matlock (1970). The SALLOP method (Simple Approach for Lateral Load on Piles) developed by Briaud (1997) was found to be a suitable method for this investigation (Briaud, 1997). This method is based on the close analogy between the cylindrical expansion of the pressuremeter and the lateral movement of the pile. The SALLOP method makes a distinction between a long-flexible pile and a short-rigid one. This distinction is based on the transfer length defined as: where E p is the elastic modulus of the pile material, I p the moment of inertia of the pile, K s the soil stiffness defined as the ratio of the line load P along the pile at a depth z and the deflection y of the pile at that same depth. In this method, the soil lateral static stiffness K s is given by: where E o is the first load modulus from the pressuremeter curve. The pile is considered long and flexible if the pile length is larger than 3l 0 , and considered short and rigid when the pile length is less than l 0 . For piles with an intermediate length (l 0 < L < 3l 0 ), a linear interpolation between the long-flexible case and the short-rigid case is used.
In the case of long-flexible piles subjected to a groundline horizontal load (H o ) and a groundline moment (M o ) (Figure 6), the shear force V at a depth z and the depth to zero-shear D v (i.e., depth where the maximum bending moment occurs) is computed using Equations (12) and (13), respectively.
Similarly, in the case of short-rigid piles, the shear force V at a depth z and the depth to zero-shear D v can be obtained as follows: For an impact load, the horizontal impact load (H o ) is usually applied to the pile at a certain height above the ground surface ( Figure 6). This load creates a moment at the groundline M 0 ¼ H 0 h impact where h imapct is the impact height.
According to the SALLOP method, the pile lateral deflection at the groundline for long flexible and short-rigid piles can be obtained using Equations (16) and (17), respectively.
Using M o equal to the product of H o by the impact height h impact in Equations (16) and (17), the static stiffness of the pile-soil system K, defined as the ratio of the load to the associated deflection is determined using Equations (18) and (19) for long-flexible and short-rigid piles, respectively:

Ultimate Lateral Capacity
In the SDOF model for each pile, a slider with a failure load of P ult was incorporated in line with the spring to account for the failure of the single pile-soil failure (Figure 3(c)). If the soil deformation exceedsy f (i.e., pile deflection at the soil failure), the load becomes equal to P ult and the spring stiffness K is excluded from the computation. In the SALLOP method, the ultimate lateral capacity of the pile, from the soil point of view, is defined as the load corresponding to a pile deflection equal to one-tenth of the pile diameter or width. At that deflection the SALLOP method gives the ultimate lateral capacity P ult as: where P L is the soil PMT limit pressure close to the surface, B the pile diameter or width, and D v the zero-shear depth given by Equations (13) and (15). If the pressuremeter test data (P L ) are not available, other in-situ or laboratory tests results, including the Standard Penetration Test (SPT) blow counts, the undrained shear strength (S u ), and the cone penetration test (CPT) point resistance (q c ) can be used to estimate P L through, for example, the empirical correlations recommended by Briaud (2013). Finally, the deflection at failurey f for an elastic perfectly plastic soil-pile model is given by: Figure 6. Horizontal load on a single pile.

Mass
The mass (M) attributed to each SDOF includes three terms: (1) the mass of pile (M p ) contributing to the lateral resistance by inertia; this mass is taken as the mass of the pile segment from the ground surface to the zero-shear depth (D v ); (2) the mass of beam (M b ) tributary to the node; this mass is taken as the mass of beam for a length of h (element length) separating two piles; and (3) the mass of soil (M s ) associated with the impact; this mass is considered to be related to the mass of the soil wedge in front of then pile mobilized during the impact (Figure 7) and estimated as: Where ρ is the soil density and B is the pile width or diameter.

Damping
In complex crash analyses, the damping term is difficult to quantify. The damping term in Equation (3) was used to approximate the unknown nonlinear energy dissipation in the pile-soil system. Although the phenomenon of soil damping under dynamic loading has been discussed in numerous studies (Wolf & Somaini, 1986;Roesset, Whitman, & Dobry, 1973;Rainer, 1974), there is no well-developed approach to quantify soil damping under impact loading. It is believed that damping in vehicle-barrier crashes originates mainly from energy dissipation either through wave propagation in the soil medium known as radiation damping, or through viscous damping associated with the rate effects in soils. Regardless of what type of damping governs the behavior, the damping value C in the simplified model of Equation (3) can be back calculated through comparison between the data from full-scale crash tests and comparison with numerical simulations on one hand and the TAMU-POST model output.
In the literature, the dynamic damping of a foundation has been correlated to such parameters as foundation dimension, foundation stiffness, and soil shear wave velocity (Wolf & Somaini, 1986). In that approach, damping increases with foundation dimension and foundation stiffness but is inversely proportional to the soil shear wave velocity. Using this established approach, the damping (C) was correlated to the pile width (B), the soil-pile stiffness (K), and the soil shear wave velocity (V SW ). The shear wave velocity (V SW ) was estimated as: where E o is the PMT first load modulus, ρ is the soil density, υ is Poisson's ratio, and α is a fitting factor. The proposed equation for the damping coefficient C is: The dimensionless factor α was obtained by optimizing the match between the TAMU-POST predictions with the crash test measurements and the LS-DYNA numerical simulations. A reasonable correlation between α and the pressuremeter limit pressure was found ( Figure 8) and Equation (25).

Energy absorbed in crushing of the vehicle
The severity of the vehicle-barrier collision depends on the impact speed, vehicle stiffness, and barrier stiffness. The energy absorbed in the crushing deformation of the vehicle is not accounted for in the TAMU-POST model. To account for the vehicle crushing, a reduction factor (κ) was applied to the vehicle velocity at impact in proportion with the energy absorbed in deforming the vehicle.
where V v and V 0 denote the reduced and initial vehicle velocity at impact, respectively. The TAMU-POST model utilizes V v (Equation (8)). The velocity reduction factor κ was determined based on the energy balance in a large number of numerical simulations. If the kinetic energy of the vehicle is E K , the energy expanded to crush the vehicle E C , and the energy expanded to deform the barrier E B , the energy balance is: The barrier energy (E B ) can be obtained from the numerical simulations and the reduction factor κcan be obtained as follows: The simulations included three classes of vehicles (medium-duty truck, pickup truck, and small car) as well as different pile group layouts. The values of κ based on the energy results of all the numeral simulations are presented in Figure 9. Since a higher value of κ is conservative in this case, it appears reasonable to select a crushing factor equal to 0.8 for reducing the impact velocity of the vehicle, thereby including the energy lost in crushing.

Crushing length
TAMU-POST predicts, among other results, the maximum deflection of the barrier D max . However, the critical measure of interest is the dynamic penetration DP. DP is defined as the distance between the initial position of the row of piles and the front edge of the vehicle flatbed at the time of maximum displacement of the barrier (Figure 2). DP can be expressed as: where DP is the dynamic penetration, D max the maximum deflection of the barrier, and L Crushed is the distance between the front bumper of the crushed vehicle at the time of maximum displacement of the barrier and the front edge of the vehicle flatbed. The sign of DP would be negative if the edge of the truck flatbed does not pass the initial position of the piles.
To determine DP, additional analyses were necessary because while D max is predicted by TAMU-POST, the crushing length L Crushed is not. Conservative estimates of L Crushed were obtained by simulating crashing vehicles against a rigid wall for each class of vehicle ( Figure 10). For these simulations, the  finite element models of a Toyota Yaris passenger car, a Chevrolet pickup truck, developed by the National Crash Analysis Center, (NCAC, 2014), and a medium-duty truck developed by the Texas A&M Transportation Institute (TTI, 2009) were used. Note that these values give an upper bound of L Crushed as the impact wall was rigid while the in-line group of piles is more flexible. Such upper bounds lead to larger dynamic penetration values, and therefore are conservative. Table 5 gives the recommended L Crushed values obtained from the simulations for a given vehicle mass (M v ) and velocity (V v ).

Evaluation of the model
The model precision was evaluated by comparing the TAMU-POST predictions to the results of a large number of computer LS-DYNA simulations covering a wide range of soil and pile properties, pile spacing, impact level (i.e., vehicle mass and velocity), and pile layouts. In addition, TAMU-POST predictions were compared to the results of two full-scale crash tests: tests nos. PU60 and M50 (Asadollahi Pajouh, Briaud, Alberson, Arrington, & Mirdamadi, 2015;Asadollahi Pajouh et al., 2014). In test no. PU60, a 2300-kg pickup truck impacted head-on at 100 km/h an in-line group of four steel piles (HSS6x6x1/2) embedded 2 m in stiff clay with a 5.2 m spacing tied together with two beams (HSS8x8x1/2). In test no. M50, a 6800-kg medium-duty truck impacted head-on at 80 km/h an in-line pile group of eight steel W14x90 embedded 3 m in loose sand and spaced at 5.2 m and tied together with two HSS8x8x1/2 beams. Figure 11 shows a comparison of the TAMU-POST

Parametric study
In order to assess the influence of the main input parameters, a parametric study was performed using TAMU-POST. These parameters included the ratio of pile spacing to pile width or diameter (s/ d), pile embedment depth (L), vehicle mass and velocity (M v and V v ), soil strength (P L ), and number of piles. This parametric study involved two different impact cases: a 2300-kg pickup truck impacting a four-pile (HSS6x6x1/2) barrier and a 6800-kg medium-duty truck impacting an eightpile (W14x90) barrier. For both barriers, the piles were linked by an HSS8x8x1/2 beam.

Effect of pile spacing and pile embedment depth
In the case of the eight-pile group impacted by the medium-duty truck, the ratio of pile spacing to pile width (s/d) was varied between 10 and 20 while the pile embedment depth ranged from 1 to 2.5 m. As shown in Figure 12, the barrier deflection increases significantly with an increase in pile spacing and a decrease in pile embedment. This impact becomes less predominant as s/d increases.

Effect of mass and velocity of the vehicle
The vehicle mass was varied from 2,000 to 7,000 kg, and the vehicle velocity was varied from 65 km/h to 100 km/h. As shown in Figure 13, the barrier deflection increases nearly linearly with an increase in vehicle mass and velocity.

Effect of soil strength
The maximum deflection of the two different barriers (i.e., four-pile and eight-pile) was studied for different soil strengths and pile embedment depths. The pressuremeter limit pressure (P L ) was varied from 0.3 MPa (very soft soil) to 2.5 MPa (very hard soil). The pile embedment depth was varied from 1 to 2 m. As shown in Figure 14, an increase in pile embedment depth significantly decreases the maximum deflection up to an optimum pile depth. This is particularly true for low strength soils. An increase in soil strength also decreases the maximum deflection. This decrease is not as pronounced for high strength soils.

Effect of number of piles
In the case of the pile group impacted by a pickup truck traveling at speed of 100 km/h, the number of piles was varied from 3 to 8 while keeping the other factors constant including the pile spacing. As shown in Figure 15, an increase in the number of piles leads to a decrease in the maximum deflection of the barrier up to 6 piles for this case.

Reliability analysis for TAMU-POST
The proposed model TAMU-POST has been presented as giving a deterministic solution. However, there are inherent uncertainties associated with the input parameters. A reliability analysis using a Monte Carlo Simulation was performed to evaluate the probability of failure for a pile group under a given vehicular impact. Probability of failure is defined here as the probability that the dynamic penetration of the vehicle will exceed 1 m as required by the ASTM standard (ASTM F2656-07). The probability of failure was quantified for two parameters: the soil limit pressure P L and the crushing factor κ. Table 6 summarizes the features of the statistical distribution considered for these two parameters.
The case considered a soil-pile barrier to contain a head-on impact by a 6800-kg medium-duty truck traveling at 100 km/h. The barrier consisted of an in-line group of eight W14x90 steel piles embedded 3 m in sand with the following properties: P L = 0.8 MPa, E 0 = 10 MPa, and ρ = 1900 kg/m 3 . The piles were spaced 5.2 m and connected by an HSS10x10x1/2 beam. The dynamic penetration was first obtained deterministically by using TAMU-POST as 0.03 m indicating that the barrier would satisfy the safety requirement (0.3 m < 1 M). The probability of failure was then calculated by using the Monte Carlo simulation as 10 −4 . Then, the pile embedment depth was varied from 1 m to 3 m to evaluate the impact of the embedment depth on the probability of failure. As shown in Figure 16, the Figure 15. Effect of number of piles on four-pile group response. Table 6. Random variables and statistical features (Briaud, 2013;Kulhawy & Trautwein, 1996) Parameter Distribution

Coefficient of variation
Model parameter Crushing factor κ Lognormal 0.05 Soil property P L (MPa) Lognormal 0.2 Figure 16. Probability density functions of dynamic penetration and probability of failurepile embedment depth ranging from 1 m to 3 m. embedment depth has a significant impact on the dynamic penetration and the probability of failure. Indeed, the probability of failure varied from 10 −4 for a 3 m embedment to 0.99 for an embedment depth of 1.5 m.

Summary and conclusions
A computationally efficient model called TAMU-POST was developed to predict the behavior of inline pile groups embedded in soil and subjected to the impact of a vehicle with a given mass and a given approach velocity. This model was based on the finite difference solution to the governing differential equation for a beam supported by piles represented by single degrees of freedom with damping, mass, spring, and a slider. This code was then calibrated against full-scale crash tests as well as a series of LS-DYNA finite element simulations. Recommendations were made on how to best obtain the input parameters.
TAMU-POST is much simpler than a computer simulation using FEM programs and much less costly to use than full-scale experiments. Yet, the predictions by TAMU-POST resulted in estimates comparable to those obtained from the full-scale crash tests and numerical simulations. The input parameters can be obtained from common laboratory and in-situ tests. These parameters are the mass and velocity of the approaching vehicle (M v and V v ), the soil strength quantified by the pressuremeter limit pressure P L or correlation thereof, the pile embedment depth (L), the ratio of pile spacing to pile width or diameter (s/d), the number of piles, and the pile and beam bending stiffness. A parametric study was performed to investigate the influence of each parameter; the results of this parametric study can assist in making more informed design decisions. The limitations of this study include the limited number of full-scale tests. Indeed, numerical simulations used to calibrate the proposed model were validated against only two full-scale experiments. Another limitation of the software is the fact that the plastic behavior of the structural elements (i.e., piles and beams) was not included to keep a simple solution for the differential equation. As such, while it has been shown that the model predicts the maximum deflection of the barrier reasonably well, it is not recommended to use this model for predicting the maximum bending moment in the piles and in the beams. Additional full-scale crash testing and numerical simulations would be required to refine the proposed simple model.