圆柱体压缩过程模拟

1、.XX学院课程实验报告实验名称 圆柱体压缩过程模拟 实验课程 材料成型计算机模拟 指导教师 . 专业班级 .姓 名 .学 号 .2021年04月23日实验一 圆柱体压缩过程模拟 1 实验目的与内容1.1 实验目的进一步熟悉AUTOCAD或PRO/E实体三维造型方法与技艺，掌握DEFORM软件的前处理、后处理的操作方法与热能，学会运用DEFORM软件分析压缩变形的变形力学问题。1.2 实验内容砧板工件锤头图1 圆柱体压缩过程模拟运用DEFORM模拟如图1所示的圆柱坯压缩过程。一压缩条件与参数锤头与砧板：尺寸20020020mm，材质DIN-D5-1U,COLD，温度室温。工件：材质DIN_CuZ

2、n40Pb2，尺寸如表1所示，温度室温。表1 实验参数序号圆柱体直径，mm圆柱体高度，mm摩擦系数，滑动摩擦锤头运动速度，mm/s压缩程度，%1100150012021001500.21203100250012041002500.2120二实验要求1运用AUTOCAD或PRO/e绘制各模具部件及棒料的三维造型，以stl格式输出；2设计模拟控制参数；3DEFORM前处理与运算参考指导书；4DEFORM后处理，观察圆柱体压缩变形过程，载荷曲线图，通过轴对称剖分观察圆柱体内部应力、应变及损伤值分布状态；5比较方案1与2、3与4、1与3和2与4的模拟结果，找出圆柱体变形后的形状差异，说明原因；6提交分

3、析报告纸质和电子版、模拟数据文件、日志文件。2 实验过程2.1工模具及工件的三维造型根据实验给定的几何尺寸，运用PRO/E分别绘制坯料、锤头和砧板的几何实体，文件名称分别为workpiece，topdie，bottomdie，输出stl格式。 2.2 压缩过程模拟2.2.1 前处理建立新问题：程序DEFORM6.1FileNew Problem Next在Problem Name栏中填写“Forging Finish进入前前处理界面；单位制度选择：点击Simulation Control按钮，点击Main，在Units栏中选中SI，并在Simulation Title栏中填入“Forging，

4、勾选“Defromation。添加对象：点击+按钮添加对象，依次为“workpiece、“topdie、“bottomdie。定义对象材料：在对象树上选择workpiece，点击General按钮设置物体类型为Plastic，填入温度20；点击Meterial，选择DIN-CuZn40Pb2完成材料属性的添加。导入毛坯：导入已经画好的workpiece的几何实体；对几何实体进展几何检查，质量符合的图形才能进展网格划分并计算，检查质量合格。实体网格化：在对象树上选择workpiece，点击 Mesh，选择Detailed Settings选项卡，将类型Type改为绝对的Absolute，尺寸比S

5、ize Ratio改为2，最小单元尺寸Min Element改为3，然后生成实体网格。工件体积补偿：在对象树上选择workpiece，点击Property，在Target Volume卡上选中Active in FEM + meshing，自动生成体积补偿。用样的方法定义上模topie，物体类型为Rigid，温度为20，勾选Primary Die。设置上模参数，进入运动参数窗口，设置参数Directiont为-Z，Speed为1mm/s。定义下模bottomdie，与上模方法一样，但不勾选Primary Die。模拟控制设置：点击Simulation Control按钮，点击Step按钮在Nu

6、mber of Simulation Steps填入模拟步数为30计算得出的Step Increment to Save设置存储增量为3在Primary Die栏中选择topie在With DieDisplacement栏中填入步长为1点击OK按钮完成模拟设置。调整对象位置关系：在工具栏点击Object Positioning按钮进入对象位置关系调整对话框改变X、Y、Z的数值来调整位置点击OK按钮完成。边界条件定义：在工具栏上点击Inter-Object按钮在对话框上选择workpiecetopdie点击Edit按钮点击Deformation卡Friction栏上选中Shear和Constan

7、t选项，填入摩擦系数方案1、3和3、4组实验的摩擦系数 点击Close按钮如此重复，依次设置其它接触关系 点击Generate all按钮点击tolerace 按钮点击OK按钮完成边界条件设置。保存k文件：在对象树上选择workpiece点击Save按钮点击保存按钮保存工件的前处理信息重复操作，依次保存各工模具的信息。2.2.2 生成库文件在工具栏上点击Database generation按钮 在Type栏选中New选项选择路径英文填入数据库文件名英文，如Forging 点击Check按钮没有错误信息那么点击Generate按钮完成模拟数据库的生成。2.2.3 退出前处理程序在工具栏上点击E

8、xi按钮，退出前处理程序界面。2.2.4 模拟运算在主控程序界面上，单击工程栏中的forging.DB文件单击Run按钮，进入运算对话框单击Start按钮开场运算单击Stop按钮停顿运算单击Summary，Preview，Message，Log按钮可以观察模拟运算情况。2.3 后处理模拟运算完毕后，在主控界面上单击forging.DB文件在Post Processor栏中单击DEFORM-3D Post按钮，进入后处理界面。1） 观察变形过程：点击播放按钮查看成型过程；2） 观察温度变化：在状态变量的下拉菜单中选择Temperature，点击播放按钮查看成型过程中温度变化情况；3） 观察最大应

9、力分布：在状态变量的下拉菜单中选择Max Stress，点击播放按钮查看成型过程中最大应力分布及其变化情况；4） 观察最大应变分布：在状态变量的下拉菜单中选择Max Strain，点击播放按钮查看成型过程中最大应变分布及其变化情况；5） 观察破坏系数分布：在状态变量的下拉菜单中选择Damage，点击播放按钮查看成型过程中可能产生破坏的情况；6） 成型过程载荷：点击Load Stroke按钮，生成变形工具加载曲线图，保存图形文件为load.png；7） 点跟踪分析：点击Point Tracking按钮，根据上图点的位置，在工件上依次点击生成跟踪点，点击Save按钮，生成跟踪信息，观察跟踪点的最大

10、应力、最大应变、温度、破坏系数，保存相应的曲线图。 分别做方案1、2、3和4组实验，在后处理中提取必要的数据、曲线图等信息。3 实验结果与分析根据以上实验过程模拟出来的数据和图形及相关分析如下：3.1观察工件变形过程下面是四组方案的工件变形情况，分别取模拟步数为第1、9、18和30为研究对象，如以下图1、2、3和4所示；坯料压缩后呈单鼓形，由于坯料上下端受到锤头和砧板的摩擦力的作用，横向流动小，而坯料中间局部是自由变形，所以压缩后呈单鼓形。图1 方案1工件变形 坯料高度为150mm,摩擦力为0图2 方案2工件变形 坯料高度为150mm，摩擦力为0.2图3 方案3工件变形 坯料高度为250mm，

11、摩擦力为0图4 方案4工件变形 坯料高度为250mm，摩擦力为0.23.2观察温度变化后处理中观察温度的变化云图如以下图5所示,点击播放按钮选择第30步查看温度变化情况。在整个挤压过程中，坯料整体的温度根本没有变化且均为室温，因为在室温20条件下且工件中心不与空气和挤压模具接触，散热和传热都也没有，故温度没有变化，即方案1、2、3和4温度变化为零。方案2 温度变化方案1 温度变化方案4 温度变化方案3 温度变化图5 四组方案的温度变化云图比照3.3观察最大应力分布后处理中观察最大应力云图如以下图6所示,点击播放按钮选择第30步查看最大应力分布情况。比较方案1与2、3与4，无摩擦时最大应力分布均

12、匀而有摩擦时最大应力分层分布且最外层应力最大；这是因为有摩擦力时在圆柱外表出现附加拉应力使其应力发生变化，环向拉应力越靠近外层越大，而径向压应力越靠近外层越小。比较1与3、2与4，得知高度较短的坯料，其最大应力较大，这是因为压下量相对越小，产生的附加应力拉应力和附加压应力就越小，应力越小。方案2 最大应力方案1 最大应力方案4 温度变化方案3 最大应力图6 四组方案的最大应力云图比照3.4观察最大应变分布后处理中观察最大应变云图如以下图7所示,点击播放按钮选择第30步查看最大应力分布情况。比较方案1与2、3与4，无摩擦时的应变分布均匀而有摩擦时的应变呈分层分布且外层应变明显高于内层，这是因为无

13、摩擦时，应变分布比较均匀，圆柱体处于一向压缩，两向拉伸状态，为自由变形应变分布均匀；又摩擦时圆柱分为三个变形区，为难变形区、自由变形区和易变形区，所以应变分层分布，其外层应力较大。比较方案1与3、2与4，压缩量相对较小时其应变较大，因为压缩相对较小时其变形小，进而产生的最大应变较小。方案1 最大应变方案2 最大应变方案4 最大应变方案3 最大应变图7 四组方案的最大应变云图比照3.5观察破坏系数分布后处理中观察破坏系数云图如以下图8所示,点击播放按钮选中第30步的破坏系数分布曲线，四组方案中的破坏洗漱根本为零，也就是说四组方案的坯料破坏很小；这是因为我们压缩量较小，变形小，坯料的不均匀变化小，

14、变形抗力小。所以几组方案的破坏情况根本一样。方案2 破坏系数方案1 破坏系数方案4 破坏系数方案3 破坏系数图8 四组方案的破坏系数云图比照3.6成型过程载荷曲线点击Load Stroke按钮，生成变形工具加载曲线图，如以下图9和载荷曲线和图10载荷曲线比照图所示。在开场极短的时间内曲线呈垂直上升，此段属于弹性变形，变形均匀，载荷均匀增加；而后阶段呈缓慢的曲线上升，此段属于塑性变形，变形不均匀，载荷不均匀增加，随着挤压过程的进展，产生加工硬化，变形抗力增加，挤压力越来越大。分析四组方案，比较方案1与2、3与4，得知在有摩擦力的情况下需要的负载大，那么所需的挤压力大，这是因为有摩擦力时，金属流动

15、越来越困难，所需的挤压力也越大；比较1与3、2与4，得知高度较短的坯料，需要的载荷大，那么所需的挤压力大。方案2 载荷曲线方案1 载荷曲线方案4 载荷曲线载荷曲线方案3 载荷曲线图9 四组方案的载荷曲线比照根据上面四方案导出的相关成型过程载荷数值并制成相应图表如下：Timesec方案一方案二方案三方案四LoadN)0 00001 162890.516632495683.8998617.362 164586.8168491.296457.2899364.913 167010.3170824.397113.7100224.74 169252.3173189.298023.85100984.95 1

16、71569.117560298823.48101829.66 173941.3178067.699631.9102695.77 176360.3180585.5102142.6103538.58 178832183157.8102977.3104406.49 181355.5185785.7103823.1105276.810 183931.5188470.7104680.9106154.611 186565.6191214.1105551.110705312 189255.4194017.4106433.910795913 192003.5196883.1107326.3108875.914

17、 194812.6199812.9108228.6109804.715 197684.4202808.9109140.7110745.516 200620.2205873.1108269.7111698.717 203622209006.6109213.1112664.218 206691.4212212.2110157.2113642.219 209830.1215492.611111811463320 213041.1218849.2112098.9115637.321 216325.8222285113087.3116654.422 219688.1225802.4114088.9117

18、685.123 223128229403.7115103.7118729.624 226651.2233091.9116132119787.825 230257.6236869.8117173.9120860.426 233950.4240740.5118229.9121947.527 237732.7244706.7119298.8123049.228 241606.8248771.7120383.6124165.929 245578.5253123.2121483125297.930 249648257405.5122591.7126445.431 123721.3127608.832 1

19、24866.6128788.233 126027.5129984.534 127201.3131197.135 128402.5132426.736 129613133673.737 130840.5134938.438 132085.6136220.639 133339.9137521.640 134617138841.441 135915.3140180.342 137233.2141538.643 138570.9142916.844 139934.4144315.245 141311.5145734.346 142697.2147174.447 144115.314863648 145

20、553.615011949 147014.4151624.850 148504.3153153.2图10 方案1、2、3和4载荷曲线比较图3.7点追踪最大应力后处理中在状态变量的下拉菜单中选择Max Stress,生成最大应力变化曲线。从以下图11中可以看出，刚开场挤压时的应力最大，然后呈垂直直线下降，而后呈缓慢曲线下降，整个锻压过程中应力呈下降趋势。比较方案1与2、3与4，有摩擦力的情况下最大应力较大；比较1与3、2与4，得知高度较短的坯料，其最大应力较大。从图中可以看出方案3曲线图波动较大，这是因为网格划分不均匀，每个节点处的受力不均，导致坯料出现横向滑移现象。方案一最大应力方案二最大应力方案2 最大应力方案1 最大应力图11四组方案的最大应力比照方案4最大应力方案3 最大应力3.8点追踪最大应变曲线后处理中在状态变量的下拉菜单中选择Max Strain,生成最大应力变化曲线。从图12中可以看出，四个方案的最大应变都是呈线性增长的趋势。比较方案1与2、3与4，相对无摩擦的情况下有摩擦力时最大应变较大。比较方案1与3、2与4，在无摩擦的情况下坯料上三点的最大应变曲线根本重合，而在有摩擦的情况下坯料上三点的最大应变曲线较分散，且中间P2点最大应变相对大一些；在压缩量一样时有摩擦的最大应变相对较小，这是因为有摩擦时金属流动较困难，变