首先是始点刚度法:

       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的更多相关文章

  1. abaqus学习笔记-abaqus与umat调用基本原理

    参考: 1.http://ivt-abaqusdoc.ivt.ntnu.no:2080/v6.14/books/sub/default.htm 2.ABAQUS 用户材料子程序开发及应用-杨曼娟 3. ...

  2. Linux下5种IO模型的小结

    概述 接触网络编程,我们时常会与各种与IO相关的概念打交道:同步(Synchronous).异步(ASynchronous).阻塞(blocking)和非阻塞(non-blocking).关于概念的区 ...

  3. Django之模型层(多表操作)

    一.创建模型 1,一对多关系 一本书只有一个出版社,一个出版社可以出版多本书,从而书与出版社之间就构成一对多关系,书是‘多’的一方,出版社是‘一’的一方,我们在建立模型的时候,把外键写在‘多’的一方, ...

  4. 巨经典论文!推荐系统经典模型Wide & Deep

    今天我们剖析的也是推荐领域的经典论文,叫做Wide & Deep Learning for Recommender Systems.它发表于2016年,作者是Google App Store的 ...

  5. iNeuOS工业互联网操作系统,三维(3D)模型在线编辑应用和实时数据统计(和值、均值、众数、方差、中位数等)

    目       录 1.      概述... 1 2.      三维(3D)模型在线编辑与应用... 2 3.      实时数据统计... 4 1.   概述 此次,iNeuOS工业互联网操作系 ...

  6. 从0开始学angularjs-笔记01

    一.angularjs简介 AngularJS 是一个为动态WEB应用设计的结构框架.它能让你使用HTML作为模板语言,通过扩展HTML的语法,让你能更清楚.简洁地构建你的应用组件.它的创新点在于,利 ...

  7. Java(JVM运行时)各种内存区域详解及扩展

    本文整理于  Java内存与垃圾回收调优 Java 堆内存 从几个sample来学习Java堆,方法区,Java栈和本地方法栈 首先来一张图让我们理清楚java运行时状态: 诚然,如上图所示:java ...

  8. WebViewJavascriptBridge

    上一篇文章介绍了通过UIWebView实现了OC与JS交互的可能性及实现的原理,并且简单的实现了一个小的示例DEMO,当然也有一部分遗留问题,使用原生实现过程比较繁琐,代码难以维护.这篇文章主要介绍下 ...

  9. Asp.Net MVC5入门学习系列④

    原文:Asp.Net MVC5入门学习系列④ 添加Model且简单的使用EF 对于EF(EntityFramework)不了解的朋友可以去百度文科或者在园子里搜一些简资源看下,假如和我一样知道EF的概 ...

随机推荐

  1. linux centos 磁盘清理

    执行df -h 与 du -sh / 所查询到的已用容量不对应 执行xfs_fsr来清理磁盘 参考 https://www.jianshu.com/p/0ded68808123

  2. Taro之百度地图显示定位点和信息

    由于没有后台数据就随机生成点来模拟了.具体代码如下. componentDidMount = () => { const { BMap, BMAP_STATUS_SUCCESS } = wind ...

  3. ES中的分词器

    基本概念: 全文搜索引擎会用某种算法对要建索引的文档进行分析, 从文档中提取出若干Token(词元), 这些算法称为Tokenizer(分词器), 这些Token会被进一步处理, 比如转成小写等, 这 ...

  4. IoU

    IoU #include <cstdio> #include <algorithm> #define re(i,a,b) for(int i=a;i<=b;i++) us ...

  5. Java框架spring 学习笔记(十五):操作MySQL数据库

    新建一个工程,添加对数据库的支持 下载mysql驱动包 mysql-connector-java-5.1.7-bin.jar,快捷键ctrl+alt+shift+s,添加jar包到工程 编写JdbcT ...

  6. Linux 进程同步和通信

    为了同步进程所以需要进程通信 管道(有名:文件形式存在,无名:仅限于父子进程间通信) 消息队列 信号量 共享存储 套接字(可用于不同机器)

  7. vue翻页器,包括上一页,下一页,跳转

    翻页组件 -- 子组件<template> <div class="pager-wrapper" ref="pager"> <di ...

  8. 【GitHub】的基本使用

    GitHub是一个常用的版本管理工具,之前安装了window版的git但是一直没使用过,今天尝试一下,去百度了一下使用方法: 引用https://www.cnblogs.com/paulwhw/p/9 ...

  9. yii2 模板twig中使用GridView::widget

    yii框架twig模板中使用GridView::widget 需要这样写: {{ grid_view_widget({ 'dataProvider': provider, 'columns': [ ' ...

  10. Codeforces Round #491 (Div. 2)

    Codeforces Round #491 (Div. 2) https://codeforces.com/contest/991 A #include<bits/stdc++.h> us ...