全的田口方法资料REV

上传人:枫** 文档编号:591998067 上传时间:2024-09-19 格式:PPT 页数:204 大小:1.65MB
返回 下载 相关 举报
全的田口方法资料REV_第1页
第1页 / 共204页
全的田口方法资料REV_第2页
第2页 / 共204页
全的田口方法资料REV_第3页
第3页 / 共204页
全的田口方法资料REV_第4页
第4页 / 共204页
全的田口方法资料REV_第5页
第5页 / 共204页
点击查看更多>>
资源描述

《全的田口方法资料REV》由会员分享,可在线阅读,更多相关《全的田口方法资料REV(204页珍藏版)》请在金锄头文库上搜索。

1、實驗設計田口方法1實驗設計簡介2為什麼需要實驗設計同樣在生產同規格的產品,為什麼有些廠商的良品率就是比較高。同樣是在生產同類型的產品,為什麼有些人的產品性能以及壽命就是比較好,而成本又比較低呢?This is the Know how相同原料相同製程為什麼良品率不一樣?相同產品相同功能更便宜的原料為什麼可以做出低成本高質量的產品?3日本工業強盛的原因日本人在多種製造業,如汽車、鋼鐵、電子和紡織方面,居於領導地位,主要是因為他們能以具競爭力的價格,生產高品質產品。美國研究後認為而他們致勝的法寶主要有二項:QFD(自顧客要求一直策劃到相應的製造管理要求)田口方法(實驗設計方法之一,簡單易用,沒有複

2、雜的統計原理)4實驗設計的想法processoutputinputsuppliercustomer客戶關心什麼,在乎什麼,抱怨什麼如何確定成為量化的產品特性?Y特性,CTQ過程中有那些的過程因子,其會影響到y特性,那些可能有相應的交互作用,從中找出CTP有那些的輸入因子會影響到y,從中找出CTI對供應應要提出那些的規格要求,尤其是CTI5Y=f(x)的思想過程產品原材料客戶供應商SupplierInputCTIxProcessCTPxOutputCTQyCustomerVOC6Y=f(x)Y=F(x)關鍵關鍵x:利用其控制利用其控制y的的平均值滿足需求,標准平均值滿足需求,標准差滿足需要差滿足

3、需要非關鍵非關鍵x:由於其影響較由於其影響較不大,用其來降低成本不大,用其來降低成本Y:優先關注客戶重視的優先關注客戶重視的特性,要能量化。特性,要能量化。有時不只一個有時不只一個Y特性特性7二階段的實驗步驟階段一階段一:篩選試驗篩選試驗決定決定y把有可能影響到y的x都要考慮並做實驗,以挑選出關鍵x階段二階段二:最佳條件最佳條件決定決定y針對已挑出的關鍵x,進行最佳條件的試驗,以決定最佳的x值。控制階段控制階段決定決定y針對關鍵少數的x參數,進行持續的控制,以spc監控其穩定性。8線外品管系統設計參數設計公差設計產品設計設計品管製程設計技術品管線上品管診斷預測測量生產製造品管服務顧客品管供應商

4、雙贏夥伴顧客需求期望滿意DOE的應用階段9田口工程簡介由田口玄一博士所提出的一套實驗方法,它在工業上較具有實際應用性,是以生產力和成本效益,而非困難的統計為依歸。廠商現在必須致力於在生產前就使複雜的產品能達到高品質。減少變異亦即要有較大的再現性和可靠性,而最終目的就是要為製造商和消費者節省更多的成本。10討論題實驗設計的目的是為了什麼?實驗設計是線上品管還是線外品管?為什麼線外品管要比線上品管早做呢?11變異和雜音雜音因素就是使機能特性,如燃料效率、換檔壓力、磨耗和轉向力等偏離目標值的因素。雜音因素可分為三類:外部雜音產品使用時,因使用條件,如溫度、濕度、灰塵等而使機能發生變異,此類條件為外部

5、雜章因素。內部雜音(劣化)產品組件的劣化。產品間雜音在既定的製造條件下,因條件變異而造成產品間發生差異。所有品質管制活動的目標就是要生產經得起各種雜音因素考驗的產品。12堅耐性所有品質管制活動的最終目標就是要生產經得起各種雜音因素考驗的產品。堅耐性堅耐性(Robustness)(Robustness)就是產品的機能特性就是產品的機能特性對雜音因素的差異不敏感,不受影響。對雜音因素的差異不敏感,不受影響。13堅耐性和雜音間的關係PROCESSINPUTOUTPUTYCUSTOMER可控因子可控因子X表示表示不可控因子不可控因子Z表示表示ZZZZXXX堅耐性設計,利用堅耐性設計,利用X使得使得Y達

6、目達目標值,且不易受標值,且不易受Z因子方影響,因子方影響,謂堅耐性設計謂堅耐性設計14討論題當你在讀書時,外界有人在吵,有人在放音樂,請問這些是信號,還是雜音?什麼是雜音,可否用比較簡易的方式來表達?什麼是堅耐性,如果以此題來表達時?15正交實驗法導入16一個瓷磚工廠的實驗在1953年,日本一個中等規模的瓷磚製造公司,花了200萬元,從西德買來一座新的隧道,窯本身有80公尺長,窯內有一部搬運平台車,上面堆疊著幾層瓷磚,沿著軌道緩慢移動,讓瓷磚承受燒烤。問題是,這些瓷磚尺寸大小的變異,他們發現外層瓷磚,有50%以上超出規格,則正好符合規格。引起瓷磚尺寸的變異,很明顯地在製程中,是一個雜音因素。

7、解決問題,使得溫度分佈更均勻,只要重新設計整個窯就可以了,但需要額外再花50萬元,投資相當大。17內部瓷磚外層瓷磚(尺寸大小有變異)上限下限尺寸大小改善前改善前外部瓷磚內部瓷磚18原材料粉碎及混合成型燒成上釉 燒成控制因素水准一(新案)水准二(現行)A:石灰石量5%1%B:某添加物粗細度細粗C:蠟石量53%43%D:蠟石種類新案組合現行組合E:原材料加料量1300公斤1200公斤浪費料回收量0%4%長石量0%5%19實驗方法一次一個因素法每次只改變一個因子,而其他因子保持固定。但它的缺點是不能保證結果的再現性,尤其是當有交互作用時。例如在進行A1和A2的比較時,必須考慮到其他因子,但目前的方法

8、無法達成。20一次一因素的實驗實驗次數ABCDEFG實驗結果1A1B1C1D1E1F1G112A2B1C1D1E1F1G123A2B2C1D1E1F1G134A2B2C2D1E1F1G145A2B2C2D2E1F1G156A2B2C2D2E2F1G167A2B2C2D2E2F2G178A2B2C2D2E2F2G2821全因子實驗法全因子實驗法這種實驗方法,所有可能的組合都必須加以深究。但相當耗費時間、金錢,例如7因子,2水準共須做128次實驗。13因子,3水準就必須做了1,594,323次實驗,如果每個實驗花3分鐘,每天8小時,一年250個工作天,共須做40年的時間。22A(64)B(32)C

9、(16)D(8)E(4)F(2)G(1)結果111111112111111231111121411111225111121161111212711112218111122291112111101112112111112121121112122131112211.1272222221128222222223正交表(Orthogonal Array)直交表(正交表)直交表用於實驗計劃,它的建構,允許每一個因素的效果,可以在數學上,獨立予以評估。可以有效降低實驗次數,進而節省時間、金錢而且又可以得到相當好的結果。24次數ABCDEFG結果123456711111111Y121112222Y23122

10、1122Y341222211Y452121212Y562122121Y672211221Y782212112Y825L8直交表A石灰石量B粗細度C蠟石量D蠟石種類E加料量F浪費回收G長石量每百件尺寸缺陷數ABCD EFG12345671234567111111115粗43現13000016211122225粗43新12004517312211225細53現13004512412222115細53新1200006521212121粗53現1200056621221211粗53新13004068722112211細43現12004042822121121細43新1300052626要素不良總數不良

11、百分比要素不良總數不良百分比A151/40012.75E112230.50A214235.5E27117.75B110726.75F15413.50B28621.5F213934.75C110125.25G113233.00C29223.00G26115.25D17619.00合計19324.12D211729.25回應表(Response Table)27最佳條件確認由於缺陷是愈小愈好,所以依此選出的最佳條件為:A1B2C2D1E2F1G2。確認實驗:將預期的缺陷數和“確認實驗”的結果做比較。但事實上廠商選得是A1B2C1D1E2F1G2,主要的原因是C(蠟石)要因的價格很貴,但改善的效果又

12、不大,所以選C1(蠟石含量為43%)28內部瓷磚外層瓷磚(尺寸大小有變異)上限下限尺寸大小改善前外部瓷磚內部瓷磚改善後29討論題從本案例中,你認為?最能提供最完整的實驗數據的是那一個方法一次一個因子法全因子法正交實驗法正交實驗法有何優點?30正交實驗配置31直交表和線點圖傳統的實驗計劃方法是由英國的在本世紀初發出來的,該方法包含多種的統計設計技巧,其需要使用比較繁複的統計技巧,所以較少使用在工業界上。田口方法:由田口玄博士所提出,它刪除許多統計設計的工作,以一種可以直接、經濟的方式一次就可以做許多因素的實驗,所以工業界上較常用。32單因素實驗所謂單因素實驗方法,即實驗過程中只允許單一因素變動,

13、其餘因素必須保持固定的實驗方法,但單因素實驗法可能之問題如下將會漏失複合因素的訊息當因素數很多時將較不真實無法保證在實際的製造條件中,可獲得實驗結果的再現性見次頁33一次一因素的實驗實驗次數ABCDEFG實驗結果1A1B1C1D1E1F1G112A2B1C1D1E1F1G123A2B2C1D1E1F1G134A2B2C2D1E1F1G145A2B2C2D2E1F1G156A2B2C2D2E2F1G167A2B2C2D2E2F2G178A2B2C2D2E2F2G2834全因素實驗全因素實驗計劃方法允許在同一時間內檢查多數因素的效果,而在做全因素實驗時所有因素的複合水準都將被檢查。35A(64)B

14、(32)C(16)D(8)E(4)F(2)G(1)結果111111112111111231111121411111225111121161111212711112218111122291112111101112112111112121121112122131112211.1272222221128222222236直交性在實驗計劃中最主要的一個特性,便是實驗結果的再現性;另外,當我們希望能在各種相異的條件,以最有效的方式比較因素水準時,都只有在直交性實驗計畫方法中才能達到利用直交表進行實驗,在實驗結果的可靠度及高再現性上,都具有高效益。不管製程條件如何變化,在不同條件下,獲得好的再現性之效果是

15、相同的。假如我們的實驗計劃均為直交,則我們在回應表中比較A1和A2時,我們將可確定A1中B效果與A2中的B效果應為相同,且當因素以直交方式變動時,其它的效應將不會混合於各因素的水准內。37次數ABCDEFG結果123456711111111Y121112222Y231221122Y341222211Y452121212Y562122121Y672211221Y782212112Y838實驗編號ABC結果1111Y12112Y23121Y34122Y4實驗編號ABC結果1111Y12122Y23212Y34221Y4假設實驗執行所需花費的成本相當高,在此情況下不管任何理由,我們希望只做四次實驗,

16、以代替全因素實驗。請問下列二表,你會選擇那一項39自由度的概念自由度實為獲取情報大小的量度,通常自由度愈大,所獲得的情報愈多例子有三個人比較身高,至少須比較多少次才可以知道結果須比較二次才可以得到結果直覺上的定義:因素的自由度為水準間所必需但不重覆的比較次數,而在數理運算上,因素的自由度可簡單的以水準數減一表示,它代表因素能夠相互獨立記述計算的數目。在實驗中因素設定的水準愈多,則自由度隨著增加,換句話說可以得到更多情報,但是相應的實驗成本會增加。40目前有三個人的身高,如果要去進行比較,最少的比較次數,而得到全部的信息。41效果A1A2效果B1B2B3二水准的情況,只須比較一次,所以自由度為一

17、。三水准的情況,須比較二次,所以自由度為二。42交互作用原先假設因素的效果不會受其它因素水準的影響,然而在實際的狀況並非如此;當一個因素的效果與其它因素水準相互影響時,因素間就有交互作用存在。一般可以繪製交互作用圖來了解其間之交互作用關係。例子:設有A, B二種藥劑,成份完全不同,且兩者都能夠使病人狀況獲得改善;單獨使用時都有功效,但合併使用,病人反而更槽。43A和B無交互作用A和B有交互作用B1B2A1A2YB1B2A1A2Y44B2A1A2B1A和B有強烈交互作用B1B2A1Y1Y2A2Y3Y4 交互作用分析表45交絡在決定是否要配置交互作用效果於一直行時,要相當謹慎,必需於交互作用極端重

18、要才可進行配置。假如A因素和B因素間並無交互作,則行3將可能配置另一個因素C,此時由於配置因素C在A*B交互作用存在的行,我們將無法再由該行估計A*B的交互作用。如果我們的判斷是錯的,且A*B相當顯著,則交互作用效果將會顯現在該行的估計值中,但是我們將無法由C因素效果中,將交互作用效果區分出來,此種現象稱之為交絡交絡。46二因素交互作用的自由度數學上之表示如下d.f.(A*B)=d.f.(A)*d.f.(B)例子,A為二水準,B為三水準,則其自由度為(1)*(2)=2效果A2A1B1B2B347直交表的自由度(二水准)表示直交表表示直交表ROW數相當於實驗總數數相當於實驗總數水準數水準數COL

19、UMN數相當於可數相當於可配置多少因子配置多少因子直交表的自由為直交表的自由為實驗執行次數減一實驗執行次數減一48直交表的自由度(三水准)表示直交表表示直交表列數相當於實驗總數列數相當於實驗總數水準數水準數行數相當於可配行數相當於可配置多少因子置多少因子直交表的自由度为直交表的自由度为实验执行次数减一实验执行次数减一49練習在二水准的直交表中,a和b有何關係,如果因子依此配置有何關係?在三水准的直交表中,a和b有何關係,如果因子依此配置有何關係?50練習試寫出直交表L8(27)可提供多少自由度,最多可以配置幾個因子。試寫出直交表L9(34)可提供多少自由度,最多可以配置幾個因子。試寫出直交表L

20、81(340)可提供多少自由度,最多可以配置幾個因子。試寫出直交表L64(421)可提供多少自由度,最多可以配置幾個因子。511231111212232124221L4(23)直交表本直交表總共須做四次實驗,總共可提供三個自由度。每一個二水准的因子需要一個自由度,所以最多只能配置三個因子。52L8(27)直交表本直交表總共須做8次實驗,總共可提供7個自由度。每一個二水准的因子需要一個自由度,所以最多只能配置7個因子。如果有二水准因子間有交互作用時,交互作用亦須配置自由度。次數ABCDEFG結果123456711111111Y121112222Y231221122Y341222211Y45212

21、1212Y562122121Y672211221Y782212112Y853直交表的運用利用自由度我們可選用最小且最合適的直交表,係依據因素數量、每個因素的水準數,以及我們所欲調查的交互作用數量等加以累加後實驗計畫的自由度來決定。例如:一實驗包含二水準因素A、B、C、D、E和交互作用A*B,A*C,請問應選用何種直交表解決此一問題54每個二水準因素具有2-1=1的自由度。每個交互作用具有1*1=1的自由度總自由度d.f.=(5個因素*1d.f.)+(2交互作用*1d.f.)=7d.f.因此,7個自由度是獲得期望資料數量所必需的自由度,而L8直交表為二水準具7個自由度的實驗計畫,因此L8直交表是

22、可以滿足此項要求的。55兩行間交互作用的配置假如我們預期兩變數存在有顯著的交互作用,則我們可能在直交表中,預先保留一直行供配置交互作用,以利清晰的估計交互作用。如果我希望避開交絡現象,則必需妥慎的配置交互作用,不應任意配置;如果不加注意,則不管是最簡單的L4直交表,或是最複雜的直交表,交互作用的追蹤分析將變得困難。可以利用的方法是三角矩陣法。練習於行(3)和(7)的交互作應配置於那裡。561234567(1)325476(2)16745(3)7654(4)123(5)32(6)1(7)L8(27)直交表的交互作用配置表57練習假設目前是要找出關鍵因子,所以都用二水准實驗,請問該用那個直交表。如

23、果溫度和時間是有交互作用,如果溫度配在L8的第三行,時間配在第五行,那麼其交互作用行應配置在那裡?熱處理熱處理原材料原材料溫溫度度時時間間淬淬火火油油種種類類硬度硬度58線點圖線點圖為田口博士在運用直交表於實驗計畫上的另一個貢獻,所謂線點圖即為三角矩陣的圖示方法,它利用圓點與直線圖型為工具,以便利完成直交中與交互作用的配置。59線點圖之使用配置因素於圓點上考慮因素間的交互作用,若交互作用存在,則配置該交互作用於聯接該兩因素圓點的直線上。若某兩因素間的交互作並未確定存在,則該兩因素的聯線上可配置其他的因素。1325467123546760交互作用的考慮主效果二因子交互三因子交互四因子交互五因子交

24、互AABABCABCDABCDEBACABDABCECADABEABDEDAEACDACDEEBCACEBDADEBEBCDCDBCECEBDEDECDE一般工業上,研究交互作用並非實際,就算交互作用存在,也是不容易對付的,所以一般高層的交互不考慮。61直交表的因素配置步驟一計算實驗總自由度。二水準因素A、B、C、D、E,交互作用B*C,C*D,故自由度=5+2=7。一個二水準m行(column)的直交表,具有m個自由度,由各種直交表中選擇一個能夠包含實驗所需自由度的直交表。62步驟二:繪出所需要的線點圖。步驟三:由標準的線點圖中選擇適當的線點圖。步驟四:將繪出的線點圖與標準線點圖相比對配合,

25、選出最合適的線點圖。步驟五:配置每個主效果和交互作用到合適的行中。63練習如果現有A,B,C,D四個因子,其AB, CD, BC是有交互作用的。請試著利用點線圖進行實驗配置。如果利用L8直交表配置不出來時該如何處理?請利用L16表配配看,是否可以滿足?64有交互作用之直交表配置計算所有因子與交互作用自由度之和,即所有要因自由度總和。選取自由度不小於要因自由之和且試驗次數最小之相同水準數直交表。選定一交互作用,將其相關之兩因子任意配置於直交表之行上,然後根據“交互作用配行表”將此交互作用配置於直交表之行上重覆上個步驟,直至所有交互作用之因子皆配置完為止。65將剩餘之因子任意配置於直交表之剩餘行上

26、。依據各要因所配置行上的數字1,2或3,決定各試驗的水準組合。依隨機順序進行全部之試驗。66線點圖的練習二水准:A,B,C,D,E,F,G,H,I,J,交互作用AB, AC, BC, BG, GH。二水准:A,B,C,D,E,F,G,H,I,交互作用AB, AC, AD, AE, EF, EG 。二水准:A,B,C,D,E,F,G,H,交互作用AB, AC, AE, BC, GH 。67如何修改標準的線點圖利用三角矩陣表使一個特殊的實驗可以利用標準的線點圖,修改為合乎使用的線點圖。此項須配合練習做法繪製符合題意的線點圖選擇合適的標準線點圖檢查三角矩陣圖,以瞭解標準線點圖可被修改的方式修改標準線

27、點圖,並且配置合適的因子效果682641211011133598715142641211011133598715143修改利用三角矩陣使一個特殊的實驗可以利用標准線點圖型,修改為合乎使用的線點圖。此圖線3已被移動,原先的直線已刪除,即1,2的交互作用是3,但12,15的交互作用也是3可以移過來。69交互作用直交表練習利用L16直交表,配置下列實驗計劃A,B,C,D,E,F,G,H,I,JAB, CD, EF, FG, GHSTEP1:繪出需要的圖 IJABCDEFGH70STEP2a:選擇合適的標准線點圖132412851510714961311STEP2b:利用L16三角矩陣表,以了解可被修

28、的方式;由於行7與10的交互作用配置於行13,因此,可將線13移動,使其聯接點7與點10 IJACDEFGHB13241285151013714961171123456789101112131415(1)32547698111013121514(2)1674510118914151213(3)765411109815141312(4)12312131415891011(5)3213121514981110(6)114151213101189(7)15141312111098(8)1234567(9)325476(10)16745(11)7654(12)123(13)32(14)1(15)L16

29、三角矩陣表72線點圖的練習二水准:A,B,C,D,E,F,G,H交互作用AB, AC, AD, AE, BC, CD, EF。二水准:A,B,C,D,E,F,G,H,I,J,交互作用AB, CD, EF, EG, GH。二水准:A,B,C,D,E交互作用AD, BE。73直交表的數據分析傳統的變異數分析係以統計檢定的方法來決定各因素對變異影響程度是否顯著。田口方法卻強調以回應分析方法(回應表及回應圖)來區分各因素平均值效果的大小。正規分析決定因素間的平均回應值。比較這些回應平均值,並選出最佳因素水準。由選出的最佳水準來估計製程平均。確認實驗的結果與估計值比較。74次數ABABCACBCABC結

30、果123456711111111Y121112222Y231221122Y341222211Y452121212Y562122121Y672211221Y782212112Y8123, 145, 246, 167一般而言,三階交互作用的存在不明顯,可以假定該交互作用不存在,所以此行可以用來配置其他主要因子。但如果交互作用明顯則會產生交絡現象。75次數ABCDEFG結果123456711111111Y121112222Y231221122Y341222211Y452121212Y562122121Y672211221Y782212112Y8飽和配置當所有的COLUMN均配置主因子時,完全忽視交互

31、作用時,則稱為飽和配置76配置練習三水准,A, B, C, AC有交互作用,請進行配置。三水准,A, B, C, D, E, AD, BE請進行配置。77配置練習在塑膠射出成型的工廠中,研究影響塑膠製品強度的三個因子分別,A,B,C。是用三水准。如果A,B,C之間沒有存交互作用,則可利用那一個直交表配置?若A與B有交互作用,則用那一個直交表來配置?78直交表的數據分析計量型79直交表的數據分析正規分析決定每個因素的平均回應值。估計每個因素及交互作用之主效果。接著比較各主效果,找出較強之主效果。完成回應表交互作用的分析只需點繪較強之主效與交互作用,因為較弱效果的因素水準,對推動力的影響極微,可以

32、忽略不計。最佳化及最佳條件的估計。確認實驗80遊艇的真空控制閥門組合的推動力目標:推動力(望大特性)交互作用:BC與CD控制因素水准一水准二A:原材料M-270M-290B:停留時間2.7秒2.2秒C:焊鎗溫度410室溫D:拼裝壓力 60磅80磅E:焊接方法音波法電阻法81次數CBBCDCDAE數據合計123456711111111383573211122224046863122112257521094122221145551005212121241488962122121282452722112212625518221211235417682決定每個因素水准的平均回應值CBDAECBCD水准

33、146.0037.5040.2542.2534.535.7538.75水准233.5042.0039.2537.254543.7540.751-212.504.501510.58283B1B2C139.7552.25C235.2531.75交互作用的計算84推動力E1E21020304050推動力D1D21020304050推動力B1B21020304050推動力A1A21020304050推動力C1C21020304050推動力C1C21020304050B2B1回應圖Tbar85最佳化條件的選定因為推動力為望大特性,從CB交互作用之回應圖知因素B與因素C之最佳水准組合為C1B2。從回應圖上

34、看出:A1的效果不錯D1的效果尚不錯E2的效果最強所以最佳條件為C1B2D1A1E286最佳水准回應值之估計為確認所定之結果的再現性,必須再估計出此最佳條件C1B2D1A1E2的推定值,並與確認實驗的結果相驗證,看是否具有再現性。87推動力E1E21020304050推動力D1D21020304050推動力B1B21020304050推動力A1A21020304050推動力C1C21020304050推動力C1C21020304050B2B1效果可加性分析TTTTTTT88確認實驗確認實的目的是為了確認結果的再現性。即在最佳條件A1B2C1D1E2之下,做了一次確認實驗,本次實驗的製程平均為。

35、CASE1:Y=58,再現性非常好。CASE2:Y=54,沒有CASE1好,但仍算好的再現性。CASE3:Y=42,再現性差,但比38好,可以先用,然後考慮再改善。CASE4:Y=30,再現性差,不可接受,必須重新考慮。CASE5:Y=65,遠較所期望的還好,可能存在某種交互作用。89確認實驗結果之說明如果確認實驗之結果不佳時,一般原因如下:可加性極差,即所選取之控制因素有極強烈的交互作用存在。最佳條件所選取之控制因素仍不足,可能遺漏了一個極顯著的控制因素。因素水準距離太小,無法測得因素水準改變所造成的效果。再現性不好的措施水准距離是否設得太小。是否漏了重要因子。有交互作用沒有考慮到。90練習

36、磨耗率控制因素水准一水准二A:活塞間隙低高B:軸承間隙低高C:襯墊厚度厚薄D:調節板型態1型態2E:閥門設計型態1型態2交互作用:AB與AC望小特性91次數ABABCACDE數據合計12345671111111166621282111222268631313122112288801684122221163651285212121273711446212212137427972211221383977822121125748105 n=16 T=96092ABCDEABAC水准169.37560.2564.62563.12551.555.12560水准250.62559.7555.37556.8

37、7568.564.875601-218.750.59.256.25179.6450回應表B1B2A164.7574A255.7545.5交互作用的計算93E1E220406080100C1C220406080100B1B220406080100D1D220406080100A1A220406080100A1A220406080100B2B1回應圖94最佳條件選定最佳水准估計確認實驗實驗結果為30實驗結果為35實驗結果為4095練習為提高某型冷氣機之EER值,考量下列控制因子皆為2水准控制因子A壓縮機規格B散熱片型式C散熱片處理冷媒銅管型式毛細管長度試設計本實驗96實驗數據分析練習ABCDEee

38、123411111112.112.132.102.1411122221.982.012.052.1412211222.182.202.222.2312222112.002.052.032.0721212121.991.892.032.0021221212.062.082.102.0422112212.262.282.192.2422121122.182.162.132.1997直交表的數據分析1231111212232124221L4(23)直交表132做四次實驗,可配置三個因子,是最小的直交表98直交表的數據分析L8(27)直交表1234567111111112111222231221122

39、4122221152121212621221217221122182212112135264。7132546799直交表的數據分析L12直交表L12是一個非常特殊的直交表,交互作用的效果平均分配到該直交表的11個縱行上。它沒有線點圖可以使用。使用之前提在於交互作用並不明顯時。它的再現性很好,是田口博士所推薦使用的。1001234567891011111111111111211111222222311222111222412122122122512212212121612221221211721221122121821212221112921122212211102221111221211221

40、212111221222112121221L12(211)直交表101三水準系列直交表每一column可提供二個自由度。每個因子需佔用一行。三水準須使用二個自由度交互作用需佔用兩行(3-1) (3-1)=4。1021234111112122231333421235223162312731328321393321L9(34)直交表13,42103L18(2137)直交表此表可配置一個2水准與七個3水准。1+27=15,但事實上L18應是提供17個自由度。但實際上此表在第一行與第二行之間存在一個“內含”的交互作用,(2-1) (3-1)=2。在第一行和第二行之間可用配置表及回應圖將交互作用給檢查出

41、來。在AT&T,L18是最普遍被使用的直交表。最常使用的直交表為:L16, L18, L8, L27, L12。10412345678111111111211122222311333333412112233512223311612331122713121323813232131913313212102113322111212113321221322113132212313214222312131522312321162313231217232111231823323231L18(2137)直交表12105衝突時的解決方式此時要將y1,y2,y3都必須做分析,並針對彼此做出最佳的妥協。不然就可能

42、會出現,y1好了,但是y2可能變不好了。y1,y2,y3x1,x2x3x4x5x6106品質特性的選取田口方法係一種工程方法,擁有製程或產品的專門知識及有效率的實驗方法,才能夠設計出來一個極有效的工業實驗,因此必須懂得此兩種型態的知識才可能成功。品質特性的選取及因素與水準的區分是屬於工程專家的工作;而各因素的配置及實驗數據的解析則屬於數據分析專家的工作。田口博士視品質特性的選擇為實驗計畫中最主要的部份,也是最困難的部份。107总体方法实际问题統計問題統計解決方法實際解決方法108參數設計內外直交表109LSLUSL正規分析的目的110LSLUSL參數設計的做法LSLUSLSTEP1STEP21

43、11參數設計參數設計的目的,在於決定產品與製程的參數值,以求得產品機能的穩定,使其在高水準下運作,而受干擾的影響程度最低。參數設計在於運用因素間非線性與線性的一些關係,找出控制因素與誤差因素間的交互作用,利用非線性減少變異,再利用線性關係提高水準,即使使用便宜的材料或在不良的環境之下,製程或產品也能達到堅耐性。112參數設計所運用的技術是S/N比(訊號雜音比),它可以表示製程或產品的水準及其誤差因素影響的程度參數設計是一種提高品質而又不影響成本的設計,一般而言,要提高品質一定要把影響這個產品的不良原因消除,才能達到,如此則必須提高成本,如果不去消除原因,而把這些原因所產生的影響設法消除,則不必

44、花什麼成本,也能提高品質,此即參數設計。113參數設計的配置參數設計的第一步,為分開列出控制因素與誤差因素,然後找出具有最小交互作用的控制因素以便研究控制因素與誤差因素之間的交互作用問題。一般而言控制因素放在直交表內側,誤差因素放在直交表外側。114參數設計設置ABCDEFG誤差因素1234567N1 N21111111121112222312211224122221152121212621221217221122182212112115ABCDEFG112212345671212111111112111222231221122412222115212121262122121722112218

45、2212112誤差因素MN116誤差因素的選擇作參數設計時,雖然誤差因素愈多愈好,如此才可獲致較多情報,但實驗將會變得很大,在費用與時間將不允許,故只能在經營能力範圍之內,選擇重要的,影響較大的才予以考慮。對策為了避免太大的實驗,最好能使用誤差因素複合成1,2或最多3個。複合時可依工程知識做取捨,假如不能確知時,應事先用直交表做實驗,一定是選重要的,影響最大的。117選擇最重要的誤差,經驗告訴我們,試驗時若對最大的誤差具有堅耐性的話,對其他的誤差也必將穩定。一般採用2水準即可,並可用兩極端條件複合。118ABCDEFG11221212123456712211111111121112222312

46、211224122221152121212621221217221122182212112誤差因素MNO119信號雜音比望小特性(不包括負值,不包括不良率0%,最佳條件最理想狀態為0)。當品質特性能夠分類,而希望愈小愈好時,如產品的收縮度,劣化度、噪音、各種公害等,其標准的信號雜音計算如下120望小特性的S/N比121練習平均排序 平均平方偏差損失/單位S/NdbA10.200.200.200.2040.047.7713.98A20.100.200.300.2050.0479.1213.28A30.200.100.150.1520.024A40.100.100.100.1010.01020A5

47、0.100.300.050.1530.034假設,修理成本為70美元,求損失函數值。計算K值以及每一方法中單位產品的損失金額,本例產品的損失金額為25美元122練習123望小特性SN值的特性S/N值是量測平均值與變異程度的指標。S/N對平均值的靈敏度大過於對變異數的靈敏度。當S/N比每增加3分貝時,則單位損失將減少 一半。X分具的增益(Gain)可用S/N比數字來表 示,其公式為:124補充Why125126127練習:射出成型(望小特性參數設計)控制因素水准一水准二A:周期時間B:成型溫度C:凹處深度D:握住壓力E:旋轉速度F:握住時間G:閥門大小128練習:射出成型(望小特性參數設計)誤差

48、因素水准一水准二H:再磨光率%I:水分含量J:周圍溫度回應特性:收縮長度(吋)百分率用L8直交表排在内側,L4直交表排在外側。計算S/N比。完成回應表及回應圖。決定最佳因素組合。129參數設計配置試驗號碼ABCDEFGHIJ111122212221S/N111111112.22.12.32.3-6.95211122220.32.52.70.3-5.35312211220.53.10.42.8-6.50412222112.01.91.82.0-5.70521212123.03.13.03.0-9.62621221212.14.21.03.1-9.12722112214.01.94.62.2-10

49、.53822121122.01.91.91.8-5.58130S/N回應表根據計算的S/N比,可以製作回應表,並繪圖如下而選擇最佳的組合。因素水準ACBDEFG水准1-6.12-7.76-7.11-8.41-7.04-6.96-8.07水准2-8.72-7.09-7.74-6.44-7.81-7.88-6.76131A1A2-9-8-7-6B1B2C1C2D1D2E1E2-9-8-7-6F1F2G1G2最佳組合是 回應圖132練習:射出成型(望小特性參數設計)當S/N比計算出來後,即可求出回應表,繪出回應圖,並可選出最佳因素組合。(亦可將成本與品質因素並同考慮選擇)。133練習:射出成型(望小

50、特性參數設計)由此可知:增益每增加3db,則MSD可減少一半。增益每增加6db,則MSD可減至四分至一。增益每增加9db,則MSD可減至八分至一。MSD減少則平均值減少,變異程度亦減少。比較S/N時愈大愈好,MSD愈小愈好。134練習:泡沫膠膜收縮度(望小特性參數設計) 回應特徵:泡膜膠膜收縮度%控制因素水準一水準二A:脫膜時間現行新B:凝固溫度現行新C:化學反應溫度現行新D:混合比現行新E:外殼預熱溫度現行新F:脫模厚度現行新G:射口大小現行新H:脫模曲綫現行新I:射出壓力現行新J:通過瓶現行新K:密度現行新135練習:泡沫膠膜收縮度(望小特性參數設計)誤差因素: M:定型度(Fixture

51、)計算S/N比完成回應表,找出對泡沫膠收縮影響較強的因素。繪出影響較大的回應圖。136L121234567891011M1M2M3db1111111111110.130.180.1416.42111112222220.250.270.2212.13112221112220.080.120.1219.34121221221120.080.110.1119.95122122121210.170.220.1315.06122212212110.090.180.1915.97212211221210.240.320.2811.08212122211120.180.190.1115.7921122212

52、2110.190.200.3010222111122120.220.250.2211221212111220.240.180.1612221121212210.080.100.11配置與數據137回應表水準1水準2水準1水準2平均差合 計 平 均A98.686.416.4314.402.03dbB87.197.914.5216.321.80dbC95.389.715.8814.950.93dbD92.292.815.3715.470.10dbE82.3102.713.7217.123.40dbF99.685.416.6014.232.37dbG90.294.815.0315.800.77dbH

53、101.683.4I92.192.9J93.391.7K91.193.9138S/N回應圖A1A251015B1B2C1C2S/ND1D251015E1E2F1F2S/N139S/N回應圖S/NG1G251015H1H2I1I2J1J251015K1K2S/N140練習:泡沫膠膜收縮度(望小特性參數設計)最佳組合的推定從回應圖中吾人選出:A1B2E2F1H1為顯著性較強的因素。C1D2G2I2J1K2為顯著性次佳的因素。故最佳組合為:A1B2C1D2E2F1G2H1I2J1K2141練習:泡沫膠膜收縮度(望小特性參數設計)在最佳組合情況下回應特性值之預測爲了避免由於高估變異的誤差,吾人在計算回

54、應特性值時僅持用影響因素較高的因素來估算。142請問倒推回來的值應是多少?時那麼其收縮率應為多少?Y2千萬要記得再做確認實驗,看其再現性是否良好。143練習:泡沫膠膜收縮度(望小特性參數設計)假如在現行狀況時吾人獲得數據如下:0.160.130.160.180.130.160.180.200.160.130.140.150.090.160.140.140.140.160.120.150.130.150.130.14n=24T=0.147144練習:泡沫膠膜收縮度(望小特性參數設計)首先計算在現行狀況時S/N比: exist=現行狀況時之S/N比值145練習:泡沫膠膜收縮度(望小特性參數設計)其

55、次計算在現行狀況下之損失(假設k=194.4)146練習:泡沫膠膜收縮度(望小特性參數設計)現行狀況與最佳狀況下S/N比之差異,可計算出兩者閒之增益。147練習:泡沫膠膜收縮度(望小特性參數設計)是以吾人可估計,在最佳狀況下之損失,並且從而估出製程改善後最佳狀態下所能節省之金額。由以上計算得知,每件可節省3美元,如每月生產10,000件,則一年可節省360,000美元。148望大特性的S/N比149O1221平均MSD排序每件損失金額S/NdbM1212V1122A14.04.04.04.040.063 3$1.8012.04A22.06.02.06.040.139 4$4.008.57A34

56、.06.04.06.050.045 213.45A41.010.01.08.050.506 5A54.04.012.012.080.035 1A61.01.21.31.510.683 6A710.020.00.55.091.013 8A81.01.01.01.011.000 7150望大特性SN值的特性S/N值是量測平均值與變異程度的指標。S/N對平均值的靈敏度大過於對變異數的靈敏度。當S/N比每增加3分貝時,則單位損失將減少1/2。X分具的增益(Gain)可用S/N比數字來表 示,其公式為:151練習塑膠表面塗鋁處理回應特性:塑膠撕裂強度控制因素水準一水準二A:分子量低高B:金屬鋁型態C:清

57、淨材料現行新D:觸媒百分率低高E:研磨壓力低高F:處理時間低高G:處理溫度低高H:處理後時間低高I:處理後溫度低高152練習塑膠表面塗鋁處理計算S/N比完成S/N回應表,並找出對撕裂強度影響較大的因素。繪出影響較大因素之回應圖。誤差因素水準一水準二M:表面粗糙度粗糙平滑B:塑膠供應商C:最初金屬溫度低高153NO.ABCDEFGHIeeMNO111122212221S/N11111111111132.545.035.050.031.821111122222262.567.555.080.036.231122211122220.030.030.037.528.741212212211220.02

58、7.527.550.028.651221221212142.555.032.560.032.861222122121122.545.047.550.030.872122112212160.080.045.062.535.382121222111245.047.527.555.031.992112221221147.562.575.080.035.9102221111221255.062.555.047.5112212121112260.045.040.042.5122211212122145.035.020.030.0配置與數據154練習塑膠表面塗鋁處理155計算上表最後一欄倒數三行之S/N數

59、值,並完成以下回應表S/N比對應表因素水準ABCDEFGHI132.432.833.631.432.830.932.2232.432.131.233.432.033.932.6由二者之差選出最佳組合156回應圖S/NA1A2303234B1B2303234E1E2303234F1F2303234157回應圖S/NH1H2303234=-10log(MSD)L=k(MSD)雜音與損失函數成正比故可由此計算雜音與損失函數成正比故可由此計算158練習塑膠表面塗鋁處理最佳狀況的決定最佳組合是:A2B1E1F2H2在最佳狀況下回應特性值的預測爲了避免由於高估變異的誤差,吾人估計在最佳組合下撕裂強度特性值

60、時,用顯著性較強的影響因素來推計。159練習塑膠表面塗鋁處理你該如何計算由現行狀況改變為最佳狀況後所產生之增益呢? 增益 =o-e160練習如果原條件的值為30,那麼請回答下列問題。請問在本次實驗所得到的最佳化值後,其損失可以降為原來的多少?請將值轉換成強度的估計值,為多少?161練習塑膠表面塗鋁處理(續)傳統分析利用每次實驗所得數據來比較傳統分析及S/N之分析完成回應表並找出影響大的因素。繪出回應圖,並找出最佳組合狀況。估計最佳組合狀況時回應特性值及opt。162NO.ABCDEFGHIeeMNO111122212221計11111111111132.545.035.050.02111112

61、2222262.567.555.080.031122211122220.030.030.037.541212212211220.027.527.550.051221221212142.555.032.560.061222122121122.545.047.550.072122112212160.080.045.062.582121222111245.047.527.555.092112221221147.562.575.080.0102221111221255.062.555.047.5112212121112260.045.040.042.5122211212122145.035.020.0

62、30.0配置與數據163計算上表的平均值並完成以下回應表平均值回應表因素水準ABCDEFGHI12選出最佳組合:164平均值回應圖A1A240505545B1B240505545E1E240505545165平均值回應圖F1F240505545H1H240505545166練習所以最佳條件為?A2B1E1H2I2估計的的最佳平均值為?和利用S/N計算的結果有不同嗎?最後請記得做確認實驗。167傳統分析和S/N分析的比較傳統分析益處可求出控制因素與誤差因素之主效果假如我們比較誤差因素時知道第一部機器比第二部機器好,就可選擇第一部機器。限制只能在變異一樣的情況下比較平均值。誤差要有常態分配。受到交

63、互作用的影響。168傳統分析和S/N分析的比較益處:S/N著眼在平均值及變異數。交互作用有,但平均分配到各效果上。把交互作用視為誤差。S/N是測試因素與誤差因素間有無交互作用的指標。S/N所顯示的數據係干擾程度的指標。S/N與成本間有密切的關聯。簡單容易導出結論來。=m2/2表示變異大小的指標。/m稱變異係數越小越好。=m2/2愈大愈好169望目特性SN值的特性S/N比(db),是變異的變化的量度。靈敏度(Sm db),是平均值變化的量度。一般其狀況如下:望目特性的因子種類:型式平均值 S/N比說明1平均值及變異度均有影響2-對平均值有影響,變異數無影響3-變異數有影響,平均數無影響4-均無影

64、響170望目特性的策略當有特定的目標時,田口導出一獨特的策,在實驗進行中有兩種型式的控制因素可運用。影響變異的因素。影響平均值但不影響變異的因素(又稱調整因素或信號因素)。用S/N比找出以上兩種型態的因素後。以影響變異的因素來縮小變異至最小的程度以信號因素來調整平均值使接近目標。171LSLUSL參數設計的做法LSLUSLSTEP1STEP2172望目特性S/N以及靈敏度 y=望目特性值。控制因素A(A1,A2,A3,A4.,An)誤差因素N(N1,N2,N3Nn)173174175練習N1N2N3TOTAVGVeSmSmdbdbA19101130101.030024.7719.99A2901

65、00110300100100.0 3000044.7719.99A39910010130010013000044.7740.00A4900100011003000100010000300000064.7719.99A59901000101030001000100300000064.7740.00A699910001001300010001300000064.7760.00176望目特性的S/N值特性平均值增加S/N增加(當變異一樣時)變異減少S/N增加(當平均值一樣時)數據變化很大時要用Smdb來計算。Sm的變化等於平均值的變化,最好用db做單位,因為算術可加性較好。177變異的衡量(S/N及

66、Ve)S/N():衡量對平均值的變異度,以正負平均值的%表示。當特性值是絕對值時,(如y0)對望目特性值,我們可計算平均值的變異度,以下式計算:Ve:衡量變異度以正負絕對值表示之。當特性可以取負值(如y0),則正負平均值%即失去意義,所以我們必須轉換成正負絕對值表示:178靈敏度分析與傳統分析Smdb對平均值的變化比較靈敏,當平均值變動很大時,如在研究開發工作上常遇到這種情形,我們應當使用Smdb。當平均值變化不大時(後續實驗),我們則使用平均值。179望目特性例題金屬印模A:滾筒高度二水准B:原料供應商二水准C:進料調整二水准J:XX1二水准I:XX2二水准H:XX3二水准控制因素控制因素誤

67、差因素誤差因素回應特性:孔邊距離回應特性:孔邊距離目標:吋目標:吋(數據已簡化,則放大數據已簡化,則放大100,0.4100=40180配置望目特性實驗J1221I1212H 1122平均Ve(S/N)dbA BC1 1113738363737.000.6733.122 1223539403336.4510.92 20.923 2124544444644.750.9233.394 2214152464245.2524.92 19.15181練習找出回應表與回應圖找出影響變異的有效因素H1H2A137.336.5A245.544.5A2A1H1H2182J1J2B141.340.5B237.84

68、4.3454035B1B2J1J2AH, BH, CHAI, BI, CIAJ, BJ, CJ依此例再計算出來。依此例再計算出來。183NoABCS/N111133.1212220.9321233.4422119.2水准ABC127.033.326.2226.320.027.1A1A22024283236B1B22024283236C1C22024283236S/N比分析184No A B C 平均值11 1 1 37.0021 2 2 36.7532 1 2 44.7542 2 1 45.25水准ABC136.8840.8841.13245.0041.0040.75A1A2202428323

69、6B1B22024283236C1C22024283236平均值分析平均值:平均值:185練習銅線鍍錫回應特性:電鍍厚度規格值:10.01.5(原數據10000倍)設計因素水准1水准2水准3A噴嘴形狀現行新無B:乾燥溫度低中高C:線速低中高D:鎔接劑低中高E:金屬溫度低中高F:冷卻溫度低中高G:線型IIIIIIH:線張力低中高M線圈2水准N線圈位置2水准噪音因子噪音因子186N0ABCDEFGHM1N1M1N2M2N1M2N1平均S/N1111111116.08.07.69.87.8513.9862112222222.03.05.05.03.757.7823113333339.810.611.

70、010.510.5026.38441211223312.411.610.413.612.0018.9855122233116.010.09.812.49.5511.0676123311226.05.05.87.25.006.92871312132311.68.48.210.09.5515.5628132321315.65.86.46.06.0024.81791331321222.014.818.823.019.6514.473102113322115.217.613.220.616.5514.315112121133210.47.410.811.810.1014.501122132211311

71、.015.812.419.614.7011.58913221231321.07.03.28.44.902.606142223121310.419.215.818.215.9012.063152231232110.014.011.213.012.0516.53116231323123.04.03.64.23.7016.868172321312310.212.411.814.012.1017.717182332123115.416.613.214.815.0020.503187ABCDEFGH115.55414.76013.72116.03213.92412.94013.34116.87214.0

72、7711.26314.65811.51816.09614.68713.139318.32316.06816.89614.42716.81917.96617.05S/N回應表回應表ABCDEFGH19.31110.5929.10812.29210.5678.41711.89211.167211.6789.9009.5589.5758.69213.8259.8503-10.99212.8179.61712.2259.2429.742望目特性原始數據回應表望目特性原始數據回應表188S/N比的計算1891012141618A1A2B1B2 B3C1C2 C3D1D2 D31012141618E1E2

73、F1F2 F3G1G2 G3H1H2 H3E319089101112A1A2B1B2 B3C1C2 C3D1D2 D31389101112E1E2F1F2 F3G1G2 G3H1H2 H313E3191練習鍍金厚度回應特性:回應特性:y=鍍金厚度鍍金厚度(望目特性望目特性)控制因素一水准二水准三水准A:金濃度0.700.751.101.15B:酸度4.204.304.40C:溫度95105115D:滾桶速度101520E:陽極大小1/41/21/1F:重量大小1/41/21/1G:液流容度1.01.52.0H:鎳濃度600650700192誤差因素誤差因素N1複合不同之誤差因素使其縮小N2複合

74、各種誤差因素使其趨大假如回應特性y值之目標值50m,y為望目特性,控制因素用L18直交表配置放內側。誤差因素複合為一因素2水准配置於外側,每次測定二個樣本,得數據如下:計算Ve,Sm,Smdb, db並填入下表。計算各因素各水准之S/N比值,完成S/N比及平均值或Smdb之回應表。繪出回應圖並討論如何縮小變異並調整平均值使接近目標值。193N0ABCDEFGHN1N2平均VeSmSmdbdb1111111118388909121122222273333381311333333575865694121122335559616751222331173757679612331122586068727

75、13121323444955588132321315054576491331321264656668102113322174798694112121133275789095122132211370768588132212313271808795142223121348565965152231232166677986162313231245535864172321312360676673182332123157657983194ABCDEFGH123B1B2B3A1A2S/N比回應表比回應表195ABCDEFGH123B1B2B3A1A2平均值回應表平均值回應表196練習某脫腊鑄造之腊膜需沾漿,

76、沾漿厚度之目標值為12mm,設下列控制因子會對沾漿厚度造成影響。A 沾漿層數B漿汁稀釋比C沾漿時間D沾漿方式E漿料廠牌F漿桶轉速G添加物多寡其中A因子為二水准,其餘均為3水准,另有二個噪音因子腊模大小腊膜清洁度皆為二水准,試定義該問題之品質特性及設計該實驗。197123456781111111111010982111222221198103113333331010884121122331314121551222331111119136123311221012121371312132314151416813232131181915169133132121616141310211332219101

77、2911212113328111381221322113121113161322123132131415171422231213121210815223123216887162313231289871723211123788718233232311313131512121122IH198參數設計的步驟-1設立目的討論:要測什麼,如何使用資料具有可加性避免(0,1)資料分類值可以變換為連續變數分類組數越多越好分類值的分析也可能發現安定性的條件。S/N比是最好的特性值(可加性的機會加大)列出因素怕少,不怕多分類為控制因素、誤差因素199參數設計的步驟-2選擇直交表控制因素配置於內側直交表,誤差因素

78、配置於外側直交表。內側直交的選擇L12,L16,L18,L27,L32較實用。推廌L12,L18模擬時宜用L32填滿各行控制因素水准間隔要大。外側直交表規模要小,但誤差引起的變化要大。實行誤差因素的複合。重要的誤差因素有兩、三個即可。有時可以不配置誤差因素200參數設計的步驟-3做模擬試驗計算S/N比(公式如次頁所附)對應分析望小:僅分析望大:僅分析望目:y0時分析db , Smdby0時分析db,原始資料最佳化: 優先,然後調節平均到目標值。不顯著因素以低成本高生產力為優先。201參數設計的步驟-4最佳條件的估計證實試驗與估計比較損失函數與目前的比較求增益第二次實驗公差設計202工廠實際範例一工廠實際範例一203工廠實際範例二工廠實際範例二204

展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 医学/心理学 > 基础医学

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号