Hi
A comment.
I think you are talking about the (scalar) energy defined as e := 0.5 \sigma :
\epsilon where \sigma := C : \epsilon.
Why do you call it a gradient? The stress is the gradient of the free energy
wrt strain.
What theory are you referring to?
cheers
Andrew
On 05.04.2012, at 19:43, Mohammad Mohsenie wrote:
> Dear All
>
> i wrote a code to compute the gradient of a functional which is :
> gradient = E:grad(u).grad(u) in which E is elasticity tensor and u is the
> answer of elastic equation . supposing isotropic material i wrote code as :
>
>
> ...
>
> for (cell ...) {
> ...
>
> fe_values.reinit (cell);
> fe_values[displacement].ge
> 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 ) .
> i know that the last thing that can be done directly in dealii is double
> contraction of two tensor as i do here but i just want to know how to get a
> scalar as gradient i.e. how do total inner product .
>
> thanks in advance,
> S.M.Mohseni
> _______________________________________________
> dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
