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                              !砌体屈服准则

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

!********************************时程分析*********************************

HB=220

DT=0.02

*DIM,ACC,,HB,1                              !定义地震位移存储数组      

*CREAT,DZ,MAC                               !创建宏

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

(1F9.5)

*END

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

/SOLU

NCNV,2,1.0E7                                !不收敛时终止但不退出

ANTYPE,TRANS                                !有应力法考虑非线性分析中的自重荷载

TRNOPT,FULL                                 !完全法求解输出所有结果

AUTOTS,ON                                  !打开自动时间步

OUTRES,ALL,ALL

TIMINT,OFF                                  !关闭时间积分效应

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                                     !*****///时程求解///*****

*ENDDO

FINISH

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

/POST1

CSYS,0

ALLSEL

NSEL,S,LOC,Z,10，10+ 0.9                              !选择截面点

ESLN,,1                                          !通过节点选择单元

NSEL,R,LOC,Z,10                                   !从中再选择截面点

SPOINT,,,,10                                      !定义力和力矩的求和位置点

*DIM,FANLI,,51                                   !定义压力剪力存储数组

*DIM,WANX,,51                                    !X弯矩

*DO,I,1,51

   SET,I

   FSUM                                          !对所选择的节点力和力矩求和

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

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

*ENDDO

*VPLOT,,FANLI

*VPLOT,,WANX

FINISH

下图是基底反力的结果：支反力随时程逐渐归零，但是找不到其中的原因，请各位

下图是弯矩的结果：