ANSYS 提取一个简易斜拉桥数值模型刚度矩阵的三种方法

/clear
/filname,cable-stayedbridge
/units,si
/vup,,z
!abbres,,GJ_1_v3.0,abbr
!对称结构取半
/prep7
!定义单元类型
et,1,beam4
et,2,link8
!link10为单向拉压，会造成影响矩阵的非线性，用link8
!定义材料属性
!tower
mp,ex,1,2.05e11
mp,prxy,1,0.25
mp,dens,1,2600
!girder
mp,ex,2,2.05e11
mp,prxy,2,0.25
mp,dens,2,2600
E_cable=2.05e11
A_cable=0.003
alpha_cable=1.2e-5
!cable-stay
mp,ex,3,E_cable
mp,prxy,3,0.1
mp,dens,3,2600
mp,alpx,3,alpha_cable
!定义实常数
r,1,2,10,0.7,1,1 !梁单元（AREA,IZZ,IYY,TKZ,TKY）
r,2,4,10,5,1,1 !索塔
r,5,A_cable !拉索（AREA,ISTRN）
/replot
h1=45 !桥塔在桥面以上高度h1
h0=5 !桥塔上拉索间距h1
*do,i,1,13
n,100+i,0,0,h1-i*h0
*enddo
l1=-60
l0=5
*do,i,1,24
n,300+i,l1+i*l0,0,0
*enddo
type,2
mat,3
real,5
i=1
en,11,101,305
en,21,101,319
type,1
mat,1
real,2
*do,i,1,12
en,100+i,100+i,100+i+1
*enddo
type,1
mat,2
real,1
*do,i,1,23
en,300+i,300+i,300+i+1
*enddo
!d,113,all
!d,301,uz,,,,,uy
!!d,305,uz,,,,,uy
!d,324,uy,,,,,rotx,roty,rotz
allsel,all
cpintf,all,0.1
```
- 方法一
```matlab
/aux2
fileaux2,cable-stayedbridge,emat
form,long
dump,all
```
-  方法二
```matlab
/output,cp,out
/debug,-1,,,1
/solu
solve
finish
/output,term
```
- 方法三
```matlab
/solu
antype,7 !substructuring分析类型
seopt,cable-stayedbridge,2
!设置文件名称和刚度矩阵类型（刚度，质量，阻尼等）
! 1为刚度矩阵，2为刚度矩阵&质量矩阵，3为刚度矩阵&质量矩阵&阻尼矩阵
nsel,all!选择所有节点
nsel,r,,,312
m,all,all!定义所有节点自由度为主自由度
solve!求解
selist,cable-stayedbridge,3!列出整体刚度矩阵
```