提取节点速度加速度（附命令流）

junzhukeer

/solu   
antype,trans
TRNOPT,full
nlgeom,on
imint,on
ALPHAD,AL
BETAD,BE  
*afun,deg   
*do,j,dt1,dt3,dt
!*do,jj,1,n_total   
*SET,dhh,(hh(jj,j/dt+1)-hh(jj,j/dt))/dt   
*SET,dpp,(pp(jj,j/dt+1)-pp(jj,j/dt))/dt   
*SET,beta(jj,1),ATAN(dhh/U)
*SET,alfa(jj,1),10-beta(jj,1)  
*SET,c_l(jj,1),cldm(alfa(jj,1),1)   
*SET,c_d(jj,1),cldm(alfa(jj,1),2)   
*SET,F_y(jj,1),FK*(c_l(jj,1)*COS(beta(jj,1))-c_d(jj,1)*sin(beta(jj,1)))
*SET,F_z(jj,1),-FK*(c_l(jj,1)*sin(beta(jj,1))+c_d(jj,1)*cos(beta(jj,1)))
F,jj,fy,F_y(jj,1)   
F,jj,fz,F_z(jj,1)   
allsel  
time,j   
DELTIM,0.01
kbc,0   
outres,all,all  
solve
FINISH
*get,hh(jj,j/dt+2),node,jj,U,y
*get,pp(jj,j/dt+2),node,jj,U,z
!*enddo
*enddo  
save
麻烦各位大侠指教一下，我的GET命令怎么好像没起作用，因为z方向没有位移和速度，都是0！

[ 本帖最后由 junzhukeer 于 2009-4-9 22:21 编辑 ]

16443

1. 数组hh和pp没有定义
2、jj是个未知数

junzhukeer

这个定义了，我是在前面定义的 ，所有的 数组和变量都定义了，就是不知道是不是循环和提取速度的时候哪里有问题，Z方向没有没有位移和速度，敬请指教