The result shown in the example problem is incorrect. The volume of the element should be the square of the global seed size of the mesh multiplied by the thickness of the individual ply, assuming there is an element through the thickness for each ply layer.
In this case, the volume would be V = 0.05 * 0.05 * 0.0072 = 0.000018 in^3.
We will have to investigate further where the reported value of 0.0115 in^3 comes from. Of course, changing the volume will also change the calculated energies, G.
As to your second question, Simulation Composite Analysis does support Puck, Christensen, and Hashin as failure initiation criteria for EBD with unidirectional materials. Only MCT, Max Stress, Max Strain, and user-specified failure criterion are supported for woven materials.
Thanks