yade-users team mailing list archive

Thread
Date

[Question #612245]: Implementation of new contact law between particles

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Yor1 <question612245@xxxxxxxxxxxxxxxxxxxxx>
Date: Wed, 05 Apr 2017 19:43:20 -0000
Reply-to: question612245@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

New question #612245 on Yade:
https://answers.launchpad.net/yade/+question/612245

Hello Yade users,

I try to implement a contact law with softening factor in the normal direction contact illustrated in the Figure 1 of the paper of Scholtes and Donze 2013 ( the paper can be downloaded from this  link : https://filex.univ-lorraine.fr/get?k=afQQ74EoCQ5ZIy35Iz1).

I implement this new contact in JointedCohesiveFrictionalPM.cpp from the line 153 to the line 163.
 JointedCohesiveFrictionalPM.cpp : https://filex.univ-lorraine.fr/get?k=KbdVLas2oIUUyKkG9nd
 JointedCohesiveFrictionalPM.hpp : https://filex.univ-lorraine.fr/get?k=Epw6QzDIm3HEjT9dH7d

After the implementation, i simulate a trixial compression test and i obtain a strange stress/strain curves in the beginning of the simulation illustrated in the figure :  https://filex.univ-lorraine.fr/get?k=dzDNhQKp3I5OXmaopWz
I don't find the fault in the implementation.

the script of the triaxial compression curves plotting : https://filex.univ-lorraine.fr/get?k=JVVacWCTBRjhLEpplth
the script of the triaxial compression : https://filex.univ-lorraine.fr/get?k=BkLLO86ERwC6a3Cc35h
the used sample : https://filex.univ-lorraine.fr/get?k=1UqqWzFs3tXJXhxZbVN

Best regards.
Jabrane.

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.