# solidDisplacementFoam - problem with weakly nonsymmetric result on symmetric geometry

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 October 3, 2012, 13:20 solidDisplacementFoam - problem with weakly nonsymmetric result on symmetric geometry #1 New Member   Join Date: Oct 2012 Location: Czech Republic Posts: 1 Rep Power: 0 Hello. 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) 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 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?

 Tags beam, elastic, soliddisplacementfoam

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Jas Main CFD Forum 10 March 30, 2013 12:26 windmechy Fidelity CFD 9 June 16, 2012 10:00 fatb0y CFX 3 January 14, 2012 19:42 sontukaku Main CFD Forum 2 December 28, 2010 00:44 Fred G. Kang Main CFD Forum 1 October 14, 1998 11:41

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

 Contact Us - CFD Online - Privacy Statement - Top