请教一个umat如何返回节点坐标的问题

CWX_0806

各位高手：大家好！
      我现在编写一个程序，需要求模型中任意节点到（0,0,0）的距离，即r=sqrt（x8*x8+y8*y8+z8*z8），采用白老师的程序如下：
c ****************** 获得单元、节点、Part的编号与节点坐标 ********************
c------------------------------------------------------------------------------
      if (etype.eq.'brick') then
c get external element id
        iext=nelmntid(i,0)
c internal node ids
c   method 1  
c        ip=mxt(j)
c        n1=ix1(j)
c        n2=ix2(j)
c        n3=ix3(j)
c        n4=ix4(j)
c        n5=ix5(j)
c        n6=ix6(j)
c        n7=ix7(j)
c        n8=ix8(j)
c   method 2
        ip=ixs(1,i)
        n1=ixs(2,i)
        n2=ixs(3,i)
        n3=ixs(4,i)
        n4=ixs(5,i)
        n5=ixs(6,i)
        n6=ixs(7,i)
        n7=ixs(8,i)
        n8=ixs(9,i)
c external part and node ids
        ipext=lqfmiv(ip)
        next1=lqfinv(n1,1)
        next2=lqfinv(n2,1)
        next3=lqfinv(n3,1)
        next4=lqfinv(n4,1)
        next5=lqfinv(n5,1)
        next6=lqfinv(n6,1)
        next7=lqfinv(n7,1)
        next8=lqfinv(n8,1)
c spot check coordinates
        x8=x(1,n8)
        y8=x(2,n8)
        z8=x(3,n8)
结果是返回不了节点坐标，程序因错误而终止。请问怎么回事？谢谢

huihui123

你的这个问题解决了么？我也碰到这样的问题，有人说把etype.eq.'brick'改成etype.eq.solid''，我改了但是还是不行，我估计是subroutine umat41 (cm,eps,...)括号里面参数的问题。

haoyi891030

楼主您好，我也想请教一下，我是轴对称模型，需要确定当前积分点所在单元的节点坐标及其位移~

liyunliang_zi

书上的是针对dyna970的，i在头文件中出现，而971里，没有i， iext=nelmntid(i,0)里的i没有这个参量吧

liyunliang_zi

我也在考虑在调试时，监控单个单元内的应力情况，可是不知道怎么选择单个单元。
直接用文件的形式把所有单元应力输出，又太慢了