有限元总复习预案.ppt

梁单元 平面三角形单元 空间问题轴对称问题动力学问题 等效节点载荷问题 每个单元的来计算、再集成到一起 分别计算集中力，分布力，集中力矩等效到每个单元上 杆单元、梁单元、平面单元 * * * * * * * * * * 位移函数： 等参单元-任意直边四边形单元（4节点） * * 平面问题有限元分析-等参单元 坐标变换函数与位移函数采用相同的形状函数 等参单元-任意直边四边形单元（4节点） 对比 节点位移个数与节点坐标个数一样，在选取单元的位移函数时，可取单元的自然坐标作为自变量，取图形（坐标）变换式的形状函数为位移函数的形状函数。 * * 平面问题有限元分析-等参单元 （i=1，2，…，4） 等参单元刚度（4节点） 根据形函数，几何矩阵为 * * 平面问题有限元分析-等参单元 在上式中， 是ξ、η的函数，因此必须用坐标变换式来转换导数关系，根据复合函数的求导法则，有 （i=1，2，…，4） 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 式中的J称为 雅可比（Jacobi）矩阵 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 求出 关键 将单元应变代入平面问题的物理方程式，就得到平面4节点等参单元的应力列阵，为 （i=1，2，…，4） 等参元刚度矩阵仍可用下式表示 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 由以上分析可知，在划分单元时，只需确定单元节点的整体坐标值，而不必画出等参元的具体形状，因为在计算中实际使用的只有单元四个节点在整体坐标系下的位移值。 等参元变换的条件为 ，因此在有限元网格划分时，要特别注意这一点。 由此式可计算出单元的刚度矩阵，此处的积分一般得不到显式形式，因此要采用数值积分法（如Gauss积分法），该积分对任何形状的四边形单元具有相同的积分表达式和相同的积分限。上式积分运算是在母单元（自然坐标）内进行的，但节点力与节点位移的方向均沿整体坐标x，y方向。所以[k]e是整体坐标内的刚度矩阵。 等参单元刚度（4节点） * * 平面问题有限元分析-等参单元 被积函数f一般很复杂，往往不能得出它的显式，在有限元计算中常采用数值积分。在单元内选出某些积分点(2～10)，算出被积函数f在这些积分点的值，然后用加权系数乘上这些积分值，再求出总和作为近似的积分值 高斯积分法用同样数目的积分点具有较高精度 一维高斯积分公式 式中,n为积分的点数，?i为积分点i的局部坐标 ，Wi为加权系数 ，f(?i)为被积函数在积分点处的函数值 高斯积分法 * 平面问题有限元分析-等参单元 * * 平面问题有限元分析-等参单元 * * 平面问题有限元分析-等参单元 * 积分点数的选取与被积函数f(?)有关，当f(?)是m次多项式时，取 当f(?)不是多项式时，需要通过一些试算来判断选取适当的n值。积分点数n(2～10)不能取得过大，否则将随着积分点的增多计算量急剧增加 高斯积分法 高斯求积公式的积分点坐标与加权系数表 n ξi Wi 2 ±0.57735 1 3 0 ±0.77460 8/9 5/9 4 高斯积分法 实际上作用在物体上的外力可以直接作用在节点上，也可以不作用在节点上，然而在有限元法中要求作用在物体上的各种外力必须用作用在节点上的力表示，这一过程称为外力的静力等效移置，所得到的节点力称为等效节点力。 处理思路 * 1 杆单元 下图示出了一维铰接杆单元，横截面积为A，长度为l，弹性模量为E，轴向分布载荷为px。单元有2个结点i，j，单元坐标为一维坐标轴x。 · · i j x l LINK px uj ui 单元结点力向量： 单元结点位移向量： * 等效节点力 单元上作用分布力px，则等效节点力计算公式仍为以下形式 当分布力集度px为常数时，有 * 集中力引起： 其中：Pj —— 沿杆轴线方向集中力大小 1.无轴向变形的平面梁单元 如图所示的梁单元为一个无轴向变形的等截面直杆，共有两个节点，节点位移包括挠度和转角，节点力包括剪力和弯矩。 2 梁单元 * 单元的节点位移向量表示为 节点力向量表示为 设单元的挠度 的位移模式取为 将单元上两个节点的坐标和相应的位移代入上式，可得单元的位移模式为 形函数表示为 * 等效节点力 如果梁上作用有集中力或集中力偶时，在划分单元时，一般将载荷作用点取为节点，在整体载荷列阵中进行叠加。如果梁上作用有横向分布载荷，则等效节点力可由下式计算 将形函数代入上式即可求得等效节点