哈工大 有限元理论.pptVIP

有限元方法（优点） 有限元可以运用于任何场问题： 没有几何形状的限制 边界条件和载荷没有限制 材料性质并不限于各向同性 具有不同行为和不同数学描述的分量可以结合。 有限元结构和被分析的物体或区域很类似 通过网格细分可以很容易地改善解的逼近度。 1. 绪论 　　　　　有限元方法（应用领域） 机械/航空航天/土木工程/自动化工程 结构分析（静力/动力分析，线形/非线性分析） 热分析/流体力学分析 电磁场分析 地质力学分析 生物医学分析 1. 绪论 有限元方法（计算机程序） 专业化小程序有 Truss Frame Plane Stress Heat Transfer 大型通用商业化程序 ANSYS ADINA ABAQUS PANTRAN SAP 1. 绪论 有限元方法（分析步骤） 第一步：问题及求解域定义 第二步：求解域离散化 第三步：确定状态变量及控制方法 第四步：单元推导 第五步：总装求解。 第六步：联立方程组求解和结果解释 　　简言之，有限元分析可分成三个阶段，前处理、处理和后处理。前处理是建立有限元模型，完成单元网格划分;后处理则是采集处理分析结果，使用户能简便提取信息，了解计算结果 1. 绪论 举例：杆系结构有限元分析的完整过程 4.1有限元分析的完整过程 E1=E2=2E7Pa A1=A2=2cm2 l1=l2=10cm P3为10N作用下二杆结构的变形。 1）离散化 该构件由两根杆件做成，因此可以自然离散成2个杆单元。 单元1; 节点位移 u1，u2； 节点力 P1，P2 （U1，U2）。 4.1有限元分析的完整过程 2）单元分析（直接法分析） 状态一：节点2固定 单元节点条件：u(1)=u1， u(2)=0 单元内应力： 材料力学：拉应力为正，压为负； 4.1有限元分析的完整过程 2）单元分析（直接法分析） 状态一：节点2固定 单元节点条件：u(1)=u1， u(2)=0 单元左端点力： 由杆单元静力平衡条件; 4.1有限元分析的完整过程 有限元：x正方向的节点力为正 2）单元分析（直接法分析） 状态二：节点1固定 单元节点条件：u(1)=0， u(2)=u2 单元左端点力： 由杆单元静力平衡条件; 4.1有限元分析的完整过程 把上述两种状态相加： 上脚标表示单元；下脚标表示节点 单元1刚度矩阵 4.1有限元分析的完整过程 4.1有限元分析的完整过程 同理可得单元2： 4.1有限元分析的完整过程 综合2单元，考虑约束： 节点1： 节点2： 节点3： 4.1有限元分析的完整过程 写成矩阵 整体刚度矩阵： 4.1有限元分析的完整过程 对称，不是偶然，满足功的互等定理。 行列式=0，奇异，物理意义无约束有位移 =0 5）建立刚度方程 简化为二阶 引入约束过程 4.1有限元分析的完整过程 6）求解节点位移 将结构参数和外载荷代入上式有 求解得（单位m） 4.1有限元分析的完整过程 9）计算支反力 具体地对于单元1，有 其中R1是节点1的支反力，P2是单元1的节点2所受的力，即单元2对该 节点的作用力，将前面求得的节点位移代入上式可得支反力大小。 4.1有限元分析的完整过程 以上是一个简单结构有限元方法求解得完整过程，对于 复杂结构，其求解过程完全相同，由于每一个步骤都具 备标准化和规范性的特征，所以可以在计算机上编程而 自动实现。 4.1有限元分析的完整过程 讨论2：由前面的步骤，我们也可以直接将各个单元的刚度矩阵按照节 点编号的对应位置来进行装配，即在未处理边界条件之前，先形成整 体刚度矩阵。 其物理意义是，表示在未处理边界条件前的基于节点描述的总体平衡 关系。在对该方程进行位移边界条件的处理后就可以求解，这样与先 处理边界条件再求系统势能的最小值所获得的方程完全相同。 4.1有限元分析的完整过程 有限元分析的基本步骤及表达式 1、物体几何区域的离散化 2、单元的研究（所有力学信息都用节点位移）来表达 3、装配集成 4、边界条件的处理并求解节点位移 5、支反力的求取以及其它力学量（应力、应变及位移三大物理量）的计算 第一章 概述 * 有限单元法基本步骤 (1) 待求解域离散化 (2) 单元分析，插值位移函数 (3) 形成单元性质的矩阵方程 (4) 形成整体系统的矩阵方程 (5) 约束处理，求解系统方程 (6)