ansys提出刚度和质量矩阵
- 格式:docx
- 大小:16.82 KB
- 文档页数:6
ansys提取质量和刚度矩阵
在ANSYS中,可以使用命令`*MATRIX, MASS, STIFF`来提取质量和刚度矩阵。
要提取质量矩阵,可以使用以下命令:
```
*MATRIX, MASS
/mesh_needed ! 控制只提取有限元网格中的质量矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后,ANSYS将会在指定的输出文件中生成质量矩阵。
要提取刚度矩阵,可以使用以下命令:
```
*MATRIX, STIFF
/mesh_needed ! 控制只提取有限元网格中的刚度矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后,ANSYS将会在指定的输出文件中生成刚度矩阵。
请注意,以上命令需要在ANSYS命令窗口中执行,并且需要
已经建立合适的有限元模型。
您还需要根据需要修改`/output`命令以指定输出文件的目录和文件名。
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中提取单元刚度矩阵的基本步骤。
通过这种方法,可以轻松地获得模型的刚度矩阵,并进一步分析结构的刚度、应
力分布和变形情况。
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中单元刚度矩阵、总刚矩阵的提取一、单元刚度矩阵的提取/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。
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单元和整体刚度矩阵的提取一、单元刚度矩阵的提取/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中导出刚度矩阵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文件,可以看到按单元写出的质量、刚度等矩阵干部教育培训工作总结[干部教育培训工作总结] 年干部教育培训工作,在县委的正确领导下,根据市委组织部提出的任务和要求,结合我县实际,以兴起学习贯彻“三个代表”重要思想新高潮为重点,全面启动“大教育、大培训”工作,取得了一定的成效,干部教育培训工作总结。
现总结报告如下:一、基本情况全县共有干部**人,其中中共党员**人,大学本科以上学历**人,大专学历**人,中专学历**人,高中及以下学历**人。
**年,以县委党校、县行政学校为主阵地,举办各类培训**期,培训在职干部**人,占在职干部总数的**.*%,培训农村党员、干部**人,其中:举办科级领导干部轮训班*期,培训**人;举办科级领导干部“三个代表”重要思想专题学习班*期,培训**人;举办科级以下公务员培训班*期,培训**人;举办企业经营管理者培训班*期,培训**人;举办专业技术人员培训班*期,培训**人;举办非中共党员干部培训班*期,培训**人;举办理论骨干培训班*期,培训**人;举办妇女干部培训班*期,培训**人;举办基层团干培训班*期,培训**人;举办农村党支部书记、村主任培训班各*期,培训**人,达到了每年培训在职干部五分之一的要求,超额完成了培训任务。
一、连续钢梁的刚度和质最矩阵弹性模量:
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。