Seeking Advice on Implementing Variable Stiffness for Particles
Hi everyone,
I am using the yade version 20240107 and python version 3.10.12. I want to modify the software code to assign distinct stiffness parameters for each particle in a simulation. I have devised two potential approaches towards achieving this:
Approach 1: In the C++ source code responsible for constructing particles, I plan to introduce a new variable that will store the individual stiffness values for each particle. This variable would then be integrated into the contact force calculation routines so that each particle's behavior reflects its unique stiffness value.
Approach 2: Alternatively, without altering the core Yade code, I could create a separate material for each particle, all sharing the same type but differing only in their stiffness parameters. Each particle would then be assigned to its corresponding material.
My question revolves around computational efficiency. Would Approach 1 likely result in less overhead due to the introduction of fewer "additional data" per particle, thus potentially having a smaller impact on computation speed compared to Approach 2 where there would be an increase in the number of material groups? Furthermore, I am open to suggestions on whether there exists a more optimal method to achieve this.
Thank you all and look forward to your valuable suggestions.
xue