Warning: mkdir(): Permission denied in /home/virtual/lib/view_data.php on line 81

Warning: fopen(upload/ip_log/ip_log_2024-03.txt): failed to open stream: No such file or directory in /home/virtual/lib/view_data.php on line 83

Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 84
Comparison of the Formulas for the Wave Forces Acting on the Perforated Caisson Breakwater

Comparison of the Formulas for the Wave Forces Acting on the Perforated Caisson Breakwater

유공케이슨 방파제에 작용하는 파력 공식의 비교

Article information

J Korean Soc Coast Ocean Eng. 2015;27(4):217-227
* Chang-Hwan Ji, Coastal & Environmental Engineering Division, Korea Institute of Ocean Science & Technology, Ansan 15627, Korea, Tel: +82-31-400-7822, Fax:+82-31-408-5823, chji@kiost.ac.kr
Received 2015 July 06; Revised 2015 August 13; Accepted 2015 August 18.

Abstract

In this study, two-dimensional physical experiment was carried out to examine the applicability of the three formulas(Takahashi and Shimosako, 1994; Tabet-Aoul and Lambert, 2003; Li, 2007), which were proposed to calculate the wave forces acting on perforated caisson breakwaters. In order to quantitatively compare the measured with the estimated values based on the wave formulas, the refined index of agreement and the coefficient of determination were calculated, by which the degree of agreement was evaluated. Among the three wave formulas, DUT formula (Li, 2007) showed the smallest deviation from the measured forces, whereas Takahashi formula (Takahashi and Shimosako, 1994) showed the largest deviation. Meanwhile, comparison of the magnitude of the measured wave forces with those from the three formulas revealed that DUT formula slightly underestimate, while the others overestimate the measured forces.

Trans Abstract

본 연구에서는 유공케이슨에 작용하는 파력을 산정하기 위해 제안된 세 가지 파력공식(Takahashi and Shimosako, 1994; Tabet-Aoul and Lambert, 2003; Li, 2007)의 적용성을 검토하기 위하여 2차원 수리실험을 수행하였다. 실험에 의한 계측값과 각 공식에 의한 계산값의 정량적 비교를 위하여 수정 일치지수 및 결정계수를 이용하였고, 이를 통해 계측값과 계산값이 어느 정도 서로 일치하는지 그 정도를 평가하였다. 세 가지 공식 중 DUT 공식(Li, 2007)이 실험값과 계산값 편차가 가장 적었으며, Takahashi 공식(Takahashi and Shimosako, 1994)은 가장 큰 편차를 나타내었다. 한편, 계측 및 계산값의 상대적인 크기를 비교한 결과 Takahashi공식과 Tabet-Aoul & Lambert 공식(Tabet-Aoul and Lambert, 2003)은 실험값에 비하여 파력을 과대산정하고, DUT 공식은 파력을 다소 과소산정하는 경향이 나타났다.

1. 서 론

우리나라 주요 항만이 대형화 되면서 최근 건설되는 방파제와 부두에 케이슨이 널리 사용되고 있다. 특히 반사파 및 파력 저감 등의 이유로 무공케이슨 보다는 유공케이슨이 주로 사용되는 실정이다. Jarlan (1961)이 해안구조물에 유공벽 개념을 도입한 이후, 유공케이슨에 대한 많은 연구들이 있었다. 하지만 이들 대부분은 유공방파제에 의한 반사파 저감에 관한 것이었으며(Fugazza and Natele, 1992; Kakuno et al., 1992; Park et al., 1993 등), 유공케이슨에 작용하는 파력에 대한 연구는 상대적으로 많지 않은 편이다.

Takahashi and Shimosako (1994)는 직립벽에 작용하는 파력에 대한 Goda의 연구(Goda, 1974)를 기초로 유공케이슨에 대한 2차원 수리모형실험을 수행하여 유공케이슨에 적용할 수 있는 파압공식과 파력보정계수를 제안하였으며, 이 결과는 현재 가장 보편적으로 사용되고 있다. 이후 Tabet-Aoul and Lambert (2003)는 수리실험 및 현장관측자료를 바탕으로 위상보정계수를 도입하고, Takahashi 공식에 기초하여 유공케이슨에 작용하는 최대 수평파력을 구하는 공식을 제안하였다. 이밖에도 2000년 이후로 중국에서는 Dalian University of Technology 연구팀에 의해 유공케이슨에 대한 다양한 연구가 진행되고 있으며(Li, 2007), 이들이 제안한 유공케이슨에 작용하는 최대 수평파력 및 연직파력을 계산하는 방법은 중국 방파제 설계 기준에 수록되어 있다(Huang et al., 2011).

본 연구에서는 유공케이슨 방파제에 작용하는 파력 계산을 위한 위 공식들의 형태 및 특징을 조사하였다. 또한 2차원 수리모형실험을 수행하여 유공케이슨 각 벽에 작용하는 파력을 독립적으로 계측하고 그 합력으로부터 케이슨 전체에 작용하는 파력을 산정하여 그 결과를 기존공식에 의한 예측값과 비교하여 각 공식들의 특징을 살펴보았다.

2. 유공케이슨 파력 공식

2.1 Takahashi 공식

Takahashi 공식은 무공직립벽의 파력을 계산하는 Goda 공식을 유공케이슨 방파제에 적용할 수 있도록 수정한 공식이다(Takahashi and Shimosako, 1994). Takahashi 공식에서 p1, p2, p3, p4, pu의 위치는 Goda 공식과 동일하며, 각각은 다음의 식을 이용하여 계산할 수 있다.

(1) p1=0.5(1+cosθ)λ1α1+λ2α*cos2γHd
(2) pu=0.5(1+cosθ)λ3α1α3γHd
(3) η0=0.75(1+cosθ)λ1Hd
(4) α*=maxα2,α1

보정계수 λ1는 입사파의 중복파 성분에 의한 파력변화, λ2는 입사파의 충격파 성분에 의한 파력변화, λ3는 양압력의 변화를 나타내는 계수이다. 정수면에서의 파압 p1을 계산할 때 Goda 공식에서의 α2 대신 α*를 사용하며, 충격파의 영향까지 고려한다. 식 (4)에서 충격파의 영향을 고려하는 인자인 αI를 계산하는 방법은 Takahashi et al. (1991)에 제시되어 있다.

Takahashi and Shimosako (1994)에 의하면 무공직립벽과는 달리 유공케이슨에서는 수평파력 및 연직파력의 최대값이 파봉이 케이슨의 전면에 위치하는 순간에 발생하지 않았다. 따라서 파의 위상에 따라 6개의 위상(Crest-I, Crest-IIa, Crest-IIb, Trough-I, Trough-II, Trough-III)을 정의하였다. 즉, 전면 유공벽에 작용하는 파력이 최대일 때 Crest-I, 후면벽에 작용하는 충격파 성분이 최대일 때 Crest-IIa, 후면벽에 작용하는 중복파 성분이 최대일 때는 Crest-IIb로 정의하였으며, 최대 수평파력 및 연직파력을 구하기 위해서는 이 세 위상을 모두 고려하여야 한다. 전면벽에 작용하는 파력이 최소일 때 Trough-I, 케이슨 전면 및 유수실 내부의 수위가 최소일 때 각각 Trough-II와 Trough-III로 정의 하였다.

Takahashi 공식의 특징은 유공케이슨 각 벽체에 각각 다른 값의 보정계수 λ1, λ2, λ3을 적용하여 파력을 계산한다는 점이다. 이 값들은 파의 위상에 따라서도 달라진다. Takahashi and Shimosako (1994)는 각 벽체에 적용할 보정계수를 Table 1에 보인 것처럼 위상별로 나타내었다. 여기서 L'은 수심 d1에서의 파장이며, B는 유수실의 폭이다. 아랫첨자 S, L, R, M, U는 유공케이슨의 각 부분을 나타내며, 그 위치는 Fig. 1에서 확인할 수 있다.

Correction coefficients for perforated caissons given by Takahashi

Fig. 1.

Pressure distributions on the perforated caisson at Crest-IIb by Takahashi formula.

2.2 Tabet-Aoul & Lambert 공식

Tabet-Aoul and Lambert (2003)는 실험자료에 근거하여 유공케이슨에 작용하는 최대 수평파력은 Takahashi and Shimosako (1994)에서 정의한 세 가지의 Crest 위상 사이에서 발생하므로, 유수실 폭과 관계가 있는 위상조정계수를 도입하여 다음과 같이 최대 수평파력 공식을 개발하였다. Fig. 2에는 Tabaet-Aoul & Lambert 공식에 의한 파압분포도를 나타내었다.

Fig. 2.

Pressure distributions on the perforated caisson by Tabet-Aoul & Lambert formula.

(5) pp1=(1+cosθ)0.21α1+B4L1+α*cos2θγHd
(6) pp3=α3pp1
(7) pp4=α4pp1
(8) pr1= 0.51+cosθ0.7-BL2α1+0.43-BL1+α*cos2θγHd
(9) pr3=α3pr1
(10) pr4=α4pr1
(11) ηp0=0.32(1+cosθ)Hd
(12) ηr0=0.751+cosθ0.7-BL2Hd

여기서, pp1, pp3, pp4는 각각 전면 유공벽의 정수면, 하단, 상단에서의 파압을, pp1, pp3, pp4는 각각 후면 무공벽의 정수면, 하단, 상단에서의 파압을 나타낸다. ηp0ηr0는 각각 유공벽과 후면벽에서의 처오름 높이이며, 그 외 다른 값들은 Takahashi 공식과 같은 값을 이용한다.

Tabaet-Aoul & Lambert 공식에서 각 벽에 작용하는 파압 분포는 Fig. 2에서와 같이 사다리꼴 분포이므로, 이를 적분하면 다음과 같은 파력공식을 얻을 수 있다.

(13) Fp=pp1+pp3d22+pp1+pp4hc*2(1-ε)
(14) Fr=pr1+pr3d22+pr1+pr4hc*2
(15) Ftot=χFp+Fr
(16) χ=1- 925BL+ 114BL2-4BL3+103BL4

이 식에서 Fp, Fr, Ftot는 각각 전면 유공벽, 후면 무공벽 및 케이슨 전체에 작용하는 최대수평파력이며, ε은 유공율, χ는 케이슨 전체에 작용하는 최대 수평파력에 대한 각 벽에 작용하는 최대 수평파력의 합의 비율로 정의되는 위상조정계수이다.

Tabaet-Aoul & Lambert 공식은 파압 계산식과 위상조정계수에 파라메터 B/L을 포함시켜 유수실 폭에 따른 최대 수평파력의 증감효과를 고려하였지만, 케이슨의 단면 형상이 일반적으로 현장에서 많이 시공되는 유수실 하부에 속채움이 있는 형태와 달라서 공식을 적용하는데 제약이 있을 수 있다.

2.3 DUT 공식

Dalian University of Technology (DUT) 연구팀에서는 2000년 이후로 유공케이슨 방파제에 대한 다양한 연구를 수행하였으며, 현재 중국의 방파제 설계기준으로 사용되고 있는 새로운 형태의 공식을 제안하였다.

(17) P1P0=0.997-1.515HsLs-0.804dLs-1.312BLs+0.25ε
(18) P2 P1=1.247+0.648dLs-0.573BLs-0.349ScLs+0.082ScLs2+0.215ε
(19) Pv2 Pv1=2.134+7.056HsLs-2.084dLs+8.273BLs-19.508BLs2+0.832ε
(20) Pv2 Pv1= 2.134-7.056HsLs+0.961dLs-0.296BLs+4.3656BLs2-0.811ScLs+0.184ScLs2+0.703ε
(21) l1l2=1.063-2.56HsLs+1.019dLs-0.88dLs2-1.432BLs+2.848BLs2+0.091ε
(23) lv1 lv0=1.01-1.917HsLs+1.063dLs-2.023BLs+2.532BLs2+0.429ε
(24) lv2 lv1=1.0
(25) Δt1 Ts=0.009+0.477BLs+0.099dLs+0.324 ε
(26) Δt2 T2=-0.237+0.304BLs+0.08dLs+0.299ScLs-0.088ScLs2+0.347 ε

여기서 HS는 입사파의 유의파고, LS는 유의파 주기의 파장, d는 수심, Sc는 정수면에서 유수실 상부 덮개까지의 거리이다. P0Pv1는 무공케이슨에 작용하는 1% 누적확률 수평 및 연직파력, P1Pv1은 유수실 상부 덮개가 없는 유공케이슨에 작용하는 1% 누적확률 수평 및 연직파력, P2Pv2는 유수실 상부 덮개가 있는 유공케이슨에 작용하는 1% 누적확률 수평 및 연직파력이다. l은 케이슨 전면 바닥에서 최대 수평파력 또는 최대 연직파력이 작용하는 지점까지의 거리, Δt는 최대 수평파력과 최대 연직파력이 작용 시점의 위상차이다.

이 공식은 파형경사(HS/LS), 파장대 수심비(d/LS), 파장대 유수실 폭 비(B/LS) 등의 영향을 고려했다는 점이 Takahashi 공식과 구별되는 특징이지만 개별 벽에 작용하는 파압분포나 파력은 계산할 수 없다.

3. 수리모형실험

3.1 실험시설 및 모형

수리모형실험은 한국해양과학기술원(KIOST)이 보유한 2차원 조파수조(길이 53 m, 높이 1.25 m, 폭 1 m)에서 수행되었다. Fig. 3에 보인 것처럼 수조의 길이방향으로 분할판을 설치하여 수조를 폭 0.6 m의 넓은수로와 0.4 m의 좁은수로로 나누고, 넓은수로에는 방파제 모형을 설치하였으며, 좁은수로에서는 통과입사파를 계측하였다.

Fig. 3.

Experimental setup (unit: mm).

케이슨 모형은 아크릴판을 이용하여, 케이슨 각 부분의 파력을 독립적으로 계측할 수 있도록 분리 제작하였다. 케이슨 모형의 높이는 47.5 cm, 폭은 70 cm이며, 상부구조물의 높이는 6 cm이다. 모형은 높이 17.5 cm의 사석 마운드 위에 설치되었으며, 방파제 전후면 마운드 수평 구간의 길이는 26 cm로 동일하였다. 유수실의 폭은 14 cm와 28 cm 두 가지로 달랐으며, 각각의 경우에 대하여 전면벽 유공부의 유공률을 20% 및 30%의 두 가지로 다르게 적용하였다. 케이슨 모형의 측면도 및 정면도를 Fig. 4에 나타내었다. Fig. 4(b)에 보인 유수실 폭이 넓은 단면형상은 Takahashi 공식이 개발된 근거가 되는 연구인 Takahashi and Shimosako (1994)의 실험에서의 형상과 거의 유사하며, 상치 구조물 부분의 높이만 다소 낮다. 한편, Fig. 4(a)Fig. 4(b)와 비교했을 때 유수실 폭만 절반이며, 나머지 제원은 모두 동일하다. 이처럼 본 연구에서 Takahashi and Shimosako (1994)의 단면 형상을 참조하여 실험 모형의 제원을 결정한 것은 2장에서 설명한 세 가지 파압 공식 중에서 Takahashi 공식이 현업 설계 실무에서 가장 폭넓게 활용되고 있기 때문에 이 공식이 개발된 기초가 되는 실험과 유사한 실험을 본 연구에서 수행함으로써 이 공식의 속성을 더 잘 이해하고 그 적용성을 검토하기 위해서였다(Oh et al., 2014).

Fig. 4.

(a) Side view of the narrow chamber caisson, (b) Side view of the wide chamber caisson, (c) Front view of the caisson with different front wall porosity (unit: mm).

3.2 계측기기의 배치 및 실험 방법

본 실험에서는 총 14개의 파고계를 사용하였으며, Fig. 3에 보인 것처럼 수조 내에 배치하였다. 또한, 파 작용 시 케이슨 벽에 작용하는 힘을 직접 계측하기 위해서 Fig. 5에 보인 것처럼 1축 하중계(load cell)를 이용한 파력 계측 시스템을 구성하여 전면벽의 무공부 및 유공부와 후면벽, 유수실 바닥판 및 상치구조물에 작용하는 파력을 독립적으로 각각 측정하였다. 그리고, 케이슨 바닥면에는 파압계 3개를 설치하여 케이슨 전체 구조물에 작용하는 연직파력을 계측하였다. 이와는 별도로 케이슨을 구성하는 각 벽체에 파압계를 부착하여 파력 뿐만 아니라 파압분포도 계측하였으나 본 연구에서는 1축 하중계를 이용하여 계측된 자료를 분석한 결과만을 제시하였다. 한편, 본 연구에 앞서 수행된 선행 연구(Oh et al., 2013)를 통해서 전면벽이 무공벽인 경우의 계측 파력은 Goda의 파압공식에 의한 산정값에 상응함을 확인하였다.

Fig. 5.

Installation of the model and arrangement of the load cells.

실험 수심은 51, 55, 59 cm의 세 가지로 다르게 하였으며, 주기 1.7, 2.1 s 일 때, 파고 11.0, 13.8, 16.5, 19.3, 22.0 cm, 주기 2.7 s 일 때, 파고 13.8, 16.5, 19.3, 22.0, 24.8 cm의 총 15개의 규칙파 조건에 대한 실험자료를 이용하여 분석을 진행하였다. Fig. 4에 보인 것처럼 유수실 폭 및 전면벽 유공률에 따라 서로 다른 4개의 모형 단면을 실험에 이용하였으므로 총 180개의 서로 다른 모형 및 실험파 조건에 대해서 분석이 이루어졌다. 각 실험파의 조파 시간은 60초였으며, 모든 파고 및 파력의 계측자료는 800 Hz로 PC의 저장장치에 기록되었다.

4. 분석 결과

어떠한 예측식에 의한 추정값과 실험에 의한 계측값이 얼마나 잘 일치하는지를 객관적인 수치로 보여주는 방법에는 여러가지가 있으며, 그 중 본 연구에서는 결정계수(coefficient of determination, R2)와 Willmott et al.(2012)에 의해 제안된 수정 일치지수(refined index of agreement, dr)를 이용하였다.

(27) R2=1-(Oi-fi)2(Oi-O¯)2
(28) dr=1-i=1nPi-Oi2i=1nOi-O¯ when i=1nPi-Oi2i=1nOi-O¯2i=1nOi-O¯i=1nPi-Oi-1  when i=1nPi-Oi>2i=1nOi-O¯

여기서, Oi는 실험에 의한 계측값, O¯ 는 계측값의 평균이며, fi는 선형회귀식에 의한 계산값, Pi는 예측식에 의한 계산값이다. 모든 관측값이 추정된 회귀식 상에만 있으면, R2 = 1.0인 반면, R2 = 0에 가까울수록 추정된 회귀식이 변수 사이의 관계를 전혀 설명해 주지 못함을 의미한다. 또한 관측값과 예측값이 완전히 일치할 때, dr = 1.0이며, 계산값이 관측값을 잘 예측하지 못할수록 dr = -1.0에 가까워진다.

4.1 Takahashi 공식과의 비교

Figs. 6-9에는 유공케이슨 방파제의 전면벽, 후면벽에 작용하는 최대 수평파력과 케이슨 전체에 작용하는 최대 수평파력 및 최대 연직파력을 Takahashi 공식을 이용한 계산값과 비교하였다. 여기서 아랫첨자 F, R, H, V는 각각 전면벽, 후면벽, 케이슨 전체에 작용하는 수평력 및 케이슨 전체에 작용하는 수직력을, 윗첨자 E와 TK는 실험에 의한 계측값과 Takahashi 공식에 의한 계산값을 나타낸다. Tables 2와 3에서도 Takahashi 공식에 의한 결과는 TK 기호로 나타내었으며, 나머지 두 가지 공식의 경우에도 각 공식을 의미하는 기호를 사용하였다.

Fig. 6.

Comparison of the maximum horizontal forces between the results by the experiments and Takahashi formula for the wave force on the front wall: (a)B=14 cm and (b)B=28 cm.

Fig. 9.

Comparison of the maximum vertical forces between the results by the experiments and Takahashi formula for the wave force on the whole caisson: (a) B=14 cm and (b)B=28 cm.

Fig. 6에는 유공케이슨 전면벽 유수실 폭이 좁은 경우와 넓은 경우에 대하여 케이슨 전면벽에 작용하는 파력을 Takahashi 공식으로 산정한 값과 실험에 의해 계측된 최대값을 비교하여 각각 나타내었다. Takahashi 공식에 의한 예측값은 실험값과 대체로 일치하는 경향을 보이긴 하였으나, 예측값이 더 큰 결과를 나타내었다. 유수실 폭에 따른 차이를 살펴보면, 유수실 폭이 넓은 경우의 적합도가 더 좋음을 확인할 수 있으며, Table 2에 제시된 수정 일치지수 및 결정계수의 값도 유수실 폭이 넓은 경우가 더 높은 수치를 보였다.

Refined index of agreement(dr) and coefficient of determination(R2)

후면벽에 작용하는 파력의 경우, 유수실 폭에 따른 차이가 더욱 분명하게 나타났다. 유수실 폭이 좁은 경우, Fig. 7(a)에 보인 것처럼 Takahashi 공식에 의한 계산값이 실험에 의해 계측된 최대값에 비해서 파력을 전체적으로 과대산정하는 경향이 나타났다. 그러나 유수실 폭이 넓은 경우, Fig. 7(b)에서 처럼 파력이 작은 경우에는 Takahashi 공식이 과대산정하지만, 파력의 크기가 증가함에 따라 Takahashi 공식과 실험값의 차이가 감소하여, 계측파력 값이 600 N/m 이상인 경우에는 오히려 과소산정하는 경향이 나타났다. 또한, 유수실 폭이 넓은 경우 FRE 의 값이 100~800 N/m로, 유수실 폭이 좁은 경우에 비하여 훨씬 넓은 범위에 분포하였다. 이 경우도 케이슨 전면벽의 경우와 마찬가지로 수정 일치지수 및 결정계수의 값 모두 유수실 폭이 넓은 경우가 더 좋은 결과를 나타내었다.

Fig. 7.

Comparison of the maximum horizontal forces between the results by the experiments and Takahashi formula for the wave force on the rear wall: (a) B=14 cm and (b)B=28 cm.

한편, Fig. 8에는 케이슨 전체에 작용하는 총 수평파력에 대한 결과를 제시하였다. 유수실 폭이 좁은 경우와 넓은 경우를 비교해 보면 전반적으로 유수실 폭이 좁은 경우의 예측값과 실험값이 더 잘 일치하였으며, 이러한 결과는 Table 2에 제시된 수정 일치지수와 결정계수의 값에서도 확인할 수 있다. Figs. 6~7에 보인 결과에서는 유수실 폭이 넓은 경우 Takahashi 공식의 예측값과 실험값의 차이가 더 적게 나타났지만, Fig. 8에 보인 총 수평파력에서는 이와는 반대경향이 나타났는데, 그 이유는 Takahashi 공식에 유수실 폭의 영향이 고려되어 있지 않기 때문인 것으로 판단된다. 즉, 전면벽 및 후면벽에서의 파력이 최대가 되는 위상에서 나머지 벽(즉, 후면벽 및 전면벽)에 작용하는 파력은 Takahashi 공식에 의한 예측값과 실험값의 차이가 크게 나타나기 때문이다. 실제로 많은 경우에 총 수평파력이 최대가 되는 위상은 Crest-IIb 인데, 이 때 전면벽의 파력은 유수실 폭 및 주기에 따라 상당한 값의 변화를 보였다.

Fig. 8.

Comparison of the maximum horizontal forces between the results by the experiments and Takahashi formula for the wave force on the whole caisson: (a) B=14 cm and (b)B=28 cm.

케이슨 전체에 작용하는 총 연직파력에 대한 결과는 Fig. 9에 나타내었다. 유수실 폭이 좁은 경우와 넓은 경우 모두 전체적으로 관측값에 비해 예측값이 큰 경향이 있었다. 또한 유수실 폭의 영향을 고려하지 않은 Takahashi 공식에서는 두 경우 모두에서 공식에 의한 계산값이 200~650 N/m 범위였지만, 실험값은 폭이 좁은 경우에서 100~550 N/m, 넓은 경우에서 200~600 N/m 으로 폭이 좁은 경우가 조금 더 작은 경향을 보이는 차이가 있었다. 이러한 이유로 유수실 폭이 넓은 경우의 수정 일치지수와 결정계수가 더 좋은 결과를 나타냄을 Table 2에서 확인할 수 있다. 또한, Table 3에 제시된 바와 같이 실험에 의한 계측값과 Takahashi 공식 산정값을 비교해 보면 평균적으로 Takahashi 공식이 실험값에 비해 30% 이상 연직파력 값을 크게 산정하는 것으로 나타났다.

Averages of normalized wave force

한편, Table 3에는 본 연구에서 검토한 세 가지 파압공식에 의한 최대파력 계산값을 실험값으로 나누어 무차원화한 값의 평균을 제시하였다. Table 3에서 Takahashi 공식이 실험값에 비해서 파력을 과대산정함을 확인할 수 있으며, 이러한 결과는 Takahashi 공식에서 제안된 파력 보정계수(λ) 값이 파력을 보수적으로 산정하고 있음을 의미하는 것이다.

4.2 Tabet-Aoul & Lambert 공식과의 비교

4.1절과 마찬가지 방법으로 Tabet-Aoul & Lambert 공식을 이용하여 유공케이슨 전면벽, 후면벽 및 전체 케이슨에 작용하는 수평파력을 실험값과 비교한 결과를 Figs. 10~12에 나타내었다. 위에서 살펴본 것처럼 전면벽 유공부의 유공률에 따른 전면벽 및 후면벽 수평파력 계측값에 차이가 거의 없었기 때문에 Figs. 6~9와는 달리 Figs. 10~12에서는 전면벽 유공률 값에 상관없이 유수실 폭의 차이에 따라서만 서로 데이터를 구별할 수 있도록 나타내었다. Figs. 10~11로부터 전면벽에 작용하는 최대 수평파력은 대체로 Tabet-Aoul & Lambert 공식에 의한 계산값보다 계측값이 더 크며, 후면벽의 경우에는 이와는 반대 경향이 나타남을 확인할 수 있다. Tabet-Aoul & Lambert 공식에 의한 전면벽 파력이 실험값보다 작게 산정된 이유는, 이 공식이 근본적으로 전면벽 전체가 슬릿으로 이루어진 형태의 구조물에 대한 실험으로부터 얻어진 식이기 때문인 것으로 추정된다. 본 실험에서 사용된 유공케이슨 모형은 전면벽 하부가 무공벽으로 구성되어있는 부분유공 형태의 구조물이므로 전면벽에서의 파력 감쇠효과가 전유공 형태에 비해서는 상대적으로 작아서 전면벽에 작용하는 최대 파력이 전유공 케이슨에 비해 더 크게 계측되었을 수 있다.

Fig. 10.

Comparison of the maximum horizontal forces between the results by the experiments and Tabet-Aoul & Lambert formula for the wave force on the front wall.

Fig. 12.

Comparison of the maximum horizontal forces between the results by the experiments and Tabet-Aoul & Lambert formula for the wave force on the whole caisson.

Fig. 11.

Comparison of the maximum horizontal forces between the results by the experiments and Tabet-Aoul & Lambert formula for the wave force on the rear wall.

전체적으로는 Fig. 12에 보인 것처럼 Tabet-Aoul & Lambert 공식에 의한 최대 총 수평파력 예측값은 실험 계측값에 비해 더 크게 평가되었으며, Table 3에 제시된 것처럼 실험값과의 차이는 Takahashi 공식보다 더 큰 편이었다. 한편, Table 2에 제시된 것처럼 최대 총 수평파력의 수정 일치 지수의 값은 0.43으로 Takahashi 공식에 비해서 더 낮았으며, 결정계수의 값은 0.69로 더 높았다. 다만, Tabet-Aoul & Lambert 공식의 경우 유수실 폭이 좁은 경우와 넓은 경우의 결과가 큰 편차를 나타내지는 않았는데, 이는 이 공식에서 유수실 폭의 영향을 고려하고 있기 때문인 것으로 보인다. 그러나 Takahashi 공식과 마찬가지로 전면벽과 후면벽을 구별하여 각각의 파력을 계산한 결과, 총 수평파력 결과와는 상당히 다른 경향을 나타내었는데, 이는 위에서 언급한 것처럼 전면벽 또는 후면벽 중 어느 한 쪽에서 파력이 최대가 되는 경우에 다른 벽체에 작용하는 파력에 대한 산정에 불확실성이 많이 내포되어 있기 때문일 것으로 판단된다.

4.3 DUT 공식과의 비교

마지막으로 DUT 공식에 의해 추정된 최대 수평파력 및 최대 연직파력과 실험에서의 계측값을 비교하여 나타낸 그림을 Figs. 13~14에 나타내었다. DUT 공식의 경우 개별 벽체에 대한 최대 수평파력 값을 따로 제공하고 있지 않기 때문에 총 수평파력에 대한 결과만을 제시하였다. 수평파력 값은 Fig. 13에서 알 수 있듯이 유수실 폭이 좁은 경우와 넓은 경우의 결과에도 큰 차이가 없을 뿐만 아니라, 실험값과 적합도 측면에서도 앞선 두 공식에 비해서 더 나은 결과를 보여주고 있다. Table 2에 제시된 일치지수와 결정계수 값도 각각 0.72와 0.71로 Takahashi 공식 및 Tabet-Aoul & Lambert 공식에 비해서 최대 수평파력을 더 정확하게 산정하고 있음을 알 수 있다.

Fig. 13.

Comparison of the maximum horizontal forces between the results by the experiments and DUT formula for the wave force on the whole caisson.

Fig. 14.

Comparison of the maximum vertical forces between the results by the experiments and DUT formula for the wave force on the whole caisson.

한편, 최대 연직파력은 Fig. 14에 보인 것처럼 수평파력에 비해서는 산포도가 더 큰 편이었다. 그렇지만 최대 수평파력과 마찬가지로 유수실 폭에 따른 차이가 거의 없었으며, Table 2에 제시된 것처럼 수정 일치지수 및 결정계수의 값이 Takahashi 공식에 의한 결과보다는 더 커서 상대적으로 실험값의 경향을 더 잘 반영하는 결과를 나타내었다. DUT 공식에 의한 최대 연직파력을 실험값으로 나누어 무차원화한 결과는 Table 3에 제시된 것처럼 1보다 작은 값을 나타내었다. 이는 실험값에 비해 파력을 작게 산정함을 의미하는데, 위에서 살펴본 것처럼 Takahashi 공식의 경우 연직파력을 과대산정한 결과와는 상반되는 경향이다.

비록, DUT 공식이 Takahashi 공식이나 Tabet-Aoul & Lambert 공식에서와 같이 각 벽별 최대파력을 제시하지 않고는 있지만 케이슨 안정성 측면에서 궁극적으로 의미있는 최대 수평파력과 최대 연직파력은 가장 적합도가 높은 결과를 제공한다는 점은 주목할 만하며, 향후 DUT 공식의 적용에 대해서 보다 적극적으로 고려할 필요가 있다고 판단된다.

5. 결 론

본 연구에서는 유공케이슨에 작용하는 수평 및 연직파력을 계산하기 위해 제안된 3가지 서로 다른 공식들의 특징을 살펴보고, 수리모형실험을 통해 계측된 결과와 비교하였다.

Takahashi 공식은 현업에서 가장 널리 활용되고 있는 공식이며, 수평파력의 경우 전면벽 유공부, 무공부 및 후면벽에 작용하는 파력을 독립적으로 평가하기 때문에 설계 시 각 벽체별 작용 파력을 개별적으로 고려할 수 있다는 장점이 있다. 그러나 전체적으로 실험값 대비 최대 수평파력은 22%, 최대 연직파력은 34% 크게 산정하는 결과를 나타내어 실험값에 비하여 파력을 상당히 보수적으로 산정하는 경향을 보였다. 이처럼 Takahashi 공식이 상당히 과도한 정도로 안전하게 파력을 평가하게 되는 이유는 이 공식의 경우 파의 작용 위상에 따라 각 벽별로 작용하는 파력을 독립적으로 산정하는데, 이 과정에서 각각의 개별 벽체에 대해서 보수적으로 평가된 파력이 합산되면서 결과적으로는 파력을 과대산정하는 효과가 더 커지는 것이라고 해석할 수 있다.

Tabet-Aoul & Lambert 공식도 유공케이슨 전면벽과 후면벽을 따로 계산한다는 점에서는 Takahashi 공식과 유사하지만, 파력 산정 시 유수실 폭의 영향을 고려할 수 있도록 공식이 제안된 점이 Takahashi 공식과 차별화되는 특징이다. 그러나 이 공식은 연직파력에 대해서는 파력 산정식을 제시하지 않는다는 한계점이 있다. 실험값과의 비교 결과 전면벽은 파력을 과소산정하며, 후면벽은 과대산정하여 매우 심하게 상반되는 결과를 보였으며, 결과적으로는 최대 수평파력을 34% 크게 산정하는 것으로 나타났다. 위에서 기술한 것처럼 이 공식은 전면벽 전체가 슬릿으로 이루어진 형태의 구조물에 대한 실험으로부터 얻어진 식이기 때문에 이러한 결과가 얻어진 것으로 판단되며, 현업에서 많이 설계가 이루어지는 유공케이슨의 경우 대부분 전면벽 하부가 무공벽이므로 Tabet-Aoul & Lambert 공식은 실무 적용성이 높지 않다고 판단된다. 다만, Takahashi 공식에 비해서는 유수실 폭에 따른 실험값과 계산값의 편차는 거의 없었는데, 이는 Tabet-Aoul & Lambert 공식에서 유수실 폭의 영향이 고려된 효과가 반영된 결과로 추정되며, 향후 더 다양한 실험 자료 분석을 통해서 이러한 결과를 확증할 필요가 있다. 그리고 이러한 경향성이 확인된다면 현업에서 단일유수실 케이슨 설계 시 주로 사용되고 있는 Takahashi 공식에 유수실 폭의 영향을 고려하는 효과를 포함시키는 것도 가능할 것이다.

마지막으로 DUT 공식은 유수실 폭 뿐만 아니라, 파고, 파장, 수심에 관련된 다양한 무차원변수를 공식에서 고려한 점이 주요한 특징이며, 앞의 두 공식과는 달리 유공케이슨의 개별 벽체에 작용하는 파력에 대한 산정값을 제시하고 있지는 않다. 그러나 최대 수평파력 및 연직파력 계산값은 세 가지 공식 중에서 실험값과 가장 잘 일치하는 결과를 나타내었다. 다만 실험값에 비해서는 수평파력은 6%, 연직파력은 12% 과소산정하는 결과를 보였는데, 이는 설계 시 케이슨 안정중량을 과소산정할 수 있는 위험성이 있으므로, 향후 이 공식의 적용성에 대해서는 추가적인 실험 및 연구를 통해서 충분한 검증이 이루어질 필요가 있다. DUT 공식의 경우 유공케이슨의 형상에 관계된 다양한 파라메터 값을 결정하는 과정 및 그 기초가 되는 수리실험에 대한 정보가 부족하여 이 공식이 유도된 과정을 명확하게 이해할 수 없다는 한계가 있으며, 향후 이에 관련된 추가적인 조사 및 연구가 필요한 실정이다. 이러한 과정 및 추가적인 실험연구를 통해 조금 더 합리적으로 DUT 공식에서 고려하고 있는 무차원변수들의 값을 재산정하게 된다면 약간 변형된 형태로서 이 공식을 설계 실무에서 활용할 수 있는 가능성도 열려있다고 할 수 있다.

Acknowledgements

본 연구는 한국해양과학기술진흥원의 첨단항만건설기술개발사업 ‘기후변화대응 항만설계기준 개선 방안연구(1단계) (PM58440)’와 한국해양과학기술원 주요사업(PE99325)의 지원을 받아 수행되었습니다.

References

1. Fugazza M., et al, Natale L.. 1992;Hydraulic design of perforated breakwaters. J. Wtrwy., Port, Coast. Oc. Engrg 118(1):1–14. 10.1061/(ASCE)0733-950X(1992)118:1(1).
2. Goda Y.. 1974;“New wave pressure formulae for composite breakwaters”. Proc. 14th Int. Conf. Coast. Engrg :1702–1720.
3. Huang Z., Li Y., et al, Liu Y.. 2011;“Hydraulic performanca and wave loadings of perforated/slotted coastal structures: A review”. Ocean Engrg 38:1031–1053. 10.1016/j.oceaneng.2011.03.002.
4. Jarlan G.E.. 1961;“A perforated vertical wall breakwater”. The Dock and Harbour Auth 41(486):394–398.
5. Kakuno S., Oda K., et al, Liu P.L.F.. 1992;“Skattering of water waves by vertical cylinders with a backwall”. Proc. 23th Int. Conf. Coast. Engrg :1258–1271.
6. Li Y.C.. 2007;“Interaction between waves and perforated-caisson breakwaters”. Proc. 4th Int. Conf. Asian Pac. Coasts :1–16.
7. Oh S.-H., Ji C.-H., Oh Y.-M., et al, Jang S.-C.. 2013;“Influence of wave chamber slab on wave pressure on first and second wall of perforated caisson breakwater”. J. Korean Soc. Civil. Engrs :2317–2328. (in Korean).
8. Oh S.-H., Ji C.-H., Oh Y.M., Jang S.-C., et al, Lee D. S.. 2014;“Comparison of maximum horizontal wave force acting on perforated caisson breakwater with single and double chamber”. J. Korean Soc. Coast. Oc. Engrg 26(5):335–341. (in Korean). 10.9765/KSCOE.2014.26.5.335.
9. Park W.S., Chun I.S., et al, Lee D.S.. 1993;“Hydraulic experiments for the reflection characteristics of perforated breakwaters”. J. Korean Soc. Coast. Oc. Engrg 33(6):198–203. (in Korean).
10. Tabet-Aoul E-H., et al, Lambert E.. 2003;“Tentative new formula for maximum horizontal wave forces acting on perforated caisson”. J. Wtrwy., Port, Coast. Oc. Engrg 129(1):34–40. 10.1061/(ASCE)0733-950X(2003)129:1(34).
11. Takahashi S., Shimosako K., et al, Sasaki H.. 1991. “Experimental study on wave forces acting on perforated wall caisson breakwaters”. Report of P.H.R.I 30(4)3–34. (in Japanese).
12. Takahashi S., et al, Shimosako K.. 1994;“Wave pressure on a perforated wall caisson”. Proc. Int. Conf. Hydro-tech. Eng. Port Harbor Const. (HYDRO-PORT ‘4) :747–764.
13. Willmott C.J., Robeson S.M., et al, Matsuura K.. 2012;“A refined index of model performance”. Int. J. Climatol 32(13):2088–2094. 10.1002/joc.2419.

Article information Continued

Table 1.

Correction coefficients for perforated caissons given by Takahashi

Crest-I Crest-IIa Crest-IIb
Slit Wall λs1 0.85 0.7 0.3
λs1 0.4 (α*≤0.75)
0.3⁄α*(α*>0.75)
0 0
Impermeable Front Wall λL1 1 0.75 0.65
λL2 0.4 (α*≤0.5)
0.2⁄α* (α*>0.5)
0 0
Wave Chamber Rear Wall λR1 0 20B30L' (B⁄L'≤0.15)
1.0 (B⁄L'>0.15)
1.4 (HDh≤0.1)
1.6-2HDh (0.1<⁄HD/h<0.3)
1.0 (HDh≥0.3)
λR2 0 0.56 (α*≤25⁄28)
0.5⁄α* (α*>25⁄28)
0
Wave Chamber Bottom Slab λM1 0 20B30L' (B⁄L'≤0.15)
1.0 (B⁄L'>0.15)
1.4 (HDh≤0.1)
1.6-2HD(0.1<HDh<0.3)
1.0 (HDh≥0.3)
λM2 0 0 0
Uplift Force λU1 1 0.75 0.65

Fig. 1.

Pressure distributions on the perforated caisson at Crest-IIb by Takahashi formula.

Fig. 2.

Pressure distributions on the perforated caisson by Tabet-Aoul & Lambert formula.

Fig. 3.

Experimental setup (unit: mm).

Fig. 4.

(a) Side view of the narrow chamber caisson, (b) Side view of the wide chamber caisson, (c) Front view of the caisson with different front wall porosity (unit: mm).

Fig. 5.

Installation of the model and arrangement of the load cells.

Fig. 6.

Comparison of the maximum horizontal forces between the results by the experiments and Takahashi formula for the wave force on the front wall: (a)B=14 cm and (b)B=28 cm.

Fig. 7.

Comparison of the maximum horizontal forces between the results by the experiments and Takahashi formula for the wave force on the rear wall: (a) B=14 cm and (b)B=28 cm.

Fig. 8.

Comparison of the maximum horizontal forces between the results by the experiments and Takahashi formula for the wave force on the whole caisson: (a) B=14 cm and (b)B=28 cm.

Fig. 9.

Comparison of the maximum vertical forces between the results by the experiments and Takahashi formula for the wave force on the whole caisson: (a) B=14 cm and (b)B=28 cm.

Table 2.

Refined index of agreement(dr) and coefficient of determination(R2)

dr R2
B = 14 cm B = 28 cm Total B = 14 cm B = 28 cm Total
TK FF 0.69 0.73 0.70 0.60 0.79 0.68
FR 0.22 0.52 0.43 0.64 0.76 0.63
FH 0.72 0.44 0.57 0.66 0.57 0.53
FV 0.62 0.80 0.71 0.47 0.69 0.55
TL FF 0.43 0.47 0.45 0.59 0.78 0.64
FR -0.30 0.44 0.16 0.64 0.76 0.61
FH 0.45 0.42 0.43 0.67 0.72 0.69
DUT FH 0.70 0.74 0.72 0.71 0.72 0.71
FV 0.75 0.86 0.81 0.50 0.82 0.68

Table 3.

Averages of normalized wave force

B = 14 cm B = 28 cm Total
TK FFTK/FFE 1.15 1.13 1.14
FRTK/FRE 1.29 1.48 1.38
FHTK/FHE 1.10 1.34 1.22
FVTK/FVE 1.47 1.18 1.34
TL FFTL/FFE 0.68 0.73 0.71
FRTL/FRE 1.51 1.54 1.53
FHTL/FHE 1.30 1.37 1.34
DUT FHDUT/FHE 0.92 0.96 0.94
FVDUT/FVE 0.81 0.94 0.88

Fig. 10.

Comparison of the maximum horizontal forces between the results by the experiments and Tabet-Aoul & Lambert formula for the wave force on the front wall.

Fig. 11.

Comparison of the maximum horizontal forces between the results by the experiments and Tabet-Aoul & Lambert formula for the wave force on the rear wall.

Fig. 12.

Comparison of the maximum horizontal forces between the results by the experiments and Tabet-Aoul & Lambert formula for the wave force on the whole caisson.

Fig. 13.

Comparison of the maximum horizontal forces between the results by the experiments and DUT formula for the wave force on the whole caisson.

Fig. 14.

Comparison of the maximum vertical forces between the results by the experiments and DUT formula for the wave force on the whole caisson.