Wie kann ich die Pariser Law Crack -Wachstumssimulation in Python umsetzen?
Posted: 19 Aug 2025, 15:07
Ich versuche, das Wachstum des Ermüdungsrisses im Python unter Verwendung von Pariser Gesetz zu modellieren:
da/dn = c (Delta k)^m
wobei:
a = Risslänge
LaD Cycls
(Stress -Nr. Bereich) (Deltak) = y. (Delta Sigma). (pi.a)^1/2
c, M = Materialkonstanten
Für einen einfachen Fall, i setze: c = 1e-12
m = 3
(Geometriefaktor) y = 1,12 < /p>
Anfangsrisslänge (0,001) m < /p>
Hier ist mein erster Versuch in Python: < /p>
>
Ist dieser iterative Ansatz (Euler vorwärts) numerisch stabil für Risswachstumsprobleme, oder sollte ich einen ODE -Solver wie scipy verwenden. )?>
da/dn = c (Delta k)^m
wobei:
a = Risslänge
LaD Cycls
(Stress -Nr. Bereich) (Deltak) = y. (Delta Sigma). (pi.a)^1/2
c, M = Materialkonstanten
Für einen einfachen Fall, i setze: c = 1e-12
m = 3
(Geometriefaktor) y = 1,12 < /p>
Anfangsrisslänge (0,001) m < /p>
Hier ist mein erster Versuch in Python: < /p>
>
Ist dieser iterative Ansatz (Euler vorwärts) numerisch stabil für Risswachstumsprobleme, oder sollte ich einen ODE -Solver wie scipy verwenden. )?>