求救求救！！！一品框架低轴往复荷载非线性分析

finish
/cLear,nostart
/fiLname,LiLiang kuangjia
/config,nres,20000 !增加荷载步
/PREP7
!1.定义单元与材料性质
et,1,soLid65,,,,,,,1
mp,ex,1,13585
mp,prxy,1,0.2
dens,1,2.7e9 !定义混凝土单元与材料
mp,ex,2,2.0e5
mp,prxy,2,0.3
!定义钢筋材料
fc=14.3
ft=1.43 !定义混凝土抗压强度及抗拉强度
tb,concr,1
tbdata,,0.5,0.95,fc,ft !定义混凝土材料及相关参数
tb,miso,1,,11 !定义混凝土应力应变曲线
tbpt,,0.0002,fc*0.19$tbpt,,0.0004,fc*0.36$tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64$tbpt,,0.001,fc*0.75$tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0,91$tbpt,,0.0016,fc*0.96$tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc$tbpt,,0.0033,fc*0.85
tb,biso,2$tbdata,,210,0 !定义纵筋应力应变关系
!定义箍筋应力应变关系
pi=acos(-1) !定义参数:
!制作边梁
!计算各区混凝土中钢筋体积配筋率
L1=0.25*Pi*4*4*3/(25*75) !拉区钢筋
L2=0.25*Pi*5*5*4/(75*25) !压区钢筋
L3=0.25*Pi*3*3*5/(75*50) !梁加密区水平
L4=0.25*Pi*3*3*5/(75*125) !梁加密区竖直
L5=0.25*Pi*3*3*16/(50*450) !非加密区水平
L6=0.25*Pi*3*3*16/(75*450) !非加密区竖直
L7=0.25*Pi*3*3*2/(75*75)!梁中部箍筋
!定义实常数
r,1,2,L1,0,90,2,L3$rmore,0,0,2,L4,90,0
r,2,2,L7,0,90
r,3,2,L2,0,90,2,L3$rmore,0,0,2,L4,90,0
r,4,2,L4,O,90,2,15$rmore,0,0,2,L6,90,0
r,5,2,L7,O,90
r,6,2,L3,O,90,2,15$rmore,0,0,2,L6,90,0
r,7,2,L1,O,90,2,15$rmore,0,0,2,L6,90,0
r,8,2,L7,O,90
r,9,2,L2,O,90,2,15$rmore,0,0,2,L6,90,0
!.创建几何模型
bLc4,,,75,125,1350
WProta,,-90$WPoff,,,25$vsbw,aLL
wpoff,,,75$vsbw,aLL$
Wpcsys,-1$Wpoff,,,125$vsbw,aLL
WPoff,,,1100$vsbw,aLL$WPcsys,-1
!定义属性
vseL,s,Loc,z,1225,1350$vseL,r,Loc,y,0,25$vatt,1,1,1
vseL,s,Loc,z,1225,1350$vseL,r,Loc,y,25,100$vatt,1,2,1
vseL,s,Loc,z,1225,1350$vseL,r,Loc,y,100,125$vatt,1,3,1
vseL,s,Loc,z,125,1225$vseL,r,Loc,y,0,25$vatt,1,4,1
vseL,s,Loc,z,125,1225$vseL,r,Loc,y,25,100$vatt,1,5,1
vseL,s,Loc,z,125,1225$vseL,r,Loc,y,100,125$vatt,1,6,1
vseL,s,Loc,z,0,125$vseL,r,Loc,y,0,25$vatt,1,7,1
vseL,s,Loc,z,0,125$vseL,r,Loc,y,25,100$vatt,1,8,1
vseL,s,Loc,z,0,125$vseL,r,Loc,y,100,125$vatt,1,9,1
!网格划分尺寸
aLLseL
LseL,s,Loc,y,0$LseL,r,Loc,z,0$Lesize,aLL,12.5
LseL,s,Loc,x,0$LseL,r,Loc,z,0$Lesize,aLL,25
LseL,s,Loc,x,0$LseL,r,Loc,y,0$Lesize,aLL,100
!建立中柱
wpcsys,-1$wpoffs,-37.5$wpoffs,,-775$wprota,,-90$bLc4,,,150,150,900!建立中柱实体
wpoffs,,,125$vsbw,aLL!切出混凝土箍筋加密区
wpoffs,,,650$vsbw,aLL!切出梁柱节点区域
wpoffs,,,-775$wprota,,-90$wpoff,,,25$vsbw,aLL!切出右侧钢筋区域
wpoff,,,100$vsbw,aLL!切出左侧钢筋区域
!计算柱子体积配筋率
z1=0.25*Pi*5*5*4/25*150!节点处左侧柱子主筋与右侧相等
z2=0.25*Pi*5*5*4/(150*100)!节点处中间柱子中间
z3=0.25*Pi*3*3*5/(150*50)!加密区柱子水平箍筋
z4=0.25*Pi*5*5*4/(100*150)!加密区柱子竖直箍筋
z5=0.25*Pi*5*5*4/(650*25)!非加密区柱子主筋左侧
z6=0.25*Pi*3*3*20/(650*50)!非加密区柱子水平箍筋
z7=0.25*Pi*3*3*20/(100*650)!非加密区柱子竖直箍筋
!定义实常数
r,11,2,z1,90,0,2,z3$rmore,0,0,2,z4,0,90
r,12,2,z2,90,0,2,z4$rmore,0,90
r,13,2,z1,90,0,2,z3$rmore,0,0,2,z4,0,90
r,14,2,z1,90,0,2,z6$rmore,0,0,2,z7,0,90
r,15,2,z5,90,0,2,z7$rmore,0,90
r,16,2,z1,90,0,2,z6$rmore,0,0,2,z7,0,90
r,17,2,z1,90,0,2,z3$rmore,0,0,2,z4,0,90
r,18,2,z2,90,0,2,z4$rmore,0,90
r,19,2,z1,90,0,2,z3$rmore,0,0,2,z3,0,90
wpcsys,-1
!定义属性
vseL,s,Loc,y,-775,-650$vseL,r,Loc,z,0,-25$vatt,1,11,1
vseL,s,Loc,y,-775,-650$vseL,r,Loc,z,-25,-125$vatt,1,12,1
vseL,s,Loc,y,-775,-650$vseL,r,Loc,z,-125,-150$vatt,1,13,1
vseL,s,Loc,y,-650,0$vseL,r,Loc,z,0,-25$vatt,1,14,1
vseL,s,Loc,y,-650,0$vseL,r,Loc,z,-25,-125$vatt,1,15,1
vseL,s,Loc,y,-650,0$vseL,r,Loc,z,-125,-150$vatt,1,16,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,0,-25$vatt,1,17,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,-25,-125$vatt,1,18,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,-125,-150$vatt,1,19,1
!网格划分属性
LseL,s,Loc,y,-775$LseL,r,Loc,z,-150$Lesize,aLL,12.5
LseL,s,Loc,x,-37.5$LseL,r,Loc,y,-775$Lesize,aLL,50
LseL,s,Loc,x,-37.5$LseL,r,Loc,z,-150$LseL,r,Loc,y,-775,0$Lesize,aLL,100
LseL,s,Loc,x,-37.5$LseL,r,Loc,z,-150$LseL,r,Loc,y,0,125$Lesize,aLL,25
!制作边柱
wpcsys,-1$wpoffs,-37.5$wpoffs,,,1500$wpoffs,,-775
wprota,,-90$bLc4,,,150,150,900 !建立边柱实体
wpoffs,,,775$vsbw,aLL!切出梁柱节点区域
wpoffs,,,-775$wprota,,-90$wpoff,,,25$vsbw,aLL!切出右侧钢筋区域
wpoff,,,100$vsbw,aLL!切出左侧钢筋区域
!计算柱体积配筋率
z11=0.25*Pi*5*5*4/25*150!节点处左侧柱子主筋与右侧主筋
z12=0.25*Pi*5*5*4/(150*100)!节点处左侧柱子中间主筋
z13=0.25*Pi*3*3*5/(150*50)!节点柱子水平箍筋
z14=0.25*Pi*3*3*5/(100*150)!节点柱子竖直箍筋
z15=0.25*Pi*5*5*4/(775*25)!柱子主筋左侧
z16=0.25*Pi*5*5*4/(775*100)!柱子中部主筋
z17=0.25*Pi*3*3*23/(775*50)!柱子水平箍筋
z18=0.25*Pi*3*3*23/(100*775)!柱子竖直箍筋
!定义实常数
r,21,2,z11,90,0,2,z13$rmore,0,0,2,z14,0,90
r,22,2,z12,90,0,2,z14$rmore,0,90
r,23,2,z11,90,0,2,z13$rmore,0,0,2,z14,0,90
r,24,2,z15,90,0,2,z17$rmore,0,0,2,z18,0,90
r,25,2,z16,90,0,2,z18$rmore,0,90
r,26,2,z15,90,0,2,z17$rmore,0,0,2,z18,0,90
wpcsys,-1
!定义属性
vseL,s,Loc,y,0,125$vseL,r,Loc,z,1475,1500$vatt,1,21,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,1375,1475$vatt,1,22,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,1350,1375$vatt,1,23,1
vseL,s,Loc,y,-775,0$vseL,r,Loc,z,1475,1500$vatt,1,24,1
vseL,s,Loc,y,-775,0$vseL,r,Loc,z,1375,1475$vatt,1,25,1
vseL,s,Loc,y,-775,0$vseL,r,Loc,z,1350,1375$vatt,1,26,1
!网格划分属性
LseL,s,Loc,y,-775$LseL,r,Loc,z,1350$Lesize,aLL,12.5
LseL,s,Loc,x,-37.5$LseL,r,Loc,y,-775$LseL,r,Loc,z,1350,1500$Lesize,aLL,50
LseL,s,Loc,x,-37.5$LseL,r,Loc,z,1350$LseL,r,Loc,y,-775,0$Lesize,aLL,100
LseL,s,Loc,x,-37.5$LseL,r,Loc,z,1350$LseL,r,Loc,y,0,125$Lesize,aLL,25
!制作边梁
wpcsys,-1$wprota,,,90$wprota,,,90$wpoffs,,,150$wpoffs,-75
bLc4,,,75,125,1350
L11=0.25*Pi*4*4*3/(25*75) !拉区钢筋
L12=0.25*Pi*5*5*4/(75*25) !压区钢筋
L13=0.25*Pi*3*3*5/(75*50) !梁加密区水平
L14=0.25*Pi*3*3*5/(75*125) !梁加密区竖直
L15=0.25*Pi*3*3*16/(50*450) !非加密区水平
L16=0.25*Pi*3*3*16/(75*450) !非加密区竖直
L17=0.25*Pi*3*3*2/(75*75)!梁中部箍筋
!定义实常数
r,31,2,L11,0,90,2,L3$rmore,0,0,2,L14,90,0
r,32,2,L17,0,90
r,33,2,L12,0,90,2,L13$rmore,0,0,2,L14,90,0
r,34,2,L14,O,90,2,L15$rmore,0,0,2,L16,90,0
r,35,2,L17,O,90
r,36,2,L13,O,90,2,L15$rmore,0,0,2,L16,90,0
r,37,2,L11,O,90,2,L15$rmore,0,0,2,L16,90,0
r,38,2,L17,O,90
r,39,2,L12,O,90,2,L15$rmore,0,0,2,L16,90,0
wprota,,-90
WPoff,,,25$vsbw,aLL
wpoff,,,75$vsbw,aLL
wprota,,90
Wpoff,,,125$vsbw,aLL
WPoff,,,1100$vsbw,aLL$WPcsys,-1
!定义属性
vseL,s,Loc,z,-275,-150$vseL,r,Loc,y,0,25$vatt,1,31,1
vseL,s,Loc,z,-275,-150$vseL,r,Loc,y,25,100$vatt,1,32,1
vseL,s,Loc,z,-275,-150$vseL,r,Loc,y,100,125$vatt,1,33,1
vseL,s,Loc,z,-1375,-275$vseL,r,Loc,y,0,25$vatt,1,34,1
vseL,s,Loc,z,-1375,-275$vseL,r,Loc,y,25,100$vatt,1,35,1
vseL,s,Loc,z,-1375,-275$vseL,r,Loc,y,100,125$vatt,1,36,1
vseL,s,Loc,z,-1500,-1375$vseL,r,Loc,y,0,25$vatt,1,37,1
vseL,s,Loc,z,-1500,-1375$vseL,r,Loc,y,25,100$vatt,1,38,1
vseL,s,Loc,z,-1500,-1375$vseL,r,Loc,y,100,125$vatt,1,39,1
!网格划分尺寸
aLLseL
LseL,s,Loc,y,0$LseL,r,Loc,z,-150$Lesize,aLL,12.5
LseL,s,Loc,x,0$LseL,r,Loc,z,-150$Lesize,aLL,25
LseL,s,Loc,x,0$LseL,r,Loc,y,0$LseL,r,Loc,z,-1500,-150
Lesize,aLL,100
!建立右边柱
wpcsys,-1$wpoffs,-37.5$wpoffs,,,-1500$wpoffs,,-775
wprota,,-90$bLc4,,,150,150,900 !建立边柱实体
wpoffs,,,775$vsbw,aLL!切出梁柱节点区域
wpoffs,,,-775$wprota,,-90$wpoff,,,25$vsbw,aLL!切出右侧钢筋区域
wpoff,,,100$vsbw,aLL!切出左侧钢筋区域
!计算柱体积配筋率
z21=0.25*Pi*5*5*4/25*150!节点处左侧柱子主筋与右侧主筋
z22=0.25*Pi*5*5*4/(150*100)!节点处左侧柱子中间主筋
z23=0.25*Pi*3*3*5/(150*50)!节点柱子水平箍筋
z24=0.25*Pi*3*3*5/(100*150)!节点柱子竖直箍筋
z25=0.25*Pi*5*5*4/(775*25)!柱子主筋左侧
z26=0.25*Pi*5*5*4/(775*100)!柱子中部主筋
z27=0.25*Pi*3*3*23/(775*50)!柱子水平箍筋
z28=0.25*Pi*3*3*23/(100*775)!柱子竖直箍筋
!定义实常数
r,41,2,z21,90,0,2,z13$rmore,0,0,2,z24,0,90
r,42,2,z22,90,0,2,z14$rmore,0,90
r,43,2,z21,90,0,2,z13$rmore,0,0,2,z24,0,90
r,44,2,z25,90,0,2,z17$rmore,0,0,2,z28,0,90
r,45,2,z26,90,0,2,z18$rmore,0,90
r,46,2,z25,90,0,2,z17$rmore,0,0,2,z28,0,90
wpcsys,-1
!定义属性
vseL,s,Loc,y,0,125$vseL,r,Loc,z,-1500,-1525$vatt,1,41,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,-1525,-1625$vatt,1,42,1
vseL,s,Loc,y,0,125$vseL,r,Loc,z,-1625,-1650$vatt,1,43,1
vseL,s,Loc,y,-775,0$vseL,r,Loc,z,-1500,-1525$vatt,1,44,1
vseL,s,Loc,y,-775,0$vseL,r,Loc,z,-1525,-1625$vatt,1,45,1
vseL,s,Loc,y,-775,0$vseL,r,Loc,z,-1625,-1650$vatt,1,46,1
!网格划分属性
LseL,s,Loc,y,-775$LseL,r,Loc,z,-1500$Lesize,aLL,12.5
LseL,s,Loc,x,-37.5$LseL,r,Loc,y,-775$LseL,r,Loc,z,-1500,-1650$Lesize,aLL,50
LseL,s,Loc,x,-37.5$LseL,r,Loc,z,-1500$LseL,r,Loc,y,-775,0$Lesize,aLL,100
LseL,s,Loc,x,-37.5$LseL,r,Loc,z,-1500$LseL,r,Loc,y,0,125$Lesize,aLL,25
!复制框架成型
aLLseL$vgen,3,aLL,,,,900,,,,0 !将框架复制多层
vsel,all
mshkey,1
vmesh,all
Nummrg,elem
Nummrg,node
Nummrg,kp
numcmp,aLL !合并图素,框架成型
!4.施加荷载和约束
/SOLU
Allsel$acel,,-9810!施加重力场
nsEl,s,loc,y,-775$d,all,all
wpoff,,,1500
Allsel$/solu
antype,trans
timint,off !打开瞬态效应
kbc,1 !荷载参数为阶跃荷载
trnopt,full !瞬态分析方法采用完全法
nropt,full !方程求解NR法
alphad,alp
betad,beta
ncnv,2,le8
cnvtol,u,,0.015
nsubst,25
outres,all,all
autots,on
lnsrch,on
neqit,50
*DIM,time_f,array,800,1
*creat,ansuit
*VREAD,time_f,saul.txt,,,1,nt
(1f10.9)
*END
/input,ansuit
*set,nt,800
*set,dt,0.02
*dim,time_f,,nt,1
*do,i,1,nt
asel,s,loc,z,0$asel,s,loc,y,1800,1925
D,nodel,all,d(t)
time,i0.02*
solve
*eddo
看看问题出在哪儿
求救！