Return to Examples Main | Next Example

Use of Random Variables

Lamina stiffness properties as random variables

For examples related to lamination theory, we need the following package from the MSC library.

The package msc`laminate` defines many lamination theory routines. Look at the help for Amat to see how the material properties and laminates get entered.

engineer can be used to get the laminate Young's modulus.

Lets assume all the material properties are independent random variables (physically impossible, but..). All the standard deviations are 10% of the mean. The stack is quasi-isotropic

Note that when the function returns a vector, errorPropagation will determine the statistics for each element of the vector

Ex | NormalDistribution[7.72621, 0.674185] |

Ey | NormalDistribution[7.72621, 0.674185] |

Gxy | NormalDistribution[2.99468, 0.254032] |

Nuxy | NormalDistribution[0.289741, 0.00836792] |

The SeriesApproximation method works remarkable well.

Ex | NormalDistribution[7.70408, 0.672427] |

Ey | NormalDistribution[7.70408, 0.672427] |

Gxy | NormalDistribution[2.98664, 0.253688] |

Nuxy | NormalDistribution[0.289523, 0.00840456] |

For this example, there is no difference between analytic and numerical derivatives. The numerical is the default so that more complex models can be inserted without any problem.

Ex | NormalDistribution[7.70408, 0.672427] |

Ey | NormalDistribution[7.70408, 0.672427] |

Gxy | NormalDistribution[2.98664, 0.253688] |

Nuxy | NormalDistribution[0.289523, 0.00840456] |

As expected, the IntervalMath method is lousy.

Ex | NormalDistribution[8.9164, 4.74047] |

Ey | NormalDistribution[8.9164, 4.74047] |

Gxy | NormalDistribution[3.56453, 1.98968] |

Nuxy | NormalDistribution[0.538001, 0.454534] |

Define a function that only returns the E1 component

Plot the CDF of the Monte Carlo simulation, along with the estimated continuous distribution.

Behavior of the mean.

An MSC derivation showed that the mean of a laminate property computed using random lamina inputs will always be less than the point value calculated from the mean of the lamina inputs.

Load the composite lamination theory package if not already loaded

Define a function that only returns the component

Mean values of the lamina stiffness

computed from the means

Mean of Ex computed from random variables with stddev of 10%

Increase standard deviation of each input to 20%

So the hypothesis holds, but its very weak for this example.

Angle as a random variable

Consider the case of a +/- 30 laminate where the angle is a random variable with a standard deviation of 1 degree

Load the composite lamination theory package if not already loaded

Define a function that only returns the component

Check the SeriesApproximation method for this case

Make sure the distribution still looks normal.

Created by Mathematica (March 7, 2004)