Merge pull request #286 from tpaviot/review/2d_fillet_example

Added 2d fillet example

Author: tpaviot

examples/core_2d_fillet.py

@@ -0,0 +1,45 @@
+#!/usr/bin/env python
+##
+##This file is part of pythonOCC.
+##
+##pythonOCC is free software: you can redistribute it and/or modify
+##it under the terms of the GNU Lesser General Public License as published by
+##the Free Software Foundation, either version 3 of the License, or
+##(at your option) any later version.
+##
+##pythonOCC is distributed in the hope that it will be useful,
+##but WITHOUT ANY WARRANTY; without even the implied warranty of
+##MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+##GNU Lesser General Public License for more details.
+##
+##You should have received a copy of the GNU Lesser General Public License
+##along with pythonOCC.  If not, see <http://www.gnu.org/licenses/>.
+
+from OCC.gp import gp_Pnt,gp_Dir,gp_Pln
+from OCC.Display.SimpleGui import init_display
+from OCC.ChFi2d import ChFi2d_AnaFilletAlgo
+from OCC.BRepBuilderAPI import BRepBuilderAPI_MakeEdge,BRepBuilderAPI_MakeWire
+
+from core_geometry_utils import make_wire
+display,start_display, add_menu,add_functionto_menu = init_display()
+
+# Defining the points
+p1 = gp_Pnt(0, 0, 0) 
+p2 = gp_Pnt(5, 5, 0)
+p3 = gp_Pnt(-5,5, 0)
+
+# Making the edges
+ed1 = BRepBuilderAPI_MakeEdge(p3,p2).Edge()
+ed2 = BRepBuilderAPI_MakeEdge(p2,p1).Edge()
+
+#Making the 2dFillet
+f = ChFi2d_AnaFilletAlgo()
+f.Init(ed1,ed2,gp_Pln())
+radius = 1.0
+f.Perform(radius)
+fillet2d = f.Result(ed1,ed2)
+
+# Create and display a wire
+w = make_wire([ed1, fillet2d, ed2])
+display.DisplayShape(w)
+start_display()

examples/core_geometry_airfoil.py

@@ -48,9 +48,10 @@ def __init__(self, chord, span, profile):
 
     def make_shape(self):
         # 1 - retrieve the data from the UIUC airfoil data page
-        foil_dat_url = 'http://www.ae.illinois.edu/m-selig/ads/coord_seligFmt/%s.dat' % self.profile
+        foil_dat_url = 'http://m-selig.ae.illinois.edu/ads/coord_seligFmt/%s.dat' % self.profile
+        print("Connecting to m-selig, retrieving foil data")
         f = urllib2.urlopen(foil_dat_url)
-
+        print("Building foil geometry")
         plan = gp_Pln(gp_Pnt(0., 0., 0.), gp_Dir(0., 0., 1.))  # Z=0 plan / XY plan
         section_pts_2d = []
 
@@ -104,8 +105,8 @@ def _Tcol_dim_1(li, _type):
     return pts
 
 if __name__ == '__main__':
+    airfoil = UiucAirfoil(50., 200., 'b737a')
     from OCC.Display.SimpleGui import init_display
     display, start_display, add_menu, add_function_to_menu = init_display()
-    airfoil = UiucAirfoil(50., 200., 'b737a')
     display.DisplayShape(airfoil.shape, update=True)
     start_display()