杆系单元悬索桥模型不收敛
最近用杆系单元建了个悬索桥的模型，算自重的时候出现很多警告There are 399 small equation solver pivot terms.
Coefficient ratio exceeds 1.0e8 - Check results.
自振特性可以计算，但明显不正确，很多索夹处的局部振型。
附上命令流，期待解答。
FINISH
/CLEAR
/PREP7
!!!!!!!!!定义单元类型
ET,1,4
ET,2,10
ET,3,21
KEYOPT,2,3,0
!!!定义质量单元自由度
KEYOPT,3,3,2
!!!!!!!!!定义材料特性
!!!钢混结合主梁&钢塔柱Q345qD
MP,EX,1,2.06E11
MP,PRXY,1,0.3
MP,DENS,1,7850
MP, ALPX, 1, 1.2E-5
MP, DAMP, 1, 0.02
!!!锚跨砼梁&塔柱下部C50
MP,EX,2,3.45E10
MP,PRXY,2,0.2
MP,DENS,2,2549
MP, ALPX, 2, 1.0E-5
MP, DAMP, 2, 0.05
!!!主缆高强钢丝1770
MP,EX,3,2.05E11
MP,PRXY,3,0.3
MP,DENS,3,8005
MP,ALPX,3,1.2E-5
MP,DAMP,3,0.02
!!!吊索钢绞线1860
MP,EX,4,1.95E11
MP,PRXY,4,0.3
MP,DENS,4,8005
MP,ALPX,4,1.2E-5
MP,DAMP,4,0.02
!!!刚臂单元无密度
MP,EX,5,2E17
MP,PRXY,5,0.2
MP,DENS,5,0
!!!!!!!!!定义实常数
!!!!!!!!!质量单元实常数（砼梁二期恒载质量实常数在后）
R,1,163814.3461!吊点1处
R,2,163974.046 !其他吊点and主塔处
R,3,2326.7 !SJ1
R,4,1832.06 !SJ2
R,5,1634.45 !SJ3
R,6,1432.41 !SJ4
R,7,1015.78 !SJ5
R,8,903.99 !SJ6
R,9,447.9 !SJ7
!!!主缆、吊索及主梁初应变变量定义，方便修改
R10=-2.84E-04
R22=0.002587594
R23=0.002600247
R24=0.002627051
R25=0.002656878
R26=0.002689586
R27=0.002725067
R28=0.002762838
R29=0.002803213
R30=0.002847252
R31=0.00289326
R32=0.002940602
R33=0.002990587
R34=0.003042942
R35=0.00306963
R36=0.002967126
R37=0.0029415
R38=0.002893092
R39=0.002847611
R40=0.002804122
R41=0.002763808
R42=0.002726541
R43=0.002691445
R44=0.002658607
R45=0.002629285
R46=0.002602979
R47=0.002579752
R48=0.002559728
R49=0.00254283
R50=0.002529217
R51=0.002518996
R52=0.002512219
R53=0.002508842
R54=1.29E-04
R55=3.08E-04
R56=3.80E-04
R57=4.57E-04
R58=5.38E-04
R59=6.15E-04
R60=7.01E-04
R61=8.11E-04
R62=8.95E-04
R63=9.65E-04
R64=1.06E-03
R65=1.16E-03
R66=6.13E-04
R67=5.42E-04
R68=9.78E-04
R69=8.76E-04
R70=7.93E-04
R71=6.92E-04
R72=5.99E-04
R73=5.25E-04
R74=4.51E-04
R75=3.67E-04
R76=2.96E-04
R77=2.30E-04
R78=1.72E-04
R79=1.22E-04
R80=7.87E-05
R81=4.44E-05
R82=1.96E-05
R83=4.90E-06
R84=0.00E+00
!!!!!!!!!主梁实常数
R,10,2.2282,357.85,2.725,3.369,39,$RMORE,-1.72E-04,R10 !Composite Beam
R,11,31.704,3705.8,43.79,3.2,39, $RMORE,,143.79 !N5-5
R,12,64.314,12637,104.12,3.9385,45, $RMORE,,314.03 !N6-6
R,13,137.94,30960,464.02,5.7045,45, $RMORE,,1536 !N7-7
R,14,143.47,32144,427.81,5.49,45, $RMORE,,1437.2 !N8-8
R,15,108.78,23020,147.67,3.8,45, $RMORE,-2.10E-05,543.07 !N9-9
R,16,29.280,2522.1,39669,3.2,33.914,$RMORE,,126.49 !S5-5
R,17,143.47,32144,427.81,3.2,39, $RMORE,-1.59E-05,1437.2 !S8-8
R,18,64.314,12637,104.12,3.9385,45, $RMORE,,314.03 !S9-9
R,19,132.43,29364,395.76,5.3885,45, $RMORE,,1333.5 !S10-10
R,20,137.66,30475,366.72,5.2,45, $RMORE,-1.59E-05,1251.6 !S11-11
R,21,108.78,23020,147.67,3.8,45, $RMORE,-2.10E-05,543.09 !S12-12
!!!!!!!!!主缆
R,22,0.1264,R22 !SECTION0主缆
R,23,0.1264,R23 !SECTION1主缆
R,24,0.1264,R24 !SECTION2主缆
R,25,0.1264,R25 !SECTION3主缆
R,26,0.1264,R26 !SECTION4主缆
R,27,0.1264,R27 !SECTION5主缆
R,28,0.1264,R28 !SECTION6主缆
R,29,0.1264,R29 !SECTION7主缆
R,30,0.1264,R30 !SECTION8主缆
R,31,0.1264,R31 !SECTION9主缆
R,32,0.1264,R32 !SECTION10主缆
R,33,0.1264,R33 !SECTION11主缆
R,34,0.1264,R34 !SECTION12主缆
R,35,0.1264,R35 !SECTION13主缆
R,36,0.1264,R36 !SECTION14主缆
R,37,0.1264,R37 !SECTION15主缆
R,38,0.1264,R38 !SECTION16主缆
R,39,0.1264,R39 !SECTION17主缆
R,40,0.1264,R40 !SECTION18主缆
R,41,0.1264,R41 !SECTION19主缆
R,42,0.1264,R42 !SECTION20主缆
R,43,0.1264,R43 !SECTION21主缆
R,44,0.1264,R44 !SECTION22主缆
R,45,0.1264,R45 !SECTION23主缆
R,46,0.1264,R46 !SECTION24主缆
R,47,0.1264,R47 !SECTION25主缆
R,48,0.1264,R48 !SECTION26主缆
R,49,0.1264,R49 !SECTION27主缆
R,50,0.1264,R50 !SECTION28主缆
R,51,0.1264,R51 !SECTION29主缆
R,52,0.1264,R52 !SECTION30主缆
R,53,0.1264,R53 !SECTION31主缆
!!!!!!!!!吊索,,,,
R,54,0.0033549,R54 !吊索1
R,55,0.0033549,R55 !吊索2
R,56,0.0033549,R56 !吊索3
R,57,0.0033549,R57 !吊索4
R,58,0.0033549,R58 !吊索5
R,59,0.0033549,R59 !吊索6
R,60,0.0033549,R60 !吊索7
R,61,0.0033549,R61 !吊索8
R,62,0.0033549,R62 !吊索9
R,63,0.0033549,R63 !吊索10
R,64,0.0033549,R64 !吊索11
R,65,0.0033549,R65 !吊索12
R,66,0.0033549,R66 !吊索13
R,67,0.0033549,R67 !吊索14
R,68,0.0033549,R68 !吊索15
R,69,0.0033549,R69 !吊索16
R,70,0.0033549,R70 !吊索17
R,71,0.0033549,R71 !吊索18
R,72,0.0033549,R72 !吊索19
R,73,0.0033549,R73 !吊索20
R,74,0.0033549,R74 !吊索21
R,75,0.0033549,R75 !吊索22
R,76,0.0033549,R76 !吊索23
R,77,0.0033549,R77 !吊索24
R,78,0.0033549,R78 !吊索25
R,79,0.0033549,R79 !吊索26
R,80,0.0033549,R80 !吊索27
R,81,0.0033549,R81 !吊索28
R,82,0.0033549,R82 !吊索29
R,83,0.0033549,R83 !吊索30
R,84,0.0033549,R84 !吊索31
!!!主塔
R,85,1.46,6.28,6.28,5.47925,5.47925$RMORE,,7.00 !Tower标准断面
R,91,2.48E+01,1.28E+02,1.23E+02,6.6,6.6$RMORE,,2.08E+02 !Tower1-1Middle
R,92,3.38E+01,1.48E+02,1.46E+02,6.60E+00,6.60E+00,$RMORE,,2.50E+02, !Tower1-1UpAndDownZ轴沿顺桥向
R,93,4.32E+01,1.55E+02,1.54E+02,6.6,6.6$RMORE,,2.66E+02 !Tower2-2
R,101,0.03072,2.19E-03,7.20E-04,6.00E-01,6.00E-01,$RMORE,3.49E-06 !THL10-10Z轴沿翼缘
R,102,0.02272,1.52E-03,2.13E-04,4.00E-01,6.00E-01,$RMORE,2.42E-06 !THL11-11_400mmZ轴沿翼缘
R,112,0.02272,1.52E-03,2.13E-04,4.00E-01,6.00E-01,90$RMORE,2.42E-06 !THL11-11_400mmZ轴沿翼缘（旋转坐标轴）
R,103,0.02672,1.86E-03,4.17E-04,5.00E-01,6.00E-01,$RMORE,2.96E-06 !THL11-11_500mmZ轴沿翼缘
R,113,0.02672,1.86E-03,4.17E-04,5.00E-01,6.00E-01,90$RMORE,2.96E-06 !THL11-11_500mmZ轴沿翼缘（旋转坐标轴）
R,104,0.012816,7.20E-05,2.13E-04,3.00E-01,3.00E-01,90$RMORE,9.70E-07 !THL12-12Z轴沿腹板
R,105,0.03072,2.19E-03,7.20E-04,6.00E-01,6.00E-01,$RMORE,3.49E-06 !THL4-4Z轴沿翼缘
R,106,0.02272,1.52E-03,2.13E-04,4.00E-01,6.00E-01,$RMORE,2.42E-06 !THL5-5Z轴沿翼缘
R,116,0.02272,1.52E-03,2.13E-04,4.00E-01,6.00E-01,90$RMORE,2.42E-06 !THL5-5Z轴沿翼缘(旋转坐标系)
R,107,0.012816,7.20E-05,2.13E-04,3.00E-01,3.00E-01,90$RMORE,9.70E-07 !THL6-6Z轴沿腹板
R,108,0.2592,8.42E-02,8.42E-02,1.48E+00,1.48E+00,$RMORE,1.21E-01 !TowerColumn
R,109,0.01872,1.19E-03,9.01E-05,6.00E-01,3.00E-01,90$RMORE,1.89E-06 !BattenPlateZ轴沿腹板
!!!刚臂
R,86,1,1,1,1,1$RMORE,,1
!!!!!!!!!建立模型
!!!主梁节点（非锚跨）
*DO,I,1,32
N,I,(I-1)*7
*ENDDO
NGEN,2,1,32,,,2.5
ALLSEL
!!!主缆节点
N,,0,-20,7.74
N,,7,-20,7.86962963
N,,14,-20,8.258518519
N,,21,-20,8.906666667
N,,28,-20,9.814074074
N,,35,-20,10.98074074
N,,42,-20,12.40666667
N,,49,-20,14.09185185
N,,56,-20,16.0362963
N,,63,-20,18.24
N,,70,-20,20.70296296
N,,77,-20,23.42518519
N,,84,-20,26.40666667
N,,91,-20,29.64740741
N,,98,-20,33.14740741
N,,105,-20,36.90666667
N,,112,-20,40.92518519
N,,119,-20,45.20296296
N,,126,-20,49.74
N,,133,-20,44.34589967
N,,140,-20,39.24712944
N,,147,-20,34.44368931
N,,154,-20,29.93557929
N,,161,-20,25.72279937
N,,168,-20,21.80534956
N,,175,-20,18.18322984
N,,182,-20,14.85644024
N,,189,-20,11.82498073
N,,196,-20,9.088851331
N,,203,-20,6.648052034
N,,210,-20,4.50258284
N,,217,-20,2.65244375
N,,230,-20,0
ALLSEL
!!!刚臂节点
NSEL,S,,,1,32
NGEN,2,70,ALL,,,,-20
ALLSEL
!!!主梁单元（非锚跨）
TYPE,1
MAT,1
REAL,10
*DO,J,1,32
E,J,J+1
*ENDDO
!!!主梁附加质量单元
TYPE,3
REAL,1
E,33
REAL,2
*DO,K,1,32
E,K
*ENDDO
!!!刚臂单元建立
TYPE,1
MAT,5
REAL,86
*DO,K,1,32
E,K,K+70
*ENDDO
!!!吊索单元
TYPE,2
MAT,4
*DO,N,65,53,-1
REAL,(119-N)
E,N,N+37
*ENDDO
*DO,N,51,34,-1
REAL,(118-N)
E,N,N+37
*ENDDO
!!!主缆单元
TYPE,2
MAT,3
*DO,N,66,35,-1
REAL,(88-N)
E,N,N-1
*ENDDO
!!!索夹质量单元(索夹7忽略之)
TYPE,3
REAL,8
*DO,M,34,37
E,M
*ENDDO
REAL,7
*DO,M,38,40
E,M
*ENDDO
REAL,6
*DO,M,41,43
E,M
*ENDDO
REAL,5
*DO,M,44,46
E,M
*ENDDO
REAL,4
*DO,M,47,49
E,M
*ENDDO
REAL,3
*DO,M,50,51
E,M
*ENDDO
REAL,3
*DO,M,53,56
E,M
*ENDDO
REAL,4
*DO,M,57,59
E,M
*ENDDO
REAL,5
*DO,M,60,62
E,M
*ENDDO
REAL,6
*DO,M,63,65
E,M
*ENDDO
!!!主塔单元
NGEN,2,51,52,,,,,-8.17
NGEN,2,1,103,,,2,2
NGEN,2,2,103,,,2,-2
NGEN,2,3,103,,,-2,-2
NGEN,2,4,103,,,-2,2
NGEN,2,4,104,107,1,,,-1.2
NGEN,12,4,108,111,1,,,-3
NGEN,2,4,152,155,1,,,-4.5
TYPE,1
MAT,5
REAL,86
E,52,103
*DO,I,104,107
E,I,103
*ENDDO
TYPE,1
MAT,1
REAL,108
*DO,I,104,107
*DO,J,1,13
E,I+(J-1)*4,I+J*4
*ENDDO
*ENDDO
NGEN,2,50,110,111,,2,,-1.5
NGEN,2,52,108,,,-2,,-1.5
NGEN,2,1,160,,,2,-2
NGEN,2,1,161,,,-2,-2
NGEN,2,1,162,,,-2,2
NGEN,11,4,160,163,,,,-3
TYPE,1
MAT,1
REAL,109
*DO,I,1,10
*DO,J,1,3
E,112+(I-1)*4+J-1,112+(I-1)*4+J-1+48
E,112+(I-1)*4+J-1,112+(I-1)*4+J-1+49
E,112+(I-1)*4+J-1,112+(I-1)*4+J-1+52
E,112+(I-1)*4+J-1,112+(I-1)*4+J-1+53
*ENDDO
*ENDDO
E,108,160
E,108,161
E,109,161
E,109,162
E,110,162
E,110,163
E,111,160
E,111,163
E,152,200
E,152,201
E,153,201
E,153,202
E,154,202
E,154,203
E,155,203
E,155,200
*DO,I,115,151,4
E,I,I+45
E,I,I+48
E,I,I+49
E,I,I+52
*ENDDO
NGEN,2,50,158,,,2,2
ALLSEL
NGEN,2,1,208,,,,,-3
NGEN,2,1,209,,,,,-3.5
Real,93
mat,2
type,1
E,208,209
E,209,210
NGEN,2,1,210,,,,,-3.9
Real,92
mat,2
type,1
E,210,211
NGEN,2,1,211,,,,,-4.35
Real,91
mat,2
type,1
E,211,212
NGEN,2,1,212,,,,,-3.9
Real,92
mat,2
type,1
E,212,213
NGEN,2,1,213,,,,,-2.5
Real,93
mat,2
type,1
E,213,214
TYPE,1
MAT,5
REAL,86
*DO,I,156,159
E,208,I
*ENDDO
!!!塔横梁
NGEN,2,160,116,,,,1.282,1.887
NGEN,2,161,116,,,,1.3,11
NSEL,U,,,1,275
NGEN,2,2,276,,,,3.718,2.263
NGEN,2,1,278,,,,3.5,1.463
NGEN,2,1,279,,,,3.25,0.911
NGEN,2,1,280,,,,3.254,0.522
NGEN,2,1,281,,,,2.996,0.154
NGEN,2,6,277,,,,3.7
NGEN,2,1,283,,,,3.5
NGEN,3,1,284,,,,3.25
NGEN,2,1,286,,,,3
Real,101
mat,1
type,1
E,277,283
E,283,284
E,283,284
E,284,285
E,285,286
E,276,278
E,278,279
E,279,280
E,280,281
Real,103
mat,1
type,1
E,277,278
REAL,113
E,278,283
REAL,103
E,283,279
REAL,113
E,279,284
Real,102
mat,1
type,1
E,284,280
E,285,281
REAL,112
E,280,285
REAL,105
MAT,1
TYPE,1
E,286,287
E,281,282
REAL,106
MAT,1
TYPE,1
E,286,281
E,286,282
ALLSEL
NSEL,U,,,1,275
NGEN,2,15,ALL,,,-4
ALLSEL
ESEL,R,REAL,,101,120
ESEL,U,REAL,,108,109
ENGEN,25,2,15,ALL
ALLSEL
!CP,NEXT,ALL,161,152,153,154,155
CP,NEXT,ALL,103,277,292
CP,NEXT,ALL,116,276
CP,NEXT,ALL,119,291
!!!镜像单元
ALLSEL
NSYM,Y,350,ALL
ENSYM,500,,350,ALL
ALLSEL
CPSGEN,2,350,ALL
NUMMRG,NODE
!!!塔横梁
ESEL,S,REAL,,101,200
ESEL,U,REAL,,108,109
REAL,116
TYPE,1
MAT,1
E,287,282
E,302,297
REAL,104
TYPE,1
MAT,1
*DO,I,278,287
E,I,I+15
*ENDDO
*DO,I,628,631
E,I,I+15
*ENDDO
*DO,I,633,636
E,I,I+15
*ENDDO
NGEN,2,400,293,,,2,,3.425
E,278,693
E,283,693
E,298,693
E,293,693
NGEN,2,400,295,,,2,,2.238
E,280,695
E,285,695
E,300,695
E,295,695
NGEN,2,400,297,,,2,,1.9
E,282,697
E,287,697
E,302,697
E,297,697
NSYM,Y,350,693,695,2
NSEL,S,,,693,695,2
ESLN,S
ENSYM,20,,350,ALL
ALLSEL
!!!镜像单元
ALLSEL
NSYM,X,1500,ALL
ENSYM,1500,,1500,ALL
cpsgen,2,1500,all
!!!砼梁
K,,219.5 !钢混结合点
K,,224.35 !5#墩左
K,,233.5 !5#墩
K,,244.2 !5#墩右
K,,290.5 !6#墩
K,,338.5 !7#墩
ALLSEL
*DO,I,1,5
L,I,I+1
*ENDDO
ALLSEL
*DO,I,1,3
LSEL,S,,,I
LATT,2,(22-I),1
ALLSEL
*ENDDO
ALLSEL
*DO,I,4,5
LSEL,S,,,I
LATT,2,(21-I),1
ALLSEL
*ENDDO
LESIZE,ALL,7
LMESH,ALL
LSYMM,X,ALL,,,10,1
ALLSEL
*DO,I,6,8
LSEL,S,,,I
LATT,2,(21-I),1
ALLSEL
*ENDDO
LSEL,S,,,9,10
LATT,2,11,1
ALLSEL
LSEL,S,,,6,10
LMESH,ALL
!!!砼跨二期恒载质量单元
R,100,66857.14286
TYPE,3
REAL,100
*DO,I,2546,2585
E,I
*ENDDO
!!!主缆自锚
TYPE,1
MAT,5
REAL,86
E,2548,66
E,2548,416
E,2568,1566
E,2568,1916
!!!!!!!!!边界条件
ALLSEL
LPLOT
KSEL,S,,,15,16
KSEL,A,,,13
KSEL,A,,,3
KSEL,A,,,5,6
DK,ALL,UZ,,,,ROTX,UY
ALLSEL
CP,NEXT,UY,209,89
CP,NEXT,UZ,209,89
CP,NEXT,UY,559,439
CP,NEXT,UZ,559,439
CP,NEXT,UY,1709,1589
CP,NEXT,UZ,1709,1589
CP,NEXT,UY,2059,1939
CP,NEXT,UZ,2059,1939
!CP,NEXT,ALL,954,66,316
!CP,NEXT,ALL,974,566,816
ALLSEL
NWPAVE,214
CSYS,4
NSEL,R,LOC,Z,0
D,ALL,ALL
ALLSEL
WPCSYS,,0
NUMMRG,ALL
FINI 以下是求解部分的命令：
/SOLU
ACEL,,,9.8
antype,0
nsubst,auto!20,30,5,on
outres,all,all
sstif,on
nlgeom,on !几何非线性 大变形
nropt,full
ematwrite,yes!大变形(大位移)预应力模态分析,在变形后的位置上结构的振动特性
allsel
solve
finish
/solu
antype,modal
!upcoord,1,on !位移清零
pstres,on
modopt,lanb,100,,,,on !振型对单位1归一
mxpand,100
psolve,eiglanb
finish
/solu
expass,on
psolve,eigexp!扩展特征向量（振型）
我用杆系单元建立的模型，应该不存在单元畸变的问题，边界条件也检查过好几次了，尝试过不同的方法，但还是不行。希望大侠帮忙看下模型
页:
[1]