abaqus简单立方体胞元周期性边界条件施加计算脚本源代码

abaqus中周期性边界条件的施加一般通过方程约束，手动设置不仅繁琐而且很容易出错。根据文献《Unit cells for micromechanical analyses of particle-reinforced composites》中简单立方体胞元周期性边界条件的施加方法，开发Python脚本，可以根据用户提供的三维数组创建网格，并施加周期性边界条件以及自动提交abaqus计算。在此提供程序的Python源代码，和大家一起学习。代码中重要语句都进行了注释，对照参考文献可以很好的理解周期性边界条件施加过程及方法，代码书写不易，希望大家多多支持，共同进步。

部分代码片段如下图所示，

用户首先需提供一个三维数组，如上图中tensor，其中1表示此处有材料，0表示此处无材料，则其形成的网格

Tensor可以手工编辑亦可以由其他程序生成，但需保证其对称性。

在abaqus中通过file→run script选择该脚本运行，脚本会输出X，Y，Z三个方向的计算结果ODB文件。