有限元及工程软件第十章.ppt

* * 第十章 热应力问题的有限元法 第一节 平面稳定温度场 结构受热 温度改变 变形与应力 有限元分析：首先分析结构的温度场，然后计算热 应力。 以平面问题为例介绍。 一、平面稳定温度场的概念及分类 平面结构受热与放热互相平衡时，结构温度不 再随时间而改变，平面内的温度只是坐标的函数， 称为平面稳定温度场。 在平面求解区域内，有 热传导系数 热源强度 一般的工程结构，本身不产生热量，热量多是由外 界传入，有 对应不同的热边界条件，微分方程的解是不同的， 对于平面结构有不同的温度分布。 第一类边界条件 平面结构的边界上保持给定的分布温度，即 边界的分布温度 第三类边界条件 在边界处与周围介质存在热交换，即 进入的热流 周围介质的温度 放热系数 上式中包含边界温度和温度梯度，称为第三类边 界条件，是一种混合边界。 第二类边界条件 确定结构边界处的温度梯度，称为第二类边界条件。 当 ， 取不同值时，上式可以转化为不同的 边界条件，因此可以统一地编制有限元程序。同时， 对于平面结构的不同边界部位，通过改变 ， ， 即可。 绝热条件：在边界处和周围介质没有热交换。 二、第一类边界条件问题 等价于 有限元求解此 变分问题 简单三角形单元，单元内温度假定 为线性分布 设单元3个顶点的温度分别为 单元节点温度列阵为 单元内各点温度为 形状函数矩阵 简单三角形单元，温度场为线性分布，则形状函数为 三角形单元满足完备性和协调性，可以收敛于真实解。 对任意单元e，单元泛函为 其中： --应变矩阵 其中： --单元“刚阵” 结构的泛函为 其中： --整个结构的“刚阵” 结构总的泛函是节点温度的二次齐次式。 即 结合边界条件，求解方程组。 三、第三类边界条件问题 等价于 边界上的弧坐标 采用简单三角形单元，单元泛函表示为 其中： --e单元所拥有的边界。 只有靠近边界的单元才具有 这一项。 设单元e的m，l节点位于边界上，计算时应以直 线lm代替部分边界。 e单元的温度为 其中： --边界“刚度”矩阵 --边界“载荷” 结构总泛函为 其中： --结构总“刚度”矩阵 --结构总“载荷” 即 已知“载荷”，求解方程组。 对于三类边界条件，按上述分析方法处理，而在 上述分析时没有考虑一类边界条件，可在求解方程组 时考虑，即：使该边界处的节点温度取为给定值。 求解方程组时，边界条件的处理： 将整个边界按三类边界处理，而对于一类边界位 置，介质温度 取为给定值，并将放热系数 取为 相当大的值。 第二节 平面热应力 设等截面杆件，原长为l，温度由 变为 ，则 热应变为 各向同性三维单元，由温升引起的热应变为 当物体各部分有同样温升时，热膨胀是均匀的， 若物体受外界约束，则处于各方向应变都相同的常 应变状态，不会产生内部应力。 当物体受热，又受外界约束时，或内部受热不均 匀，则内部会产生内应力。总应变应为热应变和弹性 应变之和。 结构由于温度变化引起的内部应力 ，称为热应力。 用位移法分析结构热应力，应按温度场的改变 计算热变形，进而计算热应力。 一、平面热变形的计算 平面问题应变和位移的关系为 应变为 弹性应变 热应变 设为各向同性问题，有 设为弹性问题，有 平面问题弹性系数矩阵 取 注意：由于位移是两个分量，温升为一个分量， 插值的形函数是一样的，但形函数矩阵是不同的。 单元e的弹性应变能为 相等 总应变和节点位移的关系 应变矩阵 对于三角形单元，应变矩阵与一般平面问题的三 角形单元一致。