Soustractif atmosphérique à la vitesse de tsiolkovski


APPLICATION DE CE DEUXIÈME AMENDEMENT À LA FUSÉE DE 22,2 t



Yüklə 1,76 Mb.
səhifə7/11
tarix11.07.2018
ölçüsü1,76 Mb.
#55211
1   2   3   4   5   6   7   8   9   10   11

APPLICATION DE CE DEUXIÈME AMENDEMENT À LA FUSÉE DE 22,2 t

Pour être exhaustif, nous devons produire ici la courbe (en bleu clair) de la Vitesse Instantanée de Tsiolkovski corrigée par ce Deuxième Amendement pour la fusée de 22,2 t précédemment évoquée.

On voit que, bien que basé sur un Premier Amendement nettement divergeant en seconde moitié de propulsion (toujours en jaune), ce deuxième amendement est plus efficace il propose une vitesse instantanée satisfaisante jusqu’au 3/4 de la phase propulsive (10 % d’erreur au lieu de 30 % avec le Premier Amendement) : il semble qu’il y ait donc une reconvergence propre à ce système de calcul en chaîne :



On peut donc imaginer qu’un Troisième Amendement améliorerait encore ce résultat, et ainsi de suite…
*

COMPARAISON DE NOTRE TRAVAIL D’AMENDEMENT AVEC LES RÉSULTATS D’UNE AUTRE SOURCE

Un texte de l’US Army, DESIGN OF AERODYNAMICALLY STABILIZED FREE ROCKETS, Military Handbook, propose un abaque donnant la perte de vitesse ΔV par seconde de phase propulsive d’après les caractéristiques propulsives et balistiques de la fusée 50 :



Dans cet abaque :
ΔVboost /tb est la perte de vitesse par seconde de propulsion (en m/s par seconde). tb est la durée de propulsion (que, pour notre part, nous avons toujours nommé T). ΔVboost /tb est donc une accélération. Dans les lignes qui suivent, nous appellerons ΔVboost ΔV tout court ;

G tb est le produit, par la durée tb de propulsion, de l’accélération initiale G de la fusée exprimée en gravité terrestre g (c’est-à-dire G = P/(Minit g) si l’on appelle P la poussée initiale du moteur, Minit la masse initiale de l’engin et g la gravité terrestre). Gtb , produit d’une accélération par un temps, est donc représentatif d’une vitesse, bien qu’il ait la dimension d’un temps. Nous y reviendrons ;

Bbo est ce que le Military Handbook appelle le Burnout Ballistic Coefficient. C’est le quotient de la Masse Balistique MBal par le produit SCx de la surface de référence par le Cx. Ce Cx doit être pris immédiatement après la fin de la propulsion. Bbo s’exprime en Kg/m² ;


Quels sont les paramètres utilisés comme entrées dans cet abaque ?
 Le coefficient balistique Bbo utilisé comme abscisse dans cet abaque ne prend pas acte de la Masse Volumique de l’air terrestre ni du coefficient ½ qui le précède généralement dans les travaux d’aérodynamique.

Ce coefficient Bbo n’est donc pas un coefficient qui pourrait être utilisé tel quel dans l’ensemble de l’Univers. Nous y voyons avec satisfaction l’indice que l’armée américaine n’est pas en état d’agir dans le Cosmos lointain 51.


Remarquons cependant que les américains ont opté pour un coefficient qui se rapproche beaucoup de notre Distance Balistique 52 : leur Coefficient Balistique Bbo vaut 1,225 / 2 fois notre Distance Balistique, ceci même s’il n’est pas une distance (il s’exprime en Kg/m2 ).

On peut donc dire que, tout comme notre Distance Balistique, il donne une représentation de la distance qu’il faut à un projectile pour que sa vitesse soit divisée par un certain nombre lors d’une course balistique dans une bulle d’air soustraite à l’action de la pesanteur 53.


 La quantité Gtb utilisée également comme entrée dans l’abaque (elle commande le passage d’une courbe à une autre) est assez proche d’une impulsion. Comme l’impulsion très particulière qu’est l’Impulsion Spécifique d’un ergol, Gtb s’exprime en seconde. Ainsi que nous l’avons dit, Gtb représente le produit par la durée de propulsion tb du quotient de la Poussée sur le poids initial de la fusée : P tb/(Minit g).

Le paramètre G est assez proche lui-même de l’inverse de notre paramètre Rpip, rapport du poids initial sur la Poussée moyenne. Souvenons-nous que cette quantité Rpip naît d’elle-même sous la plume lorsque l’on manipule le Premier Terme de notre Premier Amendement…

Intuitivement, on peut ressentir que, le quotient P/Minit représente l’accélération initiale (à l’action de la gravité près). Le produit de cette quasi accélération initiale par le temps va donc donner accès à la vitesse que la fusée obtiendrait si sa masse restait constante et égale à la Masse initiale.
Néanmoins, dans la mesure où le Rapport de Masses de la fusée considérée est assez proche de l’unité, ce produit P/Minit donne assez bien l’ordre de grandeur de la Vitesse de Fin de Propulsion.
Les auteurs de l’abaque ont opté pour la division de cet ordre de grandeur par la constante de la gravité terrestre (ce qui donne pour G un résultat sans dimension). Ce n’est pas une obligation, mais nous l’avions fait également avec notre paramètre Rpip et cela n’ôte pas à la quantité Gtb la qualité de donner une assez bonne idée de la Vitesse de Fin de Propulsion (au coefficient g près).
Voici la représentation (ci-dessous en bleu) de la quantité Gtb selon l’évolution du Rapport de Masses :

En rouge cette quantité Gtb a été (re)multiplié par g , c'est-à-dire qu’on lui donne la valeur P tb/(Minit) au lieu de P tb/(Minit g)
On doit admettre que, pour les petits Rapport de Masses (inférieurs à 1,2), ce g*Gtb 54 reste assez proche de la vitesse de Tsiolkovski sans influence de la gravité Véject Ln(R) , en fuchsia (la vitesse d’éjection Véject ayant été prise pour 2450 m/s) et donc qu’il peut constituer un paramètre valable dans l’étude des Pertes de Vitesse par Traînée lors de la phase propulsive d’une fusée (8,6 % d’erreur à R = 1,2).

En deuxième analyse, d’ailleurs, la divergence de ce paramètre g*Gtb d’avec la Vitesse de Tsiolkovski sans influence de la gravité peut fort bien constituer un élément permettant à l’abaque étudié de mieux se caler sur la réalité : par exemple on peut constater que la courbe rouge rend mieux compte de l’évolution probable de la vitesse (sous l’influence de la Traînée) que la courbe fuchsia de Tsiolkovski (qui néglige l’influence de la Traînée)…

Mais nous venons de présenter ce graphe de comparaison sans préciser comment nous l’avons obtenu. Expliquons-nous :
Nous sommes partis du libellé d’origine :

Gtb =
…et nous y avons exprimé la Poussée P sous la forme qVéject :
Gtb = =

Cette valeur de la Poussée P = qVéject est souvent diminuée de l’action de la différence de pression (Péject – Patmos) à l’aval de la tuyère d’éjection sur la section de la tuyère en ce point. Nous n’en avons pas tenu compte ici, pour simplifier. 55


Tabler sur la Vitesse d’éjection des fusées à feu est tout à fait raisonnable dans la mesure où cette Vitesse d’éjection est justement une constante pour chaque type de propergols : ainsi, pour les fusées militaires à poudre (dont traite l’abaque), cette Vitesse d’éjection est toujours de l’ordre de 2450 m/s56

Ensuite de quoi, la souvenance que la Masse d’Appui d’une fusée s’exprime également sous la forme :



Q = Minit
…nous a permis de prendre comme valeur de Gtb :
Gtb =
…valeur que nous avons modernisée en :
Gtb g = Véject pour obtenir la courbe rouge.
Cette nouvelle écriture de Gtb g nous enseigne donc que ce paramètre ne dépend que de la Vitesse d’éjection et du Rapport de Masses (et il en va de même pour Gtb)…

Remarquons que cette même courbe rouge admet une asymptote horizontale y = Véject. Elle apparaît lorsqu’on élargit l’échelle des ordonnées :



L’existence de cette asymptote distingue nettement la courbe rouge de la courbe fuchsia et disqualifie sans doute le paramètre Gtb pour ce qui est de l’étude des fusées à forts Rapports de Masses (propres aux fusées spatiales et aux fusées à eau)…

Cependant, la comparaison des deux équations dessinant les courbes fuchsia et rouge, à savoir :



VTsiolSansG = Véject Ln(R)
…et :

Gtb g = Véject
…permet de tirer une valeur assez précise de Gtb libellée à partir de la Vitesse de Tsiolkovski sans action de la gravité et de R. Cette valeur exacte est :
Gtb = VTsiolSansG

…qui est notre interprétation du paramètre Gtb utilisé dans l’abaque du Military Handbook.


Nous l’emploierons plus bas...

Analyse de l’abaque de l’US Army
L’allure générale des courbes nous indique que, lorsque le Coefficient balistique d’un corps devient faible (que le corps est doté d’un Cx plus fort), sa décélération moyenne ΔV/tb durant la phase propulsive sous l’action de la Traînée, est plus forte.

C’est normal car, au premier ordre, la décélération est proportionnelle au Cx (c’est d’ailleurs à ce premier ordre que se cantonne, on le sait, notre premier amendement 57).

Comme le Military Handbook fait figurer le Cx au numérateur du Bbo, une simple proportionnalité du ΔV/tb avec le Cx donnerait une hyperbole sur un abaque orthonormé. Ici, cependant, l’abaque admet des abscisses logarithmiques, ce qui fait qu’il est difficile de dire si les trois courbes sont des hyperboles ou non.

Notre tableur peut facilement rafraîchir notre mémoire sur ce point : dessinée sur un abaque à abscisses logarithmiques, la fonction affine y = ax + b donne la courbe bleue ci-dessous. La fonction y = c/x produit, quant à elle, la courbe rouge, alors que les fonctions y = c/x0,8 et y = c/x1,2 produisent les courbes fuchsia et orange :


À ce stade de l’analyse de l’abaque du Military Handbook, on peut donc dire qu’il est probable que ses trois courbes soient des hyperboles de formule y = c/xn.


La valeur de l’exposant n de ces courbes sera d’ailleurs, à notre sens, symptomatique de la rétroaction du SCx sur les pertes par Traînée : l’accroissement du SCx diminue en effet les vitesses atteintes et, en conséquence, diminue légèrement les pertes de vitesse par Traînée par rapport à leur accroissement strictement lié à l’accroissement du SCx

Ce qui revient à dire que l’accroissement du SCx produit des effets moins que proportionnels sur la décélération moyenne par Traînée (ou sur les pertes de vitesse par Traînée.


Cela apparaît mieux si l’on reprend l’abaque de l’US Army et que l’on exprime les même ordonnées par rapport à l’inverse de la variable Bbo, en disposant ces nouvelles abscisses de façon proportionnelle (et non plus logarithmique), et en les orientant dans le sens des SCx croissant (à Masse Balistique constante). Ce sont les courbes en traits pleins ci-dessous :

On remarque alors plus aisément que la décélération moyenne (sur l’axe des ordonnées) est un peu moins que proportionnelle au SCx (ce qu’indiquent les deux tangentes à l’origine, en trait d’axe) 58  59.

Continuons l’analyse de l’abaque de l’US Army pour relever que l’on passe d’une des courbes à une autre selon la valeur du paramètre Gtb.

Nous avons utilisé notre tableur pour traduire cette famille de courbe en une formule analytique unique qui vaudrait pour la pyro-fuséologie d’amateurs.
Des travaux rapides à l’aide de ce tableur nous ont convaincu que les trois courbes américaines sont homothétiques, assez précisément en tous cas pour les deux plus basses.

Ce sont d’ailleurs ces deux courbes les plus basses qui nous intéressent plus particulièrement puisque les produits Gtb usuels pour les amateurs fuséistes sont de 6,18 " pour le Wapiti Moyenné et de 16,11 " pour l’Isard Moyenné (tous deux nettement en dessous, donc, de la courbe la plus basse qui vaut pour un Gtb de 39,9 ").


Voilà, en couleurs, la démonstration visuelle de cette homothétie.

La saisie de la courbe basse de l’abaque américain produit la courbe rouge ; les courbes fuchsia et bleu dense sont homothétique de la courbe rouge (coefficients multiplicateurs 3,7 et 12) :


On remarque que l’homothétie est quasi parfaite entre les deux courbes basses et encore acceptable entre la courbe basse et la courbe haute.

Le Military Handbook est fort chiche de renseignements quant à l’interpolation entre les trois courbes. Nous notons cependant, pour notre compte, que la croissance du coefficient d’homothétie (1, 3,7 et 12) correspond à une évolution de la puissance 2,54ème du rapport des paramètres Gtb ; ce qui revient à dire que le coefficient d’homothétie permettant de passer d’un courbe à l’autre est [Gtb / 39,9]2,54.
D’autre part, une chance n’arrivant pas toujours seule, Excel nous donne comme régressions possibles des deux courbes basses, les libellés suivants, indiqués sur le graphe :

Nous simplifions ces libellés en y = 40000 x-0,88 et y = 11000 x-0,88 , ce qui donne les marques en cercles sur le graphe suivant :

Le triangle rouge correspond à l’abscisse Bbo et à l’ordonnée ΔV/tb d’une fusée à moteur "Isard moyenné".

Une fusée "Wapiti moyenné" donnerait un point un peu à gauche de l’abscisse 1000).

L’emplacement de ces deux points donne une bonne idée de la zone où se situent les fusées à feu d’amateurs (leur perte de vitesse par seconde étant alors de l’ordre de seulement quelques m/s)…

Il est alors aisé de déduire de ces enseignements que la formule analytique unique représentant de façon satisfaisante l’évolution du quotient ΔV/tb selon les variable Gtb et Bbo de l’abaque de l’US Army est :


ΔV/tb = (Gtb / 39,9)2,54 (11000 Bbo-0,88)

Comparons à présent cette formulation, au moins dans ses ordres de grandeur, à notre Premier Amendement (qui est évidemment l’équivalent de ΔV).


Si nous en restons aux ordres de grandeurs, nous pouvons utiliser, dans un premier temps, le seul Premier Terme de notre Premier Amendement, qui est :
½ ρ SCx (VTsiolSansG)2 {}
…avec toujours :

: R le Rapport de Masses Final à T = instant de fin de propulsion

: Q la Masse d’Appui, le Débit Massique de cette Masse d’Appui étant toujours supposé constant.

: VTsiolSansG la Vitesse Instantanée de Tsiolkovski calculée sans influence de la gravité.
(nous nous sommes affranchi du signe moins qui n’est là que pour signifier que le premier terme de notre amendement prédit des pertes de vitesses et non des gains)

Comme Q peut s’exprimer sous la forme :


Q = MBal(R-1)
… ce Premier Terme de notre Premier Amendement peut encore s’écrire :
ρ SCx (VTsiolSansG)2
…et il représente la prédiction, par le Premier Terme de notre Premier Amendement, des pertes de vitesse par Traînée pendant la phase propulsive.

En faisant apparaître la Masse Balistique, nous avons recomposé l’inverse du Coefficient balistique américain SCx/MBal qui est utilisé dans l’abaque du Military Handbook pour les abscisses.

Intuitivement, on aurait pu penser que seul le SCx intervenait dans les pertes par Traînée en phase propulsive. Mais c’aurait été oublier que plus la Masse de l’engin est forte et moins il se produit de pertes de vitesse par Traînée (puisque une plus grande partie de l’énergie de propulsion est emmagasinée en énergie cinétique dans ladite Masse au lieu d’être transformée en vitesse) 60.
La Masse Balistique présente dans le Coefficient Balistique de l’US Army l’est donc au titre de masse d’inertie. C’est d’ailleurs la même Masse Balistique qui préside aux calculs classiques de trajectoire qui suivent la phase propulsive et qui figure dans notre Distance Balistique, même si cette même masse d’inertie, de résistante est alors devenue motrice 61
Note sur une autre méthode d’accès aux ordres de grandeur :
On pourrait aussi, à titre d’approche des ordres de grandeur, s’affranchir du logarithme en le remplaçant par la régression linéaire Ln(R) = R – 1 (on sait que pour les valeurs proche de 1, ce qui est le cas de nos fusées d’amateurs, la pente du logarithme reste proche de l’unité). Cette substitution donne :
½ ρ SCx (VTsiolSansG)2 {}
La quantité (R – 1) admet alors une autre écriture ; cette écriture est :
Q /(Minit – Q) ou encore : Q /(MBal)
Il en résulte que les pertes de vitesse par Traînée en phase propulsive prédites par notre premier amendement, pour les Rapports de Masses proches de l’unité, sont peu différentes de :
ΔV ≈ ρSCx (VTsiolSansG)2


Comment concilier notre formulation et celle du Military Handbook ?
Cette question nous a occupé l’esprit quelques jours.
Un peu plus haut, nous avons considéré que le produit Gtb donnait un bon ordre de grandeur de la vitesse atteinte en fin de propulsion. Nous avons même calculé notre présentation de ce produit en référence à la Vitesse de Tsiolkovski (sans influence de la gravité) et à :
Gtb = VTsiolSansG

La formule analytique unique décrivant l’abaque du Military Handbook (que nous avons dégagée plus haut) en devient :


ΔV/tb = [ VTsiolSansG / 39,9]2,54 (11000 Bbo-0,88)

Extrayons à présent le SCx du paramètre Bbo (puisque Bbo vaut MBal / SCx) :


ΔV/tb = (SCx)0,88 [ VTsiolSansG / 39,9]2,54 (11000 MBal-0,88)

Il n’a échappé à personne que la durée de propulsion tb américaine équivaut à notre durée de propulsion T.

Les pertes de vitesse par Traînée en phase propulsive prédites par l’abaque états-unien s’élèvent donc à :
ΔV = (SCx)0,88 [ VTsiolSansG / 39,9]2,54 (11000 )

…ou encore, en introduisant la Masse Volumique de l’air :


ΔV = []2,54 (SCx)0,88 VTsiolSansG2,54

Ce qui, toutes simplifications effectuées, donne :


ΔV = 0,00233 []2,54 ρ(SCx)0,88 VTsiolSansG2,54

…qui est une bonne formulation analytique des pertes de vitesse par Traînée en phase propulsive pronostiquées par l’abaque de l’US Army.

Cette formulation offre l’avantage d’utiliser les mêmes paramètres que ceux qui figurent dans notre recherche d’amendement de la Vitesse de Tsiolkovski.
Comparons cette prédiction américaine au seul Premier Terme de notre Premier Amendement encadré plus haut, à savoir :
ρ SCx (VTsiolSansG)2
Ces deux encadrés présentent des similitudes intéressantes, à quelques différences d’exposants près (différences également intéressantes).
Nous avons relevé beaucoup plus haut que ce Premier Terme de notre Premier Amendement ne donnait pas des résultats satisfaisants pour tous les types de moteurs utilisés par les pyrofuséistes amateurs. Ainsi, pour les moteurs type "Wapiti moyenné" (à phase propulsive de 3,4’’, c à d particulièrement longue), la comparaison entre la prédiction du seul Premier Terme de notre Premier Amendement (en fuchsia) et la prédiction de l’US Army (en bleu) est décevante :

(les verticales bleu clair représentent les masses balistiques extrêmes autorisées par le règlement pour ce moteur)
Le seul Premier Terme de notre Premier Amendement pronostique des pertes de vitesses 3 à 4 fois plus fortes que celles revendiquées par l’US Army !
Cette divergence est due, rappelons-le, à la très grande durée de la phase propulsive de ce moteur "Wapiti moyenné" : les pertes de vitesse par gravité au cours de cette phase propulsive est du même ordre que la Vitesse de fin de propulsion 62, or le seul Premier Terme de notre Premier Amendement ne tient pas compte de la diminution de Traînée consécutive à la perte de vitesse par gravité…

Ceci constaté, nous avons construit plus haut dans notre texte un libellé régressé de notre Premier Amendement qui possédait la qualité de toujours rester excessif pour des Rapports de Masses inférieurs à 1,2. Ce libellé est :


½ ρ SCx {1- 1,86 Rpip + Rpip2 }
…le paramètre Rpip n’étant rien d’autre que le rapport Poids initial/Poussée, à savoir : M0g / P (inverse du paramètre G américain)…
L’utilisation de ce libellé toujours excessif en Rpip améliore grandement la prédiction des Pertes de vitesse par Traînée pour le moteur "Wapiti Moyenné" (c’est la courbe en rouge ci-dessus)

Pour être exhaustif, nous avons fait figurer également ci-dessus la courbe dessinée par l’intégralité de notre Premier Amendement (avec tous ses termes, donc) (marques en croix fuchsia), intégralité qui, ainsi que nous l’avons démontré plus haut dans ce texte, possède également la qualité d’être excessive (pour toutes les valeurs de R)…

Notons d’ailleurs que ces croix fuchsia passent sous la courbe bleu américaine, ce qui peut inciter à penser que cette courbe est placée un peu haut par les auteurs états-uniens… 63
Ces croix passent également sous notre courbe rouge que nous avons pareillement donnée pour excessive : cela peut mettre mal à l’aise un instant, mais reste néanmoins logique, dans la mesure où l’une des deux courbes (croix fuchsia et courbe rouge) peut être plus excessive que l’autre…

La divagation désagréable de la courbe fuchsia (dessinant le seul Premier Terme de notre Premier Amendement) par rapport à la courbe américaine bleue est moins nette cependant pour les fusées à moteur "Isard moyenné" :



En effet, ce moteur présente une durée de fonctionnement beaucoup plus brève (1,55’’). Cette brièveté diminue les pertes de vitesse par gravité à 15,2 m/s 64, alors que la Vitesse de Tsiolkovski (sans influence de la gravité) s’élève de 140 à 290 m/s (soit 10 à 20 fois plus) selon la Masse balistique choisie (et acceptée par le règlement).
D’ailleurs la valeur intégrale de notre Premier Amendement (croix fuchsia) fait jeu égal avec la valeur pronostiquée par le libellé en Rpip toujours excessif (courbe rouge ci-dessus) .
Cette sympathique unanimité valide au passage l’usage de ce libellé en Rpip toujours excessif (au moins pour ce moteur et probablement pour ce type de moteurs normalement bref)…

Pour plus de précision, nous allons utiliser maintenant comme libellé des pertes de vitesse par Traînée les libellés spéciaux établis précédemment pour chaque type de moteur.

Par exemple, pour le moteur "Wapiti moyenné", celui-ci, établi plus haut :

ΔV = (0,695 – 0,8 Minit) ½ ρ SCx { Ln3R}

Comme c’est le cas pour l’abaque américain ou notre formule unique qui le décrit, la Masse Balistique doit constituer les abscisses. C’est égal puisque, pour chaque Masse Balistique on a une valeur de la Masse Initiale (qui vaut MBal – Q) puisque la Masse d’Appui Q est une constante du moteur (0,05 Kg). De même, pour chaque Masse Balistique on a une valeur de R, quotient de la Masse Initiale par cette Masse Balistique.


Voici la comparaison avec l’abaque américain (pour un Cx de 0,4) :

Notre prédiction spéciale au Wapiti est en vert 65, la prédiction en Rpip toujours en rouge, celle de l’intégralité de notre Premier Amendement en croix fuchsia, celle du Military Handbook en bleu et celle du Vol de la Fusée en noir à marques rouges 66.

Les deux verticales bleu clair représentent toujours les Masse Balistiques extrêmes admises par le règlement…

La comparaison des courbes semble indiquer que la prédiction des pertes de vitesse par Traînée américaine est un peu surévaluée pour ce type de fusée (qui est d’ailleurs sous motorisée par rapport aux engins militaires dont traite le texte de l’US Army)…


Pour le moteur "Isard Moyenné", l’utilisation de la formule spécialisée dégagée plus haut :



ΔV = 0,84 ½ ρ SCx { Ln3R}
…conduit au graphe suivant :

Notre prédiction propre à l’Isard est en vert, celle utilisant le Rpip toujours en rouge, celle utilisant l’intégralité de notre Premier Amendement en croix fuchsia, celle de l’US Army en bleu et celle du Vol de la Fusée en noir à marques rouges.
Ici, la courbe américaine est plus basse que nos prévisions  "excessives" et c’est normal… Mais elle l’est peut-être trop…

À cet égard, c’est donc la courbe tirée du Vol de la Fusée qui se met en infraction (mais de très peu)…

En conclusion, cette confrontation avec l’abaque de l’US Army semble indiquer qu’avec nos prédictions moteur par moteur comme avec nos prédictions utilisant l’intégralité de notre Premier Amendement ou notre libellé simplifié en Rpip, nous aboutissons à des ordres de grandeurs raisonnables…
La proximité de notre prédiction des Pertes de vitesse par Traînée en phase propulsive avec la prédiction du Vol de la Fusée est de nouveau dévoilée dans ces graphes.


Yüklə 1,76 Mb.

Dostları ilə paylaş:
1   2   3   4   5   6   7   8   9   10   11




Verilənlər bazası müəlliflik hüququ ilə müdafiə olunur ©genderi.org 2024
rəhbərliyinə müraciət

    Ana səhifə