Abstract: The penetration of anchors in the seabed affects the subsequent towing behavior and poses a great threat to pipelines, cables and tunnels laid on the seabed. Based on the coupled Euler-Lagrangian large deformation finite element method, a model for calculating the penetration of anchors in the cohesive soil seabed is established. The variation of clay strength and gravity with depth is considered. The comparison between the numerical results and the model test results shows a good agreement. The development laws of anchor force, velocity and penetration depth during the penetration of the anchor into the seabed are revealed. The influences of the parameters of soil strength, initial bottoming speed and other parameters on the penetration of the anchor are analyzed.