在psot1中提取节点解问题

16443

在post1中插入commands后提取出节点解，为什么会出现重复的节点呢？

1. !cmsel,s,nFace
   
2. cmsel,s,Edges
   
3. nsel,r,loc,z,0
   
4. Ncount=ndinqr(0,13)
   
5. *get,Nnum,NODE,,NUM,max
   
6. *dim,x,ARRAY ,Ncount
   
7. *dim,y,ARRAY ,Ncount
   
8. *dim,z,ARRAY ,Ncount
   
9. *dim,NI,ARRAY ,Ncount
   
10. *dim,Nsolx,ARRAY ,Ncount
    
11. *dim,Nsoly,ARRAY ,Ncount
    
12. *dim,Nsolz,ARRAY ,Ncount
    
13. *vget,x(1),NODE,1,LOC,X
    
14. *vget,y(1),NODE,1,LOC,Y
    
15. *vget,z(1),NODE,1,LOC,Z
    
16. *vget,Nsolx(1),NODE,1,U,X
    
17. *vget,Nsoly(1),NODE,1,U,Y
    
18. *vget,Nsolz(1),NODE,1,U,Z
    
19. 
    
20. !***********************************
    
21. *do,i,1,Ncount,1
    
22. NI(i)=node(x(i),y(i),z(i))
    
23. *enddo
    
24. !************************************
    
25. 
    
26. 
    
27. /rgb,index,100,100,100,0
    
28. /rgb,index,0,0,0,15
    
29. /show,png
    
30. 
    
31. PLNSOL,U,sum,1.0,1.0
    
32. 
    
33. *cfopen,Nnsol,dat
    
34. 
    
35. *vwrite,
    
36. ('********************Nodal location and dispacement ********************')
    
37. *vwrite,NI(1),x(1),y(1),z(1),Nsolx(1),Nsoly(1),Nsolz(1)
    
38. (F18.8,F18.8,F18.8,F18.8,E18.5,E18.5,E18.5)
    
39. 
    
40. *cfclos
    
41. 
    
42. *status,all
    
43. nplot,1

复制代码

16443

问题解决了。
cmsel,s,Edges
nsel,r,loc,z,0
Ncount=ndinqr(0,13)
prnsol,u,comp
*GET,nod_num,NODE,0,count  
*GET,nod_min,NODE,0,NUM,MIN
*dim,nsolz,,nod_num
*dim,nsoly,,nod_num
*dim,nsolx,,nod_num
*dim,NN,,nod_num
_nxt=nod_min   

*DO,i,1,nod_num,1
nsoly(i)=uy(_nxt)  
_nxt=ndnext(_nxt)   
*ENDDO

*DO,i,1,nod_num,1
nsolz(i)=uz(_nxt)
_nxt=ndnext(_nxt)   
*ENDDO

*DO,i,1,nod_num,1
nsolx(i)=ux(_nxt)  
_nxt=ndnext(_nxt)  
*ENDDO


!***********************************
_next=nod_min
*do,i,1,nod_num
NN(i)=_next
_next=ndnext(_next)
*enddo
!************************************

*cfopen,nsol,dat
*vwrite,
('***Node number********Node Solution************')
*vwrite,NN(1),nsolx(1),nsoly(1),nsolz(1)
(F8.1,E18.5,E18.5,E18.5)
*cfclos

nplot,1
*status,all

springzhang

看了程序，想问下
文件nsol.dat是在哪个目录下放着的????

pfsr

3# springzhang

可以自己在*cfopen中设定