[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] bug in von mises stress calculation?
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] bug in von mises stress calculation? |
Date: |
Mon, 31 Aug 2009 13:48:38 +0200 |
User-agent: |
KMail/1.9.9 |
I made an error in my computation. The factor sqrt(3.0/2) is indeed correct
(for 3D problems).
Yves.
On vendredi 28 août 2009, Andriy Andreykiv - 3ME wrote:
> Dear Yves and Getfem users,
>
> Want to report a bug. I was using von mises calculation in
> getfem::mdbrick_nonlinear_elasticity<>. It was giving a wrong value even
> though all the other results of the calculation were consistent with the
> benchmark. So I checked the code and found that in 3.0.1 version of getfem
> (the version which documentation is on the website) the final expression
> for the value is
>
> VM[i] = sqrt(3.0/2)*gmm::mat_euclidean_norm(sigma); (which is correct, I
> believe)
>
> while in 3.1 version the multiplier sqrt(3.0/2) is absent.
>
> Best regards,
> Andriy
--
Yves Renard (address@hidden) tel : (33) 04.72.43.87.08
Pole de Mathematiques, INSA-Lyon fax : (33) 04.72.43.85.29
20, rue Albert Einstein
69621 Villeurbanne Cedex, FRANCE
http://math.univ-lyon1.fr/~renard
---------