Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend"
 

Greetings to all!

This post has been artificially planted here by a moderator, as part of the effort to keep the news of the solidMechanics tooblox separated from the support questions.

• The news thread: http://www.cfd-online.com/Forums/ope...am-extend.html
• Wiki page that can use your help to improve: http://openfoamwiki.net/index.php/Co...solidMechanics
• Support is given on this thread you are currently reading!

Have fun!

Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend"
 

Hi Philip,

I have OpenFOAM 1.7.x in machine but I have successfully installed feature branch OpenFOAM-1.6-ext and updated the solidMechanics folder. I am using following commands mentioned in Allwmake for compiling solvers

#!/bin/sh
set -X
wmake libso solidModels

After this command it seems compile but at last is showing some error (please see below)


SOURCE=fvPatchFields/fixedRotation/fixedRotationFvPatchVectorField.C ; g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/opt/openfoam171/src/finiteVolume/lnInclude -I/opt/openfoam171/src/meshTools/lnInclude -I/opt/openfoam171/src/VectorN/lnInclude -I/opt/openfoam171/src/triSurface/lnInclude -I/opt/openfoam171/src/finiteArea/lnInclude -IlnInclude -I. -I/opt/openfoam171/src/OpenFOAM/lnInclude -I/opt/openfoam171/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linuxGccDPOpt/fixedRotationFvPatchVectorField.o
fvPatchFields/fixedRotation/fixedRotationFvPatchVectorField.C:33:31: error: RodriguesRotation.H: No such file or directory
fvPatchFields/fixedRotation/fixedRotationFvPatchVectorField.C: In member function ‘virtual void Foam::fixedRotationFvPatchVectorField::updateCoeff s()’:
fvPatchFields/fixedRotation/fixedRotationFvPatchVectorField.C:166: error: ‘RodriguesRotation’ was not declared in this scope
make: *** [Make/linuxGccDPOpt/fixedRotationFvPatchVectorField.o] Error 1

Do you think this error is coming because I do not have OpenFOAM1.6.x installed in my machine?

Thanks in advance.

Regards,
Sangeeta

Hi Sangeeta,

The solidMechanics solvers and library use a number of features (such as GGI) that are only found in the OpenFOAM Extend (not the standard OpenCFD OpenFOAM). Therefore, I would recommend you install OpenFOAM-1.6-ext if you want to use solidMechanics (Note: OpenFOAM-1.6.x is released by OpenCFD and OpenFOAM-1.6-ext is realised by the Extend Project).

It would probably be possible to compile the solvers in standard OpenFOAM but you would have to exclude quite a bit of the functionality. So it is much easier to install OpenFOAM-1.6-ext.

Philip

OpenFOAM does include an implementation of arbitrary mesh interface (a.k.a generalisation grid interface 'GGI'), since version 2.1.0:
http://www.openfoam.org/version2.1.0/ami.php

Hi Philip,

Thanks for the reply. Do I need to first install OpenFOAM 16x for installing and compiling OpenFOAM-ext? Can I install both OpenFOAM 171 and OpenFOAM-ext in same machine?


Regards,
Sangeeta

Yes you can have separate installs of OpenFOAM-1.7.1 and OpenFOAM-1.6-ext on your machine (I think you should be able to find quite a bit of information on the forum about this). And no you don't need to install OpenFOAM-1.6.x before OpenFOAM-1.6-ext as they are separate.

Best regards,
Philip

Hi Philip,

Thanks a lot for the reply.

Sangeeta

Hello Philip
I have met the problem mechy have.
When I run It said that I haven't install git.So I run the command:
And It installed fluently.Then I go into the folder OpenFOAM-1.6-ext and run this command again:
git checkout -b feature/solidMechanics remotes/origin/feature/solidMechanics But the same error came up.
Can you tell me what should I do?The version I am using is OF-1.6-ext.

Thank you very much!
Best Regards!

lg88

Hi lg88,

first OpenFOAM-1.6-ext needs to be downloaded and installed with git:
Then navigate into the OpenFOAM-1.6-ext directory and checkout the solidMechanics branch:
Hope it helps,
Philip

Hello Philip
With this command:
can I install OF-1.6-ext in any version of Ubuntu?like Ubuntu 11.04,11.10 or 12.04?Because I always installed Ubuntu 10.04 first and then intsalled OF-1.6-ext.

Thank you very much!

regards!

lg88

Hi lg88,

Yes it should be possible to install and compile OpenFOAM-1-6-ext on any Ubuntu, assuming that all the necessary dependencies are installed.

You should be able to find lots of information on the forum for installing any version of OpenFOAM on most linux distributions.

Philip

icoFSIElasticNonLinULSolidFoam
 

Hi Philip
I have found that there is a problem of the tutorials of icoFSIElasticNonLinULSolidFoam.
There are not polyMesh.tgz (both fluid or solid).I can see the folder ployMesh but only with a file:boundary inside.
I am bot sure whether it is my problem during download.

regards!

lg88

Hi lg88,

You are right, thanks for the fix.

The problem was that ".tgz" files are excluded from the git repo so they were not added. I have fixed it by extracting the meshes to the correct place.

You should be able to get the changes using:
The code can be browsed here.

If there any other mistakes/bugs, let me know!

Best regards,
Philip

Hi Philip,

thanks for the contribution! I have tested the case hotCylinder and elasticThermalSolidFoam and compared results to that from solidDisplacementFoam. Your solver gives better results.

Is that plausible?

Can you explain whats the reason for that?

See the plot attached.


BTW: I found that in some solvers the include command in Make/options for rheologyModel.H was wrong. I had to change it from
-I../../solidModels/lnInclude
to
-I../solidModels/lnInclude

Regards, Markus.

Hi Markus,

Hmnnn that's interesting, since both solvers are essentially using the same method I would expect them to produce the same results, especially considering the mesh is orthogonal in this case (elasticThermalSolidFoam will be better on non-orthogonal grids as it has boundary non-orthogonal corrections).

If I had to guess, I would say that maybe something is slightly different between the two cases (BCs, properties, fvSchemes, fvSolutions, …).

Thanks, this is the remnants of the old way the code was organised. Which solvers so I can fix it?

Best regards,
Philip

Hi Philip
I am sorry that I can not get any changes with the command:
There is not .gz file in the polyMesh folder still.
And with the link:http://openfoam-extend.git.sourcefor...solidMechanics
I can see other meshes information ,but I can not download them.Can you tell how to download the files there?

Thank you very much!

regards!

lg88

Hi lg88,

The ".tgz" files have been removed altogether and instead the mesh files (points, faces, etc.) are located where they are meant to be in constant/polyMesh (you do not need to run blockMesh or any other mesh command).

So if you checkMesh the fluid or solid case you will see that the mesh is there. Alternatively open the case in ParaView to see the mesh.

Philip

Hi Philip
I deleted the whole folder solidMechanics careless.So I tried this command again:
but it said:How can I fix the solidMechanics ?
regards!

lg88

Hi lg88,

Use this command:
This will return the applications/solvers/solidMechanics directory to what it is like in the last commit of the branch feature/solidMechanics.

Philip

Hi Philip,

thank you for the explanation. Indeed I found a mistake in the setup.

Markus

Hi,
thank you for the development of a large deformation solid solver. I'm planning to use it for fsi in arteries and I was wondering : do you have any idea of the amount of work required to use it on OF 2.1.1 ? (and which functionalities would be lost).
Thanks
Florian

Hi Philip
I always get error when I use the following command in the floder OpenFOAM-1.6-ext:
So where can I download the code package firstly(don't use git but manual download),then compile them?

Thank you!

bryant

Hi Florian,

I don't think it should be too difficult to get it to work in standard OpenFOAM, but I haven't tried.

As regards losing functionality, I don't think anything should be lost with the large strain solver. I recommend you take the solver of interest and take the required bits out of the solidModels library (rheologyModel and boundary conditions - fvPatchFields) and try compile it. It should be possible to get it working, with a little help for the smart forum people.

Best regards,
Philip

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

Thank you for your answer, I will look at it and give you some news if I manage to make it work.

Regards,
Florian

Hi Philip
The error I met just like someone have mentioned before.When I run the following command:
I get this:
I am sorry I still can't download the code the website you give me.Because there is no download button there.I can just see tree|history or other button there.


regards!

bryant

Hi bryant,

OK I believe the git command does not work for you because you did not download OpenFOAM-1.6-ext via git, is this the case?

You can download OpenFOAM-1.6-ext via git using the command:
then moving to this OpenFOAM-1.6-ext directory and executing the command:
will checkout the solidMechanics branch.

Philip

Brilliant tutorials! However...
 

Hi there,

Thanks for the upload, a very good contribution, but I'm having trouble compiling the icoFsiElasticNonLinULSolidFoam code. Please find attached my log file for the compiling process. Any ideas how to fix this? Has anyone got the icoFsiElasticNonLinULSolidFoam code working??

Best Regards,

Steven :)

Attachment 16424

Hi Steven,

The reason for this error is due to a missing lnInclude links directory in $FOAM_SRC/tetDecompositionFiniteElement, so to solve the problem:
This will make the missing lnInclude directory.

The solver should then compile.

Best regards,
Philip

Thank you
 

Dear Philip,

Many thanks for your help! The code is working now.

Best Regards,

Steven

[FIXED] it finally seems to bee working after changing the BC:)
Hi,
after trying some of the tutorials I wanted to test your solver on my cases, especially the large deformation FSI solver. For testing purposes I set up a very simple pipe with a fluid inside (the fluid mesh being generated with snappyHexMesh and the solid via extrudeToRegionMesh).
The BC are :
fluid :
inlet : fixed velocity, zero gradient pressure
outlet : zero gradient velocity, fixed value pressure
wall : movingWallVelocity, zero pressure
and fixedValue (0 0 0) for motionU

solid :
outlet/inlet : fixedDisplacement (0 0 0)
inter/outer wall :
type solidTraction;
nonLinear updatedLagrangian;
traction uniform ( 0 0 0 );
pressure uniform 0;
value uniform (0 0 0);


I also changed the gradient scheme for DU as the code politely asked me from Gauss linear to extendedLeastSquares 0.
The problem is, it can't pass the first iteration and crash in tears with a floating point exception :

I hope I missed something simple.
Both solid and fluid domains can be seen in the attached screenshot.
By the way, what is the acceptable interpolation error ?

Thank you in advanced for you help,

regards,
Florian

Hi Florian,

For future reference, what BC did you have to change to solver your problem?

I am not sure what interpolation error is acceptable, it depends on your application and what you are interested in.
More information can be found on this solver in the conference paper by Tukovic and Jasak here.

Philip

Hi, thank you for the answer.

My fluid BC are :
motionU
p
U

Do have any any if anyone has already implemented 0D BC on the solid part ?

Regards,

Florian

Hi Florian,

I am not sure I understand, what is a 0D BC?

Philip

HI Philip,
I meant using lumped parameters to set-up the outer wall BC e.g. to represent the surrounding environment around an elastic pipe with springs and dampers.

Florian

Hi Florian,

I am not aware of such boundary conditions, but I think they should be straight-forward to implement. The BC would be based on solidTraction and then the traction_ could be calculated each increment based on the displacement using what ever spring model you like.

Philip

Hi !
Thank you for your answer, I'll look at it more closely.

Florian

Hi Philip,

first of all, many thanks to you and your group for sharing the Solid Mechanics Solvers. Must have been quite a bit of work...

I'd love to have a look at the new solvers, but unfortunately, I get build errors when using gcc-4.7.1. Do you have a gcc-47 update available (maybe in form of a seperate branch) or do I have to switch back to gcc-46, at least for now?

Thanks & best regards,
Johannes

Hi Johannes,

Hopefully you will find the solvers useful.

I have not tried gcc-47 yet so I have no experience with it.
Changing back to gcc-46 is probably the quickest way to get things working.

If you post the errors here then we can try to fix them.

Best regards,
Philip

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