41 const DimensionedField<vector, volMesh>& iF
44 fixedGradientFvPatchVectorField(p, iF),
45 traction_(p.size(),
vector::zero),
46 pressure_(p.size(), 0.0)
48 fvPatchVectorField::operator=(patchInternalField());
49 gradient() = vector::zero;
53 tractionDisplacementFvPatchVectorField::
54 tractionDisplacementFvPatchVectorField
56 const tractionDisplacementFvPatchVectorField& tdpvf,
58 const DimensionedField<vector, volMesh>& iF,
59 const fvPatchFieldMapper& mapper
62 fixedGradientFvPatchVectorField(tdpvf, p, iF, mapper),
63 traction_(tdpvf.traction_, mapper),
64 pressure_(tdpvf.pressure_, mapper)
68 tractionDisplacementFvPatchVectorField::
69 tractionDisplacementFvPatchVectorField
72 const DimensionedField<vector, volMesh>& iF,
73 const dictionary& dict
76 fixedGradientFvPatchVectorField(p, iF),
77 traction_(
"traction", dict, p.size()),
78 pressure_(
"pressure", dict, p.size())
80 fvPatchVectorField::operator=(patchInternalField());
81 gradient() = vector::zero;
85 tractionDisplacementFvPatchVectorField::
86 tractionDisplacementFvPatchVectorField
88 const tractionDisplacementFvPatchVectorField& tdpvf
91 fixedGradientFvPatchVectorField(tdpvf),
92 traction_(tdpvf.traction_),
93 pressure_(tdpvf.pressure_)
97 tractionDisplacementFvPatchVectorField::
98 tractionDisplacementFvPatchVectorField
100 const tractionDisplacementFvPatchVectorField& tdpvf,
101 const DimensionedField<vector, volMesh>& iF
104 fixedGradientFvPatchVectorField(tdpvf, iF),
105 traction_(tdpvf.traction_),
106 pressure_(tdpvf.pressure_)
112 void tractionDisplacementFvPatchVectorField::autoMap
114 const fvPatchFieldMapper& m
117 fixedGradientFvPatchVectorField::autoMap(m);
118 traction_.autoMap(m);
119 pressure_.autoMap(m);
123 void tractionDisplacementFvPatchVectorField::rmap
129 fixedGradientFvPatchVectorField::rmap(ptf, addr);
131 const tractionDisplacementFvPatchVectorField& dmptf =
132 refCast<const tractionDisplacementFvPatchVectorField>(ptf);
134 traction_.rmap(dmptf.traction_, addr);
135 pressure_.rmap(dmptf.pressure_, addr);
139 void tractionDisplacementFvPatchVectorField::updateCoeffs()
146 const dictionary& mechanicalProperties =
147 db().lookupObject<IOdictionary>(
"mechanicalProperties");
149 const dictionary& thermalProperties =
150 db().lookupObject<IOdictionary>(
"thermalProperties");
161 Switch planeStress(mechanicalProperties.lookup(
"planeStress"));
165 lambda =
nu*E/((1.0 +
nu)*(1.0 -
nu));
166 threeK = E/(1.0 -
nu);
169 scalar twoMuLambda = (2*mu + lambda).value();
173 const fvPatchField<symmTensor>& sigmaD =
178 (traction_ + pressure_*n)/
rho.value()
182 Switch thermalStress(thermalProperties.lookup(
"thermalStress"));
189 const fvPatchField<scalar>&
T =
192 gradient() += n*threeKalpha.value()*T/twoMuLambda;
195 fixedGradientFvPatchVectorField::updateCoeffs();
199 void tractionDisplacementFvPatchVectorField::write(Ostream& os)
const
201 fvPatchVectorField::write(os);
202 traction_.writeEntry(
"traction", os);
203 pressure_.writeEntry(
"pressure", os);
204 writeEntry(
"value", os);