Elliptic cylinder with nested objects

1 week 3 days ago #2919 by D Mikro
Hello NG team,

I would like to ask you for advice on creating elliptical objects. I want to create an elliptical cylinder with nested holes around its perimeter. I use it for simulation in Matlab, where I also generated an attached geo file. When I started the NG mesh generator (Netgen version 5.3), I got stuck in the Start Findpoints operation, but I don't know why? Can you describe me what I'm doing wrong?

Note: The coordinates of the second and sixth object are the same, but with opposite signs with respect to the origin of the system coordinates. Second is Ok, sixth failed.
Here is the generated geo file.
#Elliptic cylinder 
#failed: 6th and 8th hole in the eliptic body 

algebraic3d

solid body0001 = ((((((((ellipticcylinder(0, 0, 0 ; 0.0957, 0, 0 ; 0, 0.0953, 0) 
and plane(0, 0, 0 ; -0, -0, -0.318) and plane(0, 0, 0.318 ; 0, 0, 0.318))))) 

and (not (cylinder(-5.51091e-19, 0.0863, 0.296 ; 4.59243e-19, 0.1028, 0.296 ; 0.0045) 
and plane(-5.51091e-19, 0.0863, 0.296 ; -1.01033e-18, -0.0165, -0) 
and plane(4.59243e-19, 0.1028, 0.296 ; 1.01033e-18, 0.0165, 0)) 

and not (cylinder(0.0612874, 0.0610419, 0.296 ; 0.0729241, 0.0727397, 0.296 ; 0.0045) 
and plane(0.0612874, 0.0610419, 0.296 ; -0.0116367, -0.0116978, -0) 
and plane(0.0729241, 0.0727397, 0.296 ; 0.0116367, 0.0116978, 0)) 

and not (cylinder(0.0867, 5.28204e-18, 0.296 ; 0.1032, 6.29661e-18, 0.296 ; 0.0045) 
and plane(0.0867, 5.28204e-18, 0.296 ; -0.0165, -1.01457e-18, -0) 
and plane(0.1032, 6.29661e-18, 0.296 ; 0.0165, 1.01457e-18, 0)) 

and not (cylinder(0.0612874, -0.061042, 0.296 ; 0.0729241, -0.0727397, 0.296 ; 0.0045) 
and plane(0.0612874, -0.061042, 0.296 ; -0.0116367, 0.0116978, -0) 
and plane(0.0729241, -0.0727397, 0.296 ; 0.0116367, -0.0116978, 0)) 

and not (cylinder(1.11688e-17, -0.0863, 0.296 ; 1.21791e-17, -0.1028, 0.296 ; 0.0045) 
and plane(1.11688e-17, -0.0863, 0.296 ; -1.01033e-18, 0.0165, -0) 
and plane(1.21791e-17, -0.1028, 0.296 ; 1.01033e-18, -0.0165, 0)) 
#6th hole 
and not (cylinder(-0.0612874, -0.0610419, 0.296 ; -0.0729241, -0.0727397, 0.296 ; 0.0045) 
and plane(-0.0612874, -0.0610419, 0.296 ; 0.0116367, 0.0116978, 0) 
and plane(-0.0729241, -0.0727397, 0.296 ; -0.0116367, -0.0116978, 0)) 

and not (cylinder(-0.0867, -1.64041e-17, 0.296 ; -0.1032, -1.84248e-17, 0.296 ; 0.0045) 
and plane(-0.0867, -1.64041e-17, 0.296 ; 0.0165, 2.02067e-18, -0) 
and plane(-0.1032, -1.84248e-17, 0.296 ; -0.0165, -2.02067e-18, 0)) 
#8th hole
and not (cylinder(-0.0612874, 0.061042, 0.296 ; -0.0729241, 0.0727397, 0.296 ; 0.0045) 
and plane(-0.0612874, 0.061042, 0.296 ; 0.0116367, -0.0116978, -0) 
and plane(-0.0729241, 0.0727397, 0.296 ; -0.0116367, 0.0116978, 0))
)))));

solid MainNetgenObject = body0001 ;
tlo MainNetgenObject -col=[0,1,0] -material=MainNetgenObject -maxh=0.02;

Please Log in or Create an account to join the conversation.

© 2019 Netgen/NGSolve