RVB ATTRACTOR SERIES 02: GUMOWSKI-MIRA

njyqqq

Gumowski-Mira attractors were developed at the CERN research centre in 1980 by I. Gumowski and C. Mira while aiming to calculate the trajectories of sub-atomic particles. They create organic patterns resembling natural/marine forms.

njyqqq

Option Explicit
‘Script written by Elif Erdine

Call GumowskiMira()
Sub GumowskiMira()


Dim i

Dim x()

Dim y()

Dim pt()

Dim maxpoints

maxpoints = 5000


Dim
B

B = 1

ReDim Preserve x(maxpoints)

ReDim Preserve y(maxpoints)

x(i) = 0

y(i) = 5

i=0

Do While (i < maxpoints)


x(i+1) = B*y(i) + GM(x(i))

y(i+1) = GM(x(i+1)) – x(i)



ReDim Preserve pt(i)

pt(i) = Array(x(i), y(i), 0)



If IsArray(pt(i)) Then

Call Rhino.AddPoint(pt(i))

Dim plane

plane = Rhino.PlaneFromNormal(pt(i), Array(0,0,1))

Call Rhino.AddCircle(plane, 0.2)

End If

i = i+1

Loop

If IsArray(pt) Then

Dim ptcloud

ptcloud = Rhino.AddPointCloud(pt)

End If

Dim ptdel

ptdel =
Rhino.GetObjects(”select points to delete”, 1, , , True)

Call Rhino.DeleteObjects(ptdel)

End Sub
Function GM (ByVal x)

Dim A

A= 0.305

GM = A*x + 2*(1-A)*x^2 / (1+x^2)
End Function

[size=0.9em]Posted:
April 16th, 2009 under
Parametric Design.
Comments:
1



.RVB ATTRACTOR SERIES 01: ROSSLER ATTRACTORRossler attractor behaves similarly to Lorenz attractor. It’s formed by 3 non-linear ordinary differential equations. (
http://en.wikipedia.org/wiki/Rossler_map
)

njyqqq

Call Rossler()
Sub Rossler()
           
Dim x, y, z
           
Dim maxpoints
           
maxpoints = 8000
           
Dim h
           
h = 0.025
           
Dim p(2)
           
Dim ptnew()
           
Dim i
           
Do While (i<maxpoints)
                       
p(0) = p(0) + h* dx(p(0), p(1), p(2))
                       
p(1) = p(1) + h* dy(p(0), p(1), p(2))
                       
p(2) = p(2) + h* dz(p(0), p(1), p(2))

                       
Rhino.AddPoint(p)
                       
ReDim Preserve ptnew(i)
                       
ptnew(i) = p
                       
i=i+1

           
Loop


           
Call Rhino.AddCurve(ptnew, 3)
End Sub


Function dx(ByVal x, ByVal y, ByVal z)
           
dx = -y-z
End Function
Function dy(ByVal x, ByVal y, ByVal z)
           
dy = x + 0.4*y
End Function
Function dz(ByVal x, ByVal y, ByVal z)
           
dz = 0.31 + z*(x-5.3)
End Function

[size=0.9em]Posted:
April 16th, 2009 under
Parametric Design.
Comments:
none



.RVB CURLICUE FRACTALI’ve recently started using Rhino Scripting as a design tool in my research. I will be documenting here my learning process through the .rvb scripts.
This study shows the Curlicue Fractal, which generates quite intricate patterns. (http://mathworld.wolfram.com/CurlicueFractal.html)

njyqqq

Call Curlicue()
Sub Curlicue()
            Dim f
            Dim i
            Dim j
            Dim pi
            pi = Rhino.Pi
            Dim sqtwo
            sqtwo = 1.4142135623730950488016887
            Dim eulers
            eulers = 0.577215664901532860606512
            Dim golden
            golden = 1.618033988749894848204586
            Dim lntwo
            lntwo = 0.69314718055994530941
            Dim e
            e= 2.71828182845904523536028747
            Dim p(2)
            Dim ptCur()
            
            For j=0 To 50000*pi Step sqtwo*pi*2
                       
                        f = f + j
                        p(0)= p(0) + Cos(f)
                        p(1)= p(1)+ Sin(f)
                        p(2)=p(2) + Cos(f)*Sin(f)
                       
                        ‘Rhino.AddPoint(p)
                        ReDim Preserve ptCur(i)
                        ptCur(i) = p
                        i=i+1
                       
            Next
            Call Rhino.addcurve(ptCur, 3)      

End Sub

Posted: April 16th, 2009 under Parametric Design.
Comments: none

stratagem

强大！！！！！！！！

zdb888

好    厉害的rvb.

qq56

来个打包的，xiexie

crowwind

厉害。。。。

wawa

very good!!!

zdb888

赞 ， 厉害 ！

wojiaxue

貌似好深奥啊

mumu7

学的好慢啊

花血殇

赞一个 谢谢楼主