FINISH

/CLE

/PREP7

!---------------------------------單元材料-------------------------------

ET,1,SOLID65

MP,EX,1,1.851843E9 

MP,PRXY,1,0.15

MP,DENSE,1,1900  

MP,DAMP,1,0.023

TB,MKIN,1                              !砌體屈服準(zhǔn)則

TBTEMP,,STRAIN

TBDATA,,0.000548,0.0015,0.0020,0.0030,0.0040

TBTEMP,0

TBDATA,,1.01481E6,1.8883E6,2.1392E6,2.3600E6,2.2950E6

blc5,,,4,4,20

wpoffs,,,10

vsbw,all

vglue,all

VATT,1,,1                                   !M,R,T

MSHAPE,0,3D

ESIZE,0.8

VSWEEP,ALL

NSEL,S,LOC,Z,0

D,ALL,ALL                                   !底部約束

ALLSEL,ALL                                  !****///全選///****

FINISH

!********************************時(shí)程分析*********************************

HB=220

DT=0.02

*DIM,ACC,,HB,1                              !定義地震位移存儲(chǔ)數(shù)組      

*CREAT,DZ,MAC                               !創(chuàng)建宏

*VREAD,ACC(1),RAR-EL-H,TXT,,,HB,1

(1F9.5)

*END

DZ                                      !讀入X方向地震波文件(加速度，單位：0.1cm)

/SOLU

NCNV,2,1.0E7                                !不收斂時(shí)終止但不退出

ANTYPE,TRANS                                !有應(yīng)力法考慮非線性分析中的自重荷載

TRNOPT,FULL                                 !完全法求解輸出所有結(jié)果

AUTOTS,ON                                  !打開自動(dòng)時(shí)間步

OUTRES,ALL,ALL

TIMINT,OFF                                  !關(guān)閉時(shí)間積分效應(yīng)

TIME,0.001

ACEL,,,9.8                                  !施加重力

NSUBST,50

SOLVE                                      !*****///靜力求解///*****

KBC,1

TIMINT,ON

NSUBST,50

*DO,TM,1,50

  TIME,TM*DT

  ACEL,,0.001*ACC(TM)                       !單位由0.1cm化為m

  SOLVE                                     !*****///時(shí)程求解///*****

*ENDDO

FINISH

!********************************各截面反力及彎矩*********************************

/POST1

CSYS,0

ALLSEL

NSEL,S,LOC,Z,10，10+ 0.9                              !選擇截面點(diǎn)

ESLN,,1                                          !通過(guò)節(jié)點(diǎn)選擇單元

NSEL,R,LOC,Z,10                                   !從中再選擇截面點(diǎn)

SPOINT,,,,10                                      !定義力和力矩的求和位置點(diǎn)

*DIM,FANLI,,51                                   !定義壓力剪力存儲(chǔ)數(shù)組

*DIM,WANX,,51                                    !X彎矩

*DO,I,1,51

   SET,I

   FSUM                                          !對(duì)所選擇的節(jié)點(diǎn)力和力矩求和

   *GET,FANLI(I),FSUM,0,ITEM,FZ                  !提取各荷載步反力

   *GET,WANX(I),FSUM,0,ITEM,MX                   !X彎矩

*ENDDO

*VPLOT,,FANLI

*VPLOT,,WANX

FINISH

下圖是基底反力的結(jié)果：支反力隨時(shí)程逐漸歸零，但是找不到其中的原因，請(qǐng)各位

下圖是彎矩的結(jié)果：