[deal.II] Elastic energy computation

2017-01-16 Thread 'Seyed Ali Mohseni' via deal.II User Group
Hi, Is it possible to compute the elastic energy for a simple 2D cube before we solve the system? Because I checked the displacement gradients and they are 0 due to having no solution probably. Hence, I wonder how would someone be able to compute the elastic energy with the following formula:

Re: [deal.II] Re: Elastic energy computation

2017-01-16 Thread Jean-Paul Pelteret
Dear Seyed, These two things are very different, and cannot be compared with one another. The “strain” in the assembly routine stems from the weak form, which is presented in the introduction to the tutorial. If one considers a variational approach, is coincides with the variation of the

[deal.II] Re: Elastic energy computation

2017-01-16 Thread 'Seyed Ali Mohseni' via deal.II User Group
I checked the second one. It doesn't coincide with the strain tensor I obtain from the energy computation after I solve the system, e.g. - in GAUSS point 1 - ASSEMBLY STRAIN TENSOR 0 0.105662 0.105662 0.211325 ENERGY STRAIN TENSOR 0.001081 -0.000738 -0.000738 -0.008525 ... etc.

[deal.II] Re: Elastic energy computation

2017-01-16 Thread Jean-Paul Pelteret
Dear Seyed, There are two get_strain() functions in step-18. The first computes the symmetric part of the shape function gradients, and the second computes the small strain tensor. To store the values of the

[deal.II] Re: Elastic energy computation

2017-01-16 Thread 'Seyed Ali Mohseni' via deal.II User Group
Thanks, I understand now. I just deleted the post, because I misunderstood what you meant earlier. I thought you are refering to the left eps_phi_i (first) and right eps_phi_j (second). But you were talking about: template inline SymmetricTensor<2,dim> get_strain (const

[deal.II] test stage for compiling deal.II suite in Spack 0.10 (and later 1.0)

2017-01-16 Thread Denis Davydov
Dear all, Supercomputer PACKage manager currently has a test stage for the upcoming 0.10 release (*), see this comment . If you are using Spack, it's* good time to test that nothing is broken*