/prep7
block,0,350,0,300,0,190 !建垫板体一
block,350,550,0,300,0,190 !建垫板体二
vgen,6,2,,,200,, !生成体三四五六七
vgen,2,1,,,1550,, !生成体八

block,350,550,300,1900,0,190 !体九左侧构造柱
vgen,2,9,,,200,, !体十
block,580,720,300,1900,30,160 !体十一（左芯柱）
vsbv,10,11,,,keep !外壳也就是砌块是体12,芯柱仍是体11
vgen,4,12,,,200,, !中间2块空心砌块体10，13，14
vgen,2,11,,,600,, !右侧芯柱体15
vgen,2,9,,,1000,, !右侧构造柱体16

k,153,380,200,30 !纵筋和箍筋
k,154,380,2000,30
l,153,154

k,155,520,200,30
k,156,520,2000,30
l,155,156

k,157,520,200,160
k,158,520,2000,160
l,157,158

k,159,380,200,160
k,160,380,2000,160
l,159,160

l,153,155
l,155,157
l,157,159
l,159,153

lsel,s,loc,x,380,520
lsel,r,loc,y,200
lplot
lgen,19,all,,,,100 !复制一圈箍筋
allsel,all



k,233,1380,200,30 !左侧 纵筋和箍筋
k,234,1380,2000,30
l,233,234

k,235,1520,200,30
k,236,1520,2000,30
l,235,236

k,237,1520,200,160
k,238,1520,2000,160
l,237,238

k,239,1380,200,160
k,240,1380,2000,160
l,239,240

l,233,235
l,235,237
l,237,239
l,239,233
lsel,s,loc,x,1380,1520
lsel,r,loc,y,200
lplot
lgen,19,all,,,,100
allsel,all

k,313,650,200,95
k,314,650,2000,95
l,313,314 !芯柱中钢筋

k,315,1250,200,95
k,316,1250,2000,95
l,315,316


block,350,550,1900,2100,0,190 !顶梁六块17-22
block,550,750,1900,2100,0,190
block,750,950,1900,2100,0,190
block,950,1150,1900,2100,0,190
block,1150,1350,1900,2100,0,190
block,1350,1550,1900,2100,0,190

block,1550,1600,1900,2100,0,190 !垫板23
几何模型建完
save

划分单元网格
ET,1,SOLID65
ET,2,LINK8
ET,3,SOLID45

R,1 !大部分砼(顶梁外)和刚性垫块为不配筋的R1
R,2,3,0.01,90,0,3,0.01, !顶梁砼对应于配筋的R2
RMORE,0,0,3,0.001,90,0,
R,3,28.3 !直径为6的钢筋(箍筋)对应于R3
R,4,113 !直径为12的钢筋(纵筋和水平拉筋)对应于R4


MP,EX,1,3.0e4 !MAT1混凝土C30基本材料属性
MP,PRXY,1,0.2
TB,CONC,1,1,9, !C30破坏准则
TBDATA,,0.8,1,2.47,-1,,
tb,miso,1,1,12 !C30混凝土本构关系（多线性等向强化）
tbpt,,0,0
tbpt,,0.0002,6
tbpt,,0.0004,10.8
tbpt,,0.0006,15.3
tbpt,,0.0008,19.2
tbpt,,0.001,22.5
tbpt,,0.0012,25.2
tbpt,,0.0014,27.3
tbpt,,0.0016,28.8
tbpt,,0.0018,29.7
tbpt,,0.0020,30
tbpt,,0.0040,24

MP,EX,2,17600 !MAT2砌体基本材料属性
MP,PRXY,2,0.15
TB,CONC,2,1,9, !砌体破坏准则
TBDATA,,0.8,1,2.47,-1,,
TB,MISO,2,1,8 !砌体本构关系
tbpt,,0,0
tbpt,,0.0001,1.76
tbpt,,0.0004,6.49
tbpt,,0.0008,11.55
tbpt,,0.0012,15.15
tbpt,,0.0016,17.32
tbpt,,0.002,18.04
tbpt,,0.005,9.02
MP,EX,3,2e5 !MAT3 HPB235
MP,PRXY,3,0.3
TB,BKIN,3,1,2,1
TBDATA,,210,20000,,,,

MP,EX,4,2e5 !MAT4 HPB400
MP,PRXY,4,0.3
TB,BKIN,4,1,2,1
TBDATA,,360,20000,,,,

MP,EX,5,2e5 !MAT5 刚性垫板
MP,PRXY,5,0.3
TB,BKIN,5,1,2,1
TBDATA,,800,20000,,,,

划分纵向钢筋
allsel,all
lsel,s,length,,1800
LPLOT
LATT,4,4,2
lesize,all,100
LMESH,all,100

划分箍筋
allsel,all
lsel,s,length,,140
lsel,a,length,,130
lsel,u,loc,x,550,1350
lplot
LATT,3,3,2
lesize,all,,,2
LMESH,all,,,2
gplot

allsel,all
lesize,98,,,4,-7/3 !左侧构造柱
lesize,107,,,4,-13/6
lesize,103,,,4,-7/3
lesize,108,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,1,1
lesize,102,,,16
vsweep,9,52,51


lesize,218,,,4,-7/3 !右侧构造柱
lesize,227,,,4,-13/6
lesize,223,,,4,-7/3
lesize,228,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,1,1
lesize,222,,,16
vsweep,16,106,105


lesize,122,,,2 !左侧芯柱
lesize,131,,,2
mshkey,1
mshape,0,3D
vatt,1,1,1
lesize,126,,,16
vsweep,11,64,63

lesize,110,,,4,-7/3 !左侧砌块
lesize,119,,,4,-13/6
lesize,115,,,4,-7/3
lesize,120,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,2,1,1
lesize,114,,,16
vsweep,12,67,68

lesize,148,,,4,-7/3 !线网格控制
lesize,134,,,2
lesize,151,,,2
lesize,150,,,2
lesize,145,,,2
lesize,147,,,4,-13/6
lesize,146,,,4,-7/3
lesize,149,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,2,1,1
lesize,156,,,16
vsweep,10,70,58

lesize,170,,,4,-7/3 !空心砌块
lesize,158,,,2
lesize,175,,,2
lesize,174,,,2
lesize,169,,,2
lesize,171,,,4,-13/6
lesize,172,,,4,-7/3
lesize,173,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,2,1,1
lesize,180,,,16
vsweep,13,80,78

lesize,206,,,2 !右侧芯柱
lesize,215,,,2
mshkey,1
mshape,0,3D
vatt,1,1,1
lesize,210,,,16
vsweep,15,100,99

lesize,194,,,4,-7/3 !右侧砌块
lesize,195,,,4,-13/6
lesize,196,,,4,-7/3
lesize,197,,,4,-13/6
lesize,182,,,2
lesize,199,,,2
lesize,198,,,2
lesize,193,,,2
mshkey,1
mshape,0,3D
vatt,2,1,1
lesize,204,,,16
vsweep,14,90,88



lesize,392,,,4,-7/3 !顶梁
lesize,401,,,4,-13/6
lesize,397,,,4,-7/3
lesize,402,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,2,1
lesize,396,,,2
vsweep,17,112,111

lesize,404,,,4,-7/3 !顶梁
lesize,413,,,4,-13/6
lesize,409,,,4,-7/3
lesize,414,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,2,1
lesize,408,,,2
vsweep,18,118,117

lesize,416,,,4,-7/3 !顶梁
lesize,425,,,4,-13/6
lesize,421,,,4,-7/3
lesize,426,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,2,1
lesize,420,,,2
vsweep,19,124,123

lesize,428,,,4,-7/3 !顶梁
lesize,437,,,4,-13/6
lesize,433,,,4,-7/3
lesize,438,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,2,1
lesize,432,,,2
vsweep,20,130,129

lesize,440,,,4,-7/3 !顶梁
lesize,449,,,4,-13/6
lesize,445,,,4,-7/3
lesize,450,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,2,1
lesize,444,,,2
vsweep,21,136,135

lesize,452,,,4,-7/3 !顶梁
lesize,461,,,4,-13/6
lesize,457,,,4,-7/3
lesize,462,,,4,-13/6
mshkey,1
mshape,0,3D
vatt,1,2,1
lesize,456,,,2
vsweep,22,142,141

lesize,464,,,1!右侧垫板
lesize,473,,,4,-13/6
lesize,474,,,4,-13/6
lesize,469,,,1
mshkey,1
mshape,0,3D
vatt,5,1,3
lesize,468,,,2
vsweep,23,148,147

lesize,11,,,4,-13/6 !底部垫板
lesize,12,,,4,-13/6
lesize,2,,,5
lesize,7,,,5
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,6,,,3
vsweep,1,4,3

lesize,23,,,4,-13/6 !底部垫板
lesize,24,,,4,-13/6
lesize,14,,,4,-7/3
lesize,19,,,4,-7/3
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,18,,,3
vsweep,2,10,9

lesize,35,,,4,-13/6 !底部垫板
lesize,36,,,4,-13/6
lesize,26,,,4,-7/3
lesize,31,,,4,-7/3
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,30,,,3
vsweep,3,16,15

lesize,47,,,4,-13/6 !底部垫板
lesize,48,,,4,-13/6
lesize,38,,,4,-7/3
lesize,43,,,4,-7/3
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,42,,,3
vsweep,4,22,21

lesize,59,,,4,-13/6 !底部垫板
lesize,60,,,4,-13/6
lesize,50,,,4,-7/3
lesize,55,,,4,-7/3
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,54,,,3
vsweep,5,28,27

lesize,71,,,4,-13/6 !底部垫板
lesize,72,,,4,-13/6
lesize,62,,,4,-7/3
lesize,67,,,4,-7/3
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,66,,,3
vsweep,6,34,33

lesize,83,,,4,-13/6 !底部垫板
lesize,84,,,4,-13/6
lesize,74,,,4,-7/3
lesize,79,,,4,-7/3
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,78,,,3
vsweep,7,40,39

lesize,95,,,4,-13/6 !底部垫板
lesize,96,,,4,-13/6
lesize,86,,,5
lesize,91,,,5
vatt,5,1,3
mshkey,1
mshape,0,3D
lesize,90,,,3
vsweep,8,46,45

!网格划分完成
numcmp,all
nummrg,all

/solu
cnvtol,f,,0.05,,2 !力的收敛
nlgeom,on !考虑大变形
pred,on
OUTRES,ERASE
OUTRES,BASI,1
autots,on
NCNV,2, !如果求解不收敛，终止分析,但不终止程序
NEQIT,40
NSUBST,150,1000,0
kbc,0 !荷载增加方式
LNSRCH,1


NSEL,S,LOC,Y,0 !选中底面节点
D,ALL,ALL !固定底面节点所有自由度
ALLSEL,ALL,NODE !重新选取全部节点，避免意想不到错误

time,1
NSEL,S,LOC,Y,2100
NSEL,u,LOC,x,1550,1600 选择垫板
SF,ALL,PRES,1.2 !于顶面施加1.2MPa的面载
ALLSEL,ALL,NODE !重新选取全部节点，避免意想不到错误
ALLSEL,ALL !重新选中数据库中所有模型

LSWRITE,1 !创建第一个荷载步文件

time,2
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,2

time,3
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,3

time,4
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,4

time,5
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,5

time,6
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,6

time,7
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,7

time,8
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,8

time,9
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,9

time,10
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,10

time,11
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,11

time,12
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,12

time,13
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,20000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,13

time,14
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,14

time,15
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,15

time,16
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,16

time,17
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,17

time,18
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,18

time,19
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,19

time,20
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,20

time,21
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,21

time,22
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,22

time,23
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,23

time,24
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,24

time,25
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,30000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,25

time,26
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-40000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,26

time,27
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,40000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,27

time,28
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-40000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,28

time,29
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,40000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,29

time,30
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-50000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,30

time,31
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,50000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,31

time,32
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,-50000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,32

time,33
NSEL,S,LOC,x,1600
NSEL,r,LOC,y,2000
F,ALL,Fx,50000
ALLSEL,ALL,NODE
ALLSEL,ALL
LSWRITE,33

LSSOLVE,1,33,1

自己更改了子步数 150 能算到29步 在往下算就提示不收敛 请问 子步数还要在增加？谢谢阿