Yade

Implementation of new contact law between particles

Asked by Yor1

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.

Question information

Language:
English Edit question
Status:
Answered
For:
Yade Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Yor1 (jabrane-hamdi) said : 		#1

Hello,

I posted my question yesterday et nobody reacts. May be there is lack of information in my post.
If you do not find the information that you deem necessary, I can provide them.

Best regards.
Jabrane.

Revision history for this message
Bruno Chareyre (bruno-chareyre) said : 		#2

Hi,
One day is very short, and this question is not really the sort of thing you can answer right after reading it. Some plots are not as you expect them, you ask why, then the answer has to be found by downloading/reading some source files.
Don't be surprised to not see an answer yet.
My advice would be, since you are modifying an existing functor: change it step by step and tell us which change (utlimately, which line no.) is the cause of the weird things.
If you don't do this pinpointing work yourself you let it to be done by the one who dare trying to help, hence minimizing the chance to get help.
Bruno

Revision history for this message
Yor1 (jabrane-hamdi) said : 		#3

Hello Bruno,

Thank you for response !
Ok, i will give the results of the sources files changing results step by step with hope
that Yade users could help me.

Best regards.
Jabrane.

Can you help with this problem?

Provide an answer of your own, or ask Yor1 for more information if necessary.

To post a message you must log in.