I simulated a beam in 4 point bending that should fail in shear.
If i use NodeTetConstraintList with penalty or BeamParticleConstraintList with RebarConcreteInteraction, i get the same elastic response regardless of the mars type i use (i.e. marsD, marsO and marsM)
If i use BeamTetConstraintList with RebarConcreteInteraction, i get the same results with others only if i use marsD, but if i use marsM or even marsO, i get extremely higher stiffness