Generate individual defect time for particle fracture
getCrackTime(theta, a, b, stress, vickers, p, const)
polar angle
axis length (axis orthogonal to rotational axis)
rotational axis length
stress level
Vickers hardness
simulation parameter set
constant parameters for crack
numeric, the individual fracture time
The particle fracture (crack) is assumed to happen orthogonal to the major axis direction along
the maximum minor axis length. Thus the projection area can be easily computed for the purpose
of defect accumulation. The parameter set is made up of six parameters. Here only the second and
third parameters are used to simulate the defect crack
times. The failure times follow a
Weibull distribution with scale parameter \(p2*a^2/(b*\sigma*cos\theta*Hv)\) and shape parameter
\(p3\) where \(\sigma\) denotes the stress, \(a\) the minor axis length and
\(Hv\) the Vickers hardness. The angle \(\theta\) is measured between the rotational axis
and the axis of main load direction. In this way we account for the orientation of particles (spheroids)
when generating fracture times dependent on their tendency to be more or less oriented towards the
main load direction.