Antag att vi inte har någon luft i hålet så att vi inte har något luftmotstånd.
Vidare antar vi att jorden inte roterar, att den är en perfekt sfär och att den
är helt homogen (inga densitetsvariationer). Hålet vi gör är givetvis för litet
för att påverka något.
När man är nere i hålet så definierar ens position r (avstånd till centrum) en
sfär med radien r. Gravitationell påverkan från massan i det 'skal' som har
innerradie r och ytterradie r0 (jordens radie) kommer att vara 0. Detta
behändiga faktum kan bevisas om man så vill. Alltså kommer vi att påverkas
av gravitationen från en sfär med radie r som påverkar en med kraften
F = G*m*M/r^2
där G = 6.67e-11 är gravitationskonstanten
m= 80 kg, vår massa
M= sfärens massa.
Sfärens massa kan skrivas M=rho*V = rho*4*pi*r^3/3
där rho = 5515 kg/m^3 är jordens medeldensitet så att kraften är
F = 4*pi/3*G*m*rho*r.
Enligt newtons andra lag gäller
m*d^2 r /dt^2 = F = -4*pi/3*G*m*rho*r
Denna differentialekvation har sinus/cosinus funktioner som lösningar
och genom att använda r(0)=r0 (jordradien) och dr(0)/dt = 0 som begynnelse-
villkor kan vi bestämma okända konstanter. Begynnelsevillkoren motsvarar
att man vid tiden t=0 befinner sig stilla vid hålets mynning på jordytan.
Lösningen blir
r(t) = r0*cos(sqrt(4*pi/3*rho*G)*t)
som är en oscillerande funktion med perioden 2*pi/(sqrt(4*pi/3*rho*G)) vilket
blir ca 1 och en halv timme. Alltså - om du i denna modell - hoppar ned i
hålet tar det 1.5 h att komma tillbaka upp igen. Positionens och
hastighetens tidsberoende är plottade i figurerna nedan
http://surl.se/rri
Men vad händer då vi lägger till luftmotstånd? Givetvis får vi då en
retarderande kraft som kommer att vara mycket viktig. Sluthastigheten
vid fritt fall vid jordytan är ca 55 m/s - jämför det med de ca 8000 m/s som
uppnåddes i exemplet ovan.
Antag att luftens densitet är konstant genom hålet. Luftmotståndet ger
upphov till en kraft
Fl = 1/2*C*rho_l*A*v^2 = k*v^2
Vid sluthastisgheten balanseras denna kraft precis av tyngdkraften så vi får
mg = kv^2
vilket ger k = 0.26 kg/m (för m=80 kg och g=9.8). Vi låter nu k vara konstant
i hålet vilket mostavarar konstant luftdensitet, kropssform etc. Newtons
andra lag får nu formen
m*d^2 r /dt^2 = F = -4*pi/3*G*m*rho*r - sign(dr/dt)*k*(dr/dt)^2
där sign(x) är teckenfunktionen som också kan definieras x/|x|. Vi behöver
den för att få rätt riktning på kraften.
Denna differentialekvation löser vi numeriskt. Man finner då positionens och
hastighetens tidsberoende i nedstående figurer
http://surl.se/rrl
Man ser att det skulle ta ca 2.7 dygn att nå jordens centrum och där skulle
vi alltså stanna med mycket små oscillationer kring jämviktsläget.
Givetvis är detta en enkel modell men eftersom utgångspunkten är ganska
otrolig så kanske det räcker.