平面问题基本理论及有限元求解.pptVIP

平面问题定义（总结） 严格地讲，任何结构都是空间的 对于某些特殊情况，空间问题可以转化为平面问题。 1、平面应力问题 满足条件： 1）几何条件 厚度尺寸远远小于截面尺寸 2）载荷条件 载荷平行于板平面且沿厚度方向均匀分 布，而板平面不受任何外力作用 此时， 应力应变分量变为： 几何方程 物理方程 2、平面应变问题 满足条件 1）几何条件 结构呈等截面的细长形 2）载荷条件 载荷垂直于厚度方向(平行横截面)且 沿厚度均匀分布，两个端面不受力 此时， 应力应变分量变为： 物理方程 弹性力学的参量及方程汇总 1）参量 位移： 应力： 应变： 2）方程 平衡方程： 几何方程： 物理方程： 其中： 第二节 平面问题的有限元法 平面问题的有限元分析步骤（平面应力问题） 1、结构离散 离散：将一个连续的弹性体分割为一定形状和数量 的单元的组合 单元也称为网格 连续体 → 有限个单元的组合体 可用于离散的单元 三角形单元 矩形单元 不规则四边形单元 DOF 节点的自由度：节点所具有的位移分量的数量。 一个单元所有节点的自由度总和称为单元自由度 单元参数只能通过节点传递到相邻单元 单元和节点必须统一编号 2、单元分析（位移、应力、应变） 任务：形成单元刚度矩阵，建立单元特性方程 因此必须建立坐标系，如下图： 1）位移函数 分片插值→ 假设一种函数来表示单元位移分布 一般选取多项式（简单而且易求导） 对于三角形单节点单元（DOF=6） 2)单元应力和应变 将位移表达式（2-3）和（2-4）代入几何方程得： 则应力向量可以表示为： [D]、[B]均为常数矩阵，因此三角形三节点单元为常应力单元。 例4-1 图示为一平面应力单元，厚度为t，求形函数矩阵、单元应变矩阵、单元应力矩阵。 3）单元刚度矩阵 设作用在单元节点上的单元节点力列阵为： 而节点发生的虚位移为： 则节点力在虚位移上做的虚功为： 将（2-5）展开得： 单元刚度矩阵的物理意义： 在一个节点处产生单位位移而其他点为零时，在该节点上需要的外力的大小。 单元刚度矩阵的特性 （1）对称性： （弹性力学互等定理） （2）奇异性： （刚体位移） 单元刚度矩阵的特性 （3）单元刚度矩阵主对角元素恒为正值； （4）单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和为零，由此可知，单元刚度矩阵各行和各列元素的总和为零。 例4-1 图示为一平面应力单元，厚度为t，求形函数矩阵、单元应变矩阵、单元应力矩阵和单元刚度矩阵。 3、总刚度矩阵的集成 通过单元特性方程 并不能求出单元节点位移 。因为 包含单元间的作用力。 因此，必须将每个单元的特性方程相加消除内力的影响。这就是总刚度矩阵集成的目的。 1）总刚集成原理 在整个结构中，一个节点为几个单元共有 在第i个节点处的平衡方程为： 例：总刚的形成过程 2）总刚度矩阵集成过程 （1）扩阶过程（可由转换矩阵完成） （2）叠加过程： 总刚矩阵的特点： （1）对称性 节省存储容量 （2）稀疏性 可能存在大量零元素 （3）带状性 半带宽与节点的编号有关 （4）奇异性 保证刚体位移 （5）主对角元素恒为正值 4. 载荷移置移置原则 （1）刚体静力等效原则——使原荷载与移置荷载的主矢量以及对同一点的主矩也相同。 4、载荷移置 移置可能在局部产生误差，但不会影响整个结构的力学特性。 1）集中力的移置（虚功等效） 2）面力的移置 3）体力的移置 5、约束处理 1）边界位移为零 2）边界位移为已知量 6、求解线性方程组 7、计算其它物理量 8、计算结果处理 9、结果显示、打印、分析 收敛准则分析 位移函数应该满足以下几个条件 （1）包括常数项（保证刚体位移） （2）包括一次项（保证常应变） （3）保证位移的连续性（形函数性质3保证） （4）各项几何同性（x, y应该是可以互换的） 满足上述三个条件的目的是满足有限元的收敛性 （1）和（2）是收敛的必要条件 完备性条件 （3）是收敛的充分条件