We combine 2-D ice flow model with linear elastic fracture mechanics (LEFM) to model the calving behavior of Thwaites Glacier, West Antarctica. We find the combination of full-Stokes (FS) model and LEFM produces crevasses that are consistent with observations. We also find that calving is enhanced with pre-existing surface crevasses, shorter ice shelves or undercut at the ice shelf front. We conclude that the FS/LEFM combination is capable of constraining crevasse formation and iceberg calving.