算例是金屬棒被套在一段由梁單元組成的環內部，對金屬棒施加Y方向位移，對環產生力直至拔斷環。

金屬棒是離散剛體，除開U2方向位移外的其他一切自由度均為固定，環是B31梁單元組成的，環是彈塑性本構，想要模擬塑料絲被金屬拉斷的效果，左邊是我自己編寫的彈塑性vumat，右邊是abaqus內置彈塑性本構，現在的問題就是子程序與自帶本構力學響應完全不一樣，這兩個算例只改了本構，其他邊界條件，接觸等等都沒變，這是什么原因導致的呢？