BJ1111--Chapter5一维单元及一维问题分析讲述.ppt

第五章 一维单元及一维问题分析 孙 会 本章内容 一、一维单元（重点，熟悉） 二、整体坐标、局部坐标和自然坐标（熟悉） 三、热传递问题（难点，了解） 四、固体力学问题 五、 ANSYS中一维单元举例（重点，熟悉） 六、ANSYS应用（重点，掌握） 七、结果验证（重点，熟悉） 引子 如图所示的等直截面散热片，散热片左端相连的基座温度已知，周围流体温度已知，试确定散热片沿其长度方向上的温度分布？ 引子 预处理阶段 1. 计算域离散化 2. 单元分析 2.1 假设描述单元行为的近似解 2.2 单元分析（直接法、最小总势能法、加权余数法） 2.3 组装单元，构造K(G) （按照标准方法进行） 2.4 施加边界条件（约束条件、载荷条件） 求解阶段 后处理阶段 一、一维单元 线性单元 预处理阶段： 计算域的离散化： 3个单元和4个节点 假设单元行为的近似解：采用分段线性函数近似。 1、线性单元 现以其中一个典型单元为研究对象 1、线性单元 1、线性单元 单元的温度分布 2、二次单元 如何希望提高计算精度怎么办？ 增加分析中使用的线性单元的数目 通过使用高阶的插值函数 采用二次函数代替原来的线性函数，此时要求使用3个节点来定义1个单元，第3个点可取在单元的中点，例如图4.3中的节点k。 2、二次单元 2、二次单元 3、三次单元 在二次单元中，采用二次函数代替线性单元中的线性函数，相应地提供了比线性单元更为精确的结果，那么是否可以采用更高阶的函数代替线性函数或二次函数呢？ 可以。如需更高的精度，可使用更高阶的插值函数，例如三次多项式。 用三次函数代替二次函数，要求至少使用4个节点来定义1个单元，单元被分成等长的3段，4个节点的取法如图5.4所示。 3、三次单元 3、三次单元 二、整体坐标、局部坐标和自然坐标 在有限元建模过程中，大多数情况下可使用多个参考系。 整体坐标系：(1)表示每个节点的位置和每个单元的方向；(2)施加边界条件和载荷；(3)表示计算求出的解。 局部坐标系和自然坐标系：当构造几何关系或计算积分时，可提供很多便利。 1、一维线性整体坐标与局部坐标 对于一维单元，整体坐标X和局部坐标x的关系为X=Xi+x，如图所示。 相应地，形函数可用局部坐标表示： 2、一维线性自然坐标 自然坐标是局部坐标的无量纲形式，用于单元刚度矩阵或传导矩阵推导中的积分计算。 使用自然坐标容易在上限1和下限-1间积分。 2、一维线性自然坐标 3、等参单元 可用形函数Si、Sj表示其他变量，如位移u： 也可通过形函数Si、Sj进行从整体坐标X或局部坐标x到自然坐标?的变换，即： 3、等参单元 4、用一维自然坐标表示二次和三次形函数 二次自然形函数： 三次自然形函数： 三、热传递问题 如图所示的等直截面散热片，散热片左端相连的基座温度已知，周围流体温度已知，试确定散热片沿其长度方向上的温度分布？ 三、热传递问题 预处理阶段： 1. 计算域离散化：3个单元、4个节点。 2. 单元分析 2.1 假设描述单元行为的近似解 采用分段线性函数描述 2.2 对单元进行分析 三、热传递问题 内容回顾--加权余数法的基本思想： 前提：已知描述问题的控制微分方程 为控制微分方程假设一个合理的近似解。该近似解必须满足给定问题的初始条件和边界条件。 将假设的近似解代入微分方程，方程右端不为0，将产生一定的余差，这一余差称为误差。 任何余差方法都要求误差在一定的间隔内或某些节点上消失(配置法、子区域法、迦辽金法、最小二乘法)。 三、热传递问题 已知直散热片的一维热传递控制微分方程： 涉及热传导、热对流 k --热传导率； A --散热片的横截面面积 h --对流热传递系数 p --散热片周长 Tf --周围流体温度。 三、热传递问题 除基座温度已知外，在散热片末端存在3种可能： 散热片末端温度等于周围流体温度，即： 散热片末端损失热量忽略不计，则： 需考虑散热片末端损失的热量，则： 三、热传递问题 一维问题的一般形式： 三、热传递问题 加权余数法中的迦辽金公式：误差对于权函数正交，且权函数是近似解的一部分。 典型单元的温度分布： 对于节点为 i 和 j 的任意单元，令迦辽金余差相对于权函数正交： 三、热传递问题 单元的热传导矩阵： 除最后一个单元，其它所有单元传导矩阵： 如果需考虑通过散热片末端的热量损失，最后一个单元的传导