CFD Online Discussion Forums - [stressAnalysis] Thermal stress of a beam

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM CC Toolkits for Fluid-Structure Interaction (https://www.cfd-online.com/Forums/openfoam-cc-toolkits-fluid-structure-interaction/)

- - [stressAnalysis] Thermal stress of a beam (https://www.cfd-online.com/Forums/openfoam-cc-toolkits-fluid-structure-interaction/128031-thermal-stress-beam.html)

Thermal stress of a beam

Hi,
I'd like to calculate the stress applied on a beam due to the temperature increasing.
The beam is 1m long and diametere is 0.1m. The displacement is neglacted all over the surface of the beam.
The problem seems simple, but solver doesn't find a solution.
Any suggestion?

Attached are the input files (except for the mesh).

Does anyone say me if there is something wrong in the fvSolution file?

Code:



solvers

{

    "(D|T)"

    {

        solver          GAMG;

        tolerance       1e-06;

        relTol          0.9;

        smoother        GaussSeidel;

        cacheAgglomeration true;

        nCellsInCoarsestLevel 20;

        agglomerator    faceAreaPair;

        mergeLevels     1;

    }

}



stressAnalysis

{

    compactNormalStress no;

    thermalStress yes;

    nCorrectors 100000;

    D 1e-06;

    T 1e-06;

}

or fvSchemes?

Code:

d2dt2Schemes

{

    default         steadyState;

}



ddtSchemes

{

    default         Euler;

}



gradSchemes

{

    default         leastSquares;

    grad(D)         leastSquares;

    grad(T)         leastSquares;

}



divSchemes

{

    default         none;

    div(sigmaD)     Gauss linear;

}



laplacianSchemes

{

    default         none;

    laplacian(DD,D) Gauss linear corrected;

    laplacian(DT,T) Gauss linear corrected;

}



interpolationSchemes

{

    default         linear;

}



snGradSchemes

{

    default         none;

}



fluxRequired

{

    default         no;

    D               yes;

    T               no;

}