1、1,第六章 计算机辅助工程分析学习资料,第一节 有限元分析法,第二节 优化设计,第三节 计算机仿真简介,2,机械产品设计过程的一个重要环节是分析、计算，其中包括对产品几何模型进行分析、计算，通过应力变形进行结构分析，对设计方案进行分析、评价等。 传统的分析方法：采用手工计算，过程繁琐，效率低，往往需要对原型进行较大的简化并引入更多的假设，分析方法一般比较粗略，分析结果不够可靠，只能用来对产品设计方案作定性的比较分析，而不能作出定量的评价。,3,工程分析：主要通过计算机，利用数值分析方法进行辅助工程分析，是CAD中应用最早、卓有成效的领域之一。分析的关键是在三维实体建模的基础上，从产品的方案设计
2、阶段开始，按照实际使用的条件进行仿真和结构分析；按照性能要求进行设计和综合评价，以便从多个设计方案中选择最佳方案。 计算机辅助工程分析:通常包括有限元分析、优化设计、仿真（模拟分析）可靠性分析、试验模态分析等。,4,1有限元分析法：是力学与近代计算机技术相结合的产物，是一种解决工程问题的数值计算方法，1960年美国Clogh教授首次提出“有限元法（The Finite Element Method）”的概念。,第一节 有限元分析法,一有限元概述,5,2分类有限元法包括有限元建模和有限元分析两部分，目前它们已成为建立分析模型、共享数据的有效途径，是解决各种工程实际问题的便利工具和有效手段。 3应
3、用有限元法可以处理任何复杂形状、不同物理特性、多变的边界条件和任何承载情况的工程问题，广泛应用于场强（力场、电场、磁场、温度场、流体场等）分析、热传导、非线形材料的弹塑性蠕变分析等研究领域中。,6,1有限元法的基本思想先把一个原来是连续的物体剖分（离散）成有限个单元，而且他们相互连接在有限个节点上（如图所示），承受等效的结点载荷（由静力等效原则转化为节点上的等效载荷），并根据平衡条件（应用虚位移原理建立平衡条件）进行分析，然后根据变形协调条件把这些单元重新组合起来，成为一个组合体，再综合求解。由于单元的个数是有限的，结点数目也是有限的，所以称为有限元法。,二有限元法的基本思想,7,2单元类型在
4、采用有限元法对结构进行分析计算时，分析对象不同，采用单元类型（形状）也不同。常见的单元类型有：杆单元、梁单元、板单元（三角形、矩形等）、多面体单元（四面体、六面体）等。,8,有限元问题最终归结为：在满足边界条件的情况下，求解基本方程，在实际求解时，先求某些未知量，再由它们求解其他未知量。,三 有限元法的基本解法,9,1有限元求解的方法模式经典的最小势能原理的位移有限元模式（位移法）基于余能原理的应力平衡模式（应力法）基于广义势能原理的位移杂交模式（位移杂交法）基于广义余能原理的应力杂交模式（应力杂交法）基于H-W混合变分原理的混合有限元模式等有限条带法、无限元、半解析有限 元以及边界元等，从而
5、大大提高了有限元法解决实际问题的能力。,10,2有限元法求解的具体步骤1)单元剖分把连续弹性体分割成许多个有限大小的单元，并为单元和结点编号。2)单元特征分析3)总体结构合成,11,12,1前置处理前置处理是在用有限元法进行结构分析之前，按所使用的单元类型对结构进行剖分；根据要求对结点进行顺序编号；输入单元特性及结点坐标；生成网格图象并显示；为了决定它是否适用或者是否应当修改，显示的图象应带有结点和单元标号以及边界条件等信息；为了便于观察，图象应能分块显示、放大或缩小。对于三维结构的网格图象需要具备能使图象作三维旋转的功能等。以上内容一般称之为前置处理，为实现这些要求而编制的程序称为前置处理程
6、序。,四、有限元分析的前置处理和后置处理,13,2前置处理程序基本功能1）生成结点坐标并按顺序编号；2）生成网格单元，下图为一复杂曲面的网格划分图；3）修改和控制网格单元；4）引进边界条件以约束一系列结点的总体位移和转角；5）单元物理几何属性编辑，如材料特性、弹性模量、厚度、惯性矩以及泊松比等；,14,15,3后置处理后置处理，是将有限元计算分析结果进行加工处理并形象化为变形图、等值线图、应力应变彩色图、应力应变曲线以及振型图等，以便对变形、应力等进行直观分析和研究。下图所示为某受力体表面应力等值线图。为了实现上述目的而编制的程序，称为后置处理程序。,16,17,4后置处理内容： （1）数据输
7、出：将结点位移、单元应力按设计者的意图整理输出，还可从大量数据中筛选出关键的有用数据，按用户要求的格式输出规格化的数据文件。,18,（2）图形显示：图形显示和绘图可形象直观地表示有限元模型和计算结果，可帮助设计者迅速了解研究对象的特征，从而对修改模型作出判断。图形显示包括有限元网格图、结构变形图、等值线图以及振型图等。等值线有应力等值线图、位移等值线图、等高线图和温度等值线图等，其中在工程结构分析中，以应力等值线图应用最多。等值线图可在彩色屏幕上用不同的颜色加以形象化。下图所示为一曲面的彩色等高线图。,19,20,1将几何建模系统和有限元分析系统有机结合，在建模系统中将有限元的前后置处理作为线
8、框建模、表面建模、实体建模的应用层。即将几何建模的几何参数和拓扑关系等数据进行加工，自动剖分成有限元的网格，在输入有限元的其他数据，生成不同有限元分析程序所需的数据网格文件。,五有限元前、后置处理的应用分类：,21,2单独为某一有限元分析程序配置前后处理功能程序，并把二者集成为一套完整的有限元分析系统，他同时具有批处理和图形编辑功能。,22,五、有限元软件,最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大
9、学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。,23,第二节 优化设计,优化设计原则：寻求最优设计优化设计手段：计算机和应用软件优化设计理论：数学规划法。,24,一、优化设计的基本概念,优化设计要解决两个关键问题： （1）建立优化设计数学模型，它包括三个要素，即优化设计的目标函数、设计变量和约束条件； （2）选择适合的优化方法,25,1设计变量1）概念在设计中，通常用对设计性能指
10、标有影响的一组基本参数来表示某个设计方案，这组参数根据其特点又分为设计常量：可以根据客观规律或具体条件预先确定的参数，如材料的力学性能，机器的工况系数等。设计变量：在设计过程中不断变化，需要在设计过程中进行选择的基本参数，称为设计变量，如几何尺寸、速度、加速度、温度等。,26,【例5-1】设计一密闭矩形容器，其容积为3m3，容器的宽度不小于1.5m，以便于装卸车搬运，为使成本最低，要求用料最省。这项设计的设计变量是矩形容器的长l、宽w、高h。,27,2）对设计变量的要求相互独立的参数，其取值都是实数。有界连续变量，称为连续量，如温度。跳跃式的量，称为离散量，如齿轮的齿数、模数，丝杆的螺距等。对
11、离散变量，在优化设计时，常常先看作连续量，在求得连续量的优化结果后在进行圆整或标准化，以求得一个实用的最优方案。,28,一项设计，若有n个设计变量1，2，n 可以按一定次序排列，可以用n维向量来表示：=1，2，nTRn它表示了设计空间的概念，即以n个设计变量为坐标轴组成的实空间，Rn代表n维空间,29,2目标函数1）概念：目标函数是指根据特定目标建立起来的、以设计变量为自变量的一个可计算的数学函数。它是设计方案评价的标准。,30,2）目标函数的描述：优化设计的过程实际上是求目标函数极小值或极大值的过程，而求目标函数极大值的问题可转化为求目标函数极小值的问题。优化设计数学模型中通常规定求目标函数
12、的极小值。故目标函数统一描述为：min F（）= F（1，2，n ） 例5-1中密闭容器优化设计的目标函数可表示为：min F（）=F（l, w, h）=2（lh+wh+lw）,31,3）分类单目标优化：如果优化问题只有一个目标函数，则称为单目标优化.多目标优化：如果优化问题有几个目标函数，则称为多目标优化。,32,3约束条件1）概念：为产生一个可接受的设计，设计变量本身或相互间应该遵循的限制条件，称为约束条件。,33,2）表示方法：约束条件一般可表示为设计变量的不等式约束函数形式和等式约束函数形式，即gi（）= gi（1，2，n）0 gi（）= gi（1，2，n）0 （i =1，2，m）hj
13、（）= hj（1，2，n）= 0 （ j =1，2，p） 式中m，p分别表示施加于该项设计的不等式约束数目和等式约束数目。,34,3）分类按约束条件，又可分为性能约束和边界约束。（1）性能约束 又称性态约束，是针对设计对象的某种性能或指标而给出的一种约束条件。如零件的计算应力不大于许用应力，轴的扭转变形应小于许用扭转角等。一般这类约束条件总可以根据设计规范中的设计公式或通过物理学和力学的基本分析导出的约束函数来表示。,35,（2）边界约束 又称区域约束，表示设计变量的物理限制和取值范围。如例5-1 中的边界约束条件为： g1（）= g1（l，w，h）= w 1.5 0 g2（）= g2（l，w
14、，h）= l 0 g3（）= g3（l，w，h）= h 0约束条件必须是以设计变量为自变量一个有定义的函数，并且各个约束条件之间不能彼此矛盾。,36,1数学模型的描述数学模型的规范化描述形式为：min F（） =1，2，nT Rngi（）= gi（1，2，n）0 i =1，2，m hj（）= hj（1，2，n）= 0 j =1，2，p,二、优化设计的数学模型,37,2线形规划当目标函数F（）、约束条件gi（）和hj（）是设计变量的线形函数时，称该优化问题为线形规划问题.3非线形规划当目标函数F（）、约束条件gi（）和hj（）中有一个或几个是设计变量的非线形函数时，称该优化问题为非线形规划问题；
15、机械设计中，绝大多数优化设计问题的数学模型都属于非线形规划问题。,38,例5-1中的优化设计问题可建立如下的数学模型：minF（）=F（l, w, h）=2（lh+wh+lw）g1（）= g1（l，w，h）= w 1.5 0 g2（）= g2（l，w，h）= l 0g3（）= g3（l，w，h）= h 0g4（）= g3（l，w，h）= l w h = 3可见，这是一个三维非线形规划问题。为了简化问题，可根据等式约束条件消去一个设计变量：h = 3 /（ l w） 则该问题从原来的三维问题转化为二维问题。,39,4建立数学模型的一般过程1）分析设计问题，初步建立数学模型 即使是同一设计对象，如
16、果设计目标和设计条件不同，数学模型也会不同。因此，要首先弄清问题的本质，明确要达到的目标和可能的条件，选用或建立适当的数学、物理、力学模型来描述问题,40,2）确定设计变量 通常应参照以往的设计经验和实际要求，将那些对目标函数影响不大的参数取为常量，以减少设计变量的个数，有利于设计问题数学模型的简化。 3）根据工程实际，提出约束条件 约束条件是对设计变量的限制，这种限制必须要根据工程实际情况来制订，以便使设计方案切实可行。约束条件的数目越多，则可行的设计方案就越少，优化设计的难度也越大。,41,4）对照实际设计，修正数学模型 初步建立模型之后，应与设计问题加以对照，并对函数值域、数学精确度和设
17、计性质等方面进行分析，若不能正确、精确地描述设计问题，则需用逐步逼近的方法对模型加以修正 5）选择正确的计算方法，确定计算误差 如果数学模型的数学表达式比较复杂，无法求出精确解，则需采用近似的数值计算方法，并估计计算误差。,42,6）进行结果分析，审查模型灵敏度 数学模型求解后，还应进行灵敏度分析，即在优化结果的最优点处，稍稍改变某些条件，检查目标函数和约束条件的变化程度。若变化大，则说明灵敏度高，就需要重新修正模型。,43,一维搜索法 一维函数黄金分割法（0.618法） 二次插值法（近似抛物线法） 坐标轮换法 鲍威尔法（Powell法） 梯度法 牛顿法 变尺度法（DRP） 复合形法 惩罚函数法（罚函数法）,三常用优化设计方法,44,目前，优化设计软件已成为一种比较成熟的软件产品。例如OPB，MOPB，PLODM和DDDU等。这些程序库中各种算法齐全，使用方便，使用者只要会选择适当的算法，并按既定格式编写具体信息的子程序即可。在CAD/CAM中，应尽可能选用现成的优化方法软件，以节省人力、机时，尽快得到优化设计结果，满足CAD/CAM的需要。,四、优化设计软件,45,五优化设计的步骤,从设计方法来看，机械优化设计和传统的机械设计方法有本质的差别。图5-6所示为优化设计的一般步骤,
