三维实体单元的UMAT子程序

亲爱的所有,
我很好在ABAQUS虽然第一次尝试子程序。
我正试图修改一些在线可用的代码，但我在这里发布的主要UMAT不起作用。
它显示。log文件中没有错误，但作业监视器显示“编译期间出现问题- D:\Abaqus工作目录\Trial_user_subroutine\try\Umat_mine.for”

请帮我一下，这件事很紧急。
谢谢大家，
附:代码以。txt格式附上
下面是代码

**
子程序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)

包括“ABA_PARAM。公司的

字符* 80 CMNAME

尺寸应力(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)

C本地数组
C ----------------------------------------------------------------
弹性应变
eplas -塑性应变
C——塑性流动的流向
C -------------------------------------------

尺寸eelas(6)、eplas(6)、flow(6)、hard (3)

参数(0 = 0。D0, = 1。D0, 2 = 2。D0,THREE=3.D0, SIX=6.D0,
& ENUMAX =。4999d0，牛顿=10,toler = 1.0d-6)
C
C ----------------------------------------------------------------
C为各向同性弹性和各向同性米塞斯塑性
C
不能用于平面应力
C ----------------------------------------------------------------
C道具(1)- e
C道具(2)- nu
C props(3..) -生成加固数据
C
很难求出屈服应力与塑性应变的曲线
C ----------------------------------------------------
C弹性性能
EMOD =道具(1)
ENU表示=道具(2)
EBULK3 = EMOD /(两人* ENU表示)
EG2 = EMOD / (1 + ENU表示)
如= EG2 / 2
= 3 *快
以拦= (EBULK3-EG2) / 3

C弹性刚度
k1 =1 ndi吗
k2 =1 ndi吗
DDSDDE (K1, K2) =拦
最后做
DDSDDE (K1, K1) = EG2 +拦
最后做
k1 = ndi +1吗
DDSDDE (K1, K1) =
最后做

C
例5:各向同性硬化塑性恢复弹性和塑性应变并向前旋转
C也能恢复等效塑性应变
调用rotsig (statev (1)， drot, eelas, 2, ndi, nshr)
调用rotsig (statev (ntens +1)， drot, eplas, 2, ndi, nshr)
EQPLAS = STATEV (1 + 2 * nten)
C
C
C
计算预测应力和弹性应变
k1 =1吗
k2 =1吗
压力(K2) =压力(K2) + DDSDDE (K1, K2) * DSTRAN (K1)
最后做
EELAS (K1) = EELAS (K1) + DSTRAN (K1)
最后做
C计算等效冯米塞斯应力
SMISES =(压力(1)压力(2))* * 2 +(压力(2)压力(3))* * 2
& +(压力(3)压力(1))* * 2
nten, K1 = NDI + 1吗
SMISES = SMISES + 6 *压力(K1) * * 2
最后做
SMISES =√SMISES /两个)

C
例5:各向同性硬化塑性从规定的硬化曲线得到屈服应力
NVALUE = NPROPS / 2 - 1
调用uhard (syiel0, hard, eqplas, eqplart, time, dtime, temp，
& DTEMP,诺埃尔,不扩散核武器条约》,层,KSPT, KSTEP, KINC, CMNAME, NSTATV,
》& STATEV NUMFIELDV,预定义,dpr NVALUE道具(3))


If (smiss . gt .(one + toler)* syiel0) then

C主动屈服分离静液和偏应力计算流动方向
SHYDRO =压力(压力(1)+(2)+(3)压力)/ 3
NDI, K1 = 1吗
压力流(K1) = ((K1) -SHYDRO) / SMISES
最后做
k1 = ndi +1吗
流(K1) = (K1) / SMISES压力
最后做
用牛顿迭代法求解等效冯米塞斯应力和等效塑性应变增量

SYIELD = SYIEL0
DEQPL = 0
牛顿等于1吗
RHS = SMISES-EG3 * DEQPL-SYIELD
DEQPL = DEQPL + RHS /(快+硬(1))
调用UHARD (SYIELD,硬,EQPLAS + DEQPL EQPLASRT,时间,DTIME,温度,
1 DTEMP,诺埃尔,不扩散核武器条约》,层,KSPT, KSTEP, KINC, CMNAME, NSTATV,
》2 STATEV NUMFIELDV预定义,dpr NVALUE道具(3))
如果(abs (rhs) .lt.toler * syiel0)转到10
最后做
C将警告消息写入.msg文件

例5:各向同性硬化塑性
写(7,2)牛顿
2 . format (//， 30x，'*** warning - plasticity algorithm not '，
1 '收敛后'，i3，'迭代')
10继续

NDI, K1 = 1吗
压力(K1) = (K1) * SYIELD + SHYDRO流动
EPLAS (K1) = EPLAS (K1) + 3 / 2 * (K1) * DEQPL流动
EELAS (K1) = EELAS (K1) / 2 * 3 (K1) * DEQPL流动
最后做
nten, K1 = NDI + 1吗
压力(K1) = (K1) * SYIELD流动
EPLAS (K1) = EPLAS (K1) + 3 (K1) * DEQPL *流
EELAS (K1) = EELAS (K1) * 3 (K1) * DEQPL流动
最后做
EQPLAS = EQPLAS + DEQPL

社民党= DEQPL * (SYIEL0 + SYIELD) / 2
EFFG = * SYIELD / SMISES
两个* EFFG EFFG2 =
EFFG3 = 3 / 2 * EFFG2
EFFLAM = (EBULK3-EFFG2) / 3
EFFHRD = *快硬(1)/(快+硬(1))-EFFG3
k1 =1 ndi吗
k2 =1 ndi吗
DDSDDE (K1, K2) = EFFLAM
最后做
DDSDDE (K1, K1) = EFFG2 + EFFLAM
最后做
k1 = ndi +1吗
DDSDDE (K1, K1) = EFFG
最后做
k1 =1吗
k2 =1吗
Ddsdde (k2, k1)= Ddsdde (k2, k1)+ effhrd *流量(k2)*流量(k1)
最后做
最后做
ENDIF
k1 =1吗
STATEV (K1) = EELAS (K1)
STATEV (K1 + nten) = EPLAS (K1)
最后做
STATEV (1 + 2 * nten) = EQPLAS
C
返回
结束

子例程UHARD (SYIELD硬、EQPLAS EQPLASRT,时间,DTIME,临时,
1 DTEMP,诺埃尔,不扩散核武器条约》,层,KSPT, KSTEP, KINC,
2 CMNAME NSTATV STATEV NUMFIELDV,
》3预定义,dpr, NVALUE表)
包括“ABA_PARAM。公司的
字符* 80 CMNAME


尺寸(3),STATEV (NSTATV)、时间(*),
》美元预定义(NUMFIELDV), dpr(*),道具(*)
C
C使用路德维希方程:sigma_yield = sigma_0 + k * eqplas ^ n
C d(SIGMA)/d(EQPLAS) = K*N*EQPLAS^(N-1)
C用户输入:

C材料未知道具1 2和3是在abaqus中定义的，而不是在代码中，所以不需要定义道具1 2 3
C
参数(0 = 0。D0,一个= 1. D0)
C
C计算
C
syfield = props (4) + props (5)* eqplas ** props (6)
困难(1)=道具(5)* EQPLAS道具(6)* * *(道具(6)——)
(2) = 0
(3) = 0
转到10
返回
结束