A moment–curvature-based constitutive model for interactive simulation of visco-plastic rods

• Published in: Engineering with computers 2024-10, Vol.40 (5), p.2971-2983
• Main Authors: Panneerselvam, Karthikeyan, De, Suvranu
• Format: Article
• Language: English
• Subjects: CAE) and Design; Calculus of Variations and Optimal Control; Optimization; Classical Mechanics; Computer Science; Computer-Aided Engineering (CAD; Constitutive models; Control; Curvature; Degrees of freedom; Dynamical systems; Fields (physics); Kinematics; Math. Applications in Chemistry; Mathematical analysis; Mathematical and Computational Engineering; Original Article; Plastic deformation; Plastic flow; Rods; Simulation; Spline functions; Stress relaxation tests; Systems Theory; Yield criteria
• ISSN: 0177-0667; 1435-5663
• DOI: 10.1007/s00366-023-01938-0

A moment–curvature constitutive model is proposed for the dynamic simulation of visco-plastic rods subject to time-varying loads and constraints at interactive rates. Smooth spline functions are used to discretize the geometry of the rod and its kinematics with the centerline coordinates as degrees of freedom (DOF) and scalar twist as degrees of freedom (DOF). The plastic curvature is defined as a uniformly varying field in contrast to localized lumped plasticity models, suitable for simulation of spatial rods that undergo uniform plastic deformation such as a cable or surgical suture thread. The yield criterion and plastic/visco-plastic flow rule are developed for spatial rods taking advantage of the availability of smooth moment–curvature fields using the spline-based formulation. With the Bishop frame field as a reference, the material curvatures are quantified using the twist degree of freedom, enabling tracking the plastic fields with scalar twist, thereby eliminating slopes as DOF. Taking advantage of the invariant sub-blocks and the sparsity of the dynamic system matrix arising from the numerical discretization, an hierarchical (H-matrix) solution approach is utilized for efficient computation. Uniform curvature bending tests and moment relaxation tests are performed to study the convergence behavior of the model. Several real-world tests involving contact are performed to demonstrate the applicability of the model in interactive simulations.