The blade root fillets which have strong influences on the performance of propellers in the both structural and hydrodynamic points of view, are mechanical parts for smooth connection surface with a blade and a hub. A few related researches (Sabol, 1983; Kennedy, 1997) have noted that 3T-T/3 double radius section design would be suitable for reducing Stress Concentration Factor(SCF) and increasing Cavitation Inception Speed(CIS). In this paper, it is confirmed that this compound cross-section design has come close to the optimum solution in the shape optimization standpoint so that it could protect the propeller blade under the frequent and various loading cases. On that basis, we suggest the definite and simple fillet design methodology that has the cross-section with nT-T/n compound radius and elliptic shape which could sustain the given derivatives information as well as the offsets at the boundary and all inner region of the fillet surface. In addition, the result of design is presented in form of IGES file format in order to connect with NC machine seamlessly.
두 면을 부드럽게 잇는 연결 표면(smooth transition connecting surface)은 많은 공학 제품에서 발생하며 그 연결 형태와 분야에 따라 각각 챔퍼(chamfer), 블랜딩(blending), 혹은 필렛(fillet) 등으로 다양하게 불린다. 프로펠러의 날개(blade)와 허브(hub 혹은 boss, flange, disk)를 연결하는 부분에도 연결 부위가 필요한데 이를 조선 분야에서는 필렛으로 호칭한다(Carlton, 1994).
비록 프로펠러의 필렛 설계는 일체형 프로펠러의 생산 때문에 어쩔 수 없이 발생되는 부산 부위(by-product)로 여겨지기도 했지만, 최근 구조 역학적 관점에서 해당 부위에 발생하는 응력 집중을 경감(stress relief)시키기 위해 필렛 형상 단면을 규정하는 최소 반경(single radius type)을 허브와 날개가 만들어내는 단면의 최대 두께와 관련된 수치로 규정하거나 타원형(elliptic section type) 혹은 복합 반경(compound radius section type)등을 권장하는 규정이 각종 선급 (Lloyd's Register, 2007; Det Norske Veritas, 2015) 등에 등장한다.
한편, 유체역학적으로도 급속한 날개 정밀도의 개선에 따라 성능 관점에서 관심이 필렛 쪽으로 옮겨가고 있는데, Kennedy, et al. (1997)의 연구에서는 다양한 프로펠러 필렛 종류별로 캐비테이션(cavitation) 성능을 평가하였다. 그 결과, 대부분 날개앞날(blade leading edge) 부근의 날개코드(chord)의 20% 정도의 구간에서 캐비테이션 최초 발생(cavity inception)이 진행되고 앞으로 설명될 3T-T/3 단면이 그 외 제안된 다른 모든 단면들보다 우수함을 제시하고 있다. 다만 간편한 실험을 위해서 평면 허브(planar hub)와 날개형 허브단면(NACA hub section)만을 대상으로 상용 프로그램인 Surfcam을 사용하여 NURBS(Non-uniform Rational B-spline Surface)를 생성한 것이 아쉽다.
본 연구의 시작은 Kennedy가 CIS 측면에서 우수하다고 제안한 표준 3T-T/3 단면 (Kim, et al., 2007)의 필렛 설계가 응력집중인자(stress concentration factor)측면에서 얼마나 우수하며, 실제로 사용되는 일체형 프로펠러 생산에 적용 가능한 표면을 생성하는 방안에 관해서 알아보고자 함이었다. 참고로 Sabol (1983)은 그의 연구에서 표준 3T-T/3단면 필렛 설계는 Philadelphia Naval Shipyard에서 처음 제안되었다고 밝힌바 있다.
한편, 기계 역학 분야에서 특히 어께 형상 필렛(shoulder type fillet)은 다단축(stepped shaft)의 응력집중 완화를 위해 많은 연구 (Pilkey & Pilkey, 2008)가 진행되었다. 여기서, 응력집중이라 함은 구조 요소에 외력이 가해졌을 때 형상이 급격하게 변화되는 곳에서 국부적으로 매우 큰 응력이 발생하는 것을 일컫는다. 본 연구에서는 비틀림, 굽힘 혹은 인장 하중에 표준 3T-T/3 단면 형상 어깨 타입 필렛이 응력집중 완화 측면에서 얼마나 효과가 있는지 형상을 변화해 가면서 검증하고자 한다. 또한 표준 3T-T/3을 포함한 다양한 단면형상을 가지면서 날개코드 방향으로 기준 두께를 달리하는 날개 뿌리 필렛을 B-spline으로 생성할 수 있는 방안을 제시한다. 이에 관한 연구로는 Sabol (1983)과 Filkins (1991) 등이 표준 3T-T/3단면을 가지는 B-spline 필렛 표면을 생성하였지만 날개와 필렛이 닿는 부위 그리고 필렛과 허브가 접하는 부위(boundary curves)의 연속성에만 집중된 내삽입법(interpolation)설계안이었다. 본 연구에서는 표면끼리의 접촉 부분은 물론이고 표면 중간 지역에서의 연속성까지 고려한 구속된 최소자승 근사(constrained least square approximation) 기법을 사용하여 경계뿐만이 아닌 표면 어떤 부분도 주어진 단면 설계안에 최대로 가까운 표면을 도출할 수 있는 방안을 제시한다. 더욱이 특징적인 부분은 복합 표면 형상에 필요한 다양한 점들을 최적화 문제로 전환하여 빠르고 정확하게 얻어내는 방안에 관해서도 기술하고자 한다.
날개뿌리 필렛 표면 설계라 함은 보통 단면형상을 생성하고 날개코드 방향으로 기준두께의 변화량을 함수로 나타내는 과정을 거친다. 우선 Fig. 1과 같이 관련 용어를 정리해 본다.
CIS관점에서 Kennedy et al. (1997)가 권고하였고 응력집중 경감관점에서는 본 연구가 권장하는 프로펠러 표준 뿌리 필렛 단면 형상인 3T-T/3에 관해서 알아본다. 가장 기본적인 필렛 단면은 단일반경(single radius) 타입이다. 설계와 생산이 비교적 쉽지만 비교적 큰 반경을 확보하지 못하면 큰 응력집중이 발생하며, 선박의 운항 기간 동안 반복하중을 받는 프로펠러로서는 피로파괴가 우려된다. 또한 프로펠러의 본연의 추진성능을 방해할 수 있으므로 반경의 크기에는 한계가 있다. 그에 비해 3T-T/3단면은 복합반경(compound radius) 혹은 이중반경(double radius) 단면으로 불리고 nT과 T/n 두 반경으로 단면을 형성한다. 두 반경이 접하는 부분에서는 유체 및 구조역학적 성능을 고려하여 C1 연속성 이상을 확보해야 한다. 여기서 n은 스케일링 비례 상수이고 T는 날개뿌리의 최대 두께에 해당하는 기준값에 해당한다. 날개코드 방향으로 다른 기준값(T)을 가지면 변동반경 단면 타입(variable radius section type), 반대 개념은 고정반경 단면 타입(fixed radius section type)이라 호칭한다. 다음은 단면 형성과정을 차례대로 정리하고자 한다.
Step 1: 날개와 허브가 공유하는 이음매 곡선을 형성하는 과정이다. 우선 각각 식 (1)과 같이 B-spline surface으로 구축된 날개(
식(1)에서
Step 2: 이음매 곡선의 매개 변수(
Step 3: 각 단면의 기준면을 형성한다. 기준면 형성을 위해서는 해당 평면을 정의하는 그리고 상호 평행하지 않은 두 방향 벡터가 필요한데 본 연구에서는 점에 해당하는 에서의 법선방향 벡터 와 의
Step 4: nT원의 중심점과 T/n원의 중심점 그리고 과의 접점 등 단면 형성에 필요한 다양한 점들을 찾는다. nT원의 중심점은 Step 2에서 구한 점에서 Step 3에서 구한 방향으로 nT 크기만큼 이동하면 얻을 수 있다. 이점을 라 정의한다.
Step 5: nT원과 T/n원에서 매개변수인 각도를 변경해 가며 원하는 개수의 학습점(offset point)과 접선(tangential) 벡터를 획득한다. 접선 벡터를 찾는 과정도 식 (5)로부터 유도된 (12)으로 충분하다. 과정 1~5을 거치면, 필렛 표면을 만들기 위한 모든 학습점과 접선벡터를 획득할 수 있다. 이를 이용하여 구속된 최소자승 근사기법으로 결과 표면을 획득할 수 있다. 식 (12)는 Piegl and Tiller (2003)의 원을 B-spline를 통해 근사하는 법을 참고하였다. 식 (13)에서 Scale factor는 두 원의 연결점에서의 접선 벡터의 값을 일치시키기 위해서 조정된 값이다.
한편 위의 과정을 응용하면 타원형 단면(Fig. 3)도 생성할 수 있다. 타원형 단면은 이중 단면과 달리 이음부분의 무한 연속성이 보장되는 장점이 있는 반면, 2번째 대상표면(hub)을 뚫고 지날 수 있는 한계를 극복하지 못한다. 구체적인 설계 방안을 살펴보면, Step 1~3 과정은 동일하고 이중 단면 설계의 nT의 중심점에 해당하는 점을 구하기 위해서 T/n만큼 방향으로 전진하고 이점을 중심으로 식 (14)과 같이 중간점을 확보한다는 차이점만 존재한다.
필렛 단면을 완성하기 위해서는 2.2절에서 구했던 스팬(span section)방향의 정보도 필요하지만 날개코드방향의 정보도 필요하다. 날개코드 방향으로 제시되는 기준 두께를 변화시켜주는 형상 함수를 정의할 수 있다면 이를 ‘기준두께 가변 필렛 표면’(variant reference thickness surface)이라 정의한다. 날개코드 방향에 변화가 없는 필렛 단면은 흔히 상선의 고정피치프로펠러 (fixed pitch propeller, FPP)에서 사용되어 왔지만, 점차 필렛 표면 주변의 유체 압력변화에 관심이 증가되면서 변화가 생기고 있다. Sabol (1983)은 코드 중심에서 최대 두께를 가지고 L.E.와 T.E.에서 최대두께의 30%를 가지도록 식 (15)과 Fig. 4과 같이 다항식 형상함수를 제안하였다. 식 (15)의
기준 두께 가변 필렛 표면의 유체역학적 성능평가가 진행 중이어서 어떤 형상 함수가 우월한지는 아직 미정이지만 본 연구에서는 다항식을 벗어나 무한 연속성을 가지는 다음과 같은 함수를 식 (16)과 Fig. 5와 같이 제안한다. 본 식의 장점은 연속성에 관한 것 이외에도 프로펠러 압력면과 흡입면의 비대칭성을 보장하기 위해서 특정 구간에서의 기준 두께 값을 달리 설정할 수 있도록 고안된 점이다. 즉 유체성능을 극대화하는 프로펠러 필렛 설계의 설계변수로 흡입면용
2.4 구속된 최소자승 근사기법을 이용한 필렛 단면 설계 방안
주어진 옵셋 혹은 미분정보를 이용하여 곡선과 표면을 생성하는 방법 중에, 내삽법(interpolation)은 주어진 옵셋들을 정확하게 만족하는 방법 (Filkins, 1991), 즉 수치적으로 연립방정식을 푸는 형식이라면, 근사법(approximation)은 주어진 Basis function set를 통해서 가장 근접한 계수를 정해가는 과정이며 수치적으로 최적화가 필요하다. 본 연구에서는 Piegl and Tiller (1997)에서 자세히 소개된 “Weighted and Constrained Least Squares Curve Fitting”을 표면용으로 확장하여 강제적으로 만족시켜야하는 단면 양 끝단의 오프셋, 접선벡터, 곡률 그리고 되도록 만족시키려는 표면 중간의 오프셋과 접선벡터를 고려한 필렛 표면 표현법으로 구성된다. Fig. 6과 7은 단면 근사를 날개 코드를 따라 500번수행하고 단면마다 17개의 옵셋과 미분정보를 획득하여 u, v 방향으로 각각 25개, 10개의 제어점들을 확보한 필렛 표면 설계의 결과를 보여준다. Fig. 7의 곡률은 로컬 형상 반경의 역수이며 곡률 반경이 줄어들면 곡률 값은 늘어나고 해당 색상은 검정(0.001)에서 파랑(0.01), 초록(0.1)을 걸쳐 빨강(1.000)으로 변경된다.
Kennedy, et al. (1997)은 비록 평면 허브(planar hub) 혹은 날개 형상(foil shape) 허브들과 날개단면의 모양이 고정코드(constant chord), 무(無)캠버(uncambered), 꼬임 없는(untwisted) 프로펠러 날개들로 제한되어 있었지만, 다양한 날개 단면 형상을 대상으로 받음각(angle of attack)을 변경시키면서 캐비테이션 초생속도(CIS)에 미치는 영향을 실험하고 분석하였다. 결론적으로 3T-T/3 이중 단면 설계가 CIS 관점에서 상대적으로 가장 우수함을 밝히고 있으며, 이 경우 L.E.부터 0.2 코드 길이만큼 하류로 내려간 부근을 캐비테이션의 발생이 시작되는 곳으로 지목하고 있다.
프로펠러 날개 뿌리 부근에 생성되는 필렛 표면의 존재 근거는 원래 CIS를 증가시키려는 것이 아니었고, 구조적 응력 집중 감소를 위해서 시작되었다. 하지만 아직까지 3T-T/3 필렛 단면 설계가 어떤 근거로 제시되었는지 밝힌 연구는 없다. 본 연구에서는 직접적인 프로펠러-필렛-허브의 정밀한 모델링을 통해 구조해석을 수행하지는 않았지만, Sonmez (2009)의 형상 최적화 연구에서 아이디어를 얻어 어께 타입 필렛(shoulder fillet)의 단면 형상을 변화시켜 가면서 3T-T/3의 단면을 응력집중 관점에서 평가하였다.
프로펠러가 받는 응력이 비록 항복 응력을 충분히 만족한다 하더라도 비정상 하중의 빈도(number of cycles)가 높을 때 생기는 피로파괴를 고려해야 하므로 필렛 설계 시에 반드시 응력집중인자(
여기서
어떤 단면을 가지는 필렛 단면이 가장 우수한
* 수치 실험 1:
* 수치 실험 2:
Fig. 9의 Case (d)와 Fig. 11의 Case n=3 (d)가 3T-T/3 단면에 해당하는데 다양한 하중의 응력 집중인자 차원에서 수렴됨을 확인 할 수 있고, 결국 3T-T/3 이중 반경 단면이 구조 역학 측면에서도 우수한 형상임을 확인하였다.
주어진 프로펠러의 날개와 허브를 이용하여 필렛 표면을 생성하고자하는 유저가 제시된 프로그램을 통해 결과를 얻기 위한 과정을 설명한다. 참고로 날개와 허브는 IGES 파일 형식으로 B-Spline 혹은 NURBS로 정의되어도 아무런 상관이 없다. IGES 파일 포맷은 다양한 상용 프로그램이 공통으로 지원하고 있으며, 이를 통한 NC 코드의 생성에도 유리한 특징을 가지고 있다. OpenGL, MS C#등으로 구축되었으며 생성된 필렛 표면을 포함한 모든 결과를 IGES로 출력할 수 있게 하였다. 필렛 표면 설계자는 IGES 파일을 통해 2.1절에서 설명된 대상 표면을 확보하고 필렛 단면 형상을 선택한다. 기본적으로 nT-T/n의 이중 반경 단면을 지원하고 타원형과 유저-정의 단면이 가능하다. nT-T/n에서 n은 2, 3, 4 중에 선택할 수 있게 하였다(Fig. 13). 다음 단계는 2.3절 말미에 설명된 압축면과 흡입면의 최대 기준 두께를 설정한다. 마지막으로 프로펠러 날개코드 방향으로의 분할 개수를 설정하고 근사에 활용될 제어점들의 개수를 u, v방향으로 지정하면 모든 준비가 완료 된다(Fig. 13 참조). Fig. 14은 필렛 표면 설계가 종료된 후에 시각적으로 그 결과를 확인한다. 확인 후에 필렛 표면 단독으로 혹은 날개 및 허브 등의 기타 표면을 포함한 IGES 파일 포맷으로 Export가 가능하다.
본 연구에서는 프로펠러 날개의 뿌리 부분에 허브와 맞닿는 부위의 필렛 표면 설계를 위한 간단하고도 실용적인 방안을 제시한다. CIS관점에서 우수하다고 알려진 3T-T/3을 포함한 nT-T/n 단면 형태와 타원형 필렛 단면을 생성하고 날개코드방향의 함수를 설정하여 기준두께를 변화시킬 수 있게끔 하였다. 특히 nT-T/n의 이중 단면 근사를 위해서 양단 경계의 미분정보는 물론이고 곡 내부에서의 연속성도 확보할 수 있도록 하였으며, 단면 설계에 필요한 다양한 파라미터를 적절한 최적 공식화를 통해 짧은 시간에 단번에 획득할 수 있는 설계 절차를 제시했다.
CIS관점에서 우수한 단면으로 증명된 3T-T/3가 응력집중인자 차원에서도 최적점에 가까운 해임을 밝혔으며, 설계를 위해 개발된 SW는 계산에 필요한 모든 정보의 입출력은 IGES 파일 포맷으로 통일화하여 NC를 사용하는 절단 기기와의 연결이 원활토록 구축하였다.
본 연구에서는 필렛을 설계하기 위한 방안을 제시하고 있지만 동시에 설계를 위한 다양한 파라미터들도 소개하고 있다. 단면의 종류, 날개 코드방향으로의 L.E., T.E., Middle 날개코드 중심점에서의 기준 두께 등이 그 대표적인 예가 될 수 있다. 현재 진행되고 있는 연구에는 해당 파라미터의 최적 조합을 통해서 CIS 관점에서의 우수한 설계안을 도출하기 위해서 수치적으로는 고차 패널법의 활용, 실험적으로는 LCT(Large Cavitation Tank)에서의 캐비테이션 실험을 계획하고 있다.
