zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

applications/test/rigidBodyDynamics/pendulum: Simplified using sphere constructor and body lookup by name

Browse Source
This commit is contained in:
Henry Weller
2016-04-04 17:08:20 +01:00
parent 1428f11b6e
commit 19a93cc696
1 changed files with 7 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
applications/test/rigidBodyDynamics/pendulum/pendulum.C
View File

@ -32,6 +32,7 @@ Description
#include "rigidBodyModel.H" #include "rigidBodyModel.H"
#include "masslessBody.H" #include "masslessBody.H"
#include "sphere.H"
#include "joints.H" #include "joints.H"
using namespace Foam; using namespace Foam;
@ -50,9 +51,9 @@ int main(int argc, char *argv[])
// by a hinge which rotates about the z-axis // by a hinge which rotates about the z-axis
if (testMerge) if (testMerge)
{ {
label hingeID = pendulum.join pendulum.join
( (
0, pendulum.bodyID("root"),
Xt(Zero), Xt(Zero),
joint::New(new joints::Rz(pendulum)), joint::New(new joints::Rz(pendulum)),
autoPtr<rigidBody>(new masslessBody("hinge")) autoPtr<rigidBody>(new masslessBody("hinge"))
@ -60,19 +61,19 @@ int main(int argc, char *argv[])
pendulum.merge pendulum.merge
( (
hingeID, pendulum.bodyID("hinge"),
Xt(vector(0, -1, 0)), Xt(vector(0, -1, 0)),
rigidBody::New("weight", 1, Zero, 0.02*I) autoPtr<rigidBody>(new sphere("weight", 1, 0.05))
); );
} }
else else
{ {
pendulum.join pendulum.join
( (
0, pendulum.bodyID("root"),
Xt(Zero), Xt(Zero),
joint::New(new joints::Rz(pendulum)), joint::New(new joints::Rz(pendulum)),
rigidBody::New("pendulum", 1, vector(0, -1, 0), 0.02*I) rigidBody::New("pendulum", 1, vector(0, -1, 0), 1e-3*I)
); );
} }