• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于動(dòng)態(tài)成本卷積的復(fù)雜產(chǎn)品批產(chǎn)路徑優(yōu)化問(wèn)題建模與求解研究

    2024-11-23 00:00:00楊麗穎楊銳意崔新豪張思悅陳練肖依永
    關(guān)鍵詞:優(yōu)化模型路徑優(yōu)化

    摘要: 具有復(fù)雜工藝路線的裝備產(chǎn)品,其生產(chǎn)系統(tǒng)設(shè)計(jì)的待決策變量之間通常相互耦合、相互影響,對(duì)產(chǎn)品最終的批產(chǎn)成本產(chǎn)生著影響,構(gòu)成了較為復(fù)雜的組合優(yōu)化問(wèn)題。對(duì)此,提出一種基于動(dòng)態(tài)成本卷積的復(fù)雜產(chǎn)品批產(chǎn)路徑優(yōu)化模型。該模型考慮了復(fù)雜產(chǎn)品工藝路線的各節(jié)點(diǎn)均可能存在多種生產(chǎn)設(shè)計(jì)選項(xiàng),在各選項(xiàng)對(duì)應(yīng)著不同生產(chǎn)投入和批量效率的情況下,以最終交付產(chǎn)品的批產(chǎn)成本最小化為目標(biāo)輸出生產(chǎn)系統(tǒng)設(shè)計(jì)和成本卷積路徑。同時(shí),建立了混合整數(shù)規(guī)劃模型,并通過(guò)對(duì)模型中非線性成分實(shí)施可控精度線性化,將模型換化為可最優(yōu)求解的線性模型。最后,針對(duì)某連續(xù)生產(chǎn)行業(yè)設(shè)計(jì)了小、中、大規(guī)模實(shí)驗(yàn)算例,驗(yàn)證了所提模型的可行性、合理性和求解效率。

    關(guān)鍵詞: 批產(chǎn)成本; 成本卷積; 路徑優(yōu)化; 優(yōu)化模型

    中圖分類號(hào): C 94

    文獻(xiàn)標(biāo)志碼: A

    DOI:10.12305/j.issn.1001-506X.2024.06.19

    Research on modeling and solution of complex product batch production routing optimization based on dynamic cost convolution

    YANG Liying1, YANG Ruiyi2, CUI Xinhao2, ZHANG Siyue2, CHEN Lian3, XIAO Yiyong2,*

    (1. National Key Laboratory of Particle Transport and Separation Technology, Research Institute of Physical and Chemical Engineering of Nuclear Industry, Tianjin 300180, China;

    2. School of Reliability and Systems Engineering, Beihang University, Beijing 100191, China;

    3. The 714 Research Institute of China State Shipbuilding Corporation Limited, Beijing 100101, China)

    Abstract: Equipment products with complex process routes often have decision variables in their production system design that are coupled and affect with each other, affecting the final batch production cost of the product and forming a complex combinatorial optimization problem. In this regard, a complex product batch production routing optimization model based on dynamic cost convolution is proposed. The proposed model considers that there may be multiple production design options at each node of the complex product process routing, and each option corresponds to different production inputs and batch efficiency. The goal is to minimize the batch production cost of the final delivered product and output the production system design and cost convolution routing. At the same time, a mixed integer programming model is established, and by implementing controllable precision linearization on the nonlinear components in the model, the model is transformed into a linear model that can be optimally solved. Finally, small, medium, and large-scale experimental examples are designed for a certain continuous production industry to verify the feasibility, rationality, and solution efficiency of the proposed model.

    Keywords: batch production cost; cost convolution; routing optimization; optimization model

    0 引 言

    生產(chǎn)系統(tǒng)優(yōu)化設(shè)計(jì)是生產(chǎn)運(yùn)作管理領(lǐng)域的重要研究?jī)?nèi)容,圍繞生產(chǎn)效率、產(chǎn)品質(zhì)量、交付時(shí)間和降低成本等多方面的評(píng)價(jià)目標(biāo),長(zhǎng)期以來(lái)在工業(yè)界和學(xué)術(shù)界開(kāi)展著廣泛的理論研究和應(yīng)用實(shí)踐。在裝備研制與生產(chǎn)領(lǐng)域,生產(chǎn)成本是產(chǎn)品壽命周期費(fèi)用的重要組成部分,也是產(chǎn)品經(jīng)濟(jì)性設(shè)計(jì)的重要內(nèi)容。當(dāng)某新研產(chǎn)品完成設(shè)計(jì)定型進(jìn)入批產(chǎn)階段后,批產(chǎn)成本就成為生產(chǎn)系統(tǒng)設(shè)計(jì)的一項(xiàng)重要優(yōu)化目標(biāo),也是評(píng)價(jià)生產(chǎn)系統(tǒng)優(yōu)劣性的重要指標(biāo)。

    產(chǎn)品的批產(chǎn)成本受到諸多因素的影響,例如生產(chǎn)批量計(jì)劃、供應(yīng)鏈管理、生產(chǎn)設(shè)施選址、生產(chǎn)設(shè)備選型、車間布局優(yōu)化、設(shè)備維修保障策略,以及生產(chǎn)工藝流程設(shè)計(jì)等。如何集成這些因素進(jìn)行綜合優(yōu)化設(shè)計(jì),以最大程度地降低產(chǎn)品批產(chǎn)成本,已經(jīng)成為生產(chǎn)系統(tǒng)優(yōu)化設(shè)計(jì)的主要內(nèi)容。不同行業(yè)產(chǎn)品的批產(chǎn)成本影響因素通常具有較大的差異性,因此通過(guò)研究生產(chǎn)系統(tǒng)的特點(diǎn),進(jìn)而有針對(duì)性地開(kāi)展生產(chǎn)系統(tǒng)優(yōu)化設(shè)計(jì),已成為當(dāng)前批產(chǎn)成本優(yōu)化的一個(gè)研究熱點(diǎn)方向。

    多級(jí)經(jīng)濟(jì)批量問(wèn)題是生產(chǎn)批量?jī)?yōu)化方面的典型問(wèn)題[1。該問(wèn)題以批產(chǎn)成本為優(yōu)化目標(biāo),研究核心企業(yè)與多級(jí)供應(yīng)商之間進(jìn)行集成優(yōu)化的生產(chǎn)經(jīng)濟(jì)批量模型及求解算法[2。文獻(xiàn)[3]研究了帶約束的經(jīng)濟(jì)批量問(wèn)題及動(dòng)態(tài)規(guī)劃算法。文獻(xiàn)[4-5]研究了經(jīng)濟(jì)批量與學(xué)習(xí)效應(yīng)之間的影響關(guān)系。文獻(xiàn)[6]研究了單元生產(chǎn)中研制訂單和批產(chǎn)訂單的混合調(diào)度問(wèn)題,考慮了隨機(jī)事件的動(dòng)態(tài)調(diào)度算法。文獻(xiàn)[7]從降本增效、提升產(chǎn)業(yè)化能力方面對(duì)武器平臺(tái)電子產(chǎn)品的批量化生產(chǎn)過(guò)程進(jìn)行探討。文獻(xiàn)[8]研究了針對(duì)制造業(yè)生產(chǎn)車間多品種、多階段生產(chǎn)問(wèn)題的生產(chǎn)邏輯模型和優(yōu)化方案。文獻(xiàn)[9]對(duì)相關(guān)的數(shù)學(xué)模型進(jìn)行了系統(tǒng)性綜述。

    生產(chǎn)線優(yōu)化設(shè)計(jì)也是現(xiàn)有文獻(xiàn)研究批產(chǎn)成本優(yōu)化的一個(gè)重要方面。文獻(xiàn)[10]研究了面向生產(chǎn)快速反應(yīng)、高效作業(yè)和低成本綜合目標(biāo)的復(fù)雜電子產(chǎn)品柔性生產(chǎn)線建設(shè)的優(yōu)化問(wèn)題。文獻(xiàn)[11]研究了批產(chǎn)制造項(xiàng)目的進(jìn)度管理,采用啟發(fā)式算法基于資源重新分配優(yōu)化了關(guān)鍵生產(chǎn)鏈問(wèn)題。文獻(xiàn)[12]研究了低成本批產(chǎn)商業(yè)通信衛(wèi)星的任務(wù)特點(diǎn)和設(shè)計(jì)要求,提出了面向批生產(chǎn)低成本衛(wèi)星的元器件保障模式。在生產(chǎn)線流程優(yōu)化方面,日本的單元式生產(chǎn)被認(rèn)為是傳統(tǒng)福特流水生產(chǎn)線模式的一種新的替代,在降低生產(chǎn)成本和提升生產(chǎn)效率方面有顯著優(yōu)勢(shì)[13-15

    文獻(xiàn)[16]研究了生產(chǎn)設(shè)施布局問(wèn)題以及動(dòng)態(tài)設(shè)施布局問(wèn)題[17,提出了在考慮外部需求變化的情況下,如何動(dòng)態(tài)調(diào)整生產(chǎn)設(shè)施布局與工藝路線流程,以降低產(chǎn)品批產(chǎn)成本為目標(biāo)建立優(yōu)化模型。文獻(xiàn)[18]研究了工廠設(shè)施的連續(xù)選址問(wèn)題,以降低運(yùn)輸成本為優(yōu)化目標(biāo),為之建立了混合整數(shù)規(guī)劃模型。文獻(xiàn)[19]研究了面向利潤(rùn)最大化的批量生產(chǎn)計(jì)劃的非陣列排序優(yōu)化模型。文獻(xiàn)[20]研究了考慮產(chǎn)能擴(kuò)張的生產(chǎn)組合計(jì)劃優(yōu)化問(wèn)題。

    對(duì)于串行生產(chǎn)系統(tǒng),由設(shè)備故障引發(fā)的生產(chǎn)系統(tǒng)暫停對(duì)生產(chǎn)效率有較大影響,因此設(shè)備維修策略是生產(chǎn)系統(tǒng)設(shè)計(jì)的重要內(nèi)容[21-22。其中,文獻(xiàn)[21]提出了帶緩沖區(qū)的串行生產(chǎn)系統(tǒng)預(yù)防性維護(hù)決策優(yōu)化模型及啟發(fā)式求解算法,文獻(xiàn)[22]提出了一種面向串聯(lián)生產(chǎn)系統(tǒng)設(shè)備維護(hù)的考慮在線決策與緩沖分配的聯(lián)合優(yōu)化模型。文獻(xiàn)[23]提出了設(shè)備狀態(tài)與預(yù)防性維護(hù)周期相結(jié)合的預(yù)防維護(hù)策略,以單位時(shí)間總質(zhì)量成本最小化為目標(biāo)建立優(yōu)化模型。相關(guān)的綜述可參考文獻(xiàn)[23],該文獻(xiàn)分析了生產(chǎn)系統(tǒng)中設(shè)備維修管理、統(tǒng)計(jì)過(guò)程控制及經(jīng)濟(jì)生產(chǎn)批量決策的相互影響和作用機(jī)制,對(duì)相關(guān)應(yīng)用問(wèn)題研究現(xiàn)狀進(jìn)行了梳理和分析。

    復(fù)雜產(chǎn)品批產(chǎn)成本優(yōu)化領(lǐng)域還有一個(gè)被稱之為組裝路徑優(yōu)化問(wèn)題(assembly routing problem, ARP)的經(jīng)典問(wèn)題[24-26。ARP考慮了在復(fù)雜產(chǎn)品具有多層零件結(jié)構(gòu)、每種零件有多個(gè)供應(yīng)商、每個(gè)供應(yīng)商有多個(gè)供應(yīng)地點(diǎn)/倉(cāng)庫(kù)的場(chǎng)景下,如何優(yōu)化物流供應(yīng)路徑,在滿足多周期客戶交付要求的前提下實(shí)現(xiàn)生產(chǎn)成本最小化。當(dāng)ARP向后延伸到考慮多客戶的交貨地點(diǎn)和緩存?zhèn)}庫(kù)并考慮終端產(chǎn)品運(yùn)輸路徑而優(yōu)化生產(chǎn)批量時(shí),就轉(zhuǎn)化為另一種形式的生產(chǎn)路徑優(yōu)化問(wèn)題(production routing problem, PRP)[27-30、倉(cāng)儲(chǔ)路徑優(yōu)化問(wèn)題(inventory routing problem, IRP)[31,以及生產(chǎn)與分運(yùn)集成調(diào)度(integrated production and distribution scheduling, IPDS)問(wèn)題[32-33。

    現(xiàn)有文獻(xiàn)中還包括有較多的針對(duì)具體行業(yè)的生產(chǎn)系統(tǒng)優(yōu)化案例研究,包括制氫行業(yè)的生產(chǎn)路徑優(yōu)化問(wèn)題[34,以及從成本、碳排放和風(fēng)險(xiǎn)等角度建立的制氫供應(yīng)鏈多目標(biāo)優(yōu)化模型[35和相關(guān)綜述文獻(xiàn)(文獻(xiàn)[36])。文獻(xiàn)[37]研究了農(nóng)業(yè)生物行業(yè)的供應(yīng)鏈優(yōu)化模型。文獻(xiàn)[38]從生產(chǎn)規(guī)模、集成度、運(yùn)輸和供應(yīng)鏈配置等方面,研究了生物燃料生產(chǎn)行業(yè)的成本優(yōu)化問(wèn)題。文獻(xiàn)[39]從產(chǎn)量、成本、能源等7個(gè)角度,研究了食品供應(yīng)鏈的多目標(biāo)優(yōu)化模型。

    上述文獻(xiàn)關(guān)于批產(chǎn)成本優(yōu)化的研究主要集中針對(duì)典型問(wèn)題模型的局部區(qū)域,如生產(chǎn)經(jīng)濟(jì)批量、工序優(yōu)化、設(shè)施布局等。關(guān)于如何從生產(chǎn)系統(tǒng)整體設(shè)計(jì)角度優(yōu)化產(chǎn)品批產(chǎn)成本的相關(guān)理論模型研究還比較少,還缺乏面向復(fù)雜產(chǎn)品全生產(chǎn)過(guò)程和生產(chǎn)系統(tǒng)批產(chǎn)成本的優(yōu)化模型。對(duì)生產(chǎn)線的優(yōu)化研究目前還主要針對(duì)單一生產(chǎn)線,較少有針對(duì)多生產(chǎn)線耦合而組成的復(fù)雜生產(chǎn)系統(tǒng)的優(yōu)化研究。

    本文根據(jù)實(shí)際工程需要,提煉出了面向生產(chǎn)系統(tǒng)整體優(yōu)化的批產(chǎn)成本優(yōu)化問(wèn)題,建立了一種基于動(dòng)態(tài)成本卷積的產(chǎn)品批產(chǎn)路徑優(yōu)化模型。該模型考慮了產(chǎn)品批產(chǎn)過(guò)程中的重復(fù)性直接投入成本和生產(chǎn)系統(tǒng)維修消耗,以及非重復(fù)性的成本投入,如研發(fā)投入、生產(chǎn)線建設(shè)等各種一般性的成本要素,為具有復(fù)雜生產(chǎn)工藝路線的產(chǎn)品批產(chǎn)提供最優(yōu)化的成本卷積路徑。所提出的優(yōu)化模型采用了線性化設(shè)計(jì)技術(shù),可應(yīng)用商業(yè)求解器(如Cplex,Gurobi等)進(jìn)行直接求解。對(duì)某連續(xù)生產(chǎn)行業(yè)的模擬算例表明,對(duì)于中小規(guī)模問(wèn)題,該模型可在1 h計(jì)算時(shí)間內(nèi)獲得最優(yōu)解。對(duì)于大規(guī)模問(wèn)題,可在可接受時(shí)間內(nèi)獲得次優(yōu)可行解。

    1 基于制程的成本動(dòng)態(tài)卷積模型

    1.1 直接生產(chǎn)成本

    生產(chǎn)制程即生產(chǎn)制造過(guò)程,是指通過(guò)使用和消耗生產(chǎn)資源,如人工、設(shè)備、能源、輔料、外部產(chǎn)品或服務(wù)等,將輸入原材料或半成品轉(zhuǎn)換為輸出產(chǎn)品的一組有序過(guò)程,也稱為生產(chǎn)線,例如:半導(dǎo)體生產(chǎn)過(guò)程、精密件加工過(guò)程、成品組裝過(guò)程,以及專有生產(chǎn)服務(wù)支持(如特種搬運(yùn)、激光調(diào)制)等。圖1描繪了一個(gè)輸出某種產(chǎn)品或服務(wù)的生產(chǎn)制程概念模型。

    生產(chǎn)系統(tǒng)往往包含多條生產(chǎn)制程,共同組成網(wǎng)狀結(jié)構(gòu)的生產(chǎn)制程關(guān)系模型,如圖2所示。其中,制程3分別領(lǐng)用制程1和制程2的產(chǎn)品1和產(chǎn)品2,輸出產(chǎn)品3,然后被制程4領(lǐng)用,輸出產(chǎn)品4,再與生產(chǎn)制程5輸出的產(chǎn)品5一起被生產(chǎn)制程6領(lǐng)用,最后輸出終端產(chǎn)品6。

    以W表示生產(chǎn)制程的集合,O表示生產(chǎn)制程輸出產(chǎn)品的集合,pij表示制程i輸出產(chǎn)品j的每批次產(chǎn)量,sij為制程i對(duì)產(chǎn)品j的每批次用量,βij為生產(chǎn)制程i的產(chǎn)出產(chǎn)品j所占成本比例,其中i∈W, j∈O。以R表示制造資源集合,如人工、設(shè)備、外購(gòu)材料、能源等,qir表示生產(chǎn)制程i對(duì)制造資源r的每批次消耗量,vr表示制造資源r的成本單價(jià),其中i∈W, r∈R。以u(píng)ij表示生產(chǎn)制程i生產(chǎn)產(chǎn)品j的單位成本,zi表示生產(chǎn)制程i的單批次成本,那么,對(duì)每一個(gè)生產(chǎn)制程i,有如下成本卷積表達(dá)式:

    zi=∑j∈Osij·uij+∑r∈Rqir·vr,?i∈W

    pij·uijij·zi,?i∈W; j∈O(1)

    1.2 生產(chǎn)設(shè)備維修成本

    (1) 預(yù)防性維修費(fèi)用

    預(yù)防性維修是為了確保生產(chǎn)制程系統(tǒng)保持設(shè)計(jì)功能和性能,以及防止故障發(fā)生而采取的周期性檢查、潤(rùn)滑、備件更換等措施,例如各種小修、定檢、中修、大修等。預(yù)防性維修活動(dòng)產(chǎn)生的費(fèi)用稱為預(yù)防性維修費(fèi)用。預(yù)防性維修是計(jì)劃性的,其費(fèi)用可以進(jìn)行預(yù)先估算。

    預(yù)防性維修費(fèi)用分析如下:令W表示生產(chǎn)系統(tǒng)的生產(chǎn)制程集合;以Fi表示對(duì)生產(chǎn)制程i的預(yù)防性維修項(xiàng)目集合,li表示生產(chǎn)制程i的單批次生產(chǎn)時(shí)間,其中i∈W;以hik表示預(yù)防性維修項(xiàng)目的工作時(shí)間間隔,h′ik表示預(yù)防性維修項(xiàng)目的日歷時(shí)間間隔,其中i∈W, k∈Fi;以cik表示單次預(yù)防性維修活動(dòng)的費(fèi)用,其中i∈W, k∈Fi。在給定的計(jì)劃時(shí)間T內(nèi),對(duì)于生產(chǎn)制程i,若其生產(chǎn)批次數(shù)為xi,那么所產(chǎn)生的預(yù)防性維修費(fèi)用為

    c′=∑i∈W∑k∈Fimaxlixihik,Th′ikcik(2)

    式中:max函數(shù)表示按先到為主原則,取按日歷時(shí)間和按工作時(shí)間計(jì)算所得的大者為實(shí)際維修次數(shù)。通常認(rèn)為日歷時(shí)間間隔h′ik遠(yuǎn)大于工作時(shí)間間隔h′ik,或者設(shè)備被共享使用時(shí),維修次數(shù)主要按工作時(shí)間計(jì)算。因此,第i生產(chǎn)制程的預(yù)防性維修費(fèi)用估算式可簡(jiǎn)寫為

    c′=∑i∈W∑k∈Filixihikcik(3)

    (2) 故障性維修費(fèi)用

    故障性維修是指生產(chǎn)系統(tǒng)發(fā)生故障后或者出現(xiàn)故障征兆等情況下進(jìn)行的非計(jì)劃性維修。故障性維修活動(dòng)產(chǎn)生的費(fèi)用稱為故障性維修費(fèi)用。故障性維修是非計(jì)劃性的,事先未必知道在什么時(shí)候要發(fā)生何種維修活動(dòng)。但根據(jù)生產(chǎn)設(shè)備的可靠性指標(biāo)以及故障模式分布規(guī)律,也可以估算出給定期間的故障性維修期望費(fèi)用。故障性維修費(fèi)用分析如下:以Gi表示生產(chǎn)制程的設(shè)備集合,以nik表示設(shè)備的安裝數(shù)量,其中i∈W, k∈Gi;以θik表示設(shè)備的可靠性指標(biāo)平均故障間隔時(shí)間(mean time between failure,MTBF),i∈W, k∈Gi;以Eik表示設(shè)備的故障模式集合,以λikj表示故障模式的占比率,以c′ikj表示設(shè)備發(fā)生故障的修復(fù)費(fèi)用。其中,i∈W, k∈Gi, j∈Eik。

    c″=∑i∈W∑k∈Fixilinikθikc′ik(4)

    綜上,生產(chǎn)系統(tǒng)維修費(fèi)用可表示為

    c=c′+c″=∑i∈W∑k∈Fixilinikθikc′ik+∑i∈W∑k∈Filixihikcik=

    ∑i∈Wxi∑k∈Filinikθikc′ik+∑k∈Filihikcik(5)

    令ci表示生產(chǎn)制程i的一批次生產(chǎn)對(duì)應(yīng)的平均生產(chǎn)系統(tǒng)維修成本,表示為

    ci=∑k∈Filinikθikc′ik+∑k∈Filihikcik(6)

    1.3 非重復(fù)性成本分?jǐn)?/p>

    在生產(chǎn)成本中考慮非重復(fù)性投入成本,如產(chǎn)品前期研發(fā)費(fèi)用、生產(chǎn)系統(tǒng)建設(shè)費(fèi)用、技改費(fèi)用等。以bi表示生產(chǎn)制程i的非重復(fù)性成本,按總生產(chǎn)批次數(shù)均攤?cè)胫瞥躺a(chǎn)成本,即bi/xi。其中,xi是制程i的生產(chǎn)批次數(shù)。這樣,在納入非重復(fù)成本均攤和生產(chǎn)維修費(fèi)用的情況下,式(1)成本卷積公式改寫為

    zi=∑j∈Osij·uij+∑r∈Rqir·vr+ci+bi/xi,?i∈W

    pij·uijij·zi, ?i∈W;j∈O(7)

    2 基于動(dòng)態(tài)成本卷積的批產(chǎn)路徑優(yōu)化模型

    2.1 符號(hào)定義

    (1) 模型參數(shù)

    W:生產(chǎn)制程集合,i∈W, i′∈W

    O:自制產(chǎn)品(或服務(wù))集合,j∈O

    R:生產(chǎn)資源集合,如人工、設(shè)備、材料、能源等,r∈R

    pij:制程i產(chǎn)出產(chǎn)品j的每批次產(chǎn)量

    p′ij:0/1參數(shù),表示制程i是否產(chǎn)出產(chǎn)品j

    βij:制程i產(chǎn)出產(chǎn)品j的成本占比,即當(dāng)產(chǎn)出多種產(chǎn)品時(shí)不同產(chǎn)品之間的成本占比,有∑j∈Oβij=1且βij≤p′ij,M為一個(gè)大數(shù)

    qir:制程i每批次對(duì)生產(chǎn)資源r的消耗量

    sij:制程i每生產(chǎn)批次領(lǐng)用產(chǎn)品j的數(shù)量

    s′ij:0/1參數(shù),表示制程i是否領(lǐng)用產(chǎn)品j

    vr:生產(chǎn)資源r的成本單價(jià)

    dj:產(chǎn)品j的外部需求量

    ci:制程i每批次設(shè)備維修費(fèi)用,見(jiàn)式(6)

    bi:制程的固定費(fèi)用(非重復(fù)性投入的研發(fā)費(fèi)、生產(chǎn)線建設(shè)費(fèi)、技改費(fèi)等)

    M:一個(gè)大數(shù)

    (2) 決策變量

    x′i:0/1變量,表示制程i是否啟用

    πii′j:0/1變量,表示成本卷積路徑,即制程i是否領(lǐng)用制程i′產(chǎn)出的產(chǎn)品j

    θii′j:非負(fù)連續(xù)變量,表示卷積的成本值,即制程i領(lǐng)用制程i′產(chǎn)出產(chǎn)品j的成本

    xi:非負(fù)整數(shù)變量,表示制程i的生產(chǎn)批次數(shù)

    zi:非負(fù)連續(xù)變量,表示制程i的單批次成本

    uij:非負(fù)連續(xù)變量,表示制程i產(chǎn)出產(chǎn)品j的每單位成本

    2.2 數(shù)學(xué)規(guī)劃模型

    (1)目標(biāo)函數(shù)

    min Total_Cost=∑i∈W∑j∈Odj·uij(8)

    (2)約束條件

    設(shè)定x與x′的關(guān)系:

    x′i≤xi,?i∈W

    xi≤Mx′i,?i∈W (9)

    (10)

    設(shè)定各制程的生產(chǎn)批次數(shù)剛好滿足需求:

    ∑i∈Wxipij≥∑i∈Wxisij+dj,?j∈O

    ∑i∈W(xi-1)pij≤∑i∈Wxisij+dj,?j∈O (11)

    (12)

    設(shè)定產(chǎn)品的成本卷積路徑:

    πii′j≤s′ijp′ij(xi+x′i)/2,?i∈W;i′∈W;j∈O

    ∑i′∈Wπii′j≥s′ijx′i,?i∈W;j∈O (13)

    (14)

    優(yōu)化產(chǎn)品的卷積成本:

    θii′j≤Mπii′j, ?i∈W;i′∈W;j∈O

    θii′j≥sijuij-M(1-πii′j), ?i∈W;i′∈W;j∈O

    θii′j≤sijuij+M(1-πii′j), ?i∈W;i′∈W;j∈O

    uij≤Mp′ijx′ij, ?i∈W;j∈O (15)

    (16)

    (17)

    (18)

    計(jì)算產(chǎn)品的批次成本:

    zi≥∑i′∈W∑j∈Oθii′j+∑r∈Rqirvr+ci+bixi

    M(1-x′i),?i∈W

    zi≤∑i′∈W∑j∈Oθii′j+∑r∈Rqirvr+ci+bixi+

    M(1-x′i),?i∈W

    (19)

    (20)

    計(jì)算產(chǎn)品的單位批產(chǎn)成本:

    pijuijijzi,?i∈W;j∈O

    zi≤Mx′i,?i∈W(21)

    (22)

    定義變量的值域:

    xi為非負(fù)整數(shù),zi≥0;uij≥0

    πii′j,x′i∈{0,1};θii′j≥0,?i,i′∈W;j∈O(23)

    上述模型中,目標(biāo)函數(shù)為滿足外部需求的產(chǎn)品總成本最小化。約束式(9)和式(10)確定了變量xi和x′i之間的關(guān)系;約束式(11)和式(12)確保了生產(chǎn)制程之間的產(chǎn)品輸出量滿足需求量,以及同時(shí)滿足外部需求量,其中式(12)限定了產(chǎn)量批次上限;約束式(12)和式(13)確保了產(chǎn)品成本卷積路徑的合理性,其中式(13)確保了制程i能夠領(lǐng)用制程i′產(chǎn)出的產(chǎn)品j的條件是:s′ij=1且p′ij=1且xi+x′i=2,式(14)確保了產(chǎn)品領(lǐng)用路徑關(guān)系必須得到滿足;約束式(14)~式(18)基于成本卷積路徑變量πii′j確定了成本卷積的數(shù)量關(guān)系θii′j,其中式(14)設(shè)定了若不存在卷積路徑(πii′j=0),則成本卷積數(shù)量θii必須為0,式(15)和式(16)用大數(shù)條件約束法設(shè)定了當(dāng)πii′j=0時(shí),θii必須取值sijuij,即領(lǐng)用的產(chǎn)品數(shù)量乘以產(chǎn)品的單位成本,式(17)設(shè)定了當(dāng)制程i不產(chǎn)出產(chǎn)品j(即p′ij=0)或生產(chǎn)制程i未啟用(即x′i=0)時(shí),制程i產(chǎn)出的產(chǎn)品j的單位成本為0(即uij=0);約束式(19)和式(20)設(shè)定了制程i的單批次成本zi的計(jì)算方法,即“所領(lǐng)用的其他制程產(chǎn)出產(chǎn)品的成本+本制程的生產(chǎn)資源消耗成本+所產(chǎn)生的生產(chǎn)設(shè)備維修成本ci+固定費(fèi)用的分?jǐn)偝杀綛i/xi”,且以大數(shù)條件約束法確保了該約束式僅當(dāng)x′i=1時(shí)有效;約束式(21)和式(22)表示當(dāng)制程i產(chǎn)出多種產(chǎn)品時(shí),產(chǎn)品之間的成本分?jǐn)偙壤凑战o定占比βij分?jǐn)偅剑?2)表示當(dāng)制程i未啟用時(shí),zi為0;約束式(23)定義了變量的值域。

    2.3 模型求解

    上述模型中,約束式(19)和式(20)中存在非線性成分,即Bi/xi。下面利用割線法將該非線性成分進(jìn)行線性化近似化。對(duì)一般性的曲線方程 y=1/x,割線法線性化是指以一組割線y=Klx+Bl(l=1,2,…,η)將曲線 y=1/x代替為線性化方法:

    y≥Klx+Bl,?l=1,2,…,η(24)

    式中:Kl和Bl分別為第l條割線的斜率和截距;η為割線組的數(shù)量,如圖3(a)所示。

    對(duì)于其中的任意一條割線(如第l條割線),令其誤差率定義為ε(x)=(y″-y′)/y′,如圖3(b)所示。可以推導(dǎo)得到,為保證每段割線的最大誤差率不超過(guò)給定值ε,所需要的最少割線數(shù)量為

    η=ln xmax-ln xminln μ(25)

    式中:[xmin, xmax]是x的取值范圍;μ=1+2ε+2ε+ε2,ε是最大誤差率閾值。同樣,可以得到割線組的直線斜率和截距的公式如下:

    Kl=-1μ2l-1·x2min

    Bl=μ+1μl·xmin,?l=1,2,…,η(26)

    按給定參數(shù)要求預(yù)先計(jì)算割線組斜率和截距的算法偽代碼如算法1所示。

    新的線性優(yōu)化模型的完整表達(dá)式由式(8)~式(18)、式(21)~式(23)、式(28)~式(31)組成。該線性優(yōu)化模型可由通用的商業(yè)求解器(如Cplex)直接求解。

    3 模型計(jì)算實(shí)驗(yàn)

    本節(jié)設(shè)計(jì)小中大規(guī)模算例,采用AMPL(a mathematical programming language, AMPL)語(yǔ)言對(duì)模型進(jìn)行了代碼編寫,在Linux系統(tǒng)的AMPL/Cplex環(huán)境,對(duì)第2節(jié)建立的數(shù)學(xué)模型開(kāi)展計(jì)算實(shí)驗(yàn)。

    3.1 小規(guī)模實(shí)證算例

    某企業(yè)研發(fā)出了某產(chǎn)品的新型材料生產(chǎn)工藝,并進(jìn)入了新生產(chǎn)系統(tǒng)優(yōu)化設(shè)計(jì)階段。如圖4所示,該產(chǎn)品需要經(jīng)歷4個(gè)主要生產(chǎn)階段,形成3種中間產(chǎn)品(產(chǎn)品1、產(chǎn)品2和產(chǎn)品3)以及最終產(chǎn)品——產(chǎn)品5。其中,第1階段產(chǎn)出中間產(chǎn)品“產(chǎn)品1”,第2階段領(lǐng)用“產(chǎn)品1”再加工后產(chǎn)出中間產(chǎn)品“產(chǎn)品2”,第3階段領(lǐng)用“產(chǎn)品2”再加工輸出中間產(chǎn)品“產(chǎn)品3”,第4階段領(lǐng)用“產(chǎn)品3”再加工輸出中間產(chǎn)品“產(chǎn)品4”。為了優(yōu)化生產(chǎn)成本,為第1階段設(shè)計(jì)了2種生產(chǎn)方案,即生產(chǎn)制程1和生產(chǎn)制程2,為第3階段也設(shè)計(jì)了2種生產(chǎn)方案,即生產(chǎn)制程4和生產(chǎn)制程5。其中,生產(chǎn)制程5除了產(chǎn)出“產(chǎn)品3”外,還產(chǎn)出“產(chǎn)品4”(為副產(chǎn)品)。同一種產(chǎn)品的不同生產(chǎn)方案對(duì)應(yīng)著不同的投資成本、生產(chǎn)效率/能力、資源消耗和設(shè)備使用情況等。

    圖4中的小規(guī)模算例不失一般性地展示了在生產(chǎn)系統(tǒng)設(shè)計(jì)階段存在多種候選生產(chǎn)方案的情況下,如何根據(jù)生產(chǎn)方案的特點(diǎn)進(jìn)行批產(chǎn)成本卷積路徑優(yōu)化,在滿足產(chǎn)品交付要求的前提下實(shí)現(xiàn)批產(chǎn)成本優(yōu)化的一般性問(wèn)題。

    將第2.2節(jié)中的優(yōu)化模型用AMPL代碼實(shí)現(xiàn),建立模型文件ProdCost.mod,如算法2所示。

    算法 2 定義模型ProdCost.mod輸入: W,O,R,Nkb,p,p1,beta,q,s,s1,v,d,c,e,K,B,b,M

    輸出: Total_Cost

    #文件名ProdCost.mod

    set W;

    #生產(chǎn)制程集合

    set O;

    #產(chǎn)品(或服務(wù))集合

    set R;

    #生產(chǎn)資源集合

    param p{W,O};

    #制程i是否產(chǎn)出產(chǎn)品j的批次產(chǎn)量

    param p1{W,O};

    #制程i是否產(chǎn)出產(chǎn)品j

    param beta{W,O};

    #制程i輸出產(chǎn)品的成本占比

    param q{W,R};

    #制程i對(duì)資源r的每批次消耗量

    param s{W,O};

    #制程i對(duì)產(chǎn)品j的批次用量

    param s1{W,O};

    #制程i對(duì)產(chǎn)品j是否具有領(lǐng)用關(guān)系

    param v{R};

    #生產(chǎn)資源r的成本單價(jià)param d{O};

    #產(chǎn)品j的外部需求量

    param c{W};

    #制程i的生產(chǎn)資源費(fèi)用

    param e{W};

    #制程i的維修費(fèi)用

    set Nkb;

    #線性化的割線集合

    param K{Nkb};

    #割線的斜率

    param B{Nkb};

    #割線的截距

    param b{W};

    #生產(chǎn)制程i的固定費(fèi)用

    param M:=99999;

    #一個(gè)大數(shù)

    var x1{W} binary;

    #0/1變量

    var x{W}gt;=0 integer; #非負(fù)整數(shù)變量

    var z{W}gt;=0;

    #非負(fù)變量

    var z1{W}gt;=0;

    #非負(fù)變量

    var u{W,O}gt;=0;

    #非負(fù)變量

    var pi{W,W,O} binary;

    #0/1變量

    var cita{W,W,O}gt;=0;

    #非負(fù)變量

    minimize Total_Cost: sum{i in W, j in O}(d[j]*u[i,j]);

    subject to Con1a{i in W}: x1[i]lt;=x[i];

    subject to Con1b{i in W}: M*x1[i]gt;=x[i];

    subject to Con2a{j in O}:

    sum{i in W}x[i]*p[i,j]gt;=sum{i in W}x[i]*s[i,j]+

    d[j];

    subject to Con2b{j in O}:

    sum{i in W}(x[i]-1)*p[i,j]lt;=sum{i in W}x[i]*s[i,j]+d[j];

    subject to Con3a{i in W, i1 in W, j in O}:

    pi[i,i1,j]lt;=0.5*s1[i,j]*p1[i1,j]*(x1[i1]+x1[i]);

    subject to Con3b{i in W, j in O}:

    sum{i1 in W}pi[i,i1,j]gt;=s1[i,j]*x1[i];

    subject to Con4a{i in W, i1 in W, j in O}:

    cita[i,i1,j]lt;=M*pi[i,i1,j];

    subject to Con4b{i in W, i1 in W, j in O: ilt;gt;i1}:

    cita[i,i1,j]gt;=s[i,j]*u[i1,j]-M*(1-pi[i,i1,j]);

    subject to Con4c{i in W, i1 in W, j in O: ilt;gt;i1}:

    cita[i,i1,j]lt;=s[i,j]*u[i1,j]+M*(1-pi[i,i1,j]);

    subject to Con4d{i in W, j in O}:

    u[i,j]lt;=M*x1[i]*p1[i,j];

    subject to Con5a{i in W}:

    z[i]gt;=sum{i1 in W, j in O}cita[i,i1,j]+c[i]+e[i]+z1[i]-M*(1-x1[i]);

    subject to Con5b{i in W}:

    z[i]lt;=sum{i1 in W, j in O}cita[i,i1,j]+c[i]+e[i]+z1[i]+M*(1-x1[i]);

    subject to Con5d{i in W}:

    z[i]lt;=M*x1[i];

    subject to Con6{i in W, j in O}:

    p[i,j]*u[i,j]=beta[i,j]*z[i];

    subject to Con7a{i in W, k in Nkb}: z1[i]gt;=b[i]*K[k]*x[i]+b[i]*B[k]-M*(1-x1[i]);

    subject to Con7b{i in W}:

    z1[i]lt;=M*x1[i];

    令算法2輸入中各生產(chǎn)制程的批產(chǎn)成本數(shù)據(jù)如表1所示,生產(chǎn)制程與產(chǎn)品之間的領(lǐng)用數(shù)量、產(chǎn)出數(shù)量和成本占比如表2所示,并令產(chǎn)品3和產(chǎn)品5的外部需求量分別為110和120,其他產(chǎn)品的外部需求量為0。這樣,準(zhǔn)備AMPL優(yōu)化模型所需數(shù)據(jù)文件Instance4x2.dat,內(nèi)容如算法3所示。

    編寫AMPL執(zhí)行腳本程序,調(diào)研Cplex求解算例及輸出,如表3所示。通過(guò)觀測(cè)數(shù)據(jù)可知,生產(chǎn)“產(chǎn)品1”啟用了生產(chǎn)制程2,而生產(chǎn)“產(chǎn)品3”啟用了生產(chǎn)制程4,獲得最低交付成本(目標(biāo)函數(shù))為806.525。對(duì)生產(chǎn)制程基本數(shù)據(jù)進(jìn)行調(diào)整后重新計(jì)算,可以獲得與之相應(yīng)的最優(yōu)結(jié)果,表明模型達(dá)到了設(shè)計(jì)的預(yù)期效果。

    3.2 大規(guī)模求解實(shí)驗(yàn)

    下面通過(guò)隨機(jī)構(gòu)造中、大規(guī)模算例驗(yàn)證所提出模型的求解效率。構(gòu)造較復(fù)雜的生產(chǎn)系統(tǒng),令產(chǎn)品生產(chǎn)過(guò)程包含n個(gè)生產(chǎn)階段,對(duì)應(yīng)n個(gè)中間產(chǎn)品,每個(gè)中間產(chǎn)品設(shè)計(jì)有m個(gè)不同的候選生產(chǎn)制程方案,問(wèn)題規(guī)模為n與m的組合,表示為 n×m。設(shè)置計(jì)算時(shí)間上限為3 600 s。表4給出了計(jì)算結(jié)果,其中黑體數(shù)據(jù)為最優(yōu)值,R.Gap表示計(jì)算終止時(shí),當(dāng)前最好解(上界)與最大松弛解(下界)之間的相對(duì)差距值。可見(jiàn),問(wèn)題組10×2、10×3、10×5、20×2、20×3和30×2的問(wèn)題都獲得了最優(yōu)解。對(duì)于問(wèn)題組20×5、30×3和30×5,本次計(jì)算在給定的時(shí)間(3 600 s)內(nèi)未能獲得最優(yōu)解,僅獲得了當(dāng)前的最好上界值為問(wèn)題的近似優(yōu)化解。這表明增大問(wèn)題規(guī)模參數(shù)n和m都對(duì)求解效率有顯著的正相關(guān)影響。計(jì)算結(jié)果表明,所提出的優(yōu)化模型具有較高的求解效率,可在可接受時(shí)間內(nèi)求解具有10至20個(gè)生產(chǎn)制程的中規(guī)模產(chǎn)品批產(chǎn)成本卷積路徑優(yōu)化問(wèn)題,且獲得了最優(yōu)解(因?yàn)镽.Gap值為0)。對(duì)于具有30個(gè)生產(chǎn)制程的大規(guī)模優(yōu)化問(wèn)題,則需要較長(zhǎng)求解計(jì)算時(shí)間。

    4 結(jié)束語(yǔ)

    本文分析了產(chǎn)品規(guī)模生產(chǎn)階段的主要重復(fù)性生產(chǎn)成本和非重復(fù)成本,提出了產(chǎn)品批產(chǎn)成本卷積路徑優(yōu)化的基本問(wèn)題,給出了該問(wèn)題的一種混合整數(shù)規(guī)劃線性模型。該模型可支持復(fù)雜產(chǎn)品的工藝路線優(yōu)化設(shè)計(jì),在多種復(fù)雜的、可行的工藝路線方案中,選擇出最優(yōu)的批產(chǎn)成本卷積路徑,對(duì)于降低產(chǎn)品批產(chǎn)成本有決策支撐作用。所提出的優(yōu)化模型可應(yīng)用于商業(yè)求解軟件直接求解,對(duì)于生產(chǎn)制程規(guī)模在10~20之間的中規(guī)模算例,可以在較短的計(jì)算時(shí)間(10 min)內(nèi)獲得最優(yōu)解。對(duì)于超過(guò)30個(gè)生產(chǎn)制程的大規(guī)模問(wèn)題,可在較長(zhǎng)計(jì)算時(shí)間(約1 h)內(nèi)獲得最優(yōu)解或接近最優(yōu)解。該問(wèn)題的未來(lái)研究方向有兩點(diǎn):① 對(duì)于更大規(guī)模算例的快速求解問(wèn)題,可設(shè)計(jì)相應(yīng)的求解效率更高的啟發(fā)式算法,以獲取給定時(shí)間內(nèi)的可接受次優(yōu)解;② 結(jié)合不同行業(yè)的工業(yè)應(yīng)用,考慮納入更多的實(shí)際要素,如考慮需求的不確定性、成本時(shí)間屬性,以及考慮產(chǎn)品生產(chǎn)中的學(xué)習(xí)效應(yīng)等。

    參考文獻(xiàn)

    [1]DELLAERT N, JEUNET J, JONARD N. A genetic algorithm to solve the general multi-level lot-sizing problem with time-varying costs[J]. International Journal of Production Economics, 2000, 68(3): 241-257.

    [2]YOU M, XIAO Y Y, ZHANG S Y, et al. Modeling the capacitated multi-level lot-sizing problem under time-varying environments and a fix-and-optimize solution approach[J]. Entropy, 2019, 21(4): 377.

    [3]XIAO Y Y, ZHANG R Q, ZHAO Q H, et al. A variable neighborhood search with an effective local search for uncapacitated multilevel lot-sizing problems[J]. European Journal of Operational Research, 2014, 235(1): 102-114.

    [4]LOLLI F, MESSORI M, GAMBERINI R, et al. Modelling production cost with the effects of learning and forgetting[J]. IFAC-PapersOnLine, 2016, 49(12): 503-508.

    [5]JABER Y M, PELTOKORPI J. The effects of learning in production and group size on the lot-sizing problem[J]. Applied Mathematical Modelling, 2020, 81: 419-427.

    [6]郭映彤. 研制與批產(chǎn)訂單混合調(diào)度建模及優(yōu)化方法研究[D]. 上海: 上海交通大學(xué), 2016.

    GUO Y T. Modeling and optimizing methods research on hybrid scheduling of Ramp;D and batch production[D]. Shanghai: Shanghai Jiao Tong University, 2016.

    [7]于雪磊, 瞿華. 基于“龍芯”處理器軍用裝備產(chǎn)品平臺(tái)化、產(chǎn)品化研究及實(shí)踐[J]. 中國(guó)航天, 2022(S1): 62-66.

    YU X L, QU H. Research and practice of platformisation and productisation of military equipment products based on longchip processor[J]. Aerospace China, 2022(S1): 62-66.

    [8]張變亞. 多品種多階段生產(chǎn)系統(tǒng)問(wèn)題建模優(yōu)化分析[J]. 蘭州石化職業(yè)技術(shù)學(xué)院學(xué)報(bào), 2021, 21(2): 18-23.

    ZHANG B Y. Modeling optimization analysis of multiple varieties and stages production system problems[J]. Journal of Lanzhou Petrochemical Polytechnic, 2021, 21(2): 18-23.

    [9]陳洪根, 閆鑫, 張艷, 等. 生產(chǎn)系統(tǒng)維修管理-統(tǒng)計(jì)過(guò)程控制-經(jīng)濟(jì)生產(chǎn)批量集成優(yōu)化研究綜述[J]. 現(xiàn)代制造工程, 2020(10): 148-155, 12.

    CHEN H G, YAN X, ZHANG Y, et al.Review on the integrated optimization of maintenance management amp; statistical process control amp; economic production quantity on production system[J]. Modern Manufacturing Engineering, 2020(10): 148-155, 12.

    [10]王得水, 黃漫玲. 淺析復(fù)雜電子信息系統(tǒng)裝備產(chǎn)品批產(chǎn)生產(chǎn)線建設(shè)[J]. 裝備制造技術(shù), 2019(11): 209-212.

    WANG D S, HUANG M L. Brief analysis on construction of batch production line of complex electronic information system equipment product[J]. Equipment Manufacturing Technology, 2019(11): 209-212.

    [11]王中婧. 某軍工批產(chǎn)MAO制造項(xiàng)目的進(jìn)度管理研究[D]. 成都: 電子科技大學(xué), 2018.

    WANG Z J. Study on progress management in a military batch production of the MAO manufacturing project[D]. Chengdu: University of Electronic Science and Technology of China, 2018.

    [12]胡楠楠, 南方. 低成本批產(chǎn)商業(yè)通信衛(wèi)星用元器件保障模式探索[J]. 數(shù)字通信世界, 2018(8): 247-259.

    HU N N, NAN F. Exploration of component assurance models for low-cost wholesale production of commercial communications satellites[J]. Digital Communication World, 2018(8): 247-259.

    [13]YU Y, TANG J F, SUN W, et al. Reducing worker(s) by converting assembly line into a pure cell system[J]. International Journal of Production Economics, 2013, 145(2): 799-806.

    [14]KAKU I. Is SERU a sustainable manufacturing system?[J]. Procedia Manufacturing, 2017, 8: 723-730.

    [15]SUER G, ULUTAS B, KAKU I, et al. Considering product life cycle stages and worker skill level in SERU production systems[J]. Procedia Manufacturing, 2019, 39: 1097-1103.

    [16]XIE Y, ZHOU S H, XIAO Y Y, et al. A β-accurate linearization method of Euclidean distance for the facility layout problem with heterogeneous distance metrics[J]. European Journal of Operational Research. 2018, 265(1): 26-38.

    [17]XIAO Y Y, ZHANG Y, KULTUREL-KONAK S, et al. The aperiodic facility layout problem with time-varying demands and an optimal master-slave solution approach[J]. International Journal of Production Research, 2021, 59(17): 5216-5235.

    [18]YOU M, XIAO Y Y, ZHANG S Y, et al. Optimal mathematical programming for the warehouse location problem with Euclidean distance linearization[J]. Computers amp; Industrial Engineering, 2019, 136(C): 70-79.

    [19]XIAO Y Y, YUAN Y Y, ZHANG R Q, et al. Non-permutation flow shop scheduling with order acceptance and weighted tardiness[J]. Applied Mathematics and Computation, 2015, 270(1): 312-333.

    [20]ZHANG R Q, ZHANG L K, XIAO Y Y, et al. The activity-based aggregate production planning with capacity expansion in manufacturing systems[J]. Computers amp; Industrial Engineering, 2012, 62(2): 491-503.

    [21]周宏明, 高順, 張祥雷, 等. 帶緩沖串行生產(chǎn)系統(tǒng)預(yù)防性維護(hù)建模及運(yùn)行參數(shù)優(yōu)化研究[J]. 運(yùn)籌與管理, 2022, 31(1): 22-29.

    ZHOU H M, GAO S, ZHANG X L, et al.Preventive maintenance modeling and operation parameter optimization for series production systems with intermediate buffers[J]. Operations Research and Management Science, 2022, 31(1): 22-29.

    [22]陸志強(qiáng), 張之磊. 串聯(lián)生產(chǎn)系統(tǒng)維護(hù)在線決策與緩沖分配聯(lián)合優(yōu)化[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2021, 49(3): 431-439.

    LU Z Q, ZHANG Z L. Joint optimization of on-line decisionmaking for maintenance and buffer allocation for serial production system[J]. Journal of Tongji University (Natural Science), 2021, 49(3): 431-439.

    [23]方鵬, 李芳, 劉凡, 等. 考慮產(chǎn)品質(zhì)量控制的生產(chǎn)系統(tǒng)預(yù)防性維護(hù)策略聯(lián)合優(yōu)化研究[J]. 上海理工大學(xué)學(xué)報(bào), 2021, 43(5): 497-504.

    FANG P, LI F, LIU F, et al.Optimization of preventive maintenance strategy for production system considering product quality control[J]. Journal of University of Shanghai for Science and Techology, 2021, 43(5): 497-504.

    [24]ADULYASAK Y, CORDEAU J F, JANS R. The production routing problem: a review of formulations and solution algorithms[J]. Computers amp; Operations Research, 2015, 55(1): 141-152.

    [25]CHITSAZ M, CORDEAU J F, JANS R. A branch-and-cut algorithm for an assembly routing problem[J]. European Journal of Operational Research, 2020, 282(3): 896-910.

    [26]ALVAREZ A, MIRANDA P, ROHMER S. Production routing for perishable products[J]. Omega, 2022, 111: 102667.

    [27]FRIFITA S, AFSAR H M, HNAIEN F. An efficient mat-heuristic algorithm for the dynamic disassembly assembly routing problem with returns[J]. European Journal of Industrial Engineering, 2020, 16(5): 584-617.

    [28]MANOUSAKIS E G, KASAPIDIS G A, KIRANOUDIS C T, et al. An infeasible space exploring matheuristic for the production routing problem[J]. European Journal of Operational Research, 2021, 298(2): 478-495.

    [29]LI Y T, CHU F, CHU C B, et al. An efficient three-level heuristic for the large-scaled multi-product production routing problem with outsourcing[J]. European Journal of Operational Research, 2019, 272(3): 914-927.

    [30]RUSSELL-ROBERT A. Mathematical programming heuristics for the production routing problem[J]. International Journal of Production Economics, 2017, 193: 40-49.

    [31]NEVES-MOREIRA F, ALMADA-LOBO B, GUIMARES L, et al. The multi-product inventory-routing problem with pickups and deliveries: mitigating fluctuating demand via rolling horizon heuristics[J]. Transportation Research Part E: Logistics and Transportation Review, 2022, 164: 102791.

    [32]MOONS S, RAMAEKERS K, CARIS A, et al. Integrating production scheduling and vehicle routing decisions at the ope-rational decision level: a review and discussion[J]. Computers amp; Industrial Engineering, 2017, 104(2): 224-245.

    [33]HUANG M, DU B G, GUO J. A hybrid collaborative framework for integrated production scheduling and vehicle routing problem with batch manufacturing and soft time windows[J]. Computers amp; Operations Research, 2023, 159: 106346.

    [34]SHAMAKI B P, ROUX L G. Optimization for sustainable hydrogen production path[J]. Computer Aided Chemical Engineering, 2022, 51: 235-240.

    [35]ERDOGAN A, GULER M G. Optimization and analysis of a hydrogen supply chain in terms of cost, CO2 emissions, and risk: the case of Turkey[J]. International Journal of Hydrogen Energy, 2023, 48(60): 22752-22765.

    [36]RIERA J A, LIMA R M, KNIO O M. A review of hydrogen production and supply chain modeling and optimization[J]. International Journal of Hydrogen Energy, 2023, 48(37): 13731-13755.

    [37]WU J J, ZHANG J, YI W M, et al. Agri-biomass supply chain optimization in north China: model development and application[J]. Energy, 2022, 239(D): 122374.

    [38]DE-JONG S, HOEFNAGELS R, WETTERLUND E, et al. Cost optimization of biofuel production-the impact of scale, integration, transport and supply chain configurations[J]. Applied Energy, 2017, 195(6): 1055-1070.

    [39]DROFENIK J, PAHOR B, KRAVANJA Z, et al. Multi-objective scenario optimization of the food supply chain—Slovenian case study[J]. Computers amp; Chemical Engineering, 2023, 172: 108197.

    作者簡(jiǎn)介

    楊麗穎(1982—),女,高級(jí)工程師,碩士,主要研究方向?yàn)槌杀竟こ獭?/p>

    楊銳意(1999—),男,博士研究生,主要研究方向?yàn)榭煽啃韵到y(tǒng)工程。

    崔新豪(2000—),男,碩士研究生,主要研究方向?yàn)閺?fù)雜系統(tǒng)建模與優(yōu)化。

    張思悅(1996—),女,博士研究生,主要研究方向?yàn)槲锪骶W(wǎng)絡(luò)系統(tǒng)優(yōu)化。

    陳 練(1980—),男,高級(jí)工程師,碩士,主要研究方向?yàn)榻?jīng)濟(jì)性分析。

    肖依永(1973—),男,副教授,博士研究生導(dǎo)師,博士,主要研究方向?yàn)橄到y(tǒng)優(yōu)化理論與經(jīng)濟(jì)性工程。

    猜你喜歡
    優(yōu)化模型路徑優(yōu)化
    基于人工魚群算法優(yōu)化神經(jīng)網(wǎng)絡(luò)在網(wǎng)絡(luò)入侵檢測(cè)中的應(yīng)用研究
    考慮災(zāi)民感知滿意度的突發(fā)事件應(yīng)急救援人員派遣模型
    基于GEM模型的現(xiàn)代化物流產(chǎn)業(yè)集群競(jìng)爭(zhēng)力評(píng)價(jià)和路徑優(yōu)化
    信息時(shí)代數(shù)控銑削的刀具路徑優(yōu)化技術(shù)
    眾籌筑屋優(yōu)化設(shè)計(jì)方案
    經(jīng)濟(jì)發(fā)展方式轉(zhuǎn)變背景下流通體系路徑優(yōu)化策略探討
    山西省異地就醫(yī)直接結(jié)算路徑優(yōu)化研究
    CVRP物流配送路徑優(yōu)化及應(yīng)用研究
    基于優(yōu)化理論的眾籌筑屋模型
    基于意義建構(gòu)視角的企業(yè)預(yù)算管理優(yōu)化路徑探究
    久久久色成人| 亚洲美女视频黄频| 日韩欧美精品v在线| 久久久亚洲精品成人影院| 国产男女超爽视频在线观看| 亚洲经典国产精华液单| 高清视频免费观看一区二区| 国产成人一区二区在线| 欧美+日韩+精品| 十八禁网站网址无遮挡 | 视频中文字幕在线观看| 久久久久久久精品精品| 综合色丁香网| 99热网站在线观看| 久久精品国产自在天天线| 亚洲成人精品中文字幕电影| 九九在线视频观看精品| 免费观看av网站的网址| 精品国产露脸久久av麻豆| 精品人妻偷拍中文字幕| 亚洲真实伦在线观看| 亚洲精品成人久久久久久| 在线观看一区二区三区| 成人国产av品久久久| 国产精品国产三级专区第一集| 中国三级夫妇交换| 91在线精品国自产拍蜜月| 一本一本综合久久| 日日啪夜夜爽| 夜夜看夜夜爽夜夜摸| 久久久久久久久大av| 午夜福利视频精品| 在线免费十八禁| 国产一区有黄有色的免费视频| av福利片在线观看| 在线观看一区二区三区| 一区二区三区四区激情视频| 国产精品久久久久久久久免| 伦理电影大哥的女人| 婷婷色综合大香蕉| 午夜亚洲福利在线播放| 免费看av在线观看网站| 亚洲欧美成人精品一区二区| 亚洲欧美日韩无卡精品| 国产精品精品国产色婷婷| 亚洲人与动物交配视频| 国产探花在线观看一区二区| 久久精品人妻少妇| 日本猛色少妇xxxxx猛交久久| 大香蕉97超碰在线| 老师上课跳d突然被开到最大视频| 五月开心婷婷网| 内地一区二区视频在线| 亚洲国产日韩一区二区| 黄片wwwwww| 久久久久久久精品精品| 极品教师在线视频| 日韩国内少妇激情av| av播播在线观看一区| 我要看日韩黄色一级片| 99热6这里只有精品| 80岁老熟妇乱子伦牲交| 成人毛片60女人毛片免费| www.色视频.com| 久久人人爽人人爽人人片va| 大香蕉97超碰在线| 一级黄片播放器| 99久久精品一区二区三区| 边亲边吃奶的免费视频| 天天躁夜夜躁狠狠久久av| 麻豆国产97在线/欧美| av在线天堂中文字幕| 97精品久久久久久久久久精品| 国产精品精品国产色婷婷| 好男人在线观看高清免费视频| 日韩免费高清中文字幕av| 18禁裸乳无遮挡免费网站照片| 少妇高潮的动态图| 欧美+日韩+精品| 亚洲色图av天堂| 国产探花极品一区二区| 欧美高清成人免费视频www| 亚洲人成网站在线观看播放| 久久久a久久爽久久v久久| 真实男女啪啪啪动态图| 国产 精品1| 免费观看的影片在线观看| 我要看日韩黄色一级片| 精品少妇久久久久久888优播| 高清av免费在线| 精品久久久久久久久亚洲| av免费观看日本| 汤姆久久久久久久影院中文字幕| 噜噜噜噜噜久久久久久91| 国产成人91sexporn| 亚洲国产精品专区欧美| 中国国产av一级| 精品久久久久久久末码| 成人亚洲精品av一区二区| 国产精品国产av在线观看| 日韩欧美 国产精品| 视频区图区小说| 国产亚洲午夜精品一区二区久久 | 插阴视频在线观看视频| 狂野欧美激情性xxxx在线观看| 高清av免费在线| 国产成人精品婷婷| 少妇裸体淫交视频免费看高清| 欧美激情在线99| 色播亚洲综合网| 男插女下体视频免费在线播放| 国产精品国产三级国产专区5o| 91精品国产九色| 中文字幕亚洲精品专区| av黄色大香蕉| 黄色配什么色好看| 日产精品乱码卡一卡2卡三| 亚洲精品乱久久久久久| 男人爽女人下面视频在线观看| 国产日韩欧美亚洲二区| 国产成人精品福利久久| 久久久久久久久大av| 亚洲国产精品成人综合色| 久久久久性生活片| 特级一级黄色大片| 午夜福利在线在线| 国产一区亚洲一区在线观看| 成年免费大片在线观看| 可以在线观看毛片的网站| 国产午夜福利久久久久久| 91久久精品电影网| 亚洲欧美一区二区三区黑人 | 日本免费在线观看一区| 亚洲欧美精品自产自拍| 国产黄片美女视频| 国产精品久久久久久精品古装| 在线观看人妻少妇| 中文在线观看免费www的网站| 在线观看一区二区三区激情| 麻豆精品久久久久久蜜桃| 美女cb高潮喷水在线观看| 女人十人毛片免费观看3o分钟| 性插视频无遮挡在线免费观看| 18禁裸乳无遮挡免费网站照片| 亚洲欧美清纯卡通| av在线app专区| 精品熟女少妇av免费看| 日本wwww免费看| 久久97久久精品| 久久精品国产鲁丝片午夜精品| 国产探花极品一区二区| 免费电影在线观看免费观看| 色播亚洲综合网| 午夜爱爱视频在线播放| 欧美丝袜亚洲另类| 国产淫语在线视频| 久久久成人免费电影| av国产精品久久久久影院| 国产乱来视频区| 亚洲精品乱码久久久v下载方式| 国产亚洲av嫩草精品影院| 熟妇人妻不卡中文字幕| 亚洲三级黄色毛片| 国产乱人视频| 欧美精品国产亚洲| 亚洲av一区综合| 熟妇人妻不卡中文字幕| 麻豆成人av视频| 哪个播放器可以免费观看大片| 亚洲精品亚洲一区二区| 亚洲欧美一区二区三区黑人 | 国产片特级美女逼逼视频| 欧美区成人在线视频| 久久亚洲国产成人精品v| 国产成人一区二区在线| 中文字幕久久专区| 可以在线观看毛片的网站| 可以在线观看毛片的网站| 网址你懂的国产日韩在线| 尾随美女入室| 欧美一级a爱片免费观看看| 亚洲人成网站在线观看播放| 欧美97在线视频| av在线老鸭窝| 中国美白少妇内射xxxbb| 波野结衣二区三区在线| 菩萨蛮人人尽说江南好唐韦庄| 色播亚洲综合网| 亚洲在线观看片| 制服丝袜香蕉在线| 美女国产视频在线观看| 亚洲欧美一区二区三区国产| 美女高潮的动态| 日韩国内少妇激情av| 亚洲国产精品专区欧美| 久久久色成人| 国产日韩欧美在线精品| 久久99蜜桃精品久久| 国产精品精品国产色婷婷| 国产在线一区二区三区精| 美女脱内裤让男人舔精品视频| 亚洲av成人精品一二三区| 午夜爱爱视频在线播放| 国产综合懂色| 国产 一区精品| 大片电影免费在线观看免费| 国产精品蜜桃在线观看| 成人亚洲精品一区在线观看 | 精华霜和精华液先用哪个| av国产精品久久久久影院| 69av精品久久久久久| 欧美潮喷喷水| 国产精品一区二区性色av| 国产精品人妻久久久影院| 一级毛片电影观看| 少妇猛男粗大的猛烈进出视频 | 你懂的网址亚洲精品在线观看| 免费看a级黄色片| 中文精品一卡2卡3卡4更新| 国模一区二区三区四区视频| 精品人妻熟女av久视频| 高清毛片免费看| 美女被艹到高潮喷水动态| 国内精品美女久久久久久| 国产真实伦视频高清在线观看| 最近最新中文字幕免费大全7| 99九九线精品视频在线观看视频| 一区二区三区免费毛片| 欧美极品一区二区三区四区| 80岁老熟妇乱子伦牲交| 三级国产精品欧美在线观看| 精品国产三级普通话版| 久久久久久久午夜电影| 国模一区二区三区四区视频| 嘟嘟电影网在线观看| av黄色大香蕉| 亚洲av免费高清在线观看| 国产日韩欧美亚洲二区| 国产成人午夜福利电影在线观看| 中文天堂在线官网| 国产永久视频网站| 日日撸夜夜添| 色视频在线一区二区三区| 欧美极品一区二区三区四区| 特级一级黄色大片| 日韩制服骚丝袜av| 毛片一级片免费看久久久久| 国产精品99久久99久久久不卡 | 22中文网久久字幕| 亚洲精品色激情综合| av卡一久久| 97精品久久久久久久久久精品| 亚洲伊人久久精品综合| 久久人人爽av亚洲精品天堂 | 777米奇影视久久| av卡一久久| 久热久热在线精品观看| 国产伦在线观看视频一区| 91久久精品电影网| 久久久国产一区二区| 亚洲va在线va天堂va国产| 国产一区有黄有色的免费视频| av天堂中文字幕网| 日韩制服骚丝袜av| 男人爽女人下面视频在线观看| 国产老妇伦熟女老妇高清| 国产欧美日韩一区二区三区在线 | 内射极品少妇av片p| 国产老妇女一区| 亚洲av.av天堂| 午夜老司机福利剧场| 精品久久久噜噜| 精品熟女少妇av免费看| 大码成人一级视频| 老司机影院成人| 欧美高清性xxxxhd video| www.av在线官网国产| av在线观看视频网站免费| 观看美女的网站| 禁无遮挡网站| 男女那种视频在线观看| 亚洲电影在线观看av| 久久久久久久久大av| 一本久久精品| 少妇的逼水好多| 精品国产一区二区三区久久久樱花 | 蜜桃亚洲精品一区二区三区| 日韩国内少妇激情av| 久久久午夜欧美精品| 久久久精品94久久精品| 欧美成人午夜免费资源| 日日摸夜夜添夜夜添av毛片| 亚洲自拍偷在线| 色播亚洲综合网| 别揉我奶头 嗯啊视频| 久久久久久久久久成人| 边亲边吃奶的免费视频| 中文精品一卡2卡3卡4更新| 麻豆成人午夜福利视频| 新久久久久国产一级毛片| 成人亚洲精品一区在线观看 | 久久精品国产a三级三级三级| 三级国产精品欧美在线观看| 午夜激情久久久久久久| 老司机影院成人| 小蜜桃在线观看免费完整版高清| 午夜日本视频在线| 日韩强制内射视频| 欧美 日韩 精品 国产| 亚洲最大成人中文| 22中文网久久字幕| 成人国产麻豆网| 国产精品久久久久久精品电影小说 | 我的老师免费观看完整版| 久久久国产一区二区| 一个人观看的视频www高清免费观看| 99热这里只有是精品在线观看| 另类亚洲欧美激情| 自拍偷自拍亚洲精品老妇| 国产亚洲午夜精品一区二区久久 | 99热这里只有是精品在线观看| 亚洲欧美日韩东京热| av国产免费在线观看| 超碰av人人做人人爽久久| 在线免费十八禁| 久久久欧美国产精品| 我的老师免费观看完整版| av在线天堂中文字幕| 国产色爽女视频免费观看| 肉色欧美久久久久久久蜜桃 | 少妇的逼水好多| 又大又黄又爽视频免费| 国产精品国产三级国产专区5o| 你懂的网址亚洲精品在线观看| 欧美一区二区亚洲| 另类亚洲欧美激情| 七月丁香在线播放| 欧美3d第一页| 中文乱码字字幕精品一区二区三区| 人妻制服诱惑在线中文字幕| 国产男女内射视频| 成人亚洲精品av一区二区| 男女啪啪激烈高潮av片| 国产成人福利小说| 久久久久国产精品人妻一区二区| 伊人久久精品亚洲午夜| 欧美日韩视频精品一区| 男女国产视频网站| 黄色视频在线播放观看不卡| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 中国三级夫妇交换| 国产亚洲av嫩草精品影院| 小蜜桃在线观看免费完整版高清| 欧美精品人与动牲交sv欧美| 国产成年人精品一区二区| 毛片一级片免费看久久久久| 精品99又大又爽又粗少妇毛片| 欧美xxⅹ黑人| 亚洲av免费在线观看| 五月天丁香电影| 大又大粗又爽又黄少妇毛片口| 2022亚洲国产成人精品| 深爱激情五月婷婷| 伊人久久国产一区二区| 99热网站在线观看| 婷婷色av中文字幕| 日本黄大片高清| 婷婷色综合www| 99re6热这里在线精品视频| 超碰av人人做人人爽久久| 精品一区二区三卡| 久久久国产一区二区| 肉色欧美久久久久久久蜜桃 | 亚洲第一区二区三区不卡| 91精品伊人久久大香线蕉| 欧美日本视频| 久久久久久久午夜电影| 又黄又爽又刺激的免费视频.| 亚洲av男天堂| 亚洲精品一区蜜桃| 国产成人91sexporn| 韩国高清视频一区二区三区| 91久久精品电影网| 亚洲最大成人手机在线| 三级国产精品欧美在线观看| 亚洲va在线va天堂va国产| 国产老妇女一区| 国产成人免费无遮挡视频| 亚洲人成网站高清观看| 欧美日韩精品成人综合77777| 狂野欧美白嫩少妇大欣赏| 日日啪夜夜撸| 在线免费十八禁| 国产精品.久久久| 亚洲精品成人久久久久久| 少妇人妻精品综合一区二区| 夫妻性生交免费视频一级片| 搡老乐熟女国产| 超碰97精品在线观看| 亚洲最大成人中文| 亚洲国产精品专区欧美| 精品人妻熟女av久视频| 成人欧美大片| 黄片无遮挡物在线观看| 日韩免费高清中文字幕av| 午夜免费观看性视频| 久久久久久久久久人人人人人人| 麻豆成人av视频| 日本与韩国留学比较| 麻豆成人午夜福利视频| 国产在线男女| 欧美亚洲 丝袜 人妻 在线| 男插女下体视频免费在线播放| 一级毛片aaaaaa免费看小| av国产免费在线观看| 国产av国产精品国产| 亚洲图色成人| 97人妻精品一区二区三区麻豆| 精品久久久久久久末码| 亚洲精品久久久久久婷婷小说| 国产一区二区亚洲精品在线观看| 人妻系列 视频| 人妻 亚洲 视频| 国产永久视频网站| 欧美另类一区| 日韩免费高清中文字幕av| 国产成人a∨麻豆精品| 亚洲色图综合在线观看| 超碰97精品在线观看| av天堂中文字幕网| 国产成年人精品一区二区| 人妻少妇偷人精品九色| 欧美成人a在线观看| 精品久久久久久久末码| 国产91av在线免费观看| av福利片在线观看| 日本猛色少妇xxxxx猛交久久| 日本爱情动作片www.在线观看| 国产久久久一区二区三区| 国产爱豆传媒在线观看| 中国美白少妇内射xxxbb| 欧美+日韩+精品| 九九久久精品国产亚洲av麻豆| 2021少妇久久久久久久久久久| 欧美日韩亚洲高清精品| 亚洲精品亚洲一区二区| 99九九线精品视频在线观看视频| 别揉我奶头 嗯啊视频| 一区二区三区精品91| 97热精品久久久久久| 日日摸夜夜添夜夜添av毛片| 亚洲精品国产成人久久av| 亚洲精品国产av蜜桃| 日本-黄色视频高清免费观看| 99久国产av精品国产电影| 五月开心婷婷网| 国产男女内射视频| av在线亚洲专区| 久久久久久久午夜电影| 特级一级黄色大片| 少妇裸体淫交视频免费看高清| 国产探花在线观看一区二区| 国产亚洲91精品色在线| 大又大粗又爽又黄少妇毛片口| 久久久精品欧美日韩精品| 亚洲精品国产av蜜桃| 日韩大片免费观看网站| 国产中年淑女户外野战色| 日本午夜av视频| 亚洲精品中文字幕在线视频 | 久久久久久伊人网av| 免费人成在线观看视频色| 青春草亚洲视频在线观看| 免费看光身美女| 99视频精品全部免费 在线| 亚洲av男天堂| 97在线人人人人妻| 黄片无遮挡物在线观看| 涩涩av久久男人的天堂| 精品久久久久久电影网| 欧美97在线视频| 我的女老师完整版在线观看| 秋霞伦理黄片| 2018国产大陆天天弄谢| 人人妻人人爽人人添夜夜欢视频 | 亚洲国产日韩一区二区| 日本爱情动作片www.在线观看| 欧美区成人在线视频| 亚洲婷婷狠狠爱综合网| 纵有疾风起免费观看全集完整版| 黄色视频在线播放观看不卡| 成人黄色视频免费在线看| 超碰97精品在线观看| 国产亚洲5aaaaa淫片| 日本三级黄在线观看| 美女主播在线视频| 亚洲在线观看片| 美女国产视频在线观看| 国产精品国产av在线观看| 亚洲综合色惰| www.av在线官网国产| 亚洲av一区综合| 插逼视频在线观看| 婷婷色综合www| 蜜臀久久99精品久久宅男| 亚洲精品成人av观看孕妇| 99热国产这里只有精品6| 狠狠精品人妻久久久久久综合| 七月丁香在线播放| 中文字幕久久专区| 亚洲av一区综合| 九草在线视频观看| 国产精品精品国产色婷婷| 夫妻性生交免费视频一级片| 久久午夜福利片| 亚洲激情五月婷婷啪啪| 狂野欧美激情性bbbbbb| 少妇人妻 视频| 熟妇人妻不卡中文字幕| 国国产精品蜜臀av免费| 成年女人在线观看亚洲视频 | 男女无遮挡免费网站观看| 欧美一区二区亚洲| 亚洲在久久综合| 亚洲四区av| 欧美xxⅹ黑人| 香蕉精品网在线| 黄片无遮挡物在线观看| 亚洲av不卡在线观看| 日韩精品有码人妻一区| 日本三级黄在线观看| 欧美另类一区| 人人妻人人爽人人添夜夜欢视频 | 欧美激情久久久久久爽电影| 2021少妇久久久久久久久久久| 日韩人妻高清精品专区| 中文字幕亚洲精品专区| 亚洲伊人久久精品综合| 久久久欧美国产精品| 青春草亚洲视频在线观看| 精品国产露脸久久av麻豆| 久久久色成人| 亚洲一区二区三区欧美精品 | 欧美日韩在线观看h| 蜜桃亚洲精品一区二区三区| 一级片'在线观看视频| 欧美潮喷喷水| 性插视频无遮挡在线免费观看| 国产高清不卡午夜福利| 建设人人有责人人尽责人人享有的 | av在线老鸭窝| 97超视频在线观看视频| 免费看光身美女| 麻豆成人午夜福利视频| 久久精品人妻少妇| 国产久久久一区二区三区| 成人国产av品久久久| 欧美日韩国产mv在线观看视频 | 国产一区二区三区av在线| 久久久a久久爽久久v久久| 国产亚洲午夜精品一区二区久久 | 亚洲av男天堂| 天堂网av新在线| 97在线视频观看| 天天一区二区日本电影三级| 亚洲色图综合在线观看| 日本免费在线观看一区| 亚洲欧美中文字幕日韩二区| freevideosex欧美| 极品教师在线视频| 亚洲怡红院男人天堂| 最近最新中文字幕大全电影3| 亚洲国产欧美人成| 亚洲无线观看免费| 一区二区三区四区激情视频| 久久精品人妻少妇| 高清在线视频一区二区三区| 一边亲一边摸免费视频| av福利片在线观看| 在线免费十八禁| 久久久成人免费电影| 大又大粗又爽又黄少妇毛片口| 一个人观看的视频www高清免费观看| 一个人看的www免费观看视频| 狂野欧美激情性bbbbbb| 日本与韩国留学比较| 亚洲天堂国产精品一区在线| 亚州av有码| 伊人久久精品亚洲午夜| 国产精品一区www在线观看| 国产亚洲5aaaaa淫片| 干丝袜人妻中文字幕| av网站免费在线观看视频| 在现免费观看毛片| 国产爱豆传媒在线观看| 大香蕉97超碰在线| 成人免费观看视频高清| 日本一二三区视频观看| 国产成人精品久久久久久| 午夜激情福利司机影院| 亚洲国产高清在线一区二区三| 天天躁日日操中文字幕| 成人亚洲精品一区在线观看 | 国产精品麻豆人妻色哟哟久久| 国国产精品蜜臀av免费| 人妻一区二区av| 在线观看免费高清a一片| 99视频精品全部免费 在线| 少妇裸体淫交视频免费看高清| 成人美女网站在线观看视频| 日韩精品有码人妻一区| 国产精品蜜桃在线观看| 欧美亚洲 丝袜 人妻 在线| 中文字幕制服av|