Ansys基础教程.pptx

《Ansys基础教程.pptx》由会员分享，可在线阅读，更多相关《Ansys基础教程.pptx（85页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第一章 ANSYS主要功能与模块 ANSYS是世界上著名的大型通用有限元计算软件，它包括热、电、磁、流体和结构等诸多模块，具有强大的求解器和前、后处理功能，为我们解决复杂、庞大的工程项目和致力于高水平的科研攻关提供了一个优良的工作环境，更使我们从繁琐、单调的常规有限元编程中解脱出来。ANSYS本身不仅具有较为完善的分析功能，同时也为用户自己进行二次开发提供了友好的开发环境。ANSYS程序自身有着较为强大三维建模能力，仅靠ANSYS的GUI(图形界面)就可建立各种复杂的几何模型；此外，ANSYS还提供较为灵活的图形接口及数据接口。因而，利用这些功能，可以实现不同分析软件之间的模型转换。第1页/共

2、85页一、主要功能简介1.结构分析 1）静力分析-用于静态载荷.可以考虑结构的线性及非线性行为。线性结构静力分析 非线性结构静力分析 几何非线性：大变形、大应变、应力强化、旋转软化 材料非线性：塑性、粘弹性、粘塑性、超弹性、多线性弹性、蠕变、肿胀等 接触非线性：面面/点面/点点接触、柔体/柔体刚体接触、热接触 单元非线性：死/活单元、钢筋混凝土单元、非线性阻尼/弹簧元、预紧力单元等第2页/共85页 2）模态分析-计算线性结构的自振频率及振形.谱分析 是模态分析的扩展，用于计算由于随机振动引起的结构应力和应变(也叫作 响应谱或 PSD).3）谐响应分析-确定线性结构对随时间按正弦曲线变化的载荷的

3、响应.4）瞬态动力学分析-确定结构对随时间任意变化的载荷的响应.可以考虑与静力分析相同的结构非线性行为.5）谱分析 6）随机振动分析等 7）特征屈曲分析-用于计算线性屈曲载荷并确定屈曲模态形状.(结合瞬态动力学分析可以实现非线性屈曲分析.)8）专项分析:断裂分析,复合材料分析，疲劳分析第3页/共85页2.高度非线性瞬态动力分析（ANSYS/LS-DYNA）全自动接触分析，四十多种接触类型 任意拉格郎日欧拉（ALE）分析 多物质欧拉、单物质欧拉 适应网格、网格重划分、重启动 100多种非线性材料模式 多物理场耦合分析：结构、热、流体、声学 爆炸模拟，起爆效果及应力波的传播分析 侵彻穿甲仿真，鸟撞

4、及叶片包容性分析，跌落分析 失效分析，裂纹扩展分析 刚体运动、刚体柔体运动分析 实时声场分析 BEM边界元方法，边界元、有限元耦合分析 光顺质点流体动力（SPH）算法 第4页/共85页3.热分析 稳态、瞬态温度场分析 热传导、热对流、热辐射分析 相变分析 材料性质、边界条件随温度变化 4.电磁分析 静磁场分析计算直流电（DC)或永磁体产生的磁场 交变磁场分析 计算由于交流电(AC)产生的磁场 瞬态磁场分析计算随时间随机变化的电流或外界引起的磁场 电场分析用于计算电阻或电容系统的电场.典型的物理量有电流密度、电荷密度、电场及电阻热等。高频电磁场分析用于微波及RF无源组件，波导、雷达系统、同轴连接

5、器等分析。第5页/共85页5.流体动力学分析 定常/非定常分析 层流/湍流分析 自由对流/强迫对流/混合对流分析 可压缩流/不可压缩流分析 亚音速/跨音速/超音速流动分析 任意拉格郎日欧拉分析（ALE）多组份流动分析（多达6组份）牛顿流与非牛顿流体分析 内流和外流分析 共轭传热及热辐射边界 分布阻尼和风扇模型 移动壁面及自由界面分析 第6页/共85页6.声学分析 定常分析 模态分析 动力响应分析 7.压电分析 稳态、瞬态分析 模态分析 谐响应分析 8.多场耦合分析 热结构 磁热 磁结构 流体热 流体结构 热电 电磁热流体结构 第7页/共85页9.优化设计及设计灵敏度分析 单一物理场优化 耦合场

6、优化 10.二次开发功能 参数设计语言 用户可编程特性 用户自定义界面语言 外部命令 11.ANSYS土木工程专用包 ANSYS的土木工程专用包ANSYS/CivilFEM用来研究钢结构、钢筋混凝土及岩土结构的特性，如房屋建筑、桥梁、大坝、硐室与隧道、地下建筑物等的受力、变形、稳定性及地震响应等情况，从力学计算、组合分析及规范验算与设计提出了全面的解决方案，为建筑及岩土工程师提供了功能强大且方便易用的分析手段。第8页/共85页二、主要模块简介ANSYS/ANSYS/FLOTRANFLOTRANANSYS/ANSYS/EmagEmagANSYS/ANSYS/StructuralStructura

7、lANSYS/ANSYS/MultiphysicsMultiphysicsANSYS/ANSYS/LS-DYNALS-DYNAANSYS/ANSYS/MechanicalMechanicalANSYS/ANSYS/LinearPlusLinearPlusANSYS/ANSYS/ThermalThermal第9页/共85页1.ANSYS/Mechanical：该模块提供了范围广泛的工程设计分析与优化功能，这些功能包括完整的结构、热、压电及声学分析。是一个功能强大的设计校验工具，可用来确定位移、应力、作用力、温度、压力分布以及其它重要的设计标准。2.ANSYS/Structural：通过利用其先进

8、的非线性功能，该模块可进行高目标的结构分析，具体包括：几何非线性、材料非线性、单元非线性及屈曲分析。该模块可以使用户精确模拟大型复杂结构的性能。3.ANSYS/Linear plus：该模块是从ANSYS/Structural派生出来的，一个线性结构分析选项，可用于线性的静态、动态及屈曲分析，非线性分析仅包括间隙元和板/梁大变形分析。第10页/共85页4.ANSYS/Thermal：该模块同样是从ANSYS/Mechanical中派生出来的，是一个可单独运行的热分析程序，可用于稳态及瞬态热分析。5.ANSYS/Flotran：该程序是个灵活的CFD软件，可求解各种流体流动问题，具体包括：层流、

9、紊流、可压缩流及不可压缩流等。通过与ANSYS/Mechanical耦合，ANSYS/FLOTRAN是唯一一个具有设计优化能力的CFD软件，并且能提供复杂的多物理场功能。6.ANSYS/Emag：该程序是一个独立的电磁分析软件包，可模拟电磁场、静电学、电路及电流传导分析。当该程序与其它ANSYS模块联合使用时，则具有了多物理场分析功能，能够研究流场、电磁场及结构力学间的相互影响。7.ANSYS/Preppost：该模块为用户在前处理阶段提供了强大的功能，使用户能够便捷地建立有限元模型。其后处理器能够使用户检查所有ANSYS分析的计算结果。第11页/共85页8.ANSYS/ED：该模块是一个功能

10、完整的设计模拟程序，它拥有ANSYS隐式产品的全部功能，只是解题规模受到了限制（目前节点数1000）。该软件可独立运行，是理想的培训教学软件。9.ANSYS/LS-DYNA：该程序是一个显示求解软件，可解决高度非线性结构动力问题。该程序可模拟板料成形、碰撞分析、涉及大变形的冲击、非线性材料性能以及多物体接触分析，它可以加入第一类软件包中运行，也可以单独运行。PrepPost：该程序具有所有的ANSYS/LS-DYNA的前后处理功能，具体包括：实体建模、网格剖分、加载、边界条件、等值线显示、计算结果评价以及动画，但没有求解功能。：该模块是一个功能完整的设计模拟程序，它拥有ANSYS隐式产品的全部

11、功能，只是解题规模受到了限制（目前节点数16000和32000两种）。该软件可独立运行，适用与高校进行教学或科研。第12页/共85页：该模块是ANSYS的低端产品，适用与设计工程师在产品概念设计初期对产品进行基本分析，以检验设计的合理性。其分析功能包括：线性静力分析、模态分析、基本热分析、基本热力耦合分析、拓扑优化。其他功能有：CAD模型读取器、自动生成分析报告、自动生成ANSYS数据库文件、自动生成ANSYS分析模板。产品详细分类：DesignSpace for MDT DesignSpace for SolidWorks Standalone DesignSpace：（支持的CAD模型有：

12、Pro/E、UG、SAT、Parasoild）：ANSYS与CAD软件的接口产品。可以将CAD模型数据或国际标准格式CAD模型数据直接读入并进行任意ANSYS支持的分析。目前支持的CAD软件有：Pro/E、UG、CADDS。支持的国际标准格式有：IGES、SAT、Parasolid.第13页/共85页第二章第二章 ANSYS基本使用方法基本使用方法第14页/共85页一、典型分析过程1.前处理创建有限元模型 1）单元属性定义（单元类型、实常数、材料属性）2）创建或读入几何实体模型 3）有限元网格划分 4）施加约束条件、载荷条件2.施加载荷进行求解 1）定义分析选项和求解控制 2）定义载荷及载荷步

13、选项 2）求解 solve3.后处理 1）查看分析结果 2）检验结果第15页/共85页ANSYS的分析方法(续)1.建立有限元模型3.查看结果2.施加载荷求解主菜单分析的三个主要步骤可在主菜单中得到明确体现.Objective2-2.ANSYS分析步骤在GUI中的体现.第16页/共85页ANSYS的分析方法(续)ANSYS GUI中的功能排列按照一种动宾结构，以动词开始（如Create),随后是一个名词(如Circle).菜单的排列，按照由前到后、由简单到复杂的顺序，与典型分析的顺序相同.第17页/共85页二、ANSYS文件及工作文件名一些特殊的文件 数据库文件二进制Log 文件文本结果文件二

14、进制图形文件二进制 ANSYS的数据库，是指在前处理、求解及后处理过程中，ANSYS保存在内存中的数据。数据库既存储输入的数据，也存储结果数据:输入数据-必须输入的信息(模型尺寸、材料属性、载荷等).结果数据-ANSYS计算的数值(位移、应力、应变、温度等).第18页/共85页启动ANSYSObjective1-1.启动ANSYS软件.注意:以下过程默认为在 Windows NT环境下.要启动ANSYS:Windows NT 屏幕:Start Programs ANSYS 5.5 1.选择 Interactive.Procedure1.2.3.第19页/共85页启动ANSYS(续)3.选择AN

15、SYS的工作目录，ANSYS所有生成的文件都将写在此目录下。缺省为上次运行定义的目录。4.如果配置了3D显卡则选择3D.5.设定初始工作文件名，缺省为上次运行定义的工作文件名，第一次运行缺省为 file.7.选择 Run to 运行ANSYS.Note:还可以设定其他的交互选项，但通常以上几项通常要设置好.如果在第1步中选择 Run Interactive Now，将读取上一次的设置，跳过此窗口，直接运行ANSYS,2.选择ANSYS产品.6.设定ANSYS工作空间及数据库大小(参考ANSYS安装及配置手册).第20页/共85页启动ANSYS(续)当显示出这六个窗口后，就可以使用ANSYS了.

16、第21页/共85页ANSYS窗口Objective1-2.ANSYS GUI中六个窗口的总体功能输入显示提示信息，输入ANSYS命令，所有输入的命令将在此窗口显示。主菜单包含ANSYS的主要功能，分为前处理、求解、后处理等。输出显示软件的文本输出。通常在其他窗口后面，需要查看时可提到前面。应用菜单包含例如文件管理、选择、显示控制、参数设置等功能.工具条将常用的命令制成工具条，方便调用.图形显示由ANSYS创建或传递到ANSYS的图形.第22页/共85页三、前处理实体建模 参数化建模 体素库及布尔运算 拖拉、旋转、拷贝、蒙皮、倒角等 多种自动网格划分工具，自动进行单元形态、求解精度检查及修正 自

17、由/映射网格划分、智能网格划分、自适应网格划分 复杂几何体Sweep映射网格生成 六面体向四面体自动过渡网格：金字塔形 边界层网格划分 在几何模型或FE模型上加载：点载荷、分布载荷、体载荷、函数载荷 可扩展的标准梁截面形状库 第23页/共85页1.实体模型及有限元模型 现今几乎所有的有限元分析模型都用实体模型建模.类似于CAD，ANSYS以数学的方式表达结构的几何形状，用于在里面填充节点和单元，还可以在几何模型边界上方便地施加载荷.但是，几何实体模型并不参与有限元分析.所有施加在几何实体边界上的载荷或约束必须最终传递到有限元模型上（节点或单元上）进行求解.由几何模型创建有限元模型的过程叫作网格

18、划分Meshing几何实体模型有限元模型第24页/共85页ANSYS中的图元(即使想从CAD模型中传输实体模型，也应该知道如何使用ANSYS建模工具修改传入的模型.)下图示意四类图元.体(3D模型)由面围成，代表三维实体.面(表面)由线围成.代表实体表面、平面形状或壳（可以是三维曲面）.线(可以是空间曲线)以关键点为端点，代表物体的边.关键点(位于3D空间)代表物体的角点.AreasVolumeKeypointsLinesAreaObjective3-3.四类实体模型图元,以及它们之间的层次关系.第25页/共85页ANSYS中图元(续)层次关系从最低阶到最高阶，模型图元的层次关系为：关键点（K

19、eypoints）线（Lines）面（Areas）体（Volumes）提示:如果低阶的图元连在高阶图元上，则低阶图元不能删除.KeypointsLinesAreasVolumesIll justchangethis lineLinesKeypointsAreasVolumesOOPs!第26页/共85页2.布尔操作选择图形类型.将弹出 选取菜单(见下页)提示选择图形进行布尔操作.要使用布尔操作:Main Menu:Preprocessor -Modeling-Operate Procedure1.2.3.选择一种布尔操作(例如：Add)第27页/共85页四、加载、求解Objective4-1.

20、列表和分类载荷ANSYS中的载荷可分为:自由度DOF-定义节点的自由度（DOF）值(结构分析_位移、热分析_ 温度、电磁分析_磁势等)集中载荷-点载荷(结构分析_力、热分析_ 热导率、电磁分析_ magnetic current segments)面载荷-作用在表面的分布载荷(结构分析_压力、热分析_热对流、电磁分析_magnetic Maxwell surfaces等)体积载荷-作用在体积或场域内(热分析_ 体积膨胀、内生成热、电磁分析_ magnetic current density等)惯性载荷-结构质量或惯性引起的载荷(重力、角速度等)第28页/共85页加载Objective4-2a.

21、加载.可在实体模型或 FEA 模型(节点和单元)上加载.在关键点处约束实体模型沿线均布的压力在关键点加集中力在节点处约束FEA 模型沿单元边界均布的压力在节点加集中力第29页/共85页加载(续)+几何模型加载独立于有限元网格.重新划分网格或局部网格修改不影响载荷.+加载的操作更加容易,尤其是在图形中直接拾取时.直接在实体模型加载的优点:Guidelines第30页/共85页加载(续)无论采取何种加载方式,ANSYS求解前都将载荷转化到有限元模型.因此,加载到实体的载荷将自动转化到 其所属的节点或单元上。实体模型加载到实体的载荷自动转化到其所属的节点或单元上FEA 模型沿线均布的压力 均布压力转

22、化到以线为边界的各单元上第31页/共85页加载(续)输入一个压力值即为 均布载荷,两个数值 定义坡度压力说明：压力数值为正表示其方向指向表面Main Menu:Solution -Loads-Apply Pressure On Lines加载面力载荷拾取Line第32页/共85页加载(续)VALI=500VALI=500VALJ=1000VALI=1000VALJ=500500L3500L31000500L31000500坡度压力载荷沿起始关键点(I)线性变化到第二个关键点(J)。如果加载后坡度的方向相反,将两个压力数值颠倒即可。加载面力载荷（续）第33页/共85页加载(续)轴对称载荷可加载到

23、具有对称轴的3-D 结构上。3-D 轴对称结构可用一2-D 轴对称模型描述。加载轴对称载荷10”直径5”半径轴对称模型3-D 结构对称轴第34页/共85页加载(续)加载 轴对称载荷,注意以下方面:载荷数值(包括输出的反力)基于360度转角的3-D结构。在右图中，轴对称模型中的载荷是3-D结构均布面力载荷的总量。Total Force=2pr=47,124 lb.准则3-D 结构2-D 有限元模型Axis of symmetry第35页/共85页加载(续)在关键点加载位移约束:加载约束载荷Main Menu:Solution -Loads-Apply-Structural-Displacemen

24、t On Keypoints+procedure1.2.3.Expansion option 可使相同的载荷加在位于两关键点连线的所有节点上拾取keypoints例要固定一边，只要拾取关键点6、7,并设置 all DOFs=0 和 KEXPND=yes.K6K7第36页/共85页加载(续)加载约束载荷（续）在线和面上加载位移约束:Main Menu:Solution -Loads-Apply-Structural-Displacement On Lines+OR On Areas+步骤1.2.3.拾取 lines 拾取areas第37页/共85页4-1.求解结果文件结果数据数据库求解器结果输入

25、数据第38页/共85页求解时模型是否准备就绪?在求解初始化前，应进行分析数据检查，包括下面内容:统一的单位单元类型和选项材料性质参数考虑惯性时应输入材料密度热应力分析时应输入材料的热膨胀系数实常数(单元特性)单元实常数和材料类型的设置实体模型的质量特性(Preprocessor Operate Calc Geom Items)模型中不应存在的缝隙壳单元的法向节点坐标系集中、体积载荷面力方向温度场的分布和范围热膨胀分析的参考温度(与 ALPX 材料特性协调?)第39页/共85页求解过程:1.求解前保存数据库2.将Output 窗口提到最前面观看求解信息3.Main Menu:Solution -

26、Solve-Current LS.进行求解Objective4-3.描述求解过程第40页/共85页进行求解(续)没有获得结果的原因是什么?往往是求解输入的模型不完整或存在错误，典型原因有:约束不够!(通常出现的问题)。当模型中有非线性单元(如缝隙 gaps、滑块sliders、铰hinges、索cables等），整体或部分结构出现崩溃或“松脱”。材料性质参数有负值,如密度或瞬态热分析时的比热值。未约束铰接结构，如两个水平运动的梁单元在竖直方向没有约束。屈曲-当应力刚化效应为负（压）时，在载荷作用下整个结构刚度弱化。如果刚度减小到零或更小时，求解存在奇异性，因为整个结构已发生屈曲。第41页/共8

27、5页五、后处理计算报告自动生成及定制工具：自动生成符合要求格式的计算报告 结果显示菜单：图形显示、抓图、结果列表 图形：云图、等值线、矢量显示、粒子流迹显示、切片、透明及半透明显示、纹理 各种结果动画显示，可独立保存及重放 3D图形注注释功能 直接生成BMP、JPG、VRML、WMF、EMF、PNG、PS、TIFF、HPGL等格式的图形 计算结果排序、检索、列表及再组合 钢筋混凝土单元可显示单元内的钢筋、开裂情况以及压碎部位 梁、管、板、复合材料单元及结果按实际形状显示，显示横截面结果；显示梁单元弯矩图 显示优化灵敏度及优化变量曲线 提供对计算结果的加、减、积分、微分等计算 显示沿任意路径的结

28、果曲线，并可进行沿路径的数学计算第42页/共85页ANSYS 有两个后处理器:通用后处理器(即“POST1”)只能观看整个模型在某一时刻的结果(如：结果的照相“snapshot”).时间历程后处理器(即“POST26”)可观看模型在不同时间的结果。但此后处理器只能用于处理瞬态和/或动力分析结果。Objective第43页/共85页静力分析结果后处理的步骤主要包括:1.绘变形图2.变形动画3.支反力列表4.应力等值线图5.网格密度检查GuidelinesObjective介绍静力分析结果后处理的五个步骤第44页/共85页第三章 ANSYS补充说明一、坐标系1.工作平面坐标系wpcs：类似于绘图图

29、板，缺省时总与总体坐标系重合，能以网格捕捉形式显示，并可相对当前激活总体坐标系移动或旋转，其编号永远为“4”2.总体坐标系 global cs：包括三种形式 总体直角坐标系（x,y,z)编号为“0”总体柱坐标系 （r,z)编号为“1”总体球坐标系 （r,)编号为“2”3.局部坐标系 local cs：局部坐标系是在任意位置的用户定义坐标系，即不一定与总体坐标系平行或重合，可以是任意方向，编号为大于等于“11”第45页/共85页FEM坐标系4.节点坐标系 node cs：所有的力及其他方向的与节点相关的载荷都是在节点坐标系下进行的，例如力的方向等只与节点坐标系相关 节点坐标系上可以输入力和力矩；

30、位移约束；耦合及约束过程5.单元坐标系 element cs：即材料坐标系，例如弹性模量在材料为各向异性时每一方向将不同，此时则根据单元坐标系输入不同方向的 E 6.结果坐标系：结果的输出形式位移，支反力，力矩等都是与结果坐标系相关的，结果坐标系即当前激活坐标系，同节点坐标系一样，二者可以是任何一种当前激活坐标系 第46页/共85页二、CAD模型建模原则应考虑多少细节：如倒角和孔处，对分析无用时可忽略，但对分析目标有用，而且此处将会出现最大应力则不能忽略是否具有对称性：包括轴，旋转，平面或镜面，重复或平移对称等。但下列因素必须对称-几何形状；材料属性；载荷工况。此时可取一部分分析，而后叠加即可

31、应力奇异：指在有限元模型中那些应力值无限大的点处，如点载荷的集中力和力矩作用处；孤立的约束点；尖角处等。建模时最好避免之第47页/共85页三、网格划分器自由式 free：对复杂的拓扑结构无限制，形状不定映射式 mapped：拓扑结构有限制，只适用规则的体 形状，如四，六面体等，可通过global set 进行密度设置扫略 sweep：适用于柱体形状，同mapped一样可控制密度Smart size：智能尺寸是根据几何模型的形状，确定网格密度，适于free划分，可通过滑杆确定网格密度第48页/共85页网格划分原则网格划分的单元形状四方和六方的没有可比性Sweep扫略网格须上下面即对应面完全一致能

32、用mapped,sweep划分网格最好先用之，不行再用自由式free网格划分最好按线，面，体的顺序分配单元属性千万不能分配错误面尽量用四边形的网格，体尽量用六面体的网格关心应力结果的区域须进行详细网格划分仅关心位移结果的地方网格可以粗糙些第49页/共85页四、ANSYS求解器类型用于求解表征结构自由度的线性方程组 直接消去求解器波前求解器：最稳定，速度慢，小内存时用sparse求解器：速度快，非线性最适合 迭代求解器PCG：预条件共轭梯度求解器ICCG：不完全的乔里斯基共轭梯度求解器JCG：雅可比共轭梯度求解器 位置在求解器/solu中的求解选项analysis options，包括求解精度公

33、差更改第50页/共85页五、如何加快计算速度在大规模结构计算中，计算速度是一个非常重要的问题，下面就如何提高计算速度作一些建议：充分利用ANSYS MAP分网和SWEEP分网技术，尽可能获得六面体网格，这一方面减小解题规模，另一方面提高计算精度在生成四面体网格时，用四面体单元而不要用退化的四面体单元。比如95号单元有20节点，可以退化为10节点四面体单元，而92号单元为10节点单元，在此情况下用92号单元将优于95号单元选择正确的求解器。对大规模问题，建议采用PCG法。此法比波前法计算速度要快10倍以上（前提是您的计算机内存较大）。对于工程问题，可将ANSYS缺省的求解精度从1E-8改为1E-

34、4或1E-5即可第51页/共85页六、Animate菜单 Mode shape：变形模态系列Deformed shape：结构变形动画Deformed result：结构变形等值线Over time：随时间变化的变形等值线动画Over result：某一子步范围结果的顺序等值线动画Q-slice contours：变形等值线切片云图动画Q-slice vectors：变形等值线切片云图动画Isosurfaces：变形等势面云图动画Partice flow：粒子流动或带电粒子运动的动画系列第52页/共85页第四章第四章 动力学分析动力学分析第53页/共85页ANSYS动力学分析是用来确定惯性力和

35、阻尼力不可忽略时的系统动力学特性，研究固有频率，振动，减振及瞬态特性动力学控制方程：MU+CU+KU=f(t)其中 U U U 为节点位移，速度，加速度 M为质量矩阵 C 为阻尼矩阵 K 为刚度矩阵模态分析即f(t)=0的解谐响应分析的f(t),u(t)都为谐函数，如xsint瞬态动力学的f(t)为时间历程载荷第54页/共85页动力学建模原则必须定义密度和弹性模量单位制要严格统一，如使用英制单位，要定义质量密度而不是重力密度静力学关于形状和网格的东西，动力学一样遵循关心应力结果的区域须进行详细网格划分，仅关心位移结果的区域网格可以粗糙些非线性问题在完全瞬态动力学分析中允许使用，在所有其他动力学

36、分析类型中，非线性将被忽略，也就是说最初的非线性问题将一直保持不变第55页/共85页瞬态动力学时间积分步长T设置原则：T 即两个时间点间的时间增量，它决定求解的精确度，必须采用相应的值才能得到分析现象。通常在每个循环子步中，T 至少有20个时刻点应是足够的 即 T=1/20(f)f是所关心的最高响应频率，而施加阶梯载荷时，为紧紧跟随载荷的阶跃变化，T也许要小到和1/180f相近接触分析中T=1/30(fc)fc为接触频率 fc=1/2(k/m)1/2 m 为有效质量 k 为间隙刚度第56页/共85页 各分析类型的简要叙述模态分析：分析有/无阻尼系统或液固耦合自由振动，哪个方向的刚度最小，就最先

37、出现模态有预应力时，做静力分析和模态分析必须打开预应力开关，即读入预应力矩阵 位置在求解器Analysis options对话框中prestress on 开关第57页/共85页谐响应分析定义：在一个或多个同频率的正弦（简谐）载荷作用下，确定系统动力响应特性的技术输入特性：1.已知大小和频率的谐波载荷（力，压强，强迫位移等）2.同一频率的多种载荷，具有相同或不相 相同的相位输出特性：1.每一个自由度上的谐位移通常和施加的 载荷不同相 2.其他多种导出量如应力和应变等第58页/共85页谐波载荷描述公式 Fsin(t+)其中载荷值大小代表振幅 F；相角是指二个或二个以上谐波载荷间的相位差，即力为幅

38、值+相位输入，但ANSYS始终显示为实部+虚部瞬态动力学分析：把静力学改为瞬态动力学并采用多载荷步求解即为瞬态动力学分析。其关键细节为：瞬态动力学必须指定初始条件，如初始位移，速度，加速度等如有阻尼必须输入阻尼时间积分步长T 要足够小第59页/共85页区分分析类型原则如果在相对较长时间内载荷是一个常数，则选择静力分析，否则为动态分析如果激励频率小于结构最低阶固有频率的1/3，则可以进行静力分析线性分析是假设忽略载荷对结构刚度变化的影响。典型特征是小变形，弹性范围内的应变和应力，无诸如两物体接触或分离时的刚度突变等，即应力及应变为线性变化如加载引起结构刚度的显著变化必须进行非线性分析，典型因素有

39、应变超过弹性范围（塑性）；大变形；两体间接触等第60页/共85页解决应力奇异方法：在应力奇异点处网格越细化，应力值也随之增加而不收敛，如应力奇异发生在对结果感兴趣区域较近处，可用如下方法解决在尖角处增加倒角重分析代替节点力载荷为等效压力载荷“散布”位移约束至一个节点集第61页/共85页Animate菜单 Mode shape：变形模态系列Deformed shape：结构变形动画Deformed result：结构变形等值线Over time：随时间变化的变形等值线动画Over result：某一子步范围结果的顺序等值线动画Q-slice contours：变形等值线切片云图动画Q-slice

40、 vectors：变形等值线切片云图动画Isosurfaces：变形等势面云图动画Partice flow：粒子流动或带电粒子运动的动画系列第62页/共85页第五章第五章 非线性分析非线性分析第63页/共85页非线性的一些定义非线性具有bauschinger效应，即拉伸屈服后，再拉伸应力反而减小应力偏量：一般应力可分解为静水压应力和应力偏量率相关性：非线性塑性应变大小可能是加载荷快慢的函数，如果应变不须考虑时间，则称为率相关性ANSYS的非线性主要为塑性，即施加力后将产生永久变形ANSYS非线性包括材料非线性；接触非线性；几何非线性，将在后面叙述第64页/共85页非线性应力应变曲线如图一，AN

41、SYS中其屈服点和极限点将归结为一点。理想塑性曲线如图二；应变强化曲线如图三。当非线性材料施加后，可用命令 prep7material propsData tablesGraph观看曲线形式第65页/共85页非线性塑性强化准则：总体方向为各向同性强化及随动强化（即轴线漂移）ANSYS塑性强化选项：1.双线性随动强化 BKIN 双条线表示应力应变曲线，输入材料时须输入弹性模量E 和屈服应力2.双线性各向同性强化 BISO3.多线性随动强化 MKIN4.多线性随动强化 KINH如右图3和4项适于多线段应力-应变分析，应用于金属的小应变分析第66页/共85页5.多线性各向同性强化 MISO6.非线性

42、随动强化 CHAB7.非线性各向同性强化 NLIS8.各向异性 ANISO9.Drucker-prager DP：DP材料适于颗粒状材料其屈服面为圆锥形，须输入粘性值C；内部摩擦角和膨胀角r第67页/共85页10.ANAND 模型 ANAND：用于金属在热加工状态的大应变响应。注意此时材料温度假设为高于熔点的一半以上为材料非线性塑性一些选项，材料非线性可以认为是加载的时间量（蠕变等）；环境状况（如温度）；加载历史等因素影响到材料的应力应变性质变化蠕变：在常载荷作用下，不可恢复的应变随时间变化持续增加，结构将会持续变形应力松弛：如果结构承受常位移，应力将会随时间而减小应力刚化：面内应力对面外刚度

43、的影响，二者间的耦合即称应力刚化第68页/共85页 非线性单元options综合选项不协调模式（特殊形函数）：solid45的缺省选项，适用于弯曲变形选择缩减积分（BBar）：接近不可压缩材料，为体积变形统一缩减积分（URI）：接近不可压缩材料，适用于弯曲变形混合U-R形式：不可压缩或接近不可压缩材料第69页/共85页推荐使用的单元 对率无关的塑性推荐用：忽略弯曲的体积变形用plane182;solid185;但需选B-Bar选项对于小应变分析用不协调模式plane42;solid45对于大应变分析用plane182和solid185;solid95的URI选项,以上更适用于大模型,也可用vi

44、sco106;visco107;visco108单元对梁单元推荐用beam188;beam189对壳单元推荐用shell181 第70页/共85页非线性求解选项使用求解控制,不带自适应下降的全为Newtonraphson选项，推荐用。位置在求解器analysis options对话框里，多数用缺省项自动选择塑性要求时间的载荷增量T 必须足够小如果应力应变是非线性变化的，T 如果大，则结果就可能不收敛如果要应力应变曲线平滑，须打开预测器，减少迭代数第71页/共85页网格划分注意事项缩减积分单元需要更细的网格。因塑性计算发生在单元的积分点处，因此，积分点处网格密度非常重要，否则不能捕捉到进入塑性的

45、范围对弯曲分析在厚度上要有足够的单元，可从粗-细逐渐过渡第72页/共85页几何非线性用于分析大位移，大应变，大转角状态。其特性为：大位移指如果一个单元变化，其单元矩阵也变化，而且单元面的应力将影响刚度。大应变指应变不在小，而是有限或无限大的。以上这些ANSYS将根据模型自动判别第73页/共85页几何非线性求解选项求解前须将soluAnalysis options对话框中第一项大变形Large deform effects 设为on，以激活支持大应变打开大变形时，须包含应力刚化项，即stress stiffness or prestress 设为“stress stiff on”非线性分析必须打

46、开自动时间步长，并给出最大最小范围，以保证收敛（AUTOTS）线性搜索选项（LNSRCH）对收敛振荡问题有帮助。位置在solunonlinearline search一般缺省为prog chosen第74页/共85页非线性载荷方向在变形中多数情况下将保持一致的方向大转角时力将随单元改变方向加速度与集中力保持初始方向，永远不变表面压力载荷随单元变化，总是垂直变形单元的表面，施加的常压力，总压力将随表面积的改变而改变第75页/共85页几何非线性注意事项：注意载荷方向 在大应变分析中预测网格扭曲，要划分适当的网格 避免过分约束边界处的变形，以免产生应力奇异 避免使用带中间节点的单元 用适当的单元类型

47、和积分准则以解决网格自锁 时间积分步长大小合适，大转角时每一个子步须小于5度或10度 大转角分析关闭预测solunonlinearpredictor-off梁，壳单元要有足够的网格密度，没有一个单元可承受30度的弯曲度，超过30度就用一个单元描述第76页/共85页接触非线性接触分类：刚体对柔体-刚度差别较大，如金属成型柔体对柔体-表面刚度差不多，如螺栓，法兰联接等接触协调条件 两个表面须建立关系，防止表面相互渗透。即建立接触面和目标面。防止穿透标准为判断接触面的节点是否进入目标面积分点，如进入目标范围则称之接触第77页/共85页接触单元分类节点对节点：指接触的最终位置是已知的节点对面：接触位置

48、不定，并允许大滑动面对面：接触位置不定，并允许大滑动，最好用之。过程是用接触面和目标面形成接触对，但目标面ANSYS要求为刚性，称为刚性目标面第78页/共85页接触分析步骤创建有限元模型 /prep7creat设置单元选项和实常数创建目标表面单元 meshtool创建接触表面单元指定接触面和目标面/prep7creatcontact wizard施加边界条件/solu定义求解选项和载荷步 求解 solve观察结果/post1第79页/共85页选择目标面方针平面，凹面与凸面接触，平面，凹面为目标面如两表面网格划分一个粗，一个细，则网格粗的为目标面刚度大的为目标面大表面的为目标面如有一面划分成高次

49、单元，另一面为低次单元，则表面划分低次单元的为目标面第80页/共85页接触单元（conta174)关键点选项：增强的拉格朗日法penal and lagrange k(2)=0是缺省选项，推荐一般使用。对罚刚度不太敏感，但也要求给出一个侵入容差可用罚函数法 penalty method k(2)=1选项，推荐使用单元非常扭曲，大摩擦系数或增强的拉格朗日法收敛行为不好的问题罚刚度即法向刚度 normal penalty stiffness，无滑移时，摩擦问题ANSYS自动处理ANSYS面对面单元默认用k(4)高斯积分点 Gauss point来判断第81页/共85页自接触问题，用k(6)项不对称

50、接触Unsymmetrical更有效如在静态分析中，开始不连接的物体，建立接触前产生刚体位移，则可在刚体上加一些K非常小的软弹簧即可接触状态分类：（单元输出结果）远区域开接触 ID号为0 近区域开接触 ID号为1 滑动闭接触 ID号为2 粘着闭接触 ID号为3第82页/共85页接触方式k(12)contact surface behavior标准的方式 Standard粗糙接触方式：无滑动，摩擦系数无穷大 Rough不分离接触方式：不能脱开但允许滑动，法向无间隙，No separation(sliding)绑定方式：不能脱开移动，即法向无间隙 Bonded“总是”不分离方式：No separa