FreeFOAM The Cross-Platform CFD Toolkit
createFields.H
Go to the documentation of this file.
1 Info<< "Reading displacement field D\n" << endl;
3 (
4  IOobject
5  (
6  "D",
7  runTime.timeName(),
8  mesh,
9  IOobject::MUST_READ,
10  IOobject::AUTO_WRITE
11  ),
12  mesh
13 );
14 
15 Info<< "Reading/setting displacement correction field Dcorr\n" << endl;
16 volVectorField Dcorr
17 (
18  IOobject
19  (
20  "Dcorr",
21  runTime.timeName(),
22  mesh
23  ),
24  D
25 );
26 Dcorr *= 0.0;
27 
28 Info<< "Calculating stress field sigmaD\n" << endl;
29 volSymmTensorField sigmaD
30 (
31  IOobject
32  (
33  "sigmaD",
34  runTime.timeName(),
35  mesh,
36  IOobject::NO_READ,
37  IOobject::NO_WRITE
38  ),
39  mu*twoSymm(fvc::grad(D)) + (lambda*I)*tr(fvc::grad(D))
40 );
41 
42 volTensorField sigmaExp
43 (
44  IOobject
45  (
46  "sigmaExp",
47  runTime.timeName(),
48  mesh,
49  IOobject::NO_READ,
50  IOobject::NO_WRITE
51  ),
52  (lambda - mu)*fvc::grad(Dcorr) + mu*fvc::grad(Dcorr)().T()
53  + (lambda*I)*tr(fvc::grad(Dcorr))
54 );
55 
56 // ************************ vim: set sw=4 sts=4 et: ************************ //