ansys提取最大应力应变命令

CYCLES USED/ALLOWED = 0.1000E+05/ 8455（实际循环数/许用循环数）= PARTIAL USAG E（局部损伤)=1.18267

PRODUCE ALTERNATING SI (SALT) = 85.508 (应力幅值)

使用系数＝要求的循环周期数/允许的循环周期数

用ANSYS的APDL方式直接获得最大应力和最大应变的节点编号及其数值的有好多...

用ANSYS的APDL方式直接获得最大应力和最大应变的节点编号及其数值的有好多种方法：

如下面的这种，先对节点的值进行SORT，在提取最大的值即可。

最大应力节点编号及其数值的提取：

ALLSEL

NSORT,S,EQV,0,0,ALL

*GET,MAX_SEQV,SORT,0,IMAX

*GET,MAXSEQV,NODE,MAX_EQV,S,EQV

依次类推，可以最大应变节点编号及其值数值的提取：

ALLSEL

NSORT,EPTO,EQV,0,0,ALL

*GET,MAX_EPTOEQV,SORT,0,IMAX

*GET,MAXEPTOEQV,NODE,MAX_EPTOEQV,S,EQV

最大位移节点编号及其数值的提取：

ALLSEL

NSORT,U,SUM,0,0,ALL

*GET,MAX_U,SORT,0,IMAX

*GET,MAXU,NODE,MAX_U,U,SUM

还有其他方法：如把所有的节点应力应变数值提取然后找最大值

*GET,NUMALL,NODE,0,COUNT

*GET,NUMSTART,NODE,0,NUM,MIN

SMAX=0

INUSE=NUMSTART

*DO,I,1,NUMALL,1

*GET,SSUM,NODE,INUSE,S,EQV

*IF,SMAX,GE,SSUM,THEN

SMAX=SMAX

*ELSEIF,SMAX,LT,SSUM,THEN

SMAX=SSUM

*ENDIF

INUSE=NDNEXT(INUSE)

*ENDDO

还有一种是首先在窗口上显示应力应变云图，然后利用GET命令得到最大值。如提取最大应力为：

PLNSOL,S,EQV,0

*GET, PAR, PLNSOL, 0, MAX

或将所有的节点定义成一个路径，*GET, PAR, PATH, 0, MAX，

另：获取结果各个时间点的最大应力值

*DIM,MAXSEQVNODE,TABLE,63,1,1,,, !定义数组

*DO,ISET1,1,63,1 !根据结果存储的数目而定,如63

/POST1

SET,1,,1,,,,ISET1,

PLNSOL,S,EQV,0

*GET,MAXSEQV_NODEI,PLNSOL,0,MAX !获取结果,这里是最大MISES等效应力

*SET,MAXSEQV_NODE(ISET-1,0,1),(ISET-1)*4E-4 !确定对应的时间

*SET,MAXSEQV_NODE(ISET-1,1,1),MAXSEQV_NODEI !写入最大等效应力值

*ENDDO

。。。。。。。。。。。。。

ansys 怎样提取某个节点的应力和应变(网络软件)

ansys 怎样提取某个节点的应力和应变？ 时间:2010-11-14 来源:网络浏览次数: 次 1. 最简单的办法是使用NSORT，打印出结果，可以通过控制使其输出到文件 2. 使用apdl能复杂一点，下面是以前经常用的一段命令流，参考着修改一下吧*CREATE,GET_node_inf,mac, *GET,Nnod,NODE,0,COUNT !获取所选择的节点总数 *DIM,S_Xyz,ARRAY,NNOD,5 !定义1个数组存放数据 *GET,Nd,NODE,0,NUM,MIN !获取最小的节点编号 *DO,I,1,Nnod,1 S_Xyz(I,1)=Nd !将节点列表放数组第1列 S_Xyz(I,2)=NX(Nd) !节点的X坐标放数组第2列 S_Xyz(I,3)=NY(Nd) !节点的Y坐标放数组第3列 S_Xyz(I,4)=NZ(Nd) !节点的Z坐标放数组第4列 !*GET,S_Xyz(I,5),NODE,ND,S,EQV !节点的von mises值放数组第5列 *GET,S_Xyz(I,5),NODE,ND,U,SUM !节点的总变形值值放数组第5列 Nd=NDNEXT(Nd) !读出下一个节点编号 *ENDDO *END *CREATE,OUT_node_inf,mac, *CFOPEN,node_info,txt,, *VWRITE,S_Xyz(1,1),S_Xyz(1,2),S_Xyz(1,3),S_Xyz(1,4),S_Xyz(1,5) (F10.0,3F15.4,E15.5) *CFCLOS *END GET_node_inf OUT_node_inf /delete,GET_node_inf,mac /delete,OUT_node_inf,mac 另附 1.先对节点的值进行SORT，在提取最大的值即可。 最大应力节点编号及其数值的提取：

ansys质量矩阵刚度矩阵提取

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.full A simple way to dump the matrix is as follow: ------------------- /aux2 fileaux2,file,emat form,long dump,all ------------------- 2. 可以使用/DEBUG命令来得到。详细步骤参见下面的宏文件 finish /clear PI=3.1415926 w1=3 w2=10 w3=6 w4=1.2 r=.8 t=0.08 /PREP7 !* ET,1,SHELL63 R,1,t ET,2,MASS21 R,2,500,500,500,2000,2000,2000,

!* UIMP,1,EX, , ,2e11 UIMP,1,NUXY, , ,0.3, UIMP,1,DAMP, , ,0.2, UIMP,1,DENS, , ,7800, BLC4,0,0,w2,w1 ESIZE,1.5,0, AMESH,all NSEL,S,LOC,X,0.0 D,all, , , , , ,ALL, , , , , allsel,all SFA,all,1,PRES,12 FINISH /OUTPUT,cp,out,, ! 将输出信息送到cp.out文件 /debug,-1,,,1 ! 指定输出单元矩阵 /SOLU SOLVE finish /OUTPUT, TERM ! 将输出信息送到output windows中 ! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵 3. 其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。 面是一个小例题，自可明白。 /prep7 k,1 k,2,3000 l,1,2 et,1,beam3 mp,ex,1,2e5 mp,prxy,1,0.3 mp,dens,1,2e3 r,1,5000,2e7,200 lesize,all,,,10 lmesh,all

ANSYS中杆梁壳单元内力的输出

杆梁壳单元内力的输出 计算完成后，就可以对杆梁壳单元进行后处理，显示位移、应力（需打开单元开关显示开关）的操作和普通的单元的后处理是一样的。但其内力的输出和显示需要进行特殊的处理。首先要保证在定义单元类型时，打开了内力输出的选项；其次，在完成计算后需要通过定义单元表的方法来提取。 一、定义输出选项 Beam4单元说明的input summary中可以看到关于内力输出的开关选项说明： 二、定义单元表 查看对应项和序号：在单元说的后的单元输出表中查到MFOR(X,Y,Z)和MMOM(X,Y,Z)对应的项和序列号如下表，注意I和J节点上对应的序列号不同。 GUI方式如下图： 命令流过程如下： （I节点处结果） Etable,forx_i,smisc,1 ！x方向剪力（即轴力） Etable,fory_i,smisc,2 ！y方向剪力 Etable,forz_i,smisc,3 ！z方向剪力 Etable,momx_i,smisc,4 ！x方向弯矩 Etable,momy_i,smisc,5 ！y方向弯矩 Etable,momz-i,smisc,6 ！z方向弯矩

（J节点处结果） Etable,forx_j,smisc,7 ！x方向剪力（即轴力） Etable,fory_j,smisc,8 ！y方向剪力 Etable,forz_j,smisc,9 ！z方向剪力 Etable,momx_j,smisc,10 ！x方向弯矩 Etable,momy_j,smisc,11 ！y方向弯矩 Etable,momz-j,smisc,12 ！z方向弯矩 三、绘制内力图 对线单元和面单元，其绘制单元表的结果的命令是不一样的。 线单元的内力云图的绘制为：PLLS. 面单元的内力云图的绘制为：PLETAB 四、注意事项 1.弯矩图的绘制与显示比例 在绘制弯矩图时（其他内力图也一样），可以给定一个显示比例系数Fact（默认等于1），当为负值时，反向显示，数值为放大系数。同时还须注意的是，弯矩图的显示还与单元坐标系的方向有关，因此，如果同一段梁中各个单元的单元坐标系不一致（可以通过改变线的方向来改变线单元的单元坐标系），还将导致弯矩图显示不连续。壳单元的后处理中更加要注意单元坐标系的方向，而且如果面上的壳单元划分不整齐，还将导致提取出来的内力图方向不统一，没有参考价值。 2.划分数的影响 要特别注意梁单元划分数对弯矩图的影响，除了因为有限元法是一种把连续结构体离散化的近似求解方法，还因为弯矩图是以梁单元的I和J节点结果进行连线绘制的。 可通过例子中用10份，2份，100份进行对比。

Ansys中节点力提取

Ansys中节点力提取几个问题的说明 对于ansys中节点力提取的命令，一般有如下命令可以用， *GET,Par, NODE, N, RF,FX(FY/FZ/MX/MY/MZ) 这组命令是我们最开始用的，用来提取节点反力，但是有个缺陷，节点反力只在有约束位置才能提取，如果在结构中任何一个节点处提取此节点所受合力，界面操作有两种方法。 Main Menu>General Postproc>Nodal Calcs>Total Force Sum Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node 但是执行上面两个操作有个前提，需要选出对应的单元和节点，下面举例说明： 如下图：800臂架结构 由于要对连接架+塔帽进行单独的详细分析，需要提取旋转架与塔帽连接处铰点对塔帽的作用力。而且为了在详细模型中施加载荷的时候方便，提取结果的坐标系需要是X向沿着主臂的局部坐标系，见示图1。 运用Main Menu>General Postproc>Nodal Calcs>Total Force Sum 或者 Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node这两个操作可以实现。 下面就这两个操作的的结果进行比较。就如上模型，研究塔帽和旋转架连接左侧铰耳处受力提取，见示图2。 将塔帽上与此铰点相连接的两个单元选出来，选择此节点，见示图3。 读取结果文件，设置结果坐标系为要求的局部坐标系（文件中为局部坐标系11）。

执行Main Menu>General Postproc>Nodal Calcs>Total Force Sum，选取Active Rsys，结果界面如下： 执行Main Menu>General Postproc>Nodal Calcs>Sum @ Each Node，结果界面如下： 从结果界面上可以看出，了；两个结果界面中显示的结果来看，“结果显示1”与“结果显示3”的结果数据相同，而“结果显示2”与其他两个在在力的显示数据上相同，而矩的显示数据上并不相同。之所以出现这种差别，原因在于在求矩的作用数据时，矩的作用中心不相同。“结果显示2”中所显示的数据其矩的作用中心为提取载荷点（1600078），而“结果数据1”和“结果数据3”中矩的作用中心为结果坐标系的原点。 要想使上述3个结果数据显示值相同，只需要将结果显示中矩的作用中心设置到提取点上就可以了。 通过Main Menu>General Postproc>Nodal Calcs>Summation Pt>At Node将矩的作用点设置到1600078，再执行命令Main Menu>General Postproc>Nodal Calcs>Total Force Sum，结果

梁结构应力分布ANSYS分析汇总

J I A N G S U U N I V E R S I T Y 先进制造及模具设计制造实验 梁结构应力分布ANSYS分析 学院名称：机械工程学院 专业班级：研1402 学生姓名：XX 学生学号：S1403062 2015年5 月

梁结构应力分布ANSYS分析 （XX，S1403062，江苏大学） 摘要：本文比较典型地介绍了如何用有限元分析工具分析梁结构受到静力时的应力的分布状态。我们遵循对梁结构进行有限元分析的方法，建立了一个完整的有限元分析过程。首先是建立梁结构模型，然后进行网格划分，接着进行约束和加载，最后计算得出结论，输出各种图像供设计时参考。通过本论文，我们对有限元法在现代工程结构设计中的作用、使用方法有个初步的认识。 关键词：梁结构；应力状态；有限元分析；梁结构模型。 Beam structure stress distribution of ANSYS analysis （Dingrui, S1403062, Jiangsu university） Abstract: This article is typically introduced how to use the finite element analysis tool to analyze the stress of beam structure under static state distribution. We follow the beam structure finite element analysis method, established the finite element analysis of a complete process. Is good beam structure model is established first, and then to carry on the grid, then for constraint and load, calculated the final conclusion, the output of images for design reference. In this article, we have the role of the finite element method in modern engineering structural design, use method has a preliminary understanding. Key words: beam structure; Stress state; The finite element analysis; Beam structure model. 1引言 在现代机械工程设计中，梁是运用得比较多的一种结构。梁结构简单，当是受到复杂外力、力矩作用时，可以手动计算应力情况。手动计算虽然方法简单，但计算量大，不容易保证准确性。相比而言，有限元分析方法借助计算机，计算精度高，

提取刚度矩阵

============为什么折腾这个文档======== 我有一个计算线性动力学方程组的瞬态、谐响应和静力学的python程序，现希望开发一个将ANSYS组集好的总体矩阵导入该PYTHON程序中的接口。 该问题可分解为： [STEP1] [ANSYS]->[包含矩阵信息的文件] [STEP2] [包含矩阵信息的文件]->[python通用数据对象] [STEP3] [python通用数据对象]->[程序特定数据对象]->[进行计算] 因此检索了一些帖子，基本上完成了这项工作，本文是对[STEP1]和[STEP2]的整理，并且利用[STEP3]对结果进行了验证 ============主要内容================== 1，了解从ANSYS中提取总体矩阵和载荷向量的方法； 2，了解提取出来的矩阵是怎样表示的； 3，说明在Python中，如何读取这样的矩阵； 4，构造一个简单的算例，说明整个【建模】-【提取】-【读取】过程及其正确性； =========站内检索综述==================== 检索词：提取矩阵 得到21个结果，代表性的帖子有下面这9个： 编号[1] 标题：ansys中怎样提取质量，刚度，阻尼矩阵？ 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：pengweicai给出了一段网上最常见的提取代码,该程序以fortran 写成，可以利用.full文件以及一些列约定将ANSYS中的总体矩阵读入FORTRAN中。 编号[2] 标题：如何得知HBMAT命令提取的质量、刚度矩阵对应的自由度？ 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：提出了使用HBMAT命令提取稀疏矩阵时常见的问题：我们如何知道提取出来的信息是怎么储存的呢？ 编号[3] 标题：[分享]ANSYS中整体、单元刚度和质量矩阵的提取 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：在该帖子的7楼，其实已经给出了帖子[2]中问题的解答，即HBMAT 中提取出来的矩阵是Harwell-Boeing格式的，并且给出了该格式的细节，可惜是英文的，没引起多少关注。 编号[4] 标题：帮我看看提取的刚度与质量矩阵 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：这个帖子所示的矩阵并非是使用HBMAT命令提出出来的，而应该是SELIST命令列举出来的未压缩的矩阵，后续楼层的回帖给了大家一个提示，即有可能提取出来的矩阵是引入了边界条件的（即删除了被约束的行和列的）。 编号[5] 标题：提取刚度矩阵的问题 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：本帖作者的工作是基于单元刚度矩阵的，因此ANSYS中提取的单元刚度矩阵是否处于总体坐标系就成为问题。该问题并非本文内容，但仍值得关注。 编号[6] 标题：提取刚度矩阵丢失节点的问题 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：帖子[5]作者的又一帖，在这里帖子[5]的问题得到了欧阳中华老师的回答。 编号[7] 标题：提取刚度矩阵的ANSYS操作过程 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：实际上这就是使用HBMAT从ANSYS中提取总体矩阵的全过程！只是还有一些细节待确定。 编号[8] 标题：提取整体刚度矩阵、质量矩阵及阻尼矩阵的简单方法 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：给出了利用“不减缩的”子结构方法来得到总体矩阵的方法（这也是网络上常见的代码之一） 编号[9] 标题：质量矩阵、刚度矩阵如何提取？ 地址：https://www.wendangku.net/doc/e114516268.html,/forum-vi ... fromuid-159019.html 要点：16443在5楼的回帖中给出了提取刚度矩阵的三种方法 =======站外检索略述======================== 百度检索：提取矩阵 比较好的帖子有： 编号[10] 来源:百度文库 标题:怎样从ansys中提取单元刚度矩阵与质量矩阵 地址:https://www.wendangku.net/doc/e114516268.html,/view/3cf5e567f5335a8102d220d9.html 要点:这应该就是16443在帖子[9]中回复的内容了，全面的总结了在帖子[3,4,5,9]中涉及的问题。 编号[11] 来源:中华钢结构标题:ansys刚度矩阵Harwell-Boeing格式的具体含义讨论 地址:https://www.wendangku.net/doc/e114516268.html,/forum/viewthread.php?tid=184007 要点:如题，后续楼层给出了一些将矩阵读入ANSYS的APDL（好不容易读出来，又读进去干嘛呢……） 编号[12] 来源:simwe 标题:关于ANSYS（质量、刚度、阻尼）矩阵Harwell-boeing格式数据的说明 地址:https://www.wendangku.net/doc/e114516268.html,/archiver/tid-924778.html 要点:比[11]更透彻的HB格式说明! ============================================================= =======1.从ANSYS中提取总体矩阵的方法================================= ============================================================= 1,用/DEBUG命令 2,子结构法

Ansys梁分析实例

工程介绍： 某露天大型玻璃平面舞台的钢结构如图1所示，每个分格（图2中每个最小的矩形即为一个分格）x方向尺寸为1m，y方向尺寸为1m；分格的列数（x向分格）=8，分格的行数（y向分格）=5。 钢结构的主梁（图1中黄色标记单元）为高140宽120厚14的方钢管，其空间摆放形式如图3所示；次梁（图1中紫色标记单元）为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间（如不是正处于X方向正中间，偏X坐标小处布置）的次梁的两端，如图2中标记为 U R处。主梁和次梁之间是固接的。 xyz xyz 玻璃采用四点支撑与钢结构连接（采用四点支撑表明垂直作用于玻璃平面的面载荷将传递作用于玻璃所在钢结构分格四周的节点处，表现为点载荷；试对在垂直于玻璃平面方向的42 KN m的面载荷（包括玻璃自重、钢结构自重、活载 / 荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析。（每分格面载荷对于每一支撑点的载荷可等效于1KN的点载荷）。 作业提交的内容至少应包括下面几项： （1）屏幕截图显示该结构的平面布置结构，图形中应反映所使用软件的部分界面，如图2； （2）该结构每个支座的支座反力； （3）该结构节点的最大位移及其所在位置； （4）对该结构中最危险单元（杆件）进行强度校核。 图1

图2 图3 本操作中选用的单位为：（N，mm，MPa）。具体操作及分析求解： 1．更该工作文件和标题。如图1.1-1.5所示

图1.1 图1.2

图1.3 图1.4 图1.5

图1.6 2.选择单元类型。 根据题目要求，选择单元类型为beam-3D-2node-188单元。 执行Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add ，选择beam-3D-2node-188。如图2.1所示。 图2.1 3．定义材料属性 该钢结构材料为碳素结构钢Q235，则将弹性模量设置为200GPa，泊松比设置为0.3。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2.05e，在PRXY框中输入0.3。操作步骤为如图3.1；3.2所示。

ansys etable提取结果

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。 有很多不同类型的数据结果可以被存在单元表中。例如，许多单元的选项只有一个值（也就是说，每一个单元对应此选项只有一个值）单一值的选项包括：SERR, SDSG, TERR, TDSG, SENE, TENE, KENE, JHEAT, JS, VOLU和CENT. 其余的选项是多个值的（也就是说，这些值在单元中是变化的，每个节点有不同的值）。因为每个单元只能有一个值存在单元表中，多值的选项存入的是平均值（视节点数而定）。例外的是FMAG和所有的单元力选项，它们存入的是相关节点值的和。（这段话的意思是说，单值的单元选项，如单元体积，存入表中的就是这个值；而在单元不同位置有不同值的选项，如应力？，写入表中的是单元的平均值。根本原因在于一个单元只能对应表中的一个数据。） ETABLE命令中可以使用两种数据访问方法，视你想储存的数据不同而不同。一些结果只用通用名就可以访问（要素名法），而另一些结果需要一个标志名和标志数（序列数法）。要素名法用于访问常用的单元数据（也就是说，绝大部分单元类型都有的数据）。所有的单值选项以及一部分多值选项可以用要素名法访问。不同的单元值视计算方法不同和选择集不同而不同。(AVPRIN, RSYS, LAYER, SHELL,和ESEL) 尽管节点值不用单元表也可以很容易地列出和显示，你仍然可能需要利用单元表储存这些节

提取单元刚度矩阵

单元刚度矩阵的提取 刚度矩阵在有限元求解过程中扮演者非常重要的角色，以最小位能原理求解过程为例最终越是转换为含有结构刚度矩阵的能量泛函的取值问题。有限元过程中涉及到三类刚度：单元刚度矩阵，组合结构刚度矩阵和最终求解刚度矩阵。 其中单元刚度矩阵：仅与单元的自身自由度有关，同一编号的单元矩阵的维数是固定。组合结构刚度：矩阵根据求解的初始变量个数决定刚度矩阵的维数，属于单元组装后的初始刚度，维数和整个单元初始变量个数相等。最终求解刚度矩阵：代入边界条件简化后的刚度。以《Finite Element Analysis-Theory and Application With ANSYS》中的梁单元例子为例，解释刚度提取过程： 此模型的单元刚度矩阵：（学则beam3梁单元后，该单元包含两节点，每个节点具有三个自由度，因此对应单元刚度矩阵为6*6的方阵）

组合结构刚度矩阵：（该结构含有三个节点，每个节点具有三个原始自由度，因此组合结构刚度矩阵具有9*9阶的形式） 最终求解刚度矩阵：（由于边界条件的存在，该结构中，1,3点的自由度不存在，求解参数中有六个参数已知，因此对最终求解刚度矩阵为三阶方阵） 通过最终的刚度矩阵组成的方程，求解出2节点的位移解，再以这些原始解得出应力，应变，支反力的其他的解。 ansys实现过程： 提取思路如下：通过/debug提取单元刚度矩阵，通过filname.full文件提取后两者的矩阵 ansys实现过程如下： finish /clear /filname,k,1 /prep7 N,1 N,2,120 N,3,120,-108 et,1,beam3 mp,ex,1,3.0e7 mp,prxy,1,0.3 R,1,7.65,204,10 E,1,2 E,2,3 /debug,-1,,,1,,,,,

ANSYS提取单元内力的方法

ANSYS提取单元或节点内力的方法 方法1：节点荷载（List Results→Nodal Loads） 方法2：节点合力计算（Nodal Cals→Sum @ Each Node） 方法3：单元解中的节点解（List Results→Element Solution→Structural Forces & Moments）方法4：支座反力（List Results→Reaction Solu） 方法5：单元表（List Results→Elem Table Data） 上述各方法提取的结果关系如下： （1）方法1和方法2提取的结果完全相同，但结果为0的项在方法1的结果列表中不显示，而方法2的结果列表则会全部显示。 （2）方法3提取的结果是每个单元各节点在该单元中的内力，针对同一节点，将其在各个单元中的内力求和，其累加结果与方法1和2得到的结果一致。 （3）方法4提取的结果只显示有施加位移约束的节点反力，其数值大小与方法1和2得到的结果相差一个正负号，即节点内力和节点反力刚好是一对作用力与反 作用力。 （4）方法5提取的结果是单元的内力，如果单元的形函数为线性（如BEAM188单元设置“KEYOPT(3)=0”），则ANSYS会取单元中点作为积分点并将其数值代 替单元内的线性变化，因此其输出结果的绝对值等于方法3中对应单元的各节 点相应内力绝对值的平均值；如果单元的形函数为非线性（如BEAM188单元 设置“KEYOPT(3)=2”），则单元各节点的内力不同，其结果与方法3得到的结 果一致。 （5）方法1~4提取的结果都是默认基于整体坐标系的，而方法5提取的结果是基于单元坐标系的，因此提取结果的方向和正负号需特别注意。有限元中力的方向 和结构力学中的方向是有区别的，不论是什么结果坐标系，力的正方向取为对 应结果坐标的正方向，弯矩则是对应坐标轴的顺时针为正。

ansys怎样提取某个节点的应力和应变

ansys 怎样提取某个节点的应力和应变 时间:2010-11-14 来源:网络浏览次数: 次 1.最简单的办法是使用NSORT，打印出结果，可以通过控制使其输出到文件 2. 使用apdl能复杂一点，下面是以前经常用的一段命令流，参考着修改一下吧*CREATE,GET_node_inf,mac, *GET,Nnod,NODE,0,COUNT !获取所选择的节点总数 *DIM,S_Xyz,ARRAY,NNOD,5 !定义1个数组存放数据 *GET,Nd,NODE,0,NUM,MIN !获取最小的节点编号 *DO,I,1,Nnod,1 S_Xyz(I,1)=Nd !将节点列表放数组第1列S_Xyz(I,2)=NX(Nd) !节点的X坐标放数组第2列S_Xyz(I,3)=NY(Nd) !节点的Y坐标放数组第3列S_Xyz(I,4)=NZ(Nd) !节点的Z坐标放数组第4 列 !*GET,S_Xyz(I,5),NODE,ND,S,EQV !节点的 von mises值放数组第5列 *GET,S_Xyz(I,5),NODE,ND,U,SUM !节点的总变形值值放数组第5列 Nd=NDNEXT(Nd) !读出下一个节点编号 *ENDDO *END *CREATE,OUT_node_inf,mac, *CFOPEN,node_info,txt,, *VWRITE,S_Xyz(1,1),S_Xyz(1,2),S_Xyz(1,3),S_Xyz(1,4),S_Xyz(1,5) ,, *CFCLOS *END GET_node_inf OUT_node_inf /delete,GET_node_inf,mac /delete,OUT_node_inf,mac 另附 1.先对节点的值进行SORT，在提取最大的值即可。 最大应力节点编号及其数值的提取：

ANSYS_MPC连接梁壳单元实例

2010-05-21 22:12:04 作者：zhz2004 来源：机械CADl论坛浏览次数：621 网友评论0 条 近日在论坛看到些用ansys的坛友问及beam单元和shell单元、beam单元和solid单元、shell单元和solid单元的连接问题。其实解决此类问题的方法不只一种，耦合约束方程、绑定接触都是有效的方法。其中耦合约束方程适用于小变形，而绑定接触即可用于小变形，也可用于大变形的几何非线性分析。下面，我将本人所做的用MPC方法连接shell单元和b eam单元的详细步骤提供给大家，与各位共勉。 添加shell单元(略) 添加beam单元(略) 添加shell实常数 添加shell实常数：shell厚度 添加beam截面：圆钢 内经、外径及网格密度

预览网格 开始建模：转动工作平面工作平面z轴向上 建立圆面 继续： 将面拉伸成体 定义拉伸高度：0.5m 删除体，留面 显示面 删除空圆柱的顶面和底面

创建点：用于建立梁单元的第一个点。两点之间创建（正中）。 复制点：用于建立梁单元的第二个点。复制：Y方向0.5m 连接两点，用于创建梁单元。 继续 定义材料属性，有点晚^_^ 准备划分壳单元 划分壳，映射方法

准备划分梁单元 划分梁单元 选中要划分梁单元的线 完成，定义mpc接触 ＧＵＩ:ＭａｉｎＭｅｎｕ→Ｐｒｅ-ｐｒｏｃｅｓｓｏｒ→Ｍｏｄｅｌｉｎｇ→Ｃｒｅａｔｅ→ＣｏｎｔａｃｔＰａｉｒ,进入接触向导,然后按照提示与帮助说明进行选择目标面接触面等操作[4]。在创建接触对前,单击Ｏｐｔｉｏｎａｌｓｅｔｔｉｎｇ按钮弹出Ｃｏｔａｃｔｐｒｏｐｅｒｔｉｅｓ对话框,将Ｂａｓｉｃ选项卡中的Ｃｏｎｔａｃｔａｌｇｏｒｉｔｈｍ即接触算法设置为ＭＰＣａｌｇｏｒｉｔｈｍ即可。操作完成后,ＡＮＳＹＳ自动定义目标单元与接触单元类型,并生成接触对。 定义主控点 选择梁单元的下面一个关键点（当然也可以选择梁单元的最下一个node，相应选项要选pick existing node...） 选择梁单元的下面一个关键点 继续下一个： 施加集中力x方向10000n 计算结果，位移云图 显示梁截面的位移云图 显示梁单元形状

ansys动力学分析全套讲解

第一章模态分析 §模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS产品家族中的模态分析是一个线性分析。任何非线性特性，如塑性和接触（间隙）单元，即使定义了也将被忽略。ANSYS提供了七种模态提取方法，它们分别是子空间法、分块Lanczos法、PowerDynamics法、缩减法、非对称法、阻尼法和QR阻尼法。阻尼法和QR阻尼法允许在结构中存在阻尼。后面将详细介绍模态提取方法。 §模态分析中用到的命令 模态分析使用所有其它分析类型相同的命令来建模和进行分析。同样，无论进行何种类型的分析，均可从用户图形界面（GUI）上选择等效于命令的菜单选项来建模和求解问题。 后面的“模态分析实例（命令流或批处理方式）”将给出进行该实例模态分析时要输入的命令（手工或以批处理方式运行ANSYS时）。而“模态分析实例（GUI方式）” 则给出了以从ANSYS GUI中选择菜单选项方式进行同一实例分析的步骤。（要想了解如何使用命令和GUI选项建模，请参阅<>）。<>中有更详细的按字母顺序列出的ANSYS命令说明。 §模态提取方法 典型的无阻尼模态分析求解的基本方程是经典的特征值问题： 其中： =刚度矩阵, =第阶模态的振型向量（特征向量）, =第阶模态的固有频率（是特征值）, =质量矩阵。 有许多数值方法可用于求解上面的方程。ANSYS提供了7种方法模态提取方法，下面分别进行讨论。 1.分块Lanczos法 2.子空间（Subspace）法 Dynamics法

ANSYS关于接触刚度

【原创】为什么在接触分析中要自己定义接触刚度呢？ 决定接触刚度 所有的接触问题都需要定义接触刚度，两个表面之间渗量的大小取决了接触刚度，过大的接触刚度可能会引起总刚矩阵的病态，而造成收敛困难，一般来谘，应该选取足够大的接触刚度以保证接触渗透小到可以接受，但同时又应该让接触刚度足够小以使不会引起总刚矩阵的病态问题而保证收敛性。 程序会根据变形体单元的材料特性来估计一个缺省的接触刚度值，你能够用实常数FKN来为接触刚度指定一个比例因子或指定一个真正的值，比例因子一般在0.01和10之间，当避免过多的迭代次数时，应该尽量使渗透到达极小值。 为了取得一个较好的接触刚度值，又可需要一些经验，你可以按下面的步骤过行。 1、开始时取一个较低的值，低估些值要比高估些值好因为由一个较低的接触刚度导致的渗透问题要比过高的接触刚度导致的收敛性困难，要容易解决。 2、对前几个子步进行计算 3、检查渗透量和每一子步中的平衡迭代次数，如果总体收敛困难是由过大的渗透引起的（而不是由不平衡力和位移增量引起的），那么可能低估了FKN的值或者是将FTOLN的值取得大小，如果总体的收敛困难是由于不平衡力和位移增量达到收敛值需要过多的迭代次数，而不是由于过大的渗透量，那么FKN的值可能被高估。 4、按需要调查FKN或FTOLN的值，重新分析。(ANSYS公司的资料） 我的理解：接触刚度与接触面等材料属性无关，理论上接触刚度越大越好，尽量小的接触渗透。但难收敛。 通过共享实常数来判别接触对。要注意使用一个contact element 和一个target element共享实常数。 如： type,1 ! defined 1 as a contact element real,1 mat,1 !mesh type,2 !defined 2 as a target element real,1 mat,1 !mesh 在有限元分析中，接触单元通常用来描述两物体相互接触或滑动的界面。近年来，ANSYS开发了一系列的接触单元。刚开始有节点对节点单元CONTAC12和CONTAC52，接着有节点对地单元CONTAC26，然后有节点对面单元CONTAC48和CONTAC49。最近几年，我们引入一类面对

ANSYS梁单元如何提取应力

ANSYS梁单元如何提取应力 问题1：梁BEAM188的应力如何提取？最大、最小应力如何提取？ 就这个问题我上网找了一下，主要在https://www.wendangku.net/doc/e114516268.html,上找，找了几个帖子，先贴过来看看ANSYS梁单元如何提取应力 1、如何取出梁单元中的最大应力作为优化参数值? https://www.wendangku.net/doc/e114516268.html,/forum/viewthread.php?tid=133155&pid2=650360&keywords= beam188%20%E6%A2%81%20%E5%BA%94%E5%8A%9B&searchstyle=3&i ssearch=true#pid650360 问： 我用的188单元作谐响应分析,求解结束后,我想取出模型中的最大应力值作为参数,然后在接下来的优化当中用该最大应力作为状态变量,请问我应该怎么做啊,注意优化时,对应于每组参数值,最大应力点的位置都可能不同. 请高手指点一下谢谢 ---------- 以下程序段分别得到目标变量（总体积），约束变量SV的最大应力值。 /POST1 SET, NSORT,U,Y *GET,DMAX,SORT,,MAX ETABLE,VOLU,VOLU ETABLE,SMAX_I,NMISC,1 ETABLE,SMAX_J,NMISC,3 ssum *GET,VOLUME,SSUM,,ITEM,VOLU ESORT,ETAB,SMAX_I,,1 ！按照单元SMAX_I的绝对值大小进行排序 *GET,SMAXI,SORT,,MAX ESORT,ETAB,SMAX_J,,1 *GET,SMAXJ,SORT,,MAX SMAX=SMAXI>SMAXJ ！约束变量SV：SMAX=最大应力值 FINISH =============== 你这个程序段是针对beam3 吧，对 beam188好像不行。 对beam188，要求所有单元的最大、最小应力可以用命令 allsel

ANSYS动力学分析中的矩阵

1.3.1.ANSYS动力学分析中的矩阵 1.3.1.1.质量矩阵[M] 动力学分析中需要质量矩阵[M]。ANSYS有2种类型的质量矩阵：一致质量矩阵和集中质量矩阵。 一致（consistent）质量矩阵通过单元的形函数计算出来。矩阵的对角线和非对角线上一般均有元素值。一致质量矩阵是大多数单元的缺省选项。 集中（lumped）质量矩阵的质量被单元各节点平分，对角线上元素不为0，非对角线上元素均为0。通过分析选项来激活，命令LUMPM,Key。默认KEY=off，单元的一致质量矩阵，大多数单元的缺省选项。KEY=on，使用集中质量矩阵。其GUI为Main Menu>Solution>Analysis Type>Analysis Options。 一致质量矩阵对大多数分析来说，精度比集中质量矩阵高，但其计算时间较长。 若结构在一个方向的尺寸与另两个方向相比很小时，可采用简化质量矩阵（如果可能得到的话）或集中质量矩阵例如细长的梁或很薄的壳。集中质量矩阵在大规模模型以及高速动力学如波传播问题上，具有较好的节省计算时间的优势。 1.3.1. 2.阻尼矩阵[C] 1.3.1. 2.1.阻尼的分类 阻尼的作用机制比较复杂。大多简化处理。阻尼的数值主要取决于材料、运动速度和振动频率。阻尼可分类如下：粘性阻尼，滞后或固体阻尼，库仑或干摩擦阻尼。 一、粘性阻尼 粘性阻尼一般物体在液体中运动时发生。由于阻尼力与速度成正比，比例常数 c 称作阻尼常数。通过规定阻尼比ξ，Rayleigh阻尼常数α（后面将进行讨论），或定义带有阻尼矩阵的单元，可将粘性阻尼纳入考虑。 通常用阻尼比ξ（阻尼常数 c 对临界阻尼常数cc*的比值）来量化表示粘性阻尼：ξ=c/cc。其中粘性阻尼c，临界阻尼常数cc。临界阻尼定义为出现振荡和非振荡行为之间的阻尼的极值, 此处阻尼比= 1。对一个质量为m ，频率为w的单自由度弹簧质量系统，cc = 2mw 。 注意: 阻尼比ξ = 对于螺栓或铆钉链接结构常常为2%到15% 二、滞后和固体阻尼 滞后和固体阻尼是材料的固有特性，也称结构阻尼。在动力学分析中应该考虑，认识还不是很透彻，因此很难定量的确定。通过规定另一种Rayleigh 阻尼常数β（后面将进行讨论）可将滞后或固体阻尼纳入考虑。 三、库仑或干摩擦阻尼 库仑或干摩擦阻尼是物体在干表面上滑动时产生的阻尼。阻尼力与垂直于表面的力成正比。比例常数m 就是摩擦系数。动力学分析中一般不予考虑。通过规定带有摩擦性能的接触表面单元和间隙单元，可将库仑阻尼纳入考虑（可参见ANSYS 结构分析指南）。 1.3.1. 2.2.ANSYS11中的阻尼矩阵 阻尼矩阵可以用于谐响应、有阻尼模态分析和瞬态分析，以及子结构生成。可以在一个模型中指定多种阻尼形式，程序将加总指定的阻尼形式后形成阻尼矩阵[C]。ANSYS中总阻尼矩阵为： （15-20）其中

ansys提取最大应力值

ANSYS提取最大应力最大位移值 用ansys的apdl方式如何直接获得最大应力和最大位移点的节点编号？ 在后处理中，用Query Results可以找到最大、最小应力和位移的节点号，及其相应值。然后用*Get提取该节点号的各项计算值。 最大应力节点编号的提取： allsel nsort,s,eqv,0,0,all *get,max_eqv,sort,0,imax 最大位移节点编号的提取： allsel nsort,u,sum,0,0,all *get,max_u,sort,0,imax 在书上看到GUI的操作如下： （1）Main Menu>General PostProc>Element Table>Define Table, 在弹出来的对话框中User label for item输入变量名，假定为smaxe，下面Results datas item中左侧的框内选By sequence num 右侧选LS, 并在下方出现的LS,后面输入1 ，单击OK （2）General Postproc>List Results>Sorted Listing>Sort Elems

在弹出来的对话框中由上到下依次为“Descending order”“yes”“（空着）”“smaxe”，单击OK （3）Utility Menu>Parameters>Get Scalar Data 在弹出来的对话框中选择Results Data>other operations，单击OK就完成了。 先用 Nsort 将位移排序，再用 *get 取得最大值，比如，要查找 Ux 的最大值： NSORT,U,x,0 ! 降序排列 *get,ux_max,SORT,0,MAX ! 最大值 *get,ux_imax,SORT,0,IMAX ! 对应节点号 在后处理中，用Query Resulys可以找到最大、最小应力和位移的节点号，及其相应值。然后用*Get提取该节点号的各项计算值。 最大应力节点编号的提取： allsel nsort,s,eqv,0,0,all *get,max_eqv,sort,0,imax 最大位移节点编号的提取： allsel nsort,u,sum,0,0,all *get,max_u,sort,0,imax