04结构力学平面有限元2.pptVIP

哈工大 土木工程学院 * / 75 ? 单元位移函数 ? 单元协调性 位移函数中包含常数项和一次项，故满足完备性条件；边界上线性变化且由结点确定，故满足位移连续。所以矩形双线性单元是收敛的。 ? 哈工大 土木工程学院 * / 75 5-3-2 应变和应力矩阵 应变矩阵 ? 应变矩阵 ? 哈工大 土木工程学院 * / 75 ? 应力矩阵 应力矩阵 对于平面应力问题 ? 哈工大 土木工程学院 * / 75 对于平面应变问题 做如下替换 从应变矩阵和应力矩阵可见：单元内应力和应变沿 x 方向是线性变化的，沿 y 向也是线性变化的。 这是由所设位移模式所决定的。 ? 哈工大 土木工程学院 * / 75 5-3-3 单元刚度矩阵和单元等效荷载列阵 利用势能原理推导单元刚度矩阵和等效结点荷载列阵 ? 单元刚度矩阵 ? 哈工大 土木工程学院 * / 75 平面应力问题 ? 哈工大 土木工程学院 * / 75 平面应变问题 做如下替换 ? 单元等效结点荷载 ? End 哈工大 土木工程学院 * / 75 5-3-4 计算结果整理 整理方法与三角单元结果整理方法相同 ? 哈工大 土木工程学院 * / 75 5-3-5 计算实例 例题1 Page 110 ? 哈工大 土木工程学院 * / 75 §5-4 计算程序简要说明 End 哈工大 土木工程学院 * / 75 如在 ij 边 3. 在三角形单元内任一边上的形函数只与该边的两端点坐标有关且是线性的，而与其它结点坐标无关。 2 1 i、j边的直线方程为 在i、j边上形函数 根据性质1 ? 哈工大 土木工程学院 * / 75 注意应用 则有 2 1 ? 哈工大 土木工程学院 * / 75 在 ij 边的位移 由此性质可知：单元边界是线性变化的。相邻单元在公共结点有相同位移，所以能保证相邻单元的位移协调。 ? 哈工大 土木工程学院 * / 75 5-2-2 用面积座标建立单元位移场 为以后方便理解和确定复杂单元内任意一点处形函数的值，在此引入面积坐标的概念 在单元内任一点 P(x , y) 可用三个数（Li，Lj，Lk）来确定 ? 哈工大 土木工程学院 * / 75 因为 所以 可见面积坐标就是形函数 ? 哈工大 土木工程学院 * / 75 由面积坐标的定义不难发现，平行于某边直线（ik）上的所有点具有相同的坐标（Lj），并且该坐标就等于“该直线到至（ik）边的距离”与“点 j 到（ik）边的距离”之比。 ? 哈工大 土木工程学院 * / 75 三个结点的面积坐标分别为： 不难验证面积坐标与直角坐标的变换关系 ? 哈工大 土木工程学院 * / 75 5-2-3 基于势能原理的单元分析 微分算子矩阵 ? 哈工大 土木工程学院 * / 75 几何矩阵只与结点坐标有关，它的元素都是常数，所以单元内每点的应变也必为常数，所以三结点三角元称常应变三角元. ? 哈工大 土木工程学院 * / 75 应变势能 外力势能 体系势能 ? 哈工大 土木工程学院 * / 75 由最小势能原理 ? 哈工大 土木工程学院 * / 75 荷载列阵 （1）单元上任意点的集中荷载 直接加在相应位移方向上 非结点荷载等效到结点上 所以 ? 哈工大 土木工程学院 * / 75 （2）单元上体积力 所以 ? 哈工大 土木工程学院 * / 75 结果分析 单元应力 ? ?x ? xy ? yx ?y ? n ? 哈工大 土木工程学院 * / 75 应力圆半径 平均应力 最大主应力 最小主应力 主方向角 ? 2? ?1 ?2 ?y ?m ?x ? ? O ?xy ? 哈工大 土木工程学院 * / 75 5-2-4 几点结论 一.结点的选择和单元划分 1.集中力作用点、分布力突变点、支承点应选作结点。 2.不同厚度、不同材料的部分不应划在同一个单元。 3.应力变化大处单元应密集一些。结点的多少与疏密要考虑计算机的容量和计算精度。 有限元分析应注意的问题和结果整理 ? 哈工大 土木工程学院 * / 75 4.单元边界的边长之比应尽可能靠近1。 不宜 5.相邻单元的尺寸尽可能接近。 6.结点所连接的单元个数尽可能一致。 适宜 适宜 不宜 ? 哈工大 土木工程学院 * / 75 二.结点编码 尽可能使相关结点的结点编码差值最小，以缩小刚度矩阵的带宽。从而提高储存和运算效率 总刚半带宽=(相关结点最大差值+1)*结点位移数 ? 哈工大 土木工程学院 * / 75 总刚半带宽=(相关结点最大差值+1)*结点位移数 总刚半带宽=16总刚需占用的存贮空间为: 16?14?2=448 总刚半带宽=6总刚需占用的存贮空间为: 6?14?2=168 1 2 3