Changes between Version 4 and Version 5 of salome

Timestamp:

Nov 14, 2017, 3:52:07 PM (7 years ago)

Author:

brankm

Comment:

Dodan tutorial za izdelavo panela.

salome

@@ -119 +119 @@
 salome.sg.updateObjBrowser(True)
 }}}
-[[BR]][[BR]][[BR]][[BR]]
+[[BR]][[BR]][[BR]][[BR]][[BR]][[BR]]
 '''Korak 4''': Izdelava povrsine. 
 
@@ -135 +135 @@
 
 Drevesno strukturo posodobimo, ce zelimo prikazati nove geometrijske elemente v modelnem oknu.
-[[BR]][[BR]][[BR]][[BR]][[BR]][[BR]]
+[[BR]][[BR]][[BR]][[BR]][[BR]][[BR]][[BR]][[BR]][[BR]]
 '''Korak 5''': Izvlek povrsine
 
@@ -151 +151 @@
 
 Celoten postopek izdelave izvleka je v priponki (''salomePrismExample.py'').
+
+=== Primer izdelave panela v fuzijskem reaktorju ITER
+
+V tem primeru bo prikazana izdelava krivulj v 3D prostoru in izvlek krivulj, s katerim generiramo povrsino v prostoru. Koncni cilj je izdelava panela, prikazanega na sliki. V datoteki data.txt imamo podane koordinate tock (x, y in z), iz katerih bomo generirali krivujo. To krivuljo bomo nato izvlekli v prostor s pomocjo daljic, ki so prav tako definirana s tockami (xe, ye in ze). Podatki so podani v formatu 
+[[Image(prismSalomePrismExample.png, 300, right)]]
+{{{
+x1 y1 z1 xe1 ye1 ze1
+x2 y2 z2 xe2 ye2 ze2
+x3 y3 z3
+...
+xn yn zn
+}}}
+
+1. Uvoz modulov
+
+Najprej je potrebno v python uvoziti vse knjiznice, ki jih bomo uporabili v Salome okolju. Nato pa je potrebno definirati novo studijo. 
+
+{{{
+#!python
+# Uvoz knjiznic
+import salome
+import numpy as np
+import GEOM
+from salome.geom import geomBuilder
+geompy = geomBuilder.New(salome.myStudy)
+}}}
+
+2. Uvoz podatkov v Python
+
+Podatke v stolpcih v tekstovnem formatu lahko preberemo na vec nacinov, v tem primeru smo uporabili funkcijo `loadtxt()` iz knjiznice `numpy`. Podatke za krivuljo shranimo v spremenljivko tipa `np.array()` z imenoma ''coord_x'' in ''coord_y''. Nato te tocke zrcalimo preko ravnine y-z, da dobimo krivuljo se za levo stran panela, in jih shranimo v spremenljivki `coord_x_mirror` in `coord_y_mirror`. 
+
+{{{
+#!python
+# Load data
+data = np.loadtxt('/pathtofile/data.txt')
+# convert coordinates to mm
+coord_x = data[:,0]*(1000)
+coord_y = data[:,1]*(1000)
+# mirror x coordinates over yz - plane
+coord_x_mirror = coord_x[::-1]*(-1)
+coord_y_mirror = coord_y[::-1]
+}}}
+
+Nato generiramo B-zlepek krivuljo s pomocjo metode `geompy.MakeInterpol()`, ki sprejme seznam tock, na katerih izvede interpolacijo. rezultat je polinom, ki predstavlja podano krivuljo. Nato izracunano krivuljo dodamo v studijo.
+{{{
+#!python
+# Create B-Spline curve on the set of points with
+# geompy.MakeInterpol() and add curve to Salome study.
+b_spline = geompy.MakeInterpol(points_x)
+b_splineID = geompy.addToStudy(b_spline,"b_spline_ft")
+b_spline_mirror = geompy.MakeInterpol(points_y)
+b_splineID_mirror = geompy.addToStudy(b_spline_mirror,"b_spline_ft_mirrored")
+}}}
+
+Sedaj je potrebno kreirati se daljico, s katero bomo povezali obe krivulji (na sliki oznacena z rumeno barvo).
+
+{{{
+#!python
+# Create line between both curves.
+ml_point1 = geompy.MakeVertex(-coord_x[0],coord_y[0],0)
+ml_point2 = geompy.MakeVertex(coord_x[0],coord_y[0],0)
+middle_line = geompy.MakeLineTwoPnt(ml_point1,ml_point2)
+middle_lineID = geompy.addToStudy(middle_line,"middle_line")
+}}}
+
+Nato je potrebno vse tri komponente zdruziti v skupen objekt, ki ga bomo izvlekli v prostor. To storimo z metodo `geompy.MakeWire()`, ki sprejme seznam krivulj, robov in daljic, ki jih zelimo zdruziti v en objekt.
+
+{{{
+#!python
+# Create wire out of all three components
+wire = geompy.MakeWire([b_spline_mirror,middle_line,b_spline])
+wireID = geompy.addToStudy(wire,"wire_bottom")
+}}}
+
+Nato izdelamo se krivuljo za izvlek, ki je prikazana na sliki. Ponovno uvozimo tocke in jih shranimo v spremenljivke tipa `np.array()`.
+
+{{{
+#!python
+# Load data for extrusion
+dataExtrusion = np.loadtxt('/pathtofile/dataExtrusion.txt')
+# Read data for extrusion
+coord_xe = dataExtrusion[:,0]
+coord_ye = dataExtrusion[:,1]
+coord_ze = dataExtrusion[:,2]
+}}}
+
+Koordinate ponovno pretvorimo v tocke tipa `geompy` z ukazom `geompy.MakeVertex()` ter nato generiramo polyline.
+
+{{{
+#!python
+# Convert data to geompy point and store it in a list.
+points_e = []
+for point_e in range(len(coord_xe)):
+    point_e = geompy.MakeVertex(coord_xe[point_e],coord_ye[point_e],coord_ze[point_e])
+    points_e.append(point_e)
+# Create polyline for extrusion and add it to study.
+polylineVert = geompy.MakePolyline(points_e)
+polylineVertID = geompy.addToStudy(polylineVert,"b_spline_ft_vert")
+}}}
+
+Na koncu ustvarimo izvlek iz obeh danih krivulj. Izvlek dodamo v studijo in posodobimo drevesno strukturo.
+
+{{{
+#!python
+# Create panel shape
+panel_FW = geompy.MakePipe(polylineVert, wire)
+geompy.addToStudy(panel_FW,'Panel')
+salome.sg.updateObjBrowser(True)
+}}}