橋梁移動荷載響應分析示例

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