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

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

    2018-06-07 08:07:15張西林譚躍進楊志偉
    系統(tǒng)工程與電子技術 2018年6期
    關鍵詞:成本模型

    張西林, 譚躍進, 楊志偉

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

    0 引 言

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

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

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

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

    1 高端裝備研制任務仿真問題及子任務關系描述

    1.1 問題描述

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

    1.2 DSM、返工、重疊及學習曲線的量化表示

    1.2.1 DSM的量化表示

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

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

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

    (1)

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

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

    1.2.3 學習曲線

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

    (2)

    (3)

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

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

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

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

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

    OAij=min(Di×ORij,Dj)

    (4)

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

    ADi=Di+RWij

    (5)

    (6)

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

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

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

    1.4 子任務失敗的量化表示

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

    (7)

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

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

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

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

    (8)

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

    3 仿真算法設計

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

    表1 仿真過程中的關鍵變量

    圖1 單次仿真運行流程圖Fig.1 Flow chart of a single run simulation

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

    (9)

    (10)

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

    4 算例研究

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

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

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

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

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

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

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

    數(shù)據(jù)模型工期/天均值標準差成本/千美元均值標準差數(shù)據(jù)1文獻[11]1381463763文獻[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 仿真輸出的各子任務平均執(zhí)行次數(shù)

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

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

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

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

    (2) 子任務失敗對研制任務的影響

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

    (3) 風險分析

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

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

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

    (4) 任務復雜性對研制任務的影響

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

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

    圖8 復雜性對研制任務的影響Fig.8 Influence of complexity on development tasks

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

    5 結 論

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

    參考文獻:

    [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] 王濤, 蔡建峰. 航空復雜產品項目活動時間混合分布估算模型的改進[J]. 運籌與管理, 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] 劉紅旗, 方志耕, 陶良彥. 復雜裝備研制項目進度規(guī)劃GERT網絡“反問題”模型[J]. 系統(tǒng)工程與電子技術, 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的復雜產品研制進度規(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返工風險評價矩陣的項目優(yōu)化與仿真[J]. 系統(tǒng)工程理論與實踐, 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] 李洪波,徐哲.考慮活動隨機重疊和資源沖突的復雜產品開發(fā)流程仿真建模[J].系統(tǒng)工程與電子技術,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
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    溫子仁,你還是適合拍小成本
    電影(2018年12期)2018-12-23 02:18:48
    鄉(xiāng)愁的成本
    特別健康(2018年2期)2018-06-29 06:13:42
    “二孩補貼”難抵養(yǎng)娃成本
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉換方法初步研究
    獨聯(lián)體各國的勞動力成本
    亚洲 国产 在线| 后天国语完整版免费观看| av欧美777| 亚洲中文字幕日韩| 韩国高清视频一区二区三区| 国产高清国产精品国产三级| 成人手机av| 成年美女黄网站色视频大全免费| 丁香六月天网| 亚洲五月色婷婷综合| 国产成人免费无遮挡视频| 精品少妇久久久久久888优播| 后天国语完整版免费观看| 捣出白浆h1v1| 少妇粗大呻吟视频| 咕卡用的链子| 欧美日韩亚洲国产一区二区在线观看 | 欧美 日韩 精品 国产| 91老司机精品| 日本vs欧美在线观看视频| 永久免费av网站大全| 久久天堂一区二区三区四区| 天天躁夜夜躁狠狠躁躁| 男人爽女人下面视频在线观看| 制服人妻中文乱码| 亚洲人成电影观看| 成人亚洲精品一区在线观看| 日韩欧美一区视频在线观看| 日韩制服丝袜自拍偷拍| 欧美国产精品va在线观看不卡| 亚洲国产精品一区三区| av又黄又爽大尺度在线免费看| 欧美人与性动交α欧美精品济南到| 中文字幕另类日韩欧美亚洲嫩草| 精品视频人人做人人爽| 91字幕亚洲| 最近手机中文字幕大全| 看免费av毛片| 99热国产这里只有精品6| 色网站视频免费| 亚洲精品一区蜜桃| 97人妻天天添夜夜摸| 黄片播放在线免费| 一区在线观看完整版| 性高湖久久久久久久久免费观看| 亚洲精品一区蜜桃| 亚洲专区国产一区二区| 国产1区2区3区精品| 国产亚洲午夜精品一区二区久久| 青草久久国产| 国产亚洲精品久久久久5区| 国产一卡二卡三卡精品| 极品少妇高潮喷水抽搐| 亚洲精品第二区| 咕卡用的链子| 老司机深夜福利视频在线观看 | 欧美日韩成人在线一区二区| 久久鲁丝午夜福利片| 乱人伦中国视频| 国产精品久久久久久精品古装| 欧美国产精品va在线观看不卡| 精品国产国语对白av| 国产av国产精品国产| 国产一区二区 视频在线| 久久久国产一区二区| av片东京热男人的天堂| 亚洲精品日韩在线中文字幕| 国产av精品麻豆| 久久久久精品人妻al黑| 久热这里只有精品99| 国产在线观看jvid| 国产一级毛片在线| 一边摸一边抽搐一进一出视频| 日韩大片免费观看网站| 国产欧美日韩综合在线一区二区| 欧美成狂野欧美在线观看| 亚洲av国产av综合av卡| 亚洲av欧美aⅴ国产| 男女免费视频国产| 丝袜美腿诱惑在线| 亚洲精品美女久久久久99蜜臀 | 十分钟在线观看高清视频www| 免费在线观看黄色视频的| 视频在线观看一区二区三区| av网站免费在线观看视频| 91麻豆精品激情在线观看国产 | 黄色视频在线播放观看不卡| 好男人视频免费观看在线| 国产精品一区二区免费欧美 | 男女高潮啪啪啪动态图| 少妇的丰满在线观看| 亚洲国产av新网站| 亚洲国产中文字幕在线视频| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美成人综合另类久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩一本色道免费dvd| 婷婷色麻豆天堂久久| 在线观看免费日韩欧美大片| 两性夫妻黄色片| 男女下面插进去视频免费观看| 交换朋友夫妻互换小说| 1024视频免费在线观看| 1024视频免费在线观看| 久久精品国产亚洲av高清一级| 大香蕉久久网| 欧美在线黄色| 五月天丁香电影| av网站在线播放免费| 久久精品国产亚洲av高清一级| 久久久精品免费免费高清| 狂野欧美激情性bbbbbb| 欧美大码av| 欧美日韩视频精品一区| 国产真人三级小视频在线观看| 婷婷色综合www| 日韩大片免费观看网站| 亚洲一码二码三码区别大吗| 国产又爽黄色视频| 在线 av 中文字幕| 成年女人毛片免费观看观看9 | 国产亚洲欧美在线一区二区| 建设人人有责人人尽责人人享有的| 精品久久久久久电影网| 色94色欧美一区二区| av有码第一页| 女性生殖器流出的白浆| 午夜福利,免费看| 黄片小视频在线播放| 久久精品熟女亚洲av麻豆精品| 青春草视频在线免费观看| 十八禁高潮呻吟视频| 免费看av在线观看网站| 一级a爱视频在线免费观看| 高清欧美精品videossex| 中文字幕av电影在线播放| 精品一区在线观看国产| 亚洲国产精品一区三区| 日日摸夜夜添夜夜爱| 久久久久久久久免费视频了| 国产视频首页在线观看| 2018国产大陆天天弄谢| 两人在一起打扑克的视频| 日韩大码丰满熟妇| 天天操日日干夜夜撸| 久久久欧美国产精品| 久久久精品免费免费高清| 国产精品国产av在线观看| 国产精品香港三级国产av潘金莲 | 日韩伦理黄色片| 成人黄色视频免费在线看| 欧美日韩一级在线毛片| 国产欧美日韩综合在线一区二区| 亚洲欧美成人综合另类久久久| 日韩大片免费观看网站| 多毛熟女@视频| 久久人人爽av亚洲精品天堂| 黄色a级毛片大全视频| 国产视频首页在线观看| 夫妻午夜视频| 国产亚洲午夜精品一区二区久久| 国产视频首页在线观看| 久久久久久久久免费视频了| 日日爽夜夜爽网站| 一级毛片我不卡| 中文字幕精品免费在线观看视频| 国产在视频线精品| www日本在线高清视频| 欧美日韩av久久| 国产日韩欧美视频二区| 国产精品一二三区在线看| 国产无遮挡羞羞视频在线观看| 飞空精品影院首页| 成人国产一区最新在线观看 | 国产精品欧美亚洲77777| 少妇人妻久久综合中文| 成人国产av品久久久| 自拍欧美九色日韩亚洲蝌蚪91| 日韩制服骚丝袜av| 久久这里只有精品19| 亚洲精品国产一区二区精华液| 欧美成狂野欧美在线观看| 女性生殖器流出的白浆| 午夜免费成人在线视频| 2021少妇久久久久久久久久久| 一边亲一边摸免费视频| 国产精品偷伦视频观看了| 国产精品av久久久久免费| 久久久精品免费免费高清| 久久久久久免费高清国产稀缺| 午夜91福利影院| 欧美另类一区| 亚洲精品日韩在线中文字幕| 午夜福利免费观看在线| 久久精品国产亚洲av高清一级| 久久精品成人免费网站| 国产熟女午夜一区二区三区| 中文字幕色久视频| 国产亚洲午夜精品一区二区久久| 国产99久久九九免费精品| 老司机深夜福利视频在线观看 | 日本欧美国产在线视频| 亚洲专区国产一区二区| 激情视频va一区二区三区| 韩国精品一区二区三区| 欧美精品高潮呻吟av久久| 国产又色又爽无遮挡免| 女性生殖器流出的白浆| av网站免费在线观看视频| 别揉我奶头~嗯~啊~动态视频 | 99热网站在线观看| 国产视频首页在线观看| 中文字幕高清在线视频| 欧美精品亚洲一区二区| 一级a爱视频在线免费观看| 日日摸夜夜添夜夜爱| 99国产精品一区二区蜜桃av | 日本91视频免费播放| 777久久人妻少妇嫩草av网站| 高清不卡的av网站| 婷婷色综合www| 亚洲精品一卡2卡三卡4卡5卡 | 精品一区在线观看国产| 国产成人av激情在线播放| 久久国产亚洲av麻豆专区| 男人爽女人下面视频在线观看| 国产黄色视频一区二区在线观看| 韩国高清视频一区二区三区| 婷婷色av中文字幕| 国产野战对白在线观看| 国产伦理片在线播放av一区| 久久女婷五月综合色啪小说| 日韩av不卡免费在线播放| 中文欧美无线码| 一区二区av电影网| 欧美日韩亚洲国产一区二区在线观看 | 久久综合国产亚洲精品| 国产av国产精品国产| 不卡av一区二区三区| 手机成人av网站| 国产爽快片一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲国产精品999| 人妻一区二区av| 男女免费视频国产| 99热国产这里只有精品6| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩黄片免| 男人舔女人的私密视频| 两个人免费观看高清视频| 久久国产精品大桥未久av| 1024视频免费在线观看| 精品国产一区二区久久| 国产精品 欧美亚洲| 王馨瑶露胸无遮挡在线观看| 免费看av在线观看网站| 免费少妇av软件| 精品人妻1区二区| 日韩电影二区| 香蕉丝袜av| 美女主播在线视频| 久久久久视频综合| 免费在线观看日本一区| 手机成人av网站| 嫩草影视91久久| 国产午夜精品一二区理论片| 久久精品人人爽人人爽视色| 18禁黄网站禁片午夜丰满| 亚洲精品国产色婷婷电影| 亚洲av片天天在线观看| 一边摸一边抽搐一进一出视频| 国产精品久久久久久精品电影小说| 不卡av一区二区三区| 婷婷色综合www| 五月天丁香电影| 18禁观看日本| 岛国毛片在线播放| 国产成人影院久久av| 十分钟在线观看高清视频www| 啦啦啦视频在线资源免费观看| 亚洲欧美一区二区三区国产| 亚洲精品一区蜜桃| 一二三四社区在线视频社区8| 久久免费观看电影| 国产成人欧美在线观看 | 天堂8中文在线网| 一本久久精品| 亚洲国产欧美日韩在线播放| 久久久精品免费免费高清| 亚洲欧美一区二区三区黑人| 国产精品国产三级专区第一集| 中文字幕色久视频| 一本久久精品| 久久久久国产精品人妻一区二区| 曰老女人黄片| 人人妻,人人澡人人爽秒播 | 亚洲国产av影院在线观看| 各种免费的搞黄视频| 一本久久精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品久久成人aⅴ小说| 免费观看人在逋| 一区二区av电影网| 黄片播放在线免费| 久久国产精品影院| 久久久国产欧美日韩av| 国产三级黄色录像| 最新的欧美精品一区二区| 亚洲av男天堂| 一级毛片女人18水好多 | 国产三级黄色录像| 丝袜人妻中文字幕| 可以免费在线观看a视频的电影网站| 人人妻人人添人人爽欧美一区卜| 777久久人妻少妇嫩草av网站| 亚洲成人免费av在线播放| 午夜福利在线免费观看网站| 中文欧美无线码| 亚洲av电影在线观看一区二区三区| 女性被躁到高潮视频| 亚洲伊人色综图| 国产一区二区激情短视频 | 国产有黄有色有爽视频| 精品欧美一区二区三区在线| 另类精品久久| 99国产综合亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 欧美+亚洲+日韩+国产| 欧美xxⅹ黑人| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 乱人伦中国视频| 久久精品人人爽人人爽视色| 1024视频免费在线观看| 一边亲一边摸免费视频| 国产片特级美女逼逼视频| 亚洲国产欧美日韩在线播放| 亚洲,欧美精品.| 汤姆久久久久久久影院中文字幕| 久久久久精品国产欧美久久久 | 少妇人妻久久综合中文| 9热在线视频观看99| 欧美av亚洲av综合av国产av| 国产欧美亚洲国产| 中国美女看黄片| 国产成人av教育| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 天天躁狠狠躁夜夜躁狠狠躁| 五月开心婷婷网| 观看av在线不卡| 高潮久久久久久久久久久不卡| 欧美 亚洲 国产 日韩一| av在线app专区| 性色av乱码一区二区三区2| 国产精品人妻久久久影院| 国产在视频线精品| 精品亚洲成国产av| 国产成人精品在线电影| 香蕉国产在线看| 午夜91福利影院| 性色av一级| 亚洲视频免费观看视频| 老鸭窝网址在线观看| av网站在线播放免费| 久久久久网色| 日韩视频在线欧美| 男女之事视频高清在线观看 | 午夜免费鲁丝| 丰满饥渴人妻一区二区三| 久久av网站| 99精品久久久久人妻精品| 日本欧美国产在线视频| 国产午夜精品一二区理论片| 男人操女人黄网站| 国产免费视频播放在线视频| 激情视频va一区二区三区| 人妻一区二区av| 亚洲熟女精品中文字幕| 国产片内射在线| 亚洲成人免费电影在线观看 | 亚洲综合色网址| 热99国产精品久久久久久7| 人妻 亚洲 视频| 国产精品成人在线| 黄色怎么调成土黄色| 久久精品国产a三级三级三级| 成人亚洲欧美一区二区av| 高清av免费在线| 激情视频va一区二区三区| 视频区欧美日本亚洲| 在线观看免费午夜福利视频| 一区二区三区乱码不卡18| 欧美精品高潮呻吟av久久| 日韩av在线免费看完整版不卡| 国产人伦9x9x在线观看| 日韩制服丝袜自拍偷拍| 亚洲av日韩在线播放| 欧美大码av| 18禁观看日本| 免费高清在线观看日韩| 高清av免费在线| 亚洲精品国产区一区二| 精品久久久久久久毛片微露脸 | 欧美日韩综合久久久久久| 成年女人毛片免费观看观看9 | 久久影院123| 国产精品国产三级国产专区5o| 19禁男女啪啪无遮挡网站| 国产一区亚洲一区在线观看| 男人爽女人下面视频在线观看| 久久亚洲国产成人精品v| 最黄视频免费看| 欧美人与性动交α欧美软件| av视频免费观看在线观看| 制服诱惑二区| 女人爽到高潮嗷嗷叫在线视频| 欧美在线黄色| 久久精品成人免费网站| 日日爽夜夜爽网站| 一个人免费看片子| 婷婷色麻豆天堂久久| 久久青草综合色| 在线观看人妻少妇| 日韩av不卡免费在线播放| 高潮久久久久久久久久久不卡| 国产在线观看jvid| 这个男人来自地球电影免费观看| 国精品久久久久久国模美| 男女午夜视频在线观看| 日韩av在线免费看完整版不卡| 久久精品久久精品一区二区三区| 日本五十路高清| 大陆偷拍与自拍| 国产高清国产精品国产三级| 999久久久国产精品视频| 国产成人精品久久久久久| 脱女人内裤的视频| 国产亚洲av片在线观看秒播厂| 亚洲精品日韩在线中文字幕| 99久久综合免费| 精品国产一区二区三区四区第35| 嫁个100分男人电影在线观看 | 欧美人与善性xxx| 国产亚洲午夜精品一区二区久久| tube8黄色片| 性少妇av在线| 亚洲专区中文字幕在线| 中文字幕高清在线视频| 亚洲国产欧美日韩在线播放| 老司机影院毛片| 男女国产视频网站| 精品人妻1区二区| 久久免费观看电影| 大片免费播放器 马上看| 国产精品国产av在线观看| 国产麻豆69| 亚洲精品成人av观看孕妇| 国产精品.久久久| 国产精品三级大全| 久久久久久久大尺度免费视频| 大话2 男鬼变身卡| 真人做人爱边吃奶动态| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 国产亚洲一区二区精品| 欧美日韩综合久久久久久| 青青草视频在线视频观看| 97在线人人人人妻| 欧美精品一区二区免费开放| 免费在线观看视频国产中文字幕亚洲 | 亚洲午夜精品一区,二区,三区| 国产在线免费精品| 国产人伦9x9x在线观看| 波多野结衣av一区二区av| 又大又黄又爽视频免费| 国产亚洲av高清不卡| 天堂中文最新版在线下载| 国产1区2区3区精品| a级片在线免费高清观看视频| 成人国语在线视频| 黄网站色视频无遮挡免费观看| 国产成人精品久久二区二区91| 国产成人精品久久二区二区免费| 精品人妻熟女毛片av久久网站| 亚洲精品国产一区二区精华液| 欧美+亚洲+日韩+国产| 丝袜人妻中文字幕| 国产亚洲一区二区精品| 免费女性裸体啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 少妇被粗大的猛进出69影院| 国产精品亚洲av一区麻豆| 亚洲国产欧美日韩在线播放| 1024视频免费在线观看| www日本在线高清视频| 亚洲av日韩精品久久久久久密 | av国产久精品久网站免费入址| 久久精品国产亚洲av涩爱| 欧美少妇被猛烈插入视频| 丁香六月天网| 香蕉国产在线看| 丰满饥渴人妻一区二区三| 黄色视频不卡| 国产成人av激情在线播放| 日韩制服丝袜自拍偷拍| 国产一区二区三区av在线| 午夜激情av网站| 天天添夜夜摸| 黄色毛片三级朝国网站| 国产午夜精品一二区理论片| 欧美黄色淫秽网站| 亚洲图色成人| 欧美日韩视频精品一区| 精品少妇一区二区三区视频日本电影| 国产成人精品久久久久久| 久久久精品94久久精品| 亚洲三区欧美一区| 国产精品九九99| 18禁裸乳无遮挡动漫免费视频| 超碰成人久久| 国产91精品成人一区二区三区 | 日韩伦理黄色片| 制服人妻中文乱码| 一边摸一边抽搐一进一出视频| 午夜免费成人在线视频| 日韩,欧美,国产一区二区三区| 久久久亚洲精品成人影院| 国产一区二区在线观看av| 黄色a级毛片大全视频| 亚洲精品日本国产第一区| 妹子高潮喷水视频| 搡老岳熟女国产| 天天操日日干夜夜撸| 男女边吃奶边做爰视频| 亚洲成人免费电影在线观看 | 黄色片一级片一级黄色片| 国语对白做爰xxxⅹ性视频网站| 欧美在线黄色| 久久免费观看电影| 亚洲 欧美一区二区三区| videos熟女内射| 亚洲国产欧美一区二区综合| 天堂8中文在线网| av天堂在线播放| 青春草亚洲视频在线观看| 欧美97在线视频| 日韩大片免费观看网站| 亚洲国产精品国产精品| 欧美av亚洲av综合av国产av| 亚洲欧洲国产日韩| 国产欧美日韩一区二区三区在线| 视频在线观看一区二区三区| 超碰成人久久| 久久人妻熟女aⅴ| 我要看黄色一级片免费的| 夫妻性生交免费视频一级片| 免费一级毛片在线播放高清视频 | 日韩,欧美,国产一区二区三区| 国产高清国产精品国产三级| 国产精品香港三级国产av潘金莲 | 久久久久久人人人人人| 男人操女人黄网站| 亚洲三区欧美一区| 丝袜人妻中文字幕| 一级黄片播放器| 精品福利永久在线观看| 美女国产高潮福利片在线看| 亚洲国产精品一区二区三区在线| 热99国产精品久久久久久7| 51午夜福利影视在线观看| 无限看片的www在线观看| 国产有黄有色有爽视频| 咕卡用的链子| 国产有黄有色有爽视频| 夜夜骑夜夜射夜夜干| 美女脱内裤让男人舔精品视频| 久久精品国产综合久久久| 人人妻人人澡人人爽人人夜夜| 王馨瑶露胸无遮挡在线观看| 成人亚洲欧美一区二区av| 一本久久精品| 叶爱在线成人免费视频播放| 国语对白做爰xxxⅹ性视频网站| 90打野战视频偷拍视频| 欧美乱码精品一区二区三区| 美女高潮到喷水免费观看| 久久久久久久精品精品| 亚洲男人天堂网一区| 两个人看的免费小视频| 亚洲,欧美,日韩| 精品一区二区三区av网在线观看 | 母亲3免费完整高清在线观看| 久久精品亚洲熟妇少妇任你| 国产熟女欧美一区二区| 人妻人人澡人人爽人人| 国产熟女欧美一区二区| 婷婷色综合www| 国产伦人伦偷精品视频| 一本久久精品| 大型av网站在线播放| 精品福利观看| 日韩av不卡免费在线播放| 亚洲欧美日韩另类电影网站| 一区二区三区精品91| 国产亚洲一区二区精品| 成人黄色视频免费在线看| 老汉色∧v一级毛片| 亚洲人成电影免费在线| 韩国精品一区二区三区| 国产97色在线日韩免费| 久久狼人影院| 国产高清国产精品国产三级| 国产精品二区激情视频| 久久 成人 亚洲|