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
再补充一下，
（1）上面*get命令中，将imax改为imin即获得最小应力节点编号或最小位移节点编号；
（2）将imax改为max即获得响应的最大应力值或最大位移值，改为min即为最小值；
（3）将sum改为z，即为z方向最大位移值。
CYCLES USED/ALLOWED = 0.1000E+05/ 8455（实际循环数/许用循环数）= PARTIAL USA GE（局部损伤)=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）研究方法和分析流程 本次计算利用ANSYS软件来进行象鼻岭碾压混凝土拱坝全过程温控仿真计算分析。具体分析流程如下： 1）收集资料：包括工程气象水文资料、大坝体型、热力学参数、工程进度、施工措施、防洪度汛和蓄水等。 2）整理分析资料：参数拟合、分析建模方法。 3）建模：采用ANSYS软件进行建模，划分网格。 4）编写计算批处理程序：根据资料结合模型编写计算温度场的ANSYS批处理程序。 5）检查计算批处理程序：首先检查语句，然后导入计算模型检查所加荷载效果。 6）计算温度：使用ANSYS软件温度计算模块进行计算。 7）分析温度结果：主要分析各时刻的温度场分布和典型温度特征值。 8）应力计算建模：模型结构尺寸与温度分析模型相同，需要改变把温度分析材料参数改为应力分析材料参数。 9）计算应力：使用ANSYS软件温度应力计算模块和自编的二次开发软件进行计算。 10）分析应力结果：主要分析应力场分布和典型应力特征值。 11）编写报告：对计算流程和结果实施进行提炼总结，提出可行的温控指标和措施。 （2）前处理 1）建模方法选择。 有限元建模一般有两种方法：一种为通过点线面几何拓扑的方法建模，这种建模方法精确，但是比较费时。对于较大规模的建模任务花费时间太多。另一种为通过其他软件导入，如CAD，通过在其他软件中建模，然后输出为ANSYS 可以识别的文件类型，再导入ANSYS中完成建模过程，这种建模方式精度较直接建模的精度要稍低一些，但是由于要求建模的模型已经在CAD软件中完成了
初步建模，可以直接拿来稍作处理即可应用，时间花费较少。本计算选用从CAD 软件导入ANSYS中来建立模型。 2）建模范围。 建模范围可以分为全坝段建模和单坝段建模，全坝段建模可以全面反映整个坝体的温度和应力情况，但是建模难度高、计算量大；单坝段建模建模难度小，计算量也相对较小，一般情况下单坝段建模即可满足要求。 3）施工模拟层厚。 根据已建碾压混凝土坝经验，碾压层厚一般为0.3m左右，按照0.3m一层建模是最精确的，但是如果按照0.3m一层建模，计算网格数量巨大，计算时间长，对于硬件要求较高，在硬件和时间达不到要求的情况下，按照3m一层以下精度都是可以基本满足要求的。 4）分区模拟。 由于各分区混凝土水化热差别较大，对于温度计算影响较大，因此建模要尽量反映混凝土大坝内部分区变化。基岩由于对混凝土只是导热作用，且影响范围在10m左右，因此在计算时可以认为是均质体，计算热力学参数采用靠近建基面的地层参数。 5）参数选取。 参数一般选择可研阶段的材料试验报告，如果项目部未能提供这些资料，可以在征求同意的前提下，通过查阅相关书籍，尽量采取相似工程的资料。 （3）计算 1）ANSYS计算模块。 ANSYS计算温度场模块由其自带，可以直接进入模块计算。 2）化学产热模拟。 通过ANSYS中产热命令BFE模拟。 3）边界条件模拟。 ①对流边界条件通过命令SFA模拟。 ②接触散热边界条件通过命令D模拟。 4）浇筑模拟。 通过ANSYS中的生死单元功能实现，初始阶段所有单元均为死单元，死单
ANSYS热应力分析经典例题
ANSYS热应力分析例题 实例1圆简内部热应力分折： 有一无限长圆筒，其核截面结构如图13—1所示，简内壁温度为200℃，外壁温度为20℃，圆筒材料参数如表13．1所示，求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长，忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图13—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元，对圆筒进行热分析，然后将热分析单元转化为相应的结构单元，对圆筒进行结构分析；直接法是采用热应力藕合单元，对圆筒进行热力藕合分析。 /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours Plnsol,s,eqv $finish
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热应力分析实例 当一个结构加热或冷却时，会发生膨胀或收缩。如果结构各部分之间膨胀收缩程度不同，和结构的膨胀、收缩受到限制，就会产生热应力。 7.1热应力分析的分类 ANSYS提供三种进行热应力分析的方法： 在结构应力分析中直接定义节点的温度。如果所以节点的温度已知，贝U可以 通过命令直接定义节点温度。节点温度在应力分析中作为体载荷，而不是节点自由度 间接法。首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中。 直接法。使用具有温度和位移自由度的耦合单元，同时得到热分析和结构应力分析的结果。 如果节点温度已知，适合第一种方法。但节点温度一般是不知道的。对于大多数问题，推荐使用第二种方法一间接法。因为这种方法可以使用所有热分析的功能和结构分析的功能。如果热分析是瞬态的，只需要找出温度梯度最大的时间点，并将此时间点的节点温度作为荷载施加到结构应力分析中去。如果热和结构的耦合是双向的，即热分析影响结构应力分析，同时结构变形又会影响热分析（如大变形、接触等），则可以使用第三种直接法一使用耦合单元。此外只有第三种方法可以考虑其他分析领域（电磁、流体等）对热和结构的影响。 7.2间接法进行热应力分析的步骤 首先进行热分析。可以使用热分析的所有功能，包括传导、对流、辐射和表面效应单元等，进行稳态或瞬态热分析。但要注意划分单元时要充分考虑结构分析的要求。例如，在有可能有应力集中的地方的网格要密一些。如果进行瞬态分析，在后处理中要找出热梯度最大的时间点或载荷步。 表7-1热单元及相应的结构单元
重新进入前处理，将热单元转换为相应的结构单元，表7-1是热单元与结构 单元的对应表。可以使用菜单进行转换： Mai n Menu>Prep roeessor>Eleme nt Typ e>Switeh Eleme nt Type ，选择Thermal to Struetual 。 但要注意设定相应的单元选项。例如热单元的轴对称不能自动转换到结构单元中，需要手工设置一下。在命令流中，可将原热单元的编号重新定义为结构单元，并设置相应的单元选项。 设置结构分析中的材料属性（包括热膨胀系数）以及前处理细节，如节点耦 合、约束方程等。 读入热分析中的节点温度, GUI: Solution>Load Apply>Temperature>From Thermal Analysis 。输入或选择热分析的结果文件名*.rth。如果热分析是瞬态的，则还需要输入热梯度最大时的时间点或载荷步。节点温度是作为体载荷施加的，可通过Utility Men u>List>Load>Body Load>On all nodes 列表输出。 设置参考温度，Mai n Men u>Solutio n>Load Setti ng>Refere nee Temp 。 进行求解、后处理。 7.3间接法热应力分析实例 7.3.1 问题描述 图7-1冷却栅示意图
ansys各种结构单元介绍
一、单元分类 MP - ANSYS/Multiphysics DY - ANSYS/LS-Dyna3D FL - ANSYS/Flotran ME - ANSYS/Mechanical PR - ANSYS/Professional PP - ANSYS/PrepPost ST - ANSYS/Structural EM - ANSYS/Emag 3D ED - ANSYS/ED
LINK1 —二维杆单元 单元描述： LINK1单元有着广泛的工程应用，比如：桁架、连杆、弹簧等等。这种二维杆单元是杆轴方向的拉压单元，每个节点有2个自由度：沿节点坐标系x、y方向的平动。就象在铰接结构中的表现一样，本单元不承受弯矩。单元的详细特性请参考理论手册。三维杆单元的描述参见LINK8。 下图是本单元的示意图。 PLANE2 —二维6节点三角形结构实体单元 单元描述： PLANE2是与8节点PLANE82单元对应的6节点三角形单元。单元的位移特性是二次曲线，适合于模拟不规则的网格（比如由不同的CAD/CAM系统得到的网格）。 本单元由六个节点定义，每个节点有2个自由度：沿节点坐标系x、y 方向的平动。本单元可作为平面单元（平面应力或平面应变）或者作为轴对称单元使用。本单元还具有塑性、蠕变、膨胀、应力刚化、大变形、大应变等功能。详细特性请参考理论手册。 下图是本单元的示意图。
BEAM3二维弹性梁单元 BEAM3是一个轴向拉压和弯曲单元，每个节点有3个自由度：沿节点坐标系x、y方向的平动和绕z轴的转动。单元的详细特性请参考理论手册。其它的二维梁单元是塑性梁单元（BEAM23）和变截面非对称梁单元（BEAM54)。 下图是本单元的示意图。 BEAM4三维弹性梁单元 单元描述： BEAM4是一个轴向拉压、扭转和弯曲单元，每个节点有6个自由度：沿节点坐标系的x、y、z方向的平动和绕x、y、z轴的转动。本单元具有应力刚化和大变形功能。在大变形（有限转动）分析中允许使用一致切线刚度矩阵选项。本单元的详细特性请参考理论手册。变截面非对称弹性梁单元的描述参见BEAM44，三维塑性梁单元的描述参见BEAM24。
梁结构应力分布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引言 在现代机械工程设计中，梁是运用得比较多的一种结构。梁结构简单，当是受到复杂外力、力矩作用时，可以手动计算应力情况。手动计算虽然方法简单，但计算量大，不容易保证准确性。相比而言，有限元分析方法借助计算机，计算精度高，
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热应力分析命令流
/FILNAME,Double,1 ！定义工作文件名。 /TITLE,Temperature Analysis ！定义工作标题。 !* /PREP7 !定义单元。 ET,1,SOLID70 !* !定义材料属性。 MPTEMP,,,,,,,, !定义材料1。 MPTEMP,1,0 MPDATA,KXX,1,,238*3.6 !定义材料1的传热系数KXX1。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,500 !定义材料1的比热C1。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,3e-6 !定义材料1密度DENS1。 !* MPTEMP,,,,,,,, !定义材料2。 MPTEMP,1,0 MPDATA,KXX,2,,15*3.6 !定义材料2的传热系数KXX2。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,2,,100 !定义材料2的比热C2。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,2.2e-6 !定义材料2密度DENS2。 !* !建立几何模型。 BLC4,-80,-10,160,20,700 VOFFST,3,20, , !* !网格划分。 FLST,5,20,4,ORDE,2 FITEM,5,1 FITEM,5,-20 CM,_Y,LINE LSEL, , , ,P51X
CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1,10, , , , , , ,1 !定义网格大小。 !* TYPE, 1 MAT, 1 REAL, ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 TYPE, 1 MAT, 2 REAL, ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 2 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !定义网格大小完成。 !* FINISH /SOL ANTYPE,4 !定义瞬态分析类型。
ANSYS_Beam188单元应用
Beam188/189单元基于Timoshenko梁理论（一阶剪切变形理论：横向剪切应变在横截面上是常数，也就是说，变形后的横截面保持平面不发生扭曲）而开发的，并考虑了剪切变形的影响，适合于分析从细长到中等粗细的梁结构。该单元提供了无约束和有约束的横截面的翘曲选项。 Beam188是一种3D线性、二次或三次的2节点梁单元。Beam189是一种3D二次3节点梁单元。每个节点有六个或者七个自由度，包括x、y、z 方向的平动自由度和绕x、y、z 轴的转动自由度，还有一个可选择的翘曲自由度。该单元非常适合线性、大角度转动或大应变非线性问题。 beam188的应力刚化选项在任何大挠度分析中都是缺省打开的，从而可以分析弯曲、横向及扭转稳定问题(进行特征值屈曲分析或（采用弧长法或非线性稳定法）破坏研究)。 Beam188/beam189单元支持弹性、塑性，蠕变及其他非线性材料模型。这种单元还可以采用多种材料组成的截面。该单元还支持横向剪力和横向剪应变的弹性关系，但不能使用高阶理论证明剪应力的分布变化。下图是单元几何示意图：该单元的几何形状、节点位置、坐标体系和压力方向如图所示，beam188 由整体坐标系的节点i 和j 定义。 对于Beam188梁单元，当采用默认的KEYOPT(3)=0，则采用线性的形函数，沿着长度用了一个积分点，因此，单元求解量沿长度保持不变；当KEYOPT(3)=2，该单元就生成一个内插节点，并采用二次形函数,沿长度用了两个积分点，单元求解量沿长度线性变化；当KEYOPT(3)=3，该单元就生成两个内节点，并采用三次形函数,沿长度用了三个积分点，单元求解量沿长度二次变化； 当在下面情况下需要考虑高阶单元内插时，推荐二次和三次选项： 1）变截面的单元； 2）单元内存在非均布荷载（包含梯形荷载）时，三次形函数选项比二次选项提供更好的结果。（对于局部的分布荷载和非节点集中荷载情况，只有三次选项有效）； 3）单元可能承受高度不均匀变形时。（比如土木工程结构中的个别框架构件用单个单元模拟时） Beam188单元的二次和三次选项有两个限制： 1）虽然单元采用高阶内插，但是beam188的初始几何按直线处理； 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) 尽管节点值不用单元表也可以很容易地列出和显示，你仍然可能需要利用单元表储存这些节
ansys梁单元
当一个结构构件的一个方向尺寸远远大于另外两个方向的尺寸时，3D构件就可以理想化为1D构件以提高计算效率。这样的单元有两类：以承受轴向拉压作用为主的杆单元，和承受弯曲作用为主的梁单元。 ANSYS提供的单元类型中共有9种梁单元，分别为BEAM3, BEAM4, BEAM23, BEAM24, BEAM44, BEAM54, BEAM161, BEAM188, BEAM189。在结构分析中常用的是BEAM4和BEAM188或BEAM189这三中梁单元。 BEAM4单元 1.BEAM4单元是一种具有拉压弯扭能力的3D弹性单元。每节点6个自由度。 2.BEAM4单元的定义包括：几何位置的确定，单元坐标系的确定，截面特性 的输入。 BEAM4单元包含两个节点(i，j)或三个节点(i，j，k)，k为单元的方向节点；单元的截面特性用实常数(REAL)给出，主要包括截面(area)，两个 方向的截面惯性矩(IZZ)和(IYY)，两个方向的厚度(TKY和TKZ)，相对单元坐标系x轴的方向角(THETA)，扭转惯性矩(IXX)。其中惯性矩，厚度，方向角都是在单元坐标系下给出的。 3.BEAM4单元坐标系的方向确定如下：单元坐标系X轴由节点i，j连线方 向确定由i指向 j；对于两节点确定的BEAM4单元，若方向角theta=0，则单元坐标系y轴默认平行于整体坐标系的x-y平面；若单元坐标系x 轴与整体坐标系z轴平行，则单元坐标系y轴默认平行整体坐标系的y 轴，z轴由右手法则判定；若用户希望自己来控制单元绕单元坐标系x轴的转动角，则可以通过方向角theta或第三个节点k来实现，i，j，k 确定一个平面，单元坐标系的Z轴就在该平面内。 可以用下列命令查看单元坐标系及截面： /ESHAPE, 1 /PSYMB, ESYS 说明：在指定网格划分属性时，可将某一关键点作为方向点属性赋予所需划分的线，这样就生成包含3个节点的梁单元。(具体见后面) 4.单元压力荷载(pressure)的施加比较特殊。只能用SFBEAM命令来实现， 通过其他方式施加荷载都是无效的，其中LKEY为荷载方向号。 5.beam4单元应力输出：包括轴向正应力，弯曲应力，两者的合应力。 命令：PRESOL，ELEM GUI：LIST RESULT〉ELEM SOLUT〉LINEELEM RESULT
ANSYS梁单元的选择
ANSYS中有七八种梁单元，它们的特点和适用范围各不相同。了解这些单元之间的异同，有助于正确选择单元类型和得到较为理想的计算结果。 梁是一种几何上一维而空间上二维或三维的单元，主要用于模拟一个方向长度大于其它两方向的结构形式。也就是说，主要指那些细长、像柱子一样的结构，只要横截面的尺寸小于长度尺寸，就可以选用梁单元来模拟（这在一定意义上和壳单元在一个方向上比另外两个方向都薄原理相似）。通常来讲，横截面尺寸需要小于长度的1/20或1/30，这里的长度是指两支撑点间的物理意义上的距离。梁单元本身可以进行任意的网格划分，且不支配梁理论的适用性；反过来，就像刚才提到的那样，物理尺寸和特性将决定选择哪种单元更为合适。 有两种基本的梁单元理论：铁木辛格（剪切变形）理论和欧拉-伯努力理论。ANSYS 中的如下单元是基于欧拉-伯努力梁理论： 1.2D/3D elastic BEAM3/4 2.2D plastic BEAM23 3.2D/3D offset tapered,unsymmetric BEAM54/44 4.3D thin-walled,plastic BEAM24 欧拉-伯努力梁理论建立在如下假定的基础上： 1.单元形函数为Hermitian多项式，挠度是三次函数； 2.弯矩可以线性改变； 3.不考虑横截面剪切变形； 4.扭转时截面不发生翘曲； 5.只具有线性材料能力（部分单元BEAM23/24具有有限的非线性材料能力）； 6.非常有限的前后处理能力（除了BEAM44）。 ANSYS中有两种梁单元（BEAM188和BEAM189）是基于铁木辛格（剪切变形）理论，这种梁理论主要建立在如下假定基础上： 1.单元形函数为拉格朗日插值多项式，具有线性或二次的位移函数； 2.横向剪应力沿厚度方向为常数（一阶剪切变形梁单元）； 3.可以模拟自由或约束扭转效应； 4.支持丰富的模型特性（塑性和蠕变）； 5.强大的前生处理能力。 使用中需要注意： （1）铁木辛格（剪切变形）理论是基于一阶剪切变形理论的，它不能准确地求解短粗梁，因此，ANSYS在帮助里指出该类型梁的适用范围是：GAl2/EI>30，对于那些高跨比较大的梁应选用实体单元求解； （2）ANSYS中2结点的铁木辛格（剪切变形）单元BEAM188对网格密度的依赖性较强，选用时单根构件单元数应不小于5或不小于3，并且打开KEYOPT(3)，否则误差会较大。
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梁单元如何提取应力 问题1：梁BEAM188的应力如何提取？最大、最小应力如何提取？ 就这个问题我上网找了一下，主要在https://www.360docs.net/doc/713803369.html,上找，找了几个帖子，先贴过来看看ANSYS梁单元如何提取应力 1、如何取出梁单元中的最大应力作为优化参数值? https://www.360docs.net/doc/713803369.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提取最大应力值
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
ANSYS热应力分析
题目：ANSYS热应力分析 专业：材料成型及控制工程 班级：型职141 学号：14615118 姓名：武学杰 指导教师：张转转 2017年10月17日至11月13日共4周 指导教师(签字) 系主任(签字)
题目： 第一步：更改文件名
第二步：选择单元 第三步：设置材料属性 1、给定材料的导热系数40W(m·℃)。 Main Menu>Preproessor>Material Props>Material Models 第四步：建立实体模型（国际单位制） 1、创建矩形A1:X1，Y1（0，0）、X2,Y2（0.01，0.07）MainMenu>Preprocessor>Modeling>Creaate>Areas>Rectangle>By Dimensions
2、创建矩形A2：X1,Y1(0,0.05)X2,Y2(0.08,0.07) 3、显示面的编号 Utility Menu>PlotCtrls>Numbering
4、对面A1和A2进行overlap操作 Main Menu>Preocessor>Modeling>Operate>Booleans>Overlap>Areas 第五步：划分网格 1、打开Meshtool对话框； Main Menu>Preprocessor>Meshing>MeshTool 2、设定网格尺寸为0.002，网格形状为四面体映射网格； 3、Mesh。
第六步：施加载荷 1、进入Solution处理器。 Main Menu>Solution 2、设定分析类型为“steady-state” Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>On Lines 3、在外面界线上定义温度载荷60 Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>On Lines 4、在内边界线定义温度载荷0。 5、图形窗口显示线。 Utility Menu>Plot>Lines