各向同性线弹性材料本构模型umat的应用

个人专栏—ABAQUS专栏

1. Abaqus2023的用法教程——与VS2022、oneAPI 2024子程序的关联方法 Abaqus2023的用法教程——与VS2022、oneAPI 2024子程序的关联方法
2. Abaqus有限元分析——有限元网格划分基本原则 Abaqus有限元分析——有限元网格划分基本原则
3. 各向同性线弹性材料本构模型umat的应用 各向同性线弹性材料本构模型umat的应用
4. ABAQUS用户子程序二次开发——UMAT介绍（一） ABAQUS用户子程序二次开发——UMAT介绍（一）
5. Fortran语法介绍（一） Fortran语法介绍（一）
6. Fortran语法介绍（二） Fortran语法介绍（二）
7. Fortran语法介绍（三） Fortran语法介绍（三）

---

---

线弹性材料本构模型

• 对于多数材料而言，在微小变形的假设下，会满足线弹性理论，数学可以表示为：

$$\begin{array}{rl}{\sigma }_{ij}& =C_{ijkl}{\epsilon }_{kl}\\ & =\frac{E}{1+\nu }({\epsilon }_{ij}+\frac{\nu }{1-2\nu }{\epsilon }_{kk}{\delta }_{ij})\end{array}$$

$C_{ijkl}$是材料的弹性刚度，该四阶张量共有81个分量（$i、j、k、l$皆为1到3的整数）

• 考虑任意一个计算域内微小元素的力平衡、力矩平衡以及在微小变形理论中应变的计算得：
  $$\begin{array}{c}{\sigma }_{ij}={\sigma }_{ji}\\ {\epsilon }_{kl}=\frac{1}{2}(\frac{\partial u_k}{\partial x_l}+\frac{\partial u_l}{\partial x_k})\\ {\epsilon }_{kl}={\epsilon }_{lk}\\ \to C_{ijkl}=C_{jikl}=C_{ijlk}\end{array}$$
  可知刚度矩阵具有对称关系，故刚度矩阵元素中独立元素个数减少为36个。采用张量标记法Voigt标记刚度矩阵：
  $$\left[\begin{array}{c}{\sigma }_1\\ {\sigma }_2\\ {\sigma }_3\\ {\sigma }_4\\ {\sigma }_5\\ {\sigma }_6\end{array}\right]=\frac{E}{(1+\nu )(1-2\nu )}\left[\begin{array}{cccccc}1-\nu & \nu & \nu & 0& 0& 0\\ \nu & 1-\nu & \nu & 0& 0& 0\\ \nu & \nu & 1-\nu & 0& 0& 0\\ 0& 0& 0& \frac{1-2\nu }{2}& 0& 0\\ 0& 0& 0& 0& \frac{1-2\nu }{2}& 0\\ 0& 0& 0& 0& 0& \frac{1-2\nu }{2}\end{array}\right]\left[\begin{array}{c}{\epsilon }_1\\ {\epsilon }_2\\ {\epsilon }_3\\ 2{\epsilon }_4\\ 2{\epsilon }_5\\ 2{\epsilon }_6\end{array}\right]$$

• 为方便参数设定，采用参数替换

  1. 拉梅常数
     $$\lambda =\frac{E\nu }{(1+\nu )(1-2\nu )}$$
  2. 剪切模量
     $$\mu =\frac{E}{2(1+\nu )}$$
  3. 刚度矩阵
     $$C=\left[\begin{array}{cccccc}\lambda +2\mu & \lambda & \lambda & 0& 0& 0\\ \lambda & \lambda +2\mu & \lambda & 0& 0& 0\\ \lambda & \lambda & \lambda +2\mu & 0& 0& 0\\ 0& 0& 0& \mu & 0& 0\\ 0& 0& 0& 0& \mu & 0\\ 0& 0& 0& 0& 0& \mu \end{array}\right]$$

线弹性材料本构的umat编写

       SUBROUTINE UMAT(
C **** 参数列 ***********************************************C
     1 STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     2 RPL,DDSDDT,DRPLDE,DRPLDT,
     3 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     4 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     5 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC
C **** 宣告参数***********************************************C
       INCLUDE`ABA_PARAM.INC`
C------------------------------------------------------------C
       CHARACTER*80 CMNAME
       DIMENSION STRESS(NTENS),STATEV(NSTATV),
     1 DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
     2 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),
     3 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)
     4 JSTEP(4)
C------------------------------------------------------------C
       REAL*8 E, NU, MU, LAM, MU2
       INTEGER*4 I, J
C **** 材料模型主程序***********************************************C       
C **** 材料使用错误，停止分析
       IF (NTENS .EQ. 1) THEN
         WRITE(7,*) `错误：本线弹性模型不支持一维元素`
         CALL XIT
       ENDIF 
C **** 材料刚度系数归零
       DO I = 1,NTENS
         DO J = 1,NTENS
           DDSDDE(I,J) = 0.D0
         ENDDO
       ENDDO
C **** 材料参数读取与计算
       E  = PROPS(1)
       NU = PROPS(2)
       MU = E/(2.D0+2.DO*NU)
       LAM= MU*NU/(0.5D.0-NU)
C **** 剪应力和刚度矩阵
       DO J = 1,NSHR
          I = J+NDI
          DDSDDE(I,I) = MU
          STRESS(I) = STRESS(I)+ MU*DSTRAN(I)     
C **** 法向应力和刚度矩阵
       MU2 = 2.D0*MU
       DO I = 1,NDI
         DO J = 1,NDI
           DDSDDE(I,J) = LAM
           STRESS(I) = STRESS(I)+ LAM*DSTRAN(J)
         ENDDO
         DDSDDE(I,I) = LAM+ MU2
         STRESS(I) = STRESS(I) + MU2*DSTRAN(I)
       ENDDO
C------------------------------------------------------------C
       RETURN
       END SUBROUTINE UMAT

UMAT应用案例

• 假设一线弹性材料厚板，大小为$500mm\times 500mmm$，在中心挖一半径为$25mm$的圆，并在板两侧施加$0.1mm$的位移，已知材料的弹性模量$E=200GPa$，泊松比$\nu =0.3$。

---

欢迎对Abaqus感兴趣的朋友们查看：Abaqus-UMAT开发精品书籍及umat子程序学习

Abaqus非线性粘弹性模型子程序umat——广义MAXWELL粘弹性模型umat解析（朱-王-唐本构模型）

如果你喜欢以上内容，或者对Abaqus学习有兴趣，欢迎收藏关注，博主将持续更新。你的关注、收藏是我持续创作的动力！