
[Sponsors] 
solidDisplacementFoam  problem with weakly nonsymmetric result on symmetric geometry 

LinkBack  Thread Tools  Display Modes 
October 3, 2012, 13:20 
solidDisplacementFoam  problem with weakly nonsymmetric result on symmetric geometry

#1 
New Member

Hello.
I have this problem. I'd like to compute stress in basic school case  beam of prism shape witch square crosssection strained in pulling. Symbolic list of vertices in my blockMeshDict file follows: (L 0 0) (L 0 0) (L T 0) (L T 0) (L 0 T) (L 0 T) (L T T) (L T T), where L is length of the beam and T is thichness of beam and both are positive numbers. I use tractionDisplacement boundary condition type. Nonzero stress as boundary condition is only put on sides (L 0 0) (L T 0) (L T T) (L 0 T) and (L 0 0) (L T 0) (L T T) (L 0 T). On these sides the stress is put in direction outward from body of beam and perpendicular to corresponding side and constant on entire corresponding side. On rest of sides are stresses put equal zero. I derived this case from tutorial case placed in tutorials/stressAnalysis/solidDisplacementFoam/plateHole. In file constant/mechanicalProperties I changed Poisson ratio from 0.3 to 0 and Young's modulus of elasticity from 2e+11 to 1e+10. When I run solidDisplacementFoam solver, the result of displacement D is not symmetric according to plane YZ, it means plane with normal (1 0 0). Near the side with xcoordinate L has a little bit another size than near the side L, although in these areas the displacement should be of same size. What is the reason of that?? Result of stress is good  except for numerical error of constant size on entire body of beam. One solution could be to use symmetryPlane boundary condition type. In this case only half of beam is considered. More precisely: Instead of L in xcocrdinates of vertices the 0 is used and symmetryPlane condition is put on side (0 0 0)(0 T 0)(0 T T)(0 0 T) and only on side with L in xcoordinate of vertices the nonzero stress as boundary condition is used. In this case the results are even exact, without numeric errors which in previous case probably are. So is using symmetryPlane bc type the good solution for above described weak nonsymmetry? What do you think? 

Tags 
beam, elastic, soliddisplacementfoam 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Future CFD Research  Jas  Main CFD Forum  10  March 30, 2013 13:26 
problem to defined the geometry in autogrid 5  windmechy  NUMECA  9  June 16, 2012 10:00 
Problem Importing Geometry ProE to CFX  fatb0y  CFX  3  January 14, 2012 20:42 
Problem in keps turbulence code in complex geometry  sontukaku  Main CFD Forum  2  December 28, 2010 01:44 
Geometry setup for phase change problem  Fred G. Kang  Main CFD Forum  1  October 14, 1998 11:41 