有限元基础理论教程 lecture优秀培训书.pptVIP

5.4六面体等参单元 三维弹性力学问题的有限元法的基本步骤与平面问题的步骤一样，包括单元离散化、选择单元位移模式、单元分析、整体分析和方程求解。在分析三维问题时，所选择的单元主要为四面体单元和六面体单元。每个单元节点上定义有三个位移分量u、v、w。 三维问题有限元法有以下两个主要难点： 1）单元划分比较复杂 2）计算规模大 常用的三维等参单元有六面体八节点等参单元和六面体二十结点等参单元。等参单元的位移模式和坐标变化式采用相同的形函数， 用高斯积分计算出单元刚度矩阵的系数。用上节中类似的公式就可以在局部坐标下完成单元的载荷移置。 用单元的积分控制选项确定积分方式： Solid45单元采用2×2×2积分(默认方式) ，或采用单个积分点的减缩积分。 Solid95单元采用14点积分(默认方式)，或采用2×2×2的减缩积分。 当单元表面为四边形时，Solid45单元采用2×2积分，Solid95单元采用3×3积分。 当单元表面为三边形时， Solid45单元采用3点积分， Solid95单元采用6点积分。 6. 稳态热传导问题的有限元法 6.1热传导方程与换热边界 6.2加权余量法 6.3稳态温度场分析的一般有限元列式 6.4三角形单元的有限元列式 6.5温度场分析举例 6.1热传导方程与换热边界 物体内部的温度分布取决于物体内部的热量交换，以及物体与外部介质之间的热量交换，一般认为是与时间相关的。传热过程是基本的物理过程之一，涉及许多工程应用领域。 法国数学家Jean Fourier在1822年提出了导热计算的基本规律：在物体内发生纯导热时，单位时间内通过垂直于热流方向的单位面积的热流量，与等温面法线方向上的温度增量成正比，而与法向距离成反比。 6.2 加权余量法 * ANSYS提供的Solid45单元就是六面体八节点等参单元，每个节点有代表x、y、z三个方向位移的三个自由度（DOF，Degree of Freedom） 六面体八节点等参单元的形函数为， 为结点的局部坐标 ANSYS提供的Solid95单元是六面体二十节点等参单元，每个节点有代表x、y、z三个方向位移的三个自由度。 利用坐标转换，将整体坐标系下的积分转换为在局部坐标系下的积分。 在整体坐标系中的体积微元为 ， 微矢量在局部坐标系中表示为， 三维单元的雅可比矩阵为， 单元刚度矩阵在局部坐标系中的积分式为， 体力移置， 面力移置， 集中力的移置， 对于各向同性材料，热传导方程可写为以下形式， 除了热传导方程，计算物体内部的温度分布，还需要指定初始条件和边界条件。初始条件是指物体最初的温度分布情况， 物体内部的热交换采用以下的热传导方程（Fourier方程）来描述， 边界条件是指物体外表面与周围环境的热交换情况。在传热学中一般把换热边界条件分为三类。 1)给定物体边界上的温度，称为第一类边界条件。 2)给定物体边界上的热量输入或输出，称为第二类边界条件。 3）给定对流换热条件，称为第三类边界条件。 物体与其相接触的流体介质之间的对流换热系数和介质的温度为已知。 h为换热系数，W/(m2 K) Tf为流体介质的温度， Ts为物体表面的温度。 如果边界上的换热条件不随时间变化，物体内部的热源也不随时间变化，在经过一定时间的热交换后，物体内各点温度也将不随时间变化，这类问题称为稳态（Steady state）热传导问题。 三维问题的稳态热传导方程为， 对于各向同性的材料，可以得到Poisson方程， 不包含内热源，各向同性材料中的温度场满足Laplace方程 在分析稳态热传导问题时，不需要考虑物体的初始温度分布对最后的稳定温度场的影响，因此不必考虑温度场的初始条件，而只需考虑换热边界条件。计算稳态温度场实际上是求解偏微分方程的边值问题。 由于实际工程问题中的换热边界条件比较复杂，在许多场合下也很难进行测量，如何定义正确的换热边界条件是温度场计算的一个难点。 微分方程的边值问题，可以一般地表示为未知函数u满足微分方程组， 满足边界条件， （6-10） (6-11) 二维稳态热传导问题的控制方程与边界条件如下 ： 上述微分方程只有在区域比较规则、边界条件比较简单时才得到精确解。 对于复杂的工程问题，微分方程组的精确解往往很难找到，需要设法寻找近似解。所选取的近似解是一族带有待定参数的已知函数，一般表示为 为待定系数， 为已知函数，被称为试探函数。 试探函数要取自完全的函数序列，是线性独立的。 *