本发明公开了一种模拟裂纹自适应扩展的数值算法，步骤S1、根据裂纹拓展的目标问题确定出求解区域，以及将求解区域中裂纹未扩展的区域作为连续求解域，再利用克隆点技术在连续求解域内进行裂纹拓展的求解计算；步骤S2、在连续求解域内的网格上作用最大拉应力和Mohr‑Coulomb复合强度准则，以作为裂纹萌生以及扩展方向的依据，并在裂缝萌生的区域进行网格重构；步骤S3、采取空间映射技术将原有网格信息继承至新生成网格当中，并在新增裂纹面处设置接触面，通过在接触面内增加法向和切向弹簧进行裂纹面的力学特征描述。本发明实现了裂缝的动态扩展过程模拟，克服了裂缝扩展路径受制于初始网格拓扑的问题，实现了裂缝扩展方向的精确模拟。