目前在做砌体在地震下的裂缝开展，但是一直达不到想要的结果即主拉应力的剪切破坏以斜裂缝为主，以下是我的命令流，请各位帮帮忙，不胜感激！
/prep7
et,1,solid65!砌体定义
KEYOPT,1,5,1!输出每个积分点的解
KEYOPT,1,6,3!输出质心的解和积分点的解
KEYOPT,1,7,1!考虑应力松弛,有助于计算收敛
r,1
mp,ex,1,2.0e9
mp,prxy,1,0.15
mp,dens,1,2160
mp,alpx,1,5e-6
tb,miso,1,1,8
tbpt,,0.0006,1.2e6
tbpt,,0.0012,2.2e6
tbpt,,0.0018,2.7e6
tbpt,,0.0024,3.0e6
tbpt,,0.0030,3.2e6
tbpt,,0.0036,3.3e6
tbpt,,0.0042,3.3e6
tbpt,,0.0048,3.3e6
tb,concr,1,1,9
tbdata,,0.15,0.9,2.1e5,-1
  
*SET,fuconcr,1.43e7
*SET,ftcon,1.43e6
et,2,solid65!混凝土门窗过梁
KEYOPT,1,5,1!输出每个积分点的解
KEYOPT,1,6,3!输出质心的解和积分点的解
KEYOPT,1,7,1!考虑应力松弛，有助于计算收敛
r,2
mp,ex,2,3.06e10!C30
mp,prxy,2,0.2
mp,dens,2,2550
mp,alpx,2,1e-5
tb,mkin,2,1
tbtemp,,strain
tbdata,1,0.4*fuconcr/(3.06e10),0.001,0.0015,0.002,0.0035!应变
tbtemp,0.0
tbdata,1,fuconcr*0.4,fuconcr*0.75,fuconcr*15/16,fuconcr,fuconcr!应力
tb,concr,2,1,9
tbdata,,0.5,0.9,ftcon,-1
  
et,3,solid45!混凝土楼盖定义
mp,ex,3,3.06e10!C30
mp,prxy,3,0.2
mp,dens,3,2800!做法折算成自重
mp,alpx,3,1e-5
  
/VIEW,1,1,1,1
block,0,0.24,0,3,0,0.24
lesize,7,,,2,,,,,0
lesize,1,,,30,,,,,0
lesize,11,,,2,,,,,0
mat,1
type,1
real,1
vmesh,1
  
block,0.24,1.04,0,3,0,0.24
lesize,14,,,8,,,,,0
lesize,18,,,30,,,,,0
lesize,22,,,2,,,,,0
mat,1
type,1
real,1
vmesh,2
  
block,2.64,3.44,0,3,0,0.24
lesize,29,,,8,,,,,0
lesize,30,,,30,,,,,0
lesize,34,,,2,,,,,0
mat,1
type,1
real,1
vmesh,3
  
block,1.04,2.64,2.8,3,0,0.24
lesize,46,,,2,,,,,0
lesize,41,,,16,,,,,0
lesize,42,,,2,,,,,0
mat,1
type,1
real,1
vmesh,4
  
block,1.04,2.64,2.4,2.8,0,0.24!窗过梁
lesize,54,,,4,,,,,0
lesize,53,,,16,,,,,0
lesize,58,,,2,,,,,0
mat,2
type,2
real,2
vmesh,5
  
block,1.04,2.64,0,0.9,0,0.24
lesize,61,,,9,,,,,0
lesize,64,,,16,,,,,0
lesize,72,,,2,,,,,0
mat,1
type,1
real,1
vmesh,6
  
block,3.44,3.68,0,3,0,0.24
lesize,77,,,2,,,,,0
lesize,75,,,30,,,,,0
lesize,82,,,2,,,,,0
mat,1
type,1
real,1
vmesh,7
  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
  
block,3.44,3.68,0,3,0.24,3.54!*****右侧横墙*****
lesize,74,,,30,,,,,0
lesize,78,,,33,,,,,0
lesize,75,,,2,,,,,0
mat,1
type,1
real,1
vmesh,8
  
block,3.44,3.68,0,3,3.54,3.78
lesize,87,,,2,,,,,0
lesize,86,,,30,,,,,0
lesize,90,,,2,,,,,0
mat,1
type,1
real,1
vmesh,9
  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
  
block,0,0.24,0,3,0.24,3.54!*****左侧横墙*****
lesize,90,,,30,,,,,0
lesize,94,,,33,,,,,0
lesize,89,,,2,,,,,0
mat,1
type,1
real,1
vmesh,10
  
block,0,0.24,0,3,3.54,3.78
lesize,101,,,2,,,,,0
lesize,102,,,30,,,,,0
lesize,106,,,2,,,,,0
mat,1
type,1
real,1
vmesh,11
  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
  
block,0.24,1.34,0,3,3.54,3.78  !*****前纵墙*****
lesize,105,,,11,,,,,0
lesize,103,,,30,,,,,0
lesize,110,,,2,,,,,0
mat,1
type,1
real,1
vmesh,12
  
block,2.34,3.44,0,3,3.54,3.78  
lesize,114,,,11,,,,,0
lesize,120,,,30,,,,,0
lesize,124,,,2,,,,,0
mat,1
type,1
real,1
vmesh,13
  
block,1.34,2.34,2.2,3,3.54,3.78
lesize,130,,,8,,,,,0
lesize,129,,,10,,,,,0
lesize,134,,,2,,,,,0
mat,1
type,1
real,1
vmesh,14
  
block,1.34,2.34,2.1,2.2,3.54,3.78
lesize,142,,,1,,,,,0
lesize,141,,,10,,,,,0
lesize,146,,,2,,,,,0
mat,2
type,2
real,2
vmesh,15
  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
  
block,0,0.24,3.0,3.1,0,0.24!*****屋盖*****
lesize,141,,,2,,,,,0
lesize,145,,,2,,,,,0
lesize,137,,,1,,,,,0
mat,3
type,3
vmesh,16
  
block,0.24,3.44,3.0,3.1,0,0.24
lesize,148,,,32,,,,,0
lesize,157,,,2,,,,,0
lesize,149,,,1,,,,,0
mat,3
type,3
vmesh,17
  
block,3.44,3.68,3.0,3.1,0,0.24
lesize,165,,,2,,,,,0
lesize,169,,,2,,,,,0
lesize,161,,,1,,,,,0
mat,3
type,3
vmesh,18
  
block,3.44,3.68,3.0,3.1,0.24,3.54
lesize,182,,,33,,,,,0
lesize,177,,,2,,,,,0
lesize,176,,,1,,,,,0
mat,3
type,3
vmesh,19
  
block,3.44,3.68,3.0,3.1,3.54,3.78
lesize,193,,,2,,,,,0
lesize,189,,,2,,,,,0
lesize,188,,,1,,,,,0
mat,3
type,3
vmesh,20
  
block,0,0.24,3.0,3.1,0.24,3.54
lesize,206,,,33,,,,,0
lesize,201,,,2,,,,,0
lesize,200,,,1,,,,,0
mat,3
type,3
vmesh,21
  
block,0,0.24,3.0,3.1,3.54,3.78
lesize,213,,,2,,,,,0
lesize,217,,,2,,,,,0
lesize,212,,,1,,,,,0
mat,3
type,3
vmesh,22
  
block,0.24,3.44,3.0,3.1,3.54,3.78
lesize,223,,,32,,,,,0
lesize,229,,,2,,,,,0
lesize,221,,,1,,,,,0
mat,3
type,3
vmesh,23
  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
  
block,0.24,3.44,3.0,3.1,0.24,3.54
lesize,190,,,33,,,,,0
lesize,185,,,32,,,,,0
lesize,184,,,1,,,,,0
mat,3
type,3
vmesh,24
  
block,2.64,2.88,0,3,0.24,3.54
lesize,193,,,30,,,,,0
lesize,200,,,33,,,,,0
lesize,192,,,2,,,,,0
mat,1
type,1
real,1
vmesh,25
  
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
  
nsel,s,loc,y,0
d,all,all
allsel,all
  
nt=150
*dim,time,array,nt
*vread,time,time,txt,,ijk,nt
(f6.2)
*dim,acez,array,nt
*vread,acez,acez,txt,,ijk,nt
(f10.6)
  
/solu
antype,trans
trnopt,full
timint,off
time,1e-8
nsubst,2
kbc,1
acel,0,9.8,0
solve
nlgeom,on
timint,on
autots,on
OUTRES,NSOL,last
OUTRES,RSOL,last
OUTRES,ESOL,last
OUTRES,STRS,last
OUTRES,EPEL,last
ALPHAD,6.40
BETAD,0.00035
neqit,4
cnvtol,f,,0.01,2
ncnv,0
*do,i,1,nt
time,time(i)
kbc,0
nsubst,1,,,1
acel,0,0,0.03*acez(i)
solve
*enddo
finish

 

 

 看看哪里不对，以下是裂缝图，困惑的是都是水平裂缝，怎么才能出现斜裂缝或者交叉裂缝。