****mesher

  %=================================================
  % mesh of a plate with circular hole
  %=================================================

 ***mesh plate_circular_hole.geof
  **mesh_quarter_plate_hole
   *ncut_radius 45                  % number of elements in radial-direction
   *ncut_theta  30                  % number of elements in theta-direction
   *size_square 1.00                % plate edge length / 2
   *size_hole 0.2                   % radius a
  %*nb_reg 3                        % number of mesh layers with the same hole roundness
   *quad                            % quadratic mesh
   *geometrical_progression 1.05    % mesh progression in radial direction,
                                    % finer mesh near the hole
   *hole_anisotropy 1.              % aspect ratio b/a , 1 = hole is circular
   *balance_nodes                  

  % the following bset are required for calcul and post-processing
  **bset x0 *use_nset x0
  **bset y0 *use_nset y0
  **bset x1 *use_nset x1
  **bset y1 *use_nset y1
  **rename_set *nsets layer_0  hole_edge
  **bset hole_edge *use_nset hole_edge
  **set_reduced                     % reduced integration

  %=================================================
  % mesh of a plate with anisotropic hole
  %=================================================

 ***mesh plate_elliptic_hole.geof
  **mesh_quarter_plate_hole
   *ncut_radius 45               
   *ncut_theta  30              
   *size_square 1.00           
   *size_hole 0.2             
%  *nb_reg 1                 
   *quad                    
   *geometrical_progression 1.05
   *hole_anisotropy 0.5           
   *balance_nodes                  

  **bset x0 *use_nset x0
  **bset y0 *use_nset y0
  **bset x1 *use_nset x1
  **bset y1 *use_nset y1
  **rename_set *nsets layer_0  hole_edge
  **bset hole_edge *use_nset hole_edge
  **set_reduced      

****return