有限元基础理论教程lecture08.pptVIP

4.5 参数化建模 日志文件 用户操作过程，包括几何建模、网格划分、求解等，以命令流形式保存在日志文件中(*.log)。 将日志文件保存起来，在下一次使用ANSYS软件分析同一个模型时，可以先执行日志文件，得到前一次上机所完成的结果。 日志文件的缺点是保存了许多中间操作，不够简洁。 APDL语言与参数化建模 在建模过程中，比较容易修改，能逐步调试，适合建立复杂模型。 能够实现参数化建模，只要拓扑结构不变，就可以采用不同的参数值完成建模。 对ANSYS数据库进行一些直接操作、导出ANSYS模型中的数据。 将其它软件的计算结果导入ANSYS，利用ANSYS的后处理功能。 5 等参单元 5.1 四结点矩形单元 5.2 等参单元的基本概念 5.3四边形八结点等参单元 5.4等参单元的单元分析 5.5高斯积分 5.6六面体等参单元 5.1四结点矩形单元 对于包含n个结点的一维单元，采用Lagrange插值多项式直接构造位移插值函数，结点参数是待求的场函数的结点值，多项式定义如下。 为方便构造单元插值函数，引入单元自然坐标， 对于四结点矩形单元，引入自然坐标，用两个方向上的插值多项式的乘积，可以构造出矩形单元的插值函数。 5.2 等参单元的基本概念 5.3 四边形八节点等参单元 为了更好地反映物体的应力变化，适应曲线边界，在平面问题分析中经常使用四边形八节点等参单元。如图5-4所示。 结点2、3、4的形函数为： 对于结点5，可得到： 5.4等参单元的单元分析 以平面问题的四边形八结点等参单元为例，介绍构造等参单元的单元刚度矩阵的基本过程。 由于等参单元的形函数是局部坐标的函数，因此应变矩阵[B]也是局部坐标的函数。形成等参单元的单元刚度矩阵需要在整体坐标系中对局部坐标的函数进行积分，包括以下三个基本步骤： 1）计算用局部坐标表示的形函数对整体坐标x、y的偏导数； 2）将整体坐标系中的面积积分转换为在局部坐标系中的面积积分； 3）用数值积分计算出单元刚度矩阵中的元素。 1）计算形函数对整体坐标的偏导数 形函数Ni 是局部坐标的函数，同时也可以看作是整体坐标的函数，由复合函数求导法则可得： 2）将整体坐标系中的面积积分转换为在局部坐标系中的面积积分 在整体坐标系中，面积微元为x方向和y方向微矢量的叉乘的模量， 3）用数值积分计算出单元刚度矩阵中的元素 等参单元刚度矩阵的每个元素都是局部坐标的函数，不用解析的办法来计算局部坐标系中的积分，而采用数值积分方法。通常采用高斯积分方法计算单元刚度矩阵中的元素。 5.5 高斯积分 高斯积分方法预先定义了积分点和相应的加权系数，求出被积分的函数在指定积分点上的数值，加权后求和，就得到了该函数的积分。 二维、三维高斯积分可以用多重积分的方法，由一维高斯积分得到。二维积分，先令其中一个变量为常数，按照一维高斯积分公式进行内层积分 ， 进行三次分层积分，得到三维高斯积分公式， 用两个积分点进行高斯积分， 对于前面章节讲过的三角形单元，也能使用数值积分计算单元刚度矩阵的系数。 能否用高斯积分来计算？ 选择数值积分阶次的原则如下： 1）保证积分的精度，尽量满足精确积分的条件。 2）保证结构总刚度矩阵K是非奇异的。全部积分点能提供的独立关系的数目要大于系统的独立自由度数。 为了提高计算效率，可选取高斯积分的阶次低于被积函数精确积分所需要的阶次进行计算，这种积分方案被称为减缩积分。 四边形四结点、八结点的等参单元通常采用2×2的高斯积分。 将作用在单元上的外载荷同样表示为局部坐标的函数，就可以在局部坐标下完成单元的载荷移置。体力移置公式为， 在 的边上受到面力作用，面力移置公式为， 集中力移置公式为， 一维高斯积分的定义如下， 为积分点的位置 为对应的加权系数 表5-2 高斯积分的积分点坐标和加权系数 积分点数n 积分点坐标 加权系数 1 0.000000000000000 2.000000000000000 2 1.000000000000000 3 0.000000000000000 0.555555555555555 0.888888888888888 对应于(-1,1)积分域，高斯积分中所采用的积分点位置和对应的加权系数如下表所示， 再用一维高斯积分公式进行外层积分，得到二维高斯积分公式， 在任意积分域(a，b)上，用高斯积分公式计算时，先要把积分区域变换到(-1,1)区间上。 采用n个积分点的高斯积分可以达到2n-1阶的精度。如果被积分的函数是2n-1次多项式，用高斯积分可以得到精确的积分结果。计算两个简单多项式的积分，将高斯积分与精确积分结果做比较。 将积分点坐标和加权系数代入， 用三个积分点进行高斯积分， 在三角形单元