第四章平面问题有限元及程序设计教程.ppt

最近的研究课题 系泊浮体在波群作用下的运动响应 锦州大笔架山天桥修复方案研究 营口白沙湾岸滩冲刷防治措施研究 狭长型海湾污染综合治理措施研究 海域管理辅助决策系统研究 新型码头结构形式研究 桥墩冲刷防治措施研究 水质预报反演研究 港工建筑物对水动力环境影响仿真研究 吸力锚研究 眼睛流体动力研究 旅顺鱼人码头波浪、潮流、泥沙数值模拟研究 福佳大化护岸损坏修复波浪研究 若干港口泊稳条件物理模型试验和数值模型试验研究 东海狭长形海湾综合整治研究 典型海湾破坏机理研究 海事搜救系统研究 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 讲一讲实际结构的约束，以及在有限元方法中的简化。 * 补充杆系内力特征、与弹性力学平面问题的差别，讲清楚场值与点值的关系。 连续弹性体被离散，但单元个体还是一个连续、均匀、各向同性的完全弹性体。 注意规律性。 * 讲清楚这些条件背后的为什么。 * 讲清楚为什么位移的精度和应力、应变的精度 以后再讲高次位移函数。 * 以后再介绍完备、协调的概念。 * 1、虚功原理：应力在虚应变上所做的功等于外力在虚位移上所做的功。 2、这仅是推导单元刚度矩阵的一种方法，以后再讲另一种方法-变分原理。 * 保证b、c相等，即可保证单元刚度矩阵相同。错！！不一定要相等，只要两个矩阵B相乘相等即可。要保证这点，单元编号规矩要一致。 * 注意荷载的正负号 * 先介绍整体刚度方程（由结点平衡得到），然后再介绍刚度矩阵。 * 结合此例，详细介绍单元刚度矩阵每个元素应累加到整体刚度矩阵的位置。两种对应编码方式：结点编码、自由度编码。 结点荷载的累加。 * 每个元素的物理含义再提一下。 * * * * §4.3 平面问题的单元分析 单元性质分析 位移 是x、y的线性函数； 误差是△ x、 △ y的二阶小量； 应变 应力 常量； 相邻单元连续 相邻单元不连续 误差是△ x、 △ y的一阶小量； 提高精度方法： 1）减小单元尺寸； 2）采用高次位移函数，提高位移、应变和应力的精度； 收敛速度和精度估计 若单元的插值函数是完备而协调的，当单元尺寸不断缩小而趋于零时，有限元解将趋于真正解。 在有些情况下，如果用于有限元场函数近似解的多项式能精确地拟合真正解，则在有限数目的单元划分（甚至仅仅是一个单元）的条件下，也能得到精确的解答。例如真正解是二次函数，若有限元的插值函数也包含了二次的完全多项式，则有限元解就能得到精确的解答。由此我们可以得到精度与单元尺寸的关系。例如位移可以展开成Taylor级数： 这只是形式上的精度估计，并不能对有限元解的误差做出具体的估计。而后者在 实际分析工作中更有用。一般可以通过两种途径解决： 单元结点力 §4.3 平面问题的单元分析 §4.3.3 单元的 刚度矩阵 单元结点位移 单元应力向量： 给定一个虚位移： 单元虚应变： i (xi, yi) x y j (xj, yj) m (xm, ym) 虚功原理:内力虚功等于外力虚功 §4.3 平面问题的单元分析 §4.3.3 单元的 刚度矩阵 t 为单元厚度 由于虚位移是任意给定的可能位移，故： 其中， 是单元的刚度矩阵 (6 × 6)； §4.3 平面问题的单元分析 单元的刚度矩阵为： 写成分块矩阵： 其中： 平面应力： 平面应变： 写成元素矩阵： §4.3 平面问题的单元分析 单元刚度矩阵的特点： 1) 对称性： 2) 与单元尺寸无关，放大或缩小尺寸，单元刚度矩阵不变； 3）奇异性：它不存在逆阵 4）主元（对角线元素）恒正 §4.3 平面问题的单元分析 单元①的刚度矩阵为： 单元②的刚度矩阵为： 作用在单元上的荷载，既有结点荷载，也有非结点荷载，因此需要将非结点荷载转换成等效的结点荷载。 §4.3.4 等效结点荷载 等效结点荷载和原荷载在任何虚位移产生的虚功相等； 刚体： 等效原则： §4.3 平面问题的单元分析 结点荷载： 非结点荷载： 直接集成到荷载列向量； 等效成结点荷载； 原荷载与等效结点荷载在任一轴上的投影之和相等，对任一轴的力矩之和也相等。 等效结点荷载向量： 静力等效 §4.3 平面问题的单元分析 1）体积力的等效结点荷载 单元结点为 i, j, m，密度为 ρ ，任意一点的体积力向量为： 假设单元各结点发生虚位移： 则单元内任意一点的虚位移为： 根据虚功原理： §4.3 平