Advanced
Diffraction and Radiation of Waves by Array of Multiple Buoys
Diffraction and Radiation of Waves by Array of Multiple Buoys
Journal of Ocean Engineering and Technology. 2016. Jun, 30(3): 151-160
Copyright © 2016, Korean Society of Ocean Engineers
  • Received : March 04, 2016
  • Accepted : June 24, 2016
  • Published : June 30, 2016
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
About the Authors
일형 조

Abstract
The diffraction and radiation of linear waves by an array of truncated floating multiple buoys are solved using the interaction theory based on a matched eigenfunction expansion method (MEEM). The interaction processes between multiple buoys are very complex and numerous, because the scattered and radiated waves from each buoy affect the others in the array. Our primary aim is therefore to construct the rigorous wave exciting forces and hydrodynamic forces to deal with the problem of multiple interactions. This present method is applied to a square array of four buoys with two incidence angles, and the results are given for the wave excitation forces on each buoy, heave RAO for each buoy heaving independently, and wave elevations around the buoys and wave run-up. The analytical solutions are in good agreement with the numerical solutions obtained from commercial code (WAMIT).
Keywords
1. 서 론
많은 종류의 해양구조물의 하부 구조는 다수 부체가 배열된 형태로 이루어져 있다. 반 잠수식 시추선, 다수 파력발전장치로 구성된 에너지 농장, 부교, 부유식 공항 등이 대표적인 예이다. 따라서 다수 부체로 구성된 하부 구조물의 유체역학적 특성을 정확히 파악하기 위해서는 다수 부체의 상호작용 문제를 풀어야 한다. 단일 부체와 달리 다수 부체 주변에서의 파계는 자체의 회절파와 방사파 뿐만 아니라 주변의 다른 부체로부터의 산란파와 방사파가 더해지는 다중 산란(Multiple scattering)으로 매우 복잡하다. 따라서 다수 부체간의 상호작용 문제는 수학적으로 매우 복잡할 뿐만 아니라 많은 계산시간을 요구한다.
단일 부체에 대한 회절문제는 Garrett(1971) 에 의해, 방사문제는 Tung(1979) Cho and Kweon(2011) 에 의해 해석되었다. 다수 부체가 배열된 경우 일어나는 회절문제는 Kagemoto and Yue(1986) 에 의해 처음으로 다루어졌다. Kagemoto and Yue(1986) 가 제안한 상호작용 이론(Interaction theory)은 단일 부체에 대한 해석해가 존재할 때 적용이 가능하다. Murai et al.(1999) 는 다수 부체로 지지되는 대형 부유 구조물에 Kagemoto and Yue(1986) 가 제안한 상호작용 이론을 적용하였다. Linton and Evans(1990) 는 해저면에 고정된 원기둥의 다수 배열에 대한 회절문제를 풀었다. Yilmaz and Incecik(1998) Kagemoto and Yue(1986) 의 상호작용 이론을 일정한 흘수를 갖는 원기둥의 다수 배열 문제에 적용하여 파기진력과 동유체력을 구하였다. Yilmaz et al.(2001) 는 4개의 원형 부체로 지지되는 반 잠수식 해양구조물에 대한 회절과 방사문제를 풀어 구조물 주변의 수면 변위를 구하고 실험결과와 비교하였다. Siddorn and Taylor (2008) Kagemoto and Yue(1986) 의 상호작용 이론에 기반을 둔 이론을 사용하여 다수 배열된 원형 부체에 작용하는 파기진력과 운동변위를 구하였다. Child and Venugopal(2010) 는 점 흡수식 파력발전장치가 여러 개 배열되어 있는 경우에 파와 다수 부체 사이의 상호작용 문제를 해석하여 파력발전장치의 최적 배열을 연구하였다.
본 연구에서는 다수 원형 부체에 대한 회절문제와 방사문제를 풀어 부체 간의 상호작용을 고려한 파기진력과 운동변위 그리고 다수 부체 주변의 수면변위를 구하였다. 이를 위하여 고유 함수전개법(Matched eigenfunction expansion method)에 기초를 둔 Child and Venugopal(2010) 의 해석방법을 사용하였다. 먼저 단일 원기둥에 대한 회절문제와 방사문제를 풀어 해석 해를 구하고 이를 산란 방정식과 결합하여 다수 원형 부체 간의 상호작용 효과를 고려한 해석이론을 정립하였다. 해석결과의 타당성을 검증하기 위하여 상용 수치코드인 WAMIT(Wave analysis MIT)의 수치해석 결과와 비교하였다.
2. 문제의 정식화
본 연구에서는 반경이 aj 이고, 흘수가 dj NB 개의 원형 부체( j =1,2,..., NB )가 일정한 수심( h ) 위에 배열되어 있을 때 다수 부체간의 상호작용을 고려한 각 부체에 작용하는 수직방향 파기진력과 수직 운동변위를 구하였다. Fig. 1(b) 와 같이 공간상에 고정된 좌표계 ( X , Y , z )와 NB 개의 다수 부체의 중심에 고정된 좌표계( rj , θj , z ), j =1,2,..., NB 를 도입한다. j 번째 부체의 중심의 좌표는 ( xj , yj )이다. i 번째 부체의 중심에 대한 j 번째 부체의 중심의 상대 위치는 ( Rij , αij )로 표현할 수 있다( Fig. 1(b) 참조). 선형 포텐셜 이론을 가정하여 속도포텐셜을 도입한다. 속도포텐셜을 시간과 공간의 함수로 분리하면 Φ ( r , θ , z , t ) = Re { ϕ ( r , θ , z ) e iwt }로 쓸 수 있다. 속도포텐셜의 공간함수 ϕ ( r , θ , z )는 아래와 같은 경계치 문제를 만족한다.
PPT Slide
Lager Image
(a) Side view of buoy and (b) plan view of multiple buoys and coordinate systems
PPT Slide
Lager Image
여기서 g 는 중력가속도이며, wj (=− iwξj )는 j 번째 부체의 수직 운동속도로 복소수 값을 갖는다. 속도포텐셜 ϕ ( r , θ , z )은 입사포텐셜( ϕI )과 산란포텐셜( ϕS ) 그리고 방사포텐셜( ϕR )의 합으로 표현된다. 이때 산란포텐셜과 방사포텐셜은 아래와 같은 방사경계조건식을 만족한다.
PPT Slide
Lager Image
여기서 k 1 는 입사파의 파수이다.
- 2.1 입사포텐셜
입사파가 진폭( A )과 주파수( w )를 가지고 x 축의 양의 방향과 β 의 입사각도를 이루며 진행할 때 j 번째 원형 부체의 중심에 고정된 좌표계( rj , θj , z )에서 입사포텐셜을 수식으로 표현하면 다음과 같다.
PPT Slide
Lager Image
여기서 Λj (= e ik1(xjcosβ+yjsinβ) )는 공간 고정 좌표계( OXY )의 원점과 j 번째 부체의 중심( xj , yj )간의 위상차를 나타내는 함수이다. 식 (3)에 나타난 입사포텐셜을 공간함수와 계수로 분리하여 벡터식으로 표현하면 다음과 같다.
PPT Slide
Lager Image
여기서 윗첨자 T 는 전치행렬을 나타낸다. 벡터식 aj
PPT Slide
Lager Image
의 원소들은 다음과 같다.
PPT Slide
Lager Image
PPT Slide
Lager Image
여기서 Jn , In 는 각각 제1종 Bessel함수와 변형된 Bessel함수를 나타낸다.
식 (6)에서 고유값( k 10 =− ik 1 , k 1m , m =1, 2,.. )은 선형분산식( k 1 tanh k 1 h = w 2 / g , k 1m tan k 1m h =− w 2 / g , m =1, 2,.. )을 만족하며, 고유함수 f 1m ( z )은 다음과 같이 정의된다.
PPT Slide
Lager Image
식 (7)에 정의된 고유함수 f 1m ( z )는 아래와 같은 직교성(Orthogonality)을 만족한다.
PPT Slide
Lager Image
여기서 δmn n = m 일 때 1이고, n m 일 때 0의 값을 갖는 Kronecker-Delta함수이다.
- 2.2 산란포텐셜
산란포텐셜은 고정된 부체에 의한 입사파의 산란문제를 풀어 구할 수 있다. 산란문제에 고유함수전개법을 적용하기 위하여 유체영역을 Fig. 1(a) 와 같이 영역 (I)과 영역 (II)로 나눈다. 따라서 영역 (I)은 rj aj , − h z ≤ 0 , 영역 (II)은 0 ≤ rj aj , − h z ≤ − dj 으로 정의된다.
영역 (I)의 산란포텐셜은 Laplace방정식과 함께 자유표면, 해저면, 방사 경계조건식을 만족한다. i 번째 부체 주변에서의 산란포텐셜은 다음과 같다.
PPT Slide
Lager Image
여기서 제2종 변형된 Bessel함수 Kn m = 0일 때
PPT Slide
Lager Image
이며, Hn 는 제1종 Hankel함수이다.
입사포텐셜과 마찬가지로 산란포텐셜을 공간함수와 계수로 분리하여 벡터식으로 표현하면 다음과 같다.
PPT Slide
Lager Image
여기서 공간함수
PPT Slide
Lager Image
의 원소들은 다음과 같다.
PPT Slide
Lager Image
NB 개의 다수 부체들 간의 상호작용을 고려하기 위하여 i 번째 부체에 의해 산란된 파는 나머지 부체( J =1, 2,..., NB , i j )에게는 입사파로 작용한다. 산란포텐셜
PPT Slide
Lager Image
을 입사포텐셜
PPT Slide
Lager Image
로 표현하기 위하여 아래 주어진 Graf의 Bessel함수의 덧셈정리(Additional theorem)를 사용한다( Linton and Evans, 1990 ; Maniar and Newman, 1997 ; Kim, 1993 ).
PPT Slide
Lager Image
식 (12)을 이용하여 i j 인 모든 부체에 대하여
PPT Slide
Lager Image
PPT Slide
Lager Image
로 표현하면 다음과 같다.
PPT Slide
Lager Image
여기서 전달함수 Tij 의 원소는 다음과 같다.
PPT Slide
Lager Image
식 (13)을 식 (10)에 대입하면 i j 일 때 i 번째 부체에 의한 산란된 파는 j 번째 부체에게 입사파로 작용함을 나타내는 산란 포텐셜을 아래와 같이 구할 수 있다.
PPT Slide
Lager Image
i 번째 부체의 영역 (II)에서의 속도포텐셜은 변수분리법을 사용하여 다음과 같이 구할 수 있다.
PPT Slide
Lager Image
식 (15)의 공간함수
PPT Slide
Lager Image
의 원소는 다음과 같다.
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
i 번째 부체의 영역 (II)에서의 고유값이다.
- 2.3 방사포텐셜
방사포텐셜은 입사파가 없는 상태에서 부체의 수직운동으로 생성된 파형을 나타내는 속도포텐셜이다. i 번째 부체의 수직운동에 의한 영역 (I)의 방사포텐셜은 다음과 같다.
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
는 무차원화된 수직운동의 복소진폭이다.
산란포텐셜과 마찬가지로 i j 일 때 i 번째 부체의 수직운동으로 생성된 방사파는 나머지 부체( j = 1, 2,..., NB , i j )에 입사파로 작용한다. 2.2절과 동일한 방법으로 Graf의 덧셈정리 식 (12)을 이용하여 식 (17)을 다시 쓰면 다음과 같다.
PPT Slide
Lager Image
같은 방법으로 i 번째 부체의 영역 (II)에서의 방사포텐셜은 다음과 같다.
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
이며,
PPT Slide
Lager Image
는 영역 (II)의 특별해이다.
- 2.4 산란 방정식
j 번째 부체에 입사파로 작용하는 속도포텐셜은 먼 외해로부터 들어오는 순수한 입사포텐셜과 자신을 제외한 나머지 부체들에 의한 산란포텐셜과 방사포텐셜의 합으로 표현할 수 있다.
PPT Slide
Lager Image
식 (20)의 입사파의 계수
PPT Slide
Lager Image
와 식 (10)의 산란파의 계수 Ai 사이에는 회절 전달 행렬(Diffraction transfer matrix) Bj 에 의해 아래와 같은 관계식이 성립한다.
PPT Slide
Lager Image
여기서 회절 전달 행렬 Bj Ri 는 단일 원기둥에 대한 회절문제와 방사문제를 풀어 구할 수 있다(부록 A, B 참조).
영역 (II)에 대해서도 같은 방법이 적용된다.
PPT Slide
Lager Image
- 2.5 운동 방정식
j 번째 부체의 영역 (II)에서의 속도포텐셜은 회절포텐셜과 방사포텐셜의 합으로 아래와 같이 표현된다.
PPT Slide
Lager Image
j 번째 부체의 수직방향 유체력
PPT Slide
Lager Image
을 구하기 위하여 동압력을 부체 바닥면에 따라 적분한다.
PPT Slide
Lager Image
j 번째 부체에 대한 수직방향 운동방정식을 쓰면 다음과 같다.
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
j 번째 부체의 질량이며,
PPT Slide
Lager Image
j 번째 부체의 점성 감쇠계수로 자유감쇠실험(Free-decay test)을 통하여 실험적으로 구할 수 있다. 본 연구에서는 무차원화된 감쇠계수로 κj = 0.02을 사용하였다. 여기서 ωoj 는 비감쇠 공진주파수로
PPT Slide
Lager Image
이며, 부체의 질량이 부가질량보다 큰 경우( mj µj ),
PPT Slide
Lager Image
이다.
식 (25)의 오른쪽 항의 각 부체의 수직방향 유체력은 다음과 같다.
PPT Slide
Lager Image
여기서
PPT Slide
Lager Image
식 (26)에서 입사파에 의한 j 번째 부체에 작용하는 파기진력을 추출하면 다음과 같다.
PPT Slide
Lager Image
식 (26)을 식 (25)에 대입하여 정리하면
PPT Slide
Lager Image
이다. 여기서
PPT Slide
Lager Image
행렬식 (21)과 (28)을 풀어 미지수
PPT Slide
Lager Image
, Ai ( i = 1,2,.., NB )를 구하여 영역 (I)에서 완전한 해를 구하고 영역 (II)에서의 완전한 해는 식 (22)을 이용한다.
부체가 고정되어 있을 때 다수 부체 배열에 의한 각 부체의 주변의 회절파의 파형식은 식 (4)와 식 (9)의 합으로 다음과 같다.
PPT Slide
Lager Image
윗식으로부터 j 번째 부체 벽면( rj = aj )에서의 파형식인 처올림(Run-up) 파는 아래와 같다.
PPT Slide
Lager Image
3. 계산결과 및 고찰
다수 원형 부체 배열에 대한 해석해를 검증하기 위하여 반경 a = 1 m 인 4개의 동일한 부체가 정사각형으로 배열되어 있는 계산모델에 대하여 해석해를 구하고 WAMIT을 이용한 수치계산 결과와 비교하였다. 4개의 원형 부체 간의 간격과 수심은 4m이다. 계산에 사용한 z 축 방향의 고유함수 개수는 20개이며, θ 축 방향의 고유함수의 개수는 11개이다. 한편 WAMIT에서 격자수는 각 부체 당 280개이다.
먼저 해석해를 검증하기 위하여 단일 원형 부체( d / a = 2.0, h / a = 4.0)에 대하여 WAMIT의 수치계산 결과와 Fig. 3 에서 비교하였다. 여기서 x 축은 k 1 a 로 무차원화된 파장이며, y 축은 무차원화된 파기진력, 부가질량, 방사감쇠계수이다. 실선은 고유함수전개법을 이용한 해석결과이며, 동그라미는 WAMIT을 이용한 수치계산 결과이다. 고유함수 개수가 비교적 적음에도 불구하고 무차원화된 부가질량과 방사감쇠계수 그리고 파기진력 모두 서로 잘 일치하고 있음을 보여주고 있다.
PPT Slide
Lager Image
Non-dimensional wave exciting force, added mass and radiation damping coefficient for single buoy with d/a = 2.0, h/a = 4.0 (line : present solutions, symbol: WAMIT numerical solution)
Fig. 4 는 4개의 원형 부체가 Fig. 2 와 같이 배열되어 있을 때 각 부체에 작용하는 파기진력을 WAMIT의 계산결과와 비교한 그림이다. 여기서 선은 해석해의 결과이며, 기호는 WAMIT의 결과이다. 입사파의 각도( β )는 0°, 45°이다. 단일 원기둥에서와 같이 해석결과는 수치해석 결과를 잘 따라가고 있음을 보여준다. 특히, 입사각도가 0도일 때는 1, 2번과 3, 4번 부체는 입사파의 진행방향에 서로 대칭으로 놓여 있으므로 파기진력은 서로 같음을 볼 수 있다. 그러나 입사각도가 45도일 때는 2번과 3번 원형 부체만이 입사파의 진행방향에 대해 서로 대칭을 이뤄 동일한 파기진력을 주지만 파의 전면에 놓인 1번 부체에서 파기진력이 가장 크고, 후면에 놓인 4번 부체에서 상대적으로 가장 작은 파기진력이 작용한다.
PPT Slide
Lager Image
Configuration Sketch for the square array of four buoys
PPT Slide
Lager Image
Non-dimensional wave exciting forces on each buoy in four-buoy array with d/a = 2.0, s/a = 4.0, h/a = 4.0 (line : present solutions, symbol: WAMIT numerical solution)
Fig. 5 Fig. 4 와 동일한 계산 조건하에서 4개의 부체에 작용하는 파기진력들을 모두 합한 결과
PPT Slide
Lager Image
를 보여 주고 있다. 즉, 4개의 부체로 지지되는 반 잠수식 시추선에 작용하는 파기진력이라 할 수 있다. k 1 a 가 0.5보다 작은 장파 영역에서는 입사각도 β = 0°, 45°에 따른 파기진력의 차이는 거의 나타나지 않았지만 k 1 a 가 0.7보다 큰 단파 영역에서는 입사각도 β = 0°일 때의 파기진력이 β = 45° 보다 다소 커지는 것을 볼 수 있다.
PPT Slide
Lager Image
Non-dimensional wave exciting forces on a four-buoy array with d/a = 2.0, s/a = 4.0, h/a = 4.0 (line : present solutions, symbol: WAMIT numerical solution)
Fig. 6 에서는 각 부체의 수직운동 RAO(Response amplitude operator)를 살펴보았다. 4개의 원형 부체 사이의 상호작용을 고려하였을 때 전면에 놓인 부체의 수직운동 변위는 단일 부체가 있을 때와 비슷하거나 오히려 증가한다. 그러나 후면에 놓인 원형 부체의 수직운동 변위는 파의 차단 효과로 인해 다소 줄어드는 결과를 주었다. 수직운동의 RAO(| ξj / A |) 곡선의 최고점을 나타내는 k 1 a 는 부체의 공진주파수( k 1 a = 0.42)와 일치하고 있다. 공진주파수 주변에서 각 부체 간의 운동변위는 다소 차이를 보이나 나머지 주파수영역에서는 큰 차이가 나타나지 않았다.
PPT Slide
Lager Image
Heave RAO of each buoy in four-buoy array with d/a = 2.0, s/a = 4.0, h/a = 4.0
Fig. 7 은 식 (26)에 주어진 j 번째 부체에 작용하는 수직방향 유체력을 보여주고 있다. j 번째 부체에 작용하는 유체력은 파기진력과 j 번째 부체의 운동과 j 번째 부체를 제외한 나머지 부체들의 운동에 의한 동유체력의 합으로 구성되어 있다. 따라서 저주파수와 고주파수 영역에서는 수직방향 유체력은 파기진력의 거동 특성을 따라가지만 공진주파수 부근에서는 부체의 수직운동 특성을 따라감을 보여준다.
PPT Slide
Lager Image
Non-dimensional fluid forces on each buoy in four-buoy array with d/a = 2.0, s/a = 4.0, h/a = 4.0
Fig. 8 은 부체가 고정되어 있을 때 처올림 진폭비(| ηD ( aj , θj )/ A )를 보여주고 있다. 이때 무차원화된 입사파의 파장은 k 1 a =1.0이다. x 축은 원형 부체의 원주 각도이다. (a)는 입사각도 β = 0°에서의 처올림 진폭비이며 (b)는 β = 45°에서의 처올림 진폭비이다. 입사각도 β = 0°일 때는 파의 전면(90° < θ < 270°)에서의 처올림 진폭비가 파의 후면(0° < θ < 90°, 270° < θ < 360°) 보다 더 큰 것을 알 수 있다. x 축에 대한 부체의 대칭 배열로 특정한 원주 각도 범위내에서 동일한 처올림 진폭비를 보여주고 있다. 예를 들어 1번 부체의 0° < θ < 90°의 처올림 진폭비는 2번 부체의 270° < θ < 360°의 처올림 진폭비와 서로 같다. 입사각도가 β = 45°일 때 1번 원형 부체의 처올림 진폭비가 파의 전면 (180° < θ < 270°)과 후면(0° < θ < 90°) 모두에서 나머지 부체들의 처올림 진폭비 보다 크게 나타났다.
PPT Slide
Lager Image
Wave run-up on wall of each buoy of four-buoy array with d/a = 2.0, s/a = 4.0, h/a = 4.0, k1a = 1.0
Fig. 9 는 파의 입사각도( β )는 0°, 45°일 때 특정한 입사파의 파장( k 1 a = 1.0)에서 사각형으로 배열된 4개의 고정된 부체 주변에서의 회절파(| ηD |/ A )의 분포를 보여주고 있다.
PPT Slide
Lager Image
Diffracted wave eleva (|ηD(x,y)|/A) around four-buoy array with d/a = 2.0, s/a = 4.0, h/a = 4.0, k1a = 1.0
Acknowledgements
본 연구는 해양수산부의 해양청정에너지개발사업 "10MW급 부유식 파력-해상풍력 연계형 발전시스템 설계기술 개발"과제의 지원으로 수행되었으며, 연구비 지원에 감사드립니다.
References
Child B.F.M. , Venugopal V. 2010 Optimal Configuration of Wave Energy Device Arrays Ocean Engineering 37 (6) 1402 - 1417    DOI : 10.1016/j.oceaneng.2010.06.010
Cho I.H. , Kweon H.M. 2011 Extraction of Wave Energy Using the Coupled Heaving Motion of a Circular Cylinder and Linear Electric Generator Journal of Ocean Engineering and Technology 25 (6) 9 - 16
Garrett C.J.R. 1971 Wave Forces on a Circular Dock Journal of Fluid Mechanics 46 129 - 139    DOI : 10.1017/S0022112071000430
Kagemoto H. , Yue D.K.P. 1986 Interactions among Multiple Three Dimensional Bodies in Water Waves: an Exact Algebraic Method Journal of Fluid Mechanics 166 189 - 209    DOI : 10.1017/S0022112086000101
Kim M.H. 1993 Interaction of Waves with N-vertical Circular-cylinders Journal of Waterway Port Coastal and Ocean Engineering 119 (6) 671 - 689    DOI : 10.1061/(ASCE)0733-950X(1993)119:6(671)
Linton C.M. , Evans D.V. 1990 The Interaction of Waves with Arrays of Vertical Circular Cylinders Journal of Fluid Mechanics 215 549 - 569    DOI : 10.1017/S0022112090002750
Maniar H.D. , Newman J.N. 1997 Wave Diffraction by a Long Array of Cylinders Journal of Fluid Mechanics 339 309 - 330    DOI : 10.1017/S0022112097005296
Murai M. , Kagemoto H. , Fujino M. 1999 On the Hydroelastic Responses of a Very Large Floating Structure in Waves Journal of Marine Science and Technology 4 123 - 153    DOI : 10.1007/s007730050015
Siddorn P. , Taylor R.E. 2008 Diffraction and Independent Radiation by an Array of Floating Cylinders Ocean Engineering 35 (13) 1289 - 1303    DOI : 10.1016/j.oceaneng.2008.06.003
Tung C.C. 1979 Hydrodynamic Forces on Submerged Vertical Circular Cylindrical Tanks under Ground Excitation Applied Ocean Research 1 (2) 75 - 78    DOI : 10.1016/0141-1187(79)90020-8
Yilmaz O. , Incecik A. 1998 Analytical Solutions of the Diffraction Problem of a Group of Truncated Vertical Cylinders Ocean Engineering 25 (6) 385 - 394    DOI : 10.1016/S0029-8018(97)00003-6
Yilmaz O. , Incecik A. , Barltrop N. 2001 Wave Enhancement Due to Blockage in Semi-submersible and TLP Structures Ocean Engineering 28 (5) 471 - 490    DOI : 10.1016/S0029-8018(00)00023-8