A basis for developing a general approach to solve geotechnical engineering problems through dynamic finite element analysis of nonlinear porous media is presented. A new series of constitutive models named here as δ* series are developed under the general hierarchical single surface (HISS) modeling approach, to include the behavior of, saturated normally consolidated and overconsolidated clay, under drained or undrained, static and cyclic loading conditions. Algorithms for determination of material parameters for these models from laboratory models are also developed. Constitutive parameters for Sabine Clay are obtained using triaxial test results from undisturbed samples and the model is verified by back predicting the laboratory behavior of this clay. Sensitivity analyses for all the material parameters have been also carried out. Efficient and reliable algorithms for calculating strain increments for given stress increments and vice versa are developed. All the models in the δ* series are implemented in the finite element program POROUS which is based on the theory of dynamics of nonlinear porous media. Here a modular approach is used to facilitate easy modification of all the functions associated with these models (yield function, potential function, hardening function and interpolation functions). Complete test procedure for field load tests on two pile segments (3 inch and 1.72 inch) are numerically simulated and compared with field measurements. Initial stresses before pile driving are estimated using results of s self boring pressuremeter test on the site. Strain distributions just after pile driving are evaluated using the strain path method. The effective stress distribution is obtained by integrating constitutive equations for given strain paths. Corresponding total stresses and pore pressure distributions are obtained using the computer program POROUS. Consolidation after pile driving and all the static and cyclic tests followed are also simulated using the program POROUS. The predicted normalized pore pressure dissipation curve matches the field behavior. This analysis provides good predictions of shear transfer from which the pile capacity can be evaluated. Even though the variation of pore pressure during a cycle is not matched exactly, the accumulation of pore pressures are predicted well. The demonstrated ability of dynamic finite element analysis of nonlinear porous media, to simulate slow consolidation and cyclic load tests provides a basis for developing a general approach for solving geotechnical engineering problems.