CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions > OpenFOAM CC Toolkits for Fluid-Structure Interaction

[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend"

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

Like Tree134Likes

Closed Thread
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 18, 2012, 07:47
Default
  #41
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by johannes View Post
Hi Philip,

thanks for your fast response.
Luckily, my problem has been a false alarm. I don't know exactly how this happened but somehow during git checking out/compiling/cleaning I managed to kill my OF-extend binaries, so all the errors were just dependency issues. Yesterday I recompiled everything including the solidMechanics branch without a single error so everything is working now.

Best regards,
Johannes
That's good to hear.

Hopefully the solvers will do what you want!

Last edited by bigphil; December 18, 2012 at 07:47. Reason: typo
bigphil is offline  

Old   January 14, 2013, 08:56
Default empty page in source-forge for the new solidMechanics solver
  #42
New Member
 
Tian Tang
Join Date: Jun 2012
Location: Copenhagen, Denmark
Posts: 18
Rep Power: 14
tiat is on a distinguished road
hi, dear Philip,

thanks a lot for sharing these new solvers. Since i found your old post saying that the codes can be manually download from the link, however when i click into it, an error shows: "no such project". I was wondering is this because the solvers have neem moved? may I ask you about where i can found and look at the codes now? Your kind response would greatly appreciated!

Tian



Quote:
Originally Posted by bigphil View Post
Hi bryant,

You can browse and download the code manually from the source-forge here.

What error do you get when you run the git command?

Philip
tiat is offline  

Old   January 14, 2013, 10:06
Default
  #43
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by tiat View Post
hi, dear Philip,

thanks a lot for sharing these new solvers. Since i found your old post saying that the codes can be manually download from the link, however when i click into it, an error shows: "no such project". I was wondering is this because the solvers have neem moved? may I ask you about where i can found and look at the codes now? Your kind response would greatly appreciated!

Tian
Hi Tian,

I think something might be wrong with the git on sourceforge.

The OpenFOAm-Extend project is found here:
http://sourceforge.net/projects/openfoam-extend/develop

And normally it is possible to browse the git code by clicking on the git browse code but it does not seem to be working at the moment.

Hmmn...

Philip
bigphil is offline  

Old   January 14, 2013, 10:36
Default
  #44
New Member
 
Tian Tang
Join Date: Jun 2012
Location: Copenhagen, Denmark
Posts: 18
Rep Power: 14
tiat is on a distinguished road
Hi, dear Philip

thanks for the quick response. yes, you are right, there seems a problem with the git repository of openfoam-extend. I just saw the news from OpenFOAM wiki that:

Currently links to the SVN and GIT repositories of openfoam-extend on Sourceforge are not working because these were taken offline. The reason for that is that ANSYS claims three classes (none of them really essential) of the 1.6-ext are under their copyright. The matter is being worked on. All other services at openfoam-extend are still functional.

hopefully, the problem will be fixed soon so I can look at your solidMechanics solver

regards,
Tian

Quote:
Originally Posted by bigphil View Post
Hi Tian,

I think something might be wrong with the git on sourceforge.

The OpenFOAm-Extend project is found here:
http://sourceforge.net/projects/openfoam-extend/develop

And normally it is possible to browse the git code by clicking on the git browse code but it does not seem to be working at the moment.

Hmmn...

Philip
tiat is offline  

Old   January 16, 2013, 12:20
Default
  #45
Member
 
Join Date: Sep 2012
Posts: 51
Rep Power: 14
Hiroshiman is on a distinguished road
Hi,
I was wondering (and found no answer to this yet) if one could use these solid solvers on an anisotropic materials like composites.
I guess the tricky part is importing the materials properties in the good reference frame but I'm not sure of this. The goal is the simulation of the aging (thermal / time equivalent) of a composite like bended pipe.
Any help would be greatly appreciated !
Regards,

Florian
Hiroshiman is offline  

Old   January 16, 2013, 12:39
Default
  #46
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by Hiroshiman View Post
Hi,
I was wondering (and found no answer to this yet) if one could use these solid solvers on an anisotropic materials like composites.
I guess the tricky part is importing the materials properties in the good reference frame but I'm not sure of this. The goal is the simulation of the aging (thermal / time equivalent) of a composite like bended pipe.
Any help would be greatly appreciated !
Regards,

Florian
Hi Florian,

Actually, I started working on simulation of orthotropic composite materials last year and I have developed a large strain orthotropic finite volume solver (elasticOrthoNonLinULSolidFoam) where the principal fibre directions are general. It requires 9 material properties (3 Young's moduli, 3 Poisson's ratio and 3 shear moduli) and the fibre direction field must be defined.

I have a written a paper on the development of the solver and it is currently under-review.
Once the paper is (hopefully) accepted, then I will add it to the solidMechanics branch.

If you are in a rush for the solver, then you can pm me and maybe we could collaborate.

Best regards,
Philip
bigphil is offline  

Old   January 16, 2013, 14:08
Default
  #47
Member
 
Join Date: Sep 2012
Posts: 51
Rep Power: 14
Hiroshiman is on a distinguished road
Quote:
Originally Posted by bigphil View Post
Hi Florian,

Actually, I started working on simulation of orthotropic composite materials last year and I have developed a large strain orthotropic finite volume solver (elasticOrthoNonLinULSolidFoam) where the principal fibre directions are general. It requires 9 material properties (3 Young's moduli, 3 Poisson's ratio and 3 shear moduli) and the fibre direction field must be defined.

I have a written a paper on the development of the solver and it is currently under-review.
Once the paper is (hopefully) accepted, then I will add it to the solidMechanics branch.

If you are in a rush for the solver, then you can pm me and maybe we could collaborate.

Best regards,
Philip
Thank you for the answer. The thing is, I can't use an anisotropic material (5° angle between "fibers" orientation) but I can start with an isotropic material.

Also, I was wondering, is it complex to modify the elasticThermalSolidFoam into a large deformation solver with also a plastic behaviour ? This may sounds a bit messy but I'm trying to mimic the aging of human tissues (very simple at first) using thermics. The material proprieties would change with temperature in the same way the time does, on long term.
Regards,
Florian
Hiroshiman is offline  

Old   January 17, 2013, 12:54
Default
  #48
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by Hiroshiman View Post
Thank you for the answer. The thing is, I can't use an anisotropic material (5° angle between "fibers" orientation) but I can start with an isotropic material.

Also, I was wondering, is it complex to modify the elasticThermalSolidFoam into a large deformation solver with also a plastic behaviour ? This may sounds a bit messy but I'm trying to mimic the aging of human tissues (very simple at first) using thermics. The material proprieties would change with temperature in the same way the time does, on long term.
Regards,
Florian
Hi Florian,

In principle, it should not be that difficult to develop a large strain thermal-elasto-plastic solver, essentially all the work is done, you just have to merge the solvers elasticThermalSolidFoam, elasticPlasticSolidFoam and elasticNonLinTLSolidFoam.
The momentum equation will have small strain terms, large strain terms, a thermal stress term and a plastic stress term. Also make sure the traction boundary condition has all these terms.
That's the idea, let me know if you have any specific questions about it.

Best regards,
Philip
bigphil is offline  

Old   January 18, 2013, 16:58
Default
  #49
Member
 
Join Date: Sep 2012
Posts: 51
Rep Power: 14
Hiroshiman is on a distinguished road
Thank you very much for your help ! I'll try to do that soon.

Best regards,

Florian
Hiroshiman is offline  

Old   January 19, 2013, 11:32
Default
  #50
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Philip,

I have a small question with respect to your contribution, now that the git is offline.

I have found a formulation to a boundary condition, which requires a fixedCurvature type condition on the pressure, i.e. the second spatial derivative of the pressure normal to the boundary. Such a boundary condition is to my knowledge not available in OF, however, you will need that type of BC, when your are modelling solids, correct?

My question is: Will I be able to find an implementation, which I can modify, in the Solid Mechanics toolbox?

All the best

Niels
ngj is offline  

Old   January 19, 2013, 12:11
Default
  #51
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by ngj View Post
Hi Philip,

I have a small question with respect to your contribution, now that the git is offline.

I have found a formulation to a boundary condition, which requires a fixedCurvature type condition on the pressure, i.e. the second spatial derivative of the pressure normal to the boundary. Such a boundary condition is to my knowledge not available in OF, however, you will need that type of BC, when your are modelling solids, correct?

My question is: Will I be able to find an implementation, which I can modify, in the Solid Mechanics toolbox?

All the best

Niels
Hi Niels,

Hmmnn all the solidMechanics solvers are displacement based and therefore (hydrostatic) pressure is not needed - it can be calculated post processing if required.
But mixed displacement-pressure solvers should not need the second spatial derivative of pressure, actually you don't even need the first derivative as gradU/pressure is uniform for a cell.

Maybe I am not quite understanding the question, but all the solid BCs essentially just need gradU at the boundary and none use grad(gradU).

Best regards,
Philip
bigphil is offline  

Old   January 19, 2013, 12:26
Default
  #52
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Hi Philip,

I probably was not that clear. The boundary condition in question is related to a flow solver and not related to solids.

The reason, I ask this question to you, is that I remember from my classes in civil engineering many years ago that a zero moment at a boundary is related to the second derivative of the displacement. Therefore I thought you might have implemented an abstract base class for this type of boundary condition, which I could "steal" for my purpose.

Nonetheless, I can understand that you do not need those in your formulation, thus I will try to see, if I can find the time to do the work myself.

Thanks,

Niels
ngj is offline  

Old   January 19, 2013, 13:02
Default
  #53
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Hi Niels,

OK I see, unfortunately I can't help.

Good luck!
Philip
bigphil is offline  

Old   January 19, 2013, 13:03
Default
  #54
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 37
ngj will become famous soon enoughngj will become famous soon enough
Thanks anyway
ngj is offline  

Old   February 13, 2013, 16:18
Default
  #55
Member
 
Join Date: Jul 2012
Posts: 31
Rep Power: 14
kpax is on a distinguished road
hey phil,

i am currently testing the icoFsiElasticNonLinULSolidFoam - looks amazing!

just a minor suggestion: i think it's a bit misleading that the maximum number of Aitken iterations (nOuterCorr) is set in fvSolution under "PISO"...
wouldn't it make more sense if that parameter was set in couplingProperties, just like outerCorrTolerance etc.?
kpax is offline  

Old   April 19, 2013, 02:49
Default
  #56
Member
 
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13
SamerAli is on a distinguished road
Hi Phil
i am using IcoFSINonLinUlSolidFoam for simulating fluid structure interaction and it's working very well for large deformations and the updated lagrangian used is perfect.
However i have a question and i am new to OpenFoam, how can i add the temperature equation and the temperature coupling at the interface with this solver??
SamerAli is offline  

Old   April 19, 2013, 06:05
Default
  #57
Super Moderator
 
bigphil's Avatar
 
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,096
Rep Power: 34
bigphil will become famous soon enoughbigphil will become famous soon enough
Quote:
Originally Posted by SamerAli View Post
Hi Phil
i am using IcoFSINonLinUlSolidFoam for simulating fluid structure interaction and it's working very well for large deformations and the updated lagrangian used is perfect.
However i have a question and i am new to OpenFoam, how can i add the temperature equation and the temperature coupling at the interface with this solver??
Hi Samer,

That's good the solver is working for you.

To add temperature to icoFsiNonLinULSolidFoam, you need to add the energy equation for the solid and for the fluid. I think this stuff is already done for the conjugate heat transfer solvers (chtMultiRegionFoam, chtMultiRegionSimpleFoam).
So essentially you have to merge icoFsiNonLinULSolidFoam and chtMultiRegionFoam/chtMultiRegionSimpleFoam.
You can see how the interface temperature coupling is taken care of in the CHT solvers using a derived wall flux boundary condition.

Best regards,
Philip

EDIT: The wiki has a nice account of how to add temperature to the fluid here.
And you can see how temperature is added to the solid by comparing elasticSolidFoam and elasticThermalSolidFoam.
Then the last thing you need is the temperature coupling.

Last edited by bigphil; April 19, 2013 at 06:07. Reason: more info
bigphil is offline  

Old   April 19, 2013, 06:43
Default
  #58
Member
 
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13
SamerAli is on a distinguished road
One more question: if i added the temperature equation and the coupling, do i have to make changes concerning running the solver in parallel??
SamerAli is offline  

Old   April 19, 2013, 06:51
Default
  #59
Senior Member
 
Matthias Voß
Join Date: Mar 2009
Location: Berlin, Germany
Posts: 449
Rep Power: 20
mvoss is on a distinguished road
1# maybe move to new thread
2# o.f. routines will handle the parallelization
3# http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam
mvoss is offline  

Old   April 20, 2013, 15:17
Default
  #60
New Member
 
brant
Join Date: Apr 2013
Posts: 2
Rep Power: 0
brant is on a distinguished road
Hi Bigphil,

I am experienced at 3D modeling and have limited math skills..

I would like to model a hollow sphere or shell with interior resonant modes as well as shell modes, driven either from the interior or exterior... For extra credit I would like to model electric currents in the shell and use a material like iron to show magnetostrictive motions...

Is this at all possible in openFOAM? I would like to know if at least the shell with modes could be modeled before I invest a lot of time into FOAM...

What about plasma/solid interactions?

Thanks for your time,

Brant
brant is offline  

Closed Thread

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
GPU Linear Solvers for OpenFOAM gocarts OpenFOAM Announcements from Other Sources 37 August 17, 2022 15:22
[Virtualization] OpenFOAM oriented tutorial on using VMware Player - support thread wyldckat OpenFOAM Installation 2 July 11, 2012 17:01
New OpenFOAM Forum Structure jola OpenFOAM 2 October 19, 2011 07:55
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 3 September 8, 2010 07:25
OpenFOAM Debian packaging current status problems and TODOs oseen OpenFOAM Installation 9 August 26, 2007 14:50


All times are GMT -4. The time now is 05:23.