Based on the constitutive equations given by Yang and Ahmed et al.[1-5] to predict the liquefied large deformation of sand

the hardening rule of model and the calculation method of elastoplastic modulus are modified. Then

it is extended to 3-D equations and is implemented by using Fortran subroutine program and is implanted into ABAQUS finite-element software. By using ABAQUS software with the advanced model

how the main model parameters to affect the dynamic pore pressures arising and the stress-strain curves of sand

are analyzed in detail

respectively. At the same time

the sensitivities of model parameters are also analyzed. It is proved that the subroutine program can work perfectly to model the liquefied large deformation of sand; and the combination of different model parameters’ values can predict perfectly the different deformation stage of sand. These conclusions verified and given here can be used expediently to fixed the model parameters by using dynamic triaxial tests. Accordingly

a new effectively numerical method used to calculate the large liquefaction is advanced by using the advanced finite-element ABAQUS software with special constitutive model.