ForkMaster/netgen
1
0
Fork 0
mirror of https://github.com/NGSolve/netgen.git synced 2024-11-11 16:49:16 +05:00
Code Issues Packages Projects Releases Wiki Activity

move py-tutorials

Browse Source
This commit is contained in:
Joachim Schoeberl 2014-10-01 13:52:08 +00:00
parent becaacaf73
commit b950bd6ca6
3 changed files with 0 additions and 112 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
python/exportNeutral.py
View File

@ -1,26 +0,0 @@
import sys
def Export (mesh, filename):
""" export Netgen mesh to neutral format """
print ("export mesh in neutral format to file = ", filename)
f = open (filename, 'w')
points = mesh.Points()
print (len(points), file=f)
for p in points:
print (p.p[0], p.p[1], p.p[2], file=f)
volels = mesh.Elements3D();
print (len(volels), file=f)
for el in volels:
print (el.index, end=" ", file=f)
for j in el.vertices:
print (j.nr, end=" ", file=f)
print(file=f)

21
python/mesh.py
View File

@ -1,21 +0,0 @@
import sys
sys.path.append("/opt/netgen/lib")
from libmesh.meshing import *
from libcsg.csg import *
geo = CSGeometry("shaft.geo")
param = MeshingParameters()
param.maxh = 10
print (param)
m1 = GenerateMesh (geo, param)
import exportNeutral
exportNeutral.Export (m1, "shaft.mesh")

65
python/shaft.geo
View File

@ -1,65 +0,0 @@
#
## Crankshaft
#
algebraic3d
solid p1 = plane (0, 0, 0; -1, 0, 0)
and plane (10, 0, 0; 1, 0, 0)
and
( plane (35, 0, 28; 0, -1, 3)
and plane (35, 0, -28; 0, -1, -3)
and plane (35, 0, 0; 0, 1, 0)
and plane (35, -30, 0; 0, -1, 0)
or cylinder (-10, 0, 0; 20, 0, 0; 30)
or cylinder (-10, -30, 0; 20, -30, 0; 20)
);
solid p2 = plane (35, 0, 0; -1, 0, 0)
and plane (45, 0, 0; 1, 0, 0)
and
( plane (35, 0, 28; 0, -1, 3)
and plane (35, 0, -28; 0, -1, -3)
and plane (35, 0, 0; 0, 1, 0)
and plane (35, -30, 0; 0, -1, 0)
or cylinder (30, 0, 0; 50, 0, 0; 30)
or cylinder (30, -30, 0; 50, -30, 0; 20)
);
solid p3 = plane (80, 0, 0; -1, 0, 0)
and plane (90, 0, 0; 1, 0, 0)
and
( plane (0, 0, 28; 0, 1, 3)
and plane (0, 0, -28; 0, 1, -3)
and plane (0, 0, 0; 0, -1, 0)
and plane (0, 30, 0; 0, 1, 0)
or cylinder (70, 0, 0; 100, 0, 0; 30)
or cylinder (70, 30, 0; 100, 30, 0; 20)
);
solid p4 = plane (115, 0, 0; -1, 0, 0)
and plane (125, 0, 0; 1, 0, 0)
and
( plane (35, 0, 28; 0, 1, 3)
and plane (35, 0, -28; 0, 1, -3)
and plane (35, 0, 0; 0, -1, 0)
and plane (35, 30, 0; 0, 1, 0)
or cylinder (110, 0, 0; 130, 0, 0; 30)
or cylinder (110, 30, 0;130, 30, 0; 20)
);
solid sh1 = cylinder (-50, 0, 0; 10, 0, 0; 15)
and plane (-40, 0, 0; -1, 0, 0)
and plane (5, 0, 0; 1, 0, 0);
solid sh2 = cylinder (30, 0, 0; 90, 0, 0; 15)
and plane (40, 0, 0; -1, 0, 0)
and plane (85, 0, 0; 1, 0, 0);
solid sh3 = cylinder (110, 0, 0; 170, 0, 0; 15)
and plane (120, 0, 0; -1, 0, 0)
and plane (165, 0, 0; 1, 0, 0);
solid pl1 = cylinder (0, -30, 0; 50, -30, 0; 10)
and plane (5, 0, 0; -1, 0, 0)
and plane (40, 0, 0; 1, 0, 0);
solid pl2 = cylinder (80, 30, 0; 130, 30, 0; 10)
and plane (85, 0, 0; -1, 0, 0)
and plane (120, 0, 0; 1, 0, 0);
#
#
solid main = p1 or p2 or p3 or p4 or sh1 or sh2 or sh3 or pl1 or pl2;
tlo main;