For各项异性UMAT

最新推荐文章于 2024-03-07 09:00:00 发布

astragin

最新推荐文章于 2024-03-07 09:00:00 发布

阅读量933

点赞数

分类专栏： Fortran 文章标签： fortran

本文链接：https://blog.youkuaiyun.com/weixin_43725312/article/details/109135070

版权

这篇笔记介绍了如何在2020年10月17日使用Fortran编程实现一个各向异性本构模型，重点在于弹性各项异性材料的刚度矩阵，该矩阵遵循对称性原则。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

For各项异性UMAT

2020年10月17日
本次笔记想给大家呈现一个各项异性本构
弹性各项异性本构如下图所示，根据张量知识，刚度矩阵必为对称矩阵。
在这里插入图片描述

      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
     1 RPL,DDSDDT,DRPLDE,DRPLDT,
     2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
     3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
     4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)
C
      INCLUDE 'ABA_PARAM.INC'
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*********THE CONTENT OF UMAT***************C
      REAL ATEMP(3,3)
      PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.