abaqus编写本构方程vumat_基于ABAQUS的木材本构关系数值模拟方法与流程

本发明属于木材本构研究技术领域，提供了木材在复杂应力状态下应力-应变关系的数值模拟方法。

背景技术：

木材是各向异性材料，其本构的复杂主要表现为在受压作用下发生塑性变形，而在拉、剪作用下发生脆性破坏，同时各向拉、压强度均不相等。在木结构设计中，除极少数情况外(如单轴受拉、受压或受剪)，木构件受力是复杂的，一个能准确描述木材各向异性本构关系的模型是木结构有限元模型的重要组成部分。

在有限元分析中，木材通常可以被简化为正交各向异性材料。目前，ABAQUS内嵌的材料本构中，针对正交各向异性材料发生受压塑性变形的只有Hill屈服准则，但该本构模型不能同时描述木材顺纹理想塑性和横纹塑性硬化。

此外，ABAQUS内嵌的渐进损伤和失效模型尽管可以描述拉、剪作用下发生的脆性破坏，但会引起顺纹和横纹方向同时失效，这与木材的真实受力情况不符。综上所述，在目前的ABAQUS软件中，其内嵌的材料本构还不能准确描述木材在不同方向受压的塑性变形和拉、剪的脆性破坏。

技术实现要素：

本发明提供了一种基于ABAQUS的木材本构关系数值模拟方法，通过嵌入自定义的木材本构子程序VUMAT，能准确描述木材的正交各向异性，并预测木材在复杂受力状态下的破坏模式。

本发明的技术方案：

基于ABAQUS的木材本构关系数值模拟方法，步骤如下：

步骤一：建立弹性阶段本构方程，表达弹性阶段正交各向异性材料的应力-应变关系，如式(1)所示：

σ＝Dε (1)

式中：σ为应力矩阵、ε为应变矩阵、D为刚度矩阵；

展开式(1)得到式(2)：

式中，D11＝E1(1-v23v32)γ，D22＝E2(1-v13v31)γ，D33＝E3(1-v21v12)γ，D12＝E1(v21+ν23ν31)γ，D13＝E1(ν21+ν23ν31)γ，D23＝E2(ν32+ν12ν31)γ，D44＝G12，D55＝G23，D66＝G13，γ＝(1-ν21ν12-ν23ν32-v13v31-2v21v32v13)-1；

E1，