有限元和数值方法-讲稿6-7 弹性力学有限元一般原理和格式.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 改进的矩形单元 引入内部自由度和附加的二次位移模式，构造非协调元 Q4 矩形单元有剪切自锁现象（源于缺少位移二次项，不能准确模拟纯弯曲） Q6 单元中对 Q4 矩形单元进行改进、消除剪切自锁现象的基本思想： * 改进的双线性单元（Q6单元） 补救“剪切自锁”的一个办法是在矩形单元中引入常曲率变形模式，构造Q6单元 其中， 以上位移模式可精确模拟纯弯曲 四个内部自由度（不通过节点与其他单元相联系） 这种位移模式是非协调的——非协调单元 x y 2a 2b * 非协调元 Ｑ6 非协调元，有六个形函数，但是只有四个节点；二个形函数相当在边界的泡泡函数，相应的四个未知数相应于内部自由度，和其它相邻单元无关。相邻单元因此变形不协调，只在单元顶点联结。结构变形时，可以在单元间出现裂缝或叠合 这样的单元比 Q4 柔顺，因此虽然不协调，但给出的结果更好； 其结果从上面趋于真解；（协调元从下面趋于真解） 在特定的载荷下，单元边界间将有缝隙或重叠 * 非协调元 Ｑ6非协调元，有六个形函数，但是只有四个节点；二个形函数相当在边界的泡泡函数，相应的四个未知数相应于内部自由度，和其它相邻单元无关。相邻单元因此变形不协调，只在单元顶点联结。结构变形时，可以在单元间出现裂缝或迭合 x y 2a 2b 这样的单元比Ｑ４柔顺，因此虽然不协调，但给出的结果更好； 其结果从上面趋于真解；（协调元从下面趋于真解） * * * * * * * * * * * * * * * * * * * * * * * * * * * * * (1)结构离散：网格剖分(选用单元类型)； (2)结构分析：形成单元刚度矩阵、节点力向量； (3)集成总体刚度阵和节点力； (4)引入位移边界条件，修改总刚和节点力； (5)求解方程，得到节点位移； (6)计算单元应力，应变； (7)后处理，输出结果。 有限元分析过程 * 位移法有限元分析的步骤 一、结构的离散化 离散化---- 将连续体结构划分成有限个单元，并在单元的指定点设置节点，使相邻单元的有关参数具有一定的连续性。为了使该离散的计算模型有效地逼近实际的连续体，就需要合理选择单元的形状和网格，确定单元和节点的数目、节点的位置和自由度等问题。 1、单元的几何形状 选择单元的形状既要考虑单元的精度和计算量，又要考虑对边界形状的适应性。 * 位移法有限元分析的步骤 2、单元的划分和节点的布置 通常集中载荷的作用点、分布载荷集度的突变点、支承点、边界的转折点等都应取作节点 厚度或者材料突变处应布置节点 单元的大小要根据精度要求和计算机容量等因素综合考虑。对应力变化急剧的区域单元网格应划分得细些，应力变化平缓的区域单元可划分粗些 单元节点的设置应使单元的几何形状能够容易地由节点位置坐标来描述。一些在单元的端点或角点上设置节点，还可在角点之间的边上增设节点，这样可描述单元边界是曲线的情况，用二次或更高次的曲线来描述单元的几何特征。 * 位移法有限元分析的步骤 3、节点自由度数的确定 节点的自由度数应能正确地模拟单元的力学特征和易于构造单元的近似解，保证相邻单元在公共边界上的连续性。 综合起来，结构离散化应符合几何近似和物理近似的原则。有限元的计算模型在几何上应近似于真实结构。各单元的物理特性和受力特性应与真实结构在单元所在区域的情况近似。 * 位移法有限元分析的步骤 二、选择位移模式 在对典型单元进行特性分析时，为了能用节点位移表示单元内各点的位移、应变、应力，必须对单元内的位移分布作出一定的假设，也就是在每个单元内假设位移是坐标的某种简单函数，将单元内的位移场表示成单元节点位移的插值函数形式，称为位移模式。用矩阵形式写成 u —单元内位移函数列阵? ? de —单元节点位移列阵? N —称形状函数矩阵，它的元素是位置坐标的函数，称位移插值函数或形函数。 * 位移法有限元分析的步骤 三、单元的力学特性分析 1、单元的应变与单元节点位移的关系 L —微分算子矩阵，例如对于二维问题 * 位移法有限元分析的步骤 2、单元内任一点应力与单元节点位移的关系 D —弹性模量矩阵。对于二维平面应力问题 3、单元节点力与节点位移之间的关系 将单元视为独立的求解域，不计体力并考虑单元只在节点处受力时，线弹性体的单元势能泛函为 * 位移法有限元分析的步骤 由于单元处于平衡，根据最小势能原理，有 此即单元节点力与节点位移之间的关系 以上导出了单元刚度矩阵的普遍表达式。如果计及单元的体力和考虑到单元的边界有面力作用的情况，还可求得单元等效节点载荷 * 位移法有限元分析的步骤 1、