CN113111540B - 面向铝压铸的熔化、配送和保温参数集成的节能优化方法 - Google Patents
面向铝压铸的熔化、配送和保温参数集成的节能优化方法 Download PDFInfo
- Publication number
- CN113111540B CN113111540B CN202110491310.8A CN202110491310A CN113111540B CN 113111540 B CN113111540 B CN 113111540B CN 202110491310 A CN202110491310 A CN 202110491310A CN 113111540 B CN113111540 B CN 113111540B
- Authority
- CN
- China
- Prior art keywords
- melting
- heat preservation
- furnace
- melting furnace
- energy
- 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.)
- Active
Links
- 238000002844 melting Methods 0.000 title claims abstract description 189
- 230000008018 melting Effects 0.000 title claims abstract description 189
- 238000004321 preservation Methods 0.000 title claims abstract description 129
- 238000004512 die casting Methods 0.000 title claims abstract description 54
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 title claims abstract description 51
- 229910052782 aluminium Inorganic materials 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- 238000009826 distribution Methods 0.000 title claims abstract description 15
- 230000010354 integration Effects 0.000 title claims abstract description 9
- 238000004519 manufacturing process Methods 0.000 claims abstract description 33
- 238000005265 energy consumption Methods 0.000 claims abstract description 27
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 5
- 230000002068 genetic effect Effects 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims description 41
- 239000007788 liquid Substances 0.000 claims description 22
- 230000014509 gene expression Effects 0.000 claims description 9
- 238000009413 insulation Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000003647 oxidation Effects 0.000 claims description 4
- 238000007254 oxidation reaction Methods 0.000 claims description 4
- 229910045601 alloy Inorganic materials 0.000 claims description 3
- 239000000956 alloy Substances 0.000 claims description 3
- 239000000725 suspension Substances 0.000 claims description 3
- AZDRQVAHHNSJOQ-UHFFFAOYSA-N alumane Chemical compound [AlH3] AZDRQVAHHNSJOQ-UHFFFAOYSA-N 0.000 claims description 2
- 239000011810 insulating material Substances 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 14
- 238000003723 Smelting Methods 0.000 abstract description 8
- 230000003068 static effect Effects 0.000 abstract description 2
- 238000005266 casting Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 238000004134 energy conservation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000007872 degassing Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012768 molten material Substances 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22D—CASTING OF METALS; CASTING OF OTHER SUBSTANCES BY THE SAME PROCESSES OR DEVICES
- B22D17/00—Pressure die casting or injection die casting, i.e. casting in which the metal is forced into a mould under high pressure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- 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/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Mechanical Engineering (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Manufacturing & Machinery (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Vertical, Hearth, Or Arc Furnaces (AREA)
- Manufacture And Refinement Of Metals (AREA)
- Furnace Details (AREA)
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
本发明公开了一种面向铝压铸的熔化、配送和保温参数集成的节能优化方法,包括:(1)通过现场或者试验收集优化模型需要的各种静态参数和系数;(2)建立基于仿真的成本最优化模型,该模型以熔炉能耗、熔炉烧损、运输能耗、运输人力成本、保温炉能耗、生产暂停惩罚综合成本最小为目标;(3)用遗传算法搜寻近似最优的参数组合,对于一个给定生产计划,将计划参数输入到模型中,可快速地给出成本近似最优的参数组合。本发明方法能对一个给定的压铸生产计划选出成本最低的铝压铸熔化、配送和保温的参数组合,能有效降低制造成本。
Description
技术领域
本发明属于铝压铸过程节能优化技术领域,具体涉及一种面向铝压铸的熔化、配送和保温参数集成的节能优化方法。
背景技术
工业是碳排放的主体,是我国实现“2030年碳达峰”和“2060年碳中和”的重要领域,是当前实施节能减排的首要对象。作为一种汽车轻量化的使能技术,铝压铸过程具有高耗能高排放的特点,亟需开展节能减排工作。压铸过程中,熔化和保温过程耗能大(消耗了约50%的总能耗)、能效低,存在较大节能空间,是压铸过程节能减排的首要对象,值得关注。
如图1所示,一个典型的压铸车间配有n(n≥1)台熔化炉Fm={Fm1,Fm2,…,Fmn}以及k(k≥1)台保温炉Fh={Fh1,Fh2,…,Fhk}和压铸机,多台除气除渣设备及运输包,能生产多种铝合金压铸产品。在制定熔化、配送和保温的工艺参数时,每台压铸机需要的合金是由那台熔炉Fmn提供,每台压铸机的浇注质量Mpk、浇注的周期Tpk、每种产品的计划生产数量Qpk是已知的。每台熔炉的熔化速度Vmn、熔炉保温区的上限保温质量Mmupn、下限保温质量Mmlown、各保温炉从熔化炉每次取料的时刻与质量保温炉上限的保温高度Mhupk以及保温炉下限的保温高度Mhlowk是需要工艺人员制定的。熔化炉之间的关系是并列和互相不影响的,可以把每台熔炉Fmn和它供应的一个或者多个保温炉组成的集合分为一个整体,分别独立设置每个整体的运作参数。
为了防范铝液不能按时供给,当前工艺人员一般采用额定的熔化速度(最大熔化速度)、最高的熔炉上下限保温高度以及较大的保温上下限高度,且针对不同的生产计划,不会调整这些参数,保温炉加料的时间间隔也不精确,有时供料不及时或者保温炉保温量长期处于高位。当前采用的运作参数,一般能使得某些指标较好,如较高的熔化效率、较低熔化烧损、较少运输次数(使得人力成本较低)和较大的运输间隔,但是常常会导致非常高的保温能耗,造成了一定的浪费,对这些参数进行集成优选能使得整体的能耗成本有效降低,但是目前缺乏相应的方法。
针对上述问题,现有的技术比较缺乏,仅有一项相关的专利技术(公布号为CN111325408A),该专利针对的是分布式熔化模式下的铝液熔化和保温(一台熔化和保温集成炉对应一个压铸单元,没有铝液运输环节),对集中熔化模式(一台熔化炉,对应多台保温炉和多个压铸单元,有铝液配送环节)并不适用,因此需要研究新的技术来解决上述问题。
发明内容
鉴于上述,本发明提供了一种面向铝压铸的熔化、配送和保温参数集成的节能优化方法,能够针对不同的铝压铸生产计划及信息,给出最节能的熔化炉、运输和保温炉的参数方案,能有效的降低铝压铸过程的铝熔化、铝液配送和保温的成本。
一种面向铝压铸的熔化、配送和保温参数集成的节能优化方法,其针对车间内以某一台熔化炉为主的生产整体,该整体内还包括N台由熔化炉提供合金的压铸机及其对应的N台保温炉,N为大于1的自然数;
首先,获取已知的相关生产信息作为输入,包括该整体中压铸机每次的浇注质量、压铸机每次浇注的间隔时间、压铸机所需生产的产品数量、熔化炉在初始时刻的保温量、初始时刻熔化炉熔化烧嘴的状态(0代表待机,1代表正在熔化)、保温炉在初始时刻的保温量、初始时保温炉是否已取料(0代表未取料,1代表已取料)及取料时刻和取料量、初始时压铸机离下次加料的时间;
然后,建立以整体的总成本最小为目标的优化模型如下:
min{Co=Cm+∑Ctk+∑Chk+∑Cpk}
其中:Co为该整体的总成本,Cm为该整体中熔化炉的使用成本,Ctk为该整体中第k台保温炉的铝液运输成本,Chk为第k台保温炉的使用成本,Cpk为第k台保温炉因铝液缺少导致的暂停成本,k为自然数且1≤k≤N;
最后,对上述优化模型进行求解,得到该整体中熔化炉的熔化速度、熔化炉保温区的上限保温质量和下限保温质量、保温炉从熔化炉每次取料的时刻和质量、保温炉的上限保温质量和下限保温质量。
进一步地,所述优化模型的约束条件如下:
Vm∈[Vmmin,Vmmax]
Mmmin≤Mmlow<Mmup≤Mmmax
Mhmink≤Mhlowk<Mhupk≤Mhmaxk
其中:Vm为熔化炉的熔化速度,Vmmin和Vmmax分别为熔化炉的最小熔化速度和最大熔化速度,Mmup和Mmlow分别为熔化炉保温区的上限保温质量和下限保温质量,Mmmin和Mmmax分别为熔化炉保温区的最小保温质量和最大保温质量,Mhupk和Mhlowk分别为第k台保温炉的上限保温质量和下限保温质量,Mhmink和Mhmaxk分别为第k台保温炉的最小保温质量和最大保温质量。
进一步地,采用遗传算法对优化模型进行求解,求解过程中:先求得一组整体最小成本的优化参数,在该优化参数基础上运行使该整体中一台压铸机生产任务完成,然后以上述结束的状态为初始状态,求解另外一组整体最小成本的优化参数,在该优化参数基础上运行直到另一台压铸机生产任务完成,重复上述操作,直至所有压铸机生产任务完成。
进一步地,所述熔化炉使用成本Cm的表达式如下:
Cm=Em*pem+Mo*pa
Em=Emm+Emh+Ems+Emc
Emm=Pmm*Tmm
Ems=Pms*Tms
Emc=ΔEmc*Nmc
Pmh=cmh1*Mmh+cmh2
其中:Em为熔化炉的熔化能源消耗量,pem为熔化炉的熔化用能源单价,Mo为熔化炉的材料氧化烧损量,pa为熔化炉的熔化材料单价,Emm为熔化炉的熔化能耗,Emh为熔化炉保温区的铝液保温能耗,Ems为熔化炉的熔化待机能耗,Emc为熔化炉的熔化状态切换能耗,Pmm为熔化炉的熔化功率,Tmm为熔化炉的熔化时间,Pmh为熔化炉保温区的保温功率,Tmh为熔化炉保温区的保温时间,Pms为熔化炉的熔化待机功率,Tms为熔化炉的熔化待机时间,ΔEmc为熔化炉一次熔化状态切换的能量损失,Nmc为熔化炉的熔化状态切换次数,Vm为熔化炉的熔化速度,cmm1、cmm2、cmm3、cmh1、cmh2、cmo1、cmo2均为给定的系数,Mmh为熔化炉保温区的溶液质量。
进一步地,所述第k台保温炉铝液运输成本Ctk的表达式如下:
Ctk=Ntk*pet+Ntk*pht
Htk=Ntk*ct
其中:Etk为第k台保温炉的铝液运输能源消耗量,pet为铝液运输用能源的单价,Htk为第k台保温炉的运输人力使用量,pht为单位人力成本,Ntk为第k台保温炉的铝液运输次数,为第k台保温炉一次铝液运输消耗的能量,ct为人力使用量的折算系数。
进一步地,所述第k台保温炉使用成本Chk的表达式如下:
Chk=Ehk*peh
Phk=ch1*Mhk+ch2
其中:Ehk为第k台保温炉的能源消耗量,peh为保温用能源单价,Phk为第k台保温炉的保温功率,Thk为第k台保温炉的保温时间,ch1和ch2均为给定的系数,Mhk为第k台保温炉的保温铝液质量。
进一步地,所述第k台保温炉因铝液缺少导致的暂停成本Cpk表达式如下:
Cpk=Tmspk*pmsp
其中:Tmspk为第k台保温炉的暂停时长,pmsp为暂停的处罚系数。
基于上述技术方案,本发明具有以下有益技术效果:
1.本发明面向铝压铸的熔化、配送和保温参数集成的节能优化方法,能够针对不同的生产计划,快速的选出成本最低的工艺参数组合,能有效的降低压铸生产成本。
2.本发明可以植入到企业的制造执行系统,能够方便快递的进行熔炉、配送和保温炉参数优化,能有效的降低工艺人员的工作量。
附图说明
图1为压铸车间内铝熔化、铝液运输及保温的设备流水线示意图。
图2为本发明参数寻优流程示意图。
图3(a)为熔化炉的保温质量变化趋势示意图。
图3(b)为保温炉的保温质量变化趋势示意图。
图4为本发明仿真流程示意图。
图5为本发明具体实施案例的流程示意图。
图6为本发明具体实施案例中一组近似最优解的铝液变化示意图。
具体实施方式
为了更为具体地描述本发明,下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
本发明面向铝压铸的熔化、配送和保温参数集成的节能优化方法,包括如下步骤:
(1)以一个车间内某一台熔化炉(熔化炉编号为n)和多台保温炉(保温炉编号为k)组成的整体中每台保温炉给对应压铸机的每次浇注质量Mpk、每次浇注的间隔时间Tpk、每台压铸机需要生产的产品数量Qpk、熔化炉在初始时刻的保温量Mmsn、熔化炉初始时刻熔化烧嘴的状态Ssimsn(0代表待机,1代表在熔化)、各保温炉在初始时刻的保温量Mhsk、初始时是否有保温炉已经取料的标志(0为没取,1为已经取料)、已取料的时间和取料量开始时每台压铸机离下次加料的时间Tnextcsk为输入,建立以该整体成本最小为目标的优化模型,需要求解的参数组合包括该整体内熔化炉的熔化速度Vmn、熔化炉保温区的上限保温质量Mmupn、熔化炉保温区的下限保温质量Mmlown、各保温炉从熔化炉每次取料的时刻与质量j=1,2,3…,各保温炉上限的保温高度Mhupk、各保温炉下限的保温高度Mhlowk。
该优化模型的目标函数如下:
约束条件为:
Vmn∈[Vmminn,Vmmaxn]
Mmminn≤Mmlown<Mmupn≤Mmmaxn
Mhmink≤Mhlowk<Mhupk≤Mhmaxk
其中:Con为该整体的总成本,Cmn为该整体中熔化炉的使用成本,Ctk和∑Ctk分别为第k台保温炉的铝液运输成本和该整体中所有保温炉铝液运输成本之和(∑为求和符号),Cpk和∑Cpk分别为第k台保温炉的使用成本和该整体中所有保温炉的使用成本之和,Cpk和∑Cpk分别为第k台保温炉的因铝液缺少导致的暂停成本和该整体中所有保温炉因铝液缺少导致的暂停的成本之和,Fhk为车间的某一保温炉,为该整体中保温炉的集合,Vmminn,Vmmaxn为该整体中熔化炉最小和最大的熔化速度,Mmminn,Mmmaxn为该整体中熔化炉保温区的最小和最大的保温质量,Mhmink,Mhmaxk为该整体中各保温炉的最小和最大保温质量。
Cmn的计算表达如下:
Cmn=Emn*pemn+Mon*pan
Emn=Emmn+Emhn+Emsn+Emcn
Emmn=Pmmn*Tmmn
Emsn=Pmsn*Tmsn
Emcn=ΔEmcn*Nmcn
其中:Emn,Mon,pemn,pan分别为该整体中熔化炉的能源消耗量、熔化炉的材料氧化烧损量、熔化炉的熔化用能源单价、熔化炉的熔化材料单价;Emmn,Emhn,Emsn,Emcn分别为该整体中熔化炉的熔化能耗、熔化炉保温区铝液保温能耗、熔化炉熔化待机能耗、熔化炉熔化状态切换能耗;Pmmn,Pmhn,Pmsn分别为该整体中熔化炉的熔化功率、熔化炉保温区保温功率、熔化炉熔化待机功率;Tmmn,Tmhn,Tmsn分别为该整体中熔化炉的熔化时间、保温区保温时间、熔化待机时间;ΔEmcn,Nmcn分别为该整体中熔化炉的一次熔化状态切换的能量损失和切换次数,为Pmmn的系数,Mmhn为熔化炉保温区溶液质量,为Pmhn的系数,为Mon的系数。
Ctk的计算表达如下:
Ctk=Etk*pet+Htk*pht
Htk=Ntk*ct
其中:Etk,pet分别为第k台保温炉的铝液运输能源消耗量和运输用能源的单价,Htk,pht分别为第k台保温炉的运输人力使用量和单位人力成本,Ntk为第k台保温炉的铝液运输次数,为第k台保温炉的一次铝液运输消耗的能量,ct为人力使用量的折算系数。
Chk的计算表达如下:
Chk=Ehk*pehk
Cpk的计算表达如下:
Cpk=Tmspk*pmspk
其中:Tmspk,pmspk分别为第k台保温炉的暂停的时长和暂停的处罚系数。
需要说明是,上述各计算表达式的系数中,能源成本、材料成本、人力成本、处罚系数是已知的,Pmsn,ΔEmcn,ct,这些参数是静态参数,可以通过现场收集或者实验获得;Tmmn,Tmhn,Tmsn,Nmcn,Ntk,Thk,Tmspk需要根据不同的参数组合进行计算,这些参数采用仿真的方式获得;Tmhn等于熔化炉运行的时间,即从开始到所有生产任务结束,Thk为第k台保温炉运行的时间,即从开始到该保温炉对应的生产任务结束;Tmmn,Tmsn,Nmcn,Ntk,Tmspk需要通过仿真的得到,分别是熔化和保温质量曲线在时间上的积分,也需要通过仿真得到,仿真的一些设定如下:
①运输设备和运输人员数量足够多,各运输包的型号一样,运输包最大容量为Mtmax。
②保温炉正常需要加料的时刻为Mhk降到Mhlowk附近(根据实际情况,有等于、大于和小于,波动小于±Mpk)的时刻再加上半个Tpk,每次加料要加到Mhupk,当需要加料的量大于运输包的最大容量,就只加运输包的容量。
③从熔炉取料到给各保温炉加料的间隔(包括除气除渣等操作)为Ttcn,k。
④保温炉的保温质量低于设定的Mhlowk值时,也可以继续生产,当保温质量在[Mhmink,Mhmink+Mpk)时,生产暂停,直到铝液供上。
⑤当两台及以上保温炉正常取料的时间间隔小于一个设定的值ΔTc时,且它们的总加料量不超过运输包最大容量Mtmax时,运输包可以一次给多台保温炉加料,保温炉加料的时间比正常的有一些提前,否者就是运输包一次给一台保温炉加料。
⑥如果Mmsn大于或等于Mmupn时,工作开始时的熔炉熔化的初始状态为待机,否者熔化的初始状态为熔化。
⑦取料时,熔炉的当时的最大可用保温质量(Mmhn-Mmminn)小于这次需要加料的量,如果处于熔化状态,取料和加料延迟,等待直到最大可用保温质量满足这一次的需求,如果处于熔化待机状态,则这组参数直接不合格,返回不合格标志flaguq=0。
⑧取料和加料是一瞬间完成的,两个取料或者加料的最小间隔为一个仿真步长。
⑨如图3(b)所示,某一生产计划生产结束时,对应保温炉的质量需要降到Mhmink。
⑩如图3(a)所示,当所有生产结束后,熔炉的保温质量需要降到Mmminn。
仿真输入的参数如下:
①熔炉的参数(Vmn,Mmminn,Mmupn,Mnlown,Mmsn),熔炉熔化起始的状态Ssims。
②运输包的最大容量Mtmax、取加料的时间Ttcn,k、合并运输的间隔ΔTc、开始时是否有保温炉已经取料的标志(0为没取,1为已经取料),已取料的时间取料的量开始时每台压铸机离下次加料的时间Tnextcsk。
仿真输出的参数如下:
①这组参数是否合格标志flaguq,0为不合格,1为合格。
④如果所有生产结束,熔化炉的停止熔化时间Temn,使得熔化炉保温质量回到Mmminn。
⑤仿真结束时熔化炉保温质量Mmsen、熔化的状态Ssimen。
⑦各保温炉结束时的保温质量Mhek。
如图4所示,仿真的具体流程如下:
1.先定义各种需要的变量和变量初始化,具体流程如下:
2.接着先判断是否满足仿真结束的条件,具体流程如下:
3.如果满足,则返回各种需要的参数和结束仿真。
4.否则按照规则更新熔炉状态和熔炉保温质量,接着更新保温炉保温质量,然后根据规则判断是否取料,最后根据规则判断给各保温炉加料,具体流程如下:
5.重复执行2~4,直到满足仿真结束条件。
(2)采用遗传算法对步骤(1)的模型进行求解,求解的流程如图2所示,先求解一组该整体最小成本的参数,将此运行到该整体中一台压铸机生产任务完成,然后以上述结束的状态为初始状态,求解另外一组最小成本参数,将此运行到另一台压铸机生产任务完成,重复这个动作,直到所有生产任务结束。
下面以一个具体的案例说明本发明的实施过程,实施过程如图5所示,其中某一整体由一台熔炉和两个压铸单元组成,需要生产100件A产品(浇注质量3千克,周期60秒)和生产100件B产品(浇注质量2千克,周期50秒),需要选择合适的熔炉、配送和保温炉参数。
先收集相关的参数,Vmminn,Vmmaxn分别为100和600kg/h,Mmminn,Mmmaxn分别为100和500kg,Mhmink,Mhmaxk分别为(100,100)kg和(300,300)kg,为(0.5,0.6)kWh,Ttcn,k为(420,450)秒,两个保温炉的能源单价都为0.48kWh/元,两个压铸单元的惩罚系数都为200元/小时,两个保温炉的保温系数都为0.054和5.64,熔化炉的能源单价0.48kWh/元,熔化材料的当家15元/kg,熔化状态切换一次的能源损失14kWh,熔化炉待机功率27kW,熔化的三个功率系数分别为(20,1.9,500),熔化炉保温的两个系数为(0.054,5.64),氧化烧损的两个系数为(0.727,-0.643),运输用能源单价1元/kWh,人力资源单价5元/人次,人力成本折算系数1,运输包最大运输量400kg,合并加料的时间120秒,熔炉的起始量150kg,保温炉的起始量(140,160)kg,熔炉起始时待机,开始时没有取料的,距离下次加料时间正常。
将上述参数代入到优化模型中,将上述优化模型在MATLAB中进行编程,将目标函数作为优化函数,约束作为限制,用MATLAB自带的遗传算法工具箱来寻找最低成本的参数组合。如图6所示,输出的一组近似最优的参数为熔化速度570,熔化炉最低和最高保温质量为(260,340),保温炉最高和最低的保温高度为最低(100,120),最高(220,190),停止熔化时间为2212秒,最小成本250.2,A产品的取料时间和取料量为[438,2726,131]和[119,117,24],B产品取料时间和取料量[981,2726]和[86,54]。
上述对实施例的描述是为便于本技术领域的普通技术人员能理解和应用本发明。熟悉本领域技术的人员显然可以容易地对上述实施例做出各种修改,并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此,本发明不限于上述实施例,本领域技术人员根据本发明的揭示,对于本发明做出的改进和修改都应该在本发明的保护范围之内。
Claims (7)
1.一种面向铝压铸的熔化、配送和保温参数集成的节能优化方法,其针对车间内以某一台熔化炉为主的生产整体,该整体内还包括N台由熔化炉提供合金的压铸机及其对应的N台保温炉,N为大于1的自然数;其特征在于:
首先,获取已知的相关生产信息作为输入,包括该整体中压铸机每次的浇注质量、压铸机每次浇注的间隔时间、压铸机所需生产的产品数量、熔化炉在初始时刻的保温量、初始时刻熔化炉熔化烧嘴的状态、保温炉在初始时刻的保温量、初始时保温炉是否已取料及取料时刻和取料量、初始时压铸机离下次加料的时间;
然后,建立以整体的总成本最小为目标的优化模型如下:
min{Co=Cm+∑Ctk+∑Chk+∑Cpk}
其中:Co为该整体的总成本,Cm为该整体中熔化炉的使用成本,Ctk为该整体中第k台保温炉的铝液运输成本,Chk为第k台保温炉的使用成本,Cpk为第k台保温炉因铝液缺少导致的暂停成本,k为自然数且1≤k≤N;
最后,对上述优化模型进行求解,得到该整体中熔化炉的熔化速度、熔化炉保温区的上限保温质量和下限保温质量、保温炉从熔化炉每次取料的时刻和质量、保温炉的上限保温质量和下限保温质量。
2.根据权利要求1所述的节能优化方法,其特征在于:所述优化模型的约束条件如下:
Vm∈[Vmmin,Vmmax]
Mmmin≤Mmlow<Mmup≤Mmmax
Mhmink≤Mhlowk<Mhupk≤Mhmaxk
其中:Vm为熔化炉的熔化速度,Vmmin和Vmmax分别为熔化炉的最小熔化速度和最大熔化速度,Mmup和Mmlow分别为熔化炉保温区的上限保温质量和下限保温质量,Mmmin和Mmmax分别为熔化炉保温区的最小保温质量和最大保温质量,Mhupk和Mhlowk分别为第k台保温炉的上限保温质量和下限保温质量,Mhmink和Mhmaxk分别为第k台保温炉的最小保温质量和最大保温质量。
3.根据权利要求1所述的节能优化方法,其特征在于:采用遗传算法对优化模型进行求解,求解过程中:先求得一组整体最小成本的优化参数,在该优化参数基础上运行使该整体中一台压铸机生产任务完成,然后以上述结束的状态为初始状态,求解另外一组整体最小成本的优化参数,在该优化参数基础上运行直到另一台压铸机生产任务完成,重复上述操作,直至所有压铸机生产任务完成。
4.根据权利要求1所述的节能优化方法,其特征在于:所述熔化炉使用成本Cm的表达式如下:
Cm=Em*pem+Mo*pa
Em=Emm+Emh+Ems+Emc
Emm=Pmm*Tmm
Emh=∫0 TmhPmh
Ems=Pms*Tms
Emc=ΔEmc*Nmc
Pmh=cmh1*Mmh+cmh2
其中:Em为熔化炉的熔化能源消耗量,pem为熔化炉的熔化用能源单价,Mo为熔化炉的材料氧化烧损量,pa为熔化炉的熔化材料单价,Emm为熔化炉的熔化能耗,Emh为熔化炉保温区的铝液保温能耗,Ems为熔化炉的熔化待机能耗,Emc为熔化炉的熔化状态切换能耗,Pmm为熔化炉的熔化功率,Tmm为熔化炉的熔化时间,Pmh为熔化炉保温区的保温功率,Tmh为熔化炉保温区的保温时间,Pms为熔化炉的熔化待机功率,Tms为熔化炉的熔化待机时间,ΔEmc为熔化炉一次熔化状态切换的能量损失,Nmc为熔化炉的熔化状态切换次数,Vm为熔化炉的熔化速度,cmm1、cmm2、cmm3、cmh1、cmh2、cmo1、cmo2均为给定的系数,Mmh为熔化炉保温区的溶液质量。
7.根据权利要求1所述的节能优化方法,其特征在于:所述第k台保温炉因铝液缺少导致的暂停成本Cpk表达式如下:
Cpk=Tmspk*pmsp
其中:Tmspk为第k台保温炉的暂停时长,pmsp为暂停的处罚系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110491310.8A CN113111540B (zh) | 2021-05-06 | 2021-05-06 | 面向铝压铸的熔化、配送和保温参数集成的节能优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110491310.8A CN113111540B (zh) | 2021-05-06 | 2021-05-06 | 面向铝压铸的熔化、配送和保温参数集成的节能优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113111540A CN113111540A (zh) | 2021-07-13 |
CN113111540B true CN113111540B (zh) | 2022-07-12 |
Family
ID=76721284
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110491310.8A Active CN113111540B (zh) | 2021-05-06 | 2021-05-06 | 面向铝压铸的熔化、配送和保温参数集成的节能优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113111540B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5762680A (en) * | 1995-11-22 | 1998-06-09 | Norsk Hydro A.S. | Distribution system for molten magnesium |
CN107423512A (zh) * | 2017-07-28 | 2017-12-01 | 广东省智能制造研究所 | 一种大型压铸机尾板的轻量化设计方法 |
CN111325408A (zh) * | 2020-03-09 | 2020-06-23 | 浙江大学 | 一种面向铝压铸熔炉的工艺参数节能优化方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015151369A1 (ja) * | 2014-03-31 | 2015-10-08 | アイシン軽金属株式会社 | アルミニウム合金及びダイカスト鋳造方法 |
CN110059842A (zh) * | 2018-01-19 | 2019-07-26 | 武汉十傅科技有限公司 | 一种考虑熔炉和砂型大小的铸造车间生产计划优化方法 |
-
2021
- 2021-05-06 CN CN202110491310.8A patent/CN113111540B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5762680A (en) * | 1995-11-22 | 1998-06-09 | Norsk Hydro A.S. | Distribution system for molten magnesium |
CN107423512A (zh) * | 2017-07-28 | 2017-12-01 | 广东省智能制造研究所 | 一种大型压铸机尾板的轻量化设计方法 |
CN111325408A (zh) * | 2020-03-09 | 2020-06-23 | 浙江大学 | 一种面向铝压铸熔炉的工艺参数节能优化方法 |
Non-Patent Citations (3)
Title |
---|
An Internet of Things-enabled model-based approach to improving the energy efficiency of aluminum die casting processes;Weipeng Liu等;《Energy》;20200430;全文 * |
铝合金压铸行业节能评估要点剖析;夏维等;《资源节约与环保》;20160625(第06期);全文 * |
铝合金熔体制备供应能耗计算与分析;王甜;《万方》;20151231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113111540A (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105160451B (zh) | 一种含电动汽车的微电网多目标优化调度方法 | |
CN109858759A (zh) | 一种工业园区综合能源平衡调度方法 | |
CN103840457B (zh) | 考虑电动汽车充放电影响的配电网内dg优化配置方法 | |
CN102170129B (zh) | 一种基于大规模风电并网的电力系统低碳调度方法及其装置 | |
CN105811409B (zh) | 一种含电动汽车混合储能系统的微网多目标运行调度方法 | |
CN109995091B (zh) | 一种考虑预测误差的交直流混合微电网经济调度方法 | |
CN107546781B (zh) | 基于pso改进算法的微电网多目标运行优化方法 | |
CN103986194B (zh) | 一种独立型微网优化配置方法和装置 | |
CN107590744B (zh) | 考虑储能和无功补偿的主动配电网分布式电源规划方法 | |
CN105515110B (zh) | 一种电动汽车有序充电实时控制系统 | |
CN105868844A (zh) | 一种含电动汽车混合储能系统的微网多目标运行调度方法 | |
CN105139147A (zh) | 微电网系统的经济调度方法 | |
CN105207259A (zh) | 基于能量管理的并网状态下的微电网系统调度方法 | |
CN112464477A (zh) | 计及需求响应的多能耦合综合能源运行仿真方法 | |
CN105870950A (zh) | 一种智能楼宇可调资源多目标协调系统及优化控制方法 | |
CN105117797A (zh) | 基于孤岛划分的微网日前经济最优调度方法 | |
CN110165666A (zh) | 一种基于igdt的主动配电网调度方法 | |
CN110048459A (zh) | 一种基于弃风利用的塔筒电梯储能系统优化配置方法 | |
CN101989743B (zh) | 基于直流潮流的节能发电调度优化方法 | |
CN105207207B (zh) | 基于能量管理的孤网状态下的微电网系统调度方法 | |
CN110084394A (zh) | 一种考虑温控设备和电动汽车的工业用户智能用电优化方法 | |
CN112700066A (zh) | 一种电-热综合能源系统调度最佳时间尺度配合方法 | |
CN117096868A (zh) | 一种考虑多种柔性负荷和电动汽车的微电网能量调度方法 | |
CN113111540B (zh) | 面向铝压铸的熔化、配送和保温参数集成的节能优化方法 | |
CN115688448A (zh) | 一种考虑共享储能的多区域综合能源系统优化调度方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |