abaqus邓肯张模型umat
首先是始点刚度法：
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
RPL,DDSDDT,DRPLDE,DRPLDT,
STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C 始点刚度法
CHARACTER* CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),
DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
STRAN(NTENS),DSTRAN(NTENS),TIME(),PREDEF(),DPRED(),
PROPS(NPROPS),COORDS(),DROT(,),DFGRD0(,),DFGRD1(,),
JSTEP()
DIMENSION PS(),DSTRESS(NTENS)
PARAMETER (ONE=.0D0,TWO=.0D0,THREE=.0D0,SIX=.0D0) EK=PROPS()
EN=PROPS()
RF=PROPS()
C=PROPS()
FAI=PROPS()/180.0*3.1415926
UG=PROPS()
UD=PROPS()
UF=PROPS()
EKUR=PROPS()
PA=PROPS()
DFAI=PROPS()/180.0*3.1415926
S1S3O=STATEV()
S3O=STATEV()
SSS=STATEV()
CALL GETPS(STRESS,PS,NTENS) FAI=FAI-DFAI*LOG10(S3O/PA)
CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF
,SSS,S1S3O)
EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
DSTRESS=0.0
CALL GETSTRESS(DDSDDE,DSTRESS,DSTRAN,NTENS) DO I1=,NTENS
STRESS(I1)=STRESS(I1)+DSTRESS(I1)
CONTINUE
CALL GETPS(STRESS,PS,NTENS)
CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF,
SSS,S1S3O)
EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
IF(PS().GT.S3O)S3O=PS()
IF((PS()-PS()).GT.S1S3O)S1S3O=PS()-PS()
IF(S.GT.SSS)SSS=S
STATEV()=S1S3O
STATEV()=S3O
STATEV()=SSS
END SUBROUTINE GETPS(STRESS,PS,NTENS)
INCLUDE 'ABA_PARAM.INC'
DIMENSION PS(),STRESS(NTENS)
CALL SPRINC(STRESS,PS,,,)
DO I=,
DO J=I+,
IF(PS(I).GT.PS(J))THEN
PPS=PS(I)
PS(I)=PS(J)
PS(J)=PPS
END IF
CONTINUE
CONTINUE
DO K1=,
PS(K1)=-PS(K1)
CONTINUE
RETURN
END SUBROUTINE GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O
,UG,UD,UF,SSS,S1S3O)
INCLUDE 'ABA_PARAM.INC'
DIMENSION PS()
S=(-SIN(FAI))*(PS()-PS())
IF(PS().LT.(-C/TAN(FAI))) THEN
S=0.99
ELSE S=S/(*C*COS(FAI)+*PS()*SIN(FAI))
IF(S.GE.0.99)then
S=0.99
end if
END IF
EMOD=EK*PA*((S3O/PA)**EN)*((-RF*S)**) AA=UD*(PS()-PS())
AA=AA/(EK*PA*((S3O/PA)**EN))
AA=AA/(-RF*S)
ENU=UG-UF*LOG10(S3O/PA)
ENU=ENU/(-AA)/(-AA)
IF(ENU.GT.0.49)ENU=0.49
IF(ENU.LT.0.05)ENU=0.05
IF(S.LT.SSS.AND.(PS()-PS()).LT.S1S3O)THEN
EMOD=EKUR*PA*((S3O/PA)**EN)
END IF
END SUBROUTINE GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
INCLUDE 'ABA_PARAM.INC'
DIMENSION DDSDDE(NTENS,NTENS)
DO K1=,NTENS
DO K2=,NTENS
DDSDDE(K2,K1)=0.0
CONTINUE
CONTINUE
DO K1=,NDI
DO K2=,NDI
DDSDDE(K2,K1)=ELAM
CONTINUE
DDSDDE(K1,K1)=EG2+ELAM
CONTINUE
DO K1=NDI+,NTENS
DDSDDE(K1,K1)=EG
CONTINUE
RETURN
END SUBROUTINE GETSTRESS(DDSDDE,STRESS,DSTRAN,NTENS)
INCLUDE 'ABA_PARAM.INC'
DIMENSION DDSDDE(NTENS,NTENS),STRESS(NTENS),DSTRAN(NTENS)
DO K1=,NTENS
DO K2=,NTENS
STRESS(K2)=STRESS(K2)+DDSDDE(K2,K1)*DSTRAN(K1)
CONTINUE
CONTINUE
RETURN
END
abaqus邓肯张umat始点刚度法
abaqus邓肯张umat中点刚度法：
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
RPL,DDSDDT,DRPLDE,DRPLDT,
STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C 中点刚度法
CHARACTER* CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),
DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
STRAN(NTENS),DSTRAN(NTENS),TIME(),PREDEF(),DPRED(),
PROPS(NPROPS),COORDS(),DROT(,),DFGRD0(,),DFGRD1(,),
JSTEP() DIMENSION PS(),DSTRESS(NTENS),TDSTRESS(NTENS),TSTRESS(NTENS)
PARAMETER (ONE=.0D0,TWO=.0D0,THREE=.0D0,SIX=.0D0) EK=PROPS()
EN=PROPS()
RF=PROPS()
C=PROPS()
FAI=PROPS()/180.0*3.1415926
UG=PROPS()
UD=PROPS()
UF=PROPS()
EKUR=PROPS()
PA=PROPS()
DFAI=PROPS()/180.0*3.1415926
S1S3O=STATEV()
S3O=STATEV()
SSS=STATEV()
CALL GETPS(STRESS,PS,NTENS)
FAI=FAI-DFAI*LOG10(S3O/PA)
CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF
,SSS,S1S3O)
EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
TDSTRESS=0.0
CALL GETSTRESS(DDSDDE,TDSTRESS,DSTRAN,NTENS)
DO I1=,NTENS
TSTRESS(I1)=STRESS(I1)+TDSTRESS(I1)*0.5
CONTINUE CALL GETPS(TSTRESS,PS,NTENS)
CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF,
SSS,S1S3O)
EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
DSTRESS=0.0
CALL GETSTRESS(DDSDDE,DSTRESS,DSTRAN,NTENS)
DO I1=,NTENS
STRESS(I1)=STRESS(I1)+DSTRESS(I1)
CONTINUE
CALL GETPS(STRESS,PS,NTENS)
CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF,
SSS,S1S3O)
EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG) IF(PS().GT.S3O)S3O=PS()
IF((PS()-PS()).GT.S1S3O)S1S3O=PS()-PS()
IF(S.GT.SSS)SSS=S
STATEV()=S1S3O
STATEV()=S3O
STATEV()=SSS END SUBROUTINE GETPS(STRESS,PS,NTENS)
INCLUDE 'ABA_PARAM.INC'
DIMENSION PS(),STRESS(NTENS)
CALL SPRINC(STRESS,PS,,,)
DO I=,
DO J=I+,
IF(PS(I).GT.PS(J))THEN
PPS=PS(I)
PS(I)=PS(J)
PS(J)=PPS
END IF
CONTINUE
CONTINUE
DO K1=,
PS(K1)=-PS(K1)
CONTINUE
RETURN
END SUBROUTINE GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O
,UG,UD,UF,SSS,S1S3O)
INCLUDE 'ABA_PARAM.INC'
DIMENSION PS()
S=(-SIN(FAI))*(PS()-PS())
IF(PS().LT.(-C/TAN(FAI))) THEN
S=0.99
ELSE
S=S/(*C*COS(FAI)+*PS()*SIN(FAI))
IF(S.GE.0.99)S=0.99
END IF
EMOD=EK*PA*((S3O/PA)**EN)*((-RF*S)**)
AA=UD*(PS()-PS())
AA=AA/(EK*PA*((S3O/PA)**EN))
AA=AA/(-RF*S)
ENU=UG-UF*LOG10(S3O/PA)
ENU=ENU/(-AA)/(-AA)
IF(ENU.GT.0.49)ENU=0.49
IF(ENU.LT.0.05)ENU=0.05
IF(S.LT.SSS.AND.(PS()-PS()).LT.S1S3O)THEN
EMOD=EKUR*PA*((S3O/PA)**EN)
END IF
END SUBROUTINE GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
INCLUDE 'ABA_PARAM.INC'
DIMENSION DDSDDE(NTENS,NTENS)
DO K1=,NTENS
DO K2=,NTENS
DDSDDE(K2,K1)=0.0
CONTINUE
CONTINUE
DO K1=,NDI
DO K2=,NDI
DDSDDE(K2,K1)=ELAM
CONTINUE
DDSDDE(K1,K1)=EG2+ELAM
CONTINUE
DO K1=NDI+,NTENS
DDSDDE(K1,K1)=EG
CONTINUE
RETURN
END SUBROUTINE GETSTRESS(DDSDDE,STRESS,DSTRAN,NTENS)
INCLUDE 'ABA_PARAM.INC'
DIMENSION DDSDDE(NTENS,NTENS),STRESS(NTENS),DSTRAN(NTENS)
DO K1=,NTENS
DO K2=,NTENS
STRESS(K1)=STRESS(K1)+DDSDDE(K1,K2)*DSTRAN(K2)
CONTINUE
CONTINUE
RETURN
END
abaqus邓肯张umat中点刚度法
上面的是E-v模型。
abaqus邓肯张EB模型umat：
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,
RPL,DDSDDT,DRPLDE,DRPLDT,
STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME,
NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,
CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,JSTEP,KINC)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER* CMNAME
DIMENSION STRESS(NTENS),STATEV(NSTATV),
DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),
STRAN(NTENS),DSTRAN(NTENS),TIME(),PREDEF(),DPRED(),
PROPS(NPROPS),COORDS(),DROT(,),DFGRD0(,),DFGRD1(,),
JSTEP()
DIMENSION PS(),DSTRESS(NTENS)
PARAMETER (ONE=.0D0,TWO=.0D0,THREE=.0D0,SIX=.0D0) EK=PROPS()
EN=PROPS()
RF=PROPS()
C=PROPS()
FAI=PROPS()/180.0*3.1415926
DFAI=PROPS()/180.0*3.1415926
EKB=PROPS()
EM=PROPS()
EKUR=PROPS()
PA=PROPS()
S1S3O=STATEV()
S3O=STATEV()
SSS=STATEV() CALL GETPS(STRESS,PS,NTENS)
FAI=FAI-DFAI*LOG10(S3O/PA)
CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,EKB,EM
,SSS,S1S3O) EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE
CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
DSTRESS=0.0
CALL GETSTRESS(DDSDDE,DSTRESS,DSTRAN,NTENS)
DO I1=,NTENS
STRESS(I1)=STRESS(I1)+DSTRESS(I1)
CONTINUE CALL GETPS(STRESS,PS,NTENS) CALL GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,EKB,EM
,SSS,S1S3O) EBULK3=EMOD/(ONE-TWO*ENU)
EG2=EMOD/(ONE+ENU)
EG=EG2/TWO
EG3=THREE*EG
ELAM=(EBULK3-EG2)/THREE CALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG) IF(PS().GT.S3O)S3O=PS()
IF((PS()-PS()).GT.S1S3O)S1S3O=PS()-PS()
IF(S.GT.SSS)SSS=S
STATEV()=S1S3O
STATEV()=S3O
STATEV()=SSS END SUBROUTINE GETPS(STRESS,PS,NTENS)
INCLUDE 'ABA_PARAM.INC'
DIMENSION PS(),STRESS(NTENS)
CALL SPRINC(STRESS,PS,,,)
DO I=,
DO J=I+,
IF(PS(I).GT.PS(J))THEN
PPS=PS(I)
PS(I)=PS(J)
PS(J)=PPS
END IF
CONTINUE
CONTINUE
DO K1=,
PS(K1)=-PS(K1)
CONTINUE
RETURN
END SUBROUTINE GETEMOD(PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O
,EKB,EM,SSS,S1S3O) INCLUDE 'ABA_PARAM.INC'
DIMENSION PS() S=(-SIN(FAI))*(PS()-PS())
IF(PS().LT.(-C/TAN(FAI))) THEN
S=0.99
ELSE S=S/(*C*COS(FAI)+*PS()*SIN(FAI))
IF(S.GE.0.99)then
S=0.99
END IF
END IF
EMOD=EK*PA*((S3O/PA)**EN)*((-RF*S)**)
IF((S.LT.SSS).AND.((PS()-PS()).LT.S1S3O))THEN
EMOD=EKUR*PA*((S3O/PA)**EN)
END IF AA=EKB*PA*((S3O/PA)**EM)
ENU=0.5*(1.0-EMOD/./AA)
IF(ENU.GT.0.49)ENU=0.49
IF(ENU.LT.0.05)ENU=0.05
END SUBROUTINE GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)
INCLUDE 'ABA_PARAM.INC'
DIMENSION DDSDDE(NTENS,NTENS)
DO K1=,NTENS
DO K2=,NTENS
DDSDDE(K2,K1)=0.0
CONTINUE
CONTINUE
DO K1=,NDI
DO K2=,NDI
DDSDDE(K2,K1)=ELAM
CONTINUE
DDSDDE(K1,K1)=EG2+ELAM
CONTINUE
DO K1=NDI+,NTENS
DDSDDE(K1,K1)=EG
CONTINUE
RETURN
END SUBROUTINE GETSTRESS(DDSDDE,STRESS,DSTRAN,NTENS)
INCLUDE 'ABA_PARAM.INC'
DIMENSION DDSDDE(NTENS,NTENS),STRESS(NTENS),DSTRAN(NTENS)
DO K1=,NTENS
DO K2=,NTENS
STRESS(K2)=STRESS(K2)+DDSDDE(K2,K1)*DSTRAN(K1)
CONTINUE
CONTINUE
RETURN
END
abaqus邓肯张EB模型UMAT
abaqus邓肯张模型umat的更多相关文章
- abaqus学习笔记-abaqus与umat调用基本原理
参考: 1.http://ivt-abaqusdoc.ivt.ntnu.no:2080/v6.14/books/sub/default.htm 2.ABAQUS 用户材料子程序开发及应用-杨曼娟 3. ...
- Linux下5种IO模型的小结
概述 接触网络编程,我们时常会与各种与IO相关的概念打交道:同步(Synchronous).异步(ASynchronous).阻塞(blocking)和非阻塞(non-blocking).关于概念的区 ...
- Django之模型层（多表操作）
一.创建模型 1,一对多关系 一本书只有一个出版社,一个出版社可以出版多本书,从而书与出版社之间就构成一对多关系,书是‘多’的一方,出版社是‘一’的一方,我们在建立模型的时候,把外键写在‘多’的一方, ...
- 巨经典论文！推荐系统经典模型Wide & Deep
今天我们剖析的也是推荐领域的经典论文,叫做Wide & Deep Learning for Recommender Systems.它发表于2016年,作者是Google App Store的 ...
- iNeuOS工业互联网操作系统，三维（3D）模型在线编辑应用和实时数据统计（和值、均值、众数、方差、中位数等）
目 录 1. 概述... 1 2. 三维(3D)模型在线编辑与应用... 2 3. 实时数据统计... 4 1. 概述 此次,iNeuOS工业互联网操作系 ...
- 从0开始学angularjs－笔记01
一.angularjs简介 AngularJS 是一个为动态WEB应用设计的结构框架.它能让你使用HTML作为模板语言,通过扩展HTML的语法,让你能更清楚.简洁地构建你的应用组件.它的创新点在于,利 ...
- Java（JVM运行时）各种内存区域详解及扩展
本文整理于 Java内存与垃圾回收调优 Java 堆内存 从几个sample来学习Java堆,方法区,Java栈和本地方法栈 首先来一张图让我们理清楚java运行时状态: 诚然,如上图所示:java ...
- WebViewJavascriptBridge
上一篇文章介绍了通过UIWebView实现了OC与JS交互的可能性及实现的原理,并且简单的实现了一个小的示例DEMO,当然也有一部分遗留问题,使用原生实现过程比较繁琐,代码难以维护.这篇文章主要介绍下 ...
- Asp.Net MVC5入门学习系列④
原文:Asp.Net MVC5入门学习系列④ 添加Model且简单的使用EF 对于EF(EntityFramework)不了解的朋友可以去百度文科或者在园子里搜一些简资源看下,假如和我一样知道EF的概 ...
随机推荐
- card布局解决复杂操作的布局问题
一直不是很待见直接使用card布局,直到对于一些稍微复杂点的业务, 通过border布局和弹窗体的方式解决特别费劲之后,才想起了card布局, 发现card布局真是一个很好的解决办法. 那个使用起来很 ...
- Multiple dex files define Lcom/google/gson/internal/Streams$AppendableWriter$CurrentWrite;
开发中引入第三方 aar 时编译同过,运行时出现问题: Multiple dex files define Lcom/google/gson/internal/Streams$AppendableWr ...
- SpringBoot之配置
回顾 SpringBoot之基础 配置文件 ① 两种全局配置文件(文件名是固定的) 配置文件放在src/main/resources目录或者类路径/config下 application.proper ...
- python入门（十四）：面向对象（属性、方法、继承、多继承）
1.任何东西1)属性(特征:通常可以用数据来描述)(类变量和实例变量)2)可以做一些动作(方法) 类来管理对象的数据.属性:类变量和实例变量(私有变量)方法: 1)实例方法 2)类方法 ...
- 什么是云？Iaas,Paas和SaaS
周围的朋友听说我是做云相关的,总是爱问啥是云?别不是虚幻的概念吧.云计算当然不是虚幻的概念,“云”其实是互联网的一个隐喻,简单地说,云计算是通过Internet(“云”)交付计算服务——服务器.存储. ...
- JAVA字符串的处理
问题描述: 从键盘数入若干文字,最后输入的一行"end"代表结束标记. 统计该段文字中英文字母的个数 将其中的所有单词the全部改为a,输出结果 将该段文字所有的数字串找出来输出 ...
- Game Engine Architecture 7
[Game Engine Architecture 7] 1.SRT Transformations When a quaternion is combined with a translation ...
- GCC后端移植杂记
为GCC后端移植新的处理器架构的确是一件非常痛苦的事情. 拿到GCC源码时,发现它的规模比预想中庞大得多,但好在体系结构十分清晰,前端后端高度解耦.解耦的关键就在于内部的RTL中间代码(表达形式类似于 ...
- Linux - 其他命令
其他命令 目标 查找文件 find 软链接 ln 打包和压缩 tar 软件安装 apt-get 01. 查找文件 find 命令功能非常强大,通常用来在 特定的目录下 搜索 符合条件的文件 序号 命令 ...
- pycharm 配置svn
1.找到本地svn.exe 2.配置svn路径 3.输入svn路径,然后会进行密码校验 问题 本地的TortoiseSVN在bin目录没有发现svn.exe 原因: 安装TortoiseSVN的时候 ...