Dear Getfem users,
i would like to take into account the fiber direction in my model;
to do so i computed the vector associated to each point of the mesh in the order given by mf.point_of_basic_dof, stored it in a vector, used an initialized fem data
model.add_initialized_fem_data("fibre",mf_u,Param_fibre);
and included it in the model (basically, add the component T_a*fibre*fibre^T to the second Piola-Kirchhoff tensor in the reference configuration)
getfem::add_nonlinear_generic_assembly_brick(model, mim,"T_a*(Id(3)+Grad_u)*[fibre(1)*fibre(1),fibre(1)*fibre(2),fibre(1)*fibre(3);fibre(2)*fibre(1),fibre(2)*fibre(2),fibre(2)*fibre(3);fibre(3)*fibre(1),fibre(3)*fibre(2),fibre(3)*fibre(3)]:Grad_Test_u");
The results i obtained are qualitatively correct but a visible difference still remains between the reference solution and my own solution.