A three-dimensional approximate computer-based method, YSPR (Yonsei Piled Raft), was developed for analysis of behavior of piled raft foundations. The raft was modeled as a flat shell element having 6 degrees of freedom at each node and the pile was modeled as a beam-column element. The behaviors of pile head and soil were controlled by using $6{ imes}6$ stiffness matrix. To model the non-linear behavior, the soil-structure interaction between soil and pile was modeled by using nonlinear load-transfer curves (t-z, q-z and p-y curves). Comparison with previous model and FEM analysis showed that YSPR gave similar load-displacement behaviors. Comparison with field measurement also indicated that YSPR gave a reasonable result. It was concluded that YSPR could be effectively used in analysis and design of piled raft foundations.