FR 2627302 A1 19890818 FR 8801857 A 19880217 PROCEDE ET DISPOSITIF DE CALCUL DU PAS D'INTEGRATION DE TRAJECTOIRE D'OBUS La présente invention concerne un procédé et un dispositif de calcul du pas d'intégration de trajectoire d'obus. Dans le domaine de la balistique, on sait calculer les éléments de tir (hausse, azimut, event) à appliquer à un canon, relativement à une configuration de tir donnée (position canon-cible, données géographiques et météorologiques, carac- téristiques de l'obus) pour que l'obus ainsi tiré atteigne la' cible. I0 On effectue pour cela plusieurs calculs d'intégration de trajectoire, avec corrections successives à partir de conditions de tir choisies initialement, jusqu'à ce que le point de chute de l'obus soit situé dans des limites définies par rapport à la cible. La méthode d'intégration elle-même consiste clas- siquement en un algorithme d'intégration quelconque (prédicteur- correcteur, Runge-Kutta, Taylor...) avec un pas d'intégration variable en fonction du point considéré de la trajectoire, donné + par h = k. I, o k désigne une constante, et v et y les vecteurs vitesse et accélération de l'obus. Cette méthode, bien que précise, a cependant pour inconvénient de manquer de rapidité car le pas d'intégration ainsi calculé n'est pas réellement optimisé en fonction du point 2 considéré de la trajectoire, en ce sens qu'il suffirait parfois 2 d'utillser un pas d'intégration de valeur supérieure sans que la précision de calcul en soit affectée. La présente invention a pour objet une méthode de calcul d'un tel pas d'intégration optimisé, c'est-à-dire du pas d'intégration maximal compatible avec la précision -demandée, permettant de gagner sur le temps d'intégration des trajectoires sans perdre sur la précision du calcul d'intégration. Le procédé de calcul du pas d'intégration de tra- jectoires d'obus suivant l'invention consiste essentiellement - dans une étape préalable: à calculer, pour différentes configurations de tir, une trajectoire dite de référence, avec un pas d'intégration donné, relativement faible, dit pas de référence, à essayer successivement, en différents points des différentes trajectoires de référence ainsi obtenues, diffé- rentes valeurs du pas d'intégration, croissant à partir du pas de référence jusqu'à ce que l'écart avec la trajectoire de référence soit supérieur à une limite prédéterminée, dite erreur de position tolérable par pas d'intégration, auquel cas le pas d'intégration retenu, dit pas configuration de tir considérée. D'autres objets et caractéristiques de la présente invention apparaîtront plus clairement à la lecture de la description suivante d'un exemple de réalisation, faite en relation avec les dessins ci-annexés dans lesquels - la figure 1 est un diagramme destiné à illustrer le procédé de recherche du pas d'intégration limite; - les figures 2a et 2b sont des diagrammes illustrant la modé- lisation du pas d'intégration limite - la figure 3 est un organigramme qui recense les opérations à effectuer dans les moyens de caleul associés à la batterie de tir, en vue de calculer le pas d'intégration limite. Z627302 On expose tout d'abord, en relation avec la figure 1, le principe de la méthode qui permet d'optimiser le pas d'intégration. Soit un algorithme d'intégration quelconque (prédicteur-correcteur, Runge-Kutta, Taylor...) destiné à estimer la trajectoire. A chaque pas, l'état cinématique (i+1) est construit à partir - de l'état cinématique (i) (position, vitesse aérodynamique, vitesse de rotation, obliquité vectorielle); - de données météorologiques locales (température, pression, vitesse du vent) - de données géographiques (latitude du lieu de tir) - de données caractéristiques de l'obus (masse m...). L'ensemble de ces variables constitue un point P de la trajectoire. A un point P quelconque, on peut associer un pas maximal d'intégration, appelé temps limite, noté hlimite, et qui est déterminé de la façon suivante. A partir de P, on génère - une trajectoire de référence Tréf avec un pas H faible (H = 0.05 seconde par exemple) - plusieurs trajectoires avec un pas que l'on fait croître tant que l'écart E avec Tréf reste inférieur à une valeur p dite erreur de position tolérable par pas, dont la détermination sera précisée ultérieurement. Sur la figure t on a ainsi représenté à titre d'exemple trois trajectoires: trajectoire 1 (essai de h = 2H): E ú On obtient' donc dans ce cas, pour le point P considéré: hlimte = 3H. On constitue ainsi en répétant ce calcul pour divers points P une base de données constituée de couples (hlimite, X), le terme X recouvrant les variables à considérer 3 en chaque point P, à savoir état cinématique en ce point, données géographiques et météorologiques en ce point, caractéristiques de l'obus. Ce calcul est ensuite répété pour différentes confi- gurations de tir T, chacune d'elles étant caractérisée par les paramètres suivants: conditions de tir appliquées au canon, données géographiques et météorologiques, caractéristiques de l'obus (masse, vitesse d'éjection). On constitue ainsi une base de données (T, X, E hlm) him étant caleulé comme indiqué précédemment. Cette méthode est relativement. coûteuse en temps de calcul, mais elle n'a pas à être mise en oeuvre sur le terrain, dans les moyens de caleul associés à la batterie de tir; elle sera avantageusement mise en oeuvre préalablement, avec des moyens de caleul plus performants que ceux associés à la batterie de tir, dont les performances sont évidemment limitées pour des raisons d'encombrement. Cette base de données ne pouvant cependant être exhaustive, il est nécessaire, si l'on veut effectivement diminuer le temps de caleul des trajectoires sur le terrain, à l'aide des moyens de calcul associés à la batterie de tir, de déterminer un moyen rapide de calcul de la valeur hi correspondant à un point P quelconque, à partir de cette base de données et en fonction des paramètres existant au point considéré et de la configuration de tir considérée. Pour cela, une fois la base de données constituée, on procède à une modélisation de huim, sachant que le modèle recherché doit satisfaire au mieux aux contraintes de précision, de temps de caleul, et de place mémoire dans les moyens de calcul associés à la batterie de tir. On procède pour cela dans un premier temps à une recherche des variables dites explicatives, c'est-à-dire des variables en fonction desquelles la variable à modéliser hum s'exprime suivant une certaine fonction. Ceci revient géométriquement à réaliser une projection de l'ensemble des points représentant hlim dans l'espace des différentes variables "x" de l'ensemble (T, X, Ep), sur chacun p des plans (him, x), puis à sélectionner celles de ces variables "x" ayant le plus fort pouvoir explicatif. Pour tenir compte des objectifs de rapidité recherchés, on exprime hi linéairement, par domaines des lim variables explicatives. Après mise en oeuvre de cette méthode, on est ainsi amené à constater que dans une zone "Mach élevé", soit pour les valeurs de M supérieures à une certaine valeur Mo, hlm s'exprime linéairement, par domaine, en fonction de M 2, yp Pr, m, Te et Uo, soit une expression de la forme hit ai M- 2 =ï o I+a2 +al3P hlim= aO +ailM 2+ ai2 ep + ai3 Pr + a4 m +a5 Te+a6 U o M désigne le nombre de Mach (rapport de la vitesse de l'obus à la vitesse du son),- Pr la pression locale, m la masse de l'obus, Te la température et U la vitesse d'éjection de eo l'obus. A l'intérieur de la zone "Mach élevé" (M > Mo) on est encore amené à distinguer quatre zones délimitées par des -2 valeurs M1, M2 et M3 de M 2; à chacune de ces zones correspond respectivement un ensemble de coefficients a i i i ai5 i ai a2, a3, a4, 5 eta6 (1 4 i utilisés dans l'expression du modèle ci-dessus. Dans une zone "Mach faible", soit pour les valeurs de M inférieures à Mo, hm s'exprime linéairement, par domaine, en fonction de p. Ep, M, Te et WX o WX désigne la composante du vent suivant la direction horizontale, et p la pente de la trajectoire soit une expression de la forme hlim = bio + b p + bi2 ep * bi3 M hib0+b1p+ 2 p 3 + bi4 Te + bi5 WX A l'intérieur de la zone "Mach faible" (M inférieur à une valeur M4. A chacune des huit zones ainsi définies correspond un ensemble de coefficients bil1 b2 b3 b4 2bi 3, bi4, b15 (1 4 i sera décrite ultérieurement. Sur la figure 2b on a représenté graphiquement le modèle permettant d'approximer hlim dans un espace limité pour plus de clarté aux deux variables explicatives cE et M2 et dans cet espace on a considéré seulement le cas M > M (zone "Mach élevé"). Ce modèle est constitué par une o succession de plans qui se projettent chacun, comme représenté sur la figure 2a, dans chacun des plans (hlim Ep) et (hium, M-2) suivant les segments de droite résultant de la modélisation linéaire. Dans le cas d'un espace non limité à deux variables explicatives, la notion de plan serait remplacée par celle d'hyperplan. Les modèles étant ainsi déterminés, le caleul du pas optimisé d'intégration des trajectoires, effectué sur le terrain par les moyens de calcul associés à la batterie de tir se fait de la façon représentée sur l'organigramme de la figure 3. On remarquera que le nombre d'opérations rajoutées par le calcul du modèle est négligeable devant le nombre total d'opérations nécessitées par une phase d'intégration (qui demande deux calculs de la dérivée du vecteur d'état). Les opérations de test de positionnement et de calcul de modèle correspondent à celles décrites précédemment et ne seront donc pas redécrites. Les valeurs des paramètres qu'il est nécessaire de connaître au point considéré sont les valeurs des paramètres M, UX, Uy, U Z Ao, Te' Pr' U, o WX' m. Les paramètres UX, Uy, Uz désignent les composantes de la vitesse absolue de l'obus au' point considéré; A désigne la hausse du canon. Les valeurs des paramètres UX, Uy U sont nécessaires au calcul du nombre de Mach M (et de M) et de la pente p. La valeur de la hausse A est nécessaire au calcul o de l'erreur de position tolérable par pas ú p au point considéré. Le calcul de E est en effet effectué lui aussi en p ayant recours à un modèle préétabli, et lors de l'établissement de ce modèle, on constate que cette variable S peut être p modélisée en fonction de Ao, suivant trois segments de droite de la forme: = Ci +ci A p o 1 o chacun des trois segments de droite étant défini par un ensemble de coefficients ci cI (1 = 1, 3) A prenant ses O' i ' valeurs sur trois domaines [A A1 A2] [A2 A3] et [A3, A4]. Cette modélisation de Ep est obtenue à partir d'une base de données (, T) qui peut elle-même être obtenue de la façon suivante. Pour une configuration de tir T donnée et pour chaque point P, on essaye différentes valeurs de cp en appliquant à chaque fois la méthode "du temps limite" P telle que décrite plus haut en relation avec la figure 1, et ceci successivement pour tous les points de la portion de trajectoire reliant le point P considéré au point de chute de l'obus. On retient alors la valeur maximum compatible avec la précision demandée sur le point de chute de l'obus, et on renouvelle cette méthode pour différentes configurations de tir. Le procédé suivant l'invention permet de réduire notablement les temps de caleul et de se passer d'opérateurs spécialisés dans les moyens de calcul associés à la batterie de tir; l'architecture de ceux-ci s'en trouve simplifiée et leur consommation d'énergie diminuée, ce qui est fondamental pour des problèmes de portabilité et d'autonomie. A titre d'exemple pour un tir à distance moyenne (15 à.20 km) alors qu'un calcul de trajectoire traditionnel nécessite soixante pas d'intégration environ le procédé suivant l'invention permet de ramener ce nombre à une vingtaine. REVENDICATIONS 1. Dispositif de calcul de trajectoire d'obus en vue de la commande d'un canon, caractérisé en ce qu'il comporte, pour le calcul du pas d'intégration de ladite trajectoire: des moyens de calcul, pour différentes configurations de tir, d'une trajectoire dite de référence, avec un pas d'intégration donné,relativement faible, dit pas de référence; des moyens de caleul, en différents points des différentes trajectoires de référence ainsi obtenues, de l'écart entre la trajectoire de référence et différentes trajectoires obtenues en faisant croître le pas d'intégration, à partir du pas de référence, jusqu'à ce que cet écart soit supérieur à une limite prédéterminée, dite erreur de position tolérable par pas d'intégration, auquel cas le pas d'intégration retenu, dit pas d'intégration limite, est celui correspondant à la dernière trajectoire obtenue, des moyens pour effectuer une modélisation du pas d'intégration limite à partir des données ainsi obtenues, en fonction de différents paramètres définissant les points des trajectoires et les configurations de tir. 2. Dispositif selon la revendication 1, caractérisé en ce qu'il comporte des moyens de calcul du pas d'intégration limite à l'aide du modèle ainsi établi, en fonction des paramètres existant au point de la trajectoire o ce calcul est à effectuer, et de la configuration de tir considérée. 3. Dispositif selon la revendication 1, caractérisé en ce qu'il comporte en outre: des moyens de calcul du pas d'intégration limite pour différentes valeurs de l'erreur de position, et ce pour toute la portion de trajectoire reliant le point considéré au point de chute de l'obus des moyens pour déterminer, parmi ces différentes valeurs de l'erreur de position, la valeur maximum compatible avec la précision tolérée sur le point de chute de l'obus, -des moyens pour effectuer une modélisation de cette valeur maximum, dite erreur de position tolérable par pas d'inté- gration, à partir des données ainsi obtenues. 4. Dispositif selon la revendication 2, caractérisé en ce qu'il comporte, pour calculer le pas d'intégration limite - des moyens pour calculer, au point o ce caleul de la trajectoire est à effectuer, le nombre de Mach M de l'obus, la -2 valeur M, et la pente p de la trajectoire, -des moyens pour positionner ces trois valeurs par rapport à des valeurs particulières déterminées par modélisation, respectivement Me, M1, M2, M3, M4 P1, P2, P3s telles que: h ai + ai M 2 +ai + ai P Bim o a 2 p 3 r + a4 m a5 + a6 U pour des valeurs élevées du nombre de Mach de l'obus, soit M > M et: hlim = b + b Pp bi2 + bi3 M li1p+b2 p, 3 + b 4 Te + bi5 WX pour des valeurs faibles du nombre de Mach de l'obus, soit M Mo, ou suivant la position de la pente p de la trajectoire au point considéré et du nombre de Mach M de l'obus, par rapport, respectivement, aux valeurs P1, P2, p3 et M4 dans le cas M' du pas d'intégration limite au point considéré. 5. Dispositif selon les revendications 2 et 3, caractérisé en ce qu'il comporte en outre des moyens pour calculer l'erreur de position E tolérable au point considéré, p suivant le modèle: C = ci +ci A p o o les termes c. désignant un ensemble de coefficients de valeurs différentes suivant la position de A (A0 désignant la hausse du canon) par rapport à des valeurs particulières A. déterminées par modélisation de l'erreur de position tolérable par pas d'intégration.