Perform a probabilistic damage tolerance analysis on a crackR object

Description

This is the main analysis function for the sequential importance sampling approach to probabilistic damage tolerance analysis in the crackR package. User can enter an object of class crackRparameters or class crackR. If a crackRparameters object is entered, crackRinit will be run to initialize a crackR object. The sample size for an existing crackR object can be increased by utilizing an existing crackR object and specifying add=TRUE.

Usage

1
analyze(obj,add=FALSE)

Arguments

obj

Object of class Sing, Mult, CD, or an appropriate list of input parameters

add

Logical, indicating whether this run will add trials to an existing run

Details

If add=TRUE, the number of particles specified in the input object parameter Np will be added to the existing particle count and the results will be weighted according to the existing particle count. The new set of results, containing only the results of the subset run, will be included in a new component of the output: 'new.results'. The cumulative results will as before be contained in component 'results'. Note, if the input object was already initialized using crackRinit, it's class will be either Sing, Mult, or CD. If an uninitialized list of crackR parameters is submitted, this function will attempt to initialize it.

Value

state

List containing vectors representing the final state of each particle

parameters

List of the analysis parameters (not altered by the running of analyze

results

List of the results of a single run (where add=FALSE) or of a cumulative set of runs (when add=TRUE)

new.results

List of the results for the latest run in a cumulative set of runs (when add=TRUE)

Author(s)

Keith Halbert <keith.a.halbert@gmail.com>

See Also

crackRinit

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
set.seed(327)

## cp7ext data is part of this package
data(cp7ext)

## prepare a crackR object for analysis
cp7ext.init <- crackRinit(cp7ext)

## perform analysis
cp7ext.out  <- analyze(cp7ext.init)

## Plot SFPOF and PCD results
plot(cp7ext.out)