""" Ikosaederstumpf 12.01.2021 www.3d-meier.de """ import c4d import math # Variablen und Konstanten Titel = 'Ikosaederstumpf' # Name NP = 60 # Anzahl Punkte N5 = 12 # Anzahl Fuenfecke N6 = 20 # Anzahl Sechsecke Radius = 3 # Radius der Eckpunktkugeln phi = (1 + math.sqrt(5))/2.0 # Goldener Schnitt a = 3 * phi b = 2 + phi c = 2 * phi d = phi e = 2*phi + 1 Punkte = [( 0, 1, a), ( 0, 1, -a), ( 0, -1, a), ( 0, -1, -a), ( 1, b, c), ( 1, b, -c), ( 1, -b, c), ( 1, -b, -c), (-1, b, c), (-1, b, -c), (-1, -b, c), (-1, -b, -c), ( d, 2, e), ( d, 2, -e), ( d, -2, e), ( d, -2, -e), (-d, 2, e), (-d, 2, -e), (-d, -2, e), (-d, -2, -e), ( 1, a, 0), ( 1, -a, 0), (-1, a, 0), (-1, -a, 0), ( b, c, 1), ( b, -c, 1), (-b, c, 1), (-b, -c, 1), ( b, c, -1), ( b, -c, -1), (-b, c, -1), (-b, -c, -1), ( 2, e, d), ( 2, -e, d), (-2, e, d), (-2, -e, d), ( 2, e, -d), ( 2, -e, -d), (-2, e, -d), (-2, -e, -d), ( a, 0, 1), (-a, 0, 1), ( a, 0, -1), (-a, 0, -1), ( c, 1, b), (-c, 1, b), ( c, 1, -b), (-c, 1, -b), ( c, -1, b), (-c, -1, b), ( c, -1, -b), (-c, -1, -b), ( e, d, 2), (-e, d, 2), ( e, d, -2), (-e, d, -2), ( e, -d, 2), (-e, -d, 2), ( e, -d, -2), (-e, -d, -2)] Fuenfecke = [(16, 0, 12, 4, 8), (57, 49, 45, 53, 41), (39, 23, 35, 27, 31), (29, 25, 33, 21, 37), (52, 44, 48, 56, 40), (46, 54, 42, 58, 50), ( 9, 5, 13, 1, 17), (36, 20, 32, 24, 28), (30, 26, 34, 22, 38), (51, 59, 43, 55, 47), ( 7, 11, 19, 3, 15), (10, 6, 14, 2, 18)] Sechsecke = [( 8, 4, 32, 20, 22, 34), (17, 1, 3, 19, 51, 47), (28, 24, 52, 40, 42, 54), (21, 33, 6, 10, 35, 23), (31, 27, 57, 41, 43, 59), (14, 48, 44, 12, 0, 2), (21, 23, 39, 11, 7, 37), (11, 39, 31, 59, 51, 19), (43, 41, 53, 26, 30, 55), (17, 47, 55, 30, 38, 9), (36, 5, 9, 38, 22, 20), (54, 46, 13, 5, 36, 28), (50, 15, 3, 1, 13, 46), (40, 56, 25, 29, 58, 42), ( 6, 33, 25, 56, 48, 14), (49, 18, 2, 0, 16, 45), (53, 45, 16, 8, 34, 26), (12, 44, 52, 24, 32, 4), (27, 35, 10, 18, 49, 57), (29, 37, 7, 15, 50, 58)] Faktor = 100 # Skalierungsfaktor k = 0.5 # 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 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 # Zähler für Polygone zurücksetzen 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 CreateSechsecke(): obj = c4d.BaseObject(c4d.Opolygon) obj.ResizeObject(NP, N6*2) obj.SetName('Sechsecke') # 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 zz = 0 # Zähler für Polygone zurücksetzen for i in xrange(N6): obj.SetPolygon(zz, c4d.CPolygon(Sechsecke[i][0], Sechsecke[i][1], Sechsecke[i][2], Sechsecke[i][5])) zz = zz + 1 obj.SetPolygon(zz, c4d.CPolygon(Sechsecke[i][2], Sechsecke[i][3], Sechsecke[i][4], Sechsecke[i][5])) zz = zz + 1 obj.Message(c4d.MSG_UPDATE) return obj #************************************************************************ def CreateSplineFuenfecke(): obj = c4d.BaseObject(c4d.Ospline) obj.SetName("Spline-Fuenfecke") 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 CreateSplineSechsecke(): obj = c4d.BaseObject(c4d.Ospline) obj.SetName("Spline-Sechsecke") obj.ResizeObject(N6*6) zz = 0 for i in xrange(N6): for j in xrange(6): x = Punkte[Sechsecke[i][j]][0] y = Punkte[Sechsecke[i][j]][1] z = Punkte[Sechsecke[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, N6) for i in range(0, N6): obj.SetSegment(i, 6, 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 = CreateFuenfecke() plyobj2 = CreateSechsecke() splobj1 = CreateSplineFuenfecke() splobj2 = CreateSplineSechsecke() 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()