|
les_courbes_1D #------------
|
|
|
|
# Fonction déplacement noeuds pôles
|
|
f_dep COURBEPOLYLINEAIRE_1_D
|
|
Debut_des_coordonnees_des_points
|
|
Coordonnee dim= 2 0 0.
|
|
Coordonnee dim= 2 1 0.
|
|
Coordonnee dim= 2 2 1.
|
|
Coordonnee dim= 2 3 1.
|
|
Fin_des_coordonnees_des_points
|
|
|
|
# Fonction de pondération pour le loi_melange
|
|
f_ponderation COURBEPOLYLINEAIRE_1_D
|
|
Debut_des_coordonnees_des_points
|
|
Coordonnee dim= 2 0.5 1.
|
|
Coordonnee dim= 2 1.0 0.
|
|
Coordonnee dim= 2 2.0 0.
|
|
Fin_des_coordonnees_des_points
|
|
|
|
# Fonction temps_courant -> altitude (h) pour le pilotage du chargement
|
|
f_h(t) COURBE_EXPRESSION_LITTERALE_1D
|
|
f(x)= 35.5-x
|
|
fin_parametres_courbe_expression_litterale_
|
|
|
|
# fichier de points (COURBEPOLYLINEAIRE) contenant pressAir(h) et rhoAir(h) en fct de l'altitude h (table atmo)
|
|
< Tatmo.dat
|
|
|
|
# fonction pressAir(t)
|
|
f_p(h) F1_ROND_F2
|
|
courbe1= f_p_air(h)
|
|
courbe2= f_h(t)
|
|
|
|
# fonction TempAir(t)
|
|
f_T(h) F1_ROND_F2
|
|
courbe1= f_T_air(h)
|
|
courbe2= f_h(t)
|
|
|
|
|
|
les_fonctions_nD #------------
|
|
|
|
# Fonction pour le pilotage du lambda
|
|
f_lambda FONCTION_EXPRESSION_LITTERALE_nD
|
|
un_argument= temps_courant
|
|
fct= ((temps_courant <= 1.) ? 0.554 : 2.000 )
|
|
fin_parametres_fonction_expression_litterale_
|
|
|
|
# Fonction temps_courant -> altitude (h) pour le pilotage du chargement
|
|
f_altitude FONCTION_COURBE1D
|
|
un_argument= temps_courant
|
|
courbe= f_h(t)
|
|
fin_fonction_courbe1D_
|
|
|
|
# Fonction temps_courant -> pression_air(t)
|
|
f_pressAir FONCTION_COURBE1D
|
|
un_argument= temps_courant
|
|
courbe= f_p(h)
|
|
fin_fonction_courbe1D_
|
|
|
|
# Fonction temps_courant -> Température_air(t)
|
|
f_TempAir FONCTION_COURBE1D
|
|
un_argument= temps_courant
|
|
courbe= f_T(h)
|
|
fin_fonction_courbe1D_
|
|
|
|
# calcul du volume total
|
|
le_volume FONCTION_EXPRESSION_LITTERALE_nD
|
|
deb_list_var_ X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_
|
|
fct= -(X1*N_surf_1 + X2*N_surf_2 + X3*N_surf_3)/3 * 74
|
|
fin_parametres_fonction_expression_litterale_
|
|
|
|
# gravité en fonction de l'altitude : g(h) = g0*R^2/(h+R)^2
|
|
f_g FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_altitude
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ temps_courant fin_list_var_
|
|
fct= 9.80665*(6370/(f_altitude+6370))^2
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# HE - masse molaire de l'hélium : Ch = Mh/(R*Tair), évolue en fonction des incréments (altitude). On suppose que T_He = T_Air
|
|
f_Ch FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_TempAir
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ temps_courant fin_list_var_
|
|
fct= 4.0026/(8314.32*f_TempAir)
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# AIR - masse molaire de l'air : Ca = Ma/(R*T), évolue en fonction des incréments
|
|
f_Ca FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_TempAir
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ temps_courant fin_list_var_
|
|
fct= 28.964/(8314.32*f_TempAir)
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# calcul de la masse de gaz initiale : m_h = int(rho(z)dV) = rho_h^0*int(-Ch*g*z)
|
|
# et rho_h^0 = rho_a^0*Mh/Ma car p_h = p_a
|
|
# La valeur est affichée dans le fichier Maple, variable int_vol_ballon_E_quad_fct_nD_f_mh0
|
|
f_mh0 FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_rhoAir
|
|
fct_base= f_Ch
|
|
fct_base= f_g
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_
|
|
fct= f_rhoAir*4.0026/28.964 * exp(-f_Ch*f_g*X3) * (-X1*N_surf_1-X2*N_surf_2-X3*N_surf_3)/3 * 74
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# AIR - masse volumique de l'air : rhoAir = p_air*Ca = p_air*Ma/(R*T)
|
|
f_rhoAir FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_pressAir
|
|
fct_base= f_Ca
|
|
fin_fcts_interne_fonction_combinee_
|
|
fct= f_pressAir * f_Ca
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# HE - calcul de exp(-Ch*g*z)*dV pour un calcul de l'intégrale de vol.
|
|
f_integ FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_Ch
|
|
fct_base= f_g
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_
|
|
fct= exp(-f_Ch*f_g*X3) * (-X1*N_surf_1-X2*N_surf_2-X3*N_surf_3)/3 * 74
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# HE - masse volumique de l'helium : rhoHe = m0/integ(exp(-Ch*g*z)*dV)
|
|
f_rhoHe FONCTION_EXPRESSION_LITTERALE_nD
|
|
deb_list_var_ int_vol_ballon_E_quad_fct_nD_f_integ fin_list_var_
|
|
fct= 163.1457473757/int_vol_ballon_E_quad_fct_nD_f_integ
|
|
fin_parametres_fonction_expression_litterale_
|
|
|
|
# HE - p_he = rho_he / Ch
|
|
f_pressHe FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_rhoHe
|
|
fct_base= f_Ch
|
|
fin_fcts_interne_fonction_combinee_
|
|
fct= f_rhoHe/f_Ch
|
|
fin_parametres_fonction_combinee_
|
|
|
|
f_loiPressHe FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_pressHe
|
|
fct_base= f_rhoHe
|
|
fct_base= f_Ch
|
|
fct_base= f_g
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ X3 fin_list_var_
|
|
fct= f_pressHe - f_rhoHe*f_g*X3
|
|
# fct= f_pressHe*exp( -f_Ch*f_g*X3)
|
|
fin_parametres_fonction_combinee_
|
|
|
|
f_loiPressAir FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_pressAir
|
|
fct_base= f_rhoAir
|
|
fct_base= f_Ca
|
|
fct_base= f_g
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ X3 fin_list_var_
|
|
fct= f_pressAir - f_rhoAir*f_g*X3
|
|
# fct= f_pressAir*exp( -f_Ca*f_g*X3)
|
|
fin_parametres_fonction_combinee_
|
|
|
|
# supression initiale pour calcul forme et m0 : Dp = rho_0^a * (1-Mh/Ma)*g*z
|
|
f_Dp0 FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_rhoAir
|
|
fct_base= f_g
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ X3 X3_t X3_t0 temps_courant fin_list_var_
|
|
# fct= f_rhoAir*(1-4.0026/28.964)*f_g*X3
|
|
fct= 0.0813655926677*(1-4.0026/28.964)*9.7*X3
|
|
fin_parametres_fonction_combinee_
|
|
|
|
f_Dp FONC_SCAL_COMBINEES_ND
|
|
fct_base= f_loiPressHe
|
|
fct_base= f_loiPressAir
|
|
fin_fcts_interne_fonction_combinee_
|
|
deb_list_var_ X3 X3_t X3_t0 temps_courant fin_list_var_
|
|
fct= f_loiPressHe - f_loiPressAir
|
|
fin_parametres_fonction_combinee_
|
|
|
|
|
|
|