• 제목/요약/키워드: 수치해법

검색결과 374건 처리시간 0.025초

MT 법의 3차원 모델링 개관 (Review on the Three-Dimensional Magnetotelluric Modeling)

  • 김희준;남명진;송윤호;서정희
    • 지구물리와물리탐사
    • /
    • 제7권2호
    • /
    • pp.148-154
    • /
    • 2004
  • 자기지전류(MT)법의 3차원 모델링에 대해 소개한다. 3차원 MT 모델링은 MT 반응의 물리적 특성의 이해뿐만 아니라 지하의 3차원적 전기비저항 구조를 재구성하기 위한 역산법의 개발에도 필수적이다. 지난 20년 동안 3차원 모델링에 관한 여러 수치기법들이 개발되었으나 그 실용성에는 많은 한계가 있었다. 그러나 최근에는 컴퓨터의 급속한 발전과 대형 연립방정식에 대한 반복해법의 발전에 힘입어 이전에는 어려웠던 복잡한 3차원 구조에 대한 MT 반응을 효율적으로 모델링할 수 있게 되었다. 유한차분법에서는 자기 flux와 전류의 보존법칙을 만족하면서 전기장의 불연속을 표현할 수 있는 staggered 격자의 사용이 보편화되었다. 대형 연립방정식에 대한 수치해의 수렴성은 Krylov 부분공간법, 적당한 전처리 기술 및 정적 발산보정법을 채택함으로써 크게 향상된다. 변요소를 사용하는 벡터 유한요소법으로도 전기장의 불연속 문제를 해결할 수 있으며 이 방법이 가진 기하학적 유연성은 불규칙한 지표기복을 포함한 복잡한 구조를 모델화할 때 특히 유용하다.

외기대류계수에 관한 실험적 연구 (Experimental Study on Coefficient of Air Convection)

  • 전상은;김진근
    • 콘크리트학회논문집
    • /
    • 제15권2호
    • /
    • pp.305-313
    • /
    • 2003
  • 본 연구에서는 외기와의 열전달을 나타내는 외기대류계수에 관한 실험을 실시하였다. 외기대류계수에 관한 기존의 모델에서 나타났던 문제점을 해결하기 위해 실험변수로 풍속외에 거푸집의 종류(목재, 철재)를 선정하였다. 실험결과를 이용하여 외기대류계수를 산정하고자 열평형 방정식을 이용한 수치해법을 사용하였으며, 이론적인 고찰을 통해 각 거푸집별로 풍속에 따른 외기대류계수의 변화를 예측할 수 있는 모델식을 제안하였다. 제안된 모델식에 의하면, 모든 경우에 풍속에 따라 외기대류계수가 증가하는 경향을 보였으나 거푸집의 사용여부나 거푸집 재료에 따라 다른 양상을 보이는 것을 알 수 있었다. 이러한 양상의 차이는 거푸집 재료의 열특성에 의해 결정되는 것으로 외기대류계수는 거푸집이 없는 경우, 철재 거푸집을 사용한 경우, 목재 거푸집을 사용한 경우의 순으로 풍속의 영향을 받는 것으로 나타났다. 제안된 모델식을 이용하면 수화열에 의한 콘크리트 구조물의 온도해석시 보다 정확한 결과를 얻을 수 있을 것으로 사료되며, 향후 이러한 열특성계수에 대한 연구가 필요할 것으로 판단된다.

도시철도 역사 스크린 도어 개폐에 따른 냉방 기류 해석 및 효율 비교 분석 (Analysis of Cooling Air Current and Efficiency of Air Conditioning in the Underground Subway Station with Screen-Door Opening and Closing)

  • 장용준;류지민;정호성
    • 한국철도학회논문집
    • /
    • 제17권5호
    • /
    • pp.328-335
    • /
    • 2014
  • 도시철도 지하역사 냉방 기류 및 냉방 효율을 조사하기 위하여 수치해법을 이용하여 해석하고 현장 실험 결과와 비교하여 분석하였다. 해석 대상 역사로는 지하 8층의 깊이 43.6m인 서울 5호선 신금호 역사를 선정하였다. 전체 역사를 해석 영역으로 하였으며, 공조기 모드는 평상시 모드로 고정시켰다. 냉방 공조를 위하여 대합실 천정에 총 94개의 정사각형($0.6m{\times}0.6m$) 환기구를 모델하였으며, 승강장은 총 222개의 환기구가 승강장 천정에 모델되었다. 대합실에서 급기되는 공기는 $47,316m^3/h$, 배기되는 공기량은 $33,980m^3/h$이며, 승강장에서 급기되는 공기는 $33,968m^3/h$, 배기되는 공기량은 $76,190m^3/h$로 현장의 풍량을 반영하였다. 승강장에서 스크린도어(PSD)는 닫힌 경우와 열린 경우 각각을 조사하였다. 총 750만개의 격자가 사용되었으며, 전체 영역을 22개의 다중 블록으로 나누어서 계산하고, MPI를 이용하여 각각의 블록에서 계산된 결과를 교환하였다. LES 기법을 이용하여 운동량 방정식 및 에너지 방정식을 계산하였다.

유전자 알고리즘을 이용한 동적통행배정에 관한 연구 (Dynamic Traffic Assignment Using Genetic Algorithm)

  • 박경철;박창호;전경수;이성모
    • 대한공간정보학회지
    • /
    • 제8권1호
    • /
    • pp.51-63
    • /
    • 2000
  • 최근 교통문제를 해결하기 위한 방법으로 교통계획분야에 GIS나 ITS를 활용한 다양한 연구가 활발히 진행 중에 있다. 이와 함께 정보환경의 급격한 발달과 더불어 대안 경로의 선정, 또는 교통예보 서비스와 같은 온라인 상에서의 교통정보 제공이 이루어지고 있어 GIS 환경 내에서도 가로망의 교통량을 정확하게 예측할 수 있는 기능이 요구되고 있어 통행배정모형의 중요성이 증가하고 있다. 그런데, 전통적인 정적 통행배정모형은 급변하는 교통상황에 적합하지 않기 때문에 실시간 교통상황에 대한 교통흐름을 예측할 수 있는 동적 통행배정모형의 개발이 요구되고 있다. 그러나, 동적 통행배정모형은 시공간적인 변수들의 복잡성으로 인해서 그 최적해를 찾는데 많은 수학적인 어려움과 제약조건이 존재한다. 따라서, 이를 해결하기 위한 여러 가지 해법이 연구되어왔지만, 기존의 방법은 목적함수나 제약조건이 convex(하지 않은 경우에는 적용이 불가능한 단점을 가지고 있다. 본 연구에서는 인공지능방법(Artificial Intelligence Technique)의 한 분야로 활발히 연구되고 있는 유전자 알고리즘(Genetic Algorithm)을 동적 통행배정 모형에 도입하여 그 해결 방법을 제시하였다. 논문에서 사용한 동적 통행배정모형은 제약조건이 convex(하지 않은 Merchant-Nemhauser모형이고, 새로운 해결기법으로 사용된 유전자 알고리즘은 일반적인 제약조건을 처리할 수 있다고 알려진 GENOCOP III시스템이다. 새로 도입된 방법의 효율성과 유의성을 검증하기 위해 간단한 네트워크에 적용하였다. 그 결과 GENOCOP III 시스템이 계산과정의 효율성에 있어서 기존의 비선형 해법 알고리즘보다 우수한 것으로 입증되었다.연구가 진행되어야 할 것이다. 실질적으로 성감별 수정란의 대량생산이 가능할 것으로 사료되며, 농가차원에서 산업적 실용화가 될 수 있을 것으로 기대한다.twork descrition)를 통해 교통분석후의 제반 교통특성(교통량, 교통량/용량 비(比), 속도 등)을 교통망상에 표시할 수 있음으로서 의사결정에 보다 많은 도움을 줄 수 있을 것이다. 비트율의 증가와 화질 열화는 각각 최대 1.32%와 최대 0.11dB로 무시할 수 있을 정도로 작음을 확인 하였다.을 알 수 있었다. 현지관측에 비해 막대한 비용과 시간을 절약할 수 있는 위성영상해석방법을 이용한 방법은 해양수질파악이 가능할 것으로 판단되며, GIS를 이용하여 다양하고 복잡한 자료를 데이터베이스화함으로써 가시화하고, 이를 기초로 공간분석을 실시함으로써 환경요소별 공간분포에 대한 파악을 통해 수치모형실험을 이용한 각종 환경영향의 평가 및 예측을 위한 기초자료로 이용이 가능할 것으로 사료된다.염총량관리 기본계획 시 구축된 모형 매개변수를 바탕으로 분석을 수행하였다. 일차오차분석을 이용하여 수리매개변수와 수질매개변수의 수질항목별 상대적 기여도를 파악해 본 결과, 수리매개변수는 DO, BOD, 유기질소, 유기인 모든 항목에 일정 정도의 상대적 기여도를 가지고 있는 것을 알 수 있었다. 이로부터 수질 모형의 적용 시 수리 매개변수 또한 수질 매개변수의 추정 시와 같이 보다 세심한 주의를 기울여 추정할 필요가 있을 것으로 판단된다.변화와 기흉 발생과의 인과관계를 확인하고 좀 더 구체화하기 위한 연구가 필요할 것이다.게 이루어질 수 있을 것으로 기대된다.는 초과수익률이 상승하지만, 이후로는 감소하므로, 반전거래전략을 활용하는 경우 주식투자기간은 24개월이하의 중단기가 적합함을

  • PDF

사질토를 지나 풍화암에 소켓된 매입 PHC말뚝에서 지반의 허용압축지지력 산정도표 및 산정공식 개발에 관한 연구(VII) - 지반의 허용압축지지력 산정공식 - (Study(VII) on Development of Charts and Equations Predicting Bearing Capacity for Prebored PHC Piles Socketed into Weathered Rock through Sandy Soil Layers - Allowable Axial Compressive Bearing Capacity Formulae -)

  • 권오균;남문석;이원제;여규권;최용규
    • 한국지반공학회논문집
    • /
    • 제35권12호
    • /
    • pp.69-89
    • /
    • 2019
  • 사질토를 지나 풍화암에 4D 소켓된 매입 PHC말뚝에 대한 매개변수 수치해석 자료를 분석하여 지지력 표해 및 도해(I)이 제안된 바 있다(Choi et al., 2019a). 본 논문에서는 이를 활용하여 직경의 5% 침하량에서 발현되는 동원지지력을 산정할 수 있는 새로운 산정공식을 제안하였다. 제안식은 두 가지이며, 그 중에서 EQ-G1 제안식은 각각의 상대근입길이(L/D)에서 말뚝직경과 ${\bar{N}}$값(보정 N 값)을 고려하여 각 지지력 성분을 평가하며, 검증 결과는 다음과 같다. RQP는 약 71~94%로 나타났으며, 이는 일반 설계에서 사용하고 있는 지지력 산정공식으로 구한 SRF(26~37%)보다 큰 값을 나타냈다. RQP는 적정 설계의 범위로 분포하였는데, 4개 사례에서는 78, 136, 142, 180%로 나타났다. DE는 β1에 따라 달라질 수는 있는데, 본 연구에서는 0.85로 가정하였으므로 DE는 약 85%로 나타났다. 따라서 본 연구에서 제안한 EQ-G1 제안식을 사용하여 PHC말뚝 몸체의 허용압축하중(Pall)까지 활용할 수 있는 적정설계를 수행할 수 있는 것으로 판단된다. 그리고 EQ-G2-3 제안식은 D, ${\bar{N}}$ 및 L/D를 동시에 고려하여 지지력 성분들을 평가할 수 있으며, 매입 PHC말뚝에서 지반의 압축지지력 산정 시 근사해법으로 활용할 수 있다.

절편법과 유한차분법에 의한 사면안정해석 비교연구 (Slope Stability Analysis by Slice Method and Finite Difference Method- A Comparative Study -)

  • 박연준;채영수;유광호;백영식
    • 한국지반공학회논문집
    • /
    • 제15권6호
    • /
    • pp.263-272
    • /
    • 1999
  • 사면안정 문제에 널리 이용되는 절편법은 사용이 간편하고 계산시간이 짧다는 장점이 있으나, 각 절편들 사이에 작용하는 힘에 대한 가정에 따라 여러 가지 방법으로 분류되며, 그 결과 또한 다르게 나오는 단점이 있다. 이러한 방법간의 차이는 경우에 따라 매우 심하여 공학적인 결정을 어렵게 하는 경우가 드물지 않다. 본 논문에서는 이러한 문제에 대한 해결방안의 하나로 유한요소법 및 유한차분법을 이용한 사면의 안정성 해석 및 안전율 계산방법을 제시하고자 한다. 이를 위하여 유한차분법을 이용한 지반해석용 프로그램인 FLAC을 이용하여 사면의 안전율을 계산하는 routine을 작성하고, Chen의 한계평형해와 비교하여 검증하였다. Fredlund 와 Krhan이 제시한 문제에 적용하여 그 결과를 흔히 사용되는 몇 가지 절편법에 대한 해석결과와 비교한 결과 균질한 지반의 경우에는 절편법과 0.01 이내의 차이를 보였으며, 연약면이 존재하는 경우에는 절편법에 의한 해와 5% 이내의 차이를 보였다. 대부분의 지반해석용 프로그램들은 복잡한 지반조건은 물론, 지반보강재의 효과도 고려할 수 있으므로 수치해석에 의한 해법이 공학적인 결정에 많은 도움을 줄 수 있을 것으로 생각되며, 개발된 routine은 절편법이 제한적으로만 적용되는 절리암반사면이나 풍화암 사면의 안전율 계산에도 활용될 수 있을 것으로 기대된다.

  • PDF

열(熱)흐름을 동반(同伴)한 정상지하수(定常地下水)의 흐름해석(解析) 수치모형(數値模型) (A Numerical Model for Analysis of Groundwater Flow with Heat Flow in Steady-State)

  • 왕수균;조원철;이원환
    • 대한토목학회논문집
    • /
    • 제11권4호
    • /
    • pp.103-112
    • /
    • 1991
  • 본(本) 연구(硏究)는 일정균등(一定均等)한 열적(熱的) 특성(特性)을 가지고 있으며 상변화(相變化) 없는 등방(等方) 이질성(異質性)의 3차원(次元) 대수층계(帶水層系)의 열(熱)흐름과 정상상태(定常狀態)의 지하수(地下水)흐름을 모의발생(模擬發生)할 수 있는 유한차분(有限差分) 모형(模型)을 확립(確立)한 것이다. 이 모형(模型)은 대규모(大規模) 지하수(地下水) 흐름체계(體系)에서 폐기물(廢棄物)의 지하저류시(地下貯溜時) 지하수(地下水) 흐름과 발생(發生) 혹은 주입(注入)된 열(熱)의 흐름을 예측(豫測) 분석(分析)하기 위하여 확립(確立)된 것이다. 이러한 대수층계(帶水層系)의 지하수(地下水) 흐름에 작용(作用)하는 조건(條件)으로는 강우주입(降雨注入)으로 인한 수문학적(水文學的) 조건(條件)과 고정(固定)된 수리수두(水理水頭) 경계조건(境界條件) 등(等)이 포함(包含)되고, 열(熱)흐름에는 지열(地熱)의 흐름, 지표면(地表面)으로의 전도(傳導), 주입(注入)에 의한 이류(移流), 고정(固定) 수두경계(水頭境界)로 향(向)한 또는 고정수두경계(固定水頭境界)로 부터의 이류(移流) 등(等)이 포함(包含)된다. 본(本) 모형(模型)에서는 지하수(地下水)흐름과 열(熱)흐름 방정식(方程式)을 번갈아 푸는 교대반복과정(交代反復過程)을 사용(使用)하고, 두 방정식(方程式)의 계산(計算)에는 직접해법(直接解法)을 사용(使用)한다. 이동시간(移動時間)은 모형공간(模型空間)에서 입자추적(粒子追跡)으로 결정(決定)되며, 분할(分轄)된 구역내(區域內)의 지하수(地下水) 유속(流速)은 구역내(區域內)의 유속(流速)을 선형(線形)으로 보간(補間)하여 계산(計算)한다. 본(本) 모형(模型)을 경상북도(慶尙北道) 영일군(迎日郡) 송라면(松羅面) 지경리(地境里) 일대(一帶)의 지하수계(地下水系)에 적용(適用)하여 이 일대(一帶) 지하(地下) 암반층(岩盤層)의 수두분포(水頭分布), 유동로(流動路), 이동시간(移動時間) 및 지하수온분포(地下水溫分布)를 계산(計算)하여 지하수(地下水) 유동체계(流動體系)를 분석(分析)하였다.

  • PDF

건착망의 운동특성(II) -그물감이 다른 건착망의 침강특성 해석- (Characteristics on the Motion of Purse Seine(II) -An Analysis on the Sinking Characteristics of Model Purse Seine by Different Netting Material-)

  • 김석종;박정식
    • 수산해양기술연구
    • /
    • 제31권4호
    • /
    • pp.372-378
    • /
    • 1995
  • 건착망의 침강 운동특성을 구명하기 위한 기초연구로서 그물감의 재료가 다른 세 종류의 건착망을 수중 총중량(60g)이 동일하게 되도록 발돌의 양을 조절하여 제작하고, 침강특성을 해석하였다. 실험에 사용한 건착망은 그물실의 직경 및 발의 길이가 같은 폴리프로피렌계(밀도 0.91g/$cm^3$) 및 폴리에스터계(밀도 1.38g/$cm^3$)의 매듭 없는 그물감을 사용하여, 뜸줄의 길이를 420cm, 그물의 폭을 86cm가 되도록 제작하고, 이 그물들을 각각 PP, PA 및 PES 그물이라고 하였다. 회유수조의 수로상에 투망장치를 설치해서 정지상태의 수중에 투망하고, 측면에 설치한 비디오 카메라를 이용하여 촬영 녹화하였다. 그리고, 그물에 표시한 측정점의 좌표를 화상해석장치로 읽고 실험치를 구하였다. 여기에는, 건착망의 수직방향의 침강운동을 나타내는 미분방정식을 구하고 Runge-Kutta-Gill법에 의한 연립 미분방정식 해법을 이용하여 수치해석을 행하였으며, 그 결과는 다음과 같다. 1. 그물 아랫자락의 침강속도는 PP 그물이 가장 빠르고, PA 및 PES 그물 순으로 늦게 나타났다. 2. 그물감의 저항계수 $K_D$는 계산결과 $K_D=0.061({\frac{\rho}{{\rho}_w}})^4$의 관계식으로 나타낼 수 있었다. 3. 그물다발의 저항계수 $C_R$은 계산결과 $C_R=0.91({\frac{\rho}{{\rho}_w}}$의 관계식으로 나타낼 수 있었다. 4. 건착망 투망후 경과시간에 따른 그물 아랫자락의 도달수심에 대한 실험치와 계산치는 매우 잘 일치하여 meas.=0.99cal.의 관계였다.

  • PDF

FEM을 이용한 왕복동 공기압축기의 피스톤 및 커넥팅로드의 구조해석 (Transient Structural Analysis of Piston and Connecting Rods of Reciprocating Air Compressor Using FEM)

  • ;양창조;김준호;김부기
    • 해양환경안전학회지
    • /
    • 제23권4호
    • /
    • pp.393-399
    • /
    • 2017
  • 왕복동식 압축기에서 피스톤과 커넥팅로드는 중요한 부분이다. 이러한 주요부에 기계적 부하가 과도하게 가해지면 해당 기부속이 손상될 수 있으며, 교체하기도 쉽지 않고 비용도 많이 든다. 따라서 내구성과 수명에 영향을 미치는 요인을 분석할 필요가 있다. 본 연구의 주요 목적은 피스톤과 커넥팅로드의 최대 응력 집중 위치를 확인하는 것이다. 이를 위해 설계된 공기압축기의 작업 공정의 동적계산을 기반으로 피스톤 및 커넥팅로드의 응력 분석을 수행하였다. 공기압축기의 피스톤과 커넥팅로드의 3 차원 모델을 따로 설계하고, 이러한 부품들의 유한요소 해석은 수치해석적인 근사해법을 사용하였다. 피스톤은 열 경계 조건 없이 크랭크 샤프트의 각도에 따라 압력 부하를 받는다. 시뮬레이션 결과는 피스톤과 커넥팅로드의 응력 집중 위치와 그 값을 예측하고 추정할 수 있다. 그 결과 크랭크 각도 $135^{\circ}$$225^{\circ}$에서 피스톤은 190MPa, 커넥팅로드는 123MPa 이상의 최대 등가응력이 나타났으며 이는 인장 항복강도 이하의 값이다. 또한, 커넥팅로드와 피스톤에 계산 된 안전 계수는 1보다 높게 나타났다. 더욱이, 이러한 결과는 왕복동 공기압축기 제작사에 피스톤 및 커넥팅로드를 설계함에 있어서 최적화를 위한 참고 자료로 활용 될 수 있다.

선망의 침강 저항 해석 - d/l이 다른 모형망의 경우 - (An Analysis of Sinking Resistance for Purse Seine - In the Case of the Model Seine with Different d/l -)

  • 김석종;박정식
    • 수산해양기술연구
    • /
    • 제34권3호
    • /
    • pp.274-282
    • /
    • 1998
  • 그물실의 직경과 발 길이의 비율 d/l이 선망의 침강 운동에 어떤 영향을 끼칠 것인가를 해명하기 위해서 그물실의 직경은 0.45mm(PES 28 tex${\times}$2${\times}$2)의 것으로 모두 같고, 발 길이는 4.3, 5.0, 5.5, 6.0, 6.6 및 7.7mm로 각각 다르게 편망된 6종류의 무결절 그물감을 사용해서 뜸줄의 길이 450cm, 그물의 폭 85cm가 되도록 제작하고 뜸 160g과 발돌 50g을 매달아서, 각각 I,II,III,IV,V및 Ⅵ형 그물이라고 했다. 회류 수조의 수로 상에 투망 장치를 설치해서 정지 상태의 수중에 선망 그물을 투망하고 상방과 측면에 설치한 비디오 카메라를 사용하여 그물의 침강과정을 촬영 톡화했으며, 그물에 표시한 측정점의 좌표를 화상 해석 장치로 읽고 실험 값을 구했다. 그리고, Runge-Kutta-Gill법에 의한 연립 미분방정식 해법을 이용하여 수치 해석을 행하였으며, 그 결과는 다음과 같다. 1.그물 아랫자락의 침강속도는 d/l이 가장 작은 Ⅵ형 그물이 가장 빠르고, V,IV,III,II및 I형 그물 순으로 늦게 나타났다. 2.그물 벽에 대한 저항 계수는 d/l의 크기에 의존하여 $K{\arg}$=0.081$(d/l)^-0.5$의 관계식으로 나타낼 수 있었다. 3.그물 다발의 저항 계수는 d/l의 크기에 의존하지 않은 $C_R=0.91(\frac{p}{$p_u$})$의 관계식으로 나타낼 수 있었다.

  • PDF