fix(tet-based-metrics): compliance shear entries missing (+1 more)

[andrewwhitecdw]

Small fixes in utils/tet_based_metrics.py:

fix: compliance shear entries missing factor of two in strain energy

Fix: Replace:

    D_inv[3][3] = invE * (1.0 + nu)
    D_inv[4][4] = invE * (1.0 + nu)
    D_inv[5][5] = invE * (1.0 + nu)

with:

    D_inv[3][3] = 2 * invE * (1.0 + nu)
    D_inv[4][4] = 2 * invE * (1.0 + nu)
    D_inv[5][5] = 2 * invE * (1.0 + nu)

fix: voigt stress vector uses wrong xy shear component

Fix: Replace:

    voigt_stress = np.array([ts.x.x, ts.y.y, ts.z.z, ts.z.y, ts.z.x, ts.y.z])

with:

    voigt_stress = np.array([ts.x.x, ts.y.y, ts.z.z, ts.z.y, ts.z.x, ts.y.x])

Files changed

• utils/tet_based_metrics.py