36 const DimensionedField<vector, volMesh>& iF
39 fixedGradientFvPatchVectorField(p, iF),
40 traction_(p.size(),
Zero),
41 pressure_(p.size(), 0.0)
53 const DimensionedField<vector, volMesh>& iF,
54 const fvPatchFieldMapper& mapper
57 fixedGradientFvPatchVectorField(tdpvf, p, iF, mapper),
58 traction_(mapper(tdpvf.traction_)),
59 pressure_(mapper(tdpvf.pressure_))
67 const DimensionedField<vector, volMesh>& iF,
68 const dictionary& dict
71 fixedGradientFvPatchVectorField(p, iF),
72 traction_(
"traction", dict, p.size()),
73 pressure_(
"pressure", dict, p.size())
86 fixedGradientFvPatchVectorField(tdpvf),
87 traction_(tdpvf.traction_),
88 pressure_(tdpvf.pressure_)
96 const DimensionedField<vector, volMesh>& iF
99 fixedGradientFvPatchVectorField(tdpvf, iF),
100 traction_(tdpvf.traction_),
101 pressure_(tdpvf.pressure_)
109 const fvPatchFieldMapper& m
112 fixedGradientFvPatchVectorField::autoMap(m);
113 m(traction_, traction_);
114 m(pressure_, pressure_);
124 fixedGradientFvPatchVectorField::rmap(ptf, addr);
127 refCast<const tractionDisplacementFvPatchVectorField>(ptf);
129 traction_.rmap(dmptf.traction_, addr);
130 pressure_.rmap(dmptf.pressure_, addr);
141 const dictionary& mechanicalProperties =
142 db().lookupObject<IOdictionary>(
"mechanicalProperties");
144 const dictionary& thermalProperties =
145 db().lookupObject<IOdictionary>(
"thermalProperties");
148 const fvPatchField<scalar>& rho =
151 const fvPatchField<scalar>& rhoE =
154 const fvPatchField<scalar>& nu =
162 Switch planeStress(mechanicalProperties.lookup(
"planeStress"));
166 lambda = nu*E/((1.0 +
nu)*(1.0 - nu));
167 threeK = E/(1.0 -
nu);
174 const fvPatchField<symmTensor>& sigmaD =
179 (traction_ - pressure_*
n)/rho
183 Switch thermalStress(thermalProperties.lookup(
"thermalStress"));
187 const fvPatchField<scalar>& threeKalpha=
190 const fvPatchField<scalar>&
T =
193 gradient() +=
n*threeKalpha*T/twoMuLambda;
196 fixedGradientFvPatchVectorField::updateCoeffs();
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
fvPatchField< vector > fvPatchVectorField
virtual tmp< Field< Type > > snGrad() const
Return patch-normal gradient.
virtual void autoMap(const fvPatchFieldMapper &)
Map (and resize as needed) from self given a mapping object.
virtual void write(Ostream &) const
Write.
tractionDisplacementFvPatchVectorField(const fvPatch &, const DimensionedField< vector, volMesh > &)
Construct from patch and internal field.
virtual void updateCoeffs()
Update the coefficients associated with the patch field.
virtual void write(Ostream &) const
Write.
Macros for easy insertion into run-time selection tables.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
SymmTensor< scalar > symmTensor
SymmTensor of scalars.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
List< label > labelList
A List of labels.
virtual void rmap(const fvPatchVectorField &, const labelList &)
Reverse map the given fvPatchField onto this fvPatchField.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
void writeEntry(Ostream &os, const HashTable< T, Key, Hash > &ht)
const dimensionedScalar mu
Atomic mass unit.
virtual void operator=(const UList< Type > &)
dimensionedScalar lambda(laminarTransport.lookup("lambda"))
Field< vector > vectorField
Specialisation of Field<T> for vector.
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)