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

嗨,,

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

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

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

该UMAT还将以太坊和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 (stress, statev, ddsde, sse, spd, scd, rpl，
1 . dsddt, dpldt, drpldt, stran, dstran, time, dtime, temp, dtemp，
predef, dpred, cmname, ndi, nshr, ntens, nstatv, props, nprops，
3、drot、pnewdt、celent、dfgrd0、dfgrd1、noel、npt、layer、
4 . kspt, kstep, king)
C
包括“ABA_PARAM。公司的
C
字符* 8 CMNAME
C
尺寸应力(ntens)，状态(nstatv)， ddsdde (ntens, ntens)，
1 ddsdt (ntens)， drplde (ntens)， stran (ntens)， dstran (ntens)，
2 predef (1)， dpred (1)， props (nprops)， coordinates (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值呈线性变化
C模量-不能用于平面应力
C ----------------------------------------------------------------
C道具(1)- e (10)
C props (2) - nu (0)
C道具(3)- 0
C道具(4)- e (t1)
C道具(5)- nu (t1)
C道具(6)- t1
C道具(7)- α
C props (8) - t_initial
C增量开始时的弹性特性
C

写(*)“进入UMAT”

k1 =1吗
写(*)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吗
k2 =1吗
压力(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吗
STATEV (K1) = EELAS (K1)
STATEV (K1 + nten) = ETHERM (K1)
最后做
返回
结束

谢谢你的帮助。

陈

电子邮件:chen1164@purdue.edu