Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 2005-0461(Print)
ISSN : 2287-7975(Online)
Journal of Society of Korea Industrial and Systems Engineering Vol.41 No.2 pp.174-183
DOI : https://doi.org/10.11627/jkise.2018.41.2.174

Multivariate Process Capability Index Using Inverted Normal Loss Function

Hye-Jin Moon, Young-Bae Chung†
Department of Industrial and Management Engineering, Incheon National University
Corresponding Author : ybchung@inu.ac.kr
16/05/2018 27/06/2018 28/06/2018

Abstract


In the industrial fields, the process capability index has been using to evaluate the variation of quality in the process. The traditional process capability indices such as Cp, Cpk, CpM and CpM+ have been applied in the industrial fields. These traditional process capability indices are mainly applied in the univariate analysis. However, the main streams in the recent industry are the multivariate manufacturing process and the multiple quality characteristics are corrected each other. Therefore, the multivariate statistical method should be used in the process capability analysis. The multivariate process indices need to be enhanced with more useful information and extensive application in the recent industrial fields. Hence, the purpose of the study is to develop a more effective multivariate process index (MCpI) using the multivariate inverted normal loss function. The multivariate inverted normal loss function has the flexibility for the any type of the symmetrical and asymmetrical loss functions as well as the economic information. Especially, the proposed modeling method for the multivariate inverted normal loss function (MINLF) and the expected loss from MINLF in this paper can be applied to the any type of the symmetrical and asymmetrical loss functions. And this modeling method can be easily expanded from a bivariate case to a multivariate case.



역정규 손실함수를 이용한 다변량 공정능력지수

문 혜진, 정 영배†
인천대학교 산업경영공학과

초록


    1 서 론

    최근 기업들의 비즈니스 환경은 점점 까다롭고 다양 해져 가는 고객의 요구사항을 만족시켜야만 치열한 시장 경쟁에서 살아남을 수 있게 되었다. 그로 인하여 기업들 은 제품의 다양성, 고품질의 규격, 경쟁력 있는 가격 등 여러 생산조건들을 동시에 만족시키고자 기업의 역량을 강화하고 있다. 그리고 이러한 기업의 환경변화들로 인 하여 생산공정과 품질관리의 변수들은 복잡하면서도 다 양해지고 독립적이기 보다는 서로 연관성을 가지고 있다 고 할 수 있다. 그러기에 지금의 생산 공정과 품질요소들 을 관리하기 위해서는 합리적이고도 현실성 있는 유용한 통계적 척도가 요구될 것이다.

    생산현장에서 공정능력(Process Capability)이 가지는 의 미는 관리상태하에 있을 때 그 공정에서 생산되는 제품의 품질변동이 어느 정도인가를 나타내는 것으로써, 공정의 균일성(Uniformity)을 의미한다. 그리고 공정능력평가를 위해 대표적으로 사용되는 척도가 공정능력지수(Process Capability Index)이다. 공정능력지수는 공정의 변동 과 제품 의 규격한계를 정량적으로 비교하여 공정능력을 평가하 게 된다. 전통적인 공정능력지수 Cp, Cpk 는 품질의 산포 및 공정평균의 치우침만을 고려하여 평가한 반면 목표치 로부터 품질 특성치의 변동까지도 고려하여 공정능력을 평가하는 CpM이 제안되었다[2]. C p m + 는 목표치로부터 품 질 특성치의 변동을 Taguchi의 이차손실 함수의 기대값을 적용하여 평가하는 척도로써 제안되었다[1]. 그리고 더 나아가 목표치로부터 품질의 변동 에 따른 경제적 손실까 지도 공정능력 평가에 포함함으로써 보다 유용한 정보력 과 폭넓은 응용력을 지닌 다양한 공정능력지수들이 연구 되어 왔다. <Table 1>

    특히 목표치로부터 품질의 변동에 따른 경제적 손실을 잘 표현해 주는 대표적인 수리모형이 손실 함수(Loss Function) 이다. 그리고 Taguchi의 이차 손실함수는 오랫동안 다양한 산업분야에서 사용 되어 왔다. 하지만 Taguchi의 이차손실함수는 다음 과 같은 한계점을 가지고 있다. 품 질 특성치가 규격 내에 있을 때나 규격을 벗어났을 때나 일률 적으로 이차함수를 적용하고 있다는 점과 목표치 를 중심으로 좌우대칭일 경우에만 잘 설명 된다는 점이다. <Table 2>

    이러한 한계점을 극복하고자 Spiring은 대안적 손실함 수를 제안하였다[13]. Spiring의 손실함수는 목표치로터 손실에 대한 확률밀도함수의 역함수를 이용하였다. 이로 인하여 생산공정의 현실성을 반영하여 다양한 분포의 형 태와 품질특성마다 발생할 수 있는 서로 다른 경제적 손 실을 비교적 제약없이 사용할 수 있다는 장점을 가지고 있다.

    그리고 최근 생산공정에서 관리해야 할 요인들은 2개 이상으로 다변량 분석(Multivariate Analysis)이 요구된다 는 점이다. 또한 생산 및 측정 기술력의 발전과 공정의 자 동화 등으로 인하여 생산공정은 더욱더 복잡해지고 수많 은 품질요인들이 서로 연관되어 품질 변동에 중요한 영향 을 미치고 있다. 그러기에 많은 연구자들은 1개의 품질 특 성치만을 고려한 기존의 전통적 공정능력 지수들로 최신 의 생산공정을 정확하게 평가하기에는 한계점이 있다고 판단하였고, 2개 이상의 품질변수에 대한 규격과 산포를 고려한 다변량 공정능력지수(Multivariate Capability Index) 들이 연구되고 있다.<Table 3>

    따라서 본 연구에서는 최신의 산업현장에서 현실적이 면서도 보다 유용한 정보들을 가지고 있고 공정상태를 보다 민감하게 감지하고 관리할 수 있는 다변량 공정능 력지수를 제안하고자 한다. 경제적 측면의 정보력을 포 함하기 위해서 손실 함수를 이용할 것이다. 그리고 손실 함수의 모형은 공정의 다양성을 잘 반영하고 응용력이 뛰어난 Spiring의 대안적 손실함수를 채택하고자 한다. 특히, 목표치로부터 손실이 정규분포를 따를 때 정규분 포의 역함수, 즉 역정규분포의 손실함수(Inverted Normal Loss Function)를 이용하여 경제적 손실을 공정능력 평가 에 반영할 것이다. 그러면 서도 동시에 여러 변수들을 같 이 분석할 수 있는 다변량 분석기법을 이용할 것이다. 또 한 목표치를 중심으로 최대손실에 대한 분포의 유형이 좌우 대칭뿐만 아니라 비대칭일 때도 제약없이 평가 가 능한 유용성 있는 다변량 공정능력지수를 제안 하고자 한다.<Table 4>

    2 이론적 배경

    2.1 다변량 공정능력지수 MCpm

    다변량 공정능력지수는 전통적인 단변량 공정능력지수 를 벡터화된 모형으로 확장한 개념이라 할 수 있다. 즉, 정규분포일 경우를 기준으로 단변량과 다변량의 특징을 비교해 보면, 분포의 모양은 단변량일 때 종모양(ballshaped) 인 반면 다변량일 때는 타원형(2차원) 또는 타원체 (3차원)을 띤다. 확률적 표현에 있어서 단변량일 때는 간 격/폭 에 의한 길이지만 다변량일 때는 영역/범위로 나타 내는 면적 또는 부피로 표현된다.

    대표적 다변량 공정능력지수로 CpMMCpm을 비교 하고자 한다[11, 14]. <Figure 1>은 2변량 정규 분포일 때의 규격허용공차(USLi - LSLi)와 공정영역(UPLi - LPLi)을 보여주고 있다.

    다변량 공정능력지수는 규격의 허용오차영역과 공정 영역의 면적 또는 부피의 비율로써 공정능력을 측정하게 되며, CpMMCpm의 식은 각각 다음과 같다.(1)

    C p M = [ 허용 공차영역의 면적 또는 부피 수정 공정영역의 면적 또는 부피 ] 1 n = [ i = 1 n ( U S L i L S L i ) i = 1 n ( U P L i L P L i ) ] 1 n
    (1)

    M C p m = 수정 규격허용오차영역 ( R 1 ) 의 면적 또는 부피 공정영역 ( R 2 ) 의 면적 또는 부피 = 수정 규격허용오차영역 ( R 1 ) 의 면적 또는 부피 [ ( X μ ) T T 1 ( X μ ) K ( n ) ] 의 면적 또는 부피 = 수정 규격허용오차영역 ( R 1 ) 의 면적 또는 부피 | Σ T | 1 2 ( π K ) n 2 [ Γ ( n 2 + 1 ) ] 1 = 수정 규격허용오차영역 ( R 1 ) 의 면적 또는 부피 | Σ T | 1 2 ( π K ) n 2 [ Γ ( n 2 + 1 ) ] 1 = | Σ T | 1 2 ( π K ) n 2 [ Γ ( n 2 + 1 ) ] 1 [ 1 + ( μ T ) T Σ 1 ( μ T ) 1 2 ]
    (2)

    여기서 n은 품질 특성치의 수이고, ΣT는 평균제곱오차 매트릭스, Σ는 분산-공분산 매트릭스, Γ(α)는 감마분포 를 의미한다. K(n)은 변수가 n개일 때 χ2분포의 99.73 번째 분위수를 의미하며, 이는 “6σ”의 개념으로 χ2(n, 0.9973)의 값을 말한다. 이때 MCpm= 1이면 공정의 평균 과 목표치가 일치함을 의미하며 1보다 크거나 같으면 공 정능력이 양호하다고 평가할 수 있다. <Table 5>

    상기 <Figure 1>과 식 (1)과 식 (2)를 비교해 보면, CpM은 타원형의 공정영역을 사각형의 영역으로 수정하여 그 면적을 적용하였다면 MCpm은 사각 형의 규격허용오차영 역을 실제의 공정영역에 맞추어 타원형으로 수정하여 적 용하였다는 차이점 을 가지고 있다. 이러한 차이점은 결과 적으로 MCpmCpM보다 목표치로부터 공정평균의 편차 를 보다 민감하게 반영하고 있음을 증명 하였다. 따라서 본 연구에서는 MCpm의 다변량 모형을 응용하고자 한다.

    2.2 다변량 역정규 손실함수(MINLF)

    다양한 산업분야에서 널리 사용되는 손실함수는 Taguchi 의 이차손실함수이다. 이는 품질 특성치 x가 목표치 T와 일치하면 손실이 발생하지 않지만 목표 치로부터 편 차가 크면 클수록 손실이 증가 한다는 가정하에 이를 이 차식으로 근사화한 손실함수이다. <Table 6>

    반면 Spiring은 Taguchi의 이차손실함수의 대안 적 손실함 수로써 정규분포의 확률밀도함수(p.d.f)의 역함수를 활용한 손실함수를 제안하였고, 이를 역정규 손실함수(INLF : Inverted Normal Loss Function)라고 정의하였다[9, 13].

    Taguchi의 이차손실함수에서 규격을 벗어나도 손실이 계속 증가한다는 한계점을 보강한 이차 손실함수 <Figure 2>와 Spiring의 역정규 손실함수 <Figure 3>의 그래프를 비교해 보면 다음과 같다[11]. 이때 (T - Δ,T+ Δ)는 소비 자의 허용차 구간이고, 이 구간을 벗어나면 최대손실비용 A가 발생하게 된다.<Table 7>

    두 그래프를 비교해 보면, 역정규 손실함수 역시 Taguchi 의 이차손실함수의 특징인 손실이 항상 양수의 값을 가 지며 (-∞,T]구간에서는 감소하고 [T,∞)구간에서는 증 가하는 성향을 나타내고 있다. 또한 목표치 T에서 유일 한 최소값을 가진다는 손실함수의 대표적 특징을 동일하 게 나타내고 있다. 그리고 정규분포의 확률밀도함수의 역함수를 적용함으로써 품질변동 에 대한 손실을 보다 합리적으로 설명이 가능하며 다음의 <Figure 4>와 같이 목표치(T)를 중심으로 최대손실(A)에 대하여 대칭과 비 대칭 모두일 때 적용이 가능하다는 장점을 가지고 있다.<Table 8>

    실제 공정 데이터의 분포에 대한 확률밀도 함수를 f(x)라 하고, 손실함수에 이용되는 목표치에 대한 확률 밀도함수를 g(x,T)로 나타낸다. 따라서 역함수에 의한 손실함수 L(x)는 함수 g(x,T)의 역함수를 취한 식(3)과 같이 나타낼 수 있으며 그의 기대값 EL(x)]은 다음의 식(4)와 같다.

    L ( x ) = A [ 1 g ( x , T ) m ]
    (3)

    E [ L ( x ) ] = A [ 1 g ( x , T ) m ] f ( x ) d x
    (4)

    여기서 m은 함수 g (x, t)가 최대로 도달할 수 있는 값 으로 품질 특성치 x가 목표치 T 와 일치할 때 최대값을 가지게 된다. 따라서 x = T 때의 g (x, t) 값이 m의 값이 된다.

    그리고 상기 식 (3)과 식 (4)는 실제 공정 데이터의 분 포가 정규분포, 감마분포, 베타분포 등 다양성에 따라 해 당 분포에 대한 확률밀도함수 f(x)와 손실함수에 이용 되는 목표치에 대한 확률밀도함수를 g(x,T)를 적용하여 활용할 수도 있다[9, 12].

    본 연구에서는 공정데이터가 정규분포를 따른다는 가 정하에서 다음과 같이 수식을 전개하고자 한다.

    f ( x ) = 1 2 π σ e x p { ( x μ ) 2 2 σ 2 } g ( x , T ) = 1 2 π γ e x p { ( x T ) 2 2 γ 2 }

    여기서 상기 손실함수 L(x)의 식 (3)에 g(x,T)식과 m = 1 2 π γ 을 대입하면 식 (5)와 같다.

    L ( x ) = A [ 1 g ( x , T ) m ] = A [ 1 exp { ( x T ) 2 2 γ 2 } ]
    (5)

    그리고 기대손실 E[L(x)]은 다음과 같다[13].(6)

    E [ L ( x ) ] = A L ( x ) f ( x ) d x = A [ 1 γ σ 2 + γ 2 e x p { ( μ T ) 2 2 ( σ 2 + γ 2 ) } ]
    (6)

    목표치 T를 중심으로 (-∞,T]구간과 [T,∞)구간에서 비대칭 손실함수 L(x)은 다음과 같이 나타낼 수 있다.(7)

    L ( x ) = { A 1 1 exp { ( x T ) 2 2 γ 2 } ] , x < T A 2 [ 1 exp { ( x T ) 2 2 γ 2 } ] , x T
    (7)

    이때 기대손실 E[L(x)]은 다음과 같다.(8)

    E [ L ( x ) ] = T A 1 L ( x ) f ( x ) d x + T A 2 L ( x ) f ( x ) d x = T A 1 [ 1 exp { ( x T ) 2 2 γ 2 } ] f ( x ) d x + T A 2 [ 1 exp { ( x T ) 2 2 γ 2 } ] f ( x ) d x
    (8)

    단변량 모형을 벡터로 확장하면 다변량 역정규 손실함 수(MINLF : Multivariate Inverted Normal Loss Fuction)를 구할 수 있으며, 손실함수 L(X,T)와 기대손실 E[L(X,T)]은 다음과 같다.(9)(10)

    L ( X, T ) = A [ 1 g ( X, T ) m ]
    (9)

    E [ L ( X, T ) ] = A [ 1 g ( X, T ) m ] f ( X ) d x
    (10)

    본 연구범위에 따른 다변량 정규분포의 확률 밀도 함 수 f(X)와 손실함수에 이용되는 목표치에 대한 확률밀 도함수 g(X,T)는 다음과 같다.

    f ( X ) = 1 ( 2 π ) n 2 | Σ | 1 2 e x p { 1 2 ( X μ ) T Σ T 1 ( X μ ) } g ( X, T ) = 1 ( 2 π ) n 2 | Λ | 1 2 e x p { 1 2 ( X T ) T Λ T 1 ( X T ) }

    상기 f(X)와 g(X,T)를 대입하면 역정규 손실 함수 L(X,T)와 기대손실 E[L(X,T)]은 다음과 같다[6].

    L ( X, T ) = A [ 1 g ( X, T ) m ] A [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ]
    (11)

    E [ L ( X, T ) ] = A L ( X, T ) f ( X ) d X = A [ 1 | Λ 1 + Σ 1 | 1 / 2 | Σ | 1 / 2 × exp { 1 2 ( μ T Σ 1 μ + T T Λ 1 T ) } ] × exp { 1 2 ( μ T Σ 1 + T T Λ 1 ) ( Λ 1 + Σ 1 ) 1 ( Σ 1 μ + Λ 1 T ) }
    (12)

    3 연구모형

    3.1 기호

    • X : 표본 품질 특성치 벡터(X′ = xi = x1,x2,⋯,xn)

    • μ: 공정 평균벡터(μ′ = μi = μ1,μ2,⋯,μn)

    • T: 목표지 벡터(T′= Ti= T1,T2,⋯,Tn)

    • n: 품질 특성치의 수

    • i : 품질 특성치(i = 1,2,⋯,n)

    • Aj: 목표치 Ti중심으로 특성치 xi의 비대칭 구간에서 최대 손실(j = 1,2,⋯,2n)

    • Σ: 분산-공분산 매트릭스

    • Σ-1 : 분산-공분산 매트릭스 역행렬

    • ΣT: 평균제곱오차(MES) 매트릭스

    • D : 목표치로부터 공정평균의 편차

    • Λ : 척도모수(scaling parameter) 매트릭스

    • K(n): χ2(n,0.9973)의 값

    • Γ(α): 감마함수

    • L(X,T): 목표치 T에 대한 표본벡터 Xi의 손실함수

    • E[L(X,T)]: 목표치 T에 대한 표본벡터 Xi의 기대손실

    3.2 다변량 공정능력 지수 MCpI모형

    본 연구모형에 근간으로 적용하기 위해 채택한 다변량 공정능력지수 MCpm은 규격범위와 목표치로부터 공정의 산포를 비교한 다변량 공정능력 지수로써 단변량의 CpM과 동일한 개념을 가지고 있으면서 벡터로 확장한 다변 량 공정능력지수이다.

    MCpm모형에 경제적 측면을 반영할 수 있는 품질손실 함수를 적용하기 위한 연구모형의 논리적 근거는 단변량 공정능력지수 CpM으로부터 이차 손실함수의 기대값을 적 용한 C p m + 를 도출한 것과 같은 방식으로 전개하고자 한다. 즉, C p m + 는 단변량 CpM에서 목표로부터 공정 산포(σ) 대 신 Taguchi의 이차손실함수의 기대값을 응용한 공정능력 지수라 할 수 있다. 따라서 MCpm모형에 목표로부터 공 정의 산포값 대신 다변량의 역정규 손실함수를 적용할 것 이다. 역정규 손실함수를 이용함으로써 다수의 품질 특성 치마다 서로 다른 손실비용을 반영하고 그로 인한 손실함 수의 형태가 대칭 또는 비대칭 구조에 의한 분석의 제약 성을 받지 않는 모형을 연구하는데 초점을 두었다.

    따라서 본 연구는 첫째로 다수의 품질 특성치를 고려 한 다변량 공정능력분석을 목적으로 한다. 둘째는 공정 에서의 경제적 손실척도로써 역정규 손실함수를 적용한 다변량 공정능력분석에 있다. 실제 생산현장에서 유용한 정보력을 가지면서 손실함수의 형태에 따른 제약성 없이 용이한 다변량 공정능력 분석을 목적으로 한다.

    상기 <Figure 5>는 본 연구모형을 위한 모델링 방법론 을 나타내고 있다. 기존의 이론적 배경을 바탕으로 하여 다변량의 MCpm모형으로부터 목표치를 고려한 공정의 산포 대신에 다변량 정규분포의 역함수를 이용한 손실함 수의 기대값을 적용한 공정능력지수 MCpI모형을 제시 하고자 한다.

    우선 제안하고자 하는 다변량 공정능력지수 모형의 기 초모델이 되는 MCpm의 식 (2)에서 본 연구의 초점에 되 는 부분은 다음의 식과 같다.

    M C p m = 수정 규격허용오차영역의 면적 또는 부피 | Σ T | 1 2 ( π K ) n 2 [ Γ ( n 2 + 1 ) ] 1
    (13)

    상기 식 (13)에서 평균제곱오차매트릭스 T 는 아래의 구조식에서 설명해 주듯이 공정평균이 목표치로부터 얼마 나 떨어져 있는가에 대한 편차를 의미하며, 이는 손실함수 가 가지는 특징으로도 동일하게 설명이 되는 부분이다.

    Σ T = E [ ( X - T ) ( X - T ) ] = 1 n i = 1 n ( X i  - T ) ( X i  - T ) = Σ + ( μ T ) ( μ T )

    이러한 특징을 기반으로 하여 앞에서 설명한 것과 같이 CpM의 산포 대신 C p m + 는 이차손실 함수의 기대값 E[L(x)]을 적용했듯이 MCpm의 평균제곱오차 매트릭스 ΣT대신에 다변량 역정규 손실함수(MINLF)의 기대손실 E[L(X,T)]을 대입 하는 방법으로써 다변량 공정능력지수 MCpI모형 을 다음과 같이 제안한다.

    M C p I = 수정 규격허용오차영역의 면적 또는 부피 | E [ L ( X , T ) ] | 1 2 ( π K ) n 2 [ Γ ( n 2 + 1 ) ] 1
    (14)

    상기 식 (14)의 연구모형에서 중요한 입력변수는 다변 량의 기대손실 E[L(X,T)]이다.

    다변량인 MCpI의 구조식에 맞추어 분모의 기대 손실 E[L(X,T)]역시 다변량의 모형으로 확장 하여 정의하는 것이 필요하다. 또한 손실함수의 형태, 즉 대칭과 비대칭 구조를 모두 설명할 수 있는 일반화된 다변량 역정규 손 실함수(MINLF)의 수리적 모델을 정의하고자 한다.

    최대손실 A에 대하여 대칭적 손실함수라면 앞에서의 식 (11)과 식 (12)에 대입하여 구할 수 있다. 하지만 비대칭의 경우에는 우선적으로 다변량의 특성치 xi마다 목표치 Ti를 중심으로 비대칭을 이루는 구간의 전개와 해당 구간에 서의 최대손실 비용 Aj에 대하여 정의가 필요하다. 구간을 나누는 기준은 손실함수의 특징인 목표치 Ti일 때 손실값 0(zero)을 중심으로 규격 상한선(USL)과 하한선(LSL) 범위 내에서 품질 특성치 xi들 간에 연계된 구간들을 구분하여 전개하면 <Figure 6>과 같다. 그리고 각 구간에서 발생하는 최대손실비용 Aj의 값은 생산현장으로부터 규정된 손실비 용을 기반으로 한 정해진 상수값으로 가정한다. <Figure 6> 의 구간은 역정규 손실함수의 기대값 E[L(X,T)]을 구하는 정적분의 구간이 된다. 또한 <Figure 6> 방법에 의해 구한 손실함수와 기대손실은 대칭과 비대칭 모두 사용할 수 있 는 일반화된 구조식을 도출할 수 있기 때문에 응용범위가 넓어질 수 있고 사용이 용이하다는 정점을 가지고 있다.

    따라서 <Figure 6>의 구간전개방식에 따라 대칭〮 비대 칭 구조를 망라한 역정규 손실함수 L(X,T)와 기대손실 E[L(X,T)]은 다음과 같다.(15)(16)

    L ( X , T ) = { A 1 [ 1 exp { 1 1 2 ( X T ) T Λ 1 ( X T ) } ] , x 1 < T 1 , x 2 < T 2 , , x n < T n A 2 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , x 1 < T 1 , x 2 < T 2 , , x n T n A 2 n 1 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , x 1 T 1 , x 2 T 2 , , x n < T n A 2 n [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , x 1 T 1 , x 2 T 2 , , x n T n
    (15)

    E [ L ( X , T ) ] = T 1 T n A 1 L ( X , T ) f ( X ) d X + T 1 T n A 2 L ( X , T ) f ( X ) d X T 1 T n A 2 n 1 L ( X , T ) f ( X ) d X + T 1 T n A 2 n L ( X , T ) f ( X ) d X
    (16)

    본 연구의 시뮬레이션을 위한 범위와 가정은 다음과 같다.

    • (1) 품질 특성치는 망목특성이며, 연속적인 값을 가진다.

    • (2) 공정데이터는 다변량 정규분포(Multivariable Normal Distribution)을 따른다.

    • (3) 대칭일 때와 비대칭일 때 모두를 고려한다.

    • (4) 목표치로부터 공정평균에 치우침이 있는 경우와 치 우침이 없는 경우 모두 고려한다.

    • (5) 척도 매트릭스(Scaling matrix)에서 λ는 공정 파라미 터가 규격한계선에 있을 때 손실이 90%일 때를 기준 으로 한다(λi = 0.233 × 규격공차i).

    • (6) 구간별 최대손실비용 Aj은 기업으로부터 정해진 상수 값이다.

    • (7) 시뮬레이션 분석에 사용한 프로그램은 MATLAB R2015b 이다.

    본 연구모형의 시뮬레이션은 2변량(Bivariable)과 3변량 (Three-variable)일 때의 모형으로 분석 하고자 한다.

    1) 2변량 모형(Bivariable)

    X = [ x 1 x 2 ] , μ = [ μ 1 μ 2 ] , μ = [ T 1 T 2 ] , USL = [ U S L 1 U S L 2 ] , LSL = [ L S L 1 L S L 2 ] , Λ = [ λ 1 2 λ 1 λ 2 λ 2 λ 1 λ 2 2 ] , Σ = [ σ 1 2 ρ σ 1 σ 2 ρ σ 2 σ 1 σ 2 2 ] A = [ A 1 , A 2 , A 3 , A 4 ]

    2변량의 역정규 손실함수와 기대손실은 다음과 같다.

    L ( X , T ) = { A 1 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , L S L 1 < x 1 < U S L 1 , L S L 2 < x 2 < T 2 A 2 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , T 1 < x 1 < U S L 1 , L S L 2 < x 2 < T 2 A 3 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , L S L 1 < x 1 < T 1 , T 2 < x 2 < U S L 2 A 4 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , L S L 1 < x 1 < U S L 1 , L S L 2 < x 2 < U S L 2

    E [ L ( X , T ) ] = L S L 1 T 1 L S L 2 T 2 A 1 L ( X , T ) f ( X ) d X + L S L 1 T 1 L S L 2 T 2 A 2 L ( X , T ) f ( X ) d X + T 1 U S L 1 L S L 2 T 2 A 3 L ( X , T ) f ( X ) d X + T 1 U S L 1 T 2 U S L 2 A 4 L ( X , T ) f ( X ) d X = L S L 2 T 1 L S L 2 T 2 A 1 L ( x 1 , x 2 ) f ( x 1 , x 2 ) d x 2 d x 1 + L S L 1 T 1 T 2 U S L 2 A 2 L ( x 1 , x 2 ) f ( x 1 , x 2 ) d x 2 d x 1 + T 1 U S L 1 L S L 2 T 2 A 2 L ( x 1 , x 2 ) f ( x 1 , x 2 ) d x 2 d x 1 + T 1 U S L 1 T 2 U S L 2 A 4 L ( x 1 , x 2 ) f ( x 1 , x 2 ) d x 2 d x 1

    2) 3변량 모형(Three-Variable)

    X = [ x 1 x 2 x 3 ] , μ = [ μ 1 μ 2 μ 3 ] , μ = [ T 1 T 2 T 3 ] , USL = [ U S L 1 U S L 2 U S L 3 ] , LSL = [ L S L 1 L S L 2 L S L 3 ] , Λ = [ λ 1 2 λ 1 λ 2 λ 1 λ 3 λ 2 λ 1 λ 2 2 λ 2 λ 3 λ 3 λ 1 λ 3 λ 2 λ 3 2 ] , Σ = [ σ 1 2 ρ σ 1 σ 2 ρ σ 1 σ 3 ρ σ 2 σ 1 σ 2 2 ρ σ 2 σ 3 ρ σ 3 σ 1 ρ σ 3 σ 2 σ 3 2 ] A = [ A 1 , A 2 , A 3 , A 4 , A 5 , A 6 , A 7 , A 8 ]

    3변량의 역정규 손실함수와 기대손실은 다음과 같다.

    L ( X , T ) A 1 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , = { L S L 1 < x 1 < T 1 , L S L 2 < x 2 < T 2 , L S L 3 < x 3 < T 3 A 2 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , L S L 1 < x 1 < T 1 , L S L 2 < x 2 < T 2 , L S L 3 < x 3 < U S L 3 A 3 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , L S L 1 < x 1 < T 1 , T 2 < x 2 < U S L 2 , L S L 3 < x 3 < T 3 A 4 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , L S L 1 < x 1 < T 1 , T 2 < x 2 < U S L 2 , T 3 < x 3 < U S L 3 A 5 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , T 1 < x 1 < U S L 1 , L S L 2 < x 2 < T 2 , L S L 3 < x 3 < T 3 A 6 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , T 1 < x 1 < U S L 1 , L S L 2 < x 2 < T 2 , T 3 < x 3 < U S L 3 A 7 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , T 1 < x 1 < U S L 1 , T 2 < x 2 < U S L 2 , L S L 3 < x 3 < T 3 A 8 [ 1 exp { 1 2 ( X T ) T Λ 1 ( X T ) } ] , T 1 < x 1 < U S L 1 , T 2 < x 2 < U S L 2 , T 3 < x 3 < U S L 3

    E [ L ( X , T ) ] = L S L 1 T 1 L S L 2 T 2 L S L 3 T 3 A 1 L ( X , T ) f ( X ) d X + L S L 1 T 1 L S L 2 T 2 T 3 U S L 3 A 2 L ( X , T ) f ( X ) d X + L S L 1 T 1 T 2 U S L 2 L S L 3 T 3 A 3 L ( X , T ) f ( X ) d X + L S L 1 T 1 T 2 U S L 2 T 3 U S L 3 A 4 L ( X , T ) f ( X ) d X + T 1 U S L 1 L S L 2 T 2 L S L 3 T 3 A 5 L ( X , T ) f ( X ) d X + T 1 U S L 1 L S L 2 T 2 T 3 U S L 3 A 6 L ( X , T ) f ( X ) d X + T 1 U S L 1 T 2 U S L 2 L S L 3 T 3 A 7 L ( X , T ) f ( X ) d X + T 1 U S L 1 T 2 U S L 2 T 3 U S L 3 A 8 L ( X , T ) f ( X ) d X = L S L 1 T 1 L S L 2 T 2 L S L 3 T 3 A 1 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + L S L 1 T 1 L S L 2 T 2 T 3 U S L 3 A 2 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + L S L 1 T 1 T 2 U S L 2 L S L 3 T 3 A 3 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + L S L 1 T 1 T 2 U S L 2 T 3 U S L 3 A 4 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + T 1 U S L 1 L S L 2 T 2 L S L 3 T 3 A 5 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + T 1 U S L 1 L S L 2 T 2 T 3 U S L 3 A 6 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + T 1 U S L 1 T 2 U S L 2 L S L 3 T 3 A 7 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1 + T 1 U S L 1 T 2 U S L 2 T 3 U S L 3 A 8 L ( x 1 , x 2 , x 3 ) f ( x 1 , x 2 , x 3 ) d x 3 d x 2 d x 1

    4 수치 예

    다음의 수치 예는 4가지의 공정조건에 대하여 품질 특 성치들 간의 상관계수(correlation coefficient)가 변할 때 단변량과 다변량의 모형을 비교하고자 한다.

    4.1 2변량(Bivariate)

    4.1.1 대칭일 때

    (1) 치우침이 없을 때(μ = T)

    T1 = 5.5, T2 = 5.5, x1~N(5.5, 0.8), x2~N(5.5, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, A′= [1, 1, 1, 1], λ1 1.398, λ2 = 1.398, є = 0.65

    (2) 치우침이 있을 때(μT)

    T1 = 5.5, T2= 5.5, x1~N(6.5, 0.8), x2~N(3.5, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, A′= [1, 1, 1, 1], λ1= 1.398, λ2= 1.398, є = 0.65

    4.1.2 비대칭일 때

    (1) 치우침이 없을 때(μ = T)

    T1 = 5.5, T2 = 5.5, x1~N(5.5, 0.8), x2~N(5.5, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, A′= [1, 2, 3, 4], λ1 = 1.398, λ2 = 1.398, є= 0.65

    (2) 치우침이 있을 때(μT)

    T1 = 5.5, T2 = 5.5, x1~N(6.5, 0.8), x2~N(3.5, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, A′= [1, 2, 3, 4], λ1 = 1.398, λ2 = 1.398, є= 0.65

    상기 2변량의 대칭과 비대칭 경우를 비교해 보면, 치우 침이 없을 때 MCpm은 상관계수가 1에 가까워지면 공정능 력을 과대평가하고 있는 반면 MCpI는 과대평가 없이 안정 된 공정능력의 평가 결과를 보여주고 있다. 또한 MCpm은 다양한 손실비용(A)을 반영하지 못하고 대칭일 때나 비대 칭일 때 동일한 공정능력 값을 나타내지만 MCpI는 다양한 손실비용(A)을 공정능력에 반영하고 있음을 확인할 수 있 다. 그리고 MCpmMCpI모두 치우침이 있을 때 보다는 없을 때 더 나은 공정능력 결과값을 보여주고 있다.

    4.2 3변량(Three-Variate)

    4.2.1 대칭일 때

    (1) 치우침이 없을 때(μ = T)

    T1 = 5.5, T2 = 10.0, T3 = 7.5, x1~N(5.5, 0.8), x2~N(10.0, 1.5), x3~N(7.5, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, LSL3 = 2.5, USL3 = 8.5, A′= [1, 1, 1, 1, 1, 1, 1, 1], λ1 = 1.398, λ2= 1.398, λ3= 1.398, є= 0.65

    (2) 치우침이 있을 때(μT)

    T1= 5.5, T2= 10.0, T3= 7.5, x1~N(4.5, 0.8), x2~N(8.0, 1.5), x3~N(9.0, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, LSL3 = 2.5, USL3 = 8.5, A′= [1, 1, 1, 1, 1, 1, 1, 1], λ1= 1.398, λ2= 1.398, λ3= 1.398, є= 0.65

    4.2.2 비대칭일 때

    (1) 치우침이 없을 때(μ= T)

    T1= 5.5, T2= 10.0, T3= 7.5, x1~N(5.5, 0.8), x2~N(10.0, 1.5), x3~N(7.5, 1.0), LSL1 = 2.5, USL = 8.5, LSL2 = 2.5, USL2 = 8.5, LSL3 = 2.5, USL3 = 8.5, A′= [1, 2, 3, 4, 5, 6, 7, 8], λ1= 1.398, λ2= 1.398, λ3= 1.398, є= 0.65

    (2) 치우침이 있을 때(μT)

    T1 = 5.5, T2= 10.0, T3= 7.5, x1~N(4.5, 0.8), x2~N(8.0, 1.5), x3~N(9.0, 1.0), LSL1 = 2.5, USL1 = 8.5, LSL2 = 2.5, USL2 = 8.5, LSL3 = 2.5, USL3 = 8.5, A′= [1, 2, 3, 4, 5, 6, 7, 8], λ1= 1.398, λ2= 1.398, λ3= 1.398, є= 0.65

    3변량의 시뮬레이션 결과 역시 2변량일 때와 동일한 특징들을 나타내고 있다. 특히 변수가 증가한 3변량 분 석에서 MCpm은 대칭과 비대칭의 치우침이 없을 때 상관 계수가 1에 가까워지면 급격하게 과대평가하고 있으나 MCpI는 안정된 평가결과값을 보여주고 있다.

    5 결 론

    본 연구에서 제시하는 역정규 손실함수를 이용한 다 변량 공정능력지수 MCpI를 기존의 단변량 및 다변량 공 정능력지수와 비교한 결과 다음과 같다.

    2변량 및 3변량 모두 MCpm은 상관계수가 1에 가까워 지면 과대평가되는 반면 MCpI는 상관계수의 변화에 따른 결과값들이 크게 과대평가되지 않고 있다. 그리고 각각의 단변량 C p m + 와 다수의 특성들을 동시에 분석한 다변량 MCpI를 비교해 보아도 상관관계에 대하여 로버스트한 결 과들을 보여주고 있다. 또한 현실적으로 품질 특성치마다 발생하는 손실이 일률적이지 않기 때문에 서로 다른 손실 비용을 적용할 필요가 있다. 이러한 측면에서 MCpm은 경 제적 손실 정도를 반영하지 못하는 반면 MCpI는 다변량 간에 발생하는 다양한 손실비용을 반영하여 공정능력을 평가하고 있다. 그리고 손실비용의 증감에 따른 민감성도 가지고 있다.

    MCpI는 공정의 관리상태나 공정의 조건변화에 따라 민 감한 공정능력 평가결과를 보여주고 있다. 상기 시뮬레이 션에서 4가지 경우의 공정조건들을 비교해 보면, 치우침 있을 때보다는 치우침이 없을 때가 더 좋은 성향을 나타 내고, 비대칭의 손실비용이 클 때 보다는 대칭의 손실비 용이 작을수록 더 좋은 결과값을 나타내 주고 있다. 또한 본 연구에서 제시한 방법론에 따라 역정규분포의 손실함 수와 기대손실을 적용하게 되면 대칭일 때나 비대칭일 때 모두 유연성 있게 사용할 수 있다는 것과 다변량으로의 확장성이 용이하다는 장점을 가지고 있다.

    Figure

    JKISE-41-174_F1.gif

    Tolerance and Process Regions for Bivariate Case

    JKISE-41-174_F2.gif

    Modified Quadratic Loss Function

    JKISE-41-174_F3.gif

    Symmetric Inverted Normal Loss Function

    JKISE-41-174_F4.gif

    Asymmetric Inverted Normal Loss Functionv

    JKISE-41-174_F5.gif

    Modeling Method of MCpI

    JKISE-41-174_F6.gif

    Modeling Method of Asymmetric Interval for MCpI

    Table

    Comparison between Univariate and Bivariate Cases in Symmetrical Situation(μ = T)

    Comparison between Univariate and Bivariate Cases in Symmetrical Situation(μ ≠ T

    Comparison between Univariate and Bivariate Cases in Asymmetrical Situation(μ = T)

    Comparison between Univariate and Bivariate Cases in Asymmetrical Situation (μ ≠ T)

    Comparison between Univariate and Three-variate Cases in Symmetrical Situation(μ = T)

    Comparison between Univariate and Three-variate Cases in Symmetrical Situationμ ≠ T

    Comparison between Univariate and Three-variate Cases in Asymmetrical Situation(μ = T)

    Comparison between Univariate and Three-variate Cases in Asymmetrical Situationμ ≠ T

    Reference

    1. R.A. Boyles (1991) The Taguchi Capability index., J. Qual. Technol., Vol.23 (1) ; pp.17-26
    2. L.K. Chan , S.W. Chung , F.A. Spiring (1988) A New Measure of Process Capability : Cpm., J. Qual. Technol., Vol.20 (3) ; pp.162-175
    3. Y.B. Chung , Y.S. Kim (2007) Comparison and Application of Process Capability indices., Journal of Society of Korea Industrial and Systems Engineering, Vol.30 (4) ; pp.182-189
    4. Y.B. Chung , H.J. Mun (2002) A Study on Process Capability Index using Reflected Normal Loss Function., Journal of The Korean Society for Quality Management, Vol.30 (3) ; pp.66-78
    5. Y.B. Chung (2005) A Multivariate Process Capability Index using Expected Loss., Journal of Society of Korea Industrial and Systems Engineering, Vol.28 (4) ; pp.116-123
    6. D. Drain , A.M. Gough (1996) Applications of the Upside- Down Normal Loss., IEEE Trans. Semicond. Manuf., Vol.9 (1) ; pp.143-145
    7. S.J. Hashemi , S. Ahmed , F. Khan (2015) Operational loss modeling for process facilities using multivariate loss functions., Chem. Eng. Res. Des., Vol.104 ; pp.333-345
    8. F. Khan , H. Wang , M. Yang (2016) Application of loss functions in process economic risk assessment., Chem. Eng. Res. Des., Vol.111 ; pp.371-386
    9. B.P.K. Leung , F.A. Spiring (2004) Some Properties of the Family of Inverted Probability Loss Functions., Qual. Technol. Quant. Manag., Vol.1 (1) ; pp.125-147
    10. H.J. Mun , Y.B. Chung (2002) A Study on Multivariate Process Capability Index using Quality Loss Function., Journal of Society of Korea Industrial and Systems Engineering, Vol.25 (2) ; pp.1-10
    11. H. Shahriari , N.F. Hubele , F.P. Lawrence (1995) A multivariate process capability vector, Proceeding of the 4th Industrial Engineering Research Conference, ; pp.304-309
    12. F.A. Spiring , A.S. Yeung (1998) A General Class of Loss Functions with Industrial Applications., J. Qual. Technol., Vol.30 (2) ; pp.152-162
    13. F.A. Spiring (1993) The Reflected Normal Loss Function., Can. J. Stat., Vol.21 ; pp.321-330
    14. W. Taam , P. Subbaiah , J.W. Liddy (1993) A Note on Multivariate Capability Indices., J. Appl. Stat., Vol.20 ; pp.339-351