[frantov]

Hi Foamers


Im working on a problem using solids4foam but I cant get it to work.


One set of surfaces have a 'negative' pressure and another set of surfaces undergoes deformation. The rest of the domain is fixed.

What are the ideal settings for this case?

So far, I have tried a few different ones but, interestingly the set of surfaces undergoing deformation "does not deform at all".


The image attached is from:

solidProperties:
solidModel unsLinearGeometry;

mechanicalProperties:linearElastic;

The mesh was done using snappyHexMesh using 8 stls.

I have also tried

nonLinearGeometryUpdatedLagrangian;

nonLinearGeometryTotalLagrangianTotalDisplacement; ;


In all these models, the surface that has to be deformed in theory, does not deform at all.

I set up this case in Ansys and I got the deformation I was expecting, therefore Im cleary missing something in the OF setup.


Any help is appreciated, and thanks for the development of these amazing libraries!

[frantov]

Im thinking the problem is the mesh. checkMesh tells me it is a 2D mesh!. Im trying to fix that issue now, anyone has experienced something similar?



$ checkMesh


Checking geometry...
This is a 2-D mesh
Overall domain bounding box (-0.000326 -0.000451 -1e-06) (0.000326 0.000451 0.000451)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
Mesh (non-empty, non-wedge) dimensions 2
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-1.93205e-16 -1.63864e-16 -4.70151e-16) Threshold = 1e-06 OK.
Max cell openness = 2.02888e-16 OK.
Max aspect ratio = 1.99558 OK.
Minumum face area = 1.8419e-10. Maximum face area = 1.47026e-09. Face area magnitudes OK.
Min volume = 4.15348e-15. Max volume = 3.32279e-14. Total volume = 2.25227e-10. Cell volumes OK.
Mesh non-orthogonality Max: 39.3783 average: 5.31682 Threshold = 70
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.657972 OK.

Mesh OK.

[frantov]

Hi All



I was able to solve my problem. The issue was a silly error in blockMesh.

The boundary section was commented out. Hence the background mesh for snappy was ambiguous.



It actually helped me to document what was wrong, so I'll leave the post here.