劉小寧 劉 岑 陳 剛 劉 兵 楊 帆 張紅衛(wèi)

(1.武漢軟件工程職業(yè)學(xué)院 湖北 武漢：430205；2.湖北輕工職業(yè)技術(shù)學(xué)院 湖北 武漢：430070)

深冷容器是醫(yī)療、航空航天與能源等領(lǐng)域經(jīng)常使用的承壓設(shè)備，非預(yù)應(yīng)變與9%預(yù)應(yīng)變奧氏體不銹鋼S30408是制造深冷容器的常用鋼材[1]，探索奧氏體不銹鋼S30408超低溫機(jī)械性能參數(shù)的概率分布，即探索其分布規(guī)律與確定其分布參數(shù)是構(gòu)建深冷容器強(qiáng)度可靠性設(shè)計(jì)方法的基礎(chǔ)工作[2]。

奧氏體不銹鋼S30408超低溫的機(jī)械性能參數(shù)主要包括抗拉強(qiáng)度、屈服強(qiáng)度與屈強(qiáng)比(屈服強(qiáng)度與抗拉強(qiáng)度之比)，研究表明[3-4]，奧氏體不銹鋼S30408超低溫抗拉強(qiáng)度與屈強(qiáng)比，分別是基本符合正態(tài)分布的隨機(jī)變量。但是，目前對(duì)奧氏體不銹鋼S30408超低溫屈服強(qiáng)度概率分布的研究似不深入，有待于進(jìn)一步加強(qiáng)。

為了推動(dòng)深冷容器強(qiáng)度可靠性設(shè)計(jì)方法的建立，文中應(yīng)用概率論和數(shù)理統(tǒng)計(jì)知識(shí)[2,5-9]，基于非預(yù)應(yīng)變與9%預(yù)應(yīng)變奧氏體不銹鋼S30408的有限有效試驗(yàn)數(shù)據(jù)[1,10-12]，對(duì)其超低溫屈服強(qiáng)度的分布規(guī)律與分布參數(shù)進(jìn)行了探索。

1 基本理論與方法

1.1 分布規(guī)律的假設(shè)檢驗(yàn)

(1)計(jì)算有效試驗(yàn)數(shù)據(jù)的準(zhǔn)確度與精密度。當(dāng)通過(guò)試驗(yàn)測(cè)得奧氏體不銹鋼S30408超低溫屈服強(qiáng)度R的n組有效試驗(yàn)數(shù)據(jù)Ri(i=1，2，…，n)時(shí)，其準(zhǔn)確度與精密度分別為：

(1)

(2)

(2)假設(shè)超低溫屈服強(qiáng)度R的分布規(guī)律。假設(shè)超低溫屈服強(qiáng)度R是符合正態(tài)分布的隨機(jī)變量。

(3)預(yù)分組檢驗(yàn)。

①列表分組計(jì)算理論概率與理論頻數(shù)。根據(jù)有效數(shù)據(jù)個(gè)數(shù)n，把其為M個(gè)組：

M=1+3.3lgn

并取整數(shù)。對(duì)于符合正態(tài)分布的隨機(jī)變量R，其統(tǒng)計(jì)量Ri落在分組[a1，a2]，[a2，a3]，…，[aM，aM+1]內(nèi)的理論概率為：

(3)

式中，φ(·)為標(biāo)準(zhǔn)正態(tài)積分；a1=(Ri)min，aM+1=(Ri)max，(Ri)min、(Ri)max分別為Ri中的最小與最大值。

對(duì)于n個(gè)有效試驗(yàn)數(shù)據(jù)，Ri位于分組[aj，aj+1]內(nèi)的理論頻數(shù)為n×pj。

②列表試算分組的皮爾遜統(tǒng)計(jì)量并計(jì)算其之和。第j個(gè)分組的實(shí)際頻數(shù)Nj與理論頻數(shù)n×pj差異的皮爾遜統(tǒng)計(jì)量為：

(4)

M個(gè)分組實(shí)際頻數(shù)與理論頻數(shù)差異的皮爾遜統(tǒng)計(jì)量之和為：

(5)

其檢驗(yàn)條件為：若符合

(6)

則在顯著度為δ時(shí)，假設(shè)成立，否則存在分組不合理引起檢驗(yàn)不符合假設(shè)的可能，需要重新進(jìn)行合理分組，然后再一次檢驗(yàn)。

(4)重新分組檢驗(yàn)。

當(dāng)重新分組后的組數(shù)為m，每組的實(shí)際頻數(shù)與理論頻數(shù)分別為Ni與n×pi，第i個(gè)分組的實(shí)際頻數(shù)Ni與理論頻數(shù)n×pi差異的皮爾遜統(tǒng)計(jì)量為：

(8)

m個(gè)分組實(shí)際頻數(shù)與理論頻數(shù)差異的皮爾遜統(tǒng)計(jì)量之和為：

(9)

(10)

則在顯著度為δ時(shí)，假設(shè)成立，否則在現(xiàn)有條件下假設(shè)不成立。本文取顯著度δ=0.05，所用的χ2系數(shù)見(jiàn)表1[2,5-6]。

表1 χ2與t系數(shù)

1.2 分布參數(shù)的取值區(qū)間

若奧氏體不銹鋼S30408的超低溫屈服強(qiáng)度基本符合正態(tài)分布，當(dāng)雙側(cè)置信度為(1-α)時(shí)，其分布參數(shù)的取值區(qū)間按如下方法計(jì)算。

1.2.1 均值

屈服強(qiáng)度均值的取值區(qū)間為：

μ∈[μmin，μmax]

(11)

其中

(12)

(13)

式中，μ為屈服強(qiáng)度R的均值；μmin與μmax分別為μ在雙側(cè)置信度為(1-α)時(shí)的下限與上限；tn-1, 1-0.5α為t分布系數(shù)，由自由度(n-1)與0.5α查得。

α=0.02時(shí)，文中所用的t分布系數(shù)值見(jiàn)表1[2,5-6]。

1.2.2 標(biāo)準(zhǔn)差

屈服強(qiáng)度R標(biāo)準(zhǔn)差σ的取值范圍為：

σ∈[σmin，σmax]

(14)

其中

(15)

(16)

α=0.02時(shí)，文中所用的χ2分布系數(shù)值見(jiàn)表1[2,5-6]。

1.2.3 變異系數(shù)

屈服強(qiáng)度R變異系數(shù)C的取值范圍為：

C∈[Cmin，Cmax]

(17)

其中

Cmin=σmin/μmax，Cmax=σmax/μmin

(18)

式中，C為材料性能R的變異系數(shù)；Cmin與Cmax分別為C在雙側(cè)置信度為(1-α)時(shí)的下限與上限。

2 分布規(guī)律的假設(shè)檢驗(yàn)

深冷容器常用非預(yù)應(yīng)變與9%預(yù)應(yīng)變奧氏體不銹鋼S30408制造[1]，下面分別討論其超低溫屈服強(qiáng)度的分布規(guī)律。

2.1 9%預(yù)應(yīng)變奧氏體不銹鋼的超低溫屈服強(qiáng)度

在超低溫(液氮溫度即-196℃)時(shí)，文獻(xiàn)[1]獲得了9%預(yù)應(yīng)變奧氏體不銹鋼S30408屈服強(qiáng)度的43組試驗(yàn)數(shù)據(jù)，在雙側(cè)置信度為99%時(shí)，經(jīng)用文獻(xiàn)[10-12]的方法判別，其中42組是有效的，見(jiàn)表2。

將42組有效試驗(yàn)數(shù)據(jù)代入式(1)與式(2)中，可得到其平均值與精密度，一并列入表2。

假設(shè)9%預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度是基本符合正態(tài)分布的隨機(jī)變量。

表2 9%預(yù)應(yīng)變S30408鋼超低溫屈服強(qiáng)度的有效試驗(yàn)數(shù)據(jù)

表3 9%預(yù)應(yīng)變S30408鋼超低溫屈服強(qiáng)度預(yù)分組檢驗(yàn)的皮爾遜統(tǒng)計(jì)量

根據(jù)以上分析，可先將第3組與第4組、第5組與第6組分別合并，然后重新分組檢驗(yàn)。

重新分組后的組數(shù)m=4，第i個(gè)分組的實(shí)際頻數(shù)Ni、理論概率pi、理論頻數(shù)n×pi、以及反映每個(gè)分組實(shí)際頻數(shù)與理論頻數(shù)差異的皮爾遜統(tǒng)計(jì)量及其之和見(jiàn)表4。

4個(gè)組的自由度為f=4-1-2=1，取顯著度δ=0.05，

顯然，預(yù)分組檢驗(yàn)是重新分組檢驗(yàn)基礎(chǔ)，重新分組檢驗(yàn)是預(yù)分組檢驗(yàn)的補(bǔ)充與完善。

2.2 非預(yù)應(yīng)變奧氏體不銹鋼的超低溫屈服強(qiáng)度

在超低溫時(shí)，文獻(xiàn)[1]通過(guò)試驗(yàn)，也獲得了非預(yù)應(yīng)變奧氏體不銹鋼S30408屈服強(qiáng)度的60組試驗(yàn)數(shù)據(jù)，在雙側(cè)置信度為99%時(shí)，經(jīng)用文獻(xiàn)[10-12]建立的方法判別，這些數(shù)據(jù)都是有效的，見(jiàn)表5。將60

表4 9%預(yù)應(yīng)變S30408鋼超低溫屈服強(qiáng)度重新分組檢驗(yàn)的皮爾遜統(tǒng)計(jì)量

組有效試驗(yàn)數(shù)據(jù)代入式(1)與式(2)中，可得到其平均值與精密度，一并列入表5。

表5 非預(yù)應(yīng)變S30408鋼超低溫屈服強(qiáng)度的有效試驗(yàn)數(shù)據(jù)

采用與上面相同的方法，在顯著度為0.05時(shí)可以證明，非預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度是基本符合正態(tài)分布的隨機(jī)變量。

3 分布參數(shù)的取值區(qū)間

根據(jù)以上分析，在超低溫時(shí)，非預(yù)應(yīng)變與9%預(yù)應(yīng)變奧氏體不銹鋼S30408屈服強(qiáng)度分別是符合正態(tài)分布的隨機(jī)變量，因此，在雙側(cè)置信度為98%時(shí)，可根據(jù)式(11)～式(18)分別計(jì)算其分布參數(shù)的取值區(qū)間。

3.1 9%預(yù)應(yīng)變奧氏體不銹鋼的超低溫屈服強(qiáng)度

將表2中有效試驗(yàn)數(shù)據(jù)的平均值與精密度代入式(11)～式(18)，可得到9%預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度分布參數(shù)的取值區(qū)間。

均值μ的取值區(qū)間為：

μ∈[546MPa，583MPa]

(19)

標(biāo)準(zhǔn)差σ的取值范圍為：

σ∈[39.50MPa，66.52MPa]

(20)

變異系數(shù)C取值范圍為：

C∈[0.0678，0.1218]

(21)

3.2 非預(yù)應(yīng)變奧氏體不銹鋼的超低溫屈服強(qiáng)度

將表5中有效試驗(yàn)數(shù)據(jù)的平均值與精密度代入式(11)～式(18)，可得到非預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度分布參數(shù)的取值區(qū)間。

均值μ的取值區(qū)間為：

μ∈[487MPa，526MPa]

(22)

標(biāo)準(zhǔn)差σ的取值范圍為：

σ∈[51.42MPa，79.25MPa]

(23)

變異系數(shù)C的取值范圍為：

C∈[0.0978，0.1627]

(24)

3.3 討論

相對(duì)于非預(yù)應(yīng)變而言，9%預(yù)應(yīng)變是在深冷容器制造前對(duì)奧氏體不銹鋼S30408的一種處理工藝。根據(jù)文獻(xiàn)[14-15]建立的兩個(gè)正態(tài)隨機(jī)變量分布參數(shù)的比較與評(píng)價(jià)方法，分別比較式(19)與式(22)、式(21)與式(24)可知，在超低溫時(shí)，9%預(yù)應(yīng)變是顯著提高奧氏體不銹鋼S30408屈服強(qiáng)度均值的有效方法，盡管其變異系數(shù)(精度)沒(méi)有顯著變化，但存在變異系數(shù)減小即精度提高的趨勢(shì)。

4 結(jié)語(yǔ)

應(yīng)用概率論與數(shù)理統(tǒng)計(jì)方法，對(duì)分布規(guī)律的假設(shè)檢驗(yàn)方法進(jìn)行了補(bǔ)充完善。在雙側(cè)置信度為99%時(shí)，分析了試驗(yàn)數(shù)據(jù)的有效性；基于有限的有效試驗(yàn)數(shù)據(jù)及98%的雙側(cè)置信度，探索了非預(yù)應(yīng)變與9%預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度的分布規(guī)律與分布參數(shù)。

1)顯著度為0.05時(shí)，非預(yù)應(yīng)變與9%預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度基本符合正態(tài)分布。

2)雙側(cè)置信度為98%時(shí)，9%預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度的均值位于546MPa與583MPa之間，標(biāo)準(zhǔn)差位于39.50MPa與66.52MPa之間，變異系數(shù)位于0.0678與0.1218之間。

3)雙側(cè)置信度為98%時(shí)，非預(yù)應(yīng)變奧氏體不銹鋼S30408超低溫屈服強(qiáng)度的均值位于487MPa與526MPa之間，標(biāo)準(zhǔn)差位于51.42MPa與79.25MPa之間，變異系數(shù)位于0.0978與0.1627之間。

4)在變異系數(shù)(精度)沒(méi)有顯著變化時(shí)，9%預(yù)應(yīng)變是顯著提高奧氏體不銹鋼S30408超低溫屈服強(qiáng)度均值的有效方法。