PURPOSES : In this study, the applicability of the water content, suction, and suction stress in a resilient modulus prediction model for a subbase was reviewed. METHODS : To compare the applicability of water content, suction, and suction stress models for resilient modulus prediction, the suction stress was determined based on the soil water characteristic curve. The model parameters for each approach were derived from the measured resilient moduli. Finally, the relationships between the degree of saturation and resilient modulus were analyzed using the calculated model parameters. RESULTS : Prediction models of the resilient modulus based on water content and suction demonstrated high correlation with measured values, but overestimated the resilient modulus at saturation levels beyond the laboratory testing range. In contrast, the model accounting for suction stress effectively reduced this overestimation, likely owing to a decrease in suction stress as the suction increased. CONCLUSIONS : Based on the above results, the resilient modulus of subbase materials could be estimated through the change in the degree of saturation and the stress-dependent resilient modulus model using the suction stress proposed in this study.

PURPOSES : In this study, an empirical approach was established to estimate the parameters of the resilient modulus based on various geotechnical properties of subgrade soils. METHODS : Multiple regression analyses were performed to analyze the relationship between resilient modulus (k1) and deformation. The most important factors are the #200 sieve passing ratio, moisture content, and dry unit weight of the soil. The applicability of this approach was verified using selected field data and the literature. RESULTS : The correlation between the results predicted using the prediction equation of the model constant (k1) and the actual k1-value was high. The applicability of the prediction equation was considered high owing to its high suitability with the existing data. The range of values obtained using the constant prediction equation of the proposed model was also judged to be reasonable. In the comparison of the CBR value of the subgrade material of the actual design section and the predicted elastic modulus (k1), almost no relationship was observed between the CBR and the model coefficient (k1). Thus, the estimation of the elastic modulus through CBR is likely to contain errors. CONCLUSIONS : Based on these results, the parameters of the universal model can be predicted using the stress-dependent modulus model proposed in this study.

PURPOSES : In this study, the wetting band depths of road slopes estimated using numerical analyses and one-dimensional empirical equations were evaluated. METHODS : The one-dimensional empirical equations used in this study to estimate the wetting band depth were the Pradel and Raad equation, based on modifying the Green and Ampt equation, Lumb's equation, and Sun equation. The numerical analysis of a finite load slope model was carried out using the Seep/w program (2D). In particular, the effect of the initial suction, which indicated the effect of the antecedent rainfall based on the soil–water characteristic curve, was examined as one condition. The results of the wetting band depths obtained using the empirical equations were evaluated and compared with those of the numerical analysis. RESULTS : The wetting band depths obtained using one-dimensional empirical equations were greater than those from the analytical results. In the case of empirical equations, the estimated results obtained for the wetting band depth might be misleading because it has the limitation of being expressed using a one-dimensional equation with an error, owing to several assumptions for the water infiltration phenomenon. It was also found that the accuracy of the wetting band depth was closely related to the results of the soil–water characteristic curve. CONCLUSIONS : Because the wetting band depths obtained using the empirical equation may lead to overestimation, the slope stability could be evaluated as low; however, there was an advantage in terms of inducing conservative design of the road slope. In addition, it was confirmed that the estimated value of the wetting band depth obtained using the Pradel and Raad equation varied with the suction and volumetric function ratios, and further attention should be paid to these two variables.

본 연구는 우리나라 흑비둘기 최대 서식지인 울릉도에서 2016년 3월부터 2018년 2월까지 흑비둘기의 지역별 개체군 변동과 선호하는 생태적 환경에 대하여 조사하였다. 조사기간 동안 40개 정점에서 확인 된 흑비둘기의 누적개체수는 1,956개체가 조사되었고 최대개체수는 2016년 8월에 326개체가 출현하였다. 월별 개체수 변동은 10월부터 출현개체수가 감소하다가 겨울철인 12월에서 2월에는 거의 관찰되지 않았다. 지역별 개체수 비율은 2016년에는 나리분지 지역이 33.43%로 가장 높았고 2017년에는 태하, 학포 지역이 29.39%로 흑비둘기가 많이 출현하였다. 흑비둘기가 선호하는 서식지 특성은 활엽수림 (77.5%), 임분의 영급은 4영급(75%), 흉고직경(DBH)은 중경목(80%), 수관수밀도는 밀밀도를 선호하였다. 해발고는 평균 178.22±36.60 m(SD, n=40, range 10∼472)였고 10∼160 m 지역이 전체 지역의 57.5%를 차지하여 저지대 산림을 선호하는 것으로 나타났다.

PURPOSES : Recently, there has been an increase in the use of discrete randomly distributed fiber materials for reinforcing pavement foundations. However, very limited study has been made on this from the perspective of pavement engineering. Therefore, this study evaluates the performance of soil-geofibers used in pavement foundations as well as the effects of stress dependency with various mixtures. METHODS: To estimate the behavior of soil-geofiber mixtures under traffic loadings, laboratory resilient modulus data for the mechanical characteristics of geofiber mixtures were used, and they were adopted to evaluate the structural response and analyze the stress dependency through 2-D finite element analysis. As the host materials, poorly graded and uniformly graded sand were selected, and each soil was mixed with three different types of fiber, namely monofilament, fibrillated, and tape. RESULTS: The stress dependent response on resilient modulus and Poisson’s ratio were mainly considered by conducting linear and nonlinear elastic analyses on various geofiber mixtures. As a result, it was found that the response and yield function of geofiber mixed layers in pavements were affected considerably depending on the gradation of the soils and the confinement conditions. A small change was found when the particle size was homogenized. CONCLUSIONS : From this, it can be concluded that the finite element model with stress dependency is suitable for estimating the performance on geofiber mixtures. It is also noted that all the responses of geofiber mixtures were relatively sensitive to the gradation of host soils. This indicates that the effects of the nonlinearity and stress-dependency of geofiber mixtures under repetitive loadings could be substantial.

PURPOSES : This study compares the results of field-monitored unsaturated hydraulic data to verify the wetting path estimation method using the drying path on roadside soil slope. METHODS: A method proposed by Fredlund was used to estimate the wetting curves of soil water characteristic curves. Wetting curves were estimated from the drying curve acquired mainly from laboratory testing, and the estimated wetting curve was applied to the infiltration analysis to verify the estimation method. The infiltration analysis results obtained from the drying, wetting, and estimated wetting curves were compared to the field-monitored data at the same location. RESULTS : The volumetric water contents measured at 40 cm were closest to analysis results obtained from the wetting curve, whereas the volumetric water content measured at 80 cm were closest to the analysis results obtained from the drying curve. The results of the estimated wetting curves were between the results of the drying and wetting curves. The measured suction results were within the drying and wetting curves, and the estimated wetting curves were within the drying and wetting curves. Additionally, the safety factor obtained from applying the drying curves was carefully calculated, and the factor of safety of the estimated wetting curves was close to the results of the wetting curves. CONCLUSIONS : Applying a drying curve to calculate the slope stability during the rainy season will cause an overestimation of the factor of safety at roadside soil slopes. However, if the proposed estimation method of the wetting curve can be applied to the numerical analysis, the stability of the road slope can be reasonably estimated.

PURPOSES : Road subsidence occurs owing to road cavities, which cause many social and environmental problems, especially in cities. Recently, road cavities were detected by various ground radars and repair works were carried out against the detected cavities. The condition assessments related to the road cavities are necessary to understand the potential risk of the cavities. Therefore, in this study, a numerical study was performed to assess the various conditions of road cavities. METHODS : The numerical method adopted in this study is the discrete element approach, and it is suitable for analyzing the condition because it can consider the movement of the soil particles in the surrounded cavity areas. In addition, the triaxial test was modeled and performed under various cavity conditions inside the specimens. RESULTS: The conditions of different cavity locations and shapes were analyzed to identify the effect of cavity state. Three general cases of particle size distributions were formulated to identify the effect of surrounding ground conditions. As a result, the degree of decrement and volumetric strain were varied depending on the locations and shapes of the cavity. Only minor changes were observed when the particle size distributions were altered. CONCLUSIONS: The strength reduction was higher when the cavity formed was larger and located in the upper zone. Similar to the cavity shape, strength reduction and volume deformation are more influenced by the width than the length of the cavities. There is an influence from ground conditions such as the particle size distribution, especially on the wide cavity.

도로 하부에 동공이 발생하게 되면 지반재료의 거동이 크게 변화하게 된다. 최근, 동공은 GPR장치를 활용하여 탐지가 가능하며 탐지된 동공은 그 위험도에 따라 적당한 보수 대책을 마련하게 된다. 지하의 동공은 다양한 위치에 다양한 모양으로 생성될 수 있으며 이에 따른 강도 저하가 발생하는데 특히 위치 및 형상에 따라 도로함몰에 미치는 위험도가 다르다. 따라서 도로 하부에 발생한 다양한 동공의 위험도를 평가하기 위한 많은 연구가 실내시험 및 수치해석을 통해 수행되었다. Hewage and Renuka(2012)는 시료 내에 글루코스 블록을 통해 동공을 모사한 삼축압축 시험을 수행 하여 동공이 발생한 지반재료의 공학적 거동을 파악하였다. 그러나 글루코스 블록을 통하여 모사된 동공 은 물에 용해되는 과정을 거치며 초기 블록의 모양과는 다르게 형성되기에 다양한 위치, 모양을 구현하는 데 한계점이 있다. 따라서 이러한 실내 시험의 어려움 때문에 본 연구에서는 수치 해석적 방법을 이용하 여 동공의 위치와 형상이 지반재료의 강도와 변형에 미치는 영향을 알아보았다. 연구에 사용한 해석방법 은 개별요소법(Discrete Element Method)으로 시간의 흐름에 따라 개별 입자의 이동이 발생하여 동공의 변화에 따른 연속적인 해석이 가능하다. 본 연구에서는 동공의 위치 및 형상을 모사한 삼축압축 모델 테 스트를 진행하였다. 한편, 박형민 등(2016)은 DEM기반의 PFC-2D 프로그램을 이용하여 입도분포곡선을 반영한 삼축시험 모델링 해석을 수행하였으며 실내시험 데이터와 비교한 결과 그 타당성을 입증하였다. 이를 바탕으로 특 정 시료에 대한 입도분포 알고리즘을 적용하여 동공의 위치와 구속압 조건, 형상에 대한 영향을 알아보기 위한 수치해석을 진행하였다. 동공의 위치에 따른 영향을 알아보기 위하여 모델링 된 시료에 상, 중, 그리고 하부에 동공을 생성하여 해석을 수행한 결과 동공이 시료의 윗부분에 위치할수록 낮은 축차응력 수치를 보였다. 또한, 실제 동공 이 존재하는 지반 조건을 깊이 별로 모사하기 위하여 시료 내 동공은 동일한 위치에 생성한 후 지반의 깊 이를 구속압으로 표현해 각각 약 1m, 2m, 3m 깊이에 해당하는 구속압 조건에서 해석을 수행하였다. 결 과 동공이 깊은 깊이(75kPa구속압)에 위치한 경우 동공 발생으로 인해 감소한 강도가 얕은 깊이일 때 보 다 낮은 변형률에서 다시 회복되어 증가하는 현상을 확인 할 수 있었다. 또한 지반 내에서 동공은 여러 형태로 생성 될 수 있기 때문에 동공의 면적과 중심점은 동일한 조건에서 동공의 형상을 변화하여 해석을 수행하였으며 그 결과 동공의 초기 강도는 높이보다 폭에 의해 더 큰 영향을 받으며 더 큰 부피 변화를 수반하는 것을 확인 할 수 있었다. 본 연구의 결과를 토대로 도로 하부의 동공을 탐지한 경우 동공의 형상과 위치에 따라 위험도를 평가 하여 그에 따른 관리 대책을 수립 할 수 있을 것으로 기대한다.

도로의 하부지반은 크고 작은 입자들로 구성되어 있으며 크기의 분포에 따라 전체적인 강도 및 변형특 성이 다르게 나타난다. 이를 반영하여 여러 종류의 입도분포로 실험이 이뤄져 왔고 수치해석을 통한 비교 연구도 진행되어져 왔다. 본 논문에서는 수치해석 기법 중 개별요소법(Discrete Element Method)을 이 용한 방법을 활용하여 입도분포의 변화에 따른 하부지반의 전단강도 및 변형특성을 알아보고자 하였다. 입도분포를 표현하기 위하여 Fredlund & Xing(1994)과 Van Genuchen(1981)이 제안한 함수특성곡선 맞춤 식을 도입하여 입도분포에 구현하였다. Van Genuchen의 경우가 입도분포를 표현하는데 보다 더 적합하여 이를 반영하고 맞춤 변수들의 변동에 따른 입도분포 형태를 적용하여 강도 및 변형의 변화를 살펴보았다. 입도분포를 변동시키기 전 기존의 입도분포와 동일한 조건으로 수행된 실내시험결과를 비교하여, 초기 조건과 경계조건이 실내시험과 유사한 결과를 갖도록 모형화 하였다. 분포 모델은 식 (1)과 같이 표현되며, d min 은 입도분포의 최소 입경을 나타낸다. a, n, m의 값은 입도분포 식의 변수를 나타내며 <표 1>은 입도분포 변화에 따른 변수를 나타낸다. 여기서, y는 누가 통과율, x는 입자의 크기, a,n,m은 맞춤변수이다. <표 1>과 같이 n값의 변화에 따른 입도분포곡선을 산정하고, 이를 개별요소법에 적용하여 해석을 실시 하였고 구속압의 변화(100, 200, 300kPa)에 따른 강도와 변형을 확인하여 보았다. 구속압이 100kPa일 때 에는 n값의 변화에 따라 강도의 변화가 비슷하게 나타났으며, 구속압이 높아질수록 강도의 차이가 다 르게 나타나는 것을 확인할 수 있었다. 이때, n값이 커질수록 입도분포에서 세립분이 함유하는 비율이 커 지는 입도분포를 나타내었다. 또한, n값이 커질수록 강도가 크게 산정되는 것을 확인 할 수 있었다. 본 논문에서는 n값의 변화에 따른 입도분포를 변화하였으며, 이때 변화하는 강도의 차이를 확인하여 보았다. n값이 증가할수록 세립분 함유량이 증가하였고, n값의 변화에 따른 강도변화가 있음을 확인할 수 있었다. 추후 연구를 통하여, a값과 m값의 변화에 따른 입도분포를 산정하고, 개별요소법에 활용하여 a와 n, m값의 변화에 따른 강도의 차이를 확인할 수 있을 것으로 판단된다.

The natures of fatigue crack growth under Mode Ⅱ loading are studied. End notched flexure beam specimens were used. The effects of adherend thickness, rubber modification and adhesive thickness on fatigue crack growth were examined. The experimental results show that some of these parameters apparently do affect fatigue crack growth. Resistance to ModeⅡ fatigue crack growth are increased by rubber modification. The effects of adhesive thickness and rubber content on fatigue crack growth were explained by von Mises's equivalent stress using BEM analysis. For unmodified epoxy adhesives, the fatigue crack growth properties under Mode Ⅱ loadings were significantly different in all regions. For rubber-modified epoxy adhesives, they were also different in the first and second regions, but in the third region, they were similar

In this paper, the features of fatigue crack growth under ModeⅠ loading are studied. Double cantilever beam specimens were used. The effects of adherend thickness, rubber modification and adhesive thickness on fatigue crack growth were examined. The experimental results show that some of these parameters do apparently affect fatigue crack growth. Resistance to ModeⅠ fatigue crack growth are increased by rubber modification. The effects of adhesive thickness and rubber content on fatigue crack growth were explained by von Mises's equivalent stress using BEM analysis.

The crack initiation equals to fracture for bonded joint with brittle adhesive. The criterion is formulated in terms of the quasi-stress intensity factor Kp, for the maximum principle stress, that is analogous to the stress intensity factor used to characterize the stress field in the vicinity of bond terminus. Kp is evaluated using a boundary element analysis. The crack initiation at the terminus of adhesive bonded joints is estimated with the critical quasi-stress intensity factor Kp. This method presented here hardly pays attention to the crack propagation. Since there is a large influence of crack propagation on the strength of adhesive joints and structures, crack propagation must be taken into account on strength prediction of bonded joints. The quasi-stress intensity factor Kp for the maximum principle stress can use as the criteria of the crack initiation at the terminus of adhesive bonded joints having various shapes.