ANSYS有限元分析课件.ppt
《ANSYS有限元分析课件.ppt》由会员分享,可在线阅读,更多相关《ANSYS有限元分析课件.ppt(935页珍藏版)》请在咨信网上搜索。
1、-1-2024/2/25 周日周日ANSYSANSYS-2-II.ANSYS软软件及其件及其应应用用2024/2/25 周日周日2第第1章章有限元基本理论有限元基本理论第第2章章ANSYS功能简介功能简介第第3章章ANSYS基本过程基本过程第第4章章ANSYS入门与准备入门与准备第第5章章模型输入及修复模型输入及修复第第6章章坐标系坐标系第第7章章选择、组件与部件选择、组件与部件第第8章章实体建模技术实体建模技术第第9章章布尔操作布尔操作第第10章章单元属性单元属性第第11章章网格划分网格划分第第12章章加载求解技术加载求解技术第第13章章后处理技术后处理技术第第14章章结构非线性分析结构非线
2、性分析第第15章章模态分析模态分析第第16章章耦合和约束方程耦合和约束方程第第17章章APDL基础基础第第18章章子模型子模型第第19章章热分析热分析第第20章章热热-应力耦合分析应力耦合分析-3-第一章第一章有限元基本理论有限元基本理论平衡方程平衡方程几何方程几何方程物理方程物理方程边界条件边界条件物理系统物理系统有限元离散有限元离散单元的位移场单元的位移场(假定单元内位移函数假定单元内位移函数)单元节点关系单元节点关系求解区域的位移场、应力场求解区域的位移场、应力场简简单单化化-4-1.1有限元分析有限元分析(FEA)有限元分析有限元分析 是利用数学近似的方法对真实物理是利用数学近似的方法
3、对真实物理系统(几何和载荷工况)进行模拟。它利用简系统(几何和载荷工况)进行模拟。它利用简单而又相互作用的元素,即单元,用有限数量单而又相互作用的元素,即单元,用有限数量的未知量去逼近无限未知量的真实系统。的未知量去逼近无限未知量的真实系统。-5-1.2有限单元法的基本思想有限单元法的基本思想将将连连续续的的结结构构离离散散成成有有限限个个单单元元,并并在在每每一一单单元元中中设设定定有有限限个个节节点点,将将连连续续体体看看作作只只在在节节点点处处相相连连接接的的一组单元的集合体。一组单元的集合体。选选定定场场函函数数的的节节点点值值作作为为基基本本未未知知量量,并并在在每每一一单单元元中中
4、假假设设一一近近似似插插值值函函数数,以以表表示示单单元元中中场场函函数数的的分分布规律。布规律。利利用用力力学学中中的的某某种种变变分分原原理理去去建建立立用用以以求求节节点点未未知知量量的的有有限限单单元元法法方方程程,将将一一个个连连续续域域中中有有限限自自由由度度问问题化为离散域中有限自由度问题。题化为离散域中有限自由度问题。-6-1.3物理系统举例物理系统举例几何体几何体 载荷载荷 物理系统物理系统结构结构热热电磁电磁-7-1.3.1平衡方程平衡方程-8-1.3.2几何方程几何方程-9-1.3.3物理方程物理方程(本构方程本构方程)拉梅系数拉梅系数体积应变体积应变剪切模量剪切模量-1
5、0-1.3.4边界条件边界条件应力边界条件应力边界条件位移边界条件位移边界条件-11-1.4有限元模型有限元模型真实系统真实系统有限元模型有限元模型 有限元模型有限元模型 是真实系统理想化的数学抽象是真实系统理想化的数学抽象。-12-1.5自由度自由度(DOFs)自由度自由度(DOFs)用于描述一个物理场的响应特性用于描述一个物理场的响应特性。结构结构 DOFs 结构结构 位移位移 热热 温度温度 电电 电位电位 流体流体 压力压力 磁磁 磁位磁位 问题问题 自由度自由度ROTZUYROTYUXROTXUZ-13-1.6节点和单元节点和单元节点节点:空间中的坐标位置,具有一定空间中的坐标位置,
6、具有一定自由度和自由度和存在相互存在相互物理作用物理作用。单元单元:一组节点自由度间相互作用的一组节点自由度间相互作用的数值、矩阵描述(称为刚度或系数数值、矩阵描述(称为刚度或系数矩阵矩阵)。单元有线、面或实体以及二。单元有线、面或实体以及二维或三维的单元等种类。维或三维的单元等种类。有限元模型由一些简单形状的有限元模型由一些简单形状的单元单元组成,单组成,单元之间通过元之间通过节点节点连接,并承受一定连接,并承受一定载荷载荷。载荷载荷载荷载荷-14-1.6节点和单元节点和单元(续续)信息是通过单元之间的公共节点传递的。信息是通过单元之间的公共节点传递的。分离但节点重叠的分离但节点重叠的单元单
7、元A和和B之间没有之间没有信息传递(需进行信息传递(需进行节点合并处理)节点合并处理)具有公共节点具有公共节点的单元之间存的单元之间存在信息传递在信息传递.AB.AB.1 node2 nodes-15-1.6节点和单元节点和单元(续续)节点自由度是随连接该节点节点自由度是随连接该节点 单元类型单元类型 变化的。变化的。JIIJJKLILKIPOMNKJIL三维杆单元三维杆单元(铰接铰接)UX,UY,UZ三维梁单元三维梁单元UX,UY,UZ,ROTX,ROTY,ROTZ二维或轴对称实体单元二维或轴对称实体单元UX,UY三维四边形壳单元三维四边形壳单元UX,UY,UZ,ROTX,ROTY,ROTZ
8、三维实体热单元三维实体热单元TEMPJPOMNKJIL三维实体结构单元三维实体结构单元UX,UY,UZ-16-1.7单元形函数单元形函数FEA仅仅求解节点处的仅仅求解节点处的DOF值。值。单元单元形函数形函数是一种数学函数,规定了从节点是一种数学函数,规定了从节点DOF值值到单元内所有点处到单元内所有点处DOF值的计算方法。值的计算方法。因此,单元形函数提供出一种描述单元内部结果的因此,单元形函数提供出一种描述单元内部结果的“形状形状”。单元形函数描述的是给定单元的一种单元形函数描述的是给定单元的一种假定假定的特性。的特性。单元形函数与真实工作特性吻合好坏程度直接影响单元形函数与真实工作特性吻
9、合好坏程度直接影响求解精度。求解精度。-17-真实的二次曲线真实的二次曲线.节点节点单元单元 二次曲线的线性近似二次曲线的线性近似 (不理想结果不理想结果).21.7单元形函数单元形函数(续续)节点节点单元单元 DOF值二次分布值二次分布.1节点节点 单元单元 线性近似线性近似(更理想的结果更理想的结果)真实的二次曲线真实的二次曲线.3节点节点单元单元二次近似二次近似(接近于真实的二次近似拟合接近于真实的二次近似拟合)(最理想结果最理想结果).4-18-1.7单元形函数单元形函数(续续)DOF值值可可以以精精确确或或不不太太精精确确地地等等于于在在节节点点处处的的真真实实解解,但单元内的平均值
10、与实际情况吻合得很好。但单元内的平均值与实际情况吻合得很好。这些平均意义上的典型解是从单元这些平均意义上的典型解是从单元DOFs推导出来的推导出来的(如:结构应力、热梯度)。(如:结构应力、热梯度)。nodalsolutionelementsolutionUX,UY,UZ,ROTX,ROTY,ROTZ、E-19-1.7单元形函数单元形函数(续续)如如果果单单元元形形函函数数不不能能精精确确描描述述单单元元内内部部的的DOFs,就就不不能能很很好好地地得得到到导导出出数数据据,因因为为这这些些导导出出数数据据是是通通过单元形函数推导出来的。过单元形函数推导出来的。当当选选择择了了某某种种单单元元
11、类类型型时时,也也就就十十分分确确定定地地选选择择并并接接受受该种单元类型所假定的单元形函数。该种单元类型所假定的单元形函数。在在选选定定单单元元类类型型并并随随之之确确定定了了形形函函数数的的情情况况下下,必必须须确确保保分分析析时时有有足足够够数数量量的的单单元元和和节节点点来来精精确确描描述述所所要求解的问题。要求解的问题。-20-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题-21-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题(续续)就就整整个个直直杆杆来来说说,位位移移函函数数U(x)是是未未知知的的,但但对对每每一一单单元元可可以以近近似似地地假假设设一一位位移移
12、函函数数,它它在在结结点点上上等等于于结结点点位位移移。此此处处,假假设设单单元元中中的的位位移移按按线线性性分分布布,即:即:-22-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题(续续)有了位移插值函数,就可以按材料力学公式求出应变有了位移插值函数,就可以按材料力学公式求出应变和应力用节点位移表示的公式:和应力用节点位移表示的公式:-23-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题(续续)外载荷与结点的平衡方程外载荷与结点的平衡方程为第为第i个结点上承受的外载荷个结点上承受的外载荷-24-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题(续续)假定将直杆分割成假定
13、将直杆分割成3个单元,每个单元长为个单元,每个单元长为a=L/3,则,则对结点对结点2,3,4列出的平衡方程为:列出的平衡方程为:-25-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题(续续)-26-1.8直杆受自重作用的拉伸问题直杆受自重作用的拉伸问题(续续)联立求解线性代数方程组得:联立求解线性代数方程组得:-27-1.9有限单元法解题的一般步骤有限单元法解题的一般步骤结构的离散化结构的离散化选择位移模式选择位移模式建立平衡方程建立平衡方程求解节点位移求解节点位移计算单元中的应力和应变计算单元中的应力和应变-28-1.9.1结构的离散化结构的离散化将将分分析析的的结结构构物物分分割
14、割成成有有限限个个单单元元体体,使使相相邻邻的的单单元元体体仅仅在在节节点点处处相相连连接接,而而以以如如此此单单元元的的结结合合体体去代替原来的结构。去代替原来的结构。-29-1.9.2选择位移模式选择位移模式(形函数形函数)首先对单元假设一个位移差值函数,或称之为位移模首先对单元假设一个位移差值函数,或称之为位移模式,得到用节点位移表示单元体内任一点的唯一的式,得到用节点位移表示单元体内任一点的唯一的关系式关系式有了位移模式,就可利用几何关系和应力有了位移模式,就可利用几何关系和应力-应变关系表应变关系表出用单元节点位移表示单元中应变和应力的表达式出用单元节点位移表示单元中应变和应力的表达
15、式-30-1.9.3三角形单元的形函数三角形单元的形函数基基本本假假定定:假假定定单单元元内内的的位位移移可可以以用用一一个个比比较较简简单单的的函函数数来来表表示示,如如线线性性插插值值函函数数。这这在在单单元元划划分分比比较较密密的的情情况况下下是是合理可行的。合理可行的。-31-1.9.3三角形单元的形函数三角形单元的形函数(续续)将将三三角角形形单单元元的的3个个顶顶点点的的2个个方方向向位位移移代代入入位位移移函函数数可可求求出出6个个待待定定系系数数。即即可可用用节节点点的的位位移移表表示示内内部任意一点的位移:部任意一点的位移:-32-1.9.4建立平衡方程建立平衡方程可可利利用
16、用最最小小势势能能原原理理建建立立结结构构的的节节点点载载荷荷和和节节点点位位移移之间的关系式,即结构的平衡方程之间的关系式,即结构的平衡方程1.9.5 求解结点位移求解结点位移v将边界条件代入线性代数方程组将边界条件代入线性代数方程组 后,后,经解算可求得所有未知的结点位移。经解算可求得所有未知的结点位移。-33-1.9.6计算单元中的应变和应力计算单元中的应变和应力依据求得的结点位移,由依据求得的结点位移,由可求得单元中任一点的应变和应力。可求得单元中任一点的应变和应力。-34-平面问题的有限单元法平面问题的有限单元法-35-结构的离散化结构的离散化 用有限元法对结构进行应力分析时,首先要
17、将结构进用有限元法对结构进行应力分析时,首先要将结构进行离散化。即将一个连续体看成由有限个单元组成的行离散化。即将一个连续体看成由有限个单元组成的体系。弹性力学平面问题中最常见的单元是三角形单体系。弹性力学平面问题中最常见的单元是三角形单元。元。所有作用在单元上的载荷都按静力等效的原则移置到所有作用在单元上的载荷都按静力等效的原则移置到结点上,并在受几何约束的结点处设置相应的铰支座。结点上,并在受几何约束的结点处设置相应的铰支座。这样就得到了用以代替原来弹性体的有限单元计算模这样就得到了用以代替原来弹性体的有限单元计算模型。型。-36-位移模式位移模式 取一个典型的三角形单元进行力学分析。在有
18、限单元位移法取一个典型的三角形单元进行力学分析。在有限单元位移法中,假设结点上的位移是基本未知量。为了能用单元的结点中,假设结点上的位移是基本未知量。为了能用单元的结点位移表示单元中的应变和应力分量,必须假定一个位移模式,位移表示单元中的应变和应力分量,必须假定一个位移模式,也就是说根据单元的结点位移去构造单元上的位移插值函数。也就是说根据单元的结点位移去构造单元上的位移插值函数。-37-位移模式(续)位移模式(续)-38-位移插值函数位移插值函数 采用线性插值,即假定单元上的位移分量是坐标的线采用线性插值,即假定单元上的位移分量是坐标的线性函数:性函数:它们可以由结点位移确定如下:它们可以由
19、结点位移确定如下:-39-位移模式位移模式(续续)联立求解上述方程,可得:联立求解上述方程,可得:-40-位移模式位移模式(续续)其中:其中:而:而:是三角形是三角形ijm的面积。的面积。-41-位移模式位移模式(续续)于是可以得到:于是可以得到:其中:其中:同理得:同理得:-42-位移模式位移模式(续续)可以将位移模式改写为矩阵模式:可以将位移模式改写为矩阵模式:-43-单元中的应变和应力单元中的应变和应力 有了单元的位移模式,就可以借助平面问题的几何和有了单元的位移模式,就可以借助平面问题的几何和物理方程,导出用单于的结点位移表示单元中的应变物理方程,导出用单于的结点位移表示单元中的应变和
20、应力分量的公式。和应力分量的公式。由:由:-44-单元中的应变和应力单元中的应变和应力(续续)得到:得到:或简写为:或简写为:-45-单元中的应变和应力单元中的应变和应力(续续)将应变代入物理方程:将应变代入物理方程:可得:可得:即为用单元中的结点位移表示单元中应力的关系式。即为用单元中的结点位移表示单元中应力的关系式。-46-单元中的应变和应力单元中的应变和应力(续续)式中式中D为弹性矩阵,对于平面应力问题,矩阵为:为弹性矩阵,对于平面应力问题,矩阵为:-47-单元的总势能单元的总势能我们已经知道由各个单元的位移模式就形成了整个结构我们已经知道由各个单元的位移模式就形成了整个结构的位移模式。
21、按弹性力学最小势能原理,结构中最接的位移模式。按弹性力学最小势能原理,结构中最接近于真实解的位移应该是使结构总势能取得最小值的近于真实解的位移应该是使结构总势能取得最小值的那组位移函数。那组位移函数。由于在位移函数公式中,结点位移为自变量,这样就使由于在位移函数公式中,结点位移为自变量,这样就使一个泛函的极值问题变为一个多元函数的极值问题。一个泛函的极值问题变为一个多元函数的极值问题。为此我们来讨论单元的总势能关于结点位移的表达式。为此我们来讨论单元的总势能关于结点位移的表达式。每一个单元的总势能由该单元的应变能以及此单元上所每一个单元的总势能由该单元的应变能以及此单元上所有外力的势能组成。有
22、外力的势能组成。-48-单元的应变能单元的应变能 平面应力状态下,设物体厚度为平面应力状态下,设物体厚度为h,则单元中的应变能,则单元中的应变能为:为:-49-单元的应变能(续)单元的应变能(续)将将和和Bi代入上式,应用矩阵相乘的转置的逆序法则代入上式,应用矩阵相乘的转置的逆序法则,注意到弹性矩阵注意到弹性矩阵D的对称性的对称性,有:有:-50-单元的应变能(续)单元的应变能(续)因为矩阵因为矩阵B及及D的元素都是常量,所以可记:的元素都是常量,所以可记:-51-单元的应变能(续)单元的应变能(续)从而单元的应变能可写为:从而单元的应变能可写为:利用利用=Be,有:,有:-52-单元的应变能
23、(续)单元的应变能(续)注意到注意到B=Bi Bj Bm,记子矩阵,记子矩阵-53-单元上体积力的势能单元上体积力的势能物体中常见的体力为旋转离心体力和重力。在平面问题中,物体中常见的体力为旋转离心体力和重力。在平面问题中,体积力在体积力在z轴方向的分力为零,设单元体积中的体积力为:轴方向的分力为零,设单元体积中的体积力为:单元上体积力具有的势能为:单元上体积力具有的势能为:-54-单元上表面力的势能单元上表面力的势能设物体边界上一单元某边上受到表面力的作用,单位长度上设物体边界上一单元某边上受到表面力的作用,单位长度上所受到的表面力为:所受到的表面力为:则单元上表面力的势能为:则单元上表面力
24、的势能为:-55-单元节点上集中力的势能单元节点上集中力的势能如果弹性物体受到集中力如果弹性物体受到集中力Re 的作用,通常划分单元网格时的作用,通常划分单元网格时都在集中力的作用点设置结点。设某单元都在集中力的作用点设置结点。设某单元3个结点上所受个结点上所受到的集中力为:到的集中力为:于是该单元上集中力的势能是:于是该单元上集中力的势能是:-56-单元中的总势能单元中的总势能综合前面的几种情况,可以得到单元中的总势能为:综合前面的几种情况,可以得到单元中的总势能为:-57-单元中的总势能单元中的总势能分别引进单元体积力,表面力,集中力向量如下:分别引进单元体积力,表面力,集中力向量如下:-
25、58-单元中的总势能单元中的总势能则单元中的总势能可以表示为:则单元中的总势能可以表示为:-59-物体中的总势能物体中的总势能把各单元的总势能叠加起来,就可得到整个弹性体的总势能。把各单元的总势能叠加起来,就可得到整个弹性体的总势能。为了便于叠加和归并,需将单元刚度矩阵表达式为了便于叠加和归并,需将单元刚度矩阵表达式(2-18)作作适当的改写。适当的改写。假设结构离散化后共有假设结构离散化后共有n个结点,将编号为个结点,将编号为 l的结点位移记为:的结点位移记为:则结构的结点位移向量:则结构的结点位移向量:是一个是一个2n维的列向量。维的列向量。-60-物体中的总势能(续)物体中的总势能(续)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 有限元分析 课件
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【胜****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【胜****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。