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

    基于混合博弈的產(chǎn)品多學(xué)科柔性設(shè)計決策

    2017-08-31 13:18:35李廖平
    中國機械工程 2017年15期
    關(guān)鍵詞:子系統(tǒng)柔性決策

    陳 亮 李廖平

    福州大學(xué)機械工程及自動化學(xué)院,福州,350116

    基于混合博弈的產(chǎn)品多學(xué)科柔性設(shè)計決策

    陳 亮 李廖平

    福州大學(xué)機械工程及自動化學(xué)院,福州,350116

    針對產(chǎn)品多學(xué)科多目標的設(shè)計決策問題,給出了一種融合博弈論和柔性機制進行建模和求解的方法。通過混合博弈來表征學(xué)科之間及目標之間的耦合關(guān)系,提出了混合博弈模型框架及求解方法。考慮到不確定性和變化的影響,為使產(chǎn)品設(shè)計決策具有柔性,引入性能水平滿意度函數(shù)來柔性表達設(shè)計要求,并尋求設(shè)計變量的柔性范圍解來替代固定值的點解,從而增加產(chǎn)品的柔性;柔性范圍解處于可行設(shè)計空間內(nèi),使各學(xué)科性能處于滿意水平并靠近期望值,且具有較好的穩(wěn)定性。利用概率分布對相應(yīng)的性能變化進行建模,引入設(shè)計偏好指數(shù)對柔性設(shè)計方案進行評價。最后通過減速器算例的設(shè)計來驗證該方法的可行性和有效性。

    多學(xué)科設(shè)計;柔性決策;博弈;偏好指數(shù);減速器

    0 引言

    復(fù)雜產(chǎn)品設(shè)計往往涉及多個學(xué)科領(lǐng)域,需要綜合考慮多學(xué)科間的相互作用和協(xié)同效應(yīng),進行協(xié)同設(shè)計決策來獲得各方滿意的設(shè)計方案。由于多學(xué)科間存在著復(fù)雜的耦合關(guān)系,學(xué)科間相互制約,甚至存在矛盾沖突,存在信息的不確定性,使得設(shè)計方案的確定十分困難,因此,復(fù)雜產(chǎn)品設(shè)計決策中的多學(xué)科耦合分析是設(shè)計的關(guān)鍵,如何有效地處理學(xué)科間的耦合關(guān)系,獲得各學(xué)科均滿意且綜合性能好的設(shè)計方案是復(fù)雜產(chǎn)品設(shè)計開發(fā)的核心問題之一。

    多學(xué)科設(shè)計優(yōu)化(multidisciplinary design optimization ,MDO)是常用的方法之一,多學(xué)科設(shè)計優(yōu)化算法可以分為單級算法和多級算法[1]。單級算法(如單學(xué)科可行法和多學(xué)科可行法)[2]將所有的決策優(yōu)化集中到一個優(yōu)化器中,其他學(xué)科只在需要的時候提供學(xué)科分析,它是一種集中式的決策模型,問題復(fù)雜時難以求解。多級算法如并行子空間優(yōu)化(concurrent subspace optimization,CSSO)算法[3]和協(xié)同優(yōu)化(collaborative optimization, CO)算法[4]克服了單級算法的缺點,將系統(tǒng)的優(yōu)化設(shè)計問題分為兩級:一個系統(tǒng)級和并行的多個學(xué)科級。由于系統(tǒng)層和學(xué)科層之間需要不斷地迭代,當問題維數(shù)高時求解困難,計算復(fù)雜性和組織復(fù)雜性難以消除,收斂性也難以保證。另外,子系統(tǒng)間復(fù)雜的耦合關(guān)系會帶來設(shè)計信息的不確定性,因為各子系統(tǒng)只能控制整個系統(tǒng)的部分設(shè)計變量,它所需的由其他子系統(tǒng)控制的設(shè)計變量就成為其未知和不確定的耦合變量。MDO方法通過一致性約束(如協(xié)同優(yōu)化等)來考慮子系統(tǒng)間的相互耦合影響,但對耦合變量的處理過于復(fù)雜,會消耗大量的分析迭代次數(shù),忽略了學(xué)科設(shè)計間的策略互動與動態(tài)調(diào)整過程[5-6]。為此,復(fù)雜產(chǎn)品的設(shè)計需要對子系統(tǒng)間的交互作用進行建模,并對子系統(tǒng)間的耦合所導(dǎo)致的信息不確定性進行分析和處理,以減少子系統(tǒng)之間的耦合,并尋求提高系統(tǒng)質(zhì)量和穩(wěn)定性的魯棒解,提高系統(tǒng)的整體性能。

    博弈論在表達工程團隊間的交互時是有效的[7-8],可有效應(yīng)用于產(chǎn)品設(shè)計中[9-10]。利用均勻設(shè)計試驗法能確定機械產(chǎn)品多目標優(yōu)化設(shè)計博弈分析中博弈方的策略集[11]。利用博弈理論進行產(chǎn)品設(shè)計可避免人為確定各目標權(quán)重或建立特別綜合評價函數(shù)帶來的影響[12],博弈方遵循博弈過程中的理性法則進行競爭與合作,動態(tài)協(xié)調(diào)最終達到均衡點,克服了MDO方法忽略學(xué)科間策略互動與動態(tài)調(diào)整的不足。另外,通過去除系統(tǒng)層團隊,給各學(xué)科團隊賦予自治權(quán),能夠減少迭代次數(shù),降低計算和組織復(fù)雜性。

    復(fù)雜產(chǎn)品設(shè)計是一個漸進過程,存在很多不確定性因素,為了減少設(shè)計迭代反復(fù),獲得綜合性能好、抵御不確定因素干擾和適應(yīng)變化的設(shè)計解,需要在設(shè)計中融入柔性[13]。柔性是一種適應(yīng)不斷變化的能力[14],為了增加產(chǎn)品的柔性,需要在產(chǎn)品設(shè)計中考慮影響產(chǎn)品設(shè)計的變化因素,如設(shè)計要求的變化、產(chǎn)品設(shè)計參數(shù)的變化等,并采用有效的方法和機制處理這些變化因素,以保證在這些變化情況下仍保持設(shè)計解的可行性,以及性能水平的滿意度和穩(wěn)定性。針對上述問題,本文采用混合博弈對多學(xué)科交互設(shè)計決策過程進行分析,建立多學(xué)科設(shè)計決策的博弈框架,通過引入柔性機制和設(shè)計偏好指數(shù)來考慮不確定性的影響,處理各博弈方的策略選擇和進行各博弈方之間的收益均衡,在此基礎(chǔ)上,給出了一種復(fù)雜產(chǎn)品多學(xué)科柔性設(shè)計決策博弈模型,通過求解該模型獲得綜合性能高且各方滿意的解,并以減速器算例來說明該方法的可行性和有效性。

    1 多學(xué)科設(shè)計問題的博弈分析框架

    復(fù)雜產(chǎn)品設(shè)計涉及多個學(xué)科,而各個學(xué)科通常又有多個目標。本文采用博弈論的方法處理產(chǎn)品多學(xué)科多目標設(shè)計問題,將決策者設(shè)置為局外人,把各個學(xué)科或各個目標作為博弈方。通過計算設(shè)計變量對各個學(xué)科的影響因子,根據(jù)設(shè)計變量與學(xué)科的相關(guān)性進行模糊聚類劃分,將設(shè)計變量分配給與之關(guān)系較為密切的各個學(xué)科和目標。設(shè)計變量方案集合是博弈策略集,一種策略對應(yīng)博弈方的一種收益。通過博弈的理性法則,自然引導(dǎo)各博弈方的競爭與合作,最終達到均衡點,得到對應(yīng)多學(xué)科設(shè)計問題的博弈均衡解。這種基于博弈理性的均衡解具有較好的穩(wěn)定性與合理性,考慮了多個學(xué)科多個目標之間的相互影響和沖突,實現(xiàn)了綜合優(yōu)化和均衡協(xié)調(diào)。決策者這個局外人則根據(jù)自己的偏好,從博弈解中選取其滿意解。根據(jù)學(xué)科間的關(guān)系及彼此信息了解情況,學(xué)科間可以是合作博弈、非合作博弈或領(lǐng)導(dǎo)-隨從博弈。學(xué)科下的優(yōu)化目標隸屬于學(xué)科,與學(xué)科之間的關(guān)系是領(lǐng)導(dǎo)-隨從關(guān)系,用領(lǐng)導(dǎo)-隨從博弈來分析和建模,而學(xué)科下各優(yōu)化目標之間視為合作關(guān)系,采用合作博弈來分析和建模,相應(yīng)的混合博弈分析框架如圖1所示。

    圖1 博弈分析框架Fig.1 Game framework

    1.1 合作博弈(Pareto博弈)

    在合作博弈中,各博弈方彼此間共享信息,以集體理性為基礎(chǔ),以追求集體利益為目標,共同遵守某個具有合作性的約束力協(xié)議,通過協(xié)商和合作來獲取聯(lián)盟的整體最優(yōu)利益。通過合作獲取的聯(lián)盟整體收益將按照一定的規(guī)則進行分配,基本原則是保證每個博弈方所得不少于它們單獨行動(脫離聯(lián)盟,不合作)而獲得的利益,其結(jié)果對各博弈方來說不一定是最優(yōu)結(jié)果,但一定是可以接受的相對較優(yōu)的結(jié)果,即合作博弈結(jié)果是一個Pareto最優(yōu)解。

    定義合作博弈的決策模型為

    (1)

    由于各設(shè)計目標之間的數(shù)量級可能存在差異,為了避免數(shù)量級相差太大而對設(shè)計結(jié)果產(chǎn)生影響,需要對各個設(shè)計目標進行歸一化處理,將各個博弈方置于同一水平,以消除數(shù)量級的影響。歸一化處理后合作博弈的決策模型為

    (2)

    1.2 非合作博弈(Nash博弈)

    非合作博弈發(fā)生于由于組織、信息或過程等障礙而相互隔離的情形,各博弈方之間信息不共享,各自決策,以個體理性為基礎(chǔ),各博弈方選擇一種對自己有益的策略,而不考慮其他博弈方的得失,也不考慮集體利益,以自身利益最大為目標進行博弈,獲得Nash均衡解,在均衡點位置任何博弈方都不愿意改變自己的策略選擇,否則將會降低自身的得益,所以Nash均衡解是一個穩(wěn)定的收斂解,但不一定是 Pareto 最優(yōu)的。

    本文采用Nash遺傳算法進行非合作博弈求解。每個博弈方被分配一個種群和一個目標函數(shù)進行優(yōu)化。每個博弈方只能修改分配給它的設(shè)計變量,通過標準遺傳算子執(zhí)行操作。在每個時間點(即一代)后,博弈方和其他博弈方交換它們的最優(yōu)值,算法最終收斂于Nash均衡解。具體過程如圖2所示,假設(shè)有兩博弈方P1和P2,博弈方P1控制變量x,博弈方P2控制變量y。分別給每個博弈方分配一個種群,開始時,初始化兩種群,通過標準遺傳算子,兩個種群在各自種群內(nèi)優(yōu)化。需要注意的是,種群P1優(yōu)化時,x是變量,y為種群P2上代優(yōu)化后的最優(yōu)解,是一個確定的量;同理,種群P2優(yōu)化時,y是變量,x為種群P1上代優(yōu)化后的最優(yōu)解,是一個確定的量。當優(yōu)化達到終止條件時,達到Nash均衡,否則繼續(xù)優(yōu)化。與并行遺傳算法PGA中每個子種群使用相同的準則不同,Nash遺傳算法則中各子種群使用不同的準則,因此有均衡的概念。

    圖2 Nash遺傳算法Fig.2 Nash genetic algorithm

    1.3 領(lǐng)導(dǎo)-隨從博弈

    領(lǐng)導(dǎo)-隨從博弈是一種考慮了博弈方設(shè)計決策順序的博弈模型,是Nash博弈的一個特例。領(lǐng)導(dǎo)-隨從博弈分領(lǐng)導(dǎo)博弈方和隨從博弈方,兩方的地位不同,領(lǐng)導(dǎo)博弈方在博弈中取得較高的決策權(quán)。該博弈是由領(lǐng)導(dǎo)博弈方在假定隨從博弈方理性行為的情況下做出自己的設(shè)計決策,隨后隨從博弈方再根據(jù)領(lǐng)導(dǎo)博弈方的決策做出自己的決策,并且隨從博弈方的決策要與領(lǐng)導(dǎo)博弈方的決策協(xié)調(diào)一致,盡量靠近領(lǐng)導(dǎo)博弈方的決策。第i個隨從博弈方的博弈模型為

    (3)

    1.4 混合博弈

    混合博弈是指在一個博弈過程中同時擁有以上三種博弈中的兩種及兩種以上的博弈形式,圖1所示的多學(xué)科博弈框架即混合博弈,其中涉及合作博弈、非合作博弈和領(lǐng)導(dǎo)-隨從博弈三種博弈方式,其求解過程如圖3所示。

    圖3 混合博弈求解算法Fig.3 Mixed game algorithm

    (1)各學(xué)科之間先進行非合作博弈(或合作博弈、領(lǐng)導(dǎo)-隨從博弈,根據(jù)學(xué)科間的關(guān)系來確定),得到Nash均衡解(或Pareto解)。

    (2)各學(xué)科作為領(lǐng)導(dǎo)方,其下各子目標為隨從方,進行領(lǐng)導(dǎo)-隨從博弈。領(lǐng)導(dǎo)方的解為步驟(1)中的Nash均衡解(或Pareto解),由領(lǐng)導(dǎo)方的解可求得各隨從方的期望值,然后隨從方進行設(shè)計決策獲得其隨從解。

    (3)各子目標間進行合作博弈得到Pareto前沿,再從中選取與步驟(1)中的解最接近的一個解。

    2 產(chǎn)品柔性設(shè)計決策

    為了適應(yīng)不確定性和變化的條件,產(chǎn)品設(shè)計決策需要具有柔性??紤]影響產(chǎn)品設(shè)計要求的變化,應(yīng)尋求設(shè)計要求的柔性表述,替代傳統(tǒng)的設(shè)計要求被指定為剛性約束或目標(表現(xiàn)為一個精確的數(shù)值)的情形,允許設(shè)計性能要求是變化值,不同的性能水平滿意度不同,通過引入性能水平的滿意度函數(shù)S(y)來衡量不同性能水平的滿意度,其定義為滿意度S與性能水平y(tǒng)之間關(guān)系的函數(shù)??紤]到性能水平離期望值越遠,設(shè)計者對目標的滿意度呈加速下降的趨勢,本文使用二次函數(shù)來替代傳統(tǒng)的剛性精確值(如圖4為小好型(STB),根據(jù)設(shè)計問題類型,也可以是大好型(LTB)或中間好型(CTB))。函數(shù)值在0~1之間,1代表完全滿意,表示設(shè)計方案被認為是完全可接受或滿意的;函數(shù)值為0代表完全不滿意,表示設(shè)計方案被認為是不可接受的,這種方法與傳統(tǒng)優(yōu)化方法相比更加具有柔性。

    圖4 柔性設(shè)計決策Fig.4 Flexible design and decision

    為了對設(shè)計進行評價,引入設(shè)計偏好指數(shù)(design preference index ,DPI)IDP[15],其定義為設(shè)計變化范圍內(nèi)性能滿意度函數(shù)S(y)的數(shù)學(xué)期望:

    (4)

    當設(shè)計變量在設(shè)計解范圍內(nèi)變化時,引起設(shè)計性能和性能滿意度在相應(yīng)范圍內(nèi)變化,IDP的含義是衡量和評價在這種變化下對設(shè)計要求的總體滿足情況,即對本文給出的柔性設(shè)計方案進行衡量和評價的一種指標。IDP值越大,表示性能分布函數(shù)和滿意度函數(shù)的重合度越大。IDP為1時意味著在這種變化下設(shè)計能完全滿足設(shè)計要求,即在設(shè)計解范圍內(nèi)的點均是可行的,且相應(yīng)的性能水平具有高的滿意度,這是理想情況,顯然,使IDP盡量接近1是所期望的。

    (5)

    然而,IDP值為1在某些情況下是不能達到的,如設(shè)計方案②、③所示。

    基于上述分析,將設(shè)計決策模型中的目標函數(shù)定義為設(shè)計能力指數(shù)IDP值與1之差,給設(shè)計變量xi以變化范圍Δxi,并增加約束Δxi≥Δximin,這是為了確保設(shè)計解的最小柔性,若沒有這個約束,設(shè)計變量Δxi范圍將會趨近于0,設(shè)計解退化為常規(guī)的單點解,性能水平則成為確定值。增加約束Δyi≤yimax,可限定性能的最大變化,使性能保持必要的穩(wěn)定性。由此得到柔性設(shè)計決策模型:

    設(shè)計變量:x,Δxi

    設(shè)計目標:min 1-IDPi;i=1,2,…,m

    設(shè)計約束:Δxi≥Δximin,i=1,2,…,n

    Δyi≤Δyimax,i=1,2,…,m

    邊界范圍:

    lbj≤xj≤ubj

    通過求解該模型,可獲得設(shè)計變量柔性范圍解,它處于可行設(shè)計空間內(nèi),且在此范圍解集內(nèi),各學(xué)科性能處于滿意水平并盡可能靠近期望值,而且具有較好的穩(wěn)定性。

    3 實例分析

    齒輪減速器問題[16](speed reducer problem,SRP)被NASA認為是檢驗MDO方法性能的十大標準算例之一。該問題包含齒輪設(shè)計和軸設(shè)計兩部分,目標是在滿足齒輪的彎曲應(yīng)力、接觸應(yīng)力以及軸的位移和應(yīng)力等約束條件下使減速器的質(zhì)量最小。其數(shù)學(xué)表達如下:

    g7=x2x3/40-1.0≤0

    g8=5x2/x1-1.0≤0

    g9=x1/(12x2)-1.0≤0

    g10=(1.5x6+1.9)/x4-1.0≤0

    g11=(1.1x7+1.9)/x5-1.0≤0

    2.6≤x1≤3.6 0.7≤x2≤0.8

    17≤x3≤28 7.3≤x4≤8.3

    7.3≤x5≤8.3 2.9≤x6≤3.9 5.0≤x7≤5.5

    式中,x1為齒寬;x2為齒輪模數(shù);x3為小齒輪齒數(shù);x4為軸1上軸承間距;x5為軸2上軸承間距;x6為軸1直徑;x7為軸2直徑;g1為齒輪輪齒最大彎曲應(yīng)力約束;g2為齒輪輪齒最大接觸應(yīng)力約束;g3為軸1橫向最大撓度約束;g4為軸2橫向最大撓度約束;g5為軸1最大應(yīng)力約束;g6為軸2最大應(yīng)力約束;g7、g8、g9為尺寸和空間約束;g10為軸1尺寸經(jīng)驗設(shè)計約束;g11為軸2尺寸經(jīng)驗設(shè)計約束。

    在MDO中,將學(xué)科定義為系統(tǒng)中本身相對獨立、相互之間又有數(shù)據(jù)交換關(guān)系的基本模塊,又稱作子系統(tǒng)或子空間[17]。將減速器設(shè)計問題分解為兩個子系統(tǒng)和三個目標的設(shè)計優(yōu)化問題,子系統(tǒng)1負責齒輪設(shè)計,子系統(tǒng)2負責軸設(shè)計,見表1。

    3.1 單目標最優(yōu)解

    對各目標進行優(yōu)化得到相應(yīng)的單目標最優(yōu)解:

    表1 學(xué)科子系統(tǒng)劃分

    3.2 策略集劃分

    表2 各目標偏導(dǎo)數(shù)情況

    可得出各變量對三個目標的影響因子:Δ1={1,0,0};Δ2={1,0,0};Δ3={1,0,0};Δ4={0,1,0};Δ5={0,0,1};Δ6={0,1,0};Δ7={0,0,1}。然后通過計算相似度,得到模糊相似矩陣;再利用傳遞閉包法求得具有傳遞性的模糊關(guān)系矩陣;最后根據(jù)實際情況,選取合適閾值,得到相應(yīng)分類,確定各博弈方的策略集,限于篇幅,不再贅述,可參見相應(yīng)模糊聚類分析的文獻,最終獲得的策略集劃分結(jié)果為:{x1,x2,x3}隸屬于y1;{x4,x6}隸屬于y2;{x5,x7}隸屬于y3。

    3.3 多學(xué)科混合博弈模型框架的建立

    建立圖5所示的減速器混合博弈模型框架,齒輪設(shè)計子系統(tǒng)和軸設(shè)計子系統(tǒng)間進行Nash博弈, 軸設(shè)計子系統(tǒng)下的軸1設(shè)計和軸2設(shè)計子目標間進行合作博弈。為了便于比較,分為不使用設(shè)計偏好指數(shù)和使用設(shè)計偏好指數(shù)兩種情況進行分析。

    圖5 減速器博弈模型框架Fig.5 Game model of speed reducer

    3.3.1 不使用設(shè)計偏好指數(shù)

    (3)合作博弈。建立合作博弈方1(軸1)模型和合作博弈方2(軸2)模型,軸1、軸2設(shè)計處于同一數(shù)量級水平,進行歸一化處理,確定合作博弈決策模型。

    (4)求解結(jié)果?;谇笆龅幕旌喜┺那蠼馑惴?,設(shè)定ε=0.01,N=100進行求解,迭代56次得到混合博弈結(jié)果:

    x*=

    {3.5047,0.7,17.0174,7.3016,7.7558,3.3526,5.2872}

    y1=1380.8y2=346.213y3=1275.4

    3.3.2 使用設(shè)計偏好指數(shù)

    (1)滿意度函數(shù)的設(shè)定。各個學(xué)科子系統(tǒng)滿意度函數(shù)是根據(jù)設(shè)計者的經(jīng)驗和偏好來確定的,設(shè)定減速器各學(xué)科子系統(tǒng)的滿意度函數(shù)如圖6所示。

    圖6 滿意度函數(shù)Fig.6 Satisfaction function

    i=1,2,3

    S(y1)=

    (3)非合作博弈模型。

    Nash博弈方1模型:

    Δxi≥0.1xiΔyi≤0.2yi

    Nash博弈方2模型:

    (4)領(lǐng)導(dǎo)隨從博弈模型。

    領(lǐng)導(dǎo)博弈方模型:

    隨從博弈方1模型:

    隨從博弈方2模型:

    (5)合作博弈模型。

    合作博弈方1模型:

    合作博弈方2模型:

    合作博弈決策模型:

    s.t.gi(x)≤0i=1,2,…,11

    (6)求解。設(shè)定ε=0.01,N=100進行求解,迭代78次后得到混合博弈柔性結(jié)果:

    x={3.5208±0.3291,0.7024±0.0722,17.014±2.7187,

    7.4671±0.7774,7.7491±0.8267,3.3495±0.3332,

    5.2890±0.5409}

    y1=1396.2±594.58y2= 346.77±97.184

    y3=1276.5±374.66

    表3 結(jié)果對比

    4 結(jié)語

    本文基于博弈論對多學(xué)科多目標交互設(shè)計決策過程進行分析,通過混合博弈來建立學(xué)科及目標之間的耦合關(guān)系模型, 提出了混合博弈模型框架及求解方法。為了使設(shè)計具有柔性以適應(yīng)不確定性和變化,采用二次型滿意度函數(shù)來柔性表達設(shè)計要求,另外以設(shè)計變量的范圍解來替代點解,從而使得設(shè)計解具有適應(yīng)變化的柔性。用概率分布密度函數(shù)反映性能變化,并用設(shè)計偏好指數(shù)對設(shè)計方案進行評價。實例運行結(jié)果表明,本文方法是可行和有效的,柔性范圍解處于可行設(shè)計空間內(nèi),各學(xué)科目標性能處于滿意水平并靠近期望值,而且具有較好的穩(wěn)定性。對于復(fù)雜產(chǎn)品多學(xué)科多目標設(shè)計決策問題,得到綜合性能好且能適應(yīng)不確定性和變化因素影響的范圍解比追求單純的最優(yōu)解更有實際意義,所以本文方法具有較好的指導(dǎo)意義,并且獲得的柔性范圍解可作為進一步詳細設(shè)計的基礎(chǔ)。

    [1] 劉克龍,姚衛(wèi)星,余雄慶.幾種新型多學(xué)科設(shè)計優(yōu)化算法及比較[J].計算機集成制造系統(tǒng),2007,13(2):209-216. LIU Kelong , YAO Weixing, YU Xiongqing. Some New Algorithms for Multidisciplinary Design Optimization and Their Comparison [J].Computer Integrated Manufacture Systems,2007,13(2):209-216.

    [2] HAFTKA R T , SOBIESZCZANSKI-SOBIESKI J, PADULA S L . On Options for Interdisciplinary Analysis and Design Optimization[J]. Structural and Multidisciplinary Optimization,1992,4(2):65-74.

    [3] SELLAR R S, BATILL S M, RENAUD J E.Response Surface Based,Concurrent Subspace Optimization for Multidisciplinary System Design[C]//34th AIAA Aerospace Sciences Meeting and Exhibit.Reno, 1996,AIAA-1996-0714.

    [4] SOBIESKI I P, KROO I M. Collaborative Optimization Using Response Surface Estimation[J]. AIAA Journal, 2000,38(10):1931-1938.

    [5] 馬明旭,王承恩,張嘉易,等. 復(fù)雜產(chǎn)品多學(xué)科設(shè)計優(yōu)化技術(shù)[J].機械工程學(xué)報,2008,44(6):15-26. MA Mingxu, WANG Chengen, ZHANG Jiayi,et al.Multidisciplinary Design Optimization for Complex Product Review[J].Journal of Mechanical Engineering,2008,44(6):15-26.

    [6] 葉國青,姜江,陳森,等.武器裝備體系設(shè)計問題求解框架與優(yōu)化方法[J].系統(tǒng)工程與電子技術(shù), 2012,34(11):2256-2263. YE Guoqing,JIANG Jiang,CHEN Sen,et al.Design Optimization of Weapon System of Systems Based on Collaborative Multidisciplinary Decision-making[J].System Engineer and Electronics,2012,34(11):2256-2263.

    [7] LEWIS K, MISTREE F. Modeling Interactions in Multidisciplinary Design: a Game Theoretic Approach[J]. AIAA Journal,1997,35(8):1387-1392.

    [8] Shun T K. A Game-theoretic Model of Collaboration in Engineering Design[J].Journal of Mechanical Design, 2010,132(5):051005-1-9.

    [9] DHINGRA A K, RAO S S. A Cooperative Fuzzy Game Theoretic Approach to Multiple Objective Design Optimization[J].European Journal of Operational Research,1995,83(3):547-567.

    [10] LIANG Z X, YAN L, SHANG J Z.Ship Cabin Layout Design Using Game Theory[J]. Journal of Marine Science and Technology, 2008,13(4):446-454.

    [11] 陳加明,宋小文,徐浩.結(jié)合模糊聚類和合作競爭博弈的優(yōu)化方法[J].中國機械工程,2015,26(16):2179-2182. CHEN Jiaming, SONG Xiaowen, XU Hao.Optimization Method Based on Fuzzy Clustering and Cooperative Competition Game[J]. China Mechanical Engineering,2015,26(16):2179-2182.

    [12] 鄭丞,金隼,來新民,等.基于非合作博弈的公差分配優(yōu)化[J].機械工程學(xué)報,2009,45(10):159-165. ZHENG Cheng, JIN Sun, LAI Xinmin,et al.Tolerance Allocation Optimization Based on Non-cooperative Game Analysis[J].Journal of Mechanical Engineering,2009,45(10):159-165.

    [13] SIMPSON T W, ROSEN D, ALLEN J K, et al. Metrics for Assessing Design Freedom and Information Certainty in the Early Stages of Design[J]. Journal of Mechanical Design, 1998,120(4):628-635.

    [14] JOSEPH H S, DANIEL E H, DAVA J N. Flexibility in System Design and Implications for Aerospace Systems[J]. Acta Astronautica,2003,53(12):927-944.

    [15] CHEN W, YUAN C. A Probabilistic-based Design Model for Achieving Flexibility in Design[J]. Journal of Mechanical Design,1999,121(1):77-83.

    [16] PADULA S L, ALEXANDROV N, GREEN L L. MDO Test Suite at NASA Langley Research Center[C]//Proceedings of the 6th NASA/ISSMO Symposium on Multidisciplinary Analysis and Optimization.Bellevue, 1996:AIAA-1996-4028.

    [17] 王振國,陳小前,羅文彩,等.飛行器多學(xué)科設(shè)計優(yōu)化理論與應(yīng)用研究[M].北京:國防工業(yè)出版社,2006. WANG Zhenguo, CHEN Xiaoqian, LUO Wencai,et al.Research on the Theory and Application of Multidisciplinary Design Optimization of Flight Vehicles[M].Beijing: National Defend Industry Press,2006.

    [18] LI X, LIU C, LI W, et al. An Alternative Formulation of Collaborative Optimization Based on Geometric Analysis[J].Journal of Mechanical Design,2011,133(5):623-635.

    (編輯 陳 勇)

    Multidisciplinary Flexible Design and Decision Based on Mixed Game

    CHEN Liang LI Liaoping

    College of Mechanical Engineering and Automation, Fuzhou University, Fuzhou ,350116

    For the problems of multidisciplinary and multiobjective design and decision, a novel method was given using game theory and flexibility mechanism to model and solve the problems. The coupling relations among disciplines or targets were presented by mixed game, and the mixed game model framework and the corresponding solving methods were proposed. Considering the influences of uncertainty and changes, in order to make the design and decision be flexible, the satisfaction function of performance was introduced to express the flexible design requirements, and the flexible range solutions of design variables were sought rather than a point solution to increase product flexibility. The flexible range solutions were in the feasible design spaces, where every discipline’s performances were satisfactory and close to the expected values and had good stability. Probabilistic representation was used to model the corresponding design performance changes, and design preference index (DPI) was introduced to evaluate the flexible design. Finally, a gear reducer design example was given to verify the feasibility and effectiveness of the proposed method.

    multidisciplinary design; flexibility decision; game; preference index; reducer

    2016-09-01

    國家自然科學(xué)基金資助項目(50875049);福建省自然科學(xué)基金資助項目(2014J01184)

    TH122

    10.3969/j.issn.1004-132X.2017.15.014

    陳 亮,男,1963年生。福州大學(xué)機械工程及自動化學(xué)院教授、博士。主要研究方向為多學(xué)科協(xié)同設(shè)計、優(yōu)化和決策,現(xiàn)代設(shè)計理論和方法等。E-mail: chenliang@fzu.edu.cn。李廖平,男,1989年生。福州大學(xué)機械工程及自動化學(xué)院碩士研究生。

    猜你喜歡
    子系統(tǒng)柔性決策
    不對中轉(zhuǎn)子系統(tǒng)耦合動力學(xué)特性研究
    一種柔性拋光打磨頭設(shè)計
    灌注式半柔性路面研究進展(1)——半柔性混合料組成設(shè)計
    石油瀝青(2021年5期)2021-12-02 03:21:18
    為可持續(xù)決策提供依據(jù)
    高校學(xué)生管理工作中柔性管理模式應(yīng)用探索
    GSM-R基站子系統(tǒng)同步方案研究
    決策為什么失誤了
    駝峰測長設(shè)備在線監(jiān)測子系統(tǒng)的設(shè)計與應(yīng)用
    車載ATP子系統(tǒng)緊急制動限制速度計算
    柔性的思維
    棋藝(2001年7期)2001-06-14 03:21:00
    熟妇人妻久久中文字幕3abv| 亚洲七黄色美女视频| 国产精华一区二区三区| 中文在线观看免费www的网站| 久久人妻av系列| 国产欧美日韩精品一区二区| 国产探花极品一区二区| 欧美日韩黄片免| 欧美午夜高清在线| 一本久久中文字幕| 国产aⅴ精品一区二区三区波| 国产精品国产高清国产av| 亚洲精品一区av在线观看| 国产一区二区在线av高清观看| 久久性视频一级片| 欧美国产日韩亚洲一区| 亚洲精品一卡2卡三卡4卡5卡| 美女免费视频网站| 国产精品久久久久久人妻精品电影| 国产亚洲欧美98| 少妇的逼水好多| 999久久久精品免费观看国产| 亚洲激情在线av| 精品久久久久久久末码| 久久久久久久亚洲中文字幕 | 国产99白浆流出| 国产精品av视频在线免费观看| 天美传媒精品一区二区| 99久国产av精品| 国产精品亚洲美女久久久| 波多野结衣高清无吗| 亚洲熟妇熟女久久| 51国产日韩欧美| 脱女人内裤的视频| 欧美+亚洲+日韩+国产| 女生性感内裤真人,穿戴方法视频| 久久久色成人| 午夜视频国产福利| 午夜福利18| 久久久精品欧美日韩精品| 在线免费观看不下载黄p国产 | 精品久久久久久久久久久久久| 又黄又爽又免费观看的视频| 日韩欧美一区二区三区在线观看| 午夜免费激情av| 尤物成人国产欧美一区二区三区| 久久精品91蜜桃| 国产精品,欧美在线| 欧美成人一区二区免费高清观看| 国产高清三级在线| 老汉色∧v一级毛片| 国产成年人精品一区二区| 久久久久九九精品影院| 脱女人内裤的视频| 此物有八面人人有两片| 18禁美女被吸乳视频| 最新在线观看一区二区三区| 少妇熟女aⅴ在线视频| eeuss影院久久| 69人妻影院| 国产v大片淫在线免费观看| 日韩国内少妇激情av| 亚洲精品成人久久久久久| 日本a在线网址| 国产v大片淫在线免费观看| 成人高潮视频无遮挡免费网站| 精品久久久久久,| 性色av乱码一区二区三区2| 亚洲精品456在线播放app | netflix在线观看网站| 日韩精品青青久久久久久| 成年女人永久免费观看视频| 亚洲第一欧美日韩一区二区三区| 欧美不卡视频在线免费观看| 亚洲av美国av| 午夜激情欧美在线| 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 亚洲av免费高清在线观看| 中文字幕久久专区| 好男人电影高清在线观看| 淫秽高清视频在线观看| 色综合欧美亚洲国产小说| 国产毛片a区久久久久| av欧美777| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 精品人妻1区二区| 18禁在线播放成人免费| 亚洲国产中文字幕在线视频| 老汉色∧v一级毛片| 久久精品综合一区二区三区| 亚洲精品456在线播放app | 在线十欧美十亚洲十日本专区| 亚洲精品日韩av片在线观看 | 中文资源天堂在线| 国产不卡一卡二| 欧美黑人巨大hd| 黄色丝袜av网址大全| 国产亚洲精品久久久久久毛片| 黄片小视频在线播放| 亚洲最大成人手机在线| 欧美乱码精品一区二区三区| 久久亚洲精品不卡| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 少妇丰满av| 久久天躁狠狠躁夜夜2o2o| 国产又黄又爽又无遮挡在线| 日韩精品青青久久久久久| 欧美一区二区亚洲| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文字幕日韩| 亚洲精品国产精品久久久不卡| 在线免费观看不下载黄p国产 | 欧美高清成人免费视频www| 亚洲成人精品中文字幕电影| 中国美女看黄片| 午夜久久久久精精品| 国产精品久久久久久久电影 | 亚洲欧美日韩东京热| 午夜a级毛片| 日本一二三区视频观看| 少妇的逼水好多| 日韩成人在线观看一区二区三区| 高清毛片免费观看视频网站| 日韩欧美在线乱码| 香蕉av资源在线| 欧美色视频一区免费| 色综合婷婷激情| 欧美av亚洲av综合av国产av| 日本一二三区视频观看| 中文字幕高清在线视频| avwww免费| 一边摸一边抽搐一进一小说| 中文字幕久久专区| 国产一区在线观看成人免费| 一本一本综合久久| 久久精品国产综合久久久| 中文字幕人成人乱码亚洲影| 久9热在线精品视频| 99久久99久久久精品蜜桃| 3wmmmm亚洲av在线观看| 精品欧美国产一区二区三| 小蜜桃在线观看免费完整版高清| 美女免费视频网站| 一级作爱视频免费观看| 成人国产综合亚洲| 亚洲精品美女久久久久99蜜臀| 三级毛片av免费| 伊人久久精品亚洲午夜| 亚洲熟妇熟女久久| 国产成人啪精品午夜网站| 午夜福利在线观看免费完整高清在 | 琪琪午夜伦伦电影理论片6080| 国产精品亚洲av一区麻豆| 亚洲专区国产一区二区| a级毛片a级免费在线| 久久精品国产99精品国产亚洲性色| 免费无遮挡裸体视频| 亚洲av免费在线观看| 无人区码免费观看不卡| 一个人看的www免费观看视频| 欧美黄色片欧美黄色片| 少妇的丰满在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品影院6| 国产亚洲欧美98| 欧美日韩综合久久久久久 | 日本五十路高清| 99国产精品一区二区蜜桃av| 中文字幕高清在线视频| 国产精品,欧美在线| 99热6这里只有精品| 一个人看的www免费观看视频| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看| 12—13女人毛片做爰片一| 日韩欧美免费精品| 国产精品电影一区二区三区| 啦啦啦免费观看视频1| 免费观看的影片在线观看| 中文字幕av在线有码专区| eeuss影院久久| 午夜精品在线福利| 激情在线观看视频在线高清| 天天添夜夜摸| 99热这里只有是精品50| 色噜噜av男人的天堂激情| 麻豆国产av国片精品| 欧美绝顶高潮抽搐喷水| 久久精品影院6| 99久国产av精品| 三级国产精品欧美在线观看| 免费人成视频x8x8入口观看| 中文字幕av在线有码专区| 久久婷婷人人爽人人干人人爱| 精品熟女少妇八av免费久了| 三级男女做爰猛烈吃奶摸视频| 亚洲精品国产精品久久久不卡| 成熟少妇高潮喷水视频| 嫩草影视91久久| 99久久无色码亚洲精品果冻| 欧美日韩国产亚洲二区| 老汉色av国产亚洲站长工具| 成人18禁在线播放| 国产午夜精品久久久久久一区二区三区 | 色尼玛亚洲综合影院| 亚洲 国产 在线| 日本 欧美在线| 午夜视频国产福利| 亚洲国产欧洲综合997久久,| 午夜福利成人在线免费观看| 在线观看av片永久免费下载| 久久天躁狠狠躁夜夜2o2o| 在线国产一区二区在线| 天堂动漫精品| 亚洲第一电影网av| 国产成人系列免费观看| 免费在线观看亚洲国产| 身体一侧抽搐| 最近最新中文字幕大全免费视频| a级一级毛片免费在线观看| 色噜噜av男人的天堂激情| 欧美日韩一级在线毛片| 动漫黄色视频在线观看| 国产精品久久视频播放| 欧美乱码精品一区二区三区| 亚洲国产欧美人成| 九色国产91popny在线| 午夜福利18| 日韩亚洲欧美综合| 在线观看66精品国产| 乱人视频在线观看| 麻豆一二三区av精品| 丰满的人妻完整版| 婷婷丁香在线五月| 日韩精品青青久久久久久| 性色avwww在线观看| 亚洲一区二区三区不卡视频| 日本一二三区视频观看| av天堂中文字幕网| 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 国内揄拍国产精品人妻在线| 一进一出抽搐动态| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| 一区二区三区国产精品乱码| 亚洲熟妇中文字幕五十中出| 国产 一区 欧美 日韩| 美女大奶头视频| 成人av在线播放网站| 变态另类成人亚洲欧美熟女| 免费观看人在逋| 成人欧美大片| 精品国产美女av久久久久小说| 日本黄色片子视频| 亚洲av成人不卡在线观看播放网| 国产欧美日韩一区二区精品| 一区二区三区激情视频| 一边摸一边抽搐一进一小说| 婷婷亚洲欧美| 亚洲av日韩精品久久久久久密| 午夜久久久久精精品| 91av网一区二区| 国产成人av教育| 欧美三级亚洲精品| 日韩欧美精品免费久久 | 中文字幕人成人乱码亚洲影| 国产单亲对白刺激| 国产av在哪里看| 精品一区二区三区视频在线 | 国产精品美女特级片免费视频播放器| 精品不卡国产一区二区三区| 欧美黑人巨大hd| 久久精品亚洲精品国产色婷小说| 九九热线精品视视频播放| 国产高清三级在线| 91麻豆精品激情在线观看国产| 色综合欧美亚洲国产小说| 午夜福利在线在线| bbb黄色大片| 成年女人看的毛片在线观看| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av涩爱 | 亚洲人成伊人成综合网2020| 91久久精品电影网| 一本一本综合久久| 国产极品精品免费视频能看的| 欧美在线一区亚洲| 成人高潮视频无遮挡免费网站| 国产一区二区三区在线臀色熟女| 看免费av毛片| 亚洲熟妇中文字幕五十中出| 高潮久久久久久久久久久不卡| 国产免费av片在线观看野外av| 成年女人毛片免费观看观看9| 在线观看av片永久免费下载| 可以在线观看的亚洲视频| 成年女人看的毛片在线观看| 青草久久国产| 久久香蕉精品热| 亚洲av熟女| 美女被艹到高潮喷水动态| 草草在线视频免费看| 99久国产av精品| 亚洲精品一区av在线观看| 在线观看免费午夜福利视频| 日韩欧美 国产精品| 高清毛片免费观看视频网站| 亚洲成人免费电影在线观看| 十八禁人妻一区二区| 国产精品99久久99久久久不卡| 51午夜福利影视在线观看| 午夜日韩欧美国产| 噜噜噜噜噜久久久久久91| 国产欧美日韩精品一区二区| 女生性感内裤真人,穿戴方法视频| 黑人欧美特级aaaaaa片| 中文字幕人妻丝袜一区二区| 中文字幕久久专区| 有码 亚洲区| 夜夜躁狠狠躁天天躁| 免费av毛片视频| 亚洲男人的天堂狠狠| 亚洲国产精品999在线| xxxwww97欧美| 久久人人精品亚洲av| 国内精品久久久久久久电影| 人人妻,人人澡人人爽秒播| av女优亚洲男人天堂| 亚洲电影在线观看av| 久久精品91无色码中文字幕| 亚洲色图av天堂| 岛国在线免费视频观看| 老熟妇乱子伦视频在线观看| 亚洲精华国产精华精| 老司机在亚洲福利影院| 国产精品 国内视频| www日本在线高清视频| 免费看美女性在线毛片视频| 亚洲一区二区三区不卡视频| 嫩草影院精品99| 18+在线观看网站| 国产精品日韩av在线免费观看| 成人一区二区视频在线观看| 欧美一区二区国产精品久久精品| 免费av观看视频| 亚洲av日韩精品久久久久久密| 国产精品野战在线观看| 国产三级黄色录像| 少妇的逼水好多| 午夜福利免费观看在线| 韩国av一区二区三区四区| 久久久国产精品麻豆| 天堂av国产一区二区熟女人妻| 国产欧美日韩精品亚洲av| 国产成人av激情在线播放| 99精品欧美一区二区三区四区| 中文在线观看免费www的网站| 欧美午夜高清在线| 精品一区二区三区av网在线观看| 一二三四社区在线视频社区8| 国产成人av激情在线播放| 男女床上黄色一级片免费看| 99精品在免费线老司机午夜| 19禁男女啪啪无遮挡网站| 可以在线观看毛片的网站| 天堂√8在线中文| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区久久| 国内精品久久久久精免费| 中文字幕人妻丝袜一区二区| 日日摸夜夜添夜夜添小说| 床上黄色一级片| 欧美最黄视频在线播放免费| 综合色av麻豆| 一本综合久久免费| 18禁美女被吸乳视频| 欧美精品啪啪一区二区三区| 欧美黑人巨大hd| 日本成人三级电影网站| 亚洲精品久久国产高清桃花| 中文亚洲av片在线观看爽| 久久久久九九精品影院| 成人av一区二区三区在线看| 少妇人妻精品综合一区二区 | www日本黄色视频网| 色av中文字幕| 亚洲激情在线av| 男女那种视频在线观看| 老汉色av国产亚洲站长工具| 色综合婷婷激情| 久久久精品欧美日韩精品| 亚洲一区二区三区不卡视频| 丝袜美腿在线中文| 午夜两性在线视频| 精品欧美国产一区二区三| 久久久久久人人人人人| 欧美色欧美亚洲另类二区| 日韩欧美精品v在线| 99久国产av精品| 日本 av在线| 国产精品永久免费网站| 国产精品一区二区三区四区久久| 99久久精品一区二区三区| 国产精品爽爽va在线观看网站| 在线观看舔阴道视频| 超碰av人人做人人爽久久 | 舔av片在线| 亚洲无线观看免费| 99国产精品一区二区蜜桃av| 九九在线视频观看精品| 欧美日韩一级在线毛片| 亚洲乱码一区二区免费版| 亚洲欧美精品综合久久99| 国产高清激情床上av| 黄色丝袜av网址大全| 岛国视频午夜一区免费看| 99热精品在线国产| 深爱激情五月婷婷| 欧美激情久久久久久爽电影| 午夜免费激情av| 老汉色∧v一级毛片| 一进一出抽搐动态| 97人妻精品一区二区三区麻豆| 欧美另类亚洲清纯唯美| 午夜精品久久久久久毛片777| 伊人久久大香线蕉亚洲五| 色视频www国产| 午夜福利欧美成人| 国产 一区 欧美 日韩| 少妇的逼好多水| 99久久99久久久精品蜜桃| 51午夜福利影视在线观看| 国产精品99久久99久久久不卡| 男人的好看免费观看在线视频| 日本黄大片高清| 狂野欧美白嫩少妇大欣赏| 欧美又色又爽又黄视频| 欧美另类亚洲清纯唯美| 国产成年人精品一区二区| 久久久久久久久中文| 噜噜噜噜噜久久久久久91| 亚洲国产日韩欧美精品在线观看 | 国产精品亚洲一级av第二区| 操出白浆在线播放| 少妇人妻一区二区三区视频| 精品人妻一区二区三区麻豆 | 国产在视频线在精品| 88av欧美| 91在线观看av| 欧美成人a在线观看| 亚洲国产高清在线一区二区三| 午夜精品久久久久久毛片777| 女人高潮潮喷娇喘18禁视频| 欧美另类亚洲清纯唯美| 亚洲国产欧洲综合997久久,| 国产亚洲精品一区二区www| 两个人视频免费观看高清| 在线观看66精品国产| 麻豆成人午夜福利视频| 亚洲美女黄片视频| 级片在线观看| 久久久国产精品麻豆| 国产精品乱码一区二三区的特点| 日韩高清综合在线| 大型黄色视频在线免费观看| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区久久| 亚洲无线观看免费| 亚洲精品久久国产高清桃花| 又粗又爽又猛毛片免费看| 精品午夜福利视频在线观看一区| 日韩欧美三级三区| 亚洲成av人片在线播放无| 免费人成视频x8x8入口观看| 悠悠久久av| 又爽又黄无遮挡网站| 久久婷婷人人爽人人干人人爱| 免费观看精品视频网站| 国产伦人伦偷精品视频| 久久精品91蜜桃| 久久久久性生活片| 国产高清激情床上av| 看免费av毛片| 91在线精品国自产拍蜜月 | 色视频www国产| 久久精品国产亚洲av香蕉五月| 99热精品在线国产| 亚洲七黄色美女视频| 亚洲 欧美 日韩 在线 免费| 免费在线观看成人毛片| 亚洲片人在线观看| 免费观看精品视频网站| 91在线精品国自产拍蜜月 | 欧美3d第一页| www.色视频.com| 真人做人爱边吃奶动态| 免费搜索国产男女视频| 一区二区三区高清视频在线| 国产单亲对白刺激| 深夜精品福利| 久久精品国产亚洲av涩爱 | 看免费av毛片| 狠狠狠狠99中文字幕| 动漫黄色视频在线观看| 久久久国产成人精品二区| 丰满人妻一区二区三区视频av | 国产av不卡久久| 免费看光身美女| 欧美乱妇无乱码| 色在线成人网| 亚洲av成人精品一区久久| 乱人视频在线观看| 国产亚洲精品久久久久久毛片| 99热6这里只有精品| 九九久久精品国产亚洲av麻豆| 国产欧美日韩精品一区二区| 欧美在线一区亚洲| 亚洲专区中文字幕在线| 精品日产1卡2卡| 国产av不卡久久| 欧美日韩瑟瑟在线播放| 在线观看日韩欧美| 久久亚洲精品不卡| 少妇的逼水好多| 色精品久久人妻99蜜桃| 国产精品久久久久久亚洲av鲁大| 99久久成人亚洲精品观看| 淫秽高清视频在线观看| 小说图片视频综合网站| x7x7x7水蜜桃| 久久精品国产清高在天天线| 国产亚洲精品av在线| 免费看十八禁软件| 日本与韩国留学比较| 欧美大码av| 香蕉丝袜av| 国产极品精品免费视频能看的| 久久精品国产99精品国产亚洲性色| 美女高潮喷水抽搐中文字幕| 久久精品亚洲精品国产色婷小说| 成年人黄色毛片网站| 亚洲成人中文字幕在线播放| 99精品在免费线老司机午夜| 一个人观看的视频www高清免费观看| ponron亚洲| 国产久久久一区二区三区| 久久久久久久午夜电影| 亚洲av二区三区四区| 女警被强在线播放| www.色视频.com| 搡老熟女国产l中国老女人| 国产精华一区二区三区| 精华霜和精华液先用哪个| 免费观看精品视频网站| 99久久无色码亚洲精品果冻| 神马国产精品三级电影在线观看| 最近最新中文字幕大全电影3| 在线观看美女被高潮喷水网站 | 亚洲国产精品sss在线观看| 日韩大尺度精品在线看网址| 欧美zozozo另类| 亚洲精品亚洲一区二区| av黄色大香蕉| 偷拍熟女少妇极品色| 日本黄大片高清| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 在线观看66精品国产| 国产97色在线日韩免费| 一二三四社区在线视频社区8| 国产精品一区二区三区四区免费观看 | 欧美日韩综合久久久久久 | 91久久精品国产一区二区成人 | 99在线视频只有这里精品首页| 国产不卡一卡二| 中文字幕av在线有码专区| 久久精品91无色码中文字幕| 男女午夜视频在线观看| 中文字幕久久专区| 亚洲精华国产精华精| 日本与韩国留学比较| 日韩大尺度精品在线看网址| 天天躁日日操中文字幕| 亚洲精品一区av在线观看| 天堂网av新在线| 夜夜爽天天搞| 亚洲精品一区av在线观看| 最新美女视频免费是黄的| 91在线精品国自产拍蜜月 | 日韩欧美在线二视频| av在线天堂中文字幕| 亚洲一区二区三区色噜噜| 国产精品久久电影中文字幕| 国产精品自产拍在线观看55亚洲| 亚洲中文字幕一区二区三区有码在线看| 精品日产1卡2卡| 国产一区二区三区在线臀色熟女| 亚洲中文字幕一区二区三区有码在线看| 亚洲人成网站在线播| 久久久精品欧美日韩精品| 身体一侧抽搐| 成人18禁在线播放| 国产伦人伦偷精品视频| 亚洲成人免费电影在线观看| 99热6这里只有精品| 热99在线观看视频| 久久久久免费精品人妻一区二区| 麻豆国产av国片精品| 五月伊人婷婷丁香| 国产精品一区二区免费欧美| 露出奶头的视频|