matlab刚度矩阵求逆,ANSYS求结构整体刚度矩阵逆矩阵APDL命令流

虽然采用matlab等相关软件对于大型矩阵的求解非常方便，但是这里涉及到ANSYS与matlab接口的处理问题，比较繁琐，本文编制一个直接在ANSYS中提取结构整体刚度矩阵和对整体刚度求逆的APDL程序，能在ANSYS中直接操作刚度矩阵，比较实用。

利用ANSYS求结构刚度矩阵逆矩阵，模型采用单层球面网壳！

相关的APDL命令流如下所示。

finish

/clear

/filname,modal

/prep7

et,1,beam188

mp,ex,1,2.1e11

mp,prxy,1,0.28

mp,dens,1,7850

*afun,deg

!用户界面设计，输入几何参数

multipro,'start',4

*cset,1,3,f,'Enter the f',6.8

*cset,4,6,span,'Enter the span',30

*cset,7,9,Kn,'Enter the radial number',24

*cset,10,12,Nx,'Enter the node circle number',3

*cset,61,62,'please enter the geometry parameters'

multipro,'end'

!计算节点坐标位置，并定义节点

csys,2

r=(f*f+span*span/4)/(2*f)

Dalfa=atn(span/2/sqrt(r*r-(span/2)*(span/2)))/Nx