Esempio n. 1
0
 def percentThickness(self, percent):
     leadingedge = self.contour.getLeadingEdge()
     trailingedge = self.contour.getTrailingEdge()
     
     # initialise the surface matrices
     #surf = ones(self.resolution)*nan #*self.thick
     width = self.contour.getWidth()
     height = self.contour.getHeight()
     x_axis = linspace(0,width,self.resolution[0])
     y_axis = linspace(0,height,self.resolution[1])
     
     # interpolate the edges to fit the grid
     leading_edge = interpolate(leadingedge[:,1],leadingedge[:,0],y_axis,len(leadingedge[:,1]),self.resolution[1])
     trailing_edge = interpolate(trailingedge[:,1],trailingedge[:,0],y_axis,len(trailingedge[:,1]),self.resolution[1])
             
     chordlength = trailing_edge - leading_edge
     
     thick = chordlength * float(percent)
     
     self.set_basethickness(thick[0])
     
     return thick
Esempio n. 2
0
 def gen_surface(self, exponConstant):
     """ generate the 3D surface of the fin """
     
     leadingedge = self.contour.getLeadingEdge()
     trailingedge = self.contour.getTrailingEdge()
     
     # initialise the surface matrices
     #surf = ones(self.resolution)*nan #*self.thick
     width = self.contour.getWidth()
     height = self.contour.getHeight()
     x_axis = linspace(0,width,self.resolution[0])
     y_axis = linspace(0,height,self.resolution[1])
     
     # interpolate the edges to fit the grid
     leading_edge = interpolate(leadingedge[:,1],leadingedge[:,0],y_axis,len(leadingedge[:,1]),self.resolution[1])
     trailing_edge = interpolate(trailingedge[:,1],trailingedge[:,0],y_axis,len(trailingedge[:,1]),self.resolution[1])
     
     
     # rescale to grid
     dx = width / self.resolution[0]
     leading_edge_grid = (leading_edge / dx).round()
     trailing_edge_grid = (trailing_edge / dx).round()
     
     # generate the profiles and surface
     if exponConstant == 0:
         thick = self.thickness(height)
     elif exponConstant < 0:
         thick = self.percentThickness(-exponConstant)
     else:
         thick = self.exponentialThickness(exponConstant)
     
     surf = make_surface(self.resolution[0],leading_edge_grid,trailing_edge_grid,thick,self.resolution[1])
     
     #print surf
     self.surf = surf
     self.x_axis = x_axis
     self.y_axis = y_axis