急问！关于弦支穹顶稳定性分析问题（附命令流）

boomsday

小弟我正在做一个弦支穹顶的稳定性分析，预计做双非线性分析，可是现在线性分析就不对劲，主要的问题有：一、得到的特征值屈曲系数太小，没有按实际荷载进行工况组合，而是每个节点加单位力1N，特征值屈曲系数只有几十；二、环索的预应力损失很大，我用的是降温法加预应力，比如最外圈预计要加3000KN，计算结果显示只有2000KN；三、支座反力与单层网壳相比，竖向力不变，水平力反而增加了。
     我估计原因可能有以下几方面：一、我是按成形之后的结构建模的，只给环索加上了预应力，而此时其实网壳的杆件里已经有应力了，是否事先要进行找形计算？；二、索预应力损失大，弦支穹顶的变形大；三、ansys中特征值稳定性分析的一般做法是先加单位力静力分析，再进行屈曲分析，而我在静力分析中还加入了重力和索的预应力，这个是否有影响？
     小弟我做了半个多月了，还是没有什么进展，非常着急，希望论坛上的ggjj们帮我看看，提点意见啊。本人的邮箱是boomsday@emails.bjut.edu.cn，qq号是287719205，欢迎交流。
     命令流如下：
   /filname,suspend-dome
/prep7
csys,2
*do,i,1,56                    
k,i,138,(i-1)*360/56,70.5
*enddo
*do,i,57,112
k,i,138,(i-56)*360/56-180/56,72
*enddo
*do,i,113,168
k,i,138,(i-113)*360/56,73.5
*enddo
*do,i,169,224
k,i,138,(i-168)*360/56-180/56,75
*enddo
*do,i,225,280
k,i,138,(i-225)*360/56,76.5
*enddo
*do,i,281,308
k,i,138,(i-280)*360/28-360/56,78
*enddo
*do,i,309,336
k,i,138,(i-309)*360/28,79.5
*enddo
*do,i,337,364
k,i,138,(i-336)*360/28-360/56,81
*enddo
*do,i,365,392
k,i,138,(i-365)*360/28,82.5
*enddo
*do,i,393,406
k,i,138,(i-392)*360/14-360/28,84
*enddo
*do,i,407,420
k,i,138,(i-407)*360/14,85.5
*enddo
*do,i,421,434
k,i,138,(i-421)*360/14-360/28,87
*enddo
*do,i,435,441
k,i,138,(i-435)*360/7+360/28,88.5
*enddo
k,442,138,0,90
*do,i,1,5
  *do,j,1,55,1
   l,56*(i-1)+j,56*(i-1)+j+1
  *enddo
l,i*56,i*56-55
*enddo

*do,i,1,4
  *do,j,1,27,1
   l,280+28*(i-1)+j,280+28*(i-1)+j+1
  *enddo
l,280+i*28,280+i*28-27
*enddo

*do,i,1,3
  *do,j,1,13,1
   l,392+14*(i-1)+j,392+14*(i-1)+j+1
  *enddo
l,392+i*14,392+i*14-13
*enddo

*do,j,1,6,1
l,434+j,434+j+1
*enddo
l,441,435
  
!建立径杆单元
*do,i,1,4
  *do,j,1,56,1
   l,56*(i-1)+j,56*i+j
   *enddo
*enddo
*do,j,2,56,1
  l,j,55+j
  
  l,55+j,112+j
  
  l,112+j,167+j
  
  l,167+j,224+j
  
*enddo
l,1,112
l,112,113
l,113,224
l,224,225
*do,j,1,28,1
  l,223+2*j,280+j
  l,224+2*j,280+j
*enddo
*do,j,1,27,1
  l,225+2*j,280+j
  
*enddo
l,225,308
!建立第二道三层径干
*do,i,1,3
  *do,j,1,28,1
   l,280+28*(i-1)+j,280+28*i+j
   *enddo
*enddo
*do,j,1,27,1
  l,280+j,309+j
  l,309+j,336+j
  l,336+j,365+j
*enddo
l,308,309
l,309,364
l,364,365
!建立第三道径杆
*do,j,1,14,1
  l,363+2*j,392+j
l,364+2*j,392+j
*enddo
*do,j,1,13,1
  l,365+2*j,392+j
*enddo
l,365,406
!谢干
*do,i,1,2
  *do,j,1,14,1
   l,392+14*(i-1)+j,392+14*i+j
   *enddo
*enddo
*do,j,1,13,1
  l,392+j,407+j
*enddo
l,406,407
*do,j,1,13,1
  l,406+j,421+j
*enddo
l,420,421
!最后一圈
*do,j,1,7,1
  l,419+2*j,434+j
  l,420+2*j,434+j
*enddo
*do,j,1,6,1
  l,421+2*j,434+j
*enddo
l,421,441
!屋顶节点杆
*do,i,435,441,1
  l,i,442
*enddo

csys,0    !复制节点，要转换到卡式坐标系
kgen,2,407,420,,,,-3.9,36  !从网壳顶起第一道索
!环向索
*do,i,1,13,1
  l,442+i,443+i
*enddo
l,456,443
kgen,2,366,392,0,,,-3.9,91   !从网壳顶起第二道索
!环向索
*do,i,1,13,1
  l,455+2*i,457+2*i
*enddo
l,483,457
kgen,2,309,336,1,,,-3.9,175   !从网壳顶起第三道索
*do,i,1,27,1                  !环向索
  l,483+i,484+i
*enddo
l,511,484              !倒511节点
kgen,2,226,280,1,,,-3.9,286   !从网壳顶起第四道索
*do,i,1,27,1                  !环向索
  l,510+2*i,512+2*i
*enddo
l,566,512              !倒566节点
kgen,2,113,167,1,,,-3.9,454   !从网壳顶起第五道索
*do,i,1,27,1                  !环向索
  l,565+2*i,567+2*i
*enddo
l,621,567              !倒621节点

!径向索
*do,i,1,14,1            !径向索
  l,442+i,364+2*i
*enddo
*do,i,1,13,1      
  l,443+i,364+2*i
*enddo
l,443,392
*do,i,1,14,1           !径向索
  l,455+2*i,307+2*i
*enddo
*do,i,1,13,1      
  l,455+2*i,309+2*i
*enddo
l,483,309            
*do,i,1,28,1                  !径向索
   l,483+i,224+2*i
*enddo
*do,i,1,27,1
   l,484+i,224+2*i
*enddo
l,484,280
*do,i,1,28,1                  !径向索
   l,510+2*i,111+2*i
*enddo
*do,i,1,27,1
   l,510+2*i,113+2*i
*enddo
l,566,113
*do,i,1,28,1                  !径向索
   l,565+2*i,2*i
*enddo
*do,i,1,27,1
   l,567+2*i,2*i
*enddo
l,567,56

*do,i,1,14,1     !竖向撑杆
  l,442+i,406+i
*enddo
*do,i,1,14,1     !竖向撑杆
  l,455+2*i,364+2*i
*enddo
*do,i,1,28,1     !竖向撑杆
  l,483+i,308+i
*enddo
*do,i,1,28,1     !竖向撑杆
  l,510+2*i,224+2*i
*enddo
*do,i,1,28,1     !竖向撑杆
  l,565+2*i,111+2*i
*enddo
!划分网壳单元
!环向杆材料模型定义
et,1,beam188
mp,ex,1,2.1e11    !等
mp,prxy,1,0.3
mp,gxy,1,8e10  !定义剪切模量
mp,dens,1,7.85e+3
SECTYPE,1,BEAM,CTUBE,huanx1, 0
SECOFFSET, CENT
SECDATA,0.1085,0.1225,0,0,0,0,0,0,0,0
SECTYPE,2,BEAM,CTUBE,huanx2, 0
SECOFFSET, CENT
SECDATA,0.1105,0.1225,0,0,0,0,0,0,0,0
SECTYPE,3,BEAM,CTUBE,huanx3, 0
SECOFFSET,CENT
SECDATA,0.1135,0.1225,0,0,0,0,0,0,0,0  

!===============================

lsel,s,line,,1,280,1  !划分外5圈网壳单元
lesize,all,,,1
LATT,1,,1, , , ,1
lmesh,1,280,1
lsel,s,line,,281,392,1  !划分外6-8圈网壳单元
lesize,all,,,1
LATT,1,,1, , , ,2
lmesh,281,392,1
lsel,s,line,,393,441,1  !划分外9-12圈网壳单元
lesize,all,,,1
LATT,1,,1, , , ,3
lmesh,393,441,1

!=====================================================
!划分径向杆单元

mp,ex,2,2.1e11
mp,prxy,2,0.3
mp,gxy,2,8e10
mp,dens,2,7.85e+3
et,2,beam188
SECTYPE,4,BEAM,CTUBE,jingx1, 0
SECOFFSET,CENT
SECDATA,0.1105,0.1225,0,0,0,0,0,0,0,0
SECTYPE,5,BEAM,CTUBE,jingx2, 0
SECOFFSET,CENT
SECDATA,0.1135,0.1225,0,0,0,0,0,0,0,0   
!===================================
lsel,s,line,,442,1141,1  !划分外8圈径向杆单元
lesize,all,,,1
LATT,2,,2, , , ,4
lmesh,442,1141,1

!划分径向杆单元
lsel,s,line,,1142,1267,1  !划分9-13圈径向杆单元
lesize,all,,,1
LATT,2,,2, , , ,5
lmesh,1142,1267,1
!=====================================================

!划分环向锁单元
mp,ex,3,1.6e11 !定义link10单元的弹性模量
mp,prxy,3,0.3 !定义link10单元的泊松系数
mp,alpx,3,1.0e-5!定义线膨胀系数
mp,dens,3,2.78e+3
et,3,link10
r,1,7.6e-3 !定义link10单元的面积
lsel,s,line,,1268,1379,1  !划分环锁杆单元
lesize,all,,,1
LATT,3,1,3, , , ,
lmesh,1268,1379,1
mp,ex,4,2.1e11 !定义link10单元的弹性模量
mp,prxy,4,0.3 !定义link10单元的泊松系数
mp,gxy,5,8e10
mp,dens,4,7.85e+3
et,4,link8
r,2,7.6e-3 !定义link8单元的面积

!划分径向锁
lsel,s,line,,1380,1603,1  
lesize,all,,,1
LATT,4,2,4, , , ,
lmesh,1380,1603,1
!=================================================
!开始讨论撑杆的建立
mp,ex,5,2.1e11
mp,prxy,5,0.3
mp,gxy,5,8e10
mp,dens,5,2.78e+3
et,5,link8
r,3,3.1e-3         
!划分撑杆单元
lsel,s,line,,1604,1715,1  !划分撑杆单元
lesize,all,,,1
LATT,5,3,5, , , ,
lmesh,1604,1715,1
allsel

/solu
antype,0   !静力分析
pstres,on
acel,0,0,9.8
*do,i,2,56,2  !施加边界条件
  dk,i,ux
  dk,i,uy
  dk,i,uz
*enddo
egjx=2e5                       !Ey
agjx=7600  !单根钢绞线面积       !施加负预应力
ehnt=4e4 !Eh
xzxs=1.0e-5                     !线胀系数
yjl1=3800000 !定义预加力
yjl2=170000
yjl3=950000
yjl4=650000
yjl5=280000
*do,i,1268,1281,1   
  bfl,i,temp,-yjl5/(xzxs*egjx*agjx) !对5钢绞线施加温度
*enddo
*do,i,1282,1295,1   
  bfl,i,temp,-yjl4/(xzxs*egjx*agjx) !对4钢绞线施加温度
*enddo
*do,i,1296,1323,1   
  bfl,i,temp,-yjl3/(xzxs*egjx*agjx) !对3钢绞线施加温度
*enddo
*do,i,1324,1351,1   
  bfl,i,temp,-yjl2/(xzxs*egjx*agjx) !对2钢绞线施加温度
*enddo
*do,i,1352,1379,1   
  bfl,i,temp,-yjl1/(xzxs*egjx*agjx) !对1钢绞线施加温度
*enddo
*do,i,1,442,1    !施加外力得到的frequency实际上是临界载荷系数,也就是说
   fk,i,fz,-1
*enddo
solve
finish

/solu
antype,1
bucopt,lanb,6,0,0,0
mxpand,6,0,0,1,0.001,
solve