一个斜拉桥命令流（天津永和斜拉桥）

我参照天津永和斜拉桥的博士论文，建了模型，出现了众多问题。
1）刚开始出现错误，我害怕是beam4单元方向不对衬，如在x向上建模时，同一条梁我有的从左往右，有的从右往左，而且使
用了转换坐标系（因为需要镜称，从0，0，刚开始我建错了，最后各node都减去250m），但这些我都改过来了，害得我许多本
来以*do实现的，现在手动建单元，结果发现几乎没关系。
2）如果直接进行模态分析，不开预应力开关，跟实测有点差距，于是我更改质量（我害怕原文没读懂，密度设错了），较为
接近，但也有30%的差距，而且振型与实测不同。
3）直接模态分析，开预应力开关，ANS出现错误，估计是没计算出预应力文件，显示缺少一个什么文件大概就是那个意思
4）先静力分析，开预应力，成桥，结果发现不收敛
5）先静力分析，开预应力，适当调整刚度，或者调整索的初始应变（0.1倍），开大变型，但是接下来模态不能分析。因为模
态是线形的
6）将预应力调为极小（0.01甚至0.001倍），不开大变型，静力分析能成桥，结果发现实测频率前三阶与实测差几乎几个数量
级

请高手指点我的错误在哪（包括我定义的质量、刚度的单位、数量级），以及上述我的分析对吗，有预应力的桥算模态，必须
先静力分析计算预应力吗？
!改顺序，大改
finish
/clear
/TITLE,cable-stayed bridge
/prep7
!!!!!!!!!!!!!!!!!!!!!!!单元类型和属性
hgd=1
sgd=1
bgd=1e-10
xs=1
mm=1
zlm=1
suotanmo=1
jl=255
et,1,beam4
et,2,link10
!!!!!!!!!!!!!!!!!!!!!!!!!实常数
!索
!右侧 长 (并且从柱近到远)
r,1,62.05e-4,0.00250685*xs,
r,2,54.19e-4,0.0022541*xs,
r,3,62.05e-4,0.00240935*xs,
r,4,66.76e-4,0.00248915*xs,
r,5,76.19e-4,0.00252195*xs,
r,6,84.04e-4,0.00251485*xs,
r,7,95.03e-4,0.0025329*xs,
r,8,109.17e-4,0.00252995*xs,
r,9,120.17e-4,0.00257425*xs,
r,10,110.74e-4,0.0025113*xs,
r,11,170.43e-4,0.00263655*xs,
!左侧，短
r,12,62.05e-4,0.00228605*xs,
r,13,54.19e-4,0.0023187*xs,
r,14,62.05e-4,0.00241175*xs,
r,15,66.76e-4,0.0026116*xs,
r,16,76.19e-4,0.00252195*xs,
r,17,84.04e-4,0.00264935*xs,
r,18,84.04e-4,0.0024982*xs,
r,19,95.03e-4,0.00236345*xs,
r,20,156.29e-4,0.0023069*xs,
r,21,140.59e-4,0.0024902*xs,
r,22,170.43e-4,0.00231325*xs,
!乘弹模差不多,质量也是，搞定!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!主梁
R,23,8.496,141.914,4.425,2,13.6 , ,  
RMORE, ,6.3678, , , , ,
R,24,21.683,252.8265,7.0283,2,13.6, ,  
RMORE, ,11.7888, , , , ,

R,25,6.6385,152.2056,1.9611,2,13.6, ,  
RMORE, ,2.3932, , , , ,

!塔
R,26,7.2,6.264,3.9,3,3, ,   
RMORE, ,2.5527, , , , ,
R,27,5.76,5.8752,5.8752,3,3, ,  
RMORE, ,1.870, , , , ,
!  不太确定？？？？？？？？？？？？？？？？
R,28,3.5,1.5417,2.3229,3,3, ,  
RMORE, ,0.7832, , , , ,

!边梁

r,29,1,0.08333333,0.08333333,1,1,,
RMORE, ,0.3333, , , , ,
!鱼刺梁
r,30,1,0.08333333,0.08333333,1,1,,
RMORE, ,0.3333, , , , ,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!材料属性
!主梁三段
!l1
mp,ex,1,3.5e10*hgd
mp,prxy,1,0.1667
mp,dens,1,2727*zlm*0.312
!l2
mp,ex,2,3.5e10*hgd
mp,prxy,2,0.1667
mp,dens,2,2649*zlm*0.797
!l3
mp,ex,3,3.5e10*hgd
mp,prxy,3,0.1667
mp,dens,3,2762*zlm*0.24
!塔
mp,ex,4,3.5e10*sgd
mp,prxy,4,0.1667
mp,dens,4,2600
!鱼刺梁
mp,ex,5,10e15
mp,prxy,5,0
mp,dens,5,0
!边梁
mp,ex,6,3.5e10*bgd
mp,prxy,6,0
mp,dens,6,2500*mm
!索
mp,ex,7,2e11*suotanmo
mp,prxy,7,0.25
mp,dens,7,1200
!?????????????????????????
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!开始建模
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!节点
n,1,0-jl,,
n,2,21.03-jl,,
n,3,24.09-jl,,
n,4,26.9-jl,,
n,5,36.75-jl,,
n,6,48.35-jl,,
n,7,59.95-jl,,
n,8,71.53-jl,,
n,9,83.75-jl,,
n,10,94.75-jl,,
n,11,106.35-jl,,
n,12,117.95-jl,,
n,13,125-jl,,
n,14,132.05-jl,,
n,15,143.65-jl,,
n,16,155.25-jl,,
n,17,166.85-jl,,
n,18,178.45-jl,,
n,19,190.05-jl,,
n,20,201.65-jl,,
n,21,213.25-jl,,
n,22,224.85-jl,,
n,23,236.45-jl,,
n,24,248.05-jl,,
n,25,255-jl,,
n,26,19.35-jl,
n,27,25.15-jl
n,28,30.95-jl
i=28
a=-7.25
n,i+1,0-jl,a,
n,i+2,21.03-jl,a,
n,i+3,24.09-jl,a,
n,i+4,26.9-jl,a,
n,i+5,36.75-jl,a,
n,i+6,48.35-jl,a,
n,i+7,59.95-jl,a,
n,i+8,71.53-jl,a,
n,i+9,83.75-jl,a,
n,i+10,94.75-jl,a,
n,i+11,106.35-jl,a,
n,i+12,117.95-jl,a,
n,i+13,125-jl,a,
n,i+14,132.05-jl,a,
n,i+15,143.65-jl,a,
n,i+16,155.25-jl,a,
n,i+17,166.85-jl,a,
n,i+18,178.45-jl,a,
n,i+19,190.05-jl,a,
n,i+20,201.65-jl,a,
n,i+21,213.25-jl,a,
n,i+22,224.85-jl,a,
n,i+23,236.45-jl,a,
n,i+24,248.05-jl,a,
n,i+25,255-jl,a,
n,i+26,19.35-jl,a
n,i+27,25.15-jl,a
n,i+28,30.95-jl,a
i=56
a=7.25
n,i+1,0-jl,a,
n,i+2,21.03-jl,a,
n,i+3,24.09-jl,a,
n,i+4,26.9-jl,a,
n,i+5,36.75-jl,a,
n,i+6,48.35-jl,a,
n,i+7,59.95-jl,a,
n,i+8,71.53-jl,a,
n,i+9,83.75-jl,a,
n,i+10,94.75-jl,a,
n,i+11,106.35-jl,a,
n,i+12,117.95-jl,a,
n,i+13,125-jl,a,
n,i+14,132.05-jl,a,
n,i+15,143.65-jl,a,
n,i+16,155.25-jl,a,
n,i+17,166.85-jl,a,
n,i+18,178.45-jl,a,
n,i+19,190.05-jl,a,
n,i+20,201.65-jl,a,
n,i+21,213.25-jl,a,
n,i+22,224.85-jl,a,
n,i+23,236.45-jl,a,
n,i+24,248.05-jl,a,
n,i+25,255-jl,a,
n,i+26,19.35-jl,a
n,i+27,25.15-jl,a
n,i+28,30.95-jl,a

!84ge
!塔
*do,i,0,10
n,85+i,125-jl,-1*a,36+i*1.575
n,96+i,125-jl,a,36+i*1.575
*end do

!106ge
*do,j,1,20
n,106+j,125-jl,-1*a-0.1435*j,36-j*2.075
n,128+j,125-jl,a+0.1435*j,36-j*2.075
*end do
/view,1,1,1,1
/ANGLE,1,270,XM,0   
/REPLOT
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建单元
!索单元
type,2
mat,7
*do,k,1,11
real,23-k
e,29+k,96-k
e,57+k,107-k
real,12-k
e,53-k,96-k
e,81-k,107-k
*end do
!主梁单元0-19.35
type,1
mat,1
real,23
e,1,26

!主梁单元19.35-30.95
type,1
mat,2
real,24
e,26,2
e,2,3
e,3,27
e,27,4
e,4,28

!主梁单元30.95-255
type,1
mat,3
real,25
e,28,5
*do,i,5,24
e,i,i+1
*end do
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!边梁单元0-19.35
type,1
mat,6
real,29

e,29,54
e,57,82
!边梁单元19.35-30.95

p=28
e,26+p,2+p
e,2+p,3+p
e,3+p,27+p
e,27+p,4+p
e,4+p,28+p
p=56
e,26+p,2+p
e,2+p,3+p
e,3+p,27+p
e,27+p,4+p
e,4+p,28+p
!边梁单元30.95-255
e,56,33
e,84,61
*do,i,5,24
e,i+28,i+29
e,i+56,i+57
*end do


!鱼刺梁
type,1
mat,5
real,30

*do,i,1,28
e,i,i+28
e,i,i+56
*end do

!塔
!索区
type,1
mat,4
real,26
e,95,94
e,94,93
e,93,92
e,92,91
e,91,90
e,90,89
e,89,88
e,88,87
e,87,86
e,86,85
e,106,105
e,105,104
e,104,103
e,103,102
e,102,101
e,101,100
e,100,99
e,99,98
e,98,97
e,97,96
!非索区
type,1
mat,4
real,27

e,85,107
e,96,129
*do,i,107,125
e,i,i+1
e,i+22,i+23
*end do


!坐标系
nsym,x,200,all
esym,x,200,all
nummrg,all

!塔横梁
type,1
mat,4
real,28
n,1000,-130,,36
n,1001,-130,,51.75
e,1000,85
e,1000,96
e,1001,95
e,1001,106

n,2000,130,,36
n,2001,130,,51.75
e,2000,285
e,2000,296
e,2001,295
e,2001,306

!释放
cpintf,uy
cpintf,uz
cpintf,rotx
cpintf,rotz
! 给桥台加约束
allsel
d,1,uy,uz

d,27,uy,uz

d,201,uy,uz

d,227,uy,uz

nsel,,loc,z,-5.5
d,all,all
allsel
NUMCMP,all

！下面分别为几种求解方式
!求解
/solu
pstres,on
solve
finish
save

/solu
antype,modal
modopt,lanb,20
mxpand,10,0,100
sstif,on
solve
finish



/solu
antype,modal
modopt,lanb,20
mxpand,10,0,100
lumpm,on
pstres,on
upcoord,1.0,on
solve
finish
/POST1  
SET,LIST


!求解
/solu
NLGEOM,1
lumpm,on
pstres,on
solve
finish
save

有没有用实体建模的啊？？？

回复 icewangqing 的帖子

找了很久，终于找到了，非常好！在此谢过了！

看看 正在学习中

很好的学习资料！

very good thanks

这是哪篇博士论文啊