单自由度体系在矩形脉冲下的瞬态分析

楼主对王光远翻译的《结构动力学》第二版P64的例题进行模拟。算例是质量为M的质量块通过刚度为K的弹簧固定，在T1到T1+T2时间内作用三角形脉冲，命令流如下。

FINISH

/CLEAR

K=10000

M=600

F=1000

/PREP7

ET,1,MASS21,,,4             !2-D mass without rotary inertia

ET,2,COMBIN14,,1           !UX DEGREE OF FREEDOM

R,1,M

R,2,K

N,1,1,0,0

TYPE,1

REAL,1

E,1

N,2,0,0,0

TYPE,2

REAL,2

EN,2,1,2

D,2,ALL

/SOLU

T1=0.05

T2=0.05

T3=5

ANTYPE,TRANS      !TRANSIDENT ANALYSIS

OUTRES,ALL,ALL

TRNOPT,FULL 

!LOAD STEP 1

TIMINF,OFF

TIME,T1

AUTOTS,ON

F,1,FX,0

LSWRITE,1

!LOAD STEP 2

TIMINF,ON

TIME,T1+T2

DELTIM,0.00001

KBC,1

F,1,FX,F

LSWRITE,2

!LOAD STEP 3

TIMINF,ON

TIME,T1+T2+T3

DELTIM,0.00001

KBC,0

F,1,FX,0

LSWRITE,3

LSSOLVE,1,3,1

FINISH

/POST26

NSOL,2,1,U,X

PLTIME,0,T1+T2+T3/5

PLVAR,2

但是模型的结果和理论解对应不了。

根据书上的讲解，最大反应比应该出现在在T1到T1+T2时间段内，最大位移为0.133。但ANSYS给出的结果是，最大位移为0.185，出现的时间远大于T1+T2。

不知道模型哪里错了，求纠正，谢谢。