CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

solidDisplacementFoam - problem with weakly nonsymmetric result on symmetric geometry

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   October 3, 2012, 13:20
Default solidDisplacementFoam - problem with weakly nonsymmetric result on symmetric geometry
New Member
Join Date: Oct 2012
Location: Czech Republic
Posts: 1
Rep Power: 0
Sergejevicz is on a distinguished road
Send a message via ICQ to Sergejevicz

I have this problem. I'd like to compute stress in basic school case - beam of prism shape witch square cross-section 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)


(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 x-coordinate -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 x-cocrdinates 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 x-coordinate 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?
Sergejevicz is offline   Reply With Quote


beam, elastic, soliddisplacementfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

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

All times are GMT -4. The time now is 12:49.