• 지구물리와물리탐사
• /
• 제5권3호
• /
• pp.169-177
• /
• 2002

일반 반사법 탐사 모델링에서 효율적인 주파수영역 수치모델링의 실용화를 위해 무엇보다 해결해야할 과제는 파장당 격자수를 줄이는 것이다. 본 연구에서는 이에 착안하여 수치분산 및 수치이방성을 최소화시키면서 한 파장당 필요한 격자수를 줄일 수 있는 가중평균 유한요소법을 개발하였다. 강성행렬과 질량행렬은 네 개의 사각형 유한요소로 확장하였으며, 모든 격자점이 포함되도록 유한요소를 배열하여 조합하였다. 확장된 네 개의 강성행렬과 질량행렬은 가중평균계수를 주어 선형결합하는 방법으로 가중평균하였다. 가중평균계수는 확장된 25점 평균차분법을 사용하여 가중평균계수를 결정하였다. 또한, 정확도 향상은 2차원 균질모델 과 수평층 모델에서 해석해와 한 파장당 4개의 격자점을 준 가중평균 유한요소법 수치해 비교를 통하여 검증하였다. 또한, 석유탐사에서 활용성이 높은 향사구조 모델을 선정하여 이의 반응을 관찰한 결과 지층경계면외에 네 개의 사각형 유한요소들의 구성으로 인한 인위적인 파의 도달이 인식되지 않았다. 따라서, 본 연구에서 고안된 가중평균 유한요소법은 주파수영역에서 폭 넓은 수치모델링연구을 가능하게 할 것이다.

• 화약ㆍ발파
• /
• 제37권1호
• /
• pp.1-13
• /
• 2019

발파굴착을 수행하는 터널시공 현장에서, 발파진동이 반복적으로 작용할 경우 인근 보안물건에 손상을 입힐 수 있으며, 이러한 발파진동을 저감시킬 수 있는 방법으로는 전자뇌관의 사용, 분산장약, 심발공법의 변경 등이 고려될 수 있다. 본 연구에서는 수치해석 기법을 이용하여 심발공법에 따른 발파진동 저감효과를 분석하고자, 3차원 유한차분법 프로그램인 FLAC3D를 이용하여 발파공, 지연시차, 장약량 등을 수치모델에 반영하고, 발파공별로 동하중을 적용하여 수치해석을 수행하였다. 심발공법은 터널 시공현장에서 많이 사용되고 있는 V-cut 공법 및 다단평행의 심발공법을 적용하였으며, 실제 터널 시공현장에서 시험발파를 수행하여 현장의 발파진동 계측자료와 수치해석의 진동속도를 비교하였다. 수치해석과 시험발파의 발파진동은 매우 유사한 경향을 나타낸 것으로 분석되었다.

• 한국진공학회지
• /
• 제15권5호
• /
• pp.506-511
• /
• 2006

$Zn(NO_3)_2{\cdot}6H_2O$, 수산화나트륨, cyclohexylamine, 에탄올, 물이 혼합된 용액을 수열합성하여 ZnO 나노로드를 합성하고, 합성한 물질의 이산화 질소$(NO_2)$와 일산화 탄소(CO)에 대한 감응특성을 조사하였다. 혼합 용액에 첨가되는 물의 양을 변화시켜 ZnO 나노로드의 형상과 응집현상을 조절할 수 있었다. 이는 물과 cyclohexylamine의 반응에 의해 발생되는 $OH^-$ 이온의 농도변화에 의한 것으로 해석된다. 물의 함량이 낮을 때에는 뭉쳐진 성게모양의 ZnO 나노로드를, 물의 함량이 많을 때에는 잘 분산된 ZnO 나노로드를 각각 합성할 수 있었다. 잘 분산된ZnO 나노로드는 공기 중에서 50 ppm 의 CO에 노출되었을 때 주목할 만한 반응을 보이지 않는 반면, 1 ppm 의 $NO_2$에 노출되었을 때에는 저항이 1.8배 증가하였다. 이러한 선택적 반응을 보이는 ZnO 나노로드는 자동차용 자동환기 시스템의 핵심부품인 매연센서의 감응물질로 사용될 수 있다.

• 한국환경과학회:학술대회논문집
• /
• 한국환경과학회 2008년도 추계학술발표회 발표논문집
• /
• pp.428-432
• /
• 2008

대수층의 저유량이 풍부한 강변여과수 개발 예정지역의 충적충(지표면하 25$\sim$35 m 구간)에서 수리전도도와 증분산지수의 규모종속효과를 규명하기 위해 양수시험과 수렴흐름 추적자시험이 수행되었다. 양수시험과 추적자시험의 규모는 2 m 와 5 m 이었으며 양수시험은 5개 공, 추적자시험은 3개 공을 이용하여 수행되었다. 양수시험은 일정한 양수율(2,500 m$^3$/day)로 수행되었으며, 양수 시작 후 경과시간에 따른 수위변화 자료를 AQTESOLV 3.5 프로그램에 입력하여 해석하였다. 시험대수층의 수리전도도는 양수정에서 1.745$\times$10$^{-3}$ m/sec, 양수정에서 이격거리가 2 m 구간에서는 2.161$\times$10$^{-3}$ m/sec와 2.270$\times$10$^{-3}$ m/sec 이었으며, 이격거리가 5 m 구간에서는 2.452$\times$10$^{-3}$ m/sec와 2.591$\times$10$^{-3}$ m/sec로 산정되었다. 그리고, 양수정에서 회복시험 시 Theis(Recovery) 방법에 의해 해석된 수리전도도는 1.603$\times$10$^{-3}$ m/sec이었다. 양수정에서 관측정의 이격거리(d)에 따른 수리전도도(K) 증가함수는 log K = 0.0693logd-2.071와 log K = 0.08171og d-2.655로 추정되었으며, 결정계수는 각각 0.965와 0.979로서 매우 높게 나타났다. 따라서 양수정에서의 이격거리가 멀수록 수리전도도가 증가하는 규모종속을 확인하였으며, 또한 시험대수층의 수리전도도가 방사상으로 유사하게 분포하고 있음을 알 수 있었다. 수렴흐름 추적자시험의 양수율은 2,500 m$^3$/day 이었으며, 2개의 주입정에 염소이온 5 kg을 순간 주입하였다. 염소이온의 농도이력곡선을 작성하여 초기도달시간과 최고농도의 차이를 분석하였으며, 누적질량회수곡선을 통해 양수 후 경과시간에 따른 염소이온의 질량회수율을 분석하였다. 그리고, 염소이온농도 대 누적질량회수율의 이력그래프를 작성하여 누적질량회수율에 따른 염소이온농도의 증가와 감소 변화를 분석하였다. 또한, 염소이온농도의 증가/감소 구간에 대한 선형회귀분석을 수행하여 농도 증가율과 감소율의 변화를 파악하였다. 양수정에서 관측된 경과시간별 염소이온농도 자료를 CATTI 코드의 "Converging Radial Flow With Instantaneous Injection" 해석법에 적용하여 종분산지수를 추정하였다. 양수정에서 이격거리가 2 m인 경우의 종분산지수는 0.4152 m, 이격거리가 5 m인 경우의 종분산지수는 3.2665 m 이었다. 따라서 양수정에서 이격거리가 멀수록 종분산지수가 증가하는 규모종속효과를 확인하였으며, 또한 이격거리에 대한 종분산지수의 비는 각각 0.21과 0.65 정도로서 증가하였다.

• 한국전산구조공학회논문집
• /
• 제27권3호
• /
• pp.199-206
• /
• 2014

본 논문에서는 변단면 기능경사재료 보에서 중립면 탄성계수가 축방향을 따라 공간적 불확실성을 가질 경우에 대한 구조 응답변화도 산정을 위한 정식화에 대해 논한다. 기능경사재료는 두 이질재료의 체적비가 두께방향으로 연속적으로 변하며 고체화되는 과정으로 제작되는 재료로서 온도 및 응력 등에서 연속적인 변화를 가능하게 하여, 전통 복합재료에서 나타나는 층분리나 균열 발생 등이 제거되는 장점을 가지고 있다. 그러나 이론적으로 설정된 기능경사에 맞는 재료의 제작이 어려우며, 이에 따라 내재적인 불확실성을 가지고 있다. 이를 모사하기 위하여 중립면 탄성계수에서의 불확실성을 추계장으로 모델링하고, 추계적분에 의한 확률변수를 도입하여, 변위의 1, 2차모멘트를 산정할 수 있는 방법을 제시하였다. 제안된 해석 방법은 스펙트럼모사법을 적용한 몬테카를로 해석으로 검증하였다. 추계장의 상관관계거리에 따른 분산계수의 변화, 재료지수 및 기하인수가 변위의 분산계수에 미치는 영향 등을 고찰하였고, 몬테카를로 해석 대비 제안 해석법의 효율성에 대해서도 논하였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2005년도 춘계학술대회논문집
• /
• pp.372-378
• /
• 2005

Finite element analysis of concrete slab system in apartment building was executed using the tables of orthogonal arrays, and optimal design process was proposed. At first, experimental results show that sound peak components to influence the overall level and the rating of floor impact sound insulation were coincident with natural frequencies of the reinforced concrete slab. Finite element model of concrete slab was compared with experimental results, and well corresponded with an error of less than 10%. The tables of orthogonal arrays were used for finite element analysis with 8 factors. 3 related to material properties and 5 related to slab shape parameters and its results were analyzed by statistical method, ANOVA. The most effective factor among them was slab thickness, and main effect factor from slab shape parameters was different from each natural frequency. The interaction was found in the higher mode over $3^{rd}$ natural frequency. From main effect plot and interaction plot, the optimal design factor to increase the natural frequency was determined.

• 대한기계학회논문집A
• /
• 제33권9호
• /
• pp.963-968
• /
• 2009

This paper deals with shrink fit analysis of rotor by 2D cross-section, 2D axis-symmetry, and 3D FEM model. And this paper presents 2nd order approximation function of thermal expansion displacement by design variables (shape dimension, heating temperature, sleeve length, interference etc.), table of orthogonal array and RSM(response surface methodology). The possibility of the rotor with shrink fit is evaluated by thermal expansion displacement. If thermal expansion displacement is larger than interference, shrink fit enable to make the rotor. 2D axis-symmetry model and 3D model are more reasonable than 2D cross-section model, because stress and strain is different along length of shaft.

• 한국결정성장학회지
• /
• 제9권2호
• /
• pp.240-244
• /
• 1999

화학기상증착법(CVD)을 사용하여 graphite 기판 위에 단층의 SiC와 SiC/C 경사기능층을 증착시켰다. 먼저 상업적인 전산 구조 해석 프로그램을 사용하여 열충격하에 있는 재료층 내부의 온도 분포 및 열응력 분포를 계산하였고, 계산에 의해 설계, 제작된 경사기능층 시편들의 열특성을 조사하였다. 열충격 시험결과, 경사기능층 시편들이 계면에서 매우 효율적인 열응력 이완특성을 나타내는 것을 알 수 있었고, $\Delta$T=1600K의 열충격에도 변형이 없음을 알 수 있었다.

• 콘크리트학회지
• /
• 제9권5호
• /
• pp.189-196
• /
• 1997

본 연구에서는 유한요소법을 이용한 콘크리트 인장균열의 발생 및 전파양상을 모형화하기 위하여 분산균열(smeared crack) 모델의 일종인 2차원 균질화된 균열(homogenized crack)모형을 제안하였다. 제안한 모형은 인장균열면을 따라 속도 불연속계(velocity discontinuity)를 도입하고 평형방정식 빛 적합방정식을 이용하여 인장균열을 포함한 콘크리트의 구성방정식을 유도할 수 있으며 인장균열이 소성연화거동을 하는 경우. 유한요소망내에서 객관성이 유지될 수 있음을 밝히기 위하여 1차원 영역내에서 엄밀해를 유도하였다. 제안한 모형을 이용한 1차원 또는 2차원 유한요소 해석은 기존의 노치를 포함한 콘크리트 시편에 대한 실험과 상응하는 결과를 보였을 뿐만 아니라 유한요소망의 객관성이 유지되고 있음을 밝혔다. 제안한 모형은 큰 어려움없이 3차원 영역으로 확대할 수 있을 것이며 이에 대한 추가적인 연구가 계속될 것이다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1999년도 봄 학술발표회 논문집
• /
• pp.50-57
• /
• 1999

Structural analysis of large-scale structures involving large amount of computational load and data storage requires high-performance computing resources. We have previously developed PC-level distributed structural analysis algorithms based on substructuring technique where each personal computer assigned to a slave node has been involved in the computations for single substructures. Recently, it has been proved by the authors that the performance of distributed structural analysis algorithm can be further enhanced by changing substructuring schemes. Therefore a new distributed structural analysis algorithm with one PC to multiple substructures scheme is presented in this paper. The algorithm is implemented on the network of multiple personal computers and applied to structural analysis of two dimensional frame structures.