第六章节计算机辅助工程分析解析_机械cad_cam技术.pptVIP

第六章 计算机辅助工程分析 本章主要介绍CAD/CAM系统中的计算机辅助工程分析方法，有限元、优化设计和仿真技术等内容。 6.1有限元（Finite Element Method） 有限元方法基本思路：将连续体离散为有限个单元，单元间用节点相连，将受力等效到节点上，根据受力平衡条件建立方程，将所有单元重新组合，连立方程，并求解，这样的方法称为有限元。 力学分析方法类型： 1、解析法：根据数学方法，分析计算典型的受力情况，并得出经典的理论解，该方法只能用于计算受力、结构都比较简单，能用数学方法进行求解的简单而典型物理问题。 2、数值计算方法：采用数值计算原理和方法，将复杂的物理问题简化为简单的问题，使用计算机进行求解，该方法可以应用于结构、受力等复杂，又无法用经典的数学物理方法进行求解的物理问题。 有限单元方法和有限差分方法为两个工程上最为常用的数值计算方法，本教材简单介绍有限单元方法。 有限单元方法采用的单元类型有很多，根据不同的物理问题可以分为图6.2所示的集中常用单元类型： 一、弹性力学的基本知识： 1、常用物理量 (1)外力：分集中力（N）、线载荷(N/m)、面载荷(N/m2)和体积力(N/m3) (2)应力:分正应力和剪应力，单位： N/m2，见图6.3 (3)应变：分正应变(单位长度的变形量)和剪应变(变形前后的线或面之间的夹角变化量)，见图6.4 (4)位移：在力或温差等载荷作用下，物体内点的位置变化，分正位移和角位移。 2、力学方程 (1)应变与位移间的关系（几何方程）： 应变写成阵列形式： (2)应力应变关系（物理方程）： 写成矩阵形式： 二、有限元方法应用举例： 如图所示：等截面直杆，上端固定，下端自由，设杆的截面积A，长为L，单位杆长的重力q，试采用有限元方法计算杆中各点的位移。 1、将杆离散化有限个杆单元（图6.6a） 2、计算单元中任意位置的位移 取三个单元中任一单元，e为单元号，i,j为节点号，建立局部坐标系，假设单元中的位移呈线性分布，那么得到式(6.3) 式中系数a1,a2待定。 设xi=0时，u=ui，xj=l时，u=uj，代入方程6.3解得系数a1，a2 将系数a1，a2代回公式6.3，记f=u，得到： 设：N=（Ni，Nj）= ，称为形状函数，因为函数只与l、与I点的距离有关。 Δe= ui 可以得到： uj f=N Δe (6.4) 3、等效节点载荷 根据力学上的虚功原理，将单元载荷等效到单元的节点i，j上。 计算得到的各单元的节点载荷分别为：（计算过程省略） (6.5) 4、建立单元刚度矩阵 (1) 应变与位移之间关系： 几何方程：ε=du/dx，代入公式6.4，得到： (6.6) (2)应力与位移之间关系： 物理方程： σ =Eε，代入公式6.6 ，得到： (6.7) (3) 建立力与位移之间关系： 由虚功原理得到力与位移之间关系： (6.8) 式中K为单元刚度矩阵，将公式(6.6)和(6.7)中的系数代入公式(6.8)得到： (6.9) 5、建立总刚度矩阵 对于单元2、3，同理建立 (6.8)类型的方程 根据节点的力平衡，得到： 整理上面公式，得到整个问题的刚度矩阵： 记： 可以得到： 6.12 式中：K----总体刚度矩阵 Δ----总体节点位易矩阵 F----总体节点载荷矩阵 计算各矩阵中的系数，得到： 加入边界条件：u1=0，方程变成： 求解方程 三、有限元方法的基本算法与步骤 1、单元划分： 将一个连续体分割离散成许多个一定大小的单元，并将单元和节点进行编号。 2、单元特征分析 以节点位移为未知量，设一个单元位移函数（以节点位移表示单元位移）。 3、建立单元刚度矩阵 使用几何方程以节点位移表示单元应变 使用物理方程以节点位移表示单元应力 使用虚功方程以节点位移表示节点力，建立单元刚度矩阵 4、建立总刚度矩阵 通过各单元刚度矩阵形成总刚度矩阵 5、建立载荷矩阵 通过节点的力平衡方程，建立节点力矩阵，从而得到以节点位移为未知数的，总刚度矩阵为系数的线性代数方程。 6、边界条件的处理 应用位移约束或已知的边界条件 7、求解线性代数方程组 求解线性