diff --git a/applications/test/externalCoupler/Make/files b/applications/test/externalCoupler/Make/files
new file mode 100644
index 0000000000..0e55c301e2
--- /dev/null
+++ b/applications/test/externalCoupler/Make/files
@@ -0,0 +1,3 @@
+Test-externalCoupler.C
+
+EXE = $(FOAM_USER_APPBIN)/Test-externalCoupler
diff --git a/applications/test/externalCoupler/Make/options b/applications/test/externalCoupler/Make/options
new file mode 100644
index 0000000000..e2c645cc88
--- /dev/null
+++ b/applications/test/externalCoupler/Make/options
@@ -0,0 +1,7 @@
+EXE_INC = \
+ -I$(LIB_SRC)/finiteVolume/lnInclude \
+ -I$(LIB_SRC)/lumpedPointMotion/lnInclude
+
+EXE_LIBS = \
+ -lfiniteVolume \
+ -llumpedPointMotion
diff --git a/applications/test/externalCoupler/Test-externalCoupler.C b/applications/test/externalCoupler/Test-externalCoupler.C
new file mode 100644
index 0000000000..90cd38d07c
--- /dev/null
+++ b/applications/test/externalCoupler/Test-externalCoupler.C
@@ -0,0 +1,97 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see ("max", 1000);
+
+ externalCoupler coupler;
+
+ if (args.optionFound("slave"))
+ {
+ const word role = "slave";
+ Info<< "Running as " << role << " max=" << maxCount << endl;
+
+ for (label count = 0; count < maxCount; ++count)
+ {
+ // Wait for master, but stop if status=done was seen
+
+ Info<< role << ": waiting for master" << endl;
+ if (!coupler.waitForMaster())
+ {
+ Info<< role << ": stopping. status=done was detected" << endl;
+ break;
+ }
+
+ Info<< role << ": switch to master" << endl;
+ coupler.useMaster();
+ }
+ }
+ else
+ {
+ const word role = "master";
+ Info<< "Running as " << role << " max=" << maxCount << endl;
+
+ for (label count = 0; count < maxCount; ++count)
+ {
+ // Wait for slave, but stop if status=done was seen
+
+ Info<< role << ": waiting for slave" << endl;
+ if (!coupler.waitForSlave())
+ {
+ Info<< role << ": stopping. status=done was detected" << endl;
+ break;
+ }
+
+ Info<< role << ": switch to slave" << endl;
+ coupler.useSlave();
+ }
+ }
+
+ return 0;
+}
+
+
+// ************************************************************************* //
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointForces/Make/files b/applications/utilities/postProcessing/lumped/lumpedPointForces/Make/files
new file mode 100644
index 0000000000..969f66c6f4
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointForces/Make/files
@@ -0,0 +1,3 @@
+lumpedPointForces.C
+
+EXE = $(FOAM_APPBIN)/lumpedPointForces
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointForces/Make/options b/applications/utilities/postProcessing/lumped/lumpedPointForces/Make/options
new file mode 100644
index 0000000000..684551e192
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointForces/Make/options
@@ -0,0 +1,10 @@
+EXE_INC = \
+ -I$(LIB_SRC)/finiteVolume/lnInclude \
+ -I$(LIB_SRC)/fileFormats/lnInclude \
+ -I$(LIB_SRC)/meshTools/lnInclude \
+ -I$(LIB_SRC)/lumpedPointMotion/lnInclude
+
+EXE_LIBS = \
+ -lmeshTools \
+ -lfiniteVolume \
+ -llumpedPointMotion
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointForces/lumpedPointForces.C b/applications/utilities/postProcessing/lumped/lumpedPointForces/lumpedPointForces.C
new file mode 100644
index 0000000000..44ddcdeca3
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointForces/lumpedPointForces.C
@@ -0,0 +1,177 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see
+autoPtr loadField
+(
+ const volMesh::Mesh& mesh,
+ const IOobject* io
+)
+{
+ if (io && io->headerClassName() == GeoFieldType::typeName)
+ {
+ Info<< "Reading " << GeoFieldType::typeName
+ << ' ' << io->name() << endl;
+
+ return autoPtr
+ (
+ new GeoFieldType
+ (
+ IOobject
+ (
+ io->name(),
+ io->instance(),
+ io->local(),
+ io->db(),
+ IOobject::MUST_READ,
+ IOobject::AUTO_WRITE,
+ io->registerObject()
+ ),
+ mesh
+ )
+ );
+ }
+
+ return autoPtr();
+}
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+int main(int argc, char *argv[])
+{
+ argList::addNote
+ (
+ "Extract force/moment information from existing calculations based "
+ "on the lumpedPoints pressure zones."
+ );
+
+ argList::addBoolOption
+ (
+ "vtk",
+ "create visualization files of the forces"
+ );
+
+ timeSelector::addOptions(true, false);
+ argList::noFunctionObjects(); // Never use function objects
+
+ #include "addRegionOption.H"
+ #include "setRootCase.H"
+
+ const bool withVTK = args.optionFound("vtk");
+
+ #include "createTime.H"
+ instantList timeDirs = timeSelector::select0(runTime, args);
+
+ #include "createNamedMesh.H"
+
+ autoPtr movement = lumpedPointIOMovement::New
+ (
+ runTime
+ );
+
+ if (!movement.valid())
+ {
+ Info<< "no valid movement given" << endl;
+ return 1;
+ }
+
+ const labelList patchLst = lumpedPointTools::lumpedPointPatchList(mesh);
+ if (patchLst.empty())
+ {
+ Info<< "no patch list found" << endl;
+ return 2;
+ }
+
+ movement().setMapping(mesh, patchLst, lumpedPointTools::points0Field(mesh));
+
+ List forces, moments;
+ forAll(timeDirs, timei)
+ {
+ runTime.setTime(timeDirs[timei], timei);
+
+ Info<< "Time = " << runTime.timeName() << endl;
+
+ if (mesh.readUpdate())
+ {
+ Info<< " Read new mesh" << nl;
+ }
+
+ // Search for list of objects for this time
+ IOobjectList objects(mesh, runTime.timeName());
+
+ // Pressure field
+ autoPtr pField
+ = loadField(mesh, objects.lookup("p"));
+
+ // The forces per zone
+ if (movement().forcesAndMoments(mesh, forces, moments))
+ {
+ Info<<"forces per zone: " << forces << endl;
+ Info<<"moments per zone: " << moments << endl;
+
+ if (withVTK && Pstream::master())
+ {
+ const word outputName =
+ Foam::name("forces_%06d.vtp", runTime.timeIndex());
+
+ Info<<" " << outputName << endl;
+
+ movement().writeForcesAndMomentsVTP
+ (
+ outputName,
+ forces,
+ moments
+ );
+ }
+ }
+ }
+
+ Info<< "End\n" << endl;
+
+ return 0;
+}
+
+// ************************************************************************* //
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointMovement/Make/files b/applications/utilities/postProcessing/lumped/lumpedPointMovement/Make/files
new file mode 100644
index 0000000000..0efb02170d
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointMovement/Make/files
@@ -0,0 +1,3 @@
+lumpedPointMovement.C
+
+EXE = $(FOAM_APPBIN)/lumpedPointMovement
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointMovement/Make/options b/applications/utilities/postProcessing/lumped/lumpedPointMovement/Make/options
new file mode 100644
index 0000000000..684551e192
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointMovement/Make/options
@@ -0,0 +1,10 @@
+EXE_INC = \
+ -I$(LIB_SRC)/finiteVolume/lnInclude \
+ -I$(LIB_SRC)/fileFormats/lnInclude \
+ -I$(LIB_SRC)/meshTools/lnInclude \
+ -I$(LIB_SRC)/lumpedPointMotion/lnInclude
+
+EXE_LIBS = \
+ -lmeshTools \
+ -lfiniteVolume \
+ -llumpedPointMotion
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointMovement/lumpedPointMovement.C b/applications/utilities/postProcessing/lumped/lumpedPointMovement/lumpedPointMovement.C
new file mode 100644
index 0000000000..fdbc810983
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointMovement/lumpedPointMovement.C
@@ -0,0 +1,263 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see ("max", 0));
+
+ const label span =
+ Foam::max(1, args.optionLookupOrDefault("span", 1));
+
+ const scalar relax = args.optionLookupOrDefault("scale", 1);
+
+ const bool slave = args.optionFound("slave");
+ const bool removeLock = args.optionFound("removeLock");
+
+ #include "createTime.H"
+
+ autoPtr movement = lumpedPointIOMovement::New
+ (
+ runTime
+ );
+
+ if (!movement.valid())
+ {
+ Info<< "no valid movement given" << endl;
+ return 1;
+ }
+
+ List responseTable =
+ lumpedPointTools::lumpedPointStates(args[1]);
+
+ Info<< "Using response table with " << responseTable.size()
+ << " entries" << endl;
+
+ Info << "Increment input by " << span << nl;
+
+ if (maxOut)
+ {
+ Info<< "Stopping after " << maxOut << " outputs" << endl;
+ }
+
+ if (slave)
+ {
+ Info<< "Running as slave responder" << endl;
+
+ externalCoupler& coupler = movement().coupler();
+
+ label count = 0;
+ for (label index = 0; index < responseTable.size(); index += span)
+ {
+ Info<< args.executable() << ": waiting for master" << endl;
+
+ // Wait for master, but stop if status=done was seen
+ if (!coupler.waitForMaster())
+ {
+ Info<< args.executable()
+ << ": stopping status=done was detected" << endl;
+ break;
+ }
+
+ lumpedPointState state = responseTable[index].second();
+ state.relax(relax, movement().state0());
+
+ // Generate input for OpenFOAM
+ OFstream os(coupler.resolveFile(movement().inputName()));
+ if
+ (
+ movement().inputFormat()
+ == lumpedPointState::inputFormatType::PLAIN
+ )
+ {
+ state.writePlain(os);
+ }
+ else
+ {
+ os.writeEntry("time", responseTable[index].first());
+ state.writeDict(os);
+ }
+
+ Info<< args.executable()
+ << ": updated to state " << index
+ << " - switch to master"
+ << endl;
+
+ // Let OpenFOAM know that it can continue
+ coupler.useMaster();
+
+ if (maxOut && ++count >= maxOut)
+ {
+ Info<< args.executable()
+ <<": stopping after " << maxOut << " outputs" << endl;
+ break;
+ }
+ }
+
+ if (removeLock)
+ {
+ Info<< args.executable() <<": removing lock file" << endl;
+ coupler.useSlave(); // This removes the lock-file
+ }
+ }
+ else
+ {
+ runTime.setTime(instant(0, runTime.constant()), 0);
+
+ #include "createNamedPolyMesh.H"
+
+ const labelList patchLst = lumpedPointTools::lumpedPointPatchList(mesh);
+ if (patchLst.empty())
+ {
+ Info<< "no patch list found" << endl;
+ return 2;
+ }
+
+ pointIOField points0 = lumpedPointTools::points0Field(mesh);
+ movement().setBoundBox(mesh, patchLst, points0);
+
+ label index = 0;
+
+ // Initial geometry
+ movement().writeVTP("geom_init.vtp", mesh, patchLst, points0);
+
+ forAll(responseTable, i)
+ {
+ const bool output = ((i % span) == 0);
+ lumpedPointState state = responseTable[i].second();
+ state.relax(relax, movement().state0());
+
+ if (output)
+ {
+ Info<<"output [" << i << "/"
+ << responseTable.size() << "]" << endl;
+ }
+ else
+ {
+ continue;
+ }
+
+ // State/response = what comes back from FEM
+ {
+ const word outputName = Foam::name("state_%06d.vtp", index);
+ Info<<" " << outputName << endl;
+
+ state.writeVTP(outputName, movement().axis());
+ }
+
+ {
+ const word outputName = Foam::name("geom_%06d.vtp", index);
+ Info<<" " << outputName << endl;
+
+ movement().writeVTP(outputName, state, mesh, patchLst, points0);
+ }
+
+ {
+ ++index;
+
+ bool canOutput = !maxOut || (index <= maxOut);
+ if (!canOutput)
+ {
+ Info<<"stopping output after "
+ << maxOut << " outputs" << endl;
+ break;
+ }
+ }
+ }
+
+ }
+
+ Info<< args.executable() << ": End\n" << endl;
+
+ return 0;
+}
+
+
+// ************************************************************************* //
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointZones/Make/files b/applications/utilities/postProcessing/lumped/lumpedPointZones/Make/files
new file mode 100644
index 0000000000..a9c8e8f9f4
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointZones/Make/files
@@ -0,0 +1,3 @@
+lumpedPointZones.C
+
+EXE = $(FOAM_APPBIN)/lumpedPointZones
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointZones/Make/options b/applications/utilities/postProcessing/lumped/lumpedPointZones/Make/options
new file mode 100644
index 0000000000..684551e192
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointZones/Make/options
@@ -0,0 +1,10 @@
+EXE_INC = \
+ -I$(LIB_SRC)/finiteVolume/lnInclude \
+ -I$(LIB_SRC)/fileFormats/lnInclude \
+ -I$(LIB_SRC)/meshTools/lnInclude \
+ -I$(LIB_SRC)/lumpedPointMotion/lnInclude
+
+EXE_LIBS = \
+ -lmeshTools \
+ -lfiniteVolume \
+ -llumpedPointMotion
diff --git a/applications/utilities/postProcessing/lumped/lumpedPointZones/lumpedPointZones.C b/applications/utilities/postProcessing/lumped/lumpedPointZones/lumpedPointZones.C
new file mode 100644
index 0000000000..3103b01422
--- /dev/null
+++ b/applications/utilities/postProcessing/lumped/lumpedPointZones/lumpedPointZones.C
@@ -0,0 +1,107 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see movement = lumpedPointIOMovement::New
+ (
+ runTime
+ );
+
+ if (!movement.valid())
+ {
+ Info<< "no valid movement found" << endl;
+ return 1;
+ }
+
+ const labelList patchLst = lumpedPointTools::lumpedPointPatchList(mesh);
+ if (patchLst.empty())
+ {
+ Info<< "no patch list found" << endl;
+ return 2;
+ }
+
+ pointIOField points0 = lumpedPointTools::points0Field(mesh);
+ movement().setMapping(mesh, patchLst, points0);
+
+ // Initial geometry, but with zone colouring
+ movement().writeZonesVTP("lumpedPointZones.vtp", mesh, points0);
+
+ // Initial positions/rotations
+ movement().writeStateVTP("initialState.vtp");
+
+ Info<< nl
+ << "wrote 'lumpedPointZones.vtp'" << nl
+ << "wrote 'initialState.vtp'" << nl
+ << "End\n" << endl;
+
+ return 0;
+}
+
+// ************************************************************************* //
diff --git a/src/Allwmake b/src/Allwmake
index abe861514e..d897043005 100755
--- a/src/Allwmake
+++ b/src/Allwmake
@@ -76,6 +76,7 @@ wmake $targetType regionCoupled
functionObjects/Allwmake $targetType $*
+wmake $targetType lumpedPointMotion
wmake $targetType sixDoFRigidBodyMotion
wmake $targetType rigidBodyDynamics
wmake $targetType rigidBodyMeshMotion
diff --git a/src/lumpedPointMotion/Make/files b/src/lumpedPointMotion/Make/files
new file mode 100644
index 0000000000..b525c884b8
--- /dev/null
+++ b/src/lumpedPointMotion/Make/files
@@ -0,0 +1,13 @@
+externalCoupler.C
+
+lumpedPointMovement.C
+lumpedPointMovementWriter.C
+lumpedPointState.C
+lumpedPointStateWriter.C
+lumpedPointIOMovement.C
+
+lumpedPointDisplacementPointPatchVectorField.C
+
+lumpedPointTools.C
+
+LIB = $(FOAM_LIBBIN)/liblumpedPointMotion
diff --git a/src/lumpedPointMotion/Make/options b/src/lumpedPointMotion/Make/options
new file mode 100644
index 0000000000..9aa6ebf53d
--- /dev/null
+++ b/src/lumpedPointMotion/Make/options
@@ -0,0 +1,10 @@
+EXE_INC = \
+ -I$(LIB_SRC)/finiteVolume/lnInclude \
+ -I$(LIB_SRC)/dynamicMesh/lnInclude \
+ -I$(LIB_SRC)/fileFormats/lnInclude \
+ -I$(LIB_SRC)/meshTools/lnInclude
+
+LIB_LIBS = \
+ -lfiniteVolume \
+ -ldynamicMesh \
+ -lmeshTools
diff --git a/src/lumpedPointMotion/externalCoupler.C b/src/lumpedPointMotion/externalCoupler.C
new file mode 100644
index 0000000000..d5f5a17464
--- /dev/null
+++ b/src/lumpedPointMotion/externalCoupler.C
@@ -0,0 +1,309 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2015-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see ("initByExternal", false);
+
+ waitInterval_ = dict.lookupOrDefault("waitInterval", 1u);
+ if (!waitInterval_)
+ {
+ // Enforce non-zero sleep
+ waitInterval_ = 1u;
+ }
+
+ timeOut_ = dict.lookupOrDefault("timeOut", 100*waitInterval_);
+
+ log = dict.lookupOrDefault("log", false);
+
+ return true;
+}
+
+
+bool Foam::externalCoupler::initialized() const
+{
+ return runState_ != NONE;
+}
+
+
+bool Foam::externalCoupler::slaveFirst() const
+{
+ return slaveFirst_;
+}
+
+
+void Foam::externalCoupler::useMaster(const bool wait) const
+{
+ const bool wasInit = initialized();
+ runState_ = MASTER;
+
+ if (Pstream::master())
+ {
+ if (!wasInit)
+ {
+ // First time
+ Foam::mkDir(commsDir_);
+ }
+
+ const fileName lck(lockFile());
+
+ // Create lock file - only if it doesn't already exist
+ if (!Foam::isFile(lck))
+ {
+ Log << type() << ": creating lock file" << endl;
+
+ OFstream os(lck);
+ os << "status=openfoam\n";
+ os.flush();
+ }
+ }
+
+ if (wait)
+ {
+ waitForMaster();
+ }
+}
+
+
+void Foam::externalCoupler::useSlave(const bool wait) const
+{
+ const bool wasInit = initialized();
+ runState_ = SLAVE;
+
+ if (Pstream::master())
+ {
+ if (!wasInit)
+ {
+ // First time
+ Foam::mkDir(commsDir_);
+ }
+
+ Log << type() << ": removing lock file" << endl;
+
+ Foam::rm(lockFile());
+ }
+
+ if (wait)
+ {
+ waitForSlave();
+ }
+}
+
+
+bool Foam::externalCoupler::waitForMaster() const
+{
+ if (!initialized())
+ {
+ useMaster(); // was not initialized
+ }
+
+ bool isDone = false;
+ if (Pstream::master())
+ {
+ const fileName lck(lockFile());
+
+ double prevTime = -1;
+ double modTime = 0;
+
+ // Wait until file disappears (modTime == 0)
+ // But also check for status=done content in the file
+ while ((modTime = highResLastModified(lck)) > 0)
+ {
+ if (prevTime < modTime)
+ {
+ prevTime = modTime;
+ isDone = checkIsDone(lck);
+ if (isDone)
+ {
+ break;
+ }
+ }
+ sleep(waitInterval_);
+ }
+ }
+
+ // MPI barrier
+ Pstream::scatter(isDone);
+
+ return !isDone;
+}
+
+
+bool Foam::externalCoupler::waitForSlave() const
+{
+ if (!initialized())
+ {
+ useSlave(); // was not initialized
+ }
+
+ bool isDone = false;
+ if (Pstream::master())
+ {
+ const fileName lck(lockFile());
+ unsigned totalTime = 0;
+
+ Log << type() << ": waiting for lock file to appear " << lck << endl;
+
+ while (!Foam::isFile(lck))
+ {
+ sleep(waitInterval_);
+
+ if (timeOut_ && (totalTime += waitInterval_) > timeOut_)
+ {
+ FatalErrorInFunction
+ << "Wait time exceeded timeout of " << timeOut_
+ << " s" << abort(FatalError);
+ }
+
+ Log << type() << ": wait time = " << totalTime << endl;
+ }
+
+ // But check for status=done content in the file
+ isDone = checkIsDone(lck);
+
+ Log << type() << ": found lock file " << lck << endl;
+ }
+
+ // MPI barrier
+ Pstream::scatter(isDone);
+
+ return !isDone;
+}
+
+
+Foam::fileName Foam::externalCoupler::resolveFile
+(
+ const word& file
+) const
+{
+ return fileName(baseDir()/file);
+}
+
+
+void Foam::externalCoupler::shutdown() const
+{
+ if (Pstream::master() && runState_ == MASTER && Foam::isDir(commsDir_))
+ {
+ const fileName lck(lockFile());
+
+ Log << type() << ": lock file status=done" << endl;
+ OFstream os(lck);
+ os << "status=done\n";
+ os.flush();
+ }
+
+ runState_ = DONE; // Avoid re-triggering in destructor
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/externalCoupler.H b/src/lumpedPointMotion/externalCoupler.H
new file mode 100644
index 0000000000..567bb75adf
--- /dev/null
+++ b/src/lumpedPointMotion/externalCoupler.H
@@ -0,0 +1,207 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2015-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify i
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see patchLst(bf.size());
+ forAll(bf, patchi)
+ {
+ // All patches of this type
+ if (isA(bf[patchi]))
+ {
+ patchLst.append(patchi);
+ // or patchLst.append(bf[patchi].patch().index());
+ }
+ }
+
+ return patchLst.shrink();
+}
+
+
+// * * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * //
+
+const Foam::pointField&
+Foam::lumpedPointDisplacementPointPatchVectorField::points0() const
+{
+ const objectRegistry& obr = this->patch().boundaryMesh().mesh().db();
+
+ // Obtain starting locations from the motionSolver
+ return obr.lookupObject
+ (
+ "dynamicMeshDict"
+ ).points0();
+}
+
+
+const Foam::lumpedPointMovement&
+Foam::lumpedPointDisplacementPointPatchVectorField::movement() const
+{
+ const objectRegistry& obr = this->patch().boundaryMesh().mesh().db();
+ const lumpedPointIOMovement* ptr =
+ lumpedPointIOMovement::lookupInRegistry(obr);
+
+ if (ptr)
+ {
+ return *ptr; // Already exists
+ }
+
+ // create and register with this patch as the owner
+ autoPtr obj = lumpedPointIOMovement::New
+ (
+ obr,
+ this->patch().index()
+ );
+
+ return objectRegistry::store(obj);
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
+
+Foam::lumpedPointDisplacementPointPatchVectorField::
+lumpedPointDisplacementPointPatchVectorField
+(
+ const pointPatch& p,
+ const DimensionedField& iF
+)
+:
+ fixedValuePointPatchField(p, iF)
+{}
+
+
+Foam::lumpedPointDisplacementPointPatchVectorField::
+lumpedPointDisplacementPointPatchVectorField
+(
+ const pointPatch& p,
+ const DimensionedField& iF,
+ const dictionary& dict
+)
+:
+ fixedValuePointPatchField(p, iF, dict)
+{}
+
+
+Foam::lumpedPointDisplacementPointPatchVectorField::
+lumpedPointDisplacementPointPatchVectorField
+(
+ const lumpedPointDisplacementPointPatchVectorField& pf,
+ const pointPatch& p,
+ const DimensionedField& iF,
+ const pointPatchFieldMapper& mapper
+)
+:
+ fixedValuePointPatchField(pf, p, iF, mapper)
+{}
+
+
+Foam::lumpedPointDisplacementPointPatchVectorField::
+lumpedPointDisplacementPointPatchVectorField
+(
+ const lumpedPointDisplacementPointPatchVectorField& pf,
+ const DimensionedField& iF
+)
+:
+ fixedValuePointPatchField(pf, iF)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
+
+Foam::lumpedPointDisplacementPointPatchVectorField::
+~lumpedPointDisplacementPointPatchVectorField()
+{
+ // de-register movement if in use and managed by this patch
+ const lumpedPointIOMovement* ptr = lumpedPointIOMovement::lookupInRegistry
+ (
+ this->patch().boundaryMesh().mesh().db()
+ );
+
+ if (ptr && ptr->ownerId() == this->patch().index())
+ {
+ movement().coupler().shutdown();
+
+ const_cast(ptr)->checkOut();
+ }
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
+
+void Foam::lumpedPointDisplacementPointPatchVectorField::updateCoeffs()
+{
+ if (this->updated())
+ {
+ return;
+ }
+
+ const bool masterPatch = (movement().ownerId() == this->patch().index());
+ if (masterPatch)
+ {
+ if (lumpedPointIOMovement::debug)
+ {
+ Pout<<"masterPatch: " << this->patch().index() << endl;
+ }
+
+ const polyMesh& mesh = this->patch().boundaryMesh().mesh().mesh();
+
+ // need face 'zones' for calculating forces
+ // likely need bounding box for the movement
+ // -> do both now if required
+ if (!movement().hasMapping())
+ {
+ const_cast(movement()).setMapping
+ (
+ mesh,
+ // All patches of this type
+ patchIds
+ (
+ static_cast
+ (
+ this->internalField()
+ )
+ ),
+ this->points0()
+ );
+ }
+
+
+ if
+ (
+ movement().coupler().initialized()
+ || !movement().coupler().slaveFirst()
+ )
+ {
+ // Synchronized for all processes
+ List forces, moments;
+ movement().forcesAndMoments(mesh, forces, moments);
+
+ if (lumpedPointIOMovement::debug)
+ {
+ Pout<<"gatherForces: " << forces << " called from patch "
+ << this->patch().index() << endl;
+
+ if (Pstream::master())
+ {
+ Pout<<"output forces to file: "
+ << movement().locations() << " called from patch "
+ << this->patch().index() << nl
+ <<"# " << forces.size() << " force entries" << nl
+ <<"# fx fy fz" << nl
+ <<"output forces to file: "
+ << forces << " called from patch "
+ << this->patch().index() << endl;
+ }
+ }
+
+ if (Pstream::master())
+ {
+ movement().writeData(forces, moments);
+
+ // signal external source to execute
+ movement().coupler().useSlave();
+ }
+ }
+
+ // Wait for slave to provide data - includes MPI barrier
+ movement().coupler().waitForSlave();
+
+ // Read data passed back from external source - includes MPI barrier
+ const_cast(movement()).readState();
+ }
+
+ tmp tdisp = movement().displacePoints
+ (
+ this->points0(),
+ this->patch().meshPoints()
+ );
+
+ this->operator==(tdisp);
+
+ fixedValuePointPatchField::updateCoeffs();
+}
+
+
+void Foam::lumpedPointDisplacementPointPatchVectorField::write(Ostream& os)
+const
+{
+ pointPatchField::write(os);
+ writeEntry("value", os);
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointDisplacementPointPatchVectorField.H b/src/lumpedPointMotion/lumpedPointDisplacementPointPatchVectorField.H
new file mode 100644
index 0000000000..6e3a59558a
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointDisplacementPointPatchVectorField.H
@@ -0,0 +1,200 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see
+{
+ // Private data
+
+ //- Convenience typedefs
+ typedef lumpedPointDisplacementPointPatchVectorField patchType;
+ typedef DimensionedField fieldType;
+
+ // Private Member Functions
+
+protected:
+
+ //- The starting locations (obtained from the motionSolver).
+ const pointField& points0() const;
+
+ //- The auto-vivifying singleton for movement.
+ const lumpedPointMovement& movement() const;
+
+
+public:
+
+ //- Runtime type information
+ TypeName("lumpedPointDisplacement");
+
+
+ // Constructors
+
+ //- Construct from patch and internal field
+ lumpedPointDisplacementPointPatchVectorField
+ (
+ const pointPatch& p,
+ const DimensionedField& iF
+ );
+
+ //- Construct from patch, internal field and dictionary
+ lumpedPointDisplacementPointPatchVectorField
+ (
+ const pointPatch& p,
+ const DimensionedField& iF,
+ const dictionary& dict
+ );
+
+ //- Construct by mapping given patchField onto a new patch
+ lumpedPointDisplacementPointPatchVectorField
+ (
+ const lumpedPointDisplacementPointPatchVectorField& pf,
+ const pointPatch& p,
+ const DimensionedField& iF,
+ const pointPatchFieldMapper& mapper
+ );
+
+ //- Construct and return a clone
+ virtual autoPtr> clone() const
+ {
+ return autoPtr>
+ (
+ new lumpedPointDisplacementPointPatchVectorField
+ (
+ *this
+ )
+ );
+ }
+
+ //- Construct as copy setting internal field reference
+ lumpedPointDisplacementPointPatchVectorField
+ (
+ const lumpedPointDisplacementPointPatchVectorField& pf,
+ const DimensionedField& iF
+ );
+
+ //- Construct and return a clone setting internal field reference
+ virtual autoPtr> clone
+ (
+ const DimensionedField& iF
+ ) const
+ {
+ return autoPtr>
+ (
+ new lumpedPointDisplacementPointPatchVectorField
+ (
+ *this,
+ iF
+ )
+ );
+ }
+
+ //- Destructor
+ virtual ~lumpedPointDisplacementPointPatchVectorField();
+
+
+ // Member functions
+
+ //- The ids for all patches of this type
+ static labelList patchIds(const pointVectorField& pvf);
+
+ //- Update the coefficients associated with the patch field
+ virtual void updateCoeffs();
+
+ //- Write
+ virtual void write(Ostream&) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointIOMovement.C b/src/lumpedPointMotion/lumpedPointIOMovement.C
new file mode 100644
index 0000000000..dbe1fde102
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointIOMovement.C
@@ -0,0 +1,142 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see
+ (
+ lumpedPointMovement::dictionaryName
+ );
+}
+
+
+Foam::autoPtr
+Foam::lumpedPointIOMovement::New
+(
+ const objectRegistry& obr,
+ label ownerId
+)
+{
+ return autoPtr
+ (
+ new lumpedPointIOMovement
+ (
+ IOobject
+ (
+ lumpedPointMovement::dictionaryName,
+ obr.time().caseSystem(),
+ obr,
+ IOobject::MUST_READ,
+ IOobject::NO_WRITE,
+ true // register object
+ ),
+ ownerId // tag this patch as owner too
+ )
+ );
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
+
+Foam::lumpedPointIOMovement::lumpedPointIOMovement
+(
+ const IOobject& io,
+ label ownerId
+)
+:
+ lumpedPointMovement(),
+ regIOobject(io)
+{
+ bool ok =
+ (
+ readOpt() == IOobject::MUST_READ
+ || readOpt() == IOobject::MUST_READ_IF_MODIFIED
+ );
+
+ if (ok)
+ {
+ ok = readData(readStream(typeName));
+ close();
+
+ if (ok)
+ {
+ this->ownerId(ownerId);
+ }
+ }
+}
+
+
+// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
+
+Foam::lumpedPointIOMovement::~lumpedPointIOMovement()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
+
+
+bool Foam::lumpedPointIOMovement::readData(Istream& is)
+{
+ dictionary dict(is);
+
+ readDict(dict);
+
+ return is.check(FUNCTION_NAME);
+}
+
+
+bool Foam::lumpedPointIOMovement::writeData(Ostream& os) const
+{
+ os << *this;
+ return os.good();
+}
+
+
+// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
+
+Foam::Ostream& Foam::operator<<(Ostream& os, const lumpedPointIOMovement& obj)
+{
+ obj.writeDict(os);
+
+ os.check(FUNCTION_NAME);
+ return os;
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointIOMovement.H b/src/lumpedPointMotion/lumpedPointIOMovement.H
new file mode 100644
index 0000000000..f1f121b74b
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointIOMovement.H
@@ -0,0 +1,133 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see New
+ (
+ const objectRegistry& obr,
+ label ownerId = -1
+ );
+
+
+ // Constructors
+
+ //- Construct from IOobject, optionally with some owner information
+ explicit lumpedPointIOMovement
+ (
+ const IOobject& io,
+ label ownerId = -1
+ );
+
+
+ //- Destructor
+ ~lumpedPointIOMovement();
+
+
+ // Member Functions
+
+ //- readData member function used by regIOobject
+ bool readData(Istream& is);
+
+ //- writeData member function required by regIOobject
+ bool writeData(Ostream& os) const;
+
+
+ // IOstream Operators
+
+ friend Ostream& operator<<
+ (
+ Ostream& os,
+ const lumpedPointIOMovement& obj
+ );
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointMovement.C b/src/lumpedPointMotion/lumpedPointMovement.C
new file mode 100644
index 0000000000..b12a3ef941
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointMovement.C
@@ -0,0 +1,790 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see
+ Foam::lumpedPointMovement::formatNames
+ {
+ { outputFormatType::PLAIN, "plain" },
+ { outputFormatType::DICTIONARY, "dictionary" }
+ };
+
+
+const Foam::word
+Foam::lumpedPointMovement::dictionaryName("lumpedPointMovement");
+
+
+namespace Foam
+{
+//! \cond fileScope
+//- Write list content with size, bracket, content, bracket one-per-line.
+// This makes for consistent for parsing, regardless of the list length.
+template
+static void writeList(Ostream& os, const string& header, const UList& lst)
+{
+ // Header string
+ os << header.c_str() << nl;
+
+ // Write size and start delimiter
+ os << lst.size() << nl
+ << token::BEGIN_LIST << nl;
+
+ // Write contents
+ forAll(lst, i)
+ {
+ os << lst[i] << nl;
+ }
+
+ // Write end delimiter
+ os << token::END_LIST << token::END_STATEMENT << nl << nl;
+}
+//! \endcond
+
+}
+// namespace Foam
+
+
+// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
+
+Foam::IOobject Foam::lumpedPointMovement::selectIO
+(
+ const IOobject& io,
+ const fileName& f
+)
+{
+ return
+ (
+ f.size()
+ ? IOobject // construct from filePath instead
+ (
+ f,
+ io.db(),
+ io.readOpt(),
+ io.writeOpt(),
+ io.registerObject(),
+ io.globalObject()
+ )
+ : io
+ );
+}
+
+
+// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
+
+void Foam::lumpedPointMovement::calcThresholds() const
+{
+ thresholdPtr_ = new scalarField(locations_);
+
+ scalarField& thresh = *thresholdPtr_;
+
+ for (label i=1; i < thresh.size(); ++i)
+ {
+ thresh[i-1] =
+ (
+ locations_[i-1]
+ + (division_ * (locations_[i] - locations_[i-1]))
+ );
+ }
+}
+
+
+Foam::label Foam::lumpedPointMovement::threshold(scalar pos) const
+{
+ const scalarField& thresh = this->thresholds();
+ const label n = thresh.size();
+
+ // Clamp above and below
+ //
+ // ? may need special treatment to clamp values below the first threshold ?
+ // ? special treatment when 'axis' is negative ?
+
+ for (label i=0; i < n; ++i)
+ {
+ if (pos < thresh[i])
+ {
+ return i;
+ }
+ }
+
+ // everything above goes into the last zone too
+ return n-1;
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
+
+Foam::lumpedPointMovement::lumpedPointMovement()
+:
+ axis_(0,0,1),
+ locations_(),
+ division_(0),
+ relax_(1),
+ interpolationScheme_(linearInterpolationWeights::typeName),
+ ownerId_(-1),
+ boundBox_(),
+ centre_(),
+ autoCentre_(true),
+ forcesDict_(),
+ coupler_(),
+ inputName_("positions.in"),
+ outputName_("forces.out"),
+ inputFormat_(lumpedPointState::inputFormatType::DICTIONARY),
+ outputFormat_(outputFormatType::DICTIONARY),
+ state0_(),
+ state_(),
+ thresholdPtr_(0),
+ interpolatorPtr_()
+{}
+
+
+Foam::lumpedPointMovement::lumpedPointMovement
+(
+ const dictionary& dict,
+ label ownerId
+)
+:
+ axis_(0,0,1),
+ locations_(),
+ division_(0),
+ relax_(1),
+ interpolationScheme_
+ (
+ dict.lookupOrDefault
+ (
+ "interpolationScheme",
+ linearInterpolationWeights::typeName
+ )
+ ),
+ ownerId_(ownerId),
+ boundBox_(),
+ centre_(),
+ autoCentre_(true),
+ forcesDict_(),
+ coupler_(),
+ state0_(),
+ state_(),
+ thresholdPtr_(0),
+ interpolatorPtr_()
+{
+ readDict(dict);
+}
+
+
+// * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * * //
+
+Foam::lumpedPointMovement::~lumpedPointMovement()
+{
+ deleteDemandDrivenData(thresholdPtr_);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
+
+void Foam::lumpedPointMovement::readDict(const dictionary& dict)
+{
+ // assume the worst
+ deleteDemandDrivenData(thresholdPtr_);
+
+ dict.lookup("axis") >> axis_;
+
+ division_ = 0;
+ if (dict.readIfPresent("division", division_))
+ {
+ if (division_ < 0)
+ {
+ division_ = 0;
+ }
+ else if (division_ > 1)
+ {
+ division_ = 1;
+ }
+ }
+
+ dict.readIfPresent("relax", relax_);
+
+ dict.lookup("locations") >> locations_;
+
+ if (dict.readIfPresent("interpolationScheme", interpolationScheme_))
+ {
+ interpolatorPtr_.clear();
+ }
+
+ autoCentre_ = !dict.readIfPresent("centre", centre_);
+
+ // Initial state from initial locations, with zero rotation
+ tmp pts = locations_ * axis_;
+ state0_ = lumpedPointState(pts);
+ state_ = state0_;
+
+ forcesDict_.merge(dict.subOrEmptyDict("forces"));
+
+ const dictionary& commDict = dict.subDict("communication");
+ coupler_.readDict(commDict);
+
+ // TODO: calcFrequency_ = dict.lookupOrDefault("calcFrequency", 1);
+
+ commDict.lookup("inputName") >> inputName_;
+ commDict.lookup("outputName") >> outputName_;
+
+ inputFormat_ = lumpedPointState::formatNames.lookup
+ (
+ "inputFormat",
+ commDict
+ );
+
+ outputFormat_ = lumpedPointMovement::formatNames.lookup
+ (
+ "outputFormat",
+ commDict
+ );
+}
+
+
+void Foam::lumpedPointMovement::setBoundBox
+(
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+)
+{
+ boundBox_ = boundBox::invertedBox;
+
+ const polyBoundaryMesh& patches = mesh.boundaryMesh();
+ for (const label patchId : patchIds)
+ {
+ boundBox_.add(points0, patches[patchId].meshPoints());
+ }
+
+ boundBox_.reduce();
+
+ if (autoCentre_)
+ {
+ centre_ = boundBox_.midpoint();
+ centre_ -= (centre_ & axis_) * axis_;
+ if (lumpedPointIOMovement::debug)
+ {
+ Pout<<"autoCentre on " << centre_ << " boundBox: " << boundBox_ << endl;
+ }
+ }
+ else
+ {
+ // User-specified centre
+ if (lumpedPointIOMovement::debug)
+ {
+ Pout<<"centre on " << centre_ << " boundBox: " << boundBox_ << endl;
+ }
+ }
+}
+
+
+void Foam::lumpedPointMovement::setMapping
+(
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+)
+{
+ setBoundBox(mesh, patchIds, points0);
+
+ faceZones_.clear();
+
+ // Local storage location (boundary faces only)
+ const label firstFace = mesh.nInternalFaces();
+ labelList faceToZoneID(mesh.nFaces() - firstFace, -1);
+
+ // Number of faces per zone
+ labelList nFaces(thresholds().size(), 0);
+
+ const polyBoundaryMesh& patches = mesh.boundaryMesh();
+
+ for (const label patchId : patchIds)
+ {
+ const polyPatch& pp = patches[patchId];
+
+ // Classify faces into respective zones according to their centres
+ label meshFaceI = pp.start();
+
+ forAll(pp, i)
+ {
+ const face& f = mesh.faces()[meshFaceI];
+ label zoneId = this->threshold(f.centre(points0));
+
+ // localize storage location
+ faceToZoneID[meshFaceI - firstFace] = zoneId;
+ ++nFaces[zoneId];
+
+ ++meshFaceI;
+ }
+ }
+
+ if (lumpedPointIOMovement::debug)
+ {
+ Pout<<"faces per zone:" << nFaces << endl;
+ }
+
+ faceZones_.setSize(nFaces.size());
+
+ label nZonedFaces = 0;
+ forAll(nFaces, zonei)
+ {
+ label n = nFaces[zonei];
+ labelList& addr = faceZones_[zonei];
+ addr.setSize(n);
+
+ n = 0;
+ forAll(faceToZoneID, faceI)
+ {
+ if (faceToZoneID[faceI] == zonei)
+ {
+ // from local storage to global mesh face
+ label meshFaceI = faceI + firstFace;
+
+ addr[n] = meshFaceI;
+ ++n;
+ }
+ }
+
+ addr.setSize(n);
+ nZonedFaces += n;
+
+ if (lumpedPointIOMovement::debug)
+ {
+ Pout<< "Created faceZone[" << zonei << "] "
+ << returnReduce(n, sumOp()) << " faces, "
+ << thresholds()[zonei] << " threshold" << nl;
+ }
+ }
+}
+
+
+bool Foam::lumpedPointMovement::forcesAndMoments
+(
+ const polyMesh& pmesh,
+ List& forces,
+ List& moments
+) const
+{
+ forces.setSize(faceZones_.size());
+ moments.setSize(faceZones_.size());
+
+ if (faceZones_.empty())
+ {
+ WarningInFunction
+ << "Attempted to calculate forces without setMapping()"
+ << nl << endl;
+
+ forces.setSize(thresholds().size(), Zero);
+ moments.setSize(thresholds().size(), Zero);
+ return false;
+ }
+
+ // Initialize with zero
+ forces = Zero;
+ moments = Zero;
+
+ // The mass centres
+ const pointField& lumpedCentres = state().points();
+
+ const polyBoundaryMesh& patches = pmesh.boundaryMesh();
+
+ const word pName = forcesDict_.lookupOrDefault("p", "p");
+
+ scalar pRef = forcesDict_.lookupOrDefault("pRef", 0.0);
+ scalar rhoRef = forcesDict_.lookupOrDefault("rhoRef", 1.0);
+
+
+ // Calculated force per patch - cache
+ PtrMap forceOnPatches;
+
+ const volScalarField* pPtr = pmesh.lookupObjectPtr(pName);
+
+ // fvMesh and has pressure field
+ if (isA(pmesh) && pPtr)
+ {
+ const fvMesh& mesh = dynamicCast(pmesh);
+ const volScalarField& p = *pPtr;
+
+ // face centres (on patches)
+ const surfaceVectorField::Boundary& patchCf =
+ mesh.Cf().boundaryField();
+
+ // face areas (on patches)
+ const surfaceVectorField::Boundary& patchSf =
+ mesh.Sf().boundaryField();
+
+ // Pressure (on patches)
+ const volScalarField::Boundary& patchPress =
+ p.boundaryField();
+
+ // rhoRef if the pressure field is dynamic, i.e. p/rho otherwise 1
+ rhoRef = (p.dimensions() == dimPressure ? 1.0 : rhoRef);
+
+ // Scale pRef by density for incompressible simulations
+ pRef /= rhoRef;
+
+ forAll(faceZones_, zonei)
+ {
+ const labelList& zn = faceZones_[zonei];
+
+ forAll(zn, localFaceI)
+ {
+ const label meshFaceI = zn[localFaceI];
+
+ // locate which patch this face is on,
+ // and its offset within the patch
+ const label patchI = patches.whichPatch(meshFaceI);
+ if (patchI == -1)
+ {
+ // Should not happen - could also warn
+ continue;
+ }
+
+ const label patchFaceI = meshFaceI - patches[patchI].start();
+
+ if (!forceOnPatches.found(patchI))
+ {
+ // Patch faces are +ve outwards,
+ // so the forces (exerted by fluid on solid)
+ // already have the correct sign
+ forceOnPatches.set
+ (
+ patchI,
+ (
+ rhoRef
+ * patchSf[patchI] * (patchPress[patchI] - pRef)
+ ).ptr()
+ );
+ }
+
+ const vectorField& forceOnPatch = *forceOnPatches[patchI];
+
+ // Force from the patch-face into sum
+ forces[zonei] += forceOnPatch[patchFaceI];
+
+ // effective torque arm:
+ // - translated into the lumped-points coordinate system
+ // prior to determining the distance
+ const vector lever =
+ (
+ patchCf[patchI][patchFaceI] - centre_
+ - lumpedCentres[zonei]
+ );
+
+ // Moment from the patch-face into sum
+ moments[zonei] += lever ^ forceOnPatch[patchFaceI];
+ }
+ }
+ }
+ else
+ {
+ Info<<"no pressure" << endl;
+ }
+
+ Pstream::listCombineGather(forces, plusEqOp());
+ Pstream::listCombineScatter(forces);
+
+ Pstream::listCombineGather(moments, plusEqOp());
+ Pstream::listCombineScatter(moments);
+
+ return true;
+}
+
+
+// \cond fileScope
+template
+static inline T interpolatedValue
+(
+ const Foam::UList input,
+ const Foam::labelUList& indices,
+ const Foam::scalarField& weights
+)
+{
+ T result = weights[0] * input[indices[0]];
+ for (Foam::label i=1; i < indices.size(); ++i)
+ {
+ result += weights[i] * input[indices[i]];
+ }
+
+ return result;
+}
+// \endcond
+
+
+Foam::tmp
+Foam::lumpedPointMovement::displacePoints
+(
+ const pointField& points0,
+ const labelList& pointLabels
+) const
+{
+ return displacePoints(state(), points0, pointLabels);
+}
+
+
+Foam::tmp
+Foam::lumpedPointMovement::displacePoints
+(
+ const lumpedPointState& state,
+ const pointField& points0,
+ const labelList& pointLabels
+) const
+{
+ // storage for the interpolation indices/weights
+ labelList indices;
+ scalarField weights;
+ const interpolationWeights& interp = interpolator();
+
+ // The mass centres
+ const pointField& lumpedCentres = state.points();
+
+ // local-to-global transformation tensor
+ const tensorField& localToGlobal = state.rotations();
+
+ // could also verify the sizes (state vs original)
+
+ tmp tdisp(new pointField(pointLabels.size()));
+ pointField& disp = tdisp.ref();
+
+ forAll(pointLabels, ptI)
+ {
+ const point& p0 = points0[pointLabels[ptI]];
+
+ // Interpolation factor based on the axis component
+ // of the original point (location)
+ scalar where = p0 & axis_;
+
+ interp.valueWeights(where, indices, weights);
+
+ vector origin = interpolatedValue(lumpedCentres, indices, weights);
+ tensor rotTensor = interpolatedValue(localToGlobal, indices, weights);
+
+ if (indices.size() == 1)
+ {
+ // out-of-bounds, use corresponding end-point
+ where = locations_[indices[0]];
+ }
+
+ // local point:
+ // - translated into the lumped-points coordinate system
+ // - relative to the interpolation (where) location
+ const vector local = (p0 - (where * axis_)) - centre_;
+
+ // local-to-global rotation and translation
+ // - translate back out of the lumped-points coordinate system
+ //
+ // -> subtract p0 to return displacements
+ disp[ptI] = (rotTensor & local) + origin + centre_ - p0;
+ }
+
+ return tdisp;
+}
+
+
+void Foam::lumpedPointMovement::writeDict(Ostream& os) const
+{
+ os.writeEntry("axis", axis_);
+ os.writeEntry("locations", locations_);
+ os.writeEntry("division", division_);
+ // os.writeEntry("interpolationScheme", interpolationScheme_);
+}
+
+
+bool Foam::lumpedPointMovement::readState()
+{
+ lumpedPointState prev = state_;
+
+ bool status = state_.readData
+ (
+ inputFormat_,
+ coupler().resolveFile(inputName_)
+ );
+
+ state_.relax(relax_, prev);
+
+ return status;
+}
+
+
+bool Foam::lumpedPointMovement::writeData
+(
+ const UList& forces
+) const
+{
+ if (!Pstream::master())
+ {
+ return false;
+ }
+
+ const fileName output(coupler().resolveFile(outputName_));
+ OFstream os(output); // ASCII
+
+ if (outputFormat_ == outputFormatType::PLAIN)
+ {
+ os <<"# output from OpenFOAM" << nl
+ <<"# N, points, forces" << nl
+ << this->size() << nl;
+
+ const char* zeroVector = "0 0 0";
+
+ forAll(locations_, i)
+ {
+ const vector pos = locations_[i] * axis_;
+
+ os << pos.x() << ' '
+ << pos.y() << ' '
+ << pos.z() << ' ';
+
+ if (i < forces.size())
+ {
+ os << forces[i].x() << ' '
+ << forces[i].y() << ' '
+ << forces[i].z();
+ }
+ else
+ {
+ os << zeroVector;
+ }
+
+ os << nl;
+ }
+ }
+ else
+ {
+ // Make it easier for external programs to parse
+ // - exclude the usual OpenFOAM 'FoamFile' header
+ // - ensure lists have consistent format
+
+ os <<"// output from OpenFOAM" << nl << nl;
+
+ writeList(os, "points", (locations_*axis_)());
+ writeList(os, "forces", forces);
+ }
+
+ return true;
+}
+
+
+bool Foam::lumpedPointMovement::writeData
+(
+ const UList& forces,
+ const UList& moments
+) const
+{
+ if (!Pstream::master())
+ {
+ return false;
+ }
+
+ const fileName output(coupler().resolveFile(outputName_));
+ OFstream os(output); // ASCII
+
+ if (outputFormat_ == outputFormatType::PLAIN)
+ {
+ os <<"# output from OpenFOAM" << nl
+ <<"# N, points, forces, moments" << nl
+ << this->size() << nl;
+
+ const char* zeroVector = "0 0 0";
+
+ forAll(locations_, i)
+ {
+ const vector pos = locations_[i] * axis_;
+
+ os << pos.x() << ' '
+ << pos.y() << ' '
+ << pos.z() << ' ';
+
+ if (i < forces.size())
+ {
+ os << forces[i].x() << ' '
+ << forces[i].y() << ' '
+ << forces[i].z() << ' ';
+ }
+ else
+ {
+ os << zeroVector << ' ';
+ }
+
+ if (i < moments.size())
+ {
+ os << moments[i].x() << ' '
+ << moments[i].y() << ' '
+ << moments[i].z();
+ }
+ else
+ {
+ os << zeroVector;
+ }
+ os << nl;
+ }
+ }
+ else
+ {
+ // Make it easier for external programs to parse
+ // - exclude the usual OpenFOAM 'FoamFile' header
+ // - ensure lists have consistent format
+
+ os <<"// output from OpenFOAM" << nl << nl;
+
+ writeList(os, "points", (locations_*axis_)());
+ writeList(os, "forces", forces);
+ writeList(os, "moments", moments);
+ }
+
+ return true;
+}
+
+
+const Foam::interpolationWeights&
+Foam::lumpedPointMovement::interpolator() const
+{
+ if (!interpolatorPtr_.valid())
+ {
+ interpolatorPtr_ = interpolationWeights::New
+ (
+ interpolationScheme_,
+ locations_
+ );
+ }
+
+ return interpolatorPtr_();
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointMovement.H b/src/lumpedPointMotion/lumpedPointMovement.H
new file mode 100644
index 0000000000..3a6733b127
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointMovement.H
@@ -0,0 +1,399 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see formatNames;
+
+
+private:
+
+ // Private data
+
+ //- Reference axis for the locations
+ vector axis_;
+
+ //- The locations of the lumped points within the reference axis
+ // The interpolator needs scalarField, not scalarList.
+ scalarField locations_;
+
+ //- The division when calculating pressure forces (0-1)
+ scalar division_;
+
+ //- The relaxation factor when moving points (default: 1)
+ scalar relax_;
+
+ //- The interpolation type (linear|spline)
+ word interpolationScheme_;
+
+ //- Optional owner information
+ label ownerId_;
+
+ //- The bounding box (if set)
+ boundBox boundBox_;
+
+ //- The offset of patch points to compared to the locations
+ point centre_;
+
+ //- Calculate centre based on the bounding box
+ bool autoCentre_;
+
+ //- Dictionary of controls for force calculation
+ dictionary forcesDict_;
+
+ //- Communication control
+ externalCoupler coupler_;
+
+ //- File io
+ word inputName_;
+ word outputName_;
+ lumpedPointState::inputFormatType inputFormat_;
+ outputFormatType outputFormat_;
+
+
+ // Demand-driven private data
+
+ //- The initial state (positions, rotations)
+ lumpedPointState state0_;
+
+ //- The current state (positions, rotations)
+ // Eg, as response from external application
+ lumpedPointState state_;
+
+ //- Threshold locations for pressure forces
+ mutable scalarField* thresholdPtr_;
+
+ //- User-specified interpolator
+ mutable autoPtr interpolatorPtr_;
+
+ //- Pressure zones (only used from the master patch)
+ mutable List faceZones_;
+
+
+ // Private Member Functions
+
+ //- Calculate threshold locations
+ void calcThresholds() const;
+
+ //- Classify the position to be located in one of the threshold zones
+ label threshold(scalar pos) const;
+
+
+ //- Disallow default bitwise copy constructor
+ lumpedPointMovement(const lumpedPointMovement&) = delete;
+
+ //- Disallow default bitwise assignment
+ void operator=(const lumpedPointMovement&) = delete;
+
+public:
+
+ // Static data members
+
+ //- The standard dictionary name
+ static const word dictionaryName;
+
+ //- Helper: return IOobject with optionally absolute path provided
+ static IOobject selectIO(const IOobject&, const fileName&);
+
+
+ // Constructors
+
+ //- Construct null
+ lumpedPointMovement();
+
+ //- Construct from dictionary, optionally with some owner information
+ lumpedPointMovement(const dictionary& dict, label ownerId=-1);
+
+
+ //- Destructor
+ virtual ~lumpedPointMovement();
+
+
+ // Member Functions
+
+ //- Update settings from dictionary
+ void readDict(const dictionary& dict);
+
+
+ //- If no number of lumped points (locations) were specified
+ inline bool empty() const;
+
+ //- The number of lumped points (number of locations)
+ inline label size() const;
+
+ //- The normalized reference axis
+ inline const vector& axis() const;
+
+ //- Read access to the locations
+ inline const scalarField& locations() const;
+
+ //- The division (0-1) when calculating pressure forces
+ inline scalar division() const;
+
+ //- An owner Id, if needed for bookkeeping purposes
+ inline label ownerId() const;
+
+ //- Change the owner id, if needed for bookkeeping purposes
+ inline void ownerId(label id);
+
+
+ //- Threshold locations for pressure forces
+ inline const scalarField& thresholds() const;
+
+ //- Classify the position to be located in one of the threshold zones
+ inline label threshold(const point& position) const;
+
+
+ //- Communication control
+ inline const externalCoupler& coupler() const;
+
+ //- Communication control
+ inline externalCoupler& coupler();
+
+ //- The initial state (positions/rotations)
+ inline const lumpedPointState& state0() const;
+
+ //- The current state (positions/rotations)
+ inline const lumpedPointState& state() const;
+
+ //- The relaxation factor when changing states
+ inline scalar relax() const;
+
+ //- The relaxation factor when changing states
+ inline scalar& relax();
+
+ //- The input (state) file name
+ inline const word& inputName() const;
+
+ //- The output (forces) file name
+ inline const word& outputName() const;
+
+ //- The input (state) file format
+ inline lumpedPointState::inputFormatType inputFormat() const;
+
+ //- The output (forces) file format
+ inline lumpedPointMovement::outputFormatType outputFormat() const;
+
+
+ //- Define the bounding-box required to enclose the specified patches.
+ // Calculates the centre as required.
+ //
+ // \param mesh The volume mesh reference
+ // \param patchIds The patch ids to be included in the bounding box.
+ // \param points0 The initial mesh points, prior to movement
+ void setBoundBox
+ (
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+ );
+
+
+ //- Define the pressure-zones mapping for faces in the specified
+ // patches.
+ // The face centres are compared to the threshold positions,
+ // which are determined by locations along the defined axis.
+ //
+ // \param mesh The volume mesh reference
+ // \param patchIds The patch ids to be included in the mapping
+ // \param points0 The initial mesh points, prior to movement
+ void setMapping
+ (
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+ );
+
+
+ //- True if the pressure-zones mapping has already been performed
+ inline bool hasMapping() const;
+
+ //- Return the pressure-zones mapping with the associated
+ // patch face ids.
+ inline const List& zones() const;
+
+
+ //- The forces and moments acting on each pressure-zone.
+ // The zones must be previously defined via setMapping.
+ bool forcesAndMoments
+ (
+ const polyMesh& pmesh,
+ List& forces,
+ List& moments
+ ) const;
+
+
+ //- Displace points according to the current state
+ tmp displacePoints
+ (
+ const pointField& points0,
+ const labelList& pointLabels
+ ) const;
+
+
+ //- Displace points according to specified state
+ tmp displacePoints
+ (
+ const lumpedPointState& state,
+ const pointField& points0,
+ const labelList& pointLabels
+ ) const;
+
+
+ //- Interpolation weights
+ const interpolationWeights& interpolator() const;
+
+ //- Write axis, locations, division as a dictionary
+ void writeDict(Ostream& os) const;
+
+
+ //- Write points, forces
+ bool writeData
+ (
+ const UList& forces
+ ) const;
+
+ //- Write points, forces, moments
+ bool writeData
+ (
+ const UList& forces,
+ const UList& moments
+ ) const;
+
+
+ //- Read state from file, applying relaxation as requested
+ bool readState();
+
+ //- Write state as VTK PolyData format.
+ void writeStateVTP(const fileName& file) const;
+
+ //- Write forces on points as VTK PolyData format.
+ void writeForcesAndMomentsVTP
+ (
+ const fileName& file,
+ const UList& forces,
+ const UList& moments
+ ) const;
+
+
+ //- Write pressure-zones geometry, write as VTK PolyData format.
+ void writeZonesVTP
+ (
+ const fileName& file,
+ const polyMesh& mesh,
+ const pointField& points0
+ ) const;
+
+
+ //- Write displaced geometry according to the current state,
+ // write as VTK PolyData format.
+ void writeVTP
+ (
+ const fileName& file,
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+ ) const;
+
+ //- Write displaced geometry according to the specified state,
+ // write as VTK PolyData format.
+ void writeVTP
+ (
+ const fileName& file,
+ const lumpedPointState& state,
+ const polyMesh& mesh,
+ const labelUList& patchLst,
+ const pointField& points0
+ ) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "lumpedPointMovementI.H"
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointMovementI.H b/src/lumpedPointMotion/lumpedPointMovementI.H
new file mode 100644
index 0000000000..36f189d1a2
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointMovementI.H
@@ -0,0 +1,162 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see &
+Foam::lumpedPointMovement::zones() const
+{
+ return faceZones_;
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointMovementWriter.C b/src/lumpedPointMotion/lumpedPointMovementWriter.C
new file mode 100644
index 0000000000..c8fcc888b4
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointMovementWriter.C
@@ -0,0 +1,447 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see & forces,
+ const UList& moments
+) const
+{
+ OFstream fos(file);
+ std::ostream& os = fos.stdStream();
+
+ autoPtr format = vtk::newFormatter
+ (
+ os,
+ vtk::formatType::INLINE_ASCII
+ );
+
+ format().xmlHeader()
+ .beginVTKFile(vtk::fileTag::POLY_DATA, "0.1");
+
+ //
+ // The 'spine' of lumped mass points
+ //
+ const label nPoints = state().points().size();
+
+ {
+ format()
+ .openTag(vtk::fileTag::PIECE)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POINTS, nPoints)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_VERTS, nPoints)
+ .closeTag();
+
+ // 'points'
+ {
+ const uint64_t payLoad = (nPoints*3* sizeof(float));
+
+ format()
+ .tag(vtk::fileTag::POINTS)
+ .openDataArray(vtk::dataArrayAttr::POINTS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ vtk::writeList(format(), state().points());
+ format().flush();
+
+ format()
+ .endDataArray()
+ .endTag(vtk::fileTag::POINTS);
+ }
+
+ //
+ format().tag(vtk::fileTag::VERTS);
+
+ //
+ // 'connectivity'
+ //
+ {
+ const uint64_t payLoad = (nPoints*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::CONNECTIVITY)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ for (label i=0; i(vtk::dataArrayAttr::OFFSETS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ for (label i=0; i
+ }
+
+ format().tag(vtk::fileTag::POINT_DATA);
+
+ // forces
+ if (forces.size() == nPoints)
+ {
+ const uint64_t payLoad = (nPoints * 3 * sizeof(float));
+
+ format().openDataArray("forces")
+ .closeTag();
+
+ format().writeSize(payLoad);
+ vtk::writeList(format(), forces);
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ // moments
+ if (moments.size() == nPoints)
+ {
+ const uint64_t payLoad = (nPoints * 3 * sizeof(float));
+
+ format().openDataArray("moments")
+ .closeTag();
+
+ format().writeSize(payLoad);
+ vtk::writeList(format(), moments);
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::POINT_DATA);
+
+ format().endTag(vtk::fileTag::PIECE);
+ format().endTag(vtk::fileTag::POLY_DATA);
+ format().endVTKFile();
+}
+
+
+void Foam::lumpedPointMovement::writeZonesVTP
+(
+ const fileName& file,
+ const polyMesh& mesh,
+ const pointField& points0
+) const
+{
+ OFstream fos(file);
+ std::ostream& os = fos.stdStream();
+
+ autoPtr format = vtk::newFormatter
+ (
+ os,
+ vtk::formatType::INLINE_ASCII
+ );
+
+ format().xmlHeader()
+ .beginVTKFile(vtk::fileTag::POLY_DATA, "0.1");
+
+ forAll(faceZones_, zoneI)
+ {
+ uindirectPrimitivePatch pp
+ (
+ UIndirectList(mesh.faces(), faceZones_[zoneI]),
+ points0
+ );
+
+ format()
+ .openTag(vtk::fileTag::PIECE)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POINTS, pp.nPoints())
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POLYS, pp.size())
+ .closeTag();
+
+ // 'points'
+ {
+ const uint64_t payLoad = (pp.nPoints()*3* sizeof(float));
+
+ format()
+ .tag(vtk::fileTag::POINTS)
+ .openDataArray(vtk::dataArrayAttr::POINTS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ vtk::writeList(format(), pp.localPoints());
+ format().flush();
+
+ format()
+ .endDataArray()
+ .endTag(vtk::fileTag::POINTS);
+ }
+
+ //
+ format().tag(vtk::fileTag::POLYS);
+
+ //
+ // 'connectivity'
+ //
+ {
+ uint64_t payLoad = 0;
+ forAll(pp, facei)
+ {
+ const face& f = pp[facei];
+ payLoad += f.size();
+ }
+
+ format().openDataArray(vtk::dataArrayAttr::CONNECTIVITY)
+ .closeTag();
+
+ format().writeSize(payLoad * sizeof(label));
+
+ for (const face& f : pp.localFaces())
+ {
+ vtk::writeList(format(), f);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ //
+ // 'offsets' (connectivity offsets)
+ //
+ {
+ const uint64_t payLoad = (pp.size() * sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::OFFSETS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+
+ label off = 0;
+ forAll(pp, facei)
+ {
+ const face& f = pp[facei];
+ off += f.size();
+
+ format().write(off);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::POLYS);
+
+
+ format().tag(vtk::fileTag::CELL_DATA);
+
+ // zone Id
+ {
+ const uint64_t payLoad = (pp.size() * sizeof(label));
+
+ format().openDataArray("zoneId")
+ .closeTag();
+
+ format().writeSize(payLoad);
+
+ forAll(pp, facei)
+ {
+ format().write(zoneI);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::CELL_DATA);
+
+ format().endTag(vtk::fileTag::PIECE);
+ }
+
+ format().endTag(vtk::fileTag::POLY_DATA);
+ format().endVTKFile();
+}
+
+
+void Foam::lumpedPointMovement::writeVTP
+(
+ const fileName& file,
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+) const
+{
+ writeVTP(file, state(), mesh, patchIds, points0);
+}
+
+
+void Foam::lumpedPointMovement::writeVTP
+(
+ const fileName& file,
+ const lumpedPointState& state,
+ const polyMesh& mesh,
+ const labelUList& patchIds,
+ const pointField& points0
+) const
+{
+ const polyBoundaryMesh& boundaryMesh = mesh.boundaryMesh();
+
+ OFstream fos(file);
+ std::ostream& os = fos.stdStream();
+
+ autoPtr format = vtk::newFormatter
+ (
+ os,
+ vtk::formatType::INLINE_ASCII
+ );
+
+ format().xmlHeader()
+ .beginVTKFile(vtk::fileTag::POLY_DATA, "0.1");
+
+ for (const label patchId : patchIds)
+ {
+ const polyPatch& pp = boundaryMesh[patchId];
+
+ format()
+ .openTag(vtk::fileTag::PIECE)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POINTS, pp.nPoints())
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POLYS, pp.size())
+ .closeTag();
+
+ // 'points'
+ {
+ const uint64_t payLoad = (pp.nPoints()*3* sizeof(float));
+
+ format()
+ .tag(vtk::fileTag::POINTS)
+ .openDataArray(vtk::dataArrayAttr::POINTS)
+ .closeTag();
+
+ // Could be more efficient, but not often needed
+ tmp tpts = displacePoints
+ (
+ state,
+ points0,
+ pp.meshPoints()
+ ) + pointField(points0, pp.meshPoints());
+
+ const pointField& pts = tpts();
+
+ format().writeSize(payLoad);
+ vtk::writeList(format(), pts);
+ format().flush();
+
+ format()
+ .endDataArray()
+ .endTag(vtk::fileTag::POINTS);
+ }
+
+ //
+ format().tag(vtk::fileTag::POLYS);
+
+
+ //
+ // 'connectivity'
+ //
+ {
+ uint64_t payLoad = 0;
+ for (const face& f : pp)
+ {
+ payLoad += f.size();
+ }
+
+ format().openDataArray(vtk::dataArrayAttr::CONNECTIVITY)
+ .closeTag();
+
+ format().writeSize(payLoad * sizeof(label));
+
+ for (const face& f : pp.localFaces())
+ {
+ vtk::writeList(format(), f);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ //
+ // 'offsets' (connectivity offsets)
+ //
+ {
+ const uint64_t payLoad = (pp.size() * sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::OFFSETS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+
+ label off = 0;
+ forAll(pp, facei)
+ {
+ const face& f = pp[facei];
+ off += f.size();
+
+ format().write(off);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::POLYS);
+
+ format().endTag(vtk::fileTag::PIECE);
+ }
+
+ format().endTag(vtk::fileTag::POLY_DATA);
+ format().endVTKFile();
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointState.C b/src/lumpedPointMotion/lumpedPointState.C
new file mode 100644
index 0000000000..511c756697
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointState.C
@@ -0,0 +1,362 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see
+ Foam::lumpedPointState::formatNames
+ {
+ { inputFormatType::PLAIN, "plain" },
+ { inputFormatType::DICTIONARY, "dictionary" }
+ };
+
+
+// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
+
+//! \cond fileScope
+static Foam::string getLineNoComment(Foam::ISstream& is)
+{
+ Foam::string line;
+ do
+ {
+ is.getLine(line);
+ }
+ while ((line.empty() || line[0] == '#') && is.good());
+
+ return line;
+}
+//! \endcond
+
+
+// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
+
+void Foam::lumpedPointState::calcRotations() const
+{
+ rotationPtr_ = new tensorField(angles_.size());
+ forAll(angles_, itemi)
+ {
+ rotationPtr_->operator[](itemi) = EulerCoordinateRotation
+ (
+ angles_[itemi],
+ degrees_ // true=degrees, false=radians
+ ).R();
+ }
+}
+
+
+void Foam::lumpedPointState::readDict(const dictionary& dict)
+{
+ dict.lookup("points") >> points_;
+ dict.lookup("angles") >> angles_;
+ degrees_ = dict.lookupOrDefault("degrees", false);
+ deleteDemandDrivenData(rotationPtr_);
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
+
+Foam::lumpedPointState::lumpedPointState()
+:
+ points_(0),
+ angles_(0),
+ degrees_(false),
+ rotationPtr_(nullptr)
+{}
+
+
+Foam::lumpedPointState::lumpedPointState(const lumpedPointState& rhs)
+:
+ points_(rhs.points_),
+ angles_(rhs.angles_),
+ degrees_(rhs.degrees_),
+ rotationPtr_(nullptr)
+{}
+
+
+Foam::lumpedPointState::lumpedPointState
+(
+ const pointField& pts
+)
+:
+ points_(pts),
+ angles_(points_.size(), Zero),
+ degrees_(false),
+ rotationPtr_(nullptr)
+{}
+
+
+Foam::lumpedPointState::lumpedPointState
+(
+ tmp& pts
+)
+:
+ points_(pts),
+ angles_(points_.size(), Zero),
+ degrees_(false),
+ rotationPtr_(nullptr)
+{}
+
+
+Foam::lumpedPointState::lumpedPointState
+(
+ const dictionary& dict
+)
+:
+ points_(0),
+ angles_(0),
+ degrees_(false),
+ rotationPtr_(nullptr)
+{
+ readDict(dict);
+}
+
+
+// * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * * //
+
+Foam::lumpedPointState::~lumpedPointState()
+{
+ deleteDemandDrivenData(rotationPtr_);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
+
+void Foam::lumpedPointState::operator=(const lumpedPointState& rhs)
+{
+ points_ = rhs.points_;
+ angles_ = rhs.angles_;
+ degrees_ = rhs.degrees_;
+
+ deleteDemandDrivenData(rotationPtr_);
+}
+
+
+void Foam::lumpedPointState::relax
+(
+ const scalar alpha,
+ const lumpedPointState& prev
+)
+{
+ points_ = prev.points_ + alpha*(points_ - prev.points_);
+
+ scalar convert = 1.0;
+ if (degrees_ != prev.degrees_)
+ {
+ if (prev.degrees_)
+ {
+ // Was degrees, now radians
+ convert = degToRad(1);
+ }
+ else
+ {
+ // Was radians, now degrees
+ convert = radToDeg(1);
+ }
+ }
+
+ angles_ = convert*prev.angles_ + alpha*(angles_ - convert*prev.angles_);
+
+ deleteDemandDrivenData(rotationPtr_);
+}
+
+
+bool Foam::lumpedPointState::readPlain(Istream& is)
+{
+ // Assume generic input stream so we can do line-based input
+ ISstream& iss = dynamic_cast(is);
+
+ string line = getLineNoComment(iss);
+
+ label count;
+ {
+ IStringStream isstr(line);
+ isstr >> count;
+ }
+
+ points_.setSize(count);
+ angles_.setSize(count);
+
+ count = 0;
+ forAll(points_, i)
+ {
+ iss.getLine(line);
+ IStringStream isstr(line);
+
+ isstr
+ >> points_[count].x() >> points_[count].y() >> points_[count].z()
+ >> angles_[count].x() >> angles_[count].y() >> angles_[count].z();
+
+ ++count;
+ }
+
+ points_.setSize(count);
+ angles_.setSize(count);
+
+ degrees_ = false;
+ deleteDemandDrivenData(rotationPtr_);
+
+ return count;
+}
+
+
+bool Foam::lumpedPointState::readData(Istream& is)
+{
+ dictionary dict(is);
+ readDict(dict);
+
+ return points_.size();
+}
+
+
+bool Foam::lumpedPointState::writeData(Ostream& os) const
+{
+ writeDict(os);
+ return true;
+}
+
+
+void Foam::lumpedPointState::writeDict(Ostream& os) const
+{
+ os.writeEntry("points", points_);
+ os.writeEntry("angles", angles_);
+ if (degrees_)
+ {
+ os.writeKeyword("degrees") << "true;" << nl;
+ }
+}
+
+
+void Foam::lumpedPointState::writePlain(Ostream& os) const
+{
+ os <<"# input for OpenFOAM\n"
+ <<"# N, points, angles\n"
+ << points_.size() << "\n";
+
+ forAll(points_, i)
+ {
+ const vector& pt = points_[i];
+
+ os << pt.x() << ' '
+ << pt.y() << ' '
+ << pt.z() << ' ';
+
+ if (i < angles_.size())
+ {
+ os << angles_[i].x() << ' '
+ << angles_[i].y() << ' '
+ << angles_[i].z() << '\n';
+ }
+ else
+ {
+ os << "0 0 0\n";
+ }
+ }
+}
+
+
+bool Foam::lumpedPointState::readData
+(
+ const inputFormatType& fmt,
+ const fileName& file
+)
+{
+ bool ok = false;
+ if (Pstream::master())
+ {
+ IFstream is(file);
+
+ if (fmt == inputFormatType::PLAIN)
+ {
+ ok = this->readPlain(is);
+ }
+ else
+ {
+ ok = this->readData(is);
+ }
+ }
+
+ if (Pstream::parRun())
+ {
+ // Scatter master data using communication scheme
+
+ const List& comms =
+ (
+ (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+ ? Pstream::linearCommunication()
+ : Pstream::treeCommunication()
+ );
+
+ // Get my communication order
+ const Pstream::commsStruct& myComm = comms[Pstream::myProcNo()];
+
+ // Receive from up
+ if (myComm.above() != -1)
+ {
+ IPstream fromAbove
+ (
+ UPstream::commsTypes::scheduled,
+ myComm.above(),
+ 0,
+ Pstream::msgType(),
+ Pstream::worldComm
+ );
+
+ fromAbove >> points_ >> angles_ >> degrees_;
+ }
+
+ // Send to downstairs neighbours
+ forAllReverse(myComm.below(), belowI)
+ {
+ OPstream toBelow
+ (
+ UPstream::commsTypes::scheduled,
+ myComm.below()[belowI],
+ 0,
+ Pstream::msgType(),
+ Pstream::worldComm
+ );
+
+ toBelow << points_ << angles_ << degrees_;
+ }
+
+ deleteDemandDrivenData(rotationPtr_);
+
+ // MPI barrier
+ Pstream::scatter(ok);
+ }
+
+ return ok;
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointState.H b/src/lumpedPointMotion/lumpedPointState.H
new file mode 100644
index 0000000000..0d32fa8927
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointState.H
@@ -0,0 +1,194 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see formatNames;
+
+
+private:
+
+ // Private data
+
+ //- Positions of lumped points
+ pointField points_;
+
+ //- Orientation of lumped points (as Euler angles)
+ vectorField angles_;
+
+ //- Euler angles in degrees instead radians
+ bool degrees_;
+
+ //- Tensor rotation of lumped points
+ mutable tensorField* rotationPtr_;
+
+
+ // Private Member Functions
+
+ void calcRotations() const;
+
+ void readDict(const dictionary& dict);
+
+public:
+
+ // Constructors
+
+ //- Construct null
+ lumpedPointState();
+
+ //- Copy constructor
+ lumpedPointState(const lumpedPointState& rhs);
+
+ //- Construct from points with zero-rotation
+ lumpedPointState(const pointField& pts);
+
+ //- Construct from points with zero-rotation
+ lumpedPointState(tmp& pts);
+
+ //- Construct from dictionary
+ lumpedPointState(const dictionary& dict);
+
+
+ //- Destructor
+ virtual ~lumpedPointState();
+
+
+ // Member Functions
+
+ //- Has positions and consistent number of rotations?
+ inline bool valid() const;
+
+ //- If no points were specified
+ inline bool empty() const;
+
+ //- The number of points
+ inline label size() const;
+
+ //- The points corresponding to mass centres
+ inline const pointField& points() const;
+
+ //- The orientation of the points (mass centres)
+ inline const vectorField& angles() const;
+
+ //- The local-to-global transformation for each point
+ inline const tensorField& rotations() const;
+
+ //- Relax the state
+ // alpha = 1 : no relaxation
+ // alpha < 1 : relaxation
+ // alpha = 0 : do nothing
+ void relax(const scalar alpha, const lumpedPointState& prev);
+
+
+ //- Read input as dictionary content
+ bool readData(Istream& is);
+
+ //- Output as dictionary content
+ bool writeData(Ostream& os) const;
+
+ //- Output as dictionary content
+ void writeDict(Ostream& os) const;
+
+ //- Read input as plain content
+ bool readPlain(Istream& is);
+
+ //- Output as plain content
+ void writePlain(Ostream& os) const;
+
+ //- Read input from file (master process only) using specified format
+ bool readData(const inputFormatType& fmt, const fileName& file);
+
+ //- Output as VTK file for debugging/visualization
+ // The points are joined as lines, the rotation is visualized
+ // by planes, write as VTK PolyData format.
+ void writeVTP(const fileName& outputFile, const vector& axis) const;
+
+
+ // Member Operators
+
+ //- Assignment operator
+ void operator=(const lumpedPointState& rhs);
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "lumpedPointStateI.H"
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointStateI.H b/src/lumpedPointMotion/lumpedPointStateI.H
new file mode 100644
index 0000000000..af738a670a
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointStateI.H
@@ -0,0 +1,67 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see standardCorners
+{
+ {-0.1, -0.1, 0},
+ {+0.1, -0.1, 0},
+ {+0.1, +0.1, 0},
+ {-0.1, +0.1, 0}
+};
+
+
+// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
+
+void Foam::lumpedPointState::writeVTP
+(
+ const fileName& outputFile,
+ const vector& axis
+) const
+{
+ // local-to-global transformation tensors
+ const tensorField& localToGlobal = rotations();
+
+ OFstream fos(outputFile);
+ std::ostream& os = fos.stdStream();
+
+ autoPtr format = vtk::newFormatter
+ (
+ os,
+ vtk::formatType::INLINE_ASCII
+ );
+
+ format().xmlHeader()
+ .beginVTKFile(vtk::fileTag::POLY_DATA, "0.1");
+
+ //
+ // The 'spine' of lumped mass points
+ //
+ {
+ format()
+ .openTag(vtk::fileTag::PIECE)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POINTS, points_.size())
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_VERTS, points_.size())
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_LINES, 1)
+ .closeTag();
+
+ // 'points'
+ {
+ const uint64_t payLoad = (points_.size()*3* sizeof(float));
+
+ format()
+ .tag(vtk::fileTag::POINTS)
+ .openDataArray(vtk::dataArrayAttr::POINTS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ vtk::writeList(format(), points_);
+ format().flush();
+
+ format()
+ .endDataArray()
+ .endTag(vtk::fileTag::POINTS);
+ }
+
+ //
+ format().tag(vtk::fileTag::VERTS);
+
+ //
+ // 'connectivity'
+ //
+ {
+ const uint64_t payLoad = (points_.size()*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::CONNECTIVITY)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ forAll(points_, i)
+ {
+ format().write(i);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ //
+ // 'offsets' (connectivity offsets)
+ // = linear mapping onto points (with 1 offset)
+ //
+ {
+ const uint64_t payLoad = (points_.size()*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::OFFSETS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ forAll(points_, i)
+ {
+ format().write(i+1);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::VERTS);
+ //
+
+
+ //
+ format().tag(vtk::fileTag::LINES);
+
+ //
+ // 'connectivity'
+ //
+ {
+ const uint64_t payLoad = (points_.size()*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::CONNECTIVITY)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ forAll(points_, i)
+ {
+ format().write(i);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ //
+ // 'offsets' (connectivity offsets)
+ // = single line
+ //
+ {
+ const uint64_t payLoad = (1*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::OFFSETS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ format().write(points_.size());
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::LINES);
+ format().endTag(vtk::fileTag::PIECE);
+ }
+
+ // Standard corners in local axis
+ const axesRotation cornerTransform(axis);
+
+ FixedList corners;
+ forAll(standardCorners, corni)
+ {
+ corners[corni] = cornerTransform.transform(standardCorners[corni]);
+ }
+
+
+ //
+ // Planes to visualize location/rotation
+ //
+ {
+ const label nPoints = 4*points_.size(); // 4 points per quad
+ const label nPolys = points_.size();
+
+ format()
+ .openTag(vtk::fileTag::PIECE)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POINTS, nPoints)
+ .xmlAttr(vtk::fileAttr::NUMBER_OF_POLYS, nPolys)
+ .closeTag();
+
+ // 'points'
+ {
+ const uint64_t payLoad = (nPoints*3*sizeof(float));
+
+ format()
+ .tag(vtk::fileTag::POINTS)
+ .openDataArray(vtk::dataArrayAttr::POINTS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+
+ forAll(points_, posI)
+ {
+ const point& origin = points_[posI];
+ const tensor& rotTensor =
+ (
+ posI < localToGlobal.size()
+ ? localToGlobal[posI]
+ : pTraits::I
+ );
+
+
+ for (const point& cornerPt : corners)
+ {
+ // Local-to-global rotation and translation
+ const point pt = (rotTensor & cornerPt) + origin;
+
+ vtk::write(format(), pt);
+ }
+ }
+
+ format().flush();
+
+ format()
+ .endDataArray()
+ .endTag(vtk::fileTag::POINTS);
+ }
+
+ //
+ format().tag(vtk::fileTag::POLYS);
+
+ //
+ // 'connectivity' - 4 points (ie, quad)
+ //
+ {
+ const uint64_t payLoad = (4*nPolys*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::CONNECTIVITY)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ for (label i=0; i < 4*nPolys; ++i)
+ {
+ format().write(i);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ //
+ // 'offsets' (connectivity offsets)
+ // = single quad
+ //
+ {
+ const uint64_t payLoad = (nPolys*sizeof(label));
+
+ format().openDataArray(vtk::dataArrayAttr::OFFSETS)
+ .closeTag();
+
+ format().writeSize(payLoad);
+ for (label i=0; i < nPolys; ++i)
+ {
+ const label off = 4 * (i+1);
+ format().write(off);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::POLYS);
+
+#if 0
+ format().tag(vtk::fileTag::CELL_DATA);
+
+ // zone Id
+ {
+ const uint64_t payLoad = (points_.size()*sizeof(label));
+
+ format().openDataArray("zoneId")
+ .closeTag();
+
+ format().writeSize(payLoad);
+ for (label i=0; i < nPolys; ++i)
+ {
+ format().write(i);
+ }
+ format().flush();
+
+ format().endDataArray();
+ }
+
+ format().endTag(vtk::fileTag::CELL_DATA);
+#endif
+
+ format().endTag(vtk::fileTag::PIECE);
+ }
+
+ // Finally
+ // could add a 'ghost' level above to visualize extrapolated values
+ // draw as two triangles to distingush from real levels ...
+
+ format().endTag(vtk::fileTag::POLY_DATA);
+ format().endVTKFile();
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointTools.C b/src/lumpedPointMotion/lumpedPointTools.C
new file mode 100644
index 0000000000..fa4e741d01
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointTools.C
@@ -0,0 +1,161 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016-2017 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see
+ static autoPtr loadPointField
+ (
+ const pointMesh::Mesh& mesh,
+ const IOobject* io
+ )
+ {
+ if (io && io->headerClassName() == GeoFieldType::typeName)
+ {
+ Info<< "Reading " << GeoFieldType::typeName
+ << ' ' << io->name() << endl;
+
+ return autoPtr
+ (
+ new GeoFieldType
+ (
+ IOobject
+ (
+ io->name(),
+ io->instance(),
+ io->local(),
+ io->db(),
+ IOobject::MUST_READ,
+ IOobject::AUTO_WRITE,
+ io->registerObject()
+ ),
+ mesh
+ )
+ );
+ }
+
+ return autoPtr();
+ }
+}
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+Foam::List
+Foam::lumpedPointTools::lumpedPointStates(Istream& is)
+{
+ dictionary contents(is);
+ List entries(contents.lookup("response"));
+
+ DynamicList> states(entries.size());
+
+ for (const dictionary& dict : entries)
+ {
+ states.append
+ (
+ lumpedPointStateTuple
+ (
+ readScalar(dict.lookup("time")),
+ lumpedPointState(dict)
+ )
+ );
+ }
+
+ return states.shrink();
+}
+
+
+Foam::List
+Foam::lumpedPointTools::lumpedPointStates(const fileName& file)
+{
+ IFstream is(file);
+ return lumpedPointStates(is);
+}
+
+
+Foam::pointIOField
+Foam::lumpedPointTools::points0Field(const polyMesh& mesh)
+{
+ pointIOField pts
+ (
+ IOobject
+ (
+ "points",
+ mesh.time().constant(),
+ polyMesh::meshSubDir,
+ mesh,
+ IOobject::MUST_READ,
+ IOobject::NO_WRITE,
+ false // Do not re-register
+ )
+ );
+
+ return pts;
+}
+
+
+
+Foam::labelList
+Foam::lumpedPointTools::lumpedPointPatchList(const pointVectorField& pvf)
+{
+ return lumpedPointDisplacementPointPatchVectorField::patchIds(pvf);
+}
+
+
+Foam::labelList Foam::lumpedPointTools::lumpedPointPatchList
+(
+ const polyMesh& mesh
+)
+{
+ IOobjectList objects0(mesh, "0");
+
+ pointMesh pMesh(mesh);
+
+ autoPtr displacePtr = loadPointField
+ (
+ pMesh,
+ objects0.lookup("pointDisplacement")
+ );
+
+ if (!displacePtr.valid())
+ {
+ Info<< "no valid pointDisplacement" << endl;
+ return labelList();
+ }
+
+ return lumpedPointPatchList(displacePtr());
+}
+
+
+// ************************************************************************* //
diff --git a/src/lumpedPointMotion/lumpedPointTools.H b/src/lumpedPointMotion/lumpedPointTools.H
new file mode 100644
index 0000000000..d31c700d28
--- /dev/null
+++ b/src/lumpedPointMotion/lumpedPointTools.H
@@ -0,0 +1,86 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2016 OpenCFD Ltd.
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see lumpedPointStateTuple;
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace lumpedPointTools
+{
+
+//- Load a list of states from an Istream
+List lumpedPointStates(Istream& is);
+
+//- Load a list of states from a file
+List lumpedPointStates(const fileName& file);
+
+
+//- Return the 0 or constant points field
+pointIOField points0Field(const polyMesh& mesh);
+
+//- Return the patch-ids associated a "lumpedPointDisplacement" type
+labelList lumpedPointPatchList(const pointVectorField& pvf);
+
+//- Get the "pointDisplacement" at time 0 and use that to determine which
+// patches have a "lumpedPointDisplacement" type
+labelList lumpedPointPatchList(const polyMesh& mesh);
+
+} // End namespace lumpedPointTools
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //