ANSYS分析结果的后处理知识.ppt

后处理是用图形、数据的方式显示有限元计算的结果，以供用户进行判断、分析，从而检查有限元模型的合理性。这是有限元结构分析的重要环节,5.1后处理,ANSYS提供了许多方式对模型施加加载，而且借助于载荷步选项，用户可以在求解中逐步对模型施加载荷。在ANSYS的术语中，载荷(Loads)包括边界条件和外部或者内部作用力。,4.1.1载荷种类,在不同的学科中，载荷的具体含义也不尽相同，下面为不同学科中所指的载荷术语。结构分析：位移、力、压力、温度、重力；热力分析：磁场、热流速率、对流、内部热生成、无限表面等；磁场分析：磁场、磁通量、磁场段、源流密度、无限表面；电场分析：电势（电压）、电流、电荷、电荷密度、无限表面等；流体分析：流速、压力等对不同学科的载荷而言，程序中的载荷可以分为六类：（1）DOFconstraint(DOF约束）:定义节点的自由度值，也就是将某个自由度赋予一个已知值。在结构分析中该约束被指定为位移和对称边界条件；在热力分析中被指定为温度和热通量平行的边界条件。,（2）Force（集中载荷或力载荷）：施加于模型节点上的集中载荷。在结构分析中被指定为力和力矩；在热分析中为热流速率；在磁场分析中为电流段。（3）Surfaceload(表面载荷)：为施加于模型某个表面上的分布载荷。在结构分析中被指定为压力；在热分析中为对流和热通量。（4）Bodyload(体积载荷)：为施加于模型上的体积载荷或者场载荷。在结构分析中为温度；热力分析中为热生产率。（5）Inertiaload（惯性载荷）：由物体的惯性引起的载荷，如重力加速度、角速度、角加速度。主要在结构分析中使用。（6）Coupled-fieldloads（耦合场载荷）：为以上载荷的一种特殊情况，是从一种分析得到的结果作为另一种分析的载荷。例如可以施加磁场分析中计算的磁力作为结构分析中的力载荷。,1)载荷步(loadstep)载荷步就是我们平时讲的分步施加加载荷，以模拟真实的载荷配置。左图所示显示了一个需要三个载荷步的载荷历程曲线：第一个载荷步用于线性载荷，第二个载荷步用于不变载荷，第三个载荷步用于卸载。载荷值在载荷部的结束点达到全值。2)载荷子步子步(Substep):将一个载荷步分成几个子步施加载荷。3)时间的作用在所有静态和瞬态分析中，ANSYS使用时间作为跟踪参数，而不论分析是否依赖于时间。在指定载荷历程时，在每个载荷步的结束点赋予时间,4.1.1.1载荷步和子步,值。时间也可作为一个识别载荷步和载荷子步的计算器。这样计算得到的结果也将是与时间有关的函数，只不过在静力分析中，时间取为常量0；在瞬态分析中，时间作为表示真实时间历程的变量在变化；在其它分析中，时间仅作为一个计算器识别求解时所采用的不同载荷步。,从时间的概念上讲，载荷步就是作用在给定时间间隔内的一系列载荷；子步为载荷步中的时间点，并在这些点上求得中间解。,在ANSYS程序中，用户可以把载荷施加在实体模型（关键点、线、面、体等）上，也可以施加在有限元模型（结点、单元）上。如果载荷施加在几何模型上，ANSYS在求解前先将载荷转化到有限元模型上。这两种情况各有各自的优缺点。(1)施加在实体模型上优点:（1）模型载荷独立于有限元网格之外，这样就不必因为网格重新划分而重新加载；（2）通过图形拾取来加载时，因为实体较少，所以施加载荷简易。缺点:（1）不能显示所有的实体模型载荷；（2）施加关键点约束的扩展时，在两个关键点施加的约束会扩展到关键点之间的直线上所有的节点上，有时这种约束并不是实际的约束情况，因此在使用扩展约束时，在关键点上施加约束要特别小心。,4.1.2加载方式及其优缺点,(2)施加在有限元模型上优点：（1）约束可以直接施加在节点上，所以扩展约束没有影响；（2）载荷可以直接施加在节点上缺点：（1）任何对于有限元网格的修改都会使已施加的载荷无效，需要删除先前的载荷并在新网格上重新施加载荷；（2）不方便处理线载荷和面载荷，因为原来施加在一条线上的载荷需要逐个结点来拾取，原来施加在一个面上的载荷需要逐个单元来拾取，非常麻烦。实例：connect.db,任何实际结构都会受到一定的约束条件来保持其稳定性，因此给结构模型施加合适的约束条件是进行有限元分析的一个重要步骤。在结构分析涉及到的所有载荷中，惯性载荷相对于整体相对于整体笛卡儿坐标系施加于整个模型，除此之外，其它载荷既可以施加于实体图元（点、线、面），也可以施加在有限元模型上（结点、单元）。载荷可以进行施加(Apply)、删除(Delete)、运算(Operate)。施加载荷可以通过前处理器Preprocessor或求解器Solution中的Loads项完成。左图示菜单为第一种方式。弹出相应对话框后：,4.2载荷的施加,(1)选择载荷形式：如Displacement(位移),Force/Moment(力和力矩),Pressure(压力)、Temperature（温度）等；(2)选择加载的对象：如:OnKeypoints、OnLines、OnAreas、OnNodes、OnElement等；(3)指定载荷的方向和数值,在结构分析中，DOF约束中相应的自由度有平移和旋转，共有6个，即UX、UY、UZ(X、Y、Z方向平动自由度）及RTOX、RTOY、RTOZ（X、Y、Z方向的转动平动自由度），它可以施加在线上，也可以施加在结点上。位移方向与总体坐标轴正向相同时取正值，否则取负值。DOF的复杂情况如下：施加对称或者反对称DOF约束施加耦合DOF约束,4.2.1加载自由度(DOF)约束,4.2.1.1在关键点（或节点）上加载位移约束,命令：DK,KPOI,Lab,VALUE,VALUE2,KEXPND,Lab2,Lab3,Lab4,Lab5,Lab6D,NODE,Lab,VALUE,VALUE2,NEND,NINC,Lab2,Lab3,Lab4,Lab5,Lab6GUI：SolutionDefineLoadsApplyStructuralDisplacementOnKeypoints（或OnNodes）参数说明：KOPI、NODE要施加约束的关键点号、节点号LabUX、UY、UZ、ROTX、ROTY、ROTZ等符号标识；VALUE、VALUE2自由度值、第二个自由度值；NEND，NINC指定相同的约束值到NODE到NEND的节点上（缺省为NODE），其节点号增量为NINC（缺省1）；KEXPND关键点自由度扩展选项（0约束只施加在关键点处的节点上；1将关键点上的约束进行扩展(见下页）Lab2,Lab3,Lab4,Lab5,Lab6附加自由度，这些自由度取相同的值施加在这些节点上。,对话框中的KEXPND选项设为YES，可使相同的约束施加在位于两关键点连线的所有节点上。如左图所示。只要拾取关键点K1和K2，再在设置对话框中选择AllDOF，并在VLAUE框中输入0，设置EXPAND为Yes，则K1至K2之间的所有节点都将被约束，即相当于固定了这条边；KEPANDNo则只固定了K1、K2上的两个节点被约束。,K1,K2,4.2.1.2在线（或面）上加载位移约束,命令：DL,LINE,AREA，Lab,Value,Value2DA,AREA,Lab,Value,Value2参数说明：略GUI：SolutionDefineLoadsApplyStructuralDisplacementOnLines（或OnAreas）,4.2.1.2对称约束与反对称约束,在位移约束中还有两个非常重要的约束，就是施加在对称面上的对称约束和反对称约束。（1）对称约束限制对称面内所有节点的两个方向旋转自由度，同时限制了垂直于对称面的位移自由度。（2）反对称约束限制了对称面内所有节点在对称面内的两个方向位移自由度，同时限制了垂直于对称面的旋转自由度。,对称约束与反对称约束的示意图如下：,这两种约束条件应用在不同的对称模型的场合，通过设定对称约束边界条件达到简化建模的效果。施加对称约束和反对称约束的采用的命令仍然是DL、DA，其参数Lab设为SYMM（对称）或ASYM（反对称）。GUI:SolutionDefineSymmetryB.C.（或AntisymmB.C.）,在结构分析中，集中载荷包括力（FX、FY、FZ，即在X轴、Y轴、Z轴方向的集中力和力矩（MX、MY、MZ，即绕X轴、Y轴、Z轴的力矩）。它们只能施加在关键点和节点上FK，KOPI，Lab，Vlaue，Vlaue2F，NODE，Lab，Vlaue，Vlaue2，MEND，NINCGUI：|LoadsDefineLoadsApplyStructuralPressureOnOnKeypoints（或OnNodes）参数说明：KOPI、NODE关键点、节点Lab：FX，FY，FZ（力）或MX，MY，MZ（力矩）,4.2.2施加集中载荷,正值表示力的方向与坐标轴正向一致，负值表示力的方向与坐标轴正向相反,在结构分析中，表面载荷就是施加的压力。当表面载荷方向指向物体内部时取正值，否则取负值，ANSYS不仅可以将表面载荷施加在线上、面上，还可以将表面载荷施加在结点、单元以及梁上。3.2.3.1在线（或面）上施加分布载荷命令：SFL，线编号，分布载荷类型，VALI,VALJGUI：.|LoadsDefineLoadsApplyStructuralPressureOnLines,4.2.3施加分布载荷,若分布载荷为均布载荷，只需在对话框的第一个输入栏中输入相应的分布载荷值；若同时输入第二个值，则表示在这条线上，从线的起始点到线的终点，沿线的方向，承受从第一个值到第二个值线性过渡的分布载荷。注意：（1）ANSYS中的线是有方向的，相当于从起始关键点到终止关键点的一条矢量线，这在很多分析中非常重要。观察方向从实用菜单PlotCtrlsSymbol中设置Linedirrctionon（2）对于非线性的函数分布载荷，可以通过分段近似线性加载的方法，或者通过不同节点处加载不同集中力的方法进行模拟（使用数组载荷定义）,线上分布力加载的起始方向,3.2.3.2在相连的几个节点上施加分布载荷命令：SF，NList，Lab，VALUE,VALUE2GUI：.|Loads|Apply|StructualPressureOnNodes采用GUI操作，在弹出拾取对话框后，在模型上选取几个相连的节点（要施加分布载荷的节点），单击OK按钮，弹出如下所示分布载荷大小设置对话框：,注意：在节点上施加分布载荷必须选取两个以上的节点,3.2.3.3在单元上施加分布载荷命令：SFE，ELEM，LKEY，Lab，KVAL，VAL1，VAL2，VAL3，VAL4GUI：.|LoadsDefineLoadsApplyStructuralPressureOnElements,参数说明：ELEM施加分布载荷的单元号；All，则选取所有单元；P，则其后所有参数失效；LKEY载荷方向标示，默认1；Lab有效载荷符号（结构分析中LabPRES）采用GUI操作，在弹出拾取对话框后，在模型上拾取需要施加分布载荷的单元后，单击OK按钮，弹出如下所示分布载荷大小设置对话框：如果单元上的载荷是均布的，可以只在VALUE输入框中输入载荷值即可；如果不是均布载荷，则需要输入其它节点上的载荷值。各节点之间的载荷分布规律，按线性变化处理。,4.2.4施加体积载荷,在结构分析当中，体积载荷有温度和流量两种。在ANSYS程序当中，体积载荷可以施加在结点、关键点、线、面、体、单元上GUI：MainMenu|Preprocessor|Loads|Apply|TemperatureMainMenu|Solution|Loads|Apply|TemperatureGUI：MainMenu|Preprocessor|Loads|Apply|Other|FluenceMainMenu|Solution|Loads|Apply|Other|Fluence,4.2.5施加惯性力载荷,在结构分析中，惯性载荷有加速度、角加速度、角速度等。惯性载荷只有在模型具有质量的时候才有效。GUI：MainMenu|Preprocessor|Loads|Apply|GravityMainMenu|Solution|Loads|Apply|Gravity列表显示惯性载荷：UtilityMenu|List|Status|Solution|InertiaANSYS程序不允许用户删除惯性载荷，要取消惯性载荷只需将其值设为0,命令：CP，NSET，Lab，NODE1，NODE2，GUI:PreprocessorCoupling/CeqnCoupleDOFs参数说明：NSET耦合标号；Lab耦合自由度；GUI方式在弹出拾取对话框后，拾取需要耦合的若干节点，单击OK，进入耦合设置对话框，如下图所示：,ANSYS中可以设置耦合约束来模拟铰链、无摩擦滑动器、万向节、无摩擦接触面等问题。一个耦合约束是设置一组被约束在一起，有着相同大小但数值未知的自由度。,4.2.6耦合约束(CoupleDOFs),耦合编号,提示：（1）耦合设置的标号是为了区分不同耦合设置的，必须是一个自然数，每次生成新的耦合都应该输入不同的标号；（2）耦合设定的都是节点，因此耦合操作必须在划分单元后进行，不可以耦合关键点、线和面；对于需要在同一位置的所有节点之间自动生成耦合关系的情况，选择主菜单中：GUI:PreprocessorCoupling/CeqnCoincidentNodes，如左图所示。,在该对话框中选择耦合自由度，然后输入容差的值，单击OK，则完成距离小于容差设定的相同位置节点的自由度耦合。,例如：需要耦合节点301、31、32的X方向自由度，使之成为一个可以沿Y轴滑动的滑动幅，使用如下命令：CP，1，UX，30，31，32,很多力学问题需要加载预应力，例如螺钉的预紧问题等。ANSYS中提供了直接的预应力加载。,4.2.7加载预应力,4.3.1在图形中显示载荷通常情况下，施加在几何实体模型上的载荷，用相应的符号显示在几何实体模型的体、面、线或关键点上。而有限元模型上加载的载荷，则显示在节点或单元上。显示符号的设置如下：命令：/PBC/PSF/PEF/PICE/PSYMBGUI:UtilityMneuPlotcrlsSymbol选择应用菜单中的plotVolumes（Areas，Lines）等命令仅显示几何实体模型加载的情况；若选择plotElement，则仅显示有限元模型加载情况；若选择PlotMulti_Plots则显示全部加载情况。记住：通过图形显示仅仅能够看到加载和约束的基本情况。需要知道具体的加载位置和大小，应选择主菜单中的ListLoads子菜单中的相关命令。4.3.2载荷删除对于已施加在几何模型上或有限元模型上的载荷，可以进行删除操作。命令:LSCLEAR,LabGUI：defineLoadsDelete.,4.3载荷显示与控制,求解的过程大部分是由计算机自动完成的，在完成建模和加载工作后，就可以直接进行ANSYS求解了。命令：/SATTUS,SOLUSOLVEGUI:MainMenuSolutionSloveCurrentLS参数说明：/STATUS命令查看信息；参数SOLU只显示求解的控制信息和载荷信息，否则将显示包括模型信息在内的所有详细信息。SOLVE求解命令在应用ANSYS进行实际分析时，对于到多数求解过程，需要对求解过程进行控制设置，这些设置主要包括：（1）分析类型设置；（2）求解基本选项设置；,4.4求解过程控制,用当前载荷步求解,4.4.1分析类型设置,在建模之前已经作过GUI参数过滤设置，目的在于隐藏无关的菜单和命令，使操作简化。命令：/PMETH，Key，OPTIONGUI：MainMenuPreferences参数说明：Key：ON/OFF，激活/关闭求解选项设置；STAT，显示当前设置OPTION:若KeyON，结构分析时OPTION0另外，结构分析种类也很多，ANSYS默认的结构分析类型是静力分析(Static)，当要改变分析类型时，如模态分析、谐分析、瞬态分析等，需要采用以下方式进行设置：命令：ANTYPE，Antype，Status，LDSTEP，SUBSTEPGUI：MainMenuSolutionAnalysisTypeNewAnalysis,模态分析,谐波分析,瞬态分析,频谱分析,子结构分析,4.4.2求解基本选项设置,求解前使用该对话框将各个控制项设置好，就可以求解了。命令：SOLCONTROL，Key1，Key2，Key3，VtolGUI：MainMenuSolutionAnalysisTypeSolutionCtrl,求解控制对话框常用设置说明：AnalysisOptions：设置线性/非线性问题、小变形/大变形问题TimeControl：时间控制选项。这是求解控制中的一个重要部分。ANSYS在结构分析中，使用时间作为一个过程参数，而不论分析结果是否依赖于时间。这样做的好处是，在所有的情况下，以一个通用的“计数器”来反映分析历程。此外时间的单调增加与分析过程一致。显然，在瞬态分析、蠕变分析中，ANSYS的“时间”代表了实际的时间。然而在不依赖时间的静态分析或模态分析中，时间仅仅是一个识别载荷步与子步的计数器。在时间控制选项中的“Timeatendofloadstep”栏设置的就是这个时间计数器的终止时间，例如将其设置为1。注意设为0或留空，则ANSYS使用默认时间值1。在设置了终止时间后，就可以设置时间子步。所谓时间子步就是将分析过程进行认为分段，使有限元方程在迭代求解过程中分段前进，最终获得目标结果的正确值。例如在”Numberofsubsteps”栏中输入4，则表示将时间过程等分为4个部分进行分段迭代，若终止时间1，则它等同于设置每个时间子步0.25。时间控制中还有一个自动时间子步选项”Automatictimestepping”，如果,将其打开，则ANSYS将会在某一个迭代子步不收敛或很难收敛的情况下自动将该子步进行二分，该设置有助于在不明确求解子步是否设置得合适得情况下使用。另外，在动力学分析中，可以通过该选项设置使用弧长法求解。,4.5载荷步的设置与求解操作,在工程实际问题中，结构的载荷、边界条件可能存在几种工况，或载荷、边界条件周期地随时间阶段变化，因此，为了求得结构各阶段的结果，必须对各阶段的边界、载荷分别求解。为让ANSYS自动完成全部工况的计算分析，引入了载荷步的概念和载荷步文件法自动求解法。一个载荷步是指边界条件和载荷项的一次设置或一种工况。若希望使用载荷步文件法求解，则需要将每一个载荷步的加载情况依次写到相应的载荷步文件中。载荷步文件取名为“jobname.sxx”,其中”xx”表示载荷步号，然后再用一条命令来读入每个载荷步文件并开始逐个求解。4.5.1载荷步文件的建立法（1）定义一个载荷步的边界条件和加载情况；（2）写入第一个载荷步文件：SolutionLoadStepOptsWriteLSFile（3）按第二个边界条件、加载情况修改模型，即定义第二个载荷步（4）写入第二个载荷步文件。直到全部定义完毕。,4.5.2载荷步文件法求解GUI:SolutionSolveFromLSFile,起始载荷步文件号,终止载荷步文件号,载荷步文件号增量,4.5.3载荷步文件的查看或修改、删除,GUI：SolutionLoadStepOptsReadLSFile查看修改GUI:SolutionDefineLoadsOperateDeleteLSFiles,4.6简单实例,4.6.1平面问题：板中圆孔的应力集中实例1：如图所示为承受双向拉伸的板件，其中心位置有一个小圆孔，尺寸(mm)如图所示。,弹性模量E2105N/mm2，泊松比v=0.3拉伸载荷:q=20N/mm平板的厚度：t=20mm,解题思路分析：1.属于平面应力问题2.中心带孔，应使用8节点四边形单元或三角形单元3.注意单位：尺寸mm，力N，故应力N/mm24.最大变形约为0.001mm（忽略孔的影响）,最大应力在孔的顶部和底部，大小约为3.9N/mm2，即3.9MPa。依次检验有限元的分析结果。,一.相关设置设置jobname为bracket、Title为”aExampleforBracket”，且过滤参数为Structural;二.建立网格模型1.创建几何模型：在XY平面内建立一个矩形和圆，并用布尔subtract得到几何模型；2.定义单元类型：选择StructuralSolid选项下的Quad8node82单元，确定返回最后自动得到单元类型为PLANE82。它是8节点的四边形单元，是平面4节点单元PLANE42的高阶形式，更适合有曲线边界的模型。对于本问题，我们需要有厚度的平面应力单元，只需单击单元类型表中的options按钮，弹出PLANE82选项设置窗口，如图所示，在K3对应的方框中选择”Planestrsw/thk“，使得可以设置板的厚度。,3.定义实常数：PreprocessorRealConstantsAdd/Edit/Delete。弹出如下对话框,输入板的厚度20,4.定义材料特性：EX200000，PRXY0.35.定义网格尺寸：网格边长25，划自由网格。如下页所示PreprocessorMeshingSizeCntrlsManualSizeAllAreasPreprocessorMeshingMeshAreasFree或用MeshTool6.保存工作,三.加载和求解1.定义分析类型：SolutionAnalysisTypeNewAnalysis，设为Static，即结构静态分析；2.施加约束：SolutionDefineLoadsApplyStructuralDisplacementonLines，用鼠标点选模型最左侧边，并全部约束(AllDOF)3.施加载荷：板右侧边缘上有一个背离平板的20N/mm的均布线载荷，则均布压力线载荷除以板厚20mm1N/mm2。对对模型右侧边施加1的均布表面压力。SolutionDefineLoadsApplyStructuralPressureOnLines4.求解：SolutionSolveCurrentLS,网格模型图,约束、载荷模型图,三.检查计算结果，观察收敛情况,决定是否修改网格模型1.节点最大应力检查由于最大应力点的应力值有解析解，因此可以检查该点的应力值。首先显示节点编号，找出孔部对应的节点：UtilityMenuPlotNodes（如果没有显示节点编号，则UtilityMenuPlotCtrlsNumbering，打开节点编号），记下与圆顶部对应的节点编号。如图所示列出应力值：GeneralPostprocListResultsNodalSolutionStress，检查所要考察节点