기존의 건설정보 분류체계는 객체들 간의 관계성 표현 제한으로 인해 교량 구조물에 그대로 적용하기에는 한계가 있다. 본 연구에서는 이러한 문제점을 보완하기 위하여 분류체계 기술의 국제표준인 PLIB Part 42를 활용하여 이종의 도메인 간정보공유의 기본이 되는 형상 정보모델링에 특화된 교량 구성요소의 제품 분류체계를 제시하였다. 특히, 제안한 분류체계는 교량 구성요소의 기능적 특징을 고려한 부품의 의미적 정보가 포함 가능하도록 하였다. 또한 본 연구에서 제안한 분류체계를 실제 모델링에 활용하기 위한 기본 프레임워크를 제안하고 이를 활용한 교량 모델을 생성함으로써 제안한 분류체계가 실무에서 활용 가능함을 보였다.

The trawl fishery is an important fishery accounting for around 30.0% of adjacent water fishery production in Korea, and the commercial fishing is associated with high rate of fatal and non-fatal occupational injury. The hazard analysis for the fishermen's safety of offshore trawler was conducted to serve as basic data for improving the healthy and safe working environment of fishermen using fishermen's occupational accidents of the national federation of fisheries cooperatives (NFFC) from 2012 to 2014 (n=464). As a result, the occupational accident occurrence rate of large powered trawl fishery was 241.4‰ in east sea trawl fishery, 6.9 times the rate of that. In addition, death and missing rate was found to have very serious level management to 97.7‱ in east sea trawl fishery of death at 6.2 times. The accident occurred in 91.7 to 100.0% was happened at sea. The slipping, others and struck by object and son on occurred more frequently in order in the frequency of accident occurrence pattern. However, the occurrence rate of death and missing did not match the frequency of accident pattern. In other words, slipping occurred more frequently while death and missing risk was not high. The fall in the waters was low while death and missing risk was high. The results are expected to help identify and assess safety hazard occurred in offshore trawlers.

In military, empty cartridge recovery case of personal small arms is a device used to collect the empty cartridge after consumption of bullets. It has different shapes depending on the shape of small arms. However, It should be designed in a shape that wraps around the outlet for empty cartridge and it should be no restrictions to the movement of small arms. It has been used to spread military. but military is demanding design improvement due to frequent damage of empty cartridge recovery case. In this report, We improved the shape of the product in order to prevent damage to the product and verified through the FEM analysis and prototype test. According to the result of simulation, best modified modelling of this study has been reduced about 19% of the stress compared to initial modelling. In addition, modified products were confirmed durability, wearability and fixability through fire test.

The axial thrust acting on the turbocharger rotor is basically generated by the unbalance between turbine wheel gas forces and compressor wheel air forces. It has a significant influence on the friction losses, which reduces the overall efficiency and performance of high-speed turbocharger. Therefore, it’s important to calculate the thrust forces under operating conditions (surge, choke and etc.) in a turbocharger. The purpose of this paper is the development of numerical simulation methods which were verified by experimental results of axial thrust and thermally induced constraint tests of the turbocharger. The first FE model showed the relationship between thrust forces and strains by calculating the strains on specially designed thrust bearing and were compared with test results. And the second one is to identify the thermally induced strains in order to remove the thermal effects from measured strains. With these models, it’s possible to inversely predict the magnitudes of the axial thrust by directly measured strains and temperatures under operating turbocharger.

In this study, TbDyFe thin films with the thickness of 1000 Å are fabricated by DC magnetron sputtering. TbDyFe thin films are prepared by DC magnetron sputtering method. The pressure of Ar gas below 1.33 kPa and DC input power of 200 W are used for the sputtering conditions. During sputtering process the substrate holder is heated up to 150 ℃. The thin films are deposited to a thickness of 1000 Å on polyimide substrate with a thickness of 2 μm. The fabricated microstructures are observed by X-ray diffraction (XRD) and the film thickness is measured. Magnetostrictions are determined from the curvature of the thin films which are measured by the optical cantilever method. The experimental results are discussed with numerical data.

The fashion product image preference changes depending on one’s lifestyle and personal inclination. Women want to show the fashion product image preference, often through their clothing and makeup choices. Brand image includes those elements related to the brand. Advertising is the primary method for introducing brands. This study was conducted to investigate the effects of fashion product image preference on emphasis of brand image and advertising factor evaluation targeting working women in their 20s and 30s. The fashion product image preference was s sophisticated image. Preference for a sophisticated image has a significant effect on emphasis of brand image when selecting fashion products. Emphasis of brand image has a notable effect on the brand direct advertising factor evaluation. Fashion product image preference has a significant effect on both brand’s direct and indirect advertising factor evaluation. When selecting a fashion product (clothing and cosmetics), brand image importance was found to have a positive effect on a brand’s direct advertising factor evaluation. Therefore, fashion companies should take advantage of their brand logo. Companies should also pay attention to clothing and product containers used in advertising to show the brand. In addition, every company should create an advertising image that represents their overall brand, by using a combination of detailed advertising factor evaluation.

The purpose of this study is to evaluate topographic characteristics revealed in the predicted distribution areas of Aporia crataegi, according to climate change. Towards this end, this study compared the differences of topographic factors, such as altitude, mountain slope and the aspect of slope, in the distribution areas with different potential inhabitation possibilities of the Aporia crataegi. The inhabitation possibilities of the Aporia crataegi were different, according to altitude and topographic slope, and the inhabitation possibility is judged to be affected more by the topographic conditions including altitude and mountain slope than by the aspect of slope. Especially, the inhabitation possibility of the Aporia crataegi was higher in the higher altitude area, as time goes on furthermore. The reason is that the current climate environment, which is suitable for the potential inhabitation of the Aporia crataegi, is forecast to be formed with an area with high altitude. Although the difference in the aspect of slope was not statistically significant according to inhabitation possibility, the reason why the inhabitation possibility of the Aporia crataegi varies in the mainly southeast slope is conjectured to be derived from the warmer heat environmental condition to grow from a larva into an imago. The result drawn in this study is expected to be utilized as basic data to establish a policy soundly preserving and managing the habitat of biospecies in consideration of climate change and topographic conditions in the natural ecosystem field by using the already built up various biological resources information.

시멘트 페이스트의 동적인 거동을 수치적으로 해석하는 방법으로는 전산유체역학(Computational Fluid Dynamics, CFD)를 이용하는 방법과 이산요소법(Discrete Element Method, DEM)을 이용하는 방 법이 있다. 그러나 시멘트 페이스트와 같이 입자의 크기가 작은 건설재료에 대하여 이산요소법을 활용하 기 위해서는 매우 작은 입자(particle)를 적용해야하는데, 이때 필요한 해석시간은 급격하게 증가하기 때 문에 실제 적용에 어려움이 많다. 따라서 합리적인 수치해석 도구로서 이산요소법을 효율적으로 활용하기 위해서는 재료의 동적 특성에 대하여 상사성을 최대한 유지하는 조건에서 최대 입자크기를 합리적으로 결 정하는 것이 가장 현실적인 방법이라 할 수 있다. 이러한 결정조건을 활용하면, 이산요소법 또는 이산요 소법과 전산유체역학과의 결합(Coupling)을 통하여 도로하부의 동공의 발생이나 함몰의 발생, 자가다짐 콘크리트(Self-Compacting Concrete, SCC), 사면에서 흙의 전단파괴 및 3D 프린터에서 사용되는 유동 적 건설재료의 거동을 합리적으로 예측하거나 수치해석을 기반으로 조건에 맞는 재료의 조성을 결정할 있 는 가능성이 매우 높아진다. 영국의 Swansea 대학의 Feng 교수는 상사적 특성(Similarity)를 이용하여 해석시간을 단축시키기 위한 이론적인 접근 방법을 제시하기 위하여 기하학적 변수, 운동학적 변수 및 역 학적 물성 모두에 조정계수를 적용하는 완전 조정법(Exact Scaling Approach)과 특정조건에서 얻어진 역학적 물성을 상수로 취급하고 기하학적 변수와 운동학적 변수에만 조정계수를 적용하는 부분 조정법 (Partial Scaling Approach)을 제시하였다. 본 연구에서는 정적 액체 가교 모형(Static Liquid Bridge Model)에 부분 조정법을 적용하였으며, 이 산요소법을 이용하여 시멘트 페이스트에 대한 미니콘 실험(Mini-Cone Test)을 모사하고자 하였다. 해석 시간의 단축과 결과의 합리성을 판단하기 위하여 입자의 크기를 변화시키기고 이에 대하여 재료의 물성을 역학기반 이론으로 조정하여 모사를 수행하였다.

노면의 요철은 차량의 동역학적 거동에 의해 시트의 진동 또는 움직임을 유발하여 탑승자의 승차감을 저해시키는 요인이 된다. 차량의 동역학적 거동을 단순화한 모델 중 IRI(International Roughness Index)를 비롯하여 널리 적용되고 있는 Quarter-car모델은 그림 1과 같이 네 바퀴 차량의 1/4 모델로서, 현가상질량(sprung mass)의 수직 변위를 계산하는 모델이다. 그러나 실제 차량의 거동은 수직 방향 거동 외에도 네 바퀴의 거동 차이에 의해 그림 2와 같이 중심축에 대한 회전 운동을 하며, 이로 인해 시트에 앉아 있는 탑승자의 회전진동을 유발한다. 탑승자의 진동방향은 3축 병진운동과 각 축의 회전운동으로 그 림 3과 같이 정의할 수 있으며, 직진 주행 중인 차량으로 가정할 경우 z축에 대한 회전운동인 yaw 거동은 발생하지 않는다. 따라서 본 연구에서는 roll, pitch 회전거동과 연직 z방향 거동에 노출된 탑승자의 승차 감을 모사, 평가하기 위해 그림 4와 같이 시뮬레이터를 개발하였다. 또한 sin파형에 반응하는 시뮬레이터 패널 평가결과를 바탕으로 인체에 민감하게 영향을 주는 차량 거동요소의 주파수 특성을 분석하였다. 시뮬레이터 분석을 통하여 roll, pitch 및 z방향 거동의 특성에 따라 시뮬레이터에 탑승한 패널의 정성 적인 승차감 평가가 달라지는 것을 그림 5와 같이 확인하였다. roll과 pitch거동의 경우 8Hz 대역에서, z방향 거동은 4~8Hz대역에서 패널의 승차감에 가장 민감한 영향을 주는 것으로 분석되었다. 또한 양호- 불량의 경계에 해당하는 거동특성을 비교한 결과, 동일한 회전변위에 대해서 전반적으로 pitch의 영향이 더 큰 것으로 나타났다. 향후 본 연구결과를 토대로 승차감을 정량적으로 표현할 수 있는 방안을 마련할 계획이며, 도로이용자에게 보다 우수한 승차감을 제공하는데 기여할 수 있을 것으로 판단된다.

최근 서울 도심 지역에 위치한 주요 포장도로에서의 동공발생 문제가 다양한 원인으로 인하여 이전 보다 더 빈번하게 나타나는 것으로 대두되고 있는 실정이다. 이 문제는 포트홀과 싱크홀 수 증가와 포장 도로에서의 다른 문제들을 야기 시킨다. 사실, 포장도로에서의 동공발생을 일으키는 주체적인 원인은 강 우, 강설과 같은 기후 변화와 밀접한 관련이 있다. 이 연구에서는, 한국형도로포장설계법에 기반하여 개 발된 KPRP(Korea Pavement Research Program : 한국형 포장설계프로그램)를 이용하여 동상방지층 의 구조적 성능을 모니터링 하고, 유한요소해석 프로그램 중 하나인 KICTPAVE 프로그램을 이용하여 KPRP프로그램의 신뢰도를 확인하고자 한다. KPRP 프로그램은 한국형 도로포장설계법에 기반을 두고 개발 된 프로그램 이다. 사실. 이전에 사용되어지던 설계 법은 AASHTO Interim Guide Method(1972) 로, 한국형포장설계프로그램의 개발 전 국내에서 가장 보편적으로 사용된 설계법이다. 이 설계법은 주어 진 조건 하에서만 포장 단면을 설계하는 법을 도출해내기 때문에 다양한 조건에서의 설계와 구조적인 개 량이 어렵다는 지적이 나오게 되었다. 이를 보완하기 위하여 국내에서의 주어진 조건에 적합하게 도로 포장의 횡단을 설계하기 위해 개발된 한국형포장설계프로그램이 나오게 되었다. 본 연구에서는, KPRP 와 KICTPAVE를 이용하여 동상방지층의 구조적 성능을 더욱 정밀하게 연구하게 된다. 한국형포장설계프로그램은 설계조건을 반영한 후, 횡단면 조건 설정, 기상조건 반영, 교통량 반영, 포장 재료 물성치 반영 및 동상방지층 설치 여부 판단과정을 거쳐 공용성을 예측하게 된다. 설계조건 반영은 한 국형포장설계프로그램에서 제공하는 기본 DB를 활용하여 횡단 설계를 하였다. 본 연구는 효과적인 공용성 분석을 위하여 세 가지 경우로 나누어서 해석을 실시하였는데, 그 조건들은 아래 표 1에 나온 것과 같다. 표 1은 동상방지층의 구조적 성능을 한국형포장설계프로그램을 이용하여 분석해낸 결과 이다. 표 1을 보면 두께와 상관없이 모든 구조적 성능이 일치하게 되는 결과를 가져와서 KICTPAVE 프로그램을 이용하여 조 건을 5가지로 나누어 유한요소해석을 실시하였다. 각 층(동상방지층-노상토, 보조기층-동상방지층, 기층 -보조기층, 표층-기층, 표층 사이의 노드와 요소)은 5cm간격으로 0cm에서 20cm 까지 다섯 가지 경우의 조건으로 나누어서 해석을 하였다. 그 결과 보조기층과 동상방지층 사이의 변형률은 동상방지층의 두께가 두꺼워질수록 증가되는 현상을 보였다. 하지만, 이와는 반대로 노상토와 동상방지층 사이에서 발 생하는 변형률은 동상방지층의 두께가 두꺼워 질 수록 감소되는 현상을 보였다.

본 연구에서는 항온기 내의 원통형 아스팔트 혼합물 시편에 대하여 이산요소법(Discrete Element Method, DEM)을 이용한 열전도 모사연구를 수행하였다. 이를 위하여 기존 연구를 통하여 얻어진 분석적 해(Analytical Solution)를 이용한 아스팔트 혼합물 시편의 온도변화 예측결과를 이용하였으며, 이를 이 산요소법에서 얻어지는 온도의 변화와 비교하였다. 이산요소법은 개별 입자의 동적 운동을 모사하는데 적 합하여, 시간단위의 해석시간이 요구되는 열전달 해석에는 적용이 사실상 불가능하다. 그러나 본 연구에 서는 열전달 해석에 요구되는 열전달 계수에 상사적 개념을 적용하여 결과에 오차가 발생하지 않는 한도 에서 입자의 크기를 증대시켰으며, 이에 적정한 열전달 계수, 열용량 등의 계수를 조정할 수 있는 열역학 적인 근거를 확립하였다. 개별요소법을 이용한 열전달 해석은 열역학적 물성이 서로 다른 입자들을 이용 하여 무작위적으로 혼합된 입자 혼합물의 온도변화를 해석하는데 유용하게 적용될 수 있으며, 시간이나 온도의 변화에 따라서 상(Phase)가 변화하는 재료에 대하여 이를 미시적으로 고려할 수 있는 기반으로 활 용될 수 있기 때문에 큰 의미가 있다. 본 연구에서는 기존의 연구에서 결정된 열역학적 물성을 이산요소해석 소프트웨어에 적용하여 온도의 변화를 예측하였으며, 이들 예측값은 실측값 및 분석적 해에서 도출된 외부 및 내부의 온도변화와 매우 근접하는 것으로 나타났다. 운동하는 입자와 외부 유체와의 열전달에 적용하기 위해서는 보다 정밀한 상 사적 이론(Similarity Principle)과 열역학적 물성의 결정이 필요하겠으나, 본 연구에서는 역학기반의 이 산요소 압자의 크기와 물성의 결정방법을 제시하였으며, 이는 향후 확장된 연구범위에 적용할 수 있다.

지난해 2월말 경 용산역 근처 공사현장 인근 인도에서 도로함몰이 발생하여 보행자가 추락하여 부상을 입는 사고가 발생했다. 이를 계기로 서울시는 도로함몰의 조기예방을 위해 세계 최고수준의 기술을 보유 한 일본 업체와 협력하여 차량형 GPR(Ground Penetrating Rader : 전자파를 지표에 투과하여 지하의 빈 공간 형상 데이터를 수집하는 레이더 장치)을 이용한 도로하부 동공탐사 용역을 실시하였다. 서울시 동공탐사용역을 계기로 현장적응 시운전과 기술진의 분석훈련을 거쳐 실질적 탐사는 지난해 12 월부터 시작되었으며 약 4개월 동안 도로함몰 개연성이 높은 주요 간선도로 48km구간에서 105개의 공동 을 발견하였으며 금년 말까지는 약 200여개를 포함해 총 300여개의 공동이 발견될 것으로 예상했다. 이번 조사로 발견된 공동 중 위험도가 가장 높은 A급 공동이 약 60%로 61개가 발견되었고 그보다 낮은 등급의 B급 공동이 35개가 발견되어 5월 말까지 복구할 예정이며 도로함몰 가능성이 그나마 낮은 C급 공동 8개는 연구를 위해 일정 기간 관찰 후 복구할 예정이라고 밝혔다. 서울시 자치구별 공동 발견 현황을 살펴보면 송파구가 약 30%인 32개로 가장 많은 공동이 발견되었고 뒤이어 용산구 21개, 종로구 19개 순으로 나타났 다. 본 연구는 상용 유한요소해석 프로그램인 ABAQUS를 이용하여 공동 규모에 따른 아스팔트 포장도로 의 거동에 대한 분석을 실시하였다. 이를 위해 공동의 규모는 공동 위험도 평가 등급에 따른 A, B ,C 등급 으로 분류하였으며 공동의 발생 범위는 1.5m 이내로 설정하였다. 해석을 위한 단면 구성은 아스팔트 콘크 리트 층, 보조기층, 노상의 구조로 이루어져 있으며 접합상태를 의미하는 tie contact를 선택하고 각 층의 재료에 대한 입력물성은 일반적으로 아스팔트 콘크리트 포장구조체에 사용되는 범위를 적용하였다. 해석 모델은 연직하중만이 작용하는 Axisymmetric으로 구성하였으며 해석에 사용된 하중은 포장도로 설계시 사용되는 원형 등분포하중으로 설정하였다. 해석은 각 층의 단면, 재료물성과 하중 그리고 공동의 규모를 변화시키며 실시하였고 아스팔트 포장체의 표면 처짐, 변형률 등을 분석하였다.

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

This study investigates the analytical study for developing PEB(Pre-Engineered Building) connection Steel damper. The state-of the art of the steel beam-column seismic connection details were studied previously to develop proper damper for PEB connection. Finite element analysis is performed to develop the prototypes of damper under cyclic load. The study parameters of analysis are the shape, length, thickness, orientation and location of damper. As a result of finite element analysis, three prototype damper details for PEB seismic connection are derived. One is tapered plate on the lower flange of the rafter, another is C-plate on the center of panel zone and the other is brace on the beam-column connection.

본 연구에서는 개별요소해석 프로그램인 PFC3D를 이용하여 지오그리드와 재료입자 사이의 상관관계를 지오그리드 반복인장 시험 및 전단시험을 통해 확인 하고자 한다. 재료의 입도 유효입경 D50 (누과통과율 50%)을 기준으로 40mm, 30mm, 20mm, 15mm, 10mm 로 비교 하였고, 합성수지로 제작된 지오그리드 개별요소를 미시 물성치 보정을 통해 형성하였다. 반복인장에서는 변형률 0.5%를 최대 인장 변형률로 설정하여 반복인장을 500회 실시 하였으며, 전단 시험에서는 지오그리드 생성위치를 전단면으로 설정하였고 횡방향으로 변형률 10%까지 전단시켰다. 해석 결과 유효입경 20mm에서 가장 큰 보강효과를 확인 하였고, 40mm 및 30mm에서 보강효과가 감소되는 것을 확인하였다. 지오그리드의 격자 크기를 고려하였을 때, 재료의 유효입경의 크기는 격자의 내접원 직경의 60% 정도 확보 시 보강효과가 가장 좋은 것으로 확인하였다.

PURPOSES : Nowadays, cavity phenomena occur increasingly in pavement layers of downtown areas. This leads to an increment in the number of potholes, sinkholes, and other failure on the road. A loss of earth and sand from the pavement plays a key role in the occurrence of cavities, and, hence, a structural-performance evaluation of the pavement is essential. METHODS: The structural performance was evaluated via finite-element analysis using KPRP and KICTPAVE. KPRP was developed in order to formulate a Korean pavement design guide, which is based on a mechanical-empirical pavement design guide (M-EPDG). RESULTS: Installation of the anti-freezing layer yielded a fatigue crack, permanent deformation, and international roughness index (IRI) of 13%, 0.7 cm, and 3.0 m/km, respectively, as determined from the performance analysis conducted via KPRP. These values satisfy the design standards (fatigue crack: 20%, permanent deformation: 1.3 cm, IRI: 3.5 m/km). The results of FEM, using KICTPAVE, are shown in Figures 8~12 and Tables 3~5. CONCLUSIONS: The results of the performance analysis (conducted via KPRP) satisfy the design standards, even if the thickness of the anti-freezing layer is not considered. The corresponding values (i.e., 13%, 0.7 cm, and 3.0 m/km) are obtained for all conditions under which this layer is applied. Furthermore, the stress and strain on the interlayer between the sub-grade and the anti-freezing layer decrease gradually with increasing thickness of the anti-freezing layer. In contrast, the strain on the interlayer between the sub-base and the anti-freezing layer increases gradually with this increase in thickness.