有限元与数值方法-讲稿6.pptVIP

State Key Laboratory of Structural Analysis for Industrial Equipment Dalian University of Technology, China 5. 2. 平面问题的三角形单元： （1）三节点三角形单元 三节点三角形平面单元的形函数 三角形的面积坐标 三角形的面积坐标 用面积坐标表示的形函数 形函数的性质及位移插值 单元应变矩阵 单元应变矩阵 单元应力矩阵 Constant Strain Triangle Stiffness matrix for element k =BTEB tA The CST gives good results in regions of the FE model where there is little strain gradient Otherwise it does not work well. 单元刚度矩阵 Linear Strain Triangle Will this element work better for the problem? 有限元与数值方法第八1讲 5.3 弹性平面问题四边形单元 5.3 双线性矩形（Q4）平面单元 Bilinear Quadratic The element stiffness matrix is obtained the same way A big challenge with this element is that the displacement field has a bilinear approximation, which means that the strains vary linearly in the two directions. But, the linear variation does not change along the length of the element. 双线性矩形（Q4）平面单元 单元缺陷-剪切自锁 单元缺陷-剪切自锁 改进的双线性单元（Q6单元） 作业题 上述单元位移模式满足位移模式选择的基本要求： 反映了单元的刚体位移和常应变 单元在公共边界上位移连续 在矩形单元的边界上，坐标x和y的其中一个取常量，因此在边界上位移是线性分布的，由两个结点上的位移确定 与三结点三角形单元相比，四结点矩形单元的位移模式是坐标的二次函数，能够提高计算精度，但也有显著的缺点，两种单元的比较如下 不能适应曲线边界和非正交的直线边界 单元内的应力、应变是线性变化的，计算精度较高 四结点矩形单元 计算精度低 适应复杂形状， 单元大小过渡方便 三结点三角形单元 缺点 优点 单元类型 双线性矩形单元(Q4)不能表示纯弯曲 纯弯曲中，材料的真实应变应为 而Q4单元被弯曲时，顶边和底边仍是直的，单元位移和应变为 即使在纯弯载荷下，单元也产生附加的剪切变形。因此单元对于弯曲变形来说过刚 即Q4单元在模拟了弯曲变形的同时自动产生了所谓“寄生剪应变”，从而吸收了一定的应变能 寄生剪切将产生“剪切自锁”，特别是当单元纵横比较大时明显。 不建议采用细长的Q4单元 如果 ， 则产生的转角比为 1 2 3 4 x y M1 M2 a b 若纵横比增加，则单元弯曲变形将趋于零 改进的矩形单元 引入内部自由度和附加的二次位移模式，构造非协调元 Q4 矩形单元有剪切自锁现象（源于缺少位移二次项，不能准确模拟纯弯曲） Q6 单元中对 Q4 矩形单元进行改进、消除剪切自锁现象的基本思想： 补救“剪切自锁”的一个办法是在矩形单元中引入常曲率变形模式，构造Q6单元 其中， 以上位移模式可精确模拟纯弯曲 四个内部自由度（不通过节点与其他单元相联系） 这种位移模式是非协调的——非协调单元 x y 2a 2b 非协调元 Ｑ6 非协调元，有六个形函数，但是只有四个节点；二个形函数相当在边界的泡泡函数，相应的四个未知数相应于内部自由度，和其它相邻单元无关。相邻单元因此变形不协调，只在单元顶点联结。结构变形时，可以在单元间出现裂缝或叠合 这样的单元比 Q4 柔顺，因此虽然不协调，但给出的结果更好； 其结果从上面趋于真解；（协调元从下面趋于真解） 在特定的载荷下，单元边界间将有缝隙或重叠 非协调元 Ｑ6非协调元，有六个形函数，但是只有四个节点；二个形函数相当在边界的泡泡函数，相应的四个未知数相应于内部自由度，和其它相邻单元无关。相邻单元因此变形不协调，只在单元顶点联结。结构变形时，可以在单元间出现裂缝或迭合 x y 2a 2b 这样的单元比Ｑ４柔顺，因此虽然不协调，但给出的结果更好； 其结果从上面趋于真解；（协调元从下面趋于真解） 简答题 加权残数