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中单元刚度矩阵、总刚矩阵的提取
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提出刚度和质量矩阵
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单元和整体刚度矩阵的提取一、单元刚度矩阵的提取/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。
2016新编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文件,可以看到按单元写出的质量、刚度等矩阵干部教育培训工作总结[干部教育培训工作总结] 年干部教育培训工作,在县委的正确领导下,根据市委组织部提出的任务和要求,结合我县实际,以兴起学习贯彻“三个代表”重要思想新高潮为重点,全面启动“大教育、大培训”工作,取得了一定的成效,干部教育培训工作总结。
现总结报告如下:一、基本情况全县共有干部**人,其中中共党员**人,大学本科以上学历**人,大专学历**人,中专学历**人,高中及以下学历**人。
**年,以县委党校、县行政学校为主阵地,举办各类培训**期,培训在职干部**人,占在职干部总数的**.*%,培训农村党员、干部**人,其中:举办科级领导干部轮训班*期,培训**人;举办科级领导干部“三个代表”重要思想专题学习班*期,培训**人;举办科级以下公务员培训班*期,培训**人;举办企业经营管理者培训班*期,培训**人;举办专业技术人员培训班*期,培训**人;举办非中共党员干部培训班*期,培训**人;举办理论骨干培训班*期,培训**人;举办妇女干部培训班*期,培训**人;举办基层团干培训班*期,培训**人;举办农村党支部书记、村主任培训班各*期,培训**人,达到了每年培训在职干部五分之一的要求,超额完成了培训任务。
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提取刚度矩阵和质量矩阵单元刚度矩阵和单元质量矩阵提取:直接用/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,filename,out,, ! 将输出信息送到filename.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中---------------------华丽的分割线---------------整体刚度矩阵和质量矩阵提取有两种方法:(1)二次开发,用自空间法进行模态分析,从.full 文件中提取,具体见/forum/viewthread.php?tid=10860&extra=&page=1(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,filename,1 !设置文件名称和刚度矩阵类型(刚度,质量,阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,filename,3 !列出整体刚度矩阵。
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进行有限元分析时,我们通常需要提取单元的刚度矩阵。
刚度矩阵是描述单元刚度的一个重要参数,它可以用来计算单元的应力、应变和位移等参数。
提取单元刚度矩阵的方法如下:
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 -。
ANSYS动力学分析中的矩阵
ANSYS动力学分析中的矩阵在ANSYS动力学分析中,矩阵是非常重要的工具。
矩阵作为数据结构的一种形式,可以方便地表示和处理线性方程组和向量的关系。
下面我们将介绍一些在ANSYS动力学分析中常用的矩阵。
1. 刚度矩阵(Stiffness Matrix):刚度矩阵表示了结构中各个节点之间的刚度关系。
它的行和列对应于结构的自由度,每个元素表示了对应节点之间的刚度。
在动力学分析中,刚度矩阵被用于描述结构在受到外力作用时的刚度响应。
2. 质量矩阵(Mass Matrix):质量矩阵表示了结构中各个节点之间的质量关系。
它的行和列对应于结构的自由度,每个元素表示了对应节点的质量。
在动力学分析中,质量矩阵被用于描述结构的惯性特性,以及结构在外力作用下的加速度响应。
3. 全局刚度矩阵(Global Stiffness Matrix):当结构由多个单元组成时,每个单元的刚度矩阵可以通过组装操作得到全局刚度矩阵。
全局刚度矩阵描述了整个结构的刚度关系,可以用于求解结构在受到外力作用时的位移响应。
4. 结构方程矩阵(Structural Equation Matrix):结构方程矩阵是动力学分析中的一种常用的矩阵形式,它将结构的质量矩阵、刚度矩阵和外力向量组合在一起,用于求解结构在受到外力作用时的位移响应。
结构方程矩阵可以通过求解线性方程组得到结构的位移响应。
5. 阻尼矩阵(Damping Matrix):阻尼矩阵描述了结构在受到外力作用时的阻尼特性。
阻尼矩阵在动力学分析中被用来模拟结构的阻尼效应,以及计算结构在受到外力作用时的振动响应。
6. 模态矩阵(Modal Matrix):模态矩阵描述了结构在自由振动状态下的振型特性。
模态矩阵包含了结构的特征向量,它可以用来计算结构的特征值和特征向量,以及结构在外力作用下的振动响应。
7. 自由度矩阵(Degree of Freedom Matrix):自由度矩阵用于表示结构的自由度信息。
- 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。