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

科研拓展人

已于 2024-07-12 17:16:51 修改

阅读量2.2k

点赞数 57

分类专栏： Abaqus应用文章标签： ABAQUS UMAT 弹性材料本构模型 ZWT 朱王唐本构模型

于 2024-03-06 09:00:00 首次发布

本文链接：https://blog.csdn.net/Zh531445436/article/details/136473332

版权

Abaqus应用专栏收录该内容

6 篇文章 6 订阅

订阅专栏

个人专栏—ABAQUS专栏

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

文章目录

个人专栏—ABAQUS专栏

线弹性材料本构模型

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

$\begin{align*} \sigma_{ij}&=C_{ijkl}\varepsilon_{kl}\\ &=\frac{E}{1+\nu}(\varepsilon_{ij}+\frac{\nu}{1-2\nu}\varepsilon_{kk}\delta_{ij}) \end{align*}$

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

考虑任意一个计算域内微小元素的力平衡、力矩平衡以及在微小变形理论中应变的计算得：
$\begin{gather*} \sigma_{ij}=\sigma_{ji}\\ \varepsilon_{kl}=\frac{1}{2}(\frac{\partial u_k}{\partial x_l}+\frac{\partial u_l}{\partial x_k})\\ \varepsilon_{kl}=\varepsilon_{lk}\\ \rightarrow \quad C_{ijkl}=C_{jikl}=C_{ijlk} \end{gather*}$
可知刚度矩阵具有对称关系，故刚度矩阵元素中独立元素个数减少为36个。采用张量标记法Voigt标记刚度矩阵：
$\begin{bmatrix} \sigma_1\\ \sigma_2\\ \sigma_3\\ \sigma_4\\ \sigma_5\\ \sigma_6 \end{bmatrix}=\frac{E}{(1+\nu)(1-2\nu)} \begin{bmatrix} 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{bmatrix} \begin{bmatrix} \varepsilon_1\\ \varepsilon_2\\ \varepsilon_3\\ 2\varepsilon_4\\ 2\varepsilon_5\\ 2\varepsilon_6 \end{bmatrix}$
为方便参数设定，采用参数替换
1. 拉梅常数
  $\lambda=\frac{E\nu}{(1+\nu)(1-2\nu)}$
2. 剪切模量
  $\mu=\frac{E}{2(1+\nu)}$
3. 刚度矩阵
  $C=\begin{bmatrix} \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{bmatrix}$

线弹性材料本构的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