fem-model文档.pptVIP

Univ. Of Petroleum 有限元数值模拟分析 理论与方法 薛世峰 Sep.2006 数学、力学模型方程 近似数值解的概念 数值方法的分类 有限元法、有限差分法、边界元法 Galerkin 的基本概念 Part1、有限元方法概述 FEPG 系统简介 FEPG 系统入门、应用过程 FEPG 系统高级技术研讨 流固耦合模型数值软件技术 Part 2 、Galerkin有限元与FEPG自动生成技术 解析解 半解析解 数值解 级数解 FEM FDM BEM 外界输入 特定属性的对象或系统 系统响应 数学力学模型 定量描述三部分相互作用与联系 质量、动量、能量的一般规律； 系统或对象物理与力学属性的规律； 补充方程与定解条件。 模型方程一般包括： 工程力学问题主要包括三类方程： 平衡方程或运动方程； 物理性质方程； 补充方程； P 代数方程、常微分方程、偏微分方程： 三类方程的区别： 对复杂变化规律的描述程度不同； 不同的方程性质； 常微分方程、偏微分方程还需要定解条件，如： 边界条件、初始条件等。 求解的难易程度不同。 大部分工程力学涉及的数学方程为偏微分方程： 1：弹性力学平衡方程 2：线弹性应力应变关系 3：变形几何方 程 如：弹性力学问题，包括 偏微分方程分类： 与时间无关的椭圆型方程 与时间相关的抛物型、双曲型方程 小结与讨论： 工程力学中的问题相当数量属偏微分方程类型； 偏微分方程的理论解析解是很难得到的； 近似的数值解是必需的，其技术也是比较复杂的。 能否将复杂偏微分方程的求解在一定力学意 义上等效转化为对代数方程（组）的求解 ？ 数值解的理论基础： 如何将原来偏微分方程等效转换为代数方程； 数值解方法具体实施过程中若干技术问题； 数值解结果评价指标，即精度问题。 等价形式的能量方程 等价形式的积分方程 偏微分方程（组） 1、弹性力学问题的等价能量积分方程形式 变形能 外力功 fx fy Ti 2、等价形式的积分方程 以瞬态热传导方程为例，微分形式 为温度场，需在整个区域内严格满足方程关系，并与边界条件、初始条件吻合； 为权函数，一般取为一Lagrangian线性多项式： 其中 注解： 1、T(x,y)为温度场分布，需在域内近似满足； 2、T(x,y)为二阶连续可微函数； 3、在全场域内严格满足非常困难。 积分形式 分步积分公式： 区域满足 （子区域满足） i 有限元数值计算 1、连续介质区域离散为有限自由度的网格结构； 2、在每一单元上建立积分方程的近似式； 3、将各个单元上的代数式组装起来，形成总体代数方程。 简化为：结构离散、单元分析、结构总装 单元特性分析 i j k (xi, yi) (xj, yj) (xk, yk) 一次或二次近似 Physical Problem Mathematical Model Numerical model Does answer make sense? Refine analysis Happy ? YES! No! Improve mathematical model Design improvements Structural optimization Change physical problem Pre-Processor Initial State Results Dataset Boundary conditions updating Deformation Stress , strain Solid Module Fluid Module shear strain Pressure, sanding volume Flux , porosity Post-Processor Visualized Results Display Coupled code analysis chart Time Recorder T=T+DT Part 2 、Galerkin有限元与FEPG自动 生成技术 最经典最成熟的方法 从工程角度的理解把实际结构划分成 很多细小的单元 用比较简单的数学模型近似单元内部的 应力和变形 将这些单元组合，通过求解刚度矩阵得 到整体结构的反应 1、FEM 评述 对于线弹性问题，当实际结构位移场函数 连续光滑时，能够得到收敛解 对于任意复杂结构，理论上总是可以通过细 分单元的方法获得足够近似的模拟 刚度矩阵系数带状，在结构不出现软化的时 候还是对称正定的，求解方便 长期大量工程应用，积累了丰富的经验 2、有限元法的优点 非线性的数值计算是很复杂的，很难为一般工程技术人员所掌握。 为此近年来国外一些公司花费了大量的人力和投资开发诸