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的概 ...
随机推荐
- [Sw] Swoole-4.2.9 可以尝试愉快应用 Swoole 协程
大家知道 Swoole 提供了方便于服务器.网络编程的模式,简化了多进程编程. 这直接让 PHP 的运行很容易变成常驻内存的 Server 程序,执行效率上有了数倍的提升. 但是这一切还没有让人足够兴 ...
- 基于STM32F429+HAL库编写的定时器主从门控模式级联输出固定个数PWM脉冲的程序
硬件设备 42步进电机,步进电机驱动器,正点原子F429开发板 开发软件 keil5,Cube 综述 一般要精准的控制电机,就要控制单片机的引脚输出指定个数的PWM波,有多种可实现的方法 ...
- yii2.0 url美化-apache服务器
//配置内容 'urlManager' => [ 'enablePrettyUrl' => true, 'enableStrictParsing' => false, //不启用严格 ...
- angluar1+ionic详情页返回在原来的位置（缓存数据和页面高度）
因为是老项目,近期开发遇到了个需求就是从详情页按返回按钮要求返回到原来列表的页面位置,刚开始准备用的cache:true,但是存在大大的问题就是新增和编辑后返回数据都不是最新的,无法重新刷新页面rel ...
- TCP报文格式
转载自https://blog.csdn.net/mary19920410/article/details/58030147 1.TCP报文格式 TCP报头中的源端口号和目的端口号同IP数据报中的源I ...
- Celery 1
Celery是一个用Python开发的异步的分布式任务调度模块 Celery有以下优点: 简单:一但熟悉了celery的工作流程后,配置和使用还是比较简单的 高可用:当任务执行失败或执行过程中发生连接 ...
- Codeforces Round #485 (Div. 2)
Codeforces Round #485 (Div. 2) https://codeforces.com/contest/987 A #include<bits/stdc++.h> us ...
- java实验环境搭建，eclise下载与使用
一.1.官方下载地址:http://www.oracle.com/technetwork/java/javase/downloads/index.html 下载后安装,默认安装即可. 2.在 wind ...
- centos7安装zabbix server
1.参照下列网址方法,打开端口:80,3306,443,22,10050,10051(可能实际不需要打开这么多) https://www.cnblogs.com/lw-2019forlinuxpyth ...
- ubuntu 中 vim 的使用
安装 sudo apt install vim vim file_name #创建或者打开文件 vim file_name 定位到文件开头 vim file_name + 定位到文件末尾 vim f ...