Loading...
93
3 years ago

A planet orbiting a star in orbits of different eccentricities.

A long time ago I wrote a script in Unity to simulate the motion of bodies in elliptical orbits.

#Unity #Programming #Planet #Star #Orbit #Orbiting #Elliptical #Elliptic


This script calculates the position of the body as a function of time based on the formulas for mean, eccentric and true anomaly.

Since the equation for the eccentric anomaly cannot be solved algebraically, I wrote a simple numerical solver using the bisection method (RootFinder).

		
			public class Orbiter : MonoBehaviour
{
    public float period, distance, eccentricity;
    public float time = 0;
    private Transform Trs;
    private float M, E, T, r, tmp;
    private Vector2 pos = new Vector2();

    private float Ecc(float X)
    {
        return X - eccentricity * Mathf.Sin(X) - M;
    }

    void Start ()
    {
        Trs = GetComponent<Transform>();
    }
	
    void Update ()
    {
        //M
        M = ((time / period) * 2 * Mathf.PI) % (2 * Mathf.PI);
        //E
        try { E = RootFinder.Find(Ecc, 0, 2 * Mathf.PI, 1e-3f); } catch { Debug.Log("Range"); }
        //T
        tmp = Mathf.Cos(E);
        T = Mathf.Acos(eccentricity / (tmp * eccentricity - 1) - tmp / (tmp * eccentricity - 1));
        if (E > Mathf.PI) T = 2 * Mathf.PI - T;
        //R
        tmp = Mathf.Cos(T);
        r = distance / (1 + eccentricity * tmp);
        //Pos
        pos.x = r * Mathf.Sin(T);
        pos.y = r * tmp;
        Trs.localPosition = pos;
        //Time
        time += Time.deltaTime;
    }
}

public class RootFinder
{
    public static uint maxIter = 100;

    public delegate float Function(float x);

    public static float Find(Function F, float A, float B, float P)
    {
        float
            C = (A + B) / 2,
            Cv = F(C),
            Av = F(A),
            Bv = F(B);
        uint i = 0;
        if (Mathf.Abs(Av) <= P) return A;
        if (Mathf.Abs(Bv) <= P) return B;
        if (Mathf.Sign(Av) == Mathf.Sign(Bv)) throw new Exception("Bad Interval");
        while (Mathf.Abs(Cv) > P && i < maxIter)
        {
            if (Cv > 0)
            {
                if (Av > 0)
                {
                    A = C;
                    Av = Cv;
                }
                else
                {
                    B = C;
                    Bv = Cv;
                }
            }
            else
            {
                if (Av > 0)
                {
                    B = C;
                    Bv = Cv;
                }
                else
                {
                    A = C;
                    Av = Cv;
                }
            }
            C = (A + B) / 2;
            Cv = F(C);
            i++;
        }
        return C;
    }
}
		
	


0 comments

Loading...

Next up

Chibi drawing by Aikeji, modeled and animated in 3D.

Made in Blender.

#Blender #3D #3DArt #Modeling #Animation #Chibi #Aikeji

Beware of the armored frog

Kikkerstein's tank made in Blender for the War Thunder quest.

#Blender #3D #3DArt #Modeling #Animation #Kikkerstein #Kikkertank #Tank #WarThunder

It's just rain

Made in Blender using Dynamic Paint and particles.

#Blender #3D #3DArt #Animation #Rain #Storm #Window #DynamicPaint

Doki Doki Literature Club is celebrating Monika's birthday and the 8th anniversary of DDLC with exclusive stickers, avatar frames and backgrounds on Game Jolt! 

Head over to the shop to collect yours 👉 https://gamejolt.com/#shop

Spirit of the snowy forest

Made in Blender.

Anime girl model made in VRoid Studio (last image).

#Blender #3D #3DArt #Modeling #Animation #Anime #Girl #VRoid #VRoidStudio #Snow #Forest #Winter #Night #Blobcat

Play @pm-studios ' Bandit Trap Open Beta!

Grab your rocket chairs, vacuums, and toilet tornadoes! Report one more time for our final beta! Live through March 23rd 6 am CET.

Complete quests based on the game!

Play the beta: https://bit.ly/BanditTrapBeta

Lord Vader decides to delete his Minecraft world

Made in Blender.

#StarWarsArt #FanArtFriday #ArtWeeklies #Minecraft #Blender #Animation #3D #3DArt #DeathStar #Space

What's Your Favorite Video Game Console Of All Time?

PC - I've never had a console.

#GJAsks #Blender #3DArt #3D #PC #Room #Darkness

BIG NEWS: Your Bandit can wear a Game Jolt hat and shirt in Bandit Trap!

Defend your home or steal treasure in the open beta: https://bit.ly/BanditTrapBeta

It ends on March 23rd at 6 am CET!