A Phantom-Node Method with Edge-Based Strain Smoothing for Linear Elastic Fracture Mechanics

<table>The convergence in the stress intensity factor <svg height="16.375" id="M151" style="vertical-align:-3.27605pt" version="1.1" viewbox="0 0 20.262501 16.375" width="20.262501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.625,11.675)"><path d="M764 650l-7 -26q-61 -10 -91.5 -22.5t-76.5 -48.5l-235 -184q104 -152 208 -276q31 -37 54.5 -48.5t68.5 -16.5l-7 -28h-156q-138 170 -213 284q-19 29 -33.5 35t-33.5 -7l-29 -173q-15 -75 -5.5 -90.5t74.5 -20.5l-5 -28h-260l7 28q58 5 74 21.5t30 89.5l70 378
q12 68 2 83t-72 22l6 28h257l-7 -28q-62 -7 -76 -21.5t-27 -83.5l-32 -172q29 11 78 46q112 84 217 179q18 16 23.5 25.5t-1 14.5t-26.5 8l-30 4l5 28h249z" id="x1D43E"></path></g>
<g transform="matrix(.012,-0,0,-.012,13.812,15.75)"><path d="M414 650l-6 -28q-63 -4 -78.5 -17.5t-27.5 -71.5l-77 -416q-11 -58 0.5 -71.5t77.5 -17.5l-5 -28h-275l8 28q63 4 79.5 18.5t27.5 70.5l80 416q11 57 -2 71t-81 18l6 28h273z" id="x1D43C"></path></g>
</svg> versus <svg height="13.55" id="M152" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 10.0125 13.55" width="10.0125" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.625,12.725)"><path d="M499 88q-41 -46 -83.5 -73t-68.5 -27q-41 0 -19 98l52 225q10 43 8 59.5t-16 16.5q-37 0 -112 -65t-117 -135q-14 -42 -38 -179q-36 -7 -75 -20l-7 8l127 594q8 40 2.5 49t-35.5 9h-35l5 25q36 3 71.5 13t56.5 18t25 8q11 0 3 -38l-88 -417h2q77 102 195 167q44 24 79 24
q61 0 30 -130l-52 -220q-8 -33 4 -33q8 0 31 14.5t43 33.5z" id="x210E"></path></g>
</svg> (mesh size) for a sheet with an edge crack under shear.</table>

Journal of Applied Mathematics

fig20

Figure 20

Figure 20: A Phantom-Node Method with Edge-Based Strain Smoothing for Linear Elastic Fracture Mechanics 