simTimes {simLife} | R Documentation |
Simulate crack times
simTimes(S, param, vickers, stress, cores = getOption("simLife.mc", 1L), verbose = FALSE)
S |
(non-overlapping) geometry system |
param |
parameters for generating failure times |
vickers |
vickers hardness |
stress |
stress level to be applied |
cores |
optional, number of cores for mulicore parallization with |
verbose |
logical, not used yet |
The function randomly generates phase dependant failure times of defect types crack
and delam
.
The accumulation structure of defects is initialized containing the failure times of objects in ascending order.
The failure times of the defect type crack
follow a Weibull distribution, see simCrackTime
.
The failure times of the defect type delam
roughly depend on the projected area of the object, the
applied overall stress amplitude and whether the object lies in the interior of the simulation box or hits
one of the box boundaries. The argument param
consists of a distribution parameter list which contains
numeric vectors for both reinforcement objects (labled by P
) and an optional second phase (labled by F
).
If no second phase is considered the corresponding parameter set is simply ignored. Each parameter vector is made up of
six parameters in the following order: p1 probability of already materialized defects, scale factor p2,
shape factor p3, shift parameter p4 of log times, the slope p5 and p6 as the standard deviation
of the random error of log times.
list of increasing failure times
M. Baaske
## Simulation of individual defect times ## of some particle system data(AL2MC_20p_k10_F2p_S) ## generate individual failure times opt <- list("vickers"=107,"distTol"=1,"Tmax"=10^11, "inAreafactor"=1.56, "outAreafactor"=1.43, "pointsConvHull"=10, "scale"=1e+06,"pl"=0) par <- list("P"=c(0.01,6,0.5,75,-15,3), "F"=c(0,0,0,105,-12,1), "const"=NULL) ## simulate times CLT <- simTimes(S,par,vickers=opt$vickers,stress=125,cores=1L) ## times T <- unlist(sapply(CLT,`[[`,"T")) V <- unlist(sapply(CLT,`[[`,"V")) U <- unlist(sapply(CLT,`[[`,"U")) ## show estimated densities showDensity(list("Delamination"=log10(V),"Crack"=log10(U),"Time"=log10(T)),xlim=c(-2,15))