欢迎来到三一办公! | 帮助中心 三一办公31ppt.com(应用文档模板下载平台)
三一办公
全部分类
  • 办公文档>
  • PPT模板>
  • 建筑/施工/环境>
  • 毕业设计>
  • 工程图纸>
  • 教育教学>
  • 素材源码>
  • 生活休闲>
  • 临时分类>
  • ImageVerifierCode 换一换
    首页 三一办公 > 资源分类 > PPT文档下载  

    有限元法及应用知识点超全总结课件.ppt

    • 资源ID:3762171       资源大小:1.78MB        全文页数:84页
    • 资源格式: PPT        下载积分:16金币
    快捷下载 游客一键下载
    会员登录下载
    三方登录下载: 微信开放平台登录 QQ登录  
    下载资源需要16金币
    邮箱/手机:
    温馨提示:
    用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)
    支付方式: 支付宝    微信支付   
    验证码:   换一换

    加入VIP免费专享
     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

    有限元法及应用知识点超全总结课件.ppt

    有限元法及应用总结,串讲,1.,有限元的作用是什么?,?,1,)减少模型试验的数量;,?,计算机模拟允许对大量的假设情况进行快速而有效的,试验。,?,?,?,?,2,)模拟不适合在原型上试验的设计;例如:器,官移植,比如人造膝盖。,3,)节省费用,降低设计与制造、开发的成本;,4,)节省时间,缩短产品开发时间和周期;,5,)创造出更可靠、高品质的设计。,2.,有限元的基本概念,用有限个单元来描述。,有限元法:把求解区域看作由许多小的在节点,处相互连接的单元(子域)所构成,其模型给,出基本方程的分片(子域)近似解,由于单元,(子域)可以被分割成各种形状和大小不同的,尺寸,所以它能很好地适应复杂的几何形状、,复杂的材料特性和复杂的边界条件。,再加上它有成熟的大型软件系统支持,使其已,成为一种非常受欢迎的、应用极广的数值计算,方法。,?,有限元:通俗的讲就是对一个真实的系统,?,?,有限元模型与有限元分析,?,?,有限元模型:它是真实系统理想化的数学,抽象。由一些简单形状的单元组成,单元,之间通过节点连接,并承受一定载荷。,有限元分析:是利用数学近似的方法对真,实物理系统(几何和载荷工况)进行模拟。,并利用简单而又相互作用的元素,即单元,,就可以用有限数量的未知量去逼近无限未,知量的真实系统。,3.,有限单元法的特点有哪些?,?,?,?,?,?,1,)把连续体划分成有限个单元,把单元的交界结点(节点),作为离散点;,2,)不考虑微分方程,而从单元本身特点进行研究。,3,)理论基础简明,物理概念清晰,且可在不同的水平上建,立起对该法的理解。,4,)具有灵活性和适用性,适应性强。(它可以把形状不同、,性质不同的单元组集起来求解,故特别适用于求解由不同,构件组合的结构,应用范围极为广泛。它不仅能成功地处,理如应力分析中的非均匀材料、各向异性材料、非线性应,力、应变以及复杂的边界条件等问题,且随着其理论基础,和方法的逐步完善,还能成功地用来求解如热传导、流体,力学及电磁场领域的许多问题。),5,)在具体推导运算过程中,广泛采用了矩阵方法。,4.,有限元法涉及的内容有哪些?,?,?,?,?,?,?,有限元法在数学和力学领域所依据的理论;,单元的划分原则;,形状函数的选取及协调性;,有限元法所涉及的各种数值计算方法及其误,差、收敛性和稳定性;,计算机程序设计技术;,向其他各领域的推广。,5.,有限元法的分类,?,有限元法可以分为两类,即线弹性有限元,法和非线性有限元法。其中线弹性有限元,法是非线性有限元法的基础,二者不但在,分析方法和研究步骤上有类似之处,而且,后者常常要引用前者的某些结果。,线弹性有限元,?,?,线弹性有限元是以理想弹性体为研究对象的,,所考虑的变形建立在小变形假设的基础上。在,这类问题中,材料的应力与应变呈线性关系,,满足广义胡克定律;应力与应变也是线性关系,,线弹性问题可归结为求解线性方程问题,所以,只需要较少的计算时间。如果采用高效的代数,方程组求解方法,也有助于降低有限元分析的,时间。,线弹性有限元一般包括线弹性静力学分析与线,弹性动力学分析两方面。,非线性有限元,?,?,?,?,?,非线性问题与线弹性问题的区别:,1,)非线性问题的方程是非线性的,一般需要迭代,求解;,2,)非线性问题不能采用叠加原理;,3,)非线性问题不总有一致解,有时甚至没有解。,以上三方面的因素使得非线性问题的求解过程比,线弹性问题更加复杂、费用更高和更具有不可预,知性。,1,)材料非线性问题,?,?,?,?,有限元求解非线性问题可分为以下三类:,1,)材料非线性问题,材料的应力和应变是非线性的,但应力与应变却很微小,,此时应变与位移呈线性关系,这类问题属于材料的非线性,问题。由于从理论上还不能提供能普遍接受的本构关系,,所以,一般材料的应力与应变之间的非线性关系要基于试,验数据,有时非线性材料特性可用数学模型进行模拟,尽,管这些模型总有他们的局限性。,在工程实际中较为重要的材料非线性问题有:非线性弹性,(包括分段线弹性)、弹塑性、粘塑性及蠕变等。,2,)几何非线性问题,?,?,几何非线性问题是由于位移之间存在非线,性关系引起的。,当物体的位移较大时,应变与位移的关系,是非线性关系。研究这类问题一般都是假,定材料的应力和应变呈线性关系。它包括,大位移大应变及大位移小应变问题。如结,构的弹性屈曲问题属于大位移小应变问题,,橡胶部件形成过程为大应变问题。,3,)非线性边界(接触问题),?,?,?,在加工、密封、撞击等问题中,接触和摩擦的,作用不可忽视,接触边界属于高度非线性边界。,平时遇到的一些接触问题,如齿轮传动、冲压,成型、轧制成型、橡胶减振器、紧配合装配等,,当一个结构与另一个结构或外部边界相接触时,通常要考虑非线性边界条件。,实际的非线性可能同时出现上述两种或三种非,线性问题。,*6.,有限元的基础理论包括哪几部分?,?,?,?,?,1.,加权余量法,加权余量法:是指采用使余量的加权函数为零,求得微分方程近似解的方法称为加权余量法。,(,Weighted residual method WRM,),加权余量法是求解微分方程近似解的一种有效,的方法。,显然,任何独立的完全函数集都可以作为权函,数。按照对权函数的不同选择得到不同的加权,余量计算方法,主要有:配点法、子域法、最,小二乘法、力矩法和伽辽金法。其中伽辽金法,的精度最高。,2.,里兹方法,?,里兹方法:如果微分方程具有线性和自伴随的,性质,那么它不仅可以建立它的等效积分形式,,并利用加权余量法求其近似解,而且还可以建,立与之相等效的变分原理,从而得到的另一种,近似求解方法。,?,自然变分原理:原问题的微分方程和边界条件的等效,积分的伽辽金法等效于它的变分原理,即原问题的微,分方程和边界条件等效于泛函的变分为零,亦即泛函,取驻值。反之,如果泛函取驻值则等效于满足问题的,微分方程和边界条件。而泛函可以通过原问题的等效,积分的伽辽金法而得到,我们称这样得到的变分原理,为自然变分原理。,2.,里兹方法(续),?,对于具有线性、自伴随性质的微分方程在得到,与它相等效的变分原理以后,可以用来建立求,近似解,这一过程即里兹方法。它的实质是从,一族假定解中寻求满足泛函变分的“最好的”,解。显然,近似解的精度与试探函数(形函数,或试函数)的选择有关,如果知道所求解的一,般性质,那么可以通过选择反映此性质的试探,函数来改进近似解,提高近似解的精度。,3.,虚功原理,平衡方程和几何方程的等效积分“弱”形,式,?,虚功原理包含虚位移原理和虚应力原理,是虚位,移原理和虚应力原理的总称。他们都可以认为是,与某些控制方程相等效的积分“弱”形式。虚功,原理:变形体中任意满足平衡的力系在任意满足,协调条件的变形状态上作的虚功等于零,即体系,外力的虚功与内力的虚功之和等于零。,?,虚位移原理是平衡方程和力的边界条件的等效积,分的“弱”形式;,?,虚应力原理是几何方程和位移边界条件的等效积,分“弱”形式。,3.,虚功原理(续),平衡方程和几何方程的等效积分“弱”,形式,?,虚位移原理的力学意义:如果力系是平衡的,,则它们在虚位移和虚应变上所作的功的总和为,零。反之,如果力系在虚位移(及虚应变)上,所作的功的和等于零,则它们一定满足平衡方,程。所以,虚位移原理表述了力系平衡的必要,而充分条件。,?,一般而言,虚位移原理不仅可以适用于线弹性,问题,而且可以用于非线性弹性及弹塑性等非,线性问题。,?,但是否适用所有的问题呢?,3.,虚功原理(续),平衡方程和几何方程的等效积分“弱”,形式,?,虚应力原理的力学意义:如果位移是协调的,则虚应,力和虚边界约束反力在他们上面所作的功的总和为零。,反之,如果上述虚力系在他们上面所作的功的和为零,,则它们一定是满足协调的。所以,虚应力原理表述了,位移协调的必要而充分条件。,?,虚应力原理可以应用于线弹性以及非线性弹性等不同,的力学问题。,?,但是必须指出,无论是虚位移原理还是虚应力原理,,他们所依赖的几何方程和平衡方程都是基于小变形理,论的,他们不能直接应用于基于大变形理论的力学问,题。,4.,最小位能原理和最小余能原理,?,明确:最小位能原理是建立在虚位移原理基础上,的,而最小余能原理建立在虚应力原理基础上。,?,最小位能原理是指在所有可能位移中,真实位移,使系统总位能取最小值。,?,总位能是指弹性体变形位能和外力位能之和。,?,最小余能原理是指在所有的应力中,真实应力使,系统的总余能取最小值。,?,总余能是指弹性体余能和外力余能总和。,4.,最小位能原理和最小余能原理(续),?,一般而言,利用最小位能原理求得位移近似解,的弹性变形能是精确解变形能的下界,即近似,的位移场在总体上偏小,也就是说结构的计算,模型显得偏于刚硬;而利用最小余能原理求得,的应力近似解的弹性余能是精确解余能的上界,,即近似的应力解在总体上偏大,结构的计算模,型偏于柔软。,?,当分别利用这两个极值原理求解同一问题时,,我们将获得这个问题的上界和下界,可以较准,确地估计所得近似解的误差,这对工程计算具,有实际意义。,*7.,单元划分原则是什么,?,?,梁、杆单元划分的原则,?,两个节点之间的杆构成一个单元,节点可按以,下原则划分:,1,)杆件的交点一定要选为节点,(,梯子);,2,)阶梯形杆截面变化处一定取为节,点(阶梯轴);,3,)支撑点与自由端要选为节,点(悬臂梁);,4,)集中载荷作用处最好选为,节点;,5,)欲求位移的点要选为节点;,6,)单元,长度最好基本相同。,平面单元划分原则,?,1.,单元形状:常用单元形状有三角形单元、矩形单元和等,参数单元。他们的特点是单元的节点数越多,其计算精,度越高,三角形单元与等参数单元可适应任意边界。,?,2.,划分原则:,?,1,)划分单元的个数,视计算机要求的精度和计算机容量,而定,单元分得越多,块越小其精度越高,但需要的计,算机容量越大,因此,须根据实际情况而定。,?,2,)划分单元的大小,可根据部位不同有所不同,在位,移或应力变化大的部位取得单元要小;在位移或应力变,化小的部位取得单元要大,在边界比较平滑的部位,单,元可大。,平面单元划分原则(续),?,3,)划分单元的形状,一般均可取成三角形或,等参元。对于平直边界可取成矩形单元,有时,也可以将不同单元混合使用,但要注意,必须,节点与节点相连,切莫将节点与单元的边相连。,4,)单元各边的长不要相差太大,否则将影响,求解精度。,?,5,)尽量把集中力或集中力偶的作用点选为节,点。,?,6,)尽量利用对称性,以减少计算量(有限元,法的最大优点在于使用了矩阵的方法)。,*8.,有限元法分析过程,?,有限元法分析过程大体可分为:前处理、分析、,后处理三大步骤。,?,对实际的连续体经过离散化后就建立了有限元,分析模型,这一过程是有限元的前处理过程。,在这一阶段,要构造计算对象的几何模型,要,划分有限元网格,要生成有限元分析的输入数,据,这一步是有限元分析的关键。,*8.,有限元法分析过程,(续),?,有限元分析过程主要包括:单元分析、整体分析、,载荷移置、引入约束、求解约束方程等过程。这一,过程是有限元分析的核心部分,有限元理论主要体,现在这一过程中。,?,有限元法包括三类:有限元位移法、有限元力法、,有限元混合法。,?,在有限元位移法中,选节点位移作为基本未知量;,?,在有限元力法中,选节点力作为未知量;,?,在有限元混合法中,选一部分基本未知量为节点位,移,另一部分基本未知量为节点力。,*8.,有限元法分析过程,(续),?,有限元位移法计算过程的系统性、规律性强,特,别适宜于编程求解。一般除板壳问题的有限元应,用一定量的混合法外,其余全部采用有限元位移,法。因此,一般不做特别声明,有限元法指的是,有限元位移法。,?,有限元分析的后处理主要包括对计算结果的加工,处理、编辑组织和图形表示三个方面。它可以把,有限元分析得到的数据,进一步转换为设计人员,直接需要的信息,如应力分布状态、结构变形状,态等,并且绘成直观的图形,从而帮助设计人员,迅速的评价和校核设计方案。,9.,有限元法的收敛性概念与收敛条件,?,有限元法是一种数值分析方法,因此应,考虑收敛性问题。,?,有限元法的收敛性是指:当网格逐渐加,密时,有限元解答的序列收敛到精确解;,或者当单元尺寸固定时,每个单元的自,由度数越多,有限元的解答就越趋近于,精确解。,9.,有限元法的收敛性概念与收敛条件,(续),?,有限元的收敛条件包括如下四个方面:,?,1,)单元内,位移函数必须连续。多项式是单值连续函,数,因此选择多项式作为位移函数,在单元内的连续,性能够保证。,?,2,)在单元内,位移函数必须包括常应变项。每个单元,的应变状态总可以分解为不依赖于单元内各点位置的,常应变和由各点位置决定的变量应变。当单元的尺寸,足够小时,单元中各点的应变趋于相等,单元的变形,比较均匀,因而常应变就成为应变的主要部分。为反,映单元的应变状态,单元位移函数必须包括常应变项。,9.,有限元法的收敛性概念与收敛条件,(续),?,3,)在单元内,位移函数必须包括刚体位移项。一般情,况下,单元内任一点的位移包括形变位移和刚体位移,两部分。形变位移与物体形状及体积的改变相联系,,因而产生应变;刚体位移只改变物体位置,不改变物,体的形状和体积,即刚体位移是不产生变形的位移。,空间一个物体包括三个平动位移和三个转动位移,共,有六个刚体位移分量。,?,由于一个单元牵连在另一些单元上,其他单元发生变,形时必将带动单元做刚体位移,由此可见,为模拟一,个单元的真实位移,假定的单元位移函数必须包括刚,体位移项。,9.,有限元法的收敛性概念与收敛条件,(续),?,4,)位移函数在相邻单元的公共边界上必须协,调。对一般单元而言,协调性是指相邻单元在,公共节点处有相同的位移,而且沿单元边界也,有相同的位移,也就是说,要保证不发生单元,的相互脱离开裂和相互侵入重叠。要做到这一,点,就要求函数在公共边界上能由公共节点的,函数值唯一确定。对一般单元,协调性保证了,相邻单元边界位移的连续性。,?,但是,在板壳的相邻单元之间,还要求位移的,一阶导数连续,只有这样,才能保证结构的应,变能是有界量。,9.,有限元法的收敛性概念与收敛条件,(续),?,总的说来,协调性是指在相邻单元的公共边界,上满足连续性条件。,?,前三条又叫完备性条件,满足完备条件的单元,叫完备单元;第四条是协调性要求,满足协调,性的单元叫协调单元;否则称为非协调单元。,完备性要求是收敛的必要条件,四条全部满足,,构成收敛的充分必要条件。,?,在实际应用中,要使选择的位移函数全部满足,完备性和协调性要求是比较困难的,在某些情,况下可以放松对协调性的要求。,9.,有限元法的收敛性概念与收敛条件,(续),?,需要指出的是,有时非协调单元比与它对应的协调单元,还要好,其原因在于近似解的性质。假定位移函数就相,当于给单元施加了约束条件,使单元变形服从所加约束,,这样的替代结构比真实结构更刚一些。但是,这种近似,结构由于允许单元分离、重叠,使单元的刚度变软了,,或者形成了(例如板单元在单元之间的绕度连续,而转,角不连续时,刚节点变为铰接点)对于非协调单元,上,述两种影响有误差相消的可能,因此利用非协调单元有,时也会得到很好的结果。在工程实践中,非协调元必须,通过“小片试验后”才能使用。,10.,应力的单元平均或节点平均处理方法,?,?,最简单的处理应力结果的方法是取相邻单元或,围绕节点各单元应力的平均值。,?,1.,取相邻单元应力的平均值,?,这种方法最常用于,3,节点三角形单元中。这种,最简单而又相当实用的单元得到的应力解在单,元内是常数。可以将其看作是单元内应力的平,均值,或是单元形心处的应力。由于应力近似,解总是在精确解上下振荡,可以取相邻单元应,力的平均值作为此两个单元合成的较大四边形,单元形心处的应力。,10.,应力的单元平均或节点平均处理方法,?,(续),?,如,2,单元的情况下,取平均应力可以采用算术平,均,,?,即平均应力,=,(单元,1,的应力,+,单元,2,的应力),/2,。,?,也可以采用精确一些的面积加权平均,,?,即平均应力,=,单元,1,应力,单元,1,的面积,+,单元,2,应力,单元,2,面积,/,(单元,1,面积,+,单元,2,面积),?,当相邻两单元面积相差不大时,两者的结果基本,相同。在单元划分时应避免相邻两单元的面积相,差太多,从而使求解的误差相近。,10.,应力的单元平均或节点平均处理方法,?,(续),?,一般而言,,3,节点三角形单元的最佳应力点是,单元的中心点,此点的应力具有,1,阶的精度。,?,2.,取围绕节点各单元应力的平均值,?,首先计算围绕该节点(,i,)周围的相关单元在该,节点出的应力值,?,,然后以他们的平均值作为,该节点的最后应力值,?,,即,1,?,?,?,?,m,e,i,i,m,i,e,i,e,?,1,?,其中,,1m,是围绕在,i,节点周围的全部单元。取,平均值时也可进行面积加权。,*11.,有限元法求解问题的基本步骤,1.,结构离散化,对整个结构进行离散化,将其分割成若干,个单元,单元间彼此通过节点相连;,2.,求出各单元的刚度矩阵,?,K,?,(,e,),是由单元节点位移量,?,?,?,求单元,节点力向量,?,F,?,的转,移矩阵,其关系,式为,?,F,?,(,e,),?,k,(,e,),?,?,?,(,e,),(,1,?,1,),(,e,),?,K,?,(,e,),(,e,),3.,集成总体刚度矩阵,K,并写出总体平衡方程:,?,总体刚度矩阵,K,是由整体节点位移向量,?,?,?,求整体节点力向量,?,F,?,的转移矩阵,其关,系式为,?,F,?,?,K,?,?,?,,此即为总体平衡,方程。,?,4.,引入支撑条件,求出各节点的位移,?,节点的支撑条件有两种:一种是节点,n,沿,某个方向的位移为零,另一种是节点,n,沿,某个方向的位移为一给定值。,?,5.,求出各单元内的应力和应变,。,12.,单元刚度矩阵的特性,?,单元刚度矩阵无论在局部坐标系中还是在整,体坐标系中都具有相同的三个特性:,?,1,),对称性,?,由材料力学中的位移互等定理可知,对一个,构件,作用在点,j,的力引起点,i,的绕度等于有,同样大小而作用于点,i,的力引起的点,j,的绕度,,(,e,),(,e,),k,?,k,ij,ij,即,,表明单元刚度矩阵是一个对称,矩阵。,12.,单元刚度矩阵的特性(续),?,2,),奇异性,?,无逆阵的矩阵就叫做奇异矩阵,其行列,式的值为,0,,即,?,k,?,?,0,,这一点可以从例,题直接得到验证。其物理意义是引入支,撑条件之前,单元可平移。,(,e,),12.,单元刚度矩阵的特性(续),?,3,),分块性,(,e,),?,有前面所讲的内容可以看出,矩阵,?,k,?,可以,用虚线分成四块,因此可写成如下的分块形,式,,?,?,f,?,1,?,?,?,?,?,f,?,2,?,(,e,),mn,(,e,),?,?,k,?,11,?,?,?,?,k,?,21,?,k,?,12,?,?,k,?,22,?,?,(,e,),?,?,?,?,1,?,?,?,?,?,?,?,2,?,(,e,),?,式中,?,k,?,局部坐标系中单元,(e),按局部码,标记的节点,m,、,n,之间的刚度子矩阵。,*13.,求图中所示刚架中各单元在整体坐标系中,的单元刚度矩阵,(1),。,?,设两杆的长度与截面尺寸彼此相等。,(空心杆)其中,,7,2,L=200cm,D=5cm,d=4cm,E=2,10,N/cm,。,*13,:求图所示结构的节点位移向量。,?,已知节点,1,处承受外载,M,条件同前例。,1,?,5,?,10,Ncm,5,,其余,*14.,刚架结构中非节点载荷的处理的方法,?,在刚架结构以及其他较复杂的结构上,他们所受,的载荷可以直接作用在节点上,又可以不直接作,用在节点上而作用于单元节点间的其他位置上。,后一种情况下的载荷称为非节点载荷。有限元分,析时,总体刚度方程中所用到的力向量,?,F,?,是节点,力向量。因此在进行整体分析前应当进行载荷的,移植,将作用于单元上的力移植到节点上。移植,时按静力等效的原则进行。,*14.,刚架结构中非节点载荷的处理的方法,(续),?,处理非节点载荷一般可直接在整体坐标,系内进行,其过程为:,?,1,)将各杆单元看成一根两端固定的梁,,分别求出两个固定端的约束反力。其结,果可直接利用材料力学的公式求得;,?,2,)将各固定端的约束反力变号,按节点,进行集成,获得各节点的等效载荷,*15.,总体刚度矩阵的集成法,?,使用刚度矩阵获得的方法获得总体刚度矩,阵。在此将其扩展到由整体坐标系中的单,元刚度矩阵的子矩阵集成总体刚度矩阵。,步骤如下:,?,1,)对一个有,n,个节点的结构,将总体刚度,矩阵,K,划分为,n,n,各子区间,然后按节点,总码的顺序进行编号;,*15.,总体刚度矩阵的集成法(续),?,2,)将整体坐标系中单元刚度矩阵的各子,矩阵根据其下标的两个总码对号入座,,写在总体刚度矩阵相应的子区间;,?,3,)同一子区间内的子矩阵相加,成为总,体刚度矩阵中的相应的子矩阵。,16.,总体刚度矩阵的特性,?,1,)对称性:因为由此特性,在计算机中只需存储其,上三角部分;,?,2,)奇异性:物理意义仍为在无约束的情况下,整个,结构可做刚体运动;,?,3,)稀疏性:,K,中有许多零子矩阵,而且在非零子矩,阵中还有大量的零元素,这种矩阵称为稀疏矩阵。大,型结构的总体刚度矩阵一般都是稀疏矩阵;,?,4,)分块性:这个性质已经利用过,在此不再叙述。,?,了解刚度矩阵的这些特性非常有用,它可以大大减少,计算机的内存和计算工作量。,*17.,平面问题离散化时的规定,?,?,?,?,?,1,)单元之间只在节点处相连;,2,)所有的节点都为铰接点;,3,)单元之间的力通过节点传递;,4,)外载荷都要移植到节点上;,5,)在节点位移或某一分量可以不计之处,就必须,在该节点安置一个铰支座或相应的连杆支座。,?,通过以上的规定来建立平面有限元分析模型。,*18.,平面离散化的一些定性的规律,?,?,?,?,?,?,?,?,1,)结构对称性的利用,2,)对称结构的网格布局,3,)划分网格要兼顾精度和经济性,4,)不连续出的自然分割,5,)几何形状的近似与过渡圆角的处理,6,)单元形态的选择,7,)边界条件的确定,8,)单元和节点编号,19.,结构对称性的利用规律,?,一般来说,作用在对称结构上的载荷系统分为对称的、,反对称的和一般的三种情况。,?,1.,结构对称,载荷对称或反对称,?,这种情况下,对称面上的边界条件可按一下规则确定:,?,A.,在不同的对称面上,将位移分量区分为对称分量和,反对称分量;,?,B.,将载荷也按不同的对称面分别区分为对称分量和反,对称分量;,?,C.,对于同一个对称面,如载荷是对称的,则对称面上,位移的反对称分量为零,如载荷是反对称的,则对称,面上位移的对称分量为零。,结构对称,载荷一般的情况,?,如果所分析的结构对称,但载荷是不,对称的,也不是反对称的,这时可以,将这种结构系统简化成载荷为对称和,/,或反对称情况的组合,仍可以简化分,析过程,提高分析的综合效率。,?,如图,a,所示,结构对称,载荷一般,可,将其载荷分解为图,b,和图,c,的组合。图,b,为对称结构,载荷对,x,、,y,轴均为对称,,图,c,为结构对称,载荷对,x,轴反对称、,对,y,轴对称,此时可取相同的四分之一,进行研究,分别施加对称面上节点的,边界条件,进行两次分析计算,并将,计算结果迭加起来,即可得到原结构,四分之一的解答,进而得出整个结构,的解答。,对称性利用中的特殊问题,?,利用结构的对称性取某一部分建立有限元模型时,往,往会产生约束不足现象。,?,例如,若取上例中图,c,的四分之一建立有限元时,根据,上述分析,在两对称面上应加水平放置的滚动铰支座,,因此模型在垂直方向存在刚体位移。对这种约束不足,问题,利用有限元分析时,必须增加附加约束,以消,除模型的刚体位移。在本例中,垂直方向可以用刚度,很小的杆单元或边界弹簧单元连接到模型某节点上,,使得既消除了模型的刚体位移,又不致于因附加的杆,单元或边界弹簧单元刚度太大而影响结构原有的变形,状态。,20.,单元形态的选择原则,?,单元形态包括单元形状、边中节点的位置、细长比等,,在结构离散化过程中必须合理选择。一般来说,为了,保证有限元分析的精度,必须是单元的形态尽可能的,规则。,?,对于三角形单元,三条边长尽量接近,不应出现大的,钝角、大的边长。这是因为根据误差分析,应力和位,移的误差都和单元的最小内角的正弦成反比。因而,,等边三角形单元的形态最好,它与等腰直角三角形单,元的误差之比为,sin45,:sin60,=1:1.23,。但是为了适应,弹性体边界,以及单元由小到大逐渐过渡,不可能是,所有的三角形单元都接近等边三角形。实际上,常常,使用等腰直角三角形。,20.,单元形态的选择原则(续),?,对于矩形单元来说,细长比不宜过大。细长比是指单,元最大尺寸和最小尺寸之比。最优细长比在很大程度,上取决于不同方向上位移梯度的差别。梯度较大的方,向,单元尺寸要小些,梯度小的方向,单元尺寸可以,大一些;如果各方向上位移梯度大致相同,则细长比,越接近,1,,精度越高。有文献推荐,一般情况下,为了,得到较好的位移结果,细长比不应超过,7,;为了获得较,好的应力结果,细长比不应超过,3,。一般情况下,正方,形单元的形态最好。,?,对于一般的四边形单元应避免过大的边长比,过大的,边长比会导致病态的方程组。,*21.,边界条件的确定,?,确定边界条件是建立有限元模型的重要一环,合理确,定有限元模型的边界条件是成功地进行结构有限元分,析的基本要求。,?,一般情况下,建模对象的边界条件是明确的。根据力,学模型的边界条件可以很容易确定其有限元模型的边,界条件。例如电线杆插入地基的一端为固定端,桥梁,一端为固定铰支座,另一端为滚动较支座。,?,但是,在机械工程中,建模对象往往是整个结构中的,一部分,在建立有限元模型,确定其边界条件时,必,须考虑其余部分的影响。这方面主要考虑如下两类问,题。,?,1.,边界位置的确定,?,在建立连续弹性体局部区域的有限元模型时,往,往取该局部区域为隔离体,取其隔离边界条件为,零位移约束,并通过试探校正确定零位移边界条,件的位置。例如,进行齿轮齿有限元分析时,取,一个轮齿的局部区域为隔离体,如图所示,设定,PQRS,的边界条件为零位移约束,通过改变边界,深度,PQ,和边界宽度,PS,研究边界位置对齿根最大,拉应力的影响,最后确定合理的边界条件。,?,2.,边界条件的确定,?,有些分析对象的边界位置是零部件的连接部位。,在建立有限元模型时,必须研究如何给定边界,位置上的边界条件,以反映相连接结构的影响。,确定这种问题的边界条件是用简单支撑连杆替,代相连接结构的作用,使替代后结构的系统刚,度等价于原结构的系统刚度。如分析机床主轴,和传动轴时,可以利用等刚度的杆单元替代轴,承和支座的作用,使轴的分析中包含有轴承和,支座的影响。,*22.,单元和节点编号规则,?,当利用整体刚度矩阵的带状特征进行存贮和求,解方程组时,单元节点编号直接影响系统整体,刚度矩阵的半带宽,也就是影响在计算机中存,贮信息的多少、计算时间和计算费用。因而,,要求合理的节点编号使带宽极小化。半带宽的,计算公式:,?,半带宽,d=,(单元节点号的最大差值,+1,)节点,自由度,?,由此,进行网格节点编号时应使网格中单元节,点号的最大差值最小,这样才能保证半带宽最,小。试比较下图。,*22.,单元和节点编号规则(续),?,图所示网格的四种编号,方案中,单元节点标号,的最大差值分别为,5,,,3,,,5,,,9,。显然,图,2,方案要,合理。由此得出结论:,沿着短边方向按列,-,列,-,列,-,列地顺序编号比沿着长,度方向按行,-,行,-,行,-,行地,顺序要合理(半带宽小),*22.,单元和节点编号规则(续),?,然而,对于具有中间节点的单元或空间问题,,须借助于带宽极小化的优化程序来对节点重新,编号,先进的有限元程序包一般都配备有这样,的程序。,?,对单元的编号只影响整体刚度矩阵的装配时间。,由于这一时间在有限元运算时间中只占很小的,比例,因而对于单元的编号并无特殊的要求。,*23.,掌握分析三角形单元的位移模,式求解方法,?,如图所示,在局部坐标系中,三角形平,面单元的三个节点分别为,1,、,2,、,3,,其,编号按逆时针方向进行,节点坐标分别,?,x,y,?,、,?,x,y,?,为,?,x,y,?,、,1,1,2,2,3,3,24.,求解平面问题中局部坐标,?,?,系中的单元刚度矩阵,?,k,?,e,?,将几何方程和弹性方程代入虚功方程经,整理后得:局部坐标系中,?,f,?,?,e,?,?,?,k,?,?,?,?,?,e,?,?,e,?,?,e,?,?,?,k,?,式中,单元刚度矩阵,?,k,?,?,e,?,?,?,?,B,?,?,D,?,B,?,tdxdy,?,?,B,?,?,D,?,B,?,?,tdxdy,?,?,B,?,?,D,?,B,?,t,?,T,T,T,?,其中,t,三角形单元平板的厚,度,,?,三角形单元的面积,25.,平面问题中非节点载荷转换为等效,节点载荷,?,由于三角形单元复杂的力学性质,不能像,分析刚架时那样简单地利用材料力学公式,来求解,而要用虚功方程将加在结构上的,非节点载荷转换为等效节点载荷。,?,掌握以下两种常见的非节点载荷的移植结,果。,1,)作用在单元一条侧边上的集中力,设,Q,平行于,x,方向,如图,4-14,所示,则,等效节点载荷为,l,i,F,ix,?,Q,F,jx,?,Q,F,kx,?,0,F,iy,?,0,F,jy,?,0,F,ky,?,0,l,l,l,j,若,Q,平行于,y,方向,结果与此相仿。,1,)作用在单元一条侧边上的集中力,2,)作用在单元一条侧边上呈三角形,分布的载荷,?,设载荷平行于,x,方向,如图,4-15,所示,则等效,节点载荷为,qlt,qlt,F,ix,?,F,jx,?,F,kx,?,0,F,iy,?,0,F,jy,?,0,F,ky,?,0,6,3,?,若分布载荷为集度是,q,的均布载荷,则,F,ix,?,F,jx,qlt,?,3,?,其余分量为零。,2,)作用在单元一条侧边上呈三角形,分布的载荷,*26.,例:,求例,4-7,图所示结构节点的位移量。,?,已知,Q,?,100,N,q,?,0,.,1,N,/,mm,L,1,?,40,mm,L,2,?,60,mm,a,?,100,mm,t,?,10,mm,E,?,2,?,10,N,/,mm,u,?,0,5,2,2,ANSYS,软件基本知识,?,?,?,?,?,?,?,1.ANSYS,图形用户界面(,GUI,)有哪几部分组成?,2.,比较对话框中的,:“OK”,与“,Apply”,的区别;,3.,熟悉单元类型的含义;,4.ANSYS,文件及工作文件名的含义,.,5.,应用,ANSYS,软件计算,如图所示的平面桁架,长度单位为,m,,,求支座反力和各杆内力。设弹性模量为,2e+11,,泊松比,0.3,,杆件,截面面积为,0.01m,2,.,6.,给定一个简单的物理现象,能够使用,ANSYS,创建一个,2D,的有,限元模型,.,7.,熟练运用将几何模型划分网格后,进行加载与求解及结果的后,处理,.,有关软件的几个问题的处理,1.,载荷与载荷分类,ANSYS,中的载荷可分为,:,自由度,DOF,-,定义节点的自由度(,DOF,),值,(,结构分析,_,位移、,热分析,_,温度、电磁分析,_,磁势等,),?,集中载荷,-,点载荷,(,结构分析,_,力、热分析,_,热导率、电磁分析,_ magnetic current segments),?,面载荷,-,作用在表面的分布载荷,(,结构分析,_,压力、热分析,_,热,对流、电磁分析,_magnetic Maxwell surfaces,等,),?,体积载荷,-,作用在体积或场域内,(,热分析,_,体积膨胀、内生成热,、电磁分析,_,magnetic current density,等,),?,惯性载荷,-,结构质量或惯性引起的载荷,(,重力、角速度等,),2.,添加载荷应遵循的原则,?,简化假定越少越好。,?,使施加的载荷与结构的实际承载状态保持吻合;,?,如果没法做得更好,只要其它位置结果正确也是可以认为,是正确的,但是你必须忽略“不合理”边界的附近一定区,域内的应力。,?,加载时,必须十分清楚各个载荷的施加对象及定义载荷。,?,除了对称边界外,实际上不存在真正的刚性边界。,?,不要忘记泊松效应;,?,添加刚体运动约束,但不能添加过多的(其它)约束。,?,实际上,集中载荷是不存在的;,?,轴对称模型具有一些独一无二的边界特性,。,3.,求解时模型是否准备就绪,?,在求解初始化前,应进行分析数据检查,包括下面,内容,:,?,1.,统一的单位,;2.,单元类型和选项,;3.,材料性质参数,:,考虑惯性时应输入材料密度,;,热应力分析时应输入,材料的热膨胀系数,;4.,实常数,(,单元特性,);5.,单元实,常数和材料类型的设置,;6.,实体模型的质量特性,(Preprocessor Operate Calc Geom Items);7.,模型中,不应存在的缝隙,;8.,壳单元的法向,;9.,节点坐,标系,;10.,集中、体积载荷面力方向,;11.,温度场的分,布和范围,;12.,热膨胀分析的参考温度。,有必要注意:,在求解过程中,应将,OUTPUT,窗口提到最前面。,ANSYS,求解过程中的一系列信息都将显示在此,窗口中,主要信息包括,:,?,模型的质量特性,-,模型质量是精确的,-,质心和,质量矩的值有一定误差。,?,单元矩阵系数,-,当单元矩阵系数最大,/,最小值的,比率,1.0E8,时将预示模型中的材料性质、实,常数或几何模型可能存在问题。当比值过高时,,求解可能中途退出。,?,模型尺寸和求解统计信息。,?,汇总文件和大小。,4.,没有获得结果的原因是什么,?,往往是求解输入的模型不完整或存在错误,典型原因有,:,?,约束不够,!(,通常出现的问题,),。,?,当模型中有非线性单元,(,如缝隙,gaps,、滑块,sliders,、铰,hinges,、索,cables,等),整体或部分结构出现崩溃或“松,脱”。,?,材料性质参数有负值,如密度或瞬态热分析时的比热值。,?,未约束铰接结构,如两个水平运动的梁单元在竖直方向,没有约束。,?,屈曲,-,当应力刚化效应为负(压)时,在载荷作用下整,个结构刚度弱化。如果刚度减小到零或更小时,求解存,在奇异性,因为整个结构已发生屈曲。,5.,应力奇异,应力奇异,是有限元模型中由于几何构造或载荷引起弹性,理论计算应力值无限大。即使是奇异点,材料的非线性特,性不可能允许应力值出现无限增大情况,在理论上总体应,变也是有限的(许多设计准则都是根据应力制订的,例如,设计疲劳曲线,但实际上是基于应变制订的)。,在应力奇异处,:,.,单元网格越是细化,越引起计算应力无限增加,并且,不再收敛。,.,网格疏密不均匀时网格离散误差也大小不一(自适应,网格划分结果是失败的或者网格错误)。,一般应力奇异发生情形:,?,?,?,添加在节点上的集中载荷(集中力)与施,加在与该节点相连单元上的均布或变化的,面载荷(压力)等相当的话,这些节点处,就成为应力奇异点。,离散约束点导致非零反力的出现,就如同,在节点上施加一集中力,这时约束点也就,成为应力奇异点。,锐利(零半径倒角)拐角处。,不常见的应力奇异情形:,?,?,?,由于在划分单元网格时出错,模型中存在,的“裂缝”。,曲边单元中处在极不理想位置的中间点,(,ANSYS,单元形状检查会发出警告)。,严重扭曲的单元(,ANSYS,单元形状检查,会发出警告)。,实际结构中并不存在应力奇

    注意事项

    本文(有限元法及应用知识点超全总结课件.ppt)为本站会员(小飞机)主动上传,三一办公仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知三一办公(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    备案号:宁ICP备20000045号-2

    经营许可证:宁B2-20210002

    宁公网安备 64010402000987号

    三一办公
    收起
    展开