simTimes: Crack time simulation

Description

Simulate crack times

Usage

simTimes(S, param, vickers, stress, cores = getOption("simLife.mc", 1L),
  verbose = FALSE)

Arguments

(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 cores=1L (default) by mclapply which also can be set by a global option "simLife.mc"

verbose

logical, not used yet

Value

list of increasing failure times

Details

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.

Examples

Run this code

# NOT RUN {
## 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))
# }