基于Matlab的碳纖維復合材料的B基準值算法

董本正+李偉+高麗紅+郭朝龍+査萌+吉愛紅+戴振東

摘 要： B基準值是評價碳纖維復合材料性能的重要參數，對B基準值的高效精確的計算有助于評價復合材料的性能。根據航空工業標準HB7618?2009，開發出基于Matlab的 GUI模塊的B基準值的算法，該程序可以實現樣本母體檢驗、異常數據檢查、統計分布檢驗、離散系數修正與檢驗、B基準值的計算等功能，最后以Word格式導出計算結果。與標準提供范例提供的計算結果進行對比，證明程序運行準確可靠，并且其界面操作簡單，顯著地提高了B基準值計算效率。

關鍵詞： 碳纖維增強復合材料； B基準值； 統計學分析； Matlab

中圖分類號： TP319 文獻標識碼： A 文章編號： 1004?373X（2014）16?0087?05

Matlab?based algorithm for B reference value of carbon fiber reinforced polymer

DONG Ben?zheng1， LI Wei1， GAO Li?hong2， GUO Chao?long1， ZHA Meng2， JI Ai?hong1， DAI Zhen?dong1

（1. Institute of Bio?inspired Structure and Surface Engineering， NUAA， Nanjing 210016， China； 2. The First Aircraft Institute， AVIC， Xian 710089， China）

Abstract： B reference value is an important parameter to evaluate carbon fiber reinforced polymer （CRFRP）， thus the efficient and accurate computation for B reference value of CFRP is helpful to evaluate its properties. According to the aviation industry standard of HB7618?2009， the B reference value algorithm based on Matlab GUI module is developed， which can realize the functions of specimen inspection， abnormal data check， statistical distribution test， discrete coefficient correction and test， B reference value calculation. The calculated result of B reference value is exported in the Word format. In comparison with calculation results from the Aviation Industry Standard of HB7618?2009， the algorithm is accurate and reliable. It improved the calculation efficient of B reference value significantly.

Keywords： carbon fiber reinforced polymer； B reference value； statistical analysis； Matlab

0 引 言

碳纖維復合材料（Carbon Fiber Reinforced Polymer，CFRP）是以樹脂為基體，碳纖維為增強體的復合材料，作為一種先進的復合材料具有高比強度、高比模量、抗疲勞性等優點，正是這些獨特的特性使CFRP成為與鋁合金、鈦合金、鋼合金并列的四大航空材料之一[1?3]。飛機在服役的過程中會遇到各種各樣極端情況，例如低溫、高溫、潮濕環境，這都會影響CFRP的性能[4?6]。航空工業一般采用B基準值來評價CFRP的性能。為了得到CFRP在這些環境下的B基準值，就需要進行大量的實驗；根據航空工業標準HB7618?2009《聚合基復合材料性能數據表達準則》提供的計算步驟 [7]， B基準值的計算方法比較繁瑣。因而面對大量的實驗數據如何高效準確地計算B基準值對于實驗人員來說是急需解決的問題?，F在對于B基準值計算方法仍然是借助Excel的手動計算為主，這種方式無疑計算耗時長，并且容易出錯，無法滿足快速處理大量實驗數據的要求。為了提高B基準值的計算效率，本文根據Matlab在矩陣運算、圖像處理和數值處理效率方面的優勢[8]，將B基準值的計算和Matlab相結合，利用其GUI（Graphical User Interface）模塊開發B基準值計算算法。

1 基準值計算原理

1.1 B基準值

為了保證材料的結構的可靠性和安全性，就需要用到一些高級的統計學方法來確定其合理的設計許用值，一般采用B基準值[7]。B基準值是建立在統計學上的衡量材料性能的參數，在95%的置信度下，90%性能數值群的值高于此值[9]。航空工業標準HB7618?2009關于B基準值的流程復雜，本文將著重從批次相容性檢驗、檢查異常數據、分布函數檢驗、離散系數修正與檢驗、B基準值計算及導出結果幾方面進行說明。

1.2 主要計算步驟

（1） 檢驗不同批次的實驗數據是否來自同一母體；檢驗方法是K樣本的Anderson ?Darling檢驗， 其統計公式如下[7]：

[ADK=n-1n2（k-1）i=1k1nij=1lhjnFij-njHjHj（n-Hj）-nhj4] （1）

如果ADK小于錯判臨界值ADC（5%錯判風險），則不同批次的實驗數據來自不同母體。如果相容性檢驗不通過的，為了增加樣本數據來自同一母體的可能性，采用1%錯判風險繼續判斷是否來自同一母體。

（1） 檢查異常數據。采用最大賦范殘差的方法定量的檢查異常數據：

[MNR=maxxi-xs] （2）

（2） 分布函數檢驗。根據經驗威布爾計算的B基值過于保守，因此優先檢驗樣本數據是否符合正態分布，如果不符合正態分布，再依次檢驗樣本數據是否服從威布爾分布和對數正態分布。

（3） 離散系數修正與檢驗：當離散系數低于4%時，容易忽略材料差異、實驗方法等方面的實際差異性，因此要對離散系數進行修正[10]，修正方法為：

[x*i=xi-α（xi）Δ] （3）

離散系數修正結束后，對不同環境下的樣本數據進行離散系數檢驗，離散系數的檢驗公式為：

[F=i=1k（wi-w）（k-1）i=1kj=1niwij-wi（n-k）] （4）

（4） 將不同環境下的數據歸一化，根據歸一化后的實驗數據的平均值和標準差計算B基準值。求出對應的折減系數，再與不同環境下樣本的平均值相乘，從而得到B基準值，折減系數的公式為：

[Bj=x-（kj）*s] （5）

2 程序實現

2.1 基準值算法的漸進流程圖和程序界面

圖1是根據航標HB7618?2009給出的計算步驟編制的程序流程圖，該程序流程圖包括了批次相容性檢驗、檢查異常數據、分布函數檢驗、離散系數修正與檢驗、B基準值計算及導出結果等功能。

圖1 B基準值算法的漸進流程圖

2.2 變量傳遞

B基準值的算法程序包含幾個不同的模塊（見圖1），而每一個模塊又包含了多個程序和子程序。為了實現變量在不同模塊和函數間進行傳遞，本文采用handles結構體和guidata函數來實現這個功能。在B基準值的算法程序界面中的所有控件使用同一個handles結構體，handles結構體中不僅保存了圖形窗口中所有對象的句柄，而且還可以獲取或設置某個對象的屬性。在程序的計算過程中需要將每一個模塊計算后的結果賦給handles結構體，以便于將這個變量傳遞到下一個計算模塊。而每次handles結構體添加了新的元素，guidata函數便會更新handles結構體，使handle結構體中的新添加元素可以傳遞到下一個模塊。guidata函數用法guidata（hObject，handles），其中，hObject是執行回調的控件對象的句柄，handle為結構體。handles結構體和guidata函數的用法如下所示：

function pushbutton13_Bvalue_Callback（hObject， eventdata， handles）

…

handles.Bvalue（j）=roundn（mean（nn_data（：））*Bj（j），?3）；

…

guidata（hObject，handles）

2.3 設置實驗環境

B基準值的算法界面程序首先要檢驗相同環境下不同批次的樣本數據是否來自同一母體。程序界面中包括四種實驗環境：低溫干態（CTD）、常溫干態（RTD）、高溫干態（ETD）、高溫濕態（ETW）， 為了避免母體檢驗的過程中不同環境下的樣本數據相互干擾，需要對程序界面中代表著不同實驗環境的復選框（Check Box）進行設置。當其中的一種環境的復選框處于選中狀態時，其他的三種環境的復選框需要設置成未選中狀態，這樣就可以避免相互干擾。下面是當CTD處于選中狀態，其他的三種實驗環境處于未選中狀態時的代碼：

if get（handles.checkbox1_CTD，′Value′）

set（handles.checkbox2_RTD，′value′，0）；

set（handles.checkbox3_ETD，′value′，0）；

set（handles.checkbox4_ETW，′value′，0）；

set（handles.checkbox5_hebing，′value′，0）；

end

當復選框處于選中狀態時，其Value值等于1；當復選框處于未選中狀態時，其Value值等于0。

2.4 分布函數檢驗

對于樣本數據的分布函數檢驗，本文采用Matlab本身的內置函數normplot（）和weibplot（）函數進行判斷。若樣本數據服從正態分布和威布爾分布，樣本數據點分布成一條直線，否則成一條曲線。至于樣本數據的對數正態分布檢驗，程序首先對樣本數據取對數將對數分布轉化為正態分布后，再用normplot（）函數檢驗是否服從正態分布。

在B基準值的計算過程中首先檢驗的是正態分布，可見樣本數據的正態分布對于B基準值計算的重要性。為了增加正態分布檢驗的準確性，本文除了采用normplot（）函數對樣本數據進行直觀的目視正態分布檢驗外，還借助函數H=jbtest（）對樣本數據進行進一步的正態分布檢驗。若H=0，則認X服從正態分布；若H=1，則否定X服從正態分布。對CTD環境下樣本數據進行正態分布檢驗的程序代碼如下，檢測結果如圖2所示。

…

normplot（handles.tdata（：））； %正態分布概率圖

handles.H=jbtest（handles.tdata（：））； %擬合優度檢驗

…

圖2 CTD環境下樣本數據進行正態分布檢驗圖

圖2中CTD環境下的樣本數據點分布大致成一條直線，并且檢驗結果H=0，表示接受正態分布假設?？梢奀TD樣本數據服從正態分布。

2.5 生成實驗報告

面對大量的實驗數據，書寫報告往往是重復性勞動工作，此時可以利用Matlab生成一表格模板，每次自動導入數據生成報告，這就可以節約大量的時間。Matlab作為一種面向對象的語言，它支持COM技術，利用它的COM編譯器可以把Matlab開發的程序轉換成方便使用的COM組件。而ActiveX控件技術是建立在COM技術之上，由微軟公司推出的共享程序數據和功能的技術。因此，Matlab可以利用ActiveX控件技術，在Matlab中調用Microsoft Word 程序插入表格，此時的Word是組件，是服務程序，而Matlab是控制器程序。下面是生成報告的部分程序代碼：

…

try

Word = actxGetRunningServer（′Word.Application′）；

%若Word服務器已經打開，返回其句柄Word

catch

Word = actxserver（′Word.Application′）；

%創建一個Microsoft Word服務器，返回句柄Word

…

Tables=Document.Tables.Add（Selection.Range，21，5）；

%插入一個21行5列的表格

…

Document.Save %保存文檔3 計算實例

本文根據航標HB7618?2009提供的范例數據來詳細說明B基準值的計算流程，并驗證程序界面計算的可靠性，如圖3所示。范例包含了所有實驗環境下的樣本數據，見表1。

圖3 B基準值程序界面

（1） 通過Anderson?Darling方法檢驗不同環境下的不同批次的實驗數據是否來自同一母體（依次勾選對應的實驗環境），檢驗后發現ETW（高溫濕態）環境下的三批樣本數據的ADK大于臨界值（ADK=2.258 2>ADC=1.917 4），可見該環境下三批數據來自不同母體，如圖4（a）所示。如果放大為1%錯判風險下， ADK小于臨界值（ADK=2.258 2

（2） 相容性檢驗檢驗結束后，根據公式（2）計算每一試樣的MNR值（單擊異常數據檢查），如果MNR值大于臨界值，那么就可以判斷該試樣為異常值。對于能夠反映材料、工藝參數和實驗環境的差異性的異常數據需要保留。而當異常數據是由不合格的試樣、實驗設備和夾具的缺陷造成時，異常數據要刪除（單擊刪除異常數據，見圖3），否則就會導致統計結果的失真。

圖4 母體檢驗

（3） 對不同環境下實驗數據進行正態分布假設檢驗，如果不符合正態分布再依次檢驗是否服從威布爾分布和對數正態分布。

（4） 檢驗相同環境下樣本數據的離散系數是否小于4%（單擊離散系數）。由于較低的離散系數容易忽略材料和實驗方法等方面的差異性，因此對于小于4%的情況，程序會根據式（3）將離散系數修正為4%。修正后的數據平均值保持不變，原來小于平均值的數據變得更小，原來大于平均值的數據變得更大（見圖5）。然后，再檢查不同環境下離散系數是否相等，根據公式（4）計算F值小于臨界值（F=2.38<2.699），由此可以判斷不同環境下的離散系數相等。

圖5 離散系數修正圖

（5） 合并不同環境下的實驗數據再一次進行正態分布檢驗（勾選合并，見圖3）。在第（3）步中已經驗證每種環境下的實驗數據的是否符合正態分布，如圖6（a）所示，由圖6（b）可見合并后的實驗數據依然符合正態分布。然后，用每種環境下的樣本均值對數據進行歸一化，求出歸一化后合并數據的平均值、標準差和離散系數見表1。

根據公式（5）可以求出相應的環境下的折減系數Bj，并與不同環境下的樣本數據平均值相乘，從而獲得B基準值。

圖6 數據的正態分布檢驗

為了驗證程序的可靠性，采用標準中提供的范例中的計算結果進行對比說明。表2中列出了在四種不同實驗環境下的程序計算結果和航標HB7618?2009提供的結果。從表中的數據對比可以發現誤差在+0.3%以內，誤差主要由容限系數和折減系數的近似計算造成?？梢姵绦蛴嬎愕慕Y果符合工程精度（誤差±5%）的要求，此外相比于借助Excel的手動計算法方法，效率顯著提高。

4 結 語

本文根據航空工業標準HB7618?2009給出的B基準值的計算步驟，將Matlab的GUI模塊與B基準值的計算相結合開發出B基準值的計算程序。相比于借助Excle的手動計算方法，該程序界面具有以下優點：

（1） 程序可以判斷樣本是否來自同一母體，還具有異常數據處理、統計分布檢驗、離散系數檢驗與修正等功能，最后實驗結果以Word格式導出。

（2） 使用該程序進行B基準值的計算方便快捷、效率高，并且具有良好的人機界面，顯著提高了B基準值的計算效率。

表2 不同環境下計算結果對比

參考文獻

[1] 張駿華，盛祖銘，孫繼同.復合材料結構設計指南[M].北京：北京宇航出版社，1999.

[2] 李威，郭權鋒.碳纖維復合材料在航天領域的應用[J].中國光學，2011，4（3）：201?212.

[3] 陳紹杰.復合材料設計手冊[M].北京：航空工業出版社，1990.

[4] 劉建華，曹東，張曉云，等.樹脂基復合材料T300/5405的吸濕性能及濕熱環境對力學性能的影響[J].航空材料學報，2010，30（4）：75?80.

[5] 呂小軍，張琦，馬兆慶，等.濕熱老化對碳纖維/環氧樹脂基復合材料力學性能影響研究[J].材料工程，2005（11）：50?53.

[6] 鄭錫濤，李野，劉海燕，等.濕熱譜老化對復合材料層壓板強度的影響[J].航空學報，1998，19（4）：462?465.

[7] 中國航空工業總公司.HB 7618?2009 聚合物基復合材料力學性能數據表達準則[S].北京：中國航空工業總公司，2009：50?85.

[8] 陳子為.基于 Matlab GUI 掃雷游戲的設計與實現[J].現代電子技術，2009，32（24）：85?88.

[9] FAA. DOT/FAA/AR?03/19： material qualification and equivalency for polymer matrix composite material systems： updated procedure [S]. USA： FAA， 2003：55?78.

[10] 孫堅石，葉強.復合材料力學性能數據B基準值計算程序[J]. 航空制造技術，2009（z1）：19?21.

…

圖2 CTD環境下樣本數據進行正態分布檢驗圖

圖2中CTD環境下的樣本數據點分布大致成一條直線，并且檢驗結果H=0，表示接受正態分布假設。可見CTD樣本數據服從正態分布。

2.5 生成實驗報告

面對大量的實驗數據，書寫報告往往是重復性勞動工作，此時可以利用Matlab生成一表格模板，每次自動導入數據生成報告，這就可以節約大量的時間。Matlab作為一種面向對象的語言，它支持COM技術，利用它的COM編譯器可以把Matlab開發的程序轉換成方便使用的COM組件。而ActiveX控件技術是建立在COM技術之上，由微軟公司推出的共享程序數據和功能的技術。因此，Matlab可以利用ActiveX控件技術，在Matlab中調用Microsoft Word 程序插入表格，此時的Word是組件，是服務程序，而Matlab是控制器程序。下面是生成報告的部分程序代碼：

…

try

Word = actxGetRunningServer（′Word.Application′）；

%若Word服務器已經打開，返回其句柄Word

catch

Word = actxserver（′Word.Application′）；

%創建一個Microsoft Word服務器，返回句柄Word

…

Tables=Document.Tables.Add（Selection.Range，21，5）；

%插入一個21行5列的表格

…

Document.Save %保存文檔3 計算實例

本文根據航標HB7618?2009提供的范例數據來詳細說明B基準值的計算流程，并驗證程序界面計算的可靠性，如圖3所示。范例包含了所有實驗環境下的樣本數據，見表1。

圖3 B基準值程序界面

（1） 通過Anderson?Darling方法檢驗不同環境下的不同批次的實驗數據是否來自同一母體（依次勾選對應的實驗環境），檢驗后發現ETW（高溫濕態）環境下的三批樣本數據的ADK大于臨界值（ADK=2.258 2>ADC=1.917 4），可見該環境下三批數據來自不同母體，如圖4（a）所示。如果放大為1%錯判風險下， ADK小于臨界值（ADK=2.258 2

（2） 相容性檢驗檢驗結束后，根據公式（2）計算每一試樣的MNR值（單擊異常數據檢查），如果MNR值大于臨界值，那么就可以判斷該試樣為異常值。對于能夠反映材料、工藝參數和實驗環境的差異性的異常數據需要保留。而當異常數據是由不合格的試樣、實驗設備和夾具的缺陷造成時，異常數據要刪除（單擊刪除異常數據，見圖3），否則就會導致統計結果的失真。

圖4 母體檢驗

（3） 對不同環境下實驗數據進行正態分布假設檢驗，如果不符合正態分布再依次檢驗是否服從威布爾分布和對數正態分布。

（4） 檢驗相同環境下樣本數據的離散系數是否小于4%（單擊離散系數）。由于較低的離散系數容易忽略材料和實驗方法等方面的差異性，因此對于小于4%的情況，程序會根據式（3）將離散系數修正為4%。修正后的數據平均值保持不變，原來小于平均值的數據變得更小，原來大于平均值的數據變得更大（見圖5）。然后，再檢查不同環境下離散系數是否相等，根據公式（4）計算F值小于臨界值（F=2.38<2.699），由此可以判斷不同環境下的離散系數相等。

圖5 離散系數修正圖

（5） 合并不同環境下的實驗數據再一次進行正態分布檢驗（勾選合并，見圖3）。在第（3）步中已經驗證每種環境下的實驗數據的是否符合正態分布，如圖6（a）所示，由圖6（b）可見合并后的實驗數據依然符合正態分布。然后，用每種環境下的樣本均值對數據進行歸一化，求出歸一化后合并數據的平均值、標準差和離散系數見表1。

根據公式（5）可以求出相應的環境下的折減系數Bj，并與不同環境下的樣本數據平均值相乘，從而獲得B基準值。

圖6 數據的正態分布檢驗

為了驗證程序的可靠性，采用標準中提供的范例中的計算結果進行對比說明。表2中列出了在四種不同實驗環境下的程序計算結果和航標HB7618?2009提供的結果。從表中的數據對比可以發現誤差在+0.3%以內，誤差主要由容限系數和折減系數的近似計算造成?？梢姵绦蛴嬎愕慕Y果符合工程精度（誤差±5%）的要求，此外相比于借助Excel的手動計算法方法，效率顯著提高。

4 結 語

本文根據航空工業標準HB7618?2009給出的B基準值的計算步驟，將Matlab的GUI模塊與B基準值的計算相結合開發出B基準值的計算程序。相比于借助Excle的手動計算方法，該程序界面具有以下優點：

（1） 程序可以判斷樣本是否來自同一母體，還具有異常數據處理、統計分布檢驗、離散系數檢驗與修正等功能，最后實驗結果以Word格式導出。

（2） 使用該程序進行B基準值的計算方便快捷、效率高，并且具有良好的人機界面，顯著提高了B基準值的計算效率。

表2 不同環境下計算結果對比

參考文獻

[1] 張駿華，盛祖銘，孫繼同.復合材料結構設計指南[M].北京：北京宇航出版社，1999.

[2] 李威，郭權鋒.碳纖維復合材料在航天領域的應用[J].中國光學，2011，4（3）：201?212.

[3] 陳紹杰.復合材料設計手冊[M].北京：航空工業出版社，1990.

[4] 劉建華，曹東，張曉云，等.樹脂基復合材料T300/5405的吸濕性能及濕熱環境對力學性能的影響[J].航空材料學報，2010，30（4）：75?80.

[5] 呂小軍，張琦，馬兆慶，等.濕熱老化對碳纖維/環氧樹脂基復合材料力學性能影響研究[J].材料工程，2005（11）：50?53.

[6] 鄭錫濤，李野，劉海燕，等.濕熱譜老化對復合材料層壓板強度的影響[J].航空學報，1998，19（4）：462?465.

[7] 中國航空工業總公司.HB 7618?2009 聚合物基復合材料力學性能數據表達準則[S].北京：中國航空工業總公司，2009：50?85.

[8] 陳子為.基于 Matlab GUI 掃雷游戲的設計與實現[J].現代電子技術，2009，32（24）：85?88.

[9] FAA. DOT/FAA/AR?03/19： material qualification and equivalency for polymer matrix composite material systems： updated procedure [S]. USA： FAA， 2003：55?78.

[10] 孫堅石，葉強.復合材料力學性能數據B基準值計算程序[J]. 航空制造技術，2009（z1）：19?21.

…

圖2 CTD環境下樣本數據進行正態分布檢驗圖

圖2中CTD環境下的樣本數據點分布大致成一條直線，并且檢驗結果H=0，表示接受正態分布假設?？梢奀TD樣本數據服從正態分布。

2.5 生成實驗報告

面對大量的實驗數據，書寫報告往往是重復性勞動工作，此時可以利用Matlab生成一表格模板，每次自動導入數據生成報告，這就可以節約大量的時間。Matlab作為一種面向對象的語言，它支持COM技術，利用它的COM編譯器可以把Matlab開發的程序轉換成方便使用的COM組件。而ActiveX控件技術是建立在COM技術之上，由微軟公司推出的共享程序數據和功能的技術。因此，Matlab可以利用ActiveX控件技術，在Matlab中調用Microsoft Word 程序插入表格，此時的Word是組件，是服務程序，而Matlab是控制器程序。下面是生成報告的部分程序代碼：

…

try

Word = actxGetRunningServer（′Word.Application′）；

%若Word服務器已經打開，返回其句柄Word

catch

Word = actxserver（′Word.Application′）；

%創建一個Microsoft Word服務器，返回句柄Word

…

Tables=Document.Tables.Add（Selection.Range，21，5）；

%插入一個21行5列的表格

…

Document.Save %保存文檔3 計算實例

本文根據航標HB7618?2009提供的范例數據來詳細說明B基準值的計算流程，并驗證程序界面計算的可靠性，如圖3所示。范例包含了所有實驗環境下的樣本數據，見表1。

圖3 B基準值程序界面

（1） 通過Anderson?Darling方法檢驗不同環境下的不同批次的實驗數據是否來自同一母體（依次勾選對應的實驗環境），檢驗后發現ETW（高溫濕態）環境下的三批樣本數據的ADK大于臨界值（ADK=2.258 2>ADC=1.917 4），可見該環境下三批數據來自不同母體，如圖4（a）所示。如果放大為1%錯判風險下， ADK小于臨界值（ADK=2.258 2

（2） 相容性檢驗檢驗結束后，根據公式（2）計算每一試樣的MNR值（單擊異常數據檢查），如果MNR值大于臨界值，那么就可以判斷該試樣為異常值。對于能夠反映材料、工藝參數和實驗環境的差異性的異常數據需要保留。而當異常數據是由不合格的試樣、實驗設備和夾具的缺陷造成時，異常數據要刪除（單擊刪除異常數據，見圖3），否則就會導致統計結果的失真。

圖4 母體檢驗

（3） 對不同環境下實驗數據進行正態分布假設檢驗，如果不符合正態分布再依次檢驗是否服從威布爾分布和對數正態分布。

（4） 檢驗相同環境下樣本數據的離散系數是否小于4%（單擊離散系數）。由于較低的離散系數容易忽略材料和實驗方法等方面的差異性，因此對于小于4%的情況，程序會根據式（3）將離散系數修正為4%。修正后的數據平均值保持不變，原來小于平均值的數據變得更小，原來大于平均值的數據變得更大（見圖5）。然后，再檢查不同環境下離散系數是否相等，根據公式（4）計算F值小于臨界值（F=2.38<2.699），由此可以判斷不同環境下的離散系數相等。

圖5 離散系數修正圖

（5） 合并不同環境下的實驗數據再一次進行正態分布檢驗（勾選合并，見圖3）。在第（3）步中已經驗證每種環境下的實驗數據的是否符合正態分布，如圖6（a）所示，由圖6（b）可見合并后的實驗數據依然符合正態分布。然后，用每種環境下的樣本均值對數據進行歸一化，求出歸一化后合并數據的平均值、標準差和離散系數見表1。

根據公式（5）可以求出相應的環境下的折減系數Bj，并與不同環境下的樣本數據平均值相乘，從而獲得B基準值。

圖6 數據的正態分布檢驗

為了驗證程序的可靠性，采用標準中提供的范例中的計算結果進行對比說明。表2中列出了在四種不同實驗環境下的程序計算結果和航標HB7618?2009提供的結果。從表中的數據對比可以發現誤差在+0.3%以內，誤差主要由容限系數和折減系數的近似計算造成。可見程序計算的結果符合工程精度（誤差±5%）的要求，此外相比于借助Excel的手動計算法方法，效率顯著提高。

4 結 語

本文根據航空工業標準HB7618?2009給出的B基準值的計算步驟，將Matlab的GUI模塊與B基準值的計算相結合開發出B基準值的計算程序。相比于借助Excle的手動計算方法，該程序界面具有以下優點：

（1） 程序可以判斷樣本是否來自同一母體，還具有異常數據處理、統計分布檢驗、離散系數檢驗與修正等功能，最后實驗結果以Word格式導出。

（2） 使用該程序進行B基準值的計算方便快捷、效率高，并且具有良好的人機界面，顯著提高了B基準值的計算效率。

表2 不同環境下計算結果對比

參考文獻

[1] 張駿華，盛祖銘，孫繼同.復合材料結構設計指南[M].北京：北京宇航出版社，1999.

[2] 李威，郭權鋒.碳纖維復合材料在航天領域的應用[J].中國光學，2011，4（3）：201?212.

[3] 陳紹杰.復合材料設計手冊[M].北京：航空工業出版社，1990.

[4] 劉建華，曹東，張曉云，等.樹脂基復合材料T300/5405的吸濕性能及濕熱環境對力學性能的影響[J].航空材料學報，2010，30（4）：75?80.

[5] 呂小軍，張琦，馬兆慶，等.濕熱老化對碳纖維/環氧樹脂基復合材料力學性能影響研究[J].材料工程，2005（11）：50?53.

[6] 鄭錫濤，李野，劉海燕，等.濕熱譜老化對復合材料層壓板強度的影響[J].航空學報，1998，19（4）：462?465.

[7] 中國航空工業總公司.HB 7618?2009 聚合物基復合材料力學性能數據表達準則[S].北京：中國航空工業總公司，2009：50?85.

[8] 陳子為.基于 Matlab GUI 掃雷游戲的設計與實現[J].現代電子技術，2009，32（24）：85?88.

[9] FAA. DOT/FAA/AR?03/19： material qualification and equivalency for polymer matrix composite material systems： updated procedure [S]. USA： FAA， 2003：55?78.

[10] 孫堅石，葉強.復合材料力學性能數據B基準值計算程序[J]. 航空制造技術，2009（z1）：19?21.