calculateStress.H
Go to the documentation of this file.
1 if (runTime.writeTime())
2 {
3  volSymmTensorField sigma
4  (
5  IOobject
6  (
7  "sigma",
8  runTime.timeName(),
9  mesh,
10  IOobject::NO_READ,
11  IOobject::AUTO_WRITE
12  ),
13  sigmaD
14  );
15 
16  volScalarField sigmaEq
17  (
18  IOobject
19  (
20  "sigmaEq",
21  runTime.timeName(),
22  mesh,
23  IOobject::NO_READ,
24  IOobject::AUTO_WRITE
25  ),
26  sqrt((3.0/2.0)*magSqr(dev(sigma)))
27  );
28 
29  Info<< "Max sigmaEq = " << max(sigmaEq).value()
30  << endl;
31 
32  runTime.write();
33 }
Foam::volSymmTensorField
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:62
Foam::max
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
Definition: fvMeshStitchersMoving.C:102
Foam::Ostream::write
virtual Ostream & write(const char)=0
Write character.
Foam::sqrt
dimensionedScalar sqrt(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:142
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
mesh
fvMesh & mesh
Definition: createMeshesPostProcess.H:9
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:58
Foam::constant::physicoChemical::sigma
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
Foam::dev
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
Definition: dimensionedSymmTensor.C:104
Foam::magSqr
dimensioned< scalar > magSqr(const dimensioned< Type > &)
Foam::Info
messageStream Info