zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Resetting tangential overlap when slip occurs

Browse Source
This commit is contained in:
graham
2009-09-17 18:20:56 +01:00
parent bc4a71af71
commit 3d2d8f3fd8
1 changed files with 27 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

31
src/lagrangian/intermediate/submodels/Kinematic/CollisionModel/PairCollision/PairFunction/SpringSliderDashpot/SpringSliderDashpot.C
View File

@ -211,11 +211,34 @@ void Foam::SpringSliderDashpot<CloudType>::evaluatePair
scalar& etaT = etaN; scalar& etaT = etaN;
// // Tangential force
// vector fT_AB =
// -min(kT*mag(tangentialOverlap_AB), mu_*mag(fN_AB))
// *tangentialOverlap_AB/mag(tangentialOverlap_AB)
// - etaT*USlip_AB;
// Tangential force // Tangential force
vector fT_AB = vector fT_AB;
-min(kT*mag(tangentialOverlap_AB), mu_*mag(fN_AB))
*tangentialOverlap_AB/mag(tangentialOverlap_AB) if (kT*mag(tangentialOverlap_AB) > mu_*mag(fN_AB))
- etaT*USlip_AB; {
// Tangential force greater than sliding friction, particle slips
fT_AB =
-mu_*mag(fN_AB)
*tangentialOverlap_AB/mag(tangentialOverlap_AB)
- etaT*USlip_AB;
tangentialOverlap_AB = vector::zero;
tangentialOverlap_BA = vector::zero;
}
else
{
fT_AB =
-kT*mag(tangentialOverlap_AB)
*tangentialOverlap_AB/mag(tangentialOverlap_AB)
- etaT*USlip_AB;
}
pA.f() += fT_AB; pA.f() += fT_AB;
pB.f() += -fT_AB; pB.f() += -fT_AB;