搜索
基于KMAS的国际板成形Benchmark试验研究*闫康康，郭威，申国哲，胡平1吉林大学车身与模具工程研究所，汽车动态模拟国家重点实验室，长春，130025摘要：通过运用自主开发的商品化冲压成形CAE软件KMAS（King-meshAnalysisSystem）对国际板成形会议（NUMISHEET2002）上提出的若干Benchmark试验数据进行计算机仿真，把模拟结果与实验值及其他软件模拟结果进行对比。对在不同冲压工艺下，KMAS对板材的各向异性凸缘制耳、起皱、接触状况、回弹的模拟精度做出评估，验证KMAS的动力半显式算法和回弹显式卸载算法的正确性。关键词:板材成形数值模拟；KMAS软件；NUMISHEET2002；Bechmark试验中图分类号：TG385ResearchaboutInternationalBenchmarkTestofSheetMetalFormingBasedonCAESoftwareKMASYanKangkang1,GuoWei1,ShenGuozhe2,HuPing2(1.DepartmentofMaterialScienceandEngineeringofJilinUniversity,Changchun,130025;2.InstituteofAuto-bodyandDieEngineeringofJilinUnversity,Changchun,130025)Abstract:TheBenchmarktestingproposedfromNUMISHEET2002aretakenbyusingindependentlydevelopedcommercialsoftwareKMAS（King-meshAnalysisSystem）forsheetmetalforming,whichisacommercialCAEsoftwareathome.ComparedexperimentalresultswiththeresultsobtainedbyotherCAEsoftware(suchasPAM-STAMP,DYNAFORMetc.),thesimulatedresultscomingfromKMASareshowntobewithveryhighsimulationaccuracy,especiallyforthesimulatedresultsinanisotropicofmaterial,wrinkling,contacttreatmentandspring-backanalysis,thevalidityofpresentdynamicsemi-explicitalgorithmwithdynamicloadingandexplicitunloadingschemeforspring-backisprovedout.Keywords:Sheetmetalformingnumericalsimulation;KMASsoftware;NUMISHEET2002；Benchmarktesting;1.引言在汽车生产过程中，大多数结构件和覆盖件都是由板材冲压制造，而板材的成形性能和模具设计则是实现冲压工艺的关键，因而对板材的冲压成形分析（CAE）就显得至关重要。目前板材冲压成形CAE软件都以弹塑性大变形有限元技术为基础，并且已在国内外车身与模具设计中广泛应用[1]。*国家杰出青年科学基金项目（10125208）和国家自然科学基金重点项目（19832020）的联合资助1*通讯联系人：胡平，吉林大学教授，男，pinghu@jlu.edu.cn,联系电话：0431－5687162KMAS(King-meshAnalysisSystem)是由吉林大学自主开发的商业化的三重非线性有限元专用分析软件，其CAE前处理数据准备更加简捷（与PAM-STAMP和DYNAFORM比较），回弹模拟的准确性在多数情况下高于国外同类CAE软件。本文所研究的Benchmark是2002年在韩国举行的板材成形数值模拟会议(NUMISHEET)提出的标准试题,用于对全世界板成形CAE软件进行测试评估。NUMISHEET是国际上最著名的数值模拟会议之一,其Benchmark对于板材成形CAE软件的测试具有权威性[2],其测试结果是值得信赖和借鉴的。2.KMAS的基本理论KMAS由网格生成器(Mesher)，前处理器(Top)，求解器(Solver)，后处理器(View)四部分组成。可对一般（一次）成形、多阶段成形、回弹变形过程进行仿真[1]，成形方式包括工业常用的单动与双动压机两种。2.1.KMAS算法模型特点1.膜单元组弯曲理论，建立了金属薄板冲压过程的空间拟弯曲膜单元模型[3]，有效地解决了“沙漏控制”问题,实例分析结果表明能够有效地模拟弯曲和起皱问题。2.采用Mindlin板壳单元理论，四结点等参单元模型，共旋坐标理论，建立了可考虑横向剪切应变的适合于金属板材冲压成形数值模拟的BT壳有限元数学模型。3.采用具有微观物理基础的B-L各向异性屈服准则,可以有效模拟出塑性诱导各向异性的剪切“制耳”。4.在动力显式中心差分算法的基础上，采用牛顿迭代法，提出了基于速度迭代的弹塑性大变形动力半显式算法[4]，给出了迭代格式和收敛标准。采用该算法能明显提高加载过程的应力场精度，从而获得了比单纯用动力显式模拟加载过程更好的回弹结果。5.针对回弹问题，KMAS提出了一种新的求解方法，该方法采用有模具卸载方式，将动力显式中心差分方法与动力松弛方法结合起来。该方法弥补了单纯用动力松弛法计算回弹时间长、不易收敛的问题，以及采用隐式方法计算时的内存消耗大的问题。6.为提高凹模圆角处的模拟精度,KMAS引入了有限元自适应网格重划加密技术(AdaptiveRemeshing),使得在圆角处网格更细密,更加贴近实验效果。3.BenchmarkA:DeepDrawingofCylindricalCup试验研究3.1试验目的(1)研究板材在高压边力下变形的各向异性行为。(2)研究板材在低压边力下变形时的起皱现象。3.2.实验所用材料参数(右表)3.3相关实验参数铝合金板（Al6111-T4）的实物冲压实验由日本神户钢铁公司(KobeSteel)的HaruyukiKonishi(AE-01)和美国铝业公司(ALCOA)的RobertEDick共同完成;低碳钢板的实物冲压实验由英国-荷兰钢铁厂研发中心的(CorusGroupResearch,D&T)的E.H.Atzema（AE-07）完成，并提供实验参数及结果数据。其中：实验设备为双动液压机（如图1），所有工具(凹模,凸模,压边圈)由硬质工具钢(hardenedtoolsteel)-SKD11制造,工具表面热处理,工具表面粗糙度小于5微米。铝合金冲压坯料用激光切割为直径为90.0mm的圆形；低碳钢板冲压坯料用激光切割为直径为105.0mm的圆形。实验时液压机冲头速度：30mm/sec，冲头行程：40mm，润滑剂采用拉延油(FD-1500)。压边力铝材分别作50KN和10KN两组实验，钢板取70KN和10KN两组实验。143面内各向异性系数材料名称及型号板厚(mm)R0R45R90Al6111-T41.00.8940.6110.66STEELDDQ1.02.161.6112.665屈服应力(MPa)弹性模量(MPa)硬化指数强化系数(MPa)泊松比181.24705000.228539.290.34161.952213700.269547.760.3图1.:实验液压机及模具fig1:MachineandToolingSpecifications3.4.KMAS数值模拟参数KMAS数值模拟参数基本上和实验参数相同，考虑到液压机的保压性取冲头速度0.1m/s，压边速度0.5m/s。因为润滑条件良好，取摩擦系数为：0.0096（铝）和0.0426(钢)。模拟建立几何模型如图2所示。3.5.各测试点表征方案(1)冲头最终行程（40mm）时板材从中心到边沿的壁厚分布作来研究材料在拉深后的减薄程度，沿与轧制方向成0、45、90度方向截面上分别测量厚度作三次实验。（2）冲头力（kN）与冲头行程（mm）组合曲线来研究冲压过程中压力的变化。（3）板材变形后外轮廓（制耳）几何形状，用制件中心到边沿的距离沿圆周角度展开来表示制件外轮廓的展开图。角度展开起始位置如图3所示。（4）法兰面上的起皱波纹条数来研究KMAS对于起皱的预示情况。3.6.实验结果与模拟结果比较模拟结果：KMAS对于板料采用B-T薄壳[5]单元模型划分有限元网格，其中铝合金板共4104个单元，低碳钢板共10297个单元，计算采用Pentium4—2.4GHz微机，计算时间共16分38秒，消耗内存3.56MB。其中对比图中AE-01…等表示实验值，KMAS、Dynaform等表示所用软件的模拟值。3.6.1.不同材料在不同压边力下成形制件在沿不同方向厚度分布对比(表3-1---3-3)表3-1：沿轧制方向制件截面的厚度分布table3-1:Thicknessdistributioncutat0degreefrom(RD)图2.工具和板料几何尺寸(单位:毫米)Fig2.ToolsandblankgeometrydimensionR1=50.0,R2=51.25,R3=9.5,R4=7.0表3-1--3-3中横轴表示制件中心到边沿的距离，纵轴表示制件与板料轧制方向成一定角度的截面厚度分布。其中上图表示压边力（BHF）=50KN的铝合金板试验，下图表示压边力（BHF）=10KN的低碳钢板试验。144145表3-2:与轧制方向成45度截面的厚度分布table3-2:Thicknessdistributioncutat45degreefromRD图3.制件外轮廓形状试验角度取向方式Fig3.Coordinateformeasuringangle表3-3:与轧制方向成90度截面厚度分布table3-3:Thicknessdistributioncutat90degreefromRD3.6.2不同压边力下成形时冲头力与行程曲线对比表3-4:冲头力与冲头行程曲线table3-4:Punchforce(kN)vs.punchdisplacement(mm)表中横轴表示凸模行程，纵轴表示冲头压力。上下图分别表示压边力为50KN和10KN的铝合金板试验。3.6.3.高压边力下材料拉延时表现出的各向异性对到边缘距离（即凸缘制耳）沿圆称，起皱数table4-6:NumberofWrinkingwavesattheflange.7.试验结果小结：材料表现出略微的的各向异性，低压边力对于感。4.试比示意图(表3-5)表3-5:铝合金制件中心周方向展开分布图table3-5:Outerprofile(R0)asafunctionofθ.其中角度θ取值方法见图4-3所示，因为制件中心对KMAS模拟数据只测量了半个法兰面外缘形状[6]。3.6.4.低压边力下材料拉延时的起皱情况对比下图表示了铝合金板在压边力为10KN时，压料面上量的实验结果及各种软件（纵轴）模拟结果对比。表4-6压料面上起皱数量3法兰面起皱的影响相当敏BenchmarkB:UnconstrainedCylindricalBending验研究4.1.测试目的:(1)研究板材弯曲时的回弹现象。(2)研究不同行程下冲头与板材的复杂接触状态。4.2.实验所用材料参数面内各向异性系数材料名称及板厚型号(mm)R0R45R90A0.8966l6111-T41.040.6110.STEELSHS1.01.522.0962.357屈服应力(MPa)弹性模量(强化系数MPa)硬化指数(MPa)泊松比181.24705000.228539.290.34219.812221700.253631.1380.34.3.验参的实压验由韩国大学niversity)的Kim和HyungJong共同完成,即图中所示BE-01，并提供实验参数及结果数据。具表