怎样从ansys中提取单元刚度矩阵与质量矩阵
ansys提取质量和刚度矩阵
ansys提取质量和刚度矩阵
在ANSYS中，可以使用命令`*MATRIX, MASS, STIFF`来提取质量和刚度矩阵。
要提取质量矩阵，可以使用以下命令：
```
*MATRIX, MASS
/mesh_needed ! 控制只提取有限元网格中的质量矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后，ANSYS将会在指定的输出文件中生成质量矩阵。
要提取刚度矩阵，可以使用以下命令：
```
*MATRIX, STIFF
/mesh_needed ! 控制只提取有限元网格中的刚度矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后，ANSYS将会在指定的输出文件中生成刚度矩阵。
请注意，以上命令需要在ANSYS命令窗口中执行，并且需要
已经建立合适的有限元模型。
您还需要根据需要修改`/output`命令以指定输出文件的目录和文件名。
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。
(2021年整理)ANSYS中整体、单元刚度和质量矩阵的提取
ANSYS中整体、单元刚度和质量矩阵的提取编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（ANSYS中整体、单元刚度和质量矩阵的提取）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为ANSYS中整体、单元刚度和质量矩阵的提取的全部内容。
ANSYS中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。
1.1、用户程序法，需要二次开发该功能需要进行二次开发，由ansys形成的二进制文件.full提取整体刚度和质量矩阵。
基于ansys的一个用户开发程序例子编了一个程序（附件中）。
开发环境：compaqfortran6.5运行环境：win2000。
一个主文件self。
for；另一个文件matrixout.f90用于矩阵输出；binlib。
lib为ansys提供的库文件，将其引入项目中（也可直接扔进debug目录)；。
full文件由子空间迭代模态分析获得.运行编译后的可执行文件.exe其他文件见http：///f/EE24A2F87F524606。
html1.2、超单元法其原理很简单，即使用ansys的超单元即可解决问题。
定义超单元，然后列出超单元的刚度矩阵即可。
下面是一个小例题，自可明白。
/prep7k,1k，2，3000l，1，2et,1,beam3mp，ex，1，2e5mp，prxy,1，0。
3r，1,5000,2e7，200lesize,all,，,10lmesh，allfinish！——-—以上正常建立模型，不必施加约束和荷载/soluantype,7 ！substructuring分析类型seopt，matname,1 !设置文件名称和刚度矩阵类型（刚度，质量,阻尼等）nsel，all ！选择所有节点m,all,all !定义所有节点自由度为主自由度solve ！求解selist,matname,3 !列出整体刚度矩阵1.3、HBMAT命令法提取整体矩阵13.1、命令说明命令:HBMAT，fname,ext，—-,form,matrx，rhs其中：Fname--—输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ansys提取单元刚度矩阵
ansys提取单元刚度矩阵
本文将介绍如何使用ansys软件从有限元模型中提取单元刚度
矩阵。
在有限元分析中，单元刚度矩阵是一个重要的参数，它描述了单元受力后产生的应变响应。
通过提取单元刚度矩阵，可以进一步分析结构的刚度、应力分布和变形情况。
在ansys中，提取单元刚度矩阵的方法有多种，其中一个常用的方法是使用APDL命令。
首先，需要使用命令*SET,TYPE,1输入自由度约束来定义模型的自由度类型。
接下来，使用*SMAT命令来设置解决方案的类型为“刚度矩阵”，并使用*GET命令来提取单元刚度矩阵。
具体来说，以下是在ansys中提取单元刚度矩阵的基本步骤：
1. 定义模型的自由度类型
在命令窗口中输入命令*SET,TYPE,1来定义模型的自由度类型。
这将使得ansys自动将整个模型划分为自由度或受限制的自由度。
2. 设置解决方案类型为“刚度矩阵”
在命令窗口中输入命令*SMAT,OPER,SOLU来设置解决方案类型为“刚度矩阵”。
3. 提取单元刚度矩阵
在命令窗口中输入命令*GET,MAT,ALL来提取整个模型的单元刚度矩阵，或者*GET,MAT,ELEM,elemnum来提取特定单元（elemnum）的刚度矩阵。
以上就是在ansys中提取单元刚度矩阵的基本步骤。
通过这种方法，可以轻松地获得模型的刚度矩阵，并进一步分析结构的刚度、应
力分布和变形情况。
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
ansys提出刚度和质量矩阵
EDGE PLOT DATA = EDGE Z2 1 0.0000 2 0.0000 3 0.0000 0.0000 0.0000 3.0000 6.0000 0.0000 0.0000 0.0000 0.0000 6.0000 9.0000 0.0000 0.0000 0.0000 0.0000 3.0000 X1 Y1 Z1 X2 Y2
r,12,248500 r,13,192100 r,14,152200 type,1 *do,i,1,3,1 real,11+i e,i,i+1 *enddo !约束 nsel,s,node,,2,4,1 d,all,uz,0 d,all,uy,0 allsel d,1,all,0 /solu antype,7 seopt,submat,3 m,all,all solve selist,submat,3 /
-192100.00 2 0.0000000 1 344300.00 2 2760.0000 1 -152200.00 2 0.0000000
0.0000000
0.0000000
-152200.00
0.0000000
152200.00
2300.0000
0.0000000
0.0000000
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。
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中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。
1.1、用户程序法，需要二次开发该功能需要进行二次开发，由ansys形成的二进制文件.full提取整体刚度和质量矩阵。
基于ansys的一个用户开发程序例子编了一个程序（附件中）。
开发环境：compaqfortran6.5运行环境：win2000。
一个主文件self.for；另一个文件matrixout.f90用于矩阵输出；binlib.lib 为ansys提供的库文件，将其引入项目中（也可直接扔进debug目录）；.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe其他文件见/f/EE24A2F87F524606.html1.2、超单元法其原理很简单，即使用ansys的超单元即可解决问题。
定义超单元，然后列出超单元的刚度矩阵即可。
下面是一个小例题，自可明白。
/prep7k,1k,2,3000l,1,2et,1,beam3mp,ex,1,2e5mp,prxy,1,0.3r,1,5000,2e7,200lesize,all,,,10lmesh,allfinish!----以上正常建立模型，不必施加约束和荷载/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵1.3、HBMAT命令法提取整体矩阵13.1、命令说明命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
form---定义输出矩阵文件的格式，其值可取：=ASCII：ASCII码格式；=BIN：二进制格式。
ANSYS中整体、单元刚度和质量矩阵的提取
ANSYS中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。
1.1、用户程序法，需要二次开发该功能需要进行二次开发，由ansys形成的二进制文件.full提取整体刚度和质量矩阵。
基于ansys的一个用户开发程序例子编了一个程序（附件中）。
开发环境：compaqfortran6.5运行环境：win2000。
一个主文件self.for；另一个文件matrixout.f90用于矩阵输出；binlib.lib 为ansys提供的库文件，将其引入项目中（也可直接扔进debug目录）；.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe其他文件见/f/EE24A2F87F524606.html1.2、超单元法其原理很简单，即使用ansys的超单元即可解决问题。
定义超单元，然后列出超单元的刚度矩阵即可。
下面是一个小例题，自可明白。
/prep7k,1k,2,3000l,1,2et,1,beam3mp,ex,1,2e5mp,prxy,1,0.3r,1,5000,2e7,200lesize,all,,,10lmesh,allfinish!----以上正常建立模型，不必施加约束和荷载/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵1.3、HBMAT命令法提取整体矩阵13.1、命令说明命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
form---定义输出矩阵文件的格式，其值可取：=ASCII：ASCII码格式；=BIN：二进制格式。
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提取单元刚度矩阵
ansys提取单元刚度矩阵
在使用ANSYS进行有限元分析时，我们通常需要提取单元的刚度矩阵。
刚度矩阵是描述单元刚度的一个重要参数，它可以用来计算单元的应力、应变和位移等参数。
提取单元刚度矩阵的方法如下：
1. 首先，在ANSYS中建立模型并进行网格划分。
确定需要提取刚度矩阵的单元，例如在ANSYS中选择PLANE42单元。
2. 在ANSYS命令窗口中输入“/POST1”命令，进入后处理模式。
3. 输入“ET,LIST”命令，显示所有单元类型的列表。
根据需要选择需要提取刚度矩阵的单元类型。
4. 输入“SET,LIST”命令，选择需要提取刚度矩阵的单元集合。
5. 输入“MAT,LIST”命令，选择需要提取刚度矩阵的材料属性。
6. 输入“SGLSTIF,ALL”命令，提取所有单元的刚度矩阵。
也可以使用“SGLSTIF,LIST”命令，提取指定单元的刚度矩阵。
7. 输入“PRRSOL,U”命令，显示单元的未知位移。
8. 输入“PRRSOL,F”命令，显示单元的未知力。
通过以上步骤，我们就可以成功提取单元的刚度矩阵，并使用它来计算单元的应力、应变和位移等重要参数。
- 1 -。
Ansysworkbench17.0提取刚度矩阵的方法
设置文件名称和刚度矩阵类型刚度质量阻尼等nsel3756
Ansysworkbench17.0提 取 刚 度 矩 阵 的 方 法
第一步：导入模型，划分网格 第二步：定义remote point，位置设置在想要提取刚度矩阵的位置。
第三步：添加remote displacement约束，给第二步定义的点添加约束。（本例为全约束）
第四步：在static structural处插入command，内容如下： /solu antype,7 !substructuring分析类型 seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等) nsel,,,,3756 ！此处节点数（3756）由mesh中statistics 查到，nodes 值为3755，则remote point为3756、3757……m,all,all !选择所有节点 allsel !定义所有节点自由度为主自由度 solve !求解 selist,matname,3 !列出整体刚度矩阵 finish 第五步：solve 第六步：在solution中solution information查看结果（结果按照行分布）
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文件，可以看到按单元写出的质量、刚度等矩阵方法2/prep7k,1k,2,3000l,1,2et,1,beam3mp,ex,1,2e5mp,prxy,1,0.3r,1,5000,2e7,200lesize,all,,,10lmesh,allfinish!----以上正常建立模型，不必施加约束和荷载/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等) nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
下面提供三种方法：方便与其他程序进行接口编程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.full
/OUTPUT, TERM ! 将输出信息送到output windows中
! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵
3.
其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。
面是一个小例题，自可明白。
/prep7
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵
3、HBMAT命令法提取整体矩阵
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
第2行：格式（5I14），分别表示该文件的总行数（不包括文件头）、矩阵列指针的总行数、矩阵行索引的总行数、矩阵元素数值的总行数、右边项总行数。
第3行：格式（A3,11X,4I14），分别为矩阵类型、矩阵行数、矩阵列数、矩阵行索引数（对组装后的矩阵，该值等于矩阵行索引数）、单元元素数（对组装后的矩阵此值为0）。
第6行后：矩阵元素值（单列）。
矩阵类型用3个字符表示，第1个字符可取：R---实数矩阵、C---复数矩阵、P---仅矩阵结构（无元素数值）；第2个字符可取：S---对称矩阵、U---不对称矩阵、H---Hermitian矩阵、Z---病态对称矩阵；R---带状矩阵；第3个字符可取：A---组装的矩阵、E---单元矩阵（未组装）。对称矩阵只存储下三角元素，如结构刚度矩阵为对称矩阵，Harwell-Boeing格式则仅记录下三角元素。
seopt
回复 2# 的帖子
/solu
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
UIMP,1,DENS, , ,7800,
BLC4,0,0,w2,w1
ESIZE,1.5,0,
AMESH,all
NSEL,S,LOC,X,0.0
=BIN：二进制格式。
matrix---定义输出矩阵的类型，其值可取：
=STIFF：输出刚度矩阵。可用于写入了.FULL文件的任何类型的分析。
=MASS：输出质量矩阵。可用于特征值屈曲、子结构分析、模态分析。
t=0.08
/PREP7
!*
ET,1,SHELL63
R,1,t
ET,2,MASS21
R,2,500,500,500,2000,2000,2000,
!*
UIMP,1,EX, , ,2e11
UIMP,1,NUXY, , ,0.3,
UIMP,1,DAMP, , ,0.2,
Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
ext---输出矩阵文件的扩展名，缺省为.matrix。
form---定义输出矩阵文件的格式，其值可取：
=ASCII：ASCII码格式；
*do,irow,stacol,endcol-1
truerow=rowind(irow)
smatr(truerow,icol)=values(irow)
*enddo
*enddo
*do,irow,1,nrow
*do,icol,1,ncol
smatr(irow,icol)=smatr(icol,irow)
⑵ Harwell-Boeing文件格式
用HBMAT命令可输出结构刚度矩阵、质量矩阵和阻尼矩阵，其文件记录格式为大型稀疏矩阵的标准交换格式，采用索引存储方法仅记录矩阵的非零元素。文件基本格式是前面有4或5行描述数据，其后为单列矩阵元素值，说明如下：
第1行：格式（A72），为文件头的字符型解释，如刚度矩阵或质量矩阵等标题。
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件
/debug,-1,,,1 ! 指定输出单元矩阵
/SOLU
SOLVE
finish
/OUTPUT, TERM ! 将输出信息送到output windows中
这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵
selist,matname,3 !列出整体刚度矩阵
ANSYS单元和整体刚度矩阵的提取
一、单元刚度矩阵的提取
/DEBUG命令
详细说明：
finish
/clear
PI=3.1415926
w1=3
w2=10
w3=6
w4=1.2
r=.8
=DAMP：输出阻尼矩阵。仅用于有阻尼的模态分析。
rhs---右边项输出控制（右边项指用矩阵所表示方程的等号右端矢量，这里可为节点荷载向量），如rhs=YES则输出，如rhs=NO则不输出。
模态分析时，因仅LANB和QR法可生成完整的质量矩阵，因此也仅采用这两种方法时才可使用HBMAT命令得到质量矩阵文件。
命令流：
/aux2
file,hbfile,full
hbmat,hbfile,txt,,ascii,stiff,yes f Nhomakorabeanish
!从hbfile.txt读入数据，并还原为满矩阵存储
*dim,contline,,5
*vread,contline(1),hbfile,txt,,,5,,,1
A simple way to dump the matrix is as follow:
-------------------
/aux2
fileaux2,file,emat
form,long
dump,all
-------------------
2.
可以使用/DEBUG命令来得到。详细步骤参见下面的宏文件
*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
*enddo
*enddo
pointr=
rowind=
values=
rhsval=
icol=
irow=
ls0=
stacol=
endcol=
truerow=
totcrd=
ptrcrd=
indcrd=
valcrd=
rhscrd=
以上这些是我总结的，结果不同就错了。
(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)
根据Harwell-Boeing文件格式，可读取矩阵的任意行列元素的数值，也可编程还原为满矩阵存储，以便它用，很显然这种提取方式比较方便。如当生成.FULL文件后，可采用命令/AUX2$FILE,mywork,full$HBMAT,mystiff,txt,ASCII,STIFF,YES$FINISH将二进制mywork.full文件输出为ASCII码文件mystiff.txt，并输出右边项。
二、整体刚度矩阵的提取（有三种方法：用户程序法、超单元法、HBMAT命令法）
1、用户程序法：需要二次开发（略）
2、超单元法
/solu
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nrow=contline(2)
ncol=contline(3)
strline=
contline=
*if,rhscrd,eq,0,then
ls0=4
*else
ls0=5
*endif
*dim,pointr,,ptrcrd
*dim,rowind,,indcrd
*dim,values,,valcrd
D,all, , , , , ,ALL, , , , ,
allsel,all
SFA,all,1,PRES,12
FINISH
/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件