《有限元方法》课件

有限元方法概论一种将复杂问题分解为简单部分的数值分析技术世纪年代发展，现已成为工程计算基石2060广泛应用于结构、流体、热、电磁、生物等领域有限元方法发展历程1960年代1990年代飞机结构分析应用图形界面普及12341970年代2000年后商业软件出现多物理场耦合分析有限元法基本思想整体近似通过子区域拼接获得系统响应区域剖分将复杂几何分为简单单元局部分析每个子区域内进行简化求解典型问题举例固体力学结构变形与应力分析热传导温度场与热流计算电磁场电场、磁场分布模拟流体力学流场与压力分布计算有限元与传统数值方法比较方法类型几何适应性精度计算量有限元法极强高较大有限差分法一般中等中等谱方法较弱极高较小有限元法数学基础泛函分析空间与函数逼近Sobolev变分原理最小势能原理应用矩阵理论线性方程组的求解技术平衡方程与结构力学基础外力作用内力产生1集中力、分布荷载力的传递与平衡2平衡状态结构响应43静力平衡方程满足变形与应力发展变分原理简介变分原理Ritz法Galerkin法基于泛函极值寻找解通过试函数展开近似解加权残值法典型应用物理系统趋于能量最小状态最小化能量泛函确定系数与法等价（自伴算子）Ritz有限元离散化概念几何离散复杂几何划分为单元节点定义单元连接点建立场函数近似形函数构造连续解基本单元类型形函数（插值函数）原理形函数定义连接节点值的插值多项式节点特性本节点值为，其他节点为10连续性保证单元间满足位移连续阶次选择一次形函数线性变化二次形函数曲线变化节点与网格生成结构化网格非结构化网格自适应网格规则几何区域适用复杂几何形状适用根据计算结果优化网格质量控制避免尖锐角和扁平单元有限元逼近与收敛性误差来源收敛率离散化误差单元尺寸减小，误差降低数值积分误差次元素，收敛率为p Oh^pCéa引理有限元解与精确解误差上界有限元基本流程前处理几何建模与网格划分求解矩阵组装与方程求解后处理结果可视化与分析弱形式推导强形式弱形式微分方程本身加权积分形式••要求高阶导数存在降低连续性要求••边界条件显式满足自然边界条件隐含••一维问题有限元建模21节点数自由度/节点一维线性单元节点数简单拉伸杆单节点自由度EA/L刚度系数单元刚度矩阵元素典型一维单元举例轴向杆单元梁单元弹簧单元承受拉压，不承受弯曲可承受弯曲和剪切模拟弹性连接二维问题有限元建模平面应力问题平面应变问题轴对称问题薄板问题厚体截面问题围绕轴线旋转体σz=τxz=τyz=0εz=γxz=γyz=0圆柱坐标系描述三角形单元详解节点配置形函数节点（线性）或节点（二次）线性或二次多项式36数值积分自然坐标点或点高斯积分面积坐标系表示13四边形单元详解形状与节点节点（双线性）或节点（二次）48/9等参单元使用相同形函数映射几何和场变量数值积分或高斯积分点2×23×3优势比三角形单元精度更高三维问题有限元分析四面体单元六面体单元节点（线性）节点（线性）•4•8节点（二次）节点（二次）•10•20适合复杂几何计算效率高••棱柱单元节点（线性）•6过渡层应用•特殊区域填充•建立单元刚度矩阵单元力向量组装集中力分布力边界力直接施加到节点通过形函数积分转为等效节点力边界条件产生的等效力映射到全局自由度面荷载位移、温度等约束•体积力•惯性力•总刚度矩阵组装单元贡献1根据节点编号映射重叠区域叠加共享节点累加刚度稀疏矩阵形成大量零元素高效存储有限元方程求解直接法高斯消元法•分解•LU分解•Cholesky迭代法共轭梯度法•法•GMRES大规模问题适用•特殊技术稀疏矩阵存储•带宽优化•并行求解•计算结果的处理位移场应力云图矢量图变形后形状与放大系数多种应力分量可视化流场、热流方向展示误差分析离散化误差数值积分误差网格尺寸与近似阶次限制高斯点数不足导致计算误差建模误差舍入误差与条件数影响几何简化与材料模型选择收敛性与精度改进h-细化p-加密hp-自适应增加单元数量提高形函数阶次结合两种方法保持单元阶次不变保持单元数量不变区域选择性应用收敛率指数级收敛最佳收敛效率Oh^p求解非线性问题几何非线性材料非线性大变形与旋转效应塑性、超弹性、蠕变接触非线性求解算法边界条件变化与摩擦牛顿拉夫森法与弧长法-动力学有限元分析热传导问题有限元法热传导方程边界条件温度梯度与热流关系指定温度•指定热流•瞬态与稳态情况对流换热•辐射换热•协同场问题直接耦合同时求解所有场方程序贯耦合交替求解各场方程单向耦合一个场影响另一场接触问题有限元分析点到点接触面到面接触接触模型节点与节点之间接触表面间的互相作用无摩擦、库伦摩擦数值方法罚函数法和拉格朗日乘子法有限元软件介绍有限元软件流程演示几何建模内部建模或导入模型CAD网格划分自动手动网格划分与检查/设置分析材料、边界条件、求解设置求解与后处理监控求解与结果可视化工程应用案例一机械零件分析齿轮建模应力分析变形分析精细齿形与接触区网格细化齿根应力集中评估啮合精度和刚度评估工程应用案例二桥梁结构分析5000+75%单元数量安全系数复杂桥梁分析所需单元极端荷载条件下的设计预留

3.2Hz固有频率首阶固有频率需大于风振频率工程应用案例三航天器结构分析轻量化设计动力学分析拓扑优化技术发射振动环境••减重目标微重力条件•30%•结构强度保证热变形响应••碰撞仿真着陆缓冲设计•空间碎片防护•部件脱落预测•工程应用案例四生物医学骨骼模型植入物分析个性化治疗由扫描数据重建界面接触与应力分布基于患者数据定制设计CT多相材料属性分配骨整合预测模型术前规划与优化有限元建模注意事项网格质量控制避免高畸变单元边界条件合理性反映真实约束状态简化策略保留关键特征的同时减少计算量验证与确认与解析解和试验对比有限元前处理（建模）几何导入几何清理模型格式转换与修复小特征抑制与简化CAD网格生成几何分区控制网格尺寸与质量便于结构化网格划分有限元后处理（结果分析）云图展示切面视图动画效果场变量分布可视化内部结构应力分析变形过程与模态形状编写有限元代码基础数据结构设计节点、单元、材料信息存储单元函数实现形函数与刚度矩阵计算全局组装稀疏矩阵高效算法求解器开发直接法或迭代法实现开源有限元平台介绍新发展与挑战多尺度分析机器学习结合实时仿真原子到宏观跨尺度模拟神经网络加速计算模型简化技术分子动力学与有限元耦合数据驱动参数优化手术导航与训练模拟计算资源与并行技术有限元法未来展望智能辅助决策自动推荐优化方案数字孪生实时映射物理系统云端协同计算分布式资源实时调度大规模并行计算百亿自由度问题求解国内外有限元发展现状国家地区主要软件应用领域发展方向/中国航空航天、土自主软件TFEA,JIFEX木美国全领域智能化ANSYS,NASTRAN欧洲汽车、能源开源生态ABAQUS,Code_Aster日本电子、材料多尺度MSC,J-OCTA课程复习与提纲总结理论基础离散化方法变分原理与插值近似单元类型与形函数2工程应用矩阵方程3软件使用与案例分析组装与求解技术结束与问题答疑学习资源联系方式《有限元基础》周志成邮箱•professor@university.edu《有限元分析理论》王勖成•办公室工程楼区A305《》•The FiniteElement MethodZienkiewicz答疑时间周二四下午点/3-5。