Discussion :

[Tiazeet5]

Bonjour,

Je souhaite votre aide car je bloque sur un programme qui parait simple au premier coup d'oeil.
Dans le contexte de la mécanique des fluides je souhaite tracer uniquement la courbe des pertes de charges en écoulement turbulent via la relation de Prandl-Nikuradse. Je rentre donc le code suivant pour avoir mon lambda = f(re).

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
 
%Expression de la perte de charge via Prandl-Nikuradse
%1./sqrt(lambda)=2.035*log(re *sqrt(lambda)-0,8);
 
pas=500;
re=8000:pas:200000;
exp(lambda)/lambda = 2*re.^2./6,3635; % s'écrit aussi de cette façon
 
h = plot(re,lambda);
grid
xlabel 'Re'
ylabel 'Coef. de perte de charge'
title 'Perte de charge via Prandl-Nikuradse''

Le problème est que peu importe la relation que j'essai de tracer il m'affiche l'erreur suivante :

The expression to the left of the equals sign is not a valid target for an assignment.

Dois-je retravailler plus que cela la relation en remplaçant un des lambda par une constante (ce que je trouve pas folichon...) ou matlab est-il capable de solutionner cela ?

Merci d'avance pour vos réponse.

[le fab]

Salut

une idée :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
fun = @(lambda,re) exp(lambda)/lambda - 2*re.^2./6.3635;
re = 8000;
lambda = fsolve(@(x) fun(x,re),42);

cela donne lambda = 19.8028
cette valeur est éventuellement dépendante de la valeur initiale que j'ai donner au solver (42 ici, de manière totalement pifometrique)

Fabien