每人管,只好自己想办法,用了一下午看了下帮助,写了个宏,暂时可行。麻烦大虾看看,有什么毛病或者可以简化的地方没?
这个过程就是求解个质点的加速过程,有阻力,阻力和速度平方成比例,前边准备操作都没整成命令流
/PREP7
TM_START=1E-6
TM_END=1
TM_INCR=0.05
VELO=0
FINISH
/SOLU
ANTYPE,4,NEW
TRNOPT,FULL
LUMPM,0
D,1,UY,0
*DO,TM,TM_START,TM_END,TM_INCR
TIME,TM
F,1,FX,10-VELO*VELO*10
SOLVE
FINISH
/POST26
NUMVAR,200
NSOL,2,1,U,X,UX
DERIV,3,2,1,,VEL1UY
*GET,VELO,VARI,3,RTIME,TM
/SOLU
ANTYPE,,REST
*ENDDO |