KR101671946B1 - 부등피치 재생 블로워 및 이의 최적화 설계 방법 - Google Patents

부등피치 재생 블로워 및 이의 최적화 설계 방법 Download PDF

Info

Publication number
KR101671946B1
KR101671946B1 KR1020140172727A KR20140172727A KR101671946B1 KR 101671946 B1 KR101671946 B1 KR 101671946B1 KR 1020140172727 A KR1020140172727 A KR 1020140172727A KR 20140172727 A KR20140172727 A KR 20140172727A KR 101671946 B1 KR101671946 B1 KR 101671946B1
Authority
KR
South Korea
Prior art keywords
design
objective function
regenerative blower
optimal solution
reaction surface
Prior art date
Application number
KR1020140172727A
Other languages
English (en)
Other versions
KR20160067402A (ko
Inventor
이경용
최영석
김진혁
정욱희
Original Assignee
한국생산기술연구원
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 한국생산기술연구원 filed Critical 한국생산기술연구원
Priority to KR1020140172727A priority Critical patent/KR101671946B1/ko
Priority to PCT/KR2015/013040 priority patent/WO2016089103A1/ko
Priority to DE112015005494.4T priority patent/DE112015005494B4/de
Priority to US15/533,175 priority patent/US10590938B2/en
Publication of KR20160067402A publication Critical patent/KR20160067402A/ko
Application granted granted Critical
Publication of KR101671946B1 publication Critical patent/KR101671946B1/ko

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D25/00—Pumping installations or systems
    • F04D25/02—Units comprising pumps and their driving means
    • F04D25/08—Units comprising pumps and their driving means the working fluid being air, e.g. for ventilation
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D23/00—Other rotary non-positive-displacement pumps
    • F04D23/008—Regenerative pumps
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D29/00—Details, component parts, or accessories
    • F04D29/18—Rotors
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D29/00—Details, component parts, or accessories
    • F04D29/66—Combating cavitation, whirls, noise, vibration or the like; Balancing
    • F04D29/661—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps
    • F04D29/666—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps by means of rotor construction or layout, e.g. unequal distribution of blades or vanes
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04D—NON-POSITIVE-DISPLACEMENT PUMPS
    • F04D5/00—Pumps with circumferential or transverse flow
    • F04D5/002—Regenerative pumps
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T50/00—Aeronautics or air transport
    • Y02T50/60—Efficient propulsion technologies, e.g. for aircraft

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

재생 블로워가 제공된다. 본 발명의 예시적인 실시예에 따른 재생 블로워는 원주방향으로 이격 배치되는 복수개의 날개를 포함하는 임펠러를 포함하는 재생 블로워에 있어서, 복수개의 날개는 각각의 날개 사이가 증분각도(△Θi)로 배열된다.
이때,
Figure 112016062054783-pat00014

N은 전체 날개의 수(N=2보다 큰 자연수),
Am은 날개 간 간격(등분각)의 분포 크기(0도<Am<360/N도),
i는 날개의 순번(i=1, 2, 3, 4,…N),
P1 및 P2는 주기에 영향을 미치는 요소(0<P1<N, 0≤P2≤N, 단 P1, P2는 실수)

Description

부등피치 재생 블로워 및 이의 최적화 설계 방법{Uneven pitch regenerative blower and an optimal design method thereof}
본 발명은 재생 블로워 및 이의 최적화 설계 방법에 관한 것이다.
일반적으로, 재생 블로워는 산업용 고압 송풍기(링 블로워)와 같이 상대적으로 적은 유량 및 높은 압력으로 기체를 이송하는 용도에 주로 사용되는 것으로, 최근에는 연료전지 시스템의 공기 공급기, 수소 재순환장치 등으로도 그 활용범위가 확장되고 있다.
이러한 재생 블로워는 저 유량, 고 양정을 요구하는 시스템의 공기 공급용 블로워로 사용되고 오픈 채널형(Open channel type)과 사이드 채널형(Side channel type)이 있다. 재생 블로워는 원판 형태의 회전하는 임펠러의 원주 방향으로 날개가 위치하고 있으며 날개의 홈과 케이스의 유로 사이에서 내부 순환유동을 통해 압력이 상승하는 원리로 작동한다.
한편, 재생 블로워는 양정의 상승을 위해 복수개의 날개가 필수적이며 이로 인해 날개를 통과하는 주파수, 즉 고주파 소음(BPF, Blade Passing Frequency) 및 소음(overall 소음)이 발생한다. 일반적으로 재생 블로워의 소음을 줄이기 위해 효율을 높여 상대적인 성능을 향상시켜 회전 수 감소로 인한 소음 저감을 실현할 수 있으나 한계가 존재한다.
또한, 가정용 및 의료용으로 재생 블로워가 사용될 경우 소음기를 장착하여 소음을 저감시키는 방법을 사용할 수 있으나 이는 원가가 상승하고 크기가 증가하며 소음기로 인한 약 10%의 유로의 손실이 발생할 수 있다.
다만, 종래의 재생 블로워는 날개 배열을 난수로 제어하기 때문에 날개 배열에 따른 소음 및 효율을 예상하거나 조절하는 것이 어려운 문제가 있었다.
또한, 종래의 재생 블로워의 날개의 부등피치 배열은 난수발생(random number)법을 사용하나 이러한 경우 배열에 대한 근거가 부족하고 조절이 어려운 문제가 있었다.
일본 특허공보 특허 제5001975호
본 발명의 일 실시예는 날개를 부등피치로 배열하여 날개 배열에 따른 소음 및 효율을 예상하거나 조절할 수 있는 재생 블로워 및 이의 최적화 설계 방법을 제공하고자 한다.
본 발명의 일 측면에 따르면 원주방향으로 이격 배치되는 복수개의 날개를 포함하는 임펠러를 포함하는 재생 블로워에 있어서, 상기 복수개의 날개는 각각의 날개 사이가 증분각도(△Θi)로 배열되되, 상기 증분각도는
Figure 112016062054783-pat00001
이고, 상기 Am, P1, P2는 27≤≤η≤≤32이고, 77dB(A)≤≤SPL≤83.7dB(A)를 동시에 만족할 수 있는 재생 블로워.
이때, η = (Pout - Pin)Q / (σω), SPL = 10log10(P / Pref)2, η=효율, SPL=음압 레벨(Sound pressure level), (Pout -Pin)=총 압력, Q=체적 유량, σ=토크, ω=각 속도, P=음압, Pref=기준압력(2*10-5Pa), N은 전체 날개의 수(N=2보다 큰 자연수), Am은 날개 간 간격(등분각)의 분포 크기(0도<Am<360/N도), i는 날개의 순번(i=1, 2, 3, 4,…N), P1 및 P2는 주기에 영향을 미치는 요소(0≤≤P1≤≤N, 0≤≤P2≤≤N, 단 P1, P2는 실수)일 수 있다.
삭제
삭제
삭제
삭제
삭제
삭제
삭제
또한, 상기 Am은 1도≤Am≤8.23도일 수 있다.
또한, 상기 P1, P2는 1≤P1≤38, 0≤P2≤39일 수 있다.
본 발명의 다른 측면에 따르면 전술한 재생 블로워의 최적화 설계 방법으로서, 설계 변수 및 목적 함수 선택 단계, 상기 설계 변수의 상한 및 하한 값을 결정하는 설계영역 선정단계, 상기 설계 영역에서 목적함수의 최적해를 얻는 단계를 포함하되 상기 설계 영역에서 목적함수의 최적해를 얻는 단계는 상기 설계 영역에서 라틴 하이퍼 큐브 샘플링을 통해 복수개의 실험점을 결정하는 단계 및 상기 복수개의 실험점에서 각각 공력성능시험과 소음시험을 통해 상기 목적 함수 값을 얻는 단계를 포함하는 재생 블로워의 최적화 설계 방법이 제공된다.
또한, 상기 설계 영역에서 목적 함수의 최적해를 얻는 단계는 상기 최적해가 타당한지 비교하는 단계를 더 포함할 수 있다.
또한, 상기 설계 변수 및 목적 함수 선택 단계에서 상기 설계 변수는 날개 간 간격의 분포 크기인 Am, 주기에 영향을 미치는 요소 P1 및 P2를 포함하고, 상기 목적 함수는 효율 η 및 음압 레벨 SPL을 포함할 수 있다.
또한, 상기 설계 변수의 상한 및 하한값을 결정하는 설계 영역 선정 단계에서 상기 Am은 1도≤Am≤8.23도이고, 상기 P1은 1≤P1≤38이며 상기 P2는 0≤P2≤39일 수 있다.
또한, 상기 설계 영역에서 목적함수의 최적해를 얻는 단계는 반응면 기법을 사용하여 최적해를 산출할 반응면을 구성하는 단계를 포함할 수 있다.
삭제
또한, 상기 반응면 기법을 사용하면, 상기 목적 함수들의 RSA 모델의 함수 형태는, η= - 18.8659 - 17.9578Am - 10.5773P1 - 21.7493P2 + 7.3846AmP1 + 17.3858AmP2 - 0.789P1P2 + 6.2258Am2 + 11.0769P1 2 + 16.1141P2 2, SPL= 84.2304 + 4.2557Am -11.8326P1 -6.4429P2 + 8.2626AmP1 + 4.8169AmP2 + 5.9802P1P2 - 4.2959Am2 + 4.7855P1 2 + 1.2078P2 2일 수 있다.
또한, 상기 반응면 기법을 사용하여 최적해를 산출할 반응면을 구성하는 단계 이후에 다중 목적 진화 알고리즘을 이용하여, 상기 반응면 기법에 의해 얻어진 목적 함수들의 반응면들을 토대로 각 목적 함수들을 최대화시킬 수 있는 최적해를 얻을 수 있다.
또한, 상기 각 목적 함수들을 최대화시킬 수 있는 최적해를 얻는 단계 이후에 탐색 알고리즘인 SQP(sequential quadratic programming)를 이용하여 상기 최적해를 각 목적 함수들의 국부적인 검색을 통해 좀더 개선된 값을 구할 수 있다.
또한, 상기 최적해가 타당한지 비교하는 단계는 반응면 기법으로부터 구성된 각 목적 함수들의 반응면에 대한 분산 분석(ANOVA) 및 회귀 분석을 포함할 수 있다.
본 발명의 일 실시예에 따른 재생 블로워 및 이의 최적화 설계 방법은 다중 목적 최적화를 통해 설계됨으로써 효율 및 소음을 선택적으로 조절할 수 있도록 형성된다.
도 1은 본 발명의 일 실시예에 따른 재생 블로워를 도시한 개략도이다.
도 2는 본 발명의 일 실시예에 따른 재생 블로워의 임펠러를 도시한 평면도이다.
도 3은 본 발명의 일 실시예에 따른 재생 블로워의 임펠러의 변형예를 도시한 사시도이다.
도 4는 도 3의 단면을 도시한 단면도이다.
도 5는 본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법을 도시한 순서도이다.
도 6은 본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법에서 목적함수 효율 및 음압 레벨을 도시한 그래프이다.
도 7은 본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법에서 설계 변수들의 상관관계를 도시한 그래프이다.
이하, 첨부한 도면을 참고로 하여 본 발명의 실시예에 대하여 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 상세히 설명한다. 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며 여기에서 설명하는 실시예에 한정되지 않는다. 도면에서 본 발명을 명확하게 설명하기 위해서 설명과 관계없는 부분은 생략하였으며, 명세서 전체를 통하여 동일 또는 유사한 구성요소에 대해서는 동일한 참조부호를 붙였다.
이하에서는 도면을 참조하여 본 발명의 일 실시예에 따른 재생 블로워 및 이의 최적화 설계 방법을 보다 상세히 설명하도록 한다.
도 1은 본 발명의 일 실시예에 따른 재생 블로워를 도시한 사시도이다. 도 2는 본 발명의 일 실시예에 따른 재생 블로워의 날개가 부등피치로 배열된 임펠러를 도시한 사시도이다.
도 1 및 도 2를 참조하면, 본 발명의 일 실시예에 따른 재생 블로워(1)는 임펠러(70), 제1 케이싱(10), 제2 케이싱(30) 및 모터(50)를 포함할 수 있다.
도 1을 참조하면, 본 발명의 일 실시예에 따른 재생 블로워(1)는 좌우 양측으로 분할되는 한 쌍의 제1 케이싱(10)과 제2 케이싱(30)의 내부에 임펠러(70)가 회전 가능하게 설치된다. 이때, 임펠러(70)는 모터(50)의 회전축(미도시)에 설치되어 모터에 의해 회전된다.
도 3은 본 발명의 일 실시예에 따른 재생 블로워의 임펠러의 변형예를 도시한 사시도이다. 도 4는 도 3의 단면을 도시한 단면도이다.
이하에서는 본 발명의 일 실시예에 따른 재생 블로워의 임펠러를 설명한다.
본 발명의 일 실시예에 따른 재생 블로워(1)의 임펠러(70)는 원판(71), 복수개의 날개(73)를 포함할 수 있다.
한편 도 2 내지 도 4를 참조하면, 원판(71)은 재생 블로워(1)의 회전축(미도시)에 고정되기 위해 중심부에 축고정부(71a)가 형성되고 도 2에 도시된 바와 같이 일 측면 또는 도 3 및 도 4에 도시된 바와 같이 양 측면에 복수개의 날개가 원주방향으로 이격 배치될 수 있다.
다만 본 발명의 일 실시예에 따른 재생 블로워(1)는 원판의 하나의 면에 복수개의 날개가 원주방향으로 이격 배치된 재생 블로워에 대해서 이하 설명한다. 다만 이에 한정되지는 않고 도 3 및 도 4에 도시된 바와 같이 양면이 복수개의 날개가 원주방향으로 이격 배치될 수 있다.
축고정부(71a)는 재생 블로워(1)의 회전축 즉 모터의 회전축이 고정되어 회전축과 함께 원판(71)이 회전되도록 한다.
유로홈(75)은 복수개의 날개 사이에 형성되어 단면이 반원 또는 반 타원형일 수 있으나 이에 한정되지는 않는다. 유로홈(75)은 복수개의 날개 사이에 형성되므로 복수개가 서로 이격 배치된다.
한편 복수개의 날개(73)는 각각의 날개 사이의 각도(Θi)들이 동일한 등피치가 아닌 각도들이 동일하지 않은 부등피치로 배열된다.
본 발명의 일 실시예에 따른 재생 블로워는 날개 간의 각도가 식 1에 의한 증분각도(△Θi )로 설정됨으로써 부등피치로 날개를 배열할 수 있다.
Figure 112014118032127-pat00002
N은 전체 날개의 수(N=2보다 큰 자연수),
Am은 날개 간 간격(등분각)의 분포 크기(0도<Am<360/N도),
i는 날개의 순번(i=1, 2, 3, 4,…N),
P1 및 P2는 주기에 영향을 미치는 요소(0≤P1≤N, 0≤P2≤N, 단, P1, P2는 실수)
이때, 기준형상은 날개 사이의 각도가 동일하여 날개가 등 피치로 이격 배치되는 등 피치 임펠러이고 증분각도(△Θi)의 합은 360도를 만족해야 한다.
한편, 증분각도(△Θi)에 의해 임펠러(70)는 날개(73)수가 변하더라도 같은구조의 부등피치 조건을 만들 수 있을 뿐만 아니라 생성된 함수들은 (-1)i항 때문에 진동발산함수의 형태를 띄고 있어 증분각도의 평균이 전체 평균과 유사하도록 만들 수 있게 된다.
본 발명의 일 실시예에 따른 재생 블로워(1)는 날개(73)와 근접하는 격벽을 통과하는 날개의 시간 간격이 흩어져서 고주파 음이 줄어들 뿐만 아니라 음압이 여러 주파수대로 분산되어 고주파 영역에서 BPF가 감소되게 된다.
한편, 예를 들어, 전체 날개수가 N=39일 경우 여기서 전체 날개각도의 평균값은 360도 / 39 = 9.2도가 된다.
상기의 수학식들에 제시된 조건을 만족하기 위하여 날개 간 간격(등분각)의 분포 크기인 Am과 주기에 영향을 미치는 요소인 P1, P2 값을 제어하게 된다. 이러한 Am 값과 P1, P2 값을 제어하여 난수조건(Random Pitch Condition)과 비슷한 피치 조건 및 일정 간격을 가진 피치 조건을 생성할 수 있으므로 날개 배열에 대한 예상 및 조절을 용이하게 할 수 있게 된다.
도 5는 본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법을 도시한 순서도이다.
본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법에서는 다중 목적 최적화를 이용하여 날개의 간격을 부등피치로 변형시킴으로써 재생 블로워의 효율 및 소음을 동시에 조절할 수 있다.
본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법에서 최적화란 기준 형상인 등피치 임펠러와 비교할 때, 효율 및 소음을 동시에 좋게 할 수도 있고, 효율만 좋게 하거나, 또는 소음만을 좋게 할 수 있는 등 필요에 따라 조절할 수 있다는 것이다. 이를 위하여 본 발명의 일 실시예에 따르면 재생 블로워의 최적화 설계 방법은 설계 변수 및 목적 함수 선택 단계(S10), 설계 변수의 상한 및 하한 값을 결정하는 설계 영역 선정 단계(S20), 설계 영역에서 목적 함수의 최적해를 얻는 단계(S30) 및 최적해 비교 단계(S40)를 포함할 수 있다.
본 발명의 일 실시예에 따른 재생 블로워의 최적화 설계 방법에서는 재생 블로워(10)에서 설계 변수를 선정하고 설계영역 내에서 목적 함수를 최적화한다.
먼저, 설계 변수 및 목적 함수 선택 단계(S10)에서는 공력 및 소음 성능 시험을 통해서 목적 함수를 구하고, 이렇게 구한 목적 함수를 최적화하기 위해 날개의 부등피치를 결정하기 위한 설계 변수를 선정한다.
본 실시예에서, 설계 변수는 Am, P1 및 P2인데, Am은 날개 간 간격(등분각)의 분포 크기(0도<Am<360/N도)이고, P1 및 P2는 주기에 영향을 미치는 요소(0<P1<N, 0≤P2≤N, 단 P1, P2는 실수)이다.
날개(73)의 부등피치에 관련된 기하학적 매개 변수들(Am, P1 및 P2)은 재생 블로워(1)에서 효율(η) 및 음압 레벨(SPL)을 동시에 최적화하기 위한 설계 변수들로 사용될 수 있다. 이때, 설계 변수들의 범위를 확립함으로써 형성된 이동 가능한 설계 공간을 찾는 것이 중요하다.
또한, 본 발명의 일 실시예에 따른 재생 블로워(1)는 날개의 부등피치의 형상을 최적화시켜 효율 및 소음을 동시에 최적화하는 목적이 있으므로 목적 함수는 효율 η 및 음압 레벨 SPL로 설정될 수 있다.
그 후, 설계 변수의 상한 및 하한값을 결정하는 설계 영역 선정 단계(S20)에서는 최적 설계 수행을 위해 설계 변수의 범위를 한정함으로써 적절한 설계 영역을 설정한다.
최적 설계 과정에서 변경될 각 설계 변수들의 상한과 하한은 복수개의 날개를 포함하는 임펠러를 제작할 때 사용되는 드릴이나 날 등의 최소두께를 통해 결정될 수 있는데 본 발명의 발명자에 의하여 선정된 각 설계 변수들을 식 1에 대입하여 상한과 하한을 구하면 하기 표 1과 같다.
변수들 최소값 최대값
Am 1도 8.23도
P1 1 38
P2 0 39
즉, 본 발명의 일 실시예에서, 설계 변수 Am은 1도 이상 8.23도 이하이고, P1은 1이상 38이하이며, P2는 0이상 39이하이다.
그 후, 선정된 설계 영역에서의 실험 단계(S30)에서는 선정된 설계 영역에서 실험을 수행하여 예를 들어 30개의 실험 점에서의 목적함수 값을 결정한다.
이때, 30개의 실험 점은 다차원의 분포를 갖는 설계 영역에서 특정 실험 점을 샘플링하는데 유용한 라틴 하이퍼 큐브 샘플링(LHS)에 의해 결정될 수 있다. 30개의 실험 점에서의 목적함수 η 및 SPL값은 각각 공력성능시험과 소음시험을 통해 얻어질 수 있다.
실험결과를 통해 설계 영역에서 목적함수의 최적해를 얻는 단계(S40)에서는 대리모델의 일종인 반응면 기법을 사용하여 최적점을 산출할 반응면을 구성할 수 있다.
본 발명의 일 실시예에 따른 재생 블로워(10)의 다중 목적 최적화를 통하여 재생 블로워(10)의 다양한 유체역학적 성능을 높일 수 있다. 최적화의 목적은 재생 블로워의 효율(η) 및 음압 레벨(SPL)를 동시에 최적화하기 위한 것이다. η 및 SPL은 재생 블로워의 설계 최적화를 위한 목적 함수로서 다음과 같이 규정될 수 있다.
Figure 112014118032127-pat00003
Figure 112014118032127-pat00004
이때, η= 효율, SPL= 음압 레벨(Sound Pressure Level), (Pout-Pin)=총 압력, Q=체적 유량, σ=토크, ω=각 속도, P=음압, Pref=기준압력(2*10-5Pa)이다.
반응면 기법은 실제 반응함수를 근사적인 다항 함수로 모델링하기 위해 물리적인 실험 또는 수치적인 계산으로부터 얻어진 결과들을 이용하는 일련의 수리통계적인 기법이다.
반응면 기법은 한정된 수의 실험만으로도 임의 공간 내 반응을 모델링함으로써 실험의 시행횟수를 줄일 수 있다. 여기서 사용된 2차 다항식으로 구성된 반응면은 다음과 같이 표현될 수 있다.
Figure 112014118032127-pat00005
여기서, C는 회귀 분석 계수, n은 설계변수의 개수 및 x는 설계변수를 나타낸다.
이때, 회귀 분석 계수는 하기 식 5와 같다.
Figure 112014118032127-pat00006
이때, 본 발명의 일 실시예에 따른 목적 함수들의 RSA모델의 함수 형태는 정상화된 설계 변수들에 관해 아래와 같이 표현될 수 있다.
Figure 112014118032127-pat00007
Figure 112014118032127-pat00008
그 후, 상기 식 6 및 식 7를 만족하는 η 및 SPL을 구한다.
한편, 본 발명의 일 실시예에 있어서, η 및 SPL을 동시에 최적화하기 위해서 반응면 기법을 통해 얻어진 각 목적함수들의 반응면들을 토대로 각 목적 함수들을 최대화시킬 수 있는 다중 목적 진화 알고리즘을 사용할 수 있다.
다중 목적 진화 알고리즘으로서, Deb에 의해 개발된 real coded NSGA-Ⅱ 코드가 사용될 수 있다. 여기서 real coded는 NSGA-Ⅱ의 반응을 구성하기 위해 실제 설계 공간 내 교차 및 변이가 수행되는 것을 의미한다.
다중 목적 진화 알고리즘을 통하여 얻어진 최적점들은 비지배해들의 집합체인 파레토 최적해라 부른다. 이 파레토 최적해를 통해 사용하고자 하는 목적의 의도에 따라 원하는 최적점을 선택할 수 있다.
다중 목적 진화 알고리즘은 공지된 방법이므로 상세한 설명은 생략한다.
한편 라틴 하이퍼 큐브 샘플링(LHS)에 의해 얻어진 실험점들에 대한 목적함수의 값을 평가하고 평가된 목적 함수들을 바탕으로 하여 SQP(Sequential Quadratic Programming)를 사용하여 최적점을 탐색할 수 있다.
각 목적함수의 최적해들은 구배기반의 탐색 알고리즘인 SQP(sequential quadratic programming)를 사용하여 초기 NSGA-Ⅱ에 의해 예측된 해들로부터 각 목적함수들의 국부적인 검색을 통해 좀더 개선된 각각의 최적해들을 얻을 수 있다.
이때, SQP는 비선형 제약조건 내에서 비선형 목적 함수를 최적화하기 위한 방법으로 공지된 방법이므로 상세한 설명은 생략한다.
이렇게 개선된 최적해들로부터 지배해들은 버리고, 중복되는 해들은 제거하여 결과적으로 비지배해들의 집합체인 파레토 최적해를 얻을 수 있다. 파레토 최적해 중에 집단으로 분류한 단위체를 클러스터(Cluster)라 한다.
도 6은 본 발명의 일 실시예에 따른 재생 블로워의 다중목적 수치 최적설계로부터 도출된 파레토 최적해(클러스터드 최적해, COSs)의 효율 및 음압 레벨을 도시화 그래프이다.
도 6을 참조하면, 효율과 소음에 관한 목적 함수값들이 최적화됨에 따라 파레토 최적해는 S자 형상일 수 있다. 거래분석(trade-off analysis)은 2개의 목적 함수 사이에 상관관계를 보여준다.
따라서 본 발명의 일 실시예에 따른 재생 블로워(1)에서 더 높은 효율은 더 높은 소음에서 얻을 수 있고, 반대로 더 낮은 효율은 더 낮은 소음에서 얻을 수 있다.
도 6에 도시된 바와 같이, 이때 Am, P1 및 P2는 27≤η≤32이고, 77dB(A)≤SPL≤83.7dB(A)를 동시에 만족할 수 있고 하기 표 2는 이를 만족하는 Am, P1 및 P2의 값이며 도 6에서 파레토 최적해들의 그래프와 대응하는 값이다.
Figure 112014118032127-pat00009
Figure 112014118032127-pat00010
Figure 112014118032127-pat00011
Figure 112014118032127-pat00012
이때, 하기 표 3는 효율 및 소음이 동시에 최적인 집단인 Cluster A, B, C, D, E에 대한 최적 설계 변수들 Am, P1 및 P2의 값을 나타낸다. 이때 기준 형상은 효율 η 이 27.25이고 SPL은 79dB(A)이다.
설계 설계 변수들
Am P1 P2
기준 형상 0.000 0.000 0.000
Cluster A 1 23.96992 37.72269
Cluster B 1 20.31293 26.94253
Cluster C 1.975457 18.18757 23.56059
Cluster D 3.27427 15.95297 18.60822
Cluster E 6.793103 12.29705 1.858063
표 3은 최적점 A부터 E까지 이동함에 따라 설계 변수 Am은 증가하고 P1 및 P2은 감소하는 경향이다. 다만 P2가 감소하는 기울기가 P1이 감소하는 기울기보다 더 크다. 거래 분석에서 3개의 설계 변수 중 Am은 비례하고, P1 및 P2들은 반비례 관계를 보임을 확인할 수 있다.
이때, 기준 형상은 날개가 등 피치이므로 Am, P1 및 P2는 0,(도 6에서 삼각형으로 표시된 지점)이고, Cluster A는 Am = 1, P1 = 23.96992, P2 = 37.72269이고, Cluster B는 Am = 1, P1 = 20.31293, P2 = 26.94253이고, Cluster C는 Am = 1.975457, P1 = 18.18757, P2 = 23.56059이고, Cluster D는 Am = 3.27427, P1 = 15.95297, P2 = 18.60822이고, Cluster E는 Am = 6.793103, P1 = 12.29705, P2 = 1.858063이다.
도 6 및 도 7을 참조하면, 3개의 최적 설계 변수들은 기준 형상의 값에 비해 현저히 변할 수 있고, 효율 및 소음은 모든 최적점(COSs)에서 상당히 개선되어 효율 및 음압 레벨의 값을 선택할 수 있다.
따라서, 최적점 A부터 E까지 변화하는 동안 소음 및 효율은 증가하며 최적점(COSs) A는 가장 낮은 소음과 효율을 최적점(COSs) E는 가장 높은 소음과 효율을 나타냄을 알 수 있다.
본 발명의 일 실시예에 따른 최적점 비교 단계(S50)에서 반응면 기법으로부터 구성된 각 목적 함수들의 반응면에 대한 분산 분석(ANOVA) 및 회귀 분석을 통해서 구한 최적점들이 신뢰할 수 있는지 검토한다.
하기의 표 4는 분산 분석 및 회귀 분석의 결과를 나타낸다.
목적 함수들 R2 R2 adj 평균 제곱근 오차 교차 증명 오차
η 0.977 0.948 4.73×10-1 7.50×10-1
SPL 0.898 0.933 5.49×10-1 9.40×10-1
여기서 R2은 최소 제곱 표면 피팅에서 상관 계수를 나타내고, R2 adj값은 최소 제곱 표면 피팅에서의 조절된 상관 계수를 나타낼 수 있다. 이때, Ginuta는 반응면 기법에 의한 반응모델이 정확히 예측된 경우 R2 adj값은 0.9 이상 1 이하의 값을 가진다고 하였다.
평균 제곱근 오차는 실험이나 관측에서 나타나는 오차를 제곱해서 평균한 값을 의미하고, 교차 증명 오차는 예측한 오류를 계산하는 기법이다.
본 발명의 일 실시예에 따른 최적점 비교 단계(S50)에서 계산된 각 목적 함수인 효율과 소음의 R2 adj값이 각각 0.948 및 0.933이었으며, 이에 따라 반응면이 신뢰할 수 있다고 판단될 수 있다.
본 발명의 일 실시예에 따른 재생 블로워 및 이의 최적화 설계 방법은 다중 목적 최적화를 통해 날개가 부등피치로 배열되어 효율 및 소음을 선택적으로 조절할 수 있도록 형성된다.
이상에서 본 발명의 일 실시예에 대하여 설명하였으나, 본 발명의 사상은 본 명세서에 제시되는 실시 예에 제한되지 아니하며, 본 발명의 사상을 이해하는 당업자는 동일한 사상의 범위 내에서, 구성요소의 부가, 변경, 삭제, 추가 등에 의해서 다른 실시 예를 용이하게 제안할 수 있을 것이나, 이 또한 본 발명의 사상범위 내에 든다고 할 것이다.
1 : 재생 블로워 10 : 제1 케이싱
30 : 제2 케이싱 50 : 모터
70 : 임펠러 71 : 원판
71a : 축고정부 73 : 날개
75 : 유로홈

Claims (14)

  1. 원주방향으로 이격 배치되는 복수개의 날개를 포함하는 임펠러를 포함하는 재생 블로워에 있어서,
    상기 복수개의 날개는 각각의 날개 사이가 증분각도(△Θi)로 배열되되, 상기 증분각도는
    Figure 112016062054783-pat00013
    이고, 상기 Am, P1, P2는 27≤≤η≤≤32이고, 77dB(A)≤≤SPL≤83.7dB(A)를 동시에 만족할 수 있는 재생 블로워.
    이때, η = (Pout - Pin)Q / (σω), SPL = 10log10(P / Pref)2, η=효율, SPL=음압 레벨(Sound pressure level), (Pout -Pin)=총 압력, Q=체적 유량, σ=토크, ω=각 속도, P=음압, Pref=기준압력(2*10-5Pa), N은 전체 날개의 수(N=2보다 큰 자연수), Am은 날개 간 간격(등분각)의 분포 크기(0도<Am<360/N도), i는 날개의 순번(i=1, 2, 3, 4,…N), P1 및 P2는 주기에 영향을 미치는 요소(0≤P1≤N, 0≤P2≤N, 단 P1, P2는 실수)
  2. 삭제
  3. 제1 항에 있어서,
    상기 Am은 1도≤Am≤8.23도인 재생 블로워.
  4. 제1 항에 있어서,
    상기 P1, P2는 1≤P1≤38, 0≤P2≤39인 재생 블로워.
  5. 제1 항, 제3 항 및 제4 항 중 어느 한 항에 따른 재생 블로워의 최적화 설계 방법으로서,
    설계 변수 및 목적 함수 선택 단계;
    상기 설계 변수의 상한 및 하한 값을 결정하는 설계영역 선정단계;및
    상기 설계 영역에서 목적함수의 최적해를 얻는 단계를 포함하되,
    상기 설계 영역에서 목적함수의 최적해를 얻는 단계는 상기 설계 영역에서 라틴 하이퍼 큐브 샘플링을 통해 복수개의 실험점을 결정하는 단계 및
    상기 복수개의 실험점에서 각각 공력성능시험과 소음시험을 통해 상기 목적 함수 값을 얻는 단계를 포함하는 재생 블로워의 최적화 설계 방법.
  6. 제5 항에 있어서,
    상기 설계 영역에서 목적 함수의 최적해를 얻는 단계에서 상기 최적해가 타당한지 비교하는 단계를 더 포함하는 재생 블로워의 최적화 설계 방법.
  7. 제5 항에 있어서,
    상기 설계 변수 및 목적 함수 선택 단계에서 상기 설계 변수는 날개 간 간격의 분포 크기인 Am, 주기에 영향을 미치는 요소 P1 및 P2를 포함하고,
    상기 목적 함수는 효율 η 및 음압 레벨 SPL을 포함하는 재생 블로워의 최적화 설계 방법.
  8. 제5 항에 있어서,
    상기 설계 변수의 상한 및 하한값을 결정하는 설계 영역 선정 단계는 상기 Am은 1≤Am≤8.23이고, 상기 P1은 1≤P1≤38이며 상기 P2는 0≤P2≤39을 포함하는 재생 블로워의 최적화 설계 방법.
  9. 삭제
  10. 제5 항에 있어서,
    상기 설계 영역에서 목적함수의 최적해를 얻는 단계는 반응면 기법을 사용하여 최적해를 산출할 반응면을 구성하는 단계를 포함하는 재생 블로워의 최적화 설계 방법.
  11. 제10 항에 있어서,
    상기 반응면 기법을 사용하면, 상기 목적 함수들의 RSA 모델의 함수 형태는,
    η= - 18.8659 - 17.9578Am - 10.5773P1 - 21.7493P2 + 7.3846AmP1 + 17.3858AmP2 - 0.789P1P2 + 6.2258Am2 + 11.0769P1 2 + 16.1141P2 2,
    SPL= 84.2304 + 4.2557Am -11.8326P1 -6.4429P2 + 8.2626AmP1 + 4.8169AmP2 + 5.9802P1P2 - 4.2959Am2 + 4.7855P1 2 + 1.2078P2 2인 재생 블로워의 최적화 설계 방법.
  12. 제11 항에 있어서,
    상기 반응면 기법을 사용하여 최적해를 산출할 반응면을 구성하는 단계 이후에 다중 목적 진화 알고리즘을 이용하여, 상기 반응면 기법에 의해 얻어진 목적 함수들의 반응면들을 토대로 각 목적 함수들을 최대화시킬 수 있는 최적해를 얻는 재생 블로워의 최적화 설계 방법.
  13. 제12 항에 있어서,
    상기 각 목적 함수들을 최대화시킬 수 있는 최적해를 얻는 단계 이후에 탐색 알고리즘인 SQP(sequential quadratic programming)를 이용하여 상기 최적해를 각 목적 함수들의 국부적인 검색을 통해 좀더 개선된 값을 구하는 재생 블로워의 최적화 설계 방법.
  14. 제6 항에 있어서,
    상기 최적해가 타당한지 비교하는 단계는 반응면 기법으로부터 구성된 각 목적 함수들의 반응면에 대한 분산 분석(ANOVA) 및 회귀 분석을 포함하는 재생 블로워의 최적화 설계 방법.
KR1020140172727A 2014-12-04 2014-12-04 부등피치 재생 블로워 및 이의 최적화 설계 방법 KR101671946B1 (ko)

Priority Applications (4)

Application Number Priority Date Filing Date Title
KR1020140172727A KR101671946B1 (ko) 2014-12-04 2014-12-04 부등피치 재생 블로워 및 이의 최적화 설계 방법
PCT/KR2015/013040 WO2016089103A1 (ko) 2014-12-04 2015-12-02 부등피치 재생 블로워 및 이의 최적화 설계 방법
DE112015005494.4T DE112015005494B4 (de) 2014-12-04 2015-12-02 Regenerativ-Gebläsemaschine mit irregulärem Abstand und Design-Optimierungsmethode dafür
US15/533,175 US10590938B2 (en) 2014-12-04 2015-12-02 Irregular-pitch regenerative blower and optimization design method for same

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020140172727A KR101671946B1 (ko) 2014-12-04 2014-12-04 부등피치 재생 블로워 및 이의 최적화 설계 방법

Publications (2)

Publication Number Publication Date
KR20160067402A KR20160067402A (ko) 2016-06-14
KR101671946B1 true KR101671946B1 (ko) 2016-11-16

Family

ID=56091988

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020140172727A KR101671946B1 (ko) 2014-12-04 2014-12-04 부등피치 재생 블로워 및 이의 최적화 설계 방법

Country Status (4)

Country Link
US (1) US10590938B2 (ko)
KR (1) KR101671946B1 (ko)
DE (1) DE112015005494B4 (ko)
WO (1) WO2016089103A1 (ko)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107991200B (zh) * 2017-11-24 2020-02-14 北京航空航天大学 一种钛合金叶轮的疲劳寿命预测方法
CN109723658B (zh) * 2018-12-26 2020-09-11 上海洞彻智能科技有限公司 一种蒸发风机
CN112069607A (zh) * 2020-07-17 2020-12-11 北京动力机械研究所 整体叶轮成组分类编码与几何特征参数计算方法及装置
CN113048086B (zh) * 2021-03-18 2022-05-20 江苏大学 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法
KR102380405B1 (ko) * 2021-06-09 2022-03-30 김윤성 유로확장형 스웹트 프란시스 베인을 갖는 임펠러의 설계방법, 이에 의해 제조된 임펠러, 및 유로확장형 스웹트 프란시스 베인을 갖는 임펠러를 구비한 수중펌프
CN113642121A (zh) * 2021-07-26 2021-11-12 南京工业大学 一种基于响应面设计和多目标进化算法的铝合金制动钳铸造工艺参数优化方法
CN114841032B (zh) * 2021-09-29 2024-03-22 杭州汽轮控股有限公司 一种燃气轮机热部件寿命稳健性的设计方法
CN114218688B (zh) * 2021-10-28 2024-04-12 北京建筑大学 一种通风式制动盘的分段式倾斜槽叶片特征参数优化方法
CN115263774A (zh) * 2022-06-24 2022-11-01 烟台东德实业有限公司 一种转子分体式旋涡型氢气循环泵
CN116776600B (zh) * 2023-06-21 2024-04-12 安徽工业大学 基于自适应代理模型的风力机叶片优化设计方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5001975B2 (ja) * 2008-08-29 2012-08-15 ヒョンダム インダストリアル カンパニー リミテッド 燃料ポンプ用不等ピッチインペラ

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3708336C2 (de) * 1987-03-14 1996-02-15 Bosch Gmbh Robert Laufrad zum Fördern eines Mediums
JP2003336591A (ja) 2002-03-13 2003-11-28 Aisan Ind Co Ltd ウエスコ式ポンプ
JP2003278684A (ja) 2002-03-26 2003-10-02 Denso Corp 流体吸排装置
JP2006161723A (ja) 2004-12-08 2006-06-22 Denso Corp インペラおよびそれを用いた燃料ポンプ
US9599126B1 (en) * 2012-09-26 2017-03-21 Airtech Vacuum Inc. Noise abating impeller

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5001975B2 (ja) * 2008-08-29 2012-08-15 ヒョンダム インダストリアル カンパニー リミテッド 燃料ポンプ用不等ピッチインペラ

Also Published As

Publication number Publication date
WO2016089103A1 (ko) 2016-06-09
DE112015005494B4 (de) 2023-05-25
KR20160067402A (ko) 2016-06-14
US20170363091A1 (en) 2017-12-21
DE112015005494T5 (de) 2017-11-16
US10590938B2 (en) 2020-03-17

Similar Documents

Publication Publication Date Title
KR101671946B1 (ko) 부등피치 재생 블로워 및 이의 최적화 설계 방법
US9938833B2 (en) Methods, systems, and devices for designing and manufacturing flank millable components
US10352330B2 (en) Turbomachine part with a non-axisymmetric surface
KR101679189B1 (ko) 단일 유로 펌프 벌류트의 최적화 설계방법, 이에 의하여 설계된 단일 유로 펌프 벌류트 및 원심 펌프
KR101679188B1 (ko) 단일 유로 펌프 임펠러의 고효율 최적화 설계방법 및 이에 의하여 설계된 고효율 단일 유로 펌프 임펠러
KR101629005B1 (ko) 단일 유로 펌프 임펠러의 최적화 설계방법, 이에 의하여 설계된 단일 유로 펌프 임펠러 및 원심 펌프
Li et al. Multipoint and multiobjective optimization of a centrifugal compressor impeller based on genetic algorithm
Cattanei et al. Reduction of the tonal noise annoyance of axial flow fans by means of optimal blade spacing
Sugimura et al. Multi-objective optimization and design rule mining for an aerodynamically efficient and stable centrifugal impeller with a vaned diffuser
KR101881909B1 (ko) 고효율 저유체 유발 진동 단일채널펌프의 설계방법
US20180101162A1 (en) Methods, Systems, and Devices For Designing and Manufacturing Flank Millable Components
KR101742171B1 (ko) 고효율 엇회전 펌프 수차, 이의 최적화 설계방법 및 이를 포함하는 자가 발전 시스템
Benini et al. Design optimization of vaned diffusers for centrifugal compressors using genetic algorithms
KR102102190B1 (ko) 임펠러 재설계에 따라 출력 변경이 가능한 단일채널펌프의 설계안 도출방법
KR101162611B1 (ko) 축류압축기의 케이싱 그루브 최적설계방법
Schänzle et al. As good as it can be-ventilation system design by a combined scaling and discrete optimization method
RU2672990C1 (ru) Деталь газотурбинного двигателя с неосесимметричной поверхностью
CN109033661B (zh) 一种叶轮的设计方法
CN116595874A (zh) 叶轮机械性能预测模型参数优化方法及装置、存储介质
Sugimura et al. Kriging-model-based multi-objective robust optimization and trade-off-rule mining using association rule with aspiration vector
CA2954256C (fr) Procede de modelisation d&#39;une baignoire d&#39;une aube
Ha et al. An optimization method for centrifugal compressor design using the surrogate management framework
Šlachtič et al. CFD Analysis of Impeller Trimming
Schaffrath et al. Gradient-Free Aerodynamic Optimization With Structural Constraints and Surge Line Control for Radial Compressor Stage
KR20060056738A (ko) 원심다익송풍기의 요소 검출 방법 및 그 방법에 의해제조되는 원심다익송풍기

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
FPAY Annual fee payment

Payment date: 20190925

Year of fee payment: 4