Reaction forces

[daniele.barbera86]

Dear all,

I would like to discuss a bit about the reaction forces that we get from MoFEM. Most of the time I used to calculate the reaction force over a surface, with in mind the idea that I was obtaining a resultant of all the reaction over it. However, in a recent work I am trying to replicating reaction forces are presented from a single node, or point on the geometry. How this change the way of calculating reactions? 

Thank you all

[Karol Lewandowski]

Have you tried creating a block that has only 1 node? I think it should work.

If not, I suggest looking at how a particular module is reading the entities for reaction meshsets and add vertices there as well.

[daniele.barbera86]

Yes I am testing it, by creating a BLOCKSET in cubit and picking the "Node". The rest is the same as before, so the BLOCKSET name is "REACTION". Actually the reaction forces in terminal looks really good!

I also noted that in Multi-field module, we have reactions as well in the VTK file. I remember you told me to do not look at them. It still true? If they are not useful could we only output when we are for example debugging?

Thanks for the hint, maybe we cannot complete the analysis for now but the reactions are matching as well as the plastic strains.

[Karol Lewandowski]

Well done. Yes, all these options can be implemented easily.