showDensity: Plot estimated densities

In simLife: Simulation of Fatigue Lifetimes

Description

Plot the estimated densities which result from the randomly generated individual failure times.

Usage

showDensity(dv, main = "Failure time density estimation", ...)

Arguments

dv	named list of individual failure times
main	title of the plot (optional)
...	optional, additional graphics parameters

Author(s)

M. Baaske

Examples

## 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))

Example output

simLife documentation built on May 2, 2019, 6:36 a.m.