fe_values[displacement].get_function_gradients (solution_d, solution_grads);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] ; } gradient(icell) = -0.5*grad ; } but in contrast of theory that says the gradient of each cell is a scalar , here in every cell gradient is a vector (i.e. has 2 component .one in u direction and one in v as you know ) .
Andrew already asked this in a different way but I still don't know what you mean by that. If your solution is a displacement vector [u,v] then the gradient of the solution is a rank-2 tensor (a matrix).
W. ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
