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 <[email protected]> > 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, [email protected] wrote: > > ------------------------------------------------------------ > revno: 3797 > committer: Francois <[email protected]> <[email protected]> > 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 : [email protected] > 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 : [email protected] > Unsubscribe : https://launchpad.net/~yade-dev > More help : https://help.launchpad.net/ListHelp > >
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
