!参照《ANSYS结构动力分析与应用》王新敏编著中瞬态分析
!梁上移动质量的瞬态分析（移动荷载+位移耦合法）
FINISH$/CLEAR$/CONFIG,NRES,2000
/PREP7
LB=20$NE=400$NN=NE+1$NN1=NN+1$NE1=NE+1        !梁长，梁单元数，节点数，移动质量节点，总单元数
M1=1000$M2=3000$V=72/3.6                                            !质量作用梁的力，运行速度
DELTL=LB/NE                                     !梁划分单元的长度，在一个梁单元上运行时间
DELTT=DELTL/V                                 !在一个梁单元上运行时间
B=0.3$H=0.6                                      !梁宽，梁高
EM=2.1E11$AREA=B*H$IM=1/12*B*H*H*H            !弹性模量，面积，惯性矩
DENG=7800$GRA=9.81                                             !密度，加速度
K1=3E6$C1=1E4
F1=ACOS(-1)/2/LB/LB*SQRT(EM*IM/(AREA*DENG))   !自振频率（Hz）
ET,1,BEAM3
MP,EX,1,EM$MP,NUXY,1,0.3$MP,DENS,1,DENG
R,1,AREA,IM,H  
*DO,I,1,NN$N,I,(I-1)*DELTL$*ENDDO
*DO,I,1,NE$E,I,I+1$*ENDDO 

ET,2,MASS21,,,4$R,2,M1$R,3,M2                  !定义质量2D无转动
ET,3,14,,2$R,4,K1,C1                                      !定义刚度，弹簧
N,NN1,0,0$N,NN1+1,0,0
TYPE,2$REAL,2
E,NN1
TYPE,2$REAL,3
E,NN1+1
TYPE,3$REAL,4$E,NN1,NN1+1
D,1,UX,,,,,UY$D,NN,UY                          !施加梁约束和质量单元约束
D,NN1,UX,,,,,UY
F,NN1,FY,-M1*GRA                               !施加惯性力自重
F,NN1+1,FY,-M2*GRA
FINISH
!----------------------------------------------------------------------------------------------------------------
/SOLU
ANTYPE,4
OUTRES,ALL,ALL
AUTOTS,ON
!静力分析
TIMINT,OFF                           !关闭瞬态效应
TIME,0.005                            !定义时间（小时间间隔）
DELTIM,DELTT/10$KBC,0$SOLVE          !求解
TIME,0.1$SOLVE
!----------------------------------------------------------------------------------------------------------------
!质量在梁跨内的瞬态分析
TIMINT,ON                            !打开瞬态效应，设置输出控制
ALPHAD=0                            !瑞利阻尼
BETAD=0
NSUBST,1                                !1个载荷子步
DDELE,NN1,UY                        !删除质量单元的uy约束
DDELE,NN1+1,UY
*DO,I,2,NN
FDELE,I-1,ALL
TIME,(I-1)*DELTT+0.1                   !定义荷载步结束时间
D,NN1,UX,(I-1)*DELTL                   !对质量节点时间位移（根据速度施加）
D,NN1+1,UX,(I-1)*DELTL
CPDELE,ALL,ALL                        !删除所有既有耦合方程
CP,NEXT,UY,I,NN1                     !将质量节点与梁节点耦合uy自由度
SOLVE
SAVE
*ENDDO
!----------------------------------------------------------------------------------------------------------------
!质量移出梁跨的瞬态分析（自由振动）
CPDELE,ALL
FDELE,NN,ALL
D,NN1,UX,LB+1                         !将质量块移出梁跨并施加uy约束
D,NN1+1,UX,LB+1
D,NN1,UY
TIME,LB/V+5/F1                       !荷载结束时间增加5个周期（考虑自由振动）
SOLVE
SAVE
FINISH   
!----------------------------------------------------------------------------------------------------------------
/POST26
NC=NODE(LB/2,0,0)
NSOL,2,NC,U,Y
PLVAR,2
ANSYS采用位移建立车桥耦合模型，进行瞬态分析，车在桥跨内的响应没有问题，但是当车离开桥后的自由振动的响应为一条斜直线，怎么调都不对，请朋友们帮帮忙，看看自由振动的命令流哪里不对，或者说下自由振动的分析步骤也行，在下感激不尽！