yade-users team mailing list archive
Mailing list archive
Re: How to keep a constant stress boundary?
> The algorithm is : displacement = (stress offset) / (total
On one hand you compute the total stiffness on boundary N = the sum of
all sphere-box(N) contact stiffnesses.
On the other hand you have the current stress (or, say, the current
force if you prefer) applied on boundary N. From this you can define the
"offset", which is (target stress - current stress).
Then, the displacement computed in the equation above is the one "so
that the new force will be exactly the one we want after the
displacement". (assuming spheres are not moving in the meantime, which
is not true, but it will be corrected again on the next time step, so it
finally converge relatively well)
This algorithm has no name that I know.
kan a écrit :
I checked the code, but I still did not understand the physical
meaning of the code.
May I ask : what is the physical meaning behind the code? or what is
the mathematic of the algorithm?
how is it called ?
Thanks a lot.
On Fri, Sep 11, 2009 at 10:22 AM, Bruno Chareyre
<bruno.chareyre@xxxxxxxxxxx <mailto:bruno.chareyre@xxxxxxxxxxx>> wrote:
More precisely, it is in TriaxialStressController. This engine
needs to be assigned 6 boxes (one for each boundary). The boxes
are used to define the surfaces (since you need a surface to
define a stress).
If you don't have 6 boxes, you need to find a workaround.
The algorithm is : displacement = (stress offset) / (total
Maître de Conférences
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