[Maimouna]

Dear Nicolas,

lots of thanks for your last post.

[HosseinB]

Hi all,

Thanks in advance for helping with this.

Can OpenFOAM, or other open source derivatives or branches of OpenFOAM, model the simultaneous interaction of a fluid causing large plastic deformations in a solid? Imagine a water flowing around a very thin metal film causing the film to permanently deform even if the flow is stopped.

Thanks,
Hossein

[bigphil]

Originally Posted by HosseinB

Hi all,

Thanks in advance for helping with this.

Can OpenFOAM, or other open source derivatives or branches of OpenFOAM, model the simultaneous interaction of a fluid causing large plastic deformations in a solid? Imagine a water flowing around a very thin metal film causing the film to permanently deform even if the flow is stopped.

Thanks,
Hossein

Hi Hossein,

I have some work in progress code that will work; if you PM me your email address then I will send you details.

The large strain elasto-plastic solid solver is described at:
A Lagrangian Cell-Centred Finite Volume Method for Metal Forming Simulation

Best,
Philip

[mishaparem]

Hello all,

I'm also looking for a way to model deformable surfaces (cloth) in fluid (air) to see how it deforms. I do have some (more) solid elements that would keep the cloth from flapping around wildly though. I'm feeling pretty pessimistic that I can find something that will do what I need, but I'm dreading the idea of writing up the code myself, so... here's to hoping someone else had the same requirements as me and already wrote something up.

Best regards,

Misha

[bigphil]

Originally Posted by mishaparem

Hello all,

I'm also looking for a way to model deformable surfaces (cloth) in fluid (air) to see how it deforms. I do have some (more) solid elements that would keep the cloth from flapping around wildly though. I'm feeling pretty pessimistic that I can find something that will do what I need, but I'm dreading the idea of writing up the code myself, so... here's to hoping someone else had the same requirements as me and already wrote something up.

Best regards,

Misha

Hi Misha,

For modelling a cloth, it will be thin and so will probably require a shell-element approach (even a membrane-style approach); currently, there is no out-of-the-box solver that I am aware of but I can point you in the right direction if you head down the dreaded code development path.

Philip

[mishaparem]

Originally Posted by bigphil

Hi Misha,

For modelling a cloth, it will be thin and so will probably require a shell-element approach (even a membrane-style approach); currently, there is no out-of-the-box solver that I am aware of but I can point you in the right direction if you head down the dreaded code development path.

Philip

I can't say for sure that I will have the time and patience to head down the dreaded code development path, but just in case - if you can point me down the right direction, I would very much appreciate it.

Best regards,

Misha

[bigphil]

Hi Misha,

I suggest having a look at the following articles:

• J. G. Teng, S. F. Chen, and J. L. Hu. A finite volume method for deformation analysis of woven fabrics. International Journal for Numerical Methods in Engineering, 46:2061–2098, 1999.
• S. F Chen, J. L Hu, and J. G Teng. A finite-volume method for contact drape simulation of woven fabrics and garments. Finite Elements in Analysis and Design, 37:513–531, 2001.

These present finite volume methodologies for draping fabrics.

As regards the implementation in OpenFOAM, I think the finite area framework would be a good starting point.

Best,
Philip

[mishaparem]

Sweet! Thank you!

[aow]

Thanks for a nice thread documenting all of your contact work. It was very informative.

In case anyone cares, the previously provided link for Philip Cardiff's PhD thesis is no longer working. Use the following one instead.

https://researchrepository.ucd.ie/handle/10197/7919

[andrewbickerdike]

Hi Philip,

Did you manage to develop the code for the elastic membrane you mentioned?

I am looking at the deformability of red blood cells and think the elastic membrane approach is a good place to start.

Do you think it would be possible to implement such applications into the immersed boundary method?

Best wishes,
Andrew

[bigphil]

Originally Posted by andrewbickerdike

Hi Philip,

Did you manage to develop the code for the elastic membrane you mentioned?

I am looking at the deformability of red blood cells and think the elastic membrane approach is a good place to start.

Do you think it would be possible to implement such applications into the immersed boundary method?

Best wishes,
Andrew

Hi Andrew,

You can find a finite area implementation of a Kirchhoff plate and related tutorials in solids4foam:

/Users/philipc/foam/philipc-4.0/solids4foam/src/solids4FoamModels/solidModels/kirchhoffPlateSolid
/Users/philipc/foam/philipc-4.0/solids4foam/tutorials/solids/beamsPlatesShells

It is essentially a proof of concept.

Interestingly also, OpenFOAM-v2012 also included a Kirchhoff plate implementation using the finite area method, see https://www.openfoam.com/news/main-n...er-and-physics. Apart from reading this page, I haven't looked into this in any more detail.

The application of these methods to non-planar geometry may require extensions to the formulation (I have not looked into it).

Combining this with the immersed boundary method (or even the overset method) should definitely work in theory but no doubt will require for coding. Essentially it is just the fluid model that needs to know about the immersed boundary (or overset) mesh, e.g. see:

 /Users/philipc/foam/philipc-4.0/solids4foam/src/solids4FoamModels/fluidModels/pimpleOversetFluid

[andrewbickerdike]

Hi,

This is very helpful!

Thanks
Andrew

[Destouches]

Hi all,

These solidMechanics solvers look outstanding, thank you all for this thread.

Nevertheless, whatever I do on the elasticFoam/slidingFrictionBall, or the anneau_plaque case you shared (both run well on foam-extended-5.0), the VTK I extract never moves in paraview...

SimaEq rises, so contacts are registrered, but nothing moves !

Is there a trick or something ?

[bigphil]

Originally Posted by Destouches

Hi all,

These solidMechanics solvers look outstanding, thank you all for this thread.

Nevertheless, whatever I do on the elasticFoam/slidingFrictionBall, or the anneau_plaque case you shared (both run well on foam-extended-5.0), the VTK I extract never moves in paraview...

SimaEq rises, so contacts are registrered, but nothing moves !

Is there a trick or something ?

Hi Adrien,

You need to use the "Warp By Vector" filter in ParaView with the displacement field to move the mesh to the deformed position.

Philip

[Destouches]

Originally Posted by bigphil

Hi Adrien,

You need to use the "Warp By Vector" filter in ParaView with the displacement field to move the mesh to the deformed position.

Philip

Thank you very much !
It is now working perfectly.



One side question.
I wanted to fund OF foundation, ESI or extended (I essentially use the second one).


Foundation gives you the opportunity to give 5 000€, which is quite...High...And ESI didn't answer to this enquiry, like they didn't care.


According to you, what is the best way to do it ?


Cheers

[bigphil]

Originally Posted by Destouches

Thank you very much !
It is now working perfectly.



One side question.
I wanted to fund OF foundation, ESI or extended (I essentially use the second one).


Foundation gives you the opportunity to give 5 000€, which is quite...High...And ESI didn't answer to this enquiry, like they didn't care.


According to you, what is the best way to do it ?


Cheers

As you have seen, OpenFOAM.org clearly state how to contribute funds. I am not sure if OpenFOAM.com or Extend have similar options, but in general most people are happy to accept contributions You can always send a follow up email.

[Cinek_Poland]

What plasticity models does OpenFoam have ?