《大学有限单元法幻灯片》由会员分享，可在线阅读，更多相关《大学有限单元法幻灯片（96页珍藏版）》请在金锄头文库上搜索。
1、第一章 绪论,计算机辅助工程（CAE）是一种迅速发展的信息技术，是实现重大工程和工业产品的计算分析、模拟仿真与优化设计的工程软件，是支持工程科学家进行创新研究和工程师进行创新设计的、最重要的工具和手段。,1.1 有限单元法,有限单元法（以下简称为有限元法）是在当今工程分析中获得最广泛应用的数值计算方法。 将一个表示结构或连续体的求解域离散为有限个简单单元的组合，并通过它们边界上的结点相互联结成为组合体，从而将一个连续的无限自由度的问题简化为离散的有限自由度问题。,真实系统,有限元模型,有限元模型由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 节点具有一定的自由度。,齿轮有限元
2、模型,1.2 有限单元法的发展历史,1.3 有限单元法的应用,土木工程,轴承强度分析,机械工程,机械工程,三维椭圆封头开孔补强,车辆工程,航空工程,生物工程,核电工程,冶金工程,钢板精轧机热轧制有限元分析,1.4 有限元常用软件,从20世纪70年代开始，基于有限元法在结构线性分析方面已经成熟，并被工程界广泛采用，一批由专业软件公司研制的大型通用商业软件（如NASTRAN，ASKA，SAP，ANSYS，MARC，ABAQUS，JIFEX等）公开发行和被应用。它包含众多的单元型式、材料模型及分析功能，并具有网格自动划分、结果分析和显示等先后处理功能。,ANSYS 简介,1970 年成立的美国ANS
3、YS公司是世界CAE行业最著名的公司之一,长期以来一直致力于设计分析软件的开发、研制,其先进的技术及高质量的产品赢得了业界的广泛认可。 在我国,ANSYS用户也越来越多,三峡工程、二滩 电站、黄河下游特大型公路斜拉桥、国家大剧院、浦东国际机场等在结构设计时都采用了ANSYS作为分析工具。,ANSYS中国 http：/ 热分析 电磁分析 流体分析 (CFD) 耦合场分析 - 多物理场,17,结构分析用于确定结构的变形、应变、应力及反作用力等.,结构分析是有限元分析方法最常用的一个应用领域。结构这个术语是一个广义的概念，它包括土木工程结构，如桥梁和建筑物；汽车结构，如车身骨架；海洋结构，如船舶结构
4、；航空结构，如飞机机身等；同时还包括机械零部件，如活塞，传动轴等等。,ANSYS 结构分析,结构分析的类型: 静力分析 - 用于静态载荷. 可以考虑结构的线性及非线性行为，例如: 大变形、大应变、接触、塑性及蠕变等. 模态分析 - 计算线性结构的自振频率及振形. 谱分析 是模态分析的扩展，用于计算由于随机振动引起的结构应力和应变。 谐响应分析 - 确定线性结构对随时间按正弦曲线变化的载荷的响应. 瞬态动力学分析 - 确定结构对随时间任意变化的载荷的响应. 特征屈曲分析 - 用于计算线性屈曲载荷并确定屈曲模态形状. 专项分析: 断裂分析, 复合材料分析，疲劳分析,20,用于模拟冲击、碰撞、快速成
5、型等问题.,载有假人的汽车撞击刚性墙壁的有限元模拟,21,ANSYS热分析概览,热分析之后往往进行结构分析,计算由于热膨胀或收缩不均匀引起的应力. ANSYS功能: 相变 (熔化及凝固), 内热源 (例如电阻发热等) 三种热传递方式 (热传导、热对流、热辐射),ANSYS 热分析计算物体的稳态或瞬态温度分布，以及热量的获取或损失、热梯度、热通量等.,工件淬火3.06 min 时的温度分布 (NSHT3D),23,ANSYS电磁分析概览,磁场分析中考虑的物理量是磁通量密度、磁场密度、磁力、磁力矩、阻抗、电感、涡流、能耗及磁通量泄漏等. 磁场可由电流、永磁体、外加磁场等产生.,磁场分析 用于计算磁
6、场.,24,磁场分析的类型: 静磁场分析 - 计算直流电（DC)或永磁体产生的磁场. 交变磁场分析 - 计算由于交流电(AC)产生的磁场. 瞬态磁场分析- 计算随时间随机变化的电流或外界引起的磁场.,电磁接触：磁悬浮列车仿真,25,电场分析 用于计算电阻或电容系统的电场. 典型的物理量有电流密度、电荷密度、电场及电阻热等. 高频电磁场分析 用于微波及RF无源组件，波导、雷达系统、同轴连接器等分析.,26,ANSYS 流体分析 概览,流体分析 用于确定流体的流动及热行为. 流体分析分以下几类: CFD - ANSYS/FLOTRAN 提供强大的计算流体动力学分析功能，包括不可压缩或可压缩流体、层
7、流及湍流，以及多组份流等. 声学分析 - 考虑流体介质与周围固体的相互作用, 进行声波传递或水下结构的动力学分析等. 流体动力学耦合分析 - 在考虑流体约束质量的动力响应基础上，在结构动力学分析中使用流体耦合单元.,速度,压力,超音速飞行压力分布,汽车气动分析,高速导弹气动,28,ANSYS 耦合场分析概览,耦合场分析 考虑两个或多个物理场之间的相互作用。如果两个物理场之间相互影响，单独求解一个物理场是不可能得到正确结果的，因此你需要将两个物理场组合到一起求解。,其他需要耦合场分析的典型情况有： 热应力分析 流体结构相互作用 感应加热（电磁热）, 感应振荡,两根热膨胀系数不同的棒焊接在一起，图
8、示为加热后的变形.,温度分布特征是：在布管区，管板的大部分厚度上温度接近流入或流出管板处换热器的管程流体温度，只在靠近壳程流体一侧很薄的区域温度接近壳程流体温度。,温度分布,热应力分析,工况7 正常工况下应力强度分布云图,工况的产生： 设备正常运转时，为此种工况。此工况是壳程压力与温差应力强度共同作用的操作工况与管程压力单独作用是的工况进行加和。是管程、壳程、温差应力强度共同作用的结果。,弹性力学就是研究物体在外部因素作用下产生的应力、应变及其位移规律的一门学科。 理想弹性体：指符合下述四个假定的物体。 1）连续性假定。 2）完全弹性假定。 3）均匀性假定。 4）各向同性假定。,1.5 弹性力
9、学概念,求解思路：弹性力学的基本方程是以偏微分方程组来表示的。在分析应力时，一般总是从构件的连续性出发，依据无限小单元的物理数学模型建立微分方程式，然后求得解答。工程实际中，构件的形状、受力状态、边界条件都比较复杂，无法求出解析解，只能用数值计算求近似解。,有限单元法一个简单的例子,对于单元1：,轴向力公式：,分析过程: 结构离散化, 确定位移模式, 位移模式：近似的表示单元位移分量随坐标变化的简单函数。 单元力学特性分析 (1)通过几何方程建立单元应变与节点位移关系式； （2）利用物理方程导出单元应力与节点位移关系式； （3）由虚功原理推出作用于单元上节点力与节点位移之间的关系式，及单元的刚
10、度方程。,整体分析, 单元到整体的组集： (1) 所有相邻单元在公共节点处位移相等； （2）所有各节点需满足平衡条件。 解方程,输出计算结果, 其他处理,第二章 弹性力学基础,2.1 弹性力学中的几个基本概念,1-2 弹性力学中的几个基本概念,基本概念：,外力、应力、主应力。,1. 外力,体力、面力,(1) 体力, 弹性体内单位体积上所受的外力, 体力分布集度,（矢量）,X、Y、Z为体力矢量在坐标轴上的投影,单位：,N/m3,kN/m3,(2) 面力, 作用于物体表面单位面积上的外力, 面力分布集度（矢量）, 面力矢量在坐标轴上投影,单位：,1N/m2 =1Pa (帕),1MN/m2 = 10
11、6Pa = 1MPa (兆帕),(1) 一点应力的概念,内力,由于外力作用引起的物体内部相互作用力.,P,应力：由外力引起的在 P点的某一面上内力分布集度,应力分量,应力的法向分量, 正应力,应力的切向分量, 剪应力,单位:,与面力相同,MPa (兆帕),应力关于坐标连续分布的,2. 应力,用矩阵表示：,剪应力互等定理,应力符号的意义：,第1个下标 x 表示作用面；,第2个下标 y 表示作用方向.,应力正负号的规定：应力作用面的外法线方向与坐标轴正方向一致，则应力分量沿坐标轴正方向为正，沿坐标轴负方向为负；反之亦反。,2)主应力： 在过一点的所有截面中，存在着三个互相垂直的特殊截面，在这三个截
12、面上没有剪应力，只有正应力。这种没有剪应力存在的截面称为过该点的主平面，主平面上的正应力称为该点的主应力。主应力的方向总是与主平面的法线方向平行，称为该点的主方向。,弹性力学分析问题方法,静力学：平衡微分方程 几何学：几何方程 物理学：物理方程,2-2 平衡微分方程,取微元体PABC（P点附近），,Z 方向取单位长度。,设P点应力已知：,体力：X ，Y,两边同除以dx dy，并整理得：,两边同除以dx dy，并整理得：,平面问题的平衡微分方程：,（2-2）,空间平衡微分方程,应力分量与体积力分量之间的关系式,研究弹性体内各点的应变分量与位移分量之间的关系,2-3 几何方程,2-4 变形谐调方程
13、,描述六个应变分量之间所存在的关系式,2-5 物理方程,建立：平面问题中应力与应变的关系,物理方程也称：本构方程、本构关系、物性方程。,在完全弹性和各向同性的情况下，物性方程即为材料力学中的广义虎克（Hooke）定律。,其中：E为拉压弹性模量；G为剪切弹性模量；为侧向收缩系数，又称泊松比。,边界条件：,建立应力分量与外载荷间的关系。,只有当应力分量与外载荷满足一定条件时，物体边界上的点才能处于平衡状态。,边界分类,（1）位移边界,（2）应力边界,（3）混合边界,2-6 边界条件,指在物体的全部边界上，各点的位移分量均为坐标的已知函数,指在物体的全部边界上，各点所受的面力均为坐标的已知函数,指在
14、物体的部分边界上位移分量已知，而在其余的边界上面力分量是已知的,弹性力学基本方程,一、平衡方程,二、几何方程,三、物理方程,四、变形协调方程,五、边界条件(应力,位移),位移,应力,弹性力学基本方程矩阵表示,位移列阵,体积力列阵,应力列阵,应变列阵,表面外法线方向余弦矩阵,微分算子列阵,表面力列阵,已知位移列阵,1、平面应力问题,2.7 弹性力学中的几个典型问题,2.7.1 平面问题,所谓平面应力问题，是指所研究的对象在z方向上的尺寸很小(即呈平板状)，外载荷(包括体积力)都与z轴垂直，且沿z方向没有变化，在zh/2处的两个外表面(平面)上不受任何载荷。,2、平面应变问题,当物体z方向上的尺寸
15、很长，物体所受的载荷(包括体积力)又垂直于z铀且不沿长度方向(z方向)变化，即物体的内在因素和外来作用都不沿长度方向变化，那么这种问题就称为平面应变问题。,平面应变问题,（1）平衡方程,（2）几何方程,（3)物理方程,2.7.2 轴对称问题,在空间问题中，如果弹性体的几何形状、约束状态，以及其他外在因素都是对称于某一根轴(过该轴的任一平面都是对称面)，那么弹性体的所有应力、应变和位移也就都对称于这根轴。这类问题通常称为空间轴对称问题。,2.8 圣维南原理,圣维南原理一般可以这样来叙述：如果把物体的一小部分边界上的面力，变换为分布不同但静力等效的面力(即主矢量相同、对同一点的主矩也相同)，那么，
16、近处的应力分布将有显著的改变，但远处所受的影响可以不计。,作用：可以简化边界条件,圣维南原理也可以这样来叙述：如果物体一小部分边界上的面力是一个平衡力系(主矢量及主矩都等于零)，那么，这个面力就只会使得近处产生显著的应力，远处的应力可以不计。,2.9弹性力学问题的求解方法,（1）按位移求解（位移法、刚度法）,以u、v 为基本未知函数，将平衡方程和边界条件都用u、v 表示，并求出u、v ，再由几何方程、物理方程求出应变分量与应力。,（2）按应力求解（力法，柔度法）,以应力分量 为基本未知函数，将所有方程都用应力分量表示，并求出应力分量 ，再由物理方程、几何方程求出应变分量与位移。,（3）混合求解,以部
