Re: [deal.II] general question

2016-12-07 Thread Jean-Paul Pelteret
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::I))/3.0; > > The expression '1/dim' is computed in integer arithmetic and evaluates to > zero. You want to w

Re: [deal.II] general question

2016-12-06 Thread Wolfgang Bangerth
On 12/06/2016 06:49 PM, benhour.amiria...@gmail.com wrote: Dear All, I hope you have had a great time. I have a general question. I define total strain and deviatoric part of the elastic strain tensor in two different ways as below: 1- Total strain tensor: const SymmetricTensor<2, dim> total_str

[deal.II] general question

2016-12-06 Thread benhour . amirian66
Dear All, I hope you have had a great time. I have a general question. I define total strain and deviatoric part of the elastic strain tensor in two different ways as below: 1- Total strain tensor: const SymmetricTensor<2, dim> total_strain = 0.5*(Grad_u + transpose(Grad_u)); (Tensor<2, dim>