黄铜压型有限元模拟分析j.docx

1、黄铜压型有限元模拟分析j目录一、题目描述2二、题目分析2三、解题模拟、思路2四、模拟过程 1、定义工作文件名32、定义单元类型、材料参数和实常数33、创建模型54、划分网格105、创建接触对116、施加约束、载荷及求解137、显示结果17五、模拟结果分析24六、结论24七、参考文献25一、 题目描述： 如图1.1所示为金属紫铜坯料和挤压模具结构示意图，坯料与模具之间的摩擦系数为0.15。求挤压过程中坯料内部的应力场变化、应变场变化。坯料紫铜的材料参数：弹性模量：E1=1.08e11MP；泊松比：1=0.33；密度：1=8.9e3kg/m3。模具材料参数:弹性模量：E2=2.07e11MP；泊松

2、比：2=0.30；密度：2=7.9e3kg/m3。 图 1.1 所要分析模型的剖面图二、 题目分析：该问题属于状态非线性大变形静态接触分析问题，因此可以选择对其做静态分析。在分析过程中根据结构的轴对称性，只要求分析整个模型的1/4，然后将其延展成一个完整的模型。此外，需要注意模拟分析过程中的单位统一。三、 解题模拟、思路：1、 选择单元类型、定义材料属性和实常数；2、 创建具体模型；3、 网格划分同时分配材料属性，创建接触队；4、 施加约束及载荷；5、 求解、查看最终分析结果。四、 模拟过程：1、定义工作文件名和工作标题： 1.1、定义工作文件名 执行Utility Menu-FileChan

3、g Jobname-20090771，单击OK按钮。1.2、定义工作标题 执行Utility Menu-FileChange Tile-yuyuyang20090771，单击OK按钮。1.3、更改目录 执行Utility Menu-File change the working directory D/ansys。2、定义单元类型和材料属性以及实常数： 2.1、 设置计算类型 ANSYS Main Menu: Preferences select Structural OK，如图2.1. 图2.12.2、选择单元类型执行ANSYS Main MenuPreprocessor Element T

4、ypeAdd/Edit/Delete Add select Solid Brick 8node 185 单击Apply；重新在library of element types第一个列表中选择contact,在第二个列表框中选择 3d target 170，再次单击Apply，重新在library of element types第一个列表框中选择contact,在第二个列表框中选择4 nd surf 173。如图2.2所示，选择OK接受单元类型并关闭对话框。 图2.2定义单元类型对话框 2.3、定义实常数选择main menupreprocessorreal constantsadd/edit

5、/delete,出现Real constants对话框，单击add按钮，在set 1下点击add添加type 3contact 173下点OK ,在real constant set no.中输入1，单击ok;继续添加点击add,按照同样的方式定义另一组选择type 2targe170的实常数，选择编号2 。 图2.3 定义实常数2.4、定义材料属性 、执行Main MenuPreprocessor Material Props Material Models命令，出现Define Material Model Behavior对话框。 、在Material Models Available

6、列表框中依次双击StructuralLinearElasticIsotropic选项，出现Linear Isotropic Material Properties for Material Number 1对话框。在EX输入栏中输入1.08E11，在PRXY输入栏中输入0.33，单击“OK”按钮关闭对话框，如图2.4(a)。 、依次单击StructuralFriction Coefficient选项，出现Friction Coefficient for Material Number 1对话框，在MU输入栏输入0.15，单机“OK”按钮关闭对话框，如图2.4（b）。 、依次单击 Structu

7、ralDensity选项，出现 Density for Material Number 1对话框，DENS输入栏输入8.9E3，单机“OK”按钮关闭对话框，如图2.4（c）。 、在Define Material Model Behavior对话框中单击MaterialNew Model命令，出现Define Material ID对话框，在输入栏中输入2，单机“OK”按钮关闭对话框。 、在此时，将重复操作、。只需分别将EX、PRXY、MU、DENS输入栏中的数值输入2.06E11,0.30,0.15,7.9E3，分别单击“OK”按钮关闭对话框。 、在Define Material Model

8、Behavior对话框中选择MaterialExit命令，关闭该对话框。(a) （b） （c） 图2.4输入材料性能参数对话框3、创建几何模型： 运用自底向上建立模拟模型，具体操作步骤如下：3.1、定义关键点 选择ANSYS Main Menu: Preprocessor Modeling Create Keypoints In Active CS 依次输入五个点的坐标: input: 1.(0.017,0,0)、2.（0.055,0,0）、3(0.055,0.110,0)、4（0.021,0.110,0）、5（0.021,0.040,0）。 创建第六关键点，a）、先创建局部坐标到第五点，从实

9、用菜单中选择Utility Menu：WorkPlane Local Coordinate SystemsCreate Local CSAt Speciffied Loc +。 b）、在打开的创建坐标系对话框中，在Global Cartesian 文本框中输入0.021,0.040,0然后点击“OK”，得到Create Local CS At Specified Location对话框，如图3.4。 c）、在Ref number of new coord sys 文本框中输入11，在Type of coordinator system中选择Cylindrical 1，直接单击OK退出，如图3.

10、5。 d）、ANSYS Main Menu: Preprocessor Modeling Create Keypoints In Active CS输入（0.02，-125,0）OK结束。 图3.4输入K5点坐标 图3.5创建局部坐标系 、将坐标系转到全局直角坐标系下，从实用菜单中选择Utility Menu:WorkPlaneChange Active CS toGlobal Cartesian。现在创建第七关键点坐标为（0.017,0.040,0）。3.2、生成模具截面： 、ANSYS 命令菜单栏: Main Menu: Preprocessor Modeling Create Lines

11、 In Active Coord 依次连接点1,2,3,4,5,6,7OK ，生成线L 1,2,3,4,5,6如图3.6。、将线5,6搭接：Main Menu: Preprocessor Modeling OperateBooleansOverlapLines输入5,6OK，搭接完之后生成线7,8,9,10，操作完结果如图3.7（a）。 、将线9,10删除：Main MenuPreprocessorModelingDeleteLinesOnly，分别输入线9和线10，打击OK,操作完结果如图3.7（b）。 图 3.6 生成直线 （a）线搭接后的图形 （b）线删除之后的图形 图3.7、对线4和线

12、8倒圆角：Preprocessor Modeling Create Lines Line Fillet，在输入栏输入，线3和线8，圆角半径0.005，打击OK结束。如图3.8所示。 图3.8、通过线生成截面：Preprocessor Modeling Create Areas Arbitrary By Lines拾取所有的直线OK。 3.3、创建模具和坯料1/4体：、创建一条旋转直线，由关键点K 100（0,0,0）,110（0,0.110,0），生成直线L6。、选 装90：Main MenuPreprocessorModelingOperateExtrudeAreasAbout Axis拾取

13、面1OK在对话框中输入100,110OK再在对话框中ARC项中输入90，如图3.9所示。其生成体的效果如图3.10所示。 图3.9 图3.10 经旋转90后效果图、生成坯料模型体：a）、先生成坯料截面：Preprocessor Modeling Create Areas RectangleBy Dimensions在对话框中输入（0.019,0.021）、（0.042,0.102）OK，如图3.11所示。b）、接下来在按着操作即可得到坯料的1/4体，其效果如图3.12所示。 图 3.11 图 3.12 1/4体的模型体效果图4、划分网格： 选用MainMenuPreprocessor Mesh

14、ing Mesh ToolVolumes MeshTetFree,.采用自由网格划分单元。材料属性分配，将1号材料分配给坯料，将2号材料分配给模具；以及实常数，将刚性单元赋予模具，将柔性单元赋予坯料。在SmartSize前挑勾后可以自由设置网格的密度，在此将坯料的网格密度设置为4，再将模具的密度设置成默认状态6，如图3.13所示。最终划分完网格的效果图如3.14。 图 3.13 划分网格前的设置 图3.14 网格划分完后效果图5、创建接触对:在Main Menu Preprocessor Modeling Create Contact Pair下湖出现一个Contact Manager对话框如

15、图3.15所示。点击Pick Target会出现图3.16在对话框中输入3,4,5,6点击OK后返回图3.15的界面点击Next后再输入12点击OK，在点击Next会出现一个如图3.17的窗口点击Create即完成，生成的接触面效果如图3.18。 图 3.15 图3.16 图3.18 图 3.176、施加约束、加载求解: 6.1、定义分析类型: 选择Main MenuSolutionAnalysis Type命令，出现New Analysis对话框，选择Static，单击OK按钮关闭对话框。 6.2、施加约束：选择Main MenuSolutionDefine Laods ApplyStruc

16、turalDisplacement a）、再继续选择SymmetryOn Areas，现在拾取面1,9，10,15面点击OK，即对对称面施加约束结束，如图3.18所示。 b）、再继续选择On Areas出现拾取对话框，拾取面2,4,5,6,7然后点击OK，出现图3.19。 6.3、施加第一个载荷步及求解： 、Main MenuSolutionAnalysis TypeSols Controls，会出现如图3.21对话框，按图样输入，点击OK退出对话框。、Main MenuSolutionSolveCurrent LS，出现solve current load step对话框，同时出现statu

17、s comman窗口，选择file/close命令，关闭窗口。再单击solve current load step对话框中OK按钮，第一个载荷步下开始计算求解。、求解结束时，出现note提示框，单击close按钮关闭该提示框，显示迭代过程的时间跟踪图。 6.4、施加第二个载荷步及求解：、Main MenuSolutionDefine Laods ApplyStructuralDisplacementOn Areas 出现Apply U,ROT on lines对话框，先拾取13面，在选择UY在Value displacement volue文本框中输入0.030施加向下的位移，单击OK按钮关闭

18、对话框。、选择utility menuselecteverything命令，选择所有的点、线、面。、同6.3第一步一样，首先将Time attend of loadstep中输入100，其次将Automatic time stepping 中off改成on，最后再在下面三个输入栏输入100、1000、10，单击OK结束。、操作和6.3第三步一样，这样施加约束和加载过程结束。 图3.18 施加对称约束后 图3.19 施加面上约束 图3.20 施加完约束后的图样 图3.21 施加第一个载荷步时的对话框 图3.22施加完第二载荷求解后的结果图7、显示结果： 7.1.显示变形形状。执行Main Men

19、uGeneral PosprocPlot ResultsDeformed Shape，弹出如图3.23所示的对话框。选择“Def+underformed”单选按钮，单击OK按钮。生成结果如图3.24所示。 图 3.23 图 3.24 变形形状图 7.2、列出节点的总应力结果值。执行Main MenuGeneral PosprocList ResultsNodal Solution，弹出如图3.25所示的对话框。设置好后点击OK按钮。生成如图3.26所示的结果。依次可以调试出单元应力或应变值。 图 3.25 图 3.26 部分节点应力值 7.3、浏览节点上的Von Mises应力值和应变值。执行

20、Main MenuGeneral PosprocPlot ResultsContour PlotNodal Solu，弹出如图一个设置对话框。设置好后单击OK按钮，生成结果如图3.27所示。 3.27-a 节点的应力浏览图 3.27-b 节点应变浏览图7.4、浏览单元上的Von Mises应力值和应变值。执行Main MenuGeneral PosprocPlot ResultsContour PlotElement Solu，弹出如图一个设置对话框。设置好后单击OK按钮，生成结果如图3.28所示。 图 3.28-a 单元应力浏览图 图3.28-b 单元应变浏览图7.5、总的变形结果动画图的展

21、示：操作过程为Utility MenuPlotctrlsAminateDeformed Result出现对话框直接点击OK结束，最终效果由以下几张图展示如图3.29。 图 3.291图 3.292 图 3.293 图 3.294 图 3.295五、模拟结果分析：1、节点的应力应变分析：由3.27两图可知节点的最大应力为3.36E10Pa，最大应变为3147E-3m。2、单元的应力应变分析：由3.28两图可知单元的最大应力为 3.52E10Pa，最大应变为 3525E-3m。 通过以上分析可知在挤压过程当中受到最大的应力集中在角度为35的斜面部分为最大，同时在此时应变也比较大。六、结论：通过这一

22、次的有限元模拟分析，我们从中了解到如何对一个实体进行模拟的具体过程。1、在此次模拟过程当中我觉得最容易遇见的问题有以下两个问题及相应的解决方法：1.1、在建模过程中会经常出现错误，但是在Ansys软件当中没有撤销这么一步，往往做到这么不知道怎么修改。我觉得在此时的解决方法是养成做完一步保存一步的好习惯，当遇见错误的时候可以调出保存好的文档。1.2、在拾取相应的线或则面的时候不太容易直接拾取。这是应该把模型放大或则是将其旋转当最佳位置（快捷方式Ctrl+鼠标左、右键）。2、在最后模拟的分析结果的时候，通过几次不同的验证都没有解决掉的是挤压铜坯料的短的压下距离以及坯料为什么插入模具，在此过程当中通

23、过了以下几个方法：2.1、将模具分成四部分，分别在相应关键点处创建工作平面，运用布尔运算划分模具，之后在进行同样的操作。2.2、在划分网格时选用手动划分网格，将坯料和模具的网格密度的差别具体的体现出来。 3、通过认真的对待这次模拟分析，我们从对Ansys的零基础到现在可以熟练操作。其次，在Ansys中做模拟分析的过程当中体会到有些操作用菜单操作是比较费劲的，有的时候可以结合一些命令来操作，如：删线、面、体直接可以用LDELE、ADELE、TDELE操作。如果以后要完全熟练掌握Ansys的运用首先必须学会一些常用的命令流，再有就是要把两者有力的结合起来一起运用。七、参考文献： 1、王勖成，邵敏。有限单元法基本原理和数值方法（第二版）。北京：清华大学出版社，2003 2、 李尚建主编。金属塑性成形过程模拟。北京：机械工业出版社，1999 3、 俞汉青，陈金德主编。金属塑性成形原理。北京：机械工业出版社，1998 4、王从曾主编。材料性能学。北京：北京大学工业出版社，2001 5、刘惟信主编。机械最优化设计（第二版）。北京：清华大学出版社，1994