• <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| 少妇被粗大猛烈的视频| 极品人妻少妇av视频| 久久久久精品久久久久真实原创| 在线观看三级黄色| 日日摸夜夜添夜夜爱| 欧美 亚洲 国产 日韩一| av又黄又爽大尺度在线免费看| 亚洲欧洲日产国产| 午夜福利,免费看| 国产精品.久久久| av.在线天堂| 精品人妻偷拍中文字幕| 精品一区二区免费观看| 国产亚洲最大av| 久久精品久久精品一区二区三区| 午夜激情av网站| 国产一级毛片在线| 美女中出高潮动态图| 人妻一区二区av| 黑丝袜美女国产一区| 日韩在线高清观看一区二区三区| 中文字幕人妻丝袜制服| 男人操女人黄网站| 色哟哟·www| 久久精品久久精品一区二区三区| 国产 精品1| 18禁动态无遮挡网站| 久久久a久久爽久久v久久| 一区二区日韩欧美中文字幕| 国产精品一区二区在线观看99| 中文天堂在线官网| 国产精品秋霞免费鲁丝片| 大码成人一级视频| 免费观看性生交大片5| 欧美 亚洲 国产 日韩一| av不卡在线播放| 999精品在线视频| 狂野欧美激情性bbbbbb| 亚洲少妇的诱惑av| 欧美黄色片欧美黄色片| 人成视频在线观看免费观看| 中文字幕人妻熟女乱码| 国产片内射在线| 国产精品久久久久久av不卡| 激情视频va一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲三级黄色毛片| 国产精品.久久久| 日韩av不卡免费在线播放| 观看av在线不卡| 亚洲国产欧美在线一区| av在线播放精品| 成年人午夜在线观看视频| 黄色 视频免费看| 亚洲欧美色中文字幕在线| 国产毛片在线视频| 母亲3免费完整高清在线观看 | 国产熟女午夜一区二区三区| 欧美另类一区| 亚洲av免费高清在线观看| 亚洲精品中文字幕在线视频| 男人添女人高潮全过程视频| 国产男女超爽视频在线观看| 国产精品国产三级专区第一集| 日韩一卡2卡3卡4卡2021年| 欧美精品高潮呻吟av久久| 黄网站色视频无遮挡免费观看| 精品亚洲成a人片在线观看| 少妇被粗大的猛进出69影院| 欧美日韩综合久久久久久| 亚洲精品日韩在线中文字幕| 美女国产视频在线观看| 国产1区2区3区精品| 777久久人妻少妇嫩草av网站| 狠狠精品人妻久久久久久综合| 考比视频在线观看| 欧美最新免费一区二区三区| av免费观看日本| 精品一区在线观看国产| 亚洲综合精品二区| 999久久久国产精品视频| 精品人妻在线不人妻| 日韩熟女老妇一区二区性免费视频| 国产精品 国内视频| 日本wwww免费看| 久久精品久久久久久噜噜老黄| 大话2 男鬼变身卡| 午夜av观看不卡| 有码 亚洲区| 亚洲成人手机| 黄网站色视频无遮挡免费观看| 一级爰片在线观看| 人体艺术视频欧美日本| 大片电影免费在线观看免费| av女优亚洲男人天堂| 国产片特级美女逼逼视频| 国产在线视频一区二区| 精品少妇内射三级| 男的添女的下面高潮视频| 最近中文字幕2019免费版| 少妇被粗大的猛进出69影院| 宅男免费午夜| 考比视频在线观看| 欧美av亚洲av综合av国产av | 久久久久久人妻| 国产野战对白在线观看| 久久精品国产鲁丝片午夜精品| 国产在视频线精品| 18禁裸乳无遮挡动漫免费视频| 在线看a的网站| 亚洲精品久久久久久婷婷小说| 亚洲精品国产一区二区精华液| 免费av中文字幕在线| 日韩一卡2卡3卡4卡2021年| 中文天堂在线官网| 久久鲁丝午夜福利片| 老司机影院毛片| 亚洲一码二码三码区别大吗| 成年美女黄网站色视频大全免费| 国产成人精品久久二区二区91 | 久久精品国产亚洲av涩爱| 亚洲中文av在线| 亚洲欧美一区二区三区黑人 | 另类精品久久| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费观看性视频| 国产成人精品福利久久| 国产精品秋霞免费鲁丝片| 久久午夜福利片| 91成人精品电影| 熟女电影av网| 国产片内射在线| 青春草国产在线视频| 亚洲精品美女久久久久99蜜臀 | 欧美亚洲 丝袜 人妻 在线| 成人漫画全彩无遮挡| 婷婷成人精品国产| 熟女电影av网| 精品第一国产精品| 国产在线一区二区三区精| 欧美xxⅹ黑人| 亚洲情色 制服丝袜| 黄色 视频免费看| 国产成人精品福利久久| 国产男女内射视频| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜| 免费不卡的大黄色大毛片视频在线观看| 超碰成人久久| 深夜精品福利| 国产精品亚洲av一区麻豆 | 久久久久国产精品人妻一区二区| 女性被躁到高潮视频| 中文字幕av电影在线播放| 在线观看免费高清a一片| 亚洲精品国产av成人精品| 久久久久精品性色| 在线亚洲精品国产二区图片欧美| 天天影视国产精品| 99re6热这里在线精品视频| 亚洲第一区二区三区不卡| 免费人妻精品一区二区三区视频| 少妇人妻久久综合中文| 人人妻人人添人人爽欧美一区卜| 久久久久精品性色| 国产亚洲av片在线观看秒播厂| 亚洲国产看品久久| 青春草视频在线免费观看| 欧美黄色片欧美黄色片| 日韩欧美一区视频在线观看| 91精品三级在线观看| 大话2 男鬼变身卡| 观看av在线不卡| 中文精品一卡2卡3卡4更新| xxxhd国产人妻xxx| a级毛片黄视频| 成年av动漫网址| 人妻 亚洲 视频| 波野结衣二区三区在线| 丝袜人妻中文字幕| 亚洲综合色网址| 男男h啪啪无遮挡| 精品少妇一区二区三区视频日本电影 | 亚洲av中文av极速乱| 校园人妻丝袜中文字幕| 精品视频人人做人人爽| 99热国产这里只有精品6| 春色校园在线视频观看| 亚洲,欧美,日韩| 国产一区二区在线观看av| 中文乱码字字幕精品一区二区三区| 久久久欧美国产精品| 亚洲欧洲国产日韩| av免费观看日本| 91久久精品国产一区二区三区| 欧美亚洲日本最大视频资源| 永久免费av网站大全| 亚洲精品久久午夜乱码| 久久狼人影院| 色婷婷av一区二区三区视频| 建设人人有责人人尽责人人享有的| 成人免费观看视频高清| av一本久久久久| 久久久久久人妻| 伊人亚洲综合成人网| 日本av手机在线免费观看| 极品少妇高潮喷水抽搐| 免费观看a级毛片全部| 丝袜美足系列| 午夜福利,免费看| 国产成人免费观看mmmm| 丰满饥渴人妻一区二区三| videosex国产| 有码 亚洲区| 久热久热在线精品观看| 精品福利永久在线观看| 午夜日韩欧美国产| 午夜福利,免费看| 亚洲一级一片aⅴ在线观看| 亚洲三区欧美一区| 免费观看在线日韩| 久久精品久久精品一区二区三区| 国产色婷婷99| 又大又黄又爽视频免费| 久久人人爽人人片av| 国产白丝娇喘喷水9色精品| 国产精品欧美亚洲77777| 久久久久久人人人人人| 美女大奶头黄色视频| 黄色毛片三级朝国网站| 国产精品国产av在线观看| 视频区图区小说| 九色亚洲精品在线播放| 熟女电影av网| 久久韩国三级中文字幕| 欧美人与善性xxx| 一本色道久久久久久精品综合| 青春草国产在线视频| 熟女电影av网| 美女福利国产在线| 黄片小视频在线播放| 爱豆传媒免费全集在线观看| 午夜精品国产一区二区电影| 亚洲人成电影观看| 老司机影院成人| 91精品三级在线观看| 中文精品一卡2卡3卡4更新| 久久久久久人人人人人| 久久99蜜桃精品久久| av不卡在线播放| 高清视频免费观看一区二区| 午夜久久久在线观看| 久久人人97超碰香蕉20202| 日韩人妻精品一区2区三区| av卡一久久| 另类亚洲欧美激情| 熟女少妇亚洲综合色aaa.| 成年av动漫网址| 一级,二级,三级黄色视频| 精品国产国语对白av| 免费在线观看完整版高清| 性少妇av在线| 伊人久久大香线蕉亚洲五| 亚洲av欧美aⅴ国产| 午夜日韩欧美国产| 大香蕉久久网| 久久综合国产亚洲精品| 日韩精品免费视频一区二区三区| 欧美日韩一级在线毛片| 国产精品一二三区在线看| 免费观看在线日韩| 亚洲国产欧美网| 日韩视频在线欧美| 久久久久久伊人网av| 国产精品国产三级国产专区5o| 亚洲经典国产精华液单| 久久99热这里只频精品6学生| 亚洲第一av免费看| 日韩伦理黄色片| 亚洲五月色婷婷综合| 国产高清不卡午夜福利| 日本av手机在线免费观看| 日韩一区二区视频免费看| 极品人妻少妇av视频| 超碰97精品在线观看| 伊人亚洲综合成人网| 性色avwww在线观看| 亚洲国产日韩一区二区| 成人漫画全彩无遮挡| 青春草视频在线免费观看| 精品久久久久久电影网| 久久久久久人人人人人| 日本色播在线视频| 建设人人有责人人尽责人人享有的| 成人毛片60女人毛片免费| 国产成人a∨麻豆精品| 国产不卡av网站在线观看| 亚洲av福利一区| 十八禁高潮呻吟视频| 欧美日韩视频精品一区| 国产日韩欧美视频二区| 电影成人av| 99久久精品国产国产毛片| 性高湖久久久久久久久免费观看| 日韩成人av中文字幕在线观看| 午夜福利一区二区在线看| av免费在线看不卡| 99久久人妻综合| 91aial.com中文字幕在线观看| 下体分泌物呈黄色| 国产免费福利视频在线观看| 2021少妇久久久久久久久久久| 久久精品国产亚洲av涩爱| 国产视频首页在线观看| 久久这里有精品视频免费| 少妇的丰满在线观看| 18禁观看日本| 中文字幕人妻熟女乱码| 亚洲av综合色区一区| 亚洲精品日韩在线中文字幕| 国产 精品1| 欧美激情高清一区二区三区 | 精品99又大又爽又粗少妇毛片| 欧美日韩成人在线一区二区| 亚洲精品第二区| 久久久久久久久久久久大奶| 亚洲国产精品成人久久小说| 交换朋友夫妻互换小说| 99久久精品国产国产毛片| 欧美日韩综合久久久久久| 国产成人免费无遮挡视频| 国产成人一区二区在线| 9色porny在线观看| 人妻一区二区av| 激情五月婷婷亚洲| 赤兔流量卡办理| 男女边吃奶边做爰视频| 午夜福利一区二区在线看| 日本猛色少妇xxxxx猛交久久| 一个人免费看片子| 国产男人的电影天堂91| 欧美在线黄色| 久久久国产一区二区| 久久久久久久亚洲中文字幕| 中文字幕色久视频| 男女无遮挡免费网站观看| 国产精品人妻久久久影院| 精品亚洲成a人片在线观看| 久久这里只有精品19| www.自偷自拍.com| 国产成人免费观看mmmm| 新久久久久国产一级毛片| 国产国语露脸激情在线看| 999久久久国产精品视频| 最黄视频免费看| 免费看av在线观看网站| av福利片在线| 亚洲av在线观看美女高潮| 免费播放大片免费观看视频在线观看| 三上悠亚av全集在线观看| 少妇被粗大猛烈的视频| 狠狠精品人妻久久久久久综合| 日本色播在线视频| 久久久久久人妻| 成人国产麻豆网| 中文字幕制服av| 观看av在线不卡| 中文字幕亚洲精品专区| 免费播放大片免费观看视频在线观看| 久久精品国产鲁丝片午夜精品| 日韩中文字幕视频在线看片| 国产精品麻豆人妻色哟哟久久| 18禁裸乳无遮挡动漫免费视频| 波多野结衣av一区二区av| 桃花免费在线播放| 一级爰片在线观看| 大片电影免费在线观看免费| 久久毛片免费看一区二区三区| 国产视频首页在线观看| 欧美人与善性xxx| 亚洲精品久久午夜乱码| 成人亚洲精品一区在线观看| 午夜福利影视在线免费观看| 成人毛片60女人毛片免费| 亚洲国产精品成人久久小说| 这个男人来自地球电影免费观看 | 国产精品国产av在线观看| av天堂久久9| 久久av网站| 国产精品一区二区在线不卡| 久久av网站| 午夜免费鲁丝| 五月开心婷婷网| 美女福利国产在线| av.在线天堂| 我要看黄色一级片免费的| 大话2 男鬼变身卡| 精品少妇一区二区三区视频日本电影 | 老司机影院毛片| 午夜老司机福利剧场| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 国产成人精品久久久久久| 国产片特级美女逼逼视频| 精品人妻偷拍中文字幕| 精品国产超薄肉色丝袜足j| 亚洲天堂av无毛| 在线亚洲精品国产二区图片欧美| 国产成人91sexporn| 久久久国产欧美日韩av| 欧美中文综合在线视频| 亚洲精品一二三| 人人妻人人添人人爽欧美一区卜| 亚洲经典国产精华液单| 热99国产精品久久久久久7| 纵有疾风起免费观看全集完整版| 国产精品国产三级国产专区5o| 制服诱惑二区| 久久精品国产a三级三级三级| 色吧在线观看| 在线观看免费视频网站a站| 欧美精品av麻豆av| 久久久久国产网址| videossex国产| 制服丝袜香蕉在线| 一边摸一边做爽爽视频免费| 亚洲天堂av无毛| 亚洲美女搞黄在线观看| 国产精品亚洲av一区麻豆 | 精品人妻偷拍中文字幕| 超碰97精品在线观看| 国产欧美日韩综合在线一区二区| 新久久久久国产一级毛片| 午夜激情久久久久久久| 国产成人午夜福利电影在线观看| 一区二区三区四区激情视频| 国产有黄有色有爽视频| 一级a爱视频在线免费观看| 中文欧美无线码| 国产爽快片一区二区三区| 色吧在线观看| 国产精品国产av在线观看| 久久久久久久精品精品| 免费高清在线观看日韩| 国产精品嫩草影院av在线观看| 男女边摸边吃奶| 欧美亚洲日本最大视频资源| 女人精品久久久久毛片| 国产av一区二区精品久久| 亚洲国产精品一区三区| 大香蕉久久网| 色婷婷久久久亚洲欧美| 少妇人妻精品综合一区二区| 日本色播在线视频| 老女人水多毛片| 男女边吃奶边做爰视频| 日韩,欧美,国产一区二区三区| 91精品伊人久久大香线蕉| 黄频高清免费视频| 久久人人爽人人片av| 欧美日韩一级在线毛片| 久久国产精品男人的天堂亚洲| 中文字幕最新亚洲高清| 成人亚洲精品一区在线观看| 国产精品久久久久久精品电影小说| a 毛片基地| 久久99一区二区三区| 九色亚洲精品在线播放| 欧美bdsm另类| 巨乳人妻的诱惑在线观看| 久久午夜综合久久蜜桃| 亚洲av电影在线进入| 国产片特级美女逼逼视频| av在线观看视频网站免费| 久久久久久久大尺度免费视频| 极品人妻少妇av视频| 成年动漫av网址| 国产一区二区 视频在线| 女人久久www免费人成看片| 亚洲精品国产av成人精品| 亚洲激情五月婷婷啪啪| 男人舔女人的私密视频| 亚洲国产欧美在线一区| 欧美激情 高清一区二区三区| 久久97久久精品| 天堂中文最新版在线下载| 国产日韩欧美视频二区| 菩萨蛮人人尽说江南好唐韦庄| 美女午夜性视频免费| 久久人人爽av亚洲精品天堂| 亚洲人成77777在线视频| 欧美日韩av久久| 一区在线观看完整版| 伊人久久国产一区二区| 男女边吃奶边做爰视频| 五月开心婷婷网| 男人爽女人下面视频在线观看| 91国产中文字幕| av在线app专区| 精品一区在线观看国产| 五月开心婷婷网| a级毛片黄视频| 哪个播放器可以免费观看大片| 精品人妻一区二区三区麻豆| 丰满乱子伦码专区| 伦理电影大哥的女人| 亚洲国产色片| 欧美人与性动交α欧美精品济南到 | 久久久国产精品麻豆| 下体分泌物呈黄色| 99热国产这里只有精品6| 久久狼人影院| 免费观看无遮挡的男女| 成人毛片a级毛片在线播放| 亚洲av综合色区一区| 18禁裸乳无遮挡动漫免费视频| 欧美老熟妇乱子伦牲交| 人妻人人澡人人爽人人| 国产精品偷伦视频观看了| 国产野战对白在线观看| 9色porny在线观看| 免费观看av网站的网址| 亚洲av综合色区一区| 亚洲精品乱久久久久久| 91精品三级在线观看| 亚洲av免费高清在线观看| 黄片无遮挡物在线观看| 国产亚洲av片在线观看秒播厂| 亚洲av电影在线观看一区二区三区| 99精国产麻豆久久婷婷| 久久人人爽人人片av| 色婷婷av一区二区三区视频| 国产精品不卡视频一区二区| 国产在视频线精品| 制服诱惑二区| 青春草视频在线免费观看| 免费观看性生交大片5| 满18在线观看网站| 国产精品久久久久久精品古装| 日韩中文字幕视频在线看片| 老汉色av国产亚洲站长工具| a级毛片黄视频| 久久精品国产自在天天线| h视频一区二区三区| 久久精品久久精品一区二区三区| 国产女主播在线喷水免费视频网站| 久久久久久久大尺度免费视频| 2018国产大陆天天弄谢| 成年人免费黄色播放视频| www日本在线高清视频| 欧美日韩视频精品一区| 国产淫语在线视频| 精品少妇内射三级| 欧美亚洲 丝袜 人妻 在线| 日本91视频免费播放| 青春草亚洲视频在线观看| 一级毛片黄色毛片免费观看视频| 亚洲av综合色区一区| 一级片'在线观看视频| 五月开心婷婷网| 最近的中文字幕免费完整| 亚洲av成人精品一二三区| 免费看不卡的av| 久久精品国产综合久久久| freevideosex欧美| 波多野结衣av一区二区av| 校园人妻丝袜中文字幕| 成人毛片a级毛片在线播放| 日韩制服丝袜自拍偷拍| 国产人伦9x9x在线观看 | 麻豆乱淫一区二区| 爱豆传媒免费全集在线观看| av电影中文网址| 熟女少妇亚洲综合色aaa.| 亚洲 欧美一区二区三区| 80岁老熟妇乱子伦牲交| 男女无遮挡免费网站观看| 日本av免费视频播放| 久久精品夜色国产| 黄色一级大片看看| 日本猛色少妇xxxxx猛交久久| 97在线人人人人妻| 亚洲一级一片aⅴ在线观看| 免费黄网站久久成人精品| 一区二区av电影网| 婷婷色av中文字幕| 日产精品乱码卡一卡2卡三| 一级a爱视频在线免费观看| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频| 不卡视频在线观看欧美| 天堂8中文在线网| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美精品永久| 美女午夜性视频免费| 777米奇影视久久| 精品人妻一区二区三区麻豆| 精品一区二区免费观看| 免费黄网站久久成人精品| 久久久久精品久久久久真实原创| 欧美 日韩 精品 国产| 欧美人与性动交α欧美精品济南到 | 欧美激情极品国产一区二区三区| 校园人妻丝袜中文字幕| 亚洲欧美成人精品一区二区| 国产亚洲午夜精品一区二区久久| 精品国产一区二区三区四区第35| 热99久久久久精品小说推荐| 成人漫画全彩无遮挡| 国产一区二区激情短视频 | 色哟哟·www| 少妇的逼水好多| 久久ye,这里只有精品| 日韩电影二区|