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_