Loading...
93
1 year 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

To this day, I don't know how they did it

#PhotoDay

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

PC - I've never had a console.

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

My Pokémon Baking Book arrived today

#MyPokemonBakingBook #InsightEditions #Giveaway

Happy #WIPWednesday!

Are you working on a game?

Making some art?

Practicing a song?

Something else?

Tell us in the comments!

Lord Vader decides to delete his Minecraft world

Made in Blender.

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

Happy #WIPWednesday!

Are you working on a game?

Making some art?

Practicing a song?

Something else?

Tell us in the comments!

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

Happy #WIPWednesday!

Are you working on a game?

Making some art?

Practicing a song?

Something else?

Tell us in the comments!

Comparison of view transforms in Blender using Christmas candles as an example

#Blender #3D #3DArt #ViewTransform #Christmas #Candles