马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
仿真单自由度有阻尼受冲击、简谐载荷下的位移、速度曲线可总是出问题啊。
图2为冲击载荷下位移、速度曲线。载荷大小30N,作用时间10秒。我觉得10秒之后,位移该是非零直线,速度是零才对嘛?
图3为简谐载荷下位移、速度、加速度曲线。但是蓝色的位移曲线好像不对哦?平衡点该是横坐标的周期曲线才对,是不是?
/clear
/filname, singal
/title, singal free vibration under fsin
/prep7
et,1,mass21
keyopt,0,0,4
et,2,combin14
keyopt,0,1,2
r,1,1
r,2,4,2
n,1,,,,
n,2,1,,,
type,1
real,1
esys,0
e,2
type,2
real,2
esys,0
e,1,2
M,2,UX
save
finish
!计算有阻尼频率
/SOLU
D,1,all
D,2,UY
ANTYPE,MODAL
MODOPT,QR damped,1,0,1,yes,1
MXPAND,1
LUMPM,yes
SOLVE
/post1
set,list
!冲击载荷作用下
/SOLU
ANTYPE,TRANS ! TRANSIENT DYNAMIC ANALYSIS
TRNOPT,REDUC ! LINEAR TRANSIENT
DELTIM,.05 ! INTEGRATION TIME STEP SIZE,1/60个周期
KBC,1 ! STEP BOUNDARY CONDITIONS
D,1,all
D,2,UY
F,2,FX,30
OUTPR,all,all
OUTRES,NSOL,2
SOLVE ! STATIC SOLUTION AT FIRST LOAD STEP
TIME,10 ! TIME TO INCLUDE ALMOST FOUR CYCLES
F,2,FX
SOLVE
FINISH
/POST26
FILE,,rdsp
NSOL,2,2,U,X,UX ! STORE UY DISPLACEMENTS OF NODE 2 AS UY
deriv,3,2,1,,VUX_2 ! 2节点的速度,VUX_2
/GRID,1 ! TURN GRID ON
/AXLAB,Y,DISP ! Y-AXIS LABEL AS DISP
PLVAR,2 ! DISPLAY VARIABLE 2 VS TIME
plvar,2,3 !时程作图显示
prvar,2,3 !时程列表显示
finish
! 简谐载荷作用下
NT=100
DT=0.05
*dim,p,array,NT
/INPUT,hsin,txt
/SOLU
ANTYPE,TRANS !时程分析
TRNOPT,FULL
KBC,0
NSUBST,1, , ,1
OUTPR,ALL,ALL
*do,i,1,NT
TIME,i*DT
AUTOTS,0
!加荷载
F,2,fx,p(i)
LSWRITE,i
*enddo
LSSOLVE,1,NT,1
SAVE
finish
/post26
nsol,2,2,u,x,UX_2
deriv,3,2,1,,VUX_2
deriv,4,3,1,,AUX_2
plvar,2,3
prvar,2,3 |