ansys质量矩阵刚度矩阵提取教学提纲
a n s y s质量矩阵刚度矩阵提取ansys质量矩阵刚度矩阵提取看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程1.Which matrix you would like? element stiffness matrix or full stiffness matrix?element stiffness is within file.emat. full stiffness matrix is within file.fullA simple way to dump the matrix is as follow:-------------------/aux2fileaux2,file,ematform,longdump,all-------------------2.可以使用/DEBUG命令来得到。
详细步骤参见下面的宏文件finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵3.其原理很简单，即使用ansys的超单元即可解决问题。
定义超单元，然后列出超单元的刚度矩阵即可。
面是一个小例题，自可明白。
/prep7k,1k,2,3000l,1,2et,1,beam3mp,ex,1,2e5mp,prxy,1,0.3mp,dens,1,2e3r,1,5000,2e7,200lesize,all,,,10lmesh,allfinish!----以上正常建立模型，不必施加约束和荷载/soluantype,7 !substructuring分析类型seopt,matname,2 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等) nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵ANSYS单元和整体刚度矩阵的提取一、单元刚度矩阵的提取/DEBUG命令详细说明：finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵二、整体刚度矩阵的提取（有三种方法：用户程序法、超单元法、HBMAT命令法）1、用户程序法：需要二次开发（略）2、超单元法/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵3、HBMAT命令法提取整体矩阵命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
form---定义输出矩阵文件的格式，其值可取：=ASCII：ASCII码格式；=BIN：二进制格式。
matrix---定义输出矩阵的类型，其值可取：=STIFF：输出刚度矩阵。
可用于写入了.FULL文件的任何类型的分析。
=MASS：输出质量矩阵。
可用于特征值屈曲、子结构分析、模态分析。
=DAMP：输出阻尼矩阵。
仅用于有阻尼的模态分析。
rhs---右边项输出控制（右边项指用矩阵所表示方程的等号右端矢量，这里可为节点荷载向量），如rhs=YES则输出，如rhs=NO则不输出。
模态分析时，因仅LANB和QR法可生成完整的质量矩阵，因此也仅采用这两种方法时才可使用HBMAT命令得到质量矩阵文件。
⑵ Harwell-Boeing文件格式用HBMAT命令可输出结构刚度矩阵、质量矩阵和阻尼矩阵，其文件记录格式为大型稀疏矩阵的标准交换格式，采用索引存储方法仅记录矩阵的非零元素。
文件基本格式是前面有4或5行描述数据，其后为单列矩阵元素值，说明如下：第1行：格式（A72），为文件头的字符型解释，如刚度矩阵或质量矩阵等标题。
第2行：格式（5I14），分别表示该文件的总行数（不包括文件头）、矩阵列指针的总行数、矩阵行索引的总行数、矩阵元素数值的总行数、右边项总行数。
第3行：格式（A3,11X,4I14），分别为矩阵类型、矩阵行数、矩阵列数、矩阵行索引数（对组装后的矩阵，该值等于矩阵行索引数）、单元元素数（对组装后的矩阵此值为0）。
第4行：格式（2A16,2A20），分别表示列指针格式、行索引格式、系数矩阵数值格式、右边项数值格式。
第5行：格式（A3,11X,2I14），A3各列分别表示右边项格式、应用高斯起始矢量、应用eXact求解矢量；两个整数分别表示右边项列数、行索引数。
三个字符中的第1个字符可取：F---全部存贮（如节点荷载向量的全部元素）、M---与系数矩阵相同方法。
第6行后：矩阵元素值（单列）。
矩阵类型用3个字符表示，第1个字符可取：R---实数矩阵、C---复数矩阵、P---仅矩阵结构（无元素数值）；第2个字符可取：S---对称矩阵、U---不对称矩阵、H---Hermitian矩阵、Z---病态对称矩阵；R---带状矩阵；第3个字符可取：A---组装的矩阵、E---单元矩阵（未组装）。
对称矩阵只存储下三角元素，如结构刚度矩阵为对称矩阵，Harwell-Boeing格式则仅记录下三角元素。
根据Harwell-Boeing文件格式，可读取矩阵的任意行列元素的数值，也可编程还原为满矩阵存储，以便它用，很显然这种提取方式比较方便。
如当生成.FULL文件后，可采用命令/AUX2$FILE,mywork,full$HBMAT,mystiff,txt,ASCII,STIFF,YES$FINISH将二进制mywork.full文件输出为ASCII码文件mystiff.txt，并输出右边项。
命令流：/aux2file,hbfile,fullhbmat,hbfile,txt,,ascii,stiff,yesfinish!从hbfile.txt读入数据，并还原为满矩阵存储*dim,contline,,5*vread,contline(1),hbfile,txt,,,5,,,1(5f14.0)totcrd=contline(1)ptrcrd=contline(2)indcrd=contline(3)valcrd=contline(4)rhscrd=contline(5)*vread,contline(1),hbfile,txt,,,4,,,2(a3,11x,4f14.0)nrow=contline(2)ncol=contline(3)strline=contline=*if,rhscrd,eq,0,thenls0=4*elsels0=5*endif*dim,pointr,,ptrcrd*dim,rowind,,indcrd*dim,values,,valcrd*dim,rhsval,,rhscrd*vread,pointr(1),hbfile,txt,,,ptrcrd,,,ls0(f14.0)*vread,rowind(1),hbfile,txt,,,indcrd,,,ls0+ptrcrd(f14.0)*vread,values(1),hbfile,txt,,,valcrd,,,ls0+ptrcrd+indcrd(d25.15)*vread,rhsval(1),hbfile,txt,,,rhscrd,,,ls0+ptrcrd+indcrd+valcrd(d25.15)*dim,smatr,,nrow,ncol*do,icol,1,ncolstacol=pointr(icol)endcol=pointr(icol+1)*do,irow,stacol,endcol-1truerow=rowind(irow)smatr(truerow,icol)=values(irow)*enddo*enddo*do,irow,1,nrow*do,icol,1,ncolsmatr(irow,icol)=smatr(icol,irow)*enddo*enddopointr=rowind=values=rhsval=icol=irow=ls0=stacol=endcol=truerow=totcrd=ptrcrd=indcrd=valcrd=rhscrd=以上这些是我总结的，结果不同就错了。
4./soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等) nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵5.主要通过/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵单元刚度和质量矩阵的提取。
Ansys中单元刚度矩阵、总刚矩阵的提取
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵
(d25.15)
*vread,rhsval(1),hbfile,txt,,,rhscrd,,,ls0+ptrcrd+indcrd+valcrd
(d25.15)
*dim,smatr,,nrow,ncol
*do,icol,1,ncol
stacol=pointr(icol)
endcol=pointr(icol+1)
第4行：格式（2A16,2A20），分别表示列指针格式、行索引格式、系数矩阵数值格式、右边项数值格式。
第5行：格式（A3,11X,2I14），A3各列分别表示右边项格式、应用高斯起始矢量、应用eXact求解矢量；两个整数分别表示右边项列数、行索引数。三个字符中的第1个字符可取：F---全部存贮（如节点荷载向量的全部元素）、M---与系数矩阵相同方法。
3、HBMAT命令法提取整体矩阵
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
form---定义输出矩阵文件的格式，其值可取：
=ASCII：ASCII码格式；
=BIN：二进制格式。
ansys提出刚度和质量矩阵
一、连续钢梁的刚度和质最矩阵弹性模量：EX=2.0x 10n Pa 泊松比：PRXY=0.3 密度：DENS=7811Kg/m3 截面特性如右图。
AN SYS命令流：/FILNAM.BEAM MASS AND KNIFF ANALYSIS/TITLE, BEAM MASS AND KNIFF ANALYSIS/PREP7EI1BEAM3MBEX470E11MRNUXYX0.3MRDENS 丄7800SECTYPE,1,BEAM」,,3SECO 阡SET.CENTSECDATAO150・：15,030・02O02,0・(HO0O0 K10O0K210O0一、框架的刚度和质呈:矩阵的提取模型采用右图的集中质屋模型，Ml=2762Kg ,M2=2760Kg, M3=2300Kg,层间刚度分别为：248500、192100. 152200命令流如下:finish/clear/prep7 et4z combinl4keyopt,1,2,1 !—维弹簧单元（Ux平动）et2mass21 keyopt,2,3/2 !3・D mass没有转动惯最0,1041-1）*3,0 *enddo!质最常数r z2,2762r,3,276054,2300type,2*do,i,2,4,lrealje,i*enddo!层间刚度^12,248500G1*******G14J52200 type」•doj,134 realjl+i eJJ+1 *enddo 哟束nselsnode,24JL d，all,uz，0d,all,uy,0 allsel d,l,all,0 /solu antype z7 seoptsubmatB m,all z all solve selistsubmat.B使用该命令流町以得到结构的刚度矩阵和质最矩阵结果: 刚度矩阵：440600. 00 -192100. 000. 0000000 质量矩阵：2762. 00000. 00000000. 0000000•192100.00344300. 00-152200.000. 0000000-152200. 00152200. 000.00000002760. 00000.00000000.00000000.00000002300. 0000imIT12mik-ki附件1：（运行ANSYS命令的输出结果，最后面是刚度和质最矩阵）PRINT CONTENTS OF SUPERELEMENT submatPRINT OPTION = 3HEADER =8 3 2 3 34 0 4 3 01 1 0 0 10 0 1 12 11550 1101 166 103 submat 1078 0922 928 946 953 976982 1027 1072 0 934940 106004406 0 0 30 0419 0 0 0 00 0 0 0HEADER SUMMARY:NUMBER OF ROWS = 3NUMBER OF MATRICES = 2NUMBER OF EDGE PLOT LINES = 3NUMBER OF DEGREES OF FREEDOM PER NODE = 3NUMBER OE DEGREE OF FREEDOM INDICES = 4NUMBER OF NODES = 3NUMBER OF LOAD VECTORS = 1NUMBER OE TRANSFORMATIONS = 0BASE FILE NAME二submatMAXIMUM STIFFNESS二0. 44060E+06DEGREES OF FREEDOM PER NODE =DEGREES OF FREEDOM PER NODE LABELS =UX UY UZDEGREES OF FREEDOM AS GENERATED =4 7 10REORDERED DOF POSITIONSORIGINAL DOF ORDER =DEGREE OF FREEDOM INDICES =1 2 3 4 TITLE =NODES =23 4NODE X Y THXZ2 0. 0000 3. 0000 0. 00003 0. 0000 6. 0000 0. 00004 0. 0000 9. 0000 0. 0000EDGE PLOT DATA 二EDGE XI Y1Z21 0. 0000 0. 0000 0. 00002 0. 0000 3. 0000 0. 00003 0. 0000 6. 0000 0. 0000GLOBAL DOF SET 二33 65 97GLOBAL DOF SET NODES, LABELS2 UX3 UX4 7 10MASS INFORMATION:Z THXY THYZ 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000Z1 X2 Y20. 0000 0. 0000 3. 0000 0. 0000 0. 0000 6. 0000 0. 0000 0. 0000 9. 00004 UXTOTAL MASS = 7822.0CENTROID (X,Y,Z) = 0. 0000 5. 8228 MOMENT OE INTERTIA ABOUT ORIGIN:IXX = 0.31052E+06 IYY = 0. 0000IXY = 0.0000 IYZ = 0.0000 0. 00001ZZ = 0.31052E+06 IZX = 0. 0000RON 1 MATRIX440600. 00 -192100. 00 0. 0000000ROW 1 MATRIX 22762. 0000 0. 0000000 0. 0000000ROW 2 MATRIX 1-192100. 00 344300. 00 -152200. 00ROW 2 MATRIX 20. 0000000 2760. 0000 0. 0000000ROW 3 MATRIX 10. 0000000 -152200. 00 152200. 00ROW 3 MATRIX 20. 0000000 0. 0000000 2300. 0000LOAD VECTOR 10. 0000000 0. 0000000 0. 0000000。
ansys提取刚度矩阵案例(命令流及矩阵文件)
ansys提取刚度矩阵的三种方法注：本案例借鉴了王新敏老师编著的《ansys工程结构数值分析》部分内容结构刚度矩阵K(1)单元刚阵的提取-----详见page356/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLU/OUTPUT,ELEMSTIFF,TXT/DEBUG,-1,,,1SOLVE/OUTPUTFINISH用任一文本编辑器打开ELEMSTIFF.TXT文件可得到单元刚度矩阵ELEMSTIFF.TXT(2)用HBMAT提取原始刚度矩阵和节点载荷FINISH/CLEAR/FILNAME，HBFILE/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,HBFILE,FULLHBMAT,HBFILE,TXT,ASCII,STIFF,YESFINISH用文本编辑器打开HBFILE.TXT可看到用Harwell-Boeing格式记录的文件（仅含非零项的上三角阵，刚度矩阵为对称矩阵）。
矩阵有12行12列，33个非零元素HBFILE.TXT（3）用HBMAT提取结构刚度矩阵与（2）相同，但施加约束条件即可。
所生成的HBMAT.TXT的前5行为/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,文件名缺省,FULLHBMAT,文件名缺省,TXT,ASCII,STIFF,YESFINISH其意义同上，但数值有变化。
ansys提出刚度和质量矩阵
2762.0000 0.0000000 0.0000000
附件 1：（运行 ANSYS 命令的输出结果，最后面是刚度和质量矩阵） PRINT CONTENTS OF SUPERELEMENT submat PRINT OPTION = 3 HEADER = 8 4 1 0 0 at 922 982 940 0 419 0 3 0 1 0 1101 1078 928 1027 106004406 0 0 2 4 0 1 166 0 946 1072 0 0 0 3 3 0 12 103 953 0 0 0 0 3 0 1 1155 subm 976 934 3 0 0
一、连续钢梁的刚度和质量矩阵 弹性模量：EX=2.0× 1011 Pa 泊松比：PRXY=0.3 密度：DENS=7811Kg/m3 截面特性如右图。 ANSYS 命令流： /FILNAM,BEAM MASS AND KNIFF ANALYSIS /TITLE, BEAM MASS AND KNIFF ANALYSIS /PREP7 ET,1,BEAM3 MP,EX,1,2.0E11 MP,NUXY,1,0.3 MP,DENS,1,7800 SECTYPE,1,BEAM,I,,3 SECOFFSET,CENT SECDATA,0.15,0.15,0.3,0.02,0.02,0.01,0,0,0,0 K,1,0,0,0 K,2,10,0,0
GLOBAL DOF SET = 33 65 97 GLOBAL DOF SET NODES, LABELS = 2 UX 3 UX MASS INFORMATION: TOTAL MASS = 7822.0 CENTROID (X,Y,Z) = 0.0000 5.8228 MOMENT OF INTERTIA ABOUT ORIGIN: IXX = 0.31052E+06 IYY = 0.0000 IXY = 0.0000 IYZ = 0.0000 ROW 1 MATRIX 1
提取刚度矩阵
============为什么折腾这个文档========我有一个计算线性动力学方程组的瞬态、谐响应和静力学的python程序，现希望开发一个将ANSYS组集好的总体矩阵导入该PYTHON程序中的接口。
该问题可分解为：[STEP1] [ANSYS]->[包含矩阵信息的文件][STEP2] [包含矩阵信息的文件]->[python通用数据对象][STEP3] [python通用数据对象]->[程序特定数据对象]->[进行计算]因此检索了一些帖子，基本上完成了这项工作，本文是对[STEP1]和[STEP2]的整理，并且利用[STEP3]对结果进行了验证============主要内容==================1，了解从ANSYS中提取总体矩阵和载荷向量的方法；2，了解提取出来的矩阵是怎样表示的；3，说明在Python中，如何读取这样的矩阵；4，构造一个简单的算例，说明整个【建模】-【提取】-【读取】过程及其正确性；=========站内检索综述====================检索词：提取矩阵得到21个结果，代表性的帖子有下面这9个：编号[1]标题：ansys中怎样提取质量，刚度，阻尼矩阵？地址：/forum-vi ... fromuid-159019.html 要点：pengweicai给出了一段网上最常见的提取代码,该程序以fortran 写成，可以利用.full文件以及一些列约定将ANSYS中的总体矩阵读入FORTRAN中。
编号[2]标题：如何得知HBMAT命令提取的质量、刚度矩阵对应的自由度？地址：/forum-vi ... fromuid-159019.html 要点：提出了使用HBMAT命令提取稀疏矩阵时常见的问题：我们如何知道提取出来的信息是怎么储存的呢？编号[3]标题：[分享]ANSYS中整体、单元刚度和质量矩阵的提取地址：/forum-vi ... fromuid-159019.html 要点：在该帖子的7楼，其实已经给出了帖子[2]中问题的解答，即HBMAT 中提取出来的矩阵是Harwell-Boeing格式的，并且给出了该格式的细节，可惜是英文的，没引起多少关注。
apdl提取刚度矩阵
APDL提取刚度矩阵1. 引言在有限元分析中，刚度矩阵是一个重要的概念。
刚度矩阵描述了结构体系的刚度特性，是进行结构分析和求解的基础。
在ANSYS的APDL（ANSYS Parametric Design Language）中，我们可以通过一系列命令和操作来提取刚度矩阵。
本文将详细介绍如何在APDL中提取刚度矩阵，并探讨相关的技巧和注意事项。
2. APDL提取刚度矩阵的基本步骤提取刚度矩阵的基本步骤如下： 1. 创建几何模型：在APDL中，首先需要创建一个几何模型，可以使用命令或者图形界面进行建模。
2. 定义材料属性：根据实际情况，定义材料的弹性模量、泊松比等力学性质。
3. 定义单元类型：选择适当的单元类型，如梁单元、壳单元或体单元，并设置相应的单元属性。
4. 定义边界条件：定义结构的边界条件，如约束和加载。
5. 生成刚度矩阵：使用APDL提供的命令，生成刚度矩阵。
6. 导出刚度矩阵：将刚度矩阵导出到外部文件进行后续分析。
3. 创建几何模型在APDL中，可以使用命令行或者图形界面来创建几何模型。
命令行方式更加灵活，可以通过输入命令来精确控制几何体的生成。
图形界面方式则更加直观，通过鼠标操作可以创建各种几何体。
4. 定义材料属性在APDL中，可以使用*MAT命令来定义材料属性。
根据实际情况，需要输入材料的弹性模量、泊松比等力学性质。
这些参数将用于计算刚度矩阵。
5. 定义单元类型在APDL中，可以使用*ELEMENT命令来定义单元类型。
根据实际情况选择适当的单元类型，如梁单元、壳单元或体单元，并设置相应的单元属性。
不同类型的单元具有不同的自由度和刚度特性，选择合适的单元类型对于提取准确的刚度矩阵至关重要。
6. 定义边界条件在APDL中，可以使用*BOUNDARY命令来定义结构的边界条件。
边界条件包括约束和加载，用于模拟实际结构的受力和支撑情况。
正确定义边界条件对于提取准确的刚度矩阵非常重要。
ansys质量矩阵刚度矩阵提取说课材料
a n s y s质量矩阵刚度矩阵提取ansys质量矩阵刚度矩阵提取看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程1.Which matrix you would like? element stiffness matrix or full stiffness matrix?element stiffness is within file.emat. full stiffness matrix is within file.fullA simple way to dump the matrix is as follow:-------------------/aux2fileaux2,file,ematform,longdump,all-------------------2.可以使用/DEBUG命令来得到。
详细步骤参见下面的宏文件finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵3.其原理很简单，即使用ansys的超单元即可解决问题。
Ansys中单元刚度矩阵总刚矩阵的提取
Ansys中单元刚度矩阵、总刚矩阵的提取一、单元刚度矩阵的提取/DEBUG命令详细说明：finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵二、整体刚度矩阵的提取（有三种方法：用户程序法、超单元法、HBMAT命令法）1、用户程序法：需要二次开发（略）2、超单元法/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵3、HBMAT命令法提取整体矩阵命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
提取刚度矩阵
1在ANSYS建立有限元模型2建立好模型后，求解时，采用下面的命令流/soluallselantype,staticeqslv,sparsensubst,1wrfull,1ematwrite,yessolvefini求解后，得到的文件包括.emat .full 文件将这两个文件拷贝到mor4ansys.exe 所在文件夹下mor4ansys.exe的路径为：本机D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows\mor4ansys.exe在DOS系统中进入D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows 文件夹在文件夹1．进入D盘D：回车2 进入D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows 文件夹在文件夹键入cd D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows并回车然后键入：mor4ansys -f 文件名.full 文件名.emat –w例如我用ANSYS提出的文件分别为即文件名为sphere然后键入：mor4ansys -f sphere.full sphere.emat –w 回车即可得出文件然后将mminfo mmread mmwrite MATLAB_READ TTFORM 等文件拷贝到该文件夹打开matlab执行MATLAB_READ命令得出K.mat M.mat文件在打开TTFORM.m文件编辑for k=1:2316TT(fullmodel(k,1),fullmodel(k,2))=fullmodel(k,3);End其中2316为结构的自由度总数，可根据实际情况修改修改后保存，并运行，得出TT变量导入K.mat M.mat文件进行运算K1=T*M*T’M1=T*M*T’得出的K1，M1即为刚度，质量矩阵。
NO2 如何提取Ansys整体刚度阵及转换成Matlab数据
function [a, rhsval] = hb_to_msm (input_file) M、C、K：稀疏矩阵； rhsval：载荷右端项
input_file：步骤2导出的数据文件
6. 利用 matlab 程序“readAnsysSort2.m”即可读取 EMAT 文件，程序见附录 3
function NodeDofs = readAnsysSort2(file1,file2) file1：通过5导出的EMAT文件，即“DUMP.lis” file2：模型约束文件
节点号 约束自由度
4. 该方法适用于大型结构，复杂模型。 5. 该方法生成的整体刚度阵，是按优化后节点序列排序的，需要找到各节点在刚度阵中的
对应的位置，此处通过 EMAT 文件获取，其中 EMAT 文件的格式说明见附录。
FINISH /AUX2 FORM,LONG FILEAUX2,'file','emat',' ' DUMP,5, FINISH
FINISH /AUX2 FILE,'file','full',' ' HBMAT,'k','txt',' ',ASCII,STIFF,YES FINISH
HBMAT,'m','txt',' ',ASCII,MASS,YES HBMAT,'c','txt',' ',ASCII,DAMP,NO YES-列出载荷右端项；NO-不列出载荷右端项；