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

    雙貯備系統(tǒng)冷/溫/熱貯備模型的優(yōu)化選擇研究

    2023-09-27 23:34:28金海波趙欣越
    自動(dòng)化學(xué)報(bào) 2023年9期
    關(guān)鍵詞:穩(wěn)態(tài)概率狀態(tài)

    金海波 趙欣越 桑 雨

    高可靠控制系統(tǒng)在現(xiàn)代工業(yè)的諸多領(lǐng)域中越來越重要,如航空航天、化工、核能、武器、生產(chǎn)制造等領(lǐng)域[1-3].實(shí)現(xiàn)控制系統(tǒng)高可靠性的一個(gè)主要手段是采用貯備技術(shù).目前,對(duì)于一些先進(jìn)生產(chǎn)制造企業(yè),具有雙貯備設(shè)備的控制系統(tǒng)已逐漸涌現(xiàn)并凸顯其作用[4].因此,對(duì)這些雙貯備系統(tǒng)的可靠性和經(jīng)濟(jì)性等指標(biāo)進(jìn)行定量分析和優(yōu)化具有重要的理論意義和實(shí)際價(jià)值.

    貯備系統(tǒng)按照貯備設(shè)備在系統(tǒng)運(yùn)行時(shí)是否失效通常分為: 冷貯備系統(tǒng)、熱貯備系統(tǒng)和溫貯備系統(tǒng)[5-6].冷貯備系統(tǒng)是指系統(tǒng)運(yùn)行過程中,貯備設(shè)備不參與工作也不老化.對(duì)能耗要求極高的系統(tǒng)往往配成冷貯備系統(tǒng),如冶金系統(tǒng)、武器系統(tǒng)等[7].熱貯備系統(tǒng)是指運(yùn)行設(shè)備和貯備設(shè)備在相同的環(huán)境下工作,因此兩種設(shè)備的失效率相同.熱貯備系統(tǒng)主要是對(duì)系統(tǒng)切換時(shí)間要求極高的系統(tǒng),如網(wǎng)絡(luò)打印機(jī)、飛機(jī)引擎等系統(tǒng)[8-9].溫貯備系統(tǒng)是指貯備設(shè)備在系統(tǒng)工作期間參與工作,但在溫和的環(huán)境中工作其失效率小于運(yùn)行設(shè)備.對(duì)需要平衡切換時(shí)間和能耗的系統(tǒng)大都配成溫貯備系統(tǒng),如電力系統(tǒng)、存儲(chǔ)系統(tǒng)、高性能計(jì)算系統(tǒng)和飛機(jī)控制系統(tǒng)等[10-13].

    由于三種貯備系統(tǒng)在不同應(yīng)用領(lǐng)域中都有重要的作用,因此學(xué)者們對(duì)每種貯備系統(tǒng)在故障檢測(cè)、失效分析、可靠性分析、冗余分配及優(yōu)化維修等方面都做了深入研究并取得了一定研究成果[14-15].在冷貯備系統(tǒng)研究方面,Chen 等[16]對(duì)兩部件冷貯備系統(tǒng)在不同失效機(jī)理的累積作用下進(jìn)行了可靠性分析.分析過程中作者考慮了部件工作狀態(tài)和貯備狀態(tài)之間的依賴關(guān)系,評(píng)估了系統(tǒng)在不同失效機(jī)理影響下的動(dòng)態(tài)演化過程,計(jì)算了不同階段應(yīng)力持續(xù)影響下的部件損壞程度,最后建立了基于改進(jìn)的序列二值決策圖方法的系統(tǒng)可靠性模型.Zhong 等[17]對(duì)雙設(shè)備組成的冷貯備控制系統(tǒng)提出了基于半Markov理論的最優(yōu)預(yù)防維護(hù)策略.所提策略中作者使用再生點(diǎn)技術(shù)和半Markov 過程刻畫了系統(tǒng)狀態(tài)轉(zhuǎn)移概率并用Markov 更新理論建立了系統(tǒng)每個(gè)狀態(tài)逗留時(shí)間分布的線性方程組,最后以最大化系統(tǒng)壽命為目標(biāo),優(yōu)化了系統(tǒng)最優(yōu)維護(hù)周期.Wang 等[18]研究了由雙設(shè)備組成的冷貯備系統(tǒng)的更換策略,首先假設(shè)系統(tǒng)惡化過程服從廣義Polya 過程,在此假設(shè)下,作者以最小化系統(tǒng)長期平均費(fèi)用率(Cost rate)為目標(biāo)函數(shù),給出了系統(tǒng)的最優(yōu)更換策略.陳童等[19]針對(duì)裝備系統(tǒng)中多狀態(tài)工作部件存在退化失效與突發(fā)失效競爭的情況,以冷貯備系統(tǒng)為研究對(duì)象,研究隨機(jī)檢測(cè)策略.其中工作部件在各性能水平停留時(shí)間、各類維修時(shí)間等隨機(jī)時(shí)間變量以及外部沖擊的到達(dá)過程均采用相位型 (Phase-type,PH)分布進(jìn)行描述.在此基礎(chǔ)上,建立了多狀態(tài)冷貯備系統(tǒng)可靠性模型,得到了系統(tǒng)可靠性主要參數(shù)的解析表達(dá)式.在熱貯備系統(tǒng)研究方面,研究成果相對(duì)較少,原因在于熱貯備系統(tǒng)中的工作設(shè)備和貯備設(shè)備失效率相同且以并聯(lián)方式工作.因此熱貯備的可靠性可完全借鑒并聯(lián)系統(tǒng)的可靠性方法進(jìn)行分析.曹晉華等[20]利用Markov 更新過程和交替更新過程研究了兩個(gè)不同部件組成的熱貯備系統(tǒng),推導(dǎo)了系統(tǒng)首次失效時(shí)間與其后停工時(shí)間的聯(lián)合分布、修理工在任一時(shí)刻忙的概率以及在(0,t]內(nèi)系統(tǒng)失效率的分布及其均值等性能指標(biāo).Patowary 等[21]采用Markov模型結(jié)合故障樹分析(Fault tree analysis,FTA)方法研究了熱貯備微電網(wǎng)系統(tǒng)的可靠性,所提方法與傳統(tǒng)Markov 和FTA 方法相比能夠適應(yīng)不同的系統(tǒng)失效率.與冷貯備系統(tǒng)相比,溫貯備系統(tǒng)中的貯備設(shè)備在貯備期間也逐漸惡化存在失效風(fēng)險(xiǎn),但失效風(fēng)險(xiǎn)小于熱貯備系統(tǒng).從貯備設(shè)備失效率角度看,冷貯備和熱貯備屬于溫貯備的兩個(gè)特例.因此溫貯備模型是更為一般的貯備模型,這也導(dǎo)致溫貯備系統(tǒng)的可靠性分析難度急劇增加.在溫貯備系統(tǒng)研究方面,近年來涌現(xiàn)出許多重要成果.Huang 等[22]以衛(wèi)星數(shù)據(jù)處理子系統(tǒng)為背景,研究了一類特殊結(jié)構(gòu)的溫貯備系統(tǒng),該類系統(tǒng)由兩組相同數(shù)量的不同部件組成,一組是工作部件,另一組是溫備份部件.兩組部件分別由兩個(gè)電源供電.對(duì)這類貯備系統(tǒng),作者給出了系統(tǒng)可靠度模型,建立了各個(gè)部件可靠度的閉合方程組,推導(dǎo)了系統(tǒng)可靠度的解析解.然而,該研究成果局限于各部件壽命都符合指數(shù)分布的情況,對(duì)其他壽命分布不再適用.Ma 等[23]研究了兩部件溫貯備冷卻系統(tǒng),采用多階段維納過程刻畫系統(tǒng)惡化趨勢(shì),提出了基于溫度檢測(cè)數(shù)據(jù)的優(yōu)化維護(hù)模型.該模型中,同時(shí)考慮穩(wěn)態(tài)溫度控制和系統(tǒng)壽命閾值,優(yōu)化系統(tǒng)整體維護(hù)費(fèi)用.尹東亮等[24]對(duì)具有維修和保養(yǎng)兩類活動(dòng)的多狀態(tài)溫貯備系統(tǒng)進(jìn)行了研究,考慮了維修較保養(yǎng)具有更高優(yōu)先級(jí)的情況.用PH 分布構(gòu)建了系統(tǒng)狀態(tài)轉(zhuǎn)移矩陣,推導(dǎo)了系統(tǒng)穩(wěn)態(tài)可用度、系統(tǒng)故障率、平均故障間隔時(shí)間等可靠性指標(biāo).劉寶亮等[25]研究了修理設(shè)備和開關(guān)不完全可靠情形下的溫貯備可修系統(tǒng),用補(bǔ)充變量法和Laplace 變換相結(jié)合的方式推導(dǎo)了系統(tǒng)瞬時(shí)可用度.該成果實(shí)際上是Kuo 等[26]研究成果的進(jìn)一步擴(kuò)展.

    由此可見,關(guān)于三種類型貯備系統(tǒng)的可靠性問題學(xué)者們做了大量研究,同時(shí)取得了許多重要成果.然而這些成果大多是針對(duì)“用一備一”的情況進(jìn)行研究.隨著某些尖端系統(tǒng)對(duì)高可靠性的要求,“用一備二”的情形逐漸涌現(xiàn).對(duì)該類系統(tǒng)的研究也剛剛起步,其研究成果也鮮有報(bào)道.目前,我國某自動(dòng)化廠商已經(jīng)開始研發(fā)具有雙貯備設(shè)備(即“用一備二”)的控制系統(tǒng),對(duì)該種控制系統(tǒng)配成冷/溫/熱三種模型中的何種模型是研發(fā)過程中需要解決的關(guān)鍵問題之一,具有現(xiàn)實(shí)意義.原因在于最優(yōu)貯備模型可以提高系統(tǒng)穩(wěn)態(tài)可用度(即可靠性),降低維修人員忙期穩(wěn)態(tài)概率以及系統(tǒng)穩(wěn)態(tài)平均維修次數(shù),延長系統(tǒng)壽命.除此之外,在經(jīng)濟(jì)方面最優(yōu)貯備模型還能降低系統(tǒng)維修費(fèi)用,提高系統(tǒng)單位時(shí)間內(nèi)產(chǎn)生的效益.然而據(jù)作者查閱大量相關(guān)文獻(xiàn)可知,目前關(guān)于雙貯備系統(tǒng)的最優(yōu)貯備模型選擇方面研究甚少,沒有可借鑒的通用方法.因此如何建立系統(tǒng)三種貯備模型并給出確定不同條件下最優(yōu)貯備模型的優(yōu)化選擇算法是主要研究難點(diǎn).為此,本文用Markov 及半Markov 更新理論結(jié)合Laplace 及Laplace-Stieltjes 變換技術(shù)分析系統(tǒng)狀態(tài)轉(zhuǎn)移概率和系統(tǒng)再生狀態(tài)的平均逗留時(shí)間,推導(dǎo)系統(tǒng)穩(wěn)態(tài)可用度、維修人員忙期穩(wěn)態(tài)概率和系統(tǒng)穩(wěn)態(tài)維修次數(shù)的可靠性指標(biāo)以及系統(tǒng)單位時(shí)間內(nèi)凈收益的經(jīng)濟(jì)指標(biāo),給出確定不同條件下貯備模型的優(yōu)化選擇算法.

    為便于理解,模型中主要變量和符號(hào)如表1 和表2 所示.

    表1 模型中主要變量說明Table 1 Main variables involved in models

    表2 模型中主要符號(hào)說明Table 2 Main symbols involved in models

    1 系統(tǒng)冷貯備模型分析

    該模型下系統(tǒng)由三個(gè)同類型的設(shè)備組成,其中一個(gè)運(yùn)行,另外兩個(gè)冷貯備(即貯備期間設(shè)備既不失效也不老化).令隨機(jī)變量X表示設(shè)備在運(yùn)行期間的壽命,Z表示失效設(shè)備的維修時(shí)間.因?yàn)殡娮釉O(shè)備的壽命多數(shù)服從指數(shù)分布或近似服從指數(shù)分布[27],而失效設(shè)備的維修時(shí)間由多種因素決定,如維修人員的維修水平、失效設(shè)備的復(fù)雜程度以及維修工具的先進(jìn)性等因素.因此維修時(shí)間往往不服從某一特定類型的分布[28].所以設(shè)運(yùn)行設(shè)備壽命服從參數(shù)為λ的指數(shù)分布F(t,λ),失效設(shè)備維修時(shí)間服從一般分布G(t),即X~F(t,λ),Z~G(t).經(jīng)分析,此系統(tǒng)共有6 個(gè)不同狀態(tài):

    為了建立系統(tǒng)在再生狀態(tài)下的更新方程,需要將這些狀態(tài)進(jìn)行狀態(tài)劃分.根據(jù)X~F(t,λ),Z~G(t),通過分析易知,系統(tǒng)進(jìn)入狀態(tài)S0,S1,S2,S3和S4的時(shí)刻均是系統(tǒng)的再生時(shí)刻(又稱再生點(diǎn)),而進(jìn)入狀態(tài)S5的時(shí)刻是非再生時(shí)刻.因此S0,S1,S2,S3和S4是再生狀態(tài),S5是非再生狀態(tài)(又稱滑過狀態(tài))且是失效狀態(tài).狀態(tài)之間的轉(zhuǎn)移關(guān)系如圖1所示.

    圖1 冷貯備系統(tǒng)狀態(tài)轉(zhuǎn)移圖Fig.1 State transition diagram of the cold-standby system

    1.1 系統(tǒng)半Markov 核函數(shù)

    令X(t)=Sj表示時(shí)刻t系統(tǒng)處于狀態(tài)Sj,Tn表示系統(tǒng)第n次狀態(tài)轉(zhuǎn)移時(shí)刻,Zn=X(Tn+0)表示第n次轉(zhuǎn)移時(shí)刻系統(tǒng)進(jìn)入的狀態(tài),容易驗(yàn)證{Zn,Tn,n ∈N} 是狀態(tài)空間E={Sj|j=0,···, 5} 上的Markov 更新過程,{X(t),t≥0} 是半Markov過程.因此需對(duì)系統(tǒng)的半Markov 核函數(shù)Qij(t),i,j=0,···, 5進(jìn)行分析.如圖1 所示.

    1)當(dāng)系統(tǒng)處于S0時(shí),如果運(yùn)行設(shè)備失效,則其中一個(gè)冷貯備設(shè)備被立刻激活,此時(shí)系統(tǒng)轉(zhuǎn)移至S1.因此,Q01(t)可表示為

    2)當(dāng)系統(tǒng)處于S1時(shí),如果冷貯備設(shè)備激活完成,則該設(shè)備進(jìn)入運(yùn)行狀態(tài),且失效設(shè)備進(jìn)入維修狀態(tài).此時(shí)系統(tǒng)轉(zhuǎn)移至S2.因此,Q12(t)可表示為

    3)當(dāng)系統(tǒng)處于S2時(shí),此時(shí)有下述兩種情況:

    a)如果維修設(shè)備在運(yùn)行設(shè)備失效前修好,則系統(tǒng)轉(zhuǎn)移至S0.因此,Q20(t)可表示為

    b)反之,如果運(yùn)行設(shè)備在維修設(shè)備修好前失效,則該設(shè)備進(jìn)入等待維修狀態(tài)且維修設(shè)備立刻暫停維修,冷貯備設(shè)備被激活.此時(shí)系統(tǒng)轉(zhuǎn)移至S3.因此,Q23(t)可表示為

    4)當(dāng)系統(tǒng)處于S3時(shí),冷貯備設(shè)備激活完成時(shí),該設(shè)備進(jìn)入運(yùn)行狀態(tài),暫停維修的設(shè)備開始繼續(xù)維修,此時(shí)系統(tǒng)轉(zhuǎn)移至 S4.因此,Q34(t)可表示為

    5)當(dāng)系統(tǒng)處于S4時(shí),此時(shí)有下述兩種情況:

    a)如果維修設(shè)備在運(yùn)行設(shè)備失效前修好,則修好后的設(shè)備進(jìn)入冷貯備狀態(tài),等待維修的設(shè)備開始維修,此時(shí)系統(tǒng)進(jìn)入S2.因此,Q42(t)可表示為

    b)反之,如果運(yùn)行設(shè)備在維修設(shè)備修好前失效,則系統(tǒng)轉(zhuǎn)移至S5.因此,Q45(t)可表示為

    1.2 系統(tǒng)所有再生狀態(tài)的平均逗留時(shí)間

    令μi和Fi(t),i=0,1,2,3,4,分別表示系統(tǒng)在狀態(tài)Si的平均逗留時(shí)間和概率分布函數(shù).令Pi(t)表示系統(tǒng)在狀態(tài)Si的存活函數(shù),即Pi(t)=1-Fi(t)=(t).對(duì)每個(gè)狀態(tài)的平均逗留時(shí)間μi進(jìn)行分析可得

    1.3 系統(tǒng)性能指標(biāo)

    對(duì)系統(tǒng)穩(wěn)態(tài)可用度、維修人員穩(wěn)態(tài)忙期概率以及系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)三個(gè)性能指標(biāo)進(jìn)行分析.

    1.3.1 系統(tǒng)穩(wěn)態(tài)可用度

    設(shè)當(dāng)前時(shí)刻系統(tǒng)處于狀態(tài)S0,經(jīng)過時(shí)間t后,根據(jù)系統(tǒng)狀態(tài)是否發(fā)生變化,存在兩種情況: 1)如果沒有發(fā)生變化,即系統(tǒng)依然停留在S0,此種情況可表示為P {X ≥t}.此時(shí)系統(tǒng)的瞬時(shí)可用度為A0(t)=P {X ≥t}·1=P {X ≥t}.2)如果系統(tǒng)狀態(tài)發(fā)生變化,如圖1 所示系統(tǒng)只能轉(zhuǎn)移到S1.此時(shí)根據(jù)Markov 更新理論,系統(tǒng)瞬時(shí)可用度為

    其中,符號(hào) “?” 表示卷積運(yùn)算符.根據(jù)概率加法原理可得:A0(t)=Q01(t)?A1(t)+P{X ≥t}.同理可得系統(tǒng)在其余再生狀態(tài)S1,S2,S3,S4下的瞬時(shí)可用度更新方程.從而系統(tǒng)瞬時(shí)可用度的更新方程組為

    對(duì)式(13)進(jìn)行Laplace 變換得

    解式(14)可得

    根據(jù)Abel 定理,系統(tǒng)穩(wěn)態(tài)可用度為

    1.3.2 維修人員穩(wěn)態(tài)忙期概率

    與第1.3.1 節(jié)類似,首先建立維修人員在t時(shí)刻忙期的瞬時(shí)概率更新方程組

    對(duì)式(17)進(jìn)行Laplace 變換得

    解式(18)可得

    根據(jù)Abel 定理,維修人員穩(wěn)態(tài)忙期概率為

    1.3.3 系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)

    同理,建立系統(tǒng)在 (0,t] 時(shí)間內(nèi)維修次數(shù)的更新方程組

    對(duì)式(21)進(jìn)行Laplace-Stieltjes 變換得

    解式(22)可得

    其中,Γ4(s)和 Γ6(s)與式(19)中相同.

    根據(jù)Abel 定理,可得系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)為

    其中,NV=p23+p42.

    1.3.4 系統(tǒng)穩(wěn)態(tài)激活概率

    建立系統(tǒng)在t時(shí)刻的激活概率更新方程組

    對(duì)式(25)進(jìn)行Laplace 變換得

    解式(26)可得

    Γ6(s)與式(19)中相同.

    根據(jù)Abel 定理,系統(tǒng)穩(wěn)態(tài)激活概率為

    1.4 系統(tǒng)單位時(shí)間內(nèi)產(chǎn)生的經(jīng)濟(jì)效益

    系統(tǒng)單位時(shí)間內(nèi)的經(jīng)濟(jì)效益 Θ1等于系統(tǒng)的運(yùn)行收益減去設(shè)備的維修費(fèi)用、支付給維修人員的費(fèi)用、系統(tǒng)激活期間的停工費(fèi)用以及系統(tǒng)的安裝費(fèi)用.因此

    其中,c0,c1,c2,c3分別表示系統(tǒng)單位時(shí)間內(nèi)的運(yùn)行收益、設(shè)備維修費(fèi)用、支付給維修人員的費(fèi)用、停工費(fèi)用,I表示一個(gè)設(shè)備的安裝費(fèi)用.

    2 系統(tǒng)溫貯備模型分析

    該模型下系統(tǒng)由三個(gè)同類型的設(shè)備組成,其中一個(gè)設(shè)備運(yùn)行,另外兩個(gè)設(shè)備溫貯備(即貯備設(shè)備在溫和的環(huán)境下運(yùn)行,其失效率低于運(yùn)行設(shè)備).令隨機(jī)變量Xi,i=1,2,3 表示第i個(gè)設(shè)備在運(yùn)行期間的壽命,Yi表示第i個(gè)設(shè)備在貯備期間的壽命,Zi表示第i個(gè)設(shè)備失效后的維修時(shí)間,F(t,λ)表示參數(shù)為λ的指數(shù)分布,G(t)和G1(t)表示一般分布.假設(shè)運(yùn)行設(shè)備和溫貯備設(shè)備的壽命分別服從參數(shù)為λ和λ1的指數(shù)分布,即Xi~F(t,λ),Yi~F(t,λ1).運(yùn)行設(shè)備和溫貯備設(shè)備失效后的維修時(shí)間均服從一般分布,但實(shí)際系統(tǒng)中這兩種維修時(shí)間往往不同,需要分別考慮.因此,如果Zi是運(yùn)行設(shè)備失效后的維修時(shí)間,則Zi~G(t);反之,如果Zi是溫貯備設(shè)備失效后的維修時(shí)間,則Zi~G1(t).為了便于描述模型,進(jìn)一步假設(shè):

    1)X1,X2,X3,Y1,Y2相互獨(dú)立.

    2)系統(tǒng)不同狀態(tài)之間的轉(zhuǎn)移是瞬時(shí)的.

    3)設(shè)備失效后,如果沒有其他設(shè)備正在維修,則維修人員立刻對(duì)該設(shè)備進(jìn)行維修;否則,該設(shè)備進(jìn)入等待狀態(tài)直到其他設(shè)備維修完成.設(shè)備修復(fù)后,其壽命分布像新的設(shè)備一樣.

    經(jīng)分析可得系統(tǒng)共有13 個(gè)狀態(tài):

    與冷貯備模型分析過程相同,經(jīng)分析可得:S0,S1,S2,S7和S12是再生狀態(tài);S3,S4,S5,S6,S8,S9,S10和S11是非再生狀態(tài);S5,S6,S8和S11是失效狀態(tài).狀態(tài)之間的轉(zhuǎn)移關(guān)系如圖2 所示.

    圖2 溫貯備系統(tǒng)狀態(tài)轉(zhuǎn)移圖Fig.2 State transition diagram of the warm-standby system

    2.1 系統(tǒng)半Markov 核函數(shù)

    與第1.1 節(jié)類似,經(jīng)分析可得溫貯備模型下系統(tǒng)所有半Markov 核函數(shù)為(具體分析過程見附錄A)

    2.2 系統(tǒng)在所有再生狀態(tài)的平均逗留時(shí)間

    與第1.2 節(jié)類似,經(jīng)分析可得系統(tǒng)在每個(gè)再生狀態(tài)的平均逗留時(shí)間μi,i=0,1,2,7,12 為(具體分析過程見附錄B)

    2.3 系統(tǒng)性能指標(biāo)

    本節(jié)對(duì)系統(tǒng)的穩(wěn)態(tài)可用度、維修人員穩(wěn)態(tài)忙期概率和系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)三個(gè)性能指標(biāo)進(jìn)行分析.

    2.3.1 系統(tǒng)穩(wěn)態(tài)可用度

    與第1.3.1 節(jié)類似,經(jīng)分析可得系統(tǒng)穩(wěn)態(tài)可用度為

    2.3.2 維修人員穩(wěn)態(tài)忙期概率

    與第1.3.2 節(jié)類似,經(jīng)分析可得維修人員穩(wěn)態(tài)忙期概率為

    2.3.3 系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)

    與第1.3.3 節(jié)類似,經(jīng)分析可得系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)為

    2.4 系統(tǒng)單位時(shí)間內(nèi)產(chǎn)生的經(jīng)濟(jì)收益

    與第1.4 節(jié)類似,溫貯備系統(tǒng)單位時(shí)間內(nèi)產(chǎn)生的經(jīng)濟(jì)效益 Θ2可表示為

    其中,參數(shù)c4,c5,c6的含義分別與式(29)中的c0,c1,c2相同.

    3 系統(tǒng)熱貯備模型分析

    雙貯備設(shè)備下熱貯備系統(tǒng)實(shí)質(zhì)上是三個(gè)設(shè)備以并聯(lián)的方式工作.該模型中不考慮某個(gè)設(shè)備失效后負(fù)載均衡對(duì)其他運(yùn)行設(shè)備失效率造成的影響.因此三個(gè)設(shè)備的失效率和維修率均視為相同.實(shí)際上,該模型的分析方法與溫貯備系統(tǒng)模型的分析方法相同,在分析過程中只需令λ1=λ,β1=β即可.因此,該模型的分析過程略.下面只給出系統(tǒng)單位時(shí)間內(nèi)的凈收益函數(shù)

    其中,參數(shù)c7,c8,c9的含義分別與式(29)中的c0,c1,c2相同.

    4 系統(tǒng)冷/溫/熱貯備模型優(yōu)化選擇算法

    從系統(tǒng)性能指標(biāo)(即穩(wěn)態(tài)可用度、維修人員穩(wěn)態(tài)忙期概率、系統(tǒng)穩(wěn)態(tài)平均維修次數(shù))以及經(jīng)濟(jì)指標(biāo)(即系統(tǒng)單位時(shí)間內(nèi)凈收益)兩個(gè)方面,分別給出雙貯備系統(tǒng)冷/溫/熱貯備模型的優(yōu)化選擇算法.

    算法 1.以系統(tǒng)性能為目標(biāo)的雙貯備系統(tǒng)冷/溫/熱貯備模型優(yōu)化選擇算法

    算法 2.以系統(tǒng)單位時(shí)間內(nèi)凈收益為目標(biāo)的雙貯備系統(tǒng)冷/溫/熱貯備模型優(yōu)化選擇算法

    5 實(shí)例分析

    以實(shí)際的雙貯備PLC (Programmable logic controller)控制系統(tǒng)(如圖3 所示)為研究對(duì)象,進(jìn)行實(shí)例分析.從系統(tǒng)性能指標(biāo),即系統(tǒng)穩(wěn)態(tài)可用度、維修人員穩(wěn)態(tài)忙期概率、系統(tǒng)穩(wěn)態(tài)平均維修次數(shù),以及經(jīng)濟(jì)指標(biāo),即系統(tǒng)單位時(shí)間內(nèi)凈收益兩個(gè)方面,對(duì)系統(tǒng)貯備模型優(yōu)化算法的輸入?yún)?shù)進(jìn)行假設(shè).

    圖3 雙貯備冗余控制系統(tǒng)Fig.3 Redundancy control system with dual-standby device

    為了給出合理的參數(shù)假設(shè),首先對(duì)圖3 所示的系統(tǒng)進(jìn)行簡要描述.該P(yáng)LC 控制系統(tǒng)是我國某自動(dòng)化廠商自主研制的高可靠雙貯備控制系統(tǒng).系統(tǒng)中每個(gè)PLC 設(shè)備經(jīng)過加速壽命測(cè)試后得到平均工作壽命約為1800 天,因此運(yùn)行設(shè)備的平均失效率設(shè)為λ=0.00055 (個(gè)/天).溫貯備設(shè)備失效率與工作環(huán)境有關(guān),設(shè)其為λ1∈[0.00001,0.00055] (個(gè)/天).維修時(shí)間一般服從指數(shù)分布,即G(t)=1-e-βt,.其中參數(shù)β,β1分別表示運(yùn)行設(shè)備和溫貯備設(shè)備失效后的平均維修率.實(shí)際上,設(shè)備失效后需要返廠、檢測(cè)、維修、測(cè)試、現(xiàn)場安裝調(diào)試、重新運(yùn)行等環(huán)節(jié).因此根據(jù)設(shè)備歷史維修數(shù)據(jù)可得運(yùn)行設(shè)備失效后從返廠到重新運(yùn)行約為7~30 天,溫貯備設(shè)備約為5 到15 天.因此設(shè)β ∈[0.03,0.14](個(gè)/天),β1=[0.067,0.2] (個(gè)/天).該系統(tǒng)如果配成冷貯備模式,則當(dāng)運(yùn)行設(shè)備失效后,冷貯備設(shè)備能夠自動(dòng)上電并上載控制程序以及導(dǎo)入設(shè)備失效前的數(shù)據(jù),整個(gè)激活過程需要3 min 左右.通過對(duì)歷史激活數(shù)據(jù)的統(tǒng)計(jì)分析可得貯備設(shè)備激活時(shí)間服從正態(tài)分布,即

    其中,參數(shù)γ和σ分別代表均值和方差.因此設(shè)γ=0.0021 (天),σ=0.0007.

    5.1 以系統(tǒng)性能為目標(biāo)的實(shí)例分析

    5.1.1 以穩(wěn)態(tài)可用度為目標(biāo)的實(shí)例分析

    令λ=0.00055,λ1=0.00011,β=0.05,β1=0.1,γ=0.0021,σ=0.0007,系統(tǒng)性能指標(biāo)為穩(wěn)態(tài)可用度.將這些參數(shù)代入算法1,在MATLAB (2014b)環(huán)境下運(yùn)行算法1 (公式推導(dǎo)部分利用MATLAB的符號(hào)計(jì)算),其計(jì)算結(jié)果如表3 所示.

    表3 系統(tǒng)穩(wěn)態(tài)可用度Table 3 System steady-state availability

    表3 中,模型I、模型II 和模型III 分別表示冷貯備模型、溫貯備模型和熱貯備模型.由表3 可知,.因此以系統(tǒng)穩(wěn)態(tài)可用度為目標(biāo)時(shí),算法1 輸出為: 冷貯備是最優(yōu)貯備模型.由于溫貯備設(shè)備的失效率與其工作環(huán)境有關(guān),因此屬于可變參數(shù),為了研究該參數(shù)對(duì)系統(tǒng)穩(wěn)態(tài)可用度的影響,令λ1從 0.00001 變化到 0.00055,步長為 0.00001,其他參數(shù)不變.將這些參數(shù)重新代入算法1,其計(jì)算結(jié)果如圖4 所示.

    圖4 λ1 對(duì)三個(gè)模型的穩(wěn)態(tài)可用度的影響Fig.4 Impact of λ1 on steady-state available of the three models

    由圖4 可知,溫貯備系統(tǒng)的穩(wěn)態(tài)可用度隨λ1的增加而降低.?λ1∈[0.00001,0.00055],溫貯備系統(tǒng)的穩(wěn)態(tài)可用度均小于冷貯備系統(tǒng),但均大于熱貯備系統(tǒng).可見溫貯備設(shè)備的失效率對(duì)溫貯備系統(tǒng)穩(wěn)態(tài)可用度有一定影響,但對(duì)算法1 輸出結(jié)果即最優(yōu)貯備模型沒有影響.

    5.1.2 以維修人員忙期穩(wěn)態(tài)概率為目標(biāo)的實(shí)例分析

    參數(shù)λ,λ1,β,β1,γ,σ取值與第5.1.1 節(jié)相同,系統(tǒng)性能指標(biāo)為維修人員忙期穩(wěn)態(tài)概率.將這些參數(shù)代入算法1,其計(jì)算結(jié)果如表4 所示.

    表4 維修人員忙期穩(wěn)態(tài)概率Table 4 Steady-state probability of repairmen busy

    由表4 可知,.因此以維修人員穩(wěn)態(tài)忙期概率為目標(biāo)時(shí),算法1 輸出結(jié)果為: 冷貯備是最優(yōu)貯備模型.實(shí)際上,運(yùn)行設(shè)備和溫貯備設(shè)備失效后的維修率(即參數(shù)β,β1)對(duì)維修人員穩(wěn)態(tài)忙期概率有直接影響,因此,令β∈[0.03,0.14],β1=[0.067,0.2],兩個(gè)參數(shù)的變化步長均設(shè)為0.005,其他參數(shù)不變,研究β對(duì),的影響以及β,β1對(duì)的影響.其結(jié)果分別如圖5 和圖6 所示.

    圖5 運(yùn)行設(shè)備維修率對(duì)冷、熱貯備系統(tǒng)中維修人員穩(wěn)態(tài)忙期概率的影響Fig.5 Impact of repair rates for the working device on steady-state probability of repairmen busy for cold,hot-standby system

    圖6 運(yùn)行設(shè)備和溫貯備設(shè)備的維修率對(duì)溫貯備系統(tǒng)中維修人員穩(wěn)態(tài)忙期概率的影響Fig.6 Impact of repair rates for the working and warm-standby devices on steady-state probability of repairmen busy for warm-standby system

    圖7 設(shè)備維修率對(duì)冷、熱貯備系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)的影響Fig.7 Impact of repair rate on mean repair number of the cold,hot-standby systems in steady-state

    圖8 運(yùn)行設(shè)備的維修率和溫貯備設(shè)備的失效率對(duì)系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)的影響Fig.8 Impact of repair rate of the working device and failure rate of the warm-standby device on repair number of the system in steady-state

    由圖5 可知,在參數(shù)β的取值范圍內(nèi)變化時(shí),均小于.說明以維修人員穩(wěn)態(tài)忙期概率為目標(biāo)時(shí),算法1 輸出是: 冷貯備是最優(yōu)貯備模型.由圖6 可知,隨β,β1的增大呈非線性減小.從數(shù)值計(jì)算結(jié)果可知: 當(dāng)β取值相同時(shí),總成立.說明在相同的運(yùn)行設(shè)備維修率下,算法1輸出依然是: 冷貯備是最優(yōu)貯備模型.

    5.1.3 以系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)為目標(biāo)的實(shí)例分析

    參數(shù)λ,λ1,β,β1,γ,σ取值與第5.1.1 節(jié)相同,系統(tǒng)性能指標(biāo)為系統(tǒng)穩(wěn)態(tài)平均維修次數(shù).將這些參數(shù)代入算法1.其計(jì)算結(jié)果如表5 所示.

    表5 系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)Table 5 Mean repair number of the system in steady-state

    5.2 以系統(tǒng)經(jīng)濟(jì)效益為目標(biāo)的實(shí)例分析

    根據(jù)系統(tǒng)應(yīng)用案例中的歷史財(cái)務(wù)數(shù)據(jù)和歷史維修費(fèi)用記錄,給出三個(gè)模型下凈收益中費(fèi)用的合理范圍: 設(shè)c0,c4,c7∈ [1000,5000] (元/h),c1,c5,c8∈[300,500] (元/h),c2,c6,c9∈[300,700] (元/h).實(shí)際上,冷貯備系統(tǒng)激活期間的停機(jī)對(duì)于不同行業(yè)的生產(chǎn)制造企業(yè)帶來的經(jīng)濟(jì)損失或產(chǎn)生的費(fèi)用大相徑庭、難以估計(jì).但根據(jù)使用該系統(tǒng)的某生產(chǎn)企業(yè)的停機(jī)記錄,停機(jī)費(fèi)用的合理范圍為:c3∈[500,1000](元/h),系統(tǒng)安裝費(fèi)用的合理范圍為:I ∈[100,300](元/h).

    5.2.1 以系統(tǒng)單位時(shí)間內(nèi)運(yùn)行收益為研究對(duì)象

    令參數(shù)λ,λ1,β,β1,γ,σ取值與第5.1.1 節(jié)相同,令參數(shù)c1=400,c2=500,c3=800,c5=100,c6=200,c8=300,c9=500,I=200,c0,c4,c7∈[1000,5000].將以上參數(shù)代入算法2,其計(jì)算結(jié)果如下.

    對(duì)于情況a),化簡后最終不等式為

    該不等式在?c0,c4,c7∈[1000,5000] 取值范圍內(nèi)有解,即 Θ1≥max(Θ2,Θ3)成立,此時(shí)算法2輸出結(jié)果是: 冷貯備是最優(yōu)貯備模型.

    對(duì)于情況b),化簡后的最終不等式為

    該不等式在?c0,c4,c7∈[1000,5000] 取值范圍內(nèi)有解,即 Θ2≥max(Θ1,Θ3)成立,此時(shí)算法2輸出結(jié)果是: 溫貯備是最優(yōu)貯備模型.

    對(duì)于情況c),化簡后最終不等式為

    該不等式在?c0,c4,c7∈[1000,5000] 取值范圍內(nèi)有解,即 Θ3≥max(Θ1,Θ2)成立,此時(shí)算法2輸出結(jié)果是: 熱貯備是最優(yōu)貯備模型.

    5.2.2 以系統(tǒng)單位時(shí)間內(nèi)設(shè)備維修費(fèi)用為研究對(duì)象

    令c0=3000,c4=3000,c7=3000,c1,c5,c8∈[300,500],其他參數(shù)與第5.2.1 節(jié)相同.將以上參數(shù)代入算法2,其計(jì)算結(jié)果如下.

    對(duì)于情況a),化簡后的最終不等式為

    該不等式在?c1,c5,c8∈[300,500] 取值范圍內(nèi)恒成立,即 Θ1≥max(Θ2,Θ3)恒成立,此時(shí)算法2輸出結(jié)果是: 冷貯備是最優(yōu)貯備模型.

    對(duì)于情況b)和情況c),化簡后分別得最終不等式為

    然而,在?c1,c5,c8∈[300,500] 取值范圍內(nèi),以上兩個(gè)不等式均無解.即,Θ2≥max(Θ1,Θ3)和Θ3≥max(Θ1,Θ2)均不成立,此時(shí)算法2 無輸出.說明以系統(tǒng)單位時(shí)間內(nèi)設(shè)備維修費(fèi)用為研究對(duì)象時(shí),無論參數(shù)如何取值,冷貯備均是最優(yōu)貯備模型.

    5.2.3 以單位時(shí)間內(nèi)支付給維修人員的費(fèi)用為研究對(duì)象

    令c0=3000,c4=3000,c7=3000,c2,c6,c9∈[300,700],其他參數(shù)取值與第5.2.1 節(jié)相同.將以上參數(shù)代入算法2,其計(jì)算結(jié)果如下.

    對(duì)于情況a),化簡后的最終不等式為

    該不等式在?c2,c6,c9∈[300,700] 取值范圍內(nèi)恒成立,即 Θ1≥max(Θ2,Θ3)恒成立.此時(shí)算法2輸出結(jié)果是: 冷貯備是最優(yōu)貯備模型.

    對(duì)于情況b)和情況c),化簡后分別得到的最終不等式為

    然而,在?c2,c6,c9∈[300,700] 取值范圍內(nèi),以上兩個(gè)不等式均無解,即 Θ2≥max(Θ1,Θ3)和Θ3≥max(Θ1,Θ2)均不成立.此時(shí)算法2 無輸出.說明以系統(tǒng)單位時(shí)間內(nèi)支付給維修人員費(fèi)用為研究對(duì)象時(shí),無論參數(shù)如何取值,冷貯備均是最優(yōu)貯備模型.

    6 結(jié)束語

    針對(duì)選擇哪種貯備模型才能使雙貯備系統(tǒng)實(shí)現(xiàn)性能和經(jīng)濟(jì)效益最優(yōu)的問題,本文創(chuàng)新性地提出了雙貯備系統(tǒng)貯備模型優(yōu)化選擇算法.通過分析系統(tǒng)狀態(tài)及半Markov 核函數(shù)分別建立了系統(tǒng)冷/溫/熱貯備模型下的更新方程組,利用Laplace、Laplace-Stieltjes 變換技術(shù)和Abel 定理求得了系統(tǒng)穩(wěn)態(tài)可用度、維修人員穩(wěn)態(tài)忙期概率和系統(tǒng)穩(wěn)態(tài)平均維修次數(shù)的系統(tǒng)性能指標(biāo),并給出了系統(tǒng)單位時(shí)間內(nèi)凈收益的目標(biāo)函數(shù),之后通過模型對(duì)比分析給出了分別以系統(tǒng)性能指標(biāo)和經(jīng)濟(jì)指標(biāo)為目標(biāo)的系統(tǒng)貯備模型優(yōu)化選擇算法.最后以實(shí)際的國產(chǎn)雙貯備控制系統(tǒng)作為研究對(duì)象,對(duì)所提算法進(jìn)行實(shí)例分析,實(shí)例結(jié)果表明所提算法能夠有效地確定系統(tǒng)在不同條件下的最優(yōu)貯備模型.本文是在系統(tǒng)確定參數(shù)或確定參數(shù)變化范圍的情況下進(jìn)行研究的.然而對(duì)于某些實(shí)際的工業(yè)現(xiàn)場,現(xiàn)場環(huán)境復(fù)雜,系統(tǒng)部分參數(shù)無法測(cè)量,這些參數(shù)屬于不確定參數(shù).對(duì)具有不確定參數(shù)的雙貯備系統(tǒng)如何進(jìn)行分析、建模并給出貯備模型的優(yōu)化選擇算法是下一步重點(diǎn)研究的問題.

    附錄A

    對(duì)溫貯備模型下系統(tǒng)所有半Markov 核函數(shù)進(jìn)行分析如下.

    1)當(dāng)系統(tǒng)處于S0時(shí),如果其中一個(gè)溫貯備設(shè)備先于運(yùn)行設(shè)備失效,則系統(tǒng)轉(zhuǎn)移至S1.此時(shí)有

    反之,如果運(yùn)行設(shè)備先于溫貯備設(shè)備失效,則系統(tǒng)轉(zhuǎn)移至S2.此時(shí)有

    2)當(dāng)系統(tǒng)處于S1時(shí),此時(shí)有下述幾種情況:

    a)如果失效設(shè)備在運(yùn)行設(shè)備和溫貯備設(shè)備失效前已修好,則系統(tǒng)轉(zhuǎn)移至S0.此時(shí)有

    b)如果溫貯備設(shè)備在失效設(shè)備修好前失效且在該失效時(shí)刻運(yùn)行設(shè)備依然運(yùn)行,則系統(tǒng)轉(zhuǎn)移至S3,然而由于S3是非再生狀態(tài),因此當(dāng)系統(tǒng)進(jìn)入S3后會(huì)再次轉(zhuǎn)移至哪些狀態(tài)需要分別考慮.

    c)如果維修設(shè)備在運(yùn)行設(shè)備失效前已修好,則系統(tǒng)由S3轉(zhuǎn)移至S1,此時(shí)有

    d)如果運(yùn)行設(shè)備在失效設(shè)備修好前失效,則系統(tǒng)將由S3再次轉(zhuǎn)移至S5,此時(shí)有

    e)當(dāng)系統(tǒng)轉(zhuǎn)移至S5后,正在維修的設(shè)備修好后,系統(tǒng)將轉(zhuǎn)移至S7.此時(shí)有

    當(dāng)系統(tǒng)處于S1時(shí),如果運(yùn)行設(shè)備在維修設(shè)備修好前失效且在失效時(shí)刻溫貯備設(shè)備沒有失效,則系統(tǒng)轉(zhuǎn)移至S4.由于S4是非再生狀態(tài),因此系統(tǒng)由S4轉(zhuǎn)移至哪些狀態(tài)需要分別考慮.

    f)如果失效設(shè)備在運(yùn)行設(shè)備失效前已修好,則系統(tǒng)將由S4轉(zhuǎn)移至S2.此時(shí)有

    g)反之,如果運(yùn)行設(shè)備在失效設(shè)備修好前失效,則系統(tǒng)將由S4轉(zhuǎn)移至S6.此時(shí)有

    h)當(dāng)系統(tǒng)轉(zhuǎn)移至S6后,正在維修的設(shè)備修好后,系統(tǒng)將轉(zhuǎn)移至S7.此時(shí)有

    3)當(dāng)系統(tǒng)處于S2時(shí),此時(shí)有下述幾種情況:

    a)如果維修設(shè)備在運(yùn)行設(shè)備和溫貯備設(shè)備失效前已修好,則系統(tǒng)轉(zhuǎn)移至S0.此時(shí)有

    b)如果溫貯備設(shè)備在維修設(shè)備修好前失效,且在失效時(shí)刻運(yùn)行設(shè)備正常運(yùn)行,則系統(tǒng)轉(zhuǎn)移至S10,由于S10是非再生狀態(tài),因此系統(tǒng)由S10轉(zhuǎn)移至哪些狀態(tài)需要分別考慮.

    c)如果維修設(shè)備在運(yùn)行設(shè)備失效前已修好,則系統(tǒng)由S10轉(zhuǎn)移至S1.此時(shí)有

    d)反之,如果運(yùn)行設(shè)備在維修設(shè)備修好前失效,則系統(tǒng)由S10轉(zhuǎn)移至S11.此時(shí)有

    e)當(dāng)系統(tǒng)轉(zhuǎn)移至S11后,維修設(shè)備修好后,系統(tǒng)轉(zhuǎn)移至S12.此時(shí)有

    當(dāng)系統(tǒng)處于S2時(shí),如果運(yùn)行設(shè)備在維修設(shè)備修好前失效且在失效時(shí)刻溫貯備設(shè)備沒有失效,則系統(tǒng)轉(zhuǎn)移至S9.由于S9是非再生狀態(tài),因此由S9轉(zhuǎn)移至哪些狀態(tài)需要分別考慮.

    f)如果維修設(shè)備在運(yùn)行設(shè)備失效前已修好,則系統(tǒng)由S9轉(zhuǎn)移至S2.此時(shí)有

    g)如果運(yùn)行設(shè)備在維修設(shè)備修好前失效,則系統(tǒng)由S9轉(zhuǎn)移至S8.系統(tǒng)處于S8后,維修設(shè)備修好后,系統(tǒng)由S8轉(zhuǎn)移至S7.此時(shí)有

    4)當(dāng)系統(tǒng)處于S7時(shí),此時(shí)有下述幾種情況:

    a)如果正在維修的設(shè)備在運(yùn)行設(shè)備失效前已修好,則系統(tǒng)轉(zhuǎn)移至S2.此時(shí)有

    b)反之,如果運(yùn)行設(shè)備在維修設(shè)備修好前失效,則系統(tǒng)轉(zhuǎn)移至S6.系統(tǒng)處于S6后,當(dāng)維修設(shè)備修好后,系統(tǒng)由S6轉(zhuǎn)移至S7.此時(shí)有

    5)當(dāng)系統(tǒng)處于S12時(shí),此時(shí)有下述幾種情況:a)如果維修設(shè)備在運(yùn)行設(shè)備失效前已修好,則系統(tǒng)轉(zhuǎn)移至S2,此時(shí)有

    b)反之,如果運(yùn)行設(shè)備在維修設(shè)備修好前失效,則系統(tǒng)轉(zhuǎn)移至S6.系統(tǒng)處于S6后,當(dāng)維修設(shè)備修好后,系統(tǒng)由S6轉(zhuǎn)移至S7.此時(shí)有

    附錄B

    對(duì)溫貯備模型在所有再生狀態(tài)的平均逗留時(shí)間進(jìn)行分析如下.

    1)從系統(tǒng)進(jìn)入S0開始(t=0),經(jīng)過時(shí)間t后,如果運(yùn)行設(shè)備和溫貯備設(shè)備均未失效,則有

    2)從系統(tǒng)進(jìn)入S1開始(t=0),經(jīng)過時(shí)間t后,如果運(yùn)行設(shè)備和溫貯備設(shè)備沒有失效,且正在維修的設(shè)備沒有修好,則有

    3)從系統(tǒng)進(jìn)入S2開始(t=0),經(jīng)過時(shí)間t后,如果運(yùn)行設(shè)備和溫貯備設(shè)備沒有失效,且正在維修的設(shè)備沒有修好,則有

    4)從系統(tǒng)進(jìn)入S7開始(t=0),經(jīng)過時(shí)間t后,如果運(yùn)行設(shè)備沒有失效且正在維修的設(shè)備沒有修好,則有

    5)從系統(tǒng)進(jìn)入S12開始(t=0),經(jīng)過時(shí)間t后,如果運(yùn)行設(shè)備沒有失效且正在維修的設(shè)備沒有修好,則有

    猜你喜歡
    穩(wěn)態(tài)概率狀態(tài)
    可變速抽水蓄能機(jī)組穩(wěn)態(tài)運(yùn)行特性研究
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    概率與統(tǒng)計(jì)(一)
    概率與統(tǒng)計(jì)(二)
    電廠熱力系統(tǒng)穩(wěn)態(tài)仿真軟件開發(fā)
    煤氣與熱力(2021年4期)2021-06-09 06:16:54
    元中期歷史劇對(duì)社會(huì)穩(wěn)態(tài)的皈依與維護(hù)
    中華戲曲(2020年1期)2020-02-12 02:28:18
    狀態(tài)聯(lián)想
    生命的另一種狀態(tài)
    成人特级av手机在线观看| 草草在线视频免费看| 男女床上黄色一级片免费看| 一边摸一边抽搐一进一小说| 欧美日韩亚洲国产一区二区在线观看| 国产熟女xx| 精品人妻偷拍中文字幕| 精品久久久久久久久av| 嫩草影院入口| 婷婷丁香在线五月| 午夜日韩欧美国产| 国产乱人视频| 亚洲av电影在线进入| 国产欧美日韩一区二区精品| 免费人成在线观看视频色| 看十八女毛片水多多多| 一a级毛片在线观看| 日本免费a在线| 亚洲国产色片| 可以在线观看毛片的网站| 欧美三级亚洲精品| 国产精品一区二区三区四区久久| 久久久久九九精品影院| 精品一区二区三区av网在线观看| 精品一区二区三区人妻视频| 精品乱码久久久久久99久播| 日本一二三区视频观看| 国内精品久久久久精免费| 国产精品98久久久久久宅男小说| 男女做爰动态图高潮gif福利片| 69人妻影院| 性欧美人与动物交配| 国产精品不卡视频一区二区 | 久99久视频精品免费| 亚洲在线观看片| 亚洲成人免费电影在线观看| 嫁个100分男人电影在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产在线男女| 99久久精品国产亚洲精品| 中文字幕人妻熟人妻熟丝袜美| 99久久成人亚洲精品观看| 窝窝影院91人妻| 成年版毛片免费区| 国产一区二区三区视频了| 最好的美女福利视频网| 男人舔女人下体高潮全视频| 亚洲精品粉嫩美女一区| 十八禁人妻一区二区| 国产成人欧美在线观看| 男女那种视频在线观看| 9191精品国产免费久久| 亚洲国产高清在线一区二区三| 久久精品人妻少妇| 欧美又色又爽又黄视频| 国产一级毛片七仙女欲春2| 国产一级毛片七仙女欲春2| 2021天堂中文幕一二区在线观| 99久久精品热视频| 亚洲国产欧美人成| 色播亚洲综合网| 麻豆成人午夜福利视频| 亚洲欧美精品综合久久99| 看免费av毛片| 成人高潮视频无遮挡免费网站| 日本黄色视频三级网站网址| 少妇人妻一区二区三区视频| 美女黄网站色视频| 俄罗斯特黄特色一大片| 一本久久中文字幕| 国产单亲对白刺激| 久久精品国产亚洲av天美| 日韩精品中文字幕看吧| 久久精品国产清高在天天线| 久久久国产成人免费| 亚洲av成人精品一区久久| 国产精品久久久久久精品电影| 精品免费久久久久久久清纯| 别揉我奶头 嗯啊视频| 网址你懂的国产日韩在线| 男女做爰动态图高潮gif福利片| 18美女黄网站色大片免费观看| 国产精品永久免费网站| 日韩欧美 国产精品| 亚洲18禁久久av| 好男人在线观看高清免费视频| 亚洲在线观看片| 每晚都被弄得嗷嗷叫到高潮| av在线天堂中文字幕| 女人被狂操c到高潮| 亚洲最大成人中文| 黄片小视频在线播放| 久久性视频一级片| netflix在线观看网站| 午夜免费成人在线视频| 日本a在线网址| 内地一区二区视频在线| 又紧又爽又黄一区二区| 国产精品一及| 国产中年淑女户外野战色| 赤兔流量卡办理| 日本五十路高清| 欧美一区二区亚洲| 特大巨黑吊av在线直播| 亚洲人成网站在线播放欧美日韩| 久久精品综合一区二区三区| 特大巨黑吊av在线直播| 国产精品一区二区性色av| 国产黄a三级三级三级人| 久99久视频精品免费| 国产亚洲精品av在线| 午夜福利18| 亚洲七黄色美女视频| 日本三级黄在线观看| 一级黄色大片毛片| 国产亚洲欧美在线一区二区| 国产69精品久久久久777片| 在线十欧美十亚洲十日本专区| 最新中文字幕久久久久| 国产欧美日韩一区二区精品| 在线十欧美十亚洲十日本专区| 黄色女人牲交| 国产淫片久久久久久久久 | 老熟妇乱子伦视频在线观看| 久久精品影院6| 三级毛片av免费| 国产伦在线观看视频一区| 亚洲自拍偷在线| 超碰av人人做人人爽久久| 国产成人影院久久av| 亚洲狠狠婷婷综合久久图片| 精品午夜福利在线看| 久久人人爽人人爽人人片va | 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 久久久久亚洲av毛片大全| 男女床上黄色一级片免费看| 一区二区三区免费毛片| 老司机午夜福利在线观看视频| 国产av麻豆久久久久久久| 在线天堂最新版资源| 国产不卡一卡二| 我要搜黄色片| 成人性生交大片免费视频hd| 九色国产91popny在线| 国产黄片美女视频| 18禁黄网站禁片免费观看直播| a级毛片a级免费在线| 别揉我奶头~嗯~啊~动态视频| 久久久精品大字幕| 国产精品一区二区性色av| 搡女人真爽免费视频火全软件 | 亚洲美女视频黄频| 亚洲av.av天堂| 露出奶头的视频| 久9热在线精品视频| 黄色丝袜av网址大全| 国产精品国产高清国产av| 日韩欧美精品v在线| 乱人视频在线观看| 国产真实乱freesex| 午夜久久久久精精品| 精品福利观看| av国产免费在线观看| 我要搜黄色片| 狂野欧美白嫩少妇大欣赏| 午夜精品一区二区三区免费看| 又紧又爽又黄一区二区| 黄色女人牲交| 中亚洲国语对白在线视频| 成人亚洲精品av一区二区| 好看av亚洲va欧美ⅴa在| 国产三级中文精品| 欧美午夜高清在线| 白带黄色成豆腐渣| 国内揄拍国产精品人妻在线| 亚洲精品久久国产高清桃花| 日韩欧美国产在线观看| 亚洲第一区二区三区不卡| 国产精品免费一区二区三区在线| 欧美黑人巨大hd| 丰满人妻熟妇乱又伦精品不卡| h日本视频在线播放| 欧美bdsm另类| 老鸭窝网址在线观看| 免费看日本二区| 免费观看精品视频网站| 亚洲精品亚洲一区二区| 超碰av人人做人人爽久久| 简卡轻食公司| 国产又黄又爽又无遮挡在线| 亚洲精品一卡2卡三卡4卡5卡| 一区福利在线观看| 亚洲国产精品成人综合色| 蜜桃久久精品国产亚洲av| 国产精品国产高清国产av| 国产午夜精品久久久久久一区二区三区 | 老熟妇仑乱视频hdxx| 如何舔出高潮| 久久国产精品影院| 黄片小视频在线播放| 欧美色视频一区免费| 可以在线观看毛片的网站| av在线观看视频网站免费| 亚洲狠狠婷婷综合久久图片| 国产69精品久久久久777片| 亚洲一区二区三区不卡视频| 欧美又色又爽又黄视频| 亚洲国产精品sss在线观看| 成年女人永久免费观看视频| 精品欧美国产一区二区三| 三级男女做爰猛烈吃奶摸视频| 午夜日韩欧美国产| 亚洲av二区三区四区| 99视频精品全部免费 在线| 国产男靠女视频免费网站| 亚洲精品粉嫩美女一区| 99热6这里只有精品| 国产伦在线观看视频一区| 精品久久久久久久久亚洲 | 一本综合久久免费| 亚洲中文字幕日韩| 91字幕亚洲| 97超视频在线观看视频| 日本免费a在线| 国产成人av教育| 国产在线男女| 亚洲精品亚洲一区二区| 精品久久久久久久久av| 亚洲av第一区精品v没综合| 欧美中文日本在线观看视频| 淫妇啪啪啪对白视频| 成人高潮视频无遮挡免费网站| 久久欧美精品欧美久久欧美| 国产精品日韩av在线免费观看| 久久精品综合一区二区三区| 国产精品久久视频播放| 欧美成人一区二区免费高清观看| 国产久久久一区二区三区| 少妇熟女aⅴ在线视频| 欧美不卡视频在线免费观看| 男人舔女人下体高潮全视频| av视频在线观看入口| 日本在线视频免费播放| 国产精品99久久久久久久久| ponron亚洲| 精品久久久久久久久亚洲 | www.色视频.com| 精品久久久久久久末码| 国产中年淑女户外野战色| 天堂av国产一区二区熟女人妻| 国产男靠女视频免费网站| 亚洲五月天丁香| 蜜桃亚洲精品一区二区三区| 免费无遮挡裸体视频| 草草在线视频免费看| 国产精品嫩草影院av在线观看 | 国内精品久久久久久久电影| 欧美潮喷喷水| 国产真实伦视频高清在线观看 | 男女视频在线观看网站免费| 亚洲人成伊人成综合网2020| 看黄色毛片网站| 日韩免费av在线播放| 久久久久久国产a免费观看| 久久九九热精品免费| 亚洲avbb在线观看| 色综合站精品国产| 亚洲av免费高清在线观看| 精品人妻一区二区三区麻豆 | 亚洲,欧美精品.| 日韩有码中文字幕| 成人永久免费在线观看视频| 亚洲激情在线av| 国产中年淑女户外野战色| 国产精品国产高清国产av| 搡女人真爽免费视频火全软件 | 亚洲av中文字字幕乱码综合| www日本黄色视频网| 男人和女人高潮做爰伦理| 伊人久久精品亚洲午夜| 国产精品,欧美在线| 97超级碰碰碰精品色视频在线观看| 老熟妇乱子伦视频在线观看| 精品午夜福利在线看| xxxwww97欧美| 欧美精品啪啪一区二区三区| 九九久久精品国产亚洲av麻豆| 欧美成人一区二区免费高清观看| 精品一区二区三区av网在线观看| 国内少妇人妻偷人精品xxx网站| a级毛片a级免费在线| 久久午夜亚洲精品久久| 天堂av国产一区二区熟女人妻| 国产高潮美女av| 国产午夜精品论理片| 丁香六月欧美| 国产极品精品免费视频能看的| 国产欧美日韩精品亚洲av| 日韩亚洲欧美综合| 精品久久久久久成人av| 国产精品1区2区在线观看.| 国产亚洲欧美在线一区二区| 国产久久久一区二区三区| 久久精品人妻少妇| 久久久久九九精品影院| 欧美+日韩+精品| 欧美精品国产亚洲| 欧美丝袜亚洲另类 | 如何舔出高潮| 国产精品1区2区在线观看.| 亚洲av电影在线进入| 国产一级毛片七仙女欲春2| 一个人免费在线观看的高清视频| 亚洲精品乱码久久久v下载方式| 国产精品一区二区性色av| 成年女人永久免费观看视频| 欧美另类亚洲清纯唯美| 99久久久亚洲精品蜜臀av| 国产精品自产拍在线观看55亚洲| av黄色大香蕉| 久久欧美精品欧美久久欧美| 久久久久国内视频| 日日夜夜操网爽| 2021天堂中文幕一二区在线观| 午夜福利免费观看在线| 三级毛片av免费| 久久久久久久久中文| h日本视频在线播放| 一本精品99久久精品77| 日韩精品中文字幕看吧| 深夜精品福利| 亚洲五月婷婷丁香| 国产精品野战在线观看| 国产高清激情床上av| 免费观看的影片在线观看| 淫妇啪啪啪对白视频| 波多野结衣高清无吗| 深爱激情五月婷婷| 亚洲久久久久久中文字幕| av福利片在线观看| 国产精品久久久久久精品电影| 午夜老司机福利剧场| 国产综合懂色| 国产av在哪里看| 毛片女人毛片| 久久久久久久久中文| 中出人妻视频一区二区| 欧美黄色淫秽网站| 国语自产精品视频在线第100页| 欧美黄色片欧美黄色片| 老司机午夜十八禁免费视频| 黄色配什么色好看| 午夜免费成人在线视频| 97热精品久久久久久| 免费人成在线观看视频色| 久久精品国产亚洲av香蕉五月| 男人和女人高潮做爰伦理| 国产精品98久久久久久宅男小说| 欧美一区二区亚洲| 亚洲av不卡在线观看| 国产高清视频在线播放一区| 亚洲av熟女| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲电影在线观看av| 黄色视频,在线免费观看| 级片在线观看| av在线观看视频网站免费| 一本久久中文字幕| 夜夜躁狠狠躁天天躁| 神马国产精品三级电影在线观看| 90打野战视频偷拍视频| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 麻豆av噜噜一区二区三区| 国产精品影院久久| 亚洲成a人片在线一区二区| 免费看a级黄色片| 97超级碰碰碰精品色视频在线观看| 人人妻,人人澡人人爽秒播| 美女 人体艺术 gogo| 大型黄色视频在线免费观看| 99热这里只有是精品在线观看 | 无人区码免费观看不卡| av中文乱码字幕在线| 最新在线观看一区二区三区| 亚洲自拍偷在线| 免费看日本二区| 午夜福利18| 免费看日本二区| av国产免费在线观看| 欧美极品一区二区三区四区| 久久精品国产亚洲av香蕉五月| 尤物成人国产欧美一区二区三区| 亚洲自偷自拍三级| 波多野结衣巨乳人妻| a级毛片a级免费在线| 国产熟女xx| 免费观看人在逋| 99国产精品一区二区三区| 亚洲七黄色美女视频| 成人美女网站在线观看视频| 欧美日韩瑟瑟在线播放| 伊人久久精品亚洲午夜| 国产黄a三级三级三级人| 很黄的视频免费| 伊人久久精品亚洲午夜| 99在线人妻在线中文字幕| www.熟女人妻精品国产| 欧美又色又爽又黄视频| 男插女下体视频免费在线播放| 成人高潮视频无遮挡免费网站| 国产亚洲欧美在线一区二区| 99久国产av精品| 91久久精品电影网| 男女之事视频高清在线观看| 一本精品99久久精品77| 亚洲成人精品中文字幕电影| 久9热在线精品视频| 老熟妇乱子伦视频在线观看| 色精品久久人妻99蜜桃| 精品一区二区三区人妻视频| 在线免费观看不下载黄p国产 | 亚洲综合色惰| 久久人妻av系列| 国产国拍精品亚洲av在线观看| 窝窝影院91人妻| or卡值多少钱| 国产精品国产高清国产av| 午夜福利在线观看吧| 久久香蕉精品热| 给我免费播放毛片高清在线观看| 午夜福利高清视频| 国产成人欧美在线观看| 国产精品av视频在线免费观看| 亚洲不卡免费看| 宅男免费午夜| 麻豆一二三区av精品| 国产一区二区三区在线臀色熟女| 99国产精品一区二区三区| 嫩草影院入口| 亚洲专区国产一区二区| 十八禁人妻一区二区| 亚洲成人久久性| 一级黄片播放器| 一级a爱片免费观看的视频| 少妇的逼好多水| 丁香六月欧美| 天美传媒精品一区二区| 99视频精品全部免费 在线| 久9热在线精品视频| 久久欧美精品欧美久久欧美| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 午夜视频国产福利| 91在线精品国自产拍蜜月| 村上凉子中文字幕在线| 免费黄网站久久成人精品 | 国产精品不卡视频一区二区 | av在线蜜桃| aaaaa片日本免费| 中文字幕人成人乱码亚洲影| 亚洲精品乱码久久久v下载方式| 欧美午夜高清在线| 真人做人爱边吃奶动态| 亚洲av熟女| 欧美三级亚洲精品| x7x7x7水蜜桃| 一区福利在线观看| 国产探花在线观看一区二区| 我要搜黄色片| 午夜两性在线视频| 热99re8久久精品国产| 丁香欧美五月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美日韩亚洲国产一区二区在线观看| 在线免费观看的www视频| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 免费av不卡在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲乱码一区二区免费版| 中文字幕免费在线视频6| 深夜精品福利| 欧美日韩国产亚洲二区| 最后的刺客免费高清国语| 好男人电影高清在线观看| 特大巨黑吊av在线直播| 亚洲国产欧洲综合997久久,| 欧美激情国产日韩精品一区| 十八禁人妻一区二区| 久久精品人妻少妇| 蜜桃久久精品国产亚洲av| 长腿黑丝高跟| 亚洲va日本ⅴa欧美va伊人久久| 俺也久久电影网| 此物有八面人人有两片| 亚洲精品乱码久久久v下载方式| 久久人人爽人人爽人人片va | a级毛片免费高清观看在线播放| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久com| 午夜免费成人在线视频| av黄色大香蕉| 国产一区二区三区在线臀色熟女| 婷婷精品国产亚洲av| 能在线免费观看的黄片| 男女下面进入的视频免费午夜| 国产色婷婷99| 国产高清有码在线观看视频| 麻豆久久精品国产亚洲av| 一级av片app| 亚洲av美国av| 国产免费男女视频| 久久性视频一级片| 久久国产乱子伦精品免费另类| 国产又黄又爽又无遮挡在线| 极品教师在线免费播放| 国产视频一区二区在线看| x7x7x7水蜜桃| 欧美午夜高清在线| 国产精品不卡视频一区二区 | 此物有八面人人有两片| www.熟女人妻精品国产| 一个人免费在线观看的高清视频| 久久国产精品影院| 国产日本99.免费观看| 欧美午夜高清在线| 日日摸夜夜添夜夜添av毛片 | 又黄又爽又刺激的免费视频.| 久久久精品欧美日韩精品| 波野结衣二区三区在线| 少妇人妻精品综合一区二区 | 蜜桃亚洲精品一区二区三区| www.999成人在线观看| 成人永久免费在线观看视频| 国产亚洲精品久久久com| 亚洲成人精品中文字幕电影| 日本一本二区三区精品| 搡老妇女老女人老熟妇| 999久久久精品免费观看国产| 国产单亲对白刺激| 长腿黑丝高跟| 伦理电影大哥的女人| 国产美女午夜福利| 久久久成人免费电影| 国产亚洲精品综合一区在线观看| 久久久久国内视频| 国产野战对白在线观看| 精品一区二区三区av网在线观看| 日本a在线网址| 亚洲一区二区三区色噜噜| 国产单亲对白刺激| 国产av不卡久久| 老熟妇乱子伦视频在线观看| 午夜激情欧美在线| 国产精品一及| 老司机深夜福利视频在线观看| 亚洲欧美日韩高清专用| 男女下面进入的视频免费午夜| 国产伦精品一区二区三区四那| 亚洲精品影视一区二区三区av| 久久欧美精品欧美久久欧美| 国产综合懂色| 男女之事视频高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美在线二视频| 国产精品爽爽va在线观看网站| 亚洲人成网站在线播放欧美日韩| 首页视频小说图片口味搜索| 长腿黑丝高跟| 国产精品98久久久久久宅男小说| 99久久精品国产亚洲精品| 深夜精品福利| 激情在线观看视频在线高清| 国产麻豆成人av免费视频| 欧美午夜高清在线| 内射极品少妇av片p| 99精品久久久久人妻精品| 国产伦精品一区二区三区视频9| netflix在线观看网站| 久久性视频一级片| 午夜福利欧美成人| 久久久久久久久中文| 嫁个100分男人电影在线观看| 搡老岳熟女国产| 国产伦精品一区二区三区四那| 成人无遮挡网站| 亚洲成人精品中文字幕电影| 欧美潮喷喷水| 一边摸一边抽搐一进一小说| 日韩精品青青久久久久久| 久久国产精品人妻蜜桃| 日韩欧美免费精品| 一级黄色大片毛片| 可以在线观看的亚洲视频| 欧美在线黄色| 午夜免费男女啪啪视频观看 | avwww免费| 少妇的逼水好多| 青草久久国产| 日本 av在线| 91九色精品人成在线观看| 好男人在线观看高清免费视频| 搡女人真爽免费视频火全软件 | 免费看日本二区| 国产精品日韩av在线免费观看| 熟女电影av网| 蜜桃亚洲精品一区二区三区| 成人特级av手机在线观看| 亚洲中文字幕日韩| 欧美另类亚洲清纯唯美| 欧美性猛交╳xxx乱大交人| 每晚都被弄得嗷嗷叫到高潮| 99久久成人亚洲精品观看|