使用 Abaqus 来建立颗粒/短纤维改性复合材料模型 （最详细的教程）

复合材料可以发挥各种材料的优点，克服单一材料的缺陷，提高强度、弹性模量及断裂韧性等,扩大材料的应用范围。由于单一材料性质的改变比较困难，但将不同性质的材料相互复合能够形成种类丰富的材料,而表现出千差万别的材料性能，所以复合材料在新材料技术的研发过程中占有重要作用，其设计的自由性是其他材料所无法比拟的。复合材料由于其优良的综合性能，特别是其性能的可设计性被广泛应用于航空航天、国防、交通、体育等领域。
5.1 改性复合材料的本构理论与建模方法
颗粒或短纤维增强复合材料属于复合材料家族，与连续纤维增强树脂基复合材料相比，其材料可设计性更为灵活多样，这种针对材料性能设计的科研活动，根据材料体系的是否发生根本变化和是否有全新材料性能的出现，通常可分为新材料研发或原材料改性。材料的改性是材料研究领域的一项长期普遍的研究内容。针对力学性能研究，首先需要研究该材料的本构关系。研究方法一般是利用基于应变梯度塑性理论的本构方程表征局域区和内部结构的长度-尺度依赖性，此外，用应变梯度塑性理论进行正则化，大大降低了计算力学中广为接受的有限单元法计算结果对网格的依赖性。因此，应变梯度理论似乎提供一个直接的塑性流动的尺寸效应进行建模在连续介质力学方法。因此，应变梯度理论提供了一种在连续介质力学方法中模拟塑性流动尺寸效应的直接方法。
为了表示本构方程，可采用更新的拉格朗日公式，其中每个物质粒子用坐标标记。对于弹塑性固体，假设变形速率梯度的对称部分分解为弹性部分和塑性部分。进而有
(5.1.1)

弹性响应假定是各向同性的，并由胡克定律表示：
(5.1.2)

其中G为弹性剪切模量，为泊松比，, 分别为基尔霍夫应力和柯西应力在当前状态下的Jaumann速率。
利用Mises型流动规律，得到了与应变梯度相关的塑性应变率的本构方程为：

(5.1.3)

式中为柯西应力偏量部分，为等效应力，为等效塑性应变率。
与应变速率无关的材料的屈服函数变为：

(5.1.4)
与应变梯度相关的流动应力可以沿着相同的线表示，

(5.1.5)
其中，是常规的加工硬化函数，第二项表示应变梯度相关项。是描述对应变梯度依赖程度的系数。通过式(5-4)和式(5-5),一致性条件f=0变为：

(5.1.6)