simTimes: Crack time simulation

In simLife: Simulation of Fatigue Lifetimes

Description

Simulate crack times

Usage

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

Arguments

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

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.

Value

list of increasing failure times

Author(s)

M. Baaske

See Also

getCrackTime, getDelamTime

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.