Fig. 27: Comparison of envelope curves with different crack positions as depicted in Fig. 26 with l/h = 1, l = 0.01 mm and kφ = 0.01
 

The error e^num_Γ (Eq. 61) of the numerically obtained crack length in a 2D simulation with a transformed phase-field ansatz in comparison to the expected solution 0f 0.1mm  is plotted over the number of subtriangles which are utilized for the numerical integration for different given crack positions. The simulation setup is given in Fig. 26, the position of the crack is enforced with a penalty formulation.

for each file -> Columns
1: number of sub triangles per direction
2: error 

file fig_27_data_1.csv:  yc=0.0mm
file fig_27_data_2.csv:  yc=0.05mm
file fig_27_data_3.csv:  yc=0.1mm
file fig_27_data_4.csv:  yc=0.25mm
file fig_27_data_5.csv:  yc=0.5mm

