金属塑性成形原理大作业结题报告

金属塑性成形原理大作业结题报告拉深成形模拟（ABAQUS）在本次大作业的研究中，借助有限元分析成形模拟软件ABAQUS,就筒形件的拉深为例实现对拉深成型过程的模拟仿真并将该过程详细的记录下来，分析凸凹模圆角半径的大小、凸凹模间隙、有无压边圈和各部分的摩擦力等因素对拉深成形效果的影响。2014班级学号姓名目录一对象描述1二研究目标1三研究内容1四研究方案1五模拟之前的准备工作2六分析模拟过程361制作模型3611启动软件，建立工程3612绘制零件362确定材料参数并定义和赋予截面属性663定义装配件764定义分析步和输出变量8641定义分析步8642定义场输出和历程输出865定义接触和摩擦966施加载荷和边界条件967网格划分1068生成作业11七后处理数据图表分析1171观察毛坯在拉深过程中的形状变化，分析其各主要位置的受力状态1172拉深变形起皱的情况，找到破裂的危险界面位置1373材料、磨具几何尺寸、压边力等参数对拉深成形能力的影响1674主要四部分的摩擦对拉深成形效果的影响32八工艺参数对成形结果的影响规律37九与其他计算方法的比较38十学习总结38十一参考文献39一对象描述拉深是指将平板毛坯或杯形毛坯在凸模作用下拉入凹模型腔形成开口空心零件的成形工艺方法。拉深也称为压延，是钣金成形的基础性工艺。拉深有多种形式。按零件的外形可分为筒形件、锥形件、半球形件、阶梯型件、盒形件和复杂形状零件拉深等等。而在各种拉深成形工艺中，圆筒形件拉深是最基本的拉深方法。故在本次大作业的研究中，我将借助有限元分析成形模拟软件ABAQUS,就筒形件的拉深为例来实现对拉深成形过程的仿真模拟。并将该过程详细的记录下来，加深对课程知识的理解和掌握，提高自己动手解决工程问题的能力，形成具有借鉴和学习意义的作业报告。二研究目标利用有限元分析软件ABAQUS模拟筒形件的拉深变形。研究起皱的条件和防皱的措施；找到破裂的危险截面的位置；探究改变凸凹模圆角半径、凸凹模间隙、压边力和板料相对厚度等条件对拉深性的影响。三研究内容1观察毛坯在拉深过程中的形状变化，分析其各主要位置的受力状态。2拉深变形起皱的情况，找到破裂的危险界面位置。3改变凸、凹模圆角半径、凸凹膜间隙、板料相对厚度、有无压边圈、弹性压边圈压边力、拉深变形速度、板料的材料观察改变不同的条件对拉深性的影响。4研究压边圈与板料的摩擦、凹模上表面与板料的摩擦、凸模圆角处与板料的摩擦、凹模圆角处与板料摩擦对拉深成形的影响。四研究方案1完成拉深过程的模拟，根据毛坯的变形情况和应力云图并结合已学过的应力应变知识完成研究内容第一项。2对比存在压边圈和不存在压边圈的拉深效果的区别并结合研究内容第三项来完成研究内容第二项。3保持其他条件不变，改变凹模圆角半径、凸模圆角半径、凸凹膜间隙、板料相对厚度，观察改变不同的条件对拉深性的影响，并记录下研究结果来完成研究内容第四项。4保持其他摩擦条件不变，单独改变压边圈与板料的摩擦、凹模上表面与板料的摩擦等摩擦条件来探究摩擦对拉深成形的影响。五模拟之前的准备工作1通过查阅文献查得本模拟实验的材料深冲钢DC06的力学性能如下表表51DC06材料属性参数（力学）屈服强度/MPA抗拉强度/PA伸长率/数值16025053DC06极限拉深系数LIM0432拉深件的成品尺寸设计高度；根据高度的值取修边余量；H252M当量高度；27H中径尺寸；凸模圆角半径。41D9PR毛坯直径D的计算公式如下22RHR2DR8741D取毛配直径为进行拉深。74M故可以一次拉深成型。LI1/050433根据毛坯尺寸和成品的设计尺寸计算磨具和压边圈的尺寸参数如下其中，材料的参数和磨具板料的尺寸如下所示材料、磨具、参数的初始取值1）毛坯的尺寸直径、厚度；材料深冲钢DC06；074DM1T2）压边圈的尺寸内径、外径；解析刚体；125285DM3）凸模的尺寸长度、圆角半径、直径；PL9PR40PDM解析刚体；4）凹模的尺寸长度、圆角半径、内径50DLM105D、外径；解析刚体；142285）凸凹模之间的双边间隙值；Z6）采用弹性压边圈，最小压边集中力的数值计算公式经过计算得；S00D28T1765BDQR287QKN7）弹性模量；泊松比；密度；EGPA335/CMG在以上的各种计算结果中，凸凹模的圆角半径、内径均为实验室量得的实际数据；凹模的外径为通过毛坯的尺寸经过计算而得；压边圈和压边圈上的集中力是通过书本上的公式经过严格计算所得。4ABAQUS的单位制ABAQUS没有统一的单位各单位之间只要达到统一即可六分析模拟过程61制作模型611启动软件，建立工程图61进入界面612绘制零件如下图所示依次建立各部件，下面以DIE的创建为例，展示模型的建立过程。图63模型建立121从工具栏的MODULE列表中选择PART项进入部件PART模块，点击PARTCREATE来创建新的部件。PART的定义如下图所示122运用直线、圆角工具实现草图的绘制。图63草图绘制123绘制完草图后点击鼠标中键，旋转角度选择360，完成DIE的特征生成。之后选择刚体参考点（刚体必须选择参考点）。方法TOOLSREFERENCEPOINT选择刚体上一点即可。图64旋转体设置REVOLEDIRECTION用来控制旋转的方向。若想要更改草图，在左边的模型树下选择DIE一直到SECTIONSKETCH子单元，右击选择EDIT即可。更改完成后右击FEATURES选择REGENERATE更新模型。图65凹模旋转成形本模型中拟建立的模型有四个DIE凹模、PUNCH凸模、BLANK板料、HOLDER压边圈。其中DIE、PUNCH、HOLDER为离散刚体（DISCRETERIGID）、BLANK为可变形体（DEFORMABLE）。图66建立刚体参考点不定义刚体的参考点程序运行时会报错。对刚体的操作实际上就是对参考点的操作，后面的定义边界条件时就是对参考点的约束，选择刚体的整体来约束即达不到想要的效果又会造成程序执行的冗余。62确定材料参数并定义和赋予截面属性点击创建材料。在此模块下定义了材料四个属性如下密度（DENSITY）、弹性（ELASTIC）、塑性（PLASTIC）、热膨胀率（EXPANSION）。本次模拟使用深冲钢DC06此种材料的力学性能定义如下表所示1DENSITY、ELASTIC、EXPANSION的设定表61泊松比、杨氏模量、密度、热膨胀系数的取值MASSDENSITY785E009POISSONSRATIO03YOUNGSMODULUS210000EXPANSION135E005（表格中的数据均为ABAQUS中设定的原始数据，此处没有考虑单位）2PLASTIC的设定表62DC06塑性应力应变值YIELDSTRESS160170180190200PLASTICSTRAIN001699033980509706796YIELDSTRESS210220230240250PLASTICSTRAIN08495101951189413593153（表格中的数据均为ABAQUS中设定的原始数据，此处没有考虑单位）对于动力问题必须设置MASSDENSITY。图67材料属性的赋值63定义装配件图68装配体效果图从工具栏的MODULE列表中选择进入装配体模块。定义装配件时为了完成各PART之间的顺利装配，在菜单下选择共轴约束为了正确地定义各PART之间的上下相对位置应用了菜单下的工具。使用工具移动各装配部件完成装配。定义完成的装配体如上图所示64定义分析步和输出变量641定义分析步程序中定义了五个分析步（原则上只定义两个分析步即可，但通过查阅书籍了解到将分析步细化有利于计算过程的收敛，故将分析步定义为六步）。第一步实现压边圈压向板料，第二步实现凸模下行知道与板料接触。第二步以后继续凸模的下行实现板料的拉深过程，直到筒形件的拉深过程结束。在分析步中可以选择生成自适应网格，自适应网格方法是指计算中，在某些变化较为剧烈的区域，如大变形、激波面、接触间断面和滑移面等，网格在迭代过程不断调节，将网格细化，做到网格点分布与物理解的耦合，从而提高解的精度和分辨率的一种技术。自适应网格希望在物理解变动较大的区域网格自动密集，而在物理解变化平缓区域网格相对稀疏，这样在保持计算高效率的同时得到高精度的解。具体的设置步骤如下选择点击创建自适应网格。开始时，我使用自适应网格，但是在后处理的时候发现开始定义研究的集合在运行过程中会发生变化，不利于分析，故之后舍弃了这种方法。表63分析步的分配第一分析步第二分析步第三分析步第四分析步第五分析步第六分析步0101025025025025642定义场输出和历程输出在场输出变量FIELDOUTPUT中分别选择应力、应变、位移、力选项卡下面的若干项和厚度选项卡下的（用于检测板厚的变化）。65定义接触和摩擦定义四个摩擦因数分别为0、01、025、035。分别定义7对接触分别为凸模圆柱外表面与板料上表面的接触（初始摩擦系数0）、凸模圆角外表面与板料上表面的接触（初始摩擦系数01）、凸模底外表面与板料上表面的接触（初始摩擦系数0）、凹模面与板料下表面的接触（初始摩擦系数01）、凹模圆角处与板料下表面的接触（初始摩擦系数01）、凹模内腔表面与板料下表面的接触（初始摩擦系数0）、压边圈下表面与板料上表面的接触（初始摩擦系数01）。定义完成的凸模圆角外表面与板料上表面的接触效果图如下图69定义接触66施加载荷和边界条件在弹性压边圈的模型中，程序中第二分析步将一集中力（大于30FKN最小压边力）施加于压边圈上，实现防皱的效果。在第一个分析步内实现压边圈与板料的接触，第二个分析步中实现凸模与板料的接触，在剩下的分析步中实现凸模的下行，完成拉深成形过程。第一分析步的速度是400MM/S、第二分析步的速度是500MM/S以后各步的速度是300MM/S。边界条件约束凹模的全部自由度将其固定；对于弹性压边圈放开其Y轴方向的移动自由度（必须定义刚体的转动惯量），对于固定压边圈则固定其全部的自由度；凸模Y轴方向的移动自由度不限制，给它一个向下的速度使其下行；不限制板料的任何自由度。67网格划分网格划分的粗细程度会影响最后的分析结果，网格过细会使运行的速度大大降低，网格果粗达不到需要的精度。网格的划分情况如下所示（以BLANK为例）图610设置网格属性图611网格划分ABAQUS中必须对解析刚体划分网格。在后来的实践中发现使用扫略网格会大大提高之后的后处理分析速度。故在之后的分析中部分采用了扫略网格。68生成作业完成各项设置之后点提交作业；点击实时监控JOB的执行情况；运行结束后点击查看运行之后的结果。图612提交作业七后处理数据图表分析71观察毛坯在拉深过程中的形状变化，分析其各主要位置的受力状态对变形过程进行分析并结合金属塑性成形理论知识可得出如下结论1）拉深过程中凸缘区材料不断拉入凹模型腔形成筒壁。凸缘区是拉深材料的主要变形区。凸缘材料的变形以压缩为主。2）凸模圆角的材料在拉深一开始就受到凸模的冲压，在拉深过程中受到双向拉应力和一定的弯曲作用。3）凹模圆角处材料受到径向拉深变形，同时产生塑性弯曲，进一步从凹模圆角区拉向筒壁时又要被较直，受到反向弯曲。4）圆筒侧壁受轴向拉伸。5）筒底部分的材料受平面拉伸作用，由于凸模圆角摩擦力的作用变形不大。6）由筒形件的拉深成形模拟过程可以看出毛坯主要受到四个摩擦力的作用，凸缘材料与压边圈之间、凸缘材料与凹模面之间、凹模圆角区与凹模圆角之间和凸模圆角区材料与凸模圆角之间的摩擦。毛坯在拉深过程中逐渐从圆形毛坯转变成了筒形件，变形结果如图71图71拉深成形效果图应力云图如下图所示图72拉深成形过程图72拉深变形起皱的情况，找到破裂的危险界面位置对于破裂的危险界面的探究根据模拟结果我认为应该分成两种情况来进行讨论1）凸缘部分材料面积较小时，模拟结果的应力云图如图72所示最大应力部位发生在凹模圆角处，凸缘部分的材料由凸缘区流入凹模圆角时在凹模圆角处发生弯曲和反弯曲变形使该处的应力增大，又由于凸缘部分的材料面积较小，板料所受的凹模面的摩擦较小，板料能较容易地流入凹模型腔形成筒壁，所以最大应力的部位发生在凹模圆角处。2）凸缘部分材料面积较大时，模拟结果的应力云图如图73所示最大应力部位发生在凸模圆角处，在凸模下行的过程中板料首先与凸模圆角发生接触，受到凸模的顶压和弯曲变形的影响厚度有一定的减薄，在拉深的过程中进一步收到双向拉深，厚度又有所减薄。又由于凸缘部分的材料面积较大，板料所受的凹模面的摩擦较大，板料能较难流入凹模型腔形成筒壁，通过凸模圆角传递的拉深力增大，所以最大应力的部位发生在凹模圆角处。由应力与应变的关系可知此处的也应变较大容易发生破坏。图72凸缘面积小时最大应力处图73凸缘面积大时最大应力处经测量图72凹模圆角处材料减薄019MM,凸模圆角处减薄012MM。图73凹模圆角处材料减薄011MM,凸模圆角处减薄023MM。由此可以验证以上的两条结论的正确性。再者当凸模圆角或凹模圆角较大时，板料也容易起內皱。凹模圆角较大时的情况如下（凹模圆角半径为9MM、板料相对厚度为1/74）图74凹模圆角较大起皱图凸模圆角较大时的情况如下图75凸模圆角较大起皱图使凸凹模的圆角半径限制在合理的范围内时可以限制起皱。加大板料的相对厚度也可以防止起皱。下图为凹模圆角半径为9MM/板料相对厚度为3/74时的情况，可以看出起皱的情况得到缓解。图76加大板料厚度变形图所以为了防止起皱可以选择板料的相对厚度大的材料（在允许的情况下）压边力对起皱的影响在后边会单独讨论，此处不再赘述。73材料、磨具几何尺寸、压边力等参数对拉深成形能力的影响在研究多个要素对结果的影响时常常采用控制因素（变量）的方法，把多因素的问题变成多个单因素的问题，而只改变其中的某一个因素，从而研究这个因素对事物影响，分别加以研究，最后再综合解决。在问题3的研究中，使用控制变量法，保证在其他影响因素不变的条件下对某一单一要素的取值情况进行调整，研究其对拉深变形的影响状况。但是在选择要研究的单元体时存在一定的困难，原因在于采用手工选择的方法精度不高，不能保证每次的选择都完全精确的显示在同一个单元上，误差较大。通过查阅书籍，我采用了如下的方法对单元进行研究显示各单元的编号，通过对比相同编号的单元应力状况，分析凸模圆角处的应力变化状态，得到相应的规律。显示单元编号的方法在VISUALIZATION板块选择选项下的选项，选择其下四个选项卡中的选项卡然后再选中然后单击。设置完的板料变化如图77（其中板料上方的蓝色的类似圆圈的物体就是单元的编号）图77显示网格编号731凹模圆角半径的大小对拉深成形的影响选取节点417，输出其在、下的应力应变曲线。5DR9D13DR图78应力图图79应变图表71最大应力与等效应变凹模圆角半径大小（MM）最大应力等效应变5208777（节点417）950905E039197321（节点564）663885E031319341（节点120）567641E03图710凸模圆角过大起皱图不难看出在凸模与板料接触时节点417所受的应力急剧增大，达到一定程度以后缓慢变化。在应变曲线中的最大应变远大于和时，且5DR9DR13D后两者的最大应变相差不大。随着凹模圆角半径的增大，最大应力应变呈下降的趋势。所以可以的出结论凹模圆角半径越大对拉深成形越有利。当圆角半径过大时板料悬空部分太多容易起內皱。732凸模圆角半径的大小对拉深成形的影响表72最大应力与等效应变凸模圆角半径大小（MM）最大应力等效应变520096（节点149）741671E03919941（节点3）717631E0313200595（节点572）732921E03选取节点149，输出其在、下的应力应变曲线。5DR9D1DR图711应力图图712最大应力变化图图713应变图由上图可以知随着凸模圆角半径的增大，板料上的最大应力变化不大所以可以认为在凸模圆角半径在的范围内其大小对拉深性的影响不明显。513M此处得出的结论与理论上不符，于是又补充了和的模拟实验，3DR15D模拟实验的应力结果下图所示从中仍然看出凸模尺寸增大，应力的变化范围还是不明显。结果还是与理论上有出入。通过查阅书籍发现我取的凸模圆角半径的值都在对应于板厚1MM的合理范围内，故可的出结论在凸模圆角半径在对应板厚的合理范围内时，凸模圆角半径的大小对板料的拉深性没有太大的影响。图714应变图733凸凹模间隙对拉深成形的影响使用工具建立通过原点的坐标参考面，再使用PARTITION对板料进行拆分操作，这样划分网格的时候可以在参考平面的两侧建立均匀的网格，方便拾取点来对不同位置的直径进行测量。图715剖分图网格划分之后的效果图图716剖分后网格划分图模型跑完之后的效果图（正对的是所建立的参考面附近的网格）图717剖分后模型运行结果图718C3时轻微起皱图在略低于零件口部的位置出现了轻微的內皱现象表73最大应力与等效应变凸凹模间隙C最大应力零件底部直径零件中部直径零件口部直径080MM232425400320411843414854100215958419988419578416707125MM19941424715421350416065150MM19664426445421298416337图719零件各部分直径与凸凹模间隙的关系图由以上两个图可以看出最大屈服应力随着凸凹模间隙的增大而减小，原因是随着凸凹模间隙的增大，凸凹模侧壁对板料的挤压作用相应减小，最大屈服应力降低。凸凹模间隙为08时，间隙小于板料的厚度，板料拉深过程中会发生变薄现象，在侧壁会发生应力集中。间隙过大易起皱。在凸凹模的间隙小于板料的厚度时，在凹模圆角处及以下部分发生减薄现象使零件的锥度增大。在凸凹模的间隙大于板料的厚度时，随着间隙值的增大，零件的外形尺寸变化的越来越大，锥度呈现逐渐增大的趋势。由零件各部分直径与凸凹模间隙的关系图可以直观的看出，当凸凹模间隙等于毛料的板厚时，零件的锥度最小。下图为各间隙值时零件的外形轮廓的对比08C125C10C150C图720不同间隙变形对比图734有无压边圈对拉深成形的影响在设定的板料和磨具尺寸的条件下去掉压边圈将会事板料形成不能去除的死皱，且在褶皱处产生应力集中的现象。起皱情况如下图所示图721无压边圈起皱图735改变拉深的材料对拉深成形的影响网上难以搜到深冲钢DC06和2324T351的屈服应力与塑性应变等塑性指标，但是找到了2324T351的屈服轨迹和DC06的抗拉强度和屈服极限的值。采用如下的方法近似求出2324T351和DC06的塑性指标。1根据2324T351的屈服轨迹在MATLAB中取一系列点模拟其塑性指标，模拟的曲线结果如下图所示，模拟结果与实际结果相类似。图7222324T2351真实应力应变曲线2根据深冲钢的弹性模量、抗拉强度、屈服极限和断后延伸率模拟其屈服轨迹。近似用直线代替深冲钢屈服轨迹中的屈服阶段。在模拟的直线段上取10个点来近似的代替其塑性指标。图723DC06真实应力应变曲线表74最大应力与等效应变拉深材料最大应力等效应变DC06（深冲钢）19664（节点572）638702E0345号钢920648（节点85）56968E032324T35153553（节点632）654799E03图724应变图图725应力图由图可以看出加载过程相同，不同材料同一个单元的应变状态没有多大的差异。但是就应力水平来说DC06的性能高于2324T351的性能高于45号钢的性能。原因是DC06的屈服极限远小于2324T351的屈服极限而2324T351的屈服极限又远小于45号钢的屈服极限。在同样的加载条件下，DC06能以很快的速度进入屈服阶段其次是2324T351故应力水平呈现为图示的状态。736拉深变形速度对拉深成形的影响图726应变图图727应力图表75最大应力与等效应变拉深变形的速度最大主应力等效应变300MM/S196945（节点631）635275E03500MM/S198724（节点625）657921E03800MM/S196289（节点612）619986E03由上面的表格可以看出三种加载速度所得的最大应力和等效应变差别甚微且达到最大应力的节点彼此之间相互靠近。在这种意义上来讲加载速度对拉深过程的影响不大。但就上边的应力应变图来说，他们达到最大应力和应变的时间却不相同。也就是说如果加载速度增大到凸缘材料不能及时流入凹模型腔的水平，那么在凸模下行的过程中板料就会通过减薄来实现筒壁的增长，这样对拉深过程是不利的。所以加载速度应该控制在合理的水平。737板料相对厚度对成形结果的影响图728应变图图729应力图表76最大应力与等效应变板料相对厚度最大应力等效应变1/74203727（节点498）74292E032/74197797（节点42）713772E033/7419281（节点1976）595194E03由以上可知板料的相对厚度越大，板料的拉深性能越好。最大主应力和等效应变越小。738压边力的大小对成形结果的影响图730应变图图731应力图由表中数据和应力应变图像可以看出压边力越大，最大应力和等效应变越大。原因是压边力越大引起的压边圈和板料之间和凹模面与板料之间的摩擦力越大，结果使筒壁承受更大的拉应力。由应力图还可以看出在加载的初始阶段不同压边力的应力曲线相互重合。原因是开始时压住板料所需的力较小，三种压边力均够用，而在加载的末尾阶段板料所需的压边力越来越大，小的压边力渐渐不能够满足要求，渐渐出现差别。图732采用弹性压边圈起皱图另外压边力越大起皱的程度越小，但压边力过大会导致凸缘材料被压死如以下三幅图所示图73330KN压边力图734KN压边力650图735KN压边力10通过以上的结果演示可以得到如下的结论压边力的大小对拉深成形时的起皱现象有显著的影响。压边力应该适当大小，但最小不能小于最小压边力，最大也应该有限制（保证凸缘材料顺利流入凹模型腔）739改变壳单元高斯积分点的数模对成形结果的影响表77最大应力与应变高斯积分点数目最大应力等效应变5236916（节点445）130691（445）7236916（节点445）130691（445）9236916（节点445）130691（445）11236916（节点445）130691（445）由上表可知，将拉深变形板料的高斯积分点从5依次增加到7、9、11计算的最大应力和等效应变没有变化。因此，对于拉深成形来讲，将高斯积分点设置成5个就足以达到足够的计算精度。7310回弹的探究表78回弹前后的尺寸对比尺寸零件底部直径零件中部直径零件口部直径540032041184341485411400418412578416707由上表格可以看出，零件的口部回弹量最大，而零件底部由于与凸模接触的时间较长，回弹量较少。74主要四部分的摩擦对拉深成形效果的影响各部分的摩擦条件对拉深成形的影响拉深成形过程中摩擦系数的初始设定值表79各部分摩擦值的设定板料与凹模圆角处的接触摩擦系数01板料与凸模圆角处的接触摩擦系数01板料与压边圈接触处的接触摩擦系数01板料与凹模上表面接触处的接触摩擦系数01板料与凹模内腔的接触摩擦系数0板料与凸模外圆柱面的接触摩擦系数0板料与凸模底面的接触摩擦系数0741改变板料与凹模圆角处的接触摩擦系数对拉深成形的影响前三种摩擦与材料的流动方向相反，阻碍材料向凹模型腔流动，增加了筒壁承受拉应力的负担，对拉深不利。凸模圆角处的材料受拉有向外流动的趋势，摩擦力指向内，这部分摩擦力对拉深有利，一方面可以减缓凸模圆角区材料的变薄，另一方面可以抵消一部分筒壁拉应力。下图为摩擦系数为0、01、025时应变曲线。由此曲线可知凹模圆角处的摩擦力越大应变值越大即对拉深越不利，故应该想办法减小凹模圆角处的摩擦。图736应变图下图为摩擦系数为0、01、025、035时应变曲线。由此曲线可看出当摩擦系数为035时节点1204的应变值远大于摩擦系数为0、01、025时的应变值。原因是凹模圆角处的摩擦力过大导致在凸缘区的材料向凹模型腔流动时被卡死，筒底部分的材料通过减薄而被单方面拉长，使网格过度变形。严重时会导致材料被拉裂。故应尽量减小凹模圆角处的摩擦力图737应变图图738应力图表710最大应力与等效应变摩擦系数最大屈服应力等效应变01835（502）703477E0301197221（32）711426E03025200661（100）732329E03035250（1337）305131由应力曲线可以看出当时间增长到一定程度后摩擦系数为035的曲线反而呈现下降的趋势。原因是筒底部分的材料过分变形导致应力应变的规律超出了胡克定律的范围呈现非线性的缘故。742改变板料与凸模圆角处的接触摩擦系数对拉深成形的影响图739应力图图740应变图表711最大应力与等效应变摩擦系数最大屈服应力等效应变0202749（42）737844E0301201255（451）730256E03025200743（572）723875E0303520024（110）698594E03由表中的数据和应力应变图可以看出，凸模圆角处的摩擦对拉深是有利的。743改变板料与压边圈接触处的接触摩擦系数对拉深成形的影响图741应变图图742应力图表712最大应力与等效应变摩擦系数最大屈服应力等效应变0204227（483）773388E0301202235（568）752325E03025200875（557）738509E03035200462（616）713206E03从表中数据和应力应变图中可以看出，压边圈与板料之间的摩擦系数越小同一节点的应力应变数值越小，及拉深性越好。744改变板料与凹模上表面接触处的接触摩擦系数对拉深成形的影响图743应变图图744应力图表713最大应力与等效应变摩擦系数最大屈服应力等效应变0201326（503）748709E0301197221（32）711426E03025201807（509）737633E03035199193（149）741354E03由表中数据和应力应变图看不出凹模面对板料的摩擦与拉深性的明确关系，但是凹模面对板料的摩擦对材料流入凹模型腔起阻碍的作用，很明显是不利的。八工艺参数对成形结果的影响规律改变各参数对拉深性的影响拉深性是指材料对拉深成形的适应能力，通常指不发生断裂的最大可能性。影响拉深性的主要因素通常有板料的相对厚度、凸模圆角半径、凹模0T/DPR圆角半径、凸凹模间隙等方面。通过对后处理结果的分析并结合所学的有关DR拉深的知识可以得出以下规律1板料的相对厚度越大，拉深性越好。因为相对厚度大时凸缘材料抗压失稳刚度大，不易失稳起皱，提高了拉深性能。2凸模圆角是拉深时断裂的位置所在，凸模圆角半径的值对拉深性有着PR显著的影响。当较小时，凸模圆角处材料变形大，板厚减薄严重，使危险断PR面强度降低。当较大时，凸模与板料的接触面积减小，接触压力增大，双向拉应力引起筒底板料在大范围内减薄，使拉深性降低。另外，当值较大时，PR凸模与凹模之间悬空部分材料增多，容易引起內皱。生产中，的值一般在下式给定的范