############################## # hausuebung 05: rollercoaster and UN building ############################## import rhinoscriptsyntax as rs import random, time, sys ### import math ### sys.path.append("P:/") sys.path.append("/Users/anny.more/Documents/studium/5. semester/DM2") import DM_lib as dm reload(dm) ############################## rs.UnitSystem(4) # km = 5, meters = 4, cm = 3 etc rs.ShowGrid(None, 0) # grid > 0 = off rs.ShowGridAxes(None, 1) # y/y/z axen display > 0/1 = off/on rs.ViewDisplayMode(rs.CurrentView(), "wireframe") rs.Command("cplane w t enter", 0) # cPlane World Top dm.PointRadius(displayModeX=0, rad=3, styl=3) dm.printDisplay(state=1) rs.EnableRedraw(0) ######################################################################################### # UN BUILDING ######################################################################################### ###_________________________________________# ### basic settings for grid to fit UN_slab # ### # ! no need 2 change ! floors = H = dm.H = 40 # default=40 / incl roof slabs = L = dm.L = 11 # default=11 depth = D = dm.D = 4 # default= 4 / division in building_depth floorHeight = fH = dm.fH = 4.0 # default= 4.0 / 4.0*(H-1) = 156 meters # ############# # get gridCoords L*D*H = 1760 UnoGridCoords = dm.getUnoGridCoords() # get from DM_lib as dm #################___________________________# dm.newEmptyLayer("UNO::setUp::points", [150,150,200] ) dm.newEmptyLayer("UNO::setUp::geo", [150, 20,150] ) dm.newEmptyLayer("UNO::setUp::crv", [250, 10, 10]) dm.newEmptyLayer("UNO::setUp", [120,120,140] ) ################ lengthVec = lVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) )) ## rs.AddPoint( dm.getUnoCoord(1, 0, 0) ) depthVec = dVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 1, 0), dm.getUnoCoord(0, 0, 0) )) ## rs.AddPoint( dm.getUnoCoord(0, 1, 0) ) ################ if 1: if 0: rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(10, 0, 0) ) rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord( 0, 3, 0) ) rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord( 0, 0, 39) ) rs.CurveArrows(rs.AllObjects()[0], arrow_style=2) rs.CurveArrows(rs.AllObjects()[1], arrow_style=2) rs.CurveArrows(rs.AllObjects()[2], arrow_style=2) rs.ObjectColor(rs.AllObjects()[2], [0, 222, 0] ) rs.ObjectColor(rs.AllObjects()[1], [222, 0, 0] ) rs.ObjectColor(rs.AllObjects()[0], [0, 0, 222] ) rs.ObjectPrintWidth( rs.AllObjects()[0:3], 1.0 ) p0 = dm.getUnoCoord( 0, 0, 0) p1 = dm.getUnoCoord( 10, 0, 0) p2 = dm.getUnoCoord( 10, 3, 0) p3 = dm.getUnoCoord( 0, 3, 0) unoBaseCoords = [p0, p1, p2, p3, p0] unoBaseCrv = rs.AddCurve( unoBaseCoords, 1) unoTopCoords = [dm.getUnoCoord( 0, 0, 39), dm.getUnoCoord( 10, 0, 39), dm.getUnoCoord( 10, 3, 39), dm.getUnoCoord( 0, 3, 39), dm.getUnoCoord( 0, 0, 39), ] unoTopCrv = rs.AddCurve( unoTopCoords, 1) rs.ZoomExtents() if 1: ############################## rs.CurrentLayer("UNO::setUp::points") rs.AddPoints( UnoGridCoords ) if 0: for i in range(L): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(i, 0, 0) ), [0, 222, 0] ) for i in range(D): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, i, 0) ), [222, 0, 0] ) for i in range(H): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, 0, i) ), [0, 0, 222] ) ############################## if 0: rs.CurrentLayer("UNO::setUp::crv") dm.UN_slab(showAll=1) for i in range(8): p0 = dm.getUnoCoord( random.randint(0,L-1), random.randint(0,D-1), random.randint(0,H-1) ) p1 = dm.getUnoCoord( random.randint(0,L-1), random.randint(0,D-1), random.randint(0,H-1) ) rs.AddCurve( [p0, p1] ) dm.textDots( [p0, p1] ) rs.ZoomExtents() if 0: coords = UnoGridCoords[:] random.shuffle( coords ) rs.AddCurve( coords[0:], random.randint(1, 3) ) ################# ### abbreveation: def gC (l=1, d=2, h=3): cor = dm.getUnoCoord(l, d, h) return cor ## rs.AddPoint( gC(1,0,1) ) ################# #_____________________________here you go: dm.newEmptyLayer("myPROJ", [100,100,100]) # #cor = dm.getUnoCoord(5, 0, 1) #rs.AddPoint( cor ) #rs.AddPoint( dm.getUnoCoord(5, 0, 2) ) #rs.AddPoint( dm.getUnoCoord(5, 0, 5) ) # #for z in range(40): # rs.AddPoint( dm.getUnoCoord(5, 0, z) ) # if z%3==0: # rs.AddPoint( dm.getUnoCoord(6, 2, z) ) if 0: rs.CurrentLayer("UNO::setUp::crv") p0 = gC (l=10, d=0, h=39) rs.AddPoint( p0 ) p1 = rs.VectorAdd(p0, rs.VectorScale( dVec, 100.0 ) ) p2 = gC (5, 3, 39) p3 = gC (8, 0, 20) p4 = rs.VectorAdd(p3, rs.VectorScale( dVec, -100.0 ) ) p5 = gC (0, 0, 0) coords = [p0, p1, p2, p3, p4, p5] #rs.AddCurve( coords, 1 ) crvDiv = 44 crv = rs.AddCurve( coords, 2 ) coords = rs.DivideCurve( crv, crvDiv, create_points=0 ) #rs.AddPoints( coords ) for cor in coords: pass #rs.AddCircle( cor, 5.5 ) #rs.CurvePerpFrame() ##CurvePerpFrame(curve_id, parameter) dom = rs.CurveDomain( crv ) print "dom =", dom dom = rs.CurveDomain( crv )[1] print "dom =", dom #pX = rs.EvaluateCurve( crv, dom*0.5 ) #planeX = rs.CurvePerpFrame( crv, dom*0.5 ) #rs.AddCircle( planeX, 5.5) #paraX = rs.CurveClosestPoint( crv, coords[10] ) #planeX = rs.CurvePerpFrame( crv, paraX ) #rs.AddCircle( planeX, 5.5) circDiv = 16 allCoordsLists = [] for cor in coords: paraX = rs.CurveClosestPoint( crv, cor ) planeX = rs.CurvePerpFrame( crv, paraX ) circ = rs.AddCircle( planeX, 5.5) circCoords = rs.DivideCurve( circ, circDiv, 0) rs.DeleteObject( circ ) allCoordsLists.append( circCoords ) i = 0 for i in range(circDiv): coords = [] for list in allCoordsLists: pX = list[i] coords.append( pX ) rs.ObjectColor( rs.AddCurve( coords, 2 ), [0, 100, 100] ) #rs.ObjectColor( rs.AddCurve( coords, 1 ), [220, 100, 220] ) #____________________________:here you end ################################### dm.zA( proz=0.95 ) rs.EnableRedraw(1) ### 4 the MACs #dm.printDisplay(1) rs.CurrentLayer("Default") ######################################################################################### # ROLLERCOASTER ######################################################################################### # create main helix curve ("rollercoaster") def create_rollercoaster_helix(turns=10, height=80, base_radius=15, max_radius_variation=5): rs.CurrentLayer("UNO::setUp::crv") base_point = dm.getUnoCoord(5, 2, 0) helix_points = [] for i in range(turns * 50): angle = i * (2 * math.pi / 50) radius = base_radius + max_radius_variation * math.sin(i * 0.1) x = base_point[0] + radius * math.cos(angle) y = base_point[1] + radius * math.sin(angle) z = base_point[2] + height * i / (turns * 50) helix_points.append([x, y, z]) helix_curve = rs.AddCurve(helix_points) rs.ObjectColor(helix_curve, [100, 100, 100]) return helix_curve def add_perpendicular_structures(curve, num_divisions=15): rs.CurrentLayer("UNO::setUp::geo") allCoordsLists = [] points = rs.DivideCurve(curve, num_divisions, create_points=True) for pt in points: paraX = rs.CurveClosestPoint(curve, pt) planeX = rs.CurvePerpFrame(curve, paraX) # add a circle circ = rs.AddCircle(planeX, 5.0) # add a rectangle # rect = rs.AddRectangle(planeX, 10, 10) # color the circle rs.ObjectColor(circ, [110, 200, 110]) # divide circle into points circle_points = rs.DivideCurve(circ, 32, create_points=True) allCoordsLists.append(circle_points) # add lines for support #rs.AddLine(pt, rs.VectorAdd(pt, rs.VectorScale(rs.VectorUnitize(planeX[0]), 20))) #rs.AddLine(pt, rs.VectorAdd(pt, rs.VectorScale(rs.VectorUnitize(planeX[1]), 20))) #rs.AddLine(pt, rs.VectorAdd(pt, rs.VectorScale(rs.VectorUnitize(planeX[2]), 20))) # connecting points along the "rollercoaster" with curves for i in range(32): coords = [] for coordList in allCoordsLists: coords.append(coordList[i]) rs.ObjectColor(rs.AddCurve(coords, 3), [100, 100, 190]) rs.ObjectColor(rs.AddInterpCurve(coords, degree=3, knotstyle=3), [190, 100, 190]) # generate the rollercoaster helix and add structures rollercoaster_helix = create_rollercoaster_helix(turns=10, height=155, base_radius=50, max_radius_variation=10) add_perpendicular_structures(rollercoaster_helix) # drawing functions dm.zA(proz=0.95) rs.EnableRedraw(1)