有限元法在内燃机上的应用.ppt

* * 矩形单元应变 位移函数：4个节点，4个系数 矩形单元应变。因为形函数为2次函数，求导后几何矩阵[B]元素为一次函数。所以矩形单元应变是一次函数 矩形单元应力 矩形单元应力（物理方程） 矩形单元应力。因为形函数为2次函数，求导后几何矩阵[B]元素为一次函数。所以矩形单元应力是一次函数 三角形单元应变、应力都是常量，矩形单元应变和应力都是一次函数，精度高于三角形单元，可以采用较粗的网格就可以得到较高的精度，而且求解量小 平面问题虚功方程 平面问题虚功方程（作为定律直接应用），t为平面厚度 代入矩形单元的应变 得到 矩形单元刚阵（1） 矩形单元刚阵 考虑到虚位移可以是任意量，提到积分号外 考虑到节点位移也是任意的 矩形单元刚阵（2） 考虑到几何矩阵元素为常量，应力也是常量 对比以下公式，矩形单元刚阵已经就得 矩形单元刚阵需要积分，几何矩阵[B]中元素为一次函数，积分方便，这也是单元位移函数采用多项式的原因 矩形单元应力 矩形单元应力（物理方程） 矩形单元应力。因为形函数为2次函数，求导后几何矩阵[B]元素为一次函数。所以矩形单元应力是一次函数 三角形单元应变、应力都是常量，矩形单元应变和应力都是一次函数，精度高于三角形单元，可以采用较粗的网格就可以得到较高的精度，而且求解量小 矩形单元有限元法（1） 1、如右图平面问题，划分网格， 得到节点1、2、3、4、5、6和单元1、2 2、求出单元1和单元2的单元刚阵（8x8矩阵） 3、将单元1和2的单元刚阵扩充相加，得总体刚阵。（12X12矩阵） 4、施加边界条件 A、节点1、2固定，位移为0，不用再求解，在总体刚阵中消去对应的行列（考虑到节点1、2都有x、y2个方向，所以消去4行4列，得实用总刚阵（8X8矩阵） B、由节点平衡方程，得知节点4x方向的节点力为-P，其他节点力=0， 可列出节点力和节点位移的线性代数方程组（8个方程） 矩形单元有限元法（2） 求解线性代数方程组，得到各节点位移（实际上只有节点3、4、5、6的x、y向位移，共8个，节点1、2位移为0） 节点位移已知后，平板上任一点的位移有位移函数求得(近似） 知道位移后，平板上任一点应变有几何方程求得，对位移求导，常应变。 平板上任一点应力有物理方程求得 所有量都已经求得 矩形单元有限元法（3） 可以只用1个矩形单元，4个节点，和三角形单元节点数相同，但，精度高于三角形单元 矩形单元的局限性 矩形单元不适于划分曲边边界处，需要借用三角形单元 任意四边形单元适于划分曲边边界，不需要借用三角形单元 单元的变形协调性要求（1） 相邻的两个单元变形后在公共边处既不能开裂，也不能重叠，这称为单元变形的协调性，这样的单元称为协调单元 如果单元的位移函数在公共边上退化为一次函数（或低于一次），则单元变形协调 1、三角形单元位移函数本身就是一次函数，所以单元变形协调 2、矩形单元位移函数为2次函数，但是在如下坐标系下，最右边 方程为x=a，在该边上位移函数退化为一次函数 单元的变形协调性要求（2） 3、任意四边形单元，公共边的位移函数仍然为2次函数，单元变形后，在相邻的公共边处会开裂或重叠，任意四边形变形不协调 4、任意四边形采用自然坐标系可以解决单元变形协调性问题， 内容略 第四章 体单元简介 三种形状体单元 四面体单元(棱锥） 五面体单元（楔形单元），只用于四面体和六面体的过度 六面体单元（砖形单元） 内容略 4节点四面体单元位移函数 4个节点，4个未知系数 位移函数为一次函数，类似3节点三角形单元，可以推知其单元的应变和单元应力都是常量，即4节点四面体单元是常应变、常应力单元，精度较差 类似可得到4面体单元刚阵。（略） 四面体单元法可以求解所有关于体的应力、应变和位移问题（略） 高阶体单元 10节点四面体单元 8节点六面体单元 20节点六面体单元 四面体与六面体单元的选择 4节点四面体单元常应变、常应力，精度差。细化网格才能提高精度，但会使得求解量增加 10节点四面体单元应变、应力不是常量，精度高，但是求解量大。每个节点有3个方程 8节点六面体单元不需要细化网格即可获得较高精度，求解量不大。 四面体单元的最大优点是FEA软件自动划分网格容易，而六面体单元较难自动划分网格，投入的人力大大增加 建议:静态应力计算可用10节点的4面体网格。 动态应力计算和迭代计算量大的采用六面体网格 第五章 平面温度应力问题简介 平面温度场的控制方程 控制方程 平面温度场温度边界条件 1、第一类边界条件 2、第二类边界条件 3、第三类边界