[轻木砼混合结构]模态分析做不出！！！

sturdy8353

命令流如下：三层框架[混凝土+木结构]，二三层木结构，木结构剪力墙等效为等效桁架模型，即刚性杆+非线性弹簧做非线性动力分析，但是先做模态分析却始终做不出来，模态分析是对线性结构而言，则忽略非线性成分。不知道从何下手了。各位高手帮忙分析一下。谢谢。

/prep7
et,1,BEAM188   !框架梁柱
et,2,SHELL143   !混凝土板
et,3,LINK8    !刚性杆
et,4,COMBIN39   !非线性弹簧
et,5,MASS21,,,2   !集中质量
mp,ex,1,3.15e4  !C35混凝土线弹性参数
mp,prxy,1,0.2
mp,dens,1,2500e3
TB,BKIN,1,1,2,1
TBTEMP,0
TBDATA,,16.7,3.15e3
mp,ex,2,2e8     !刚性杆材料属性
mp,prxy,2,0.2
TB,BKIN,2,1,2,1
TBTEMP,0
TBDATA,,10,2e7
sectype,1,beam,rect  !框架柱截面
secdata,270,270
sectype,2,beam,rect  !X向框架梁截面
secdata,200,300
sectype,3,beam,rect  !Y向框架梁截面
secdata,250,250
  
r,1,80    !楼板厚度
r,2,62500 !刚性杆截面面积
r,3,0.0918,256.67,0.348,262.75,1.1718,660.08
rmore,2.9298,1179.17,4.2666,1467.58,5.676,1749.58,7.3242,1980.25
rmore,9.6678,2191.75,12.7626,2441.75,15.6558,2608.33,19.446,2819.83
rmore,23.7852,2986.42,26.55,3082.58,32.355,3146.67,38.5068,3101.83
rmore,40.851,3044.17,42.114,3012.08,45.9414,2903.17,50.4456,5698.08
rmore,51.288,2666
r,4,0.0918,256.67,0.4026,961.67,0.7692,1730
rmore,1.3548,2563.33,2.1606,3589.17,3.0396,4742.5,4.23,6280.83
rmore,5.2554,7626.67,6.3174,9036.67,8.1114,9869.17,11.5908,10766.67
rmore,16.5342,12305,22.98,13330,28.821,14163.33,34.3686,13970.83
rmore,40.7778,14035,46.8936,13715,48.816,13458.33,58.557,13394.17
rmore,60.8826,13137.5
r,5,0.0552,705,0.5124,4294.17,0.9702,6345
rmore,1.2816,7498.33,1.6296,8651.67,2.1972,10382.5,2.8566,12112.5
rmore,3.4242,13394.17,4.3944,15252.5,4.9074,16150,6.0972,18072.5
rmore,7.5624,19930.83,9.5946,20892.5,10.272,20956.67,13.2384,23071.67
rmore,16.443,22943.37,19.4826,21597.5,22.119,20507.5,26.6604,19802.5
rmore,31.7874,17431.67
r,6,,,875  !二层顶集中质量
r,7,,,275  !三层顶集中质量
k,5000,2915,1750,20000   !设置BEAM188单元方向
k,1,0,0,0                !底层关键点
k,2,5830,0,0
kgen,2,1,2,,,3500        !整体关键点
kgen,2,1,4,,,,3000
kgen,2,1,4,,,,5400
kgen,2,1,4,,,,7800
/view,1,1,1,1            !变换视角
/ANGLE,1,270,xm,0
/REPLOT
allsel
lsel,none
l,1,5   !框架柱模型
l,2,6
l,3,7
l,4,8
latt,1,,1,,5000,,1
lesize,all,1500
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,5,6   !X向框架梁模型
l,7,8
latt,1,,1,,5000,,2
lesize,all,3000
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,5,7   !Y向框架梁模型
l,6,8
latt,1,,1,,5000,,3
lesize,all,2000
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,5,9  !竖向刚性杆模型
l,6,10
l,7,11
l,8,12
lgen,2,9,12,1,,,2400,4
latt,2,2,3
lesize,all,1200
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,9,10  !二三层顶X向刚性杆模型
l,11,12
l,13,14
l,15,16
latt,2,2,3
lesize,all,3000
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,9,11  !二三层顶Y向刚性杆模型
l,10,12
l,13,15
l,14,16
latt,2,2,3
lesize,all,2000
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,9,12  !二三层顶楼面交叉刚性杆模型
l,10,11
l,13,16
l,14,15
latt,2,2,3
lesize,all,3500
lmesh,all
lsel,u,,,all
allsel
lsel,none
l,6,12  !侧面交叉刚性杆
l,8,10
l,5,11
l,7,9
l,9,15
l,10,16
l,11,13
l,12,14
latt,2,2,3
lesize,all,2000
lmesh,all
lsel,u,,,all
allsel
nsel,none
a,5,6,8,7
lesize,all,,,6
amesh,all
type,2
real,1
mat,1
asel,u,,,all

kmesh,9,12,1
type,5
real,6
kmesh,13,16,1
type,5
real,7
l,5,10
l,6,9
l,7,12
l,8,11
latt,,3,4
lesize,all,,,1
lmesh,all
lsel,u,,,all
l,9,14  
l,10,13
l,11,16
l,12,15
latt,,4,4
lesize,all,,,1
lmesh,all
lsel,u,,,all

[ 本帖最后由 rodge 于 2008-1-27 11:22 编辑 ]

sturdy8353

自己已经做出来了