CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Stress analysis--BC's on non-planar surface

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 24, 2014, 21:08
Default Stress analysis--BC's on non-planar surface
  #1
New Member
 
Join Date: Oct 2011
Location: Sydney, NSW
Posts: 20
Rep Power: 5
the.drizzle is on a distinguished road
Hi there,

I'm having trouble figuring out a way to prescribe a stress condition on a non planar boundary. For example, if we consider the Plate With A Hole example from the documents and wish to prescribe a pressure condition of 10kPa pushing outward on the hole instead of the (planar) plate boundary, how can I actually do this? In particular, my problem arises from the fact that I cannot see how to employ the tractionDisplacement boundary condition to do this:

Code:
right 
{ 
    type            tractionDisplacement; 
    traction        uniform ( 10000 0 0 ); 
    pressure        uniform 0; 
    value           uniform (0 0 0); 
}
the "traction" is clearly expected to have the a vector value prescribed, but I'm not seeing something like "normal to the surface" available.

Thanks in advance for all help!
the.drizzle is offline   Reply With Quote

Old   June 25, 2014, 10:22
Default pressure
  #2
New Member
 
Jeff
Join Date: May 2014
Location: Nantes, France
Posts: 10
Rep Power: 2
MrJeff is on a distinguished road
Hi the.drizzle,

I am a newbie on OpenFOAM because I use it for 1 month. I use sometimes the solver solidEquilibriumDisplacementFoam and I think traction and pressure are in Pa, traction runs in the x,y, and z directions and pressure runs in the normal direction on the plane.

Try this:
hole { type tractionDisplacement; traction uniform ( 0 0 0 ); pressure uniform 10000; value uniform (0 0 0); }
Good luck! I wait your results.

Jeff.
MrJeff is offline   Reply With Quote

Old   June 25, 2014, 20:07
Default
  #3
New Member
 
Join Date: Oct 2011
Location: Sydney, NSW
Posts: 20
Rep Power: 5
the.drizzle is on a distinguished road
Hi Jeff,

Thanks for the reply! I don't have my OF-equipped computer handy, but I don't think that will work--I think I already tried exactly that approach already. I think. Maybe not though, I will try again later today to be sure.

It would make sense if it did work though, as the we have from the tutorial:

Quote:
The other patches are traction boundary conditions, set by a specialist traction boundary type. The traction boundary conditions are specified by a linear combination of: (1) a boundary traction vector under keyword traction; (2) a pressure that produces a traction normal to the boundary surface that is defined as negative when pointing out of the surface, under keyword pressure. The up and hole patches are zero traction so the boundary traction and pressure are set to zero. For the right patch the traction should be (1e4,0, 0) Pa and the pressure should be 0 Pa.
except even from this, I'm not quite seeing how this works as it should also work (by this reasoning) to define the BC on the right wall as:

Code:
right
{ 
    type            tractionDisplacement; 
    traction        uniform ( 0 0 0 ); 
    pressure        uniform 10000; 
    value           uniform (0 0 0); 
}
which I'm quite sure won't work...
the.drizzle is offline   Reply With Quote

Old   June 26, 2014, 04:02
Default 3 ways to model your exemple and not the same results
  #4
New Member
 
Jeff
Join Date: May 2014
Location: Nantes, France
Posts: 10
Rep Power: 2
MrJeff is on a distinguished road
Hi the.drizzle,

Good to hear from you!

As to model the traction of the right planar surface, I try three differents BC details with the exemple of the plateHole.

Ont the first one, we applied a traction of 10000 Pa on the positive direction.
Code:
right
    {
        type            tractionDisplacement;
        traction        uniform ( 10000 0 0 );
        pressure        uniform 0;
        value           uniform (0 0 0);
    }


The second one:
Code:
right
    {
        type            tractionDisplacement;
        traction        uniform ( 0 0 0 );
        pressure        uniform 10000;
        value           uniform (0 0 0);
    }
The same results from the first example but sigmaXX(1st ex) = - sigmaXX(2nd ex)
I thought the pressure in on the opposite direction. That's why I try a last example on the other direction:

Code:
right
    {
        type            tractionDisplacement;
        traction        uniform ( 0 0 0 );
        pressure        uniform -10000;
        value           uniform (0 0 0);
    }
This time, the results aren't the same as the 1st example but there are closed and not in the opposite values.

What's wrong ? I don't know ...

Please find attached the pics from my tests on OpenFOAM.

See you,

Jeff.
Attached Images
File Type: jpg sigmaXX1st.jpg (14.1 KB, 5 views)
File Type: jpg sigmaXX2nd.jpg (14.2 KB, 4 views)
File Type: jpg sigmaXX3rd.jpg (14.3 KB, 4 views)
MrJeff is offline   Reply With Quote

Old   June 26, 2014, 04:08
Default A mistake from me :D
  #5
New Member
 
Jeff
Join Date: May 2014
Location: Nantes, France
Posts: 10
Rep Power: 2
MrJeff is on a distinguished road
In the third example, I stop my atteinment at the first increment. That's why I haven't the same results as the first exemple .

Please see the real result of the sigmaXX on the third one. It is the same as the first one.

So, the pressure is like a traction if you apply it on the opposite direction.

Now, I'm waiting your results on the hole.
Attached Images
File Type: jpg sigmaXX3rd.jpg (14.6 KB, 5 views)
MrJeff is offline   Reply With Quote

Reply

Tags
non-planar bend, tractiondisplacement

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Cylindrical Quadric Surface for Post-process analysis Steady State FLUENT 0 November 7, 2013 10:03
Plane stress and plane strain cases for stress analysis Sargam05 OpenFOAM 0 May 3, 2013 18:16
Calculating the Wall Shear Stress Gradient over surface 123catty456 Main CFD Forum 1 October 1, 2012 22:27
Layers don't fully surround surface EVBUCF OpenFOAM Native Meshers: snappyHexMesh and Others 14 August 20, 2012 04:31
stress analysis rt Main CFD Forum 0 June 7, 2007 05:01


All times are GMT -4. The time now is 02:58.