Verfasst am: 24.09.2013, 07:37
Titel: Pressungsverteilung mit surf plotten
Guten Morgen zusammen,
ich habe eine Flächenpressung im Kontakt zweier Körper, deren Verteilung würde ich gerne mit surf plotten.
Leider habe ich anscheinend irgendwo einen Denkfehler, sodass mein plot nicht funktioniert.
Ich habe die Länge der Kontaktfläche l_k, die Breite der Kontaktfläche b und die Pressungsverteilung (durch eine Formel).
Die Pressung nimmt dabei vom Rand b=-b_max (p=0) bis zur Mitte zu und hat sein maximum bei b=0 und nimmt zum Rand b=b_max wieder ab.
Mein plot sieht jedoch anders (komisch) aus (siehe Anhang). Kann mir jemand helfen ?
ohne ein entsprechendes code beispiel werden dir nur leute mit hellseherischen fähigkeiten helfen können da die aber alle hände voll bei 9-irgendwas zu tun haben könntest du das beispiel ja nachreichen
für neulinge bietet sich der tread richtig fragen aus meiner signatur an
grüße winkow
_________________
k=0;
b_press = zeros(1,1);
for i = 1:50
k=0;
for j = -4:1:4
k = k +1;
if b(i,1) == 0
b_press(i,k) = 0;
else
b_press(i,k) = j/4 * b(i,1);
end end end
l_press = zeros(1,1);
for i = 1:50 for j = 1:1:9
l_press(i,j) = (j-1)* l_k1/8;
end end
p_m = zeros(1,1);
for i = 1:50
k=0;
for j = 0:1:4
k=k+1;
if Fn_p(i,1) == 0
p_m(i,k) = 0;
else
p_m(i,k) = 2/(pi()*l_Si_k1*b(i,1)) * Fn_p(i,1) *(j/4)*pi()/4;
end end end
for i = 1:50
k=5;
for j = 3:-1:0
k=k+1;
if Fn_p(i,1) == 0
p_m(i,k) = 0;
else
p_m(i,k) = 2/(pi()*l_Si_k1*b(i,1)) * Fn_p(i,1) *(j/4)*pi()/4;
end end end surf(b_press,l_press,p_m)
Ich schreibe es später so um, dass es lauffähig ist.
Aber schonmal danke für deine Hilfe.
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 25.09.2013, 07:33
Titel:
Ich habe die benötigten Variablen jetzt einfach so eingefügt, dass das Programm läuft.
Leider ist es jetzt nicht nachzuvollziehen, wo diese Werte herkommen, aber der Code, der diese berechnet wäre zu groß um jetzt hier zu posten
Code:
b = [0.0952944844944899;0.204933363122546;0.223402573281324;
0.229299485871436;0.231710198022262;0.236224283244833;
0.230464791501484;0.218596602861475;0.216162571551865;
0.204933363122546;];
k=0;
b_press = zeros(1,1);
for i = 1:N_Scheiben_k1
k=0;
for j = -4:1:4
k = k +1;
if b(i,1) == 0
b_press(i,k) = 0;
else
b_press(i,k) = j/4 * b(i,1);
end end end
l_press = zeros(1,1);
for i = 1:N_Scheiben_k1
for j = 1:1:9
l_press(i,j) = (j-1)* l_k1/8;
end end
p_m = zeros(1,1);
for i = 1:N_Scheiben_k1
k=0;
for j = 0:1:4
k=k+1;
if Fn_p(i,1) == 0
p_m(i,k) = 0;
else
p_m(i,k) = 2/(pi()*l_Si_k1*b(i,1)) * Fn_p(i,1) *(j/4)*pi()/4;
end end end
for i = 1:N_Scheiben_k1
k=5;
for j = 3:-1:0
k=k+1;
if Fn_p(i,1) == 0
p_m(i,k) = 0;
else
p_m(i,k) = 2/(pi()*l_Si_k1*b(i,1)) * Fn_p(i,1) *(j/4)*pi()/4;
end end end
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
da die aber alle hände voll bei 9-irgendwas zu tun haben könntest du das beispiel ja nachreichen 
RSS