土与上部结构共同作用命令流！瞬态分析不收敛

mylife2000

！以下为本人建立的土与上部结构共同作用命令流。瞬态分析不收敛。望各位指点。计算时候将tim、dizhenbo文件考到计算目录下。
/title,6layers_frame_base_on_board
/prep7
ET,1,BEAM4
R,1,0.36,1.08e-2,1.08e-2,0.6,0.6!一到二层柱实常数
R,2,0.25,5.2e-3,5.2e-3,0.5,0.5!三到六层柱实常数
R,3,0.18,1.35e-3,0.54e-2,0.3,0.6!标准层主梁实常数
R,4,0.195,1.46e-3,0.69e-2,0.3,0.65!顶层主梁实常数
R,5,0.125,6.51e-4,2.6e-3,0.25,0.5!次梁实常数
R,6,0.56,2.29e-2,2.99e-2,0.7,0.8!基础梁实常数
ET,2,SHELL63
R,7,0.11!板
ET,3,solid45
!定义材料属性
mp,ex,1,3e10
mp,nuxy,1,0.2
mp,dens,1,2700
mp,ex,2,6.5e8
mp,nuxy,2,0.25
mp,dens,2,1500
tb,dp,2
tbdata,1,9500,30,30
!定义一层关键点
k,1,0,0,0
k,2,3.6,0,0
k,3,3.6*2,0,0
k,4,3.6*3,0,0
k,5,3.6*4,0,0
kgen,3,1,5,,,3.6
!层数及层高在此定义
!可在此变化模型
x=6
y=3.9
!节点层间复制
kgen,6,1,15,,,,y
kgen,2,1,5,2,,,-y
kgen,2,11,15,2,,,-y
!生成各层梁网格
*do,i,1,x
l,1+15*(i-1),2+15*(i-1)
l,2+15*(i-1),3+15*(i-1)
l,3+15*(i-1),4+15*(i-1)
l,4+15*(i-1),5+15*(i-1)
l,6+15*(i-1),7+15*(i-1)
l,7+15*(i-1),8+15*(i-1)
l,8+15*(i-1),9+15*(i-1)
l,9+15*(i-1),10+15*(i-1)
l,11+15*(i-1),12+15*(i-1)
l,12+15*(i-1),13+15*(i-1)
l,13+15*(i-1),14+15*(i-1)
l,14+15*(i-1),15+15*(i-1)
l,1+15*(i-1),6+15*(i-1)
l,6+15*(i-1),11+15*(i-1)
l,2+15*(i-1),7+15*(i-1)
l,7+15*(i-1),12+15*(i-1)
l,3+15*(i-1),8+15*(i-1)
l,8+15*(i-1),13+15*(i-1)
l,4+15*(i-1),9+15*(i-1)
l,9+15*(i-1),14+15*(i-1)
l,5+15*(i-1),10+15*(i-1)
l,10+15*(i-1),15+15*(i-1)
*enddo
!创建二层及以上层柱网格
*do,i,1,x-1
l,1+15*(i-1),1+15*i
l,3+15*(i-1),3+15*i
l,5+15*(i-1),5+15*i
l,11+15*(i-1),11+15*i
l,13+15*(i-1),13+15*i
l,15+15*(i-1),15+15*i
*enddo
!创建一层柱网格
l,1,1+15*x
l,3,2+15*x
l,5,3+15*x
l,11,4+15*x
l,13,5+15*x
l,15,6+15*x
!创建各层楼板
*do,i,1,x
a,1+15*(i-1),2+15*(i-1),7+15*(i-1),6+15*(i-1)
a,2+15*(i-1),3+15*(i-1),8+15*(i-1),7+15*(i-1)
a,3+15*(i-1),4+15*(i-1),9+15*(i-1),8+15*(i-1)
a,4+15*(i-1),5+15*(i-1),10+15*(i-1),9+15*(i-1)
a,6+15*(i-1),7+15*(i-1),12+15*(i-1),11+15*(i-1)
a,7+15*(i-1),8+15*(i-1),13+15*(i-1),12+15*(i-1)
a,8+15*(i-1),9+15*(i-1),14+15*(i-1),13+15*(i-1)
a,9+15*(i-1),10+15*(i-1),15+15*(i-1),14+15*(i-1)
*enddo
!标准层主梁网格划分为6个单元
*do,i,1,x-1
lsel,s,,,1+22*(i-1),4+22*(i-1)
lsel,a,,,9+22*(i-1),12+22*(i-1)
lsel,a,,,13+22*(i-1),14+22*(i-1)
lsel,a,,,17+22*(i-1),18+22*(i-1)
lsel,a,,,21+22*(i-1),22+22*(i-1)
latt,1,3,1
lesize,all,,,6
lmesh,all
*enddo
!顶层主梁网格划分为6个单元
*do,i,x,x
lsel,s,,,1+22*(i-1),4+22*(i-1)
lsel,a,,,9+22*(i-1),12+22*(i-1)
lsel,a,,,13+22*(i-1),14+22*(i-1)
lsel,a,,,17+22*(i-1),18+22*(i-1)
lsel,a,,,21+22*(i-1),22+22*(i-1)
latt,1,4,1
lesize,all,,,6
lmesh,all
*enddo
!次梁网格划分为6个单元
*do,i,1,x
lsel,s,,,5+22*(i-1),8+22*(i-1)
lsel,a,,,15+22*(i-1),16+22*(i-1)
lsel,a,,,19+22*(i-1),20+22*(i-1)
latt,1,5,1
lesize,all,,,6
lmesh,all
*enddo
!楼板网格划分
asel,all
aatt,1,7,2
amesh,all
allsel,all
!二层柱网格划分为13个单元
lsel,s,,,1+22*x,6+22*x
latt,1,1,1
lesize,all,,,13
lmesh,all
!一层柱网格划分
lsel,s,,,1+6*(x-1)+22*x,6+6*(x-1)+22*x
latt,1,1,1
lesize,all,,,13
lmesh,all
!三到四层柱网格划分
*do,i,1,2
lsel,s,,,1+6*i+22*x,6+6*i+22*x
latt,1,2,1
lesize,all,,,13
lmesh,all
*enddo
!其余层柱网格划分
!若增加一层可将2变为3其余类推
*do,i,1,2
lsel,s,,,1+6*(i+2)+22*x,6+6*(i+2)+22*x
latt,1,2,1
lesize,all,,,13
lmesh,all
*enddo
nummrg,all, , , ,low
numcmp,all
allsel,all
/replot
!正负零以下模型如下
!定义基础埋深
d=1.8
block,-0.3,0.3,-0.3,0.3,-y-d+0.6,-y
vgen,3,1,,,7.2
vgen,2,1,3,,,7.2
block,0.3,6.9,-0.3,0.3,-y-d,-y-d+0.6
vgen,2,7,,,7.2
vgen,2,7,8,,,7.2
block,-0.3,0.3,0.3,6.9,-y-d,-y-d+0.6
vgen,3,11,,,7.2
block,0.3,6.9,0.3,6.9,-y-d,-y-d+0.6
vgen,2,14,,,7.2
block,-0.3,0.3,-0.3,0.3,-y-d,-y-d+0.6
vgen,3,16,,,7.2
vgen,2,16,18,,,7.2
block,-0.9,15.3,-0.9,8.1,-y-d-0.6,-y-d
block,-0.9,-0.3,-0.9,8.1,-y-d,-y
vgen,2,23,,,15.6
block,-0.3,14.7,-0.9,-0.3,-y-d,-y
vgen,2,25,,,,8.4
block,-0.3,14.7,0.3,6.9,-y-d+0.6,-y
block,0.3,6.9,-0.3,0.3,-y-d+0.6,-y
vgen,2,28,,,7.2
vgen,2,28,29,,,7.2
vglue,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,24,25,26,27,28,29,30,31
nummrg,all
numcmp,all
vsel,s,,,1,8
vsel,a,,,19,22
aslv,s
lsla,s
lsel,r,loc,z,-y
lesize,all,,,2
vsel,s,,,1,8
vsel,a,,,19,22
aslv,s
lsla,s
lsel,r,loc,z,-y-d
lesize,all,,,2
vsel,s,,,1,8
vsel,a,,,19,22
aslv,s
lsla,s
lsel,r,loc,z,-y-d+0.6
lesize,all,,,2
lsel,s,ndiv,,0
lesize,all,0.6
allsel,all
vsel,s,,,1,8
vsel,a,,,19,31
vatt,1,,3
vmesh,all
vsel,s,,,9,18
vatt,2,,3
vmesh,all
extopt,esize,4,0
extopt,aclear,0
type,3
mat,2
asel,s,loc,z,-y-d-0.6
vext,all,,,,,-4
extopt,esize,3,0
extopt,aclear,0
type,3
mat,2
asel,s,loc,y,-0.9
vext,all,,,,-3.6
extopt,esize,3,0
extopt,aclear,0
type,3
mat,2
asel,s,loc,y,8.1
vext,all,,,,3.6
extopt,esize,3,0
extopt,aclear,0
type,3
mat,2
asel,s,loc,x,-0.9
vext,all,,,-3.6
extopt,esize,3,0
extopt,aclear,0
type,3
mat,2
asel,s,loc,x,15.3
vext,all,,,3.6
allsel,all
cp,1,ux,2023,2406
cp,2,uy,2023,2406
cp,3,uz,2023,2406
cp,4,ux,2062,2487
cp,5,uy,2062,2487
cp,6,uz,2062,2487
cp,7,ux,2036,2433
cp,8,uy,2036,2433
cp,9,uz,2036,2433
cp,10,ux,2075,2514
cp,11,uy,2075,2514
cp,12,uz,2075,2514
cp,13,ux,2088,2541
cp,14,uy,2088,2541
cp,15,uz,2088,2541
cp,16,ux,2049,2460
cp,17,uy,2049,2460
cp,18,uz,2049,2460
allsel,all
/sol
asel,s,loc,x,-4.5
asel,a,loc,x,18.9
da,all,all
asel,s,loc,y,-4.5
asel,a,loc,y,11.7
da,all,all
asel,s,loc,z,-y-d-0.6-4
da,all,all
allsel,all
nsel,s,,,2023
nsel,a,,,2062
nsel,a,,,2036
nsel,a,,,2075
nsel,a,,,2049
nsel,a,,,2088
d,all,rotx
d,all,roty
d,all,rotz
allsel,all
!模态分析
antype,modal
modopt,subsp,10
solve
finish
/post1
set,list
!到该步模态分析完毕且给出了前10阶振型圆频率
/sol
!施加静力面荷载
!首先施加各标准层面荷载
asel,s,,,1,8*(x-1)
sfa,all,,pres,-5.4e3
!屋面面荷载
asel,s,,,8*(x-1)+1,8*x
sfa,all,,pres,-8.2e3
allsel,all
!施加楼层梁上荷载
esel,s,real,,3
sfbeam,all,1,pres,1.1e4
!施加屋顶女儿墙栏板荷载
esel,s,real,,4
sfbeam,all,1,pres,3.5e3
allsel,all
!施加地震波
*create,byb,mac
*dim,aa,,1000
*dim,bb,,1000
*vread,aa(1),dizhenbo,txt
(f8.3)
*vread,bb(1),time,txt
(f8.2)
*end
byb
!瞬态分析
k=1
antype,trans
trnopt,full
!根据计算结果应改正
alphad,0.056425
betad,0.044246
upcoord,1,on
nsubst,1
kbc,0
outres,all,1
*do,i,1,10,1
time,bb(i)
acel,,k*aa(i),9.8
lswrite,i
*enddo
lssolve,1,10,1
finish

[ 本帖最后由 xuruikl 于 2008-12-1 17:21 编辑 ]

mylife2000

请赐教!调了一个月了，是在没招了

jiegou2011

mylife2000 发表于 2008-12-1 17:13

                               
登录/注册后可看大图

请赐教!调了一个月了，是在没招了

刚看到你这个问题，我想知道你这问题怎么解决的？