""" Johnson Körper 21 Tunnel 31.01.2021 www.3d-meier.de """ import c4d import math # Variablen und Konstanten Titel = 'Johnson Körper 21 Tunnel' # Name NP = 40 # Anzahl Punkte N3 = 20 # Anzahl Dreiecke N4 = 15 # Anzahl Vierecke N5 = 10 # 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(5 * (5 + 2 * math.sqrt(5))) / 10 C3 = (1 + math.sqrt(5)) / 4 C4 = math.sqrt(10 * (5 + math.sqrt(5))) / 10 C5 = math.sqrt(2 * (5 + math.sqrt(5))) / 4 C6 = math.sqrt(10 * (25 + 11 * math.sqrt(5))) / 20 C7 = (3 + math.sqrt(5)) / 4 C8 = ((5 + math.sqrt(10 * (5 + math.sqrt(5)))) / 10)+0.5 C9 = math.sqrt(5 * (5 + 2 * math.sqrt(5))) / 5 C10 = math.sqrt(5 + 2 * math.sqrt(5)) / 2 C11 = (1 + math.sqrt(5)) / 2 C12 = ((5 + 2*math.sqrt(5*(5 + 2*math.sqrt(5)))) / 10)+0.5 C13 = C12 - 1 C14 = C8 - 1 Punkte = [( C11, 1.0, 0.0), ( C11, 0.0, 0.0), (-C11, 1.0, 0.0), (-C11, 0.0, 0.0), ( C7, 1.0, C5), ( C7, 0.0, C5), ( C7, 1.0, -C5), ( C7, 0.0, -C5), ( -C7, 1.0, C5), ( -C7, 0.0, C5), ( -C7, 1.0, -C5), ( -C7, 0.0, -C5), ( 0.5, 1.0, C10), ( 0.5, 0.0, C10), ( 0.5, 1.0, -C10), ( 0.5, 0.0, -C10), (-0.5, 1.0, C10), (-0.5, 0.0, C10), (-0.5, 1.0, -C10), (-0.5, 0.0, -C10), ( C7, C8, -C1), ( -C7, C8, -C1), ( C3, C8, C6), ( -C3, C8, C6), ( 0.0, C8, -C9), ( C3, C12, C0), ( -C3, C12, C0), ( 0.5, C12, -C2), (-0.5, C12, -C2), ( 0.0, C12, C4), ( C3, C13, C0), ( -C3, C13, C0), ( 0.5, C13, -C2), (-0.5, C13, -C2), ( 0.0, C13, C4), ( C7, C14, -C1), ( -C7, C14, -C1), ( C3, C14, C6), ( -C3, C14, C6), ( 0.0, C14, -C9)] Fuenfecke = [( 4, 22, 25, 20, 0), (16, 23, 29, 22, 12), ( 2, 21, 26, 23, 8), (18, 24, 28, 21, 10), ( 6, 20, 27, 24, 14), (36, 33, 39, 19, 11), (35, 7, 15, 39, 32), (37, 5, 1, 35, 30), (38, 17, 13, 37, 34), (36, 3, 9, 38, 31)] Vierecke = [( 4, 0, 1, 5), (16, 12, 13, 17), ( 2, 8, 9, 3), (18, 10, 11, 19), ( 6, 14, 15, 7), ( 4, 5, 13, 12), (16, 17, 9, 8), ( 2, 3, 11, 10), (18, 19, 15, 14), ( 6, 7, 1, 0), (25, 30, 32, 27), (27, 32, 33, 28), (28, 33, 31, 26), (26, 31, 34, 29), (29, 34, 30, 25)] Dreiecke = [( 4, 12, 22, 22), (16, 8, 23, 23), ( 2, 10, 21, 21), (18, 14, 24, 24), ( 6, 0, 20, 20), (25, 27, 20, 20), (29, 25, 22, 22), (26, 29, 23, 23), (28, 26, 21, 21), (27, 28, 24, 24), ( 3, 36, 11, 11), (17, 38, 9, 9), ( 5, 37, 13, 13), ( 7, 35, 1, 1), (19, 39, 15, 15), (34, 37, 30, 30), (31, 38, 34, 34), (33, 36, 31, 31), (32, 39, 33, 33), (30, 35, 32, 32)] 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 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-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 main(): nullobj1 = CreateNullobjekt1() nullobj2 = CreateNullobjekt2() # Ecken nullobj3 = CreateNullobjekt3() # Kanten nullobj4 = CreateNullobjekt4() # Polygone plyobj1 = CreateDreiecke() plyobj2 = CreateVierecke() plyobj3 = CreateFuenfecke() splobj1 = CreateSplineDreiecke() splobj2 = CreateSplineVierecke() splobj3 = 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(splobj3, nullobj3, None, True) doc.InsertObject(splobj2, nullobj3, None, True) doc.InsertObject(splobj1, nullobj3, 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()