陸曉珩,高延兵,瞿建農(nóng),任衛(wèi)禮
(1.上海市嘉定區(qū)計(jì)量質(zhì)量檢測(cè)所,上海 201818;2.新大洲本田摩托有限公司,上海 201708)
工程中的大型鋼筋水泥預(yù)制件在澆注前,內(nèi)置鋼筋線索的張拉力,是由電動(dòng)高壓油泵工作時(shí)產(chǎn)生的壓力(單位MPa),同步帶動(dòng)千斤頂緩慢工作起到收緊作用(產(chǎn)生的力值單位kN)。然而,這兩者之間的關(guān)系僅僅按照某些特定關(guān)鍵檢測(cè)點(diǎn)的測(cè)試數(shù)據(jù)是不能滿足生產(chǎn)工藝要求的,必須進(jìn)行適當(dāng)?shù)募夹g(shù)處理。
用數(shù)學(xué)上的回歸分析數(shù)理統(tǒng)計(jì)方法,對(duì)一系列測(cè)量數(shù)據(jù)進(jìn)行處理,可以得出比較符合事物內(nèi)部規(guī)律的數(shù)學(xué)表達(dá)式,能夠比較科學(xué)客觀地解決各變量之間的內(nèi)在關(guān)系。
在澆注大型鋼筋水泥預(yù)制件的實(shí)際工作中,操作人員根據(jù)計(jì)量技術(shù)機(jī)構(gòu)提供的檢測(cè)報(bào)告中的回歸方程式,事先準(zhǔn)備好對(duì)照表,以控制電動(dòng)高壓油泵壓力表產(chǎn)生的壓力值來決定千斤頂?shù)睦Α?/p>
依據(jù)文獻(xiàn)[1],將千斤頂置于經(jīng)檢定符合0.3級(jí)的材料試驗(yàn)機(jī)工作臺(tái)上,在保證其兩者的加力軸線基本重合的前提下,同時(shí)啟動(dòng)電動(dòng)高壓油泵和試驗(yàn)機(jī)。當(dāng)緩慢地給試驗(yàn)機(jī)加力后,千斤頂隨即同步帶動(dòng)電動(dòng)高壓油泵上的壓力表。至測(cè)試點(diǎn)時(shí),由兩人同時(shí)分別讀取試驗(yàn)機(jī)(其本身的示值誤差應(yīng)予修正)和電動(dòng)高壓油泵上壓力表的數(shù)據(jù)。某一組特定的關(guān)鍵檢測(cè)點(diǎn)測(cè)量數(shù)據(jù)如表1所示。
表1 一組特定關(guān)鍵檢測(cè)點(diǎn)測(cè)量數(shù)據(jù)
從表1數(shù)據(jù)可以直觀地發(fā)現(xiàn),兩個(gè)測(cè)量參數(shù)之間基本成線性關(guān)系。參考文獻(xiàn)[2],令一元線性回歸方程為
為便于數(shù)據(jù)處理,兩參數(shù)之間的關(guān)系以表2的形式表達(dá)。
表2 兩測(cè)量參數(shù)關(guān)系對(duì)應(yīng)
式(1)中的兩系數(shù)a和b分別為
于是系數(shù)取值有:
即所求的一元線性回歸方程為
如精確到千分位,則為
該一元線性回歸方程的可靠性,可以由以下兩種方法進(jìn)行檢驗(yàn)。
通常采用F檢驗(yàn)法。為此應(yīng)計(jì)算如下參數(shù):
回歸平方和 U=b×lxy=0.175 834 427×10 920.38=1 920.178 767,剩余平方和 Q=lyy-b×lxy=1 920.624-1920.178767=0.4452326。
在回歸方程中,每一個(gè)平方和都有一個(gè)與其對(duì)稱的自由度。在一元線性回歸中,總的離差平方和的自由度 υS=υU+υQ=n-1。而 υU的自由度對(duì)應(yīng)于自變量的個(gè)數(shù),在一元線性回歸中υU=1,自然υU=n-2。
查F分布表[2],當(dāng)置信概率為99%時(shí),F(xiàn)0.01(1,8)=11.26;
F=34502>>F0.01(1,8)=11.26,可以確定該一元線性回歸方程高度顯著,且可信賴程度為99%以上。
相關(guān)指數(shù)
于是相關(guān)系數(shù)
相關(guān)系數(shù)R是一個(gè)≤1的常數(shù),當(dāng)它越接近1時(shí),說明回歸直線的效果越好、可靠性越大?,F(xiàn)R=0.999884085,說明其可靠性是值得信賴的。
據(jù)此回歸方程式y(tǒng)=0.231470506+0.175834427x,可以非常容易地計(jì)算出任意千斤頂?shù)睦λ鶎?duì)應(yīng)的電動(dòng)高壓油泵壓力表的壓力值。比如:按圖紙?jiān)O(shè)計(jì)的要求,鋼筋線索的張拉力應(yīng)為200 kN,其生產(chǎn)工藝上的諸多要求力值所對(duì)應(yīng)的壓力表的壓力值如表3所示。
表3 鋼筋線索生產(chǎn)工藝上的壓力值
壓力表的讀數(shù)按現(xiàn)行有效的文獻(xiàn)[3]之5.3.3.2“壓力表的示值應(yīng)按分度值的1/5估讀”要求,該規(guī)格壓力表的讀數(shù)應(yīng)是十分位的偶數(shù)倍。于是,最后一行中“壓力表壓力值”中的數(shù)據(jù)按特定的法則決定進(jìn)位還是舍去。比如:
3.748 ÷2=1.874,按正常的“四舍五入”法則1.874≈1.9,再將 1.9×2=3.8,故 3.748=3.8;
7.265 ÷2=3.6325,按正常的“四舍五入”法則3.6325≈3.6,再將 3.6×2=7.2,故 7.265=7.2…。
如果將普通壓力表改用精密壓力表,由于后者按分度值的1/10估讀,可知測(cè)量數(shù)據(jù)的讀數(shù)精度必定有所提高。由此,回歸方程的準(zhǔn)確性、可靠性和表3中的對(duì)應(yīng)值必將更加完美。
應(yīng)用一元線性回歸數(shù)據(jù)處理方法,可順利解決兩個(gè)參數(shù)之間的關(guān)系問題,且具有準(zhǔn)確、可靠、靈活應(yīng)用等諸多優(yōu)點(diǎn)。
[1]JJG 621—2005液壓千斤頂檢定規(guī)程[S].北京:中國計(jì)量出版社,2005.
[2]費(fèi)業(yè)泰.誤差理論與數(shù)據(jù)處理[M].北京:機(jī)械工業(yè)出版社,1983.
[3]JJG 52—1999彈簧管式一般壓力表、壓力真空表和真空表檢定規(guī)程[S].北京:中國計(jì)量出版社,1999.