桥梁移动荷载响应分析示例

桥梁移动荷载响应分析示例

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