I ran a transient analysis and want to see the average velocity for a named nodal selection over the full time of the transient.


The below commands will produced the desired data.

*get,dset,active,,set,nset
*do,jj,1,dset
*if,jj,eq,1,then
set,first
cmsel,s,name ! change "name" to the named selection's name
*get,n_mx,node,,num,max
*get,n_count,node,,count
*dim,ndat,array,n_mx,4 ! nodal data with col 4 mask, col 1/2/3 vx/vy/vz
*vget,ndat(1,4),node,1,nsel
*dim,avdat,array,dset,4 ! substep data with col 4 time, col 1/2/3 vx avg/vy avg/vz avg
*dim,tavdat,table,dset,3 ! table copy of avdat for line plot with col 0 time and 1/2/3 same as avdat
*else
set,next
*endif
*get,avdat(jj,4),active,,set,time
!
nzz=0
*do,qq,1,n_count
nzz=ndnext(nzz)
*get,ndat(qq,1),node,qq,v,x
*get,ndat(qq,2),node,qq,v,y
*get,ndat(qq,3),node,qq,v,z
*enddo
!
!*vmask,ndat(1,4)
!*vget,ndat(1,1),node,1,v,x
!*vmask,ndat(1,4)
!*vget,ndat(1,2),node,1,v,y
!*vmask,ndat(1,4)
!*vget,ndat(1,3),node,1,v,z
!
*vmask,ndat(1,4)
*vscfun,avdat(jj,1),mean,ndat(1,1)
*vmask,ndat(1,4)
*vscfun,avdat(jj,2),mean,ndat(1,2)
*vmask,ndat(1,4)
*vscfun,avdat(jj,3),mean,ndat(1,3)
*vfill,ndat(1,1),ramp,0,0
*vfill,ndat(1,2),ramp,0,0
*vfill,ndat(1,3),ramp,0,0
*enddo
!
nzz=ndat(1)=
n_mx=
dset=
jj=
qq=
*vfun,tavdat(1,0),copy,avdat(1,4)
*vfun,tavdat(1,1),copy,avdat(1,1)
*vfun,tavdat(1,2),copy,avdat(1,2)
*vfun,tavdat(1,3),copy,avdat(1,3)
avdat(1)=
parsav,all
*vplo,tavdat(1,0),tavdat(1,1)
*vplo,tavdat(1,0),tavdat(1,2)
*vplo,tavdat(1,0),tavdat(1,3)





Show Form
No comments yet. Be the first to add a comment!