《有限元分析基础》PPT课件

1 有限元分析基础 2 ANSYS应用 2 内容结构 第一章概述 第六章空间问题的有限单元法 第七章轴对称旋转单元 第五章等参元 第四章平面结构问题的有限单元法 第三章杆系结构静力分析的有限单元法 第二章结构几何构造分析 3 1 1有限单元法的概念 1 2有限单元法基本步骤 1 3工程实例 第一章概述 4 1 1有限单元法的概念 基本思想 借助于数学和力学知识 利用计算机技术而解决工程技术问题 FiniteElementMethod FEMFiniteElementAnalysis 第一章概述 5 第一章概述 三大类型 按其推导方法分 1 直接刚度法 简称直接法 根据单元的物理意义 建立有关场变量表示的单元性质方程 2 变分法直接从求解泛函的极值问题入手 把泛函的极植问题规划成线性代数方程组 然后求其近似解的一种计算方法 3 加权余量法直接从控制方程中得到有限单元方程 是一种近似解法 6 1 2有限单元法基本步骤 1 待求解域离散化 2 选择插值函数 3 形成单元性质的矩阵方程 4 形成整体系统的矩阵方程 5 约束处理 求解系统方程 6 其它参数计算 第一章概述 7 图1 2工程问题有限单元法分析流程 第一章概述 8 1 3工程实例 a 铲运机举升工况测试 b 铲运机工作装置插入工况有限元分析 图1 3WJD 1 5型电动铲运机 第一章概述 9 a KOMATSU液压挖掘机 b 某液压挖掘机动臂限元分析图1 4液压挖掘机 第一章概述 10 图1 5驾驶室受侧向力应力云图图1 6接触问题结构件应力云图 第一章概述 11 图1 7液压管路速度场分布云图图1 8磨片热应力云图 图1 9支架自由振动云图 第一章概述 12 第二章结构几何构造分析 2 1结构几何构造的必要性 2 2结构计算基本知识 2 3结构几何构造分析的自由度与约束 13 2 1结构几何构造的必要性 结构是用来承受和传递载荷的 如果不计材料的应变 在其受到任意载荷作用时其形状和位置没有发生刚体位移时 称之为几何不变结构或几何稳定结构 反之则称为几何可变结构或几何不稳定结构 几何可变结构不能承受和传递载荷 对结构进行几何构造分析也是能够对工程结构作有限单元法分析的必要条件 第二章结构几何构造分析 14 a 结构本身可变 b 缺少必要的约束条件 c 约束汇交于一点图2 1几何可变结构 第二章结构几何构造分析 15 2 2结构计算基本知识 2 2 1结构计算简图 实际结构总是很复杂的 完全按照结构的实际情况进行力学分析是不可能的 也是不必要的 因此在对实际结构进行力学计算之前 必须将其作合理的简化 使之成为既反映实际结构的受力状态与特点 又便于计算的几何图形 这种被抽象化了的简单的理想图形称之为结构的计算简图 有时也称为结构的力学模型 结构计算所常用的结点和支座的简化形式 1 结点 铰结点 刚结点 混合结点 2 支座 活动铰支座 固定铰支座 固定支座 定向支座 第二章结构几何构造分析 16 2 2 2结构的分类与基本特征 按结构在空间的位置分结构可分为平面结构和空间结构两大类 2 按结构元件的几何特征分 杆系结构 梁 拱 桁架 刚架 桁构结构等 板壳结构 实体结构实体结构的长 宽 高三个尺寸都很大 具有同一量级 混合结构 第二章结构几何构造分析 17 2 3结构几何构造分析的自由度与约束 1 自由度指结构在所在空间运动时 可以独立改变的几何参数的数目 也就是确定该结构位置时所需的独立参数的数目 2 约束指减少结构自由度的装置 即限制结构结构运动的装置 a 支座链杆的约束b 铰的约束 单铰 复铰 完全铰与不完全铰 第二章结构几何构造分析 18 3 1结构离散与向量表示 第三章杆系结构静力分析的有限单元法 3 2位移函数及单元的刚度矩阵 3 3坐标变换及单元刚度矩阵 3 4整体刚度矩阵 3 5约束处理及求解 3 6计算示例 3 7ANSYS桁架结构计算示例 3 8ANSYS刚架结构计算示例 19 3 1结构离散与向量表示 工程上许多由金属构件所组成的结构 如塔式桁构支承架 起重机起重臂架 钢结构桥梁 钢结构建筑等可以归结为杆系结构 杆系结构按各杆轴线及外力作用线在空间的位置分为平面杆系和空间杆系结构 杆系结构可以由杆单元 梁单元组成 a Liebherr塔式起重机 b Liebherr履带式起重机 c 钢结构桥梁 d 埃菲尔铁塔图3 1杆系结构 第三章杆系结构静力分析的有限单元法 20 3 1 1结构离散化 由于杆系结构本身是由真实杆件联接而成 故离散化比较简单 一般将杆件或者杆件的一段 一根杆又分为几个单元 作为一个单元 杆件与杆件相连接的交点称为结点 杆系结构的离散化的要点可参考如下 a 杆件的转折点 汇交点 自由端 集中载荷作用点 支承点以及沿杆长截面突变处等均可设置成结点 这些结点都是根据结构本身特点来确定的 b 结构中两个结点间的每一个等截面直杆可以设置为一个单元 变换为作用在结点上的等效结点载荷 第三章杆系结构静力分析的有限单元法 21 c 变截面杆件可分段处理成多个单元 取各段中点处的截面近似作为该单元的截面 各单元仍按等截面杆进行计算 d 对曲杆组成的结构 可用多段折线代替 每端折线为一个单元 如若提高计算精度 也可以在杆件中间增加结点 e 在有限元法计算中 载荷作用到结点上 当结构有非结点载荷作用时 应该按照静力等效的原则将其等效结点荷载 第三章杆系结构静力分析的有限单元法 a 结点载荷处理方式 b 等效结点载荷处理方式图3 2杆系结构离散化示意图 22 3 1 2坐标系 图3 3坐标系示意图 为了建立结构的平衡条件 对结构进行整体分析 尚需要建立一个对每个单元都适用的统一坐标系 即结构坐标系或称之为整体坐标系 总体坐标系 第三章杆系结构静力分析的有限单元法 23 3 1 3向量表示 在有限单元法中力学向量的规定为 当线位移及相应力与坐标轴方向一致时为正 反之为负 转角位移和力矩 按右手法则定出的矢量方向若与坐标轴正向相一致时为正 对于任意方向的力学向量 应分解为沿坐标轴方向的分量 刚架结构示意图 b 结点位移和结点力分向量图3 4平面刚架分析示意图 第三章杆系结构静力分析的有限单元法 24 结点位移列向量为 单元e结点位移列向量为 结点力向量为 单元e结点力列向量为 第三章杆系结构静力分析的有限单元法 25 3 2位移函数及单元的刚度矩阵 3 2 1轴向拉压杆单元的位移的函数 有限单元法分析中 虽然对不同结构可能会采取不同的单元类型 采用的单元的位移模式不同 但是构建的位移函数的数学模型的性能 能否真实反映真实结构的位移分布规律等 直接影响计算结果的真实性 计算精度及解的收敛性 为了保证解的收敛性 选用的位移函数应当满足下列要求 a 单元位移函数的项数 至少应等于单元的自由度数 它的阶数至少包含常数项和一次项 至于高次项要选取多少项 则应视单元的类型而定 第三章杆系结构静力分析的有限单元法 26 由单元结点位移 确定待定系数项当时 当时 所以用结点位移表示其中 分别表示当 时 时的单元内的轴向位移状态 故称为轴向位移形函数 第三章杆系结构静力分析的有限单元法 b 单元的刚体位移状态和应变状态应当全部包含在位移函数中 c 单元的位移函数应保证在单元内连续 以及相邻单元之间的位移协调性 27 3 2 2梁单元平面弯曲的位移函数梁单元平面弯曲仅考虑结点的四个位移分量 由材料力学知 各截面的转角 故梁单元平面弯曲的位移表达式可分为仅包含四个待定系数 的多项式单元结点位移条件当时 当时 第三章杆系结构静力分析的有限单元法 28 称为形函数矩阵 第三章杆系结构静力分析的有限单元法 29 3 2 3单元的应力应变在弹性范围内 并且不考虑剪力的影响时 平面刚架单元内任一点的轴向线应变由两部分组成 即轴向应变与弯曲应变之和 其轴向应变与平面桁架轴向应变相同 轴向应变为弯曲应变为y为梁单元任意截面上任意点至中性轴 x轴 的距离 得出平面刚架单元应变 图3 5弯曲应变计算示意图 则 平面刚架梁单元的应变转换矩阵 第三章杆系结构静力分析的有限单元法 30 3 2 4平面刚架梁单元的刚度矩阵梁单元的i j结点发生虚位移为 单元内相应的虚应变应为 由虚功原理有 由于结点虚位移的任意性 故上式可写成 上式称为局部坐标下的平面刚架单元的刚度方程 简称为单刚 第三章杆系结构静力分析的有限单元法 31 横截面积A横截面对形心轴z的静矩S横截面对主惯性轴z的惯性矩I得到四个33子块所组成的局部坐标系下的平面刚架梁单元的单元刚度矩阵 第三章杆系结构静力分析的有限单元法 32 平面桁架的单元刚度矩阵为 空间桁架单元每个结点有3个位移分量 其单元结点位移列向量 空间桁架局部坐标下的单元刚度矩阵是6 6的 第三章杆系结构静力分析的有限单元法 33 空间刚架单元每个结点有6个位移分量 其单元结点位移列向量 空间刚架局部坐标下的单元刚度矩阵是12 12的 a 杆单元i端产生单位位移 b 杆单元j端产生单位位移图3 6平面桁架单元刚度系数的物理意义 a 梁单元i端产生单位位移 b 梁单元j端产生单位位移 第三章杆系结构静力分析的有限单元法 34 c 梁单元i端产生单位角位移 d 梁单元j端产生单位角位移图3 7平面刚架单元刚度系数的物理意义 3 2 5单元的刚度矩阵的性质a 单元刚度矩阵仅与单元的几何特征和材料性质有关 仅与单元的横截面积A 惯性矩I 单元长度l 单元的弹性模量E有关 b 单元刚度矩阵是一个对称阵 在单元刚度矩阵对角线两侧对称位置上的两个元素数值相等 即 根据是反力互等定理 c 单元刚度矩阵是一个奇异阵 d 单元刚度矩阵可以分块矩阵的形式表示 具有确定的物理意义 第三章杆系结构静力分析的有限单元法 35 3 3坐标变换及单元刚度矩阵 3 3 1坐标变换在整体坐标系中单元结点力向量和结点位移列向量可分别表示成 a 向量转换分析 b 向量转换图3 8向量转换示意图 第三章杆系结构静力分析的有限单元法 36 对于梁单元如图3 8 b 所示 则有 可简写为 第三章杆系结构静力分析的有限单元法 37 同理 式中 平面刚架梁单元的从局部坐标系向整体坐标系的转换矩阵 3 3 2整体坐标系下的单元刚度矩阵 式中 整体坐标下的单元刚度矩阵 和一样 为对称阵 奇异阵 第三章杆系结构静力分析的有限单元法 38 3 4整体刚度矩阵 3 4 1整体刚度矩阵的建立整体刚度矩阵也称之为结构刚度矩阵或总体刚度矩阵 简称总刚 整体刚度矩阵的求解是建立在结构平衡条件的基础之上 因此研究对象以整体坐标系为依据 图3 9载荷向量示意图 如右图所示刚架结构 其结点载荷列向量分别为 第三章杆系结构静力分析的有限单元法 39 结构载荷列向量 结点位移列向量 建立结点平衡条件方程式如右表 第三章杆系结构静力分析的有限单元法 40 用分块矩阵的形式 建立杆端内力与结点位移的关系式 第三章杆系结构静力分析的有限单元法 41 第三章杆系结构静力分析的有限单元法 42 第三章杆系结构静力分析的有限单元法 43 单元刚度矩阵由2 2的子矩阵组成 每个子矩阵是3 3的方阵 的上角标表示单元编号 下角标表示单元j端单位位移所引起的i端相应力 将杆端内力与结点位移关系式代入结点的平衡条件方程式中 经整理得 简写为 称之为结构原始平衡方程 其中 为整体刚度矩阵 第三章杆系结构静力分析的有限单元法 44 3 4 2整体刚度矩阵的集成整体刚度矩阵是由在整体坐标系下 矩阵按照结点编号的顺序组成的行和列的原则 将全部单元刚度矩阵扩展成n n方阵后对号入座叠加得到 对于单元1 对于单元2 对于单元3 单元刚度矩阵集成得出整体刚度矩阵 第三章杆系结构静力分析的有限单元法 45 3 4 3整体刚度矩阵的性质整体刚度矩阵中位于主对角线上的子块 称为主子块 其余为副子块 a 中主子块由结点i的各相关单元的主子块扩展之后叠加求得 即b 当结点i j为单元e的相关结点时 中副子块为该单元e相应的副子块 即 c 当结点i j为非相关结点时 中副子块为零子块 即 d 仅与各单元的几何特性 材料特性 即A I l E等因素有关 e 为对称方阵 f 为奇异矩阵 其逆矩阵不存在 因为建立整体刚度矩阵时没有考虑结构的边界约束条件 第三章杆系结构静力分析的有限单元法 46 g 为稀疏矩阵 整体刚度矩阵中的非零元素分布区域的宽度与结点编号有关 非零元素分布在以对角线为中心的带状区域内 称为带状分布规律 见图3 10 a 在包括对角线元素在内的区域中 每行所具有的元素个数叫做把半带宽 以d表示 最大半带宽等于相邻结点号的最大差值加1与结点自由度数的乘积 结点号差越大半带宽也就越大 计算机以半带宽方式存储 见图3 10 b 半带宽越窄 计算机的存储量就越少 而且可以大幅度减少求解方程所需的运算次数 其效果对大型结构显得尤为突出 图3 10整体刚度矩阵存储方法h 整体刚度矩阵稀疏阵 故整体刚度矩阵不能求逆 必须作约束处理方能正确地将结点位移求出 进而求出结构的应力场 a 带状分布规律 b 带状存储 第三章杆系结构静力分析的有限单元法 47 3 5约束处理及求解 3 5 1约束处理的必要性建立结构原始平衡方程式时 并未考虑支承条件 约束 也就是说 将原始结构处理成一个自由悬空的 存在刚体位移的几何可变结构 整体刚度矩阵是奇异矩阵 因此 无法求解 可以参照第2章的原则 结合实际工程结构引入支承条件 即对结构原始平衡方程式做约束处理 约束处理后的方程称为基本平衡方程 统一记为 3 5 2约束处理方法约束处理常用方法有填0置1法和乘大数法 采用这两种方法不会破坏整体刚度矩阵的对称性 稀疏性及带状分布等特性 第三章杆系结构静力分析的有限单元法 48 下面以图3 11所示刚架结构为例 解释如何进行约束处理 对于下图所示刚架结构 设结点位移列向量为设结点载荷列向量为 固定支座 b 支座强迫位移已知图3 11结构约束 第三章杆系结构静力分析的有限单元法 49 其原始平衡方程式为 按照每个结点的位移分量将上式展开为 第三章杆系结构静力分析的有限单元法 50 对于如图3 11 a 所示 结构约束 支座 位移全部为零 此时做约束处理时 采用填0置1法比较适宜 对于如图3 11 b 所示 某约束 支座 位移为给定的强迫值 此时做约束处理时 采用乘大数法比较适宜 1 填0置1法如右图所示结点1 3处为固定支座 可知将整体刚度矩阵中与之相对应的主对角元素全部置换成1 相应行和列上的其它元素均改为0 同时 所在同一行上的载荷分量替换成0 则有 第三章杆系结构静力分析的有限单元法 51 则 第三章杆系结构静力分析的有限单元法 也可简便地采用划行划列的办法 在整体刚度矩阵中将与约束位移为0的行和列划掉 包括相关的所在行的位移和载荷向量 52 处理后得基本平衡方程 2 乘大数法右图所示刚架 结点1为固定支座 结点3处在方向的约束为已知强迫位移 即将整体刚度矩阵中与之相对应的主对角元素全部乘以一个大数N 一般取 同时 将相应同一行上的载荷分量替换成N乘以其主对角刚度系数和给定的强迫位移 包括零位移 第三章杆系结构静力分析的有限单元法 53 得到 由于N足够大 可以近似认为 则得出 同时得到 求出位移之后 即可以求出结构的应力场 第三章杆系结构静力分析的有限单元法 54 第三章杆系结构静力分析的有限单元法 用有限单元法计算空间刚架结构 在原理上及推导过程与计算平面刚架结构相同 在此不再重复 但应注意到 由于空间的每一结点一般具有六个自由度 故计算较之复杂些 3 6计算示例设两杆的杆长和截面尺寸相同 杆件长m 图3 12刚架受力简图 55 结构离散化后将结构划分为4个结点 3个单元 截面积 惯性矩 2 求结点载荷首先须求局部坐标系中固定端内力 a 单元1作为两端固定梁反力示意图 b 单元2作为两端固定梁反力示意图图3 13内力示意图 第三章杆系结构静力分析的有限单元法 56 单元1 单元2 在局部坐标系下单元载荷列向量 单元1 单元2 单元3 第三章杆系结构静力分析的有限单元法 57 为了求出在整体坐标下的载荷列向量 先求单元得坐标转换矩阵 单元1 2 单元3 第三章杆系结构静力分析的有限单元法 58 求各单元在整体坐标下的等效结点载荷 第三章杆系结构静力分析的有限单元法 59 求刚架的等效结点载荷 第三章杆系结构静力分析的有限单元法 60 因为无结点载荷作用 总结点载荷即为等效结点载荷 3 求单元刚度矩阵由于单元1 2 3的尺寸相同 材料弹性模量相同 故 梁单元的局部坐标下的刚度矩阵表达式 第三章杆系结构静力分析的有限单元法 61 则 4 求整体坐标系中的 单元1 单元2 单元3 第三章杆系结构静力分析的有限单元法 62 5 求结构整体刚度矩阵 利用刚度集成法 6 建立原始平衡方程式 第三章杆系结构静力分析的有限单元法 63 7 引入约束条件解方程组 由于1 3 4为固定端 修改整体刚度矩阵中的1 3 6 12行与列 以及载荷列向量中的相应的行 既约束处理 建立基本平衡方程 即 得到 8 求各杆的杆端力 单元3结点位移列向量 第三章杆系结构静力分析的有限单元法 64 单元1杆端内力计算 单元2杆端内力计算 单元3杆端力计算 第三章杆系结构静力分析的有限单元法 65 9 作内力图 a 刚架轴力图 b 刚架剪力图 c 刚架轴弯矩图图3 14刚架内力图 第三章杆系结构静力分析的有限单元法 66 ANSYS的基本过程 一个典型的ANSYS分析过程可分为以下3个步骤 ANSYS应用分析 第三章杆系结构静力分析的有限单元法 67 1 前处理 前处理 指定工程名称和分析标题 定义单位 定义单元类型 定义单元常数 创建横截面 定义材料特性 创建有限元模型 定义分析类型求解控制 加载 68 1 1指定工程名称和分析标题 更改工程名 定义分析标题 69 1 2定义单位 使用 UNITS命令可以设置系统单位 没有相应的GUI USER 用户自定义单位 是缺省设置 SI 国际单位制 m kg s BFT 以英尺为基础的单位制 ft slug s F CGS cm g s MPA mm mg s BIN 以英寸为基础的单位制in lbm s F 70 1 3定义单元类型 71 1 4定义单元常数 单元实常数是由单元类型的特性决定的 如梁单元的横截面特性 并不是所有的单元类型都需要实常数 同类型的不同单元也可以有不同的实常数 指定单元的实常数号 72 1 5创建横截面 创建梁的横截面 73 1 6定义材料特性 定义材料特性 指定单元材料号 74 75 1 7定义分析类型求解控制 定义分析类型 求解控制 基本设置 瞬态设置 求解选项 非线性设置 求解终止的高级控制 76 包括 自由度约束 力 表面分布载荷 体积载荷 惯性载荷 耦合场载荷 载荷步 仅指可求得解的载荷设置 子步 是指在一个载荷步中每次增加的步长 主要是为了在瞬态分析和非线性分析中提高分析精度和收敛性 子步也称作时间步 代表一段时间 1 8加载 77 2 求解 求解当前载荷步 求解某载荷步 78 3 通用后处理器 画出分析的结果用列表的形式列出分析的结果查询某些结点或者单元处的应力值以及其它分析选项 79 DeformedShape 表示画出变形后的形状 有如下选项 3 1画出分析的结果 80 3 2画出节点的结果 81 位移 转角 3 3求解自由度结果 82 正应力和剪应力 主应力 应力强度 平均等效应力 3 4求解应力结果 83 正应变和剪应变 主应变 应变强度 平均等效应变 3 5求解总应变结果 84 求解能量 弹性应变 蠕变 其它应变 正应变和剪应变 主应变 应变强度 平均等效应变 3 6其它求解结果 85 3 7图形输出选项 只画出变形后的图形 画出变形前后的图形 画出变形后的图形和变形前的边界图 86 4 时间历程后处理器 适用于 瞬态分析的后处理 87 3 7ANSYS桁架结构计算示例 1m 1m 材料为Q235 1 选择单元类型 运行Preprocessor ElementType Add Edit Delete 在结点8上施加竖直向下的集中载荷F 60000N 约束为结点1处约束X Y方向自由度 结点5处约束Y方向自由度 图3 15桁架结构示意图图3 16桁架各单元横截面图 第三章杆系结构静力分析的有限单元法 图3 17单元类型对话框 88 图3 18单元类型库对话框 2 设置材料属性 运行Preprocessor MaterialProps MaterialModels 图3 19选择材料属性对话框图3 20设置材料1属性对话 3 设置单元截面形式 选择菜单Preprocessor Sections Beam CommonSections 图3 21梁截面设置对话框 第三章杆系结构静力分析的有限单元法 89 4 定义实常数 运行RealConstants Add Edit Delete 图3 22设置LINK1单元的实常数 5 建立模型 首先生成结点 运行主菜单Preprocessor Modeling Create Nodes InActiveCS 再生成单元 运行主菜单Preprocessor Modeling Create Elements AutoNumbered ThruNodes穿越结点命令 图3 23创建结点对话框 第三章杆系结构静力分析的有限单元法 90 图3 24通过结点建立单元 图3 25桁架的有限元模型 6 施加约束 运行主菜单Solution DefineLoads Apply Structural Displacement OnNodes 图3 26结点施加约束对话框 第三章杆系结构静力分析的有限单元法 91 7 施加载荷运行主菜单Solution DefineLoads Apply Structural Force Moment OnNodes 图3 27结点施加载荷对话框 8 求解运行主菜单Solution Solve CurrentLS 分析当前的负载步骤命令 弹出如图3 28所示对话框 单击OK 开始运行分析 分析完毕后 在信息窗口中提示计算完成 单击Close将其关闭 9 后处理运行主菜单GeneralPostproc PlotResults ContourPlot NodalSolu命令 运行DOFSolution Displacementvectorsum 出现桁架轴向应力云图 图3 29云图显示对话框 第三章杆系结构静力分析的有限单元法 图3 28求解对话框 92 图3 30位移云图 选择Stress vonMisesstress 则出现桁架位移云图 图3 31云图显示对话框 图3 32轴向应力云图 桁架的位移云图可知 最大位移发生在桁架的中部 最大位移为m 桁架的轴向应力云图可知 最大应力发生在2单元 最大应力45 9MPa 第三章杆系结构静力分析的有限单元法 93 3 8ANSYS刚架结构计算示例 图3 33刚架示意图 约束形式为 A D点施加全约束 在BC梁中点处受到竖直向下集中载荷的作用F1 20000N AB柱的中点处受水平向右的集中载荷F2 10000N AB 2m BC 2m 材料为钢材 弹性模量E 2 1 1011Pa 泊松比 0 3 1 选择分析范畴 图3 34选择分析范畴对话框 在主菜单中单击Preferences菜单 弹出PreferencesforGUIFiltering窗口 选择Structural 然后单击OK按钮 第三章杆系结构静力分析的有限单元法 94 2 选择单元类型运行Preprocessor ElementType Add Edit Delete 弹出ElementTypes对话框 选择BEAM188单元 图3 35单元类型对话框 图3 36单元类型库对话框 3 设置单元截面形式运行Preprocessor Section Beam CommonSections 弹出BeamTool对话框 W1选项栏中填写0 1 W2选项栏中填写0 2 t1 t4中填写0 008 设置完毕单击OK按钮 图3 37梁截面设置对话框 第三章杆系结构静力分析的有限单元法 95 4 设置材料属性运行Preprocessor MaterialProps MaterialModels 弹出DefineMaterialModelBehavior对话框 双击Isotropic选项 弹出LinearIsotropicPropertiesforMaterialNumber1对话框 在EX选项栏中设置数值2 1e11 在PRXY选项栏中设置数值0 3 设置完毕单击OK按钮 图3 38选择材料属性对话框 图3 39设置材料属性对话框 第三章杆系结构静力分析的有限单元法 96 5 建立模型设置材料属性对话框运行Preprocessor Modeling Create Keypoints InActiveCS 创建关键点1 在NPT选项栏中设置数值1 表示设置的关键点号为1 在X Y Z栏中设置数值0 0 0 表示关键点1的坐标为 0 0 0 同理设置关键点2 3 4 坐标分别为 0 2 0 2 2 0 2 0 0 运行Preprocessor Modeling Create Lines Lines StraightLine 弹出CreateStraightLine对话框 分别拾取点1 2 3 4 2 3 并经过布尔运算将两直线相加 图3 40创建关键点对话框 图3 41创建直线对话框 第三章杆系结构静力分析的有限单元法 97 图3 42刚架模型 6 划分网格选择刚架的单元属性 运行Preprocessor Meshing MeshAttributes Pickedlines 弹出LineAttributes对话框 拾取刚架后弹出MeshingAttributes对话框 采取默认设置 点击OK 图3 43划分网格拾取线对话框 图3 44设置网格单元属性 第三章杆系结构静力分析的有限单元法 98 选择Preprocessor Meshing SizeCntrls ManualSize Lines PickedLines 选择刚架 弹出对话框 在NDIV一栏中输入30 单击OK 最后在MeshTool中自由划分网格 图3 45定义单元尺寸拾取线对话框 图3 46设置线上单元尺寸对话框 7 施加约束运行Solution defineLoads Apply Structure Displacement OnKeypoints 选择关键点1 选择ALLDOF 同理对关键点4进行全约束 第三章杆系结构静力分析的有限单元法 99 图3 47对关键点施加全约束 8 施加载荷 将图形结点显示 运行PlotCtrls Numbering 激活NodeNumbers后面的选框 使它变成on形式 选择菜单Solution DefineLoads Apply Structure Force Moment OnNodes 拾取结点17 施加集中载荷Fy 20000N 同理 在结点7上施加集中载荷Fx 10000N 图3 48编号显示设置对话框 第三章杆系结构静力分析的有限单元法 100 图3 49拾取点对话框图3 50施加约束 载荷后的刚架有限元模型 9 求解选择Solution Solve CurrentLS 弹出如右图所示对话框 单击OK按钮 开始计算 计算结束会弹出计算完毕对话框 单击Close关闭对话框 计算完毕 图3 51求解对话框 第三章杆系结构静力分析的有限单元法 101 10 后处理显示位移云图 运行GeneralPostproc PlotResults ContourPlot NodalSolu 弹出如下图所示对话框 运行DOFSolution Displacementvectorsum显示刚架位移云图 图3 52云图显示对话框图3 53刚架的位移云图 第三章杆系结构静力分析的有限单元法 102 4 1平面应力问题 第四章平面结构问题的有限单元法 4 2平面应变问题 4 3平面问题的离散化 4 4平面三结点三角形单元 4 5ANSYS平面结构计算示例 103 严格地说 任何弹性体都是处于三维受力状态 因而都是空间问题 但是在一定条件下 许多空间问题都可以简化成平面问题 平面问题可以分为两类 平面应力问题和平面应变问题 图4 1平面问题应力状态 第四章平面结构问题的有限单元法 104 4 1平面应力问题 图4 2 a 平面应力问题 如图所示的深梁结构 其厚度方向的尺寸远比其它两个方向的尺寸小得多 可视为一薄板 它只承受作用在其平面内的载荷 且沿厚度方向不变 计算时以中性面为研究对象 其力学特点是 平面应力问题的应力应变转换矩阵即弹性矩阵为 第四章平面结构问题的有限单元法 105 图4 2 b 平面应变问题 4 2平面应变问题 图示为一圆形涵洞的横截面 其长度方向上的尺寸远比其它两个方向上的尺寸大得多 同样 载荷作用在xy坐标面内 且沿z轴方向均匀分布 其力学特点是 但一般情况下 平面应变问题的弹性矩阵只需将式 4 1 中的E换成 换成 即可 第四章平面结构问题的有限单元法 106 无论是平面应力问题还是平面应变问题的应力 与应变 之间的关系均为 其中 为初应变 式中 4 3平面问题的离散化 a 三结点三角形单元 b 四结点正方形单元 c 四结点矩形单元 d 四结点四边形单元 图4 3平面问题单元的主要类型 第四章平面结构问题的有限单元法 107 图4 4 a 表示的是带有椭圆孔的平板 在均匀压力作用下的应力集中问题 图4 5 b 是利用结构的对称性 采用三结点三角形单元而离散后的力学模型 各单元之间以结点相连 a 均匀受力板力学模型 b 力学模型离散化 图4 4平面问题有限单元法的计算力学模型 第四章平面结构问题的有限单元法 108 4 4平面三结点三角形单元 4 1 1位移函数 图4 5三角形单元 如果把弹性体离散成为有限个单元体 而且单元很小时 就很容易利用其结点的位移 构造出单元的位移插值函数 即位移函数 位移函数矩阵形式 第四章平面结构问题的有限单元法 109 简写为 由于位移函数适用于单元中的任意一点 所以带入3个结点的坐标后 得出结点处位移函数为 简写为 第四章平面结构问题的有限单元法 110 解出 其中 是三角形单元的面积 当三角形单元结点i j m按逆时针次序排列时 则有 4 4 2形函数矩阵 第四章平面结构问题的有限单元法 111 其中记号 表示将i j m进行轮换后 可得出另外两组带脚标的a b c的公式 单元位移函数为结点位移的插值函数 即 第四章平面结构问题的有限单元法 4 9 112 令 在式 4 10 中表示的称为形函数 于是位移函数表达式用形函数表示为 4 10 4 11 写成矩阵形式 4 12 第四章平面结构问题的有限单元法 113 由几何方程知 将式 4 9 代入式 4 13 中 并求偏导数 得 4 13 4 4 3单元的应力与应变 第四章平面结构问题的有限单元法 114 简写为 4 14 由于 B 是常量 单元内各点应变分量也都是常量 这是由于采用了线性位移函数的缘故 这种单元称为常应变三角形单元 4 15 第四章平面结构问题的有限单元法 115 由弹性力学的物理方程可知 其应力与应变有如下关系 4 16 将式 4 14 代入式 4 16 得 4 17 式中 4 18 S 称为应力转换矩阵 对平面应力问题 其子矩阵为 4 19 由式 4 17 看出 应力分量也是一个常量 在一个三角形单元中各点应力相同 一般用形心一点表示 其应变也可同样表示 第四章平面结构问题的有限单元法 116 用虚功原理来建立结点力和结点位移间的关系式 从而得出三角形单元的刚度矩阵 a 实际力系 b 虚设位移 图4 6弹性体虚功原理的应用 4 4 4三角形单元刚度矩阵 第四章平面结构问题的有限单元法 117 结点力列向量和应力列向量分别为 结点虚位移列向量和虚应变列向量为 用虚功原理建立三角形单元的虚功方程为 由式 4 12 式知 代入式 4 20 得 4 20 第四章平面结构问题的有限单元法 118 由于虚位移是任意的 等号两边可左乘 得 4 21 三角形单元的刚度矩阵可写成 4 22 用分块矩阵形式表示 4 23 第四章平面结构问题的有限单元法 119 结构的平衡条件可用所有结点的平衡条件表示 假定i结点为结构中的任一公共结点 则该结点平衡条件为 i结点的结点力列向量 围绕i结点所有单元的结点力的向量和 i结点的载荷列向量 4 4 5整体刚度矩阵 第四章平面结构问题的有限单元法 120 每个结点由两个平衡方程组成 若结构共有n个结点 则有2n个平衡方程 整个结构的平衡条件由式 4 24 求和得到 即 i 1 2 n 4 26 4 27 其中 K 为结构整体刚度矩阵 为结构的结点位移列向量 4 28 第四章平面结构问题的有限单元法 121 将式 4 26 式 4 27 代入式 4 25 中得 4 29 整体刚度矩阵也可按结点写成分块矩阵的形式 4 30 同杆系结构一样 整体刚度方程经过约束处理后 即可求出结点位移 进而求出所希望的应力场 第四章平面结构问题的有限单元法 122 4 5ANSYS平面结构计算示例 4 5 1问题描述 如图4 7所示长方形板ABCD 板厚0 04m 孔半径r 0 2m E 210GPa 泊松比 0 3 约束条件 在长方形底边AD约束全部自由度 BC边施加垂直向下均布载荷g 10000000N m 图4 7长方形板结构 4 5 2ANSYS求解操作过程 打开Ansys软件 在Ansys环境下做如下操作 第四章平面结构问题的有限单元法 123 图4 8单元类型对话框 1 选择单元类型运行Preprocessor ElementType Add Edit Delete 弹出ElementTypes对话框 如图4 8所示 单击Add 弹出LibraryofElementTypes窗口 如图4 9所示 依次选择StructuralSolid Quad8node82 单击OK 图4 9单元类型库对话框 第四章平面结构问题的有限单元法 124 在ElementTypes对话框中 如图4 10所示 单击Options 弹出如图4 11所示对话框 设置K3选项栏为Planestrsw thk 设置K5选项栏为Nodalstress 设置K6选项栏为Noextraoutput 表示单元是应用于平面应力问题 且单元是有厚度的 图4 10单元类型对话框 图4 11PLANE82单元选项设置对话框 第四章平面结构问题的有限单元法 125 2 定义实常数运行Preprocessor RealConstants Add Edit Delete 弹出如图4 12所示对话框 点击Add 弹出如图4 13所示对话框 点击OK 弹出如图4 14所示对话框 在THK选项栏中设置板厚度为0 04m 设置完毕单击OK按钮 图4 12实常数对话框 图4 13选择要设置实常数的单元类型 图4 14PLANE82实常数设置 第四章平面结构问题的有限单元法 126 3 设置材料属性运行Preprocessor MaterialProps MaterialModels 弹出如图4 15所示对话框 依次双击Structural Linear Elastic Isotropic 弹出图4 16所示对话框 在EX选项栏中设置数值2 1e11 在PRXY选项栏中设置数值0 3 设置完毕单击OK按钮 图4 15选择材料属性对话框 图4 16设置材料属性对话框 第四章平面结构问题的有限单元法 127 4 建立模型运行Preprocessor Modeling Create Area Rectangle By2Corners 弹出如图4 17所示对话框 设置参数 WPX选项栏中填写0 WPY选项栏中填写0 Width选项栏中填写1 5 Height选项栏中填写1 单击OK 继续运行Preprocessor Modeling Create Area Circle SolidCircle 得到如图4 18所示对话框 在WPX选项栏中填写0 75 WPY选项栏中填写0 5 在Radius选项栏中填写0 2 设置完毕点击OK按钮 图4 17建立矩形对话框 图4 18创建实心圆对话框 第四章平面结构问题的有限单元法 128 进行布尔运算 Preprocessor Modeling Operate Booleans Subtract Areas 先选矩形面单击OK 再单击圆面 单击OK 得到如图4 19所示图形 图4 19长方形板模型 5 划分网格运行Meshing SizeCntrls ManualSize Areas AllAreas 弹出如图4 20所示对话框 在SIZE选项栏中填写0 05 点击OK按钮 图4 20设置网格尺寸对话框 第四章平面结构问题的有限单元法 129 运行Mesh MeshTool 弹出如图4 21所示对话框 在Shape选项栏后面 选择Tri和Free 单击Mesh 划分网格 网格划分如图4 22所示 图4 21网格划分对话框 图4 22划分网格后的有限元模型 第四章平面结构问题的有限单元法 130 图4 23施加全约束 6 施加约束选择菜单Solution DefineLoads Apply Structure Displacement OnLines 选择长方形底边 弹出图4 23所示对话框 选择AllDOF 单击OK 7 施加载荷选择菜单Solution DefineLoads Apply Structure Pressure OnLines 弹出如图4 24所示对话框 拾取长方形上边 单击OK按钮 弹出如图4 25所示对话框 在VALUE选项栏中填写10000000 设置完毕点击OK完成设置 第四章平面结构问题的有限单元法 131 图4 24拾取要施加载荷的边 图4 25施加载荷对话框 8 求解运行Solution Solve CurrentLS 弹出如图4 26所示对话框 单击OK按钮 开始计算 计算结束会弹出计算完毕对话框 单击Close关闭对话框 计算完毕 图4 26求解当前步载荷对话框 第四章平面结构问题的有限单元法 132 9 后处理运行GeneralPostproc PlotResults ContourPlot NodalSolu 弹出如图4 27所示对话框 运行DOFSolution Displacementvectorsum和Stress vonMisesstress 分别显示长方形面板的位移云图和应力云图 结果显示如图4 28和图4 29所示 图4 27云图显示对话框 第四章平面结构问题的有限单元法 133 图4 28位移变形云图 图4 29等效应力云图 4 5 3结论 从图4 28长方形面板的位移云图可知 最大位移发生在圆孔的上部 最大位移为0 750 10 4m 从图4 29长方形面板的应力云图可知 最大应力发生在圆孔的两侧 最大应力为32 9MPa 第四章平面结构问题的有限单元法 134 5 1坐标变换与平面四结点等参元 第五章等参元 5 2平面八结点等参单元 5 3单元刚度矩阵 5 4ANSYS等参元计算示例 135 5 1坐标变换与平面四结点等参元 图5 1 a 为一个任意四边形单元 称为实际单元 在实际单元内以对边的中点连线建立起一个局部坐标系 通过坐标转换把实际单元 映射 为如图5 1 b 所示的一个正方形 此坐标系称为单元的自然坐标系或等参数坐标系 正方形称为基本单元 基本单元内任一点P 与实际单元内的一点P x y 唯一对应 a 直角坐标系与实际单元 b 自然坐标系与基本单元 图5 1四结点等参单元 第五章等参元 136 实际单元与基本单元的对应关系可写为 或 其中 用同样的形状函数来插值单元内任意一点 x y 的位移 5 1 5 2 第五章等参元 137 为此单元的结点位移列向量 为形状函数矩阵 这里采用了同样的形状函数式 5 2 用同样的结点插值函数表示出单元的几何坐标x y与u v 这种单元称为等参单元 5 3 即 第五章等参元 138 5 2平面八结点等参元 类似地可以推广到具有更多结点的单元 如图5 2所示 a 直角坐标系与实际单元 b 自然坐标系与基本单元 图5 2八结点等参单元 该基本单元的位移函数可取为 5 4 第五章等参元 139 其中在顶角结点与边中点上的形函数分别为 5 6 5 5 第五章等参元 140 5 3单元刚度矩阵 首先给出单元内的应变列向量 对平面问题 应有 5 7 按坐标变换关系式 5 1 有 第五章等参元 141 写成矩阵表达式为 5 8 由式 5 8 可解出 其中 称为坐标变换的雅可比 Jacabian 矩阵 其中 5 9 第五章等参元 142 合写成矩阵形式有 5 10 将式 5 3 代入式 5 7 中 则有 为应变转换矩阵 按结点分块表示 有 第五章等参元 143 而i 1 2 3 4 5 11 将式 5 9 代入式 5 11 即可得出此单元的应变转换矩阵 进而求出 同上 单元内的应力可表示为 单元刚度矩阵由虚功原理求得 即 第五章等参元 144 5 12 上述积分在自然坐标系内进行 得 刚度矩阵 5 13 一般参数单元的计算都采用数值积分求式 5 13 的近似值 同时 为了减少计算点的数目和便于编写程序 多采用高斯数值积分方法 二维积分法的高斯求积公式为 5 14 第五章等参元 145 式中 为对应的坐标位置 值 为权重系数 L M为沿 方向的积分点数目 5 4ANSYS平面结构计算示例 5 4 1问题描述 一个长方形面板 如图5 3所示 其高AB 1m 宽BC 1 5m 板厚b 0 04 孔半径R 0 2m 长方形面板的弹性模量E 210GPa 泊松比 0 3 约束条件 在长方形底边约束全部自由度 载荷 BC边施加垂直向下均布载荷F 10000000N m 图5 3长方形板结构 第五章等参元 146 5 4 2ANSYS求解操作过程 打开Ansys软件 在Ansys环境下做如下操作 1 选择单元类型运行Preprocessor ElementType Add Edit Delete 弹出ElementType对话框 如图5 4所示 单击Add 弹出LibraryofElementType窗口 如图5 5所示 选择PLANE82 图5 4单元类型对话框 图5 5单元类型库对话框 第五章等参元 147 在ElementTypes对话框中单击Options对话框 弹出如图5 6所示对话框 设置K3选项栏为Planestrsw thk 设置K5选项栏为Nodalstress 设置K6选项栏为Noextraoutput 表示单元是应用于平面应力问题 且单元是有厚度的 图5 6PLANE82单元类型选项对话框 运行Preprocessor RealConstants Add Edit Delete 弹出如图5 7所示对话框 点击Add 弹出如图5 8所示对话框 点击OK 弹出如图5 9所示对话框 在THK选项栏中设置板厚度为0 04m 设置完毕单击OK按钮完成设置 第五章等参元 148 2 设置材料属性运行Preprocessor MaterialProps Materi