第五章—结构单元(1).pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 有限单元法 崔向阳 机械与运载工程学院 湖南大学 College of Mechanical Vehicle Engineering Hunan University 崔向阳 结构单元(1)-梁问题 第五章: * 引言 梁板壳结构在工程中广泛应用。 特点： 梁——两个方向的尺度比其它方向小得多 板壳——一个方向的尺度比其它方向小得多 梁板壳结构在分析时可以根据其特性进行一定的简化。简化后仍然包括三大类基本方程和两类边界条件，只是表达形式一般与通用表达式有所不同。 板 壳 梁 * 引言 原则上，可以使用2D、3D实体单元分析梁（2D）和板壳结构（3D）问题，但存在一定的困难。 为了获得一定的计算精度，单元划分时必须保持单元在各个方向上尺度相近，这样导致单元总数过分庞大，计算效率过低。 关于梁板壳结构，通常是根据结构的特点在应变和应力方面引入一定的假定，对问题进行简化，从而构造适合梁板壳结构分析的单元。 结构单元是梁板壳单元的总称。 * 引言 基于Kirchhoff假设的经典梁单元 ——不考虑剪切变形的薄梁单元 考虑剪切变形的梁单元 经典梁理论基础上引入剪切变形 ——截面转动和挠度仍然相关（C1） Timoshenko梁单元 ——挠度和截面转动独立插值（C0） * 薄梁公式 平衡方程 (静态) (动态) 单元位移 对于薄梁： * 薄梁单元 设单元坐标位移模式为 位移模式 由单元两端点的节点位移条件，解出a1、a2、a3、a4。再代入该式，可将位移模式写为以下形式： * 薄梁单元 应变矩阵B 单元弯曲应变?b与节点位移de的关系。 由材料力学知，梁单元上任一点的应变和该点挠度之间关系为： 所以，可得单元弯曲应变和节点位移之间关系 因为 * 薄梁单元 刚度矩阵 * 薄梁等参单元 自然坐标系: 位移模式： * 薄梁等参单元 在 x= -a 或 x = -1 在 x= a 或 x = 1 * 薄梁等参单元 * 所以 薄梁等参单元 * 积分得到 薄梁等参单元 * 积分得到 薄梁等参单元 * 假设梁单元沿x方向有均布外力fy，两端点分别有集中力fs1， fs2和集中力矩ms1,ms2，则总节点力向量为： 薄梁等参单元 * 剪切梁 只有当梁的高度远小于跨度时（一般l/h5为细长梁，其计算误差满足工程精度 ），才能忽略横向剪切变形的影响。而高梁的情况下，梁内的横向剪切力将产生剪切变形并引起梁的附加挠度，使原来垂直于中面的截面变形后不再和中面垂直，但仍假定截面保持平面。 薄梁单元基础上引入剪切变形的梁单元（C1） 挠度和转角独立插值的Timoshenko梁单元（C0） 修正系数 * 剪切梁 薄梁单元基础上引入剪切变形的梁单元 对于弯曲引起的挠度，采用薄梁的三次多项式插值函数： 参见经典梁的表达形式 * 剪切梁 位移插值函数 节点位移条件 对于剪切引起的附加挠度，采用线性插值函数： 剪应变 有兴趣的同学可以根据 推导单元刚度矩阵 此处不做讲解！！！！！ * 剪切梁 挠度和转角独立插值的Timoshenko梁单元 采用线性插值函数，可以得到： 剪应变 梁的曲率 * 剪切梁 令 * 剪切梁 当 h/l 趋于0时（即梁很薄）时，希望剪应变为零 要使该式在梁单元内恒成立，不仅常数项为零，还必须一次项为零，因此要求 从而 这意味着梁不能发生弯曲，与真实情况相违背，这种现象称为剪切锁死(shear locking)。 * 剪切梁 在剪切应变的表达式中， 和 的函数表达式不是相同的阶次，无法恒满足细长梁的约束条件 换句话说，在梁很薄的情况下，不适当地夸大了剪切应变能的量级造成了剪切锁死现象。 克服“剪切锁死”，可以采用的方案 减缩积分（reduced integration） 假设剪切应变（assumed shear strains） * 剪切梁 在计算 的积分时，不采用精确积分，而用一点积分（单元的中心）来计算，这样相当于将原来的线性变化关系改为常数(中点平均值)，使得 和 保持同阶，就有可能做到使细长梁的约束条件 恒得到满足。 如此，考虑剪切变形的Timoshenko梁单元也可以用于细长梁的分析。 * 剪切梁 精确积分 减缩积分 * 薄梁例题 P=1000 N 0.5 m 0.06 m 0.1 m E=69 GPa ?=0.33 * 薄梁例题 第一步: 求单元刚度矩阵 * 薄梁例题 第二步: 施加边界条件 问题的平衡方程表示为： * 薄梁例题 第三步：求解