• 제목/요약/키워드: 2-dimensional finite element method

검색결과 901건 처리시간 0.029초

전기비저항 모니터링 방법을 이용한 충적층 지하수위 변동 감지 (Detection of Groundwater Table Changes in Alluvium Using Electrical Resistivity Monitoring Method)

  • 김형수
    • 지질공학
    • /
    • 제7권2호
    • /
    • pp.139-149
    • /
    • 1997
  • 전기비저항 모니터링을 이용하여 충적층 지하수위의 변동을 감지 할 수 있는지의 여부가 연구되었다. 연구는 유한요소법을 이용한 전기비저항 수치 모델링 시험과 20일간의 장기 양수 시험이 수행된 금강 정동리 지역의 현장 자료를 이용하여 수행되었다. 수치 모델링 시험 및 현장 자료 획득은 전극 기본 간격이 4m인 21개의 고정전극을 이용하였으며, 전극배열은 "변형 웨너 배열"과 쌍극자 배열을 이용하였다. 현장 장기 양수 시험에서 관측된 수위 변화를 단순화하여 지하 2차원 비저항 모델로 사용한 사용한 결과, 겉보기 전기비저항의 변화 양상은 수위 변화보다 둔감한 변화 양상을 보여 주며, 수위 변화가 있는 부분뿐만 아니라하부 결정질 기반암의 경계 부위에 해당하는 심도에서도 겉보기 비저항 장치가 변화하는 현상이 관측되었다. 이러한 수치 결과와 부합되는 결과가 실제 현장에서 획득된 겉보기 전기 비저항 단면에서도 나타났다. 수치 계산결과 및 현장 자료는 전기 비저한 모니터링이 지하수위 변동을 유무 감지하는데 유용하게 사용될 수 있음을 보여주었다. 그러나 겉보기 전기비저항 단면이 수위 변동량을 정량적으로 규명할 만큼의 뚜렷한 변화를 보여주지는 못하기 때문에 전기비저항 방법만으로 정략적인 수위변동을 감지하는데 한계가 있는 것으로 판단되었다. 이러한 지하 수위 변동량 등을 명확히 규명하기 위해서는 고해상도 탄성파 또는 지하레이다 탐사처럼 지하의 모습을 영상화 할 수 있는 지구물리 탐사법이 병행되는 것이 바람직하다.

  • PDF

Umbrella Arch 공법이 적용된 터널의 3차원 유한요소 해석에 관한 연구 (A Study on the Three Dimensional Finite Element Analysis for the Tunnel Reinforced by Umbrella Arch Method)

  • 김창용;배규진;문현구;최용기
    • 터널과지하공간
    • /
    • 제8권3호
    • /
    • pp.209-225
    • /
    • 1998
  • 최근터널 굴착 보조공법중의 하나인 강관 보강형 다단 그라우팀 공법(Umbrella Arch Method, UAM)은 지반을 보강하고 터널 막장의 안정성을 증진시키기 위해서 많은 현장에서 사용되고 있다. 이러한 UAM은 터널 보강목적의 forepoling과 차수목적으 grouting이 한 공정으로 구성되어 있다는 잇점 때문에, 최근 국내 지하철, 도로터널 및 전력구터널 등에서 많은 적용 사례를 찾아 볼 수 있다. 그러나 이 공법은 주로 현장 시공을 통해서 얻어진 경험적인 방법에 의해서 설계와 시공이 이루어지고 있기 때문에 본 공법에 대한 보다 정량적이고 체계적인 설계인자 평가 작업이 필요하다. 따라서, 본 연구에서는 수치해석 방법에 의한 체계적이고 정량적인 효과확인 과정을 제안하였고, 몇몇 설계인자에 대한 매개변수 변환연구를 수행하였다. 이를 위해서 먼저, UAM의 지반보강기구에 있어서 관련된 강관, 그라우트재 및 강지보재등의 역할을 밝히고자 하였고, 두 번째로 매개변수 변환연구를 통해 UAM의 설계 제요소들에 대한 영향을 평가하기 위해 1)지반조건별, 2) 토피고별, 3) 강관배치형상별, 4) 그라우트 영역별, 5)강관자체 특성별 해석을 수행하여 각 항목별로 상호 비교.분석하였다.

  • PDF

수평 자기쌍극자 송신원을 이용한 2.5차원 전자탐사 모델링 (2.5 Dimensional EM Modeling considering Horizontal Magnetic Dipole Source)

  • 권형석;송윤호;손정술;서정희
    • 지구물리와물리탐사
    • /
    • 제5권2호
    • /
    • pp.84-92
    • /
    • 2002
  • 이 연구에서는 최근에 개발되어 적용성에 대한 시험이 진행되고 있는 $1\;kHz\~1\;MHz$ 주파수대역의 수평 자기쌍극자를 송신원으로 하여 전기장과 자기장을 측정하는 탐사법의 모델링 알고리듬을 개발하고자 하였다. 이를 위해 공간주파수 영역에서 수평 자기쌍극자를 송신원으로 하는 이차장 Maxwell 방정식을 구성하고 유한요소법을 이용하여 2.5차원 모델링 알고리듬을 구현하였다. 또한, 검증된 층서구조 반응결과와 비교하여 그 정확도를 검증하였다. 개발된 알고리듬을 이용하여 균질 반무한공간 내에 전도성 및 비전도성 고립이상체가 존재하는 2차원 모형에 대한 전기장, 자기장과 임피던스를 계산하였다. 계산된 전기장, 자기장과 임피던스의 거동에 대한 고찰 결과 고립이상체와 같이 2차원 구조에서 전기장과 임피던스가 자기장에 비해 이상체의 탐지에 효과적임을 보였다. 특히, 전기장을 자기장으로 나눈 값인 임피던스는 전기장의 공간 분해능을 가지면서 보다 안정된 결과를 보여 주었다. 따라서 자기장만을 측정하는 기존의 전자탐사법에 비해 전기장과 자기장을 모두 측정하는 탐사법이 천부에 존재하는 이상체의 탐지에 보다 효과적임을 알 수 있었다.

수평 아치형 터널 갱구부 비탈면의 역학적 거동 (Mechanical Behavior of Tunnel Portal in Horizontal Arch Slope)

  • 양문상;이상덕
    • 한국터널지하공간학회 논문집
    • /
    • 제2권2호
    • /
    • pp.50-61
    • /
    • 2000
  • 터널의 갱구부는 터널에서 3차원 거동을 가장 뚜렷하게 보이는 부분이며 대체로 풍화가 진행되어 거의 토사화되어 있는 취약한 암반에 시공됨에 따라 비탈면은 잠재적인 활동체가 되게 된다. 터널 갱구부 비탈면에서 안정에 위해되는 요인이 발생했을 경우 파괴가 급속히 진행되어 대형사고로 이어질수 있으므로 특별한 안정화 대책이 필요하여 경제적으로 큰 부담이 된다. 이를 해결하기 위한 방안으로 갱구부 비탈면을 수평 아치형으로 굴착하게 되면 굴착량이 줄어서 경제적이고 자연훼손이 줄어들며 보다 안전하고 아름다운 형상을 기대할 수 있다. 그러나 인식부족으로 인하여 이와 같은 시도는 거의 이루어지지 않고 있어서 면벽식 갱구부가 보편화 되어 있는 실정이다. 본 연구에서는 아치형 갱구부 비탈면의 이점을 알아보기 위하여 3차원 수치해석을 ㄷ형과 아치형 비탈면에 대하여 각각 수행하여 그 결과를 서로 비교하였으며 터널 굴진과 암반강성에 따른 갱구부 비탈면의 안정성은 비탈면에서 발생되는 최대전단변형율과 막장면에서 터널굴진에 따른 소성접근도를 분석하여 검토하였다.

  • PDF

2차원 전기비저항 탐사를 이용한 경주 개곡리 지역의 단층조사 (Investigation of fault in the Kyungju Kaekok-ri area by 2-D Electrical Resistivity Survey)

  • 이치섭;김희준;공영세;이정모;장태우
    • 지구물리와물리탐사
    • /
    • 제4권4호
    • /
    • pp.124-132
    • /
    • 2001
  • 경주 개곡리 지역에서 단층 파쇄대와 지하 내부구조를 파악하기 위해 전기비저항 탐사를 실시하였다. 지층의 전기전도도 분포는 공극수의 함량, 지하수, 공극률 및 점토광물 등의 많은 요인에 의해 지배되고, 이러한 특징은 지층의 풍화상태와 정도, 단층 파쇄대를 간접적으로 설명하는 증거가 된다. 따라서 전기비저항 탐사는 지표에 드러나지 않아 그 연장성이 확실하지 않은 단층 파쇄대를 탐지하기에는 좋은 방법이다. 단층구조를 해석하기 위한 2차원 전기비저항 역산 기법으로 평활화제약 최소자승법을 이용하였고, 모형반응의 계산은 지형의 기복을 고려하기 위해 유한요소법을 사용하였다. 본 연구에 이용된 야외자료는 울산단층에 인접한 경주 외동읍 개곡리 부근의 제 4기 단층 파쇄대로 추정되는 곳에서 쌍극자배열을 이용하여 획득하였다. 이들 자료를 역산하여 전기비저항 2차원 단면도를 작성한 결과, 단층 파쇄대의 특징을 잘 보여주는 것으로 나타났다. 본 연구의 결과, 경주시 외동읍 개곡리 제 4기 단층은 노두에서부터 총 50m까지 단층의 연장이 확인되었으며 그 연장 방향은 $N80^{\circ}W$임이 밝혀졌다. 본 연구를 통해 단층 조사에 있어 2차원 전기비저항 탐사해석이 유용하게 적용될 수 있음을 보여주었다.

  • PDF

응력특이를 갖는 축방향 부재의 웨이블렛 급수해석 (Wavelet Series Analysis of Axial Members with Stress Singularities)

  • 우광성;장영민;이동우;이상윤
    • 한국전산구조공학회논문집
    • /
    • 제23권1호
    • /
    • pp.1-8
    • /
    • 2010
  • 푸리에 급수는 사인 곡선처럼 일정한 진폭으로 진동하는 정규파(wave)를 사용한다. 그래서 푸리에 급수에서 사용하는 함수는 진동수의 크기가 시간에 따라 변하지 않기 때문에 국부적인 영역에서 급작스런 진동이나 불연속성을 갖는 신호를 표현하기에는 한계가 있다. 그러나 이러한 푸리에 해석의 단점을 여러개의 적절한 웨이블렛의 선형조합에 의해 보완할 수 있는 것이 웨이블렛 급수해석이다. 시간에 집중되어진 궤적의 작은 잔파(wavelet)를 사용함으로써 시간과 주기의 폭을 변화시킬 수 있기 때문에 유동적이고, 특이(singular)형상을 지닌 신호들을 보다 효율적으로 표현할 수 있다. 이 연구의 주요 목적은 웨이블렛 급수해석이라고 불리는 방법을 2계 편미분방정식으로 표현되는 1차원 축방향 부재에 웨이블렛 이론을 적용함과 동시에 유한요소법과 같은 수치해석법과의 비교를 통해 성능평가를 위해 제안되었다. 여러 형태의 웨이블렛 함수의 검토 후에 HAT 함수가 웨이블렛 및 스케일링 함수로 채택되었다. 등분포하중을 받는 경우의 축방향 부재해석에서 제안된 방법은 유한요소법과 같이 효율적임을 보이며, 특히 응력특이점에서는 더 정확한 값을 보였으며, 계산시간도 절약되는 장점을 얻을 수 있었다.

H-Pile과 Plastic Sheet Pile을 결합한 토류벽체에 대한 수치해석적 연구 (Numerical study for Application of H-Pile Connection Plastic Sheet Pile Retaining Wall (HCS))

  • 이규남;임희대
    • 지질공학
    • /
    • 제27권3호
    • /
    • pp.331-343
    • /
    • 2017
  • 본 연구에서는 기존의 H-Pile+토류판 또는 H-Pile+토류판+차수그라우팅 공법의 안정성, 시공성 및 경제성을 개선하기 위해 H-Pile에 Plastic Sheet Pile(P.S.P)과 연성벽체인 P.S.P의 간격유지 및 보강기능을 위한 간격재(각형강관)를 결합한 토류벽체 System인 HCS공법을 개발하고, HCS공법을 구성하는 각 부재의 거동을 3차원 유한요소해석에 의해 규명하는 연구이다. HCS공법의 거동을 수치해석적으로 규명하기 위해 Plastic Sheet Pile 규격 3종류, H-Pile 규격 2종류 및 설치간격 3종류, 간격재 규격 1종류 및 설치간격 4종류에 대해 광범위한 3차원 유한요소해석을 실시하였다. 수치해석결과 $P.S.P-460{\times}131.5{\times}7t$ (PS7)와 H-Pile $250{\times}250{\times}9{\times}14$ (H250), $P.S.P473{\times}133.5{\times}9t$ (PS9)와 H-Pile $300{\times}200{\times}9{\times}14$ (H300)의 조합에서 상대적으로 유사한 응력비(=발생응력/허용응력)를 갖는 것으로 검토되어 이 제품의 조합이 경제적인 것으로 확인되었으며, P.S.P+H-Pile+간격재 복합체의 강성이 증가할수록 벽체의 수평변위와 상부지반의 연직변위가 감소하였다. 특히, H-Pile과 P.S.P의 강성차이로 인한 Arching 현상으로 P.S.P의 토압의 상당부분이 H-Pile로 응력(토압) 전이가 발생하여 P.S.P의 응력 및 변위는 미소하게 나타났다. 본 연구를 통하여 HCS공법을 구성하는 각각의 부재들의 거동을 확인할 수 있었으며, 확인된 연구결과를 통해 향후 HCS공법을 합리적이고 안정하며 경제적으로 적용하는 데 활용 가능하리라 판단된다.

T300/924C 탄소섬유/에폭시 복합재 적층판의 이차원 압축 강도의 크기효과 및 좌굴방지장치의 영향 (Two Dimensional Size Effect on the Compressive Strength of T300/924C Carbon/Epoxy Composite Plates Considering Influence of an Anti-buckling Device)

  • 공창덕;방조혁;이정환
    • 한국복합재료학회:학술대회논문집
    • /
    • 한국복합재료학회 2002년도 추계학술발표대회 논문집
    • /
    • pp.88-91
    • /
    • 2002
  • The two dimensional size effect of specimen gauge section (length x width) was investigated on the compressive behavior of a T300/924 [45/-45/0/90]3s, carbon fiber-epoxy laminate. A modified ICSTM compression test fixture was used together with an anti-buckling device to test 3mm thick specimens with a 30$\times$30, 50$\times$50, 70$\times$70, and 90mm$\times$90mm gauge length by width section. In all cases failure was sudden and occurred mainly within the gauge length. Post failure examination suggests that $0^{\circ}$ fiber microbuckling is the critical damage mechanism that causes final failure. This is the matrix dominated failure mode and its triggering depends very much on initial fiber waviness. It is suggested that manufacturing process and quality may play a significant role in determining the compressive strength. When the anti-buckling device was used on specimens, it was showed that the compressive strength with the device was slightly greater than that without the device due to surface friction between the specimen and the device by pretoque in bolts of the device. In the analysis result on influence of the anti-buckling device using the finite element method, it was found that the compressive strength with the anti-buckling device by loaded bolts was about 7% higher than actual compressive strength. Additionally, compressive tests on specimen with an open hole were performed. The local stress concentration arising from the hole dominates the strength of the laminate rather than the stresses in the bulk of the material. It is observed that the remote failure stress decreases with increasing hole size and specimen width but is generally well above the value one might predict from the elastic stress concentration factor. This suggests that the material is not ideally brittle and some stress relief occurs around the hole. X-ray radiography reveals that damage in the form of fiber microbuckling and delamination initiates at the edge of the hole at approximately 80% of the failure load and extends stably under increasing load before becoming unstable at a critical length of 2-3mm (depends on specimen geometry). This damage growth and failure are analysed by a linear cohesive zone model. Using the independently measured laminate parameters of unnotched compressive strength and in-plane fracture toughness the model predicts successfully the notched strength as a function of hole size and width.

  • PDF

블록 이동용 러그의 안전사용하중 결정에 관한 구조 평가법 (Structural Evaluation Method to Determination Safe Working Load of Block Handling Lugs)

  • 권오현;박주신;서정관
    • 해양환경안전학회지
    • /
    • 제29권4호
    • /
    • pp.363-371
    • /
    • 2023
  • 한 척의 선박을 건조하기 위해서는 다양한 크기의 블록(block)들을 이동 및 탑재해야 한다. 이러한 과정에서 블록의 체결 방법 및 각 조선소 설비 특성에 맞는 다양한 기능에 부합하는 러그를 사용하고 있다. 블록 구조의 중량 및 형태에 따라서 러그의 크기와 형상이 다양하며, 샤클(shackle)이 체결되는 홀 주변에 부족한 강성을 보완하기 위하여 덧판(doubling pad)을 용접하여 구조를 보강한다. 리프팅(lifting) 조건별 러그의 설계를 하는 방법은 보 이론(beam theory)에 의한 수계산 방법과 유한요소해석 모델링을 이용한 구조해석을 수행하고 있다. 해석적 방법의 경우, 요소의 종류와 모델링 방법에 따라서 결과 차이가 발생하여 표준화된 평가법의 정립이 필요한 상황이다. 이러한 모호한 방법론 적용 시 블록의 이동 및 반전(turn-over) 과정 중에서 심각한 안전 문제를 유발할 가능성이 있다. 본 연구에서는 러그의 실제 탑재공정에 따른 구조 응답을 평가할 수 있는 모델링 조건, 평가법을 확정하고자 다양한 변수의 영향을 수치 구조해석을 통하여 비교 및 분석하였다. 러그 홀(hole) 주변 덧판부와 용접 비드(bead)를 표현한 모델링 기법이 가장 실제적인 거동 결과를 주고 있다. 실제 러그와 동일한 조건(용접부 비드만 주재료와 연결)의 모델링에 등가하중을 적용한 결과는 MPC 하중 적용 결과보다 낮은 최종강도를 나타낸다. 더불어 해석 시간 단축을 위해서 2차원 쉘(shell) 요소를 적용한 경우, 덧판 두께를 85% 수준으로 감소시켜서 안전사용하중을 예측할 수 있음을 확인하였다. 논문에서 검토한 다양한 변수의 영향들 결과는 러그 설계 및 안전사용하중 예측에 근거 자료로 활용될 것으로 기대된다.

Memory Organization for a Fuzzy Controller.

  • Jee, K.D.S.;Poluzzi, R.;Russo, B.
    • 한국지능시스템학회:학술대회논문집
    • /
    • 한국퍼지및지능시스템학회 1993년도 Fifth International Fuzzy Systems Association World Congress 93
    • /
    • pp.1041-1043
    • /
    • 1993
  • Fuzzy logic based Control Theory has gained much interest in the industrial world, thanks to its ability to formalize and solve in a very natural way many problems that are very difficult to quantify at an analytical level. This paper shows a solution for treating membership function inside hardware circuits. The proposed hardware structure optimizes the memoried size by using particular form of the vectorial representation. The process of memorizing fuzzy sets, i.e. their membership function, has always been one of the more problematic issues for the hardware implementation, due to the quite large memory space that is needed. To simplify such an implementation, it is commonly [1,2,8,9,10,11] used to limit the membership functions either to those having triangular or trapezoidal shape, or pre-definite shape. These kinds of functions are able to cover a large spectrum of applications with a limited usage of memory, since they can be memorized by specifying very few parameters ( ight, base, critical points, etc.). This however results in a loss of computational power due to computation on the medium points. A solution to this problem is obtained by discretizing the universe of discourse U, i.e. by fixing a finite number of points and memorizing the value of the membership functions on such points [3,10,14,15]. Such a solution provides a satisfying computational speed, a very high precision of definitions and gives the users the opportunity to choose membership functions of any shape. However, a significant memory waste can as well be registered. It is indeed possible that for each of the given fuzzy sets many elements of the universe of discourse have a membership value equal to zero. It has also been noticed that almost in all cases common points among fuzzy sets, i.e. points with non null membership values are very few. More specifically, in many applications, for each element u of U, there exists at most three fuzzy sets for which the membership value is ot null [3,5,6,7,12,13]. Our proposal is based on such hypotheses. Moreover, we use a technique that even though it does not restrict the shapes of membership functions, it reduces strongly the computational time for the membership values and optimizes the function memorization. In figure 1 it is represented a term set whose characteristics are common for fuzzy controllers and to which we will refer in the following. The above term set has a universe of discourse with 128 elements (so to have a good resolution), 8 fuzzy sets that describe the term set, 32 levels of discretization for the membership values. Clearly, the number of bits necessary for the given specifications are 5 for 32 truth levels, 3 for 8 membership functions and 7 for 128 levels of resolution. The memory depth is given by the dimension of the universe of the discourse (128 in our case) and it will be represented by the memory rows. The length of a world of memory is defined by: Length = nem (dm(m)+dm(fm) Where: fm is the maximum number of non null values in every element of the universe of the discourse, dm(m) is the dimension of the values of the membership function m, dm(fm) is the dimension of the word to represent the index of the highest membership function. In our case then Length=24. The memory dimension is therefore 128*24 bits. If we had chosen to memorize all values of the membership functions we would have needed to memorize on each memory row the membership value of each element. Fuzzy sets word dimension is 8*5 bits. Therefore, the dimension of the memory would have been 128*40 bits. Coherently with our hypothesis, in fig. 1 each element of universe of the discourse has a non null membership value on at most three fuzzy sets. Focusing on the elements 32,64,96 of the universe of discourse, they will be memorized as follows: The computation of the rule weights is done by comparing those bits that represent the index of the membership function, with the word of the program memor . The output bus of the Program Memory (μCOD), is given as input a comparator (Combinatory Net). If the index is equal to the bus value then one of the non null weight derives from the rule and it is produced as output, otherwise the output is zero (fig. 2). It is clear, that the memory dimension of the antecedent is in this way reduced since only non null values are memorized. Moreover, the time performance of the system is equivalent to the performance of a system using vectorial memorization of all weights. The dimensioning of the word is influenced by some parameters of the input variable. The most important parameter is the maximum number membership functions (nfm) having a non null value in each element of the universe of discourse. From our study in the field of fuzzy system, we see that typically nfm 3 and there are at most 16 membership function. At any rate, such a value can be increased up to the physical dimensional limit of the antecedent memory. A less important role n the optimization process of the word dimension is played by the number of membership functions defined for each linguistic term. The table below shows the request word dimension as a function of such parameters and compares our proposed method with the method of vectorial memorization[10]. Summing up, the characteristics of our method are: Users are not restricted to membership functions with specific shapes. The number of the fuzzy sets and the resolution of the vertical axis have a very small influence in increasing memory space. Weight computations are done by combinatorial network and therefore the time performance of the system is equivalent to the one of the vectorial method. The number of non null membership values on any element of the universe of discourse is limited. Such a constraint is usually non very restrictive since many controllers obtain a good precision with only three non null weights. The method here briefly described has been adopted by our group in the design of an optimized version of the coprocessor described in [10].

  • PDF