my apology, JxW was missed in my code , now added .
for (q_piont ...) {strain[q] = symmetrize(solution_grads[q]); stress[q] = (2*d_mu*strain[q] + d_landa*trace(strain[q])* unit_symmetric_tensor<dim>()); grad += stress[q]*strain[q]*fe_values.JxW(q) ; }
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
