using monte carlo simulation to estimate the percolation threshold of composite filled with conductive fillers

I want to create a 3D rectangular space in Matlab and randomly generate a randomly-oriented cylinder (with define length and width) within that 3D rectangular space. I want to keep adding these randomly oriented cylinder until global connectivity (percolation threshold) is established throughout the 3D rectangular space that I created initially.

My understanding is that I would need to: 1) define a 3D rectangular space 2) define the dimensions of the cylinder 3) randomly generate the cylinder in the 3D rectangular space 4) check for global connectivity, if not then continue to generate more cylinder. If yes, calculate the volume fraction of cylinder require for global connectivity 5) if possible, plot the 3D rectangular space filled with cylinder for visualisation

Please help me with the matlab coding.

I would appreciate any help I can get. Does not even have to be a matlab code for my problem. Any hints, tips, references to matlab theory that would help me. I just started with matlab and need some guidance. Any ideas?

Answer 1

[Image Analyst]

Try cylinder() and rand().