Contact law: issue if fragile=False and unpMax=0
Hi,
I am working with the contact law Law2_ScGeom6D_
I would like to use non brittle contact with respect to the shear force thus attribute fragile=False.
However I understand that according to:
if (phys->unpMax && phys->unp<
return false;
contact will be lost only if unpMax!=0. But if we want to have a non brittle contact in the shear direction and no normal plasticity then we would assign unpMax=0 and the contact will live foreover instead of breaking when (-Fn)> phys->normalAdh
Consequently there are different possibilites:
1/ change nothing and just assign unpMax=
2/ change:
if (phys->unpMax && phys->unp<
return false;
into
if (phys->
return false;
and if one will want to normal plasticity for ever he will have to put unpMax=
3/ but better than 2 I would prefer to change:
if (phys->fragile && (-Fn)> phys->normalAdh
// BREAK due to tension
return false;
} else {
if ((-Fn)> phys->normalAdh
Fn=-
phys->unp = un+phys-
if (phys->unpMax && phys->unp<
return false;
}
into
if ((phys->fragile || !phys->unpMax) && (-Fn)> phys->normalAdh
// BREAK due to tension
return false;
} else {
if ((-Fn)> phys->normalAdh
Fn=-
phys->unp = un+phys-
if (phys->
return false;
}
Any preferences, or advice?
Cheers,
Luc
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Luc Sibille
- Solved:
- Last query:
- Last reply: