search

ContactEngineering / Adhesion

Cohesive zone models and intermolecular interactions for contact calculations
https://contactengineering.github.io/Adhesion/
MIT License
5 stars 2 forks source link

Adhesive Plastic Simulation Cone not working with PowerLaw or Exponential interactions #78

Closed sitangshugk95 closed 1 year ago

sitangshugk95 commented 1 year ago

I tried to use a PowerLaw interaction instead of LJ82 (to see if they both give the same results), but I get the following error: image

The only changes made from the example file are encircled below: image

Is this an error in the source code, or is there something wrong in how I am implementing it?

sitangshugk95 commented 1 year ago

@sannant

sitangshugk95 commented 1 year ago

Error also when Exponential is used (albeit different error message):

image

sannant commented 1 year ago

The exponential and Powerlaw potential do not include repulsion and are not meant to be used in combination with linearize_core (which is used to implement plasticity).

Ideally, we should check whether the potential has a repulsion and throw a more verbose error in linearize_core, but that is not a priority.