CN106156446B - 基于数值模拟的塑性成形过程能耗分析方法 - Google Patents

基于数值模拟的塑性成形过程能耗分析方法 Download PDF

Info

Publication number
CN106156446B
CN106156446B CN201610613885.1A CN201610613885A CN106156446B CN 106156446 B CN106156446 B CN 106156446B CN 201610613885 A CN201610613885 A CN 201610613885A CN 106156446 B CN106156446 B CN 106156446B
Authority
CN
China
Prior art keywords
energy consumption
forming
plastic forming
plastic
forming process
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610613885.1A
Other languages
English (en)
Other versions
CN106156446A (zh
Inventor
温彤
秦剑
朱滕薛
郑杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Original Assignee
Chongqing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University filed Critical Chongqing University
Priority to CN201610613885.1A priority Critical patent/CN106156446B/zh
Publication of CN106156446A publication Critical patent/CN106156446A/zh
Application granted granted Critical
Publication of CN106156446B publication Critical patent/CN106156446B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Forging (AREA)

Abstract

一种基于数值模拟的塑性成形过程能耗分析方法，包括以下步骤：(1)分析制件的成形工艺，按体积成形、板料成形以及坯料是否需要加热分类，确定成形道次；(2)在三维造型软件上构建制件坯料和模具的几何模型，对热成形件计算坯料体积；(3)将制件的材料力学性能参数、边界条件与成形工艺参数输入数值模拟软件，模拟各道次的成形过程；(4)根据模拟结果计算各道次的成形能耗，包括塑性变形能耗和加热能耗；(5)将所有道次的成形能耗累加，得到塑性成形过程的总能耗；(6)改变输入参数并重复步骤(4)‑(5)，得到输入参数调整后的总能耗。利用本发明可以快速、准确地计算成形过程的能耗，同时分析不同因素对能耗的影响。

Description

基于数值模拟的塑性成形过程能耗分析方法
技术领域
本发明属于金属材料的加工领域，尤其涉及一种基于数值模拟的塑性成形过程能耗分析方法。
背景技术
金属塑性成形是利用压力机等锻压加工设备或数控机床等机械加工设备，结合相应的模具、夹具等工装，将一定的外力施加在坯料上使其产生塑性变形，从而得到一定形状和尺寸并满足相应力学性能要求的制件的一种金属加工方法。根据坯料形式以及具体实现方式的不同，金属塑性成形包括体积成形(如锻造、挤压、拉拔等)和板料成形(如冲压、渐进成形以及胀形等)两大类型；按照成形时坯料是否预热，又分为不加热的冷成形(即室温下的成形)，和加热状态下的温成形(加热到室温与金属的再结晶温度之间)与热成形(加热到再结晶温度以上)等不同类型。金属塑性成形具有生产效率与材料利用率高、加工质量好、能大幅改善材料的组织性能等其他加工方法难以具备的综合优势，因而在制造领域的应用极其广泛。据不完全统计，汽车上多于60％零件，飞机上多于40％的零件要用到塑性成形的方法进行生产。
在全球能源短缺和环境污染不断加剧的情况下，既要保证经济的快速发展，又要满足节能减排的可持续发展战略要求，这已经是制造业的普遍共识。作为一种关系民生的基础工业技术，金属塑性成形也存在着能耗相对较高的缺点，要实现“低碳化”，塑性成形行业必须大大降低生产过程的能量消耗。迄今为止，对金属塑性成形的节能减排已有大量研究，从改进工艺、节约原料，采用净近成形的冷、温精锻技术，到加热与成形设备的节能改造，研发与推广新型低排放的节能设备，以及合理配置生产线辅助设施、从而发挥其最高效能等等，都是塑性成形从业者以及相关的成形设备制造商研究的课题。
对塑性成形产品而言，生产过程的能量消耗主要包括塑性成形阶段的能耗以及后续处理(如热处理、表面清理、切削加工等)的能耗。其中，塑性成形阶段的能耗是制件整个生产过程能耗的重要部分，它通常包括两大部分：(1)使坯料产生塑性变形的输入外力所做的机械功，以及(2)热成形时对坯料进行加热的能耗。由于塑性成形的具体工艺方法繁多，同一制件可以选择不同的成形方法、或者采用不同的工序组合进行生产，目前尚无能够对各种工艺方法的能耗进行快速准确分析的系统手段。迄今为止，虽然已经提出了许多针对各类塑性成形过程的变形力计算公式，但主要用于选择成形设备的规格(吨位)。与此同时，数值模拟技术在塑性成形过程的应用已经十分成熟，但目前主要是用于分析材料成形过程的应力、应变等力学参量以及流动变形的状态，很少涉及成形过程能耗的计算与分析。
综上所述，为了能够计算出各类零件在具体塑性成形过程中的能耗，找出影响能耗的主要因素，从而为实现节能减排提供必要的参考，有必要建立一套可行、可信的塑性成形过程的能耗分析方法。
发明内容
本发明提供了一种基于数值模拟的塑性成形过程能耗分析方法，可以针对各种金属塑性成形制件的具体工艺过程，快速而准确地计算出成形阶段的总能量消耗。
为了实现上述目的，本发明包括以下步骤：
(1)分析待加工制件的塑性成形工艺，按照体积成形、板料成形，以及坯料是否需要加热进行分类，确定成形道次；
(2)在三维造型软件上构建制件坯料和塑性成形模具的几何模型，计算需要加热成形的制件的坯料体积V；
(3)将制件的材料力学性能参数、边界条件与成形工艺参数输入数值模拟软件，对各道次的塑性成形工序进行模拟；
(4)根据数值模拟结果，计算各道次塑性成形的能耗En，所述En包括塑性变形能耗An和热成形的加热能耗Qn两部分。
第n道次塑性成形的塑性变形能耗An采用如下公式计算：
An＝S*K1
其中，S为第n道次塑性成形的材料塑性变形能，对于普通锻造和冲压成形，S等于数值模拟得到的载荷-行程曲线与横坐标所包围的面积；对于板料的渐进成形，S由数值模拟软件输出的外力所做的功确定；K1为考虑塑性成形过程能量损失的修正系数。
对于不需要加热的冷成形，加热能耗Qn为零；对于加热条件下的温成形和热成形，第n道次塑性成形的加热能耗Qn为将坯料从室温加热到相应成形温度的能耗，采用如下公式计算：
Qn＝c*V*ρ*(T1-T0)*K2
其中，c为材料的比热容，ρ为材料密度，T1为坯料加热的最高温度，T0为室温，K2为考虑加热过程损耗的修正系数。
第n道次塑性成形的能耗En，等于该道次塑性成形的塑性变形能耗An与加热能耗Qn之和：
En＝An+Qn
(5)将所有道次的塑性成形能耗En进行累加，得到制件整个塑性成形过程的总能耗E：
E＝ΣEn
(6)改变步骤(3)的输入参数，重复步骤(4)-(5)，得到调整输入参数后的总能耗E，比较不同输入参数下的总能耗E，找出对总能耗E影响较大的因素。
本发明具有以下有益效果：利用各类塑性成形制件成形过程的数值模拟结果，可以快速、准确地计算出整个成形过程的能量消耗，同时分析不同因素对能耗的影响程度。
附图说明
下面结合附图并分别以一个铝合金锻件的热锻成形和一个锥台壳体零件的渐进成形为例，对本发明进一步说明。
图1为本发明基于数值模拟的塑性成形过程能耗分析方法的流程图；
图2为热锻成形铝合金锻件的坯料示意图；
图3为利用DEFORM软件模拟得到的铝合金锻件热塑性成形过程的“载荷－行程”曲线；
图4为渐进成形的锥台壳体零件的示意图；
图5为锥台壳体渐进成形采用的等高线加工轨迹；
图6为锥台壳体采用等高线加工轨迹进行渐进成形时，利用Abaqus软件模拟得到的工具头在x、y、z三个方向的“作用力－时间”曲线；
图7为锥台壳体渐进成形采用的螺旋线加工轨迹；
图8为锥台壳体采用螺旋线加工轨迹进行渐进成形时，利用Abaqus软件模拟得到的工具头在x、y、z三个方向的“作用力－时间”曲线；
图9为利用Abaqus软件模拟得到的锥台壳体渐进成形过程的“能耗－时间”变化曲线，分别采用等高线和螺旋线两种加工轨迹。
具体实施方式
1、铝合金锻件的热锻成形
该铝合金锻件的材料为Al 7050，采用热锻成形。锻件坯料在水平面的外形尺寸(长*宽)为1200mm*300mm，属于特大型锻件。
(1)分析零件的塑性成形工艺，该锻件采用热模锻压力机上的模锻成形，属于体积成形范围且需要加热成形，成形道次为1次。
(2)在三维造型软件上构建锻件坯料和成形模具的几何模型，坯料体积V＝3.64X107(mm3)；
(3)将材料力学性能参数、边界条件以及成形工艺参数输入体积塑性成形的数值模拟软件DEFORM，对铝合金锻件的热锻成形过程进行模拟；
(4)计算铝合金锻件热锻成形过程的塑性变形能耗An和加热能耗Qn
该铝合金锻件热锻成形道次为1，该道次塑性变形能耗A1利用DEFORM数值模拟得到的“载荷－行程”曲线计算。具体做法是，将“载荷－行程”曲线导入Origin软件中，通过菜单栏中Analysis–mathematics-integrate求积分功能，得到“载荷－行程”曲线与横坐标所包围的面积，即为该道次塑性成形的材料塑性变形能S。对于该铝合金锻件，Origin计算出area＝8490171369.731，即S＝8490(KJ)。
该铝合金锻件的锻造采用热模锻压力机，根据现有资料，该类设备的能源利用率为20％左右，因此可确定塑性成形过程能量损耗的修正系数K2＝5。于是，该铝合金锻件热锻成形的塑性变形能耗A1为
A1＝S*K1＝8490*5＝42450(KJ)
7050铝合金的密度为2.81g/cm3，比热容c为860J/(kg*K)，始锻温度T1＝400℃，室温T0＝20℃。采用天然气炉加热，根据资料可确定考虑加热损耗的修正系数K2＝10。于是加热能耗Q1为：
Q1＝c*V*ρ*(T1-T0)*K2＝334264(KJ)
该铝合金锻件热锻成形的总能耗E，即为塑性变形能耗A1与加热能耗Q1之和：
E＝A+Q＝A1+Q1＝376714(KJ)
2、锥台壳体的渐进成形
该锥台壳体零件的小端直径80mm，大端直径160mm，高度40mm，壁厚1mm，材料为Al6061。采用不同于传统冲压的渐进成形方式制造，该方法基于“分层制造”的思想，利用一个简单的柱状工具头在数控程序的控制下对板坯料进行逐渐的成形，直到完成整个制件的加工。工具头可以采用不同的加工路径轨迹，本实施例中，工具头分别采用等高线和螺旋线两种轨迹，进给量△Z＝2mm。
具体能耗分析步骤如下：
(1)分析零件的成形工艺，该零件为钣金件且采用渐进成形，坯料不需要加热，成形道次为1次；
(2)在三维造型软件上构建锥台壳体零件的坯料和成形模具的几何模型；
(3)将成形工艺参数、材料力学性能参数以及边界条件输入数值模拟软件Abaqus，对渐进成形过程进行模拟，工具头首先采用等高线轨迹；
(4)根据Abaqus对渐进成形过程的模拟结果，利用软件的History Output模块输出整个成形过程中外力所做的功，于是得到采用等高线轨迹时消耗的功S＝3.553KJ。考虑塑性成形过程能量损耗的修正系数K1为5，则整个渐进成形过程塑性变形的能耗A1为
A1＝S*K1＝17.765(KJ)
(5)将工具头的加工轨迹改变为螺旋线形式，其余输入参数不变，重复上述数值模拟步骤，得到采用螺旋线轨迹时所消耗的功为S＝3.566KJ，于是采用螺旋线轨迹时
A1＝S*K1＝17.83(KJ)
可见工具头加工轨迹对该锥台壳体零件渐进成形过程的能耗影响较小。

Claims (4)

1.一种基于数值模拟的塑性成形过程能耗分析方法，其特征在于，包括以下步骤：
(1)分析待加工制件的塑性成形工艺，按照体积成形、板料成形，以及坯料是否需要加热进行分类，确定成形道次；
(2)在三维造型软件上构建制件坯料和塑性成形模具的几何模型，计算需要加热成形的制件的坯料体积V；
(3)将制件的材料力学性能参数、边界条件与成形工艺参数输入数值模拟软件，对各道次的塑性成形工序进行模拟；
(4)根据数值模拟结果，计算各道次塑性成形的能耗En；
(5)将所有道次的塑性成形能耗En进行累加，得到制件整个塑性成形过程的总能耗E：
E＝ΣEn
(6)改变步骤(3)的输入参数，重复步骤(4)-(5)，得到调整输入参数后的总能耗E，比较不同输入参数下的总能耗E，找出对总能耗E影响较大的因素。
2.如权利要求1所述的基于数值模拟的塑性成形过程能耗分析方法，其特征在于，第n道次塑性成形的能耗En，等于该道次塑性成形的塑性变形能耗An与加热能耗Qn之和：
En＝An+Qn。
3.如权利要求1所述的基于数值模拟的塑性成形过程能耗分析方法，其特征在于，第n道次塑性成形的塑性变形能耗An采用如下公式计算：
An＝S*K1
其中，S为第n道次塑性成形的材料塑性变形能，对于普通锻造和冲压成形，S等于数值模拟得到的载荷-行程曲线与横坐标所包围的面积；对于板料的渐进成形，S由数值模拟软件输出的外力所做的功确定；K1为考虑塑性成形过程能量损失的修正系数。
4.如权利要求1所述的基于数值模拟的塑性成形过程能耗分析方法，其特征在于，对于不需要加热的冷成形，加热能耗Qn为零；对于加热条件下的温成形和热成形，第n道次塑性成形的加热能耗Qn为将坯料从室温加热到相应成形温度的能耗，采用如下公式计算：
Qn＝c*V*ρ*(T1-T0)*K2
其中，c为材料的比热容，ρ为材料密度，T1为坯料加热的最高温度，T0为室温，K2为考虑加热过程损耗的修正系数。
CN201610613885.1A 2016-07-28 2016-07-28 基于数值模拟的塑性成形过程能耗分析方法 Expired - Fee Related CN106156446B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610613885.1A CN106156446B (zh) 2016-07-28 2016-07-28 基于数值模拟的塑性成形过程能耗分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610613885.1A CN106156446B (zh) 2016-07-28 2016-07-28 基于数值模拟的塑性成形过程能耗分析方法

Publications (2)

Publication Number Publication Date
CN106156446A CN106156446A (zh) 2016-11-23
CN106156446B true CN106156446B (zh) 2019-04-12

Family

ID=57327965

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610613885.1A Expired - Fee Related CN106156446B (zh) 2016-07-28 2016-07-28 基于数值模拟的塑性成形过程能耗分析方法

Country Status (1)

Country Link
CN (1) CN106156446B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108491588B (zh) * 2018-03-02 2021-10-22 武汉理工大学 推力滚子轴承套圈冷空间包络成形方法
CN108436474B (zh) * 2018-04-26 2020-02-07 重庆大学 一种结合钻孔与塑性变形的弯曲孔加工方法
CN111595697B (zh) * 2020-04-10 2023-08-25 唐山钢铁集团有限责任公司 一种用于评判冲压时材料抗剪切成形能力的方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100585603C (zh) * 2007-10-18 2010-01-27 中国兵器工业第五九研究所 一种塑性成形数值模拟方法
CN103020460B (zh) * 2012-12-20 2015-08-05 重庆大学 一种降低钢筋连接能耗的方法
CN104077439B (zh) * 2014-06-19 2017-04-26 西北工业大学 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法

Also Published As

Publication number Publication date
CN106156446A (zh) 2016-11-23

Similar Documents

Publication Publication Date Title
Lu et al. Model predictive control of incremental sheet forming for geometric accuracy improvement
CN104573237B (zh) 一种基于摩擦磨损cae分析的模具优化方法
CN103381441B (zh) 一种薄壁钛合金封头的热冲压成型加工方法
CN104077439B (zh) 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法
CN106156446B (zh) 基于数值模拟的塑性成形过程能耗分析方法
CN104200037B (zh) 一种小双曲度钣金件的成形模具设计方法
CN104765912A (zh) 一种铝板冲压工艺的稳健性优化方法
CN110750931A (zh) 一种高效的型材挤压模具设计方法
CN109127945A (zh) 轻量化车身覆盖件冲压成形精度的调控方法
CN105528484A (zh) 工件成型力的cae仿真分析方法以及成型装置优化方法
Sun et al. Numerical simulation of extrusion process and die structure optimization for a complex magnesium doorframe
CN104636565B (zh) 基于可加工性分析的镁合金一次模锻成形工艺优化方法
CN106503389A (zh) 一种快速确定热冲压生产线工艺参数的数值方法
Ji et al. Optimization of the extrusion die and microstructure analysis for a hollow aluminum alloy profile
Al-Ghamdi et al. Forming forces in incremental forming of a geometry with corner feature: investigation into the effect of forming parameters using response surface approach
CN103341943A (zh) 塑料板料高压热风渐进成形装置及其成形方法
Zhou et al. The multi-objective optimization design of a new closed extrusion forging technology for a steering knuckle with long rod and fork
Li et al. Accurate trimming line optimization of multi-station progressive die for complex automotive structural parts
CN109570310A (zh) 一种新型钛合金曲面成型方法
CN203371715U (zh) 塑料板料高压热风渐进成形装置
CN111889601B (zh) 增大难变形合金模锻件表层变形的波浪形坯料设计方法
CN116653347A (zh) 一种基于五次多项式优化的伺服冲压工艺曲线控制方法
Kong et al. Numerical and experimental investigation of preform design in non-axisymmetric warm forming
CN207508182U (zh) 一种环套类零件批量锻造冲孔模具
CN108804816A (zh) 一种轨道车辆用抗侧滚扭杆轴镦粗工艺仿真分析方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190412

Termination date: 20200728