CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   bug in solidDisplacementFoam and other structural codes (http://www.cfd-online.com/Forums/openfoam-bugs/114058-bug-soliddisplacementfoam-other-structural-codes.html)

mkraposhin March 4, 2013 12:48

bug in solidDisplacementFoam and other structural codes
 
Dear Foamers, after inspecting solvers like solidDisplacementFoam (and similar) i found, that increase of pressure on surface leads to expansion of material!

This error have roots in patch description: surface normal pressure is added to traction (but it must be subtracted), let's see to tractionDisplacementFvPatchVectorField.C

Code:

    gradient() =
    (
        (traction_ + pressure_*n)/rho
      + twoMuLambda*fvPatchField<vector>::snGrad() - (n & sigmaD)
    )/twoMuLambda;

But it must be:

Code:

    gradient() =
    (
        (traction_ - pressure_*n)/rho
      + twoMuLambda*fvPatchField<vector>::snGrad() - (n & sigmaD)
    )/twoMuLambda;

I think, that this error must be eliminated in next versions

----
And of course, excuse me my English


All times are GMT -4. The time now is 07:22.