FINISH
/CLEAR
/input,model-GAI,txt !--------load structure model-----
*create,readload
n=4095
*dim,ZK1,,n,1
*dim,ZK2,,n,1
*cfopen,ZK1,txt
*vread,ZK1(1,1),ZK1,txt
(F16.13)
*cfclos
*cfopen,ZK2,txt
*vread,ZK2(1,1),ZK2,txt
(F16.13)
*cfclos
*end
/wait,8
/input,readload
/solu
!m=6
!*dim,FN,,m,1 !---excitions node
!FN(1)=1010,1021,585,785,1032,1043 !------##################
!*do,i,1,m
!DDELE,FN(i),ux !-deleting Translational degrees of freedom
!*enddo
finish
!---------------Large mass method--------------------------
/prep7
*get,netype,etype,0,num,max !--max number of element type--
!*status,ntype
*get,nreal,rcon,0,num,max !--max number of element type--
!*status,nreal
ET,netype+1,MASS21
KEYOPT,netype+1,1,0
KEYOPT,netype+1,2,0
KEYOPT,netype+1,3,0
!----------------------------------------------------
LM=1 !----large mass #####################
/prep7
r,nreal+1,LM,LM,LM,
type,netype+1$real,nreal+1 !------##################
*do,i,1,m
e,FN(i)
*enddo
finish
allsel,all
/config,nres,100000
/SOLU
allsel
ANTYPE,TRANS
TRNOPT,FULL !---指定为完全瞬态分析
!DELTIM,dt,dt/10,dt*10
dmpr,0.05
!OUTRES,nsol,last !--输出基本项,每一项都输出
!OUTRES,esol,last
TIMINT,on !----------打开时间积分效应
!OUTPR,nsol,last
!OUTPR,esol,last
KBC,1 !----设置加载方式为线性加载
!------------------------------------------------
tt=40.95 !-----------地震持时
dt=0.01 !-----------时间步长
!-----------------------------定义和读入时程曲线
*do,iii,1,tt/dt
t=dt*iii
g=1
k=-1
D,1010,Ux,LM*g*ZK1(iii)
D,1021,Ux,LM*g*ZK1(iii)
D,585,Ux,LM*g*ZK1(iii)
D,785,Ux,LM*g*ZK1(iii)
D,1032,Ux,LM*g*ZK1(iii)
D,1043,Ux,LM*g*ZK1(iii)
D,1010,Uy,LM*g*ZK2(iii)
D,1021,Uy,LM*g*ZK2(iii)
D,585,Uy,LM*g*ZK2(iii)
D,785,Uy,LM*k*ZK2(iii)
D,1032,Uy,LM*k*ZK2(iii)
D,1043,Uy,LM*k*ZK2(iii)
time,t
solve
*enddo
FINISH
/CLEAR
/input,model-GAI,txt !--------load structure model-----
*create,readload
n=4095
*dim,ZK1,,n,1
*dim,ZK2,,n,1
*cfopen,ZK1,txt
*vread,ZK1(1,1),ZK1,txt
(F16.13)
*cfclos
*cfopen,ZK2,txt
*vread,ZK2(1,1),ZK2,txt
(F16.13)
*cfclos
*end
/wait,8
/input,readload
/solu
!m=6
!*dim,FN,,m,1 !---excitions node
!FN(1)=1010,1021,585,785,1032,1043 !------##################
!*do,i,1,m
!DDELE,FN(i),ux !-deleting Translational degrees of freedom
!*enddo
finish
!---------------Large mass method--------------------------
/prep7
*get,netype,etype,0,num,max !--max number of element type--
!*status,ntype
*get,nreal,rcon,0,num,max !--max number of element type--
!*status,nreal
ET,netype+1,MASS21
KEYOPT,netype+1,1,0
KEYOPT,netype+1,2,0
KEYOPT,netype+1,3,0
!----------------------------------------------------
LM=1 !----large mass #####################
/prep7
r,nreal+1,LM,LM,LM,
type,netype+1$real,nreal+1 !------##################
*do,i,1,m
e,FN(i)
*enddo
finish
allsel,all
/config,nres,100000
/SOLU
allsel
ANTYPE,TRANS
TRNOPT,FULL !---指定为完全瞬态分析
!DELTIM,dt,dt/10,dt*10
dmpr,0.05
!OUTRES,nsol,last !--输出基本项,每一项都输出
!OUTRES,esol,last
TIMINT,on !----------打开时间积分效应
!OUTPR,nsol,last
!OUTPR,esol,last
KBC,1 !----设置加载方式为线性加载
!------------------------------------------------
tt=40.95 !-----------地震持时
dt=0.01 !-----------时间步长
!-----------------------------定义和读入时程曲线
*do,iii,1,tt/dt
t=dt*iii
g=1
k=-1
D,1010,Ux,LM*g*ZK1(iii)
D,1021,Ux,LM*g*ZK1(iii)
D,585,Ux,LM*g*ZK1(iii)
D,785,Ux,LM*g*ZK1(iii)
D,1032,Ux,LM*g*ZK1(iii)
D,1043,Ux,LM*g*ZK1(iii)
D,1010,Uy,LM*g*ZK2(iii)
D,1021,Uy,LM*g*ZK2(iii)
D,585,Uy,LM*g*ZK2(iii)
D,785,Uy,LM*k*ZK2(iii)
D,1032,Uy,LM*k*ZK2(iii)
D,1043,Uy,LM*k*ZK2(iii)
time,t
solve
*enddo
FINISH
