The shape memory effect observed in many alloys arises due to stress induced transformation between variants of the martensitic phase. It is difficult to study this process in detail using continuum approaches and particle based methods are eminently more suitable. In this work, we study detwinning, which is the transformation between martensite variants due to applied stress using a novel discrete particle model. The approach uses a novel multibody interparticle interaction defined directly using a non-convex free energy potential pertinent to such material behavior. This model is able to describe simultaneously occurring multiple microscale events during the detwinning process: nucleation and propagation of ledges along twin boundary. Due to these underlying microscale events the plateau region of stress–strain response shows a jerky nature. The effect of temperature and morphology on the stress–strain behavior of the self-accommodated martensite microstructure is studied in detail. From the simulations, we identify the morphological features affecting the transformation stress for detwinning. The critical parameters are found to be the length and number of twin and macro-twin boundaries and the number of mobile ledges.