材料成形过程数值模拟.ppt

1,材料成形过程数值模拟,2,课程简介,课时数：40课时学分数：2.5先修课程线性代数、工程力学、传质传热学、材料科学基础、材料成型原理、热处理原理及工艺、铸造工艺及设备、焊接工艺及设备等,3,目的：系统掌握材料成形过程数值模拟的基本原理和一般步骤了解材料成形时各加工方式的特点、基本规律、数值分析原理与技术要点培养学习者分析解决实际工程问题的能力。要求：掌握材料成形数值模拟的基本知识和理论掌握各种材料成形过程的相关理论、数值方法、实现过程和实际应用了解材料加工过程中建模与数值分析技术现状和发展趋势,课程目的和要求,4,第一章 绪论,基本概念CAE将一个成形过程（或过程的某一方面）定义为由一组控制方程加上边界条件构成的定解问题，从而获得对成形过程的定量认识。材料成形包括:金属铸造成形、塑性成形、连接成形高分子材料黏流态注射成形,5,1.2 材料成形数值模拟的工程意义及应用现状,工程意义制定和优化材料成形方案与模具设计方案解决工模具调试或产品试成形过程中的技术问题解决成形制品批量生产中的质量控制问题应用现状铸造成形、塑性成形、焊接成形、黏流态成形,6,1.3 材料成形数值模拟的发展趋势,（1）模拟分析由宏观进入微观（2）加大多物理场的耦合分析（3）拓宽数值模拟在特种成形中的应用（4）强化基础性研究（5）关注反向模拟技术应用,7,第二章 有限元与有限差分法基础,CAE的工具：有限元法（FEM）、有限差分法（FDM）、边界元法（BEM）、有限体积法（FVM）、无网格法等等在材料成形的CAE中主要使用的是有限元法和有限差分法,8,“有限元法”的基本思想早在20世纪40年代初期就有人提出，但真正用于工程中则是电子计算机出现以后。“有限元法”这一名称是1960年美国的克拉夫（Clough，R.W.）在一篇题为“平面应力分析的有限元法”论文中首先使用。此后，有限元法的应用得到蓬勃发展。到20世纪80年代初期国际上较大型的结构分析有限元通用程序多达几百种，从而为工程应用提供了方便条件。由于有限元通用程序使用方便，计算精度高，其计算结果已成为各类工业产品设计和性能分析的可靠依据。,9,有限元法最初用于飞机结构的强度设计，由于它在理论上的通用性，因而它可用于解决工程中的许多问题。目前，它可以解决几乎所有的连续介质和场的问题，包括热传导、电磁场、流体动力学、地质力学、原子工程和生物医学等方面的问题。机械设计中，从齿轮、轴、轴承等通用零部件到机床、汽车、飞机等复杂结构的应力和变形分析（包括热应力和热变形分析）。有限元法不仅可以解决工程中的线性问题、非线性问题，而且对于各种不同性质的固体材料，如各向同性和各向异性材料，粘弹性和粘塑性材料以及流体均能求解；对于工程中最有普遍意义的非稳态问题也能求解。,10,2.1 有限元法基础,基本思想：将一个连续求解域（对象）离散（剖分）成有限个形状简单的子域（单元）利用有限个节点将各子域连接起来在给定的初始条件和边界条件下进行综合计算求解，从而获得对复杂工程问题的近似数值解,11,为什么要离散？,1.无法得到复杂实际问题的解析解2.将域划分成一些微小而形状规则的单元后，便于在一个单元内得到近似解3.域中所有单元的解可视为该复杂问题的近似解,12,有限元分析的过程,1.连续体离散化 2.单元分析 3.整体分析 4.确定约束条件 5.方程求解 6.结果分析与讨论,13,1.连续体离散化,连续体：是指所求解的对象（如物体或结构）。离散化（划分网格或网络化）：是将所求解的对象划分为有限个具有规则形状的微小块体，把每个微小块体称为单元，相邻两个单元之间只通过若干点互相连接，每个连接点称为节点。相邻单元只在节点处连接，载荷也只通过节点在各单元之间传递，这些有限个单元的集合体，即原来的连续体。*单元划分后，给每个单元及节点进行编号；*选定坐标系，计算各个节点坐标；*确定各个单元的形态和性态参数以及边界条件等。,14,单元的划分基本上是任意的，一个结构体可以有多种划分结果。但应遵循以下划分原则：(1)分析清楚所讨论对象的性质，例如，是桁架结构还是结构物，是平面问题还是空间问题等等。(2)单元的几何形状取决于结构特点和受力情况，单元的几何尺寸(大小)要按照要求确定。一般来说，单元几何形体各边的长度比不能相差太大。(3)有限元模型的网格划分越密，其计算结果越精确，但计算工作量就越大。因此，在保证计算精度的前提下，单元网格数量应尽量少。(4)在进行网格疏密布局时，应力集中或变形较大的部位，单元网格应取小一些，网格应划分得密一些，而其他部分则可疏一些。,15,(5)在设计对象的厚度或者弹性系数有突变的情况下，应该取相应的突变线作为网格的边界线；(6)相邻单元的边界必须相容，不能从一单元的边或者面的内部产生另一个单元的顶点。(7)网格划分后，要将全部单元和节点按顺序编号，不允许有错漏或者重复。(8)划分的单元集合成整体后，应精确逼近原设计对象。原设计对象的各个顶点都应该取成单元的顶点。所有网格的表面顶点都应该在原设计对象的表面上。所有原设计对象的边和面都应被单元的边和面所逼近。,16,有限元分析模型图例,将悬臂梁划分为许多三角形单元三角形单元的三个顶点都是节点载荷直接施加在节点上,悬臂梁及其有限元模型,17,2.单元分析,连续体离散化后，即可对单元体进行特性分析，简称为单元分析。单元分析工作主要有两项：(1)选择单元位移模式(位移函数)用节点位移来表示单元体内任一点的位移、应变和应力，就需搞清各单元中的位移分布。一般是假定单元位移是坐标的某种简单函数，用其模拟内位移的分布规律，这种函数就称为位移模式或位移函数。通常采用的函数形式多为多项式。根据所选定的位移模式，就可以导出用节点位移来表示单元体内任一点位移的关系式。,18,2.单元分析(2),(2)分析单元的特性，建立单元刚度矩阵 进行单元力学特性分析，将作用在单元上的所有力（表面力、体积力、集中力）等效地移置为节点载荷；采用有关的力学原理建立单元的平衡方程，求得单元内节点位移与节点力之间的关系矩阵单元刚度矩阵。,19,3.整体分析,把各个单元的刚度矩阵集成为总体刚度矩阵，以及将各单元的节点力向量集成总的力向量，求得整体平衡方程。集成过程所依据的原理是节点变形协调条件和平衡条件。,20,4.确定约束条件,由上述所形成的整体平衡方程是一组线性代数方程，在求解之前，必修根据具体情况分析，确定求解对象问题的边界约束条件，并对这些方程进行适当修正。,21,5.有限元方程求解,应用有限元法求解机械结构应力类问题时，根据未知量和分析有三种基本解法：,位移法 力法 混合法,22,(1)位移法以节点位移作为基本未知量，通过选择适当的位移函数，进行单元的力学特性分析。在节点处建立单元刚度方程，再组合成整体刚度矩阵，求解出节点位移后，进而由节点位移求解出应力。位移法优点是比较简单，规律性强，易于编写计算机程序。所以得到广泛应用，其缺点是精度稍低。(2)力法以节点力作为基本未知量，在节点处建立位移连续方程，求解出节点力后，再求解节点位移和单元应力。力法的特点是计算精度高。(3)混合法取一部分节点位移和一部分节点力作为基本未知量，建立平衡方程进行求解。,23,单元特性的推导方法,单元刚度矩阵的推导是有限元分析的基本步骤之一。目前，建立单元刚度矩阵的方法主要有以下四种：,直接刚度法 虚功原理法 能量变分法 加权残数法,24,1.直接刚度法 直接刚度法是直接应用物理概念来建立单元的有限元方程和分析单元特性的一种方法。这一方法仅能适用于简单形状的单元，如梁单元。但它可以帮助理解有限元法的物理概念。,图1所示是xoy平面中的一简支梁简图，现以它为例，来说明用直接刚度法建立单元刚度矩阵的思想和过程。,图1平面简支梁元及其计算模型,25,梁在横向外载荷（可以是集中力或分布力或力矩等）作用下产生弯曲变形，在水平载荷作用下产生线位移。对于该平面简支梁问题：梁上任一点受有三个力的作用：水平力Fx，剪切力Fy，弯矩Mz。相应的位移为：水平线位移u，挠度v，转角 z。,由上图可见：,水平线位移和水平力向右为正，挠度和剪切力向上为正，转角和弯矩逆时针方向为正。,通常规定：,26,为使问题简化，可把图示的梁看作是一个梁单元。如图1所示，当令左支承点为节点 i，右支承点为节点 j 时，则该单元的节点位移和节点力可以分别表示为：,称为单元的节点位移列阵。,称为单元的节点力列阵；若 F 为外载荷，则称为载荷列阵。,（1-1）,（1-2）,写成矩阵形式为,q(e)=ui,vi,zi,vj,uj,zjT,ui,vi,zi,vj,uj,zj,F(e)=Fxi,Fyi,Mzi,Fxj,Fyj,MzjT,Fxi,Fyi,Mzi,Fxj,Fyj,Mzj,27,显然，梁的节点力和节点位移是有联系的。在弹性小变形范围内，这种关系是线性的，可用下式表示,或,（1-3b）,（1-3a）,28,上式(1-3b)称为单元有限元方程，或称为单元刚度方程，它代表了单元的载荷与位移之间（或力与变形之间）的联系；式中，K(e)称为单元刚度矩阵，它是单元的特性矩阵。,对于图1所示的平面梁单元问题，利用材料力学中的杆件受力与变形间的关系及叠加原理，可以直接计算出单元刚度矩阵K(e)中的各系数 kst(s,t=i,j)的数值,29,2.虚功原理法,下面以平面问题中的三角形单元为例，说明利用虚功原理法来建立单元刚度矩阵的步骤。如前所述，将一个连续的弹性体分割为一定形状和数量的单元，从而使连续体转换为有限个单元组成的组合体。单元与单元之间仅通过节点连结，除此之外再无其他连结。也就是说，一个单元上的只能通过节点传递到相邻单元。,从分析对象的组合体中任取一个三角形单元：设其编号为 e，三个节点的编号为i、j、m，在定义的坐标系 xoy 中，节点坐标分别为(x j,y j)、(xi,y i)、(xm,ym)，如图2所示。,图2三节点三角形单元,30,由弹性力学平面问题的特点可知，单元每个节点有两个位移分量，即每个单元有6个自由度，相应有6个节点载荷，写成矩阵形式，即,单元节点载荷矩阵：F(e)=Fxi,Fyi,Fxj,Fyj,Fxm,FymT,单元节点位移矩阵：q(e)=ui,vi,uj,vj,um,vmT,图2三节点三角形单元,31,(1)设定位移函数,按照有限元法的基本思想：首先需设定一种函数来近似表达单元内部的实际位移分布，称为位移函数，或位移模式。,三节点三角形单元有6个自由度，可以确定 6个待定系数，故三角形单元的位移函数为,（1-4）,式(1-4)为线性多项式，称为线性位移函数，相应的单元称为线性单元。,u=u(x,y)=1+2x+3yv=v(x,y)=4+5x+6y,32,上式(5-5)也可用矩阵形式表示，即,式中，d为单元内任意点的位移列阵。,（1-5）,33,由于节点 i、j、m 在单元上，它们的位移自然也就满足位移函数式(1-4)。设三个节点的位移值分别为(ui,vi)、(uj,vj)、(um,vm)，将节点位移和节点坐标代入式(1-4)，得,34,（1-6）,式中,（1-7）,由上可知，共有6个方程，可以求出6个待定系数。解方程，求得各待定系数和节点位移之间的表达式为,为三角形单元的面积。其中：,35,（1-8）,将式(1-7)及式(1-8)、式(1-9)代入式(1-6)中，得到,（1-9）,（1-10）,36,式中，矩阵N 称为单元的形函数矩阵；为单元节点位移列阵。其中，为单元的形函数，它们反映单元内位移的 分布形态，是x,y 坐标的连续函数，且有,（1-11）,式(1-10)又可以写成,（1-12）,上式清楚地表示了单元内任意点位移可由节点位移插值求出。,37,(2)利用几何方程由位移函数求应变,根据弹性力学的几何方程，线应变 剪切应变 则应变列阵可以写成,式中，B称为单元应变矩阵，它是仅与单元几何尺寸有关的常量矩阵，即,（1-13）,38,（1-14）,上述方程(1-13)称为单元应变方程，它的意义在于：单元内任意点的应变分量亦可用基本未知量即节点位移分量来表示。,39,(3)利用广义虎克定律求出单元应力方程,根据广义虎克定律，对于平面应力问题,上式（1-15）也可写成,（1-15）,（1-16）,式中，为应力列阵；D 称为弹性力学平面问题的弹性矩阵，并有,40,则有如下单元应力方程,由式(1-18)可求单元内任意点的应力分量，它也可用基本未知量即节点位移分量来表示。,（1-17）,（1-18）,41,(4)由虚功原理求单元刚度矩阵,根据虚功原理，当弹性结构受到外载荷作用处于平衡状态时，在任意给出的微小的虚位移上，外力在虚位移上所做的虚功 AF等于结构内应力在虚应变上所存储的虚变形势能 A，即,设处于平衡状态的弹性结构内任一单元发生一个微小的虚位移，则单元各节点的虚位移 为,（1-20）,（1-21）,（1-19）,则单元内部必定产生相应的虚应变，故单元内任一点的虚应变 为,42,显然，虚应变和虚位移之间关系为,设节点力为,则外力虚功为,（1-24）,（1-22）,（1-23）,单元内的虚变形势能为,43,根据虚功原理,因为,（1-26）,（1-25）,代人式(1-25)，则有,式中，均与坐标 x,y 无关，故可以从积分符号中提出，可得：,44,其中，单元刚度矩阵,（1-27）,式(1-27)称为单元有限元方程，或称单元刚度方程，其中 是单元刚度矩阵。,（1-28）,因为三角形单元是常应变单元，其应变矩阵B、弹性矩阵D均为常量，而，所以式(1-28)可以写成,（1-29）,45,式中，t 为三角形单元的厚度；为三角形单元的面积。,对于图2所示的三角形单元，将D 及B代入式(1-28)，可以得到单元刚度为,（1-30）,式中:K为66阶矩阵，其中每个子矩阵为22阶矩阵，由下式给出,（1-31）,46,按照力学的一般说法，任何一个实际状态的弹性体的总位能是这个系统从实际状态运动到某一参考状态(通常取弹性体外载荷为零时状态为参考状态)时它的所有作用力所做的功。弹性体的总位能 是一个函数的函数，即泛函，位移是泛函的容许函数。,从能量原理考虑，变形弹性体受外力作用处于平衡状态时，在很多可能的变形状态中，使总位能最小的就是弹性体的真正变形，这就是最小位能原理。用变分法求能量泛函的极值方法就是能量变分原理。能量变分原理除了可解机械结构位移场问题以外，还扩展到求解热传导、电磁场、流体力学等连续性问题。,3.能量变分原理法,47,该方法是将假设的场变量的函数(称为试函数)引入问题的控制方程式及边界条件，利用最小二乘法等方法使残差最小，便得到近似的场变量函数形式。该方法的优点是不需要建立要解决问题的泛函式，所以，即使没有泛函表达式也能解题。,4.加权残数法,48,有限元解的收敛性有限元解是近似解近似解是否收敛于真实解、近似解收敛速度、近似解的稳定性近似解的收敛条件:1.完备性准则（必要条件）试探函数（插值函数）的次数（m）不小于场函数的最高可导阶数2.协调性准则（充分条件）试探函数在m-1次连续可导。,49,有限元分析的误差,有限元分析误差,建模误差,计算误差,离散误差,物理离散误差,几何离散误差,边界条件误差,单元形状误差,舍入误差,截断误差,插值函数与真实函数之间的差异1.减小单元特征尺寸，称为h法2.提高插值函数的阶次，称为p法,单元组合体与求解对象几何形状的差异1.网格局部加密2.选用边或面上带有节点的单元,边界条件的复杂性1.准确测定，完善模型2.细分边界网格,单元严重畸变而退化细分局部网格或者控制调整关键区域的网格,数据储存,计算方法、解题性质、解题规模,注意网格的划分 选择合适的解算方法 控制解题的规模,减少运算次数，降低解题规模,选择合适的解算方法，控制解题规模,50,材料成形中的非线性问题,1.材料非线性材料本构方程非线性 弹塑性、刚弹性、刚黏塑性、黏弹塑性2.几何非线性3.边界非线性,51,2.2 有限差分法基础,一种直接将微分问题转变成代数问题的近似数值解法。基本思想数值微分法是把连续的定解区域划分成差分网格，用有限个节点代替原连续求解域。把原方程和定解条件中的微商用差商来近似把原微分方程和定解条件近似地用代数方程组代替，即有限差分方程,52,差分网格通常为矩形在边界不规则或者形状复杂时精度降低,有限元网格,有限差分网格,53,差分概念,自变量x的解析函数 y=f(x)，则有：dx,dy自变量和函数的微分 函数对自变量的一阶导数 函数对自变量的一阶差商,差商,54,差分方向,向前差分向后差分中心差分,55,差商的截断误差,将函数f(x+x)按Taylor级数展开向前向后中心,56,二阶中心差商通常采用向前差商的向后差商截断误差与(x)2 同一数量级一阶向前差商 一阶向后差商 一阶计算精度一阶中心差商 二阶中心差商 二阶计算精度,57,我们在弹性体上，用相隔等间距h而平行于坐标轴的两组平行线织成正方形网格，x=y=h，如图。,设f=f(x,y)为弹性体内的某一个连续函数。该函数在平行于x轴的一根网线上，如在-上,它只随x坐标的改变而变化。在邻近结点处，函数f可展为泰勒级数如下：,58,我们将只考虑离开结点充分近的那些结点，即（x-x0）充分小。于是可不计（x-x0）的三次及更高次幂的各项，则上式简写为：,在结点，x=x0-h,在结点1,x=x0+h，代入(b)得：,59,联立(c),(d),解得差分公式：,同理，在网线-上可得到差分公式,60,从而可导出其它的差分公式如下：,61,相隔2h的两结点处的函数值来表示中间结点处的一阶导数值，可称为中点导数公式。,以相邻三结点处的函数值来表示一个端点处的一阶导数值，可称为端点导数公式。,中点导数公式与端点导数公式相比，精度较高。因为前者反映了结点两边的函数变化，而后者却只反映了结点一边的函数变化。,62,边界元法简介,边界元法（boundary element method）一种结合有限元法和边界积分法发展起来的一种新数值方法只在定义域的边界上划分单元，用满足控制方程的函数去逼近边界条件。适用于应力（薄板）、流体力学、声场、电磁场等问题,63,边界元法基本思想,以微分控制方程的基本解为权函数，利用加权余量法将区域积分转化为边界积分，并结合求解域边界的离散，构建基于边界单元的代数方程组，然后进行计算求解以定义在边界上的边界积分方程为控制方程，通过对边界元插值离散，化为代数方程组求解降低了问题的维数，从而显著降低了自由度数边界的离散比区域的离散方便得多，可用较简单的单元准确地模拟边界形状，最终得到阶数较低的线性代数方程组,64,加权余量法简介,一种直接从所需求解的微分方程及边界条件出发，寻求边值问题近似解的数学方法。从静力发展到了动力、稳定、材料非线性和几何非线性等各方面。在求解域上建立一个试函数 试函数由完备函数集的子集构成。已被采用过的试函数有幂级数、三角级数、样条函数、贝赛尔函数、切比雪夫和勒让德多项式等等。试函数与真实解之间的偏差，即余量（内部和边界）引入权函数，定义消除余量的条件加权余量法就是一种定义近似解与真解之间余量，并设法使其最小的方法。,65,设问题的控制微分方程为:在V域内 在S边界上 式中:L、B分别为微分方程和边界条件中的微分算子；f、g 为与未知函数u无关的已知函数域值；u为问题待求的未知函数。,66,当利用加权余量法求近似解时，首先在求解域上建立一个试函数，一般具有如下形式：,式中：待定系数，也可称为广义坐标；,取自完备函数集的线性无关的基函数。,由于 一 般只是待求函数u的近似解，因此代入后将得不到满足，若记：,在V域内,在S边界上,显然 反映了试函数与真实解之间的偏差，它们分别称做内部和边界余量。,67,若在域V内引入内部权函数，在边界S上引入边界权函数则可建立n个消除余量的条件，一般可表示为：,不同的权函数 和 反映了不同的消除余量的准则。从上式可以得到求解待定系数矩阵C的代数方程组。一经解得待定系数，由式(5.1.3)即可得所需求解边值问题的近似解。,68,由于试函数 的不同，余量 和 可有如下三种情况，依此加权余量法可分为：1内部法试函数满足边界条件，也即 此时消除余量的条件成为:2边界法试函数满足控制方程，也即 此时消除余量的条件为:,69,3混合法 试函数不满足控制方程和边界条件，此时消除余量的条件为:显然，混合法对于试函数的选取最方便，但在相同精度条件下，工作量最大。对内部法和边界法必须使基函数事先满足一定条件，这对复杂结构分析往往有一定困难，但试函数一经建立，其工作量较小。,70,边界元法特点,1 前处理工作量小 2 解算规模小 3 求解奇异性问题时计算精度高 4 在载荷集中和半无限域等问题上有优势相对于有限元法，边界元法发展较慢，,71,有限元法解决应力集中问题,在应力分析中对于应力集中区域必须划分很多的单元，从而增加了求解方程的阶数，计算费用也就随之增加用位移型有限元法求解出的应力的精度低于位移的精度，对于一个比较复杂的问题必须划分很多单元，相应的数据输人量就很大，同时，在输出的大量信息中，又有许多并不是人们所需要的。,72,1.使问题的维数降低一维原为三维空间的可降为二维空间，原为二维空间的问题可降为一维。2.只需将边界离散而需将区域离散化所划分的单元数目远小于有限元，这样它减少了方程组的方程个数和求解问题所需的数据，不但减少了准备工作，而且节约了计算时间。,边界元法与有限元法比较,73,3.由于直接建立在问题控制微分方程和边界条件上的，不需要事先寻找任何泛函。有限元法是以变分问题为基础，如果泛函不存在就难于使用。可以求解经典区域法无法求解的无限域类问题。4.由于边界元法引入基本解，具有解析与离散相结合的特点，因而具有较高的精度。,74,边界元法的缺点,对于非均值和非线性问题求解比较困难用边界元法解非线性问题时，遇到同非线性项相对应的区域积分，这种积分在奇异点附近有强烈的奇异性，使求解遇到困难。,75,有限体积法,基本思想将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积。将待解的微分方程对每一个控制体积积分，便得出一组离散方程。其中的未知数是网格点上的因变量的数值。有限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足 属于加权余量法中的子域法,76,加权余量法子域法(Subdomain Method)此法首先将求解域V划分成n个子域，在每个子 域内令权函数等于1，而在子域之外取权函数为零，也即:,如果在各个子域里分别选取试函数，那么它的求解在形式上将类似于有限元法。,77,有限体积法特点,即使在粗网格情况下，也显示出准确的积分守恒 就离散方法而言，有限体积法可视作有限元法和有限差分法的中间物 有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布；有限元法必须假定值在网格点之间的变化规律（既插值函数），并将其作为近似解。有限体积法只寻求结点值；有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。,78,2.4 应用数值方法模拟的若干注意事项,1 简化模型简化细节（如倒角等等）；利用对称结构2 选择单元节点尽量少；曲（线）面边界用带（边）面节点的单元；轴对称结构采用轴对称单元3 划分网格细分网格；网格自适应（remeshing）4 建立初始条件和边界条件5定义材料参数,79,有限元常用软件,从20世纪70年代开始，基于有限元法在结构线性分析方面已经成熟，并被工程界广泛采用，一批由专业软件公司研制的大型通用商业软件（如NASTRAN，ASKA，SAP，ANSYS，MARC，ABAQUS，JIFEX等）公开发行和被应用。它包含众多的单元型式、材料模型及分析功能，并具有网格自动划分、结果分析和显示等先后处理功能。,80,81,几种有限元程序及其应用范围,82,几种有限元程序及其应用范围（续）,83,ANSYS 简介,1970 年成立的美国ANSYS公司是世界CAE行业最著名的公司之一,长期以来一直致力于设计分析软件的开发、研制,其先进的技术及高质量的产品赢得了业界的广泛认可。在我国,ANSYS用户也越来越多,三峡工程、二滩 电站、黄河下游特大型公路斜拉桥、国家大剧院、浦东国际机场等在结构设计时都采用了ANSYS作为分析工具。,84,ANSYS的分析方法,1.建立有限元模型,3.查看结果,2.施加载荷求解,主菜单,分析的三个主要步骤可在主菜单中得到明确体现.,Objective,2-2.ANSYS分析步骤在GUI中的体现.,85,ANSYS的分析方法(续),ANSYS GUI中的功能排列按照一种动宾结构，以动词开始（如Create),随后是一个名词(如Circle).,菜单的排列，按照由前到后、由简单到复杂的顺序，与典型分析的顺序相同.,86,ANSYS文件及工作文件名,一些特殊的文件 数据库文件jobname.db二进制Log 文件jobname.log文本结果文件jobname.rxx二进制图形文件jobname.grph二进制 ANSYS的数据库，是指在前处理、求解及后处理过程中，ANSYS保存在内存中的数据。数据库既存储输入的数据，也存储结果数据:输入数据-必须输入的信息(模型尺寸、材料属性、载荷等).结果数据-ANSYS计算的数值(位移、应力、应变、温度等).,87,启动ANSYS,Objective,1-1.启动ANSYS软件.,注意:以下过程默认为在 Windows NT环境下.,要启动ANSYS:,Windows NT 屏幕:Start Programs ANSYS X.x,1.选择 Interactive.,88,启动ANSYS(续),3.选择ANSYS的工作目录，ANSYS所有生成的文件都将写在此目录下。缺省为上次运行定义的目录。4.如果配置了3D显卡则选择3D.,5.设定初始工作文件名，缺省为上次运行定义的工作文件名，第一次运行缺省为 file.,7.选择 Run to 运行ANSYS.,Note:还可以设定其他的交互选项，但通常以上几项通常要设置好.如果在第1步中选择 Run Interactive Now，将读取上一次的设置，跳过此窗口，直接运行ANSYS,2.选择ANSYS产品.,6.设定ANSYS工作空间及数据库大小(参考ANSYS安装及配置手册).,89,ANSYS窗口,Objective,1-2.ANSYS GUI中六个窗口的总体功能,输入显示提示信息，输入ANSYS命令，所有输入的命令将在此窗口显示。,主菜单包含ANSYS的主要功能，分为前处理、求解、后处理等。,输出显示软件的文本输出。通常在其他窗口后面，需要查看时可提到前面。,应用菜单包含例如文件管理、选择、显示控制、参数设置等功能.,工具条将常用的命令制成工具条，方便调用.,图形显示由ANSYS创建或传递到ANSYS的图形.,90,一、主要功能简介,1.结构分析 1）静力分析-用于静态载荷.可以考虑结构的线性及非线性行为。线性结构静力分析 非线性结构静力分析 几何非线性：大变形、大应变、应力强化、旋转软化 材料非线性：塑性、粘弹性、粘塑性、超弹性、多线性弹性、蠕变、肿胀等 接触非线性：面面/点面/点点接触、柔体/柔体刚体接触、热接触 单元非线性：死/活单元、钢筋混凝土单元、非线性阻尼/弹簧元、预紧力单元等,91,2）模态分析-计算线性结构的自振频率及振形.谱分析 是模态分析的扩展，用于计算由于随机振动引起的结构应力和应变(也叫作 响应谱或 PSD).3）谐响应分析-确定线性结构对随时间按正弦曲线变化的载荷的响应.4）瞬态动力学分析-确定结构对随时间任意变化的载荷的响应.可以考虑与静力分析相同的结构非线性行为.5）谱分析 6）随机振动分析等 7）特征屈曲分析-用于计算线性屈曲载荷并确定屈曲模态形状.(结合瞬态动力学分析可以实现非线性屈曲分析.)8）专项分析:断裂分析,复合材料分析，疲劳分析,92,2.高度非线性瞬态动力分析（ANSYS/LS-DYNA）全自动接触分析，四十多种接触类型 任意拉格郎日欧拉（ALE）分析 多物质欧拉、单物质欧拉 适应网格、网格重划分、重启动 100多种非线性材料模式 多物理场耦合分析：结构、热、流体、声学 爆炸模拟，起爆效果及应力波的传播分析 侵彻穿甲仿真，鸟撞及叶片包容性分析，跌落分析 失效分析，裂纹扩展分析 刚体运动、刚体柔体运动分析 实时声场分析 BEM边界元方法，边界元、有限元耦合分析 光顺质点流体动力（SPH）算法,93,3.热分析 稳态、瞬态温度场分析 热传导、热对流、热辐射分析 相变分析 材料性质、边界条件随温度变化 4.电磁分析 静磁场分析计算直流电（DC)或永磁体产生的磁场 交变磁场分析 计算由于交流电(AC)产生的磁场 瞬态磁场分析计算随时间随机变化的电流或外界引起的磁场 电场分析用于计算电阻或电容系统的电场.典型的物理量有电流密度、电荷密度、电场及电阻热等。高频电磁场分析用于微波及RF无源组件，波导、雷达系统、同轴连接器等分析。,94,5.流体动力学分析 定常/非定常分析 层流/湍流分析 自由对流/强迫对流/混合对流分析 可压缩流/不可压缩流分析 亚音速/跨音速/超音速流动分析 任意拉格郎日欧拉分析（ALE）多组份流动分析（多达6组份）牛顿流与非牛顿流体分析 内流和外流分析 共轭传热及热辐射边界 分布阻尼和风扇模型 移动壁面及自由界面分析,95,6.声学分析 定常分析 模态分析 动力响应分析 7.压电分析 稳态、瞬态分析 模态分析 谐响应分析 8.多场耦合分析 热结构 磁热 磁结构 流体热 流体结构 热电 电磁热流体结构,96,9.优化设计及设计灵敏度分析 单一物理场优化 耦合场优化 10.二次开发功能 参数设计语言 用户可编程特性 用户自定义界面语言 外部命令 11.ANSYS土木工程专用包 ANSYS的土木工程专用包ANSYS/CivilFEM用来研究钢结构、钢筋混凝土及岩土结构的特性，如房屋建筑、桥梁、大坝、硐室与隧道、地下建筑物等的受力、变形、稳定性及地震响应等情况，从力学计算、组合分析及规范验算与设计提出了全面的解决方案，为建筑及岩土工程师提供了功能强大且方便易用的分析手段。,97,ANSYS主要模块简介,98,1.ANSYS/Mechanical：该模块提供了范围广泛的工程设计分析与优化功能，这些功能包括完整的结构、热、压电及声学分析。是一个功能强大的设计校验工具，可用来确定位移、应力、作用力、温度、压力分布以及其它重要的设计标准。2.ANSYS/Structural：通过利用其先进的非线性功能，该模块可进行高目标的结构分析，具体包括：几何非线性、材料非线性、单元非线性及屈曲分析。该模块可以使用户精确模拟大型复杂结构的性能。3.ANSYS/Linear plus：该模块是从ANSYS/Structural派生出来的，一个线性结构分析选项，可用于线性的静态、动态及屈曲分析，非线性分析仅包括间隙元和板/梁大变形分析。,99,4.ANSYS/Thermal：该模块同样是从ANSYS/Mechanical中派生出来的，是一个可单独运行的热分析程序，可用于稳态及瞬态热分析。5.ANSYS/Flotran：该程序是个灵活的CFD软件，可求解各种流体流动问题，具体包括：层流、紊流、可压缩流及不可压缩流等。通过与ANSYS/Mechanical耦合，ANSYS/FLOTRAN是唯一一个具有设计优化能力的CFD软件，并且能提供复杂的多物理场功能。6.ANSYS/Emag：该程序是一个独立的电磁分析软件包，可模拟电磁场、静电学、电路及电流传导分析。当该程序与其它ANSYS模块联合使用时，则具有了多物理场分析功能，能够研究流场、电磁场及结构力学间的相互影响。7.ANSYS/Preppost：该模块为用户在前处理阶段提供了强大的功能，使用户能够便捷地建立有限元模型。其后处理器能够使用户检查所有ANSYS分析的计算结果。,100,8.ANSYS/ED：该模块是一个功能完整的设计模拟程序，它拥有ANSYS隐式产品的全部功能，只是解题规模受到了限制（目前节点数1000）。该软件可独立运行，是理想的培训教学软件。9.ANSYS/LS-DYNA：该程序是一个显示求解软件，可解决高度非线性结构动力问题。该程序可模拟板料成形、碰撞分析、涉及大变形的冲击、非线性材料性能以及多物体接触分析，它可以加入第一类软件包中运行，也可以单独运行。10.ANSYS/LS-DYNA PrepPost：该程序具有所有的ANSYS/LS-DYNA的前后处理功能，具体包括：实体建模、网格剖分、加载、边界条件、等值线显示、计算结果评价以及动画，但没有求解功能。11.ANSYS/University：该模块是一个功能完整的设计模拟程序，它拥有ANSYS隐式产品的全部功能，只是解题规模受到了限制（目前节点数16000和32000两种）。该软件可独立运行，适用与高校进行教学或科研。,101,12.ANSYS/DesignSpace：该模块是ANSYS的低端产品，适用与设计工程师在产品概念设计初期对产品进行基本分析，以检验设计的合理性。其分析功能包括：线性静力分析、模态分析、基本热分析、基本热力耦合分析、拓扑优化。其他功能有：CAD模型读取器、自动生成分析报告、自动生成ANSYS数据库文件、自动生成ANSYS分析模板。产品详细分类：DesignSpace for MDT DesignSpace for SolidWorks Standalone DesignSpace：（支持的CAD模型有：Pro/E、UG、SAT、Parasoild）13.ANSYS/Connection：ANSYS与CAD软件的接口产品。可以将CAD模型数据或国际标准格式CAD模型数据直接读入并进行任意ANSYS支持的分析。目前支持的CAD软件有：Pro/E、UG、CADDS。支持的国际标准格式有：IGES、SAT、Parasolid.,