incorrect computation of potential energy for frictional interactions

Bug #831910 reported by Dion Weatherley
10
This bug affects 1 person
Affects Status Importance Assigned to Milestone
ESyS-Particle
Fix Released
Medium
SteffenAbe

Bug Description

Hi esys developers and users,

As I see the function getPotentialEnergy for the linear contact model, it only accounts for the normal potential energy. I'm thinking is there any tangential potential energy if the sliding does not occur (the tangential force is smaller than the threshold)? Since the tangential part is also considered as a spring with stiffness k_s, it will contain some potential energy, right? The energy is something like 0.5*f_s*f_s/k_s.

Best regards!
Ning

Related branches

Changed in esys-particle:
importance: Undecided → Medium
SteffenAbe (s-abe)
Changed in esys-particle:
assignee: nobody → SteffenAbe (s-abe)
SteffenAbe (s-abe)
Changed in esys-particle:
status: New → Confirmed
Revision history for this message
SteffenAbe (s-abe) wrote :

Added missing tangential part as Ning suggested. The current implementations is still only an approximation (using total current tangential force vs. incremental updates), but tests suggest an accuracy within a few percent (or better).
DOES NOT FIX WRONG POTENTIAL ENERGY IN CHertzianViscoElasticFrictionInteraction !

Changed in esys-particle:
status: Confirmed → Fix Committed
Vince Boros (v-boros)
Changed in esys-particle:
status: Fix Committed → Fix Released
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Related questions

Remote bug watches

Bug watches keep track of this bug in other bug trackers.