47 int main(
int argc,
char *argv[])
51 "Transient segregated finite-volume solver of linear-elastic," 52 " small-strain deformation of a solid body, with optional thermal" 53 " diffusion and thermal stresses" 62 #include "createControls.H" 63 #include "createFields.H" 67 Info<<
"\nCalculating displacement field\n" <<
endl;
76 scalar initialResidual = 0;
107 initialResidual = DEqn.solve().max().initialResidual();
111 divSigmaExp =
fvc::div(DEqn.flux());
135 #include "calculateStress.H"
GeometricField< tensor, fvPatchField, volMesh > volTensorField
constexpr char nl
The newline '\n' character (0x0a)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
dimensionedSymmTensor twoSymm(const dimensionedSymmTensor &dt)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
tmp< GeometricField< Type, faPatchField, areaMesh > > d2dt2(const dimensioned< Type > dt, const faMesh &mesh)
static const Identity< scalar > I
dimensionedScalar lambda("lambda", dimTime/sqr(dimLength), laminarTransport)
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
dimensionedScalar tr(const dimensionedSphericalTensor &dt)
Info<< "Reading field D\"<< endl;volVectorField D(IOobject("D", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);autoPtr< volScalarField > Tptr
const dimensionedScalar mu
Atomic mass unit.
fvMatrix< vector > fvVectorMatrix
messageStream Info
Information stream (stdout output on master, null elsewhere)
const dimensionedScalar & D
Execute application functionObjects to post-process existing results.
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)