""" Johnson Körperr J76 11.01.2021 www.3d-meier.de http://dmccooey.com/polyhedra/DiminishedRhombicosidodecahedron.html """ import c4d import math # Variablen und Konstanten Titel = 'Johnson Körper J76' # Name NP = 55 # Anzahl Punkte N3 = 15 # Anzahl Dreiecke N4 = 25 # Anzahl Vierecke N5 = 11 # Anzahl Fuenfecke N10 = 1 # Anzahl Zehnecke Radius = 2 # Radius der Eckpunktkugeln C0 = (1 + math.sqrt(5)) / 4 C1 = (3 + math.sqrt(5)) / 4 C2 = (1 + math.sqrt(5)) / 2 C3 = (5 + math.sqrt(5)) / 4 C4 = (2 + math.sqrt(5)) / 2 Punkte = [( 0.5, 0.5, C4), ( 0.5, 0.5, -C4), ( 0.5, -0.5, C4), ( 0.5, -0.5, -C4), (-0.5, 0.5, C4), (-0.5, 0.5, -C4), (-0.5, -0.5, C4), (-0.5, -0.5, -C4), ( C4, 0.5, 0.5), ( C4, 0.5, -0.5), ( C4, -0.5, 0.5), ( C4, -0.5, -0.5), ( -C4, 0.5, 0.5), ( -C4, 0.5, -0.5), ( -C4, -0.5, 0.5), ( -C4, -0.5, -0.5), ( 0.5, C4, -0.5), ( 0.5, -C4, 0.5), ( 0.5, -C4, -0.5), (-0.5, C4, -0.5), (-0.5, -C4, 0.5), (-0.5, -C4, -0.5), ( 0.0, C1, -C3), ( 0.0, -C1, C3), ( 0.0, -C1, -C3), ( C3, 0.0, C1), ( C3, 0.0, -C1), ( -C3, 0.0, C1), ( -C3, 0.0, -C1), ( C1, C3, 0.0), ( C1, -C3, 0.0), ( -C1, C3, 0.0), ( -C1, -C3, 0.0), ( C1, C0, C2), ( C1, C0, -C2), ( C1, -C0, C2), ( C1, -C0, -C2), ( -C1, C0, C2), ( -C1, C0, -C2), ( -C1, -C0, C2), ( -C1, -C0, -C2), ( C2, C1, C0), ( C2, C1, -C0), ( C2, -C1, C0), ( C2, -C1, -C0), ( -C2, C1, C0), ( -C2, C1, -C0), ( -C2, -C1, C0), ( -C2, -C1, -C0), ( C0, C2, -C1), ( C0, -C2, C1), ( C0, -C2, -C1), ( -C0, C2, -C1), ( -C0, -C2, C1), ( -C0, -C2, -C1)] Zehnecke = [( 0, 33, 41, 29, 16, 19, 31, 45, 37, 4 )] Fuenfecke = [( 22, 52, 19, 16, 49 ), ( 23, 53, 20, 17, 50 ), ( 24, 51, 18, 21, 54 ), ( 25, 33, 0, 2, 35 ), ( 26, 36, 3, 1, 34 ), ( 27, 39, 6, 4, 37 ), ( 28, 38, 5, 7, 40 ), ( 29, 41, 8, 9, 42 ), ( 30, 44, 11, 10, 43 ), ( 31, 46, 13, 12, 45 ), ( 32, 47, 14, 15, 48 )] Vierecke = [( 1, 22, 49, 34 ), ( 2, 23, 50, 35 ), ( 3, 36, 51, 24 ), ( 5, 38, 52, 22 ), ( 6, 39, 53, 23 ), ( 7, 24, 54, 40 ), ( 8, 41, 33, 25 ), ( 9, 26, 34, 42 ), ( 10, 25, 35, 43 ), ( 11, 44, 36, 26 ), ( 12, 27, 37, 45 ), ( 13, 46, 38, 28 ), ( 14, 47, 39, 27 ), ( 15, 28, 40, 48 ), ( 16, 29, 42, 49 ), ( 17, 30, 43, 50 ), ( 18, 51, 44, 30 ), ( 19, 52, 46, 31 ), ( 20, 53, 47, 32 ), ( 21, 32, 48, 54 ), ( 0, 4, 6, 2 ), ( 1, 3, 7, 5 ), ( 8, 10, 11, 9 ), ( 12, 13, 15, 14 ), ( 17, 20, 21, 18 )] Dreiecke = [( 22, 1, 5, 5 ), ( 23, 2, 6, 6 ), ( 24, 7, 3, 3 ), ( 25, 10, 8, 8 ), ( 26, 9, 11, 11 ), ( 27, 12, 14, 14 ), ( 28, 15, 13, 13 ), ( 30, 17, 18, 18 ), ( 32, 21, 20, 20 ), ( 34, 49, 42, 42 ), ( 35, 50, 43, 43 ), ( 36, 44, 51, 51 ), ( 38, 46, 52, 52 ), ( 39, 47, 53, 53 ), ( 40, 54, 48, 48 )] 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 CreateVierecke(): obj = c4d.BaseObject(c4d.Opolygon) obj.ResizeObject(NP, N4) obj.SetName('Vierecke') # 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)) # Vierecke setzen for i in xrange(N4): obj.SetPolygon(i, c4d.CPolygon(Vierecke[i][0], Vierecke[i][1], Vierecke[i][2], Vierecke[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 CreateZehnecke(): obj = c4d.BaseObject(c4d.Opolygon) obj.ResizeObject(NP, N10*5) obj.SetName('Zehnecke') # 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 # Zähler für Polygone zurücksetzen for i in xrange(N10): obj.SetPolygon(zz, c4d.CPolygon(Zehnecke[i][0], Zehnecke[i][1], Zehnecke[i][9], Zehnecke[i][9])) zz = zz + 1 obj.SetPolygon(zz, c4d.CPolygon(Zehnecke[i][1], Zehnecke[i][2], Zehnecke[i][8], Zehnecke[i][9])) zz = zz + 1 obj.SetPolygon(zz, c4d.CPolygon(Zehnecke[i][2], Zehnecke[i][3], Zehnecke[i][7], Zehnecke[i][8])) zz = zz + 1 obj.SetPolygon(zz, c4d.CPolygon(Zehnecke[i][3], Zehnecke[i][4], Zehnecke[i][6], Zehnecke[i][7])) zz = zz + 1 obj.SetPolygon(zz, c4d.CPolygon(Zehnecke[i][4], Zehnecke[i][5], Zehnecke[i][6], Zehnecke[i][6])) 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 CreateSplineVierecke(): obj = c4d.BaseObject(c4d.Ospline) obj.SetName("Spline-Vierecke") obj.ResizeObject(N4*4) zz = 0 for i in xrange(N4): for j in xrange(4): x = Punkte[Vierecke[i][j]][0] y = Punkte[Vierecke[i][j]][1] z = Punkte[Vierecke[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, N4) for i in range(0, N4): obj.SetSegment(i, 4, 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 CreateSplineZehnecke(): obj = c4d.BaseObject(c4d.Ospline) obj.SetName("Spline-Zehnecke") obj.ResizeObject(N10*10) zz = 0 for i in xrange(N10): for j in xrange(10): x = Punkte[Zehnecke[i][j]][0] y = Punkte[Zehnecke[i][j]][1] z = Punkte[Zehnecke[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, N10) for i in range(0, N10): obj.SetSegment(i, 10, 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 = CreateVierecke() plyobj3 = CreateFuenfecke() plyobj4 = CreateZehnecke() splobj1 = CreateSplineDreiecke() splobj2 = CreateSplineVierecke() splobj3 = CreateSplineFuenfecke() splobj4 = CreateSplineZehnecke() 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(splobj4, nullobj3, None, True) doc.InsertObject(splobj3, nullobj3, None, True) doc.InsertObject(splobj2, nullobj3, None, True) doc.InsertObject(splobj1, nullobj3, None, True) doc.InsertObject(plyobj4, nullobj4, None, True) doc.InsertObject(plyobj3, nullobj4, 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()