본 연구에서는 Oh et al. (2014)의 수리실험 자료 중 이중유수실 케이슨에 대한 실험자료 및 Fig. 3(c)에 보인 것처럼 높은 마운드 조건에서 이루어진 추가 실험자료를 이용하여 이중유수실 케이슨 각 벽체에 작용하는 파력을 산정하였다. Oh et al. (2014)에 제시된 것처럼 파 작용 시 케이슨 벽에 작용하는 힘을 직접 계측하기 위해서 1축 하중계(load cell)를 이용한 파력 계측 시스템을 구성하여 전면벽의 무공부 및 유공부, 중간벽, 후면벽, 그리고 상치구조물에 작용하는 파력이 독립적으로 각각 측정되었다. Oh et al. (2014)에서는 Fig. 3(b)와 같이 마운드가 낮은 조건에서 단일 및 이중 유수실 유공케이슨에 작용하는 수평파력을 상호 비교하는 연구가 이루어졌으며, 본 연구에서의 초점은 이중 유수실 케이슨의 파력 산정식을 Takahashi 파압 공식과 유사한 형태로 제안하는 데 있다.

실험은 길이 53 m, 높이 1.25 m, 폭 1 m인 2차원 조파수조에서 수행되었으며, 실험에 사용한 케이슨 방파제 모형은 전체 높이 53.5 cm, 폭 59 cm, 길이 70 cm로서 투명 아크릴로 제작되어 사석 마운드 위에 설치되었다. Fig. 2에는 실험 모식도를, Fig. 3에는 방파제 모형의 정면도 및 측면도를 제시하였다. Fig. 3에 보인 케이슨 모형은 Takahashi and Shimosako(1994)의 실험 모형과 매우 비슷한 형상이지만, 유수실이 2개인 이중 유수실 구조인 점은 차이가 있다. 방파제 모형은 높이가 서로 다른 2개의 마운드 위에 각각 설치되어 별도로 실험이 수행되었다. 높이가 17.5 cm인 낮은 마운드의 케이슨 전후면으로부터의 길이는 26 cm였으며, 높이가 26 cm인 마운드의 길이는 70 cm였다. 이 두 가지 마운드 형상은 Takahashi and Shimosako (1994)의 마운드 제원을 고려하여 결정되었다. 한편, 유수실 뒤쪽에는 6 cm 높이의 상부구조물이 위치하였다. 케이슨 전면벽 유공부 및 중간벽의 유공률(ε)은 20% 및 30%의 두 가지로 변화시켰다. Fig. 3(a)에 보인 정면도는 전면벽 유공률이 30%인 경우에 해당한다. 한편, 유공케이슨 설치 위치에서의 수심은 51, 55, 59 cm의 세 가지로 다르게 변화시켰다. 이러한 조건에 대하여 주기(T)는 0.77~2.70 s, 파고(H)는 2.8~24.8 cm 범위에서 서로 다른 37개 규칙파를 조파하여 계측을 수행하였으며, 수리모형실험에 대한 보다 상세한 설명은 Oh et al. (2014)에 제시되어 있다.

이렇게 취득된 실험자료 중에서 본 연구에서는 T = 1.7~2.7 s 범위에 해당하는 실험자료만을 분석에 활용하였다. 실험파 주기가 이 범위일 때 유수실 폭(B)과 파장(L)의 비를 나타내는 파라메터인 B/L의 값은 0.05~0.08이 되며, 이는 유공케이슨 방파제의 통상적인 설계파 주기에 따른 값의 범위에 해당한다.

Takahashi and Shimosako (1994)는 단일유수실 케이슨에 작용하는 파의 위상을 정압(positive pressure)이 작용하는 경우 3 가지(Crest I, Crest IIa, Crest IIb), 부압(negative pressure)이 작용하는 경우 3 가지(Trough I, Trough II, Trough III)의 6가지 위상(phase)으로 각각 구분하였다. 유수실이 2개인 이중유수실 케이슨의 경우 유수실을 구성하는 벽체가 3개(전면벽, 중간벽, 후면벽) 존재하므로 케이슨에 정압이 작용하는 경우에 대해서는 보다 세분하여 아래와 같이 5개의 위상을 고려할 수 있다.

여기에서 Crest F, Crest Ri 및 Crest Rp는 각각 Takahashi 파압공식의 Crest I, Crest IIa 및 Crest IIb에 대응되며, Crest Mi와 Crest Mp는 본 연구에서 새롭게 정의된 위상이다. Fig. 4에는 위 다섯 개 위상에서 케이슨에 작용하는 파 형상을 개념적으로 제시하였다. 한편, 이중유수실 케이슨에 부압이 작용하는 경우에 대해서도 비슷한 방법으로 더 세분하여 파의 작용 위상을 고려할 수 있지만 정압이 작용하는 경우에 비해서는 위상차에 따른 파력 변화의 차이가 크지 않기 때문에 본 연구에서는 별도로 세분화된 위상을 정의하지 않았다.

Fig. 5에는 이중유수실 유공케이슨을 구성하는 서로 다른 연직벽체에 작용하는 수평파력 계측 사례를 제시하였다. Fig. 5는 T = 1.70 s, H = 16.5 cm인 파가 작용했을 때 대략 한 주기 정도에 해당하는 시간 동안의 수평파력 시계열을 보여주며, 여기서 F_S는 전면벽 유공부, F_L은 전면벽 무공부, F_M은 중간벽, F_R은 후면벽에 작용하는 수평파력을 나타낸다. 따라서 전면벽에 작용하는 전체 파력은 F_S와 F_L을 더한 값이 된다. 한편, F_H는 이 4개 벽체에 작용하는 파력을 모두 더한 값으로서 유공케이슨에 작용하는 총 수평파력을 의미한다. 각각의 시계열 그림에서 ○ 기호로 표시된 곳은 제시된 시간 동안 각 벽체에서 발생한 최대파력이 발생했을 때의 파의 위상에 해당한다. 즉, F_S와 F_L의 시계열에서는 Crest F에 해당하며, F_M 및 F_R 시계열에서는 각각 Crest Mp 및 Crest Rp에 해당하는 위상이다. 한편, F_M과 F_R의 그림에는 ◇ 기호로 표시된 곳이 있는데, 이는 파가 중간벽 및 후면벽에 부딪칠 때의 위상인 Crest Mi와 Crest Ri에 해당한다. Fig. 5에 보인 그림에서는 중간벽 및 후면벽에서 모두 Crest Mi와 Crest Ri 위상에 비해서 Crest Mp와 Crest Rp 위상에서 계측된 파력값이 각각 더 컸는데, 파 조건에 따라서는 이와는 반대의 경향이 나타나기도 한다. 한편, Fig. 5에 보인 실험파 조건에서는 총파력(F_H)의 최대값은 Crest Ri (F_R의 시계열 그림에서 ◇ 기호로 표시된 위상)에서 나타나게 됨을 확인할 수 있다.

Fig. 6에는 Takahashi 파압공식(Takahashi and Shimosako, 1994)에서 λ_S1 파라메터를 도출할 때와 동일한 방법을 적용하여 마운드 높이가 낮은 경우 Crest F 위상에서 전면벽 유공부에 작용하는 무차원 평균파력강도(∫SE/∫SG)를 파고 대 수심의 비(H⁄h)에 대해서 나타낸 것이다. 여기서 ∫SE는 전면벽 유공부에 작용하는 평균파력강도, 즉 유공부에 작용하는 파력 계측값을 유공부의 면적으로 나눈 값이며, ∫SG는 비슷한 방법으로 Goda의 파압공식에 의해 계산된 평균파력강도를 의미한다. 따라서 ∫SE/∫SG는 전면벽 유공부를 만약 무공벽이라고 가정했을 경우 Goda 공식에 의해 계산되는 파력에 비해서 실제로 작용하는 파력이 어느 정도의 값을 가지는지를 나타내는 지표가 된다. 즉, 이 값이 만약 1이라고 한다면 이것은 전면벽 유공부에 대하여 계측된 파력이 Goda 파압공식에 의한 무공벽 작용 파력값과 동일하다는 의미가 된다. 따라서 이러한 그림으로부터 전면벽 유공부에 대한 무차원 평균파력강도 ∫SE/∫SG의 값을 어떠한 상수 또는 함수 형태로 나타내게 된다면 그것을 Crest F 위상에서의 전면벽 유공부에 대한 파력저감계수인 λ_S1의 값으로 사용할 수 있게 된다.

Fig. 6에서 무차원 평균파력강도(∫SE/∫SG)의 값은 대략 0.5~0.9 범위에 분포하며, 전체 데이터의 평균값은 0.72가 된다. 한편, Takahashi and Shimosako (1994)의 논문에서 이와 대응되는 그림을 살펴보면 전면벽 유공부에 대한 무차원 평균파력강도의 값이 본 실험과 유사하게 0.6~0.9 범위에 분포함을 확인할 수 있다. Takahashi and Shimosako (1994)는 이로부터 전면벽 유공부 파력저감계수 λ_S1의 값을 0.85로 제시하였다. 이 값은 그들의 실험자료 평균값과 최대값 사이에 해당하게 되는데 이 값을 제시하게 된 이유나 근거는 설명되어 있지 않다. 아마도 방파제의 안정성과 직결되는 파력 산정의 중요성을 감안하여 다소 보수적인 관점에서 전체데이터의 평균값보다는 조금 더 큰 값을 주관적으로 채택하였을 것으로 추정된다.

이처럼 파력을 보수적으로 평가하는 관점은 유지하되, 보다 통계적으로 의미있고 일관성 있는 기준을 제시하기 위하여 본 연구에서는 ∫SE/∫SG의 평균(µ)에 표준편차(σ)를 더한 값으로서 파력저감계수 λ_S1의 값을 산정하였다. Fig. 6에 보인 붉은 색 실선이 이에 해당하는 값을 나타낸 것이며, 통계적으로 이는 전체 데이터에서 상위 16%에 해당하는 값을 의미하게 된다. 한편, 단일유수실 케이슨에 대하여 제안된 Takahashi 파압공식에서는 λ_S1의 값으로 0.85를 제시하였으며, 본 연구에서 산정한 파력보정계수 값과 거의 같다. 본 연구에서의 실험은 이중유수실 케이슨에 대한 것이지만, Crest F 위상에서는 파가 처음으로 구조물에 부딪치는 시점이므로 단일유수실과 이중유수실 케이슨에 작용하는 파력이 크게 다르지 않을 것으로 추정할 수 있다.

한편, Crest F 위상에서 전면벽 유공부에 작용하는 충격파성분 (식(1)에서 두 번째 항에 해당하는 성분)에 의한 파력보정계수를 나타내는 λ_S2 의 값을 산정하는 과정 역시 Takahashi and Shimosako (1994)에 제시된 방법을 따랐다. 즉, 마운드가 높은 조건에서 수행된 실험자료에 대하여 유공벽 전면부에 대한 파력 계측값 ∫SE에서 식 (1)에 λ₁ = 0.84, α^*=0을 대입하여 계산한 값(∫SC)을 뺀 후 물의 단위중량(w₀) 및 파고(H)로 나눈 값을 계산하고, 이를 Fig. 7에 보인 것처럼 α_I의 값에 대하여 그래프로 나타내었다. α^*는 Goda 파압공식에서 쇄파의 영향을 나타내는 α₂와 충격파적 영향을 나타내는 α_I 의 큰 값을 채택하도록 정의되어 있으며, 마운드가 높은 경우에는 α^*의 값은 대체로 α_I 의 값에 의해 결정되므로 Fig. 7의 가로축을 α_I에 대하여 나타내었으며, 이는 Takahashi and Shimosako (1994)에서도 동일하게 적용된 방법이다.

Fig. 7로부터 (∫SE-∫SC)/w0H의 값은 α_I 가 증가함에 따라서 근소하게 증가하는 경향이 있음을 확인할 수 있으며, 따라서 Takahashi 파압공식과 마찬가지로 원점을 지나는 선형 회귀식을 도출하였다. Takahashi and Shimosako (1994)에서는 α_I≥ 0.75 범위에서는 λ_S2= 0.4, α^* < 0.75 범위에서는 λ_S2= 0.3/α^*로 제시하고 있는데, 역시 이 경우에도 위 값이 도출된 근거에 대한 명확한 설명은 제시되지 않았으며 실험자료를 바탕으로 경험적인 판단에 따라 결정된 것으로 추정된다. 본 연구에서는 Fig. 6에서와 마찬가지로 다소 보수적인 관점에서 회귀분석 결과로부터 84% 확률을 가지는 신뢰구간(confidence interval)의 상한값에 대응되는 회귀식을 계산하였으며 그 결과를 Fig. 7에 나타내었다. Fig. 7의 회귀직선의 기울기는 0.69이며, 따라서 λ_S2 의 값으로 0.69를 제시할 수 있다. 이 값은 Takahashi 파압공식에 의한 λ_S2 값에 비해서는 더 크며, 따라서 충격파적 성분의 파 작용 시의 파력을 Takahashi 공식보다 더 크게 산정하는 결과가 된다. 이는 본 연구에서 마운드가 높은 조건에서 수행된 실험에서의 계측 파력이 Takahashi and Shimosako (1994)의 실험보다 더 컸음을 시사하는 것이다.

Table 2에는 λ_S1과 λ_S2의 산정 결과를 제시하였으며, Figs. 6~7에 보인 것과 동일한 분석 방법을 통해서 도출된 전면벽무공부의 파력저감계수 λL1 및 λL2의 값도 나타내었다. Crest F 위상에서는 파가 아직 중간벽 및 후면벽에는 도달하지 않은 상태이기 때문에 중간벽과 후면벽에 대한 파력보정계수의 값은 자연스럽게 0이 된다. 한편, Table 2에는 이후에 설명하게 될 Crest Mi에서부터 Crest Rp까지의 위상에 대한 각 벽체 별 파력보정계수 값을 모두 수록하였다.

Fig. 8에는 Cest Mi 위상에서 두 개의 유수실 사이에 위치한 중간벽에 작용하는 계측파력(∫ME)을 이에 대응되는 Goda 파압공식에 의한 계산값(∫MG)으로 나눈 결과를 제시하였다. 여기서 ∫MG는 Goda 파압공식에서 α₂= 0으로 설정하였을 때의 값이다. Fig. 8은 근본적으로 Fig. 6과 동일한 형태의 그림이며, 파 작용 위상이 Crest Mi이고 검토하는 대상 벽체가 중간벽이라는 점만이 다르다. Fig. 8로부터 무차원 평균파력강도 ∫ME/∫MG의 값은 대략 0.1에서 0.7 정도에 분포하였으며, 이는 Fig. 7과 비교해 봤을 때 상대적으로 더 작은 값이다. 즉, Crest Mi 위상에서 중간벽에 작용하는 평균파력강도가 Crest F 위상에서 전면벽 유공부에 작용하는 평균파력 강도에 비해서 더 작다는 것을 알 수 있다. 이 경우에도 4.1절에서 설명한 것처럼 전체 데이터의 평균에 표준편차를 더한 값인 0.48을 λ_M1 의 값으로서 Table 2에 제시하였다. Table 2에는 Crest Mi 위상에서의 λ_M1의 값뿐만 아니라 이 위상에서의 λ_S1, λL1 및 λ_A1 산정 결과도 제시되어 있는데, 이 값들은 유공부 전면벽 및 무공벽, 그리고 첫 번째 유수실의 바닥판에 작용하는 파력에 대해서 각각 Fig. 8에 보인 것과 동일한 분석을 실시하여 얻어진 값이다.

유수실 중간벽에는 실험파 조건에 따라서 매우 큰 충격파 성분에 의한 파력이 작용할 수 있으며, Takahashi and Shimosako (1994)에서 후면벽에 작용하는 충격파 성분을 고려하기 위하여 적용한 방법을 따라 Crest Mi 위상에서의 충격파 성분에 대한 파력보정계수 λ_M2의 값을 산정할 수 있다. 즉, Crest-Mi 위상에서 중간벽에 작용하는 파력의 계측 시계열로부터 0.2초 동안의 이동평균을 시행한 시계열을 뺀 값을 계산함으로써 충격파적 성분(∫MEi)을 분리해내었다. 그리고 이 충격파적 성분을 w₀H로 나누어 무차원화한 값을 Fig. 9에 보인 것처럼 α_I′의 함수로써 나타내었다. 여기서 α_I′는 Takahashi and Shimosako (1994)에 제시된 것처럼 유수실의 수심(d′)과 폭(B_M′) 및 유수실 내 파장(L′)을 고려하여 계산된 값이다. Fig. 9에서 ∫MEi/w0H의 값은 대부분 0~0.4 범위에 분포하였다. 위에서와 마찬가지로 평균에 표준편차를 더한 값으로써 파력저감계수( )를 산정하였으며, 이 값은 0.26이었다. 그런데, 이 경우에는 파력을 무차원화한 방법이 의 경우와는 다르기 때문에 의 값은 Table 2에 보인 것처럼 = 와 같이 표현된다.

한편, Crest Mi 위상에서는 전면벽 유공부와 무공부 및 첫 번째 유수실 바닥에도 파력이 작용하게 되며, Fig. 6 및 Fig. 8 과 마찬가지로 해당 벽체에 대한 계측 파력을 Goda 파압 공식에 의해 계산된 파력으로 무차원화함으로써 파력저감계수를 각각 산정할 수 있다. Table 2에는 이렇게 얻어진 , 및 의 값을 제시하였다.

중간벽에 중복파 성분의 파력이 작용하게 되는 위상인 Crest Mp에서의 파력저감계수 역시 위에서 제시된 방법에 따라서 산정하였다. Fig. 10에는 그 한 예로서 중간벽에 작용하는 파력에 대한 결과를 제시하였다. 즉 Crest Mp 위상에서 중간벽에 대하여 계측된 파력(∫MG)을 Goda 파압공식에 의한 값(∫ME)으로 나눈 값을 그림으로 나타내었다. ∫ME/∫MG의 값은 대체로 0.3~0.6의 범위에 분포하며, Fig. 10에 붉은 선으로 표시된 것처럼 이 위상에서의 의 값은 0.51로 산정되었다. Crest Mp 위상에서 다른 벽체에 작용하는 파력저감계수의 값 역시 Table 2에 수록되었으며, 이 위상에서는 충격파 성분의 파력은 고려하지 않으므로 λ₂의 값은 모든 벽체에서 0이 된다.

Crest Ri 및 Crest Rp 위상에서의 파력보정계수 산정 방법은 4.2절과 4.3절에 제시된 방법과 동일하다. 다만, 이 두 가지 위상에서는 모든 벽체가 파력을 받기 때문에 λ₁의 값은 Table 2에 제시된 것처럼 각 벽체 별로 모두 제시된다. 한편, λ₂의 값은 후면벽에 충격파적 파력이 작용하는 Crest Ri 위상에 대해서만 산정하게 된다.

Fig. 11과 Fig. 12에는 Crest Ri 위상에서 유수실 후면벽에 작용하는 파력에 대한 실험 결과를 나타내었다. Fig. 11은 이 위상에서 λ_R1의 값을 산정하기 위한 그림이며 의 값은 이 위상에서 0.2에서부터 1.6까지 매우 큰 범위에서 변화한다. 여기서 ∫RE는 Goda 파압공식에서 α₂= 0으로 설정하였을 때의 값이며, 이는 Takahashi 파압공식에서 Crest IIa (Fig. 4에서의 Crest Ri) 위상의 파력저감계수 λ_R1을 계산할 때의 과정과 동일하다. 전체 데이터의 평균에 표준편차를 더한 값은 1.04였으며, Fig. 9에 붉은 색 선으로 표시되었다. 또한, Fig. 12는 λ_R2의 값에 대응되는 그림으로서 ∫REi/w0H의 값은 많은 경우 0 근처에 머무르며, 이러한 데이터들은 해당 파 조건에서 충격파가 거의 작용하지 않았음을 의미하는 것이다. 전체 데이터의 평균에 표준편차를 더한 값은 0.46이었으며, 따라서 Table 2에 보인 것처럼 로 표현된다.

마지막으로, Crest Rp에서의 파력저감계수를 산정하는 사례로 유수실 후면벽에 작용하는 파력을 무차원화한 그림을 Fig. 13에 나타내었다. 이 위상에서 ∫RE/∫RG의 값은 대략 0.5에서 1.5사이의 범위에서 분포하며 Fig. 13에 붉은 선으로 표시한 것처럼 의 값은 1.07로 나타났다.

Crest Ri 및 Rp 위상에서 나머지 벽체들에 작용하는 파력보정계수 값들도 이와 동일한 과정을 통해서 Table 2에 제시하였다.