A novel modeling technique based on the coupled Eulerian-Lagrangian (CEL) method is provided to solve the geotechnical problems with large deformations. The technique is intended to solve the update problem of soil mechanical properties during spudcan penetration in normally consolidated clay soil. In the CEL model, the normal method of assigning an increasing shear strength profile with depth (NA) is defective due to its Eulerian framework. In this paper, a new technique is proposed to update soil material properties by introducing thermo-mechanical coupled analysis (TMCA) to the CEL models. During establishment of the CEL models, the optimal penetration velocity and minimum mesh size are determined through parametric studies. Reasonability and accuracy are then verified through comparison of the preliminary results with the soil flow configuration and penetration resistance ( F v) of a centrifuge test, and the results of the proposed method are compared with those of the remeshing and interpolation technique with small strain (RITSS) method. To achieve a CEL model with satisfactory accuracy, the NA and TMCA methods implemented in the CEL models and the RITSS method are first adopted in weightless soil. Comparison of the findings with those obtained in previous studies shows that the TMCA method can update material properties and predict F v. The TMCA method is then applied to soils with self-weight and different shear strength profiles. Results show that the proposed method is capable of accurately modeling the large deformation problem of spudcan penetration in non-homogeneous clay.