CAD LOGIC at DIA Dessau

Group Valentina and Irina_Tensed Riemann Surface [diagrams and renders]

2009-09-06T21:40:00.000-07:00

Group Valentina and Irina_booklet_Tensed Riemann Surface

2009-09-05T07:09:00.000-07:00

Open publication - Free publishing - More pavilion

Booklet, Alexandra & Ban Xu Group

2009-08-17T09:15:00.001-07:00

Open publication - Free publishing - More slitting

Booklet, Mircea and Qin

2009-08-14T10:41:00.000-07:00

Open publication - Free publishing

Group_Alex_Shi_ CAD logic_DIASS_Booklet

2009-08-07T03:29:00.000-07:00

Open publication - Free publishing - More cad logic

booklet , savanna +miruna

2009-07-31T01:12:00.000-07:00

Open publication - Free publishing - More study

Andreea Nica_Masha Nazarenko_Slitting_ Booklet

2009-07-30T07:38:00.000-07:00

Open publication - Free publishing - More slitting

Group Alic' & Alec'_Quinconce_Booklet

2009-07-29T12:16:00.000-07:00

Open publication - Free publishing - More slitting

Group Mauro and Cosmos_booklet_motion_slitting

2009-07-29T06:53:00.000-07:00

Open publication - Free publishing - More slitting

Group Mauro and Cosmos_motion_physical model

2009-07-29T05:13:00.000-07:00

Group Grisha and Liu, slitting, Tensile structure, Animation

2009-07-27T10:11:00.000-07:00

HIGH FASHION: Group Anastasia Globa & Wan Da_Booklet

2009-07-27T09:17:00.000-07:00

Open publication - Free publishing - More wan da

Tensile Surface - Booklet

2009-07-27T06:20:00.000-07:00

Open publication - Free publishing - More tensile

Group Grisha and Liu, slitting, Tensile structure, Booklet

2009-07-27T05:47:00.000-07:00

Open publication - Free publishing - More dia

HIGH FASHION: Group Anastasia Globa & Wan Da_Animation

2009-07-26T02:59:00.000-07:00

Tudor&Claudia_SLITTING_Animation

2009-07-24T06:26:00.000-07:00

Esteban & Diana_Slitting_V cuts+Tension_Final Presentation

2009-07-23T19:35:00.000-07:00

Open publication - Free publishing - More dia

HIGH FASHION: Group Anastasia Globa & Wan Da_Final

2009-07-22T03:13:00.000-07:00

Group_Esteban & Diana_Slitting_V-Cuts + Tension_Surface Relaxation

2009-07-20T14:50:00.000-07:00

Weft: Tension of fibers in V-direction.

Warp: Tension of fibers in U-direction.

Group Esteban & Diana_Slitting_V cuts&Tension_Surface deformation

2009-07-20T14:19:00.000-07:00

Student work 2009 - Exhibition

2009-07-13T10:43:00.000-07:00

Group Alex and Shi / Digital model_NURBS surface relaxation

2009-07-13T08:30:00.000-07:00

Group Alex and Shi / Digital model_NURBS surface relaxation

2009-07-13T07:02:00.000-07:00

Group Alexandra and Ban Xu, Slitting, Final model, Logic and Renders

2009-07-11T12:36:00.001-07:00

Experiment No.4, Slitting, Group Alexandra and Ban Xu

2009-07-11T12:24:00.000-07:00

Group Mauro and Cosmos / Digital model

2009-06-24T08:58:00.000-07:00

Below is our first approach to a digital model. To reach this aim we follow the next steps:

1.representation of deformation on physical model in a 3D surface.

2.projection of slitting lines.

3.representing slit thickness with lines.

4.extrusion of lines along slitting lines.

Group Mauro and Cosmos / Combination of angles

2009-06-24T08:38:00.000-07:00

This stage involves combining the 3 slit angles to create our model.

Group Joshua and Deborah - Digital models 3

2009-06-24T06:37:00.001-07:00

3D Model Experiments

2009-06-23T20:38:00.000-07:00

Deborah & Joshua - Model 2

2009-06-23T19:30:00.000-07:00

Group Joshua and Deborah - Digital models 1

2009-06-23T19:13:00.000-07:00

Group Joshua and Deborah - Digital models - Animation 2

2009-06-23T18:27:00.000-07:00

Group Joshua and Deborah - Digital models - Slitting experiments

2009-06-23T18:09:00.000-07:00

Group Grisha and Liu, 3D model test02

2009-06-23T17:14:00.001-07:00

Group:Shi xinyu&Alexander Kalachev _ computer test 02

2009-06-22T15:14:00.000-07:00

Group Grisha and Liu, 3D model test

2009-06-21T15:24:00.000-07:00

Group:Shi xinyu&Alexander Kalachev _ computer test 01

2009-06-21T04:29:00.000-07:00

Nurbs Relaxation Rhino Script

2009-06-21T04:01:00.000-07:00



'Version december 3rd 2004

'Routines developed and copyrighted by Gelfling '04 aka. David Rutten
'Parts of these routines were developed while I was a intern at Prat-SA, Toulouse, France
'Surface relaxation using FDM (force density method) and VR(vibration reduction)
'(This method was written specifically for Rhinoceros 3.x and it's scripting engine.
'However only the interface section of the source code *requires* RhinoScript methods...)

'Internal mesh definitions differ from those required by the Rhinoceros 3.x Scripting engine.
'Instead of faces these routines use edges. Translation from one definition to the other
'will have to be performed prior to and after relaxation. Since only the coordinates
 'of mesh vertices (or nodes) will be altered, this should be fairly straightforward.

'Furthermore in these routines vector objects play an important role. Vectors are arrays
'of 2 3D point definitions:
' Dim vecExample(1)
' vecExample(0) = Array(0,0,0) 'Starting point of vector
' vecExample(1) = Array(1,0,1) 'Ending point of vector
'Vector variables may be bigger than these dimensions but not smaller. They always need at
'least an x, y and z component.

Option Explicit

'This function performs a number of relaxation steps on the nodes of a mesh-definition.
'arrNodes is an array of 3D-points identifying the node coordinates in euclidian space
'arrConnections is an array of arrays containing information about connections between nodes
' Every index in the array is another array containing 2 numeric values:
' 0 = the index of the node where the connection links to
' 1 = the tension of the connection.
'  Note that arrNodes and arrConnections must have the same size (every index describes a single  node)
' Also note that if connections link to non-existent nodes an overflow error will occur.
' Then note that every node needs at least one connection, although logically 2 is the absolute  minimum.
' Nodes can be connected to themselves, but this makes no sense so try to avoid it.
'arrConstrains is an array that describes the constrains a node can have in space.
' Every entry describes a single node. It can have the following values:
' Non-UUID string = the point is free to move in space if there is no object that has an ID   which is identical
'     to the supplied string. Using zero length strings to make  sure nodes are free.
' UUID-String = the point is locked to an object (can be curve, surface or polysurface) by  ClosestPoint
'    (ClosestPoint routines are not included in this routine. Instead  native Rhino 3.x script
'    methods are called. If you intend to run this script on alternate  platform you'll have to
'    translate these calls into native methods. Also if other platforms  use different object
'    identifyers (instead of strings), you should make further  alterations).
'Damping = a real number controlling the damping of the relaxation routine. Numbers higher than one
' will accelerate the relaxation, likely causing instable results. This however can be  combatted by
' limiting the maximum distance nodes are allowed to travel. Numbers between zero and one will  deccelerate
' the relaxation thus preventing bouncing solutions, but more steps will be required to reach a  stable solution
' within tolerance. Negative numbers should be avoided.
'Limit = a real number controlling the distance every node is allowed to move during each iteration.
'Summary will be filled with an array of data regarding the details of the entire relaxation process
' it will be filled with the following values:
' 0 = The sum-total of distances of all node translations
' 1 = The biggest resultant vector in the set
'The return value of the function is an array identical to arrNodes but with different coordinates
Private Function RelaxMeshNodes(ByRef arrNodes, ByRef arrConnections, ByRef arrConstraints, ByVal Damping, ByVal Limit, ByRef Summary)

        Dim i, j
 Dim arrV(), resV
 Dim newNodes
 Dim ptCP

 Dim maxResVec, curResVec
 Dim sumNodeTrans

 maxResVec = 0.0
 sumNodeTrans = 0.0

 newNodes = arrNodes
 For i = 0 To UBound(newNodes)
  If arrConstraints(i) <> "FIXED" Then
   Erase arrV
   ReDim arrV(UBound(arrConnections(i)))
   For j = 0 To UBound(arrV)
    arrV(j) = MultiplyVector(Array(newNodes(i),  newNodes(arrConnections(i)(j)(0))), _
           arrConnections(i)(j)(1))
   Next
   resV = ResultantVector(arrV)
   resV = MultiplyVector(resV, Damping)
   resV = LimitVector(resV, Limit)

   If Rhino.IsObject(arrConstraints(i)) Then
    If Rhino.IsCurve(arrConstraints(i)) Then
     ptCP = Rhino.EvaluateCurve(arrConstraints(i),  Rhino.CurveClosestPoint(arrConstraints(i), resV(1)))
    ElseIf Rhino.IsBREP(arrConstraints(i)) Then
     ptCP = Rhino.BrepClosestPoint(arrConstraints(i), resV(1))(0)
    End If

    If IsArray(ptCP) Then
     newNodes(i) = ptCP
      curResVec = PointPointDistance(newNodes(i),  arrNodes(i))
      If curResVec > maxResVec Then maxResVec = curResVec
      sumNodeTrans = sumNodeTrans + curResVec
    Else
     newNodes(i) = resV(1)
      curResVec = PointPointDistance(newNodes(i),  arrNodes(i))
      If curResVec > maxResVec Then maxResVec = curResVec
      sumNodeTrans = sumNodeTrans + curResVec
    End If
   Else
    newNodes(i) = resV(1)
      curResVec = PointPointDistance(newNodes(i), arrNodes(i))
     If curResVec > maxResVec Then maxResVec = curResVec
     sumNodeTrans = sumNodeTrans + curResVec
   End If
  End If
 Next

 Summary = Array(sumNodeTrans, maxResVec)
 RelaxMeshNodes = newNodes
End Function

Private Function PointPointDistance(ByRef arrPoint1, ByRef arrPoint2)
 PointPointDistance = (arrPoint1(0)-arrPoint2(0)) * (arrPoint1(0)-arrPoint2(0)) + _
      (arrPoint1(1)-arrPoint2(1)) *  (arrPoint1(1)-arrPoint2(1)) + _
      (arrPoint1(2)-arrPoint2(2)) *  (arrPoint1(2)-arrPoint2(2))
 PointPointDistance = Sqr(PointPointDistance)
End Function

'This function calculates the resultant vector of an array of vectors.
'The starting point of the first vector in the array will be used as grip
'for the resultant vector, if an invalid vector array is passed an error will occur
Private Function ResultantVector(ByRef arrVectors)
 Dim i
 Dim divX, divY, divZ
 Dim resX, resY, resZ
 Dim ptT(2), ptH(2)

 resX = 0.0
 resY = 0.0
 resZ = 0.0
 For i = 0 To UBound(arrVectors)
  divX = arrVectors(i)(1)(0) - arrVectors(i)(0)(0)
  divY = arrVectors(i)(1)(1) - arrVectors(i)(0)(1)
  divZ = arrVectors(i)(1)(2) - arrVectors(i)(0)(2)
  resX = resX + divX
  resY = resY + divY
  resZ = resZ + divZ
 Next
 ptT(0) = arrVectors(0)(0)(0)
 ptT(1) = arrVectors(0)(0)(1)
 ptT(2) = arrVectors(0)(0)(2)
 ptH(0) = ptT(0) + resX
 ptH(1) = ptT(1) + resY
 ptH(2) = ptT(2) + resZ

 ResultantVector = Array(ptT, ptH)
End Function

'This function will limit a vectorlength to a certain value
'vectors smaller than the specified length will remain untouched.
'Note that it is unhealthy to use negative numbers or zero for dblMaxLength
Private Function LimitVector(ByVal vecIn, ByVal dblLimitation)
 Dim vecOut(1)
 Dim l

 vecOut(0) = vecIn(0)
 vecOut(1) = vecIn(1)
 l = PointPointDistance(vecIn(0), vecIn(1))

 If l > dblLimitation Then
  vecOut(1)(0) = vecOut(0)(0) + ((vecOut(1)(0)-vecOut(0)(0)) / l) * dblLimitation
  vecOut(1)(1) = vecOut(0)(1) + ((vecOut(1)(1)-vecOut(0)(1)) / l) * dblLimitation
  vecOut(1)(2) = vecOut(0)(2) + ((vecOut(1)(2)-vecOut(0)(2)) / l) * dblLimitation
 End If
 LimitVector = vecOut
End Function

'This function will multiply a vectorlength by a specified factor.
Private Function MultiplyVector(ByVal vecIn, ByVal dblFactor)
 Dim vecOut(1)
 vecOut(0) = vecIn(0)
 vecOut(1) = vecIn(1)
 vecOut(1)(0) = vecOut(0)(0) + (vecOut(1)(0) - vecOut(0)(0)) * dblFactor
 vecOut(1)(1) = vecOut(0)(1) + (vecOut(1)(1) - vecOut(0)(1)) * dblFactor
 vecOut(1)(2) = vecOut(0)(2) + (vecOut(1)(2) - vecOut(0)(2)) * dblFactor
 MultiplyVector = vecOut
End Function

'----------------------------------------------------------------------------------------------
'----------------------------------------------------------------------------------------------
'----------------------------------------------------------------------------------------------
'---------------------Thus ends the mesh relaxation source code...-----------------------------
'--------------Below you will find an subroutine which implements the above  routine-----------
'------------and also adds interface code.  (This subroutine requires Rhinoceros  3.x)---------
'----------------------------------------------------------------------------------------------
'----------------------------------------------------------------------------------------------
'----------------------------------------------------------------------------------------------

'This function translates a Rhino Nurbs Surface definition into an FDM Mesh definition
'idSurface = The ID of the surface to convert. If the ID is invalid, Null will be returned
'   If the surface is a non-compatible FDM object, Null will be returned
'   Surface that are closed or periodic in the v-direction should not be used.  Trims will be ignored.
'arrNodes = Points to a Variant variable that will recieve the Mesh Node data
'arrConnections = Points to a Variant variable that will recieve the Mesh Connectivity data
'arrConstraints = Points to a Variant variable that will recieve the Mesh Constraint data
'WarpTension = A double representing the tension of the fabric in Warp direction (U)
'WeftTension = A double representing the tension of the fabric in Weft direction (V)
'EdgeTension = An Array containing tension factor for all 4 surface edges 0=umin; 1=umax; 2=vmin;  3=vmax
'EdgeLinks = An Array containing constraint information for all four edges. Values can be:
' "FIXED" = Fixed edge
' objID = Linked edge
'  = Free edge. String that are not "FIXED" and do not represent an object are  automatically treated as free.
'     Use "FREE" to skip a few steps in the relaxation loop
'TipLinks = An Array containing constraint information for all four tips. Values can be identical to  EdgeLinks
'   0=umin,vmin; 1=umax,vmin; 2=umin,vmax; 3=umax,vmax
'If the function succeeds True will be returned, if the function fails Null will be returned.
Private Function NURBS2FDM(ByVal idSurface, ByRef arrNodes, ByRef arrFaces, _
         ByRef arrConnections, ByRef arrConstraints, _
         ByVal WarpTension, ByVal WeftTension, ByVal  EdgeTension, _
         ByVal EdgeLinks, ByVal TipLinks)
 NURBS2FDM = Null

 Dim u, v, N, M
 Dim arrN(), arrF(), arrC(), arrL()
 Dim rhPoints
 Dim uCount, vCount
 Dim uDegree, vDegree
 Dim strChoice

 rhPoints = Rhino.SurfacePoints(idSurface)
 uCount = Rhino.SurfacePointCount(idSurface)(0)
 vCount = Rhino.SurfacePointCount(idSurface)(1)
 uDegree = Rhino.SurfaceDegree(idSurface, 0)
 vDegree = Rhino.SurfaceDegree(idSurface, 1)

 If Rhino.IsSurfacePeriodic(idSurface, 0) Then 'periodic in u-direction
   'Set mesh nodes
  ReDim arrN(vCount * (uCount-uDegree) - 1)
  N = 0
  For u = 0 To uCount - uDegree - 1
  For v = 0 To vCount - 1
   arrN(N) = rhPoints((u*vCount)+v)
   N = N+1
  Next
  Next

  ReDim arrF((vCount-1)*(uCount-uDegree)-1)
  ReDim arrC(UBound(arrN))
  ReDim arrL(UBound(arrN))
  'Set mesh faces
  N = 0
  For u = 0 To uCount-uDegree-1
  For v = 0 To vCount-2
   arrF(N) = Array((u*vCount)+v, _
       ((u*vCount)+v+vCount) Mod  ((uCount-uDegree)*vCount), _
       ((u*vCount)+v+vCount+1) Mod  ((uCount-uDegree)*vCount), _
       (u*vCount)+v+1)
   N = N+1
  Next
  Next

  'Set fabric
  For u = 0 To uCount-uDegree-1
  For v = 1 To vCount-2
   N = (u*vCount)+v
   arrC(N) = Array(Array(N+1, WeftTension), _
       Array(N-1, WeftTension), _
       Array(((N-vCount)+((uCount-uDegree)*vCount))  Mod ((uCount-uDegree)*vCount), WarpTension), _
       Array((N+vCount) Mod  ((uCount-uDegree)*vCount), WarpTension))
   arrL(N) = "FREE"
  Next
  Next

  'Set u-border edges
  For u = 0 To (uCount-uDegree-1)*vCount Step vCount
   N = u
   arrC(N) = Array(Array((N+vCount) Mod ((uCount-uDegree)*vCount),  EdgeTension(0)), _
       Array((N-vCount+((uCount-uDegree)*vCount))  Mod ((uCount-uDegree)*vCount), EdgeTension(0)), _
       Array(N+1, WeftTension))
   arrL(N) = EdgeLinks(0)

   N = u+vCount-1
   arrC(N) = Array(Array((N+vCount+((uCount-uDegree)*vCount)) Mod  ((uCount-uDegree)*vCount), EdgeTension(1)), _
       Array((N-vCount+((uCount-uDegree)*vCount))  Mod ((uCount-uDegree)*vCount), EdgeTension(1)), _
       Array(N-1, WeftTension))
   arrL(N) = EdgeLinks(1)
  Next
 ElseIf Rhino.IsSurfaceClosed(idSurface, 0) Then  'closed in u direction
  'Set mesh nodes
  ReDim arrN(vCount * (uCount-1) - 1)
  N = 0
  For u = 0 To uCount - 2
  For v = 0 To vCount - 1
   arrN(N) = rhPoints((u*vCount)+v)
   N = N+1
  Next
  Next

  ReDim arrF((vCount-1)*(uCount-1)-1)
  ReDim arrC(UBound(arrN))
  ReDim arrL(UBound(arrN))
  'Set mesh faces
  N = 0
  For u = 0 To uCount-2
  For v = 0 To vCount-2
   arrF(N) = Array((u*vCount)+v, _
       ((u*vCount)+v+vCount) Mod  ((uCount-1)*vCount), _
       ((u*vCount)+v+vCount+1) Mod  ((uCount-1)*vCount), _
       (u*vCount)+v+1)
   N = N+1
  Next
  Next

  'Set fabric
  For u = 0 To uCount-2
  For v = 1 To vCount-2
   N = (u*vCount)+v
   arrC(N) = Array(Array(N+1, WeftTension), _
       Array(N-1, WeftTension), _
       Array(((N-vCount)+((uCount-1)*vCount)) Mod  ((uCount-1)*vCount), WarpTension), _
       Array((N+vCount) Mod ((uCount-1)*vCount),  WarpTension))
   arrL(N) = "FREE"
  Next
   Next

  'Set u-border edges
  For u = 0 To (uCount-2)*vCount Step vCount
   N = u
   arrC(N) = Array(Array((N+vCount) Mod ((uCount-1)*vCount), EdgeTension(0)), _
       Array((N-vCount+((uCount-1)*vCount)) Mod  ((uCount-1)*vCount), EdgeTension(0)), _
       Array(N+1, WeftTension))
   arrL(N) = EdgeLinks(0)

   N = u+vCount-1
   arrC(N) = Array(Array((N+vCount+((uCount-1)*vCount)) Mod ((uCount-1)*vCount),  EdgeTension(1)), _
       Array((N-vCount+((uCount-1)*vCount)) Mod  ((uCount-1)*vCount), EdgeTension(1)), _
       Array(N-1, WeftTension))
   arrL(N) = EdgeLinks(1)
  Next
 Else 'rectangular patch
  'Set mesh nodes
  ReDim arrN(UBound(rhPoints))
  For N = 0 To UBound(rhPoints)
   arrN(N) = rhPoints(N)
  Next
  'Set mesh faces
  ReDim arrF((uCount-1)*(vCount-1)-1)
  N = 0
  For u = 0 To uCount-2
  For v = 0 To vCount-2
   arrF(N) = Array((u*vCount)+v, (u*vCount)+v+vCount, (u*vCount)+v+1+vCount,  (u*vCount)+v+1)
   N = N+1
  Next
  Next
  'Set connections and link-constraints
  ReDim arrC(UBound(arrN))
  ReDim arrL(UBound(arrN))

  For N = 0 To UBound(arrN)
   arrC(N) = Array(Array(N+1, WeftTension), Array(N-1, WeftTension), _
       Array(N-vCount, WarpTension), Array(N+vCount,  WarpTension))
   arrL(N) = "FREE"
  Next
   'Set u-border edge cables
  For N = 0 To vCount*(uCount-1) Step vCount
   arrC(N) = Array(Array(N+vCount, EdgeTension(0)), Array(N-vCount,  EdgeTension(0)), Array(N+1, WeftTension))
   arrL(N) = EdgeLinks(0)
   M = N+vCount-1
   arrC(M) = Array(Array(M+vCount, EdgeTension(1)), Array(M-vCount,  EdgeTension(1)), Array(M-1, WeftTension))
   arrL(M) = EdgeLinks(1)
  Next
  'Set v-border edge cables
  For N = 0 To vCount-1
   arrC(N) = Array(Array(N+vCount, WarpTension), Array(N+1, EdgeTension(2)),  Array(N-1, EdgeTension(2)))
   arrL(N) = EdgeLinks(2)
   M = N + (uCount-1)*vCount
    arrC(M) = Array(Array(M-vCount, WarpTension), Array(M-1, EdgeTension(3)),  Array(M+1, EdgeTension(3)))
   arrL(M) = EdgeLinks(3)
  Next
  'Set tips
  arrC(0) = Array(Array(1,EdgeTension(2)), Array(uCount, EdgeTension(0)))
  arrL(0) = TipLinks(0)
  arrC(vCount-1) = Array(Array(vCount-2, EdgeTension(2)), Array(2*vCount-1,  EdgeTension(1)))
  arrL(vCount-1) = TipLinks(2)
  arrC(vCount*(uCount-1)) = Array(Array(vCount*(uCount-1)+1, EdgeTension(3)),  Array(vCount*(uCount-2), EdgeTension(0)))
  arrL(vCount*(uCount-1)) = TipLinks(1)
  arrC(vCount*uCount-1) = Array(Array(vCount*uCount-2, EdgeTension(3)),  Array(vCount*(uCount-1)-1, EdgeTension(1)))
  arrL(vCount*uCount-1) = TipLinks(3)
 End If

 arrNodes = arrN
 arrFaces = arrF
 arrConnections = arrC
 arrConstraints = arrL
 NURBS2FDM = True
End Function

Private Function RecreateFDMSurface(idSurface, arrNodes)
 Dim newNodes()
 Dim u, v, N, index
 Dim Nu, Nv, Du

 Nu = Rhino.SurfacePointCount(idSurface)(0)
 Nv = Rhino.SurfacePointCount(idSurface)(1)
 Du = Rhino.SurfaceDegree(idSurface, 0)
 ReDim newNodes(Nu*Nv-1)

 If Rhino.IsSurfacePeriodic(idSurface, 0) Then
  N = 0
  For u = 0 To Nu-1
  For v = 0 To Nv-1
   index = (u*Nv+v) Mod ((Nu-Du)*Nv)
   newNodes(N) = arrNodes(index)
   N = N+1
  Next
  Next
 ElseIf Rhino.IsSurfaceClosed(idSurface, 0) Then
  N = 0
  For u = 0 To Nu-1
  For v = 0 To Nv-1
   index = (u*Nv+v) Mod ((Nu-1)*Nv)
   newNodes(N) = arrNodes(index)
   N = N+1
  Next
  Next
 Else
  For N = 0 To UBound(arrNodes)
   newNodes(N) = arrNodes(N)
  Next
 End If
 RecreateFDMSurface = Rhino.AddNurbsSurface(Rhino.SurfacePointCount(idSurface), _
             newNodes, _
             Rhino.SurfaceKnots(idSurface)(0), _
             Rhino.SurfaceKnots(idSurface)(1), _
             Rhino.SurfaceDegree(idSurface), _
             Rhino.SurfaceWeights(idSurface))
End Function

'This function displays a message at the command line, blinks it with a red background colour and  then offers a set of options.
'Might be handy to use for feedback...
Private Function FlashMessage(strMessage, arrOptions)
 Dim rgbCurrent, rgbFlash
 Dim lngFlashPause
 lngFlashPause = 175
 rgbCurrent = Rhino.AppearanceColor(12)
 rgbFlash = RGB(255,0,0)

 Rhino.Prompt strMessage
 Rhino.AppearanceColor 12, rgbFlash
 Rhino.Sleep lngFlashPause
 Rhino.AppearanceColor 12, rgbCurrent
 Rhino.Sleep lngFlashPause
 Rhino.AppearanceColor 12, rgbFlash
 Rhino.Sleep lngFlashPause
 Rhino.AppearanceColor 12, rgbCurrent

FlashMessage = Rhino.GetString(strMessage, "", arrOptions)

End Function

'This function creates 4 textdots at the surface corners and prompts the user to pick some.
'You have to remove the objects again by calling the function CLEARPREVIEW afterwards.
'The return value is an array that contains the indices of the edges or NULL on error/nopick
'0=umin,vmin 1=umax,vmin 2=umin,vmax 3=umax,vmax
Private Function GetTips(ByVal idSurface, ByVal strPrompt)
 GetTips = Null
 Dim arrDots(3), selDots
 Dim srfDomain, i, j
 Dim arrAll(), S

 If Rhino.IsSurfacePeriodic(idSurface, 0) Or Rhino.IsSurfacePeriodic(idSurface, 1) Or _
 Rhino.IsSurfaceClosed(idSurface, 0) Or Rhino.IsSurfaceClosed(idSurface, 1) Then
  FlashMessage "The surface has no corners... You'll need to fix or link the edges  instead.", Array("OK")

Else

  srfDomain = Array(Rhino.SurfaceDomain(idSurface,0), Rhino.SurfaceDomain(idSurface,1))
  arrDots(0) = Rhino.AddTextDot("A", Rhino.EvaluateSurface(idSurface,  Array(srfDomain(0)(0), srfDomain(1)(0))))
  arrDots(1) = Rhino.AddTextDot("B", Rhino.EvaluateSurface(idSurface,  Array(srfDomain(0)(1), srfDomain(1)(0))))
  arrDots(2) = Rhino.AddTextDot("C", Rhino.EvaluateSurface(idSurface,  Array(srfDomain(0)(0), srfDomain(1)(1))))
  arrDots(3) = Rhino.AddTextDot("D", Rhino.EvaluateSurface(idSurface,  Array(srfDomain(0)(1), srfDomain(1)(1))))
  Rhino.ObjectName arrDots, "GelflingRelaxationPreviewObjects"
 
  selDots = Rhino.GetObjects(strPrompt, 0, False, False, True, arrDots)
  If IsNull(selDots) Then Exit Function
  S = 0
  For i = 0 To UBound(selDots)
   For j = 0 To 3
    If selDots(i) = arrDots(j) Then
     ReDim Preserve arrAll(S)
     arrAll(S) = j
     S = S+1
    End If
   Next
  Next
  GetTips = arrAll
 End If
End Function

'This function is a wrapper for the Rhino.ExtractIsoCurve method. Differences:
'- it returns only a single ID instead of an array.
'-if the extracted edge is very short a point will be drawn instead.
Private Function ExtractIsoCurveWrapper(ByVal idSurface, ByVal arrParam, ByVal Direction)
 Dim varResult
 varResult = Rhino.ExtractIsoCurve(idSurface, arrParam, Direction)
 If IsNull(varResult) Then
  ExtractIsoCurveWrapper = Rhino.AddPoint(Rhino.EvaluateSurface(idSurface, arrParam))
  If IsNull(ExtractIsoCurveWrapper) Then ExtractIsoCurveWrapper = "null-object"
 ElseIf Rhino.CurveLength(varResult(0)) <= Rhino.UnitAbsoluteTolerance Then   Rhino.DeleteObject varResult(0)   ExtractIsoCurveWrapper = Rhino.AddPoint(Rhino.EvaluateSurface(idSurface, arrParam))   If IsNull(ExtractIsoCurveWrapper) Then ExtractIsoCurveWrapper = "null-object"  Else   ExtractIsoCurveWrapper = varResult(0)  End If End Function  'This function copies the edges of a surface and prompts the user to pick some. 'The return value is an array with edge indices or NULL on error/nopick '0=umin 1=umax 2=vmin 3=vmax Private Function GetEdges(ByVal idSurface, ByVal strPrompt)  GetEdges = Null  Dim arrEdges(3), selEdges  Dim srfDomain(1), midDomain(1)  Dim i, j  Dim arrAll(), S      srfDomain(0) = Rhino.SurfaceDomain(idSurface,0)  midDomain(0) = (srfDomain(0)(0) + srfDomain(0)(1))/2  srfDomain(1) = Rhino.SurfaceDomain(idSurface,1)  midDomain(1) = (srfDomain(1)(0) + srfDomain(1)(1))/2  If Rhino.IsSurfacePeriodic(idSurface, 1) Or Rhino.IsSurfaceClosed(idSurface, 1) Then    arrEdges(2) = ExtractIsoCurveWrapper(idSurface, Array(srfDomain(0)(0), midDomain(1)),  1)   arrEdges(3) = ExtractIsoCurveWrapper(idSurface, Array(srfDomain(0)(1), midDomain(1)),  1)  ElseIf Rhino.IsSurfacePeriodic(idSurface, 0) Or Rhino.IsSurfaceClosed(idSurface, 0) Then   arrEdges(0) = ExtractIsoCurveWrapper(idSurface, Array(midDomain(0), srfDomain(1)(0)),  0)   arrEdges(1) = ExtractIsoCurveWrapper(idSurface, Array(midDomain(0), srfDomain(1)(1)),  0)  Else   arrEdges(0) = ExtractIsoCurveWrapper(idSurface, Array(midDomain(0), srfDomain(1)(0)),  0)   arrEdges(1) = ExtractIsoCurveWrapper(idSurface, Array(midDomain(0), srfDomain(1)(1)),  0)   arrEdges(2) = ExtractIsoCurveWrapper(idSurface, Array(srfDomain(0)(0), midDomain(1)),  1)   arrEdges(3) = ExtractIsoCurveWrapper(idSurface, Array(srfDomain(0)(1), midDomain(1)),  1)  End If    For i = 0 To 3   If IsNull(arrEdges(i)) Or IsEmpty(arrEdges(i)) Then arrEdges(i) = "null-object"  Next  Rhino.ObjectName arrEdges, "GelflingRelaxationPreviewObjects"     selEdges = Rhino.GetObjects(strPrompt, 0, False, False, True, arrEdges)  If IsNull(selEdges) Then Exit Function  S = 0  For i = 0 To UBound(selEdges)   For j = 0 To 3    If selEdges(i) = arrEdges(j) Then     ReDim Preserve arrAll(S)     arrAll(S) = j     S = S+1    End If   Next  Next  GetEdges = arrAll End Function   'This function draws a set of curves on a surface, either in U/warp or V/weft 'You have to call the CLEARPREVIEW function to remove the curves again afterwards. Private Sub DrawWarpWeftPreview(ByVal strSurface, ByVal intDirection, ByVal dblDensity)  Dim uDomain, vDomain  Dim i, u, v  Dim crvAdd     uDomain = Rhino.SurfaceDomain(strSurface, 0)  vDomain = Rhino.SurfaceDomain(strSurface, 1)   Rhino.EnableRedraw False  For i = 0 To 1 Step 1/dblDensity   u = uDomain(0) + (uDomain(1)-uDomain(0))*i   v = vDomain(0) + (vDomain(1)-vDomain(0))*i   crvAdd = Rhino.ExtractIsoCurve(strSurface, Array(u,v), intDirection)   Rhino.ObjectName crvAdd, "GelflingRelaxationPreviewObjects"  Next  Rhino.EnableRedraw True End Sub  'This function removes all objects from the document that have been created by functions in this  script. 'Basically it removes all objects that are named "GelflingRelaxationPreviewObjects" Private Sub ClearPreview()  Dim allPreviewObjects  allPreviewObjects = Rhino.ObjectsByName("GelflingRelaxationPreviewObjects", False)  If IsNull(allPreviewObjects) Then Exit Sub  Rhino.DeleteObjects allPreviewObjects End Sub   'This function checks a surface to see whether it is suitablefor FDM relaxation. 'If it turns out not to be, then it will prompt the user for action. 'If the surface could not be or was not fixed NULL will be returned. 'If the surface is acceptable or was fixed the new ID will be returned Private Function CheckNURBSSurface(ByVal idSurface)  CheckNURBSSurface = Null  Dim strResult  Dim idNewSurface    idNewSurface = idSurface     If Not Rhino.IsSurface(idNewSurface) Then Exit Function  If (Rhino.IsSurfaceClosed(idSurface,0) Or Rhino.IsSurfacePeriodic(idSurface,0)) And _  (Rhino.IsSurfaceClosed(idSurface,1) Or Rhino.IsSurfacePeriodic(idSurface,1)) Then   strResult = FlashMessage("The surface you selected has no free edges. Relaxation is  not possible.", Array("OK"))    Exit Function  End If     If Rhino.SurfacePointCount(idNewSurface)(0) <= 3 Or Rhino.SurfacePointCount(idSurface)(1) <=3  Then   strResult = FlashMessage("The surface you selected does not have sufficient grips.",  Array("OK"))   Exit Function  End If   '  If Rhino.IsSurfaceRational(idSurface) Then '   strResult = FlashMessage("The surface you selected is rational. Weighted grips will  not translate correctly, what would you like to do?", _ '          Array("Unweight", "Ignore",  "Abort")) '   If IsNull(strResult) Then Exit Function '   Select Case UCase(strResult) '   Case "UNWEIGHT" '    idNewSurface = Rhino.AddNurbsSurface(Rhino.SurfacePointCount(idNewSurface), _ '              Rhino.SurfacePoints(idNewSurface), _  '              Rhino.SurfaceKnots(idNewSurface)(0), _  '               Rhino.SurfaceKnots(idNewSurface)(1), _ '              Rhino.SurfaceDegree(idNewSurface)) '    Rhino.DeleteObject idSurface '   Case "IGNORE" '    'Do not take any action '   Case Else '    Exit Function '   End Select '  End If    If Rhino.IsSurfaceTrimmed(idNewSurface) Then   strResult = FlashMessage("The surface you selected is trimmed. Trims are not  supported, what would you like to do?", Array("Untrim", "Ignore", "Abort"))    If IsNull(strResult) Then Exit Function   Select Case UCase(strResult)   Case "UNTRIM"    idNewSurface = Rhino.AddNurbsSurface(Rhino.SurfacePointCount(idNewSurface), _               Rhino.SurfacePoints(idNewSurface), _                Rhino.SurfaceKnots(idNewSurface)(0), _               Rhino.SurfaceKnots(idNewSurface)(1), _               Rhino.SurfaceDegree(idNewSurface), _               Rhino.SurfaceWeights(idNewSurface))    Rhino.DeleteObject idSurface    Case "IGNORE"    'Do not take any action   Case Else    Exit Function   End Select  End If    If Rhino.IsSurfacePeriodic(idNewSurface, 1) Then   strResult = FlashMessage("Only u-periodic surfaces are supported, what would you like  to do?", Array("Swap_UV_Directions", "Abort"))   If IsNull(strResult) Then Exit Function   If UCase(strResult) = "ABORT" Then Exit Function   Rhino.EnableRedraw vbFalse    Rhino.UnselectAllObjects    Rhino.SelectObject idNewSurface    Rhino.Command "-_Dir _SwapUV _Enter", vbFalse   Rhino.EnableRedraw vbTrue  End If   If Rhino.IsSurfaceClosed(idNewSurface, 1) Then   strResult = FlashMessage("Only u-continues surfaces are supported, what would you  like to do?", Array("Swap_UV_Directions", "Abort"))    If IsNull(strResult) Then Exit Function   If UCase(strResult) = "ABORT" Then Exit Function   Rhino.EnableRedraw vbFalse    Rhino.UnselectAllObjects    Rhino.SelectObject idNewSurface    Rhino.Command "-_Dir _SwapUV _Enter", vbFalse   Rhino.EnableRedraw vbTrue  End If    CheckNURBSSurface = idNewSurface End Function   'This is it. The main sub... Private Sub RelaxNURBS_Surface()  Dim arrNodes, arrFaces, arrConnections, arrConstraints  Dim EdgeTension(3), WarpTension, WeftTension  Dim EdgeLink(3), TipLink(3)  Dim strSurfaceID, newSurfaceID  Dim MeshID  Dim DampingFactor, Limitation, PreflightLoops, AccuracyMark, LimitationDecay  Dim Summary  Dim srfBBox, Diagonal  Dim RunAsDeveloper, DestinationFolder, FramePrefix, FrameFormat, FrameDim(1)   'Set to TRUE to save an animation sequence of the relaxation process.  RunAsDeveloper = False    strSurfaceID = Rhino.GetObject("Select a surface for relaxation", 8, True, True)  If IsNull(strSurfaceID) Then Exit Sub  strSurfaceID = CheckNURBSSurface(strSurfaceID)  If IsNull(strSurfaceID) Then   Rhino.Print "Relaxation process aborted due to improper geometry..."   Exit Sub  End If   Rhino.SurfaceIsoCurveDensity strSurfaceID, -1  If Not IsNull(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "LastOperation")) Then   EdgeTension(0) = CDbl(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation",  "EdgeTension0"))   EdgeTension(1) = CDbl(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation",  "EdgeTension1"))   EdgeTension(2) = CDbl(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation",  "EdgeTension2"))   EdgeTension(3) = CDbl(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation",  "EdgeTension3"))    WarpTension = CDbl(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation",  "WarpTension"))   WeftTension = CDbl(Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation",  "WeftTension"))   EdgeLink(0) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "EdgeLink0")   EdgeLink(1) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "EdgeLink1")   EdgeLink(2) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "EdgeLink2")   EdgeLink(3) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "EdgeLink3")   TipLink(0) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "TipLink0")   TipLink(1) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "TipLink1")   TipLink(2) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "TipLink2")   TipLink(3) = Rhino.GetObjectData(strSurfaceID, "Gelfling_Relaxation", "TipLink3")  Else   WarpTension = 1.0   WeftTension = 1.0   For i = 0 To 3    EdgeTension(i) = 10.0    EdgeLink(i) = "FREE"    TipLink(i) = "FIXED"   Next  If Rhino.IsSurfacePeriodic(strSurfaceID, 0) Or Rhino.IsSurfaceClosed(strSurfaceID, 0) Then     For i = 0 To 3     EdgeLink(i) = "FIXED"    Next   End If  End If    Rhino.Prompt "Loading geometry... please wait"  Summary = NURBS2FDM(strSurfaceID, arrNodes, arrFaces, arrConnections, arrConstraints, _       WarpTension, WeftTension, EdgeTension, EdgeLink,  TipLink)  If IsNull(Summary) Then   MsgBox "Error in building FDM mesh.", vbOKOnly Or vbCritical, "Relaxation error"   Exit Sub  End If     Dim arrOptions, strResult  Dim intObject, i, N  Dim idParent  Do   If RunAsDeveloper Then  arrOptions = Array("Tension_Factors", "Geometry_Links", "Summary", "Relax", "Relax_advanced", "Quit")    Else    arrOptions = Array("Tension_Factors", "Geometry_Links", "Relax", "Quit")   End If   strResult = Rhino.GetString("Relaxation settings", "Relax", arrOptions)   If IsNull(strResult) Then Exit Sub   Select Case Left(UCase(strResult),1)  Case "M"     If UCase(strResult) = "MELLON" Then RunAsDeveloper = True   Case "T"    Do     arrOptions = Array("Warp", "Weft", "Edges", "Return")     strResult = Rhino.GetString("Fabric tension factor properties",  "Return", arrOptions)      If IsNull(strResult) Then strResult = "Return"     Select Case Left(UCase(strResult), 2)     Case "WA"      DrawWarpWeftPreview strSurfaceID, 0, 30      strResult = Rhino.GetReal("Specify a new fabric warp  tension", WarpTension, 0.01, 1000)      If Not IsNull(strResult) Then       WarpTension = CDbl(strResult)      End If      ClearPreview     Case "WE"      DrawWarpWeftPreview strSurfaceID, 1, 30      strResult = Rhino.GetReal("Specify a new fabric weft  tension", WeftTension, 0.01, 1000)  If Not IsNull(strResult) Then        WeftTension = CDbl(strResult)      End If      ClearPreview     Case "ED"      intObject = GetEdges(strSurfaceID, "Select edges to set  tension")      If IsArray(intObject) Then       strResult = Rhino.GetReal("Specify a new edge cable  tension", EdgeTension(intObject(0)), 0.01, 1000)       If Not IsNull(strResult) Then        For i = 0 To UBound(intObject)          EdgeTension(intObject(i)) =  CDbl(strResult)        Next       End If      End If      ClearPreview     Case Else      Exit Do     End Select    Loop   Case "G"    Do     arrOptions = Array("Corners", "Edges", "Return")     strResult = Rhino.GetString("Surface geometry linking options",  "Return", arrOptions)     If IsNull(strResult) Then strResult = "Return"     Select Case Left(UCase(strResult), 1)     Case "C"      intObject = GetTips(strSurfaceID, "Select corners to  restrain")      If IsArray(intObject) Then       arrOptions = Array("Fixed", "Linked")       strResult = Rhino.GetString("Pick a node restraint  type", "Fixed", arrOptions)       If Not IsNull(strResult) Then        Select Case Left(UCase(strResult),2)        Case "LI"         idParent = Rhino.GetObject("Select a  guide object", 4+8+16, False, True)          If Not IsNull(idParent) Then          For i = 0 To  UBound(intObject)           TipLink(intObject(i))  = idParent          Next         End If        Case Else         For i = 0 To UBound(intObject)          TipLink(intObject(i)) =  "FIXED"         Next        End Select       End If      End If      ClearPreview      Rhino.UnSelectAllObjects     Case "E"      intObject = GetEdges(strSurfaceID, "Select edges to  restrain")      If IsArray(intObject) Then       arrOptions = Array("Fixed", "Free", "Linked")       strResult = Rhino.GetString("Edge restraint type?",  "Fixed", arrOptions)        If Not IsNull(strResult) Then        Select Case Left(UCase(strResult),2)        Case "FR"         For i = 0 To UBound(intObject)          EdgeLink(intObject(i)) = ""         Next        Case "LI"         idParent = Rhino.GetObject("Select a  guide object", 4+8+16, False, False)           If Not IsNull(idParent) Then          For i = 0 To  UBound(intObject)            EdgeLink(intObject(i)) = idParent           Next         End If        Case Else         For i = 0 To UBound(intObject)          EdgeLink(intObject(i)) =  "FIXED"         Next        End Select       End If      End If      ClearPreview      Rhino.UnSelectAllObjects     Case Else      Exit Do     End Select    Loop   Case "S"    N = 0    For i = 0 To UBound(arrConnections)     N = N + UBound(arrConnections(i))+1    Next    MsgBox "NURBS surface relaxation data for object; " & vbNewLine & _        "[" & strSurfaceID & "]" & vbNewLine & vbNewLine & _        "Mesh node count; " & UBound(arrNodes)+1 & vbNewLine & _        "Mesh face count; " & UBound(arrFaces)+1 & vbNewLine & _        "Mesh link count; " & N & vbNewLine & _        "Fabric warp-tension; " & WarpTension & "%" & vbNewLine & _        "Fabric weft-tension; " & WeftTension & "%" & vbNewLine & _        "Fabric edge-tension(0); " & EdgeTension(0) & "% (" & EdgeLink(0)  & ")" & vbNewLine & _        "Fabric edge-tension(1); " & EdgeTension(1) & "% (" & EdgeLink(1)  & ")" & vbNewLine & _        "Fabric edge-tension(2); " & EdgeTension(2) & "% (" & EdgeLink(2)  & ")" & vbNewLine & _        "Fabric edge-tension(3); " & EdgeTension(3) & "% (" & EdgeLink(3)  & ")" & vbNewLine & _        "Fabric corner-link(A); " & TipLink(0) & vbNewLine & _        "Fabric corner-link(B); " & TipLink(1) & vbNewLine & _        "Fabric corner-link(C); " & TipLink(2) & vbNewLine & _        "Fabric corner-link(D); " & TipLink(3) & vbNewLine & vbNewLine & _        "end of summary...", vbOKOnly, "FDM-mesh summary"   Case "R"    srfBBox = Rhino.BoundingBox(strSurfaceID)    Diagonal = PointPointDistance(srfBBox(0), srfBBox(6))        If Len(strResult) = 14 Then     AccuracyMark = Rhino.GetReal("Specify relaxation tolerance",  Rhino.UnitAbsoluteTolerance, Rhino.UnitAbsoluteTolerance/100)     DampingFactor = Rhino.GetReal("Specify a global damping factor", 1.0,  0.001, 1000.0)     Limitation = Rhino.GetReal("Specify an initial translation  limitation", Round(Diagonal * 0.05,2), AccuracyMark, Diagonal*10)     LimitationDecay = Rhino.GetReal("Specify a translation limitation  decay factor per iteration", 0.95, 0.1, 0.99)     PreflightLoops = Rhino.GetInteger("Specify the number of  preflight-iterations", 20, 1, 1000)     If IsNull(AccuracyMark) Then AccuracyMark =  Rhino.UnitAbsoluteTolerance     If IsNull(DampingFactor) Then DampingFactor = 1.0     If IsNull(Limitation) Then Limitation = Diagonal * 0.05      If IsNull(LimitationDecay) Then LimitationDecay = 0.95     If IsNull(PreflightLoops) Then PreflightLoops = 20    Else     DampingFactor = 1.0     Limitation = Diagonal * 0.05     PreflightLoops = 20     AccuracyMark = Rhino.UnitAbsoluteTolerance     LimitationDecay = 0.95    End If    Exit Do   Case "Q"    Exit Sub   Case Else    Rhino.Print "Unknown command entered... nothing done."   End Select     Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "LastOperation", CStr(Now)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeTension0",  CStr(EdgeTension(0))   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeTension1",  CStr(EdgeTension(1))   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeTension2",  CStr(EdgeTension(2))   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeTension3",  CStr(EdgeTension(3))   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "WarpTension",  CStr(WarpTension)    Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "WeftTension",  CStr(WeftTension)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeLink0", EdgeLink(0)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeLink1", EdgeLink(1)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeLink2", EdgeLink(2)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "EdgeLink3", EdgeLink(3)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "TipLink0", TipLink(0)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "TipLink1", TipLink(1)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "TipLink2", TipLink(2)   Rhino.SetObjectData strSurfaceID, "Gelfling_Relaxation", "TipLink3", TipLink(3)  Loop    Rhino.UnSelectAllObjects  Rhino.HideObject strSurfaceID    Summary = NURBS2FDM(strSurfaceID, arrNodes, arrFaces, arrConnections, arrConstraints, _       WarpTension, WeftTension, EdgeTension, EdgeLink,  TipLink)  If IsNull(Summary) Then   MsgBox "The surface could not be translated into an FDM-mesh."& vbNewLine & _    "The fact that you're seeing this warning means there" & vbNewLine & _    "is a bug in the script. Please mail your surface file" & vbNewLine & _    "(*.3dm fileformat) to david@rutten.as", vbOKOnly Or vbCritical, "Relaxation  routine bug"   Exit Sub  End If    Rhino.Command "-_AdvancedDisplay _Wires _ShowMeshWires=Yes _Enter _Enter", vbFalse  MeshID = Rhino.AddMesh(arrNodes, arrFaces)      DestinationFolder = Null  If RunAsDeveloper Then   DestinationFolder = Rhino.BrowseForFolder(,"Select a folder to save the animation  sequence. Press 'Cancel' to not save animation.", "Animation settings")   If Not IsNull(DestinationFolder) Then    FramePrefix = Rhino.StringBox("Enter a name for the sequence...",  "RelaxationFrame", "Animation settings")    If IsNull(FramePrefix) Then  DestinationFolder = Null     Else     FrameFormat = Rhino.ListBox(Array("bmp", "tga", "pcx", "jpg", "png",  "tif"), "Select a file format", "Animation settings")     If IsNull(FrameFormat) Then      DestinationFolder = Null     Else      Summary = Rhino.PropertyListBox(Array("Animation width",  "Animation height"), _               Array("1500", "1000"), "Enter frame dimensions", "Animation settings")      If IsNull(Summary) Then       DestinationFolder = Null      Else       FrameDim(0) = Summary(0)       FrameDim(1) = Summary(1)      End If     End If    End If     If Not IsNull(DestinationFolder) Then     Rhino.Command "-_TestViewCaptureToFile " & Chr(34) &  DestinationFolder & FramePrefix & "00000." & FrameFormat & Chr(34) & _            " _Width=" & FrameDim(0) & " _Height=" &  FrameDim(1) &" _DrawCplane=No _DrawWorldAxes=No _Enter", vbFalse    End If   End If  End If    'Analyse current FDM mesh properties  Rhino.Prompt "Performing test relaxation run... please wait"  RelaxMeshNodes arrNodes, arrConnections, arrConstraints, DampingFactor, Limitation, Summary  Rhino.Print "Initial cumulative node translation: " & Round(Summary(0)) &  Rhino.UnitSystemName(False, True, True) & _     "   Initial largest node vector: " & Round(Summary(1), 2) & " N"    'Perform preflight relaxation  For i = 1 To PreflightLoops   Rhino.Prompt "Performing preflight relaxation. Step " & i & " of " & PreflightLoops &  "..."   arrNodes = RelaxMeshNodes(arrNodes, arrConnections, arrConstraints, DampingFactor,  Limitation, Summary)   Rhino.EnableRedraw vbFalse    Rhino.DeleteObject MeshID    MeshID = Rhino.AddMesh(arrNodes, arrFaces)   Rhino.EnableRedraw vbTrue           If Not IsNull(DestinationFolder) Then    Rhino.Command "-_TestViewCaptureToFile " & Chr(34) & DestinationFolder &  FramePrefix & Right("00000" & i, 5) & "." & FrameFormat & Chr(34) & _            " _Width=" & FrameDim(0) & " _Height=" &  FrameDim(1) & " _DrawCplane=No _DrawWorldAxes=No _Enter", vbFalse   End If  Next    'Perform damped relaxation  Do   arrNodes = RelaxMeshNodes(arrNodes, arrConnections, arrConstraints, DampingFactor,  Limitation, Summary)   Rhino.Prompt "Iteration " & i & _       "   Damping: " & Round(DampingFactor*100, 2) & "%" & _       "   Limitation: " & Round(Limitation, 2) &  Rhino.UnitSystemName(False, True, True)    Limitation = Limitation * LimitationDecay   'DampingFactor = DampingFactor * 0.98     Rhino.EnableRedraw vbFalse    Rhino.DeleteObject MeshID    MeshID = Rhino.AddMesh(arrNodes, arrFaces)   Rhino.EnableRedraw vbTrue          If Not IsNull(DestinationFolder) Then    Rhino.Command "-_TestViewCaptureToFile " & Chr(34) & DestinationFolder &  FramePrefix & Right("00000" & i, 5) & "." & FrameFormat & Chr(34) & _            " _Width=" & FrameDim(0) & " _Height=" &  FrameDim(1) &" _DrawCplane=No _DrawWorldAxes=No _Enter", vbFalse   End If     If Limitation < i =" i+1" newsurfaceid =" RecreateFDMSurface(strSurfaceID,">

Nurbs Relaxation Rhino Script

2009-06-21T03:56:00.000-07:00

HIGH FASHION_Translation_1st step: Group Anastasia Globa & Wan Da

2009-06-20T13:38:00.000-07:00

Open publication - Free publishing - More wan da

Esteban & Diana_Slitting_V cuts&Tension_Model

2009-06-19T06:13:00.000-07:00

Open publication - Free publishing - More model

Alice & Alec', Surface Deformations, Computer Speculations

2009-06-18T03:27:00.001-07:00

Surface Deformations,
Computer Speculations, Python Scripting/Maya

Explanations will follow...

Shi xinyu&Alexander Kalachev _ slitting experment

2009-06-16T12:08:00.000-07:00

Group Alexandra and Ban Xu, Basic Digital Model, Screenshots of the process

2009-06-16T05:16:00.001-07:00

Zhu lei & Mark_ slitting experment

2009-06-16T04:51:00.000-07:00

Group Alexandra and Ban Xu, Screen Shots of the third Experiment

2009-06-15T17:12:00.000-07:00

Group Alexandra and Ban Xu, Third experiment, 3Exercises, Photos of the models

2009-06-15T16:33:00.000-07:00

Group Alexandra and Ban Xu, Slitting, Third Experiment, 3 Exercises, Physical Model and Diagrams,

2009-06-15T16:16:00.000-07:00

Group Andreea, Masha, Physical models, Slitting, V-cut studies

2009-06-15T10:29:00.000-07:00

Summarizing of paper slitting technics: Group Wanda & Anastasia

2009-06-15T02:45:00.000-07:00

4th and 5th experiment

2009-06-14T13:37:00.000-07:00

Slitting, Physical Model, Diagrams, Fabric

2009-06-14T13:16:00.000-07:00

group savanna and miruna - progress -

2009-06-13T19:49:00.001-07:00

Group Mauro and Cosmos - Observations of slitting experiments

2009-06-13T08:37:00.000-07:00

We made a series of observations on our last slitting experiments to understand the behave of the material through different angles of slitting.

Christian and Miléne | Membranes | Paper +/- Space

2009-06-13T02:10:00.000-07:00

Group Joshua and Deborah - Physical models - Slitting experiments

2009-06-12T14:50:00.002-07:00

Group Joshua and Deborah - Physical models - Slitting experiments

2009-06-12T14:50:00.001-07:00

Group Joshua and Deborah - Physical models - Slitting experiments

2009-06-12T14:48:00.002-07:00

2009-06-12T14:48:00.001-07:00

Group Joshua and Deborah - Physical models - Slitting experiments

2009-06-12T14:44:00.001-07:00

Alice & Alec', slitting, summary of investigation

2009-06-12T07:59:00.000-07:00

We try to summarize what we have explored so far, and we are moving on to another stage to see how to combine grasshopper (or other animation software) in order to model an object, which may be a dress, a lamp or a furniture, with these specific slitting effects.

Group Grisha and Liu, slitting, Tensile experiment01

2009-06-08T13:57:00.000-07:00

Based on the tensile material, we firstly test vertical cut, put and fix it to
the specific point of the frame. Also, the point is durable to move vertically,
therefore we can how the structure is changed by the movements of each point.

Group Valentina and Irina, 2nd Tensile Structure Experiment

2009-06-07T07:20:00.000-07:00

For this experiment we've studied the Riemann Surfaces to see what their qualities are.
We have generated a series of graphics by using Mathematica 7.

We have looked specifically into the Riemann Surfaces of Inverses of Sums of Two Trigonometric Functions.
After analising the solutions we have picked up one in order to make a physical model out of it and start experimenting and deforming the original mathematical model.

The chosen solution is composed out of 2 tangent trigonometrical functions. The physical model was made out of a continuous cylindrical surface on which we applied 2 V cuts.

Group Valentina and Irina, Tensile Structure_Streight Cuts

2009-06-06T20:28:00.000-07:00

Group Valentina and Irina, Inflatable Surface Experiment

2009-06-06T20:20:00.000-07:00

Group Valentina and Irina, 1st slitting experiments

2009-06-06T19:47:00.000-07:00

Slitting, Physical Model, Diagram, Fabric

2009-05-30T05:53:00.000-07:00

Group Shi, Alexander, Slitting, Physical models, experiments02

2009-05-28T10:16:00.000-07:00

Group Shi, Alexander, Slitting, Physical models, experiments01

2009-05-28T08:14:00.000-07:00

3rd experiment

2009-05-28T06:33:00.000-07:00

The purpose of this experiment was to research the techniques which could create the sensation of a continuous surface. Linked x-cuts are being created on either 2 sides or all four. The resulting v-shaped fabric elements are then rotated on the frame.
Further experiments concerning location of the cuts and proportions of the frame will follow.

Flip-Flop: Group Wanda & Anastasia

2009-05-27T16:58:00.000-07:00

Group Andreea, Masha, Physical models, Slitting, experiments

2009-05-27T10:54:00.000-07:00

Group Cosmos and Mauro - Slitting experiments

2009-05-27T10:34:00.000-07:00

Continuing with the idea of motion, we made a serie of experimentations to explore this condition on our model. In these series we variated the angle and density of slitting made on paper; after that we pushed the two ends towards each other to see how they behave. We made the same kind of experiment also with other materials.

Esteban, Diana,Slitting, V cuts, Physical Model, Stockings

2009-05-27T09:35:00.000-07:00

This is or first attempt experimenting with V cuts in a stocking and affecting it by tension to get minimal surfaces.

Group Joshua and Deborah - Physical models - Slitting experiments

2009-05-27T09:34:00.001-07:00

Experiments with vertical slitting in different materials, plastic , cardboard and recicled cotton.

Open publication - Free publishing - More physical model

2009-05-27T03:02:00.001-07:00

Mauro, Cosmos, SLITTING, 2nd Experimentation, Paper

We started our process of experimentation creating patterns of cuts on a standard size of paper and observing how they behave when subjected to pushing and stretching. Then we decide to create a model which didn´t need to be manipulated but was capable of rendering its own surface and form. To achieve this, we made two actions: fold and cut. After further experiments, we picked the model which presented some qualities of a surface in motion with the aim of exploiting this quality further in a more rigorous fashion.

Groupe A&A ... Work in Progress

2009-05-26T16:04:00.000-07:00

The four files you sent to laser cutting.
The little square (19x19cm) represents separed experiments
The big squares (40x40cm) represents composed experiemnts.

Photos of the physicals models will follow.

SLITTING_alexandraand__ban _xu__group

2009-05-26T15:59:00.000-07:00

Idea of slitting, connections and photos of phisical model

Group Liu and Grisha

2009-05-26T14:03:00.000-07:00

Slitting 01

Slitting 02

Group Joshua and Deborah - Slitting

2009-05-26T13:45:00.000-07:00

Open publication - Free publishing - More surfaces

Slitting 01

2009-05-26T11:24:00.000-07:00

Slitting, Physical Model, Diagram, Paper, Fabric

2009-05-26T06:12:00.000-07:00

group_alexandraand__ban Xu

2009-05-26T06:02:00.000-07:00

First experiments

-difrent materials : cardboard and textile material

-difrent cuts : X and I shape

Slitting, Physical model, Fabric

2009-05-26T05:54:00.000-07:00

Group Wanda & Anastasia

2009-05-26T05:48:00.001-07:00

1st and 2nd experiments

2009-05-26T05:24:00.000-07:00

The slitting experiments are dealing both with ornamental as well as with structural aspects. A number of experiments on how to get simple overlaying slitting patterns by using elastic materials have been accomplished.

The next steps are experiments within 3d structures, as well as types of structures that can emerge(dealing with cutting and rotating, creating different states or types of tension within the material).

More to follow.

Variations/Declinations

2009-05-20T03:27:00.000-07:00

2009-05-19T15:02:00.000-07:00

sinic frostwork ,Wan Da

2009-04-02T08:36:00.000-07:00

MATA4

2009-03-06T18:20:00.001-08:00

Fianally i managed to get it posted... hard job with a slow internet... :)

Making of MATA4 from Tudor Cosmatu on Vimeo.

Elements 4

2009-02-17T07:08:00.001-08:00

Elements 4

2009-02-17T07:02:00.000-08:00

"making of"

2009-02-16T07:11:00.000-08:00

The Round Table!

2009-02-14T15:05:00.000-08:00

MERCURIAL / Anthony Adelmann, Deanna Lauw, Diana Pérez, Wan Da

2009-02-14T06:43:00.000-08:00

Urban Chair/Du Xushen & Yang Dongyue

2009-02-14T05:28:00.000-08:00

2009-02-13T17:30:00.000-08:00

tex & relax / Mircea Mogan, Andreea Nica, Alexandra Virlan

2009-02-13T15:49:00.000-08:00

GIA-Lamp [model_photos] / Grisha Zotov, Irina Bogdan, Anastasia Globa

2009-02-11T09:49:00.000-08:00

ELEMENTS-4

2009-02-07T12:27:00.000-08:00

CAD LOGIC at DIA Dessau

Group Valentina and Irina_Tensed Riemann Surface [diagrams and renders]

Group Valentina and Irina_booklet_Tensed Riemann Surface

Booklet, Alexandra & Ban Xu Group

Booklet, Mircea and Qin

Group_Alex_Shi_ CAD logic_DIASS_Booklet

booklet , savanna +miruna

Andreea Nica_Masha Nazarenko_Slitting_ Booklet

Group Alic' & Alec'_Quinconce_Booklet

Group Mauro and Cosmos_booklet_motion_slitting

Group Mauro and Cosmos_motion_physical model

Group Grisha and Liu, slitting, Tensile structure, Animation

HIGH FASHION: Group Anastasia Globa & Wan Da_Booklet

Tensile Surface - Booklet

Group Grisha and Liu, slitting, Tensile structure, Booklet

HIGH FASHION: Group Anastasia Globa & Wan Da_Animation

Tudor&Claudia_SLITTING_Animation

Esteban & Diana_Slitting_V cuts+Tension_Final Presentation

HIGH FASHION: Group Anastasia Globa & Wan Da_Final

Group_Esteban & Diana_Slitting_V-Cuts + Tension_Surface Relaxation

Group Esteban & Diana_Slitting_V cuts&Tension_Surface deformation

Student work 2009 - Exhibition

Group Alex and Shi / Digital model_NURBS surface relaxation

Group Alex and Shi / Digital model_NURBS surface relaxation

Group Alexandra and Ban Xu, Slitting, Final model, Logic and Renders

Experiment No.4, Slitting, Group Alexandra and Ban Xu

Group Mauro and Cosmos / Digital model

Group Mauro and Cosmos / Combination of angles

Group Joshua and Deborah - Digital models 3

3D Model Experiments

Deborah & Joshua - Model 2

Group Joshua and Deborah - Digital models 1

Group Joshua and Deborah - Digital models - Animation 2

Group Joshua and Deborah - Digital models - Slitting experiments

Group Grisha and Liu, 3D model test02

Group:Shi xinyu&Alexander Kalachev _ computer test 02

Group Grisha and Liu, 3D model test

Group:Shi xinyu&Alexander Kalachev _ computer test 01

Nurbs Relaxation Rhino Script

Nurbs Relaxation Rhino Script

HIGH FASHION_Translation_1st step: Group Anastasia Globa & Wan Da

Esteban & Diana_Slitting_V cuts&Tension_Model

Alice & Alec', Surface Deformations, Computer Speculations

Shi xinyu&Alexander Kalachev _ slitting experment

Group Alexandra and Ban Xu, Basic Digital Model, Screenshots of the process

Zhu lei & Mark_ slitting experment

Group Alexandra and Ban Xu, Screen Shots of the third Experiment

Group Alexandra and Ban Xu, Third experiment, 3Exercises, Photos of the models

Group Alexandra and Ban Xu, Slitting, Third Experiment, 3 Exercises, Physical Model and Diagrams,

Group Andreea, Masha, Physical models, Slitting, V-cut studies

Summarizing of paper slitting technics: Group Wanda & Anastasia

4th and 5th experiment

Slitting, Physical Model, Diagrams, Fabric

group savanna and miruna - progress -

Group Mauro and Cosmos - Observations of slitting experiments

Christian and Miléne | Membranes | Paper +/- Space

Group Joshua and Deborah - Physical models - Slitting experiments

Group Joshua and Deborah - Physical models - Slitting experiments

Group Joshua and Deborah - Physical models - Slitting experiments

Group Joshua and Deborah - Physical models - Slitting experiments

Alice & Alec', slitting, summary of investigation

Group Grisha and Liu, slitting, Tensile experiment01

Group Valentina and Irina, 2nd Tensile Structure Experiment

Group Valentina and Irina, Tensile Structure_Streight Cuts

Group Valentina and Irina, Inflatable Surface Experiment

Group Valentina and Irina, 1st slitting experiments

Slitting, Physical Model, Diagram, Fabric

Group Shi, Alexander, Slitting, Physical models, experiments02

Group Shi, Alexander, Slitting, Physical models, experiments01

3rd experiment

Flip-Flop: Group Wanda & Anastasia

Group Andreea, Masha, Physical models, Slitting, experiments

Group Cosmos and Mauro - Slitting experiments

Esteban, Diana,Slitting, V cuts, Physical Model, Stockings

Group Joshua and Deborah - Physical models - Slitting experiments

Groupe A&A ... Work in Progress

SLITTING___alexandra__and__ban _xu__group

Group Liu and Grisha

Group Joshua and Deborah - Slitting

Slitting 01

SLITTING_alexandraand__ban _xu__group

group_alexandraand__ban Xu