Journal Archive

차량용 현가장치는 차체의 무게를 지지하고 차량의 진동을 흡수하여 승차감을 향상시키고 차륜과 지면과의 적절한 접지력을 유지하여 주행안정성을 확보하는 등 다양한 기능을 수행하며 차량의 각종 성능과 직접 관련된 매우 중요한 부분이다. 현가장치는 주로 충격흡수기, 코일 스프링, 각종 링크 등으로 구성되어 있으며, 이 중에서 특히 충격흡수기는 위에서 언급한 다양한 기능을 수행하는 가장 핵심적인 부품이다. 이와 같이 충격흡수기가 차량에서 가지는 중요성 때문에 충격흡수기에 대한 정확한 수학적 모델링과 성능에 영향을 미치는 각종 파라미터들에 대한 연구는 지금까지도 꾸준히 연구^1-5)되어 왔다. 그러나 일반적으로 충격흡수기를 포함한 기계부품에는 제조과정에서 제품의 성능을 좌우하는 주요 설계변수의 변동이 존재할 수 있다. 이러한 설계변수의 변동 즉 불확실성은 치수공차, 형상공차, 재질의 물성치 변화 등의 원인에 기인하며 제품의 성능에도 중요한 영향을 미칠 수 있다. 설계변수의 불확실성이 제품의 성능에 미치는 영향에 대한 연구는 제품의 성능을 확률적으로 예측 가능하게 하여 제품 설계과정에서 설계의 신뢰성을 높일 수 있는 중요한 연구 분야이지만 현재까지는 아직 이러한 분야에 대한 연구^6),7)가 부족한 실정이다. 본 연구는 차량용 충격흡수기를 대상으로 하여 설계변수의 불확실성에 근거한 변동이 충격흡수기의 성능에 어떠한 영향을 미치는가에 대하여 통계적인 접근법을 이용하여 연구를 실시하였다.

본 연구의 대상인 충격흡수기의 구조와 작동원리를 Fig. 1과 Fig. 2에 나타내었다. 그림에서와 같이 충격흡수기는 피스톤에 설치된 피스톤밸브 조립체(piston valve assembly)와 본체 아래 부분에 설치된 바디밸브 조립체(body valve assembly)에 감쇠력을 발생시키는 주요 밸브들이 설치되어 있다. 또한 충격흡수기는 인장실, 압축실 그리고 저장실과 같은 3개의 공간(chamber)으로 구분되어 있으며 피스톤의 왕복운동에 의한 각 챔버에서의 압력 변동에 따라 작동유체가 밸브들을 통과하여 이동하고 이때 적절한 감쇠력이 발생하도록 설계되어 있다. 이러한 유체 유동에 대한 작동원리와 이론식은 이전의 연구에서 제시되었다.^8-9) 차량용 충격흡수기의 성능특성을 해석하기 위한 해석 모델을 Fig. 3에 나타내었다. 해석모델은 지멘스사의 AMESim 소프트웨어를 사용하여 구성하였고 그림에 주요 설계변수를 표시하였다.

Fig. 1 
Schematic diagram of automotive shock absorber

Fig. 2 
Detailed configuration and piston valve assembly of automotive shock absorber

Fig. 3 
Simulation model of shock absorber for the statistical design

그림에서와 같이 시뮬레이션 모델은 본체, 피스톤밸브 조립체, 바디밸브 조립체, 압축실, 인장실, 저장실 등으로 구성되어 있고 본체에 sin파형의 변위 입력신호가 가해지는 구조로 되어 있다.

특히, 저장실은 공기를 포함하는 공간이므로 공압 모델을 이용하여 모델링하였다. 이와 같은 시뮬레이션 모델의 설계변수들을 Table 1에 나타내었고 Fig. 4에 Table 1의 설계변수를 이용한 충격흡수기의 성능지표인 감쇠력-속도 및 감쇠력-변위 선도(F-V, F-S curve)의 계산 결과를 나타내었다. 그림에서 A, B 그리고 C점은 각각 압축행정에서의 최대감쇠력, 인장행정에서의 최대감쇠력, 블로오프 포인터(blow-off point)에서의 감쇠력으로 충격흡수기의 성능특성을 나타내는 주요한 성능지표이다.

Fig. 4 
Simulation results of F-V and F-S curve

설계변수의 불확실성에 의한 출력변수의 영향을 통계적 기법에 의하여 연구하기 위하여 Fig. 7과 같이 LHS에 의해 추출된 설계변수를 Fig. 3의 시뮬레이션 모델에 입력하여 300회의 시뮬레이션을 수행하였다. 이러한 과정을 거친 시뮬레이션 결과를 분석하기 위하여 변수들 간의 상관계수를 검토하였다. 상관계수(coefficient of correlation)는 변수 간의 관계의 강도를 수치로 표시한 것으로 –1에서 1의 범위를 가진다. 여기서 상관계수가 1에 가까울수록 양의 상관관계가 강하다고 하고, –1에 가까울수록 음의 상관관계가 강하다고 한다. 이때 음의 값은 한 변수가 증가할 때 다른 변수는 감소함을 의미하고 양의 값은 한 변수가 증가할 때 다른 변수도 같이 증가함을 의미한다. 또한 0에 가까운 값은 두 변수 사이에 상관관계가 거의 없음을 의미한다. 상관계수는 다음과 같은 피어슨 상관계수(Pearson correlation coefficient)를 사용하였다.

여기서 x¯와 y¯는 표본평균을 나타내고, x-x¯와 y-y¯는 평균과의 편차를 나타낸다. Fig. 8에 피어슨의 상관계수의 원리를 이용한 설계변수와 출력변수 간의 상관관계 해석 결과를 나타내었다. 이러한 상관관계 해석을 통하여 설계변수가 출력변수에 미치는 영향을 전체적이고 체계적으로 파악할 수 있음을 그림을 통하여 확인할 수 있다. 그림에서와 같이 압축행정에서의 최대감쇠력 R1에 가장 큰 영향을 미치는 변수는 X4로서 상관계수는 – 0.774를 나타내어 강한 음의 상관관계를 나타냄을 알 수 있다. 또한 인장행정에서의 최대감쇠력 R2는 X2, 그리고 블로오프 포인트에서의 감쇠력 R3은 X3과 가장 큰 상관관계를 가지며 각각의 상관계수는 – 0.838과 0.985를 나타내었다.

Fig. 8 
Correlation between design and response variables

다음으로 반응표면(response surface)법을 이용하여 통계적으로 추출한 설계변수 변동이 출력변수의 변동에 미치는 영향을 확인하였다. 반응표면법¹³⁾은 통계적 기법에 의해 추출된 값들을 유한요소 모델이나 AMESim과 같은 동적 해석모델에 입력하여 계산결과를 3차원 표면의 형식으로 표시한 것으로, 설계변수의 변동에 대한 출력변수의 변동을 근사적인 평면함수로 표시한 것이다. 반응표면은 최소자승법을 이용하여 회귀모델(regression model)로 나타내어진다. 이러한 반응표면은 최적화 계산 등에서 계산의 간략화를 위한 근사모델로도 사용된다. 이와 같은 반응표면법을 적용하기 위하여 Fig. 8의 상관계수 해석 결과를 참조하여 출력변수 R1, R2, R3에 가장 큰 영향을 미치는 변수들을 선정하였고 반응표면법으로 추정한 결과를 Fig. 9에 나타내었다. Fig. 9(a)에서와 같이 압축행정에서의 최대감쇠력 R1은 설계변수 X3에 의해 증가하는 경향을 보이고 X4에 의해 감소하는 경향을 보임을 알 수 있다. 이것은 Fig. 8의 결과와 잘 일치함을 알 수 있다. Fig. 9(b)와 (c)에서도 Fig. 8과 일치하는 결과를 보임을 알 수 있다. 반응표면법에 의한 Fig. 9(a)의 응답을 나타내는 결과 식을 식 (2)에 나타내었다.

Fig. 9 
Simulation results of surface response

또한 Fig. 10과 Table 4에 출력변수들의 확률밀도 그래프와 통계적 특성을 나타내었다. Table 4를 참조하면 압축행정에서의 최대감쇠력 R1은 932N에서 1,114N의 범위, 인장행정에서의 최대감쇠력 R2는 1,816N에서 2,526N의 범위, 블로오프 포인트에서의 감쇠력 R3은 312N에서 429N의 범위에 95.4%의 신뢰수준으로 분포하고 있음을 확인할 수 있다.

Fig. 10 
Histogram of response variables

또한 편포도(skewness)는 인장행정에서의 최대감쇠력 R2가 가장 크게 나타남을 알 수 있어 Fig. 10(b)에서와 같이 히스토그램이 평균을 중심으로 오른쪽으로 가장 많이 편향되어 있음을 알 수 있다.

이상과 같이 통계적인 해석 방법을 이용한 차량용 충격흡수기의 성능 특성에 대한 연구를 통하여 다음과 같은 결과를 얻을 수 있었다.

1) LHS 방법에 의한 설계변수의 표본추출 방법을 통하여 정규분포를 이루는 설계변수 표본을 추출하였고 이를 이용하여 충격흡수기 설계의 핵심지표인 F-V와 F-S 특성에 대한 통계적 특성을 파악할 수 있었다.

2) 설계변수에 대한 출력변수의 상관관계를 분석한 결과, 압축행정에서의 최대감쇠력 R1에 가장 큰 영향을 미치는 변수는 바디밸브 조립체의 블로오프 밸브 유량-압력 기울기 X4로서 상관계수 –0.774의 강한 음의 상관관계를 나타내었다. 또한 인장행정에서의 최대감쇠력 R2는 피스톤밸브 조립체의 블로오프밸브 유량-압력 기울기 특성 X2, 그리고 블로오프 포인트에서의 감쇠력 R3은 바디밸브 조립체의 블로오프밸브 크래킹 압력 특성 X3에 의해 가장 큰 영향을 받음을 알았고 각각의 상관계수는 –0.838과 0.985를 나타내었다.

3) 반응표면법에 의하여 출력변수의 변동을 회귀모델의 형태로 표시한 결과, 상관관계 해석결과와 잘 일치함을 확인할 수 있었다.