第六章计算机辅助工程分析_机械CAD_CAM技术.pptxVIP

第六章 计算机辅助工程分析;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、求解线性代数方程组 求解线性代数方程，得到节点位移的未知数 8、结果 得到节点位移的未知数后，运用几何方程和物理方程，得出单元应力应变，输出结果或用图线来表示。;四、有限单元方法的前处理、后处理 有限单元方法的前处理、后处理在有限单元方法中的应用和作用。;图6.11显示正确的和错误的前处理后的单元网格图;前处理时的注意事项： a、单元一般均匀划分，但在几何突变的地方加密。 b、节点编