FINISH
/CLEAR,NOSTART
/PREP7 ! 建模
ET,1,BEAM3 ! 梁单元
MP,EX,1,207E9 ! 弹性模量
MP,DENS,1,2600 ! 密度
MP,NUXY,1,0.3 ! 泊松比
R,1,0.1,0.0001/12,0.1 ! 横截面积、惯性距、截面高
NELEM=100 ! 划分单元数
LBEAM=16 ! 梁长度
deltL=LBEAM/NELEM ! 单元长度
V=120 ! 移动荷载速度
deltT=deltL/V*3.6 ! 单元时间
NNODE=NELEM+1 ! 节点数
fpw=40 ! 移动简谐力频率
fpc=100 ! 移动简谐力振幅
*DO,i,0,NNODE,1 ! 定义节点
N,i,i*0.16,0 ! 节点间的间距为0.3
*ENDDO
*DO,i,1,NELEM ! 定义单元
E,i,i+1
* ENDDO
D,1,ALL ! 定义约束
D,NNODE,ALL
FINISH
/ESHAPE,1
/VIEW,1,1,1,1
EPLOT
/SOLU
ANTYPE,trans ! 分析类型
OUTREA,all,all ! 输出控制
*DO,ii,1,NNODE,1
time,ii*deltT ! 每个荷载步的最终时间为xtim
KBC,1 ! 加载方式
NSUBST,5 ! 每个荷载步包含五个荷载子步
FDELE,all,all ! 删除所有的集中力荷载
NSEL,s,loc,x,ii*deltL ! 选择集中力作用节点
F,ALL,FY,-fpc*cos(fpw*ii*deltT) ! 施加集中荷载
ALLS ! 选择所有
SOLVE ! 求解
*ENDDO
FINISH
/POST26
NSOL,2,50,U,Y,UY_2 ! 定义第二个变量为UY_2,值为50号节点Y方向的位移
XVAR,1 ! 定义时间变量为坐标横轴
PLVAR,2 ! 输出位移图
DERIV,3,2,1,,,,,1 ! 将位移对时间求导，得到速度，存为变量3
PLVAR,3 ! 输出速度图
FINISH
