better transfer of sets from gmsh to cgx

Author: mkraska

Linear/Crack1/README.md

@@ -42,7 +42,7 @@ From the difference of the strain energy for two different crack lengths the ene
 release rate can be calculated.
 
 In order to determine the energy release rate, run the model twice with
-different settings of the crack length a in `ct.geo`
+different settings of the crack length `a` in `ct.geo`
 
 For each run, record the crack length a and the strain energy from the file s`olve.dat`.
 
@@ -69,10 +69,9 @@ the element size at the three points in the crack plane.
 <img src="mesh.png" width="400" title="Mesh for B=20, a=20.5, t=2. Image generated by Save as png in Gmsh">
 
 The faces for symmetry constraints and load application are defined as so called physical groups in Gmsh.
-This has the advantage that properly named element sets are written to the output `gmsh.inp`.
+This has the advantage that properly named node sets are written to the output `gmsh.inp`.
 
-The mesh then is read into CGX and the surface element sets are extended to node sets (for the constraints) and to
-element face sets (for the load).
+The mesh then is read into CGX, the surface elements from Gmsh are deleted and the node set for the load is extended into a face set.
 
 <img src="Refs/sets.png" width="400" title="Sets: Load (yellow), Ysym (red), Zsym (blue)">
 

Linear/Crack1/ct.fbd

@@ -1,16 +1,8 @@
 sys gmsh ct.geo
 read gmsh.inp
 
-# Symmetry
-comp Symy do
-comp Symz do
-# force
-comp Load do
 # remove all surface elements
-seta tozap e all
-setr tozap e +C3D10
-zap tozap
-del se0
+zap +CPS6
 
 send all abq
 send Symy abq nam

Linear/Crack1/ct.geo

@@ -55,5 +55,6 @@ Physical Surface("Symy")={30};
 Physical Surface("Symz")={14};
 Physical Surface("Load")={66,70};
 Physical Volume("Part")={1};
+Mesh.SaveGroupsOfNodes = 1;
 Save "gmsh.inp";
 Save "gmsh.png";