[deal.II] Re: Elastic energy computation

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

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

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

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

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