• 한국지반공학회지:지반
• /
• 제10권2호
• /
• pp.41-56
• /
• 1994

지하수위 하에서 터널이 시공될 경우, 터널막장은 시공 중 용수에 의한 영향을 받게 되며 지보 System은 시공 후 지하수 흐름이 문제시 될 수 있다. 븐 연구는 터널막장 및 라이닝에 대해 배수조건에 따른 지하수 흐름을 고려한 적절한 해석 및 설계방법을 제시하고 있다. 첫째, 시공 완료 후 터널 라이닝이 배수조건에 따라 받게 되는 지하수의 영향을 라이닝에 작용하는 응력 및 변위로써 검토하고 각 배수조건별로 적절한 해석 및 설계방법을 제안하였다. 둘째, 시공 중 굴착에 의한 지하수의 흐름이 문제가 되는 터널막장에 대해서 지하수 흐름의 3 차원 모델링을 수행하였으며, 그 결과를 막장의 안정성에 대한 이론적 검토에 반영하여 침투력이 막장의 안정성 미치는 영향에 대하여 검토하였다. 또한 막장을 평면변형률조건(plane strain condition)으로 모델링 하여 침투 고려시와 지하수의 미고려시에 막장면에 작용하는 응력 및 변위를 산출하여 침투력이 막장면에 미치는 영향을 검토하였다. 연구 결과, 터널 시공시 지하굴착 및 배수에 의한 지하수위 저하가 크지 않은 구간에서 터널막 장은 정상류 흐름에 의한 침투력을 고려한 적절한 안정대책이 요구되며, 터널 라이닝 또한 침투를 고려한 배수개념이나 내부 라이닝이 정수압을 견뎌야 하는 비배수 개념에 의해 설계 및 시공이 이루어져야 함을 보여준다.

• 한국강구조학회 논문집
• /
• 제23권4호
• /
• pp.429-438
• /
• 2011

본 연구에서는 스트라치 시스템의 긴장설치과정 및 극한하중 해석을 수행하기 위한 명시적 해석법을 제안하였다. 스트라치 시스템은 Stressed-Arch에서 유래한 용어로 슬리브와 갭이 도입된 유동하현재 내부의 긴장재에 초기장력을 도입함으로써 갭이 점차 닫히게 되며, 이에 따라 상현재에 곡률이 도입되면서 전체 구조물이 상승하여, 최종적인 아치형태의 구조물을 형성하는 독창적인 구조시스템이다. 스트라치 시스템의 초기장력 도입과정을 긴장설치(stress-erection) 과정이라 하며, 초기곡률의 도입에 따라 유동 상현재에는 과도한 초기변형이 발생하여 소성거동에 의한 강체회전이 발생하는 불안정 구조물이 된다. 본 연구에서는 이러한 스트라치 시스템의 불안정 거동특성을 해석하기 위해서 강성행렬을 사용하지 않는 명시적 동적이완법을 사용하여 비선형 평형방정식의 해를 구하였고, 대변위 및 단면의 재료적 특성을 반영할 수 있는 필라맨트 보요소를 사용하여 연속된 상현재의 비선형 거동특성을 분석하였다. 필라맨트 보요소의 단면은 다수의 1차원 필라맨트로 구성되며, 각각의 필라맨트에 대해서 다양한 재료모델을 적용할 수 있다. 본 연구에서는 비선형 재료모델인 Ramberg-Osgood모델 및 Bi-linear 탄소성 모델을 적용하여 긴장설치 및 극한하중 해석을 수행하였고, 그 결과를 이전의 실험적 연구결과와 비교 분석하였다. 본 연구의 해석결과는 이전의 실험적 연구결과와 유사하였으며, 명시적 해석법의 특성상 효율적으로 후좌굴거동 특성까지 해석할 수 있었다.

• 대한치과교정학회지
• /
• 제34권1호
• /
• pp.93-107
• /
• 2004

Enlow의 구조적 대응체 분석법은 다수의 평균에 의한 정상치와 비교하는 분석법과는 달리 개개인에서 두 개안면의 구성양상이 발달되어온 해부학적, 발육학적 특징의 복합체를 설명하는 방법이다. 본 연구에서는 악교정 수술을 시행한 골격성 III급 부정교합 환자 80명(남자 40명, 여자 40명)을 대상으로 구조적 대응체 분석법을 이용하여 수술 전후를 비교하고, 정상교합자와 비교하여 다음과 같은 결과를 얻었다. 1. 수술에 의해서 주로 하악지 전후방 길이(B3)가 감소되어 상하악의 전후방적 부조화가 해소되었다. 2. 수술에 의해서 하악지의 배열(R3, R4)이 후방으로 재배열되었으며 교합평면(R5)은 후하방으로 회전되었다. 3. 수술후와 정상교합자의 비교결과 수술후에 골격성 III급 양상은 해소되었으나 두개저수준(R1, R2)에서는 차이를 보이고 있다. 4. 수술방법에 따른 비교결과 양악수술은 상악골의 길이(A4), 하악지의 배열(R3, R4), 교합평면의 각도(R5)에서 하악수술과 차이를 보였으나 두 그룹간의 차이는 크지 않았다. 한국인 골격성 III급 부정교합자의 특징을 살펴본 지난 연구에서 골격적 원인은 두개저의 후상방회전, 상악골의 후방위치, 하악지의 전방경사, 하악골체부의 길이증가였는데 본 연구 결과 실제 수술에서 개선되는 요소는 주로 하악지의 전후방 길이 감소, 상악골의 전진이기 때문에 상대적인 개선으로 골격성 III급 양상이 해소되었다고 볼 수 있다. 이러한 결과는 실제적인 수술부위가 상악골과 하악골의 일부에 국한되기 때문이다.

• 대한의용생체공학회:의공학회지
• /
• 제28권3호
• /
• pp.397-409
• /
• 2007

Owing to needs of biomechanical comprehension and analysis to obtain various medical treatment designs which are related with the spine in order to cure and diagnose LBP patients, the FE modeling and nonlinear analysis of lumbosacrum including a partial ilium and iliolumbar ligaments, were carried out. First, we investigated whether the geometrical configuration of vertebrae displayed by DICOM slice files is regular and normal condition. After constructing spinal vertebrae including a partial ilium, a sacrum and five lumbars (from L1 to L5)with anatomical shape reconstructed using softwares such as image modeler and CAD modeler, we added iliolumbar ligaments, lumbar ligaments, discs and facet joints, etc.. And also, we assigned material property and discretized the model using proper finite element types, thus it was completely modeled through the above procedure. For the verification of each segment, average sagittal ROM, average coronal ROM and average transversal ROM under various loading conditions(${\pm}10Nm$), average vertical displacement under compression(400N), ALL(Anterior Longitudinal Ligament) and PLL(Posterior Longitudinal Ligament) force at L12 level, strains of seven ligaments on sagittal plane at L45 level and maximal strain of disc fibers according to various loading conditions at L45 level, etc., they were compared with experimental results. For the verification of multilevel-lumbosacrum spine including partial ilium and iliolumbar ligaments, the cases with and without iliolumbar ligaments were compared with ROM of experiment. The results were obtained from analysis of the verified FE model as follows: I) Iliolumbar ligaments played a stabilizing role as mainly posterior iliolumbar ligaments under flexion and as both posterior and anterior iliolumbar ligaments of one side under lateral bending. 2) The iliolumbar ligaments decreased total ROM of 1-8% in total model according to various motion conditions, which changed facet contact forces of L5S level by approximately 0.8-1.4 times and disc forces of L5S level by approximately 0.8-1.5 times more than casewithout ilioligaments, under various loading conditions. 3) The force of lower discs such as L45 and L5S was bigger than upper discs under flexion, left and right bending and left and right twisting, except extension. 4) It was predicted that strains of posterior ligaments among iliolumbar ligaments would produce the maximum 16% under flexion and the maximum 10% under twisting. 5) It's expected that this present model applies to the development and design of artificial disc, since it was comparatively in agreement with the experimental datum.

• Wind and Structures
• /
• 제11권4호
• /
• pp.275-289
• /
• 2008

Large eddy simulations have been performed within and over different types of urban building arrays. This paper adopted three dimensionless parameters, building frontal area density (${\lambda}_f$) the variation degree of building height (${\sigma}_h$), and the staggered degree of building range ($r_s$), to study the systematic influence of building spacing, height and layout on wind and turbulent characteristics. The following results have been achieved: (1) As ${\lambda}_f$ decrease from 0.25 to 0.18, the mean flow patterns transfer from "skimming" flow to "wake interference" flow, and as ${\lambda}_f$ decrease from 0.06 to 0.04, the mean flow patterns transfer from "wake interference" flow to "isolated roughness" flow. With increasing ${\lambda}_f$, wind velocity within arrays increases, and the vortexes in front of low buildings would break, even disappear, whereas the vortexes in front of tall buildings would strengthen and expand. Tall buildings have greater disturbance on wind than low buildings do. (2) All the wind velocity profiles and the upstream profile converge at the height of 2.5H approximately. The decay of wind velocity within the building canopy was in positive correlation with ${\lambda}_f$ and $r_s$. If the height of building arrays is variable, Macdonald's wind velocity model should be modified through introducing ${\sigma}_h$, because wind velocity decreases at the upper layers of the canopy and increases at the lower layers of the canopy. (3) The maximum of turbulence kinetic energy (TKE) always locates at 1.2 times as high as the buildings. TKE within the canopy decreases with increasing ${\lambda}_f$ and $r_s$ but the maximum of TKE are very close though ${\sigma}_h$ varies. (4) Wind velocity profile follows the logarithmic law approximately above the building canopy. The Zero-plane displacement $z_d$ heighten with increasing ${\lambda}_f$, whereas the maximum of and Roughness length $z_0$ occurs when ${\lambda}_f$ is about 0.14. $z_d$ and $z_0$ heighten linearly with ${\sigma}_h$ and $r_s$, If ${\sigma}_h$ is large enough, $z_d$ may become higher than the average height of buildings.

• 대한토목학회논문집
• /
• 제10권2호
• /
• pp.23-38
• /
• 1990

본(本) 연구(硏究)에서는 다재하조건(多載荷條件), 허용응력(許容應力), 좌굴응력(座屈應力), 변위(變位), 고유진동수(固有振動數) 제약(制約)을 고려(考慮)한 트러스 구조물(構造物)의 형상(形狀)을 효율적(效率的)으로 최적화(最適化)하기 위해서 Two-Levels 분할(分割) 최적화(最適化) 기법(技法)을 택(擇)하고 Level 1에서는 허용방향법(許容方向法)에 의(依)한 단면(斷面) 최적화(最適化), Level 2에서 Powell 1법(法)의 일방향(一方向) 탐사법(探査法)에 의(依)한 목적함수(目的函數)만이 최소(最小)가 되도록 형상(形狀)을 최적화(最適化)하였다. 본 연구(硏究)의 알고리즘을 트러스의 구조모형(構造模型)에 적용(適用)하여 얻어진 연구(硏究) 결과(結果)를 요약(要約)하면 다음과 같다. 1. 본(本) 연구(硏究)의 알고리즘은 트러스의 형태(形態), 재하조건(載荷條件), 정적제약조건(靜的制約條件), 고유진동수(固有振動數) 제약조건(制約條件) 등(等)에 구애받지 않고 효율적(效率的)으로 수렴(收斂)함을 알 수 있다. 2. 트러스 구조물(構造物)의 최적형상(最適形狀)은 고려(考慮)되는 제약조건식(制約條件式)에 따라 대단(大端)히 상이(相異)함을 알 수 있다. 3. 동일(同一)한 설계조건하(設計條件下)에서 트러스의 기하학적형태(幾何學的形態)를 고정(固定)시키고 단면(斷面)만을 최적화(最適化)한 경우 보다 본(本) 연구(硏究)의 알고리즘에 의(依)하여 트러스의 형상(形狀)까지도 최적화(最適化)한 경우에는 트러스의 초기(初期)의 기하형태(幾何形態)와 설계조건(設計條件)에 따라 다소(多少) 차이(差異)가 있겠지만 중량(重量)을 상당(相當)히 감소(減少)시킬 수 있다는 사실(事實)을 알 수 있었다.

• 한국암반공학회:학술대회논문집
• /
• 한국암반공학회 2000년도 암반공학문제의 수치해석(Numerical Analysis in Rock Engineering Problems)
• /
• pp.67-75
• /
• 2000

암석의 파괴인성 측정법으로서 최근에 연구되고 있는 Brazilian 파괴인성 시험을 PFC code를 이용하여 수치해석적으로 시뮬레이션하였다. 이로부터 가압면 형태, 시험편과 모델링 입자의 크기 관계, 가압면 각도 범위, 하중속도 등의 영향 변수에 대한 분석 및 적정 기준에 대해 검 토하였다. 균등한 하중 전달을 위해 가압면이 평면인 Brazilian 시험편을 도입하였는데, 시험편 중앙에서 초기 인장균열을 생성하고 파괴인성 산정에 필요한 안정한 인장균열 진전 시점의 임계하중이 나타나는 하중-변위 곡선을 얻을 수 있는 가압면 각도는 $20^{\circ}C$ -$40^{\circ}C$ 범위로 나타났다. 가압면 각도가 $20^{\circ}C$ 인 경우에, 디스크 시험편의 구성입자 반경이 1 mm 이하이고 하중속도가 0.01 mm/s 이하인 조건에서는 거의 일정한 파괴인성값을 얻을 수가 있었다. 가압면 각도가 $20^{\circ}C$이상이고 하중속도가 0.01 mm/s 이하는 파괴인성 측정에서 중요한 기본 조건인 시험편 중 에서의 인장균열 생성 및 이 인장균열의 안정적인 진전 제어를 위한 조건인 것으로 볼 수가 있다.

• 터널과지하공간
• /
• 제10권3호
• /
• pp.320-328
• /
• 2000

암석의 파괴인성 측정법으로서 최근에 연구되고 있는 Brazilian 파괴인성 시험을 PFC code를 이용하여 수치 해석적으로 시뮬레이션하였다. 이로부터 가압면 형태, 시험편과 모델링 입자의 크기 관계, 가압면 각도 범위, 하중속도 등의 영향 변수에 대한 분석 및 적정 기준에 대해 검토하였다. 균등한 하중 전달을 위해 가압면이 평면인 Brazilian 시험편을 도입하였는데, 시험편 중앙에서 초기 인장균열을 생성하고 파괴인성 산정에 필요한 안정한 인장균열 진전 시점의 임계하중이 나타나는 하중-변위 곡선을 얻을 수 있는 가압면 각도는 20$^{\circ}$~40$^{\circ}$범위로 나타났다. 가압면 각도가 20$^{\circ}$인 경우에, 디스크 시험편의 구성입자 반경이 1mm 이하이고 하중속도가 0.01mm/s 이하인 조건에서는 거의 일정한 파괴인성값을 얻을 수가 있었다. 가압면 각도가 20$^{\circ}$이상이고 하중속도가 0.01mm/s 이하는 파괴인성 측정에서 중요한 기본 조건인 시험편 중앙에서의 인장균열 생성 및 이 인장균열의 안정적인 진전 제어를 위한 조건인 것으로 볼 수가 있다.

• 한국전산구조공학회논문집
• /
• 제25권5호
• /
• pp.413-420
• /
• 2012

본 논문은 2차원 선형탄성 직접 경계요소법에서 저매개변수 요소를 사용할 때 Kernel의 적분방법에 대하여 논의하였다. 일반적으로 등매개변수 요소의 경우 형상함수로 통칭되는 해의 기저함수와 요소의 적분을 위해 사용되는 사상함수를 동일하게 사용한다. 그러나 본 논문에서는 사상함수의 차수를 낮게 취하여 순수기저절점을 도입하고 그때 직접 경계요소의 Kernel을 적분하기 위한 방법이 모색되었다. 일반적으로 경계요소법의 적분 Kernel의 경우 Log수치적분과 코쉬주치(Cauchy principal value) 등을 통해 해결하는데, 본 논문에서는 대수적 조작을 통해 적분값의 정확도를 높일 수 있도록 새로운 수식을 유도하였다. 본 연구에서 저매개변수 기반의 직접 경계요소에 대한 강건성과 정확도를 검증하기 위해 2차원 타원형 편미분방정식으로 표현되는 평면응력과 평면변형문제에 대해 적용하였다. 적용 예제로는 단순연결영역(simple connected region)의 대표적 문제인 캔틸레버보와 다중연결영역(multiple connected region)의 대표적인 문제인 개구부가 있는 사각평면에 대해 각각 수치해석을 수행한 결과 대폭적인 자유도의 감소에 비해 정확도 측면에는 기존의 방법과 차이가 없음을 볼 수 있었다. 본 논문에서 제시된 방법은 기저함수 고차화 저매개변수 직접 경계요소법(subparametric high order boundary element)과 이에 기초를 둔 저매개변수 고차 이중경계요소법(subparametric high order dual boundary element)의 초석이 될 수 있을 것이다.

• 한국지구과학회지
• /
• 제43권1호
• /
• pp.151-164
• /
• 2022

역사지진과 계기지진 기록에 따르면 한반도 남동부는 우리나라에서 지진활성도가 가장 높게 평가되는 곳으로, 최근에 양산단층대와 울산단층대를 따라 제4기 단층이 다수 보고되어 고지진학적 연구가 활발하게 이루어지고 있다. 특히 울산단층대의 중부지역에 해당하는 경북 경주시 외동읍 말방리 일원은 울산단층대 내에서 가장 많은 활성단층이 보고된 지역이다. 따라서 이 지역에 대한 고지진학적 특성을 이해하기 위하여 먼저 LiDAR 영상 및 항공사진을 이용한 지형 및 선형구조 분석을 실시하여 단층에 의한 기복으로 추정되는 지형인자를 확인하고, 야외답사와 물리탐사를 통해 단층을 추적하여 기 보고된 말방단층 지점에서 약 300 m 북서쪽에 위치한 곳에서 길이 20 m, 너비 5 m, 깊이 5 m의 굴착조사를 실시하였다. 굴착단면을 통해 분석된 제4기 퇴적층의 특징을 바탕으로 단층의 기하학적·운동학적 특성을 해석하여 고지진학적 특성을 규명하고자 하였다. 이번 굴착단면에서 확인된 역단층의 기하를 보이는 단층의 자세는 N26°W/33°NE로 울산단층대를 따라 분포하는 기 보고된 단층들과 유사하다. 약 40 cm의 단일 겉보기 변위가 인지되었으나 단층조선의 부재로 실변위는 산출할 수 없었다. 선행연구에서 제안된 극저온구조층의 연대결과 값을 토대로 단층의 최후기 운동시기는 후기 뷔름빙기 이전으로 추정하였다. 기 보고된 연구결과와 본 굴착단면에서 획득한 단층기하를 종합하여 이 지역에 발달하는 단층계를 인편상구조로 해석하였고, 단층특성을 반영한 모델을 제시하였다. 말방리 일원에서 수 회의 굴착조사를 비롯한 다수의 선행연구가 수행되었음에도 불구하고 구체적인 단층변수에 대한 정보가 미진하고 각 지점들 간의 상관성이 명확하게 규명되지 않은 것은 역단층의 복잡한 운동학적 특성에 기인한 것으로 판단된다. 추후 고지진학적 연구가 추가적으로 수행된다면 상기의 문제점들을 해결하여 종합적인 단층의 형태와 운동사가 규명될 수 있을 것으로 판단된다.