|
[Sponsors] |
[stressAnalysis] solidDisplacementFoam displacement |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 1, 2015, 11:04 |
solidDisplacementFoam displacement
|
#1 |
New Member
Anon
Join Date: Mar 2010
Posts: 22
Rep Power: 16 |
Hi Everyone!
I'm running, what I would consider a rather simple analysis where in I am simulating the change in temperature of a large pipe (300k to 270.75k). I would have expected the displacement of the pipe to show a contraction, but the view seems to indicate an expansion of the material. What gives? Did I do something wrong in my setup? Here's my various inputs: T: Code:
dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { left { type symmetryPlane; } right { type fixedValue; value uniform 270.75; } down { type symmetryPlane; } up { type fixedValue; value uniform 270.75; } hole { type fixedValue; value uniform 270.75; } frontAndBack { type empty; } } Code:
dimensions [0 1 0 0 0 0 0]; internalField uniform (0 0 0); boundaryField { left { type symmetryPlane; } right { type tractionDisplacement; traction uniform ( 0 0 0 ); pressure uniform 0; value uniform (0 0 0); } down { type symmetryPlane; } up { type tractionDisplacement; traction uniform ( 0 0 0 ); pressure uniform 0; value uniform (0 0 0); } hole { type tractionDisplacement; traction uniform ( 0 0 0 ); pressure uniform 0; value uniform (0 0 0); } frontAndBack { type empty; } } Code:
rho { type uniform; value 7860; } nu { type uniform; value 0.3; } E { type uniform; value 201e+9; } planeStress yes; Code:
C { type uniform; value 0.000472; } k { type uniform; value 51.4; } alpha { type uniform; value 12.9e-06; } thermalStress yes; 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 yes; } Code:
solvers { "(D|T)" { solver GAMG; tolerance 1e-06; relTol 0.1; smoother GaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 20; agglomerator faceAreaPair; mergeLevels 1; } } stressAnalysis { compactNormalStress no; thermalStress yes; nCorrectors 100000; D 1e-06; T 1e-06; } |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
snappyhexmesh remove blockmesh geometry | philipp1 | OpenFOAM Running, Solving & CFD | 2 | December 12, 2014 10:58 |
[snappyHexMesh] determining displacement for added points | CFDnewbie147 | OpenFOAM Meshing & Mesh Conversion | 1 | October 22, 2013 09:53 |
Expression with time function | kestasposkus | CFX | 28 | July 24, 2011 06:47 |
FSI - Specified Mesh Displacement | Vinzent | CFX | 2 | September 17, 2010 07:09 |
how to use Total Mesh Displacement for making a subject translating and rotating? | snailstb | CFX | 3 | September 16, 2009 04:06 |