精品版有限元综合课程设计Word文档格式.docx

精品版有限元综合课程设计Word文档格式.docx

《精品版有限元综合课程设计Word文档格式.docx》由会员分享，可在线阅读，更多相关《精品版有限元综合课程设计Word文档格式.docx（21页珍藏版）》请在冰豆网上搜索。

5.3加载完成11

5.4求解11

6.后处理12

6.1查看变形图12

6.2查看应力16

6.3查看三维立体图19

7.结论21

8.参考文献22

1.用AtuoCAD绘制加农炮身管并标注尺寸

1.1绘制炮身结构

由提供的100mm加农炮身管结构图，应用AutoCAD2012软件采用第一角投影法，严格按照GB14665-2012《机械工程CAD制图》标准规定制图，炮身零件工程图见附录附录一。

1.2100mm加农炮身管说明

（1）图纸名称：

100mm加农炮身管图；

（2）绘图比例：

1:

20；

（3）数量：

1；

（4）绘图方法：

半剖式。

1.3在绘制过程中应注意的问题

（1）绘图前先设置图限、图层；

（2）不同类型的对象绘制在不同的图层上，利于以后修改；

（3）设置线宽显示比例，使粗细线显示协调；

（4）内螺纹的大径为细实线，小径为粗实线；

（5）标注前先设置标注样式，包括非圆直径、角度、标准标注等；

（6）将图形缩小到A4图纸后，应注意放大线宽。

2.利用AutoCAD绘制加农炮身管截面图

建立有限元模型有两种方法，输入法和创建法。

输入法是直接输入由其他CAD软件创建好的实体模型，创建法是在ANSYS中从无到有地创建实体模型，两者并不是完全分开的。

但由于ANSYS在创建模型是过于繁琐，本次课程设计采用输入法建立有限元模型。

利用AUTOCAD2012,以y轴为火炮身管的对称轴绘制截面图的右侧部分，并利用CAD中面域功能将所绘制的图形面域化。

将上述绘制好的图形以“SAT”文件格式导出,保存并命名。

3.问题分析

3.1100mm火炮身管静态受力分析

身管的材料参数为：

弹性模量E=207GPa,泊松比u=0.3。

该模型后端施加X轴方向固定位移约束以及从坡膛到身管尾部的Y轴位移约束，在身管的膛底和坡膛之间内侧表面垂直内表面方向上施加183MPa的均布载荷。

模型如图3-1所示。

图3-1

3.2确定问题的范围

在对问题有了初步了解以后，进一步要弄明白问题涉及的领域，即确定问题的范围。

对问题分析可知本问题属于结构力学领域中的静态分析问题且在材料变形的线性变形的范围内。

4.建立有限元模型

建立模型包括设定分析作业和标题；

定义单元类型和实常数；

定义材料属性；

导入几何模型；

划分有限元网格。

其具体步骤为：

4.1设定分析作业名、标题以及导图模型

在进行一个新的有限元分析时，通常需要修改数据库名，并在图形输出窗口中定义一个标题来说明当前进行的工作内容。

另外，对于不同的分析范畴（结构分析，热分析，流体分析，电磁场分析等），ANSYS所用的主菜单的内容不尽相同，为此，需要在分析开始时选定分析内容的范畴，以便ANSYS显示出与其相对应的菜单选项。

a.定义工作标题和更改储存路径:

打开ANSYS14.0MechanicalAPDLProductLauncher,出现如下如图4-1所示的对话框，定义工作标题和更改储存路径。

图4-1

b.导入模型及修改文件名：

执行File-Import—SAT…—浏览上述模型，如图4-2所示；

执行File-ChangeTile,如图4-3所示；

模型图如下图4-4所示。

图4-2

图4-3

图4-4

c.定义分析类型：

执行MainMenu-Preferences,选择结构分析，如图4-5所示。

图4-5

4.2定义单元类型

在进行有限元分析时，首先应该根据分析问题的几何结构，分析类型和所分析的问题精度要求等，选定合适具体分析的单元类型。

本次选用四节点四边形板单元PLANE182。

PLANE182不仅可用于计算平面应力问题，还可以用于分析平面应变和轴对称问题。

单元类型对话框如图4-6所示。

图4-6

a.定义单元类型：

执行MainMenu-Preprocessor-ElementType-Add，弹出LibraryofElementTypes对话框，选择固体下4节点182单元，如图4-7所示。

图4-7

b.单元行为方式选择：

单击option按钮，打开如图4-8所示的对话框，在Elementbehavior项选择Axisymmetri（轴对称）。

图4-8

4.3定义材料属性

本次选用的单元类型不需要定义实常数，故直接定义材料属性。

a.依次单击Structural-Linear-Elastic-Isotropic,如图4-9所示。

展开材料属性的树形结构，打开材料的弹性模量EX和泊松比PRXY的定义对话框，如图4-10所示，在对话框的EX文本框中输入弹性模量2.07e11，在PRXY文本框中输入泊松比0.3。

图4-9

图4-10

b.依次Structural-Density,打开材料密度定义对话框如图4-11所示。

图4-11

4.4网格划分

a.执行MainMenu-Preprocessor-meshing-MeshTool命令，打开MeshTool（网格工具），如图4-12所示。

图4-12

b.自定义网格大小：

点击lines域的set按钮，出现对话框如图4-13所示，要求选择定义单元个数的线，在模型图上选择炮膛内壁的线；

接着在对话框中点击Apply,出现如图4-14所示的对话框，在单元数目该项上输入60，点击Apply，再次选择坡膛内壁的线，点击Apply，输入40，其余地方选取30个单元，最后点击OK。

图4-13

图4-14

c.在Mesh栏中选择选择默认选项，点击Mesh按钮，出现如图4-13的对话框，在对话框中选择PickAll，ANSYS将按照对线的控制进行网格的划分，划分后的结果如图4-15所示。

图4-15

5.加载以及求解

5.1添加约束位置

执行Mainmenu-Solution-DefineLoad-apply-structural-Displacement-onLines出现选择对话框，选择最下端的线，单击OK按钮出现ApplyU,ROTonLines对话框，选择UY，添加炮管在Y方向的约束，如图5-1所示。

图5-1

5.2施加载荷

执行Mainmenu-Solution-DefineLoad-apply-structural-Pressure-onlines，打开选择线对话框，选择身管的坡膛到身管尾尾内侧表面线，单击OK按钮，出现ApplyPRESonlines对话框,如图5-3所示，在LoadPRESvalue文本框中输入均布载荷，大小为1.83e8，单击OK按钮施加均布载荷。

图5-3

5.3加载完成

加载以及施加约束后的模型图如图5-4所示。

图5-4

5.4求解

执行Mainmenu-Solution-Solve-CurrentLS,出现对话框，依次点击OK进行求解，直至出现如图5-5所示的对话框，求解完成！

图5-5

6.后处理

求解完成后，就可以利用ANSYS软件生成的结果文件（对于静力分析，就是Jobname.RST）进行后处理。

静力分析中通常通过POST1后处理器就可以处理和显示大多数兴趣的结果数据。

6.1查看变形图

关键的变形为径向变形，在高膛压下膛底至坡膛处径向产生过大变形，可能导致身管破裂。

a.执行ANSYSMainMenu-GeneralPostproc-PlotResults-ContourPlot-NodalSolu出现ContourNodalSolutionData（等值线显示节点数据）对话框，如图6-1所示。

图6-1

b.Itemtobecontourd（等值线显示结果项）域中选择DOFsolution（自由度解）选项，在列表中选择X-Componentofdisplacement（X向位移）选项，此时，X向位移即为径向位移。

选择Deformedshapewithundeformededge（变形后和未变形轮廓线）选项，如图6-2所示。

图6-2

c.单击OK按钮，在图形窗口中显示变形图，包含变形前的轮廓线，X向位移图如图6-2所示。

图中下方的色谱表面不同颜色对应的数值。

其局部放大图如图6-3所示。

图6-2X向位移

图6-3X向位移局部放大

d.用同样的方法显示Y向位移，其位移图和局部放大图如图6-4、6-5所示。

图6-4Y向位移

图6-5Y向位移局部放大

6.2查看应力

火炮身管膛内压力垂直于内表面，主要应力也是径向应力。

a.与查看应变类似，执行：

MainMenu-GeneralPostproc-PlotResults-ContourPlot-Nodal命令对话框，打开ContourNodalSolutionNata（等值线显示节点解数据）对话框。

b.在Itemtobecontourd（等值线显示结果项）域中选择DOFsolution（自由度解）选项，在列表中选择X-Componentofstress（X向应力）选项。

c.选择Deformedshapewithundeformededge（变形后和未变形轮廓线）选项，如图6-6所示：

点击OK,即可查看其X向应力图。

X向应力图及其局部放大图如6-7、6-8所示。

图6-6

图6-7X向应力

图6-8X向应力局部放大

d.用同样的方法查看模型的Y向应力图及其局部放大图如图6-10、6-11所示。

图6-10Y向应力

图6-11Y向应力局部放大

6.3查看三维立体图

a.从应用菜单中选择UtilityMenu：

PlotCtrls-Style-SymmetricExpansion-2DAxi–Symmetric命令，打开2DAxi-SymmetricExpansion对话框，如图6-12所示。

图6-12

b.在Selectexpansionanount栏中选择Fullexpansion,点击OK,得到如图6-13所示的整体图。

图6-13整体图

c.在Selectexpansionanount栏中选择3/4expansion，点击OK，得到如图6-14所示的3/4整体图。

图6-143/4整体图

7.结论

通过上述位移图分析可得，身管X方向的最大位移为0.011732mm,通过上述应力图分析可得，最大应力为298MPa,小于身管的屈服极限620MPa.，身管满足强度要求。

ANSYS有限元软件包是一个多用途的有限元法计算机设计程序，可以用来求解结构、流体、电力、电磁场及碰撞等问题，但是建模功能方面不是很强大。

因此，一般都是在AUTOCAD中建模，然后再导入ANSYS分析。

8.参考文献

[1]徐次达,华伯浩.固体力学有限元理论、方法及程序.北京：

水利电力出版社，1989.

[2]赵均海、汪梦甫.弹性力学及有限元.第二版.武汉：

武汉理工大学出版社，2008.

[3]丁茹、李丽丽.工程制图AtuoCAD实训教程.北京：

北京邮电大学出版社，2011.

[4]刘勇.大学计算机基础.第二版.北京：

清华大学出版社，2013.

[5]姜弘道,陈和群.有限单元法的程序设计.北京：

水利电力出版社，1983.

[6]胡仁喜，徐东什，李亚东.ANSYS13.0机械与结构有限元分析从入门到精通.北京：

机械工业出版社，2011.