Bonjour a tous,

Dans le cadre de mon TIPE, je cherche a modéliser la trajectoire d'un avion en  chute libre, en calculant d'abord la vitesse en tous points. Pour cela, j'ai utilisé le module odeint de scipy, mais ce dernier me renvoi

"Excess work done on this call (perhaps wrong Dfun type).

Run with full_output = 1 to get quantitative information."

J'obtient quand même des valeurs mais qui deviennent très vite délirantes...

 Voila le programme :

####################

#Imports

####################

import matplotlib.pyplot as plt

import numpy as np

from scipy.integrate import odeint

m=16576

rho=1.225

S=75.21

Cx=0.0305

Cz=0.45

g=9.8

vx_0=220

vz_0=0

x_0=0

z_0=1000

####################

#Function

####################

def phiVitesse (A,t):

""" entre la vitesse en x et en z dans A et donne (x'',z'')"""

dxdt=A[0]

dzdt=A[1]

d2xdt2=(-1/(2*m))*rho*S*Cx*(dxdt**2+dzdt**2)

d2zdt2=-g+(1/(2*m))*rho*S*Cz*(dxdt**2+dzdt**2)

return [d2xdt2,d2zdt2]

####################

#Main

####################

#vitesse initialle

v0=[vx_0,vz_0]

#liste de points

t=np.linspace(0,100,1000)

#solution

A=odeint(phiVitesse,v0,t)

#séparation de dx et dz dans A

dxdt=A[:,0]

dzdt=A[:,1]

plt.clf()

plt.plot(t,dxdt,'red')

plt.plot(t,dzdt,"black")

plt.show()

Avez vous une idée de ce qui cloche?

Merci d'avance !

-
Edité par FrançoisFournier2 12 août 2018 à 16:15:37

Bonjour,

Le message qui suit est une réponse automatique activée par un membre de l'équipe. Les réponses automatiques leur permettent d'éviter d'avoir à répéter de nombreuses fois la même chose, ce qui leur fait gagner du temps et leur permet de s'occuper des sujets qui méritent plus d'attention.
Nous sommes néanmoins ouverts et si vous avez une question ou une remarque, n'hésitez pas à contacter la personne en question par Message Privé.
Pour plus d'informations, nous vous invitons à lire les règles générales du forum

Merci d'utiliser la mise en forme de code

Les forums d'Openclassrooms disposent d'une fonctionnalité permettant de colorer et mettre en forme les codes source afin de les rendre plus lisibles et faciles à manipuler par les intervenants. Pour cela, il faut utiliser le bouton de l'éditeur, choisir un des langages proposés et coller votre code dans la zone prévue. Si vous utilisez l'éditeur de messages en mode Markdown, il faut utiliser les balises <pre class="brush: python;">Votre code ici</pre>.