首页 理论教育多尺度方法的算法流程与有效性验证

多尺度方法的算法流程与有效性验证

【摘要】:在平均化嵌套有限元多尺度计算过程中,宏观尺度应力并不是直接从宏观本构方程计算得到,需要通过细观尺度有限元计算所得细观应力进行平均化处理后得到。有关调用UMAT与ABAQUS内核语言Python并嵌入ABAQUS平台实现嵌套多尺度算法的更多阐述,有兴趣的读者可参考文献[24]。

嵌套多尺度算法可通过用户材料子程序UMAT与调用ABAQUS内核语言Python并嵌入ABAQUS平台实现。在平均化嵌套有限元多尺度计算过程中,宏观尺度应力并不是直接从宏观本构方程计算得到,需要通过细观尺度有限元计算所得细观应力进行平均化处理后得到。因此,可以根据UMAT完成平均化嵌套有限元多尺度计算的部分内容:通过UMAT提前宏观应变增量n+1G作为细观RVE有限元计算的边界条件,并在完成细观RVE有限元计算后对细观应力进行平均化处理,给出新的宏观等效雅可比矩阵DDSDDE及新的宏观等效应力n+1σ。

在平均化嵌套有限元多尺度计算过程中,在已知宏观应变增量的情况下需要进行细观RVE有限元有限元计算以获得细观应力分布进而获得宏观等效应力,也就是说,有限元软件需要完成自我调用,而据作者所知尚未有任何商用有限元软件能够实现此功能,因此本书通过有限元分析ABAQUS软件所提供的二次开发功能,调用ABAQUS脚本接口所提供的库函数,绕过ABAQUS/CAE界面,直接操纵ABAQUS内核完成细观RVE的建模、自平衡问题非线性分析、后处理分析结果等操作,并在用户子程序接口UMAT中通过RVE的细观迭代计算结果进行后处理操作,从而获得宏观等效应力和宏观切线模量,完成嵌套多尺度计算。有关调用UMAT与ABAQUS内核语言Python并嵌入ABAQUS平台实现嵌套多尺度算法的更多阐述,有兴趣的读者可参考文献[24]。