Wie kann ich die Pariser Law Crack -Wachstumssimulation in Python umsetzen? ⇐ Python

[Anonymous]

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. )?>