surfaceFeatures: New version of surfaceFeatureExtract with simplied controls
Surfaces are specified as a list and the controls applied to each, e.g. in the
rhoPimpleFoam/RAS/annularThermalMixer tutorial:
surfaces
(
"AMI.obj"
"shaft.obj"
"wall.obj"
"statorBlades.obj"
"rotorBlades.obj"
);
includedAngle 150; // Identifes a feature when angle
// between faces < includedAngle
trimFeatures
{
minElem 10; // minimum edges within a feature
}
writeObj yes; // writes out _edgeMesh.obj files to view features
If different controls are required for different surfaces multiple
sub-dictionaries can be used:
AMIsurfaces
{
surfaces
(
"AMI.obj"
);
includedAngle 140; // Identifes a feature when angle
// between faces < includedAngle
trimFeatures
{
minElem 8; // minimum edges within a feature
}
writeObj yes; // writes out _edgeMesh.obj files to view features
}
otherSurfaces
{
surfaces
(
"shaft.obj"
"wall.obj"
"statorBlades.obj"
"rotorBlades.obj"
);
includedAngle 150; // Identifes a feature when angle
// between faces < includedAngle
trimFeatures
{
minElem 10; // minimum edges within a feature
}
writeObj yes; // writes out _edgeMesh.obj files to view features
}
Existing feature edge files corresponding to particular surfaces can be specified using
the "files" association list:
surfaces
(
"AMI.obj"
"shaft.obj"
"wall.obj"
"statorBlades.obj"
"rotorBlades.obj"
);
files
(
"AMI.obj" "constant/triSurface/AMI.obj.eMesh";
);
includedAngle 150; // Identifes a feature when angle
// between faces < includedAngle
trimFeatures
{
minElem 10; // minimum edges within a feature
}
writeObj yes; // writes out _edgeMesh.obj files to view features
This commit is contained in:
Henry Weller
@ -0,0 +1,3 @@
|
||||
surfaceFeatures.C
|
||||
|
||||
EXE = $(FOAM_APPBIN)/surfaceFeatures
|
||||
11
applications/utilities/surface/surfaceFeatures/Make/options
Normal file
11
applications/utilities/surface/surfaceFeatures/Make/options
Normal file
@ -0,0 +1,11 @@
|
||||
EXE_INC = \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/triSurface/lnInclude \
|
||||
-I$(LIB_SRC)/surfMesh/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude
|
||||
|
||||
EXE_LIBS = \
|
||||
-lmeshTools \
|
||||
-ltriSurface \
|
||||
-lsampling
|
||||
729
applications/utilities/surface/surfaceFeatures/surfaceFeatures.C
Normal file
729
applications/utilities/surface/surfaceFeatures/surfaceFeatures.C
Normal file
@ -0,0 +1,729 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2018 OpenFOAM Foundation
|
||||
\\/ 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 <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "argList.H"
|
||||
#include "Time.H"
|
||||
#include "triSurfaceMesh.H"
|
||||
#include "featureEdgeMesh.H"
|
||||
#include "extendedFeatureEdgeMesh.H"
|
||||
#include "surfaceFeatures.H"
|
||||
#include "triSurfaceFields.H"
|
||||
#include "vtkSurfaceWriter.H"
|
||||
#include "IOdictionary.H"
|
||||
|
||||
using namespace Foam;
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
autoPtr<surfaceFeatures> extractFromFile
|
||||
(
|
||||
const fileName& featureEdgeFile,
|
||||
const triSurface& surf,
|
||||
const Switch& geometricTestOnly
|
||||
)
|
||||
{
|
||||
edgeMesh eMesh(featureEdgeFile);
|
||||
|
||||
// Sometimes duplicate edges are present. Remove them.
|
||||
eMesh.mergeEdges();
|
||||
|
||||
Info<< nl << "Reading existing feature edges from file "
|
||||
<< featureEdgeFile << nl
|
||||
<< "Selecting edges purely based on geometric tests: "
|
||||
<< geometricTestOnly.asText() << endl;
|
||||
|
||||
return autoPtr<surfaceFeatures>
|
||||
(
|
||||
new surfaceFeatures
|
||||
(
|
||||
surf,
|
||||
eMesh.points(),
|
||||
eMesh.edges(),
|
||||
1e-6,
|
||||
geometricTestOnly
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
autoPtr<surfaceFeatures> extractFromSurface
|
||||
(
|
||||
const triSurface& surf,
|
||||
const Switch& geometricTestOnly,
|
||||
const scalar includedAngle
|
||||
)
|
||||
{
|
||||
Info<< nl
|
||||
<< "Constructing feature set from included angle "
|
||||
<< includedAngle << nl
|
||||
<< "Selecting edges purely based on geometric tests: "
|
||||
<< geometricTestOnly.asText() << endl;
|
||||
|
||||
return autoPtr<surfaceFeatures>
|
||||
(
|
||||
new surfaceFeatures
|
||||
(
|
||||
surf,
|
||||
includedAngle,
|
||||
0,
|
||||
0,
|
||||
geometricTestOnly
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
autoPtr<surfaceFeatures> surfaceFeatureSet
|
||||
(
|
||||
const fileName& surfaceFileName,
|
||||
const triSurface& surf,
|
||||
const dictionary& dict,
|
||||
const scalar includedAngle
|
||||
)
|
||||
{
|
||||
const Switch geometricTestOnly = dict.lookupOrDefault<Switch>
|
||||
(
|
||||
"geometricTestOnly",
|
||||
"no"
|
||||
);
|
||||
|
||||
if (dict.found("files"))
|
||||
{
|
||||
HashTable<fileName, fileName> fileNames(dict.lookup("files"));
|
||||
|
||||
if (fileNames.found(surfaceFileName))
|
||||
{
|
||||
return extractFromFile
|
||||
(
|
||||
fileNames[surfaceFileName],
|
||||
surf,
|
||||
geometricTestOnly
|
||||
);
|
||||
}
|
||||
else
|
||||
{
|
||||
return extractFromSurface
|
||||
(
|
||||
surf,
|
||||
geometricTestOnly,
|
||||
includedAngle
|
||||
);
|
||||
}
|
||||
}
|
||||
else
|
||||
{
|
||||
return extractFromSurface
|
||||
(
|
||||
surf,
|
||||
geometricTestOnly,
|
||||
includedAngle
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void extractFeatures
|
||||
(
|
||||
const fileName& surfaceFileName,
|
||||
const Time& runTime,
|
||||
const dictionary& dict
|
||||
)
|
||||
{
|
||||
const fileName sFeatFileName = surfaceFileName.lessExt().name();
|
||||
|
||||
Info<< "Surface : " << surfaceFileName << nl << endl;
|
||||
|
||||
const Switch writeVTK =
|
||||
dict.lookupOrDefault<Switch>("writeVTK", "off");
|
||||
const Switch writeObj =
|
||||
dict.lookupOrDefault<Switch>("writeObj", "off");
|
||||
const Switch verboseObj =
|
||||
dict.lookupOrDefault<Switch>("verboseObj", "off");
|
||||
|
||||
const Switch curvature =
|
||||
dict.lookupOrDefault<Switch>("curvature", "off");
|
||||
const Switch featureProximity =
|
||||
dict.lookupOrDefault<Switch>("featureProximity", "off");
|
||||
const Switch closeness =
|
||||
dict.lookupOrDefault<Switch>("closeness", "off");
|
||||
|
||||
|
||||
Info<< nl << "Feature line extraction is only valid on closed manifold "
|
||||
<< "surfaces." << endl;
|
||||
|
||||
// Read
|
||||
// ~~~~
|
||||
|
||||
triSurface surf(runTime.constantPath()/"triSurface"/surfaceFileName);
|
||||
|
||||
Info<< "Statistics:" << endl;
|
||||
surf.writeStats(Info);
|
||||
Info<< endl;
|
||||
|
||||
const faceList faces(surf.faces());
|
||||
|
||||
// Either construct features from surface & featureAngle or read set.
|
||||
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
||||
|
||||
const scalar includedAngle = readScalar(dict.lookup("includedAngle"));
|
||||
|
||||
autoPtr<surfaceFeatures> set
|
||||
(
|
||||
surfaceFeatureSet
|
||||
(
|
||||
surfaceFileName,
|
||||
surf,
|
||||
dict,
|
||||
includedAngle
|
||||
)
|
||||
);
|
||||
|
||||
// Trim set
|
||||
// ~~~~~~~~
|
||||
|
||||
if (dict.isDict("trimFeatures"))
|
||||
{
|
||||
dictionary trimDict = dict.subDict("trimFeatures");
|
||||
|
||||
scalar minLen =
|
||||
trimDict.lookupOrAddDefault<scalar>("minLen", -great);
|
||||
|
||||
label minElem = trimDict.lookupOrAddDefault<label>("minElem", 0);
|
||||
|
||||
// Trim away small groups of features
|
||||
if (minElem > 0 || minLen > 0)
|
||||
{
|
||||
Info<< "Removing features of length < "
|
||||
<< minLen << endl;
|
||||
Info<< "Removing features with number of edges < "
|
||||
<< minElem << endl;
|
||||
|
||||
set().trimFeatures(minLen, minElem, includedAngle);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// Subset
|
||||
// ~~~~~~
|
||||
|
||||
// Convert to marked edges, points
|
||||
List<surfaceFeatures::edgeStatus> edgeStat(set().toStatus());
|
||||
|
||||
if (dict.isDict("subsetFeatures"))
|
||||
{
|
||||
const dictionary& subsetDict = dict.subDict
|
||||
(
|
||||
"subsetFeatures"
|
||||
);
|
||||
|
||||
if (subsetDict.found("insideBox"))
|
||||
{
|
||||
treeBoundBox bb(subsetDict.lookup("insideBox")());
|
||||
|
||||
Info<< "Selecting edges inside bb " << bb;
|
||||
if (writeObj)
|
||||
{
|
||||
Info << " see insideBox.obj";
|
||||
bb.writeOBJ("insideBox.obj");
|
||||
}
|
||||
Info<< endl;
|
||||
|
||||
selectBox(surf, bb, true, edgeStat);
|
||||
}
|
||||
else if (subsetDict.found("outsideBox"))
|
||||
{
|
||||
treeBoundBox bb(subsetDict.lookup("outsideBox")());
|
||||
|
||||
Info<< "Removing all edges inside bb " << bb;
|
||||
if (writeObj)
|
||||
{
|
||||
Info<< " see outsideBox.obj" << endl;
|
||||
bb.writeOBJ("outsideBox.obj");
|
||||
}
|
||||
Info<< endl;
|
||||
|
||||
selectBox(surf, bb, false, edgeStat);
|
||||
}
|
||||
|
||||
const Switch nonManifoldEdges =
|
||||
subsetDict.lookupOrDefault<Switch>("nonManifoldEdges", "yes");
|
||||
|
||||
if (!nonManifoldEdges)
|
||||
{
|
||||
Info<< "Removing all non-manifold edges"
|
||||
<< " (edges with > 2 connected faces) unless they"
|
||||
<< " cross multiple regions" << endl;
|
||||
|
||||
selectManifoldEdges(surf, 1e-5, includedAngle, edgeStat);
|
||||
}
|
||||
|
||||
const Switch openEdges =
|
||||
subsetDict.lookupOrDefault<Switch>("openEdges", "yes");
|
||||
|
||||
if (!openEdges)
|
||||
{
|
||||
Info<< "Removing all open edges"
|
||||
<< " (edges with 1 connected face)" << endl;
|
||||
|
||||
forAll(edgeStat, edgei)
|
||||
{
|
||||
if (surf.edgeFaces()[edgei].size() == 1)
|
||||
{
|
||||
edgeStat[edgei] = surfaceFeatures::NONE;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
if (subsetDict.found("plane"))
|
||||
{
|
||||
const plane cutPlane(subsetDict.lookup("plane")());
|
||||
|
||||
selectCutEdges(surf, cutPlane, edgeStat);
|
||||
|
||||
Info<< "Only edges that intersect the plane with normal "
|
||||
<< cutPlane.normal()
|
||||
<< " and base point " << cutPlane.refPoint()
|
||||
<< " will be included as feature edges."<< endl;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
surfaceFeatures newSet(surf);
|
||||
newSet.setFromStatus(edgeStat, includedAngle);
|
||||
|
||||
Info<< nl
|
||||
<< "Initial feature set:" << nl
|
||||
<< " feature points : " << newSet.featurePoints().size() << nl
|
||||
<< " feature edges : " << newSet.featureEdges().size() << nl
|
||||
<< " of which" << nl
|
||||
<< " region edges : " << newSet.nRegionEdges() << nl
|
||||
<< " external edges : " << newSet.nExternalEdges() << nl
|
||||
<< " internal edges : " << newSet.nInternalEdges() << nl
|
||||
<< endl;
|
||||
|
||||
boolList surfBaffleRegions(surf.patches().size(), false);
|
||||
|
||||
wordList surfBaffleNames;
|
||||
dict.readIfPresent("baffles", surfBaffleNames);
|
||||
|
||||
forAll(surf.patches(), pI)
|
||||
{
|
||||
const word& name = surf.patches()[pI].name();
|
||||
|
||||
if (findIndex(surfBaffleNames, name) != -1)
|
||||
{
|
||||
Info<< "Adding baffle region " << name << endl;
|
||||
surfBaffleRegions[pI] = true;
|
||||
}
|
||||
}
|
||||
|
||||
// Extracting and writing a extendedFeatureEdgeMesh
|
||||
extendedFeatureEdgeMesh feMesh
|
||||
(
|
||||
newSet,
|
||||
runTime,
|
||||
sFeatFileName + ".extendedFeatureEdgeMesh",
|
||||
surfBaffleRegions
|
||||
);
|
||||
|
||||
|
||||
if (dict.isDict("addFeatures"))
|
||||
{
|
||||
const word addFeName = dict.subDict("addFeatures")["name"];
|
||||
Info<< "Adding (without merging) features from " << addFeName
|
||||
<< nl << endl;
|
||||
|
||||
extendedFeatureEdgeMesh addFeMesh
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
addFeName,
|
||||
runTime.time().constant(),
|
||||
"extendedFeatureEdgeMesh",
|
||||
runTime.time(),
|
||||
IOobject::MUST_READ,
|
||||
IOobject::NO_WRITE
|
||||
)
|
||||
);
|
||||
Info<< "Read " << addFeMesh.name() << nl;
|
||||
addFeMesh.writeStats(Info);
|
||||
|
||||
feMesh.add(addFeMesh);
|
||||
}
|
||||
|
||||
|
||||
Info<< nl
|
||||
<< "Final feature set:" << nl;
|
||||
feMesh.writeStats(Info);
|
||||
|
||||
Info<< nl << "Writing extendedFeatureEdgeMesh to "
|
||||
<< feMesh.objectPath() << endl;
|
||||
|
||||
mkDir(feMesh.path());
|
||||
|
||||
if (writeObj)
|
||||
{
|
||||
feMesh.writeObj
|
||||
(
|
||||
feMesh.path()/surfaceFileName.lessExt().name(),
|
||||
verboseObj
|
||||
);
|
||||
}
|
||||
|
||||
feMesh.write();
|
||||
|
||||
// Write a featureEdgeMesh for backwards compatibility
|
||||
featureEdgeMesh bfeMesh
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
surfaceFileName.lessExt().name() + ".eMesh",
|
||||
runTime.constant(),
|
||||
"triSurface",
|
||||
runTime,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE,
|
||||
false
|
||||
),
|
||||
feMesh.points(),
|
||||
feMesh.edges()
|
||||
);
|
||||
|
||||
Info<< nl << "Writing featureEdgeMesh to "
|
||||
<< bfeMesh.objectPath() << endl;
|
||||
|
||||
bfeMesh.regIOobject::write();
|
||||
|
||||
|
||||
// Find distance between close features
|
||||
if (closeness)
|
||||
{
|
||||
Info<< nl << "Extracting internal and external closeness of "
|
||||
<< "surface." << endl;
|
||||
|
||||
// Searchable triSurface
|
||||
const triSurfaceMesh searchSurf
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
sFeatFileName + ".closeness",
|
||||
runTime.constant(),
|
||||
"triSurface",
|
||||
runTime
|
||||
),
|
||||
surf
|
||||
);
|
||||
|
||||
{
|
||||
Pair<tmp<triSurfaceScalarField>> closenessFields
|
||||
(
|
||||
searchSurf.extractCloseness()
|
||||
);
|
||||
|
||||
closenessFields.first()->write();
|
||||
closenessFields.second()->write();
|
||||
|
||||
if (writeVTK)
|
||||
{
|
||||
const faceList faces(searchSurf.faces());
|
||||
|
||||
vtkSurfaceWriter().write
|
||||
(
|
||||
runTime.constantPath()/"triSurface",// outputDir
|
||||
searchSurf.objectRegistry::name(), // surfaceName
|
||||
searchSurf.points(),
|
||||
faces,
|
||||
"internalCloseness", // fieldName
|
||||
closenessFields.first(),
|
||||
false, // isNodeValues
|
||||
true // verbose
|
||||
);
|
||||
|
||||
vtkSurfaceWriter().write
|
||||
(
|
||||
runTime.constantPath()/"triSurface",// outputDir
|
||||
searchSurf.objectRegistry::name(), // surfaceName
|
||||
searchSurf.points(),
|
||||
faces,
|
||||
"externalCloseness", // fieldName
|
||||
closenessFields.second(),
|
||||
false, // isNodeValues
|
||||
true // verbose
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
{
|
||||
Pair<tmp<triSurfacePointScalarField >> closenessFields
|
||||
(
|
||||
searchSurf.extractPointCloseness()
|
||||
);
|
||||
|
||||
closenessFields.first()->write();
|
||||
closenessFields.second()->write();
|
||||
|
||||
if (writeVTK)
|
||||
{
|
||||
const faceList faces(searchSurf.faces());
|
||||
const Map<label>& meshPointMap = searchSurf.meshPointMap();
|
||||
|
||||
const triSurfacePointScalarField&
|
||||
internalClosenessPointField = closenessFields.first();
|
||||
|
||||
const triSurfacePointScalarField&
|
||||
externalClosenessPointField = closenessFields.second();
|
||||
|
||||
scalarField internalCloseness(searchSurf.nPoints(), great);
|
||||
scalarField externalCloseness(searchSurf.nPoints(), great);
|
||||
|
||||
forAll(meshPointMap, pi)
|
||||
{
|
||||
internalCloseness[pi] =
|
||||
internalClosenessPointField[meshPointMap[pi]];
|
||||
|
||||
externalCloseness[pi] =
|
||||
externalClosenessPointField[meshPointMap[pi]];
|
||||
}
|
||||
|
||||
vtkSurfaceWriter().write
|
||||
(
|
||||
runTime.constantPath()/"triSurface",// outputDir
|
||||
searchSurf.objectRegistry::name(), // surfaceName
|
||||
searchSurf.points(),
|
||||
faces,
|
||||
"internalPointCloseness", // fieldName
|
||||
internalCloseness,
|
||||
true, // isNodeValues
|
||||
true // verbose
|
||||
);
|
||||
|
||||
vtkSurfaceWriter().write
|
||||
(
|
||||
runTime.constantPath()/"triSurface",// outputDir
|
||||
searchSurf.objectRegistry::name(), // surfaceName
|
||||
searchSurf.points(),
|
||||
faces,
|
||||
"externalPointCloseness", // fieldName
|
||||
externalCloseness,
|
||||
true, // isNodeValues
|
||||
true // verbose
|
||||
);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
if (curvature)
|
||||
{
|
||||
Info<< nl << "Extracting curvature of surface at the points."
|
||||
<< endl;
|
||||
|
||||
triSurfacePointScalarField k
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
sFeatFileName + ".curvature",
|
||||
runTime.constant(),
|
||||
"triSurface",
|
||||
runTime
|
||||
),
|
||||
surf,
|
||||
dimLength,
|
||||
surf.curvature()
|
||||
);
|
||||
|
||||
k.write();
|
||||
|
||||
if (writeVTK)
|
||||
{
|
||||
vtkSurfaceWriter().write
|
||||
(
|
||||
runTime.constantPath()/"triSurface",// outputDir
|
||||
sFeatFileName, // surfaceName
|
||||
surf.points(),
|
||||
faces,
|
||||
"curvature", // fieldName
|
||||
k,
|
||||
true, // isNodeValues
|
||||
true // verbose
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
if (featureProximity)
|
||||
{
|
||||
Info<< nl << "Extracting proximity of close feature points and "
|
||||
<< "edges to the surface" << endl;
|
||||
|
||||
const scalar searchDistance =
|
||||
readScalar(dict.lookup("maxFeatureProximity"));
|
||||
|
||||
scalarField featureProximity(surf.size(), searchDistance);
|
||||
|
||||
forAll(surf, fi)
|
||||
{
|
||||
const triPointRef& tri = surf[fi].tri(surf.points());
|
||||
const point& triCentre = tri.circumCentre();
|
||||
|
||||
const scalar radiusSqr = min
|
||||
(
|
||||
sqr(4*tri.circumRadius()),
|
||||
sqr(searchDistance)
|
||||
);
|
||||
|
||||
pointIndexHitList hitList;
|
||||
|
||||
feMesh.allNearestFeatureEdges(triCentre, radiusSqr, hitList);
|
||||
featureProximity[fi] = min
|
||||
(
|
||||
feMesh.minDisconnectedDist(hitList),
|
||||
featureProximity[fi]
|
||||
);
|
||||
|
||||
feMesh.allNearestFeaturePoints(triCentre, radiusSqr, hitList);
|
||||
featureProximity[fi] = min
|
||||
(
|
||||
minDist(hitList),
|
||||
featureProximity[fi]
|
||||
);
|
||||
}
|
||||
|
||||
triSurfaceScalarField featureProximityField
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
sFeatFileName + ".featureProximity",
|
||||
runTime.constant(),
|
||||
"triSurface",
|
||||
runTime,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
surf,
|
||||
dimLength,
|
||||
featureProximity
|
||||
);
|
||||
|
||||
featureProximityField.write();
|
||||
|
||||
if (writeVTK)
|
||||
{
|
||||
vtkSurfaceWriter().write
|
||||
(
|
||||
runTime.constantPath()/"triSurface",// outputDir
|
||||
sFeatFileName, // surfaceName
|
||||
surf.points(),
|
||||
faces,
|
||||
"featureProximity", // fieldName
|
||||
featureProximity,
|
||||
false, // isNodeValues
|
||||
true // verbose
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
Info<< endl;
|
||||
}
|
||||
|
||||
|
||||
void extractFeatures
|
||||
(
|
||||
const fileNameList& surfaceFileNames,
|
||||
const Time& runTime,
|
||||
const dictionary& dict
|
||||
)
|
||||
{
|
||||
forAll(surfaceFileNames, i)
|
||||
{
|
||||
extractFeatures(surfaceFileNames[i], runTime, dict);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
argList::addNote
|
||||
(
|
||||
"extract and write surface features to file"
|
||||
);
|
||||
argList::noParallel();
|
||||
|
||||
#include "addDictOption.H"
|
||||
|
||||
#include "setRootCase.H"
|
||||
#include "createTime.H"
|
||||
|
||||
const word dictName("surfaceFeaturesDict");
|
||||
#include "setSystemRunTimeDictionaryIO.H"
|
||||
|
||||
Info<< "Reading " << dictName << nl << endl;
|
||||
|
||||
const IOdictionary dict(dictIO);
|
||||
|
||||
if (dict.found("surfaces"))
|
||||
{
|
||||
extractFeatures
|
||||
(
|
||||
fileNameList(dict.lookup("surfaces")),
|
||||
runTime,
|
||||
dict
|
||||
);
|
||||
}
|
||||
else
|
||||
{
|
||||
forAllConstIter(dictionary, dict, iter)
|
||||
{
|
||||
if (!iter().isDict())
|
||||
{
|
||||
continue;
|
||||
}
|
||||
|
||||
extractFeatures
|
||||
(
|
||||
fileNameList(iter().dict().lookup("surfaces")),
|
||||
runTime,
|
||||
iter().dict()
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
||||
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
||||
<< nl << endl;
|
||||
|
||||
Info<< "End\n" << endl;
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,118 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object surfaceFeatureExtractDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
surface1
|
||||
{
|
||||
surfaces
|
||||
(
|
||||
"surface1.stl"
|
||||
);
|
||||
|
||||
// Mark edges whose adjacent surface normals are at an angle less
|
||||
// than includedAngle as features
|
||||
// - 0 : selects no edges
|
||||
// - 180: selects all edges
|
||||
includedAngle 120;
|
||||
|
||||
// Do not mark region edges
|
||||
geometricTestOnly yes;
|
||||
|
||||
// Write options
|
||||
|
||||
// Write features to obj format for postprocessing
|
||||
writeObj yes;
|
||||
verboseObj no;
|
||||
}
|
||||
|
||||
|
||||
surface2
|
||||
{
|
||||
surfaces
|
||||
(
|
||||
"surface2.nas"
|
||||
);
|
||||
|
||||
// Load from an existing feature edge file
|
||||
files
|
||||
(
|
||||
"surface2.nas" "constant/triSurface/featureEdges.nas";
|
||||
);
|
||||
|
||||
trimFeatures
|
||||
{
|
||||
// Remove features with fewer than the specified number of edges
|
||||
minElem 0;
|
||||
|
||||
// Remove features shorter than the specified cumulative length
|
||||
minLen 0.0;
|
||||
}
|
||||
|
||||
subsetFeatures
|
||||
{
|
||||
// Use a plane to select feature edges
|
||||
// (normal)(basePoint)
|
||||
// Keep only edges that intersect the plane will be included
|
||||
plane (1 0 0)(0 0 0);
|
||||
|
||||
// Select feature edges using a box
|
||||
// (minPt)(maxPt)
|
||||
// Keep edges inside the box:
|
||||
insideBox (0 0 0)(1 1 1);
|
||||
// Keep edges outside the box:
|
||||
outsideBox (0 0 0)(1 1 1);
|
||||
|
||||
// Keep nonManifold edges (edges with >2 connected faces where
|
||||
// the faces form more than two different normal planes)
|
||||
nonManifoldEdges yes;
|
||||
|
||||
// Keep open edges (edges with 1 connected face)
|
||||
openEdges yes;
|
||||
}
|
||||
|
||||
addFeatures
|
||||
{
|
||||
// Add (without merging) another extendedFeatureEdgeMesh
|
||||
name axZ.extendedFeatureEdgeMesh;
|
||||
|
||||
// Optionally flip features (invert all normals, making
|
||||
// convex<->concave etc)
|
||||
//flip false;
|
||||
}
|
||||
|
||||
// Output the curvature of the surface
|
||||
curvature no;
|
||||
|
||||
// Output the proximity of feature points and edges to each other
|
||||
featureProximity no;
|
||||
|
||||
// The maximum search distance to use when looking for other feature
|
||||
// points and edges
|
||||
maxFeatureProximity 1;
|
||||
|
||||
// Out put the closeness of surface elements to other surface elements.
|
||||
closeness no;
|
||||
|
||||
// Write features to obj format for postprocessing
|
||||
writeObj yes;
|
||||
verboseObj no;
|
||||
|
||||
// Write surface proximity and curvature fields to vtk format
|
||||
// for postprocessing
|
||||
writeVTK no;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -1381,7 +1381,8 @@ void Foam::extendedEdgeMesh::flipNormals()
|
||||
|
||||
void Foam::extendedEdgeMesh::writeObj
|
||||
(
|
||||
const fileName& prefix
|
||||
const fileName& prefix,
|
||||
const bool verbose
|
||||
) const
|
||||
{
|
||||
Info<< nl << "Writing extendedEdgeMesh components to " << prefix
|
||||
@ -1389,6 +1390,8 @@ void Foam::extendedEdgeMesh::writeObj
|
||||
|
||||
edgeMesh::write(prefix + "_edgeMesh.obj");
|
||||
|
||||
if (!verbose) return;
|
||||
|
||||
OBJstream convexFtPtStr(prefix + "_convexFeaturePts.obj");
|
||||
Info<< "Writing convex feature points to " << convexFtPtStr.name() << endl;
|
||||
|
||||
@ -1511,19 +1514,15 @@ void Foam::extendedEdgeMesh::writeStats(Ostream& os) const
|
||||
os << incrIndent;
|
||||
os << indent << "convex feature points : "
|
||||
<< setw(8) << concaveStart_-convexStart_
|
||||
//<< setw(8) << convexStart_
|
||||
<< nl;
|
||||
os << indent << "concave feature points : "
|
||||
<< setw(8) << mixedStart_-concaveStart_
|
||||
//<< setw(8) << concaveStart_
|
||||
<< nl;
|
||||
os << indent << "mixed feature points : "
|
||||
<< setw(8) << nonFeatureStart_-mixedStart_
|
||||
//<< setw(8) << mixedStart_
|
||||
<< nl;
|
||||
os << indent << "other (non-feature) points : "
|
||||
<< setw(8) << points().size()-nonFeatureStart_
|
||||
//<< setw(8) << nonFeatureStart_
|
||||
<< nl;
|
||||
os << decrIndent;
|
||||
|
||||
@ -1531,23 +1530,18 @@ void Foam::extendedEdgeMesh::writeStats(Ostream& os) const
|
||||
os << incrIndent;
|
||||
os << indent << "external (convex angle) edges : "
|
||||
<< setw(8) << internalStart_-externalStart_
|
||||
//<< setw(8) << externalStart_
|
||||
<< nl;
|
||||
os << indent << "internal (concave angle) edges : "
|
||||
<< setw(8) << flatStart_-internalStart_
|
||||
//<< setw(8) << internalStart_
|
||||
<< nl;
|
||||
os << indent << "flat region edges : "
|
||||
<< setw(8) << openStart_-flatStart_
|
||||
//<< setw(8) << flatStart_
|
||||
<< nl;
|
||||
os << indent << "open edges : "
|
||||
<< setw(8) << multipleStart_-openStart_
|
||||
//<< setw(8) << openStart_
|
||||
<< nl;
|
||||
os << indent << "multiply connected edges : "
|
||||
<< setw(8) << edges().size()-multipleStart_
|
||||
//<< setw(8) << multipleStart_
|
||||
<< nl;
|
||||
os << decrIndent;
|
||||
}
|
||||
|
||||
@ -520,7 +520,11 @@ public:
|
||||
// Write
|
||||
|
||||
//- Write all components of the extendedEdgeMesh as obj files
|
||||
void writeObj(const fileName& prefix) const;
|
||||
void writeObj
|
||||
(
|
||||
const fileName& prefix,
|
||||
const bool verbose = true
|
||||
) const;
|
||||
|
||||
//- Dump some information
|
||||
virtual void writeStats(Ostream& os) const;
|
||||
|
||||
@ -1,31 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: dev |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object surfaceFeatureExtractDictDefaults;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
extractionMethod extractFromSurface;
|
||||
|
||||
extractFromSurfaceCoeffs
|
||||
{
|
||||
includedAngle 150; // Identifes a feature when angle
|
||||
} // between faces < includedAngle
|
||||
|
||||
trimFeatures
|
||||
{
|
||||
minElem 10; // minimum edges within a feature
|
||||
}
|
||||
|
||||
writeObj yes; // writes out .obj files to view features
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
@ -10,33 +10,26 @@ FoamFile
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object surfaceFeatureExtractDict;
|
||||
object surfaceFeaturesDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
AMI.obj
|
||||
surfaces
|
||||
(
|
||||
"AMI.obj"
|
||||
"shaft.obj"
|
||||
"wall.obj"
|
||||
"statorBlades.obj"
|
||||
"rotorBlades.obj"
|
||||
);
|
||||
|
||||
includedAngle 150; // Identifes a feature when angle
|
||||
// between faces < includedAngle
|
||||
trimFeatures
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
minElem 10; // minimum edges within a feature
|
||||
}
|
||||
|
||||
shaft.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
wall.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
statorBlades.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
|
||||
rotorBlades.obj
|
||||
{
|
||||
#include "surfaceFeatureExtractDictDefaults"
|
||||
}
|
||||
writeObj yes; // writes out _edgeMesh.obj files to view features
|
||||
|
||||
// ************************************************************************* //
|
||||
Reference in New Issue
Block a user