Using command SHPE

[dichtstoff]

hello
i have questions with the command SPHE for creating a sphere by cgx.
I’m working on a sphere model and i used before for my surfaces
only default coons shape of surfaces, BLEND instead of SPHE.

Buckling of sphere

Hello Dichtstoff, you used as mentioned only the default coons shape of surfaces. Please regard the command SHPE which is used in the “SPHERE” example on Page 180 of the latest cgx manual. This command temporally creates a NURBS and relates the given surface to it which in this case defines a sphere.

Is it possible to create surfaces with SPHE with mouse and interactive commands q ?
(what is the abbreviations q like qsurf,qshp, ect.?)
or only manual per script?

if i use the SPHE command, i get a limit with the amount of elements,
and the mesh is getting erroneous!?
what is the reaon for?
these happens for meshing bodies and also with surfaces,
if i copy and rotate them?
what is the reason for these how can i avoid these?
all errors in the corner of the sphere segments!?

if i have a fine mesh with high densitiy i get on my surfaces
grooves where i have the lines !? the grove i get,
if i create solid elements with bodies, but also with shell
and expanding them to solid elements !?
what is the reason for these how can i avoid these?

quarter sphere with solid elements created with bodies with 16 elements
erroneous mesh:

quarter sphere with solid elements created with bodies with 12 elements
perfect mesh:

quarter sphere with solid elements created with bodies with 96 elements
mesh with grooves at the lines: and erroneous mesh:

and here with surface and QU elements:
with creating quarter sphere no errors, copy and rotoate them for half sphere
with erroneous mesh: the erroneuous mesh with a minimum amount of elemens 88
with less, it works fine:

seta part s A001 A002 A003
copy all all rot x 90

seta part s A001 A002 A003
copy all all rot x 90
copy all all rot x 180

creating-sphere-solid.fbd

 # written by dichtstoff
 # buckling Spherical shell
 # sigma critical according Theory of elastic stability
 # thickness of sphere
 # valu sphth 1.93
 # radius of the sphere
 # valu sphrd 200
 # valu points of sphere 
 # valu sphv1 45.000000000
 # valu sphv2 35.264389683
 # amount of element with line division
 # valu lidiv 96
 ######################################
 valu sphth 1.93
 valu sphrd 200 
 valu sphth        * sphth  +0.5
 valu sphthoffseto + sphrd sphth 
 valu sphthoffseti - sphrd sphth
 valu sphrdn       * sphrd  -1.0
 valu lidiv 96

 PNT base     0.00000       0.00000       0.000000000
 PNT D001     0.00000       0.00000       sphthoffseto
 PNT D002     0.00000       0.00000       sphthoffseti
 PNT D003     0.00000       sphthoffseto  0.000000000
 PNT D004     0.00000       sphthoffseti  0.000000000
 PNT D005     sphthoffseto  0.00000       0.000000000
 PNT D006     sphthoffseti  0.00000       0.000000000
 PNT roty     0.00000       sphrdn        0.000000000
 PNT rotx      sphrdn       0.00000       0.000000000

 LINE ! D001 D002 1
 LINE ! D003 D004 1
 LINE ! D005 D006 1

 seta L001 l L001
 seta L002 l L002
 seta L003 l L003

 copy L001 L004 rot y 45
 copy L002 L005 rot x 45
 copy L003 L006 rot z 45

 seta roty p roty
 seta rotx p rotx

 move roty rot z  45
 move rotx rot z -45

 copy L006 L007 rot roty rotx -35.264389683

 LINE ! D001 D007 base lidiv
 LINE ! D001 D009 base lidiv
 LINE ! D002 D008 base lidiv
 LINE ! D002 D00A base lidiv
 LINE ! D003 D009 base lidiv
 LINE ! D003 D00B base lidiv
 LINE ! D004 D00A base lidiv
 LINE ! D004 D00C base lidiv
 LINE ! D005 D00B base lidiv
 LINE ! D005 D007 base lidiv
 LINE ! D006 D00C base lidiv
 LINE ! D006 D008 base lidiv
 LINE ! D007 D00D base lidiv
 LINE ! D008 D00E base lidiv
 LINE ! D009 D00D base lidiv
 LINE ! D00A D00E base lidiv
 LINE ! D00B D00D base lidiv
 LINE ! D00C D00E base lidiv

 SHPE SPHi sph base sphthoffseti
 SHPE SPHo sph base sphthoffseto

 GSUR A001 + BLEND - L006 + L00O + L007 - L00P
 GSUR A002 + BLEND - L005 + L00M + L007 - L00N 
 GSUR A003 + BLEND - L004 + L00K + L007 - L00L  
 GSUR A004 + BLEND - L003 + L00G + L006 - L00I 
 GSUR A005 + BLEND - L003 + L00H + L004 - L00J 
 GSUR A006 + BLEND - L002 + L00C + L005 - L00E 
 GSUR A007 + BLEND - L002 + L00D + L006 - L00F
 GSUR A008 + BLEND - L001 + L008 + L004 - L00A 
 GSUR A009 + BLEND - L001 + L009 + L005 - L00B  
 GSUR A010 + SPHo - L00G + L00H + L00K - L00O 
 GSUR A011 + SPHi - L00I + L00J + L00L - L00P
 GSUR A012 + SPHi - L00E + L00F + L00P - L00N 
 GSUR A013 + SPHo - L00C + L00D + L00O - L00M
 GSUR A014 + SPHi - L00A + L00B + L00N - L00L 
 GSUR A015 + SPHo - L008 + L009 + L00M - L00K
 GBOD B001 NORM + A001 - A005 - A003 - A011 + A004 + A010
 GBOD B002 NORM - A001 + A006 + A002 + A013 - A007 - A012 
 GBOD B003 NORM - A002 + A008 + A003 + A015 - A009 - A014

 #div all mult 2

 elty all HE8
 mesh all

 prnt se
 plot sa all
 plot e all

written by cgx_2.22

[dichtstoff]

I have checked the file with cgx 2.17, and there is no groove and no erroneous mesh !?
So it could be a bug or smiular?
quarter mesh with line division 100 and perfect mesh:

wbr dichtstoff

[dichtstoff]

calculix_2.XX_4win / bug

cgx_2.22

cgx_2.21

cgx_2.20

cgx_2.19 no bug

[dichtstoff]

so , one step more clear,
cgx_2.19_4win
is not working with SHPE,
but the others do:
cgx_2.20_4win
cgx_2.21_4win
cgx_2.22_4win
result with node 180:

cgx 2.19
SHPE
180,9.507439188655e-01,2.811210614975e-01,-1.186187337043e-01
BLEND
180,9.507439188655e-01,2.811210614975e-01,-1.186187337043e-01

cgx 2.20
SHPE
180,9.523752493016e-01,2.816033422185e-01,-1.188224700458e-01
BLEND
180,9.507439188655e-01,2.811210614975e-01,-1.186187337043e-01

cgx 2.21
SHPE
180,9.523752493016e-01,2.816033422185e-01,-1.188224700458e-01
BLEND
180,9.507439188655e-01,2.811210614975e-01,-1.186187337043e-01

cgx 2.22
SHPE
180,9.523752493016e-01,2.816033422185e-01,-1.188224700458e-01
BLEND
180,9.507439188655e-01,2.811210614975e-01,-1.186187337043e-01

[dichtstoff]

and the problem (BUG) with SHPE still there:

test-sphere.fbd

 # written by dichtstoff
 # buckling Spherical shell
 # sigma critical according Theory of elastic stability
 # thickness of sphere
 # valu sphth 1.93
 # radius of the sphere
 # valu sphrd 200
 # valu points of sphere 
 # valu sphv1 45.000000000
 # valu sphv2 35.264389683
 # amount of element with line division
 # valu lidiv 100
 ######################################
 valu sphth 1.93
 valu sphrd 200 
 valu sphth        * sphth  +0.5
 valu sphthoffseto + sphrd sphth 
 valu sphthoffseti - sphrd sphth
 valu sphrdn       * sphrd  -1.0
 valu lidiv 10

 PNT base     0.00000       0.00000       0.000000000
 PNT D001     0.00000       0.00000       sphthoffseto
 PNT D002     0.00000       0.00000       sphthoffseti
 PNT D003     0.00000       sphthoffseto  0.000000000
 PNT D004     0.00000       sphthoffseti  0.000000000
 PNT D005     sphthoffseto  0.00000       0.000000000
 PNT D006     sphthoffseti  0.00000       0.000000000
 PNT roty     0.00000       sphrdn        0.000000000
 PNT rotx      sphrdn       0.00000       0.000000000

 LINE ! D001 D002 1
 LINE ! D003 D004 1
 LINE ! D005 D006 1

 seta L001 l L001
 seta L002 l L002
 seta L003 l L003

 copy L001 L004 rot y 45
 copy L002 L005 rot x 45
 copy L003 L006 rot z 45

 seta roty p roty
 seta rotx p rotx

 move roty rot z  45
 move rotx rot z -45

 copy L006 L007 rot roty rotx -35.264389683

 LINE ! D001 D007 base 1
 LINE ! D001 D009 base 1
 LINE ! D002 D008 base 1
 LINE ! D002 D00A base 1
 LINE ! D003 D009 base 1
 LINE ! D003 D00B base 1
 LINE ! D004 D00A base 1
 LINE ! D004 D00C base 1
 LINE ! D005 D00B base 1
 LINE ! D005 D007 base 1
 LINE ! D006 D00C base 1
 LINE ! D006 D008 base 1
 LINE ! D007 D00D base 1
 LINE ! D008 D00E base 1
 LINE ! D009 D00D base 1
 LINE ! D00A D00E base 1
 LINE ! D00B D00D base 1
 LINE ! D00C D00E base 1

 seta linedivisonset l L008 L009 L00A L00B L00C L00D
 seta linedivisonset l L00E L00F L00G L00H L00I L00J
 seta linedivisonset l L00K L00L L00M L00N L00O L00P

 SHPE SPHi sph base sphthoffseti
 SHPE SPHo sph base sphthoffseto

 GSUR A001 + BLEND - L006 + L00O + L007 - L00P
 GSUR A002 + BLEND - L005 + L00M + L007 - L00N 
 GSUR A003 + BLEND - L004 + L00K + L007 - L00L  
 GSUR A004 + BLEND - L003 + L00G + L006 - L00I 
 GSUR A005 + BLEND - L003 + L00H + L004 - L00J 
 GSUR A006 + BLEND - L002 + L00C + L005 - L00E 
 GSUR A007 + BLEND - L002 + L00D + L006 - L00F
 GSUR A008 + BLEND - L001 + L008 + L004 - L00A 
 GSUR A009 + BLEND - L001 + L009 + L005 - L00B  
 GSUR A010 + SPHo - L00G + L00H + L00K - L00O 
 GSUR A011 + SPHi - L00I + L00J + L00L - L00P
 GSUR A012 + SPHi - L00E + L00F + L00P - L00N 
 GSUR A013 + SPHo - L00C + L00D + L00O - L00M
 GSUR A014 + SPHi - L00A + L00B + L00N - L00L 
 GSUR A015 + SPHo - L008 + L009 + L00M - L00K
 #GSUR A010 + BLEND - L00G + L00H + L00K - L00O 
 #GSUR A011 + BLEND - L00I + L00J + L00L - L00P
 #GSUR A012 + BLEND - L00E + L00F + L00P - L00N 
 #GSUR A013 + BLEND - L00C + L00D + L00O - L00M
 #GSUR A014 + BLEND - L00A + L00B + L00N - L00L 
 #GSUR A015 + BLEND - L008 + L009 + L00M - L00K
 GBOD B001 NORM - A005 - A001 - A003 - A011 + A004 + A010
 GBOD B002 NORM - A001 + A006 + A002 + A013 - A007 - A012 
 GBOD B003 NORM - A002 + A008 + A003 + A015 - A009 - A014

 div  linedivisonset mult lidiv

 div all mult 2

 elty all HE20
 mesh all

 seta nodes n all
 enq  nodes fix-x   rec 0 _ _ 1
 enq  nodes fix-y   rec _ 0 _ 1
 enq  nodes fix-z   rec _ _ 0 1
 enq  nodes monitor rec 0 0 _ 1 

 prnt se
 plot n fix-x r
 plus n fix-y b
 plus n fix-z g
 send all abq
 send fix-x abq nam
 send fix-y abq nam
 send fix-z abq nam
 send monitor abq
 plus e all

[dichtstoff]

i have seen maybe one more issue:
the geometry of the elements are maybe not correct created with cgx?:

node: value: x: y: z:
7192 3.621672e-02 -198.236000 -12.962500 12.962500
7193 3.592476e-02 -200.158000 -13.088200 13.088200
7196 3.607717e-02 -199.197000 -13.025400 13.025300

node: value: x: y: z:
7199 3.470693e-02 -196.970000 -25.892500 12.871400
7200 3.442284e-02 -198.880000 -26.143500 12.996200
7203 3.458847e-02 -197.869000 -25.735900 12.932500

if you look at the first node set it makes sense
7192 - 7196 -7193 with y coordination :
-12.962500 -13.025400 -13.088200

7199 - 7203 - 7200 with y coordination:
-25.892500 -25.735900 -26.143500
value should be in the range between node 7199 and 7200 !?

wbr

[klaus.wittig]

@dichtstoff

Is it possible to create surfaces with SPHE with mouse and interactive commands q ?
(what is the abbreviations q like qsurf,qshp, ect.?)
or only manual per script?

with qshp you can only create a plane (not all the other types) but you can relate any existing shape to an existing surface. So when you had created the surface (for example with qsur) you create a shape with SHPE and you use qshp to assign it to the surface (cumbersome, I know)

The mesher you used is the regular mesher. It requires the same mesh density on opposite edges, otherwise you get a transition (happens also if you use surfs with 3 or 5 edges since cgx internally uses 4 edges for the regular mesh). In your case this transition forces two edges of the deformed elements to be on the same edge of the surfaces as it seems by looking at the picture (which results in a bad shape). So only certain divisions work well as you also mentioned. You may try the unstructured mesher by using the element tr6u.

The only explanation for the grove which comes to my mind is a tiny geometrical difference between the lines and the shape. I have to look into that.

[klaus.wittig]

@dichtstoff

and the problem (BUG) with SHPE still there:

.. I could reproduce the bad elements with your fbd file. It seems to be a problem of the nurbs representation of a sphere. Close to the poles of a nurbs sphere the algorithm seems to have problems. Then this is the reason why very fine meshes fail at that location.

[dichtstoff]

thank you for taking time.

i have updated the file:

test-sphere.fbd

now you can reproduce the problem with these line:

valu lidiv 10 = no problem
valu lidiv 12 = erroneous mesh
valu lidiv 100 = ……

i have three bodies:
GBOD B001 NORM - A005 - A001 - A003 - A011 + A004 + A010
GBOD B002 NORM - A001 + A006 + A002 + A013 - A007 - A012
GBOD B003 NORM - A002 + A008 + A003 + A015 - A009 - A014
and i think the problem exists only with bodie B001!?
Update:
and if you move your geometry, the erroneous mesh looks different:

move all tra 100 100 100

elty all HE20
mesh all

[dichtstoff]

the issue is with the surfaces
if you mesh the two surface A010 and A011 you get an erroneous mesh:
all other surfaces works fine !?

[klaus.wittig]

At least I understand the problem and can show you a way how to avoid that problem until I improved the code. You just have to rotate the shapes a bit to avoid having the poles of the sphere in the region of your mesh. I added some lines to your code.

SHPE SPHi sph base sphthoffseti
SHPE SPHo sph base sphthoffseto

#Add this lines after you created the spheres:

seta SPH sh SPHi SPHo
comp SPH do
move SPH rot y 20

[klaus.wittig]

I will later care about the small radial mismatch between the nodes in the surfaces and the nodes on the lines.

[dichtstoff]

i have create a new half sphere.
I created all nodes, and then created manual all lines, surfaces
and all 12 bodies. now i have the problem with the points of
the sphere, because i get always a erroneous mesh.
so i moved the points in one body, meshed all other bodies,
except the one with the points, then i moved them out ,
meshed the last body and add them with merg n together
you can mesh with higher density than ccx can
handle:

 # written by dichtstoff
 # buckling Spherical shell
 # sigma critical according Theory of elastic stability
 # thickness of sphere
 # valu sphth 1.93
 # radius of the sphere
 # valu sphrd 200
 # valu points of sphere 
 # valu sphv1 45.000000000
 # valu sphv2 35.264389683
 # amount of element with line division
 # valu lidiv 72
 ######################################
 valu sphth 1.93
 valu sphrd 200 
 valu sphth        * sphth  +0.5
 valu sphthoffseto + sphrd sphth 
 valu sphthoffseti - sphrd sphth
 valu sphrdn       * sphrd  -1.0
 valu lidiv 100

 PNT base     0.00000       0.00000       0.000000000
 PNT D001     0.00000       0.00000       sphthoffseto
 PNT D002     0.00000       0.00000       sphthoffseti
 PNT D003     0.00000       sphthoffseto  0.000000000
 PNT D004     0.00000       sphthoffseti  0.000000000
 PNT D005     sphthoffseto  0.00000       0.000000000
 PNT D006     sphthoffseti  0.00000       0.000000000
 PNT roty     0.00000       500           0.000000000
 PNT rotx     500           0.00000       0.000000000

 LINE ! D001 D002 1
 LINE ! D003 D004 1
 LINE ! D005 D006 1

 seta L001 l L001
 seta L002 l L002
 seta L003 l L003

 copy L001 L004 rot y 45
 copy L002 L005 rot x 45
 copy L003 L006 rot z 45

 seta roty p roty
 seta rotx p rotx

 move roty rot z  45
 move rotx rot z -45

 copy L006 L007 rot roty rotx 35.264389683

 seta part l L003 L004 L006 L007
 copy part new rot z  -90
 copy part new rot z -180
 del se part new

 seta part l L006 L007 
 copy part part rot z 90
 del se part

 seta ldone l L001 L002 L003 L004 L005 L006
 seta ldone l L007 L008 L009 L00A L00B L00C
 seta ldone l L00D L00E L00F L00G L00H L00I

 LINE ! D001 D007 base 1
 LINE ! D001 D009 base 1
 LINE ! D001 D00H base 1
 LINE ! D001 D00P base 1
 LINE ! D002 D008 base 1
 LINE ! D002 D00I base 1
 LINE ! D002 D00Q base 1
 LINE ! D002 D00A base 1
 LINE ! D007 D005 base 1
 LINE ! D007 D00D base 1
 LINE ! D007 D00L base 1
 LINE ! D008 D00E base 1
 LINE ! D008 D00M base 1
 LINE ! D008 D006 base 1
 LINE ! D009 D00X base 1
 LINE ! D009 D003 base 1
 LINE ! D009 D00D base 1
 LINE ! D00A D00Y base 1
 LINE ! D00A D004 base 1
 LINE ! D00A D00E base 1
 LINE ! D00B D00D base 1
 LINE ! D00B D003 base 1
 LINE ! D00B D005 base 1
 LINE ! D00C D00E base 1
 LINE ! D00C D004 base 1
 LINE ! D00C D006 base 1
 LINE ! D00H D00L base 1
 LINE ! D00H D00T base 1
 LINE ! D00H D00F base 1
 LINE ! D00I D00M base 1
 LINE ! D00I D00U base 1
 LINE ! D00I D00G base 1
 LINE ! D00J D00L base 1
 LINE ! D00J D005 base 1
 LINE ! D00J D00F base 1
 LINE ! D00P D00T base 1
 LINE ! D00P D00N base 1
 LINE ! D00P D00X base 1
 LINE ! D00K D00M base 1
 LINE ! D00K D006 base 1
 LINE ! D00K D00G base 1
 LINE ! D00Q D00U base 1
 LINE ! D00Q D00O base 1
 LINE ! D00Q D00Y base 1
 LINE ! D00S D00U base 1
 LINE ! D00S D00G base 1
 LINE ! D00S D00O base 1
 LINE ! D00R D00T base 1
 LINE ! D00R D00F base 1
 LINE ! D00R D00N base 1
 LINE ! D00V D00X base 1
 LINE ! D00V D00N base 1
 LINE ! D00V D003 base 1
 LINE ! D00W D00Y base 1
 LINE ! D00W D00O base 1
 LINE ! D00W D004 base 1

 GSUR A001 + BLEND - L001 + L00I + L004 - L00M 
 GSUR A002 + BLEND - L001 + L00K + L009 - L00N 
 GSUR A003 + BLEND - L001 + L00L + L00D - L00O 
 GSUR A004 + BLEND - L001 + L00J + L005 - L00P
 GSUR A005 + BLEND - L00D + L01H + L00F - L01N 
 GSUR A006 + BLEND + L00E + L01Q - L00F - L01T 
 GSUR A007 + BLEND - L009 + L019 + L00F - L01C
 GSUR A008 + BLEND - L009 + L018 + L00B - L01B 
 GSUR A009 + BLEND - L00A + L01E + L00B - L01K 
 GSUR A010 + BLEND - L004 + L00S + L00B - L00U
 GSUR A011 + BLEND - L004 + L00R + L007 - L00T 
 GSUR A012 + BLEND - L005 + L00Y + L007 - L011 
 GSUR A013 + BLEND - L006 + L012 + L007 - L015
 GSUR A014 + BLEND - L005 + L00W + L00H - L00Z 
 GSUR A015 + BLEND - L00D + L01J + L00H - L01P 
 GSUR A016 + BLEND - L00G + L01W + L00H - L01Z
 GSUR A017 + BLEND + L008 - L01R - L00E + L01U 
 GSUR A018 + BLEND - L008 - L01A + L009 + L01D 
 GSUR A019 + BLEND - L008 - L01G + L00A + L01M 
 GSUR A020 + BLEND - L002 - L013 + L006 + L016 
 GSUR A021 + BLEND - L002 - L00X + L005 + L010 
 GSUR A022 + BLEND - L002 - L01Y + L00G + L021
 GSUR A023 + BLEND - L00C - L01X + L00G + L020 
 GSUR A024 + BLEND + L00C - L01S - L00E + L01V 
 GSUR A025 + BLEND - L00C - L01I + L00D + L01O 
 GSUR A026 + BLEND - L003 - L01F + L00A + L01L 
 GSUR A027 + BLEND - L003 - L00Q + L004 + L00V 
 GSUR A028 + BLEND - L003 - L014 + L006 + L017

 SHPE SPHi sph base sphthoffseti
 SHPE SPHo sph base sphthoffseto

 #### outside
 GSUR A101 + SPHo - L00I + L00K + L018 - L00S 
 GSUR A102 + SPHo - L00I + L00J + L00Y - L00R
 GSUR A103 + SPHo - L00K + L00L + L01H - L019 
 GSUR A104 + SPHo - L00J + L00L + L01J - L00W
 GSUR A105 + SPHo - L00Q + L00S - L01E + L01F 
 GSUR A106 + SPHo - L00Q + L00R - L012 + L014 
 GSUR A107 + SPHo - L00X + L00Y - L012 + L013 
 GSUR A108 + SPHo - L00W + L00X - L01Y + L01W 
 GSUR A109 + SPHo - L01I + L01J - L01W + L01X 
 GSUR A110 + SPHo - L01H + L01I - L01V + L01T 
 GSUR A111 + SPHo - L019 + L01A - L01U + L01T 
 GSUR A112 + SPHo - L018 + L01A - L01G + L01E

 seta outside A101 A102 A103 A104 A105 A106
 seta outside A107 A108 A109 A110 A111 A112

 #### inside
 GSUR A201 + SPHi - L00M + L00N + L01B - L00U
 GSUR A202 + SPHi - L00M + L00P + L011 - L00T 
 GSUR A203 + SPHi - L00N + L00O + L01N - L01C 
 GSUR A204 + SPHi - L00O + L00P + L00Z - L01P
 GSUR A205 + SPHi - L00U + L00V - L01L + L01K 
 GSUR A206 + SPHi - L00T + L00V - L017 + L015 
 GSUR A207 + SPHi - L010 + L011 - L015 + L016 
 GSUR A208 + SPHi - L00Z + L010 - L021 + L01Z 
 GSUR A209 + SPHi - L01O + L01P - L01Z + L020 
 GSUR A210 + SPHi - L01N + L01O - L01S + L01Q 
 GSUR A211 + SPHi - L01C + L01D - L01R + L01Q 
 GSUR A212 + SPHi - L01B + L01D - L01M + L01K 

 seta inside A201 A202 A203 A204 A205 A206
 seta inside A207 A208 A209 A210 A211 A212

 GBOD B001 NORM + A001 - A008 + A010 - A201 - A002 + A101 
 GBOD B002 NORM - A012 + A001 - A011 - A102 + A004 + A202
 GBOD B003 NORM + A002 - A005 + A007 - A203 - A003 + A103  
 GBOD B004 NORM + A003 - A014 + A015 - A204 - A004 - A104
 GBOD B005 NORM + A009 - A027 - A010 + A205 + A026 + A105 
 GBOD B006 NORM + A011 - A028 - A013 - A206 + A027 - A106 
 GBOD B007 NORM - A020 + A012 + A013 + A107 - A021 - A207 
 GBOD B008 NORM + A014 - A022 - A016 - A208 + A021 + A108 
 GBOD B009 NORM - A023 + A015 + A016 + A109 - A025 - A209 
 GBOD B010 NORM + A005 + A024 + A006 - A210 + A025 + A110  
 GBOD B011 NORM - A006 - A018 - A007 + A211 - A017 - A111
 GBOD B012 NORM + A008 - A019 - A009 - A212 + A018 + A112

 seta linedive l all
 setr linedive se ldone
 div  linedive mult lidiv

 seta SPH sh SPHi SPHo
 comp SPH do
 move SPH rot z 22.5

 seta correct b B001 B002 B003 B004 B005
 seta correct b B007 B008 B009 B010 B011 B012
 seta bug b B006

 div all mult 2

 elty all HE20R
 mesh correct

 move SPH rot z 45
 mesh bug

 merg n all

 seta nodes n all
 enq  nodes boundary   rec _ _ 0 1
 enq  nodes monitor    rec 0 0 _ 1

 prnt se
 plot n boundary r
 send all abq
 send boundary abq nam
 send monitor abq
 #plot n all
 #plus pa SPH b
 #plus ba all

[klaus.wittig]

Meanwhile I could improve the accuracy and solve the bad elements problem. It will be available with the next release. If somebody needs the fixes before we update the code in the internet he may contact me.