Xa 0.6[m] "Crack length"
da 0.01[m] "Crack length increment"
Wp 1.5[m] "Plate width"
Hp 1.5*Wp "Half length of plate"
Th 10[mm] "Plate thickness"
Q0 20[MPa] "Load intensity"
AP 1.4E-11 "Paris' law coefficient (in MN/m^1.5 units)"
mP 3.1 "Paris' law exponent"
ar pi*Xa/(2*Wp) "Intermediate crack length parameter"
K1r Q0/1[Pa]*sqrt(pi*Xa/1[m])*sqrt(tan(ar)/ar)/cos(ar)*(0.752+2.02*(Xa/Wp)+0.37*(1-sin(ar))^3) "Reference stress intensity factor"