A macroelement is developed for soil-structure interaction analyses of piles in liquefiable soils, which captures efficiently the fundamental mechanisms of saturated granular soil behavior. The mechanical model comprises a nonlinear Winkler-type model that accounts for soil resistance acting along the circumference of the pile, and a coupled viscous damper that simulates changes in radiation damping with increasing material non-linearity. Three-dimensional (3D) finite element (FE) simulations are conducted for a pile in radially homogeneous soil to identify the critical parameters governing the response. The identified parameters, i.e., hydraulic conductivity, loading rate of dynamic loading, dilation angle and liquefaction potential are then expressed in dimensionless form. Next, the macroelement parameters are calibrated as a function of the soil properties and the effective stress. A semi-empirical approach that accounts for the effects of soil-structure interaction on pore pressure generation in the vicinity of pile is used to detect the onset of liquefaction. The predictions are compared with field data obtained using blast induced liquefaction and centrifuge tests and found to be in good agreement. Finally, the macroelement formulation is extended to account for coupling in both lateral directions. FEM simulations indicate that response assuming no coupling between the two horizontal directions for biaxial loading tends to overestimate the soil resistance and fails to capture features like 'apparent negative stiffness', 'strain hardening' and 'rounded corners'.