ANSYS高级培训手册(官方培训)ppt课件

欢迎各位来参加ANSYS进阶培训,主讲人：黄志新肖金花ANSYS公司北京办事处,培训安排,第一天单元库及常用单元、材料库、高级有限元模型技术第二天建模、加载、后处理的高级技术第三天APDL参数化分析技术、优化设计第四天非线性分析第五天ANSYSWorkbench的高级使用,单元库及常用单元,ANSYS单元类型,实体单元梁/管单元壳/膜单元杆/索单元弹簧元接触单元表面效应单元质量单元mesh200,.,.,.,.,.,.,单元类型,常用单元的形状,点(质量),线(弹簧，梁，杆),面(薄壳,二维实体,轴对称实体),线性,二次,体(三维实体),线性,二次,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,.,在单元手册中，ANSYS单元库有200多种单元类型，其中许多单元具有好几种可选择特性来胜任不同的功能。,具体单元名称,单元图示,ANSYS单元名称,单元特性,(类别,编号),单元类型,在结构分析中，结构的应力状态决定单元类型的选择。选择维数最低的单元去获得预期的结果(尽量做到能选择点而不选择线，能选择线而不选择平面，能选择平面而不选择壳，能选择壳而不选择三维实体)。对于复杂结构，应当考虑建立两个或者更多的不同复杂程度的模型。可以建立简单模型，对结构承载状态或采用不同分析选项作实验性探讨。,单元类型,主要单元类型举例,线单元:Beam(梁)单元是用于螺栓(杆)，薄壁管件，C形截面构件，角钢或者狭长薄膜构件(只有膜应力和弯应力的情况)等模型。Spar(杆)单元是用于弹簧，螺杆，预应力螺杆和薄膜桁架等模型。Spring单元是用于弹簧，螺杆，或细长构件，或通过刚度等效替代复杂结构等模型。,单元类型,主要单元类型举例,X-Y平面单元:在整体笛卡尔X-Y平面内（模型必须建在此面内），有几种类型的ANSYS单元可以选用。其中任何一种单元类型只允许有平面应力、平面应变、轴对称、或者谐结构特性。,O,K,N,J,M,P,L,I,I,J,K,L,O,P,N,M,TriangularOption,Y,(orAxial),X(orRadial),单元类型,平面应力假定在Z方向上的应力为零，主要有以下特点：当Z方向上的几何尺寸远远小于X和Y方向上的尺寸才有效。所有的载荷均作用在XY平面内。在Z方向上存在应变。运动只在XY平面内发生。允许具有任意厚度(Z方向上)。,平面应力分析是用来分析诸如承受面内载荷的平板、承受压力或远离中心载荷的薄圆盘等结构。,单元类型,主要单元类型举例,平面应变假定在Z方向的应变为零，主要具有以下特点：当Z方向上的几何尺寸远远大于X和Y方向上的尺寸才有效。所有的载荷均作用在XY平面内。在Z方向上存在应力。运动只在XY平面内发生。,平面应变分析是用于分析那种一个方向的尺寸（指定为总体Z方向）远远大于其它两个方向的尺寸，并且垂直于Z轴的横截面是不变的。,单元类型,主要单元类型举例,轴对称假定三维实体模型是由XY面内的横截面绕Y轴旋转360o形成的（管，锥体，圆板,圆顶盖，圆盘等）。对称轴必须和整体Y轴重合。不允许有负X坐标。Y方向是轴向，X方向是径向，Z方向是周向。周向位移是零；周向应变和应力十分明显。只能承受轴向载荷（所有载荷）。,Hoop,单元类型,主要单元类型举例,谐单元将轴对称结构承受的非轴对称载荷分解成傅立叶级数。傅立叶级数的每一部分独立进行求解，然后根据再合并到一起。谐单元较常用于单一受扭或受弯的分析求解，其中受扭和受弯对应于傅立叶级数的第1和第2项。,单元坐标系(显示的是KEYOPT(1)=0情形),O,K,J,M,P,L,I,I,J,K,L,O,P,N,M,TriangularOption,y,x,N,单元类型,主要单元类型举例,谐单元-举例：假定一承受剪力，弯矩，和/或者扭矩的轴。,轴上的扭矩以傅立叶级数的一项施加到轴上。这时，除了扭矩外，事实上是一般的轴对称问题。,弯矩和横向剪力可以分别作为傅立叶级数的其它两项施加到轴上。,谐单元还可以用于实际当中的任意循环分布载荷，这可能需要分解成50-100项傅立叶级数才能得到满意的结果。,M,V,T,单元类型,主要单元类型举例,主要单元类型举例,壳单元:Shell(壳)单元用于薄面板或曲面模型。壳单元分析应用的基本原则是每块面板的主尺寸不低于其厚度的10倍。,单元类型,三维实体单元:用于那些由于几何、材料、载荷或分析结果要求考虑的细节等原因造成无法采用更简单单元进行建模的结构。四面体模型在用CAD建模往往比使用专业的FEA分析建模更容易，也偶尔得到使用。,K,R,L,Q,O,P,M,N,J,I,Tetrahedronmesh,Brickmesh,单元类型,主要单元类型举例,专用单元:专用单元包括接触单元-用于构件间存在接触面的结构建模，如涡轮盘和叶片，螺栓头部和法兰，电触头，以及O-圈等等。做好接触分析要求有这方面的知识和经验。,单元类型,主要单元类型举例,其它可供选择的单元类型,线性单元/二次单元/p单元:一旦你决定采用平面、三维壳或者三维实体单元，还需要进一步决定采用线性单元、二次单元或P单元。线性单元和高阶单元之间明显的差别是线性单元只存在“角节点”，而高阶单元还存在“中节点”。下面还提到一些差别。,线性单元内的位移按线性变化，因此(大多数时)单个单元上的应力状态是不变的。二次单元内的位移是二阶变化的，因此单个单元上的应力状态是线性变化的。p单元内的位移是从2阶到8阶变化的，而且具有求解收敛自动控制功能，自动分析各位置上应当采用的阶数。,单元类型,我们有必要讲述一下ANSYS中各线性概念之间的区别。线性分析是指不包含任何非线性影响(如：大变形，塑性，或者接触)。线性方程求解器是指方程组解就是结构的自由度解。即使是非线性分析，这些方程还是线性的(但必须进行多次求解)。线性单元假定单元内的自由度按线性变化(跟二次单元,三次单元,或p单元相比)。,单元类型,其它可供选择的单元类型,线性单元/二次单元/p单元(续):在许多情况下，同线性单元相比，采用更高阶类型的单元进行少量的计算就可以得到更好的计算结果。下面是根据不同分析目的进行单元选择的情况。,单元类型,其它可供选择的单元类型,线性单元/二次单元/p单元(续):在进行单元选择时应考虑的其它因素。线性单元的扭曲变形可能引起精度损失。更高阶单元对这种扭曲变形不敏感。就求解的精度的差别讲，线性单元和二次单元网格之间的差别远没有平面单元和三维实体单元网格之间的差别那么惊人之大。所以经常使用线性壳单元。,高度扭曲的二次情形(非平行对边),单元类型,其它可供选择的单元类型,线性单元/二次单元/p单元(续):大多数二次单元允许忽略部分或所有边的中节点-但是，在没有中节点的边上，你只能得到线性结果。如果所有中节点均不存在，该单元就变成了线性单元，计算精度也随之降低(由于转化成线性单元的二次单元和块单元具有“不相容的位移模式”，并引起单元弯曲)。,单元类型,其它可供选择的单元类型,线性单元/二次单元/p单元(续):更高阶的单元模拟曲面的精度就越高。,低阶单元,更高阶单元,单元类型,其它可供选择的单元类型,线性单元/二次单元/p单元(续):采用越来越高阶的单元，给曲线结构划分越来越稀疏的单元网格，ANSYS开始向你发出警告，甚至发出由于单元扭曲变形超过单元允许范围而引起网格划分失败的信息。其原因是，由于模型表面单元的弯曲程度过大，使部分中节点偏离了自身位置，最终决定了你能划分单元网格的稀疏程度。同其它软件一样，ANSYS程序允许用更高阶的直边单元划分网格(降低了实际几何模型的精度，特别是对于p单元而言，通常极不理想)，也允许用不带中节点的更高阶单元划分单元网格(即降低了几何模型的精度，又降低了单元精度，所以在通常情况下更不理想)。所以，一般建议采用尽可能稀疏的单元网格，而又不至于出现形状检查警告。,单元类型,其它可供选择的单元类型,线性单元/二次单元/p单元(续):不能将接触单元同具有中节点的单元连起来(仅对于节点-节点和节点-面接触单元而言-对于面-面接触单元则是允许的)。类似地，在热分析问题中，不能将辐射link单元或者非线性对流表面添加到具有中节点的单元上。,在非线性材料特性区域内，二次单元并不比线性单元更有效。,单元类型,其它可供选择的单元类型,四边形单元/三角形单元，块单元/四面体单元:针对平面或者三维壳体分析模型而言，四边形单元和三角形单元是有差别的，下表列出了这些差异。,单元类型,其它可供选择的单元类型,四边形单元/三角形单元，块单元/四面体单元(续):全部采用三角形单元网格是很少见的。给面进行单元网格划分的实质问题是，你是否允许模型中存在一些三角形单元网格。实际上，各处存在三角形单元会相当麻烦，但是应当仔细思考下列问题:如果采用更高阶单元，三角形单元的计算精度接近于二次单元。所以，全部采用二次单元网格也是不可能的。如果你采用线性单元，三角形单元就十分糟糕-但是，不这样会使四边形单元网格扭曲。除了多数不重要的结构外，任何四边形单元网格(结构的或者非结构的)不得不包含部分形状糟糕的三角形单元网格。所以，还是不可能全部采用四边形单元。,单元类型,其它可供选择的单元类型,四边形单元/三角形单元，块单元/四面体单元(续):对三维实体分析模型而言，块单元和四面体单元是有差别的，下表列出了这些差异。,单元类型,其它可供选择的单元类型,四边形单元/三角形单元，块单元/四面体单元(续):建立三维实体模型需要作出下列选择：使用四面体单元划分网格采用简便方法建立实体模型。选用二次单元或者p单元。或者使用块单元划分单元网格选用块单元网格建立实体模型。通常需要花费更多时间和精力。划分子区域连接处理延伸采用任何块单元。,单元类型,其它可供选择的单元类型,四边形单元/三角形单元，块单元/四面体单元(续):为什么使用四面体单元划分单元网格会有这么大的困难呢？过去，有限元模型全部采用线性四面体单元网格(这种模型十分“粗糙”)。现在，使用二次单元和p单元的有限元模型变得相当理想了。四面体单元模型的自由度几乎是同等精度的块单元单元模型的3到10倍。迄今求解器技术取得了很大突破，大多数分析者还是没有高性能的计算机来求解无关紧要的四面体单元模型。,块单元网格:125个单元216个节点,四面体单元网格:679个单元1230个节点,单元类型,其它可供选择的单元类型,四边形单元/三角形单元，块单元/四面体单元(续):还有其它一些因素帮你作出选择：做接触分析，使用四面体单元划分网格时还需要进行一些处理，消除接触面上的中节点(只针对节点-节点接触单元和节点-面接触单元，而面-面接触单元则不需要)。长或薄结构划分成块单元网格可能更理想。,potentialcontactregion,单元类型,其它可供选择的单元类型,ANSYS分析采用的单位制,除电磁分析以外，你不必为ANSYS设置单位系统。简单地确定你将采用的单位制，然后保证所有输入数据均采用该种单位制就可以(即，ANSYS不能自动进行单位转换)。你确定的单位制将影响尺寸、实常数、材料特性和载荷等的输入值。,单元公式,不同材料行为、不同结构行为选用不同单元公式,高频电磁反射计算,滤波器,单元公式,传统位移方法困难：剪切锁定、体积锁定Solid45KEYOPT(1)=1由于剪切锁定而很少使用非协调模式(附加形函数)Solid45缺省选项，弯曲变形选择缩减积分(B-Bar)几乎不可压缩材料，体积变形一致缩减积分(URI)几乎不可压缩材料，弯曲变形混合U-P公式不可压缩材料，超弹性,为何有如此多的不同单元公式?普通非线性求解非常费时，采用不同的单元技术可更加有效地解决各种类型的非线性问题。不同材料行为(弹性、塑性、超弹性)和不同的结构行为(体积变形、弯曲)需要选择不同的单元公式。,单元公式,单元手册中对每一种单元的定义、特点、适用范围、输入、输出做了详细说明。应该习惯于随时查看单元手册。手册的综述部分应该耐心阅读,单元公式,传统的基于位移的单元有两个问题:剪切锁定和体积锁定:剪切锁定导致弯曲行为过分刚化(寄生剪切应力)。当细的构件承受弯曲时，这是一种几何特性。体积锁定导致过度刚化响应。当泊松比接近或等于0.5时，这是一种材料特性。重点讨论用不同单元公式解决这两个问题的方法。主要讨论连续(实体)单元。由于非线性分析花费计算机时间太多，所以有些单元公式也提供了更有效地解决非线性问题的方法。,传统单元公式,剪切锁定,在弯曲问题中，完全积分低阶单元呈现“过分刚硬”。在弯曲中这种公式包括实际上并不存在的剪切应变，称为寄生剪切。(从纯弯曲中的梁理论可知剪切应变gxy=0),微体积纯弯曲变形中，平直断面保持平直，上下两边变成圆弧，gxy=0。,完全积分低阶单元变形中，上下两边保持直线，不再保持直角，gxy不等于0。,剪切锁定的实例,当长厚比增加时，模型更容易剪切锁定.因为寄生的剪切应变/应力，所以产生的位移被低估。下面的例子是弯曲中的梁。这种情况下剪切应力接近于零，如SXY等高线图中所示，发生了剪切锁定。,体积锁定,材料行为是几乎或完全不可压缩时(泊松比接近或等于0.5)，在完全积分单元中发生体积锁定。超弹材料或塑性流动可发生不可压缩(后面讨论)。单元中产生的伪压应力导致单元对不会引起任何体积变化的变形“过度刚化”。体积锁定也会引起收敛问题。各种应力状态都会发生体积锁定，包括平面应变、轴对称及3-D应力。对平面应力问题不会发生体积锁定，因为平面外应变用于满足体积不可压缩条件。,体积锁定,泊松比接近或等于0.5引起数值上的困难:由于泊松比接近0.5,体积模量无穷大，体积应变接近零。反过来说，很小的体积应变（可能是误差）将会引起极大的静水压力（伪压力）。,由于体积应变由位移的导数计算出，所以其值不如位移精确。体积应变中任何小的误差在静水压力中被放大，这反过来又会影响位移计算。导致不会引起任何体积改变的位移无法产生，网格会锁定。,体积锁定,体积锁定实例,体积锁定可通过压应力“棋盘状”模式(相邻单元间变化显著)检测出。可用单元等值线绘图(PLESOL)绘制静水压力(HPRES)等值线来验证此行为。如怀疑存在体积锁定，可试细分高静水压力区域的网格或改变单元类型。,单元公式,下面的各部分介绍用以克服剪切和体积锁定的单元技术。非协调模式(特殊形状)：形函数，剪切锁定、体积锁定选择缩减积分(B-Bar)：积分方案，体积锁定一致缩减积分(URI)：积分方案，剪切锁定、体积锁定混合U-P公式：特殊自由度，体积锁定,作为一个简单的解释，剪切锁定和体积锁定是由于系统的过度约束。利用不同的单元公式通过放松约束或引入附加的方程求解这些约束来解决这个问题。不幸地是,没有现成的单元公式能最有效地解决锁定问题.因此在下面部分将从正反两方面来讨论每个公式。,18X单元,目前在18x单元中有四个不同的单元技术:B-Bar,URI,增强应变和混合U-P。它们用于处理剪切和体积锁定:高阶18x单元(PLANE183,SOLID186-187)通常用URI。缺省时低阶18x单元(PLANE182,SOLID185)用B-Bar。B-Bar和增强应变不能用于高阶单元。混合U-P技术独立于其它技术,所以可以和B-Bar,增强应变或URI联合使用。,单元公式的选择,单元选项允许用户选择合适的单元公式。MainMenuPreprocessorElementTypeAdd/Edit/Delete“Options”buttonindialogbox若用命令,KEYOPT(1)用于PLANE182的B-bar,URI和增强应变KEYOPT(2)用于SOLID185的B-bar,URI和增强应变KEYOPT(6)用于所有实体/平面18x单元的混合U-P。,增强应变,低阶完全积分单元的形函数可被表示常曲率状态的模式所增强，这些增加的模式作为内在的自由度，因其导致网格的缝隙和重叠而被称为非协调模式。,非协调模式,无非协调模式,F,2F,F,F,2F,F,F,2F,F,F,2F,F,记住增强应变为弯曲和几乎不可压缩应用而设计增强应变不能用于完全不可压缩分析，但对PLANE182和SOLID185可以与混合U-P公式结合使用，在下节讨论。增强应变有上述优点，但更耗费计算机时间前面幻灯片提到的附加内部DOF被凝聚在单元层次，但仍额外消耗计算机时间(和更大的*.esav文件)。只有低阶四边形PLANE182和六面体SOLID185支持增强应变。如果单元扭曲，则增强应变在弯曲中将不利，尤其是梯形单元。,增强应变,选择缩减积分,选择缩减积分(又名B-bar方法,持续膨胀单元)用低一阶的积分方法对体积项积分。应力状态可分解为静水压力(p)和偏差应力(s)两项。上面的方程中,ev是体积应变，ed是偏差应变.k是体积模量,G是剪切模量。,应变通过下式和位移相关:而计算B时,对体积项和偏差项使用不同的积分阶数。,Bv以一个积分点计算(缩减积分),另一方面,Bd以2x2积分点计算(完全积分),选择缩减积分,如前一幻灯片所示，B的体积项和偏差项不是以同一积分阶数计算，只有体积项用缩减积分，这就是该方法称为选择缩减积分的原因。因为B在体积项上平均，因此也称为B-bar法。体积项Bv缩减积分的事实使Bv因为没有被完全积分而软化，这样允许求解几乎不可压缩行为和克服体积锁定。然而，因为偏差项Bd不变，仍然存在寄生剪切应变，所以这个公式仍然容易剪切锁定。具有选择缩减积分的单元有：plane182,solid185,选择缩减积分（体积）,选择缩减积分总结,总之,选择缩减积分在体积变形占优势的问题中对几乎不可压缩材料行为(如塑性,超弹性)有用。单独的B-Bar法对完全不可压缩问题不适用，但可以和混合U-P单元(以后讨论)结合用于完全不可压缩材料。B-Bar法不能用于弯曲占优势的模型。某些单元支持选择缩减积分:可用于平面应变、轴对称和3D应力状态。体积锁定对平面应力不是问题,所以在这种情况下不需要B-Bar法。缺省时PLANE182和SOLID185用B-Bar法(KEYOPT(1)=0)。能用于各种本构模型。,一致缩减积分,一致缩减积分(URI)采用比数值精确积分所需要的阶数低一阶的积分公式这和选择缩减积分类似，但体积和偏差项都用缩减积分。这个公式更灵活，可帮助消除剪切和体积锁定。体积项的缩减积分可以求解几乎不可压缩问题。偏差项的缩减积分防止弯曲问题中的剪切锁定。然而URI可能会引起应变能为零的变形模式，这被称为零能量或沙漏模式。,沙漏模式,沙漏模式是由于变形而引起零应变能的变形模式。如右图所示两例，在只有一个积分点的低阶单元中，此单个积分点未获得任何单元应变能。这可导致出现不切实际的行为。,沙漏模式通常只是低阶URI单元中的问题。只要在每一个方向上有多于一个的单元,高阶URI单元的零能量模式就不会传播。为控制沙漏模式ANSYS使用一个小的沙漏刚度来控制变形的零能量模式。ANSYS为沙漏刚度提供了缺省值。大部分情况下可直接使用缺省值，但也可以用一个实常数缩放因子改变沙漏刚度。任何情况下都应该监控由沙漏模式产生的“虚假能量”，可以用单元表格项AENE来存储“虚假能量”。最好使“虚假能量”与总能量的比值(AENE/SENE)小于5%。,沙漏模式,有URI公式的ANSYS低阶单元包括:Plane182、Solid185、Solid45和Shell181。如果模型中发生沙漏模式，推荐采取的步骤按优先顺序排列如下所示：去掉点载荷和点约束细化网格采用其它可选单元类型增大沙漏刚度缩放因子有URI公式的ANSYS高阶（二次）单元包括:Plane82（采用2x2高斯积分规则）、Solid95（采用2x2x2高斯积分）。只有一个零能量模式，并且只要模型中有不止一个单元，零能量模式就不会传播。推荐大部分应用采用这些单元，因其一般无沙漏模式困难。,一致缩减积分,另一方面，用户在使用URI时需要注意一些事情:低阶URI单元容易沙漏，需要检查。低阶URI单元太柔软，尤其在弯曲占优势的问题中，因此需要细化网格以使位移不被高估。低阶和高阶URI单元的积分公式都比完全积分低一阶。这意味着对低阶单元应力在1点求值，对高阶单元在2x2或2x2x2点求值。因此，需要更多单元来捕捉应力梯度。URI不能用于完全不可压缩分析。,一致缩减积分,缺省时大多数ANSYS高阶结构单元（PLANE82,PLANE183,SOLID186)）用URI，这是因为高阶单元不易沙漏且有许多优点，所以很具吸引力。SOLID95采用修正的14-点积分格式，但当KEYOPT(11)=1时采用URI缺省时大多数低阶单元不采用URI。对SOLID45和SOLID185（KEYOPT(2)=1）或PLANE182（KEYOPT(1)=1）时URI被激活对PLANE42,URI不可用，建议采用支持URI的PLANE182除非特殊需要（如与LS-DYNA单元兼容）,对低阶单元鼓励用户采用B-bar或增强应变代替URI。,一致缩减积分,混合U-P公式,混合U-P单元(又名杂交单元或Herrmann单元)通过内插(并求解)静水压力做为附加自由度来处理体积锁定。单独的内插函数用于位移和静水压力DOF。由于压力可单独求解，所以静水压力的精度和体积应变、体积模量或泊松比无关.ANSYS中有两种方法实现混合u-p对几乎不可压缩用基于惩罚的混合U-P对几乎和完全不可压缩用Lagrange乘子法,基于惩罚的混合U-P,基于惩罚的混合U-P的基本方法是通过体积约束方程把静水压力(p)自由度在单元层次凝聚掉。这样,刚度矩阵仍基于位移而不必担心附加自由度。该公式用于超弹材料(Mooney-Rivlin)的HYPER56,58,74和158也用于支持率相关和率无关塑性(Anand,等向强化)的VISCO106-108该公式可用于几乎不可压缩分析。注意，根据是采用超弹性还是塑性，用户必须选择适当的HYPER或VISCO单元类型。,Lagrange乘子混合U-P,对几乎和完全不可压缩分析采用18x单元，用一个称之为Lagrange乘子法的特殊单元公式。不像基于惩罚的混合U-P公式,Lagrange乘子法将P作为独立自由度来求解。静水压力自由度和内部结点相联系，内部结点由ANSYS自动生成且对于用户是透明的，是不能访问的。该公式用于18x系列单元(KEYOPT(6)0)(PLANE182-183,SOLID185-187)ANSYS将根据材料自动采用适当的公式，因此对用户是透明的。,混合U-P总结,总之,对几乎和完全不可压缩材料,ANSYS提供了丰富的应用混合U-P公式的单元技术库。对几乎不可压缩超弹材料，用HYPER56,58,74,158或混合U-P18x系列单元。对几乎不可压缩弹塑材料，用18x系列的混合U-P公式或VISCO106-108单元。对完全不可压缩超弹材料，用18x单元的混合U-P公式。前面部分中讨论过，18x单元中的混合U-P公式可以和其它单元公式结合。混合U-P本身能解决体积锁定问题对18x单元,可将混合U-P(KEYOPT(6)0)和B-bar,URI或增强应变公式结合。,单元公式,非协调模式：弯曲、体积变形（几乎不可压缩）选择缩减积分(B-Bar)：体积变形（几乎不可压缩）一致缩减积分(URI)：弯曲、体积变形（几乎不可压缩）混合U-P公式：体积变形（完全不可压缩）,实体单元推荐,传统单元容易剪切和体积锁定，ANSYS中有很多单元技术解决这两个问题。通常根据模型选择单元技术，包括弯曲/体积变形和材料行为。只要可能，对非线性问题建议采用18x单元，因为:最新的单元技术和18x单元结合，包括B-bar,URI,增强应变和混合U-P。18x系列的单元技术和材料技术分开。这些单元具有丰富的本构模型，这也有助于缩小单元选择的范围。,实体单元推荐,对高阶单元,缺省时采用URI。用户仅需考虑的是如果材料是完全不可压缩的，应该采用混合U-P。低阶单元选择的一些指南如下:,实体单元推荐,线性分析和小应变非线性分析任何具有附加位移形式的低阶四边形/六面体单元(对PLANE42,SOLID45在非退化形式中缺省)。这些单元对剪切锁定和几乎不可压缩材料行为都有用。任何二阶单元，尤其是需要四面体网格的CAD几何图形的SOLID92(或SOLID187)。高阶四边形/六面体单元如PLANE183或SOLID186采用URI,URI对克服剪切锁定和几乎不可压缩行为也有用。,实体单元推荐,有限应变非线性分析对大应变的应用，首选低阶四边形/六面体单元(不会出现中间结点逆位问题)。先用B-Bar法;如果剪切锁定成为问题，用户可以切换到增强应变。高阶单元(缺省时用URI)也可接受。对18x单元，对几乎或完全不可压缩分析可以采用混合U-PKEYOPT(6)与其它技术的结合。对大应变，需要细化网格和预测大应变区域以确保整个求解过程保持好的单元质量。,壳单元-概述,当结构的总体厚度相对于典型长度很小时可使用壳单元，长度比厚度大20倍以上的问题可决定使用壳单元。ANSYS中的壳单元根据要求解的问题类型采用不同的公式，三个基本的壳公式包括:薄膜理论,“薄”壳理论和“厚”壳理论。,壳单元-概述,薄膜理论Shell41采用薄膜理论。Shell41忽略弯曲和横向剪切，只包含薄膜效应。经典Love-Kirchhoff理论Shell63是“薄”壳单元。Shell63包含弯曲和薄膜效应但忽略横向剪切变形。Reissner/Mindlin理论Shell43,143,181,91,93和99是“厚”壳单元。其包含弯曲、薄膜和横向剪切效应。横向剪切被表示为整个厚度上的常剪切应变,这种一阶近似只适用于“中等厚度”壳体。,平面变形中的壳单元,平面内壳的响应可认为是平面应力状态，因此对于壳单元不会出现体积锁定问题。(当绝对不可压缩，泊松比=0.5时Shell181支持超弹性)对于薄膜现象，壳单元的平面公式与平面实体单元的公式相似(非协调模式)。Shell41,43,63和181对于平面内变形支持非协调模式。Shell181也支持具有沙漏控制的一致缩减积分(缺省选项)。,壳单元推荐,线性分析如壳的厚度非常小采用Shell63，Shell63单元不包含横向剪切效应。如横向剪切变形重要，对于均匀材料采用Shell43,Shell93或Shell143，对于复合材料采用Shell91或Shell99。注意具有一致缩减积分(缺省)的单元Shell181对大模型较快，但需要较细的网格。,壳单元推荐,非线性分析等向强化塑性和超弹性采用Shell181。其优势包括：较小的.esav文件，较少的CPU时间，压力载荷刚度效果，可以导入初始应力，厚度变化。随动强化塑性，蠕变采用Shell143,Shell43和Shell93。Shell143适用于小应变塑性，Shell93是弯曲的壳(高阶)。,梁单元-概述,梁单元可用于分析主要受侧向或横向载荷的结构，长度对横截面的比率超过20:1可作为梁单元应用的原则。典型的梁应用包括：机器主轴，房屋构架，桥梁等。ANSYS中可用的两个梁单元公式为：Euler/Bernoulli梁Beam3和Beam4包括弯曲、轴向和扭转变形。横向剪切变形不包括于单元公式中。Timoshenko梁Beam188和Beam189在单元公式中包括弯曲、轴向、扭转和横向剪切变形。,梁单元推荐,线性分析对于线性模型采用Beam3,Beam4,Beam188或Beam189。Beam3和Beam4采用Hermitian多项式作为形函数并且在弯曲中具有三次响应。Beam188采用线性多项式作为形函数，Beam189采用二次多项式作为形函数。渐变截面采用beam44,beam54注意Beam188和Beam189需要更细化的网格，然而，它们有许多优秀的前后处理特色。,梁单元推荐,非线性分析采用Beam188和Beam189模拟各向同性强化塑性、大应变、屈曲(特征值和非线性屈曲)和/或大转动问题。用作壳加强单元，Beam188与Shell181完全兼容，并且Beam189与Shell93完全兼容。,材料库及常用非线性材料模型,定义材料性质时，首先给出弹性材料性质(EX、PRXY等)。然后给出非线性材料性质。,EX,屈服点,T3,T2,T1,材料属性定义,各向同性材料：EX必须输入泊松比（PRXY或NUXY）默认为0.3GXYEX/(2(1+NUXY)正交各向异性材料所有参数必须输入（EX，EY，EZ，(PRXY，PRYZ，PRXZ，orNUXY，NUYZ，NUXZ)，GXY，GYZ，andGXZ)，无默认值一般各向异性材料直接输入弹性（或柔度）矩阵,线弹性属性定义,PRXY和NUXY的意义、区别：PRXY为主泊松比，指的是在单轴作用下，X方向的单位拉（或压）应变所引起的Y方向的压（或拉）应变。NUXY为次泊松比，它代表了与PRXY成正交方向的泊松比，指的是在单轴作用下，Y方向的单位拉（或压）应变所引起的X方向的压（或拉）应变。对于正交各向异性材料，需要根据材料数据的来源确定数据的输入方式。但是对于各向同性材料来说，选择PR*或NU*来输入泊松比是没有任何区别的。,泊松比的意义,ANSYS材料库：ansys90matlib用户自定义材料库-练习,材料库的运用,非线性材料属性,弹塑性：多种屈服准则：Mises、Hill、广义Hill、Drucker-Prager、Mohr-Coulomb多种强化方式：随动、各向同性、混合双线性、多线性粘塑性：高温金属蠕变：数十种蠕变模型，显式P3=0;P4=0P1是在Xg=0处的值,按P1=2(200)+200=600计算sfe,eflat,4,pres,600,200,0,0,.表面效应单元,面5P1数值的压力矢量.方向=P2,P3,P4现在代表矢量的方向余弦，对数值无影响。例如:sfe,eflat,5,pres,1000,-1,-1,0定义了X-Y面内45度方向的压力,.表面效应单元,内容5(续):压力矢量的幅值也依赖KEYOPT(11).KEYOPT(11)=0(缺省)和1给单元的投影面（垂直于荷载）施加压力。对螺钉荷载有用(或销钉荷载).例如:sfe,ecurv,5,pres,1000,0,-1,0在曲面上定义了一螺栓荷载。如下图位置1轮廓线所示。,.表面效应单元,KEYOPT(11)=2对整个面施加压力。对斜面有用(如屋顶)或风载荷。如:sfe,eslope,5,pres,1000,0,-1,0对整个斜面定义了完全相同的压力。如下图所示,.表面效应单元,练习螺栓扭转,选择求解器,求解器的功能是求解关于结构自由度的联立线性方程组-这个过程可能需要花费几分钟(1,000个自由度)到几个小时或者几天(100,000-1,000,000自由度)，基本上取决于你所用计算机的速度。对于简单分析，可能需要一、两次求解。对于复杂的瞬态或非线性分析，可能需要进行几十次、几百次、或者甚至几千次求解。ANSYS提供了三个求解器用于一般求解：波前求解器（Frontalsolver）和PCG求解器（PCGsolver）（预条件共扼梯度,）。稀疏矩阵求解器（Sparsesolver）也可以使用，主要用于非线性问题。,波前(Wavefront)求解器,波前求解器经常发出“主对角值”或“主元”为小或负的警告或错误信息，指出求解发生奇异。任何一条信息都指出某个特定的自由度从你的约束中忽略掉。,PCG求解器,PCG求解器不检验求解的奇异问题。存在奇异的情况下，它仍可以计算求解，或者结果不收敛，但仍然进行所有的PCG迭代计算，并输出错误信息。,高级后处理技术,单元表的应用,单元表可以获取其他后处理功能不能获得的结果数据。可以利用单元表来对结果数据进行数学计算。查阅单元说明来了解后处理中可以得到的数据。梁弯矩图的绘制,.单元表,一个单元表的标准形成过程由下列几个步骤构成:首先选择需要的单元把结果数据填到表中列出或绘出数据执行所想要的算术运算GUI方式下的GeneralPostprocElementTable.中给出了单元表的功能,.单元表,注意每个单元的每个数据项只有一个值:即该单元的平均或有效值。例如,如果把节点的位移UX添加到单元表,单元节点的平均UX就被存储到每个单元。可以添加3个标量数(如单元x方向的压力，y方向的压力,z方向的压力)，合并他们可获得一矢量图(必须三个分量的单元表都输入）。查看PLVECT命令获取细节,.单元表,执行操作有点象电子数据表，可以在单元表的列间进行算术运算。一些有用的操作:两列相加(SADD)乘(SMULT)比较和保存两列的最大值和最小值(SMAXorSMIN)每列求和(SSUM对选择单元总体积的计算很有用)等等,.单元表,例如:,.单元表,若一个想要的操作无效，可以把单元表中的项上载到一数组参数，使用数组操作。使用*VGET命令（或UtilityMenuParametersGetArrayData)接着使用一个数组操作:*VOPER,*VSCFUN,等例如，象平均和标准偏差这些统计量在单元表水平上无法计算，但可以在数组上计算。,.单元表,更新单元表如果从不同的载荷步骤或从不同的结果文件读进结果，单元表不会自动更新。,执行ETABLE,REFL(或在单元表数据对话框中点击Update按钮)来“填充”或更形单元表。列头部的“current”或“previous”标示依赖数据项的状况。,.单元表,既然单元表不能自动更新，那么就很容易比较两个载荷步的结果或两个不同的分析,如下所示:读载荷步1,把结果项作为A1保存到单元表。接着读进载荷步2，把结果项作为A2保存到单元表。现在可以比较A1和A2，把最大值和最小值保存到第三列，或进行其他想要的操作。这项设置可以被扩展来比较相同模型的不同结果文件。,.单元表,单元表的保存和删除:通过保存ANSYS数据库来保存单元表:在POST1中保存数据库(ToolbarSAVE_DB或SAVE命令).或退出ANSYS时使用“Saveeverything”选项来保存(ToolbarQUIT或/EXIT,ALL命令).删除全部单元表:ElementTableEraseTable或ETABLE,ERASE删除单元表中的一项（列）:使用单元表数据对话框中的Delete按钮.或ETABLE,name,ERASE,查询拾取,查询拾取允许您在模型上“探测”任意拾取位置的应力,位移或其它的结果量.您还可以很快地为查询量的最大值和最小值定位.仅能通过GUI方式操作(无命令):GeneralPostprocQueryResultsNodalorElementorSubgridSolu.选择某个结果量,按OK,PowerGraphicsOFF,PowerGraphicsON,查询拾取,拾取模型中的任一点,以查看该点的结果值.Min和Max将显示最大和最小点的值.使用Reset清除所有值并重新开始拾取查询.注意:实体的编号,位置以及结果值都将显示在拾取菜单中.,自动生成文本注释,rib.inp,路径操作,查看结果的另一种方法是通过路径操作,这一方法允许您:在通过模型的任意一条路径上绘图输出结果数据沿某一路径进行数学运算,包括积分和微分显示一“路径图”观察结果量沿路径的变化情况此方法仅对包含2-D或3-D实体单元或壳单元的模型有效.后面结合例题讲解：定义两条路径分别绘制板的弯矩sx（序列号4），和z向位移。,.路径操作,产生路径图的三个步骤:定义一个路径将数据映射到路径上绘图输出数据1.定义一个路径需要以下信息:定义路径的点(2到1000个).您可以使用工作平面内的节点或特定位置.路径的曲率由激活的坐标系(CSYS)确定.路径名.,.路径操作,1.定义一个路径(续)首先激活需要的坐标系(CSYS).GeneralPostprocPathOperationsDefinePathByNodesorOnWorkingPlane拾取节点或工作平面上的特定位置以形成期望的路径,按OK选取一个路径名.在许多情况下,nSets和nDiv的空上最好为缺省值.,.路径操作,2.将数据映射到路径上GeneralPostprocPathOperationsMapontoPath(或PDEF命令)选定需要的量,诸如SX.为选定的量加入一个用于绘图和列表的标签.如果需要,您可以显示这一路径.GeneralPostprocPathOperationsPlotPaths(或键入命令/PBC,PATH,1续之以NPLOT或EPLOT命令),.路径操作,3.绘图输出数据您既可以采用曲线图绘出路径上的量:PLPATH或GeneralPostprocPathOperationsOnGraph.或沿路径的几何形状:PLPAGM或GeneralPostprocPathOperationsOnGeometry.,.路径操作,ANSYS允许您定义多条路径,您只需为每条路径指定一个唯一的路径名.一次只能有一条路径被激活.除绘图和列表外,还有许多其它的路径功能,包括:应力线性化在压力容器工业中用于将沿某一路径上的应力分解为膜应力及弯曲应力分量.计算功能在断裂力学中用于计算J-积分和应力集中因子.在热分析中用于计算越过某一路径的散失或获得的热量.点积和叉积在电磁分析的矢量操作中有广泛应用.,载荷工况组合,只要您求解多载荷步,每一载荷步的结果将以独立的序列存放在结果文件中(由载荷步号识别).载荷工况组合是两个结果序列之间的操作,这些序列被称为载荷工况.操作发生在数据库中的一个载荷工况和结果文件中的第二个载荷工况之间.操作的结果组合的工况存放回数据库.,.载荷工况组合,典型步骤:1.建立载荷工况2.将某一载荷工况读入数据库3.执行期望的操作,.载荷工况组合,建立载荷工况一个载荷工况可简单地充当一个结果系列的指示器.它需要如下两条信息:唯一的ID号它代表的结果序列(载荷步和载荷子步号)使用LCDEF命令或GeneralPostprocLoadCaseCreateLoadCase,.载荷工况组合,将某一载荷工况读入数据库(内存)简单地采用载荷工况号识别结果序列,使用LCASE命令或GeneralPostprocLoadCaseReadLoadCase.或在后处理中使用一个标准的“读结果”选择(SET命令).,.载荷工况组合,执行期望的操作许多有效的操作如这里的菜单所示使用LCOPER命令或GeneralPostprocLoadCaseAdd,Subtract,等.切记操作的结果存放在数据库(内存)中.,.载荷工况组合,有两种有用的选项用于存储组合的载荷工况:写一个载荷工况文件将载荷工况添加到结果文件写一个载荷工况文件(RAPPND或GeneralPostprocLoadCaseWriteLoadCase)产生一个与结果文件相似,但比它小得多的文件.添加选项(LCWRITE或GeneralPostprocWriteResults)允许您将组合的载荷工况添加到结果文件并用一给定的载荷步号和时间值识别.,切片显示,PlotCtrls-Style-Hidden-LineOptions类型:CappedZ-buffer:显示切面后面的实体。Q-SliceZ-buffer:显示切面及实体轮廓。Section:只显示切面。Z-buffered:默认，显示全部实体。切面：垂直于视线的平面（通过焦点）焦点由PlotCtrls-ViewSettings-FocusPoint设定。工作平面,图形、动画的保存,保存图形：PlotCtrls-CaptureImage,保存为BMP图形格式保存动画PlotCtrls-DeviceOption-Animationmode-Animation-动画图-Animation-SaveAnimation-文件名（AVI格式）,多窗口操作,PlotCtrls-Multiwindowlayout-Multi-plotControlPlot-Multi-Plots,后处理练习,分析结果评价,结果验证,验证分析的结果，在任何有限元分析中无疑是最为重要的步骤.在开始任何分析以前，应该至少对分析的结果有粗略的估计（来自经验、试验、标准考题等）.如果结果与预期的不一样，应该研究差别的原因.,识别无效的结果,应该知道所分析的对象的一些基本的行为:重力方向总是竖直向下的离心力总是沿径向向外的物体受热一般要膨胀没有一种材料能抵抗1000MPa的应力弯曲载荷造成的应力使一侧受压，另一侧受拉如果只有一个载荷施加在结构上，检验结果比较容易.如果有多个载荷，可单独施加一个或几个载荷分别检验，然后施加所有载荷检验分析结果.,计算出的几何项,在输出窗口中输出的质量特性，可能会揭示在几何模型、材料属性（密度）或实常数方面存在的错误.,变形/温度/应力,检验求解的自由度及应力.确认施加在模型上的载荷环境是合理的.确认模型的运动行为与预期的相符-无刚体平动、无刚体转动、无裂缝等.确认位移和应力的分布与期望的相符，或者利用物理学或数学可以解释.,反作用力或节点力,模型所有的反作用力应该与施加的点力、压力和惯性力平衡.,在所有约束节点的竖直方向的反作用力.,在所有约束节点水平方向的反作用力必须与水平方向的载荷平衡所有约束节点的反作用力矩必须与施加的载荷