|
1
|
les_courbes_1D #------------
|
|
2
|
|
|
3
|
# Fonction déplacement noeuds pôles
|
|
4
|
f_dep COURBEPOLYLINEAIRE_1_D
|
|
5
|
Debut_des_coordonnees_des_points
|
|
6
|
Coordonnee dim= 2 0 0.
|
|
7
|
Coordonnee dim= 2 1 0.
|
|
8
|
Coordonnee dim= 2 2 1.
|
|
9
|
Coordonnee dim= 2 3 1.
|
|
10
|
Fin_des_coordonnees_des_points
|
|
11
|
|
|
12
|
# Fonction de pondération pour le loi_melange
|
|
13
|
f_ponderation COURBEPOLYLINEAIRE_1_D
|
|
14
|
Debut_des_coordonnees_des_points
|
|
15
|
Coordonnee dim= 2 0.5 1.
|
|
16
|
Coordonnee dim= 2 1.0 0.
|
|
17
|
Coordonnee dim= 2 2.0 0.
|
|
18
|
Fin_des_coordonnees_des_points
|
|
19
|
|
|
20
|
# Fonction temps_courant -> altitude (h) pour le pilotage du chargement
|
|
21
|
f_h(t) COURBE_EXPRESSION_LITTERALE_1D
|
|
22
|
f(x)= 35.5-x
|
|
23
|
fin_parametres_courbe_expression_litterale_
|
|
24
|
|
|
25
|
# fichier de points (COURBEPOLYLINEAIRE) contenant pressAir(h) et rhoAir(h) en fct de l'altitude h (table atmo)
|
|
26
|
< Tatmo.dat
|
|
27
|
|
|
28
|
# fonction pressAir(t)
|
|
29
|
f_p(h) F1_ROND_F2
|
|
30
|
courbe1= f_p_air(h)
|
|
31
|
courbe2= f_h(t)
|
|
32
|
|
|
33
|
# fonction TempAir(t)
|
|
34
|
f_T(h) F1_ROND_F2
|
|
35
|
courbe1= f_T_air(h)
|
|
36
|
courbe2= f_h(t)
|
|
37
|
|
|
38
|
|
|
39
|
les_fonctions_nD #------------
|
|
40
|
|
|
41
|
# Fonction pour le pilotage du lambda
|
|
42
|
f_lambda FONCTION_EXPRESSION_LITTERALE_nD
|
|
43
|
un_argument= temps_courant
|
|
44
|
fct= ((temps_courant <= 1.) ? 0.554 : 2.000 )
|
|
45
|
fin_parametres_fonction_expression_litterale_
|
|
46
|
|
|
47
|
# Fonction temps_courant -> altitude (h) pour le pilotage du chargement
|
|
48
|
f_altitude FONCTION_COURBE1D
|
|
49
|
un_argument= temps_courant
|
|
50
|
courbe= f_h(t)
|
|
51
|
fin_fonction_courbe1D_
|
|
52
|
|
|
53
|
# Fonction temps_courant -> pression_air(t)
|
|
54
|
f_pressAir FONCTION_COURBE1D
|
|
55
|
un_argument= temps_courant
|
|
56
|
courbe= f_p(h)
|
|
57
|
fin_fonction_courbe1D_
|
|
58
|
|
|
59
|
# Fonction temps_courant -> Température_air(t)
|
|
60
|
f_TempAir FONCTION_COURBE1D
|
|
61
|
un_argument= temps_courant
|
|
62
|
courbe= f_T(h)
|
|
63
|
fin_fonction_courbe1D_
|
|
64
|
|
|
65
|
# calcul du volume total
|
|
66
|
le_volume FONCTION_EXPRESSION_LITTERALE_nD
|
|
67
|
deb_list_var_ X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_
|
|
68
|
fct= -(X1*N_surf_1 + X2*N_surf_2 + X3*N_surf_3)/3 * 74
|
|
69
|
fin_parametres_fonction_expression_litterale_
|
|
70
|
|
|
71
|
# gravité en fonction de l'altitude : g(h) = g0*R^2/(h+R)^2
|
|
72
|
f_g FONC_SCAL_COMBINEES_ND
|
|
73
|
fct_base= f_altitude
|
|
74
|
fin_fcts_interne_fonction_combinee_
|
|
75
|
deb_list_var_ temps_courant fin_list_var_
|
|
76
|
fct= 9.80665*(6370/(f_altitude+6370))^2
|
|
77
|
fin_parametres_fonction_combinee_
|
|
78
|
|
|
79
|
# 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
|
|
80
|
f_Ch FONC_SCAL_COMBINEES_ND
|
|
81
|
fct_base= f_TempAir
|
|
82
|
fin_fcts_interne_fonction_combinee_
|
|
83
|
deb_list_var_ temps_courant fin_list_var_
|
|
84
|
fct= 4.0026/(8314.32*f_TempAir)
|
|
85
|
fin_parametres_fonction_combinee_
|
|
86
|
|
|
87
|
# AIR - masse molaire de l'air : Ca = Ma/(R*T), évolue en fonction des incréments
|
|
88
|
f_Ca FONC_SCAL_COMBINEES_ND
|
|
89
|
fct_base= f_TempAir
|
|
90
|
fin_fcts_interne_fonction_combinee_
|
|
91
|
deb_list_var_ temps_courant fin_list_var_
|
|
92
|
fct= 28.964/(8314.32*f_TempAir)
|
|
93
|
fin_parametres_fonction_combinee_
|
|
94
|
|
|
95
|
# calcul de la masse de gaz initiale : m_h = int(rho(z)dV) = rho_h^0*int(-Ch*g*z)
|
|
96
|
# et rho_h^0 = rho_a^0*Mh/Ma car p_h = p_a
|
|
97
|
# La valeur est affichée dans le fichier Maple, variable int_vol_ballon_E_quad_fct_nD_f_mh0
|
|
98
|
f_mh0 FONC_SCAL_COMBINEES_ND
|
|
99
|
fct_base= f_rhoAir
|
|
100
|
fct_base= f_Ch
|
|
101
|
fct_base= f_g
|
|
102
|
fin_fcts_interne_fonction_combinee_
|
|
103
|
deb_list_var_ X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_
|
|
104
|
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
|
|
105
|
fin_parametres_fonction_combinee_
|
|
106
|
|
|
107
|
# AIR - masse volumique de l'air : rhoAir = p_air*Ca = p_air*Ma/(R*T)
|
|
108
|
f_rhoAir FONC_SCAL_COMBINEES_ND
|
|
109
|
fct_base= f_pressAir
|
|
110
|
fct_base= f_Ca
|
|
111
|
fin_fcts_interne_fonction_combinee_
|
|
112
|
fct= f_pressAir * f_Ca
|
|
113
|
fin_parametres_fonction_combinee_
|
|
114
|
|
|
115
|
# HE - calcul de exp(-Ch*g*z)*dV pour un calcul de l'intégrale de vol.
|
|
116
|
f_integ FONC_SCAL_COMBINEES_ND
|
|
117
|
fct_base= f_Ch
|
|
118
|
fct_base= f_g
|
|
119
|
fin_fcts_interne_fonction_combinee_
|
|
120
|
deb_list_var_ X1 X2 X3 N_surf_1 N_surf_2 N_surf_3 fin_list_var_
|
|
121
|
fct= exp(-f_Ch*f_g*X3) * (-X1*N_surf_1-X2*N_surf_2-X3*N_surf_3)/3 * 74
|
|
122
|
fin_parametres_fonction_combinee_
|
|
123
|
|
|
124
|
# HE - masse volumique de l'helium : rhoHe = m0/integ(exp(-Ch*g*z)*dV)
|
|
125
|
f_rhoHe FONCTION_EXPRESSION_LITTERALE_nD
|
|
126
|
deb_list_var_ int_vol_ballon_E_quad_fct_nD_f_integ fin_list_var_
|
|
127
|
fct= 163.1457473757/int_vol_ballon_E_quad_fct_nD_f_integ
|
|
128
|
fin_parametres_fonction_expression_litterale_
|
|
129
|
|
|
130
|
# HE - p_he = rho_he / Ch
|
|
131
|
f_pressHe FONC_SCAL_COMBINEES_ND
|
|
132
|
fct_base= f_rhoHe
|
|
133
|
fct_base= f_Ch
|
|
134
|
fin_fcts_interne_fonction_combinee_
|
|
135
|
fct= f_rhoHe/f_Ch
|
|
136
|
fin_parametres_fonction_combinee_
|
|
137
|
|
|
138
|
f_loiPressHe FONC_SCAL_COMBINEES_ND
|
|
139
|
fct_base= f_pressHe
|
|
140
|
fct_base= f_rhoHe
|
|
141
|
fct_base= f_Ch
|
|
142
|
fct_base= f_g
|
|
143
|
fin_fcts_interne_fonction_combinee_
|
|
144
|
deb_list_var_ X3 fin_list_var_
|
|
145
|
fct= f_pressHe - f_rhoHe*f_g*X3
|
|
146
|
# fct= f_pressHe*exp( -f_Ch*f_g*X3)
|
|
147
|
fin_parametres_fonction_combinee_
|
|
148
|
|
|
149
|
f_loiPressAir FONC_SCAL_COMBINEES_ND
|
|
150
|
fct_base= f_pressAir
|
|
151
|
fct_base= f_rhoAir
|
|
152
|
fct_base= f_Ca
|
|
153
|
fct_base= f_g
|
|
154
|
fin_fcts_interne_fonction_combinee_
|
|
155
|
deb_list_var_ X3 fin_list_var_
|
|
156
|
fct= f_pressAir - f_rhoAir*f_g*X3
|
|
157
|
# fct= f_pressAir*exp( -f_Ca*f_g*X3)
|
|
158
|
fin_parametres_fonction_combinee_
|
|
159
|
|
|
160
|
# supression initiale pour calcul forme et m0 : Dp = rho_0^a * (1-Mh/Ma)*g*z
|
|
161
|
f_Dp0 FONC_SCAL_COMBINEES_ND
|
|
162
|
fct_base= f_rhoAir
|
|
163
|
fct_base= f_g
|
|
164
|
fin_fcts_interne_fonction_combinee_
|
|
165
|
deb_list_var_ X3 X3_t X3_t0 temps_courant fin_list_var_
|
|
166
|
# fct= f_rhoAir*(1-4.0026/28.964)*f_g*X3
|
|
167
|
fct= 0.0813655926677*(1-4.0026/28.964)*9.7*X3
|
|
168
|
fin_parametres_fonction_combinee_
|
|
169
|
|
|
170
|
f_Dp FONC_SCAL_COMBINEES_ND
|
|
171
|
fct_base= f_loiPressHe
|
|
172
|
fct_base= f_loiPressAir
|
|
173
|
fin_fcts_interne_fonction_combinee_
|
|
174
|
deb_list_var_ X3 X3_t X3_t0 temps_courant fin_list_var_
|
|
175
|
fct= f_loiPressHe - f_loiPressAir
|
|
176
|
fin_parametres_fonction_combinee_
|
|
177
|
|
|
178
|
|
|
179
|
|