UG 有限元分析 说一说UG有限元分析的具体步骤

建模以后按以下步骤进行：

1、进入结构分析模块（UG4.0是设计仿真,需要自己安装FEA软件）

2、选择分析解算器（3.0自带structure P.E.可以用）

3、加载荷

4、加约束（固定约束最好有一个）

5、赋材料

6、网格化

7、模型检查

8、解算

解算以后进入导航器看result,可以看到位移(displacement),应力（strss)、应变（strain)等的云图,其它后处理不做介绍

首先讨论了如何在UG软件中完成客车车身的数值模型及如何将此数值模型进行简化转化成客车车身有限元模型，接着在ANSYS软件中对设计的客车车身骨架结构进行了静态弯曲工况、扭转工况和弯扭工况三种工况下，车身结构的强度和刚度的分析，并对该车进行了动态分析。

基于UG软件的客车车身曲面设计，客车车身曲面不同于轿车车身曲面，其曲面最复杂的地方集中于车头和车尾，侧围和顶盖的曲面相对而言较为简单。所以对于客车车身外表面最方便易性的构造方法是直接由车身的二维轮廓线出发，在计算机上绘制出车身的主要轮廓线，再由这些轮廓线出发构造车身外表面模型。由此，我们定出了9根车身外表面轮廓线，通过这些轮廓线可确定车身外表面的基本形状。如侧围曲面可由。1曲线沿c2曲线平行扫掠构成顶盖曲面由c6,c4和c8曲线沿0曲线扫掠而成后围曲面由0和c9曲线沿c8曲线扫掠而成前围曲面较为复杂，除需要车身外表面主要轮廓线c5和c6曲线外，还需根据车身的造型特点，再另外构造3根曲线，才能生成前围曲面。

侧围主视向轮廓线(客车左右侧对称，可任选一根)侧围俯视向轮廓线(一般中间是直线，两端向前后围缩一偏移顶盖侧视向轮廓线顶盖主视向轮廓线(一般顶盖为大圆弧，两端为与侧围主视向轮廓的上部相切的倒圆弧线)前围与侧围相交处轮廓线

(客车左右侧对称c6:前围与顶盖相交处轮廓线c7:后围与侧围相交处轮廓线(客车左右侧对称)c8:后围与顶盖相交处轮廓线

c9:后围侧视向轮廓线。

为确保轮廓线的光顺性，使用UG软件的曲线分析功能，对这9根车身外表面轮廓线的曲率进行分析、编辑和调整。

3种分析方法:<br />

1.基于UG软件的车身骨架设计

由于客车车身骨架截面在各个不同的空间位置上其形状和大小都保持不变，故用UG构造客车车身骨架时，可采用曲面扫描法，求出骨架杆件截面的空间运动轨迹(即车身骨架杆件外表面中心线)，将该截面沿其空间运动轨迹扫掠即可得车身骨架的实体模型。又因为客车车身骨架是一个空间多层次的杆件结构，分为底架，前围、后围、左侧围、右侧围和顶盖六大部分，在具体设计时，先根据六大片的设计参数进行布局设计，一般是先进行底架布局设计，确定底盘各总成的具体布置位置后，再根据底架设计中的一些关键参数进行前、后围、左右侧围及顶盖的设计然后利用在UG上已建立好的车身表面数字模型和骨架六大片布局设计参数求取车身骨<br />

架与车身表面数值模型的截交线即车身骨架杆件外表面中心线，构造出车身六大片的线框模型。根据客车车身结构需要，选取合适的骨架构件截面，如矩形、槽形、L形(角钢)和工字型等，由此截面沿车身六大片的线框模型扫掠构造出车身六大片骨架实体模型。最后利用UG的装配模块，进行整车装配，生成车身骨架图。

2模型的简化

因为建立车身有限元模型时，既要如实的反映客车车身实际结构的重要力学特性，又要尽量采用较少的单元和简单的单元形态，以保证较高的计算精度及缩小解题规模。在有限元模型中，我们一般人为的用一根通过截面形心的直线来代替具有一定横截面尺寸的实际构件。所以在利用ANSYS软件的数据接口程序导入在UG中完成的客车骨架结构图时，只需导入车身骨架线框图并对其进行以下简化:1、略去蒙皮和某些非承载构件2、将车身中的各微曲梁进行直化处理，侧围和顶盖中一些曲率较小的构件近视的看作由直梁单元分段组成3、对于两个靠得很近而又不重合的交叉连接点则可考虑简化为一个节点来处理。4、对于邻接构件在空间交接的轴线不重合，出现了两个离得很近的节点，在力学特性上它们的变形很接近，把它们简化成一对主从节点，这样就避免了可能出现的总刚度阵的病态，同时也可提高结构分析的效率。5、对于空间叠交的两焊接梁，若其中心线的距离a较大，平移其中一梁中心线将引起不可忽略的误差，则可于模型中加一个长度为a的梁(该梁截面、材料特性参数为两梁中较大者)来连接两梁。例如，底横梁与车架纵梁之间就存在着很大的“偏心”，横梁置于车架纵梁之上，两轴线相差距离为0.5(h+H)。为了使模型根接近实际，将底横梁于车架纵两连接处分量各节点考虑，并设其间有一刚臂连接6、对于两同向焊接梁，因其焊接处强度近似于材料内部强度，故可将其视为一根梁来简化

7、对线梁单元采取刚度补偿的方法来降低误差。以线单元表示梁，要满足梁相交的空间拓扑关系，须将其中的某些梁单元线延长至相交，这样处理将大大降低梁单元的刚度，使得位移解偏大而应力解偏小，同时增加了额外的重量。采用刚度补偿的方法来降低误差，经补偿前后结果比较后，己验证了该方法简单有效。以梁单元xoy平面内弯曲((1轴为x轴，2轴为z轴)为例，说明该补偿方法。采用二节点Hermite单元的有限元求解方程Ka=p的单元刚度矩阵K“和位移矢量1其中，l为梁单元沿1轴的长度，。，为单元节点1处的挠度，乓为单元节点l处的转角，由于模型中的梁单元比实际的延长了△l，故可通过改变E或者Iz来抵消该变化，使K“基本不变。8、确定单元长度l。用有限元法分析梁弯曲问题时，于二节点Hermite单元中，试探函数(形函数)采用3阶完全多项式，位移解的误差是o(l小若梁单元长度过长，则会引起较大的位移误差。在分析车身梁单元模型时，经FEA验证当梁单元长度15400mm时，其解已收敛到足够的精度。梁单元长度l也不应划分得过小，若梁单元长度Z过小(接近于截面尺寸)，主从自由度的原理将不再适用，模型单元简化为梁单元也就不合理。各相邻梁单元长度1相差也不应过大，理论和实践已证明，l相差过大将引起较大的刚度壁，这易导致刚阵病态而得不到方程组的解。根据以上模型的简化原则，样车车身骨架被划分为3044个长度不等，截面形状各异的单元和5929个节点。

3载荷处理

在车身计算模型中，载荷可按如下方式处理:1、对于车身骨架的自重，在ANSYS软件前处理程序中输入骨架材料的密度和重力加速度，程序便根据所输入的单元截面形状、实常数自动将单元载荷因子的信息计入总载荷，进行计算

2、安放在车身或车架上的汽车总成、设备重力，如发动机总成、备胎、蓄电瓶、油箱等，可作为集中载荷，按安放点的实际位置及各位置所分担的重力，作用于相应的节点上。

3、载重力，如乘员及座椅的重力，可作为集中载荷，按支点跨距分配于相应梁的结点上。车上有站立乘员者，可按每平方实际站立人数，作为均布载荷作用于地板上并传到底架梁单元上。由于在有限元法中认为内力或外力均由结点来传递，在整体刚度方程中的载荷项均为结点载荷。因此，当梁单元受有均布载荷或其他非节点载荷时，必须将其向结点移置，即将非结点载荷换算成作用在结点上的效果相当的集中载荷(称等效结点载荷)。非结点载荷移置方法如下

有非结点载荷作用的单元的两端位移完全约束住，再根据材料力学中求支反力的方法，求得梁单元两端的反力，称固端力，记作仇}02、将固端力反号，并进行坐标变换，即得整体坐标系中的等效结点载荷，可将它直接送入结构整体刚度方程的载荷向量中去进行计算。在ANSYS软件中，如果先在车身有限元模型上加载再进行网格划分能直接将非结点载荷转换成等效结点载荷。<br />

4边界约束条件

钢板d簧除了作d性元件外，还起导向作用，因此其在各个方向上均有刚度，且其在其他方向上的刚度要比垂直方向上的刚度大得多，故用刚性梁一柔性梁结构模拟钢板d簧。在约束处理中忽略轮胎的变形。悬架d簧刚度K用水平柔梁的垂直弯曲刚度来等效对于刚性梁，为使其受力时垂直位移远小于水平柔梁的垂直位移，取其轴向刚度为6.0x106N/mm。刚性梁截面取为正方形，面积由式A=KxLIE计算。

5.强度分析工况<br />

客车的使用工况很复杂，有弯曲工况、扭转工况、转弯工况和加速工况等。理论分析、室内试验和使用实践都表明，直接关系到车身结构强度的主要是弯曲和扭转两种工况。<br />

I、弯曲工况<br />

客车在平坦路面上以较高车速行驶时，路面的反作用力使车身承受对称的垂直载荷。它使车身产生弯曲变形，其大小取决于作用在车身各处的静载荷及垂直加速度。在ANSYS中通过约束四车轮六个方向的自由度来模拟计算客车在平坦路面上，以较高车速满载行驶产生对称垂直动载荷时，车身的刚度和强度。2、扭转工况<br />

扭转工况是车身变形最严重的工况，一般都是当汽车以低速通过崎岖不平路面时发生的。此种扭转工况下的动载，在时间上变化得很缓慢，当然此时惯性载荷也很小，所以，车身的扭转特性可以近似的看作是静态的，许多试验结果也都证实了这一点，即静扭试验下的骨架强度可以反映出实际强度。也就是说，静扭时骨架上的大应力点，就可用来判定动载时的大应力点。文中将讨论两种扭转工况，右前轮悬空工况和左后轮悬空工况。通过约束左后轮X,Y,Z方向的平动自由度和Z方向的转动自由度，左前轮和右后轮Z方向的平动自由度，来模拟车身右前轮悬空，左后轮陷入坑中的扭转工况。通过约束右前轮X,Y,Z方向的平动自由度和Z方向的转动自由度，左前轮和右后轮Z方向的平动自由度，来模拟车身左前轮悬空、右后轮陷入坑中的扭转工况。<br />

4.2.2刚度分析工况<br />

车身结构的刚度是指车身结构反映出的载荷与变形之间关系的特性。刚度不足，会引起车身的门框、窗框等开口处的变形大，以至车门卡死、玻璃砸碎、密封不严导致漏雨、渗水及内饰脱落等问题，还会造成车身振动频率低、发生结构共振，破坏车身表面的保护层和车身的密封性，从而削弱抗腐蚀能力。车身刚度包括扭转刚度和弯曲刚度两部分，理论分析和许多试验结果都表明，客车车身的弯曲变性很小，故只需考虑其弯扭工况下的扭转刚度。我们用整车总长之间车身对角线相对扭角、左右上大梁的相对扭角状况、底架两纵梁的相对扭角状况来表达车身的扭转变形。<br />

4.2.3动态特性研究<br />

用模态综合法来研究整车振动特性和动载荷时，车身结构的模态频率是最重要的参数之一。用它能够预测车身与其它部件如悬挂系统、路面、发动机及传动系等系统之间的动态干扰的可能性，通过合理的设计可以避开共振频率，一般希望车身结构整体一阶模态频率越高越好。<br />

<br />

4.3.1强度计算结果及分析<br />

1、弯曲工况<br />

弯曲工况下，车身的弯曲应力如图4.3所示。弯曲应力集中的区域有:底架主纵梁与前后钢板d簧支撑梁位置处(50-90Mpa)车顶中部与侧窗上沿的过渡连接区(30-40Mpa)中门立柱上半部的附近区域(10-30MPa)前门立柱上半部的附近区域(10-40MPa)。其中应力最大的地方是底架主纵梁与后钢板d簧支承梁位置处，应力值为90MPao<br />

<br />

2、右前轮悬空工况<br />

右前轮悬空工况下，车身X方向的应力分布如图4.4所示。应力集中的区域有:底架主纵梁与前后钢板d簧支撑梁位置处(60-123Mpa)车顶中部与侧窗上沿的过渡连接区(40-60Mpa):中门立柱上半部的附近区域(60-70MPa)。其中应力最大的地方是底架主纵梁与后钢板d簧支承梁位置处，应力值为123Mpao<br />

3、左后轮悬空工况<br />

左后轮悬空工况下，车身X方向的应力分布如图4.5所示。应力集中的区域有:底架主纵梁与前后钢板d簧支撑梁位置处(80一125Mpa)车顶中部与侧窗上沿的过渡连接区(60一90Mpa)中门立柱上半部的附近区域(90-177MPa)。其中应力最大的地方是中门上门梁位置处，应力值为177MPa<br />

<br />

通过上述三种工况的计算，我们知道弯曲工况下车身骨架的应力水平较小，应力值大于50Mpa的单元数目为30个，仅占单元总数的0.9%右前轮悬空工况下，车身骨架的应力水平要比弯曲工况下的应力水平高很多。由于发动机后置，左后轮悬空工况(弯扭联合工况)是客车行驶过程中最恶劣的工况。考虑到客车行驶过程中的动载荷、疲劳及材料缺陷引起的应力集中等问题，取安全系数为1.5，则骨架材料Q215A3钢的许用屈服应力[cr]-153MPa，底架材料09SiV低合金结构钢的许用屈服应力<br />

叶卜220MPa。可以看出，在弯扭工况下，中门上门梁位置处的应力超过了许用应力，需要对门梁的截面尺寸进行优化。另外，从整个结构来看，应力分布是不均匀的，且大小相差几个数量级。这无疑将造成材料的浪费，增大整个车身的重量。因此，从应力角度分析，可以通过优化方钢厚度来合理经济的使用材料。多梁相交处的应力值特别大，去掉一些可取掉的单元后，交点处的应力值将大大降低。<br />

4.3.2刚度计算结果及分析<br />

1、右前轮悬空工况<br />

车身右前角区域从车顶至车架各部分均有较大的位移，而且越靠近角<br />

部位移越大，垂直方向向下的最大位移为11.868mm。车身变形如图4.6<br />

<br />

4.3.3模态计算结果及分析<br />

模态分析主要是计算车身固有频率和振型。整体车身空间框架模型的6阶固有频率如表4.9所示，前六阶振型车身的变形如图4.10-4.15所<br />

<br />

图4.15车身骨架第六阶振型图<br />

车身骨架的动态优化设计要求车骨架的模态频率错开载荷激振频率。同时为防止第一阶弯曲模态和第一阶扭转模态的祸合效应，要求这两种固有频率错开3Hz以上。虽然由于客车模型略去了蒙皮的影响，略去了非承载构件，所计算的车身固有频率比实际的要低，但是该车前六阶固有频率集中在5-13Hz，而路面激励频率又往往低于20Hz，且第一阶弯曲模态和第一阶扭转模态的固有频率仅错开了2Hz左右，因而在客车行驶过程中产生局部振动的构件受此激励将在客车内部形成噪声源，影响到乘客的乘座舒适性。<br />

4.4结论<br />

从原模型计算结果可以看出，该车车身骨架的高应力区共有3个部位:中门立柱附近区域车顶中部与侧窗上沿的过渡连接区和底架主纵梁与前后钢板d簧支撑梁位置处。2、由计算结果可知，该车在弯曲工况下，骨架的变形和应力均较小，表明该车在静载下满足强度和刚度要求在左后轮悬空工况下，除了车身中门门上梁中间部位应力超过了许用应力，车身骨架的其他单元应力都未超过许用应力。而左后轮悬空工况是车身变形最严重的工况，实际上由于该车是城市公交车，不可能出现如此严重的扭转工况，因此该车车身结构是能够满足强度使用要求。<br />

3、由计算结果可知，总体上车身骨架的变形量相对较小，对于车身刚度而言，从整体结构考虑，门窗对角线变形大小尤为重要。从整理的弯扭工况下车身骨架各节点变形数据中可看出，弯扭工况下各门窗对角线位移均较小，因此该车车身结构是能够满足刚度使用要求的。<br />

4、由车身模态分析可知，车身骨架前六阶的固有频率都低于20Hz，而路面激励频率又往往低于20Hz，这会造成车身骨架发生共振，造成车内噪声过大，因此进行车身结构的动态优化设计，提高车身的固有频率很有必要。<br />

<br />

5车身结构的优化设计<br />

5.1优化设计的基本概念一般的工程问题都有许多可行的设计方案，如何根据设计任务和要求从众多的可行性方案中，寻求一个最好的方案，是设计工作者的首要任务。实践证明，结构的优化设计是保证产品具有优良的性能，减轻结构自重或体积，降低工程造价的一种有效方法。优化方法的出现可追溯到Newton,Lagrange和Cauchy时代，由Newton,Leibnitz和Weirstrass等奠定了变分学的基础Lagrange创立了包含特定乘子的约束问题优化方法，并将其命名为Lagrange乘子法Cauchy最早应用最速下降法来求解无约束极小化问题。尽管如此在20世纪中以前，优化法的进展甚小。直到后来，高速计算机的出现，才使优化程序成为可能，促使了各种新方法的进一步发展。五十年代以前，用于解决最优化问题的数学方法，仅限于古典微分法和变分法。无约束优化数值方法领域中的主要进展只是在60年代才在英国形成，数学规划方法被首次用于结构最优化，并成为优化设计中求优方法的理论基础，线性规划和非线性规划是其主要内容。1947年，Dantzig提出求解线性规划问题的单纯形法1957年，Bellman对动态规划问题提出了最优化理论。60年代初，Zoutendijk和Rosen对非线性规划右很大贡献。Canon,Fiacco和Mclomick的研究使很多非线性规划问题能用无约束优化方法予以解决。几何规划是60年代由Duffin,Zener和Peterson发展起来的。概括来讲，优化设计工作包括以下两部分内容:1、将设计问题的物理模型转变为数学模型，建立数学模型时要选取设计变量，列出目标函数，给出约束条件。2、采用适当的优化方法，求解数学模型，可归结为在给定的条件下求目<br />

标函数的极值和最优化值的问题。机械最优化设计，就是在给定的载荷或环境条件下，在对机械产品的性能、几何尺寸关系或其他因素的限制范围内，选取设计变量，建立目标函数并使其获得最优化值的一种设计方法。实际的工程优化设计按其原理不同区分为数学规划法和准则法两个分支，按其优化层次不同可分为总体方案优化和设计参数优化。<br />

5.2ANS丫S软件中的设计优化<br />

ANSYS程序提供了分析一评估一修正的循环过程对设计方案进行优化，对初始设计进行分析，根据设计要求对分析结果进行评估，然后对设计进行修正。重复执行这一循环过程直到所有设计都满足要求，得到最优设计方案。<br />

5.2.1优化方法<br />

ANSYS提供了零阶方法和一阶方法两种优化方法。绝大多数的优化问题都可以使用这两种方法。零阶方法(直接法)是一个很完善的处理方法，其中有两个重要的概<br />

念:目标函数和状态变量的逼近方法，由约束的优化问题转换为无约束的优化问题。该方法使用所有因变量(状态变量和目标函数)的逼近，而不用他们的导数，用因变量的近似值工作，而不用实际函数目标函数近似为最小值，而不是用实际的目标函数状态变量近似为使用设计约束，而不用实际状态变量，可以很有效的处理大多数的工程问题。所有变量至少要适应所有的全部现有设计集，以形成近似式:<br />

一阶方法(间接法)基于目标函数对设计变量的敏感程度，使用因变量的一阶导数来决定搜索方向并获得优化结果，因为没有近似，所以精度很高，尤其是在因变量变化大，设计空间也相对较大时，更加适合于精确的优化分析。每次迭代涉及多次分析(对分析文件的多次循环)，以确定适当的搜索方向，因此分析时间较长。当零阶方法不够精确，而精度又非常重要时，要用一阶方法进行优化。<br />

5.2.2优化工具<br />

ANSYS程序还提供了一系列的优化工具以提高优化过程的效率。优化工具是搜索和处理设计空间的技术。下面是常用的优化工具:单步运行:实现一次循环并求出一个FEA解。可以通过一系列的单次循环，每次求解前设定不同的设计变量来研究目标函数与设计变量的变化关系。随机搜索法:进行多次循环，每次循环设计变量随机变化。可以指定最大循环次数和期望和理解的数目。主要用来研究整个设计空间，并为以后的优化分析提供合理解。往往作为零阶方法的先期处理。等步长搜索法:以一个参考设计序列为起点，生成几个设计序列。按照单一步长在每次计算后将设计变量在变化范围内加以改变，用于设计空间内完成扫描分析。对于目标函数和状态变量的整体变化评估可以用本工具实现。<br />

乘子计算法:是一个统计工具，用二阶技术生成设计空间上极值点上的设计序列数值。主要用来计算目标函数和状态变量的关系和相互影响。最优梯度法:对用户指定的参考设计序列，计算目标函数和状态变量对设计变量的梯度，可以确定局部的设计敏感性。<br />

5.2.3优化变量<br />

设计变量、状态变量和目标函数总称为优化变量。设计变量为自变量，优化结果的取得就是通过改变设计变量的数值来实现的。状态变量是约束设计的数值，是“因变量”，是设计变量的函数，状态变量可能会有上下限，也可能只有单方面的限制，即只有上限或下限。目标函数是设计最小化或最大化的数值，是设计变量的函数。目标函数值由最佳合理设计到当前设计的变化应小于目标函数允差。一个合理的设计是指满足所有给定的约束条件(设计变量的约束和状态变量的约束)的设计。如果其中任一约束条件不满足，设计就被认为是不合理的。而最优设计是既满足所有的约束条件又能得到最小目标函数值得设计。(如果所有的设计序列都是不合理的，那么最优设计是最接近合理的设计，而不考虑目标函数的数值)<br />

5.3车身骨架的优化设计<br />

5.3.1参数化优化模型<br />

进行车身骨架的优化设计首先必须要建立车身骨架的参数化模型，我们采用了车身骨架的早期静态有限元模型，作为其参数化模型的原型。由于该模型的建立没有参数化，所以必须重新划分单元，简化模型，使骨架单元数控制在4000个以下，模型的简化过程中保持计算偏差在8%以内，然后提取简化模型的节点、单元、形参、单元类型等模型信息，通过这些信息生成优化分析文件。车身骨架是一个高次超静定的复杂空间杆系结构，各杆件截面形状并不相同，承受的载荷也非常复杂，如果将所有杆件截面参数都选取为设计变量，这是很不现实的。根据前面车身的静力分析得出的计算结果，我们知道扭转工况是车身承受的应力和扭转最严重的工况，该车的刚度基本上达到要求，而强度不足，所以选择扭转工况下，车身骨架应力最高区，中「〕立柱附近区域、顶盖中部区域和车身骨架应力相对较小的地方，后围、<br />

侧围搁梁区域的杆件的截面尺寸参数作为设计变量。选择车身骨架的应力作为状态变量，以车身应力最大的五个点作为应力控制点，保证车身骨架的最大应力值小于材料的许用应力。选取车身重量作为目标函数，通过改变设计变量，在满足车身应力强度的条件下，对车身进行轻量化。由于车身形状比较复杂，精确计算车身<br />

重量比较困难，因此可以通过有限元分析计算单元的重量，然后逐个单元叠加来得到整体车身的重量。<br />

5.3.2计算结果<br />

采用ANSYS软件提供的零阶方法进行了30次迭代优化计算，车身总质量由以前的2169kg减少到2131kg根据市场型材的规格及厂方实际生产条件，对主要杆件优化后的截面尺寸进行了尺寸处理，具体参数见表<br />

对弯扭工况下的车身，取优化后各杆件的截面尺寸，重新计算车身的弯曲应力，车身骨架在弯扭工况下的车身SX方向的应力分布如图5.3所

---