<p>The computational accuracy of the finite element method has a direct relationship with the size of meshes. Excessively large element sizes degrade computational accuracy while overly small sizes increase both the number of elements and degrees of freedom. Adaptive remeshing method can effectively resolve the contradiction between the computational accuracy and the computational cost. In this paper, an adaptive remeshing algorithm based on a posteriori error estimation is developed by using Python-ABAQUS and the whole process of adaptive solution is implemented. By applying the Python scripts to the finite element analysis of several types of elastoplastic problems in three dimensions, it is numerically verified the effectiveness and adaptability of the resulting adaptive method based on the remeshing techniques for solving those stress concentration problems. Adaptive remeshing method significantly improves the computational accuracy through numerical comparison between adaptive remeshing method and finite element method in Sec. 3.1. Through the Python scripts, we can realize the automatic processing from geometric modeling, mesh generation to adaptive solution and it is convenient to repeatedly modify the models considered and the corresponding parameters.</p>