sinQb

'spider evolution'

2010-11-11T13:06:00.001+01:00

during the course of my investigation of simplistic model ideas with large number of different outcomes, I found myself working on a circular stack of triangles which are defined by few simple parameters. each parameter was originally defined by the user (the coding was left within the script but deactivated), however, as I was going through different settings to get shapes which are very different from each other for testing purposes, I decided to fully randomise the script. I have settled for the boundary condition of 20 units randomly picked to speed up the process of the shape search. within a matter of several minutes I have generated some 30 random shapes which were eventually dubbed 'spider evolution'. I'm now planning to expand the script or modify it for other purposes.

.:script:.

Option Explicit

'Script written by shane gregoran

'Script version 30 August 2010 00:20:56 updated 28 October 2010

Call p04TriangularSomethings()

Sub p04TriangularSomethings()

'for this purpose all random numbers have been set to max of 20

'Dim nPts : nPts = Rhino.GetReal("how many segments", 12)

Dim nPts : nPts = Floor(rnd*20)

If nPts < 3 Then

nPts = 3

End If

Dim cp : cp = Rhino.WorldXYPlane

cp(0)(2) = Floor(rnd*10)

'Dim iR : iR = Rhino.GetReal("length of short triangulation",10)

Dim iR : iR = Floor(rnd*20)

If iR < 1 Then

iR = 1

End If

'Dim R : R = Rhino.GetReal("length of segments",15)

Dim R : R = Floor(rnd*20)

If R < 1 Then

R = 1

End If

Dim iCirc, iLength, iPts

ReDim iPts(nPts)

iCirc = Rhino.AddCircle (cp,iR)

If Rhino.IsCurve(iCirc) Then

iLength = Rhino.CurveLength(iCirc)

End If

Dim i, S, iSc

For i = 0 To (nPts-1)

iSc = iLength/nPts

If i=0 Then

S=0

Else

S = iSc*i

End If

iPts(i) = Rhino.CurveArcLengthPoint(iCirc, S)

Dim oCirc, oLength, oPts

ReDim oPts(nPts)

cp(0)(2) = Floor(rnd*10)

oCirc = Rhino.AddCircle (cp,R)

If Rhino.IsCurve(oCirc) Then

oLength = Rhino.CurveLength(oCirc)

End If

Dim outS, oSc

For i = 0 To (nPts-1)

oSc = oLength/nPts

If i=0 Then

outS=0

Else

outS = oSc*i

End If

oPts(i) = Rhino.CurveArcLengthPoint(oCirc, outS)

'Dim apxH : apxH = Rhino.GetReal("height of apex guide curvature points", 6)

Dim apxH : apxH = Floor(rnd*20)

Dim Origin, Direction, Normal, hPlane

Origin = array(0,0,apxH)

If IsArray(Origin) Then

Direction = array(0,0,apxH+1)

If IsArray(Direction) Then

Normal = VectorCreate(Direction, Origin)

Normal = VectorUnitize(Normal)

hPlane = Rhino.PlaneFromNormal(Origin, Normal)

End If

'Dim apxR : apxR = Rhino.GetReal("distance of apex point from base", 5)

Dim apxR : apxR = Floor(rnd*20)

If apxR < 1 Then

apxR = 1

End If

Dim hCirc, hLength, hPts

ReDim hPts(nPts-1)

hCirc = Rhino.AddCircle (hPlane,apxR)

If Rhino.IsCurve(hCirc) Then

hLength = Rhino.CurveLength(hCirc)

End If

Dim hS, hSc

For i = 0 To (nPts-1)

hSc = hLength/nPts

If i=0 Then

hS=hSc/2

Else

hS = (hSc/2)+(hSc*i)

End If

hPts(i) = Rhino.CurveArcLengthPoint(hCirc, hS)

Dim wo : wo = array (0,0,(apxH/2))

Dim crvInnerPTS, crvOutterPTS, crvInnerPTSstring, crvOutterPTSstring

ReDim crvOutterPTS(nPts-1), crvInnerPTS(nPts-1), crvInnerPTSstring(nPts-1), crvOutterPTSstring(nPts-1)

For i = 0 To (nPts-1)

crvOutterPTS(i) = array(wo, hPts(i), oPts(i))

crvOutterPTSstring(i) = Rhino.AddCurve (crvOutterPTS(i))

For i = 0 To (nPts-1)

If i = 0 Then

crvInnerPTS(i) = array(wo, hPts(nPts-2), iPts(i))

Else

If i = 1 Then

crvInnerPTS(i) = array(wo, hPts(nPts-1), iPts(i))

Else

crvInnerPTS(i) = array(wo, hPts(i-2), iPts(i))

End If

crvInnerPTSstring(i) = Rhino.AddCurve (crvInnerPTS(i))

Dim capCurve, capCurveString

ReDim capCurve(nPts-1), capCurveString(nPts-1)

For i = 0 To (nPts-1)

If i = (nPts-1) Then

capCurve(i) = array(oPts(i), hPts(i), iPts(1))

Else

If i = (nPts-2) Then

capCurve(i) = array(oPts(i), hPts(i), iPts(0))

Else

capCurve(i) = array(oPts(i),hPts(i), iPts(i+2))

End If

capCurveString(i) = Rhino.AddCurve (capCurve(i))

Dim LoftA, rr, LoftB, SurfA, rrr

ReDim LoftA(nPts-1), LoftB(nPts-1), SurfA(nPts-1)

For rr = 0 To nPts-1

If rr = nPts-2 Then

rrr = 0

Else

If rr = nPts-1 Then

rrr = 1

Else

rrr = rr + 2

End If

LoftA(rr) = Array(crvOutterPTSstring(rr), crvInnerPTSstring(rrr), capCurveString(rr))

rhino.DeleteObject iCirc

rhino.DeleteObject oCirc

rhino.DeleteObject hCirc

Rhino.AddEdgeSrf LoftA(rr)

End Sub

p4: gravity field visualization experiment

2010-10-29T23:44:00.000+02:00

Gfield script is an attempt to visualize the units of space-time as they are distorted by a gravity type field. given that each 'cube' represents a unit of space-time, greatly enlarged for this purpose, the distortions of each said unit will continue until the geometry can no longer represent a 'cube' type space, at which point the 'cube' will vanish. at a later time I plan to add additions to the script to uphold a change in geometry from a standard eight point geometric shape to one that can support seven or less points as well as adding mass value which will increase the strength of field generated by the attractor point for each unit of space-time that has collapsed onto the attractor point. the following video is a top view of the change in the quantity of space-time units and how they are distorted over time.

I will, given time, make another 'perspective view' video to further represent the three dimensionality of the changes taking place. script is available for download.

project 3 - HField part 3 - functionality

2010-05-06T23:23:00.002+02:00

this is after first test run with the dissipating field around a single point based attractor. this process will be refined with much larger array then the one used in this run 3*3*3. i am also going to do a complete breakdown in sections per layer for better visual representation of the way the field generates compressions of unit spaces around itself. i am also going to add a weight unit to all points and attempt to create a generation effect. stay tuned

script so far can be downloaded but it's still a work in progress
.:DOWNLOAD SCRIPT:.

project 3 - HField part 2 - derailment

2010-05-06T01:01:00.002+02:00

as i continued to modify the bellow posted images of a box space array during my application of the modified Force of Gravity equation in a field manner, a computational matrix ended up yielding some interesting results. the 5*5*5 array was turned into the new format of space, with each of the original box unit corners being scaled dependent on the distance from the attractor point. i will conduct some more test before the script is released to conduct further experiments.

as it is but it's still a work in progress
.:DOWNLOAD SCRIPT:.

project 3 - HField part 1

2010-05-03T14:57:00.002+02:00

i have been working on getting the basic point based systems with box geometry set up in such a way that it can easily be adjusted for the sample size. since numerous calculations need to be done for differential field strength, i have decided to split this project into 2, possibly 3 parts. basic system geometry is self adjusting and can easily be modified to cut the computational times - i will be working on that in the last 3rd part if i get sufficient time but eventually i plan on doing it.

these are the images of the space geometry which i'm working on easiest way to apply modifications through a field.

.:DOWNLOAD SCRIPT:.

5^3 sample

project 2 - scripting extended - BojaShell

2010-04-20T06:28:00.005+02:00

for this project i used creation of curves with user defined parameters to surface the spine of shell. the script is not user friendly until i start managing the function calls to speed up script building. until such time i am doing everything long way to understand the inner workings of scripting to the full.

'calculation x
Dim M, n, Mx
ReDim M(z), Mx(z)
For n = 0 To z
M(n) = Rhino.Distance(cPt, Pt(n))
Mx(n) = Rhino.Distance(cPt, Ptx(n))
Next

for script details and other info read more.
.::DOWNLOAD SCRIPT::.
the script will be modified for user input if i get the time, there is just too many interesting possibilities to be making pretty scripts that are fairly basic with the shapes they create. for now if you choose to play with the script do not forget to keep the creator tagged somewhere.

project 3 - BojaShell

Option Explicit
'Script written by shane gregoran
'Script copyrighted by sinQb
'Script version 17. april 2010 15:35:44

Call BojaShell()
Sub BojaShell()
Dim cPt : cPt = Array(0,0,0) 'can be user defined later
Dim cPtz : cPtz = Array(0,0,3) 'can be user defined later

'user input for the relative size will be added later
Dim Ax : Ax =20
Dim A : A = 20
Dim Ay : Ay = A*2
Dim z : z = 3 'separation on Ay
Dim z1 : z1 = Ay/z

'starting point for the arcs
Dim cirC : cirC = Array(Ax,0,-Ax)
Dim PtStrt : PtStrt = Array(Ax,0,0)

'far side end points along y
Dim Pt(), b, Ptx
ReDim Pt(z), Ptx(z)
For b = 0 To z
Pt(b) = Array (Ax*2,b*z1,0)
Ptx(b) = Array (b*z1,Ay,0)
Next

'first arc displacement - can be user defined - center point of arc
Dim distSPt1 : distSPt1 = Rhino.Distance(cPt, cirC)
Dim P : P = distSPt1*1.05 'height of central arcs

'calculation x
Dim M, n, Mx
ReDim M(z), Mx(z)
For n = 0 To z
M(n) = Rhino.Distance(cPt, Pt(n))
Mx(n) = Rhino.Distance(cPt, Ptx(n))
Next

'calculation y
Dim Gy, hy, Gx
ReDim Gy(z), Gx(z)
For hy = 0 To z
Gy(hy) = Sqr((M(hy)/2)^2 - Ax^2)
Gx(hy) = Sqr((Mx(hy)/2)^2 - A^2)
Next

'mid calculation for arcs
Dim arc1(), d, arc1x()
ReDim arc1(z), arc1x(z)
For d = 0 To z
arc1(d) = Array(Ay/2,Gy(d),-Ax+P)
arc1x(d) = Array(Gx(d),Ax*2/2,-Ax+P)
Next

'arcs from center point
Dim arrArc(), c, arrArcX()
ReDim arrArc(z), arrArcX(z)
For c = 0 To z
arrArc(c) = Rhino.AddArc3Pt(cPt,Pt(c),arc1(c))
arrArcX(c) = Rhino.AddArc3Pt(cPt, Ptx(c),arc1x(c))
Next

'making cross arcs
Dim CsA(), CsB(), Yax(), u, uu, CsAx(), CsBx()
Dim i : i = z-1
ReDim Yax(i), CsB(i), CsBx(i)
For u = 0 To i
Yax(u) = (z1/2)+(u*z1)
Next
'last digit 4 can be changed for hight
For uu = 0 To i
CsB(uu) = Array(Ax*2,Yax(uu),4)
CsBx(uu) = Array(Yax(uu),Ay,4)
Next
Dim j, jj
ReDim CsA(i), CsAx(i)
For j = 0 To i
jj = j+1
CsA(j) = Rhino.AddArc3Pt(Pt(j),Pt(jj),CsB(j))
CsAx(j) = Rhino.AddArc3Pt(Ptx(j),Ptx(jj),CsBx(j))
Next

'lets loft the thing
Dim LoftA(), r, rr, LoftB(), SurfA()
ReDim LoftA(i), LoftB(i), SurfA(i)
For r = 0 To i
rr = r + 1
LoftA(r) = Array(arrArc(r), arrArc(rr), CsA(r))
LoftB(r) = Array(arrArcX(r), arrArcX(rr), CsAx(r))
'this bit is for string set if i choose to use it in future
SurfA(r) = Rhino.AddEdgeSrf (Array(arrArc(r), arrArc(rr), CsA(r)))
Rhino.AddEdgeSrf LoftB(r)
Next

End Sub

revenge of an Entity

2010-04-15T19:34:00.001+02:00

project 1 - scripting in the first - BojaForm

2010-04-13T19:56:00.005+02:00

rendered image of the first simple VB script in Rhino. surface is created via lofted set of curves and populated with user defined spheres on x and y axis with a chosen radius.

'Get the number of points on Y axis

nY = Rhino.GetInteger("Number on Y axis", 2, 2)

If IsNull(nY) Then Exit Sub

nY = nY - 1

for script details read more.

as you can see the script itself is very simple and i have separated all sections for clearer understanding as i am still in the process of learning inner workings of scripting. for those of you who are also getting to grips with what each piece of code does, i hope you will find this useful.

Project 1 - BojaForm

Option Explicit
'Script written by Shane Gregoran
'Script copyrighted by sinQb
'Script version 1.2. april 2010 20:28:14

Call Project1BojaForm()
Sub Project1BojaForm()

'Collect data on curves to loft
Const rhObjectCurve = 4
Dim arrObjects
arrObjects = Rhino.GetObjects("select curves to loft", rhObjectCurve )
If IsArray(arrObjects) And UBound(arrObjects) > 0 Then
Rhino.AddLoftSrf arrObjects
End If

'Parametrs used for making a uv grid on surface
Dim nX, nY, strSurface
Dim u, v, a, b, arrUV(1), arrPt, R

'Get the surface object
strSurface = Rhino.GetObject("Select surface", 8)
If IsNull(strSurface) Then Exit Sub

'Get the number of points on Y axis
nY = Rhino.GetInteger("Number on Y axis", 2, 2)
If IsNull(nY) Then Print "Set to 2"
nY = 2
nY = nY - 1
End If

'Get the number of points on X axis
nX = Rhino.GetInteger("Number on X axis", 2, 2)
If IsNull(nX) Then Print "Set to 2"
nX = 2
nX = nX - 1
End If

'Choose R value of Spheres
R = Rhino.GetInteger("Radius of Spheres", 2, 1)
If IsNull(R) Then Print "Set to 1"
R=1
End If

'Get the domain of the surface
u = Rhino.SurfaceDomain(strSurface, 0)
v = Rhino.SurfaceDomain(strSurface, 1)
If Not IsArray(U) Or Not IsArray(V) Then Exit Sub

'Add the points
For a = 0 To nY
arrUV(0) = u(0) + (((u(1) - u(0)) / nY) * a)
For b = 0 To nX
arrUV(1) = v(0) + (((v(1) - v(0)) / nX) * b)
arrPt = Rhino.EvaluateSurface(strSurface, arrUV)

'replace points with Spheres of R radius
If (Rhino.IsPointOnSurface(strSurface, arrPt) = True) Then
Call Rhino.AddSphere (arrPt, R)

End If
Next
Next

End Sub

fractal invasion of an entity

2010-04-10T23:36:00.007+02:00

an interesting combination of fractal pattern breaking up an entity. made by Hana Dutina

construct

2010-04-08T02:41:00.001+02:00

in Abstract form