yxm04thu 发表于 2009-8-17 14:35:26

钢管混凝土柱,老是不收敛,望高手指教

将命令流和部分注释贴下,供大家参考,请高手指教,为什么一开始计算就出现严重的不收敛,算不下去。

FINISH
/clear
/PREP7
!!!单位制:N,mm,s
!!!-------------材料-------------!!!
!!!混凝土!!!
!!!---C80---!!!
MPTEMP,,,,,,,,
MP,EX,1,38000
MP,PRXY,1,0.2
MP,DENS,1,2.5E-9
TB,CONCR,1
TBDATA,,0.5,1.0,3.1,-1!定义砼C1,C2,R1,Ra
TB,MISO,1
TBPT,DEFI,200E-6,7.6
TBPT,DEFI,500E-6,17.2
TBPT,DEFI,1000E-6,32.5
TBPT,DEFI,1400E-6,42.1
TBPT,DEFI,1800E-6,48.3
TBPT,DEFI,2000E-6,49.8
TBPT,DEFI,2150E-6,50.2
TBPT,DEFI,2750E-6,42.4
TBPT,DEFI,3500E-6,29.1
TBPT,DEFI,4000E-6,22.9
TBPT,DEFI,5000E-6,15.4
TBPT,DEFI,6000E-6,11.3
TBPT,DEFI,7000E-6,8.8
!!!钢材!!!
!!!---Q235钢板---!!!
MP,EX,2,2E5
MP,PRXY,2,0.3
MP,DENS,2,7.85E-9
TB,BKIN,2
TBDATA,,235,2E3
!!!---HRB235钢筋---!!!
MP,EX,3,2E5
MP,PRXY,3,0.3
MP,DENS,3,7.85E-9
TB,BKIN,3
TBDATA,,235,2E3
!!!---刚体---!!!
MP,EX,4,1E8
MP,PRXY,4,0.2
MP,DENS,4,1E-12
!!!-------------单元-------------!!!
ET,1,SOLID65
KEYOPT,1,7,1
ET,2,SHELL181
ET,3,LINK8
!!!-------------实常数-------------!!!
R,1,6,6,6,6
R,2,5,5,5,5
R,3,3,3,3,3
R,4,30.8
R,5,,,,,,
!!!-------------几何模型-------------!!!
!!!---混凝土---!!!
K,1,228,156
K,2,228,-156
K,3,-228,-156
K,4,-228,156
K,5,228,156,1856
K,6,228,-156,1856
K,7,-228,-156,1856
K,8,-228,156,1856
V,1,2,3,4,5,6,7,8
!!!---刚体---!!!
K,9,228,156,1886
K,10,228,-156,1886
K,11,-228,-156,1886
K,12,-228,156,1886
V,5,6,7,8,9,10,11,12
!!!---切分纵隔板---!!!
WPOFFS,168
WPROTA,,,90
VSBW,ALL
WPOFFS,,,-336   
VSBW,ALL
WPROTA,,90
WPOFFS,,,96
VSBW,ALL
WPOFFS,,,-192   
VSBW,ALL
WPOFFS,,,96
VSBW,ALL
WPOFFS,,168
WPROTA,,90
VSBW,ALL
WPOFFS,,,114
VSBW,ALL
WPOFFS,,,-228   
VSBW,ALL
!!!---切分横隔板---!!!
ALLS
WPOFFS,-337
WPROTA,,,90
VSBW,ALL
WPOFFS,,,-333
VSBW,ALL
WPOFFS,,,-333
VSBW,ALL
WPOFFS,,,-333
VSBW,ALL
WPOFFS,,,-84
VSBW,ALL
WPOFFS,,,-256
VSBW,ALL
WPOFFS,,,-100
VSBW,ALL
!!!-------------单元属性-------------!!!
!!!---刚性层---!!!
FLST,5,24,6,ORDE,23
FITEM,5,3   
FITEM,5,7   
FITEM,5,9   
FITEM,5,12
FITEM,5,14
FITEM,5,16
FITEM,5,19
FITEM,5,21
FITEM,5,24
FITEM,5,-25
FITEM,5,27
FITEM,5,-28
FITEM,5,36
FITEM,5,41
FITEM,5,-43
FITEM,5,45
FITEM,5,47
FITEM,5,-48
FITEM,5,50
FITEM,5,-51
FITEM,5,53
FITEM,5,55
FITEM,5,-56
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,       4,       5,   1,       0
MSHAPE,0,3D
MSHKEY,1
ESIZE,60
VMESH,ALL
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1

!!!---混凝土---!!!
FLST,5,192,6,ORDE,13
FITEM,5,18
FITEM,5,20
FITEM,5,22
FITEM,5,30
FITEM,5,-35
FITEM,5,38
FITEM,5,-39
FITEM,5,57
FITEM,5,-96
FITEM,5,99
FITEM,5,-144
FITEM,5,146
FITEM,5,-240
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT,       1,       5,   1,       0
MSHAPE,0,3D
MSHKEY,1
ESIZE,60
VMESH,ALL
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!!!---钢管---!!!
ALLS
FLST,5,160,5,ORDE,156   
FITEM,5,81
FITEM,5,83
FITEM,5,91
FITEM,5,126
FITEM,5,-128
FITEM,5,133
FITEM,5,139
FITEM,5,147
FITEM,5,150
FITEM,5,158
FITEM,5,212
FITEM,5,214
FITEM,5,216
FITEM,5,220
FITEM,5,224
FITEM,5,230
FITEM,5,-231
FITEM,5,234
FITEM,5,236
FITEM,5,-238
FITEM,5,242
FITEM,5,244
FITEM,5,-246
FITEM,5,253
FITEM,5,258
FITEM,5,264
FITEM,5,270
FITEM,5,274
FITEM,5,277
FITEM,5,-278
FITEM,5,283
FITEM,5,291
FITEM,5,296
FITEM,5,298
FITEM,5,302
FITEM,5,308
FITEM,5,312
FITEM,5,315
FITEM,5,317
FITEM,5,320
FITEM,5,330
FITEM,5,-331
FITEM,5,354
FITEM,5,361
FITEM,5,-362
FITEM,5,365
FITEM,5,368
FITEM,5,371
FITEM,5,374
FITEM,5,-375
FITEM,5,378
FITEM,5,381
FITEM,5,383
FITEM,5,388
FITEM,5,391
FITEM,5,397
FITEM,5,400
FITEM,5,405
FITEM,5,407
FITEM,5,415
FITEM,5,417
FITEM,5,-418
FITEM,5,421
FITEM,5,424
FITEM,5,427
FITEM,5,-428
FITEM,5,430
FITEM,5,434
FITEM,5,443
FITEM,5,447
FITEM,5,454
FITEM,5,459
FITEM,5,465
FITEM,5,473
FITEM,5,479
FITEM,5,482
FITEM,5,487
FITEM,5,489
FITEM,5,493
FITEM,5,521
FITEM,5,-522
FITEM,5,525
FITEM,5,-526
FITEM,5,529
FITEM,5,532
FITEM,5,535
FITEM,5,538
FITEM,5,-539
FITEM,5,542
FITEM,5,545
FITEM,5,547
FITEM,5,552
FITEM,5,555
FITEM,5,561
FITEM,5,564
FITEM,5,569
FITEM,5,571
FITEM,5,575
FITEM,5,581
FITEM,5,584
FITEM,5,-585
FITEM,5,589
FITEM,5,-590
FITEM,5,595
FITEM,5,601
FITEM,5,-602
FITEM,5,606
FITEM,5,608
FITEM,5,614
FITEM,5,620
FITEM,5,622
FITEM,5,629
FITEM,5,634
FITEM,5,637
FITEM,5,643
FITEM,5,646
FITEM,5,651
FITEM,5,653
FITEM,5,681
FITEM,5,684
FITEM,5,-686
FITEM,5,689
FITEM,5,-690
FITEM,5,693
FITEM,5,696
FITEM,5,699
FITEM,5,702
FITEM,5,-703
FITEM,5,706
FITEM,5,709
FITEM,5,711
FITEM,5,716
FITEM,5,719
FITEM,5,725
FITEM,5,728
FITEM,5,733
FITEM,5,735
FITEM,5,742
FITEM,5,745
FITEM,5,750
FITEM,5,754
FITEM,5,761
FITEM,5,765
FITEM,5,773
FITEM,5,776
FITEM,5,783
FITEM,5,-784
FITEM,5,788
FITEM,5,791
FITEM,5,793
FITEM,5,798
FITEM,5,801
FITEM,5,807
FITEM,5,810
FITEM,5,815
FITEM,5,817
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT,       2,       1,   2,       0,   
MSHAPE,0,2D
MSHKEY,1
ESIZE,60
AMESH,ALL
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*

!!!---纵隔板---!!!
FLST,5,64,5,ORDE,64
FITEM,5,218
FITEM,5,222
FITEM,5,226
FITEM,5,248
FITEM,5,260
FITEM,5,280
FITEM,5,288
FITEM,5,290
FITEM,5,297
FITEM,5,307
FITEM,5,324
FITEM,5,367
FITEM,5,370
FITEM,5,373
FITEM,5,380
FITEM,5,385
FITEM,5,393
FITEM,5,396
FITEM,5,399
FITEM,5,426
FITEM,5,436
FITEM,5,450
FITEM,5,-451
FITEM,5,461
FITEM,5,468
FITEM,5,475
FITEM,5,478
FITEM,5,481
FITEM,5,531
FITEM,5,534
FITEM,5,537
FITEM,5,544
FITEM,5,549
FITEM,5,557
FITEM,5,560
FITEM,5,563
FITEM,5,578
FITEM,5,583
FITEM,5,598
FITEM,5,610
FITEM,5,617
FITEM,5,626
FITEM,5,631
FITEM,5,639
FITEM,5,642
FITEM,5,645
FITEM,5,695
FITEM,5,698
FITEM,5,701
FITEM,5,708
FITEM,5,713
FITEM,5,721
FITEM,5,724
FITEM,5,727
FITEM,5,744
FITEM,5,752
FITEM,5,756
FITEM,5,760
FITEM,5,764
FITEM,5,790
FITEM,5,795
FITEM,5,803
FITEM,5,806
FITEM,5,809
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT,       2,       2,   2,       0,   
MSHAPE,0,2D
MSHKEY,1
ESIZE,60
AMESH,ALL
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
!!!---横隔板---!!!
FLST,5,80,5,ORDE,80
FITEM,5,7   
FITEM,5,-8
FITEM,5,15
FITEM,5,20
FITEM,5,30
FITEM,5,38
FITEM,5,43
FITEM,5,48
FITEM,5,-49
FITEM,5,59
FITEM,5,63
FITEM,5,75
FITEM,5,78
FITEM,5,80
FITEM,5,90
FITEM,5,101
FITEM,5,118
FITEM,5,120
FITEM,5,131
FITEM,5,-132
FITEM,5,161
FITEM,5,185
FITEM,5,191
FITEM,5,211
FITEM,5,213
FITEM,5,219
FITEM,5,223
FITEM,5,227
FITEM,5,229
FITEM,5,241
FITEM,5,243
FITEM,5,249
FITEM,5,255
FITEM,5,-256
FITEM,5,261
FITEM,5,-262
FITEM,5,269
FITEM,5,272
FITEM,5,-273
FITEM,5,281
FITEM,5,-282
FITEM,5,293
FITEM,5,300
FITEM,5,-301
FITEM,5,311
FITEM,5,314
FITEM,5,316
FITEM,5,319
FITEM,5,329
FITEM,5,336
FITEM,5,340
FITEM,5,343
FITEM,5,-344
FITEM,5,352
FITEM,5,-353
FITEM,5,355
FITEM,5,412
FITEM,5,-413
FITEM,5,416
FITEM,5,422
FITEM,5,429
FITEM,5,440
FITEM,5,446
FITEM,5,449
FITEM,5,453
FITEM,5,455
FITEM,5,462
FITEM,5,-463
FITEM,5,467
FITEM,5,500
FITEM,5,504
FITEM,5,512
FITEM,5,515
FITEM,5,577
FITEM,5,586
FITEM,5,596
FITEM,5,603
FITEM,5,607
FITEM,5,618
FITEM,5,621
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT,       2,       3,   2,       0,   
MSHAPE,0,2D
MSHKEY,1
ESIZE,60
AMESH,ALL
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
FLST,5,32,5,ORDE,32
FITEM,5,12
FITEM,5,27
FITEM,5,56
FITEM,5,68
FITEM,5,89
FITEM,5,123
FITEM,5,202
FITEM,5,235
FITEM,5,263
FITEM,5,341
FITEM,5,348
FITEM,5,414
FITEM,5,452
FITEM,5,470
FITEM,5,495
FITEM,5,503
FITEM,5,513
FITEM,5,579
FITEM,5,594
FITEM,5,600
FITEM,5,612
FITEM,5,619
FITEM,5,660
FITEM,5,664
FITEM,5,672
FITEM,5,675
FITEM,5,747
FITEM,5,753
FITEM,5,763
FITEM,5,767
FITEM,5,775
FITEM,5,778
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT,       2,       1,   2,       0,   
MSHAPE,0,2D
MSHKEY,1
ESIZE,60
AMESH,ALL
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*

!!!-------------加载求解-------------!!!
/SOLU
ANTYPE,0
NLGEOM,ON
NROPT,FULL
CNVTOL,F, ,0.02,2,0.01,
CNVTOL,M, ,0.02,2,0.01,
AUTOTS,1 !打开自动时间步控制
LNSRCH,1 !打开线性搜索
NCNV,2 !如果不收敛时结束而不退出
NEQIT,100 !每一子步中方程的迭代次数限值
TIME,1
DELTIM,0.005,1e-10,0.1,
OUTRES,ALL,1
ALLS
NSEL,S,LOC,Z,0
CM,FIXED,NODE
D,FIXED,ALL
ALLS
FLST,5,7,1,ORDE,7   
FITEM,5,11
FITEM,5,16
FITEM,5,19
FITEM,5,24
FITEM,5,44
FITEM,5,-45
FITEM,5,71
NSEL,S, , ,P51X
CM,PUSH,NODE
ALLS
D,PUSH,UX,100
ALLS
SOLVE
页: [1]
查看完整版本: 钢管混凝土柱,老是不收敛,望高手指教