Hallo liebe Wissenden,
ich hab das problem, dass ich mich nicht wirklich mit MuPAD auskenne.
Ich möchte gern folgendes Integral numerisch berechnen:
(AS*( (BS+ln((X1/Rth)+1))^Cs)*X1*
( sqrt((Ra^2)-(X1^2))-sqrt((Ri^2)-(X1^2)) ))
wobei as bs rth cs ra und ri bekannt sind.
Sprich das integral läuft von 0 bis ri übder dy
Kann mir jemand sagen, wie man das mit MuPAD oder auch mit einem anderen Programm numerisch lösen kann !?
Ich bin schon am verzweifeln.
Hallo Schnibbi,
nachdem der Ausdruck
(AS*( (BS+ln((X1/Rth)+1))^Cs)*X1*
( sqrt((Ra^2)-(X1^2))-sqrt((Ri^2)-(X1^2)) ))
nicht von y abhängt, sollte für ein Integral
übder dy
kein Computer vonnöten sein, weder für eine analytische Lösung noch für eine numerisch.
Viele Grüße
Stefan
Sorry, hab ich mich vertan, natürlich über d X1 wäre ja sinnlos 
Is mir beim abschreiben wohl etwas verrückt.
Hallo,
mit MuPAD kenne ich mich zwar auch nicht aus, könnte aber versuchen, das Integral mit meinem selbstgeschriebenen VB-Programm auszurechnen. Dazu bräuchte ich allerdings die Werte für AS, BS, Rth, Cs, Ra und Ri.
Gruß Walter
Kriegst du:
Ze=0,0173;
Ra=5,5;
Ri=4,5;
AS=500
BS=0,0075
CS=0,2
Also wenn das klappt, musst du mir unbedinggt mehr über dein VB Programm erzählen.
Sorry Rth vergessen: 15,5
[Bei dieser Antwort wurde das Vollzitat nachträglich automatisiert entfernt]
Hallo,
ich habe mit meinem Programm das Integral
(500*((0.0075+log((x/15.5)+1))^0.2)*x*(sqr((5.5^2)-(x^2))-sqr((4.5^2)-(x^2))))
problemlos ausrechnen können und bin zu dem folgenden Ergebnis gelangt:
5.203,166162 E² für das Intervall von 0 bis 4,5.
Mit dem Wert Ze = 0,0173 habe ich jedoch nichts anfangen können, denn Ze sehe ich in Deiner Ausgangsgleichung nicht.
Wie Du siehst, habe ich in der Gleichung kleinere Anpassungen vornehmen müssen:
Dein X1 habe ich durch x ersetzt, weil ich in meinem Programm nur das kleine x als unabhängige Variable verwenden kann, und sqrt habe ich durch sqr ersetzt, weil VB nur sqr für die Quadratwurzel kennt.
Zu meinem Programm:
Für die numerische Integration wird die Fläche unter einer Kurve - wie Du sicherlich weißt - in beliebig kleine Teilflächen unterteilt. Die Gesamtfläche ergibt sich dann durch die Addition der Fläche dieser Teilflächen. Dabei gilt, daß das Ergebnis um so genauer wird, je kleiner man die Schrittweite (delta-x) wählt.
Für mein Programm, das nach dem Trapezverfahren arbeitet, habe ich eine Schrittweite von 1/1000 voreingestellt und habe damit auch Dein Integral ausgerechnet.
Es würde mich freuen, wenn ich Dir helfen konnte.
Gruß Walter
[Bei dieser Antwort wurde das Vollzitat nachträglich automatisiert entfernt]
Hallo zusammen
Also wenn das klappt, musst du mir unbedinggt mehr über dein
VB Programm erzählen.
Wenn Du Excel hast dann überprüfe es selber 
Function Simpson(a, b As Double, n As Integer, funktion As String)
'(C)oded 1998 by hw
Dim h As Double, t As Double, sum1 As Double, sum2 As Double
Dim i As Integer, m As Integer
ReDim fx(n) As Double
h = (b - a) / n
For i = 0 To n: fx(i) = Eval(a + i \* h, funktion): Next
If (n Mod 2 0) And (n \> 1) Then
t = fx(n - 3) + 3 \* (fx(n - 2) + fx(n - 1)) + fx(n)
Simpson = 3 \* h \* t / 8
m = n - 3
Else
Simpson = 0
m = n
End If
If m \> 1 Then
sum1 = 0: sum2 = 0: i = 1
While i
Und in MuPad
int(f,x=0..4.5)
0.2
int(500 x (ln(0.06451612903 x + 1) + 0.0075)
2 0.5 2 0.5
((30.25 - x ) - (20.25 - x ) ), x = 0..4.5)
float(last(1))
5203.13231
Gruß HW
Ok, ich merk schon ihr seid echte Profis. Muss zugeben in numerik bzw. coding bin ich nun wirklich nicht fit.
Also bei der ganzen Geschichte geht es darum das maximale Biegemoment bei einem Biegevorgang zu berechnen.
Allerdings scheinen mir da mehr als 5000 Nm recht viel…
Weiterhin soll es möglich sein eine maximale Randspannung auszurechnen, allerdings über eine andere Formel.
Vielleicht kann mir bei dieser Geschichte jemand weiterhelfen.
meine Emailadresse wäre [email protected]
Vielen Dank aber schoneinmal für eure Bemühungen
Im link also nochmal die Formel die Berechnet werden soll
http://www.funpic.de/fotoalbum/foto,231381,0.htm
Und daraus ergibt sich angeblich folgende Formel:
http://www.funpic.de/fotoalbum/foto,231780,0.htm
allerdings sind meine werte immer zu hoch für diese Randspannung.