由于从实验及理论角度研究岩样单轴拉伸条件下的破坏全过程及尺寸效应难度都很大。因此采用拉格朗日元法来研究这些问题。在峰值强度之前后,岩石材料的本构模型分别取为线弹性及拉破坏线性应变软化模型。为了使拉伸塑性区不出现在试样的端部,在试样的两侧面中部预制了2个凹槽。数值模拟结果表明,全程拉应力-拉应变曲线分为峰前和峰后阶段。在接近峰值的峰前阶段,由于两凹槽附近具有明显的拉应力集中现象,拉伸塑性区最先出现在两凹槽附近。随着轴向拉应变的增加,发生拉伸破坏的单元的数目增加,新发生拉伸破坏的单元越来越接近试样的中心,直到两块拉伸塑性区在应变软化阶段贯通。两凹槽连线上各单元拉应力的分布呈现3个阶段,"澡盆型"("U 型")阶段,"双峰型"("M 型")阶段及"单峰型"("П型")阶段。"澡盆型"阶段对应于全程拉应力-拉应变曲线的弹性阶段。"双峰型"阶段及"单峰型"阶段对应于全程拉应力-拉应变曲线的非弹性阶段(包括峰值强度之前的一小段,即应变硬化阶段及峰后的应变软化阶段)。增加试样高度及降低试样宽度,拉应力-拉应变曲线的软化段变得越来越陡峭,因而试样越容易发生失稳破坏。由于试样宽度较大时,试样内部的单元并非处于单向拉应力状态,因此,增加试样宽度,全程拉应力-拉应变曲线的峰值强度增加。当试样宽度较小时,从出现塑性区,到塑性区贯通所需要的时间步较小,或应变范围较窄。这说明试样的脆性较强,前兆不明显。前兆不明显的脆性破坏对应常见的是洞室岩爆、冲击地压及地震等灾害。