橋梁移動荷載響應(yīng)分析示例

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        ! 劃分單元數(shù)
LBEAM=16        ! 梁長度
deltL=LBEAM/NELEM     ! 單元長度
V=120          ! 移動荷載速度
deltT=deltL/V*3.6      ! 單元時(shí)間
NNODE=NELEM+1     ! 節(jié)點(diǎn)數(shù)
fpw=40          ! 移動簡諧力頻率
fpc=100          ! 移動簡諧力振幅
*DO,i,0,NNODE,1        ! 定義節(jié)點(diǎn)
N,i,i*0.16,0        ! 節(jié)點(diǎn)間的間距為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        ! 每個(gè)荷載步的最終時(shí)間為xtim
KBC,1          ! 加載方式
NSUBST,5        ! 每個(gè)荷載步包含五個(gè)荷載子步
FDELE,all,all        ! 刪除所有的集中力荷載
NSEL,s,loc,x,ii*deltL      ! 選擇集中力作用節(jié)點(diǎn)
F,ALL,FY,-fpc*cos(fpw*ii*deltT)    ! 施加集中荷載
ALLS          ! 選擇所有
SOLVE          ! 求解
*ENDDO
FINISH
/POST26
NSOL,2,50,U,Y,UY_2   ! 定義第二個(gè)變量為UY_2,值為50號節(jié)點(diǎn)Y方向的位移
XVAR,1         ! 定義時(shí)間變量為坐標(biāo)橫軸
PLVAR,2         ! 輸出位移圖
DERIV,3,2,1,,,,,1     ! 將位移對時(shí)間求導(dǎo)，得到速度，存為變量3
PLVAR,3         ! 輸出速度圖
FINISH