How to use gridap topology optimization while using 2 governing equations

I want to use gridap topology optimization for fracture optimization. so there is two governing equations one is for stress strain relation and the other is about fracture energy. how can i use gridap in this case? as i shoulf define matrix A and other matrices in tutorial 18 gridap.