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

    多重不確定因素影響的高端裝備研制任務(wù)仿真建模

    2018-06-07 08:07:15張西林譚躍進(jìn)楊志偉
    關(guān)鍵詞:成本模型

    張西林, 譚躍進(jìn), 楊志偉

    (1. 國防科技大學(xué)系統(tǒng)工程學(xué)院, 湖南 長沙 410073; 2. 江蘇師范大學(xué)商學(xué)院, 江蘇 徐州 221116)

    0 引 言

    高端裝備是指技術(shù)含量高、涉及學(xué)科多、資金投入大、具有高附加值、服役壽命長,一般需要跨領(lǐng)域、跨行業(yè)、跨區(qū)域的多個(gè)企業(yè)共同協(xié)作才能完成的一類技術(shù)裝備,如航空航天裝備、海洋工程裝備、軌道交通裝備、衛(wèi)星、智能制造裝備等。高端裝備研制任務(wù)通常包含許多子任務(wù),子任務(wù)之間存在串行、并行、重疊和耦合等關(guān)系,而且任務(wù)網(wǎng)絡(luò)復(fù)雜、隨機(jī)性大、影響因素多、返工迭代多。各個(gè)子任務(wù)的成功概率、完成時(shí)間、消耗成本、返工迭代次數(shù)等都存在著一定的不確定性,需要綜合考慮這些因素,以便對(duì)其工期、成本、風(fēng)險(xiǎn)等進(jìn)行更準(zhǔn)確的預(yù)測、評(píng)估和控制。

    計(jì)劃評(píng)審技術(shù)(program evaluation and review technique,PERT)是美國研制北極星導(dǎo)彈時(shí)提出并發(fā)展起來的一種網(wǎng)絡(luò)計(jì)劃技術(shù),主要用于項(xiàng)目的工期優(yōu)化、關(guān)鍵路線識(shí)別、敏感性分析等。文獻(xiàn)[1-2]利用PERT網(wǎng)絡(luò)模型研究了項(xiàng)目管理和裝備研制的工期估算問題。但是PERT是一個(gè)確定性的網(wǎng)絡(luò),而且不允許存在回路,難以對(duì)高端裝備研制任務(wù)的迭代、重疊等進(jìn)行描述。圖示評(píng)審技術(shù)(graphical evaluation and review technique,GERT)方法可以將研制任務(wù)或項(xiàng)目的各個(gè)活動(dòng)的參數(shù)設(shè)為隨機(jī)分布,而且網(wǎng)絡(luò)中允許存在自環(huán)和回路。文獻(xiàn)[3-6]利用GERT網(wǎng)絡(luò)模型從不同角度研究了復(fù)雜裝備研制問題。但是GERT無法描述高端裝備研制任務(wù)中的子任務(wù)重疊和學(xué)習(xí)曲線。

    文獻(xiàn)[7]將設(shè)計(jì)結(jié)構(gòu)矩陣(design structure matrix,DSM)應(yīng)用到復(fù)雜系統(tǒng)分析與管理,DSM在產(chǎn)品的流程優(yōu)化、協(xié)同設(shè)計(jì)、風(fēng)險(xiǎn)分析等領(lǐng)域得到了廣泛的應(yīng)用[8-10]。DSM用方陣表示項(xiàng)目中活動(dòng)之間的串行、并行、耦合和重疊等關(guān)系,能夠有效解決流程分析建模、仿真與優(yōu)化等問題。文獻(xiàn)[9-10]結(jié)合不同實(shí)例分別研究了DSM的應(yīng)用,指出研究DSM需要結(jié)合圖論、網(wǎng)絡(luò)分析、復(fù)雜網(wǎng)絡(luò)等理論。文獻(xiàn)[11]首次構(gòu)建了基于DSM的產(chǎn)品研發(fā)流程架構(gòu)仿真模型,奠定了基于DSM的研發(fā)任務(wù)仿真的基礎(chǔ),該模型可以分析不同流程下的進(jìn)度、風(fēng)險(xiǎn),但是未考慮資源約束。文獻(xiàn)[12-14]通過模型分析了不同項(xiàng)目屬性對(duì)重疊策略的影響,研究了如何確定產(chǎn)品開發(fā)過程中的最優(yōu)重疊策略。文獻(xiàn)[15]基于演化DSM和靈敏度DSM模型衡量重疊活動(dòng)之間的交互強(qiáng)度。文獻(xiàn)[16]在分析上下游活動(dòng)依賴性的基礎(chǔ)上,提出了不同重疊策略下量化產(chǎn)品設(shè)計(jì)演化過程的評(píng)估模型。高端裝備研制過程中不可避免地會(huì)涉及到迭代過程,迭代具有積極的作用,但同時(shí)會(huì)增加研制任務(wù)的工期和成本[17-18]。當(dāng)產(chǎn)品研制過程中允許活動(dòng)返工時(shí),產(chǎn)品研發(fā)工期和方差難以估計(jì),文獻(xiàn)[19]提出了一種研發(fā)工期和方差的估計(jì)和分析方法。文獻(xiàn)[20]研究了結(jié)構(gòu)復(fù)雜度對(duì)成本和工期的影響。文獻(xiàn)[21-22]綜合考慮返工的不確定性、返工風(fēng)險(xiǎn)等因素,采用遺傳算法優(yōu)化DSM。文獻(xiàn)[23]考慮了復(fù)雜產(chǎn)品研發(fā)過程中活動(dòng)的返工迭代、隨機(jī)重疊、學(xué)習(xí)效應(yīng)等多種不確定關(guān)系,建立了其仿真模型。文獻(xiàn)[24]構(gòu)造了一種基于裝配特征對(duì)的裝配建模與仿真集成框架。文獻(xiàn)[25]在DSM分析的基礎(chǔ)上,基于仿真方法評(píng)估產(chǎn)品架構(gòu)對(duì)可演化性的影響。文獻(xiàn)[26]綜合利用DSM和網(wǎng)絡(luò)分析方法研究了項(xiàng)目實(shí)施過程中任務(wù)之間的交互。文獻(xiàn)[27]構(gòu)建了能夠反映產(chǎn)品研發(fā)過程變化的DSM網(wǎng)絡(luò),并提出了基于多層次建模與仿真的產(chǎn)品開發(fā)流程管理方法。文獻(xiàn)[28]開發(fā)了一個(gè)離散事件仿真模型,用來分析產(chǎn)品設(shè)計(jì)過程中返工迭代和變更傳播的影響。文獻(xiàn)[29]在資源受限條件下,考慮節(jié)點(diǎn)之間的多種關(guān)系,實(shí)現(xiàn)了基于DSM的項(xiàng)目網(wǎng)絡(luò)仿真。文獻(xiàn)[30]考慮返工和學(xué)習(xí)效應(yīng)等因素,基于離散事件仿真估計(jì)資源受限研發(fā)項(xiàng)目的工期。

    當(dāng)前文獻(xiàn)在進(jìn)行研制任務(wù)仿真時(shí),主要考慮了其子任務(wù)的順序、并行、重疊、返工迭代等關(guān)系,沒有考慮子任務(wù)的自環(huán)和失敗。在高端裝備研制任務(wù)中,很多子任務(wù)都需要高新技術(shù)突破,完成時(shí)間、能否成功等都存在一定的不確定性,需要進(jìn)一步考慮其子任務(wù)的自環(huán)和失敗。本文在當(dāng)前研究的基礎(chǔ)上,考慮了子任務(wù)的自環(huán)、失敗等因素,對(duì)高端裝備研制任務(wù)進(jìn)行仿真建模,使得仿真模型能夠更真實(shí)地反映研制任務(wù)的實(shí)際過程。

    1 高端裝備研制任務(wù)仿真問題及子任務(wù)關(guān)系描述

    1.1 問題描述

    一般來說,高端裝備研制任務(wù)可參考的歷史數(shù)據(jù)、信息較少,沒有類似的成功案例可供借鑒,其研制過程的不確定性更大。由于受到多種因素的影響,各個(gè)子任務(wù)完工所需的工期和成本會(huì)在一定范圍內(nèi)波動(dòng)。某子任務(wù)完工引發(fā)的后續(xù)子任務(wù)及其執(zhí)行方式是不確定的,例如,某子任務(wù)完工可能會(huì)引發(fā)自身返工或其他若干子任務(wù)的返工,也可能引發(fā)其緊后子任務(wù)的執(zhí)行,而且其后續(xù)子任務(wù)可能是重疊或非重疊執(zhí)行。各個(gè)子任務(wù)會(huì)存在一定的隨機(jī)失敗風(fēng)險(xiǎn)。在整個(gè)研制任務(wù)執(zhí)行過程中,由于資源是有限的,子任務(wù)之間可能會(huì)存在資源沖突。為了實(shí)現(xiàn)整個(gè)研制任務(wù)的工期最短、成本最小、失敗率最低等目標(biāo),需要對(duì)整個(gè)研制任務(wù)過程進(jìn)行分析、優(yōu)化,由于問題的復(fù)雜性高和眾多的不確定性,難以用數(shù)學(xué)模型求解分析。因此,本文用一定的分布或概率描述各個(gè)子任務(wù)工期、成本及子任務(wù)之間關(guān)系的不確定性,采用離散事件仿真建模的方法對(duì)研制任務(wù)進(jìn)行仿真和優(yōu)化。通過仿真,能夠更直觀地表示研制任務(wù)執(zhí)行過程中子任務(wù)的工期和成本的不確定性、子任務(wù)之間的各種關(guān)系以及執(zhí)行過程中的返工迭代、重疊、失敗、學(xué)習(xí)效應(yīng)等,從而能夠有效把握整個(gè)研制任務(wù)的工期、成本等,并可以在此基礎(chǔ)上對(duì)各影響因素進(jìn)行深入分析。

    1.2 DSM、返工、重疊及學(xué)習(xí)曲線的量化表示

    1.2.1 DSM的量化表示

    DSM的元素表示該元素對(duì)應(yīng)的列任務(wù)向?qū)?yīng)的行任務(wù)輸出或者支持信息。給定n個(gè)子任務(wù)Ti(i=1,2,3,…,n)構(gòu)成的研制任務(wù)過程,則DSM中各行各列對(duì)應(yīng)的元素表示子任務(wù)與子任務(wù)之間的關(guān)系,矩陣中的元素Aij表示子任務(wù)Tj輸出信息給子任務(wù)Ti。DSM各元素的取值定義為:DSMij=0表示子任務(wù)i與子任務(wù)j之間不存在緊前緊后關(guān)系;DSMij=1表示子任務(wù)i是子任務(wù)j的緊后任務(wù);DSMij=p,p∈(0,1)表示子任務(wù)i與子任務(wù)j重疊執(zhí)行的概率是p。

    1.2.2 返工概率矩陣和返工影響矩陣

    返工概率矩陣描述了返工迭代的不確定性,以一定的返工概率進(jìn)行表示。其元素RPij表示子任務(wù)j引發(fā)子任務(wù)i發(fā)生返工的概率,其中i,j=1,2,3,…,n。某子任務(wù)每完成一次返工迭代,其存在錯(cuò)誤的概率降低、產(chǎn)生的新信息減少,再次發(fā)生返工的概率也會(huì)隨之降低。一般來說,兩個(gè)子任務(wù)之間的關(guān)系越復(fù)雜,返工子任務(wù)越難以按照返工要求完成,再次返工的概率也就越大。因此,其返工概率降低系數(shù)就越大;反之,其返工概率降低系數(shù)就越小。第k次由子任務(wù)j引發(fā)子任務(wù)i發(fā)生返工的概率可表示為

    (1)

    式中,RPij(0)為子任務(wù)j引發(fā)子任務(wù)i返工的初始返工概率;RPLij為子任務(wù)j引發(fā)子任務(wù)i返工的概率降低系數(shù)。

    返工影響矩陣表示當(dāng)發(fā)生返工時(shí),對(duì)返工子任務(wù)的影響程度。其元素RIij表示子任務(wù)j引發(fā)子任務(wù)i發(fā)生返工時(shí),子任務(wù)i返工的工作量占其原工作量的比例,其中i,j=1,2,3,…,n。

    1.2.3 學(xué)習(xí)曲線

    學(xué)習(xí)曲線表示同一子任務(wù)每重復(fù)執(zhí)行一次,完工所需工作量的減小程度。學(xué)習(xí)曲線的大小與子任務(wù)的復(fù)雜性有關(guān),一般來說,子任務(wù)的復(fù)雜性越高,熟練該子任務(wù)的速度越慢,重復(fù)執(zhí)行時(shí)所需的工作量越大,其學(xué)習(xí)曲線值就越大;反之,其學(xué)習(xí)曲線值就越小。LCi表示子任務(wù)i每次返工所需的工作量占其上次完工所需工作量的比例,其中i,j=1,2,3,…,n。子任務(wù)i第k次返工所需的工期和成本分別可表示為

    (2)

    (3)

    式中,Di、Ci分別為子任務(wù)i的工期、成本。

    由于每次返工需要子任務(wù)之間的信息交流、作業(yè)準(zhǔn)備等活動(dòng),所需的工作量不能無限小,假設(shè)某子任務(wù)返工所需的工作量降低到一定程度后,以后每次返工所需工作量是固定的。

    1.2.4 重疊比例矩陣和重疊影響矩陣

    子任務(wù)重疊是把緊后子任務(wù)提前開始與緊前子任務(wù)進(jìn)行部分并行執(zhí)行,因此子任務(wù)重疊條件判斷時(shí),只有重疊概率以及資源約束兩個(gè)條件都滿足,才能重疊執(zhí)行,而且緊前子任務(wù)開始執(zhí)行之前緊后子任務(wù)不能開始執(zhí)行。

    重疊比例矩陣表示子任務(wù)發(fā)生重疊時(shí)允許的重疊程度。ORij表示當(dāng)子任務(wù)i與子任務(wù)j發(fā)生重疊時(shí),允許的重疊時(shí)間與子任務(wù)i工期的比,其中i,j=1,2,3,…,n。子任務(wù)i與子任務(wù)j的重疊量可表示為

    OAij=min(Di×ORij,Dj)

    (4)

    重疊影響矩陣表示由于重疊執(zhí)行引發(fā)的重疊返工程度。OIij表示當(dāng)子任務(wù)i與子任務(wù)j發(fā)生重疊時(shí),需要的重疊返工量與重疊量的比,其中i,j=1,2,3,…,n。子任務(wù)i與子任務(wù)j的重疊返工量RWij=OAij×OIij,相應(yīng)的子任務(wù)i的實(shí)際工期和實(shí)際成本分別為

    ADi=Di+RWij

    (5)

    (6)

    1.3 子任務(wù)自環(huán)的量化表示

    自環(huán)是子任務(wù)自身的返工。當(dāng)某子任務(wù)完工后,進(jìn)行驗(yàn)收,確定子任務(wù)是否達(dá)到要求,如果達(dá)到要求,繼續(xù)進(jìn)行后續(xù)子任務(wù);如果沒有達(dá)到要求,需要進(jìn)行自身的返工。自環(huán)的發(fā)生具有隨機(jī)性,發(fā)生自環(huán)的概率和返工影響分別在返工概率矩陣和返工影響矩陣的對(duì)角線位置進(jìn)行表示。

    在返工概率矩陣中,對(duì)角線元素RPii表示該元素發(fā)生自環(huán)的概率,其中:0≤RPii<1。RPii取0,表示不會(huì)發(fā)生自環(huán);RPii取(0,1)之間的數(shù)值,表示以一定的概率RPii發(fā)生自環(huán)。在返工影響矩陣中,對(duì)角線元素RIii表示該元素發(fā)生自環(huán)時(shí),所引起的子任務(wù)i的額外工作量的消耗量占該子任務(wù)初始工作量的比例,其中0≤RIii<1。在RPii=0的對(duì)應(yīng)矩陣位置,RIii取0,表示不會(huì)發(fā)生自環(huán);在0

    1.4 子任務(wù)失敗的量化表示

    高端裝備研制任務(wù)的很多子任務(wù)需要高新技術(shù)突破,不能確保子任務(wù)一定能夠成功完成,因此,子任務(wù)會(huì)存在著失敗的可能。子任務(wù)失敗往往會(huì)給整個(gè)研制任務(wù)帶來毀滅性的風(fēng)險(xiǎn),因此,在仿真過程中,需要考慮子任務(wù)失敗。本文假設(shè)如果某個(gè)子任務(wù)失敗,后續(xù)子任務(wù)無法繼續(xù)進(jìn)行,造成整個(gè)研制任務(wù)的失敗,本次仿真以“任務(wù)失敗”的形式結(jié)束。子任務(wù)每執(zhí)行一次,會(huì)取得一定的技術(shù)突破,進(jìn)行該子任務(wù)返工時(shí),其失敗概率降低,因此子任務(wù)的失敗概率會(huì)隨著該子任務(wù)返工次數(shù)的增加而降低。一般來說,子任務(wù)的復(fù)雜性越高,在重復(fù)執(zhí)行時(shí)其失敗概率的降低程度越小,其失敗率降低系數(shù)越大;反之,其失敗率降低系數(shù)越小。子任務(wù)i第k次執(zhí)行時(shí),其失敗的概率可以表示為

    (7)

    式中,FPi(0)為子任務(wù)i的初始失敗率;FPLi為子任務(wù)i的失敗率降低系數(shù)。

    2 返工優(yōu)先的資源沖突處理規(guī)則

    在高端裝備研制任務(wù)執(zhí)行過程中,當(dāng)多個(gè)子任務(wù)并行執(zhí)行時(shí),可能會(huì)存在資源沖突,需要進(jìn)行資源占用的優(yōu)先級(jí)排序。在仿真模型中,假設(shè)存在一定數(shù)量的可循環(huán)使用資源(例如人員、機(jī)器設(shè)備等),子任務(wù)執(zhí)行過程中會(huì)占用一定數(shù)量的資源,子任務(wù)完工時(shí),釋放其占用的資源。

    由于資源的限制,子任務(wù)之間可能會(huì)存在資源的競爭,為了縮短研制任務(wù)工期和降低研制任務(wù)成本,需要按照一定的規(guī)則進(jìn)行資源占用的優(yōu)先級(jí)排序。本文在傳統(tǒng)方法的基礎(chǔ)上提出了“改進(jìn)的加權(quán)資源利用率與優(yōu)先關(guān)系”(improved weighted resource utilization ration and precedence,IWRURP),同時(shí)考慮子任務(wù)的緊后子任務(wù)數(shù)量、可能引發(fā)返工子任務(wù)的數(shù)量、資源占用量、工期、是否屬于返工等因素。一般來說,某子任務(wù)的緊后子任務(wù)數(shù)量和可能引發(fā)返工的子任務(wù)的數(shù)量越多,需要其信息的子任務(wù)數(shù)量就越多,該子任務(wù)需要優(yōu)先執(zhí)行。某子任務(wù)的工期越長,成為關(guān)鍵子任務(wù)的可能性越大,為了縮短總工期,其所需的資源需要優(yōu)先滿足。當(dāng)返工子任務(wù)與首次執(zhí)行子任務(wù)產(chǎn)生資源沖突時(shí),首次執(zhí)行子任務(wù)通常為當(dāng)前完工子任務(wù)的緊后子任務(wù),如果先執(zhí)行緊后子任務(wù),返工子任務(wù)的完工可能會(huì)引起緊后子任務(wù)的再次返工,因此返工子任務(wù)的優(yōu)先級(jí)應(yīng)高于首次執(zhí)行的子任務(wù)。IWRURP的計(jì)算公式為

    (8)

    式中,ωi(0<ωi<1,i=1,2,3)為各因素的權(quán)重,且ω1+ω2+ω3=1;j為子任務(wù)的編號(hào);Nj為子任務(wù)j的緊后子任務(wù)數(shù)與可能引發(fā)的返工子任務(wù)數(shù)之和;rj為子任務(wù)j的資源需求量;R為資源總量;L1、L2為優(yōu)先級(jí)因子。

    3 仿真算法設(shè)計(jì)

    基于上述模型和分析過程,采用DSM對(duì)高端裝備研制任務(wù)網(wǎng)絡(luò)進(jìn)行描述,采用離散事件仿真模型進(jìn)行仿真,使用Matlab編程實(shí)現(xiàn)。模型輸入信息為子任務(wù)工期、成本、所需資源、DSM等。某一時(shí)點(diǎn)上的子任務(wù)完工觸發(fā)狀態(tài)的轉(zhuǎn)變,狀態(tài)包括各子任務(wù)完工、正在執(zhí)行、排隊(duì)等待等。研制任務(wù)仿真的關(guān)鍵變量如表1所示,單次仿真運(yùn)行流程圖如圖1所示。

    表1 仿真過程中的關(guān)鍵變量

    圖1 單次仿真運(yùn)行流程圖Fig.1 Flow chart of a single run simulation

    每次仿真運(yùn)行,可以得到相應(yīng)的工期T和成本C,為了得到穩(wěn)定的時(shí)間和成本分布,需要多次運(yùn)行仿真。仿真不斷運(yùn)行,直到m次運(yùn)行與m-b次運(yùn)行的統(tǒng)計(jì)量的差距穩(wěn)定到一定精度α,可認(rèn)為仿真輸出結(jié)果穩(wěn)定。仿真運(yùn)行穩(wěn)定時(shí),工期T需要滿足如下表達(dá)式:

    (9)

    (10)

    對(duì)于成本C,采用同樣的方法求出其輸出結(jié)果穩(wěn)定時(shí)相應(yīng)的運(yùn)行次數(shù)。實(shí)際仿真次數(shù)取工期T和成本C仿真輸出結(jié)果穩(wěn)定時(shí)的運(yùn)行次數(shù)的最大值,即m=max(nT,nC)。

    4 算例研究

    以某無人機(jī)研制任務(wù)為例進(jìn)行算例分析,該研制任務(wù)包括14項(xiàng)子任務(wù),各子任務(wù)的工期、成本、所需資源、學(xué)習(xí)曲線及失敗概率數(shù)據(jù)如表2所示。假設(shè)各子任務(wù)的工期和成本均服從三角分布,用Tria(Do,Dm,Dp)描述各子任務(wù)工期的分布,其中Do為最樂觀工期,Dm為最可能工期,Dp為最悲觀工期;用Tria(Co,Cm,Cp)描述各子任務(wù)成本的分布,其中C0為最樂觀成本,Cm為最可能成本,Cp為最悲觀成本。在每次仿真運(yùn)行時(shí),用蒙特卡羅方法分別抽取各子任務(wù)的工期和成本。該無人機(jī)的研制任務(wù)的DSM如圖2所示,返工概率及返工影響矩陣如圖3所示,重疊比例及重疊影響矩陣如圖4所示。

    文獻(xiàn)[11]中為任務(wù)可中斷的仿真模型。任務(wù)不可中斷是指某子任務(wù)開始執(zhí)行,必須一直執(zhí)行到該子任務(wù)完工,中間不允許中斷;任務(wù)可中斷是指某子任務(wù)開始執(zhí)行后,如果有必要,在執(zhí)行過程中可以中斷一段時(shí)間,然后再繼續(xù)執(zhí)行。有些情況下,為了降低研制成本和返工次數(shù),任務(wù)可中斷是有必要的,例如圖5中:子任務(wù)T3是子任務(wù)T4和T5的緊前子任務(wù),子任務(wù)T3完工后,子任務(wù)T4和T5并行執(zhí)行,子任務(wù)T5的工期大于子任務(wù)T4的工期,子任務(wù)T4完工引發(fā)子任務(wù)T3返工。這時(shí),為了避免子任務(wù)T3完工時(shí)引發(fā)子任務(wù)T5返工,子任務(wù)T5可以暫時(shí)中斷執(zhí)行,等子任務(wù)T3完工后再繼續(xù)執(zhí)行。任務(wù)可中斷相對(duì)于任務(wù)不可中斷,可以降低成本、減少返工次數(shù),但是任務(wù)可中斷會(huì)減少子任務(wù)的并行執(zhí)行,從而導(dǎo)致整個(gè)研制任務(wù)工期增大。

    表2 研制任務(wù)的部分輸入數(shù)據(jù)

    圖2 研制任務(wù)的DSMFig.2 DSM of the development tasks

    圖3 返工概率及返工影響矩陣Fig.3 Rework probability and rework impact matrix

    圖4 重疊比例及重疊影響矩陣Fig.4 Overlap rate and overlap impact matrix

    考慮自環(huán)和失敗,在子任務(wù)可中斷條件下,研制任務(wù)仿真模型的某次運(yùn)行過程的甘特圖如圖5所示。圖中子任務(wù)T2、T8、T13在發(fā)生了自環(huán);子任務(wù)T3、T4、T5等發(fā)生了返工,例如子任務(wù)T12完工引發(fā)了T8返工,T8返工引發(fā)了T10和T11返工,T10返工又引發(fā)了T11和T12返工;子任務(wù)T5發(fā)生了中斷執(zhí)行;T6與T5、T11與T10、T12與T10重疊執(zhí)行。

    經(jīng)過多次仿真運(yùn)行測試,仿真運(yùn)行5 000次,相對(duì)精度α可達(dá)到0.000 1。為便于與文獻(xiàn)中的仿真模型進(jìn)行比較,本文首先使用了文獻(xiàn)中的數(shù)據(jù)進(jìn)行仿真,仿真結(jié)果與具有代表性的文獻(xiàn)[11,23]的仿真結(jié)果比較如表3所示,“本文(不可中斷)”是指使用本文構(gòu)建的不可中斷仿真模型;“本

    文(可中斷)”是指使用本文構(gòu)建的可中斷仿真模型;“數(shù)據(jù)1”是文獻(xiàn)中的數(shù)據(jù)(不考慮子任務(wù)發(fā)生自環(huán)和失敗);“數(shù)據(jù)2”是考慮子任務(wù)發(fā)生自環(huán)和失敗的數(shù)據(jù)。使用本文構(gòu)建的仿真模型,分別使用“數(shù)據(jù)1”和“數(shù)據(jù)2”仿真輸出的子任務(wù)執(zhí)行次數(shù)統(tǒng)計(jì)如表4所示。

    圖5 某次研制任務(wù)仿真運(yùn)行的甘特圖Fig.5 Gantt chart for the simulation of a development task

    數(shù)據(jù)模型工期/天均值標(biāo)準(zhǔn)差成本/千美元均值標(biāo)準(zhǔn)差數(shù)據(jù)1文獻(xiàn)[11]1381463763文獻(xiàn)[23]13118.5920451本文(不可中斷)135.17.8660.850.4本文(可中斷)136.18.0658.748.5數(shù)據(jù)2本文(不可中斷)144.310.0697.260.3本文(可中斷)145.510.4695.058.7

    表4 仿真輸出的各子任務(wù)平均執(zhí)行次數(shù)

    本文仿真模型允許子任務(wù)重疊,而文獻(xiàn)[11]沒有考慮子任務(wù)重疊,通過允許子任務(wù)重疊可以縮短整個(gè)研制任務(wù)工期,因此本文仿真得到工期均值小于文獻(xiàn)[11]的工期均值;但是子任務(wù)重疊會(huì)引發(fā)重疊返工,相應(yīng)會(huì)增加子任務(wù)的完成時(shí)間和成本,因此本文仿真得到的成本比文獻(xiàn)[11]的成本高。文獻(xiàn)[23]對(duì)返工進(jìn)行了限制,規(guī)定當(dāng)返工改進(jìn)程度小于一定值時(shí),則不再進(jìn)行返工,而本文規(guī)定只要引發(fā)返工就一定執(zhí)行,因此本文使用“數(shù)據(jù)1”時(shí)仿真得到的工期均值大于文獻(xiàn)[23]的工期均值。本文在仿真模型中定義了更多的運(yùn)行規(guī)則,因而本文仿真輸出數(shù)據(jù)的標(biāo)準(zhǔn)差小于文獻(xiàn)中的標(biāo)準(zhǔn)差。仿真輸出結(jié)果比較表明,本文的仿真模型能夠真實(shí)反映研制任務(wù)的運(yùn)行情況。

    高端裝備的研制成功具有重大意義,一般來說研制成功越早,越可能帶來重大價(jià)值。由表3中的仿真數(shù)據(jù)可知,任務(wù)可中斷雖然可以降低成本,但降低幅度很小,而且會(huì)增加工期,因此,對(duì)于高端裝備研制任務(wù)來說,更適用于采用任務(wù)不可中斷模式。以下使用任務(wù)不可中斷、“數(shù)據(jù)2”的仿真輸出數(shù)據(jù)進(jìn)行分析:

    (1) 自環(huán)對(duì)研制任務(wù)的影響

    對(duì)仿真輸出的研制任務(wù)成功時(shí)的數(shù)據(jù)進(jìn)行分析,由于不包括研制任務(wù)失敗時(shí)的數(shù)據(jù),相當(dāng)于消除了失敗因素的影響,單獨(dú)考慮自環(huán)因素。研制任務(wù)成功時(shí)工期均值為144.3天、成本均值為697.3千美元,分別比不考慮子任務(wù)自環(huán)和失敗時(shí)的均值增加了6.8%和5.5%,子任務(wù)總的執(zhí)行次數(shù)由21.23次增加到26.87次,增加了26.6%。數(shù)據(jù)分析表明,考慮子任務(wù)自環(huán)時(shí),研制任務(wù)工期、成本及返工次數(shù)均比不考慮自環(huán)時(shí)有所增加。

    (2) 子任務(wù)失敗對(duì)研制任務(wù)的影響

    由仿真輸出數(shù)據(jù)可知,研制任務(wù)失敗次數(shù)為490次,失敗率為9.8%,失敗時(shí)消耗的工期均值為49.9天、費(fèi)用均值為215.5千美元。如果把研制任務(wù)失敗消耗的時(shí)間和費(fèi)用分?jǐn)偟窖兄迫蝿?wù)成功的工期和成本,研制任務(wù)工期均值為149.8天、成本均值為720.6千美元,比不考慮子任務(wù)返工和失敗時(shí)分別增加了10.9%和9.0%,比單獨(dú)考慮自環(huán)時(shí)分別增加了3.8%和3.4%。考慮子任務(wù)失敗時(shí),雖然每個(gè)子任務(wù)失敗的概率非常小,但是整個(gè)研制任務(wù)的失敗風(fēng)險(xiǎn)已經(jīng)比較高;如果對(duì)研制任務(wù)失敗時(shí)消耗的時(shí)間和成本進(jìn)行分?jǐn)?會(huì)使研制任務(wù)的平均工期和成本均有較大程度的增加。

    (3) 風(fēng)險(xiǎn)分析

    風(fēng)險(xiǎn)的種類很多,這里討論任務(wù)失敗、超期、超支等風(fēng)險(xiǎn)。文獻(xiàn)[11]計(jì)劃工期和成本分別為130天和630千美元,本文考慮了自環(huán),因此需要更長的工期和更多的成本,假設(shè)計(jì)劃工期和計(jì)劃成本分別為145天和700千美元。根據(jù)仿真結(jié)果,研制任務(wù)成功時(shí)工期、費(fèi)用的頻次直方圖和累積曲線分別如圖6、圖7所示。研制任務(wù)成功時(shí),實(shí)際工期超過計(jì)劃工期的概率為43.7%,實(shí)際成本超過計(jì)劃成本的概率為42.2%,實(shí)際工期和實(shí)際成本均未超過計(jì)劃要求的概率僅為41.4%。雖然單個(gè)子任務(wù)的失敗風(fēng)險(xiǎn)很低甚至無失敗風(fēng)險(xiǎn),但是仿真結(jié)果表明的整個(gè)研制任務(wù)成功的概率僅為90.2%,研制任務(wù)成功且在計(jì)劃工期和成本內(nèi)完成的概率僅為37.3%,研制任務(wù)成功但會(huì)發(fā)生超期、超支的概率為52.9%。因此,考慮自環(huán)和失敗的情況下,高端裝備研制任務(wù)具有更高的風(fēng)險(xiǎn)。

    圖6 工期的頻次直方圖和累積曲線Fig.6 Frequency histogram and cumulative curve of the duration

    圖7 成本的頻次直方圖和累積曲線Fig.7 Frequency histogram and cumulative curve of the cost

    (4) 任務(wù)復(fù)雜性對(duì)研制任務(wù)的影響

    高端裝備研制任務(wù)的復(fù)雜性與整個(gè)任務(wù)網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)、子任務(wù)的復(fù)雜性、子任務(wù)之間關(guān)系的復(fù)雜性等因素有關(guān)。把子任務(wù)(子任務(wù)之間關(guān)系)的復(fù)雜性水平劃分為極度復(fù)雜(完全不可知)、高度復(fù)雜、中度復(fù)雜、輕度復(fù)雜、簡單化(完全確定)等五個(gè)等級(jí),不同的復(fù)雜性水平對(duì)應(yīng)不同的失敗率降低系數(shù)、學(xué)習(xí)曲線、返工概率降低系數(shù)等仿真參數(shù)。為了明顯體現(xiàn)復(fù)雜性對(duì)研制任務(wù)的影響效果,我們假定每次仿真時(shí)所有子任務(wù)的復(fù)雜性水平相同。根據(jù)子任務(wù)的復(fù)雜性水平、子任務(wù)之間關(guān)系的復(fù)雜性水平,分別給失敗率降低系數(shù)、學(xué)習(xí)曲線、返工概率降低系數(shù)等參數(shù)賦予[0,1]之間的數(shù)值。

    根據(jù)仿真數(shù)據(jù)繪制任務(wù)復(fù)雜性對(duì)研制任務(wù)影響圖如圖8所示。

    圖8 復(fù)雜性對(duì)研制任務(wù)的影響Fig.8 Influence of complexity on development tasks

    X軸表示子任務(wù)的復(fù)雜性水平,Y軸表示子任務(wù)之間關(guān)系的復(fù)雜性水平,Z軸分別表示研制任務(wù)的工期、成本、失敗率等。可見,隨著子任務(wù)復(fù)雜性水平的增高,整個(gè)研制任務(wù)的工期、成本均會(huì)增大,研制任務(wù)的失敗率也會(huì)增大。隨著子任務(wù)之間關(guān)系復(fù)雜性水平的增高,整個(gè)研制任務(wù)工期、成本均會(huì)增大,失敗率也會(huì)有一定程度的增大,但相比子任務(wù)復(fù)雜性水平影響下失敗率的增大幅度較小。同時(shí),當(dāng)子任務(wù)的復(fù)雜性水平和子任務(wù)之間關(guān)系的復(fù)雜性水平增高時(shí),整個(gè)研制任務(wù)工期、成本的不確定性程度變大,工期、成本的方差均會(huì)增大。

    5 結(jié) 論

    高端裝備研制任務(wù)中,子任務(wù)自環(huán)、返工迭代、重疊執(zhí)行、失敗以及資源沖突等因素都會(huì)對(duì)研制任務(wù)的工期、成本等產(chǎn)生影響。本文構(gòu)建了考慮多重不確定影響因素的高端裝備研制任務(wù)仿真模型,假定子任務(wù)自環(huán)、返工迭代、重疊執(zhí)行、失敗等均以隨機(jī)概率發(fā)生,并考慮了學(xué)習(xí)效應(yīng),更能有效反應(yīng)高端裝備研制任務(wù)的真實(shí)運(yùn)行情況。在此基礎(chǔ)上,根據(jù)算例仿真輸出數(shù)據(jù)對(duì)該研制任務(wù)的工期、成本、風(fēng)險(xiǎn)、復(fù)雜性影響等進(jìn)行了分析。本文提出的仿真模型可用于估計(jì)工期、成本、返工次數(shù)、失敗率、各子任務(wù)花費(fèi)的時(shí)間和成本、復(fù)雜性對(duì)研制任務(wù)的影響、資源對(duì)研制任務(wù)的影響等信息,從而能夠幫助高端裝備研制任務(wù)管理者更好地預(yù)測、評(píng)估和控制研制任務(wù)過程。

    參考文獻(xiàn):

    [1] HAHN E D. Mixture densities for project management activity times: a robust approach to PERT[J]. European Journal of Operational Research, 2008, 188(2): 450-459.

    [2] 王濤, 蔡建峰. 航空復(fù)雜產(chǎn)品項(xiàng)目活動(dòng)時(shí)間混合分布估算模型的改進(jìn)[J]. 運(yùn)籌與管理, 2012, 21(4): 214-219.

    WANG T, CAI J F. Improved mixture densities for the activity time estimation in aviation complex product project[J]. Operations Research and Management Science,2012,21(4):214-219.

    [3] WU Y, PAN X, KANG R, et al. Multi-parameters uncertainty analysis of logistic support process based on GERT[J]. Journal of Systems Engineering and Electronics, 2014, 25(6): 1011-1019.

    [4] NELSON R G, AZARON A, AREF S. The use of a GERT based method to model concurrent product development processes[J]. European Journal of Operational Research,2016,250(2):566-578.

    [5] 劉紅旗, 方志耕, 陶良彥. 復(fù)雜裝備研制項(xiàng)目進(jìn)度規(guī)劃GERT網(wǎng)絡(luò)“反問題”模型[J]. 系統(tǒng)工程與電子技術(shù), 2015, 37(12): 2758-2763.

    LIU H Q, FANG Z G, TAO Y L. Complex equipment development project planning GERT network “inverse problem” model[J]. Systems Engineering and Electronics, 2015, 37(12): 2758-2763.

    [6] 陶良彥,劉思峰,方志耕,等.基于多層次GERT的復(fù)雜產(chǎn)品研制進(jìn)度規(guī)劃“超沖突均衡”博弈模型[J].控制與決策,2014,29(11): 2002-2010.

    TAO Y L, LIU S F, FANG Z G, et al. “Super-conflict equalization” game model for complex product development scheduling based on multi-level GERT[J].Control and Decision,2014,29(11):2002-2010.

    [7] STEWARD D V. The design structure system: a method for managing the design of complex systems[J]. IEEE Trans.on Engineering Management, 1981,28(3): 71-74.

    [8] FANG C, MARLE F. Dealing with project complexity by matrix-based propagation modelling for project risk analysis[J]. Journal of Engineering Design, 2013, 24(4): 239-256.

    [9] BROWNING T R. Applying the design structure matrix to system decomposition and integration problems: a review and new directions[J]. IEEE Trans.on Engineering Management, 2001, 48(3): 292-306.

    [10] BROWNING T R. Design structure matrix extensions and innovations: a survey and new opportunities[J]. IEEE Trans.on Engineering Management, 2016, 63(1): 27-52.

    [11] BROWNING T R, EPPINGER S D. Modeling impacts of process architecture on cost and schedule risk in product development[J]. IEEE Trans.on Engineering Management, 2002, 49(4): 428-443.

    [12] LIN J, CHAI K H, WONG Y S, et al. A dynamic model for managing overlapped iterative product development[J]. European Journal of Operational Research, 2008, 185(1): 378-392.

    [13] LIN J, CHAI K H, BROMBACHER A C, et al. Optimal overlapping and functional interaction in product development[J]. European Journal of Operational Research,2009,196(3):1158-1169.

    [14] LIN J, QIAN Y, CUI W, et al. Overlapping and communication policies in product development[J]. European Journal of Operational Research, 2010, 201(3): 737-750.

    [15] YANG Q, YAO T, LU T, et al. An overlapping-based design structure matrix for measuring interaction strength and clustering analysis in product development project[J]. IEEE Trans.on Engineering Management, 2014, 61(1): 159-170.

    [16] ZHANG D, BHUIYAN N. A study of the evolution of uncertainty in product development as a basis for overlapping[J]. IEEE Trans.on Engineering Management, 2015, 62(1): 39-50.

    [17] WYNN D C, ECKERT C M. Perspectives on iteration in design and development[J]. Research in Engineering Design, 2017, 28(2): 153-184.

    [18] MEIER C, BROWNING T R, YASSINE A A, et al. The cost of speed: work policies for crashing and overlapping in product development projects[J]. IEEE Trans.on Engineering Management, 2015, 62(2): 237-255.

    [19] NASR W, YASSINE A, KASM O A. An analytical approach to estimate the expected duration and variance for iterative product development projects[J]. Research in Engineering Design, 2016, 27(1): 55-71.

    [20] REBENTISCH E, SCHUH G, RIESENER M, et al. Assessment of changes in technical systems and their effects on cost and duration based on structural complexity[J]. Procedia Cirp, 2016, 55: 35-40.

    [21] 楊青, 呂杰峰. 基于DSM返工風(fēng)險(xiǎn)評(píng)價(jià)矩陣的項(xiàng)目優(yōu)化與仿真[J]. 系統(tǒng)工程理論與實(shí)踐, 2010, 30(9): 1665-1671.

    YANG Q, Lü J F. Project optimization and simulation based on DSM rework risk evaluation matrix[J]. Systems Engineering-Theory & Practice, 2010, 30(9): 1665-1671.

    [22] COOK I, COATES G. Optimising the time-based design structure matrix using a divide and hybridise algorithm[J]. Journal of Engineering Design, 2016, 27(4/6): 306-332.

    [23] 李洪波,徐哲.考慮活動(dòng)隨機(jī)重疊和資源沖突的復(fù)雜產(chǎn)品開發(fā)流程仿真建模[J].系統(tǒng)工程與電子技術(shù),2012,34(7): 1412-1418.

    LI H B, XU Z. Simulation modeling of complex product deve-lopment process considering stochastic activities overlap and resource conflict[J]. Systems Engineering and Electronics, 2012, 34(7): 1412-1418.

    [24] ZHANG J, XU Z, LI Y, et al. Framework for the integration of assembly modeling and simulation based on assembly feature pair[J]. International Journal of Advanced Manufacturing Technology, 2015, 78(5/8): 765-780.

    [25] LUO J. A simulation-based method to evaluate the impact of product architecture on product evolvability[J]. Research in Engineering Design, 2015, 26(4): 355-371.

    [26] COLLINS S T, YASSINE A A, BORGATTI S P. Evaluating product development systems using network analysis[J]. Systems Engineering, 2009, 12(1): 55-68.

    [27] KARNIEL A, REICH Y. Multi-level modelling and simulation of new product development processes[J]. Journal of Engineering Design, 2013, 24(3): 185-210.

    [28] MAIER J F, WYNN D C, BIEDERMANN W, et al. Simulating progressive iteration, rework and change propagation to prioritise design tasks[J]. Research in Engineering Design, 2014, 25(4): 283-307.

    [29] CHO S H, EPPINGER S D. A simulation-based process model for managing complex design projects[J]. IEEE Trans.on Engineering Management, 2005, 52(3): 316-328.

    [30] ABDELSALAM H M, RASMY M H, MOHAMED H G. A simulation-based time reduction approach for resource constrained design structure matrix[J]. International Journal of Modeling & Optimization, 2014, 4(1): 51-55.

    猜你喜歡
    成本模型
    一半模型
    2021年最新酒駕成本清單
    河南電力(2021年5期)2021-05-29 02:10:00
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    溫子仁,你還是適合拍小成本
    電影(2018年12期)2018-12-23 02:18:48
    鄉(xiāng)愁的成本
    特別健康(2018年2期)2018-06-29 06:13:42
    “二孩補(bǔ)貼”難抵養(yǎng)娃成本
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    獨(dú)聯(lián)體各國的勞動(dòng)力成本
    波多野结衣av一区二区av| 丝袜人妻中文字幕| 一边亲一边摸免费视频| 天天躁日日躁夜夜躁夜夜| 国产爽快片一区二区三区| 国产精品国产av在线观看| 在线看a的网站| 久久久久精品久久久久真实原创| 国产成人免费无遮挡视频| 亚洲欧美中文字幕日韩二区| 国产成人系列免费观看| 国产伦人伦偷精品视频| 日韩大片免费观看网站| 最近手机中文字幕大全| 亚洲少妇的诱惑av| 日韩熟女老妇一区二区性免费视频| 亚洲美女黄色视频免费看| 伦理电影免费视频| 国产乱来视频区| av视频免费观看在线观看| 国产成人精品久久二区二区91 | 精品久久久久久电影网| 无限看片的www在线观看| 黑人巨大精品欧美一区二区蜜桃| 欧美最新免费一区二区三区| 看非洲黑人一级黄片| 人人澡人人妻人| 亚洲熟女毛片儿| 另类亚洲欧美激情| 亚洲精品,欧美精品| 久久久久精品国产欧美久久久 | 2018国产大陆天天弄谢| 一二三四在线观看免费中文在| 国产精品香港三级国产av潘金莲 | xxx大片免费视频| 午夜91福利影院| 自线自在国产av| 91精品三级在线观看| 亚洲精品,欧美精品| 久久久久人妻精品一区果冻| 少妇 在线观看| 午夜免费观看性视频| 高清黄色对白视频在线免费看| 久久久久视频综合| av女优亚洲男人天堂| 免费少妇av软件| av.在线天堂| 99久久精品国产亚洲精品| 色播在线永久视频| 免费黄频网站在线观看国产| 中文字幕另类日韩欧美亚洲嫩草| 一级片'在线观看视频| 2021少妇久久久久久久久久久| 男女午夜视频在线观看| 亚洲欧美精品自产自拍| 夫妻午夜视频| 丰满饥渴人妻一区二区三| 18禁国产床啪视频网站| 欧美日韩成人在线一区二区| 最近中文字幕高清免费大全6| 免费看av在线观看网站| 国产在线免费精品| 69精品国产乱码久久久| 国产免费视频播放在线视频| 女人高潮潮喷娇喘18禁视频| 满18在线观看网站| 国产亚洲午夜精品一区二区久久| 一区二区三区激情视频| 久久精品国产亚洲av涩爱| 中文欧美无线码| 青青草视频在线视频观看| 久久久久久免费高清国产稀缺| 久久av网站| 超色免费av| 成人漫画全彩无遮挡| 丰满饥渴人妻一区二区三| 99re6热这里在线精品视频| 超色免费av| 欧美日韩视频高清一区二区三区二| 精品国产超薄肉色丝袜足j| 国产毛片在线视频| 操出白浆在线播放| 日本爱情动作片www.在线观看| 国产1区2区3区精品| 国产成人a∨麻豆精品| 丁香六月欧美| 免费在线观看视频国产中文字幕亚洲 | 国产精品99久久99久久久不卡 | 波多野结衣一区麻豆| 亚洲精品乱久久久久久| 色视频在线一区二区三区| 在线 av 中文字幕| 久久久久久人人人人人| 国产xxxxx性猛交| 日韩,欧美,国产一区二区三区| 男的添女的下面高潮视频| 免费在线观看黄色视频的| 爱豆传媒免费全集在线观看| 日日爽夜夜爽网站| 免费黄色在线免费观看| 午夜福利视频在线观看免费| 免费女性裸体啪啪无遮挡网站| 成年女人毛片免费观看观看9 | 成人影院久久| 久久人人爽人人片av| 亚洲男人天堂网一区| 男的添女的下面高潮视频| 欧美另类一区| 美女视频免费永久观看网站| 51午夜福利影视在线观看| 国产精品免费大片| 欧美 日韩 精品 国产| 亚洲成色77777| 日韩不卡一区二区三区视频在线| 成年人午夜在线观看视频| 国产精品熟女久久久久浪| 亚洲精品久久成人aⅴ小说| 十八禁人妻一区二区| 日本av手机在线免费观看| 在线亚洲精品国产二区图片欧美| 黄频高清免费视频| 成人漫画全彩无遮挡| 纯流量卡能插随身wifi吗| 亚洲成人av在线免费| 中文字幕色久视频| 99久久综合免费| 亚洲欧美色中文字幕在线| 男的添女的下面高潮视频| 国产亚洲最大av| 飞空精品影院首页| 男女下面插进去视频免费观看| 激情五月婷婷亚洲| 中文字幕亚洲精品专区| 亚洲av欧美aⅴ国产| 激情五月婷婷亚洲| 国产乱人偷精品视频| 国产成人精品久久二区二区91 | 久久久久久久国产电影| 成年美女黄网站色视频大全免费| 黄片播放在线免费| 免费黄网站久久成人精品| 国产免费现黄频在线看| 久久久国产一区二区| 亚洲精品成人av观看孕妇| 日本黄色日本黄色录像| 2018国产大陆天天弄谢| 熟妇人妻不卡中文字幕| 久久国产精品大桥未久av| 色综合欧美亚洲国产小说| 久久久久国产精品人妻一区二区| 波野结衣二区三区在线| 国产精品欧美亚洲77777| 亚洲欧美激情在线| 成人亚洲欧美一区二区av| 午夜激情久久久久久久| 女人被躁到高潮嗷嗷叫费观| 电影成人av| 久久99热这里只频精品6学生| 新久久久久国产一级毛片| 亚洲av欧美aⅴ国产| 欧美日韩视频高清一区二区三区二| www.自偷自拍.com| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜爱| 一个人免费看片子| 国产黄色免费在线视频| 国产1区2区3区精品| 日本欧美视频一区| 婷婷色综合大香蕉| 国产 精品1| 久久亚洲国产成人精品v| 别揉我奶头~嗯~啊~动态视频 | 亚洲美女搞黄在线观看| 国产女主播在线喷水免费视频网站| av天堂久久9| www.自偷自拍.com| 男女边摸边吃奶| 久久热在线av| 久久久国产欧美日韩av| 一级,二级,三级黄色视频| 综合色丁香网| 亚洲成人免费av在线播放| 老鸭窝网址在线观看| 人人妻人人添人人爽欧美一区卜| 极品少妇高潮喷水抽搐| 天天躁夜夜躁狠狠久久av| 久久青草综合色| 亚洲成人免费av在线播放| 亚洲精品久久久久久婷婷小说| 嫩草影院入口| a级片在线免费高清观看视频| 99精品久久久久人妻精品| 又黄又粗又硬又大视频| 啦啦啦 在线观看视频| 大码成人一级视频| 久久久久久人妻| 精品免费久久久久久久清纯 | 啦啦啦啦在线视频资源| 人人妻人人澡人人爽人人夜夜| √禁漫天堂资源中文www| 飞空精品影院首页| 麻豆av在线久日| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 哪个播放器可以免费观看大片| 夜夜骑夜夜射夜夜干| 宅男免费午夜| av又黄又爽大尺度在线免费看| 韩国高清视频一区二区三区| 久久精品国产亚洲av高清一级| 国产精品女同一区二区软件| 国产色婷婷99| 美女主播在线视频| 久久韩国三级中文字幕| 国产黄色视频一区二区在线观看| 亚洲图色成人| 18禁国产床啪视频网站| 精品福利永久在线观看| 高清av免费在线| av女优亚洲男人天堂| 精品一区二区三卡| 十八禁网站网址无遮挡| 汤姆久久久久久久影院中文字幕| 少妇的丰满在线观看| 亚洲四区av| 女人高潮潮喷娇喘18禁视频| 叶爱在线成人免费视频播放| 日韩一卡2卡3卡4卡2021年| 91精品伊人久久大香线蕉| 中文字幕亚洲精品专区| 国产一区二区三区综合在线观看| 亚洲av中文av极速乱| 欧美少妇被猛烈插入视频| 侵犯人妻中文字幕一二三四区| 久久久欧美国产精品| 亚洲激情五月婷婷啪啪| 激情五月婷婷亚洲| 十八禁高潮呻吟视频| 超色免费av| 无限看片的www在线观看| 中文字幕人妻丝袜一区二区 | 亚洲精华国产精华液的使用体验| 国产男人的电影天堂91| 欧美 日韩 精品 国产| 亚洲av欧美aⅴ国产| 嫩草影院入口| 黄色一级大片看看| 成人影院久久| 精品人妻在线不人妻| 精品国产一区二区久久| 日韩中文字幕欧美一区二区 | 又粗又硬又长又爽又黄的视频| 纯流量卡能插随身wifi吗| 1024视频免费在线观看| 久久久久精品性色| 一级毛片我不卡| a级片在线免费高清观看视频| 国产成人欧美在线观看 | 国产成人精品在线电影| 亚洲 欧美一区二区三区| 久久久精品免费免费高清| e午夜精品久久久久久久| 国产精品免费视频内射| 国产福利在线免费观看视频| 丰满少妇做爰视频| 午夜福利免费观看在线| 看免费av毛片| 亚洲国产精品999| 可以免费在线观看a视频的电影网站 | 18在线观看网站| 99久国产av精品国产电影| 黄色一级大片看看| 久久综合国产亚洲精品| 精品一区二区三卡| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线| 香蕉国产在线看| a 毛片基地| 美女脱内裤让男人舔精品视频| 午夜免费鲁丝| 久久亚洲国产成人精品v| 国产亚洲欧美精品永久| 乱人伦中国视频| 久久性视频一级片| 黄色视频在线播放观看不卡| 国产成人精品久久久久久| 久久国产亚洲av麻豆专区| 国产一卡二卡三卡精品 | 乱人伦中国视频| 亚洲色图 男人天堂 中文字幕| 美女国产高潮福利片在线看| 在线看a的网站| 欧美日韩国产mv在线观看视频| 久久鲁丝午夜福利片| 免费av中文字幕在线| 丝袜喷水一区| 99热网站在线观看| 亚洲伊人色综图| 青草久久国产| 热re99久久国产66热| 亚洲精品一区蜜桃| 国产精品av久久久久免费| 各种免费的搞黄视频| 99久国产av精品国产电影| 另类亚洲欧美激情| 69精品国产乱码久久久| 国产精品秋霞免费鲁丝片| 国产1区2区3区精品| 欧美亚洲日本最大视频资源| 日日啪夜夜爽| 亚洲第一区二区三区不卡| 一级片'在线观看视频| 日韩一区二区三区影片| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 免费观看性生交大片5| 久久精品国产a三级三级三级| 满18在线观看网站| 久久综合国产亚洲精品| 免费黄网站久久成人精品| 精品国产露脸久久av麻豆| 美女扒开内裤让男人捅视频| 韩国高清视频一区二区三区| 五月开心婷婷网| 日韩中文字幕视频在线看片| 久久久精品区二区三区| 99国产综合亚洲精品| 伦理电影免费视频| 伊人久久国产一区二区| 欧美变态另类bdsm刘玥| 晚上一个人看的免费电影| 最近最新中文字幕大全免费视频 | 免费人妻精品一区二区三区视频| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三卡| av福利片在线| 超碰成人久久| 伦理电影大哥的女人| 夫妻午夜视频| 18禁观看日本| 天天操日日干夜夜撸| 久久狼人影院| 午夜老司机福利片| 黄色视频在线播放观看不卡| 高清视频免费观看一区二区| 亚洲精品乱久久久久久| 亚洲欧美中文字幕日韩二区| 国产不卡av网站在线观看| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 午夜免费观看性视频| 日韩不卡一区二区三区视频在线| 五月天丁香电影| 成年人免费黄色播放视频| 各种免费的搞黄视频| 天堂俺去俺来也www色官网| 日韩熟女老妇一区二区性免费视频| 亚洲国产中文字幕在线视频| 欧美 日韩 精品 国产| 国产福利在线免费观看视频| 丰满迷人的少妇在线观看| videos熟女内射| 日韩av在线免费看完整版不卡| 国产成人a∨麻豆精品| 人妻 亚洲 视频| av又黄又爽大尺度在线免费看| 丝瓜视频免费看黄片| 又大又爽又粗| netflix在线观看网站| 国产成人av激情在线播放| 精品少妇一区二区三区视频日本电影 | 狠狠精品人妻久久久久久综合| 欧美精品人与动牲交sv欧美| 热re99久久国产66热| 国产成人午夜福利电影在线观看| 日本av免费视频播放| 超碰成人久久| 日韩一区二区视频免费看| 亚洲天堂av无毛| av电影中文网址| 亚洲,欧美精品.| 色婷婷av一区二区三区视频| 狂野欧美激情性bbbbbb| videosex国产| 国产亚洲欧美精品永久| www.自偷自拍.com| 交换朋友夫妻互换小说| av卡一久久| 男女无遮挡免费网站观看| 中文字幕人妻熟女乱码| 国产老妇伦熟女老妇高清| 国产欧美亚洲国产| www.精华液| 韩国精品一区二区三区| 麻豆av在线久日| 国产在视频线精品| 黄色一级大片看看| 制服人妻中文乱码| 国产又爽黄色视频| 婷婷色综合大香蕉| 日韩一区二区三区影片| 国产伦人伦偷精品视频| 9热在线视频观看99| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 亚洲精品久久午夜乱码| 波野结衣二区三区在线| 免费日韩欧美在线观看| 久久免费观看电影| 另类精品久久| 久久这里只有精品19| 精品亚洲乱码少妇综合久久| 在线观看三级黄色| kizo精华| 最近中文字幕高清免费大全6| 国产在线免费精品| 亚洲精品国产av成人精品| 男女免费视频国产| 韩国av在线不卡| 亚洲熟女精品中文字幕| 人人妻,人人澡人人爽秒播 | 国产人伦9x9x在线观看| 亚洲精品美女久久久久99蜜臀 | 熟妇人妻不卡中文字幕| 日韩大片免费观看网站| 午夜福利乱码中文字幕| 久久人人爽av亚洲精品天堂| videosex国产| 国产精品嫩草影院av在线观看| 一级片免费观看大全| 亚洲欧美激情在线| 日本一区二区免费在线视频| 亚洲专区中文字幕在线 | 亚洲精品久久成人aⅴ小说| 国产视频首页在线观看| 九草在线视频观看| 免费在线观看视频国产中文字幕亚洲 | 国产一区有黄有色的免费视频| 日本猛色少妇xxxxx猛交久久| 国产一卡二卡三卡精品 | 亚洲成人国产一区在线观看 | 在线亚洲精品国产二区图片欧美| 在线观看免费午夜福利视频| 欧美av亚洲av综合av国产av | 日韩,欧美,国产一区二区三区| 色播在线永久视频| 午夜福利网站1000一区二区三区| 国产精品一二三区在线看| 一区二区三区精品91| 久久久久久久国产电影| 国产一区亚洲一区在线观看| 亚洲成人国产一区在线观看 | 99久久人妻综合| 最近最新中文字幕免费大全7| 嫩草影视91久久| 亚洲av成人精品一二三区| 无遮挡黄片免费观看| 夫妻午夜视频| 19禁男女啪啪无遮挡网站| 日韩欧美精品免费久久| a级毛片在线看网站| 亚洲精品中文字幕在线视频| 色婷婷av一区二区三区视频| 久热这里只有精品99| 黄色怎么调成土黄色| 天天添夜夜摸| 最近2019中文字幕mv第一页| 日本午夜av视频| 精品亚洲乱码少妇综合久久| √禁漫天堂资源中文www| 99re6热这里在线精品视频| netflix在线观看网站| 精品福利永久在线观看| 国产精品久久久久久久久免| 成人国语在线视频| 久久毛片免费看一区二区三区| av在线播放精品| 天天添夜夜摸| 少妇人妻精品综合一区二区| 日韩伦理黄色片| 在线观看一区二区三区激情| 国产日韩一区二区三区精品不卡| 老司机影院毛片| 午夜激情av网站| 国产熟女午夜一区二区三区| 国产亚洲av片在线观看秒播厂| 久久久精品94久久精品| 19禁男女啪啪无遮挡网站| 超碰成人久久| 国产精品无大码| 国产日韩欧美亚洲二区| 国产视频首页在线观看| 国产又色又爽无遮挡免| 中文字幕最新亚洲高清| 久久久久久久久免费视频了| 在线观看人妻少妇| 在线观看免费午夜福利视频| 一边摸一边做爽爽视频免费| 亚洲精品视频女| 一二三四在线观看免费中文在| 中文天堂在线官网| 中文乱码字字幕精品一区二区三区| 国产av码专区亚洲av| 咕卡用的链子| 国产免费现黄频在线看| 亚洲综合精品二区| 在线天堂最新版资源| 最近手机中文字幕大全| 美女脱内裤让男人舔精品视频| 久久 成人 亚洲| 久久天堂一区二区三区四区| 黑人猛操日本美女一级片| 午夜福利乱码中文字幕| 色综合欧美亚洲国产小说| 美女午夜性视频免费| av网站在线播放免费| 美女中出高潮动态图| 国产精品免费大片| 久久狼人影院| 午夜久久久在线观看| 美女中出高潮动态图| 色综合欧美亚洲国产小说| 男女无遮挡免费网站观看| av网站在线播放免费| 日韩电影二区| 国产亚洲欧美精品永久| 亚洲国产毛片av蜜桃av| 777米奇影视久久| 日韩一区二区视频免费看| svipshipincom国产片| 国产一区二区三区综合在线观看| 在现免费观看毛片| 最黄视频免费看| 精品视频人人做人人爽| 卡戴珊不雅视频在线播放| 美女大奶头黄色视频| 成年av动漫网址| 免费av中文字幕在线| 久久女婷五月综合色啪小说| 韩国av在线不卡| 七月丁香在线播放| 高清视频免费观看一区二区| 99久国产av精品国产电影| 亚洲色图 男人天堂 中文字幕| 国产日韩欧美在线精品| av国产精品久久久久影院| 永久免费av网站大全| 亚洲国产成人一精品久久久| 国产亚洲午夜精品一区二区久久| 久久女婷五月综合色啪小说| 国产日韩欧美视频二区| 日日撸夜夜添| 亚洲成av片中文字幕在线观看| 久久精品亚洲熟妇少妇任你| 欧美日韩一级在线毛片| 久热这里只有精品99| 好男人视频免费观看在线| 国产精品国产三级专区第一集| 丁香六月天网| 久久久久久久精品精品| 精品福利永久在线观看| 99久久人妻综合| 制服人妻中文乱码| 成年美女黄网站色视频大全免费| 欧美日韩亚洲高清精品| 夫妻午夜视频| 精品国产一区二区三区久久久樱花| 久久久精品区二区三区| 色综合欧美亚洲国产小说| 五月开心婷婷网| 亚洲精品第二区| 亚洲一级一片aⅴ在线观看| 精品少妇黑人巨大在线播放| 亚洲国产精品成人久久小说| 一区二区av电影网| 亚洲视频免费观看视频| 超碰成人久久| 日本一区二区免费在线视频| 免费看不卡的av| 亚洲欧美清纯卡通| 99精品久久久久人妻精品| 80岁老熟妇乱子伦牲交| 日韩欧美精品免费久久| 国产 一区精品| 精品一区二区三卡| 99久久综合免费| 青草久久国产| 久久婷婷青草| 丰满少妇做爰视频| 一区在线观看完整版| 成人午夜精彩视频在线观看| av电影中文网址| 美女主播在线视频| 国产又色又爽无遮挡免| 不卡av一区二区三区| 视频区图区小说| 亚洲第一青青草原| a级毛片黄视频| 亚洲欧美成人综合另类久久久| 亚洲精华国产精华液的使用体验| 不卡av一区二区三区| 少妇人妻 视频| 国产成人一区二区在线| 久久精品亚洲熟妇少妇任你| 亚洲美女视频黄频| 深夜精品福利| 九色亚洲精品在线播放| 中国三级夫妇交换| 五月开心婷婷网| 国产亚洲一区二区精品| 亚洲av福利一区| 亚洲人成77777在线视频| 亚洲精品中文字幕在线视频| svipshipincom国产片| 国产伦人伦偷精品视频|