Abstract
An algorithm designed for the determination of equilibrium shocks that appear in quasi-static evolution problems associated to elastic nonmonotonous stress-strain laws is presented in the context of one-dimensional media. Two basic procedures are involved in the proposed method: (i) enhancement of the finite element in order to describe the weak discontinuities in any point of its interior and (ii) implementation of a return mapping algorithm for the determination of the shocks, which have to satisfy the inequality constraints imposed by a maximally dissipative hypothesis. A rigorous proof of the unconditional stability property of the algorithm is also given. The present study is applied to the theoretical model presented by Abeyaratne and Knowles in the context of one-dimensional extensional deformations of bars. The numerical results are in complete agreement with the analytical ones.
Sign-in/Register to access full text options 
Free) 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call
