2D Y-code and some 3D Y-codes about FDEM can be found in vgest.net. I am simulating the fracture patterns of single crushing process of a quartz sand, but the web does not offer 3D Y-codes about breakage. I have emailed nearly all of the researchers whose paper I have read relies in 3D FDEM codes, unfortunately nobody is willing to give them to me for a research purpose. 

Thanks for reading my question.