本科毕业设计论文有限元节课论文Word文件下载.docx
《本科毕业设计论文有限元节课论文Word文件下载.docx》由会员分享,可在线阅读,更多相关《本科毕业设计论文有限元节课论文Word文件下载.docx(11页珍藏版)》请在冰豆网上搜索。
由于单元体能按照不同的链接方式组合,且单元体又有不同的形状,因此可以模拟成不同几何形状的求解小区域,然后对单元进行力学分析。
最后再整合分析,这种化整为零,化单体为多体,最后在集零为整的方法就是有限元方法的主要思想。
1.有限元方法的发展历史与进程
有限元方法起源于需要解决市政工程和航空工程方面复杂的弹性结构分析问题。
它的开发可以追溯到A.Hrennikoff(1941)和R.Courant(1942)的工作。
虽然这些先驱者使用这些方法,并且引人注目的不同,但他们都共享一个基本的特性:
把连续域的网格离散化进入一组离散的子域里。
Hrennikoff的工作是采用格子使域离散,而与之类似,为了求解起源于汽缸扭转的问题的二阶椭圆的偏微分方程式(PDEs),RichardCourant的方法是把域划分成有限的三角形子域。
对于由Rayleigh,Ritz和Galerkin开发的偏微分方程式(PDEs),RichardCourant的贡献是改进,绘制了大量的早期结果。
针对机身和结构分析的有限元方法的开发最早开始于1950年代中期,并且用于市政工程的有限元方法许多是1960年代在伯克利开始启动(见伯克利早期有限元研究)。
在1973年Strang和Fix出版的《有限元方法的分析》里,提供的方法采用了严格的数学基础,并且已经在广泛变化的工程学科,即电磁和流体力学里,针对物理系统的数字建模,归纳成为应用数学的分枝。
在结构力学里,有限元方法的开发常常是基于能量理论,即虚功原理或最小总潜能原理,对于结构工程师来说,早就强烈要求提供综合的,直觉的和物理的依据。
到20世纪80年代初期国际上较大型的结构分析有限元通用程序多达几百种,从而为工程应用提供了方便条件。
由于有限元通用程序使用方便,计算精度高,其计算结果已成为各类工业产品设计和性能分析的可靠依据。
到目前为止,有限元分析方法已经成为各工程行业必不可少的可靠分析方法。
2有限元法的基本原理
在工程分析和科学研究中,常常会遇到大量的由常微分方程、偏微分方程及相应的边界条件描述的场问题,如位移场、应力场和温度场等问题。
目前求解这类场问题的方法主要有两种:
用解析法求得精确解;
用数值解法求其近似解。
其中,能用解析法求出精确解的只能是方程性质比较简单且几何边界相当规则的少数问题。
而对于绝大多数问题,则很少能得出解析解。
这就需要研究它的数值解法,以求出近似解。
目前,工程中实用的数值解法主要有三种:
有限差分法;
有限元法;
边界元法;
其中,以有限元法通用性最好,解题效率高,工程应用最广。
目前它已成为机械产品动、静、热特性分析的重要手段,它的程序包是机械产品计算机辅助设计方法库中不可缺少的内容之一。
2.1基本思想
弹性力学解法的问题弹性力学解法的问题在于:
不论是应力函数解法、扭转函数解法、挠曲函数解法、还是基于最小势能原理的瑞利-李兹等方法,其困难在于如何给出一个在全求解区给出一个在全求解区域上均成立的试探函数。
在有限单元法里,这个问题通过定义分片插值的位移或应力函数得到了巧妙的解决。
有限元方法基于变分原理和加权余量法,把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
2.2求解步骤
建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的出发点。
区域单元剖分,根据求解区域的形状及实际问题的物理特点,将区域剖分为若干相互连接、不重叠的单元。
区域单元划分是采用有限元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插值条件的插值函数作为单元基函数。
有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则。
单元分析:
将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;
再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程。
总体合成:
在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程。
边界条件的处理:
一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。
对于自然边界条件,一般在积分表达式中可自动得到满足。
对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足。
解有限元方程:
根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。
3简单悬架梁实例应用分析
实验设置:
固定钢直尺一端,在另一端施加不同的载荷,计算施加载荷点的变形量。
实验采用实际测量数据ansys有限元分析数据与力学理论计算数据相比较的方法来相互验证。
3.1力学模型简化
此实验可以简化为简单悬臂梁受力模型
3.2实际测量数据:
(实验操作图1)(实验操作图2)
当在直尺长度为275mm处施加载荷时
变形量/
砝码质量/
当在直尺长度为225mm处施加载荷时
当在直尺长度为175mm处施加载荷时
钢条下降量/
3.3理论值数据
钢条的变形量以挠度来计,第一种情况即距离固定端处施加荷载,弹性模量,惯性矩,则其挠度,计算可得:
第一种情况即距离固定端275mm处施加载荷
挠度/
15.8
22.5
第二种情况即距离固定端225mm处施加荷载
11.9
13.2
第三种情况即距离固定端175mm处施加荷载
5.3
6.3
3.4ANSYS软件分析
步骤1:
画出模型图
用PROE画出模型图:
如下图所示
步骤2:
打开ANSYS,选择StaticStructural(静力结构分析)模块,界面如图所示。
步骤3:
选择材料,本例子选用的是不锈钢材料;
如图所示。
步骤4:
导入模型图:
注1:
在导入PROE图时要将其后缀删掉,不然导入不能生成。
如下图
变为
注2:
模型图中不能有填充面,必须要有一定的高度,不然导入不能分析,如下图,(PROE最小厚度为0.067mm)。
步骤5:
双击Model,进入结构分析界面。
如下图所示;
步骤6:
给模型赋予材料属性。
本例子为不锈钢;
步骤7:
网格划分。
本例子为自动网格划分,只是将RelevanceCenter(关联中心)的Coarse(粗糙)变成Medium(中等)。
网格划分如下图所示。
施加载荷与约束;
施加约束,本例子是在一端固定约束。
施加载荷,本例子是秤砣的重力为载荷。
在另一端施加力载荷。
如下图所示:
步骤8:
后处理;
查看模型的总变形。
可以看到200g直径为秤砣,长300mm,宽25mm,厚1.15mm的钢直尺。
最大的变形量为24.859mm。
导入模型,重复以上步骤;
得到下面数据图。
1.圆心距离固定端275mm,200g直径50mm为秤砣应力,长300mm,宽25mm,厚1.15mm的钢直尺。
2.圆心距离固定端275mm,250g直径为50mm秤砣,长300mm,宽25mm,厚1.15mm的钢直尺。
最大的变形量为31.074mm。
3.圆心距离固定端225mm,200g直径为50mm秤砣,长300mm,宽25mm,厚1.15mm的钢直尺。
最大的变形量为13.959mm。
4.4.圆心距离固定端225mm,250g直径为50mm秤砣,长300mm,宽25mm,厚1.15mm的钢直尺。
最大的变形量为17.448mm。
5.5.圆心距离固定端175mm,200g直径为50mm秤砣,长300mm,宽25mm,厚1.15mm的钢直尺。
最大的变形量为6.8176mm。
6.圆心距离固定端175mm,250g直径为50mm秤砣,长300mm,宽25mm,厚1.15mm的钢直尺。
最大的变形量为8.552mm。
最后ANSYS项目管理界面。
(因为钢直尺的材料一样,可以材料共享,不用每次导入输入材料属性)
3.5数据比对分析
第一种情况为施加载荷点出为固定端275mm处
变形量(mm)
载荷(g)
200
250
与理论值对比(%)
实验数据
23
32
68.58
70.19
理论数据
100
Ansys分析数据
24.895
31.074
64.25
72.40
第二种情况为施加载荷点出为固定端225mm处
16.3
18.3
72.98
71.94
13.959
17.448
85.25
75.65
第一种情况为施加载荷点出为固定端175mm处
7
8.7
75.44
72.19
6.817
8.552
77.73
73.67
通过实验数据的对比,可以知道,同样的情况下,施加相同载荷时,ansys分析数据与实验测量数据相接近,并且都与理论计算数据有相同量级的偏差,所以ansys软件分析模拟的结果与实际结果一致,可以作为工程中对实际情况的可靠分析与仿真的参考依据。
4心得与展望
有限元方法是求解复杂工程问