Finished AB5
This commit is contained in:
@ -40,31 +40,74 @@ public class Simulation5 {
|
||||
CodeDraw cd = new CodeDraw();
|
||||
|
||||
// create solar system with 12 bodies
|
||||
NamedBody sun = new NamedBody("Sun", 1.989E30, new Vector3(0.0, 0.0, 0.0), new Vector3(0.0, 0.0, 0.0));
|
||||
NamedBody earth = new NamedBody("Earth", 5.972E24, new Vector3(-6.13135922534815E10, -1.383789852227691E11, 2.719682263474911E7), new Vector3(26832.720535473603, -11948.23168764519, 1.9948243075997851));
|
||||
NamedBody moon = new NamedBody("Moon", 7.349E22, new Vector3(-6.132484773775896E10, -1.387394951280871E11, 1.701046736294776E7), new Vector3(27916.62329282941, -12020.39526008238, -94.89703264508708));
|
||||
NamedBody mars = new NamedBody("Mars", 6.41712E23, new Vector3(-1.7923193702925848E11, 1.726665823982123E11, 7.991673845249474E9), new Vector3(-15925.78496403673, -15381.16179928219, 68.67560910598857));
|
||||
NamedBody deimos = new NamedBody("Deimos", 1.8E20, new Vector3(-1.792255010450533E11, 1.726891122683271E11, 7.990659337380297E9), new Vector3(-17100.476719804457, -15020.348656808, 631.2927851249581));
|
||||
NamedBody phobos = new NamedBody("Phobos", 1.08E20, new Vector3(-1.792253482539647E11, 1.72661109673625E11, 7.987848354800322E9), new Vector3(-14738.203714241401, -13671.17675223948, -411.0012490555253));
|
||||
NamedBody mercury = new NamedBody("Mercury", 3.301E23, new Vector3(-5.167375560011926E10, -4.217574885682655E10, 1.14808913958168E9), new Vector3(21580.25398577148, -34951.03632847389, -4835.225596525241));
|
||||
NamedBody venus = new NamedBody("Venus", 4.86747E24, new Vector3(-3.123150865740532E10, 1.0395568504115701E11, 3.173401325838074E9), new Vector3(-33748.180519629335, -10014.25141045021, 1809.94488874165));
|
||||
NamedBody vesta = new NamedBody("Vesta", 2.5908E20, new Vector3(-3.337493557929893E11, -4.7147908276077385E10, 4.1923010146878105E10), new Vector3(4440.54247538484, -19718.49074006637, 48.06573124543601));
|
||||
NamedBody pallas = new NamedBody("Pallas", 2.14E20, new Vector3(4.3452066613895575E11, -2.057319365171432E11, 1.0549957423213101E11), new Vector3(5058.947582097117, 11184.45711782372, -8183.524138259704));
|
||||
NamedBody hygiea = new NamedBody("Hygiea", 8.32E19, new Vector3(-3.983943433707043E11, 2.325833000024021E11, -2.233667695713672E10), new Vector3(-6931.864585548552, -15686.8108598699, -690.5791992347208));
|
||||
NamedBody ceres = new NamedBody("Ceres", 9.394E20, new Vector3(3.781372641419032E11, 1.96718960466285E11, -6.366459168068592E10), new Vector3(-8555.324226752316, 14718.33755980907, 2040.230135060142));
|
||||
|
||||
NamedBody sun = new NamedBody(SolSystem4.SUN_NAMED);
|
||||
NamedBody earth = new NamedBody(SolSystem4.EARTH_NAMED);
|
||||
NamedBody moon = new NamedBody(SolSystem4.MOON_NAMED);
|
||||
NamedBody mars = new NamedBody(SolSystem4.MARS_NAMED);
|
||||
NamedBody deimos = new NamedBody(SolSystem4.DEIMOS_NAMED);
|
||||
NamedBody phobos = new NamedBody(SolSystem4.PHOBOS_NAMED);
|
||||
NamedBody mercury = new NamedBody(SolSystem4.MERCURY_NAMED);
|
||||
NamedBody venus = new NamedBody(SolSystem4.VENUS_NAMED);
|
||||
NamedBody vesta = new NamedBody(SolSystem4.VESTA_NAMED);
|
||||
NamedBody pallas = new NamedBody(SolSystem4.PALLAS_NAMED);
|
||||
NamedBody hygiea = new NamedBody(SolSystem4.HYGIEA_NAMED);
|
||||
NamedBody ceres = new NamedBody(SolSystem4.CERES_NAMED);
|
||||
// create some additional bodies
|
||||
Body[] bodies = new Body[NUMBER_OF_BODIES];
|
||||
|
||||
Random random = new Random(2022);
|
||||
|
||||
for (int i = 0; i < bodies.length; i++) {
|
||||
bodies[i] = new Body(Math.abs(random.nextGaussian()) * OVERALL_SYSTEM_MASS / bodies.length,
|
||||
bodies[i] = new Body(
|
||||
Math.abs(random.nextGaussian()) * OVERALL_SYSTEM_MASS / bodies.length,
|
||||
new Vector3(0.2 * random.nextGaussian() * AU, 0.2 * random.nextGaussian() * AU, 0.2 * random.nextGaussian() * AU),
|
||||
new Vector3(0 + random.nextGaussian() * 5e3, 0 + random.nextGaussian() * 5e3, 0 + random.nextGaussian() * 5e3));
|
||||
new Vector3(0 + random.nextGaussian() * 5e3, 0 + random.nextGaussian() * 5e3, 0 + random.nextGaussian() * 5e3)
|
||||
);
|
||||
}
|
||||
|
||||
//TODO: implementation of this method according to 'Aufgabenblatt5.md'.
|
||||
// Add both, NamedBody- and Body-objects, to your simulation.
|
||||
MassiveForceHashMap forceOnBody = new MassiveForceHashMap();
|
||||
forceOnBody.put(sun, new Vector3());
|
||||
forceOnBody.put(earth, new Vector3());
|
||||
forceOnBody.put(moon, new Vector3());
|
||||
forceOnBody.put(mars, new Vector3());
|
||||
forceOnBody.put(deimos, new Vector3());
|
||||
forceOnBody.put(phobos, new Vector3());
|
||||
forceOnBody.put(mercury, new Vector3());
|
||||
forceOnBody.put(venus, new Vector3());
|
||||
forceOnBody.put(vesta, new Vector3());
|
||||
forceOnBody.put(pallas, new Vector3());
|
||||
forceOnBody.put(hygiea, new Vector3());
|
||||
forceOnBody.put(ceres, new Vector3());
|
||||
|
||||
for (Body b : bodies) {
|
||||
forceOnBody.put(b, new Vector3());
|
||||
}
|
||||
|
||||
long seconds = 0;
|
||||
while (true) {
|
||||
seconds++;
|
||||
|
||||
for (Massive b1 : forceOnBody.keyList()) {
|
||||
Vector3 force = new Vector3();
|
||||
for (Massive b2 : forceOnBody.keyList()) {
|
||||
if (b1 != b2) {
|
||||
force = force.plus(b1.gravitationalForce(b2));
|
||||
}
|
||||
}
|
||||
forceOnBody.put(b1, force);
|
||||
}
|
||||
|
||||
for (Massive body : forceOnBody.keyList()) {
|
||||
body.move(forceOnBody.get(body));
|
||||
}
|
||||
|
||||
if ((seconds % 3600) == 0) {
|
||||
cd.clear(Color.BLACK);
|
||||
for (Massive body : forceOnBody.keyList()) {
|
||||
body.draw(cd);
|
||||
}
|
||||
cd.show();
|
||||
}
|
||||
}
|
||||
}
|
||||
}
|
||||
|
Reference in New Issue
Block a user