Besoin d'un programme matlab

Fermé
amonatia - 1 juin 2010 à 14:46
Bonjour, j'ai besoin d'un programme matlab qui résoud le systéme d'équations suivant pour avoir les rapports R0/R1, R1/R2, R2/R3

sigma1r(R1) = -P1*(R1.^2/(R1.^2-R0.^2))*[1-(R0/R1)^2]

sigma1o(R1) = -P1*(R1.^2/(R1.^2-R0.^2))*[1+(R0/R1)^2]


sigma2r(R1) = P1*(R1^2/(R2^2-R1^2))*[1-(R2/R1)^2] - P2*(R2^2/(R2^2-R1^2))*[1-(R1/R1)^2]

sigma2o(R1) = P1*(R1^2/(R2^2-R1^2))*[1+(R2/R1)^2] - P2*(R2^2/(R2^2-R1^2))*[1+(R1/R1)^2]


sigma2r(R2) = P1*(R1^2/(R2^2-R1^2))*[1-(R2/R2)^2] - P2*(R2^2/(R2^2-R1^2))*[1-(R1/R2)^2]

sigma2o(R2) = P1*(R1^2/(R2^2-R1^2))*[1+(R2/R2)^2] - P2*(R2^2/(R2^2-R1^2))*[1+(R1/R2)^2]


sigma3r(R2) = P2*(R2^2/(R3^2-R2^2))*[1-(R3/R2)^2]

sigma3o(R2) = P2*(R2^2/(R3^2-R2^2))*[1+(R3/R2)^2]

avec

v1=0.22;
E1=370000;
v2=0.29;
E2=66900;
v3=0.35;
E3=211000;

A = [(1/E2)*(sigma2o(R1) - v2*sigma2r(R1))] - [(1/E1)*(sigma1o(R1) - v1*sigma1r(R1))]

A'= [(1/E3)*(sigma3o(R2) - v3*sigma3r(R2))] - [(1/E2)*(sigma2o(R2) - v2*sigma2r(R2))]


P1=(A' - A*M)/[(M*N - (1-v2)/E2)*(R1.^2/(R2.^2-R1.^2))]

P2=((A' - A*M)*N)/[(M*N - (1-v2)/E2)*(R1^2/(R2.^2-R1.^2))]



N=((1-v2)/E2)*(R1.^2/(R2.^2-R1.^2)) + ((1+v2)/E2)*(R2.^2/(R2.^2-R1.^2)) + ((1-v1)/E1)*(R1.^2/(R1.^2-R0.^2)) + ((1+v1)/E1)*(R0.^2/(R1.^2-R0.^2))

M=(E2*(1-v3)/(2*E3))*((R2.^2-R1.^2)/(R3.^2-R1.^2)) + (E2*(1+v3)/(2*E3))*((R3.^2*(R2.^2-R1.^2))/(R2.^2*(R3.^2-R1.^2))) + (1-v2)/2 + (R1.^2/R2.^2)*((1+v2)/2)