yade-dev team mailing list archive

[François <francois.kneib@xxxxxxxxx>]

Hi Bruno,

You're right, when I wrote the code I assumed that there is 1 material per
body, which is totally wrong. I made the code this way because I didn't
want to change Kn calculation etc. Two solutions : the first one is yours,
the second is to change stiffness calculation depending on the case. I
prefer yours, because there is no sense to use frictMat with E,ν and
viscElMat with kn,ks in the same simulation.
I'll change it and add an error message, then the code will crash. What do
you think about that ?



2014/1/15 Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>

>  Hi François,
> I think there is a problem in the Ip2_FrictMat_FrictMat_FrictPhys  [1].
> The functor is supposed to define the IP, not to modify the materials. A
> problem is that it is not thread safe. Many interactions could share the
> same material and modify it at the same time.
> Besides, I don't get logic. If mat1 and mat2 are modified according to
> rad1 and rad2, they will be changed all the time when the materials are
> shared by many particles of different sizes. Is that correct?
> Couldn't we assume that the user combining FrictPhys and ViscElMat will
> have to define young and poisson?
>
> Cheers
>
> Bruno
>
>
> [1]
> https://github.com/yade/trunk/commit/6f7318d4c47830275498450a1a80874a6da3490a#diff-6311c0b0529c49d12986e8bfa90017b2R36
>
> On 08/01/14 16:16, noreply@xxxxxxxxxxxxx wrote:
>
> ------------------------------------------------------------
> revno: 3797
> committer: Francois <francois.kneib@xxxxxxxxx> <francois.kneib@xxxxxxxxx>
> timestamp: Wed 2014-01-08 15:36:59 +0100
> message:
>   Allow the contact between (coh)frictMat and viscElMat as frictPhys. Almost everything was done by inheritance, just had to convert stiffnesses to modulus and modulus to stiffnesses to ensure material compatibility.
>   Note that for the moment the timeStepper cannot handle this kind of simulations -> will be fixed soon.
> modified:
>   pkg/dem/Ip2_FrictMat_FrictMat_FrictPhys.cpp
>   pkg/dem/ViscoelasticPM.cpp
>   pkg/dem/ViscoelasticPM.hpp
>
>
> --
> lp:yadehttps://code.launchpad.net/~yade-pkg/yade/git-trunk
>
> Your team Yade developers is subscribed to branch lp:yade.
> To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription
>
>
>
> _______________________________________________
> 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
>
>
>
> --
> _______________
> Bruno Chareyre
> Associate Professor
> ENSE³ - Grenoble INP
> Lab. 3SR
> BP 53
> 38041 Grenoble cedex 9
> Tél : +33 4 56 52 86 21
> Fax : +33 4 76 82 70 43
> ________________
>
>
> _______________________________________________
> 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
>
>