yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

Small mistake indeed, my bad.

"My" comment was :

// ΔεE=(Δl/l₀)(K/A) = Δσ=σ₁-σ ↔ Δl=(σ₁-σ)l₀(A/K)

It should have been (without the "₀")
// ΔεE=(Δl/l)(K/A) = Δσ=σ₁-σ ↔ Δl=(σ₁-σ)l(A/K)

And yes, I changed the meaning of grow like in the rest of the code (no 
need to multiply then divide).
Grow is a strain increment, no unit.

If it is not clear, remark that :
eps2-eps1 = log(h2/h)-log(h1/h) = log (h2/h1) =~ (h2-h1)/h1

The last equality explains the comment.

Btw, this equation might be wrong for sheared faces (L.159) :
Vector3r 
cellArea=Vector3r(cellSize[1]*cellSize[2],cellSize[0]*cellSize[2],cellSize[0]*cellSize[1]);



Václav Šmilauer a écrit :
> > > It it actually quite clear now: instead of Δl=(σ₁-σ)l₀(A/K) you compute
> > > Δl/l₀=(σ₁-σ)(A/K). In that case, cellGrow should not be called cellGrow
> > > but something like cellRelGrow.
> > >       
> > No, I compute deltaEpsilon = deltaSigma/stifness (with epsilon=log(...) ).
> > The commentary was still correct I think (that's why I didn't change it).
> >     
>
> I don't get what you computed then. Clearly the commentary wasn't
> correct, since it ("Δl=(σ₁-σ)l₀(A/K)") contains l₀ you removed from the
> multiplication. That's why said you actually computed Δl/l₀, putting the
> l₀ to the other side, thus changing meaning of cellGrow.
>
> v
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-dev
> Post to     : yade-dev@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-dev
> More help   : https://help.launchpad.net/ListHelp
>   


--

_______________
Chareyre Bruno
Maître de Conférences

Grenoble INP
Laboratoire 3SR - bureau E145
BP 53 - 38041, Grenoble cedex 9 - France
Tél : 33 4 56 52 86 21
Fax : 33 4 76 82 70 43
________________