之前也传过一个例子，但是那个里面纤维不能在模型边界上。 现在做了一些改进。
2D的程序包括赋材料，定义周期性边界条件。运行之后可以直接提交计算。  程序里面有比较详细的注释。
3D的程序只是建模过程。
直接从‘run script’运行就可以了。
附上两张运行之后的截图：

大神啊                 

谢谢楼主的无私，            

感谢感谢感谢

顺便问下，python读取ODB中的单元应力是什么函数？

里面有个dispersion的参数是什么意思啊，求楼主qq

感谢，学习中~~

请问你这个到时候如何在界面嵌入内聚力单元？

胡佳成 发表于 2018-10-31 10:21
请问你这个到时候如何在界面嵌入内聚力单元？

我这个模型是最早学python的时候做的玩的。 这里面的随机纤维的生成算法不好，用这种算法生成的纤维的体积含量不能太高。我后来弄了一个随机移动纤维的算法可以实现大约90%的体积含量。
界面这个很好办，因为纤维的圆心都生成了，你只要给每个纤维画一个半径比它小一点点的同心圆，再做一下布尔运算就可以了。

你的意思是指建立有限厚度的内聚力单元吗？一般是零厚度的，不过有限厚度也未尝不可。