A hybrid FEM-NN optimization method to learn the physics-constrained constitutive relations from full-field data | Synapse