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的概 ...
随机推荐
- redis多实例和高可用
实验:多实例节点 10.220.5.65(三个实例) 高可用节点:10.220.5.67 创建目录,将三个实例文件都放在这个目录中 [root@ ~]# yum install redis -y [r ...
- centos下删除MYSQL 和重新安装MYSQL
centos下彻底删除MYSQL 和重新安装MYSQL 因centos系统自带的mysql版本比较低5.1,所以想卸载重新安装较新版本,下面是过程 1 删除Mysql yum remove mysql ...
- 《DOM Scripting》学习笔记-——第七章 动态创建html内容
本章内容: 1.动态创建html内容的“老”技巧:document.write()和innerHTML属性 2.DOM方法:createElement(),creatTextNode(),append ...
- 微信小程序---转发分享功能
1. 转发---onShareAppMessage 2.不带参数 //用户点击右上角分享 onShareAppMessage: function (res) { return { title: 'xx ...
- 使用MagickNet编辑图片
ImageMagick是一个免费的创建.编辑.合成图片的软件.它可以读取.转换.写入多种格式的图片.图片切割.颜色替换.各种效果的应用,图片的旋转.组合,文本,直线,多边形,椭圆,曲线 ...
- linux 安装 mysql
二进制安装mysql 1,下载安装包,使用国内站点速度会比较快,如清华站点https://mirrors.tuna.tsinghua.edu.cn 登陆linux系统后,使用wget 进行下载mysq ...
- ORACLE升级11g以上之前版本的wm_concat()函数失效
先执行: create or replace type string_sum_obj as object ( --聚合函数的实质就是一个对象 sum_string ), static function ...
- 从零开始学spring cloud(八) -------- Eureka 高可用机制
一.Eureka高可用机制介绍 Eureka服务器没有后端存储,但注册表中的服务实例都必须发送心跳以使其注册保持最新(因此可以在内存中完成). 客户端还有一个Eureka注册的内存缓存(因此,他们不必 ...
- 对palindrome的常用判断
判断String是否为palindrome:Two Pointers(left & right) 同时边扫边check 当前两边的char是否相同 code public boolean is ...
- redis学习笔记01 — 基本介绍、安装配置及常用命令
redis--NoSQL的一种 为了解决高并发.高可用.高扩展.大数据存储等一系列问题而产生的数据库解决方案,就是NoSQL NoSQL,非关系型数据库,全名:Not Only Sql,它不能代替关系 ...