yade-users team mailing list archive

[Jan Stránský <honzik.stransky@xxxxxxxxx>]

by "adding the dissipation(Frictional, Damping etc) also to my total energy
at every 50 steps" do you mean this?

def myFunction():
        pe_normal  = Mindlin.normElastEnergy() # normal PE
        pe_shear  = Mindlin.shearEnergy # shear PE
        norm_damp_dissip = Mindlin.normDampDissip # normal diss
        shear_damp_dissip = Mindlin.shearDampDissip # shear diss
        fric_dissip = Mindlin.frictionDissipation # loss due to friction
        E_tracker         = dict(O.energy.items())
        E_kin_trans = E_tracker['kinTrans']
        E_kin_rotation = E_tracker['kinRot']
        etotal =  pe_normal + pe_shear + E_kin_trans + E_kin_rotation +
fric_dissip + norm_damp_dissip + shear_damp_dissip # total energy

if so, by my last post I meant that Mindlin.normDampDissip,
Mindlin.shearDampDissip, Mindlin.frictionDissipation
might be lost during save/load process. Because it is only guess, I asked
you to check the values of these variables before save and how much energy
is not conserved after load. In the ideal case, the values would be equal..

cheers
Jan




2014-12-07 19:51 GMT+01:00 ipemath <question258679@xxxxxxxxxxxxxxxxxxxxx>:

> Question #258679 on Yade changed:
> https://answers.launchpad.net/yade/+question/258679
>
>     Status: Answered => Open
>
> ipemath is still having a problem:
> Hi Jan Stransky,
>
> Actually I f you see the python script I have posted, I am adding the
> dissipation(Frictional, Damping etc) also to my total energy at every 50
> steps. Since I am accounting for the dissipation in the total energy, it
> should be conserved isnt?
> Can I mail the scripts to somebody who can help me???
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>