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。