Merge branch 'master' of ssh://git.code.sf.net/p/foam-extend/foam-extend-3.2
[foam-extend-3.2.git] / applications / utilities / mesh / manipulation / deformedTetFemGeom / deformedTetFemGeom.C
blob3e61e70b542bf70db8ffb854c08f082506fc6653
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | foam-extend: Open Source CFD
4    \\    /   O peration     | Version:     3.2
5     \\  /    A nd           | Web:         http://www.foam-extend.org
6      \\/     M anipulation  | For copyright notice see file Copyright
7 -------------------------------------------------------------------------------
8 License
9     This file is part of foam-extend.
11     foam-extend is free software: you can redistribute it and/or modify it
12     under the terms of the GNU General Public License as published by the
13     Free Software Foundation, either version 3 of the License, or (at your
14     option) any later version.
16     foam-extend is distributed in the hope that it will be useful, but
17     WITHOUT ANY WARRANTY; without even the implied warranty of
18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
19     General Public License for more details.
21     You should have received a copy of the GNU General Public License
22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
24 Application
25     deformedTetFemGeom
27 Description
28     Deforms a polyMesh using a displacement field U from the tetFem method
29     and a scaling factor supplied as an argument.
31 \*---------------------------------------------------------------------------*/
33 #include "argList.H"
34 #include "polyMesh.H"
35 #include "pointFields.H"
36 #include "tetFemMatrix.H"
37 #include "tetPointFields.H"
38 #include "IStringStream.H"
40 using namespace Foam;
42 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
44 int main(int argc, char *argv[])
46     argList::validArgs.append("scaling factor");
48 #   include "setRootCase.H"
50     scalar scaleFactor(readScalar(IStringStream(args.additionalArgs()[0])()));
52 #   include "createTime.H"
53 #   include "createPolyMesh.H"
55     tetPolyMesh tetMesh(mesh);
57     // Get times list
58     instantList Times = runTime.times();
60     pointField zeroPoints(mesh.points());
62     runTime.setTime(Times[0], 0);
64     for (int i = 1; i<Times.size(); i++)
65     {
66         runTime.setTime(Times[i], i);
68         Info<< "Time = " << runTime.timeName() << endl;
70         IOobject Uheader
71         (
72             "U",
73             runTime.timeName(),
74             mesh,
75             IOobject::MUST_READ
76         );
78         // Check U exists
79         if (Uheader.headerOk())
80         {
81             Info<< "    Reading U" << endl;
82             tetPointVectorField U(Uheader, tetMesh);
84             vectorField Uslice
85             (
86                 vectorField::subField
87                 (
88                     U.internalField(),
89                     zeroPoints.size()
90                 )
91             );
93             pointField newPoints = zeroPoints + scaleFactor*Uslice;
95             mesh.movePoints(newPoints);
97             mesh.write();
98         }
99         else
100         {
101             Info<< "    No U" << endl;
102         }
104         Info<< endl;
105     }
107     Info<< "End\n" << endl;
109     return(0);
113 // ************************************************************************* //