COMP: avoid ambiguous construct from tmp - solvers/ stressAnalysis

This commit is contained in:
Mark Olesen
2010-12-21 09:53:39 +01:00
parent 078c427594
commit 091baa91d9
5 changed files with 12 additions and 8 deletions

View File

@ -102,7 +102,7 @@ int main(int argc, char *argv[])
} }
{ {
volTensorField gradD = fvc::grad(D); volTensorField gradD(fvc::grad(D));
sigmaD = mu*twoSymm(gradD) + (lambda*I)*tr(gradD); sigmaD = mu*twoSymm(gradD) + (lambda*I)*tr(gradD);
if (compactNormalStress) if (compactNormalStress)

View File

@ -168,7 +168,7 @@ void tractionDisplacementFvPatchVectorField::updateCoeffs()
scalar twoMuLambda = (2*mu + lambda).value(); scalar twoMuLambda = (2*mu + lambda).value();
vectorField n = patch().nf(); vectorField n(patch().nf());
const fvPatchField<symmTensor>& sigmaD = const fvPatchField<symmTensor>& sigmaD =
patch().lookupPatchField<volSymmTensorField, symmTensor>("sigmaD"); patch().lookupPatchField<volSymmTensorField, symmTensor>("sigmaD");
@ -207,7 +207,11 @@ void tractionDisplacementFvPatchVectorField::write(Ostream& os) const
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
makePatchTypeField(fvPatchVectorField, tractionDisplacementFvPatchVectorField); makeNonTemplatedPatchTypeField
(
fvPatchVectorField,
tractionDisplacementFvPatchVectorField
);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

View File

@ -1,10 +1,10 @@
if (!(runTime.timeIndex() % 5)) if (!(runTime.timeIndex() % 5))
{ {
volScalarField kineticEnergy = magSqr(Dcorr); volScalarField kineticEnergy(magSqr(Dcorr);
dimensionedScalar intKineticEnergy = fvc::domainIntegrate(kineticEnergy); dimensionedScalar intKineticEnergy = fvc::domainIntegrate(kineticEnergy);
Info<< "kineticEnergy = " << intKineticEnergy.value() << endl; Info<< "kineticEnergy = " << intKineticEnergy.value() << endl;
volScalarField kineticPower = (Dcorr - Dcorr.oldTime()) & Dcorr; volScalarField kineticPower((Dcorr - Dcorr.oldTime()) & Dcorr);
dimensionedScalar intKineticPower = fvc::domainIntegrate(kineticPower); dimensionedScalar intKineticPower = fvc::domainIntegrate(kineticPower);
Info<< "kineticPower = " << intKineticPower.value() << endl; Info<< "kineticPower = " << intKineticPower.value() << endl;

View File

@ -69,7 +69,7 @@ int main(int argc, char *argv[])
D += accFac*Dcorr; D += accFac*Dcorr;
{ {
volTensorField gradDcorr = fvc::grad(Dcorr); volTensorField gradDcorr(fvc::grad(Dcorr));
sigmaExp = sigmaExp =
(lambda - mu)*gradDcorr + mu*gradDcorr.T() (lambda - mu)*gradDcorr + mu*gradDcorr.T()

View File

@ -165,7 +165,7 @@ void tractionDisplacementCorrectionFvPatchVectorField::updateCoeffs()
lambda = nu*E/((1.0 + nu)*(1.0 - nu)); lambda = nu*E/((1.0 + nu)*(1.0 - nu));
} }
vectorField n = patch().nf(); vectorField n(patch().nf());
const fvPatchField<symmTensor>& sigmaD = const fvPatchField<symmTensor>& sigmaD =
patch().lookupPatchField<volSymmTensorField, symmTensor>("sigmaD"); patch().lookupPatchField<volSymmTensorField, symmTensor>("sigmaD");
@ -194,7 +194,7 @@ void tractionDisplacementCorrectionFvPatchVectorField::write(Ostream& os) const
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
makePatchTypeField makeNonTemplatedPatchTypeField
( (
fvPatchVectorField, fvPatchVectorField,
tractionDisplacementCorrectionFvPatchVectorField tractionDisplacementCorrectionFvPatchVectorField