ANSYS beam188单元结果SMISC和NMISC结果的提取
- 格式:doc
- 大小:459.50 KB
- 文档页数:3
ansys提取质量和刚度矩阵
在ANSYS中,可以使用命令`*MATRIX, MASS, STIFF`来提取质量和刚度矩阵。
要提取质量矩阵,可以使用以下命令:
```
*MATRIX, MASS
/mesh_needed ! 控制只提取有限元网格中的质量矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后,ANSYS将会在指定的输出文件中生成质量矩阵。
要提取刚度矩阵,可以使用以下命令:
```
*MATRIX, STIFF
/mesh_needed ! 控制只提取有限元网格中的刚度矩阵
/output, CDH ! 控制输出文件的目录和文件名
/solve ! 解析模型
```
执行上述命令后,ANSYS将会在指定的输出文件中生成刚度矩阵。
请注意,以上命令需要在ANSYS命令窗口中执行,并且需要
已经建立合适的有限元模型。
您还需要根据需要修改`/output`命令以指定输出文件的目录和文件名。
ANSYS计算结果无难事,APDL经典命令让你的模型舞起来1、让你的ANSYS模型'舞'起来ANSYS计算结果的动画可采用ANTIME、ANMODE、ANCNTR、ANHARM等自动生成动画(AVI格式),使结果展示更加生动直观,相信使用ANSYS的都会制作。
然而,几何模型或有限元模型则无动画显示功能,有时为展示模型本身,会从多个角度截取图片。
那么,模型能否也可制作动画呢?答案是肯定的。
利用ANSYS的图形存储命令/SEG可以实现此功能,让你的模型动起来。
具体过程详见命令流中及其注释,动画上传总是失败,自己生成不要观看吧。
Finish$/clear$/prep7!简单的创建几何模型以减少篇幅blc4,0,0,4,2,5cyl4,2,4,1,,2,,4!关闭图例信息/plopts,info,off!以下开始制作模型动画!删除当前储存的图形/seg,dele/seg,multi,jhdh,1 !独立存储且不覆盖,文件名为jhdh/auto,1 !自动计算与图形区合适显示方式!正视/view,1,0,0,1$vplot!侧视/view,1,1$vplot!俯视/view,1,,1$vplot!D视图/view,1,1,1,1$vplot!循环36次,每次改变10度视角*do,i,1,36$/ang,1,10,ys,1$/replot$*enddo!关闭图形存储操作,保存为jhdh.avi文件/seg,off$/anfile,save,jhdh,avi其实比较简单,一旦进入模型动画制作过程,所有的xPLOT(x=KLAVNE)绘制的图形都将进入动画序列,按显示过程形成一部连续的动画。
2、用一个命令解决ANSYS数据列表分页早年初学ANSYS时,经常用到xLIST(如NLIST、ELIST、KLIST、LLIST、ALIST、VLIST等命令)和PRxSOL(如PRNSOL、PRESOL、PRRSOL、PRETAB、PRPATH)等列表命令,并希望将这些内容保存到TXT文件中,然后再导入EXCEL中处理。
BEAM188中文说明BEAM188 —3-D 线性有限应变梁(基于Ansys 5.61的help)MP ME ST PR PP ED元素描述BEAM188 适用于分析细长的梁。
元素是基于Timoshenko 梁理论的。
具有扭切变形效果。
BEAM188 是一个二节点的三维线性梁。
BEAM188 在每个节点上有6或7个自由度,(自由度)数目的变化是由KEYOPT(1)来控制的。
当 KEYOPT(1) = 0时 (默认), 每节点有6个自由度。
分别是沿x,y,z的位移及绕其的转动。
当 KEYOPT(1) = 1时,会添加第七个自由度 (翘曲量) 。
此元素能很好的应用于线性(分析),大偏转,大应力的非线性(分析)。
BEAM188包含应力刚度,在默认情况下,在某些分析中由NLGEOM来打开。
在进行弯曲( flexural),侧向弯曲( lateral), 和扭转稳定性( torsional stability)分析时,应力刚度应该是被打开的。
BEAM188 能够采用SECTYPE, SECDATA, SECOFFSET, SECWRITE,和SECREAD来定义任何截面(形状)。
. 弹性(elasticity),蠕变( creep),和塑性( plasticity)模型都是允许的 (不考虑次截面形状)。
图1. BEAM188 3-D 线性有限应变梁输入数据(元素的)几何形状,节点为止,即元素坐标系图示于BEAM188。
BEAM188在模型坐标系中是由节点 I 和节点 J 来定义的。
节点 K 是必需的元素方向点定义。
有关方向点的相关信息详见Generating a Beam Mesh With Orientation Nodes 在ANSYS Modeling and Meshing Guide中。
于LMESH和LATT命令说明中可见节点 K 的自动定义的详细说明。
在空间中这是一个没有量纲的元素。
截面形状是用SECTYPE和SECDATA命令(详见ANSYS Commands Reference )来独立定宓摹C恳桓鼋孛嫘巫淳囟ㄒ桓?ID 号(SECNUM)。
A NSYS梁单元如何正确提取应力这几个月来一直没花时间做论文,以至于现在手忙脚乱,而且发现很多基础的知识都忘光了,抑或是以前压根就没好好学过,我发现自己是个不折不扣的混日子的,只有在最后关头才会花时间做事情,而且常常是一阵忙,忙完就忘光~~元旦前赶回老家参加妹的婚礼,回来后一直不敢去见导师,因为导师之前要求我们在元旦之前把论文初稿给他,而我和同门都没给,原因很简单,就是没花时间在论文上,大部分时间都是混,除了上网、看电影外就是和同学一起玩~~总值及其堕落~~~这两天总算花了点时间做论文,当然主要还是为了能回去过年,这次我的口号是:“不做完论文就不回家!”然而很可能我过年前做不完,所以不得不预先面对可能不能回家过年的打算!好了,也不多说废话了,转入正题!这几天论文进展缓慢除了没用心外主要还是不知道到底要提取哪些数据,如何提取,模型是否完善,其实都很糊涂,今天拿第3工况大概地进行静力分析,但随后就是提取数据,然而基础不扎实,以至于不知道如何提取数据。
问题1:梁BEAM188的应力如何提取?最大、最小应力如何提取?就这个问题我上网找了一下,主要在上找,找了几个帖子,先贴过来看看ANSYS梁单元如何提取应力1、问:我用的188单元作谐响应分析,求解结束后,我想取出模型中的最大应力值作为参数,然后在接下来的优化当中用该最大应力作为状态变量,请问我应该怎么做啊,注意优化时,对应于每组参数值,最大应力点的位置都可能不同.请高手指点一下谢谢----------以下程序段分别得到目标变量(总体积),约束变量SV的最大应力值。
/POST1SET,NSORT,U,Y*GET,DMA X,SORT,,MA XETA BLE,VOLU,VOLUETA BLE,SMA X_I,NMISC,1ETA BLE,SMA X_J,NMISC,3ssum*GET,VOLUME,SSUM,,ITEM,VOLUESORT,ETA B,SMA X_I,,1 !按照单元SMA X_I的绝对值大小进行排序*GET,SMA XI,SORT,,MA XESORT,ETA B,SMA X_J,,1*GET,SMA XJ,SORT,,MA XSMA X=SMA XI>SMA XJ !约束变量SV:SMA X=最大应力值FINISH===============你这个程序段是针对beam3 吧,对beam188好像不行。
ansys仿真结果导出格式
ANSYS是一款强大的有限元分析软件,广泛应用于各种工程领域。
在完成仿真分析后,用户需要将结果导出为不同格式的文件,以便于后续的数据处理和分析。
本篇文章将详细介绍ANSYS仿真结果导出格式及操作方法。
首先,用户需要了解ANSYS仿真结果导出格式。
ANSYS支持多种导出格式,如CSV、XY、TEXT、EXCEL、HTML和MATLAB等。
不同的导出格式对应着不同的数据处理和分析需求,用户可以根据实际需求选择合适的导出格式。
接着,用户需要导出结果至文件。
在ANSYS软件中,用户可以通过以下步骤实现结果导出:
1.在ANSYS主界面,选择“File” -> “Export” -> “Results”。
2.在“Export Results”对话框中,选择“Format”为所需导出格式,如CSV格式。
3.设置导出路径和文件名,点击“Save”按钮。
在导出结果至文件的过程中,用户可以进行高级导出设置,如自定义数据筛选、数据时间步筛选和数据单位转换等。
这些高级设置可以帮助用户更好地满足特定需求。
在完成结果导出后,用户可以通过各种可视化工具对导出的数据进行分析,如在Excel中绘制图表、在Matlab中进行数据处理等。
这有助于用户更好地理解仿真结果,并为后续的设计优化提供有力支持。
总之,ANSYS仿真结果导出格式对于用户进行数据处理和分析具有重要意
义。
用户应熟悉各种导出格式,并根据实际需求选择合适的导出方式。
(完整版)ansysetable提取结果ANSYS etable提取结果etable提取运算结果ansys求解完结束后,会把计算结果保存在结果表中,行是单元的编号,列是计算信息,如节点位移、节点力、应力、应变等。
ETABLE这个命令就是把表中的需要的信息取出来,赋值给某个自定义的向量,再通过*GET命令可以指定某个具体的单元,就可以把该单元的对应计算结果提取出来了。
ETABLE, Lab, Item, Comp命令选项解释:Lab自定义的表名,用于后续命令或输出的标题,最多可使用8个字母,不可与预定义的表名称重复。
默认的表名是Item和Comp项的前四个字母组合而成的8个字母。
如果与用户之前定义的表名相同,本次结果将被包括在同一表中。
最多可定义200个不同的表名。
以下表名是ANSYS预定义的,不可用作用户自定表名:REFL, STAT, 和ERAS. Lab = REFL以ETABLE的最新选项重写所有ETABLE命令预定义的表,但保留字段将被忽略,这个命令在载荷步改变后重写表时很方便。
Lab = STAT将显示储存的表的值。
Lab = ERAS将删除整个表。
Item选项名称。
常用的选项名称见后表。
某些选项需要栏目名。
Item = ERAS将删除表中的某一栏。
Comp选项的栏目名(如果需要的话)。
常用的栏目名见后表。
说明:定义单元值的表以便后续处理。
单元值表可以被认为是工作表,其行代表所有被选择的单元。
其列代表通过ETABLE命令输入表中的单元值。
每一列数据有一个用户定义的名称,用于列表和显示。
将数据输入单元表后,你不仅可以列出和显示你的数据,还可以对数据进行许多操作,例如列相加或列相乘[SADD, SMULT],为安全计算定义允许的应力[SALLOW],或者将一列数据和另一列相乘[SMULT]。
更多的细节请看ANSYS Basic Analysis Guide。
有很多不同类型的数据结果可以被存在单元表中。
beam188单元自定义截面的方法及预览实现第一部分定义cross section 文件:1,2种方法,一种可以在autocad种画好图形,存为".sat"格式,然后再通过ansys的i mport端口导入,一种是直接在ansys里画好。
2,对导入或画好的图形进行网格划分,用plane82网格划分好后(一定完成网格的划分,注意不是预划分),选择precessor-section-beam:首先选择write from area,定义为一个.sect文件3,然后选择read section mesh 定义section name 和section library file.选择完成。
4,绘制beam188单元,完成后,定义单元横截面属性:meshing-mesh attribute-default attributes 在section number处,选择刚才定义的横截面文件(如果提示没有定义的横截面文件,选择read section mesh 定义section name 和section library file.再重新读取一次,再meshing-mesh attribute-default attributes 在section number处,选择刚才定义的横截面文件)。
到现在,beam188单元的自定义截面就定义成功了。
第二部分下面就是怎么实现三维预览的问题了,必须要先给beam188单元划分完网格之后,才能实现三维预览,实际上,从现实三维图形的命令中,我们也能够知道,所谓的三维预览,实际上只是针对element才起作用的,element肯定是在划分完网格后才有的有限元单元。
因此,三维预览的方法:1,给beam188单元划分网格(一定完成网格的划分,注意不是预划分)。
2,选择plotctrl-style-size and shape,在displace of element 后的方格中打上勾。
Beam1883 维线性有限应变梁单元Beam188 单元描述Beam188 单元适合于分析从细长到中等粗短的梁结构,该单元基于铁木辛哥梁结构理论,并考虑了剪切变形的影响。
Beam188 是三维线性(2 节点)或者二次梁单元。
每个节点有六个或者七个自由度,自由度的个数取决于KEYOPT(1)的值。
当KEYOPT(1)=0(缺省)时,每个节点有六个自由度;节点坐标系的x、y、z 方向的平动和绕x、y、z 轴的转动。
当KEYOPT(1)=1 时,每个节点有七个自由度,这时引入了第七个自由度(横截面的翘曲)。
这个单元非常适合线性、大角度转动和/并非线性大应变问题。
当NLGEOM 打开的时候,beam188 的应力刚化,在任何分析中都是缺省项。
应力强化选项使本单元能分析弯曲、横向及扭转稳定问题(用弧长法)分析特征值屈曲和塌陷)。
Beam188/beam189 可以采用sectype、secdata、secoffset、secwrite 及secread 定义横截面。
本单元支持弹性、蠕变及素性模型(不考虑横截面子模型)。
这种单元类型的截面可以是不同材料组成的组和截面。
Beam188 从6.0 版本开始忽略任何实参数,参考seccontrols 命令来定义横向剪切刚度和附加质量。
单元坐标系统(/psymb,esys)与beam188 单元无关。
下图是单元几何示意图:BEAM188 输入数据该单元的几何形状、节点位置、坐标体系如图“BEAM Geometry”所示,beam188 由整体坐标系的节点i 和j 定义。
节点K 是定义单元方向的所选方式,有关方向节点和梁的网格划分的信息可以参见ANSYS Modeling and Meshing Guide中的Generating a Beam Mesh With Orientation Nodes。
参考lmesh 和latt 命令描述可以得到k 节点自动生成的详细资料。
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--—输出矩阵的路径和文件名,缺省为当前工作路径和当前工作文件名。
BEAM188BEAM188 — 3-D 线性有限应变梁(基于Ansys 5.61的help)MP ME ST PR PP ED元素描述BEAM188 适用于分析细长的梁。
元素是基于Timoshenko 梁理论的。
具有扭切变形效果。
BEAM188 是一个二节点的三维线性梁。
BEAM188 在每个节点上有6或7个自由度,(自由度)数目的变化是由KEYOPT(1)来控制的。
当 KEYOPT(1) = 0时 (默认), 每节点有6个自由度。
分别是沿x,y,z的位移及绕其的转动。
当 KEYOPT(1) = 1时,会添加第七个自由度 (翘曲量) 。
此元素能很好的应用于线性(分析),大偏转,大应力的非线性(分析)。
BEAM188包含应力刚度,在默认情况下,在某些分析中由NLGEOM来打开。
在进行弯曲( flexural),侧向弯曲( lateral), 和扭转稳定性( torsional stability)分析时,应力刚度应该是被打开的。
BEAM188 能够采用SECTYPE, SECDATA, SECOFFSET, SECWRITE,和SECREAD来定义任何截面(形状)。
. 弹性(elasticity),蠕变( creep),和塑性( plasticity)模型都是允许的 (不考虑次截面形状)。
图1. BEAM188 3-D 线性有限应变梁输入数据(元素的)几何形状,节点为止,即元素坐标系图示于BEAM188。
BEAM188在模型坐标系中是由节点 I 和节点 J 来定义的。
节点 K 是必需的元素方向点定义。
有关方向点的相关信息详见Generating a Beam Mesh With Orientation Nodes 在ANSYS Modeling and Meshing Guide中。
于LMESH和LATT命令说明中可见节点 K 的自动定义的详细说明。
在空间中这是一个没有量纲的元素。
截面形状是用SECTYPE和SECDATA命令(详见ANSYS Commands Reference )来独立定宓摹C恳桓鼋孛嫘巫淳 囟ㄒ桓?ID 号(SECNUM)。
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中整体、单元刚度和质量矩阵的提取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中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。
1.1、用户程序法,需要二次开发该功能需要进行二次开发,由ansys形成的二进制文件.full提取整体刚度和质量矩阵。
基于ansys的一个用户开发程序例子编了一个程序(附件中). 开发环境:compaqfortran6。
5运行环境:win2000。
一个主文件self。
for;另一个文件matrixout。
f90用于矩阵输出;binlib.lib为ansys提供的库文件,将其引入项目中(也可直接扔进debug目录);.full文件由子空间迭代模态分析获得。
运行编译后的可执行文件.exe其他文件见http://d。
/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中单元刚度矩阵、总刚矩阵的提取一、单元刚度矩阵的提取/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、整体刚度和质量矩阵的提取。
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:二进制格式。