""" Johnson Körperr J48 22.02.2021 www.3d-meier.de http://dmccooey.com/polyhedra/GyroelongatedPentagonalBirotunda.html """ import c4d import math # Variablen und Konstanten Titel = 'Johnson Körper J48' # Name NP = 40 # Anzahl Punkte N3 = 40 # Anzahl Dreiecke N5 = 12 # Anzahl Fuenfecke Radius = 2 # Radius der Eckpunktkugeln C0 = math.sqrt(10 * (5 - math.sqrt(5))) / 20 C1 = math.sqrt(10 * (5 + math.sqrt(5))) / 20 C2 = math.sqrt(2 * (math.sqrt(2 * (25 + 11 * math.sqrt(5))) - 2 * math.sqrt(5) - 4)) / 4 C3 = math.sqrt(5 * (5 + 2 * math.sqrt(5))) / 10 C4 = (1 + math.sqrt(5)) / 4 C5 = math.sqrt(10 * (5 + math.sqrt(5))) / 10 C6 = math.sqrt(2 * (5 + math.sqrt(5))) / 4 C7 = math.sqrt(10 * (25 + 11 * math.sqrt(5))) / 20 C8 = (2 * math.sqrt(10 * (5 + math.sqrt(5))) + 5 * math.sqrt(2 * (math.sqrt(2 * (25 + 11 * math.sqrt(5))) - 2 * math.sqrt(5) - 4))) / 20 C9 = (3 + math.sqrt(5)) / 4 C10 = math.sqrt(5 * (5 + 2 * math.sqrt(5))) / 5 C11 = math.sqrt(5 + 2 * math.sqrt(5)) / 2 C12 = (1 + math.sqrt(5)) / 2 C13 = (4 * math.sqrt(5 * (5 + 2 * math.sqrt(5))) + 5 * math.sqrt(2 * (math.sqrt(2 * (25 + 11 * math.sqrt(5))) - 2 * math.sqrt(5) - 4))) / 20 Punkte = [( 0.0, C13, C5), ( -C5, -C13, 0.0), ( C4, C13, C0), ( -C4, C13, C0), ( -C0, -C13, C4), ( -C0, -C13, -C4), ( 0.5, C13, -C3), (-0.5, C13, -C3), ( C3, -C13, 0.5), ( C3, -C13, -0.5), ( 0.0, C8, -C10), ( C10, -C8, 0.0), ( C9, C8, -C1), ( -C9, C8, -C1), ( C1, -C8, C9), ( C1, -C8, -C9), ( C4, C8, C7), ( -C4, C8, C7), ( -C7, -C8, C4), ( -C7, -C8, -C4), ( C12, C2, 0.0), (-C12, C2, 0.0), ( 0.0, -C2, C12), ( 0.0, -C2, -C12), ( 0.5, C2, C11), ( 0.5, C2, -C11), (-0.5, C2, C11), (-0.5, C2, -C11), ( C11, -C2, 0.5), ( C11, -C2, -0.5), (-C11, -C2, 0.5), (-C11, -C2, -0.5), ( C9, C2, C6), ( C9, C2, -C6), ( -C9, C2, C6), ( -C9, C2, -C6), ( C6, -C2, C9), ( C6, -C2, -C9), ( -C6, -C2, C9), ( -C6, -C2, -C9)] Dreiecke = [( 24, 16, 32, 32 ), ( 34, 17, 26, 26 ), ( 35, 13, 21, 21 ), ( 25, 10, 27, 27 ), ( 20, 12, 33, 33 ), ( 31, 19, 39, 39 ), ( 38, 18, 30, 30 ), ( 36, 14, 22, 22 ), ( 29, 11, 28, 28 ), ( 23, 15, 37, 37 ), ( 24, 32, 36, 36 ), ( 34, 26, 38, 38 ), ( 35, 21, 31, 31 ), ( 25, 27, 23, 23 ), ( 20, 33, 29, 29 ), ( 31, 39, 35, 35 ), ( 38, 30, 34, 34 ), ( 36, 22, 24, 24 ), ( 29, 28, 20, 20 ), ( 23, 37, 25, 25 ), ( 24, 22, 26, 26 ), ( 34, 30, 21, 21 ), ( 35, 39, 27, 27 ), ( 25, 37, 33, 33 ), ( 20, 28, 32, 32 ), ( 31, 21, 30, 30 ), ( 38, 26, 22, 22 ), ( 36, 32, 28, 28 ), ( 29, 33, 37, 37 ), ( 23, 27, 39, 39 ), ( 0, 2, 16, 16 ), ( 3, 0, 17, 17 ), ( 7, 3, 13, 13 ), ( 6, 7, 10, 10 ), ( 2, 6, 12, 12 ), ( 1, 5, 19, 19 ), ( 4, 1, 18, 18 ), ( 8, 4, 14, 14 ), ( 9, 8, 11, 11 ), ( 5, 9, 15, 15 )] Fuenfecke = [( 0, 3, 7, 6, 2 ), ( 1, 4, 8, 9, 5 ), ( 24, 26, 17, 0, 16 ), ( 34, 21, 13, 3, 17 ), ( 35, 27, 10, 7, 13 ), ( 25, 33, 12, 6, 10 ), ( 20, 32, 16, 2, 12 ), ( 31, 30, 18, 1, 19 ), ( 38, 22, 14, 4, 18 ), ( 36, 28, 11, 8, 14 ), ( 29, 37, 15, 9, 11 ), ( 23, 39, 19, 5, 15 )] Faktor = 100 # Skalierungsfaktor k = 1 # Korrekturfaktor fuer Kantenlaenge 1 #************************************************************************ def CreateNullobjekt1(): obj = c4d.BaseObject(c4d.Onull) obj.SetName(Titel) obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateNullobjekt2(): obj = c4d.BaseObject(c4d.Onull) obj.SetName('Ecken') obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateNullobjekt3(): obj = c4d.BaseObject(c4d.Onull) obj.SetName('Kanten') obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateNullobjekt4(): obj = c4d.BaseObject(c4d.Onull) obj.SetName('Polygone') obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateDreiecke(): obj = c4d.BaseObject(c4d.Opolygon) obj.ResizeObject(NP,N3) obj.SetName('Dreiecke') # Punkte uebergeben for i in xrange(NP): obj.SetPoint(i, c4d.Vector(Punkte[i][0]*Faktor*k, Punkte[i][1]*Faktor*k, Punkte[i][2]*Faktor*k)) # Dreiecke setzen for i in xrange(N3): obj.SetPolygon( i, c4d.CPolygon(Dreiecke[i][0], Dreiecke[i][1], Dreiecke[i][2], Dreiecke[i][3])) obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateFuenfecke(): obj = c4d.BaseObject(c4d.Opolygon) obj.ResizeObject(NP, N5*2) obj.SetName('Fünfecke') # Punkte uebergeben for i in xrange(NP): obj.SetPoint(i, c4d.Vector(Punkte[i][0]*Faktor*k, Punkte[i][1]*Faktor*k, Punkte[i][2]*Faktor*k)) zz = 0 # Zaehler fuer Polygone zuruecksetzen for i in xrange(N5): obj.SetPolygon(zz, c4d.CPolygon(Fuenfecke[i][0], Fuenfecke[i][1], Fuenfecke[i][2], Fuenfecke[i][4])) zz = zz + 1 obj.SetPolygon(zz, c4d.CPolygon(Fuenfecke[i][2], Fuenfecke[i][3], Fuenfecke[i][4], Fuenfecke[i][4])) zz = zz + 1 obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateSplineDreiecke(): obj = c4d.BaseObject(c4d.Ospline) obj.SetName("Spline-Dreiecke") obj.ResizeObject(N3*3) zz = 0 for i in xrange(N3): for j in xrange(3): x = Punkte[Dreiecke[i][j]][0] y = Punkte[Dreiecke[i][j]][1] z = Punkte[Dreiecke[i][j]][2] obj.SetPoint(zz, c4d.Vector(x*Faktor*k, y*Faktor*k, z*Faktor*k)) zz = zz + 1 # Segmente erzeugen obj.MakeVariableTag(c4d.Tsegment, N3) for i in range(0, N3): obj.SetSegment(i, 3, True) # Spline schliessen obj[c4d.SPLINEOBJECT_CLOSED] = True# obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateSplineFuenfecke(): obj = c4d.BaseObject(c4d.Ospline) obj.SetName("Spline-Fünfecke") obj.ResizeObject(N5*5) zz = 0 for i in xrange(N5): for j in xrange(5): x = Punkte[Fuenfecke[i][j]][0] y = Punkte[Fuenfecke[i][j]][1] z = Punkte[Fuenfecke[i][j]][2] obj.SetPoint(zz, c4d.Vector(x*Faktor*k, y*Faktor*k, z*Faktor*k)) zz = zz + 1 # Segmente erzeugen obj.MakeVariableTag(c4d.Tsegment, N5) for i in range(0, N5): obj.SetSegment(i, 5, True) # Spline schliessen obj[c4d.SPLINEOBJECT_CLOSED] = True obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def main(): nullobj1 = CreateNullobjekt1() nullobj2 = CreateNullobjekt2() # Ecken nullobj3 = CreateNullobjekt3() # Kanten nullobj4 = CreateNullobjekt4() # Polygone plyobj1 = CreateDreiecke() plyobj2 = CreateFuenfecke() splobj1 = CreateSplineDreiecke() splobj2 = CreateSplineFuenfecke() doc.InsertObject(nullobj1, None, None, True) doc.InsertObject(nullobj4, nullobj1, None, True) doc.InsertObject(nullobj3, nullobj1, None, True) doc.InsertObject(nullobj2, nullobj1, None, True) doc.InsertObject(splobj2, nullobj3, None, True) doc.InsertObject(splobj1, nullobj3, None, True) doc.InsertObject(plyobj2, nullobj4, None, True) doc.InsertObject(plyobj1, nullobj4, None, True) # Kugeln auf Eckpunkte setzen for i in range(0, NP): obj = c4d.BaseObject(c4d.Osphere) obj[c4d.PRIM_SPHERE_RAD] = Radius obj.SetName(str(i)) x = Punkte[i][0] y = Punkte[i][1] z = Punkte[i][2] obj.SetAbsPos(c4d.Vector(x*Faktor*k, y*Faktor*k, z*Faktor*k)) doc.InsertObject(obj, nullobj2, None, True) c4d.EventAdd() if __name__=='__main__': main()