中华钢结构论坛 China Structure Forum

标题: 网壳顶稳定性分析为何不收敛的问题请高手赐教 [打印本页]

作者: 5350553 时间: 2010-6-3 13:44 标题: 网壳顶稳定性分析为何不收敛的问题请高手赐教

/clear
/FILENAME,m04

/TITLE,NONLINER Buckling STABILITY ANALYSIS OF SPHERICAL RETICULATED-SHELL(FULL)
/prep7

ALLSEL,ALL

ET,1,BEAM189
MP,DENS,1,7850
MP,EX,1,2e11
MP,PRXY,1,0.3

LOCAL,11,2,0,0,0
CSCIR,11,1
K,1,46,0,60
K,37,46,180,60
K,72,46,355,60
kFILL,1,37
kFILL,37,72

kGEN,2,1,72,1,0,0,4.23,72,0

kGEN,2,1,72,1,0,0,8.71,144,0

kGEN,2,1,72,1,0,0,13.32,216,0

kGEN,2,1,72,1,0,0,17.95,288,0

kGEN,2,1,72,1,0,0,22.51,360,0

kGEN,2,1,72,1,0,0,26.88,432,0

*DO,m,1,7,1
J=72*(m-1)+1
*DO,I,J,J+70,1
L,I,I+1
*ENDDO
L,J+71,J
*ENDDO

*DO,m,1,6,1
J=72*(m-1)+1
*DO,I,J,J+71,1
L,I,I+72
*ENDDO
*ENDDO
LDIV,505,0.388, ,2,0
LDIV,523,0.388, ,2,0
LDIV,541,0.388, ,2,0
LDIV,559,0.388, ,2,0

LSTR,    439,    463
LSTR,    475,    499
LSTR,    445,    493
LSTR,    457,    481

L,362,289
L,432,289
L,289,218
L,289,288
L,218,145
L,288,145
L,145,74
L,145,144
L,74,505
L,144,505

L,378,307
L,380,307
L,307,234
L,307,236
L,234,163
L,236,163
L,163,90
L,163,92
L,90,506
L,92,506

L,396,325
L,398,325
L,325,252
L,325,254
L,252,181
L,254,181
L,181,108
L,181,110
L,108,507
L,110,507

L,414,343
L,416,343
L,343,270
L,343,272
L,270,199
L,272,199
L,199,126
L,199,128
L,126,508
L,128,508
LOVLAP,941,942,943,944
LDELE,866, 936,2 ,1

底圈梁
ALLSEL,ALL

LSEL,S,LINE,,1,72,,0

CM,Line1,LINE
!R2-R5环梁

ALLSEL,ALL

LSEL,S,LINE,,73,360,,0

CM,Line2,LINE
!R6环梁
ALLSEL,ALL

LSEL,S,LINE,,361,432,,0

CM,Line3,LINE
!R7中心圈环梁
ALLSEL,ALL

LSEL,S,LINE,,433,504,,0

CM,Line4,LINE

!中心加强筋

ALLSEL,ALL
LSEL,S,LINE,,985,996,,0
CM,Line5,LINE

!拱梁
ALLSEL,ALL

LSEL,S,LINE,,505,864,,0
LSEL,A,LINE,,865,935,2,0
LSEL,A,LINE,,937,940,,0
CM,Line6,LINE

!斜撑
ALLSEL,ALL

LSEL,S,LINE,,945,984,,0
CM,Line7,LINE

!底圈环梁
SECTYPE,1,BEAM,MESH,beam007
SECOFFSET,USER,0.04,0.297,
SECREAD,'BEAM007','SECT',' ',MESH

!R2-R5环梁
SECTYPE, 2, BEAM, CHAN, beam02, 0
SECOFFSET, CENT
SECDATA,0.077,0.077,0.220,0.0115,0.0115,0.007,0,0,0,0
!R6环梁
SECTYPE, 3, BEAM, I, beam03, 0
SECOFFSET, CENT
SECDATA,0.175,0.175,0.244,0.011,0.011,0.007,0,0,0,0

!R7中间环梁
SECTYPE, 4, BEAM, I, beam04, 0
SECOFFSET, CENT
SECDATA,0.3,0.3,0.328,0.014,0.014,0.014,0,0,0,0

!拱肋
SECTYPE, 5, BEAM, I, beam05, 0
SECOFFSET, CENT
SECDATA,0.15,0.15,0.3,0.009,0.009,0.0065,0,0,0,0

!斜撑
SECTYPE, 6, BEAM, L, beam06, 0
SECOFFSET, CENT
SECDATA,0.1,0.1,0.01,0.01,0,0,0,0,0,0
ALLSEL,ALL

k,513,0,0,0

LESIZE,ALl, , ,6, , , , ,1

!底圈环梁
CMSEL,S,Line1

LATT, , ,1, ,    , ,1

LMESH,All
ALLSEL,ALL

!R2-R5环梁
CMSEL,S,Line2

LATT, , ,1, ,    , ,2

LMESH,All
ALLSEL,ALL
!R6环梁
CMSEL,S,Line3

LATT, , ,1, ,    , ,3

LMESH,All
ALLSEL,ALL

!R7中间环梁

CMSEL,S,Line4

LATT, , ,1, ,    , ,4

LMESH,All
ALLSEL,ALL

!中间加强筋
CMSEL,S,Line5

LATT, , ,1, ,  513  , ,5

LMESH,All
ALLSEL,ALL
!拱肋

CMSEL,S,Line6

LATT, , ,1, ,  513  , ,5

LMESH,All
ALLSEL,ALL

!斜撑

CMSEL,S,Line7

LATT, , ,1, ,  513  , ,6

LMESH,All
ALLSEL,ALL

ALLSEL,ALL

KSEL,S,KP,,1,72,,0

CM,KP1,KP
!底圈梁R1

/SOl
CSYS,0

CMSEL,S,Line1

DL,all, ,all,
!DL,all, ,UY,
!DL,all, ,UZ,

ALLSEL,ALL

FK,ALL,FZ,-1

ALLSEL,ALL
/SOLU
ANTYPE,0
EQSLV,SPAR
PSTRES,ON
SOLVE
FINISH

/SOL
ANTYPE,1
BUCOPT,SUBSP,1,0,0
SUBOPT,0,0,0,0,0,ALL
OUTPR,BASIC,ALL,
MXPAND,1,0,0,1,0.001,
SOLVE
FINISH
/POST1
*get,freq1,mode,1,freq
FINISH
save

/PREP7

UPGEOM,0.153,1,1,'m04',rst

save,pert,db
resu,pert,db

SAVE
FINISH



/solu
outres,1
antype,0
nlgeom,on

outres,all,all
arclen,on,4


nsubst,30
ALLSEL,ALL

FK,ALL,FZ,-freq1*1.2

ALLSEL,ALL

/GST                                  !将竖向力值放大，放大系数为第一阶模态的主频
solve                                                 !
finish

欢迎光临中华钢结构论坛 China Structure Forum (http://bbs.okok.org/)