Thread Previous • Date Previous • Date Next • Thread Next |
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. Bruno kan a écrit :
Thanks, Bruno,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.yongfengOn 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 box-spheres stiffness) Bruno
-- _______________ 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 ________________
Thread Previous • Date Previous • Date Next • Thread Next |