FreeFOAM The Cross-Platform CFD Toolkit
calculateStress.H
Go to the documentation of this file.
1  if (runTime.outputTime())
2  {
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  if (thermalStress)
17  {
18  const volScalarField& T = Tptr();
19  sigma = sigma - I*(rho*threeKalpha*T);
20  }
21 
22  volScalarField sigmaEq
23  (
24  IOobject
25  (
26  "sigmaEq",
27  runTime.timeName(),
28  mesh,
29  IOobject::NO_READ,
30  IOobject::AUTO_WRITE
31  ),
32  sqrt((3.0/2.0)*magSqr(dev(sigma)))
33  );
34 
35  Info<< "Max sigmaEq = " << max(sigmaEq).value()
36  << endl;
37 
38  runTime.write();
39  }
40 
41 // ************************ vim: set sw=4 sts=4 et: ************************ //