Ejemplo n.º 1
0
        static internal double minLOSCalc(bodydef body)
        {
            var sat = GraphWindow.numSats;

            if (sat >= 3)
            {
                return((body.eqr * occmod(body)) / (Math.Cos(0.5 * (2f * Math.PI / sat))) - body.eqr);
            }
            return(0);
        }
Ejemplo n.º 2
0
 static internal double occmod(bodydef body)
 {
     if (!GraphWindow.occlusionModifiers)
     {
         return(1);
     }
     if (body.atm > 0)
     {
         return(GraphWindow.atmOcclusion);
     }
     else
     {
         return(GraphWindow.vacOcclusion);
     }
 }
Ejemplo n.º 3
0
        static double burnCalc(orbitdef s, orbitdef c, bodydef b)
        {
            double sta = 0;
            double cta = 0;

            if (c.Ap == s.Ap)
            {
                cta = 180;
            }
            double sr  = s.SMA * (1 - Math.Pow(s.e, 2)) / (1 + (s.e * Math.Cos(sta)));
            double sdv = Math.Sqrt(b.GM * ((2 / sr) - (1 / s.SMA)));

            double cr  = c.SMA * (1 - Math.Pow(c.e, 2)) / (1 + (c.e * Math.Cos(cta)));
            double cdv = Math.Sqrt(b.GM * ((2 / sr) - (1 / c.SMA)));

            return(Math.Round(100 * Math.Abs(sdv - cdv)) / 100);
        }
Ejemplo n.º 4
0
        public static void Update()
        {
            body = new bodydef(PlanetSelection.selectedBody);

            var sataltitude = GraphWindow.orbitAltitude;

            if (sataltitude > body.SOIAlt())
            {
                disabled = true;
            }
            if (body.geoAlt > 0 && body.geoSMA < body.SOI)
            {
                synchrorbit = body.geoAlt.ToString("N") + " m";
            }
            else
            {
                synchrorbit = "n/a";
            }

            minLOS = minLOSCalc(body);
            if (minLOS < body.SOIAlt() && minLOS > 0)
            {
                losorbit        = minLOS.ToString("N") + " m";
                losOrbitWarning = (minLOS < body.atm);

                constellationWarning = sataltitude * 1.01 < minLOS;
            }
            else
            {
                losorbit = "n/a";
            }


            double satcount = GraphWindow.numSats;
            double satratio = ((satcount + 1) / satcount);

            satelliteorbit = new orbitdef(sataltitude, sataltitude, body);
            carrierorbit   = new orbitdef(sataltitude, sataltitude, body);


            if (satcount > 0 && sataltitude > 0)
            {
                carrierorbit.modifyAp(newMAfromT(satelliteorbit.T * satratio, body));

                if (carrierorbit.Ap > body.SOIAlt() || GraphWindow.flipOrbit == true)
                {
                    satratio     = ((satcount - 1) / satcount);
                    carrierorbit = new orbitdef(sataltitude, sataltitude, body);
                    carrierorbit.modifyPe(newMAfromT(satelliteorbit.T * satratio, body));
                }

                period = satelliteorbit.oph;
                if (!periodEntry)
                {
                    periodHour = satelliteorbit.op_p(orbitdef.timePos.hours).ToString();
                    periodMin  = satelliteorbit.op_p(orbitdef.timePos.min).ToString();
                    periodSec  = satelliteorbit.op_p(orbitdef.timePos.sec).ToString();
                }
                /* if (carrierorbit.Ap > 0) */
                carrierAp        = carrierorbit.Ap.ToString("N1") + " m";
                carrierApWarning = (carrierorbit.Ap < body.atm || carrierorbit.Ap < 0);

                /*if (carrierorbit.Pe > 0)*/
                carrierPe = carrierorbit.Pe.ToString("N1") + " m";


                carrierPeWarning = (carrierorbit.Pe < body.atm);
                carrierPEUrgent  = (carrierorbit.Pe <= 0);

                carrierT = carrierorbit.oph;

                dBurnDV = burnCalc(satelliteorbit, carrierorbit, body);
                burnDV  = dBurnDV.ToString("N2") + " m/s";
            }


            if (period == "")
            {
                header[0] = "";
                header[1] = "";
            }
            else
            {
                string s = body.body.displayName.Substring(0, body.body.displayName.Length - 2);
                header[0]  = s;
                header[0] += " " + satcount + "-satellite constellation";
                header[1]  = "Ap " + carrierAp;
                header[1] += "     Pe " + carrierPe;
                header[1] += "     Δv " + burnDV;
            }
        }
Ejemplo n.º 5
0
 static public double newMAfromT(double T, bodydef body)
 {
     return(2 * Cbrt((body.GM * Math.Pow(T, 2f)) / 39.4784176));  // 39.4784176 = 4π^2
 }
 public orbitdef(double Ap, double Pe, bodydef body)
 {
     this.Ap   = Ap;
     this.Pe   = Pe;
     this.body = body;
 }