Journal of the Korean Geosynthetics Society. 30 December 2016. 1-8
https://doi.org/10.12814/jkgss.2016.15.4.001

ABSTRACT


MAIN

  • 1. 서 론

  • 2. 폐석회의 지반공학적 특성에 관한 연구 동향

  •   2.1 압밀특성

  •   2.2 강도특성

  • 3. 폐석회의 지반공학적 특성분석

  •   3.1 물리적특성

  •   3.2 압밀특성 분석

  •   3.3 강도특성 분석

  • 4. 폐석회 매립지반의 압밀침하 특성

  •   4.1 이론식에 의한 압밀침하량 산정

  •   4.2 수치해석에 의한 압밀침하량 산정

  • 5. 결 론

1. 서 론

우리나라에 산재되어 있는 폐석회는 소다회 생산과정에서 발생되어 생산위주의 공장운영과 관리소홀 등으로 석회산업의 부산물인 부산석회 즉 폐석회가 적당한 대비책이 없이 방치되어온 실정이다. 우선 폐석회의 발생현황을 살펴보면 일반적으로 시멘트공장, 생석회공장, 비료공장, 화학공장 등에서 발생한다. 그러나 발생량에서 화학공장에서 발생되는 폐석회의 양이 95%이상을 차지한다. 국내의 사업장에서 발생되는 폐석회는 주로 소다회 생산 공정에서 배출되는 폐석회(CaCO3)로 적당한 처리기술의 미비로 공장주변에 방치되어온 실정이다. 외국의 경우 폐석회의 처리방안을 살펴보면 대부분 내륙매립과 석회비료로 사용되고, 건축재료, 지반안정제, 염분토양의 중화재, 폐기물 매립장의 복토재 등으로 많은 연구가 진행되고 있다. 국내의 화학공장이나 비료공장에서 발생된 폐석회는 재활용이 극히 미미한 수준이고 대부분이 인근부지에 야적된 상태로 주변미관의 파괴와 토지이용계획에 문제점으로 대두되고 있다. 최근 친환경 및 지속가능한 개발이라는 패러다임에 따라 폐자원의 처리와 그 활용 방안에 있어서 다양한 기술과 방법에 있어서 많은 변화를 보이고 있다. 최근에는 화학공장에서 발생된 폐석회를 매립 후 주변과의 조화와 효율적인 토지이용을 위하여 최적의 처리시스템을 이용하여 이를 공원이나 체육부지로 개발하는 매립공사가 활발히 진행 중에 있다. 따라서, 폐석회 매립지가 안정화되는 동안 매립지의 침하 및 지반계수 등의 지반 공학적 특성에 관한 보다 많은 연구가 필요하다. 매립지를 건설부지로 재활용할 경우에는 폐기물 공학적 특성에 보다 많은 이해와 지반을 개량하는 공법들에 대한 연구들도 병행되어야 한다. 따라서, 이 연구에서는 폐석회로 구성된 지반에서 폐석회의 압밀특성과 강도특성을 정확히 파악하고, 그 결과를 이용하여 폐석회 지반의 침하량을 산정하고자 하였다.

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F1.jpg
Fig. 1

Creep deformation of clay for elapsed time(Ladd et. al., 1977)

2. 폐석회의 지반공학적 특성에 관한 연구 동향

2.1 압밀특성

현재 폐석회에 관련된 많은 연구가 진행되고 있다. 폐석회와 관련한 연구는 대부분이 폐석회 재활용 및 특성에 관련된 내용이 대부분으로 폐기물 매립지반에 대한 연구는 진행되고 있으나 폐석회 매립지반의 침하특성에 대한 연구는 아직 미미한 상황이다. 이 연구에서는 폐석회 매립지반의 압밀 침하특성과 유사한 특성을 나타내는 점성토 지반의 압밀침하특성에 대한 연구동향 및 폐석회 특성에 대한 기존의 연구에 대하여 고찰하였다. 점성토지반의 압밀거동은 1923년 Terzaghi(1948)의 압밀이론 발표이후 현재까지 많은 논문이 발표되고 연구되었다. 이차압밀거동을 설명하는 이론은 Mesri and Choi(1985)는 일차압밀이 끝난 시점에서의 변형은 유일하다고 보고하였다. 또한, 이차압축지수에 대한 압축지수의 비는 주어진 점토에 대하여 일정하다고 하였다. 국내에서 Shin et al.(1996)은 부산석회의 토질역학적 특성에 관한 연구를 통하여 도로의 기층재로의 활용가능성을 제시하였다. 또한, 석회안정처리공법에 근간을 두고 부산석회를 준설토와 혼합하여 연약지반 개량재로 활용하는 방안을 제안하였다. 현장 실험을 통해 부산석회와 점토 혼합토는 혼합비 20%에서 가장 큰 강도를 보였으며, 시간이 경과함에 따라 강도는 증가하는 것을 연구하였다. 그러나 현재까지도 폐석회 매립지반의 장기 침하 양상이나 메커니즘을 명확히 규명하지 못하고 있다. 따라서, 폐석회의 지반공학적 특성에 대한 이해가 필요하고, 기존 압밀이론의 활용 및 보완할 수 있는 연구가 필요하다.

2차 압밀은 1차 압밀 동안의 체적변형의 연속이며, 상당히 느린 속도로 진행된다. 또한 과잉간극수압이 소산된 이후에도 일정한 유효응력이 계속 유지된다는 점에서 1차압밀과는 크게 다르다. 이러한 2차 압밀의 발생은 각각의 입자의 변형에 대하여 입자들의 상대적인 변위로 정의되는 것이 일반적이다. 이러한 이유로 즉시침하량이나 1차 압밀 침하량이 큰 흙일수록, 또는 높은 예민비를 보이는 점토일수록 상대적으로 큰 2차 압밀 특성을 나타내는 것으로 알려져 있다(Ladd et al., 1977; Jamiolkowski et al., 1985; Kim, 2005).

2.2 강도특성

비배수전단강도(http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC23E3.gif)는 흙의 고유특성은 아니나 점토의 연경도를 나타내는데 이용되며 시험형태, 경계조건, 재하율, 초기응력상태 등에 영향을 받는다(Bjerrum, 1967). 비배수전단강도가 응력 의존성이므로 일반적으로 비배수전단강도가 측정되는 연직 유효 상재압력(http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC249F.gif)로 표현된다. Skempton et al. (1957)은 정규압밀 점토의 현장 베인 시험에서 결정된 비배수전단강도와 소성지수의 관계를 식 (1)과 같이 제안하였다. 높은 액성지수의 예민 점토는 식 (1)의 비배수전단강도와 소성지수의 관계가 잘 맞지 않으며 강도증가율은 액성지수와의 관계가 더 바람직하다고 보고된다. 정규압밀점토의 삼축압축 시험에서 평가된 비배수 전단강도로부터 산정된 강도증가율과 소성지수의 관계는 Fig. 2와 같다. 또한, Hansbo(1981)는 액성한계를 이용하여 식 (2)와 같이 제안하였다.

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F2.jpg
Fig. 2

Undrained shear strength for palstic index in normally Consolidation Clay (Bjerrum, 1967)

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC252D.gif                      (1)

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC25AB.gif                           (2)

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC25DB.gif                               (3)

여기서, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC260B.gif : 비배수전단강도, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC262B.gif  : 소성지수, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC265B.gif : 액성한계, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC269A.gif : 초기유효응력, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC26DA.gif : 강도증가율, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC270A.gif : 내부마찰각(°)

3. 폐석회의 지반공학적 특성분석

3.1 물리적특성

이 연구는 폐석회 매립지반에 대하여 지층의 연속적인 지반상태를 나타내기 위하여 매립중의 현장에서 시추조사를 실시하였으며, 자연상태의 불교란시료를 채취하여 한국공업규격에 명시된 표준실험방법(KS F)에 의거해 물리적 특성과 역학적 특성에 대한 연구를 수행하였다. 불교란 시료는 4~6개소의 샘플을 채취하였으며, 매립상태를 초기간극비로 가정하여, 느슨하고 단단한 정도에 따라 구분하였다. 이 연구에서 수행한 실내시험 중 물성실험은 5가지의 시험을 수행하였다. 시험결과는 압밀 시험에서 도출한 초기간극비를 기준으로 구분하였다. Table 1은 실내시험을 통한 폐석회의 물리적 특성을 나타내며, Fig. 3은 소성도 표에 의한 폐석회의 분류를 나타내고 있다. 또한, 추가적으로 폐석회의 압밀특성에 영향을 미칠것으로 예상되는 유기물 함량 실험을 수행하였다. 유기물 함유량 실험 결과 폐석회의 유기물함유량은 18.85%로 도출되었다.

Table 1. Result of physical experimental for waste lime
http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Table_KGSS_15_4_1_T1.jpg
http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F3.jpg
Fig. 3

Plastic chart of waste lime

3.2 압밀특성 분석

이 연구에서는 폐석회 지반의 압밀특성을 정확하게 파악하여 지반을 효율적으로 개량하고자 연구대상 폐석회 매립지반에서 채취한 불교란 시료를 대상으로 압밀시험을 실시하였다. 압밀시험은 KSF 2316에 따라 실시하였으며, 압밀시험에 의한 선행압밀압력과 압축지수, 이차압축지수, 초기간극비는 연약지반의 침하특성을 살펴보는데 매우 중요한 자료가 된다. Table 2는 폐석회의 표준압밀시험의 결과를 나타낸 것이다. 이 연구에서는 초기간극비가 이차압축지수에 어떠한 영향을 미치는지에 중점을 두었다. 2차압축비(http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2749.gif)는 Mesri(1973)가 제시한 함수비와 이차압축비와의 관계를 도시한 Fig. 4에 나타낸 바와 같다. 폐석회의 특성은 대부분 유기질 실트질 점토, 보스톤청색 점토의 특성을 가지는 것으로 판단된다.

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F4.jpg
Fig. 4

Relation of water content and secondary compression ratio(Ladd et al., 1977)

Table 2. Result of consolidation test for waste lime
http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Table_KGSS_15_4_1_T2.jpg

3.3 강도특성 분석

이 연구에서는 폐석회의 강도특성을 분석하기 위하여 삼축압축 시험기를 사용하였다.강도 특성 실험은 비배수압밀(CU)실험을 수행하였으며, 실험 방법은 시료를 우선 압밀시킨 후에, 비배수상태에서 압축하여 파괴하였다. 하중재하시 시료내부에서 과잉간극수압이 등압이 되도록 충분히 느린 속도로 재하였다. 재하중에 주응력과 시료의 압축량을 측정하였다. 내부마찰각은 초기간극비가 증가할수록 약간 감소하는 경향을 나타내었으며, 18.50~22.64°의 분포를 나타내고 있다. Fig. 5는 초기간극비에 따른 전응력 경로를 나타낸 것이다. 응력경로는 모어-쿨롱(Mohr-Coulomb) 응력원의 꼭지점을 연결한 것이다. 이렇게 연결한 선은 파괴포락선이 된다. 압밀비배수시험의 전응력 경로는 기울기 1:1이 되는 직선임을 확인할 수 있다. 전응력 경로의 직선 수는 구속응력의 변화횟수와 같게 되므로 실험시 구속응력을 3번 변화하였으므로 직선의 개수는 3개가 된다. 정규압밀점토의 경우 파괴포락선은 원점을 지나게 되나, 이 연구에서 수행한 폐석회의 파괴포락선은 Y축의 절편을 가지고 있으나 기본 물성 실험결과와 압밀실험 결과를 종합하여 볼 때 이 연구대상 지역에 매립되어 있는 폐석회는 정규압밀 상태임을 확인할 수 있다. 강도증가율은 Skempton의 제안식(식 (1)), Hansbo의 제안식(식 (2))과 비배수압밀(CU) 삼축압축시험결과 마찰각을 이용하는 방법(http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2779.gif방법, 식 (3))으로부터 산정하였다. 각 방법에 의한 폐석회 매립지반의 평균 강도증가율은 0.151(Skempton 식), 0.346(Hansbo 식), 0.264(http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC27A9.gif 방법)로 나타났다. 3가지 방법 중에서 Hansbo의 제안식에 의한 방법이 가장 크게 평가되는 것으로 나타났다. 이는 매립된 폐석회 대부분 액성한계가 50%이상의 고소성을 나타내기 때문이며, 일반적으로 Hansbo의 방법은 저소성에 대해서는 과소평가되는 것으로 알려져 있다.

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F5.jpg
Fig. 5

Stree path of waste lime

4. 폐석회 매립지반의 압밀침하 특성

4.1 이론식에 의한 압밀침하량 산정

이 연구에서는 실내에서 수행한 물성실험과 압밀시험의 결과를 이용하여 연구 대상 지반의 침하량을 산정하였다. 침하량 산정은 1차 압밀침하량과 2차 압밀침하량을 각각 계산하여 총침하량을 산정하였다. 침하량 산정에 필요한 물성은 초기간극비에 따라 구분하여 적용하였으며, 상재하중은 매립 완료 후 성토높이가 1.0~5.0m일 때를 고려하여 적용하였다. 성토체의 단위중량은 일반적으로 성토체로 사용되는 풍화토의 단위중량(18kN/m3)을 적용하였다. Table 3은 압밀도와 침하량 계산시 적용한 폐석회의 압밀특성이다. 침하량 산정결과 1차 압밀침하량은 성토높이의 영향을 받는 것을 알 수 있다. 2차 압밀침하량이 1차 압밀침하량에 비하여 크게 산정된 것은 폐석회에 함유된 유기물 함량(18.85%)이 높기 때문으로 판단된다. Table 4는 이론식 (4)과 식 (5)에 의하여 산정된 총 압밀침하량을 나타낸 것이다.

Table 3. Characteristics consolidation applied to settlement estimate ( : sweeling index)
http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Table_KGSS_15_4_1_T3.jpg
Table 4. Total consolidation settlement for theoretical estimation (unit : mm)
http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Table_KGSS_15_4_1_T4.jpg

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2846.gif                      (4)

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2895.gif                          (5)

여기서, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC28B6.gif : 압밀침하량,  http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC28D6.gif  : 압축지수, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2906.gif : 지층두께, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2926.gif : 초기간극비, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2946.gif : 초기응력

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2986.gif : 응력증가량, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC29B6.gif : 2차 압축지수, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC29F5.gif : 1차압밀 완료 후 연약층 두께

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2A44.gif : 1차압밀 완료 후 간극비, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2A84.gif : 2차압밀 시작시간, http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/PIC2AE3.gif : 2차압밀 종료시간

4.2 수치해석에 의한 압밀침하량 산정

이 연구에서 사용한 압밀해석 프로그램(K-Embank)은 윈도우용의 종합적인 연약지반 설계프로그램을 사용하였다. 해석시 배수조건은 일면배수조건으로 해석을 수행하였다. 해석지반 모델링 하는 과정에 있어서 폐석회가 매립되고 있는 실제 현장의 단면을 모델링 하였으며, 심도 21.4m, 폭은 300m의 단면을 구성하였다. 압밀해석 프로그램에서 총하중법 및 개별하중법의 2종류를 선택할 수 있다. 총하중법(Total Load Method)은 각 단계별 지반조건은 초기 원지반 상태로 규정하고, 1단계 하중 재하시 원지반조건에 대해 침하 및 시간 계산, 2단계 하중 재하시 원지반조건에 대해 침하 및 시간 계산, 1단계 침하곡선과 2단계 침하곡선이 2단계 하중 재하 시작일에서 동일한 발생 침하량으로 교차하도록 침하곡선을 작성하는 방법으로 이 연구에서는 총하중법을 적용하였다. 수치해석에 적용한 물성값은 폐석회 물리적 공학적 시험결과를 적용하였다. Fig. 6~Fig. 7은 해석 대상지반을 모델링한 것이다. 모델링시 성토높이를 1.0 ~5.0m로 변화시키면서 모델링을 실시하였다. 이는 실제 매립완료 이후 복토높이를 가정하여 적용하였다. 해석 프로그램에 의한 총 침하량 산정결과는 Fig. 8에 나타내었다. Fig. 9와 같이 이론적 방법에 의하여 산정된 압밀 침하량과 압밀해석 프로그램에 의하여 산정된 압밀 침하량을 비교 분석하였다. 이론에 의한 결과와 압밀해석 프로그램에 의한 총 침하량 결과가 동일하게 산정되었다.

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F6.jpg
Fig. 6

Modeling of 1-layer waste lime

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F7.jpg
Fig. 7

Modeling of Multi layer waste lime

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F8.jpg
Fig. 8

Numerical analysis for total settlement

http://static.apub.kr/journalsite/sites/kgss/2016-015-04/N0150150401/images/Figure_KGSS_15_4_1_F9.jpg
Fig. 9

Comaprison of theoretical and numerical analysis for total settlement

5. 결 론

이 연구는 소다회 생산과정에서 발생한 폐석회의 지반 공학적특성 분석과 압밀침하량을 산정하기 위하여 시추조사, 실내기본물성시험, 그리고 압밀시험과 삼축압축시험을 수행하였다. 압밀 시험 산정 결과를 바탕으로 이론적인 방법으로 압밀침하량을 산정하였으며, 또한 압밀해석 프로그램을 사용하여 폐석회 매립지반의 침하량을 산정하였다. 이와 같은 연구결과를 바탕으로 도출된 결론은 다음과 같다.

(1)압밀비배수 삼축압축시험 결과 폐석회의 파괴포락선은 Y축의 절편을 가지고 있으나 압밀 시험 결과와 기본물성시험 결과를 종합하여 볼 때 폐석회는 정규압밀 상태임을 확인할 수 있었다. 이를 통하여 폐석회는 압밀점토와 유사한 특성을 가지는 것으로 사료된다.

(2)실내실험 결과를 이용한 폐석회 매립지반의 총 압밀침하량을 산정하면, 초기간극비가 1.5이하일 때 성토높이에 따라 24.6cm~76.1cm이며, 초기간극비가 1.5이상에서 1.8이하 범위일 때 성토높이에 따라 63.2cm~173.2cm이고, 초기간극비가 1.8이상일 때 성토높이에 따라 67.6cm ~198.3cm로 산정되었다. 총 압밀침하량은 성토높이에 따라 증가하였고, 초기간극비가 증가하였을 경우 총 압밀침하량은 증가하는 것으로 분석되었다.

(3)성토높이가 증가할수록 침하량이 증가하는 것은 이론적 방법과 프로그램해석에 의한 방법이 동일한 결과를 나타냈다. 초기간극비의 변화에 따른 압밀 침하량은 이론에 의한 결과와 프로그램해석결과 모두 초기간극비가 증가할수록 침하량은 증가하는 것으로 분석되었다. 총 압밀침하량 또한 이론에 의한 결과와 압밀해석 프로그램에 의한 결과가 동일하게 산정되었다.

(4)이 연구결과를 바탕으로 현재 연구 대상지역은 폐석회 매립이 진행중에 있으므로 향후 매립이 완료된 후 압밀침하에 대한 계측을 실시하여 장기적인 압밀침하특성에 대한 연구가 진행되어야 할 것으로 판단된다.

References

1
Kim, Y. T. (2005), “Axisymmetric Nonlinear Consolidation Analysis for Drainage-installed Deposit Considering Secondary Compression”, Journal of Korea Geotechnical Society, Vol.21, No.3, pp.133-140.
2
Shin, E, C., Oh Y. I. and Kim S. H. (1996), “Development of technology soft ground improvement method using waste lime”, Proceeding of Korean Society of Civil Engineers, pp.487-490.
3
Bjerrum, “Engineering geology of Norwegian normally con-solidated marine clays as related to settlements of building”, Geotechnique, Vol.17, No.2, pp.81-118.
4
Hansbo, S. (1981), “Consolidation of fine-grained soils by prefabricated drains”, Proc.10th International Conference Soil Mechanics Eng. Stockholm, Vol.3, paper 12/22.
5
Jamiolkowski, M., Ladd, C. C., Germaine, J. T. and Lancellotta, R. (1985), “New developments in field and laboratory testing of soils”, Proceedings of 11th ICSMFE, San Francisco, Vol.1, pp.57-153.
6
Ladd, C. C., Foott, R., Ishihara, K., Schlosser, F. and  Poulos, H. G. (1977), “Stress-deformation and strength characteristics”, General report, Proceedings, 9th International Conference on Soil Mechanics and Foundation Engineering, Tokyo. Vol.2, pp.421-494.
7
Mesri, G. (1973), “Coeffieient of secondary compression” ASCE, Journal of Soil Mechanics and Foundation Division, Vol.90, N0.SM5, pp.133-156.
8
Mesri, G. and Choi, Y. K. (1985). “Settlement analysis of embankments on soft clays”, ASCE, Journal of Geotechnical Engineering
9
Mesri, G. and Choi, Y. K. (1985). “The uniqueness of end-of primary(EOP) void ration effective stress relationship”, Proceedings, 11th International Confercnce on Soil Mechanics and Foundation Engineering, San Francisco, Vol.2, pp.587-590.
10
Skempton, A. W. and Bjerrum, L. (1957), “A contribution to the settlement analysis of foundation on clay”, Geothechnique, Vol.7, pp.168-178.
11
Terzaghi, K. (1948), Theoretical Soil Mechanics, Jhon Wiley & Sons, New York.
페이지 상단으로 이동하기