前言

记录自己学习过程，本文针对《非线性本构关系在ABAQUS中的实现》第二章“动态回弹模量“本构的学习，书中自带程序可能由于转成PDF的原因有一些排版上的混乱，下文的UMAT子程序，有VS code编写，和书中例子基本相同。

umat子程序

      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)
CINCLUDE 'ABA_PARAM.INC'
CCHARACTER*80 CMNAMEDIMENSION 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    局部变量DIMENSION SDEV(NTENS)REAL THETA, THETA0, TAUOC, TAUOC0, DNU, DK, DN, DM, PALOGICAL SINGULAR
C     DNU:Poisson's ratio => PROPS(1)
C     DK: Material parameter, k1=> PROPS(2)
C     DN: Material parameter,k2=> PROPS(3)
C     DM: Material parameter,k3=> PROPS(4)
C     PA: atmo=> PROPS(5)
C     DL: Material parameter,k4=> PROPS(6)DNU=PROPS(1)DK=PROPS(2)DN=PROPS(3)DM=PROPS(4)PA=PROPS(5)DL=PROPS(6)C     计算弹性常数ALPHA=DNU/(1-2.0*DNU)ALPHA_BAR=3.0*ALPHA+1DLL=1.0/DL
C     偏应力与等效偏应力THETA=STRESS(1)+STRESS(2)+STRESS(3)DO I=1,NDISDEV(I)=STRESS(I)-THETA*0.333333ENDDODO I=NDI+1,NTENSSDEV(I)=STRESS(I)ENDDOTAUOC=0.0DO I=1,NTENSTAUOC=TAUOC+SDEV(I)**2ENDDOTAUOC=SQRT(TAUOC/3.0)THETA0=THETATAUOC0=TAUOCIF(THETA .GT. -0.01) THETA0=-0.01IF(TAUOC .LT. 0.00001) TAUOC0=0.00001C     更新恢复模量（Resilient modulus）RMR=DK*PA*(-THETA0/PA)**DN*(DL*TAUOC0/PA+1)**DMSTATEV(1)=RMRSINGULAR=.TRUE.DO I=1,NTENSIF (STRESS(I) .LT. -0.1) SINGULAR=.FALSE.ENDDOIF (SINGULAR) THENA=ALPHAB=0.0C=0.0D=0.0ELSEXX=TAUOC+PA*DLLYY=1.0/(XX*(1.0-DN)-DM*TAUOC)A=XX*YY*ALPHA_BAR*DN/3.0+ALPHAB=ALPHA_BAR*DN*XX*YY/THETAC=THETA*DM*YY/(9*TAUOC)END IFDDSDDE=0.0DO I=1,NDIDDSDDE(I,I)=1.0ENDDODO I=NDI+1,NTENSDDSDDE(I,I)=0.5ENDDODO I=1,NDIDO J=1,NDIDDSDDE(I,J)=DDSDDE(I,J)+A+B*SDEV(I)+C*SDEV(J)&     +D*SDEV(I)*SDEV(J)ENDDOENDDODO I=1,NDIDO J=NDI+1,NTENSDDSDDE(I,J)=DDSDDE(I,J)+C*SDEV(J)+D*SDEV(I)*SDEV(J)DDSDDE(J,I)=DDSDDE(J,I)+C*SDEV(J)+D*SDEV(I)*SDEV(J)ENDDOENDDODO I=NDI+1,NTENSDO J=NDI+1,NTENSDDSDDE(I,J)=DDSDDE(I,J)+D*SDEV(I)*SDEV(J)ENDDOENDDO
C     应力更新DO I=1,NTENSDO J=1,NTENSSTRESS(I)=STRESS(I)+DDSDDE(I,J)*DSTRAN(J)ENDDOENDDORETURNEND

DeBug

1）log文件中修改umat的语法错误

2）没有语法错误依旧出现：
Error in job Job-Lug: Abaqus/Standard Analysis exited with an error - Please see the message file for possible error messages if the file exists.
Job Job-load aborted due to errors
解决办法：1边界条件设置一定要从初始步开始
2.载荷施加设置分析步施加较小载荷提供缓冲
3.分析步的分析时间设置initial为更小的值（我这里0.05）

计算结果

ABAQUS材料子程序学习（20年12月1日）相关推荐

ABAQUS材料子程序学习（20年12月2日）
@ABAQUS材料子程序学习(20年12月2日) 前言继续记录自己学习过程,本文针对<非线性本构关系在ABAQUS中的实现>第三章"黏弹性"本构的学习,UMAT子程序 ...
ABAQUS材料子程序学习（线性各向同性硬化塑性）
ABAQUS材料子程序学习(线性各向同性硬化塑性) 前言塑性力学增量形式实现 umat子程序参数计算结果前言记录自己学习abaqus软件umat子程序的t过程,本文主要参考了<非线性本 ...
Python学习_2015年12月14日
2015年12月14日 1.Python代码中有成片代码需要加缩进?(比如少加了一个if,后面的代码块都要缩进) 解决方法:选择要缩进的一片代码,然后按Tab,所有行一起缩进.如果要取消缩进,Shif ...
Meta分析和网状Meta分析速成班（ 2019年12月28-29日上海）
适合对象零基础的研究生.临床医务人员.医药院校的科研工作者:迫切希望发表SCI文章晋升的临床工作者.本次速成班特别适合初学者和无基础的学员. 课程背景本培训班以1篇meta分析和网状meta分析为 ...
清华大学教评中心＜jpxzx2019@mail.tsinghua.edu.cn＞ 12月24日9:00-11:10试讲评审材料（庞观和朱峰老师）
➤01 邮件通知卓老师好!感谢您担任线上试讲评审专家! 1.评审时间: 2020年12月24日(周四)9:00-11:10 2.试讲教师及内容: 12月24日上午上午试讲时间 9:00-10 ...
分享Silverlight/WPF/Windows Phone/HTML5一周学习导读(11月28日-12月4日)
分享Silverlight/WPF/Windows Phone/HTML5一周学习导读(11月28日-12月4日) 本周Silverlight学习资源更新 Silverlight HttpUtil 封 ...
2017年12月6日学习笔记（JAVA面试题）
前段时间都在泡图书馆,看<java核心技术卷1>.感觉好多东西会用了,但是并不知道底层实现原理,还有好多没见过的知识点概念.觉得果然还是要多读书. 接下来准备边读书,边把各种面试题提及的知 ...
12月19日 OpenCV 实战基础学习笔记——特征匹配
文章目录前言一.特征匹配 1.Brute-force 蛮力匹配 2.1 对 1 匹配 3.k 对最佳匹配二.答题卡识别前言本文为12月19日 OpenCV 实战基础学习笔记--特征匹配,分为 ...
20计算机二级考试时间,开考省份全部发布｜20年12月计算机二级考试报名时间
20年12月计算机二级考试通知[最新] 一.福建省 [考试时间]12月5日至6日 [报名时间和方式] 2.网上报名时间:2020年11月2日9时-11月13日24时止. 3.准考证打印时间:2020年 ...

ABAQUS材料子程序学习（20年12月1日）

ABAQUS材料子程序学习（20年12月1日）

前言

umat子程序

DeBug

计算结果

ABAQUS材料子程序学习（20年12月1日）相关推荐

最新文章

热门文章