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

Mathematical representation of fixedDisplacementZeroShear boundary condition

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

Like Tree1Likes
  • 1 Post By Sargam05

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 22, 2013, 01:11
Default Mathematical representation of fixedDisplacementZeroShear boundary condition
  #1
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hello all,

I am using fixedDisplacementZeroShear (in OpenFOAM-ext) for uniaxial compression loading problem. In the problem, there is a 3D solid cube geometry which is fixed at the bottom surface, where I use fixedDisplacementZeroShear boundary condition and traction force on the top surface.

As I found fixedDisplacementZeroShear boundary condition is type of Directionmixed boundary condition which has fixed value and fixed gradient i.e. combination of Neumann and Dirichlet boundary conditions. I am using this boundary condition in the bottom surface of solid cube. I am wondering about mathematical representation of fixedDisplacementZeroShear boundary condition. Please find the mathematical representation of fixedDisplacementZeroShear boundary condition as attachment.
if l = width= length of the cube .

Am I representing it correctly? I'll be very grateful if one can give me advice.

Thank you
Attached Files
File Type: docx fixedDisplacementZeroShearBC.docx (25.5 KB, 66 views)
amuzeshi likes this.
Sargam05 is offline   Reply With Quote

Old   May 22, 2013, 06:06
Default
  #2
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,086
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Sangeeta,

If your boundary condition is this:
Code:
type fixedDisplacementZeroShear;
value uniform (4 5 6);
And the patch has the unit surface normal n,
then the normal displacement of the patch Un is set to n & U, where U is (4 5 6) and '&' is a dot product.
The shear components on the patch are then set such that the shear traction is zero (so essentially zero gradient).

For example, if your patch is flat and has the unit normal (1 0 0), then the x-displacement Ux (normal displacement) is set to (1 0 0) & ( 4 5 6) which is equal to 4.
The shear gradients are then zero, dUy/dx == dUz/dx == 0.

This means on the patch we know that: Ux = 4, sigma_xy = sigma_xz = 0.
And we don't know anything about the other sigma components.

Hope it helps,
Philip
bigphil is offline   Reply With Quote

Old   May 22, 2013, 09:33
Default
  #3
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hi Philip,

Thank you so much for the explanation. Actually I am using following code at bottom surface:
type fixedDisplacementZeroShear;
I am not mentioning the value uniform (0 0 0), but still getting the answer. Does it is okay to not to mention value uniform in case of fixed surface (i.e fixed displacement without shear)?

Best regards,
Sangeeta
Sargam05 is offline   Reply With Quote

Old   May 22, 2013, 10:13
Default
  #4
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,086
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Sangeeta,

Hmnn actually I think fixedDisplacementZeroShear was meant to give an error if the value is not specified - I suppose it assumes a value of ( 0 0 0 ) if the value is not specified so it is probably OK.
To be safe, I would specify "value uniform ( 0 0 0 );".

Best regards,
Philip
bigphil is offline   Reply With Quote

Old   May 23, 2013, 11:19
Default
  #5
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hi Philip,

Thank you so very much for the information!

Best regards,
Sangeeta
Sargam05 is offline   Reply With Quote

Old   May 23, 2013, 11:35
Default
  #6
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hi Philip,

I have one more concern about boundary condition.

As you mentioned in the http://www.cfd-online.com/Forums/ope...structure.html post:


"fixedDisplacementZeroShear which applies a displacement in the normal direction and enforces zero shear stress in the tangential direction (equivalent to directionMixed). "

It seems that
fixedDisplacementZeroShear boundary condition works as symmetry boundary condition works? When I am using Symmetry and fixedDisplacementZeroShear boundary conditions to run microstructure, fixedDisplacementZeroShear is working well and giving converge solution but symmetry condition is giving diverge solution. I do not know why this is happening. I also checked my geometry and it is okay. As I know symmetry boundary condition assumes mirror symmetry of the geometry and shear stresses are zero in the plane of symmetry. What is the different between fixedDisplacementZeroShear and symmetry conditions?


Best regards,
Sangeeta
Sargam05 is offline   Reply With Quote

Old   May 23, 2013, 11:45
Default
  #7
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,086
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Sangeeta,

Yes fixedDisplacementZeroShear should work as a symmetry plane.
One difference between them is that the shear gradients are always zero for symmetry plane and never changed,
but they shear gradient are constantly being changed in fixedDisplacementZeroShear (in updateCoeffs functions) to enforce the traction to be zero -> this must help convergence in your case.

Actually on a related note, I just realised non-orthogonal correction has not be added to fixedDisplacementZeroShear, this will improve results and possibly convergence.
I will added the correction and send it to you when I have some free time.

Philip
bigphil is offline   Reply With Quote

Old   May 23, 2013, 12:15
Default
  #8
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hi Philip,

Thank you for the quick reply and help! I appreciate it!

Best regards,
Sangeeta
Sargam05 is offline   Reply With Quote

Old   May 24, 2013, 00:49
Default
  #9
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hi Philip,

I have run some more problems with symmetry and fixedDisplacementZeroShear boundary conditions to check the difference between these two. For simple dense geometries both conditions give same answers but not for microstructure (as I mentioned previously). As you mentioned that symmetry condition always consider shear gradient to be zero that mean it also enforces traction to be zero.

As you mentioned:
"shear gradient are constantly being changed in fixedDisplacementZeroShear (in updateCoeffs functions) to enforce the traction to be zero"

Does it mean this boundary condition calculates new shear stress values (with some value not zero) in each step and then use the new calculated value to enforce the traction to be zero? As I understand if shear gradient have some value then it will change traction force with some value not zero (based on formula Traction = sigma . n).

Best regards,
Sangeeta
Sargam05 is offline   Reply With Quote

Old   May 24, 2013, 05:58
Default
  #10
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,086
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Sangeeta,

probably the easiest way to check is to look at the tractions on the surface.

I have attached a post processing utility called surfaceTractions which reads in the stress tensor sigma and calculate three tractions fields for visualisation in ParaView: totalTraction, shearTraction and normalTraction.

Check the shear tractions for both of your cases and see which is zero (or closest to zero).

Philip
Attached Files
File Type: zip surfaceTractions.zip (2.8 KB, 29 views)
bigphil is offline   Reply With Quote

Old   May 25, 2013, 18:50
Default
  #11
Member
 
Sangeeta
Join Date: Jul 2012
Location: Kingston, Canada
Posts: 70
Rep Power: 13
Sargam05 is on a distinguished road
Hi Philip,

Thank you for the valuable suggestions. I have run cases for both the boundary conditions and checked it for shear traction values. I am emailing these cases to you because files are larger in size and I cannot attach these here.

Best regards,
Sangeeta
Sargam05 is offline   Reply With Quote

Old   December 31, 2021, 04:00
Default fixedDisplacementZeroShear
  #12
New Member
 
Nikhil Chitnavis
Join Date: Sep 2019
Location: India
Posts: 16
Rep Power: 6
Nikhil2511 is on a distinguished road
I am trying to run a case of flexible rod using fsiFoam with a initial displacement given to rod is 10 mm. For giving this initial displacement I am using "fixedDisplacementZeroShear" BC for the solid but getting the error of "genericFvPatchField". I am attaching the screenshot of the test case and the error, Please help me.


boundaryField
{
FLAP
{
type fixedDisplacementZeroShear;
value uniform (0.01 0 0);

}
FIXED1
{
type fixedDisplacement;
value uniform (0 0 0);
}
FIXED2
{
type fixedDisplacement;
value uniform (0 0 0);
}
frontAndBackPlanes
{
type empty;
}
}
]



Thank you
Attached Images
File Type: png Screenshot from 2021-12-30 21-08-36.png (44.0 KB, 8 views)
File Type: png Screenshot from 2021-12-31 13-22-19.png (174.9 KB, 8 views)
Nikhil2511 is offline   Reply With Quote

Old   January 5, 2022, 13:58
Default
  #13
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,086
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by Nikhil2511 View Post
I am trying to run a case of flexible rod using fsiFoam with a initial displacement given to rod is 10 mm. For giving this initial displacement I am using "fixedDisplacementZeroShear" BC for the solid but getting the error of "genericFvPatchField". I am attaching the screenshot of the test case and the error, Please help me.


boundaryField
{
FLAP
{
type fixedDisplacementZeroShear;
value uniform (0.01 0 0);

}
FIXED1
{
type fixedDisplacement;
value uniform (0 0 0);
}
FIXED2
{
type fixedDisplacement;
value uniform (0 0 0);
}
frontAndBackPlanes
{
type empty;
}
}
]



Thank you

I guess fixedDisplacementZeroShear is only available in solids4foam, and not with fsiFoam from extend-bazaar.

Philip
bigphil is offline   Reply With Quote

Old   January 7, 2022, 13:46
Default
  #14
New Member
 
Nikhil Chitnavis
Join Date: Sep 2019
Location: India
Posts: 16
Rep Power: 6
Nikhil2511 is on a distinguished road
Hi Philip,
Thankyou for the reply, Can you please suggest how to give initial displacement to the flexible beam using solids4foam.

I am trying to run the case for flexible beam fixed at both ends but before running this FSI case initial displacement is given to the beam. I am attaching the screenshot of the geometry.


Thank you
Attached Images
File Type: png Screenshot from 2022-01-07 23-10-22.png (55.2 KB, 4 views)
Nikhil2511 is offline   Reply With Quote

Old   January 11, 2022, 07:55
Default
  #15
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,086
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by Nikhil2511 View Post
Hi Philip,
Thankyou for the reply, Can you please suggest how to give initial displacement to the flexible beam using solids4foam.

I am trying to run the case for flexible beam fixed at both ends but before running this FSI case initial displacement is given to the beam. I am attaching the screenshot of the geometry.


Thank you
The initial displacement field can be set in 0/D (or 0/DD, depending on the solid model) using setFields or some other method.
bigphil is offline   Reply With Quote

Reply

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
Domain Imbalance HMR CFX 5 October 10, 2016 06:57
Radiation interface hinca CFX 15 January 26, 2014 18:11
Mathematical expression for buoyantPressure boundary condition amwitt OpenFOAM Running, Solving & CFD 3 June 25, 2013 12:08
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00
Setting outlet Pressure boundary condition using CAFFA code Mukund Pondkule Main CFD Forum 0 March 16, 2011 04:23


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