abaqus如何在非等温情况下使用UMAT获得DSTRAN

嗨,,

我最近遇到了一个关于依赖温度的UMAT情况的大问题。

我使用UMAT进行非等温弹性，如下所示。这里没有显示DDSDDE计算部分。

你可以看到我在代码的开头输出了ETHERM(K1)，EELAS(K1)，STRAN(K1)，DSTRAN(K1)。

该UMAT还将ether和EELAS保存为状态变量数组。

作业设置是无加载和使用静态的一般步骤。但本质上它是一个自由膨胀，热膨胀系数为6e-6，温度从0到80，这只反映在UMAT中。Abaqus完成工作后，这些结果都是正确的，SDV(1-6)=EELAS=0, SDV(7-9)=ETHERM=4.8e-4。DSTRAN也是4.8e-4。

我的问题是abaqus如何知道这里的DSTRAN，因为我没有应用任何机械加载，也没有告诉abaqus这是一个热膨胀分析。我认为DSTRAN应该是由特定的剥离或热膨胀。

有人能帮我一下吗?我已经为此挣扎了一段时间。

子程序umat(应力，statev, ddsdde, sse, spd, scd, rpl，
1 ddsddt, drplde, drpldt, stran, dstran, time, dtime, temp, dtemp，
2 predef, dpred, cmname, ndi, nshr, ntens, nstatv, props, nprops，
3坐标，drot, pnewdt, celent, dfgrd0, dfgrd1, noel, npt, layer，
4 kspt, kstep, kinc)
C
包括“ABA_PARAM。公司的
C
字符* 8 CMNAME
C
维度应力(ntens)， statev (nstatv)， ddsdde (ntens, ntens)，
1 ddsddt (ntens)， drplde (ntens)， stran (ntens)， dstran (ntens)，
2 predef (1)， dpred (1)， props (nprops)，坐标(3)，drot (3,3)，
3 dfgrd0 (3,3)， dfgrd1 (3,3)

维eelas (6)， ether (6)， dtherm (6)， deldse (6,6)
C
参数(0 = 0。D0, = 1。D0, 2 = 2。D0, 3 = 3。D0, 6 = 6. D0)
C ----------------------------------------------------------------
C umat为线性变化的各向同性热弹性
C模-不能用于平面应力
C ----------------------------------------------------------------
C道具(1)- e (t0)
C道具(2)- nu (t0)
C道具(3)- t0
C道具(4)- e (t1)
C道具(5)- nu (t1)
C道具(6)- t1
C道具(7)-阿尔法
C props (8) - t_initial
C增量开始时的弹性特性
C

写(*)“进入UMAT”

k1 =1 ntens
写(*)ETHERM (K1) EELAS (K1) STRAN (K1) DSTRAN (K1)
最后做

C
C计算热膨胀
C

NDI, K1 = 1吗
ETHERM (K1) =道具(7)* (TEMP-PROPS (8))
DTHERM (K1) =道具(7)* DTEMP

最后做
nten, K1 = NDI + 1吗
ETHERM (K1) = 0
DTHERM (K1) = 0
最后做
C
C计算应力、弹性应变和热应变
C
k1 =1 ntens
k2 =1 ntens
压力(K2) =压力(K2) + DDSDDE (K1, K2) * (DSTRAN (K1) -DTHERM (K1)) + DELDSE (K1, K2) * (STRAN (K1) -ETHERM (K1))
最后做
ETHERM (K1) = ETHERM (K1) + DTHERM (K1)
EELAS (K1) = STRAN (K1) + DSTRAN (K1) -ETHERM (K1)
最后做

C
C存储弹性和热应变的状态变量数组
C
k1 =1 ntens
STATEV (K1) = EELAS (K1)
STATEV (K1 + nten) = ETHERM (K1)
最后做
返回
结束

谢谢你的帮助。

陈

电子邮件:chen1164@purdue.edu