The invention discloses a viscoelastic-plastic numerical
simulation method, which comprises the following steps of: receiving physical modeling and solving
control data of a viscoelastic-plastic material, constructing a convex cone planning model for uniformly describing elastoplasticity,
viscoelasticity and friction multi-
physical field coupling behaviors according to the physical modeling and solving
control data, expressing a plastic yield criterion and a
viscous flow law as second-order cone constraints, converting the model into a second-order cone complementary equation; converting the equation into an implicit immobile point iterative equation by adopting a
modular system matrix splitting iterative
algorithm, performing matrix splitting on a
system stiffness matrix, and constructing an iterative solution format with calculation complexity at a quadratic level; and performing iterative calculation to convergence based on the format and the input data, and outputting numerical results of a
displacement field, a
stress field and a strain field. According to the method, multi-field
coupling behaviors are unified to a convex cone planning framework, solving is conducted through an efficient iterative
algorithm, the calculation efficiency and stability of large-scale viscoelastic-plastic numerical
simulation are improved, and the problem of calculation
bottleneck caused by existing
commercial software multi-field decoupling and general
solver explicit inversion is solved.