数值模拟：第四讲平面问题(一)——基本知识与离散化精要.pptVIP

4 弹性力学基本知识与平面问题离散化 (3) 弹性力学问题的求解策略 1）解析法—精确解 a. 应力解法 b. 位移解法 2）能量法（变分法）—近似解 3）数值法—近似解 a.有限差分法 b.有限元法 c.边界元法 弹性体受力变形后任意一点都产生位置的变化，形成一个位移矢量，其在坐标轴上的投影（位移分量）用 u，v，w 表示。 4.2 弹性力学平面问题基础 4.2.1 平面应力问题___高速旋转的薄圆盘、土建中的剪力墙 4.2.2 平面应变问题——水坝，火炮身管 4.2.3 平面问题基本方程和边界条件 实体模型离散化为有限元模型 4.3.4 边界条件离散化的原则 当结构离散化后，原先连续的边界被多个单元的边界所形成的折线所代替，这时边界位移和边界外载荷等边界条件也需作相应的离散化处理，并规定如下： （1） 原结构某段边界在某个方向上受到约束或可认为某个方向上的位移忽略不计，则离散此段边界后需在相关边界节点的该方向上加上一个连杆支座约束，即单向约束。 （2）原结构某段边界的位移为零，则离散此段边界后需在相关边界节点上加上一个固定铰链支座约束，即双向约束。 （3） 所有作用在结构上的外载荷，包括集中力、面力和体力(或重力)，都要按静力等效或虚功等效的原则移到有关单元的就近节点上而成为节点等效载荷 载荷离散化 （1）单元分类 在弹性力学平面问题中用以进行结构离散化的单元形式有许多种，这样就构成了多种单元的形式。 1）从单元的几何形态来分 三角形、矩形、任意四边形等； 2）从单元的边界形状来分 直边和曲边等； 3）从单元的节点选来分 顶点、边上点和形内点等。 （2）单元举例（二维） 其中最简单的单元是仅取顶点作节点的三节点直边三角形单元，这种单元又称为简单三角形单元。 （3）单元大小的确定 由有限元误差分析可知，在收敛的前提下，应力误差与单元的尺寸成正比，位移的误差与单元的尺寸平方成正比。单元划分愈大，位移的误差愈显著。 从有限元本身来看，单元划分得越细，节点布置得越多，计算结果越精确，但随之而来的是计算量越多，计算时间越长，计算费用越高，对计算机存贮量的要求也越大，而且输入输出的数据越多，数据准备和结果整理的工作量也就越大。 大量有限元计算结果表明，当单元网格划分细密程度超过了一定的限度，不但不能提高计算结果的精度，有时反而会使精度有所下降。因为过分加密网格，会使计算量激增从而导致累积计算误差增大。 合理地布置单元的大小——对于应力和位移状态需要详细了解的重要部位，以及如结构上的凹槽、裂纹等应力和位移变化比较剧烈的部位，单元应分得细一些，网格加密一些。 对于次要部位，或应力和位移变化平缓的部位，以及边界比较平滑的部位，单元可放大一些，网格可稀疏一些。但单元由小到大应逐渐过渡，相邻单元的尺寸不宜相差悬殊，否则会引起较大的计算误差。 两端受弯矩的板梁离散化例子 带凹槽结构的单元划分例子 （5）单元的过渡 在单元由大到小过渡时，必须特别注意任一单元的节点必须同时也是相邻单元的节点，而不能是相邻单元边上的内点。 （6）曲边的单元划分与单元的大小过渡 需要指出的是：矩形单元不仅不大适合边界曲折的结构，也不便在不同的部位采用大小不同的单元，通常采有的三角形单元和矩形单元混合的方式实现单元大小的过渡。 （7）结构厚度或弹性常数突变时的单元划分 当结构厚度或弹性常数发生突变时，在结构相应部位应力也有突变。 由后续单元分析的结论可知，由于通常假定三角形单元的厚度及材料常数在一个单元内均为常量，从而得出三角形单元中的应力也是常量，只有通过不同的单元反映应力的突变。 单元划分时应将厚度或弹性系数突变线作为单元的分界线。突变线不能穿过单元，否则不能正确反映这些实际存在的应力突变。 连杆小头 右图表示一连杆小头，铜质衬套和钢质连杆具有不同的弹性系数，小头和杆身的工字形断面又各有不同的厚度，故在突变处都划成单元的边界线。 （8）载荷突变部分部位的单元划分 如果结构受有集度突变的分布载荷或集中载荷，应把受力部分的单元划细些，并在载荷突变或集中之处布置节点，以使应力的突变得到一定程度的反映。 （10）特殊部位的单元划分 在工程设计中，人们常对某些特定断面上的应力或特殊点的位移表示关注，则在单元划分时应使这些特殊点的位移能立即算得，那些特定断面上的应力分