• 제목/요약/키워드: finite field method

검색결과 2,149건 처리시간 0.044초

궤도 진동기반의 침목플로팅궤도 침목방진패드 스프링강성 추정 기법 연구 (Estimation Method of Resilience Pads Spring Stiffness for Sleeper Floating Tracks based on Track Vibration)

  • 최정열;박상욱;정지승
    • 문화기술의 융합
    • /
    • 제9권6호
    • /
    • pp.1057-1063
    • /
    • 2023
  • 본 연구대상인 도시철도 침목플로팅궤도(STEDEF)는 구조물로 전달되는 진동을 저감시키기 위한 방진궤도이다. 현재 침목플로팅궤도의 침목방진패드 교체주기(정적 스프링강성 변화율, 25±2%)는 하중기반(궤도충격계수와 궤도지지강성)으로 설정되어 운영중인 실정이다. 그러나 대부분의 선행연구는 침목방진패드의 피로수명평가와 스프링강성 증가에 따른 궤도충격계수 및 궤도지지강성의 증가 등 하중기반의 구조적 안전성 측면의 연구가 진행되었다. 따라서 본 연구에서는 분석 구간별 도상 진동가속도를 측정하고 700만회 피로시험결과를 이용하여 구간별 침목방진패드 스프링강성을 산출하고자 한다. 구간별 산출한 침목방진패드 스프링강성을 해석제원으로 설정하여 도상 진동가속도를 해석적으로 도출하였다. 구간별 해석 도상 진동가속도가 현장측정 도상 진동가속도 범위 이내로 나타나 해석모델링의 적정성이 검증되었다. 도출된 스프링강성 변화에 따른 진동가속도 선도(g-k curve)를 이용하여 측정 도상 진동가속도로 침목방진패드 스프링강성을 추정하고자 한다. 따라서 측정 도상 진동가속도를 이용한 운행선로의 침목방진패드 스프링강성을 추정할 수 있는 기법을 제시하고자 한다.

Effect of modifying the thickness of the plate at the level of the overlap length in the presence of bonding defects on the strength of an adhesive joint

  • Attout Boualem;Sidi Mohamed Medjdoub;Madani Kouider;Kaddouri Nadia;Elajrami Mohamed;Belhouari Mohamed;Amin Houari;Salah Amroune;R.D.S.G. Campilho
    • Advances in aircraft and spacecraft science
    • /
    • 제11권1호
    • /
    • pp.83-103
    • /
    • 2024
  • Adhesive bonding is currently widely used in many industrial fields, particularly in the aeronautics sector. Despite its advantages over mechanical joints such as riveting and welding, adhesive bonding is mostly used for secondary structures due to its low peel strength; especially if it is simultaneously exposed to temperature and humidity; and often presence of bonding defects. In fact, during joint preparation, several types of defects can be introduced into the adhesive layer such as air bubbles, cavities, or cracks, which induce stress concentrations potentially leading to premature failure. Indeed, the presence of defects in the adhesive joint has a significant effect on adhesive stresses, which emphasizes the need for a good surface treatment. The research in this field is aimed at minimizing the stresses in the adhesive joint at its free edges by geometric modifications of the ovelapping part and/or by changing the nature of the substrates. In this study, the finite element method is used to describe the mechanical behavior of bonded joints. Thus, a three-dimensional model is made to analyze the effect of defects in the adhesive joint at areas of high stress concentrations. The analysis consists of estimating the different stresses in an adhesive joint between two 2024-T3 aluminum plates. Two types of single lap joints(SLJ) were analyzed: a standard SLJ and another modified by removing 0.2 mm of material from the thickness of one plate along the overlap length, taking into account several factors such as the applied load, shape, size and position of the defect. The obtained results clearly show that the presence of a bonding defect significantly affects stresses in the adhesive joint, which become important if the joint is subjected to a higher applied load. On the other hand, the geometric modification made to the plate considerably reduces the various stresses in the adhesive joint even in the presence of a bonding defect.

이중관으로 연결된 자연치와 임플랜트의 악골 내응력분포에 관한 3차원 유한요소분석 (A Study on the Stress Distribution of Tooth/Implant Connected with Konus Telescope Denture Using 3-Dimensional Finite Element Method)

  • 이수옥;최대균;권긍록;우이형
    • 대한치과보철학회지
    • /
    • 제46권4호
    • /
    • pp.381-395
    • /
    • 2008
  • 서론: 임플랜트와 자연치가 혼재하는 증례에 있어서 보철물의 연결여부와 그 결과에 대해 논쟁이 많이 있어왔다. 본 연구에서는 임플랜트와 자연치를 이중관으로 연결했을 경우에 치조골 내에서 임플랜트 주변 조직의 하중에 대한 응력의 양상을 분석하여 자연 치와 임플랜트를 연결하는 보철물의 방식으로 이중관의 사용가능성을 비교 분석하고자 하였다. 연구재료 및 방법: 실험군으로서는 한 개의 자연치아와 한 개의 임플랜트를 이중관으로 연결한 모델을 사용하고, 대조군으로 두 개의 임플랜트를 이중관으로 연결한 모델을 사용하여 응력분포의 차이를 비교분석하였다. 임플랜트와 자연치 각각 통상적인 금합금을 사용하여 내,외관을 모델링한 후 각각의 외관을 연결한 구조를 형성한다. 내관과 외관은 $6^{\circ}$의 경사도 (taper)로 형성하고 내외관 사이의 공간을 인정하여, 하중을 가했을 경우 내외관 사이에 합금의 마찰력과 변위를 계산하여, 그 결과가 골내 임플랜트의 응력 분포에 영향을 줄 수 있도록 실험모델을 형성하였다. 하중은 일반적인 저작력의 크기로 임플랜트 부위와 자연치 부위에 각각 부여하여 응력의 차이를 비교하였으며, 최종적으로 대조군과의 차이를 비교하고자 하였다. 하중 조건은 임플랜트와 자연치아 사이, 임플랜트 후방, 전방부 자연치아의 사이에 하중을 가하였으며, Load case 1, 2, 3로 분류하였다. 실험을 위해서 ANSYS Version 8.1을 사용하여 3차원 유한 요소 분석을 시행하였다. 연구성적: 지대주에 가해지는 응력 - 최대 응력은 실험군이 대조군에 비해 각각 2.4 (Load Case1), 1.02 (Load Case2), 0.46 (Load Case3)배 높게 나타났으며,Load Case3의 경우 대조군의 최대응력이 더욱 높게 나타났는데, 이것은 전치부 치아의 경우 하중을 효과적으로 분산하기 때문인 것으로 보인다. 악골내에 가해지는 응력 - 악골내에서 최대 응력은 실험군이 대조군에 비해 각각 1.69 (Load Case1), 1.26 (Load Case2), 1.93 (Load Case3) 배 크게 나타났다 결론: 임플랜트와 자연치아가 혼재된 경우, 이중관으로 상부 보철물을 연결한 경우, 임플랜트만으로 구성된 대조군에 비해, 임플랜트 및 임플랜트 식립 부위의 악골에 가해지는 최대 응력은 1.02 에서 2.4배 가량 높게 나타났다. 자연치아의 움직임이 임플랜트의 경부 및 악골 내에 응력을 집중시킨 것으로 보이며, 자연치아가 있는 악골부위에서는 응력이 거의 발생하지 않았다. 본 실험의 경우 자연치아와 임플랜트에서 모든 하중을 받는 것으로 계산하였고, 실제 가철성 의치에서 나타나는 점막의 영향을 배제하였기 때문에, 실제 임상에서는 후방 임플랜트에 미치는 최대 응력의 크기는 2.4배 이하일 것으로 추정할 수 있다. 앞으로, 실제 임상에 적용하였을 경우, 임플랜트에 미치는 영향 및 자연치아에 미치는 결과에 대한 연구가 필요할 것이다.

DRASTIC 모델 및 지하수 수치모사 연계 적용에 의한 부여읍 일대의 지하수 오염 취약성 평가 (A Joint Application of DRASTIC and Numerical Groundwater Flow Model for The Assessment of Groundwater Vulnerability of Buyeo-Eup Area)

  • 이현주;박은규;김강주;박기훈
    • 한국지하수토양환경학회지:지하수토양환경
    • /
    • 제13권1호
    • /
    • pp.77-91
    • /
    • 2008
  • 본 연구에서는 지하수 오염 취약성 평가 방법 중 가장 활용도가 높은 DRASTIC 기법 및 지하수 수치모사 기법을 병행하여 충청남도 부여군 부여읍 일대의 지하수 오염 취약성 변화를 평가하였다. 특히, DRASTIC 기법에서 사용되는 7가지 수리인자들 중 취약성에 대한 민감도가 상대적으로 높고 비교적 장기간 자료의 획득이 용이하지 않은 '지하수위(depth to water)' 인자에 대하여 수치모사 기법을 적용하여 시-자료에 대한 동적 보간을 실시하였으며, 이를 통하여 시간에 따른 지하수위 변화가 지하수 오염 취약성에 미치는 영향을 검토하였다. 연구의 대상 지역은 부여군 부여읍 일대로 충적 대수층을 대상으로 기 실시된 지하수 수위자료를 보정대상으로 하였으며 연구지역의 지하수위 분포양상을 모사하기 위하여 유한차분모델인 MODFLOW를 이용하였다. 지하수 수치모사에는 정류 및 부정류 모사가 모두 실시되었다. 그 결과 기존 지구통계기법을 통한 지하수위 보간에서 가장 큰 가중치를 갖는 지형적 요인뿐만 아니라 유역 암상이나 지질특성 등과 같은 제반 수리지질학적 요인까지 복합적으로 반영하는 지하수위 자료를 얻을 수 있었다. 기존 연구에서 수위자료를 보완할 때 주로 이용되는 기법인 크리깅은 수치모사와는 달리 강우나 계절적 영향, 암상이나 지질등과 같은 요인들은 반영하지 못한다. 또한 수치모사를 통하여 얻어진 수위자료는 부정류 모사가 가능함으로 인하여 풍수기와 갈수기 그리고 평수기의 지하수위 변화 등도 파악할 수 있다는 장점이 있다. 모사를 통하여 얻어진 지하수위 자료와 공동크리깅 방법을 통해 얻어진 지하수위 자료를 입력자료로 하여 연구지역에 지하수 오염 취약성을 비교 평가한 결과 평수기의 경우 유사한 경향의 지하수 오염 취약성을 보였다. 또한 공동크리깅을 통해 분석이 어려운 분기별 지하수 오염 취약성 평가 결과 강우량이 많은 여름철을 포함한 분기와 강우량이 가장 적은 겨울철을 포함한 분기가 뚜렷한 취약성 차이를 보이는 것으로 분석되었다. 본 연구에서 제시하는 바와 같이 지하수 오염 취약성 평가에 있어 가장 큰 민감도를 갖는 수위자료의 수가 공간적 혹은 시간적으로 부족한 경우 수치모사 적용을 통해 자료의 한계성을 보완할 수 있으며 시간적인 보간이 이루어질 수 있어 유용할 수 있을 것이라 판단된다.

프로토타입 R/F 흉부 디지털 단층영상합성장치 시스템에서 잘림 아티팩트 감소를 위한 가중 정규화 접근법에 대한 연구 (Truncation Artifact Reduction Using Weighted Normalization Method in Prototype R/F Chest Digital Tomosynthesis (CDT) System)

  • 손준영;최성훈;이동훈;김희중
    • 한국방사선학회논문지
    • /
    • 제13권1호
    • /
    • pp.111-118
    • /
    • 2019
  • 흉부 디지털 단층영상합성장치는 기존 DR의 낮은 깊이 해상도, CT의 높은 피폭선량 문제를 해결할 수 있는 획기적인 영상장치로 대두되고 있다. 그러나 제한된 스캔 각도로 인해 프로젝션이 X 선 소스 동작 방향으로 흉부를 완전히 포함 할 수 없어 재구성 된 슬라이스의 위, 아래 방향 경계를 따라 강도의 불연속성이 발생하게 되는데 이러한 현상을 잘림 아티팩트 (Truncation artifact)라고 한다. 이 연구의 목적은 가중 정규화 접근법을 사용하여 잘림 아티팩트를 줄이고 리스템에서 개발한 프로토 타입 흉부 디지털 단층영상합성장치 시스템에 대한 이 접근법의 성능을 평가하는 것이다. 이 시스템의 source-to-image distance는 1100 mm 이고 X 선원의 회전 중심은 검출기 표면에서 100mm 위로 설정되었다. LUNGMAN 팬텀을 사용하여 ${\pm}20^{\circ}$의 투영 뷰를 $1^{\circ}$ 간격으로 41장을 얻은 후, filtered back projection 알고리즘으로 재구성했다. 정량적 평가를 위하여 시뮬레이션을 이용하여 기준영상을 재구성 후 peak signal to noise ratio와 structure similarity index 값을 평가하였으며 실제 실험 데이터를 이용하여 mean value of specific direction 값을 평가하였다. 시뮬레이션 결과로 아티팩트 보정 전 일반적인 filtered back projection 알고리즘으로 재구성 한 영상과 비교하여 peak signal to noise ratio값과 structure similarity index값 모두 각각 증가하였으며, 실제 실험 재구성 영상의 mean value of specific direction 결과는 아티팩트의 영향이 감소됨을 확인할 수 있었다. 결론적으로, 가중 정규화 방법은 잘림 아티팩트를 줄임으로써 진단의 어려움을 발생시키는 가능성을 개선시킬 수 있는 방법으로 사료된다.

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

개조된 케이슨 플로팅 도크의 구조 보강에 대한 연구 (A Study on the Structural Reinforcement of the Modified Caisson Floating Dock)

  • 김홍조;서광철;박주신
    • 해양환경안전학회지
    • /
    • 제27권1호
    • /
    • pp.172-178
    • /
    • 2021
  • 선박 수리시장은 선박에 의한 환경오염 방지 강화, 선박구조에 대한 안전기준 강화 등의 영향으로 유지 및 보수에 관한 관심이 꾸준히 증대되고 있다. 이러한 영향을 반영하여 서남해에 있는 수리 조선사들에 외국 선사들의 수리 요청 접수 건수가 증가하고 있다. 그러나, 서남해권 수리 조선사들은 영세한 중소업체가 대부분이라서 수리조선 업체의 통합적 시너지 효과로 이어지기가 쉽지 않고, 집적화가 되어있지 않아서 인프라 공동활용이 어려워서 수리조선업 활성화에 걸림돌로 작용하고 있다. 수리조선업을 운영하기 위해서는 플로팅 도크가 필수적으로 필요로 하며, 대부분 노후화된 케이슨 도크를 해외로부터 수입한 후, 개/보수를 통하여 운용하고 있다. 그러나, 사용 수명이 30년 이상이고, 구조물 검사 기준이 없어서 안전분야에 취약성을 갖고 있다. 본 연구에서는 개조된 케이슨 도크의 구조 안전성을 평가하고, 도출된 문제점을 해결하기 위하여 추가적인 구조 보강안을 찾기 위하여 유한요소해석 프로그램인 ANSYS를 활용하였다. 플로팅 도크의 경우, 선급 규정이 있지만 구조강도 관련해서는 규정이 미흡하여 적용성이 떨어지고 있는 실정이다. 이러한 부족한 평가 영역에 대해서는 상세 구조해석을 통하여 보완하였다. 보강안은 수리조선소 작업의 특성을 고려하여 폰툰 갑판 상부 보강과 선측 탱크 보강으로 결정하였다. 결정안에 대한 구조해석을 통하여 선측 보강안을 최종안을 선정하였고, 실제 구조물을 제작하여 보강안을 반영하였다. 도출된 주요 결과들은 유사 설비의 구조 강도 개선을 위한 참고 자료로 활용 가능하며, 개/보수 시 이러한 방법을 활용하면 빨리 최적 해를 찾을 수 있을 것으로 기대된다.

대기 누출 방사성물질 선원 위치 추적을 위한 3차원 궤적모델 개발 (Development of Three-Dimensional Trajectory Model for Detecting Source Region of the Radioactive Materials Released into the Atmosphere)

  • 서경석;박기현;민병일;김소라;양병모
    • Journal of Radiation Protection and Research
    • /
    • 제41권1호
    • /
    • pp.31-39
    • /
    • 2016
  • 연구배경: 우리나라를 포함한 중국, 대만, 북한, 일본 등에서 원전, 재처리시설과 같은 원자력시설의 증가에 따라 주변국 핵활동 분석의 종합적 대책이 필요하다. 우리나라와 포괄적핵실험금지조약기구(Comprehensive Nuclear-Test-Ban Treaty Organization, CTBTO)는 동북아시아 지역에서 핵종 탐지소를 운영 중으로, 핵종탐지 장비에서 특이 값 측정시 모니터링 자료의 분석과 더불어 배출원 탐색모델을 이용하여 핵종의 기원이 어디인지 추정하고 평가하는 것은 주변국 핵활동에 대한 감시 및 안전성 확보 측면에서 중요하다. 재료 및 방법: 주변국의 은밀한 핵활동 시 방사성핵종의 기원을 추정하기 위하여 3차원 전진/후진형 궤적모델을 개발하였다. 개발된 궤적모델은 궤적 미분방정식을 유한차분법을 이용한 방법으로 주어진 바람자료를 이용하여 방사성핵종의 방출지점으로부터 입자의 궤적을 순차적으로 찾아가는 전진형 모델과 시간 역산으로 방출기원을 추정하는 후진형 모델로 구성되었다. 결과 및 논의: 개발된 궤적모델의 검증을 위하여 체르노빌 사고 당시 측정된 농도자료를 이용하였다. 검증결과 관측지점의 농도가 높게 측정된 지점과 방출기원에서 가까운 지역으로부터 시간 역산의 방출지점을 추정한 결과의 정확도가 높았다. 3차원 궤적모델은 방출시간, 방출높이, 방출간격 등의 변수에 의해 계산결과가 달라지는 불확도를 내포하고 있는데, 이러한 궤적모델의 불확도를 최소화하기 위해 한국원자력연구원에서 개발한 대기확산모델(long-range accident dose assessment system, LADAS)를 이용하여 fields of regards (FOR) 기법에 의해 오염물 방출영역을 추정한바 신뢰성 있는 결과를 얻었다. 결론: 본 연구를 통하여 개발된 배출원 탐색모델은 주변국의 은밀한 핵활동 시 핵종 탐지장비와 연계하여 방사성핵종의 방출지역과 기원을 파악하여 우리나라의 핵종탐지 능력을 향상하고 핵활동 및 방사선 안전 분야에서 주도적 역할을 할 수 있을 것으로 생각된다.

옥천화강암과 보은화강암 관입에 의한 접촉변성작용 (The Contact Metamorphism Due to the Intrusion of the Ogcheon and Boeun granites)

  • 오창환;김창숙;박영도
    • 암석학회지
    • /
    • 제6권2호
    • /
    • pp.133-149
    • /
    • 1997
  • 옥천화강암 주변의 변성니질암에서는 관입체에 가까워 지면서 변성정도가 흑운모대에서 홍주석대를 거쳐 규선석대로 증가한다. 그리고 보은화강암 남서부 주변 변성염기성암에서는 관입체에 가까워지면서 녹색편암상-각섬암상 점이대에서 각섬암상대를 거쳐 상부각섬암상대로 변성정도가 증가한다. 보은화강암 근접부인 도이리 부근에서 규선석과 홍주석을 포함하는 광물군이 화강암 경계부로부터 500m 이내에 위치한 변성니질암에서 관찰된다. 이러한 연구결과는 보은화강암과 옥천화강암이 주변지역에 저압형의 접촉 변성작용을 야기시켰음을 지시한다. 옥천화강암 주변에서 석류석을 포함하는 홍주석대 시료로부터 계산된 온도-압력 조건은 $540{\pm}40^{\circ}C, 2.8{\pm}0.9$ kb이며 보은화강암 주변 상부각섬암상대의 변성염기성암으로부터 추정된 접촉변성 작용시의 온도는 $698{\pm}28^{\circ}C$이다. 보은화강암 주변의 변성염기성에서 나타나는 각섬석과 사장석의 성분은 매우 넓은 범위를 갖는데 이는 각섬석과 사장석의 불혼합구간과 불완전한 변성반응에 의해 잔류된 불안정한 성분에 기인한다. 변성정도가 증가하면서 안정된 각섬석과 사장석의 성분범위는 불혼합구간의 소멸에 의해 줄어든다. 보은화강암과 옥천화강암에 의한 접촉변성 작용시 주변암에 야기된 열적변화는 2차원의 유한차분법을 이용한 CONTACT2 프로그램을 이용하여 이론적으로 계산되었다. 이론적인 계산을 위해 보은화강암의 경우 화강암의 변의 길이가 20 km인 삼각형을 그리고, 옥천화강암의 경우 지름이 10 km인 원을 관입형태로 가정하였다. 이론적인 계산결과를 화강암 주변의 변성암에서 계산된 결과와 비교하여 보면 예상되는 관입온도는 옥천화강암 경우는 $800^{\circ}C$ 그리고 보은화강암의 경우는 $1000^{\circ}C$ 보다 높았을 것으로 추정된다. 하지만 옥천 및 화강암 관입시기에 연구지역의 지온 상승률이 일반 대륙지각에서의 지온 상승률보다 높았을 경우 관입온돈느 $800^{\circ}C$$1000^{\circ}C$ 보다 낮을 수 있다.

  • PDF