Curva del eclipse central

para tiempos dados

Vamos a calcular la banda de totalidad para un eclipse total o anular en un instante dado.
Consideremos el eclipse total del 31 de Julio de 1981 ( ya tratado por otros procedimientos)
Los elementos Besselianos del eclipse son:

T0X0X1X2X3
40,2152090,5479364-0,0000457-0,00000803
Y0Y1Y2Y3
0,551598-0,0891738-0,00014110,00000124
d0d1d2M0M1
18,30329-0,010036-0,000004238,4210915,001874
L10L11L12tanf1tanf2
0,5429860,0001113-0,0000110,00460610,0045832
L20L21L22E. BESSELIANOS
31/7/81
-0,0031220,0001108-0,0000115

Vamos a calcular la longitud y latitud de la zona de totalidad a las 2h 18m de TU tomando como diferencia TD-TU=52seg. El tiempo t (TD) medido en horas, desde el instante de referencia T0 vale:
t=2h 18m 52s-4h=-1,6855555 horas.
Para este instante los elementos besselianos valen:
X=X0+X1*t+X2*t 2+X3*t 3=-0,708459596
Y=Y0+Y1*t+Y2*t 2+Y3*t 3=0,701498572
Para la declinación solar
d=d0+d1*t+d2*t 2=18º,320109487
M=M0+M1*t=213º,134598
L2=L20+L21*t+L22*t 2=-0,003141432
Las variaciones horarias en X e Y se obtienen derivando:
X'=X1+2*X2*t+3*X3*t 2=0,548022017
Y'=Y1+2*Y2*t+3*Y3*t 2=-0,088687567
Entonces calculemos:
[esm14a.gif]
b=Y'-p*X*send=-0,030380954
c=X'+p*Y*send=0,605755733
y1=w*Y=0,703624257
b1=w*send=0,315278134
b2=0,99664719*w*cosd=0,94899317
Ahora hay que calcular B. Si B no existe no hay eclipse central a la hora indicada:
[esm14b.gif]
Encontremos la latitud geocéntrica y el ángulo horario de las expresiones:
[esm14c.gif]

tag H=4,17040005 ---> H=256º,5159
Entonces la longitud geográfica y la latitud de la sombra se obtienen:
[esm14d.gif]

Para calcular la duración del eclipse total o anular se aplica la expresión:
L'2=L2-B*tan f2=-0,00359237
El eclipse es Total si L'2<0
El eclipse es Anular si L'2>0
a=c-p*B*cos d=0,592146
n2=a2+b2=0,351560589 --->n=0,5929254
[esm14e.gif]
del camino de eclipse de Sol total o anular.
La fórmula para calcular la anchura no es rigurosa, pero a menos que el Sol este muy bajo proporciona una aproximación muy buena.
La anchura se mide perpendicular a la dirección de la sombra por la superficie terrestre.
La altura del Sol sobre el horizonte se encuentra mediante la fórmula usual:
[esm14f.gif]
Calculemos:
L'1=L10+L11*t+L12*t 2-B*tan f1=0,5425149
La relación entre el diámetro aparente de la Luna y del Sol vale:
[esm14g.gif]



El fichero Excel calcula para este y todos los eclipses totales o anulares entre 1951 y 2200 para diferentes instantes la posición sombra, duración, anchura de la banda etc.


Hay programas en Basic que calculan lo mismo.

This HTML document was created by GT_HTML 6.0d 10/07/98 13:51.