你在这里
粘弹性UMAT
星期二，2012-07-03 01:20 -supandanane
大家好，
我正在寻找基于麦克斯韦或开尔文元素的粘弹性UMAT代码。如果有人能帮我的话。
问候
»
- supandanane的博客
- 登录或注册发表评manbetx体育论
- 15373年读
评论
粘弹性UMAT: Maxwell和Kelvin
麦克斯韦模型可以在我的博士论文中找到
玻璃体SiO2的镦粗与粘弹性:实验、解释与模拟
http://opus.kobv.de/tuberlin/volltexte/2006/1179/
的
开尔文-福格特模型可以用类似的方法推导出来。代码是
在下面。我把推导过程写成pdf文档;但这
这里的东西不允许我上传。公开你的电邮地址，
然后我就可以寄给你了。
祝你好运
弗兰克
--------------------------------------------------------
子程序UMAT(压力、STATEV DDSDDE SSE,社民党,谈到
1 RPL DDSDDT DRPLDE DRPLDT,
2 STRAN DSTRAN时间、DTIME TEMP, DTEMP, PREDEF, dpr, CMNAME,
3 NDI NSHR nten, NSTATV道具、NPROPS,坐标,DROT, PNEWDT,
4 CELENT DFGRD0 DFGRD1,诺埃尔,不扩散核武器条约》,层,KSPT, KSTEP KINC)
C
包括“ABA_PARAM。公司的
C
字符* 80 CMNAME
维应力(nten), STATEV (NSTATV),
1 DDSDDE (nten, nten),
2 DDSDDT nten, DRPLDE (nten),
3 STRAN nten, DSTRAN (nten), (2), PREDEF (1), dpr (1),
4道具(NPROPS)坐标(3),DROT(3、3),DFGRD0(3、3),DFGRD1(3、3)
C
C
整数i, j
REAL Y n
REAL，是，是
实a b c d f g
C
Real, dimension (ntens, ntens):: ddsde
Real, dimension (ntens, ntens):: DDSDS
实，维度(NTENS):: SIGold
C
Y = props (1)
n = PROPS(2)
eta =道具(3)!剪切粘度
C
= Y*n/((1.0d0+n)*(1.0d0-2.0d0*n))
mu = Y/(2.0d0*(1.0d0+n))
thetalambda = 3.0d0*eta/Y
= eta/
C
A = lambda*(1.0d0 + thetalambda/DTIME) +
1 2.0d0*mu*(1.0d0 + thetamu/DTIME)
B = lambda*(1.0d0 + thetalambda/DTIME)
C = lambda + 2.0d0*mu
D =
F = 2.0d0*mu*(1.0d0 + thetamu/DTIME)
G = 2.0d0*mu
C
i = 1,NTENS
j = 1,NTENS
DDSDDE(i,j) = 0.0d0
DDSDE(i,j) = 0.0d0
DDSDS(i,j) = 0.0d0
最后做
最后做
C
C定义与正应力相关的张量系数
C
i = 1,NDI
j = 1,NDI
DDSDDE(i,j) = B
DDSDE(i,j) = D
DDSDS(i,j) = 0.0d0
最后做
DDSDDE(i,i) = A
DDSDE(i,i) = C
DDSDS(i,i) = -1.0d0
最后做
C
C定义与剪应力相关的张量系数
C
i = NDI+1,NTENS
DDSDDE(i,i) = F/2.0d0
DDSDE(i,i) = G/2.0d0
DDSDS(i,i) = - 1.0d0
最后做
C
在数组中存储当前的应力状态
C
i= 1,NTENS
SIGold(i) = STRESS(i)
最后做
C
C更新应力
C
i = 1,NTENS
j = 1,NTENS
应力(i) =应力(i) + DDSDDE(i,j)*DSTRAN(j) +
1 DDSDE(i,j)*STRAN(j) +
2 DDSDS (i, j) * SIGold (j)
最后做
最后做
C
返回
结束
------------------------------------------
-鲁尔大学
波鸿
德国
粘弹性UMAT
里希特博士
我知道这个论坛帖子已经有很多年了，但是我非常有兴趣看看你的开尔文-福格特模型的PDF推导。我是加州州立理工大学的一名研究生，您的博士工作对我的论文以及我的同行的工作都是非常重要的。
还能从你那里拿到这个吗?
谢谢你！
奥斯汀C。
阿库明[在]calpoly[点]edu
你好弗兰克!谢谢你
你好弗兰克!谢谢你在这里分享你的UMAT代码。我刚开始研究粘弹性材料，需要在UMAT中实现非线性各向异性粘弹性材料。您的代码似乎是一个很好的开始。你上面提到的pdf写起来进行推导。我想知道你是否可以分享它和你的博士论文，上面的ilnk不管用?我的邮箱是sergkuznet@hotmail.com
谢谢你！
UMAT表示粘弹性
我的论文已经开始了
https://depositonce.tu-berlin.de/handle/11303/1574
开尔文-福格特模型的UMAT实现:
http://万博manbetx平台m.limpotrade.com/node/19884
读也
https://polymerfem.com/forum/finite-element-modeling/constitutive-models..。
祝你好运
弗兰克
好谢谢!
好谢谢!