粗糙节理岩石真三轴压缩数值仿真

软件：LSDYNA，具体方法不限，FEM、SPG、DEM等方法均可，网格类型不限。
问题描述：模型为50*50*100mm的长方体三维岩石模型，模型含有贯穿整体的倾斜粗糙节理（裂隙），节理为2种粗糙度X4种倾角，包括有孔模型和无孔模型。模型几何具体见附图。
无孔的粗糙节理模型，需要进行真三轴加载、真三轴卸载的数值仿真，即2种计算工况*2种粗糙度*4种倾角=16个具体的计算工况。
有孔模型需要设置锚杆，施加锚固预紧力，然后进行预加载——>单向卸载同时轴向加载的数值仿真，即1种计算工况*2种粗糙度*4种倾角=8个具体的计算工况。
可使用任何材料模型，提供完整岩石试件的常规三轴力学曲线用于调整材料模型参数。
任务内容包括建模和计算。
效果要求：要求结果能得到与实验相近应力应变曲线和破坏模式，利用删除失效单元或者塑性应变分布来模拟裂隙，最终交付运算结果和运算K文件。