三個模型，一個是2D Contour integral法，一個是3D?Contour integral法，還有一個3D XFEM方法求解模型，輸出應力，應變及應力強度因子K1，并對比結果。