求助 钢框架的端板连接的有限元分析

hnnh

下面是我的命令流，求高手帮我看看 为什么就solve不了的呢 ？

/prep7
!定义参数
hc=300          !柱截面高度
bc=300          !柱截面宽度
tcf=15        !柱翼缘厚度
tcw=10        !柱腹板厚度
lc=1800         !柱构件伸出长度
hb=244            !梁截面高度
bb=210           !梁截面宽度
tbf=11           !梁翼缘厚度
tbw=7          !梁腹板厚度
lb=1500         !梁构件伸出长度

hd=400      !端板高度
bd=250     !端板宽度
td=10       !端板厚度
miu=0.4                   !板间抗滑移系数
lbt=tcf+td            !螺栓杆长度
dbt=25                     !螺栓杆直径
dbth=30                    !螺栓头和螺母直径
lbth=12.5                  !螺栓头厚度
dh0=dbt+2                  !螺栓孔直径
preten=155000              !螺栓施工预拉力
!定义材料和单元属性
!定义材料1，Q345钢材用于梁,柱，端板
et,1,solid45
r,1
mp,ex,1,2e5
mp,prxy,1,0.3
mp,dens,1,7.85e-9
!定义材料用于高强螺栓
et,2,solid92
r,2
mp,ex,2,2e7
mp,prxy,2,0.3
!梁
block,0,lb,-(hb/2-tbf),hb/2-tbf,0,tbw/2
wpoffs,,hb/2-tbf,,
block,0,lb,0,tbf,0,bb/2
vgen,2,2,,,,-hb+tbf,, ,0
! other side duanban
wpoffs,lb,-(hb/2-tbf),,
block,0,tcf,-hb/2,hb/2,0,bb/2

!端板
wpoffs,-lb,,,
block,0,-td,-hd/2,hd/2,0,bd/2

!柱
wpoffs,-td,,,
block,0,-tcf,-lc/2,lc/2,0,bc/2
wpoffs,-tcf,,
block,0,-(hc-2*tcf),-lc/2,lc/2,0,tcw/2
vgen,2,6,,,-bc+tcf,,, ,0

!打孔
wpoffs,tcf+td,hd/6,20
wprota,,,90,
cylind,dh0/2,,0,-lbt,0,360
vgen,2,9,,,,-hd/3,,
vgen,2,10,,,,-85,,
vgen,2,9,,,,85,,
vsel,s,volu,,5,6,1
cm,bb,volu
vsel,s,volu,,9,12,1
cm,vv,volu
allsel,all
vsbv,bb,vv
allsel,all
numcmp,volu
vglue,1,2,3,4,7
numcmp,volu
vglue,2,3,4
numcmp,volu
!建螺栓
wpoffs,,,-lbt/2
cylind,dbt/2,,-lbt/2,lbt/2,0,360
cylind,dbth/2,,-lbt/2-lbth,-lbt/2,0,360
cylind,dbth/2,,lbt/2,lbt/2+lbth,0,360
vglue,9,10,11
numcmp,volu
vgen,2,9,11,1,,-hd/3,,
vglue,12,13,14
numcmp,volu
vgen,2,12,14,1,,-85,,
vglue,15,16,17
numcmp,volu
vgen,2,9,11,1,,85,,
vglue,18,19,20
numcmp,volu


!mesh螺栓
vsel,s,volu,,9,20,1
VATT,2,2,2
LSEL,s,radius,,dbth/2
LESIZE,all, , ,3, , , , ,1
LSEL,s,radius,,dbt/2
LESIZE,all, , ,3, , , , ,1
LSEL,s,length,,lbt
LESIZE,all, , ,2, , , , ,1  
LSEL,s,length,,lbth
LESIZE,all, , ,1, , , , ,1
VMESH,all

!mesh其他
allsel,all
vsel,u,volu,,9,20,1
VATT,1,1,1,0

LSEL,s,length,,hc-2*tcf
LESIZE,all, , ,10, , , , ,1
LSEL,s,length,,bc/2
LESIZE,all, , ,8, , , , ,1
LSEL,s,length,,(bc-tcw)/2
LESIZE,all, , ,7, , , , ,1
LSEL,s,length,,hb
LESIZE,all, , ,8, , , , ,1
LSEL,s,length,,tcf
LESIZE,all, , ,1, , , , ,1
LSEL,s,length,,tcw/2
LESIZE,all, , ,1, , , , ,1
LSEL,s,radius,,dh0/2
LESIZE,all, , ,3, , , , ,1
LSEL,s,radius,,dbth/2
LESIZE,all, , ,3, , , , ,1
LSEL,s,length,,tbf
LESIZE,all, , ,1, , , , ,1
LSEL,s,length,,tbw/2
LESIZE,all, , ,1, , , , ,1
LSEL,s,length,,bb/2
LESIZE,all, , ,8, , , , ,1
LSEL,s,length,,(bb-tbw)/2
LESIZE,all, , ,7, , , , ,1
LSEL,s,length,,hb-2*tbf
LESIZE,all, , ,8, , , , ,1
LSEL,s,length,,hd
LESIZE,all, , ,12, , , , ,1
LSEL,s,length,,td
LESIZE,all, , ,2, , , , ,1
LSEL,s,length,,bd
LESIZE,all, , ,10, , , , ,1
LSEL,s,length,,lc
LESIZE,all, , ,40, , , , ,1
LSEL,s,length,,lb
LESIZE,all, , ,30, , , , ,1

MSHAPE,1,3D
VMESH,all
!!!!建立螺栓与孔的接触!!!!!!!!
allsel,all
MP,MU,3,miu
ET,3,170 !目标单元为targe170
ET,4,174 !接触单元为targe174
!目标面
allsel,all
ASEL,S,,,65,72,1
ASEL,a,,,75,82,1
TYPE,3
mat,3
real,3   
NSLA,S,1
ESLN,S,0
ESURF,ALL
!接触面
allsel,all
ASEL,S,,,61,62,1
ASEL,a,,,21,22,1
ASEL,a,,,91,92,1
ASEL,a,,,103,104,1
TYPE,4
mat,3
real,3   
NSLA,S,1
ESLN,S,0
ESURF,ALL


!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!施加loadstep 1,边界约束条件和螺栓预拉力
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
allsel,all
finish
wpcsys,-1,0
/solu
antype,0
nlgeom,1
autots,on
nsubst,200,0,0
NEQIT,50   
CNVTOL,F, ,0.03,2, ,
stiff,1
lnsrch,1
pred,1
ALLSEL,ALL
SLOAD,ALL,9,LOCK,FORC,preten, 1,2
ASEL,S,LOC,Z,   
DA,all,SYMM
allsel,all
asel,s,loc,y,-lc/2
da,all,all
allsel,all
SFA,54,,PRES,3000
allsel,all
solve
save

xuruikl

solve不了
提示什么呢？