2D Plane strain contact model in CCX

1

hi everyone,
I am trying a 2D model from gmsh to CCX, but it seems something wrong with it.
Do anyone has experiences with meshing a 2D model in gmsh, and put this mesh to calculate plane strain contact model in CacluiX(even a simple geometry)?

2

Hi @caeplayer, could you give more details about this? Maybe share the *.geo and *inp files to look at it.

3

In https://github.com/mkraska/CalculiX-Examples/tree/master/Contact/Eyebar there is an example with plane strain and contact. Note that both plane strain and plane stress are actually expanded into a layer of 3D elements with appropriate constraints. I remember that plane stress caused problems in the example because contact may be lost due to thickening of the slave surface (when using node to surface contact).

In plane strain, this should not be a problem. Also, surface to surface contact should not exhibit this problem.

4

hello,@ mkraska, thanks for your suggustion, precisely,my problem is which meshing type 2D or 3D I should choose,and how to change the meshing type to qu4e or element type to CPE4 in gmsh output .inp file to simulate Contact 2D problem no matter plane strain or plane stress.

plan strain
plan strain1307×541 11.9 KB

5

hi@ jbr, my case geometry is too complicated to draw and mesh in cgx ,so I use the 2D step file and Gmsh to mesh it. Eventhough, my 3D simulation is successful,I still do not know how to use Gmsh to make an output file that could simulate 2D plane stress or strain in CCX without cgx pre(I have tried some ways but is still fail).

6

hi,
the answer is actually on your diagram as shown above. if you want to do 2D analysis of plane stress/strain, you need to chose meshing 2D in GMSH then export to Abaqus input files as you noted. but seems it has plane stress by default (e.g CPS3,CPS4) so you need to change manually to plane strain element in CalculiX input decks by simple using find/replace inside text editor.

although 2D plane stress/strain element (triangle & quadrangle) are expanded to 3D element (wedge & hexahedral) it has different due to solver implementation in element nodes displacement restraint for out of plane direction.

1 Like