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

a problem about boundary condition in stressanalysis

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 8, 2011, 05:49
Default a problem about boundary condition in stressanalysis
  #1
Member
 
chenkai
Join Date: May 2010
Location: munich
Posts: 44
Rep Power: 7
xck1986 is on a distinguished road
Hi everyone,

I want to simulate the torsion of a cylinder with OpenFOAM. (see the attached picture).
But I don't know how to set the boundary condition of the outside for the cylinder, because the direction of the stress always changes.
Does anyone have idea how to set this case?
Is it possible to set the moment as boundary condition in OpenFOAM?

Thanks a lot in advance.
Attached Images
File Type: png torsion.png (10.3 KB, 9 views)
xck1986 is offline   Reply With Quote

Old   February 8, 2011, 12:28
Default
  #2
Member
 
MSR CHANDRA MURTHY
Join Date: Mar 2009
Posts: 32
Rep Power: 8
chandramurthy is on a distinguished road
It may be possible with GroovyBC.
chandramurthy is offline   Reply With Quote

Old   February 10, 2011, 06:34
Default
  #3
Member
 
chenkai
Join Date: May 2010
Location: munich
Posts: 44
Rep Power: 7
xck1986 is on a distinguished road
Quote:
Originally Posted by chandramurthy View Post
It may be possible with GroovyBC.
Hi,
thanks a lot for your reply.
But from the two tuturials of stressAnalysis in OpenFOAM, it seems that, in the time 0 folder, only the Displacement (D) boundary condition can be specified; which use the following boundary type:

type tractionDisplacement
traction uniform (0 0 0) //value for traction(stress)
pressure uniform 0
value uniform (0 0 0) //value for Dispalcement

But I want to define the traction uniform (0 0 0) not as a constant but as a variate. That means the vector depends on the position of the node. e.g.traction uniform ( 100*x/(x2+y2) 100*y/(x2+y2) 0 ).
And when I use the GroovyBC, as far as I know, it seems like I can only specify a variate for the Displacement but not for the traction.

Could please give some suggestion?
Or is it possible to have a force or stress or moment as a boundary condition, which can be also defined in the 0 folder?

Thanks again
xck1986 is offline   Reply With Quote

Old   February 10, 2011, 09:42
Default
  #4
Member
 
MSR CHANDRA MURTHY
Join Date: Mar 2009
Posts: 32
Rep Power: 8
chandramurthy is on a distinguished road
Quote:
Originally Posted by xck1986 View Post
But I want to define the traction uniform (0 0 0) not as a constant but as a variate. That means the vector depends on the position of the node. e.g.traction uniform ( 100*x/(x2+y2) 100*y/(x2+y2) 0 ).
Thanks again
the internal field is described as a function of few variables. i think this is definitely possible with http://openfoamwiki.net/index.php/Contrib_funkySetFields.

have a nice day
chandramurthy is offline   Reply With Quote

Old   February 10, 2011, 12:59
Default
  #5
Member
 
chenkai
Join Date: May 2010
Location: munich
Posts: 44
Rep Power: 7
xck1986 is on a distinguished road
Quote:
Originally Posted by chandramurthy View Post
the internal field is described as a function of few variables. i think this is definitely possible with http://openfoamwiki.net/index.php/Contrib_funkySetFields.

have a nice day
Hi,
Thanks a lot for your help.
But it also seems like that a field musst be set in funkySetFields dictory.
And I want to change the vector traction, which is not a field defined in time 0 folder, but a vector defined in the file of D(Displacement).
So it's still a problem for me.
But thanks again!!!
xck1986 is offline   Reply With Quote

Old   February 11, 2011, 03:06
Default
  #6
Member
 
MSR CHANDRA MURTHY
Join Date: Mar 2009
Posts: 32
Rep Power: 8
chandramurthy is on a distinguished road
Quote:
Originally Posted by xck1986 View Post
Hi,
And I want to change the vector traction, which is not a field defined in time 0 folder, but a vector defined in the file of D(Displacement).
So it's still a problem for me.
create a volVectorField say Ts inside createfields.H which similar to U of icoFoam/createfields.H and rebuild the solver. create a new field file 0/Ts
and solve it

hope it will help the situation
chandramurthy is offline   Reply With Quote

Reply

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
Transient Simulation: Boundary Condition Problem Shafiul CFX 7 January 11, 2011 17:40
boundary condition problem shahab.ehsanfar CFX 2 December 31, 2010 09:32
problem about periodic boundary condition in Fluent winnawinna FLUENT 0 December 29, 2010 00:32
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Problem setting boundary condition from file Antonis CFX 1 September 11, 2006 16:53


All times are GMT -4. The time now is 21:11.