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

BUG: sixDoFRigidBodyMotionRestraint. Stabilisation of divisions.

Browse Source
This commit is contained in:
graham
2010-04-28 10:35:13 +01:00
parent ad075512c1
commit 19f9f9b5e4
3 changed files with 5 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyMotion/sixDoFRigidBodyMotionRestraint/linearAxialAngularSpring/linearAxialAngularSpring.C
View File

@ -101,10 +101,10 @@ Foam::sixDoFRigidBodyMotionRestraints::linearAxialAngularSpring::restrain
// Removing any axis component from oldDir and newDir and normalising // Removing any axis component from oldDir and newDir and normalising
oldDir -= (axis_ & oldDir)*axis_; oldDir -= (axis_ & oldDir)*axis_;
oldDir /= mag(oldDir); oldDir /= (mag(oldDir) + VSMALL);
newDir -= (axis_ & newDir)*axis_; newDir -= (axis_ & newDir)*axis_;
newDir /= mag(newDir); newDir /= (mag(newDir) + VSMALL);
scalar theta = mag(acos(min(oldDir & newDir, 1.0))); scalar theta = mag(acos(min(oldDir & newDir, 1.0)));

2
src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyMotion/sixDoFRigidBodyMotionRestraint/linearSpring/linearSpring.C
View File

@ -86,7 +86,7 @@ void Foam::sixDoFRigidBodyMotionRestraints::linearSpring::restrain
scalar magR = mag(r); scalar magR = mag(r);
// r is now the r unit vector // r is now the r unit vector
r /= magR; r /= (magR + VSMALL);
vector v = motion.currentVelocity(restraintPosition); vector v = motion.currentVelocity(restraintPosition);

4
src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyMotion/sixDoFRigidBodyMotionRestraint/tabulatedAxialAngularSpring/tabulatedAxialAngularSpring.C
View File

@ -103,10 +103,10 @@ Foam::sixDoFRigidBodyMotionRestraints::tabulatedAxialAngularSpring::restrain
// Removing any axis component from oldDir and newDir and normalising // Removing any axis component from oldDir and newDir and normalising
oldDir -= (axis_ & oldDir)*axis_; oldDir -= (axis_ & oldDir)*axis_;
oldDir /= mag(oldDir); oldDir /= (mag(oldDir) + VSMALL);
newDir -= (axis_ & newDir)*axis_; newDir -= (axis_ & newDir)*axis_;
newDir /= mag(newDir); newDir /= (mag(newDir) + VSMALL);
scalar theta = mag(acos(min(oldDir & newDir, 1.0))); scalar theta = mag(acos(min(oldDir & newDir, 1.0)));