Re: [Yade-dev] friction dissipation

2010-05-18 Thread chiara modenese
Hello Bruno, I send again to the list a couple of mails I addressed to you. Chiara On 16 May 2010 01:23, chiara modenese c.moden...@gmail.com wrote: On 6 May 2010 23:11, chiara modenese c.moden...@gmail.com wrote: Hi Bruno, could you explain me in formula how did you get the plastic

Re: [Yade-dev] friction dissipation

2010-05-18 Thread Bruno Chareyre
I have one big guess. How can this formulation being incremental using trialForce? Let's take a simple case: normal force is set to a constant value. Then say that trialForce will increase. In this case shearForce is the same since normal force is as well (like maintaining

Re: [Yade-dev] friction dissipation

2010-05-18 Thread chiara modenese
On 18 May 2010 16:57, Bruno Chareyre bruno.chare...@hmg.inpg.fr wrote: No. It is the increment of displacement. Example : 1) time t : us_total=1, us_elastic=0.1, us_plastic=0.9, fs=0.1*ks 2) assume an increment dus between t and t+dt 3) trialFs = previousFs + dus*ks = us_elastic*ks +

Re: [Yade-dev] friction dissipation

2010-05-18 Thread Bruno Chareyre
If you agree could you commit it, please? Thanks. Oh, sorry, I forgot that point... :-/ I'm sending you the sources, because I have also diffs for the OPenMP plastic accumulator and I need to compile before commiting (I'm already compiling another trunk and it will take time). I will commit

Re: [Yade-dev] friction dissipation

2010-05-18 Thread chiara modenese
On 18 May 2010 17:41, Bruno Chareyre bruno.chare...@hmg.inpg.fr wrote: If you agree could you commit it, please? Thanks. Oh, sorry, I forgot that point... :-/ I'm sending you the sources, because I have also diffs for the OPenMP plastic accumulator and I need to compile before commiting

Re: [Yade-dev] friction dissipation

2010-05-15 Thread chiara modenese
On 6 May 2010 23:11, chiara modenese c.moden...@gmail.com wrote: Hi Bruno, could you explain me in formula how did you get the plastic dissipation for the friction component? I see plasticDissipation += (shearDisp+(1/currentContactPhysics-ks)*(shearForce-prevForce)).dot(shearForce) Is

[Yade-dev] friction dissipation

2010-05-06 Thread chiara modenese
Hi Bruno, could you explain me in formula how did you get the plastic dissipation for the friction component? I see plasticDissipation += (shearDisp+(1/currentContactPhysics-ks)*(shearForce-prevForce)).dot(shearForce) Is shearDisp the total shear displacement, right? I would compute the