def expand(self,iface,iedge,nedges,sign=1,angle=np.pi): """ expand face iface along edge iedge Parameters ---------- iface : int iedge : int nedges : int angle : float """ # add a new node to the graph self.nnode = self.nnode + 1 # get the selected edge (iedge) of the selected face (iface) ed = self.lfaces[iface].subshapes('Edge')[iedge] # construct the equivalent complex number of the 2 extremities of the # edge points = ed.poly() z0 = points[0][0] + 1j*points[0][1] z1 = points[1][0] + 1j*points[1][1] # construct the new face with nedges sides lz = [z0, z1] for k in range(nedges-2): dz = lz[-2]-lz[-1] lz.append(lz[-1]+ np.abs(dz)*np.exp(1j*(np.angle(dz)+sign*(nedges-2)*np.pi/nedges))) lpts = [] for k in range(nedges): lpts.append((np.real(lz[k]),np.imag(lz[k]),0.0)) lpts.append(lpts[0]) new_face = cm.plane(cm.polygon(lpts)) # control homogeneity of face normals if new_face.normal()[2]==-1: new_face = cm.plane(cm.polygon(lpts[::-1])) # append new face to list of faces self.lfaces.append(new_face) # new face number in the graph node_num = self.nnode - 1 self.add_node(node_num,normal=new_face.normal()) self.pos[node_num] = new_face.center()[0:2] self.add_edge(iface, node_num, angle=angle, iedge=iedge) # construct Shell from Face self.shell = cm.Shell(self.lfaces)
def test_polygon(self): w1 = cm.polygon([(0.0, 0.0, 0.0), (1.0, 0.0, 0.0), (1.5, 1.0, 0.0), (0.5, 1.5, 0.0), (-0.5, -0.5, 0.0)]) self.assert_(close(5.472, w1.length(), 0.001))
def test_polygon(self): w1 = cm.polygon([(0.0, 0.0, 0.0), (1.0, 0.0, 0.0), (1.5, 1.0, 0.0), (0.5, 1.5, 0.0), (-0.5, -0.5, 0.0)]) self.assert_(close(5.472, w1.length(), 0.001))
def wire_polygon(): w1 = cm.polygon([(0.0, 0.0, 0.0), (1.0, 0.0, 0.0), (1.5, 1.0, 0.0), (0.5, 1.5, 0.0), (-0.5, -0.5, 0.0)]) save_top(w1, 'wire_polygon.png', (0.0, 0.0, 0.0))
def __init__(self, **kwargs): """ Parameters ---------- pt : np.array (N x 2) points in the 0xy plane N : number of nodes of the first polygon l : edge length of the first polygon """ self.folded = False pt = kwargs.pop('pt', []) N = kwargs.pop('N', 3) self.l = kwargs.pop('l', 1) if pt == []: al = (N - 2) * np.pi / (2 * N) self.r = np.sin(al) / np.sin(2 * np.pi / N) self.s = np.sqrt(self.r**2 - self.l**2/4.) t = np.linspace(0, N, N+1) u = 2 * np.pi * t / N pt = self.r * np.c_[np.cos(u), np.sin(u)] N = pt.shape[0] # conversion to 3D self.pt = np.c_[pt[:, 0], pt[:, 1], np.zeros(N)] w0 = cm.polygon(self.pt) face0 = cm.plane(w0) self.lfaces = [face0] self.nnode = 1 nx.Graph.__init__(self) self.add_node(0,normal=face0.normal()) self.add_node(0) self.pos = dict() self.pos[0] = face0.center()[0:2]
def derived_revol(): e1 = cm.arc(1.0, -math.pi / 2, math.pi / 2) e1.translate((3.0, 0.0, 0.0)) w1 = cm.polygon([(3.0, 1.0, 0.0), (2.0, 1.0, 0.0), (2.0, -1.0, 0.0), (3.0, -1.0, 0.0)]) f1 = cm.plane(cm.wire([e1, w1])) f1.rotatex(math.pi / 2) s1 = cm.revol(f1, (0.0, 0.0, 0.0), (0.0, 0.0, 1.0), 2 * math.pi) save_iso(s1, 'derived_revol.png')
def logging_face_fillet(): w1 = cm.polygon([(-1.0, -1.0, 0.0), (1.0, -1.0, 0.0), (1.0, 1.0, 0.0), (-1.0, 1.0, 0.0), (-1.0, -1.0, 0.0)]) f1 = cm.plane(w1) f1.fillet(0.25, f1.nearest('vertex', [(1.0, 1.0, 0.0), (-1.0, -1.0, 0.0)])) save_top(f1, 'logging_face_fillet.png')
def derived_revol(): e1 = cm.arc(1.0, -math.pi / 2, math.pi / 2) e1.translate((3.0, 0.0, 0.0)) w1 = cm.polygon([(3.0, 1.0, 0.0), (2.0, 1.0, 0.0), (2.0, -1.0, 0.0), (3.0, -1.0, 0.0)]) f1 = cm.plane(cm.wire([e1, w1])) f1.rotatex(math.pi / 2) s1 = cm.revol(f1, (0.0, 0.0, 0.0), (0.0, 0.0, 1.0), 2 * math.pi) save_iso(s1, 'derived_revol.png')
def logging_face_fillet(): w1 = cm.polygon([(-1.0, -1.0, 0.0), (1.0, -1.0, 0.0), (1.0, 1.0, 0.0), (-1.0, 1.0, 0.0), (-1.0, -1.0, 0.0)]) f1 = cm.plane(w1) f1.fillet(0.25, f1.nearest('vertex', [(1.0, 1.0, 0.0), (-1.0, -1.0, 0.0)])) save_top(f1, 'logging_face_fillet.png')
def wire_polygon(): w1 = cm.polygon([(0.0, 0.0, 0.0), (1.0, 0.0, 0.0), (1.5, 1.0, 0.0), (0.5, 1.5, 0.0), (-0.5, -0.5, 0.0)]) save_top(w1, 'wire_polygon.png', (0.0, 0.0, 0.0))