一个混凝土剪力墙的命令流讨论

zyq309

各位高手，我最近在做一个剪力墙的分析，出了几个问题，下面是小弟的命令流。请各位仁兄帮我看看，怎么改进较好。
我发现当我nsel,s,loc,y,2.31后，显示我选择了节点有两层，而我nsel,s,loc,y,0.6时显示我选择的节点是一层，这才是正确的嘛！为什么会显示两层节点呢？因为后面我要nummrg,all，两层节点是要不得的。然后就是我在最左边端柱上端加了一个分布力6e5N/m*0.18m=108KN,而这时我剪力墙上的应力只有0.029Pa。这也太小了吧！我检查了我的单位，长度是m，力是N，质量是Kg，弹性模量是N/mm2没错啊！在此谢谢各位！请大家给我指点。


/FILNAME,SHEAR-WALL
/TITLE,THE ANALYSIS OF A SHEAR-WALL
/PREP7
ET,1,SOLID65
ET,3,SHELL181
!关闭压碎,keyopt(1)=0,keyopt(7)=1
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
!定义混凝土实常数
r,1,2,0.0132,0,0,2,0.0115  
rmore,90,0,2,0.0132,0,90   !1～2层约束边缘端柱的x，y，z方向的配筋
r,2,2,0.0205,0,0,2,0.0209  
rmore,90,0,2,0.0205,0,90   !1～10层暗柱的x，y，z方向的配筋
r,3,2,0.0052,0,0,2,0.0052  
rmore,90,0,2,0.0026,0,90   !1～10层剪力墙身的x，y，z方向的配筋
r,4,2,0.0057,0,0,2,0.0035  
rmore,90,0,2,0.0057,0,90   !3～10层构造边缘端柱的x，y，z方向的配筋
!连梁单元
R,5,0.0034,0.0034,0.0034,0.0034
R,6,0.0022,0.0022,0.0022,0.0022
R,7,0.0028,0.0028,0.0028,0.0028
R,8,0.0018,0.0018,0.0018,0.0018
R,9,0.0026,0.0026,0.0026,0.0026
R,10,0.0016,0.0016,0.0016,0.0016
MP,EX,3,2.1E11
MP,PRXY,3,0.3
MP,DENS,3,7850
!混凝土基本材料属性
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.3585e10
MPDATA,PRXY,1,,0.2
fc=14.3e6
ft=1.43e6
!屈服准则
tb,miso,1,1,11
TBTEMP,0
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,concr,1
tbdata,,0.5,1,ft,-1
TBPLOT,miso,1  !检查1曲线
!钢材基本属性
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,210e9
MPDATA,PRXY,2,,0.3
!屈服准则
TB,BKIN,2,1,2,1
TBTEMP,0
TBDATA,,215e6,10.3e6,,,,
TBPLOT,bkin,2 !检查2曲线
!建立模型
!连梁
wpcsys,0
wpoffs,0.874,0.48,0.09
blc4,0,0,0.832,0.12
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.12
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
asel,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.12
wprota,0,90,0
wpoffs,-0.874,-0.48,-0.09
asel,all
aplot
wpoffs,0.874,1.08,0.09
blc4,0,0,0.832,0.12
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.12
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,16,20,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.12
wprota,0,90,0
wpoffs,-0.874,-1.08,-0.09
asel,all
aplot
wpoffs,0.874,1.71,0.09
blc4,0,0,0.832,0.09
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.09
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,31,35,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.09
wprota,0,90,0
wpoffs,-0.874,-1.71,-0.09
asel,all
aplot
wpoffs,0.874,2.31,0.09
blc4,0,0,0.832,0.09
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.09
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,46,50,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.09
wprota,0,90,0
wpoffs,-0.874,-2.31,-0.09
asel,all
aplot
wpoffs,0.874,2.91,0.09
blc4,0,0,0.832,0.09
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.09
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,all
asel,s,area,,61
asel,a,area,,62,65,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.09
wprota,0,90,0
wpoffs,-0.874,-2.91,-0.09
asel,all
aplot
wpoffs,0.874,3.52,0.09
blc4,0,0,0.832,0.08
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.08
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,76,80,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.08
wprota,0,90,0
wpoffs,-0.874,-3.52,-0.09
asel,all
aplot
wpoffs,0.874,4.12,0.09
blc4,0,0,0.832,0.08
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.08
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,91,95,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.08
wprota,0,90,0
wpoffs,-0.874,-4.12,-0.09
asel,all
aplot
wpoffs,0.874,4.72,0.09
blc4,0,0,0.832,0.08
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.08
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,106,110,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.08
wprota,0,90,0
wpoffs,-0.874,-4.72,-0.09
asel,all
aplot
wpoffs,0.874,5.32,0.09
blc4,0,0,0.832,0.08
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.08
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,121,125,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.08
wprota,0,90,0
wpoffs,-0.874,-5.32,-0.09
asel,all
aplot
wpoffs,0.874,5.92,0.09
blc4,0,0,0.832,0.08
wprota,0,-90,0
blc4,0,0,0.832,0.02
blc4,0,0,0.832,-0.02
wpoffs,0,0,0.08
blc4,0,0,0.832, 0.02
blc4,0,0,0.832,-0.02
asel,s,area,,136,140,1
aglue,all
nummrg,all
numcmp,all
wpoffs,0.116,,,
wprota,0,0,90
asbw,all
wpoffs,,,0.6
asbw,all
aglue,all
nummrg,all
numcmp,all
wpoffs,,,-0.6
wprota,0,0,-90
wpoffs,-0.116,,,
wpoffs,0,0,-0.08
wprota,0,90,0
wpoffs,-0.874,-5.92,-0.09
asel,all
aplot


!剪力墙
block,0,0.18,0,6,0,0.18
block,0.18,0.99,0,6,0.06,0.12
vsel,all
vglue,1,2
nummrg,all
numcmp,all
wpcsys,0
wpoffs,0.874,0,0
wprota,0,-90,0
wprota,0,0,90
vsbw,1         
wprota,0,0,-90
wprota,0,90,0
wpoffs,-0.874,0,0
Block,1.59,2.40,0,6,0.06,0.12
Block,2.40,2.58,0,6,0,0.18
Vglue,1,5
nummrg,all
numcmp,all
wpoffs,1.706,0,0
wprota,0,-90,0
wprota,0,0,90
vsbw,1
nummrg,all
numcmp,all
!得到6个体积，从左往右是V1,V2,V3,V5,V6,V4。
wprota,0,0,-90
wprota,0,90,0
wpoffs,-1.706,0,0
wpoffs,0,0,0.06
wprota,-90,0,0
vsbw,1
wpoffs,0,0,0.06
vsbw,8
nummrg,all
numcmp,all
wpoffs,2.4,0,0
vsbw,4
wpoffs,0,0,-0.06
vsbw,10
nummrg,all
numcmp,all
!                         V7                 V10
!得到10个体积，从左往右是 V8,V2,V3,    V5,V6,V4。
!                         V1                 V9

asel,s,area,,2,5,1
asel,a,area,,7,10,1
asel,a,area,,12,15,1
asel,a,area,,17,20,1
asel,a,area,,22,25,1
asel,a,area,,27,30,1
aatt,3,5,3
asel,all
asel,s,area,,1,26,5
aatt,3,6,3
asel,all
asel,s,area,,32,35,1
asel,a,area,,37,40,1
asel,a,area,,42,45,1
asel,a,area,,47,50,1
asel,a,area,,52,55,1
asel,a,area,,57,60,1
asel,a,area,,62,65,1
asel,a,area,,67,70,1
asel,a,area,,72,75,1
aatt,3,7,3
asel,all
asel,s,area,,31,71,5
aatt,3,8,3
asel,all
asel,s,area,,77,80,1
asel,a,area,,82,85,1
asel,a,area,,87,90,1
asel,a,area,,92,95,1
asel,a,area,,97,100,1
asel,a,area,,102,105,1
asel,a,area,,107,110,1
asel,a,area,,112,115,1
asel,a,area,,117,120,1
asel,a,area,,122,125,1
asel,a,area,,127,130,1
asel,a,area,,132,135,1
asel,a,area,,137,140,1
asel,a,area,,142,145,1
asel,a,area,,147,150,1
aatt,3,9,3
asel,all
asel,s,area,,76,146,5
aatt,3,10,3
asel,all
lsel,all
lsel,s,line,,1,380,1
lesize,all,0.01
mshkey,1
mshape,0
asel,s,area,,1,150,1
et,3,shell181
amesh,all

vsel,s,volu,,2
vsel,a,volu,,6
vatt,1,3,1
vsel,all
vsel,s,volu,,3
vsel,a,volu,,5
vatt,1,2,1
vsel,all
vsel,s,volu,,1
vsel,a,volu,,4
vsel,a,volu,,7,10,1
vatt,1,1,1
vsel,all
vsel,s,volu,,3
vsel,a,volu,,5
esize,0.01
mshkey,1
mshape,0,3D
vmesh,all
vsel,all
eplot
vsel,s,volu,,2
vsel,a,volu,,6
vsel,a,volu,,1
vsel,a,volu,,4
vsel,a,volu,,7,10,1
esize,0.1
mshkey,1
mshape,0,3D
vmesh,all
vsel,all
eplot


nsel,s,loc,y,0.48
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,0.6
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,1.08
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,1.2
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,1.71
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,1.8
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,2.31
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,2.4
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,2.91
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,3.0
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,3.52
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,3.6
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,4.12
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,4.2
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,4.72
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,4.8
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,5.32
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,5.4
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,5.92
nplot
nummrg,all
numcmp,all
nsel,all
nsel,s,loc,y,6.0
nplot
nummrg,all
numcmp,all
nsel,all
allsel,all
eplot
finish

/solu
nsel,s,loc,y,0
d,all,all,0
nsel,all
nsel,s,loc,y,6
nsel,a,loc,x,0,0.99
nsel,a,loc,x,1.59,2.58
d,all,uz,0
allsel,all
f,
allsel,all

antype,0
solcontrol,on
nlgeom,on        !打开几何非线性的影响
kbc,0
nropt,full
cnvtol,f,,0.05,2,1   !定义收敛条件,使用缺省的VALUE ,力的二范数收敛准则!，5％
nsubst,20         !定义子步数   
outres,all,all       !输出每一子步的结果  
autots,1  
lnsrch,1           !打开线性搜索  
ncnv,2            !如果不收敛时结束而不退出  
neqit,10           !每一子步中方程的迭代次数限值
autots,on          !打开自动时间步控制，若50内不收敛，采用二分法
pred,on           !打开预测器  
solve
finish