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的概 ...
随机推荐
- linux重装rabbitmq的问题
一.卸载 [root@zabbix_server lib]# rpm -qa|grep rabbitmq rabbitmq-server--.noarch [root@zabbix_server li ...
- css第二天
二丶 1.字体属性font: 字体名称(font-family)字体大小(font-size):pc中通常,字体大小表示为12px,14px.移动设备中通常表示为0.57rem.字体粗细(font-w ...
- jQuery判断鼠标滚动方向
var scrolltop = new Array(); var index = 0; scrolltop[0] = 0; $(document).scroll(function(){ index + ...
- java、asp.net 通用分页码函数
<script type="text/javascript"> $(document).ready(function(){ ajaxGetPaging(1); }); ...
- 混合现实开发教程unity2017
共52节,MP4格式,英字,大小1GB 扫码时备注或说明中留下邮箱 付款后如未回复请至https://shop135452397.taobao.com/ 联系店主
- JVM学习01：内存结构
JVM学习01:内存结构 写在前面:本系列分享主要参考资料是 周志明老师的<深入理解Java虚拟机>第二版. 内存结构知识要点Xmind梳理 案例分析 分析1 package com.h ...
- webservice学习01：wsdl文档结构
webservice学习01:wsdl文档结构 wsdl文档结构 WSDL文档示例 <wsdl:definitions xmlns:xsd="http://www.w3.org/200 ...
- IIS7如何实现访问HTTP跳转到HTTPS访问
感谢原文作者,为方便后期查阅转载,原文链接:https://www.cnblogs.com/xiefengdaxia123/p/8542737.html 通常情况下我们是用的都是http的路径,对于h ...
- Python爬虫——Request模块
# 使用 Requests 发送网络请求# 1.导入 Requests 模块import requests# 2.尝试获取某个网页 # HTTP 请求类型r = requests.get('https ...
- dubbo知识体系
Dubbo负载均衡 Random随机/ RoundRobin轮询/ LeastActive最少Dubbo缓存Dubbo集群容错 Dubbo回调 服务分组: group Dubbo异步调用 <du ...