hexagon logo

Subroutine for Paris Law

I need to change fatigue behaviour/crack propagation based on grain size. Is there any subroutine on parameter equation for PARIS Law?