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

    基于Matlab的碳纖維復(fù)合材料的B基準(zhǔn)值算法

    2014-08-20 18:09:58董本正李偉高麗紅郭朝龍査萌吉愛(ài)
    現(xiàn)代電子技術(shù) 2014年16期

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

    摘 要: B基準(zhǔn)值是評(píng)價(jià)碳纖維復(fù)合材料性能的重要參數(shù),對(duì)B基準(zhǔn)值的高效精確的計(jì)算有助于評(píng)價(jià)復(fù)合材料的性能。根據(jù)航空工業(yè)標(biāo)準(zhǔn)HB7618?2009,開發(fā)出基于Matlab的 GUI模塊的B基準(zhǔn)值的算法,該程序可以實(shí)現(xiàn)樣本母體檢驗(yàn)、異常數(shù)據(jù)檢查、統(tǒng)計(jì)分布檢驗(yàn)、離散系數(shù)修正與檢驗(yàn)、B基準(zhǔn)值的計(jì)算等功能,最后以Word格式導(dǎo)出計(jì)算結(jié)果。與標(biāo)準(zhǔn)提供范例提供的計(jì)算結(jié)果進(jìn)行對(duì)比,證明程序運(yùn)行準(zhǔn)確可靠,并且其界面操作簡(jiǎn)單,顯著地提高了B基準(zhǔn)值計(jì)算效率。

    關(guān)鍵詞: 碳纖維增強(qiáng)復(fù)合材料; B基準(zhǔn)值; 統(tǒng)計(jì)學(xué)分析; Matlab

    中圖分類號(hào): TP319 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 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 引 言

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

    1 基準(zhǔn)值計(jì)算原理

    1.1 B基準(zhǔn)值

    為了保證材料的結(jié)構(gòu)的可靠性和安全性,就需要用到一些高級(jí)的統(tǒng)計(jì)學(xué)方法來(lái)確定其合理的設(shè)計(jì)許用值,一般采用B基準(zhǔn)值[7]。B基準(zhǔn)值是建立在統(tǒng)計(jì)學(xué)上的衡量材料性能的參數(shù),在95%的置信度下,90%性能數(shù)值群的值高于此值[9]。航空工業(yè)標(biāo)準(zhǔn)HB7618?2009關(guān)于B基準(zhǔn)值的流程復(fù)雜,本文將著重從批次相容性檢驗(yàn)、檢查異常數(shù)據(jù)、分布函數(shù)檢驗(yàn)、離散系數(shù)修正與檢驗(yàn)、B基準(zhǔn)值計(jì)算及導(dǎo)出結(jié)果幾方面進(jìn)行說(shuō)明。

    1.2 主要計(jì)算步驟

    (1) 檢驗(yàn)不同批次的實(shí)驗(yàn)數(shù)據(jù)是否來(lái)自同一母體;檢驗(yàn)方法是K樣本的Anderson ?Darling檢驗(yàn), 其統(tǒng)計(jì)公式如下[7]:

    [ADK=n-1n2(k-1)i=1k1nij=1lhjnFij-njHjHj(n-Hj)-nhj4] (1)

    如果ADK小于錯(cuò)判臨界值A(chǔ)DC(5%錯(cuò)判風(fēng)險(xiǎn)),則不同批次的實(shí)驗(yàn)數(shù)據(jù)來(lái)自不同母體。如果相容性檢驗(yàn)不通過(guò)的,為了增加樣本數(shù)據(jù)來(lái)自同一母體的可能性,采用1%錯(cuò)判風(fēng)險(xiǎn)繼續(xù)判斷是否來(lái)自同一母體。

    (1) 檢查異常數(shù)據(jù)。采用最大賦范殘差的方法定量的檢查異常數(shù)據(jù):

    [MNR=maxxi-xs] (2)

    (2) 分布函數(shù)檢驗(yàn)。根據(jù)經(jīng)驗(yàn)威布爾計(jì)算的B基值過(guò)于保守,因此優(yōu)先檢驗(yàn)樣本數(shù)據(jù)是否符合正態(tài)分布,如果不符合正態(tài)分布,再依次檢驗(yàn)樣本數(shù)據(jù)是否服從威布爾分布和對(duì)數(shù)正態(tài)分布。

    (3) 離散系數(shù)修正與檢驗(yàn):當(dāng)離散系數(shù)低于4%時(shí),容易忽略材料差異、實(shí)驗(yàn)方法等方面的實(shí)際差異性,因此要對(duì)離散系數(shù)進(jìn)行修正[10],修正方法為:

    [x*i=xi-α(xi)Δ] (3)

    離散系數(shù)修正結(jié)束后,對(duì)不同環(huán)境下的樣本數(shù)據(jù)進(jìn)行離散系數(shù)檢驗(yàn),離散系數(shù)的檢驗(yàn)公式為:

    [F=i=1k(wi-w)(k-1)i=1kj=1niwij-wi(n-k)] (4)

    (4) 將不同環(huán)境下的數(shù)據(jù)歸一化,根據(jù)歸一化后的實(shí)驗(yàn)數(shù)據(jù)的平均值和標(biāo)準(zhǔn)差計(jì)算B基準(zhǔn)值。求出對(duì)應(yīng)的折減系數(shù),再與不同環(huán)境下樣本的平均值相乘,從而得到B基準(zhǔn)值,折減系數(shù)的公式為:

    [Bj=x-(kj)*s] (5)

    2 程序?qū)崿F(xiàn)

    2.1 基準(zhǔn)值算法的漸進(jìn)流程圖和程序界面

    圖1是根據(jù)航標(biāo)HB7618?2009給出的計(jì)算步驟編制的程序流程圖,該程序流程圖包括了批次相容性檢驗(yàn)、檢查異常數(shù)據(jù)、分布函數(shù)檢驗(yàn)、離散系數(shù)修正與檢驗(yàn)、B基準(zhǔn)值計(jì)算及導(dǎo)出結(jié)果等功能。

    圖1 B基準(zhǔn)值算法的漸進(jìn)流程圖

    2.2 變量傳遞

    B基準(zhǔn)值的算法程序包含幾個(gè)不同的模塊(見(jiàn)圖1),而每一個(gè)模塊又包含了多個(gè)程序和子程序。為了實(shí)現(xiàn)變量在不同模塊和函數(shù)間進(jìn)行傳遞,本文采用handles結(jié)構(gòu)體和guidata函數(shù)來(lái)實(shí)現(xiàn)這個(gè)功能。在B基準(zhǔn)值的算法程序界面中的所有控件使用同一個(gè)handles結(jié)構(gòu)體,handles結(jié)構(gòu)體中不僅保存了圖形窗口中所有對(duì)象的句柄,而且還可以獲取或設(shè)置某個(gè)對(duì)象的屬性。在程序的計(jì)算過(guò)程中需要將每一個(gè)模塊計(jì)算后的結(jié)果賦給handles結(jié)構(gòu)體,以便于將這個(gè)變量傳遞到下一個(gè)計(jì)算模塊。而每次handles結(jié)構(gòu)體添加了新的元素,guidata函數(shù)便會(huì)更新handles結(jié)構(gòu)體,使handle結(jié)構(gòu)體中的新添加元素可以傳遞到下一個(gè)模塊。guidata函數(shù)用法guidata(hObject,handles),其中,hObject是執(zhí)行回調(diào)的控件對(duì)象的句柄,handle為結(jié)構(gòu)體。handles結(jié)構(gòu)體和guidata函數(shù)的用法如下所示:

    function pushbutton13_Bvalue_Callback(hObject, eventdata, handles)

    handles.Bvalue(j)=roundn(mean(nn_data(:))*Bj(j),?3);

    guidata(hObject,handles)

    2.3 設(shè)置實(shí)驗(yàn)環(huán)境

    B基準(zhǔn)值的算法界面程序首先要檢驗(yàn)相同環(huán)境下不同批次的樣本數(shù)據(jù)是否來(lái)自同一母體。程序界面中包括四種實(shí)驗(yàn)環(huán)境:低溫干態(tài)(CTD)、常溫干態(tài)(RTD)、高溫干態(tài)(ETD)、高溫濕態(tài)(ETW), 為了避免母體檢驗(yàn)的過(guò)程中不同環(huán)境下的樣本數(shù)據(jù)相互干擾,需要對(duì)程序界面中代表著不同實(shí)驗(yàn)環(huán)境的復(fù)選框(Check Box)進(jìn)行設(shè)置。當(dāng)其中的一種環(huán)境的復(fù)選框處于選中狀態(tài)時(shí),其他的三種環(huán)境的復(fù)選框需要設(shè)置成未選中狀態(tài),這樣就可以避免相互干擾。下面是當(dāng)CTD處于選中狀態(tài),其他的三種實(shí)驗(yàn)環(huán)境處于未選中狀態(tài)時(shí)的代碼:

    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

    當(dāng)復(fù)選框處于選中狀態(tài)時(shí),其Value值等于1;當(dāng)復(fù)選框處于未選中狀態(tài)時(shí),其Value值等于0。

    2.4 分布函數(shù)檢驗(yàn)

    對(duì)于樣本數(shù)據(jù)的分布函數(shù)檢驗(yàn),本文采用Matlab本身的內(nèi)置函數(shù)normplot()和weibplot()函數(shù)進(jìn)行判斷。若樣本數(shù)據(jù)服從正態(tài)分布和威布爾分布,樣本數(shù)據(jù)點(diǎn)分布成一條直線,否則成一條曲線。至于樣本數(shù)據(jù)的對(duì)數(shù)正態(tài)分布檢驗(yàn),程序首先對(duì)樣本數(shù)據(jù)取對(duì)數(shù)將對(duì)數(shù)分布轉(zhuǎn)化為正態(tài)分布后,再用normplot()函數(shù)檢驗(yàn)是否服從正態(tài)分布。

    在B基準(zhǔn)值的計(jì)算過(guò)程中首先檢驗(yàn)的是正態(tài)分布,可見(jiàn)樣本數(shù)據(jù)的正態(tài)分布對(duì)于B基準(zhǔn)值計(jì)算的重要性。為了增加正態(tài)分布檢驗(yàn)的準(zhǔn)確性,本文除了采用normplot()函數(shù)對(duì)樣本數(shù)據(jù)進(jìn)行直觀的目視正態(tài)分布檢驗(yàn)外,還借助函數(shù)H=jbtest()對(duì)樣本數(shù)據(jù)進(jìn)行進(jìn)一步的正態(tài)分布檢驗(yàn)。若H=0,則認(rèn)X服從正態(tài)分布;若H=1,則否定X服從正態(tài)分布。對(duì)CTD環(huán)境下樣本數(shù)據(jù)進(jìn)行正態(tài)分布檢驗(yàn)的程序代碼如下,檢測(cè)結(jié)果如圖2所示。

    normplot(handles.tdata(:)); %正態(tài)分布概率圖

    handles.H=jbtest(handles.tdata(:)); %擬合優(yōu)度檢驗(yàn)

    圖2 CTD環(huán)境下樣本數(shù)據(jù)進(jìn)行正態(tài)分布檢驗(yàn)圖

    圖2中CTD環(huán)境下的樣本數(shù)據(jù)點(diǎn)分布大致成一條直線,并且檢驗(yàn)結(jié)果H=0,表示接受正態(tài)分布假設(shè)??梢?jiàn)CTD樣本數(shù)據(jù)服從正態(tài)分布。

    2.5 生成實(shí)驗(yàn)報(bào)告

    面對(duì)大量的實(shí)驗(yàn)數(shù)據(jù),書寫報(bào)告往往是重復(fù)性勞動(dòng)工作,此時(shí)可以利用Matlab生成一表格模板,每次自動(dòng)導(dǎo)入數(shù)據(jù)生成報(bào)告,這就可以節(jié)約大量的時(shí)間。Matlab作為一種面向?qū)ο蟮恼Z(yǔ)言,它支持COM技術(shù),利用它的COM編譯器可以把Matlab開發(fā)的程序轉(zhuǎn)換成方便使用的COM組件。而ActiveX控件技術(shù)是建立在COM技術(shù)之上,由微軟公司推出的共享程序數(shù)據(jù)和功能的技術(shù)。因此,Matlab可以利用ActiveX控件技術(shù),在Matlab中調(diào)用Microsoft Word 程序插入表格,此時(shí)的Word是組件,是服務(wù)程序,而Matlab是控制器程序。下面是生成報(bào)告的部分程序代碼:

    try

    Word = actxGetRunningServer(′Word.Application′);

    %若Word服務(wù)器已經(jīng)打開,返回其句柄Word

    catch

    Word = actxserver(′Word.Application′);

    %創(chuàng)建一個(gè)Microsoft Word服務(wù)器,返回句柄Word

    Tables=Document.Tables.Add(Selection.Range,21,5);

    %插入一個(gè)21行5列的表格

    Document.Save %保存文檔3 計(jì)算實(shí)例

    本文根據(jù)航標(biāo)HB7618?2009提供的范例數(shù)據(jù)來(lái)詳細(xì)說(shuō)明B基準(zhǔn)值的計(jì)算流程,并驗(yàn)證程序界面計(jì)算的可靠性,如圖3所示。范例包含了所有實(shí)驗(yàn)環(huán)境下的樣本數(shù)據(jù),見(jiàn)表1。

    圖3 B基準(zhǔn)值程序界面

    (1) 通過(guò)Anderson?Darling方法檢驗(yàn)不同環(huán)境下的不同批次的實(shí)驗(yàn)數(shù)據(jù)是否來(lái)自同一母體(依次勾選對(duì)應(yīng)的實(shí)驗(yàn)環(huán)境),檢驗(yàn)后發(fā)現(xiàn)ETW(高溫濕態(tài))環(huán)境下的三批樣本數(shù)據(jù)的ADK大于臨界值(ADK=2.258 2>ADC=1.917 4),可見(jiàn)該環(huán)境下三批數(shù)據(jù)來(lái)自不同母體,如圖4(a)所示。如果放大為1%錯(cuò)判風(fēng)險(xiǎn)下, ADK小于臨界值(ADK=2.258 2

    (2) 相容性檢驗(yàn)檢驗(yàn)結(jié)束后,根據(jù)公式(2)計(jì)算每一試樣的MNR值(單擊異常數(shù)據(jù)檢查),如果MNR值大于臨界值,那么就可以判斷該試樣為異常值。對(duì)于能夠反映材料、工藝參數(shù)和實(shí)驗(yàn)環(huán)境的差異性的異常數(shù)據(jù)需要保留。而當(dāng)異常數(shù)據(jù)是由不合格的試樣、實(shí)驗(yàn)設(shè)備和夾具的缺陷造成時(shí),異常數(shù)據(jù)要?jiǎng)h除(單擊刪除異常數(shù)據(jù),見(jiàn)圖3),否則就會(huì)導(dǎo)致統(tǒng)計(jì)結(jié)果的失真。

    圖4 母體檢驗(yàn)

    (3) 對(duì)不同環(huán)境下實(shí)驗(yàn)數(shù)據(jù)進(jìn)行正態(tài)分布假設(shè)檢驗(yàn),如果不符合正態(tài)分布再依次檢驗(yàn)是否服從威布爾分布和對(duì)數(shù)正態(tài)分布。

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

    圖5 離散系數(shù)修正圖

    (5) 合并不同環(huán)境下的實(shí)驗(yàn)數(shù)據(jù)再一次進(jìn)行正態(tài)分布檢驗(yàn)(勾選合并,見(jiàn)圖3)。在第(3)步中已經(jīng)驗(yàn)證每種環(huán)境下的實(shí)驗(yàn)數(shù)據(jù)的是否符合正態(tài)分布,如圖6(a)所示,由圖6(b)可見(jiàn)合并后的實(shí)驗(yàn)數(shù)據(jù)依然符合正態(tài)分布。然后,用每種環(huán)境下的樣本均值對(duì)數(shù)據(jù)進(jìn)行歸一化,求出歸一化后合并數(shù)據(jù)的平均值、標(biāo)準(zhǔn)差和離散系數(shù)見(jiàn)表1。

    根據(jù)公式(5)可以求出相應(yīng)的環(huán)境下的折減系數(shù)Bj,并與不同環(huán)境下的樣本數(shù)據(jù)平均值相乘,從而獲得B基準(zhǔn)值。

    圖6 數(shù)據(jù)的正態(tài)分布檢驗(yàn)

    為了驗(yàn)證程序的可靠性,采用標(biāo)準(zhǔn)中提供的范例中的計(jì)算結(jié)果進(jìn)行對(duì)比說(shuō)明。表2中列出了在四種不同實(shí)驗(yàn)環(huán)境下的程序計(jì)算結(jié)果和航標(biāo)HB7618?2009提供的結(jié)果。從表中的數(shù)據(jù)對(duì)比可以發(fā)現(xiàn)誤差在+0.3%以內(nèi),誤差主要由容限系數(shù)和折減系數(shù)的近似計(jì)算造成??梢?jiàn)程序計(jì)算的結(jié)果符合工程精度(誤差±5%)的要求,此外相比于借助Excel的手動(dòng)計(jì)算法方法,效率顯著提高。

    4 結(jié) 語(yǔ)

    本文根據(jù)航空工業(yè)標(biāo)準(zhǔn)HB7618?2009給出的B基準(zhǔn)值的計(jì)算步驟,將Matlab的GUI模塊與B基準(zhǔn)值的計(jì)算相結(jié)合開發(fā)出B基準(zhǔn)值的計(jì)算程序。相比于借助Excle的手動(dòng)計(jì)算方法,該程序界面具有以下優(yōu)點(diǎn):

    (1) 程序可以判斷樣本是否來(lái)自同一母體,還具有異常數(shù)據(jù)處理、統(tǒng)計(jì)分布檢驗(yàn)、離散系數(shù)檢驗(yàn)與修正等功能,最后實(shí)驗(yàn)結(jié)果以Word格式導(dǎo)出。

    (2) 使用該程序進(jìn)行B基準(zhǔn)值的計(jì)算方便快捷、效率高,并且具有良好的人機(jī)界面,顯著提高了B基準(zhǔn)值的計(jì)算效率。

    表2 不同環(huán)境下計(jì)算結(jié)果對(duì)比

    參考文獻(xiàn)

    [1] 張駿華,盛祖銘,孫繼同.復(fù)合材料結(jié)構(gòu)設(shè)計(jì)指南[M].北京:北京宇航出版社,1999.

    [2] 李威,郭權(quán)鋒.碳纖維復(fù)合材料在航天領(lǐng)域的應(yīng)用[J].中國(guó)光學(xué),2011,4(3):201?212.

    [3] 陳紹杰.復(fù)合材料設(shè)計(jì)手冊(cè)[M].北京:航空工業(yè)出版社,1990.

    [4] 劉建華,曹東,張曉云,等.樹脂基復(fù)合材料T300/5405的吸濕性能及濕熱環(huán)境對(duì)力學(xué)性能的影響[J].航空材料學(xué)報(bào),2010,30(4):75?80.

    [5] 呂小軍,張琦,馬兆慶,等.濕熱老化對(duì)碳纖維/環(huán)氧樹脂基復(fù)合材料力學(xué)性能影響研究[J].材料工程,2005(11):50?53.

    [6] 鄭錫濤,李野,劉海燕,等.濕熱譜老化對(duì)復(fù)合材料層壓板強(qiáng)度的影響[J].航空學(xué)報(bào),1998,19(4):462?465.

    [7] 中國(guó)航空工業(yè)總公司.HB 7618?2009 聚合物基復(fù)合材料力學(xué)性能數(shù)據(jù)表達(dá)準(zhǔn)則[S].北京:中國(guó)航空工業(yè)總公司,2009:50?85.

    [8] 陳子為.基于 Matlab GUI 掃雷游戲的設(shè)計(jì)與實(shí)現(xiàn)[J].現(xiàn)代電子技術(shù),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] 孫堅(jiān)石,葉強(qiáng).復(fù)合材料力學(xué)性能數(shù)據(jù)B基準(zhǔn)值計(jì)算程序[J]. 航空制造技術(shù),2009(z1):19?21.

    圖2 CTD環(huán)境下樣本數(shù)據(jù)進(jìn)行正態(tài)分布檢驗(yàn)圖

    圖2中CTD環(huán)境下的樣本數(shù)據(jù)點(diǎn)分布大致成一條直線,并且檢驗(yàn)結(jié)果H=0,表示接受正態(tài)分布假設(shè)??梢?jiàn)CTD樣本數(shù)據(jù)服從正態(tài)分布。

    2.5 生成實(shí)驗(yàn)報(bào)告

    面對(duì)大量的實(shí)驗(yàn)數(shù)據(jù),書寫報(bào)告往往是重復(fù)性勞動(dòng)工作,此時(shí)可以利用Matlab生成一表格模板,每次自動(dòng)導(dǎo)入數(shù)據(jù)生成報(bào)告,這就可以節(jié)約大量的時(shí)間。Matlab作為一種面向?qū)ο蟮恼Z(yǔ)言,它支持COM技術(shù),利用它的COM編譯器可以把Matlab開發(fā)的程序轉(zhuǎn)換成方便使用的COM組件。而ActiveX控件技術(shù)是建立在COM技術(shù)之上,由微軟公司推出的共享程序數(shù)據(jù)和功能的技術(shù)。因此,Matlab可以利用ActiveX控件技術(shù),在Matlab中調(diào)用Microsoft Word 程序插入表格,此時(shí)的Word是組件,是服務(wù)程序,而Matlab是控制器程序。下面是生成報(bào)告的部分程序代碼:

    try

    Word = actxGetRunningServer(′Word.Application′);

    %若Word服務(wù)器已經(jīng)打開,返回其句柄Word

    catch

    Word = actxserver(′Word.Application′);

    %創(chuàng)建一個(gè)Microsoft Word服務(wù)器,返回句柄Word

    Tables=Document.Tables.Add(Selection.Range,21,5);

    %插入一個(gè)21行5列的表格

    Document.Save %保存文檔3 計(jì)算實(shí)例

    本文根據(jù)航標(biāo)HB7618?2009提供的范例數(shù)據(jù)來(lái)詳細(xì)說(shuō)明B基準(zhǔn)值的計(jì)算流程,并驗(yàn)證程序界面計(jì)算的可靠性,如圖3所示。范例包含了所有實(shí)驗(yàn)環(huán)境下的樣本數(shù)據(jù),見(jiàn)表1。

    圖3 B基準(zhǔn)值程序界面

    (1) 通過(guò)Anderson?Darling方法檢驗(yàn)不同環(huán)境下的不同批次的實(shí)驗(yàn)數(shù)據(jù)是否來(lái)自同一母體(依次勾選對(duì)應(yīng)的實(shí)驗(yàn)環(huán)境),檢驗(yàn)后發(fā)現(xiàn)ETW(高溫濕態(tài))環(huán)境下的三批樣本數(shù)據(jù)的ADK大于臨界值(ADK=2.258 2>ADC=1.917 4),可見(jiàn)該環(huán)境下三批數(shù)據(jù)來(lái)自不同母體,如圖4(a)所示。如果放大為1%錯(cuò)判風(fēng)險(xiǎn)下, ADK小于臨界值(ADK=2.258 2

    (2) 相容性檢驗(yàn)檢驗(yàn)結(jié)束后,根據(jù)公式(2)計(jì)算每一試樣的MNR值(單擊異常數(shù)據(jù)檢查),如果MNR值大于臨界值,那么就可以判斷該試樣為異常值。對(duì)于能夠反映材料、工藝參數(shù)和實(shí)驗(yàn)環(huán)境的差異性的異常數(shù)據(jù)需要保留。而當(dāng)異常數(shù)據(jù)是由不合格的試樣、實(shí)驗(yàn)設(shè)備和夾具的缺陷造成時(shí),異常數(shù)據(jù)要?jiǎng)h除(單擊刪除異常數(shù)據(jù),見(jiàn)圖3),否則就會(huì)導(dǎo)致統(tǒng)計(jì)結(jié)果的失真。

    圖4 母體檢驗(yàn)

    (3) 對(duì)不同環(huán)境下實(shí)驗(yàn)數(shù)據(jù)進(jìn)行正態(tài)分布假設(shè)檢驗(yàn),如果不符合正態(tài)分布再依次檢驗(yàn)是否服從威布爾分布和對(duì)數(shù)正態(tài)分布。

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

    圖5 離散系數(shù)修正圖

    (5) 合并不同環(huán)境下的實(shí)驗(yàn)數(shù)據(jù)再一次進(jìn)行正態(tài)分布檢驗(yàn)(勾選合并,見(jiàn)圖3)。在第(3)步中已經(jīng)驗(yàn)證每種環(huán)境下的實(shí)驗(yàn)數(shù)據(jù)的是否符合正態(tài)分布,如圖6(a)所示,由圖6(b)可見(jiàn)合并后的實(shí)驗(yàn)數(shù)據(jù)依然符合正態(tài)分布。然后,用每種環(huán)境下的樣本均值對(duì)數(shù)據(jù)進(jìn)行歸一化,求出歸一化后合并數(shù)據(jù)的平均值、標(biāo)準(zhǔn)差和離散系數(shù)見(jiàn)表1。

    根據(jù)公式(5)可以求出相應(yīng)的環(huán)境下的折減系數(shù)Bj,并與不同環(huán)境下的樣本數(shù)據(jù)平均值相乘,從而獲得B基準(zhǔn)值。

    圖6 數(shù)據(jù)的正態(tài)分布檢驗(yàn)

    為了驗(yàn)證程序的可靠性,采用標(biāo)準(zhǔn)中提供的范例中的計(jì)算結(jié)果進(jìn)行對(duì)比說(shuō)明。表2中列出了在四種不同實(shí)驗(yàn)環(huán)境下的程序計(jì)算結(jié)果和航標(biāo)HB7618?2009提供的結(jié)果。從表中的數(shù)據(jù)對(duì)比可以發(fā)現(xiàn)誤差在+0.3%以內(nèi),誤差主要由容限系數(shù)和折減系數(shù)的近似計(jì)算造成。可見(jiàn)程序計(jì)算的結(jié)果符合工程精度(誤差±5%)的要求,此外相比于借助Excel的手動(dòng)計(jì)算法方法,效率顯著提高。

    4 結(jié) 語(yǔ)

    本文根據(jù)航空工業(yè)標(biāo)準(zhǔn)HB7618?2009給出的B基準(zhǔn)值的計(jì)算步驟,將Matlab的GUI模塊與B基準(zhǔn)值的計(jì)算相結(jié)合開發(fā)出B基準(zhǔn)值的計(jì)算程序。相比于借助Excle的手動(dòng)計(jì)算方法,該程序界面具有以下優(yōu)點(diǎn):

    (1) 程序可以判斷樣本是否來(lái)自同一母體,還具有異常數(shù)據(jù)處理、統(tǒng)計(jì)分布檢驗(yàn)、離散系數(shù)檢驗(yàn)與修正等功能,最后實(shí)驗(yàn)結(jié)果以Word格式導(dǎo)出。

    (2) 使用該程序進(jìn)行B基準(zhǔn)值的計(jì)算方便快捷、效率高,并且具有良好的人機(jī)界面,顯著提高了B基準(zhǔn)值的計(jì)算效率。

    表2 不同環(huán)境下計(jì)算結(jié)果對(duì)比

    參考文獻(xiàn)

    [1] 張駿華,盛祖銘,孫繼同.復(fù)合材料結(jié)構(gòu)設(shè)計(jì)指南[M].北京:北京宇航出版社,1999.

    [2] 李威,郭權(quán)鋒.碳纖維復(fù)合材料在航天領(lǐng)域的應(yīng)用[J].中國(guó)光學(xué),2011,4(3):201?212.

    [3] 陳紹杰.復(fù)合材料設(shè)計(jì)手冊(cè)[M].北京:航空工業(yè)出版社,1990.

    [4] 劉建華,曹東,張曉云,等.樹脂基復(fù)合材料T300/5405的吸濕性能及濕熱環(huán)境對(duì)力學(xué)性能的影響[J].航空材料學(xué)報(bào),2010,30(4):75?80.

    [5] 呂小軍,張琦,馬兆慶,等.濕熱老化對(duì)碳纖維/環(huán)氧樹脂基復(fù)合材料力學(xué)性能影響研究[J].材料工程,2005(11):50?53.

    [6] 鄭錫濤,李野,劉海燕,等.濕熱譜老化對(duì)復(fù)合材料層壓板強(qiáng)度的影響[J].航空學(xué)報(bào),1998,19(4):462?465.

    [7] 中國(guó)航空工業(yè)總公司.HB 7618?2009 聚合物基復(fù)合材料力學(xué)性能數(shù)據(jù)表達(dá)準(zhǔn)則[S].北京:中國(guó)航空工業(yè)總公司,2009:50?85.

    [8] 陳子為.基于 Matlab GUI 掃雷游戲的設(shè)計(jì)與實(shí)現(xiàn)[J].現(xiàn)代電子技術(shù),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] 孫堅(jiān)石,葉強(qiáng).復(fù)合材料力學(xué)性能數(shù)據(jù)B基準(zhǔn)值計(jì)算程序[J]. 航空制造技術(shù),2009(z1):19?21.

    圖2 CTD環(huán)境下樣本數(shù)據(jù)進(jìn)行正態(tài)分布檢驗(yàn)圖

    圖2中CTD環(huán)境下的樣本數(shù)據(jù)點(diǎn)分布大致成一條直線,并且檢驗(yàn)結(jié)果H=0,表示接受正態(tài)分布假設(shè)??梢?jiàn)CTD樣本數(shù)據(jù)服從正態(tài)分布。

    2.5 生成實(shí)驗(yàn)報(bào)告

    面對(duì)大量的實(shí)驗(yàn)數(shù)據(jù),書寫報(bào)告往往是重復(fù)性勞動(dòng)工作,此時(shí)可以利用Matlab生成一表格模板,每次自動(dòng)導(dǎo)入數(shù)據(jù)生成報(bào)告,這就可以節(jié)約大量的時(shí)間。Matlab作為一種面向?qū)ο蟮恼Z(yǔ)言,它支持COM技術(shù),利用它的COM編譯器可以把Matlab開發(fā)的程序轉(zhuǎn)換成方便使用的COM組件。而ActiveX控件技術(shù)是建立在COM技術(shù)之上,由微軟公司推出的共享程序數(shù)據(jù)和功能的技術(shù)。因此,Matlab可以利用ActiveX控件技術(shù),在Matlab中調(diào)用Microsoft Word 程序插入表格,此時(shí)的Word是組件,是服務(wù)程序,而Matlab是控制器程序。下面是生成報(bào)告的部分程序代碼:

    try

    Word = actxGetRunningServer(′Word.Application′);

    %若Word服務(wù)器已經(jīng)打開,返回其句柄Word

    catch

    Word = actxserver(′Word.Application′);

    %創(chuàng)建一個(gè)Microsoft Word服務(wù)器,返回句柄Word

    Tables=Document.Tables.Add(Selection.Range,21,5);

    %插入一個(gè)21行5列的表格

    Document.Save %保存文檔3 計(jì)算實(shí)例

    本文根據(jù)航標(biāo)HB7618?2009提供的范例數(shù)據(jù)來(lái)詳細(xì)說(shuō)明B基準(zhǔn)值的計(jì)算流程,并驗(yàn)證程序界面計(jì)算的可靠性,如圖3所示。范例包含了所有實(shí)驗(yàn)環(huán)境下的樣本數(shù)據(jù),見(jiàn)表1。

    圖3 B基準(zhǔn)值程序界面

    (1) 通過(guò)Anderson?Darling方法檢驗(yàn)不同環(huán)境下的不同批次的實(shí)驗(yàn)數(shù)據(jù)是否來(lái)自同一母體(依次勾選對(duì)應(yīng)的實(shí)驗(yàn)環(huán)境),檢驗(yàn)后發(fā)現(xiàn)ETW(高溫濕態(tài))環(huán)境下的三批樣本數(shù)據(jù)的ADK大于臨界值(ADK=2.258 2>ADC=1.917 4),可見(jiàn)該環(huán)境下三批數(shù)據(jù)來(lái)自不同母體,如圖4(a)所示。如果放大為1%錯(cuò)判風(fēng)險(xiǎn)下, ADK小于臨界值(ADK=2.258 2

    (2) 相容性檢驗(yàn)檢驗(yàn)結(jié)束后,根據(jù)公式(2)計(jì)算每一試樣的MNR值(單擊異常數(shù)據(jù)檢查),如果MNR值大于臨界值,那么就可以判斷該試樣為異常值。對(duì)于能夠反映材料、工藝參數(shù)和實(shí)驗(yàn)環(huán)境的差異性的異常數(shù)據(jù)需要保留。而當(dāng)異常數(shù)據(jù)是由不合格的試樣、實(shí)驗(yàn)設(shè)備和夾具的缺陷造成時(shí),異常數(shù)據(jù)要?jiǎng)h除(單擊刪除異常數(shù)據(jù),見(jiàn)圖3),否則就會(huì)導(dǎo)致統(tǒng)計(jì)結(jié)果的失真。

    圖4 母體檢驗(yàn)

    (3) 對(duì)不同環(huán)境下實(shí)驗(yàn)數(shù)據(jù)進(jìn)行正態(tài)分布假設(shè)檢驗(yàn),如果不符合正態(tài)分布再依次檢驗(yàn)是否服從威布爾分布和對(duì)數(shù)正態(tài)分布。

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

    圖5 離散系數(shù)修正圖

    (5) 合并不同環(huán)境下的實(shí)驗(yàn)數(shù)據(jù)再一次進(jìn)行正態(tài)分布檢驗(yàn)(勾選合并,見(jiàn)圖3)。在第(3)步中已經(jīng)驗(yàn)證每種環(huán)境下的實(shí)驗(yàn)數(shù)據(jù)的是否符合正態(tài)分布,如圖6(a)所示,由圖6(b)可見(jiàn)合并后的實(shí)驗(yàn)數(shù)據(jù)依然符合正態(tài)分布。然后,用每種環(huán)境下的樣本均值對(duì)數(shù)據(jù)進(jìn)行歸一化,求出歸一化后合并數(shù)據(jù)的平均值、標(biāo)準(zhǔn)差和離散系數(shù)見(jiàn)表1。

    根據(jù)公式(5)可以求出相應(yīng)的環(huán)境下的折減系數(shù)Bj,并與不同環(huán)境下的樣本數(shù)據(jù)平均值相乘,從而獲得B基準(zhǔn)值。

    圖6 數(shù)據(jù)的正態(tài)分布檢驗(yàn)

    為了驗(yàn)證程序的可靠性,采用標(biāo)準(zhǔn)中提供的范例中的計(jì)算結(jié)果進(jìn)行對(duì)比說(shuō)明。表2中列出了在四種不同實(shí)驗(yàn)環(huán)境下的程序計(jì)算結(jié)果和航標(biāo)HB7618?2009提供的結(jié)果。從表中的數(shù)據(jù)對(duì)比可以發(fā)現(xiàn)誤差在+0.3%以內(nèi),誤差主要由容限系數(shù)和折減系數(shù)的近似計(jì)算造成??梢?jiàn)程序計(jì)算的結(jié)果符合工程精度(誤差±5%)的要求,此外相比于借助Excel的手動(dòng)計(jì)算法方法,效率顯著提高。

    4 結(jié) 語(yǔ)

    本文根據(jù)航空工業(yè)標(biāo)準(zhǔn)HB7618?2009給出的B基準(zhǔn)值的計(jì)算步驟,將Matlab的GUI模塊與B基準(zhǔn)值的計(jì)算相結(jié)合開發(fā)出B基準(zhǔn)值的計(jì)算程序。相比于借助Excle的手動(dòng)計(jì)算方法,該程序界面具有以下優(yōu)點(diǎn):

    (1) 程序可以判斷樣本是否來(lái)自同一母體,還具有異常數(shù)據(jù)處理、統(tǒng)計(jì)分布檢驗(yàn)、離散系數(shù)檢驗(yàn)與修正等功能,最后實(shí)驗(yàn)結(jié)果以Word格式導(dǎo)出。

    (2) 使用該程序進(jìn)行B基準(zhǔn)值的計(jì)算方便快捷、效率高,并且具有良好的人機(jī)界面,顯著提高了B基準(zhǔn)值的計(jì)算效率。

    表2 不同環(huán)境下計(jì)算結(jié)果對(duì)比

    參考文獻(xiàn)

    [1] 張駿華,盛祖銘,孫繼同.復(fù)合材料結(jié)構(gòu)設(shè)計(jì)指南[M].北京:北京宇航出版社,1999.

    [2] 李威,郭權(quán)鋒.碳纖維復(fù)合材料在航天領(lǐng)域的應(yīng)用[J].中國(guó)光學(xué),2011,4(3):201?212.

    [3] 陳紹杰.復(fù)合材料設(shè)計(jì)手冊(cè)[M].北京:航空工業(yè)出版社,1990.

    [4] 劉建華,曹東,張曉云,等.樹脂基復(fù)合材料T300/5405的吸濕性能及濕熱環(huán)境對(duì)力學(xué)性能的影響[J].航空材料學(xué)報(bào),2010,30(4):75?80.

    [5] 呂小軍,張琦,馬兆慶,等.濕熱老化對(duì)碳纖維/環(huán)氧樹脂基復(fù)合材料力學(xué)性能影響研究[J].材料工程,2005(11):50?53.

    [6] 鄭錫濤,李野,劉海燕,等.濕熱譜老化對(duì)復(fù)合材料層壓板強(qiáng)度的影響[J].航空學(xué)報(bào),1998,19(4):462?465.

    [7] 中國(guó)航空工業(yè)總公司.HB 7618?2009 聚合物基復(fù)合材料力學(xué)性能數(shù)據(jù)表達(dá)準(zhǔn)則[S].北京:中國(guó)航空工業(yè)總公司,2009:50?85.

    [8] 陳子為.基于 Matlab GUI 掃雷游戲的設(shè)計(jì)與實(shí)現(xiàn)[J].現(xiàn)代電子技術(shù),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] 孫堅(jiān)石,葉強(qiáng).復(fù)合材料力學(xué)性能數(shù)據(jù)B基準(zhǔn)值計(jì)算程序[J]. 航空制造技術(shù),2009(z1):19?21.

    制服诱惑二区| 日本一区二区免费在线视频| 18禁黄网站禁片午夜丰满| 国产片内射在线| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看影片大全网站| 亚洲国产欧洲综合997久久, | 成人一区二区视频在线观看| 久久香蕉精品热| 国产精品香港三级国产av潘金莲| 人成视频在线观看免费观看| 免费在线观看成人毛片| 中文资源天堂在线| 国产亚洲精品第一综合不卡| 90打野战视频偷拍视频| 99热6这里只有精品| 人妻丰满熟妇av一区二区三区| 亚洲成人免费电影在线观看| 91麻豆精品激情在线观看国产| 国产精品精品国产色婷婷| 久久人妻福利社区极品人妻图片| 中文字幕人妻熟女乱码| 久久国产精品人妻蜜桃| 人人妻人人看人人澡| 中文亚洲av片在线观看爽| 久久婷婷人人爽人人干人人爱| 亚洲av第一区精品v没综合| 国产熟女午夜一区二区三区| 日日爽夜夜爽网站| 伊人久久大香线蕉亚洲五| www.999成人在线观看| 丝袜人妻中文字幕| 国产成人av激情在线播放| 精品不卡国产一区二区三区| 在线永久观看黄色视频| 两性夫妻黄色片| 制服人妻中文乱码| 精品少妇一区二区三区视频日本电影| 久久精品国产清高在天天线| 成人手机av| 国产又黄又爽又无遮挡在线| 又紧又爽又黄一区二区| 在线免费观看的www视频| 视频在线观看一区二区三区| 制服人妻中文乱码| 麻豆久久精品国产亚洲av| 亚洲第一青青草原| 别揉我奶头~嗯~啊~动态视频| 欧美一级a爱片免费观看看 | 色综合站精品国产| 美女 人体艺术 gogo| 免费一级毛片在线播放高清视频| 国产av在哪里看| 女性生殖器流出的白浆| 在线看三级毛片| 国产成人精品久久二区二区免费| 国产av在哪里看| 欧美性猛交╳xxx乱大交人| 操出白浆在线播放| 国产亚洲精品综合一区在线观看 | 琪琪午夜伦伦电影理论片6080| 欧美日韩福利视频一区二区| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 90打野战视频偷拍视频| 欧美成人免费av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲第一av免费看| 俺也久久电影网| 丰满人妻熟妇乱又伦精品不卡| 国语自产精品视频在线第100页| 男人舔女人下体高潮全视频| 国产色视频综合| 精品第一国产精品| 国产一级毛片七仙女欲春2 | 国产乱人伦免费视频| 国产精品香港三级国产av潘金莲| 91国产中文字幕| 夜夜躁狠狠躁天天躁| 国产一区二区三区视频了| 一二三四社区在线视频社区8| 亚洲色图 男人天堂 中文字幕| 黄色a级毛片大全视频| 亚洲精品国产区一区二| 老汉色av国产亚洲站长工具| 精品熟女少妇八av免费久了| 亚洲专区字幕在线| 两个人免费观看高清视频| 色精品久久人妻99蜜桃| 脱女人内裤的视频| 亚洲 国产 在线| 精品久久久久久久毛片微露脸| 精品无人区乱码1区二区| 亚洲一区二区三区不卡视频| 久久亚洲精品不卡| 十八禁人妻一区二区| 黑人欧美特级aaaaaa片| 国产精品影院久久| 亚洲一码二码三码区别大吗| 国产aⅴ精品一区二区三区波| 亚洲一区二区三区色噜噜| 国产午夜精品久久久久久| 一进一出抽搐动态| 91九色精品人成在线观看| 国产成人系列免费观看| 欧美黑人巨大hd| 国产又黄又爽又无遮挡在线| 婷婷精品国产亚洲av在线| 久久中文字幕一级| 丰满人妻熟妇乱又伦精品不卡| 日本黄色视频三级网站网址| 亚洲av电影在线进入| 国产视频内射| 侵犯人妻中文字幕一二三四区| 免费电影在线观看免费观看| 很黄的视频免费| 亚洲专区字幕在线| 在线免费观看的www视频| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人一区二区三| 美女大奶头视频| 男女视频在线观看网站免费 | 精品国产超薄肉色丝袜足j| 天堂√8在线中文| 日本一区二区免费在线视频| 日韩精品中文字幕看吧| 国产精品二区激情视频| 免费人成视频x8x8入口观看| 色播亚洲综合网| 999精品在线视频| 丁香欧美五月| 麻豆成人av在线观看| 老熟妇仑乱视频hdxx| 精品国产美女av久久久久小说| 亚洲精品美女久久久久99蜜臀| 久久精品国产亚洲av高清一级| 欧美色欧美亚洲另类二区| 亚洲av中文字字幕乱码综合 | 国产黄色小视频在线观看| 国产色视频综合| 男女那种视频在线观看| 动漫黄色视频在线观看| 国产一卡二卡三卡精品| 两个人视频免费观看高清| 亚洲欧美激情综合另类| 国产又黄又爽又无遮挡在线| 在线十欧美十亚洲十日本专区| 国产亚洲精品久久久久久毛片| 天天躁夜夜躁狠狠躁躁| 美女午夜性视频免费| 国产黄a三级三级三级人| 色播亚洲综合网| 99久久精品国产亚洲精品| 国内精品久久久久精免费| 日韩精品中文字幕看吧| 老汉色∧v一级毛片| 欧美在线黄色| 91成年电影在线观看| 91国产中文字幕| 婷婷丁香在线五月| 国产伦在线观看视频一区| 麻豆久久精品国产亚洲av| 亚洲精品色激情综合| 欧美成人一区二区免费高清观看 | 观看免费一级毛片| 久久热在线av| 在线播放国产精品三级| 黑人操中国人逼视频| 欧美日韩黄片免| 亚洲精品久久成人aⅴ小说| 香蕉丝袜av| 欧美日韩福利视频一区二区| 伦理电影免费视频| 欧美日韩瑟瑟在线播放| 欧美+亚洲+日韩+国产| 美女高潮喷水抽搐中文字幕| 亚洲黑人精品在线| 午夜久久久在线观看| 国产亚洲精品综合一区在线观看 | 婷婷精品国产亚洲av在线| 长腿黑丝高跟| 亚洲精品国产一区二区精华液| 身体一侧抽搐| 久久久久久久精品吃奶| 在线国产一区二区在线| 亚洲,欧美精品.| 三级毛片av免费| 曰老女人黄片| 亚洲无线在线观看| svipshipincom国产片| 一区二区三区高清视频在线| 麻豆av在线久日| 中文字幕av电影在线播放| 国产又爽黄色视频| 午夜激情福利司机影院| 黑人巨大精品欧美一区二区mp4| videosex国产| 九色国产91popny在线| 亚洲av成人不卡在线观看播放网| 丝袜人妻中文字幕| www日本在线高清视频| 99精品欧美一区二区三区四区| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费| av超薄肉色丝袜交足视频| 麻豆一二三区av精品| 久久亚洲精品不卡| 欧美黄色淫秽网站| 亚洲全国av大片| 免费人成视频x8x8入口观看| 91九色精品人成在线观看| 国产精品,欧美在线| av片东京热男人的天堂| 亚洲国产欧美日韩在线播放| 久久天躁狠狠躁夜夜2o2o| 午夜日韩欧美国产| 十八禁网站免费在线| 亚洲欧美精品综合久久99| 国产成人欧美| 俺也久久电影网| av电影中文网址| 国产私拍福利视频在线观看| 色av中文字幕| 婷婷精品国产亚洲av在线| 十八禁人妻一区二区| 国产伦一二天堂av在线观看| 色综合站精品国产| 中文字幕最新亚洲高清| 精品国产一区二区三区四区第35| aaaaa片日本免费| 日韩中文字幕欧美一区二区| 精品免费久久久久久久清纯| 国产成人精品久久二区二区免费| 亚洲国产欧洲综合997久久, | 一二三四社区在线视频社区8| 91麻豆av在线| 亚洲真实伦在线观看| 午夜福利在线观看吧| 真人做人爱边吃奶动态| 十八禁网站免费在线| 亚洲avbb在线观看| 久久人人精品亚洲av| 日韩av在线大香蕉| 亚洲中文日韩欧美视频| 亚洲精品美女久久av网站| 欧美另类亚洲清纯唯美| 亚洲国产精品999在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲熟女毛片儿| 51午夜福利影视在线观看| 窝窝影院91人妻| 亚洲片人在线观看| 国内精品久久久久久久电影| 在线观看免费视频日本深夜| 少妇裸体淫交视频免费看高清 | 天天一区二区日本电影三级| 国产精品爽爽va在线观看网站 | 日本 av在线| 午夜激情福利司机影院| 国产精品综合久久久久久久免费| 日本 欧美在线| 我的亚洲天堂| 嫩草影视91久久| 欧美成人一区二区免费高清观看 | 亚洲人成网站高清观看| 久久亚洲精品不卡| 欧美日韩亚洲国产一区二区在线观看| 午夜福利18| 欧美不卡视频在线免费观看 | 亚洲成av人片免费观看| 麻豆成人av在线观看| 国产伦人伦偷精品视频| 91在线观看av| 国内揄拍国产精品人妻在线 | 天堂动漫精品| 国产高清有码在线观看视频 | 亚洲av第一区精品v没综合| 国产97色在线日韩免费| 欧美精品啪啪一区二区三区| 男人舔女人下体高潮全视频| 夜夜爽天天搞| 法律面前人人平等表现在哪些方面| 亚洲va日本ⅴa欧美va伊人久久| 丁香欧美五月| 人成视频在线观看免费观看| 精品久久久久久久久久久久久 | 19禁男女啪啪无遮挡网站| 热99re8久久精品国产| 美国免费a级毛片| 久久久国产欧美日韩av| 欧美黑人巨大hd| 国产精品一区二区免费欧美| 亚洲国产高清在线一区二区三 | 香蕉国产在线看| 中文字幕高清在线视频| 无遮挡黄片免费观看| 日本a在线网址| xxx96com| 日本免费一区二区三区高清不卡| 日日干狠狠操夜夜爽| 国产人伦9x9x在线观看| 亚洲成人免费电影在线观看| 欧美性猛交黑人性爽| 性欧美人与动物交配| 日本熟妇午夜| 国产精品精品国产色婷婷| 国产精品一区二区免费欧美| 国产男靠女视频免费网站| 欧美日韩一级在线毛片| 日韩欧美 国产精品| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 亚洲性夜色夜夜综合| 成人18禁高潮啪啪吃奶动态图| 国产男靠女视频免费网站| 亚洲黑人精品在线| 国产激情欧美一区二区| 久久久久久久久久黄片| 亚洲在线自拍视频| 国产精品一区二区三区四区久久 | 美女大奶头视频| 欧美色欧美亚洲另类二区| 精品电影一区二区在线| 黄色片一级片一级黄色片| 妹子高潮喷水视频| 久久午夜亚洲精品久久| 精品一区二区三区视频在线观看免费| 香蕉国产在线看| 亚洲中文字幕日韩| 国产精华一区二区三区| 国产一级毛片七仙女欲春2 | 亚洲欧美日韩无卡精品| 在线av久久热| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| 午夜福利高清视频| 18禁黄网站禁片午夜丰满| 12—13女人毛片做爰片一| 两性夫妻黄色片| 日本黄色视频三级网站网址| tocl精华| 久久久久久九九精品二区国产 | 久久 成人 亚洲| 老汉色av国产亚洲站长工具| 最新美女视频免费是黄的| 波多野结衣高清无吗| 1024视频免费在线观看| 在线十欧美十亚洲十日本专区| 制服丝袜大香蕉在线| av电影中文网址| 俄罗斯特黄特色一大片| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻,人人澡人人爽秒播| 精品熟女少妇八av免费久了| www.自偷自拍.com| 国产区一区二久久| 成人三级黄色视频| 午夜福利成人在线免费观看| 久久狼人影院| 日韩欧美在线二视频| 国产高清激情床上av| 欧洲精品卡2卡3卡4卡5卡区| 国产主播在线观看一区二区| 国产av一区在线观看免费| 国内揄拍国产精品人妻在线 | 国产av在哪里看| av片东京热男人的天堂| 美女免费视频网站| 法律面前人人平等表现在哪些方面| 一区福利在线观看| 午夜久久久久精精品| 国产欧美日韩一区二区三| 亚洲天堂国产精品一区在线| 男女床上黄色一级片免费看| 成人一区二区视频在线观看| 国产精品爽爽va在线观看网站 | 欧美激情极品国产一区二区三区| 人人妻人人看人人澡| 亚洲人成77777在线视频| 亚洲国产欧美网| 欧美性猛交╳xxx乱大交人| 香蕉久久夜色| 亚洲av熟女| 成年版毛片免费区| 日本a在线网址| 欧美日韩瑟瑟在线播放| 1024视频免费在线观看| 亚洲中文字幕一区二区三区有码在线看 | 黄片小视频在线播放| 观看免费一级毛片| 亚洲精品粉嫩美女一区| 亚洲午夜理论影院| 国产亚洲精品久久久久5区| 少妇熟女aⅴ在线视频| 日韩av在线大香蕉| 美国免费a级毛片| 欧美午夜高清在线| 久久天躁狠狠躁夜夜2o2o| 国产一区二区三区在线臀色熟女| www.熟女人妻精品国产| av天堂在线播放| 午夜福利在线观看吧| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 亚洲 欧美一区二区三区| 欧美黄色片欧美黄色片| 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 国产黄色小视频在线观看| 人妻丰满熟妇av一区二区三区| 一区二区日韩欧美中文字幕| xxxwww97欧美| 12—13女人毛片做爰片一| 一区二区三区精品91| 免费一级毛片在线播放高清视频| 国产三级在线视频| 欧美日韩福利视频一区二区| 亚洲成人久久性| 最近在线观看免费完整版| 岛国在线观看网站| 国产一区二区在线av高清观看| 精品国产乱子伦一区二区三区| 亚洲精品久久成人aⅴ小说| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| a级毛片在线看网站| 老司机在亚洲福利影院| 欧美激情高清一区二区三区| 国内揄拍国产精品人妻在线 | 99热这里只有精品一区 | 日韩高清综合在线| 侵犯人妻中文字幕一二三四区| 亚洲七黄色美女视频| 国产99久久九九免费精品| 国产蜜桃级精品一区二区三区| 特大巨黑吊av在线直播 | 一级作爱视频免费观看| 法律面前人人平等表现在哪些方面| 亚洲 欧美 日韩 在线 免费| 日本免费a在线| 黄色a级毛片大全视频| 日韩欧美一区二区三区在线观看| 久久久国产成人精品二区| 伊人久久大香线蕉亚洲五| 1024手机看黄色片| 欧美成人午夜精品| 极品教师在线免费播放| 又黄又粗又硬又大视频| 视频区欧美日本亚洲| xxx96com| 亚洲国产精品999在线| 欧美成狂野欧美在线观看| 熟女电影av网| 丝袜人妻中文字幕| 精品一区二区三区视频在线观看免费| 久久久久久久午夜电影| 欧美一级a爱片免费观看看 | 久久国产精品男人的天堂亚洲| 国产成人av教育| 亚洲国产欧美日韩在线播放| 国产亚洲av嫩草精品影院| 久热爱精品视频在线9| 国产亚洲精品久久久久5区| 亚洲中文日韩欧美视频| 久久精品影院6| 最近最新中文字幕大全电影3 | 欧美黑人巨大hd| 女人爽到高潮嗷嗷叫在线视频| 亚洲,欧美精品.| 一卡2卡三卡四卡精品乱码亚洲| 亚洲真实伦在线观看| 操出白浆在线播放| 欧美三级亚洲精品| 一级片免费观看大全| 老司机深夜福利视频在线观看| a级毛片a级免费在线| 中文在线观看免费www的网站 | 深夜精品福利| 老司机福利观看| 亚洲国产精品sss在线观看| 18禁黄网站禁片午夜丰满| 每晚都被弄得嗷嗷叫到高潮| 欧美乱妇无乱码| 午夜两性在线视频| 嫁个100分男人电影在线观看| 午夜两性在线视频| 久久久久九九精品影院| 国产精品一区二区精品视频观看| 一夜夜www| 99精品在免费线老司机午夜| 日本熟妇午夜| 亚洲一区二区三区色噜噜| av福利片在线| 欧美日本亚洲视频在线播放| bbb黄色大片| 神马国产精品三级电影在线观看 | 别揉我奶头~嗯~啊~动态视频| 国产午夜福利久久久久久| 给我免费播放毛片高清在线观看| 精品不卡国产一区二区三区| 高清毛片免费观看视频网站| 黄色成人免费大全| 亚洲五月天丁香| 欧美日韩精品网址| 久久亚洲真实| 天天一区二区日本电影三级| 免费女性裸体啪啪无遮挡网站| 香蕉久久夜色| av在线天堂中文字幕| 可以在线观看毛片的网站| 变态另类丝袜制服| 亚洲欧美日韩高清在线视频| 亚洲第一青青草原| 色尼玛亚洲综合影院| 一区二区日韩欧美中文字幕| 91麻豆av在线| 久久国产精品男人的天堂亚洲| 久久久国产成人精品二区| 久久国产乱子伦精品免费另类| 久久这里只有精品19| 午夜福利高清视频| 欧美日韩福利视频一区二区| 999久久久精品免费观看国产| 亚洲成av人片免费观看| 亚洲电影在线观看av| 成人手机av| 精品福利观看| 啪啪无遮挡十八禁网站| 久久久久国产精品人妻aⅴ院| 亚洲av第一区精品v没综合| 欧美性长视频在线观看| 两个人视频免费观看高清| 搡老熟女国产l中国老女人| 日本免费a在线| www日本在线高清视频| 国产精品久久久久久人妻精品电影| 色老头精品视频在线观看| 久久精品91无色码中文字幕| 亚洲全国av大片| www.精华液| 成人亚洲精品一区在线观看| 一级黄色大片毛片| 亚洲成a人片在线一区二区| 久久中文看片网| 欧美乱色亚洲激情| 禁无遮挡网站| www日本黄色视频网| 国产99久久九九免费精品| 99国产精品一区二区蜜桃av| 日韩免费av在线播放| svipshipincom国产片| 国产精品爽爽va在线观看网站 | 国产免费av片在线观看野外av| 日韩欧美国产在线观看| 亚洲中文av在线| 婷婷亚洲欧美| 国产麻豆成人av免费视频| 免费av毛片视频| 夜夜躁狠狠躁天天躁| 日韩三级视频一区二区三区| 亚洲午夜理论影院| 一本久久中文字幕| 一级a爱片免费观看的视频| 99精品在免费线老司机午夜| 国产精品 国内视频| 精品人妻1区二区| 亚洲人成网站高清观看| 成人午夜高清在线视频 | 国产亚洲欧美精品永久| 国产黄a三级三级三级人| 亚洲色图av天堂| 国内精品久久久久久久电影| 岛国视频午夜一区免费看| 在线十欧美十亚洲十日本专区| av福利片在线| 精品国产美女av久久久久小说| 国产精品亚洲美女久久久| 最近在线观看免费完整版| 欧美黑人精品巨大| 中文资源天堂在线| 精品欧美国产一区二区三| av欧美777| 欧美日本视频| 欧美色视频一区免费| 国产高清激情床上av| 亚洲av中文字字幕乱码综合 | 欧美绝顶高潮抽搐喷水| 国产主播在线观看一区二区| 亚洲aⅴ乱码一区二区在线播放 | 中文字幕av电影在线播放| 精品国产乱码久久久久久男人| 麻豆久久精品国产亚洲av| 国产久久久一区二区三区| 国产精品日韩av在线免费观看| 亚洲 欧美 日韩 在线 免费| 久久草成人影院| 亚洲中文av在线| 亚洲成人免费电影在线观看| 中出人妻视频一区二区| 欧美激情极品国产一区二区三区| 在线观看66精品国产| 免费看日本二区| 亚洲精品粉嫩美女一区| 亚洲人成网站在线播放欧美日韩| 国产v大片淫在线免费观看| 女性生殖器流出的白浆| 一区二区三区激情视频| 美女扒开内裤让男人捅视频| 欧美性猛交╳xxx乱大交人| 免费在线观看影片大全网站| 国产成人啪精品午夜网站| 亚洲九九香蕉| 午夜福利高清视频| 精品第一国产精品| 啪啪无遮挡十八禁网站|