Apprendre Maple

[alex]

bonjour je suis en 1 ère année et je fais un exposé sur le trébuchet, malgrès une bonne dizaine d'heure passé devant mon ordinateur à scruter toutes les pages du net pour trouver ma réponse en vain. Je fais appel à vous j'espere que vous pourrez m'aider!

Donc voici mon system d'équadiff: je cherche th(t) et phi(t)


eq1:={-0.30*0.19*45*(phi2+2*th2)*Cos(phi(t)) + 0.30*0.19*45*phi2*(phi2+2*th1)*Sin(phi(t)) + 1/6*(6*0.19^2*45*phi2 + 6*0.30^2*45*th2 + 6*0.19^2*45*th2 + 6*1.21^2*0.45*th2 + 2*0.19^2*1*th2 - 2*0.30*1.21*4*th2 + 2*1.21^2*1*th2 + 6*9.8*0.30*45*Sin(th(t)) - 6*9.8*1.21*0.45*Sin(th(t)) + 3*9.8*0.30*1*Sin(th(t)) -
> 3*9.8*1.21*1*Sin(th(t)) - 6*9.8*0.19*45*Sin(phi(t) + th(t)))=0, 0.19*45*(0.19*phi2 + 0.19*th2 - 0.30*th2*Cos(phi(t)) - 0.30*(th1^2)*Sin(phi(t)) -
> 9.8*Sin(phi(t) + th(t))) = 0}: eq1;

avec phi2=phi'' phi1=phi' th2...
et phi(0)=/4 et th(0)=3/4 et th'(0)=phi(0)=0

est-il possible de le résoudre avec maple juste de façon numérique?
Mes essais n'ont pas été concluant car les fonction dsolve ou fsolve ne sont que pour un systeme d'équations différentielles linéaires ordinaires(ODE).

Auriez vous une solution ?

[zozo]

Bonsoir,
déjà commence par remplacer Sin par sin, Cos par cos, th2 par diff(th(t),t$2), etc...
et phi(0)=/4??? phi(0)=0??? (il faudrait peut-être se mettre d'accord)
A+

[prof_simplet]

[alex]

ça je l'avez déja fait en faite c'était pour simplifier l'équation



eq1:={-2.5650*((diff(phi(t), `$`(t, 2)))+2*(diff(th(t), `$`(t, 2))))*Cos(phi(t))+2.5650*(diff(phi(t), `$`(t, 2)))*((diff(phi(t), `$`(t, 2)))+2*(diff(th(t), t)))*Sin(phi(t))+1.624500000*(diff(phi(t), `$`(t, 2)))+6.349411667*(diff(th(t), `$`(t, 2)))+122.5049000*Sin(th(t))-83.79000000*Sin(phi(t)+th(t)) = 0,
1.6245*(diff(phi(t), `$`(t, 2)))+1.6245*(diff(th(t), `$`(t, 2)))-2.5650*(diff(th(t), `$`(t, 2)))*Cos(phi(t))-2.5650*(diff(th(t), t))^2*Sin(phi(t))-83.790*Sin(phi(t)+th(t)) = 0}

après je fais:

vars:={phi(t),th(t)}:vars; et puis

c_init:=(phi(0)=3*3.14/4,th(0)=3.14/4,th1(0)=0,phi1(0)=0) ; et puis

sol:=dsolve({eq1 union phi(0)=3*3.14/4,th(0)=3.14/4,th1(0)=0,phi1(0)=0}, numeric ,vars);

et la un message d'erreur me dit:

Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {th(t), phi(t), t, phi(t)+th(t)}


je n'arrive pas à interpréter ce message

[zozo]

Ca ne fonctionnera pas si tu laisses th1(0)=0,phi1(0)=0 dans les conditions initiales.

[alex]

en faite à l'origine j'ai posé:

phi2:=diff(diff(phi(t),t),t);

th2:= diff(diff(th(t),t),t);

phi1:= diff(phi(t),t);

th1:=diff(th(t),t);

et il me prend bien mes conditions initiales


c_init:=(phi(0)=3*3.14/4,th(0)=3.14/4,th1(0)=0,phi1(0)=0)

c_init := phi(0) = 2.355000000, th(0) = .7850000000, diff(th(t), t)(0) = 0, diff(phi(t), t)(0) = 0

[alex]

mais si ce n'est pas solvable par maple dite le moi merci!

[ALS]

Bonjour, Maple ne donne rien pour résultat (patientez un peu pour la dernière instruction).

[alex]

et de façon numérique vous pensez que ça pourrait marcher?

[ALS]

Ca peut valoir le coup d'essayer...