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  rho*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:61
Foam::max
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
Foam::constant::physicoChemical::sigma
const dimensionedScalar & sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
runTime
engineTime & runTime
Definition: createEngineTime.H:13
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
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:57
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18
Foam::dev
dimensionedSymmTensor dev(const dimensionedSymmTensor &dt)
Definition: dimensionedSymmTensor.C:104
Foam::magSqr
dimensioned< scalar > magSqr(const dimensioned< Type > &)
Foam::Info
messageStream Info