Poids/Scilab

Le programme Scilab suivant a été utilisé pour l'article Poids, pour déterminer les maximum et minimum de g et l'angle maximum de g avec la verticale.

Catégories :

Logiciel de statistiques - Statistiques - Logiciel scientifique libre - Logiciel de calcul numérique - Physique

Page(s) en rapport avec ce sujet :

La répartition de l'énergie solaire à la surface de la Terre...... l'unité de temps est le jour, c'est-à-dire la durée de la rotation de la Terre sur son axe, mais cette durée.... Evolution envisageable de Mars suivant un processus de terraformation..... Etude de la raideur d'un ressort et d'une fonction par Scilab... (source : oboulo)
Le probl`eme posé consistait, `a partir d'une terre quelconque A, B, C, ou D, ...... Scilab c, la commande relative `a la recherche du plus court chemin entre..... Le flot est dit complet si tout chemin allant de x1 ` a x n contient au..... Le résultat suivant caractérise un couplage de poids maximal en termes... (source : cours.ensem.inpl-nancy)
Une dalle FSR 1 axe réagit à la répartition de poids entre la gauche et la.... Développer sous Scilab (logiciel libre) un programme de calcul de ... (source : bepdf)

Le programme Scilab suivant a été utilisé pour l'article Poids, pour déterminer les maximum et minimum de g et l'angle maximum de g avec la verticale. Il est donné à des fins de contrôle (rectification d'une envisageable erreur).

On considère ici un modèle particulièrement simplifié de la Terre :

Terre sphérique de rayon constant R = 6 380 km ;
vitesse de rotation de 2π radians par un jour, un jour faisant 24 h de 3 600 s chacune ;
l'axe y est l'axe de rotation de la Terre, l'axe x est dans le plan de l'équateur.

// ********** Constantes
 
R = 6.38e6; // rayon de la Terre en m
omega = 2*%pi/(24*3600); //vitesse de rotation de la Terre
M = 5.98e24; // masse de la Terre en kg
G = 6.67428e-11; // constante gravitationnelle en m3kg-1s-2
g1 = G*M/R^2; // partie constante de l'accélération de la gravité
 
// ********** Fonctions
 
// rayon du point par rapport à l'axe de rotation
// à la latitude a
 
function [y] = rayon(a) 
y = R*cos(a);
endfunction
 
// accélération centrifuge
// à la latitude a
 
function [y] = accelerationcentrifuge(a) 
y = omega^2*rayon(a);
endfunction
 
// ********** Programme principal
 
alpha=0:0.0001:%pi/2; // latitude de 0 à 90°
 
// initialisation des valeurs extrêmes
anglemax=0;
latitudeanglemax=0;
gmin=10;
latitudegmin=0;
gmax=9;
latitudegmax=0;
 
for latitude=alpha
  g2=accelerationcentrifuge(latitude);
  g1x=-g1*cos(latitude); // accélération selon x due à l'attraction terrestre seule
  gx=g1x+g2;
  gy=-g1*sin(latitude);
  norme=sqrt(gx^2 + gy^2); // accélération de la gravité totale
  angle=acos((gx*g1x + gy^2)/(g1*norme)); // angle entre l'accélération de la gravité
    // et la verticale (à partir du produit scalaire)
// détermination des extrema
  if angle>anglemax then
    anglemax=angle; // écart maximum à l'axe point-centre de la Terre
    latitudeanglemax=latitude; // latitude où l'on a cet angle maximal
  end
  if norme>gmax then
    gmax=norme; // norme maximale de g
    latitudegmax=latitude; // latitude où l'on a cette norme maximale
  end
  if norme<gmin then
    gmin=norme;// norme minimale de g
    latitudegmin=latitude; // latitude où l'on a cette norme minimale
  end
end
// affichage des résultats
gmax
latitudegmax*180/%pi
gmin
latitudegmin*180/%pi
anglemax
latitudeanglemax*180/%pi

Les résultats obtenus sont les suivants :

gmax  =

   9.8053759
ans  =

   89.994481
gmin  =

   9.7716353
ans  =

   0.
anglemax  =

   0.0017235
ans  =

   44.948539

On remarque que la gravité maximale est inférieure à la gravité mesurée aux pôles (9, 83 m·s^-2), certainement à cause de l'hypothèse de rayon constant.

Recherche sur Amazon (livres) :

Ce texte est issu de l'encyclopédie Wikipedia. Vous pouvez consulter sa version originale dans cette encyclopédie à l'adresse http://fr.wikipedia.org/wiki/Poids/Scilab.
Voir la liste des contributeurs.
La version présentée ici à été extraite depuis cette source le 07/04/2010.
Ce texte est disponible sous les termes de la licence de documentation libre GNU (GFDL).
La liste des définitions proposées en tête de page est une sélection parmi les résultats obtenus à l'aide de la commande "define:" de Google.
Cette page fait partie du projet Wikibis.