CADCAM技术基础-计算机辅助工程分析-文档资料

1、12010. 09CAD/CAMCAD/CAM技术基础技术基础CAD/CAM Technology BaseCAD/CAM Technology Base2第五章第五章 计算机辅助工程分析计算机辅助工程分析 Computer Aided EngineeringComputer Aided Engineering3引引 例例 有限元法、边界元法及结构优化设计技术等计算力学方法 4引引 例例过程概 念 设 计产 品 设 计模 具 设 计开 模试 模生 产修 改5引引 例例计算机辅助工程概 念 设 计产 品 设 计模 具 设 计开 模试 模生 产6内内 容容 75.1.1 有限元分析方法概述有限元分

2、析方法概述 l Finite Element Analysis l是一种用于求解各类工程实际问题的数值计算方法 ，广泛应用于解决力学问题和场问题 。l简称为FEA8有限元分析的思路和方法：有限元分析的思路和方法： 5.1.1 有限元分析方法概述有限元分析方法概述 9l将某个工程结构离散为由各种单元（每种单元可以是一维、二维或三维的情况）组成的计算模型，这一步称作单元剖分。l离散后单元与单元之间利用单元的节点相互连结起来。l单元节点的设置、性质、数目等应根据问题的性质、描述变形形态的需要和计算精度而定(一般情况，单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大)。 5.1.1 有

3、限元分析方法概述有限元分析方法概述 10l 选择位移模式：可选择节点位移作为基本未知量；可选择节点力作为基本未知量；亦可取一部分节点力和一部分节点位移作为基本未知量。l 分析单元的力学性质：根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式。l 计算等效节点力：用等效的节点力来替代所有作用在单元上的力。5.1.1 有限元分析方法概述有限元分析方法概述 115.1.1 有限元分析方法概述有限元分析方法概述 l利用结构力的平衡条件和边界条件把各个单元按原来的结构重新联接起来，形成整体的有限元方程 (5.1)K是整体结构的刚度矩阵；q是节点位移列阵；F是载荷列

4、阵。FKq 125.1.1 有限元分析方法概述有限元分析方法概述 l解有限元方程式得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。l有限元法的基本思想是“一分一合”，“分”是为了进行单元分析，“合”则是为了对整体结构进行综合分析。135.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法l直接方法。l虚功原理法。l能量变分原理方法14 图中所示是xoy平面中的简支梁弯曲简图，EI为梁的抗弯刚度。现在，以它为例用直接方法建立单元的刚度矩阵。图5.1 平面简支梁和它的计算模型5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法15

5、 梁在横向外载荷(可以是集中力或力矩或分布载荷等)作用下产生弯曲变形，对于平面弯曲问题，每个点(包括支承点)处的位移有两个，即挠度和倾角；相应地也有两个节点力，即与挠度对应的剪力和与倾角对应的弯矩。我们规定挠度和剪力向上为正，倾角和弯矩逆时针方向为正。 为使问题简化，把图示的梁看成是一个单元。当令左支承点为节点i，右支承点为节点j时，则节点位移和节点力可以分别写成Vi、zi、Vj、zj和Fyi、Mzi、Fyj、Mzj。也可写成矩阵形式 zjjziivvq zjyjziyiMFMFF 5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法16 梁的节点力和节点位移是有联

6、系的。在弹性小位移范围内，这种联系是线性的。zjjziizjyjziyivvkkkkkkkkkkkkkkkkMFMF44434241343332312423222114131211即 qKF它代表了单元的载荷与位移之间(或力与变形之间)的联系，称为单元的有限元方程。式中 K称为单元刚阵，它是单元的特性矩阵。（5.2）5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法17 单元刚度矩阵中任一元素kij表示j号节点的单位位移对i号节点力的贡献。如 K中第1列各元素就分别代表当在i节点处挠度方向产生单位位移i=1时，它们对其他各位移(包括vi)方向上引起的节点力Fyi，

7、Mzi，Fyj，Mzj的贡献。由功的互等定理有kij = kji，所以单元刚度矩阵是对称的。对于梁单元平面弯曲问题，可以计算出各系数kij的数值。5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法18 假设0, 1zjjziivv则梁变形情况如图所示：1zi0zjjivv图5.2 梁变形图图5.3 梁变形图5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法19 这时，平面弯曲梁单元的刚度矩阵或单元特性矩阵为： 2222346266126122646612612lllllllllllllEIK（5.3）5.1.2 有限元分析方法中单元特

8、性的导出方法有限元分析方法中单元特性的导出方法20 解题基本步骤：l设定位移函数l由位移函数求应变l根据虎克定律，通过应变求应力l由虚功原理求单元的刚度矩阵5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法21 以平面问题中的三角形单元为例，说明其方法步骤。图5.4 三角形单元5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法22 Tyxvyxuyxd,当单元很小时，单元内一点的位移可以通过节点的位移数值来表示。假设单元内位移为x, y的线性函数，即：yaxaayxu321,yaxaayxv654,设三节点三角形单元内的位移函数为：5

9、.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法23 写成矩阵形式： aSaaaaaayxyxvud65432110000001（5.4）5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法24 同理，单元的三个节点i，j，k上的位移也可用它来表示： acaaaaaayxyxyxyxyxyxvuvuvuqkkkkjjjjiiiikkjjii654321100000011000000110000001即 qca15.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法25 其中 kjikjikjikjikjikji

10、cccbbbaaacccbbbaaaAc000000000000000000211式中A是三角形面积5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法26 kjjkjkjikkjjiiyyxxyyxxyxyxyxA1112jkkjiyxyxakiikjyxyxaijjikyxyxakjiyybikjyybjikyybjkixxckijxxcijkxxc（5.5）5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法27 kkjjiikjikjikjikjikjikjivuvuvucccbbbaaacccbbbaaayxyxAvu00000

11、000000000000010000001215.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法28 上式相乘后，可写成kkjjiikkjjiivNvNvNyxvuNuNuNyxu,也可简写成 qNd 此式即为单元内某点的位移用节点位移插值表示的多项式。称N为形状函数。（5.6）（5.6a）5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法29 由弹性力学知 qBvuvuvubcbcbccccbbbAkkjjiikkjjiikjikji00000021（5.7）5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导

12、出方法30 根据虎克定律，对于平面问题，有 qBDD其中 D对于平面问题为 2100010112ED（5.8）（5.9）5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法31 根据虚功原理，当结构受载荷作用处于平衡状态时，在任意给出的节点虚位移下，外力(节点力)F及内力所做的虚功之和应等于零，即0AAF5.1.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法32 经推导得平面应力问题三角形单元刚度矩阵为：srsrsrsrsrsrsrsrrsbbcccbbcbcbbccbbAEtK21212121142kjiskjir,;,（5.16）5.1

13、.2 有限元分析方法中单元特性的导出方法有限元分析方法中单元特性的导出方法335.1.3 有限元法的解题步骤1 单元剖分和插值函数的确定根据构件的几何特性、载荷情况及所要求的变形点，建立由各种单元所组成的计算模型。再按单元的性质和精度要求，写出表示单元内任意点的位移函数u(x,y,z)，v(x,y,z)， w(x,y,z)或d =S(x,y,z) a。 qNqCSd1上式是用节点位移表示单元体内任意点位移的插值函数式。345.1.3 有限元法的解题步骤2 单元特性分析根据位移插值函数，由弹性力学中给出应变和位移关系，计算出应变。由物理关系，得应变与应力的关系式。由虚位移原理，可得单元的有限元方

14、程，或力与位移之间的关系式，即 qKF式中，K是单元特性，即刚度矩阵355.1.3 有限元法的解题步骤3 单元组集把各单元按节点组集成与原结构相似的整体结构，得到整体结构的节点与节点位移的关系：KqF 式中，K是整体结构的刚度矩阵；F是总的载荷列阵；q是整体结构所有节点的位移列阵。4 解有限元方程可采用不同的计算方法解有限元方程，得出各节点的移。5 计算应力36 5.2.1 机械优化设计问题的数学模型1 设计变量一个设计方案可以用一组基本参数的数值来表示。对某个具体的优化设计问题，并不是要求对所有的基本参数都用优化方法进行修改调整。一些参数，可以根据已有的经验预先取为定值。这样，对这个设计方案

15、来说，它们就成为设计常数。而除此之外的基本参数，则需要在优化设计过程中不断进行修改、调整，这些基本参数称为设计变量，又叫做优化参数。37 5.2.1 机械优化设计问题的数学模型1 设计变量设计变量的全体实际上是一组变量，可用一个列向量表示Tnxxxx 21称作设计变量向量。向量中分量的次序完全是任意的可以根据使用的方便任意选取。一旦规定了这样一种向量的组成，则其中任意一个特定的向量都可以是一个“设计”。 38 5.2.1 机械优化设计问题的数学模型2 约束条件一个可行设计必须满足某些设计限制条件，这些限制条件称作约束条件，简称约束。在工程问题中，根据约束的性质可以把它们分成性能约束和侧面约束两

16、大类，针对性能要求而提出的限制条件称作性能约束。不是针对性能要求，只是对设计变量的取值范围加以限制的约束称作侧面约束。39 5.2.1 机械优化设计问题的数学模型3 目标函数倘若需要优化的性质可以表示成设计变量的一个可计算函数，这个用来使设计得以优化的函数就称作目标函数。建立目标函数是整个优化设计过程中重要的问题。当对某一个性能有特定的要求，而这个要求又很难满足时，则若针对这一性能进行优化将会取得满意的效果。但在某些设计问题中，可能存在两个或两个以上需要优化的指标，这将是多目标函数的问题。40 5.2.1 机械优化设计问题的数学模型4 优化问题的数学模型优化问题的数学模型是实际优化设计问题的数

17、学抽象。在明确设计变量、约束条件、目标函数之后，优化设计问题就可以表示成一般数学形式。在实际优化问题中，对目标函数一般有两种要求形式：目标函数极小化f(x)min或目标函数级大化f(x)max。由于求f(x)的极大化与求-f(x)极小化等价，所以优化问题的数学表达一般采用目标函数极小化形式。41 5.2.2 机械优化设计问题的基本解法1 解析解法与数值解法解析解法就是把所研究的对象用数学方程(数学模型)描述出来，然后再用数学解析方法(如微分方法)求出优化解。但是，在很多情况下，优化设计的数学描述比较复杂因而不便于其至不可能用解析方法求解。数值解法不仅可用于求解复杂函数的优化解，也可以用于处理没

18、有数学解析表达式的优化设计问题。因此，它是实际问题中常用的方法。42 5.2.2 机械优化设计问题的基本解法2 优化准则法与数学规划法优化准则法是从一个初始设计xk出发，着眼于在每次迭代中满足的优化条件，按着迭代公式kkkxCx1得到一个改进的xk+1，而无需再考虑目标函数和约束条件的信息状态数学规划法是从一个初始设计xk出发，对结构进行分析，但是按照如下迭代公式kkkxxx1得到一个改进的设计xk+1（5.20）（5.21）43 5.2.2 机械优化设计问题的基本解法3 迭代终止条件准则1：当相邻两设计点的移动距离已达到充分小时，即：11kkxx准则2：当函数值的下降量已达到充分小时，即：3

19、1)(kkxfxf准则3：当某次迭代点的目标函数梯度已达到充分小时，即： 5kxf采用哪种收敛准则，可视具体问题而定。44 例题：平面四连杆机构的优化设计。 图5.5所示是一个曲柄摇杆机构。图中x1,x2,x3,x4分别是曲柄AB、连杆BC、摇杆CD和机架AD的长度。是曲柄输入角，0是摇杆输出的起始位置角。这里，规定0为摇杆的右极限位置角为0时的曲柄起始位置角，它们可以由x1,x2,x3,x4确定。通常规定曲柄长度x1=1，而在这里x4是给定的，并设x4=5，所以只有x2和x3是设计变量。图5.5 曲柄摇杆机构45 例题：平面四连杆机构的优化设计。在给定最大和最小传动角的前提下，当曲柄从0位置

20、转到0+90时，要求摇杆的输出角最优地实现一个给定的运动规律f0()。要求 200032 f对于这样的设计问题，可以取机构的期望输出角=f0()和实际输出角j=fj()的平方误差积分准则作为目标函数，使 2200dxfj最小。46 例题：平面四连杆机构的优化设计。当把输入角取s个点进行数值计算时，目标函数可以化简为 sijiixxfxf0243,最小相应的约束条件有：（1）曲柄与机架共线位置时的传动角最大传动角 max135 最小传动角 min45对本问题可以计算出47 例题：平面四连杆机构的优化设计。322322max236arccosxxxx322322min216arccosxxxx所以有01645cos2036135cos2322322322322xxxxxxxx（2）曲柄存在条件321441