calculateStress.H
Go to the documentation of this file.
1 if (runTime.writeTime())
2 {
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 }
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:62
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
virtual Ostream & write(const char)=0
Write character.
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
fvMesh & mesh
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:58
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
dimensioned< scalar > magSqr(const dimensioned< Type > &)
messageStream Info