求助！！三维环焊缝单元排序
:( :( :'(我给三维环焊缝（多道焊-单元类型为2）的单元进行排序，用的是以下程序，但是结果提示第一道焊的第一个单元杀不死，提示说是此单元没有排列进去。高手帮我看看这个程序怎么改动一下？还有提示：specified rang of 0 to 0 is not permitted!这么修改阿请指教。。。:'( :'(CSYS,0
ESEL,S,TYPE,,2
*GET,EMAX,ELEM,,NUM,MAX
ESEL,S,TYPE,,2
*GET,NSE,ELEM,,COUNT
*DIM,NE,,NSE
*DIM,NEX,,NSE
*DIM,NEY,,NSE
*DIM,NEZ,,NSE
*DIM,NER,,NSE
*DIM,NER1,,NSE
*DIM,NEA,,NSE
*DIM,NEORDER,,NSE
MINE=0
*DO,I1,1,NSE
ESEL,U,ELEM,,MINE
*GET,NSEL,ELEM,,COUNT
II=0
*DO,I,1,EMAX
*IF,ESEL(I),EQ,1,THEN
II=II+1
NE(II)=I
*ENDIF
*ENDDO
*DO,I,1,NSEL
*GET,NEX(I),ELEM,NE(I),CENT,X
*GET,NEY(I),ELEM,NE(I),CENT,Y
*GET,NEZ(I),ELEM,NE(I),CENT,Z
*ENDDO
*DO,I,1,NSEL
NER(I)=SQRT(NEX(I)**2+NEZ(I)**2)
NER1(I)=SQRT((NER(I)-r2)**2+NEY(I)**2)
*IF,NEZ(I),LE,1E-6,THEN
NEA(I)=ACOS(NEX(I)/NER(I))
*ELSE
NEA(I)=ACOS(NEX(I)/NER(I))+PI
*ENDIF
*ENDDO
MINR=1E20
MINY=1E20
MINA=1E20
*DO,I,1,NSEL
*IF,NER1(I),LT,MINR,THEN
MINR=NER1(I)
MINY=NEY(I)
MINA=NEA(I)
MINE=NE(I)
*ELSE
*IF,NER1(I),LT,MINR+0.002,AND,NER1(I),GE,MINR,THEN
*IF,NEY(I),LT,MINY,THEN
MINR=NER1(I)
MINY=NEY(I)
MINA=NEA(I)
MINE=NE(I)
*ELSE
*IF,NEY(I),LT,MINY+0.002,AND,NEY(I),GT,MINY,THEN
*IF,NEA(I),LT,MINA,THEN
MINR=NER1(I)
MINY=NEY(I)
MINA=NEA(I)
MINE=NE(I)
*ENDIF
*ENDIF
*ENDIF
*ENDIF
*ENDIF
*ENDDO
NEORDER(I1)=MINE
*ENDDO
*DO,I,1,NSE
EKILL,NEORDER(I)
ESEL,S,LIVE
同命相怜呀
我也遇到了类似的问题*** ERROR *** CP = 102.641 TIME= 15:35:00
Element 0 is not defined.
*** WARNING *** CP = 102.672 TIME= 15:35:00
Specified range of 0 to 0 is not permitted.
The D command is ignored.
*** WARNING *** CP = 123.750 TIME= 15:36:28
The mesh of area 3 contains PLANE13 triangles, which are much too stiff
in bending.Use quadratic (6- or 8-node) elements if possible.
*** WARNING *** CP = 134.891 TIME= 15:37:39
There are 1 small equation solver pivot terms.May be a transitional
effect caused by the Newton-Raphson logic.
哪位高人帮忙指点一下谢谢了 哪位高人帮忙指点一下谢谢了 必须结合模型才能对单元进行排序阿 单元的选取是不是有问题啊 新手 不会的飘过 本帖最后由 dollfish000 于 2017-8-1 16:54 编辑
CSYS,0 改为CSYS,1 其实我也是第一次做环焊缝，不会，交流加260013078
页:
[1]