关于框架结构的整体分析命令流！求指正
我一个钢筋混凝土框架结构的整体抗震分析，为一七层的框架结构，我的反应谱分析结果最大位移只有5mm, 是不是太小了？我也不知道哪出了问题 ，麻烦各位高手指点[ 本帖最后由 xuruikl 于 2007-11-19 18:44 编辑 ] 这是命令流
/prep7
/title,seismic analysis of frame
et,1,beam188
KEYOPT,1,7,1
KEYOPT,1,3,2
KEYOPT,1,8,1
et,2,beam188
KEYOPT,2,7,1
KEYOPT,1,3,2
KEYOPT,2,8,1
et,3,shell63
!下部梁
mp,ex,1,3.34e10
mp,prxy,1,0.2
mp,dens,1,2540
!下部柱子
mp,ex,2,5.79e10
mp,prxy,2,0.2
mp,dens,2,2550
!板
mp,ex,3,2.00e10
mp,prxy,3,0.2
mp,dens,3,3000
!梁板上附加质量
ET,6,MASS21
KEYOPT,6,1,0
KEYOPT,6,2,0
KEYOPT,6,3,2
!上部梁
mp,ex,4,3.18e10
mp,prxy,4,0.1667
mp,dens,4,2540
!上部柱子
mp,ex,5,5.73e10
mp,prxy,5,0.1667
mp,dens,5,2550
r,1
r,2
r,3,0.1
r,4
r,5
r,6,340
!梁截面
SECTYPE, 1, BEAM, RECT, beam2, 0
SECOFFSET, USER, 0, 0.325
SECDATA,0.4,0.65,0,0,0,0,0,0,0,0
!柱截面
SECTYPE, 2, BEAM, RECT, column, 0
SECOFFSET, CENT
SECDATA,0.8,0.8,0,0,0,0,0,0,0,0
!柱子关键点
k,1,-13.2,0,0
k,2,-6.6,0,0
k,3,0,0,0
k,4,6.6,0,0
k,5,13.2,0,0
k,6,-13.2,6.6,0
k,7,-6.6,6.6,0
k,8,0,6.6,0
k,9,6.6,6.6,0
k,10,13.2,6.6,0
k,11,-13.2,13.2,0
k,12,-6.6,13.2,0
k,13,0,13.2,0
k,14,6.6,13.2,0
k,15,13.2,13.2,0
k,16,-6.6,19.8,0
k,17,0,19.8,0
k,18,6.6,19.8,0
!生成各层柱关键点
FLST,3,18,3,ORDE,2
FITEM,3,1
FITEM,3,-18
KGEN,2,P51X, , , , ,3.9,100,0
FLST,3,18,3,ORDE,2
FITEM,3,101
FITEM,3,-118
KGEN,7,P51X, , , , ,3,100,0
!定义二层主梁
*do,ii,101,104
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,116,117
LSTR,ii,ii+1
*enddo
*do,ii,101,105
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,112,114
LSTR,ii,ii+4
*enddo
!定义三层主梁
*do,ii,201,204
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,216,217
LSTR,ii,ii+1
*enddo
*do,ii,201,205
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,212,214
LSTR,ii,ii+4
*enddo
!定义四层主梁
*do,ii,301,304
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,316,317
LSTR,ii,ii+1
*enddo
*do,ii,301,305
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,312,314
LSTR,ii,ii+4
*enddo
!定义五层主梁
*do,ii,401,404
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,416,417
LSTR,ii,ii+1
*enddo
*do,ii,401,405
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,412,414
LSTR,ii,ii+4
*enddo
!定义六层主梁
*do,ii,501,504
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,516,517
LSTR,ii,ii+1
*enddo
*do,ii,501,505
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,512,514
LSTR,ii,ii+4
*enddo
!定义七层主梁
*do,ii,601,604
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,616,617
LSTR,ii,ii+1
*enddo
*do,ii,601,605
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,612,614
LSTR,ii,ii+4
*enddo
!定义顶层主梁
*do,ii,701,704
LSTR,ii,ii+1
LSTR,ii+5,ii+6
LSTR,ii+10,ii+11
*enddo
*do,ii,716,717
LSTR,ii,ii+1
*enddo
*do,ii,701,705
LSTR,ii,ii+5
LSTR,ii+5,ii+10
*enddo
*do,ii,712,714
LSTR,ii,ii+4
*enddo
!定义柱
*do,ii,1,18
l,ii,ii+100
l,ii+100,ii+200
l,ii+200,ii+300
l,ii+300,ii+400
l,ii+400,ii+500
l,ii+500,ii+600
l,ii+600,ii+700
*enddo
!定义二层楼板
*do,ii,1,2
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,4,5
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,7,8
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,10,11
al,ii,ii+13,ii+1,ii+11
*enddo
al,6,26,13,25
al,9,27,14,26
!定义三层楼板
*do,ii,28,29
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,31,32
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,34,35
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,37,38
al,ii,ii+13,ii+1,ii+11
*enddo
al,33,53,40,52
al,36,54,41,53
!定义四层楼板
*do,ii,55,56
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,58,59
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,61,62
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,64,65
al,ii,ii+13,ii+1,ii+11
*enddo
al,60,80,67,79
al,63,81,68,80
!定义五层楼板
*do,ii,82,83
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,85,86
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,88,89
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,91,92
al,ii,ii+13,ii+1,ii+11
*enddo
al,87,107,94,106
al,90,108,95,107
!定义六层楼板
*do,ii,109,110
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,112,113
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,115,116
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,118,119
al,ii,ii+13,ii+1,ii+11
*enddo
al,114,134,121,133
al,117,135,122,134
!定义七层楼板
*do,ii,136,137
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,139,140
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,142,143
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,145,146
al,ii,ii+13,ii+1,ii+11
*enddo
al,141,161,148,160
al,144,162,149,161
!定义顶层层楼板
*do,ii,163,164
al,ii,ii+16,ii+1,ii+14
*enddo
*do,ii,166,167
al,ii,ii+15,ii+1,ii+13
*enddo
*do,ii,169,170
al,ii,ii+14,ii+1,ii+12
*enddo
*do,ii,172,173
al,ii,ii+13,ii+1,ii+11
*enddo
al,168,188,175,187
al,171,189,176,188
!划分下部主梁
TYPE, 1
MAT, 1
REAL, 1
ESYS, 0
SECNUM, 1
FLST,5,108,4,ORDE,2
FITEM,5,1
FITEM,5,-108
LSEL,S, , ,P51X
CM,bbeam,LINE
cmsel,,bbeam
lesize,bbeam,,,6,,,,,1
lmesh,bbeam
!划分上部主梁
TYPE, 1
MAT, 4
REAL, 4
ESYS, 0
SECNUM, 1
FLST,5,81,4,ORDE,2
FITEM,5,109
FITEM,5,-189
LSEL,S, , ,P51X
CM,tbeam,LINE
cmsel,,tbeam
lesize,tbeam,,,6,,,,,1
lmesh,tbeam
!划分下部柱
TYPE, 2
MAT, 2
REAL, 2
ESYS, 0
SECNUM, 2
FLST,5,72,4,ORDE,2
FITEM,5,190
FITEM,5,-261
LSEL,S, , ,P51X
CM,bcolumn,LINE
cmsel,,bcolumn
lesize,bcolumn,,,6,,,,,1
lmesh,bcolumn
!划分上部柱
TYPE, 2
MAT, 5
REAL, 5
ESYS, 0
SECNUM, 2
FLST,5,54,4,ORDE,2
FITEM,5,262
FITEM,5,-315
LSEL,S, , ,P51X
CM,tcolumn,LINE
cmsel,,tcolumn
lesize,tcolumn,,,6,,,,,1
lmesh,tcolumn
!划分楼板
TYPE, 3
MAT, 3
REAL, 3
ESYS, 0
FLST,5,70,5,ORDE,2
FITEM,5,1
FITEM,5,-70
LSEL,S, , ,P51X
CM,floor,area
cmsel,,floor
esize,1.1
mshape,0,2d
mshkey,1
amesh,floor
!梁上恒荷
lsel,s,type,,1
nsll,s,1
NPLOT
*get,node,node,,count
*get,minnd,node,,num,min
*get,maxnd,node,,num,max
TYPE, 6
REAL, 6
*do,ii,minnd,maxnd
e,ii
*enddo
NUMMRG,NODE, , , ,LOW
NUMCMP,NODE
nsel,all
nplot
finish
/SOLU
ANTYPE,0
nsel,s,loc,z,0
d,all,all
nsel,all
acel,,,9.8
SOLVE
finish
/SOLU
ANTYPE,MODAL
MODOPT,SUBSP,20
SUBOPT,,8
SOLVE
SAVE
FINISH
/solu
antype,spectr
spopt,sprs,,yes
DMPRAT,0.05
sed,,1
svtype,2
FREQ,0.100,0.125,0.143,0.167,0.200,0.250,0.333,0.400,
FREQ,0.444,0.500,0.571,0.667,0.800,1.000,1.333,2.000,
FREQ,2.857,10.000,1000.000
SV,,0.055,0.086,0.102,0.118,0.133,0.149,0.165,0.172,
SV,,0.176,0.180,0.184,0.2112,0.249,0.305,0.395,0.569,
SV,,0.784,0.784,0.353
solve
finish
/solu
antype,modal
expass,on
mxpand,20,,,yes,,
solve
finish
/solu
antype,spectr
srss,,disp
Outres ,all ,all
solve
finish
/POST1
SET,LIST
/INPUT,'seismic','mcom','D:\boss\',, 0
PLDISP,1
!反应谱结果分析，主要是层间位移和层间位移角
nsel,s,loc,x,0
nsel,r,loc,y,13.2
esln,s,1
etable,uy,u,y
pretab,uy
*get,node,node,,count
*get,minnd,node,,num,min
*get,maxnd,node,,num,max
*dim,AA,,8,1
FLOOR=7
*dim,CANGLE,TABLE,FLOOR,2
*dim,CDDISP,TABLE,8,2
*dim,HEIGH,,FLOOR
AA(1)=0
heigh(1)=3.9
*DO,I,2,FLOOR
AA(I)=Uy(MINND)
cc=nz(ndnext(MINND))
heigh(i)=cc-nz(minnd)
minnd=ndnext(MINND)
*enddo
AA(8)=uy(952)
*do,jj,1,8
cddisp(jj,0)=jj
cddisp(jj,1)=AA(jj)
*enddo
*do,ii,2,8
dd=abs(AA(ii)-AA(ii-1))
cangle(ii-1,1)=dd/heigh(ii-1)
cangle(ii-1,0)=ii-1
*enddo
save
nsel,s,loc,x,0
nsel,r,loc,y,0
esln,s,1
etable,uy,u,y
pretab,uy
*get,node,node,,count
*get,minnd,node,,num,min
*get,maxnd,node,,num,max
*dim,AA1,,8,1
AA1(1)=0
FLOOR=7
*dim,HEIGH1,,FLOOR
HEIGH1(1)=3.9
*DO,I,2,FLOOR
AA1(I)=Uy(MINND)
ff=nz(ndnext(MINND))
heigh1(i)=ff-nz(minnd)
minnd=ndnext(MINND)
*enddo
AA1(8)=uy(940)
*do,jj,1,8
cddisp(jj,2)=AA1(jj)
*enddo
*do,ii,2,8
dd=abs(AA1(ii)-AA1(ii-1))
cangle(ii-1,2)=dd/heigh1(ii-1)
*vplot,cddisp(1,0),cddisp(1,1),2
*vplot,cangle(1,0),cangle(1,1),2
!时程分析
finish
/config,nres,20000
*DIM,aa,,1501
*dim,accel,,1501
*CREATE,ansuitmp
*VREAD,aa,'record','txt',' ',1, , , , , ,
(f8.3)
*END
/INPUT,ansuitmp
*vabs,1,1
*vscfun,mmax,max,aa
a=0.35/mmax
*do,i,1,1501
accel(i)=aa(i)*a
*enddo
/solu
allsel
antype,trans
nsubst,1
auto,off
kbc,0
dmprat,0.05
outres,all,all
*do,i,1,1501
time,i*0.02
acel,,accel(i)
solve
*enddo 无数次的看怎么还没有高手伸援手呢？qq:121819959 谢谢了
页:
[1]