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中单元刚度矩阵、总刚矩阵的提取

/solu
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提取单元刚度矩阵

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提出刚度和质量矩阵

filnambeammasskniffanalysistitlebeammasskniffanalysisprep7et1beam3mpex120e11mpnuxy103mpdens17800sectype1beami3secoffsetcentsecdata015015030020020010000k21000一框架的刚度和质量矩阵的提取模型采用右图的集中质量模型m12762kgm22760kgm32300kg层间刚度分别为
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质量矩阵刚度矩阵提取教学提纲

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

一、连续钢梁的刚度和质最矩阵弹性模量：
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
/PREP7
EI1BEAM3
MBEX470E11
MRNUXYX0.3
MRDENS 丄7800
SECTYPE,1,BEAM」,,3
SECO 阡SET.CENT
SECDATAO150・：15,030・02O02,0・(HO0O0 K10O0
K210O0
一、框架的刚度和质呈:矩阵的提取模型采用右图的集中
质屋模型，Ml=2762Kg ,
M2=2760Kg, M3=2300Kg,层间刚度分别为：248500、192100. 152200
命令流如下:
finish
/clear
/prep7 et4z combinl4
keyopt,1,2,1 !—
维
弹簧单元（Ux平动）
et2mass21 keyopt,2,3/2 !3・D mass没有转动惯最
0,1041-1）*3,0 *enddo
!质最常数
r z2,2762
r,3,2760
54,2300
type,2
*do,i,2,4,l
realj
e,i
*enddo
!层间刚度^12,248500
G1*******
G14J52200 type」•doj,134 realjl+i eJJ+1 *enddo 哟束nselsnode,24JL d，all,uz，0
d,all,uy,0 allsel d,l,all,0 /solu antype z7 seoptsubmatB m,all z all solve selistsubmat.B
使用该命令流町以得到结构的刚度矩阵和质最矩阵结果: 刚度矩阵：
440600. 00 -
192100. 00
0. 0000000 质量矩阵：2762. 0000
0. 0000000
0. 0000000
•192100.
00
344300. 00
-152200.00
0. 0000000
-152200. 00
152200. 00
0.0000000
2760. 0000
0.0000000
0.0000000
0.0000000
2300. 0000
im
IT1
2
mi
k-
ki
附件1：（运行ANSYS命令的输出结果，最后面是刚度和质最矩阵）
PRINT CONTENTS OF SUPERELEMENT submat
PRINT OPTION = 3
HEADER =
8 3 2 3 3
4 0 4 3 0
1 1 0 0 1
0 0 1 12 1155
0 1101 166 103 subm
at 1078 0
922 928 946 953 976
982 1027 1072 0 934
940 106004406 0 0 3
0 0
419 0 0 0 0
0 0 0 0
HEADER SUMMARY:
NUMBER OF ROWS = 3
NUMBER OF MATRICES = 2
NUMBER OF EDGE PLOT LINES = 3
NUMBER OF DEGREES OF FREEDOM PER NODE = 3
NUMBER OE DEGREE OF FREEDOM INDICES = 4
NUMBER OF NODES = 3
NUMBER OF LOAD VECTORS = 1
NUMBER OE TRANSFORMATIONS = 0
BASE FILE NAME二submat
MAXIMUM STIFFNESS二0. 44060E+06
DEGREES OF FREEDOM PER NODE =
DEGREES OF FREEDOM PER NODE LABELS =
UX UY UZ
DEGREES OF FREEDOM AS GENERATED =
4 7 10
REORDERED DOF POSITIONS
ORIGINAL DOF ORDER =
DEGREE OF FREEDOM INDICES =
1 2 3 4 TITLE =
NODES =
2
3 4
NODE X Y THXZ
2 0. 0000 3. 0000 0. 0000
3 0. 0000 6. 0000 0. 0000
4 0. 0000 9. 0000 0. 0000
EDGE PLOT DATA 二
EDGE XI Y1
Z2
1 0. 0000 0. 0000 0. 0000
2 0. 0000 3. 0000 0. 0000
3 0. 0000 6. 0000 0. 0000
GLOBAL DOF SET 二
33 65 97
GLOBAL DOF SET NODES, LABELS
2 UX
3 UX
4 7 10
MASS INFORMATION:
Z THXY THYZ 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000
Z1 X2 Y2
0. 0000 0. 0000 3. 0000 0. 0000 0. 0000 6. 0000 0. 0000 0. 0000 9. 0000
4 UX
TOTAL MASS = 7822.0
CENTROID (X,Y,Z) = 0. 0000 5. 8228 MOMENT OE INTERTIA ABOUT ORIGIN:
IXX = 0.31052E+06 IYY = 0. 0000
IXY = 0.0000 IYZ = 0.0000 0. 0000
1ZZ = 0.31052E+06 IZX = 0. 0000
RON 1 MATRIX
440600. 00 -192100. 00 0. 0000000
ROW 1 MATRIX 2
2762. 0000 0. 0000000 0. 0000000
ROW 2 MATRIX 1
-192100. 00 344300. 00 -152200. 00
ROW 2 MATRIX 2
0. 0000000 2760. 0000 0. 0000000
ROW 3 MATRIX 1
0. 0000000 -152200. 00 152200. 00
ROW 3 MATRIX 2
0. 0000000 0. 0000000 2300. 0000
LOAD VECTOR 1
0. 0000000 0. 0000000 0. 0000000。