Dear Benhour, > > 2-Deviatoric strain > > dev_T_elastic = deviator(T_elastic); > > > > dev_T_elastic = T_elastic - (1/dim) * trace(T_elastic) * (Tensor<2, > > dim>(StandardTensors<dim>::I))/3.0; > > The expression '1/dim' is computed in integer arithmetic and evaluates to > zero. You want to write this as '1.0/dim'. That should do the trick. Want > to > try it out? >
I suspect that the trailing "/3.0" should not be there, because this factor is the one supposedly taken up by "1.0/dim". See, for example, equation 4.13 in this link <http://people.ee.ethz.ch/~luethim/pdf/script/pdg/chapter4.pdf>. Regards, Jean-Paul -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.