UMAT中递归变量的传递(存储)问题

junl06

我在编制UMAT的时候,需要用到递归,递归的次数由增量数(KINC)决定,但是最后子程序老是通不过,说出现The executable C:\ABAQUS\6.5-1\exec\standard.exe aborted with system error code 253. ....的错误,
子程序的框架如下:
SUBROUTINE UMAT(STRESS.....)
.......
DO k1=1,NTENS
   DSTRES(k1)=0.0
   STRESS(K1)=0.0
END DO
s1(kinc)=....
.......
end
recursive real function s1(kinc) result(ans1)
  if (kinc<=0)then
      ans1=0.0
  return
else if (kinc==1)then
     ans1=......
return
else
     ans1=8*s1(kinc-1)......
return
end if
     ans1=s1(kinc)
return
end      
请高手指点
还有就是上面的子程序在 UMAT中,若没有:
DO k1=1,NTENS
   DSTRES(k1)=0.0
   STRESS(K1)=0.0
END DO
则程序不能运行,但是这样是否使得每次调用UMAT的时候应力和应力增量都变成0了?前面的计算没有保存呢?不知道为啥?
s1是递归函数

junl06

1# junl06
呼唤版主回答问题啊，^_^

junl06

2# junl06
aogudusi版主帮忙解答下啊，谢了

junl06

3# junl06
咋没人回复的呢？高手跑哪去了

奔奔

你那个do

每次把应力都变成了0

怎么能对

错误肯定还有 我没仔细看

奔奔

建议通过 6通道号 把 结果输出看看  你自己每步迭代的结果

奔奔

我刚看了下

DSTRESS  是在input 文件中定义的

STRESS 只要定义初值就可以

junl06

7# 奔奔 应力增量 DSTRES怎么会从input给出呢？他应该实在UMAT中计算得到的啊，还有就是KINC怎么和递归的子程序进行传递呢？谢谢

junl06

8# junl06 是不是版主都找简单点的去回答了啊，咋这个一个都不回答呢？^_^