工程仿真课程设计报告书Word格式文档下载.docx
《工程仿真课程设计报告书Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《工程仿真课程设计报告书Word格式文档下载.docx(23页珍藏版)》请在冰豆网上搜索。
2.1问题阐述…………………………………………………………………3
2.2理论分析……………………………………………………………………………4
2.3有限元分析…………………………………………………………………………5
2.4本章小结……………………………………………………………………………8
3.实体模型分析…………………………………………………………………………9
3.1模型材料与几何参数……………………………………………………9
3.2有限元模型的建立………………………………………………………9
3.3结果与分析……………………………………………………………………11
4.心得体会………………………………………………………………………18
5.参考文献
1前言
1.1引言
有限元法(finiteelementmethod)是一种高效能、常用的数值计算方法。
科学计算领域,常常需要求解各类微分方程,而许多微分方程的解析解一般很难得到,使用有限元法将微分方程离散化后,可以编制程序,使用计算机辅助求解。
有限元法在早期是以变分原理为基础发展起来的,所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中(这类场与泛函的极值问题有着紧密的联系)。
自从1969年以来,某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程,因而有限元法可应用于以任何微分方程所描述的各类物理场中,而不再要求这类物理场和泛函的极值问题有所联系。
基本思想:
由解给定的泊松方程化为求解泛函的极值问题。
1.2结构静力学分析
通过进行结构静力学的分析,我们能更加熟悉已学过的结构求解,结构受力的情况,结构静力学的分析不仅仅关系到我们的课程作业,也与我们日常生活,我们的工作密切相关,熟练的掌握结构静力学分析能帮助我们解决许多生活、工作中遇到的静力学问题。
1.3课程设计的目的和要求
课程设计能更好的让我们把所学的各科知识联系起来,把书本上的知识与生活实际所面对的问题联系起来,能有效的提高我们解决实际问题的能力。
通过做课程设计,能有效的提高学生自主学习的能力,提高学生查阅书刊杂志,自主选择查找相关资料的能力,进行独自分析解答的能力。
本次课程设计有两个问题,第一个为桁架的静力学分析,第二个为ANSYS设计模型要求采用理论计算和ANSYS软件两种方法分别进行分析,这是进一步让我们熟悉使用ANSYS软件,并对有限元理论进行更好的总结,让我们学习更好。
2平面桁架问题分析
2.1问题阐述
平面桁架问题在各种工程结构中被广泛使用的一种特殊的刚体系,它们是一些短而直的钢杆彼此以端部连接而成的几何不变的结构。
当刚杆之间的连接能近似的看做是铰链约束时,这种杆系结构称为桁架。
已知有一个桁架,结构形状如下图所示。
施加在桁架上的力F=10kN。
桁架每根长1m,杆的截面积为0.01㎡,泊松比μ=0.3,弹性模量E=200GPa。
用两种方法分别计算出桁架的内力并进行比较。
图2-1桁架计算模型图
图2-2桁架受力分析图
2.2理论分析
求支座反力,分别对A,B点取矩:
(2-1)
解得
(2-2)
解得
以节点A为研究对象:
:
(2-3)
(2-4)
(压力)
以节点C为研究对象:
(2-5)
解得
(拉力)
(压力)
以节点D为研究对象:
(2-6)
(2-7)
以节点E为研究对象:
(2-8)
(2-9)
(拉力)
以节点F为研究对象:
(2-10)
(2-11)
以节点G为研究对象:
(2-12)
(2-13)
2.3有限元分析
(1)选择单元类型:
通过单击Preference>
ElementType>
Add/Edit/Delete定义单元类型为LINK3Dfinitstn180。
(2)定义几何特性:
通过单击Preprocessor>
RealConstants>
Add/Edit/Delete定义单元的几何特性,即:
面积为0.01㎡。
(3)定义材料特性:
MaterialProps>
MaterialModels,以及Structural>
Linear>
Elastic>
Isotropic定义弹性模量为2E11,泊松比为0.3。
(4)定义节点:
Modeling>
Create>
Nodes>
InActiveCS,定义13个节点,它们分别是1(0,0,0),2(5,0,0),3(0,1,0),4(1,1,0),5(1,2,0),6(2,1,0),7(2,2,0),8(3,1,0),9(3,2,0),10(4,1,0),11(4,2,0),12(5,1,0),13(2.5,1,0)。
(5)定义单元:
通过节点生成单元,方法如下,单击Preprocessor>
Elements>
AutoNumbered>
ThruNodes,弹出“节点选择”对话框,依次选择1、3,单击Apply,如此生成1号单元。
同理单击Apply生成剩余的22个单元。
(6)施加位移约束:
依次单击Solution>
DefineLoads>
Apply>
Structural>
Displacement>
OnNodes,在1号节点施加AllDOF约束,在2号节点施加UY=0的约束。
(7)施加集中力载荷:
Force/Moment>
OnNodes在6号和8号节点施加F=-10000N的荷载。
(8)求解结果:
Solve>
CurrentLS,求出结果。
图2-3桁架有限元分析加载图
图2-4桁架总体变形图
图2-5桁架整体应变图
表2-1桁架杆件内力表
Element
MF(N)
1
-10000.
2
10000.
3
0.17618E-09
4
-14142.
5
0.23717E-09
6
-20000.
7
10000.
8
14142.
9
0.10842E-09
10
-0.54210E-10
11
20000.
12
13
0.81315E-10
14
15
-0.27105E-10
16
17
18
0.0000
19
20
21
22
23
退出系统
ANSYSUtilityMenu:
File>
Exit…>
SaveEverything>
OK
2.4本章小结
通过对比手动理论计算的结果和ANSYS软件分析的结果,我们可以发现,在理论计算杆力为0的地方,ANSYS分析结果为极小的数值,两者相差并不大,这说明了在面对平面桁架结构这种问题时,理论计算的结果还是大体正确的,但是在面对更为复杂的问题时,我们就要求助于ANSYS等工程类软件的应用了。
在这个平面桁架问题上,当我们运用ANSYS软件时,能够更加快速,更加精确地计算出问题的结果,并且误差小,而手动理论计算的误差更大,计算过程中也有可能出错,所以,熟练运用ANSYS软件解决工程问题,是我们必须学会的一项技能。
3实体模型分析
3.1模型材料与几何参数
现已知一扇门,高度为2m,宽度为0.9m,厚度为0.05m,要求对于门铰链的位置进行优化分析,那么我们可以设定三组不同位置的门铰链位置数据进行对比,设置门的弹性模量为2e10,泊松比为0.3,进行相关计算。
3.2有限元模型的建立
(1)定义工作名、工作标题、过滤参数:
定义工作名:
Utilitymenu>
Jobname
工作标题:
ChangeTitle(XX)
(2)设置解题类型:
启动ANSYS,单击Preference交互界面对话框,选中Structural,来对后面的分析进行菜单及相应的图形界面过滤,单击OK按钮。
(3)选择单元类型:
MainMenu>
Preprocessor>
Add/Edit/Delete>
Add>
Solid>
10node187
(4)设置材料属性:
定义材料的弹性模量和泊松比:
MainMenu>
MaterialModels>
Isotropic定义弹性模量EX为2e10,泊松比PRXY为0.3
(5)实体建模,建立薄板模型:
Volumes>
Block>
ByDiemnsions:
X方向0.9m,Y方向0.05m,Z方向2m。
点击OK。
形成薄板模型。
图3-1实体模型图
(6)网格划分:
Meshing>
MeshTool>
Lines>
Set:
逐个点击薄板的各个边>
OK。
设置单元长度为0.1m,点击Mesh,选中薄板,点击Apply。
图3-2门框网格划分图
(7)施加约束:
Solution>
DefinesLoads>
OnNodes。
选择四个约束点,点击AllDOF,点击OK。
图3-3施加约束图
(8)施加载荷:
Solution>
Pressure>
OnAreas。
选择施加载荷的面,施加150N的力,点击OK。
(9)求解:
3.3结果与分析
ANSYSMainMenu:
GeneralPostproc>
PlotResult>
DeformedShape>
Def+undeformed>
Apply。
获得整体变形图。
ContourPlot>
NodalSolu>
DOFSolution。
获得Y方向位移图。
Stress。
获得门框应力分布图。
第一组:
施加约束于薄板两个顶端
图3-4门框变形图
图3-5Y方向位移图
图3-6门框应力分布图
第二组:
施加约束于薄板边缘距两端五分之一处
图3-7门框变形图
图3-8Y方向位移图
图3-9门框应力分布图
第三组:
施加约束于薄板距边缘两端十分之三处
图3-10门框变形图
图3-11Y方向位移图
图3-12门框应力分布图
表3-1第一组的最大最小应力值表
MINIMUMVALUES
NODE
87
67
VALUE
-107.25
-1063.9
-4647.8
3.1690
2.7677
MAXIMUMVALUES
6299.2
1493.5
145.18
6154.0
5602.9
表3-2第二组的最大最小应力值表
199
963
64
-142.43
-723.51
-2873.7
0.30130
0.26715
197
961
3058.8
581.17
74.517
3285.0
2964.9
表3-3第三组的最大最小应力值表
203
959
177
-74.383
-242.14
-1589.3
0.38053
0.37202
201
175
2461.2
287.30
9.7723
2825.3
2562.5
结论:
对比上面三组图形,可以看出,第一组在门框的两端安装门铰链时,会使门整体所受应力过于集中,从而造成整体结构的不稳定,因此应舍弃,而第三组的门铰链位于距门框两端十分之三处,则使应力在门框中部集中,使得门框易从中部受损,因此舍弃,第二组在距门框两端五分之一处安装门铰链,其受力与应变最为均匀,不会发生大的形变,也就不会轻易受损,所以是最合适的门铰链位置。
综上,经过ANSYS的计算与分析,可以明确指出第二组为最优方案,也就是铰链位置相对分散更易于对门框结构的保护。
4心得体会
这一次的工程仿真课程设计,不仅使我加深了对于使用有限单元法的理解,同时也让我能更熟练地运用ANSYS软件,进行相关的计算分析。
在没有接触到ANSYS这款软件之前,每当我们要进行有关力学的计算时,只能用纸和笔,花大量的时间进行缓慢计算,耗时耗精力,并且,人工进行的计算经常会因为一点不小心而失误,造成整个题目的计算错误,而重新查找错误的过程又要花去大量时间,当然,这还只是一些较简单的结构方面的计算,如果是更为接近实际的结构的计算,凭借人工计算,对于我们来说几乎是不可能的。
然而,有了ANSYS软件,这一切就变得方便许多,尤其是一些桁架,圆盘类的受力的计算,只需在软件中建立相关模型,加上条件,就能获得相应结果。
刚开始接触ANSYS时,由于对有限元,单元,节点,形函数等基本概念没有清楚的了解,学起来特别吃力。
通过这个我认识到,应该把以前的书都看一下。
把《材料力学》,《结构力学》等都进行了系统的复习。
之后才开始使用ANSYS软件。
ANSYS不是一门单独的学科,而是和很多专业知识联系在一起的。
只有把他们有机地结合在一起,才可以使自己的能力得到真正的提高。
比如这一次的的平面桁架分析,要求用20根左右的桁架搭建起一个结构,先要用结构力学相关知识进行手动计算求解,之后才要用ANSYS进行分析,我自己复习了有关桁架内力的计算知识,搭建了桁架结构,进行计算,之后再将桁架模型输入ANSYS中,加上材料性质,载荷,求出结果。
而第二项作业,则是对于我们建模能力的极大考验,因为题目要求的是实际生活中存在的结构,需要我们自己先把实体结构简化成符合要求的相关模型,再加上题目中要求的荷载,输入软件进行计算,光是建立模型就让我们动了好一番脑筋,却也给予了我们很多收获。
总的来说,通过这一次的课程设计的学习,使我更加认识到ANSYS软件的重要性,并且进一步熟悉了软件的使用,获益良多。
5参考文献
〔1〕王铎孙毅程靳编《理论力学》2009高等教育
〔2〕胡于进王璋奇著《有限元分析及应用》2009清华大学
〔3〕赵经文王红钰编《结构有限元分析(第2版)》2001科学
〔4〕傅永华编《有限元分析基础》2003XX大学
〔5〕王新敏李义强编《ANSYS结构分析单元与应用》2001人民交通