“…In the p-y method the pile is assumed to behave as an Euler-Bernoulli beam with the soil modelled as a series of discretely spaced springs, each connected to one of the pile segments into which the pile is discretised. The springs model the soil response to loading through p-y curves ( p is the unit resistance per unit pile length offered by the springs, and y is the pile deflection), which are developed empirically by adjusting the curves until they match actual load-displacement results Briaud et al, 1984;Yan & Byrne, 1992;Brown et al, 1994;Gabr et al, 1994;Briaud, 1997;Wu et al, 1998;Bransby, 1999;Ashour & Norris, 2000). However, the p-y method often fails to predict pile response (Anderson et al, 2003;Kim et al, 2004), for it is not capable of capturing the complex three-dimensional interaction between the pile and the soil.…”