一个简单的斜拉桥模型 命令流

!最近在学习ansys ，发一个简单的斜拉桥模型命令流吧
!20130312_ansys斜拉桥最小弯曲能量法
!该方法要将索 塔 梁 的轴向刚度增至足够大（通常将梁 塔 索的面积增大10e4倍，同时为了保证梁塔索的自重不变
!还要将其密度降低10e4倍）进行一次成桥计算即可，这样得到的成桥主梁和塔中弯矩总体上是很均匀且很小的索力总体上也上合理!的，但可能在尾索出现不合理的现象，这里采用一种简单方法，通过调整边墩的刚度和在边墩处主梁施加压重来处理
!另外，这是个简单的例子，许多方面都没有考虑，
!如果有什么错误的话，请大家指出，我也是新手，也有许多不懂的地方
finish  
/CLE
/PREP7
!定义主梁
ET,1,BEAM3
MP,EX,1,3.65E10   !E = 3.65E10 PA
MP,DENS,1,0.29441   !DENS=2944.1 KG/M3 缩小10000倍
R,1,0.15365E6,0.98662E1   !面积增大10000倍
!定义塔
ET,2,BEAM3
MP,EX,2,3.5E10
MP,DENS,2,0.27 !DENS=2700 KG/M3 缩小10000倍
R,2,0.20540E6, 0.86365E2
R,3,0.25840E6, 0.10584E3
R,4,0.44424E6, 0.25260E3

!定义索
ET,3,BEAM3
MP,EX,3,2.0E11
MP,DENS,3,0.8 !DENS=8000 KG/M3 缩小10000倍
R,7,0.83896E2, 0   !抗弯惯性矩=0
R,8,0.97751E2, 0
R,9,0.11622E3, 0
R,10,0.12546E3,0
R,11,0.14393E3,0
R,12,0.16240E3,0
R,13,0.24091E3,0
!定义墩
ET,4,BEAM3
MP,EX,4,3.0E10
MP,DENS,4,2500  !DENS=2500 KG/M3 不用缩小10000倍
MP,EX,5,3.0E10
MP,DENS,5,0
R,5,0.28275E2, 0.15905E2
R,6,0.35343E2, 0.19880E2
R,14,0.16E1,0.15E-4  !抗弯惯性矩=0取很小可认为是铰接
!定义节点
!主塔部分
N,1,-5.2,6.0
N,2,0   ,6.0
N,3, 5.2,6.0
*DO,I,0,2,1
N,I+4,0,26.00+7*I
*ENDDO
N,7,0,47.310
N,8,0,49.680
N,9,0,78.310
N,10,0,101.110
N,101,-1.4,101.110  !生成索在塔上的锚固点
N,102, 1.4,101.110  !生成索在塔上的锚固点
N,11,0,103.610
N,111,-1.41,103.610  !生成索在梁上的锚固点
N,112, 1.41,103.610  !生成索在梁上的锚固点
N,12,0,105.810
N,121,-1.425,105.810  !生成索在梁上的锚固点
N,122, 1.425,105.810  !生成索在梁上的锚固点
N,13,0,107.610
N,131,-1.445,107.610  !生成索在梁上的锚固点
N,132, 1.445,107.610  !生成索在梁上的锚固点
*DO,I,0,9,1
N,14+I,0,109.110+1.2*I
*ENDDO
N,141,-1.460,109.110  !生成索在梁上的锚固点
N,142, 1.460,109.110  !生成索在梁上的锚固点
N,151,-1.480,110.310  !生成索在梁上的锚固点
N,152, 1.480,110.310  !生成索在梁上的锚固点
N,161,-1.495,111.510  !生成索在梁上的锚固点
N,162, 1.495,111.510  !生成索在梁上的锚固点
N,171,-1.510,112.710  !生成索在梁上的锚固点
N,172, 1.510,112.710  !生成索在梁上的锚固点
N,181,-1.525,113.910  !生成索在梁上的锚固点
N,182, 1.525,113.910  !生成索在梁上的锚固点
N,191,-1.520,115.110  !生成索在梁上的锚固点
N,192, 1.520,115.110  !生成索在梁上的锚固点
N,201,-1.535,116.310  !生成索在梁上的锚固点
N,202, 1.535,116.310  !生成索在梁上的锚固点
N,211,-1.545,117.510  !生成索在梁上的锚固点
N,212, 1.545,117.510  !生成索在梁上的锚固点
N,221,-1.555,118.710  !生成索在梁上的锚固点
N,222, 1.555,118.710  !生成索在梁上的锚固点
N,231,-1.565,119.910  !生成索在梁上的锚固点
N,232, 1.565,119.910  !生成索在梁上的锚固点

N,24,0,121.410
N,241,-1.575,121.410  !生成索在梁上的锚固点
N,242, 1.575,121.410  !生成索在梁上的锚固点
N,25,0,122.910
N,251,-1.580,122.910  !生成索在梁上的锚固点
N,252, 1.580,122.910  !生成索在梁上的锚固点
!主梁部分
*DO,I,0,15,1
!生成主梁上的结点
N,26+2*I,-13.2-8*I,49.680
N,26+2*I+1,13.2+8*I,49.680
!生成索在梁上的锚固点
N,300+26+2*I,-13.2-8*I+0.5,49.680-1.297
N,300+26+2*I+1,13.2+8*I-0.5,49.680-1.297
*ENDDO
!墩
N,58,-137.0,49.680
N,59,137.0,49.680
N,60,-137.0,6.0
N,61,137.0,6.0
!形成单元
TYPE,3
MAT,4
REAL,5
E,1,4
REAL,6
E,2,4
REAL,5
E,3,4
TYPE,2
MAT,2
REAL,4
*DO,I,0,3,1
E,4+I,5+I
*ENDDO
REAL,3
E,8,9
E,9,10
REAL,2
*DO,I,0,14,1
E,10+I,11+I
*ENDDO
TYPE,1
MAT,1
REAL,1
E,26,8
E,8,27
K=10
V=26
*DO,I,0,15,1
TYPE,2
MAT,3
*IF,I,LT,3,THEN
  REAL,7
*ENDIF
*IF,I,GE,3,AND,I,LT,6,THEN
  REAL,8
*ENDIF
*IF,I,GE,6,AND,I,LT,9,THEN
  REAL,9
*ENDIF
*IF,I,GE,9,AND,I,LT,12,THEN
  REAL,10
*ENDIF
*IF,I,GE,12,AND,I,LT,14,THEN
  REAL,11
*ENDIF
*IF,I,EQ,14,THEN
  REAL,12
*ENDIF
*IF,I,EQ,15,THEN
  REAL,13
*ENDIF
E,300+V,10*K+1   !形成索
E,10*K+2,300+V+1
CERIG,K,10*K+1   !将索的锚定点与主塔中心刚接
CERIG,K,10*K+2
CERIG,V,300+V    !将索的锚固点与主梁中心刚接
CERIG,V+1,300+V+1
K=K+1
TYPE,1
REAL,1
MAT,1
E,V,V+2
E,V+1,V+3

V=V+2
*ENDDO
TYPE,3
REAL,14
MAT,5
E,60,58
E,61,59  
FINISH
!加荷载和边界
/SOLU
D,1,ALL,,,3
D,60,ALL
D,61,ALL
F,58,FY,-1000000 !边墩处主梁施加压重来防止索力不合理
F,59,FY,-1000000
ACEL,,10  !加自重
SOLVE
FINISH
/POST1
PLDISP,1
!定义单元表
ETABLE,UX,U,X
ETABLE,UY,U,Y
ETABLE,ROTZ,ROT,Z
ETABLE,FX1,SMISC,1  !I端轴力
ETABLE,FX2,SMISC,7
ETABLE,FY1,SMISC,2   !I端剪力
ETABLE,FY2,SMISC,8
ETABLE,M1,SMISC,6    !I端弯矩
ETABLE,M2,SMISC,12
ETABLE,SMAX1,NMISC,1  !I端最大应力
ETABLE,SMAX2,NMISC,3  !J端最大应力
ETABLE,SMMIN1,NMISC,2  !I端最小应力
ETABLE,SMIN2,NMISC,4  !J端最小应力
PRETAB,UY   !列表显示UY