• 제목/요약/키워드: model averaging

검색결과 300건 처리시간 0.033초

Stochastic finite element based seismic analysis of framed structures with open-storey

  • Manjuprasad, M.;Gopalakrishnan, S.;Rao, K. Balaji
    • Structural Engineering and Mechanics
    • /
    • 제15권4호
    • /
    • pp.381-394
    • /
    • 2003
  • While constructing multistorey buildings with reinforced concrete framed structures it is a common practice to provide parking space for vehicles at the ground floor level. This floor will generally consist of open frames without any infilled walls and is called an open-storey. From a post disaster damage survey carried out, it was noticed that during the January 26, 2001 Bhuj (Gujarat, India) earthquake, a large number of reinforced concrete framed buildings with open-storey at ground floor level, suffered extensive damage and in some cases catastrophic collapse. This has brought into sharp focus the need to carry out systematic studies on the seismic vulnerability of such buildings. Determination of vulnerability requires realistic structural response estimations taking into account the stochasticity in the loading and the system parameters. The stochastic finite element method can be effectively used to model the random fields while carrying out such studies. This paper presents the details of stochastic finite element analysis of a five-storey three-bay reinforced concrete framed structure with open-storey subjected to standard seismic excitation. In the present study, only the stochasticity in the system parameters is considered. The stochastic finite element method used for carrying out the analysis is based on perturbation technique. Each random field representing the stochastic geometry/material property is discretised into correlated random variables using spatial averaging technique. The uncertainties in geometry and material properties are modelled using the first two moments of the corresponding parameters. In evaluating the stochastic response, the cross-sectional area and Young' modulus are considered as independent random fields. To study the influence of correlation length of random fields, different correlation lengths are considered for random field discretisation. The spatial expectations and covariances for displacement response at any time instant are obtained as the output. The effect of open-storey is modelled by suitably considering the stiffness of infilled walls in the upper storey using cross bracing. In order to account for changes in soil conditions during strong motion earthquakes, both fixed and hinged supports are considered. The results of the stochastic finite element based seismic analysis of reinforced concrete framed structures reported in this paper demonstrate the importance of considering the effect of open-storey with appropriate support conditions to estimate the realistic response of buildings subjected to earthquakes.

Generation of Ionospheric Delay in Time Comparison for a Specific GEO Satellite by Using Bernese Software

  • Jeong, Kwang Seob;Lee, Young Kyu;Yang, Sung Hoon;Hwang, Sang-wook;Kim, Sanhae;Song, Kyu-Ha;Lee, Wonjin;Ko, Jae Heon
    • Journal of Positioning, Navigation, and Timing
    • /
    • 제6권3호
    • /
    • pp.125-133
    • /
    • 2017
  • Time comparison is necessary for the verification and synchronization of the clock. Two-way satellite time and frequency (TWSTFT) is a method for time comparison over long distances. This method includes errors such as atmospheric effects, satellite motion, and environmental conditions. Ionospheric delay is one of the significant time comparison error in case of the carrier-phase TWSTFT (TWCP). Global Ionosphere Map (GIM) from Center for Orbit Determination in Europe (CODE) is used to compare with Bernese. Thin shell model of the ionosphere is used for the calculation of the Ionosphere Pierce Point (IPP) between stations and a GEO satellite. Korea Research Institute of Standards and Science (KRISS) and Koganei (KGNI) stations are used, and the analysis is conducted at 29 January 2017. Vertical Total Electron Content (VTEC) which is generated by Bernese at the latitude and longitude of the receiver by processing a Receiver Independent Exchange (RINEX) observation file that is generated from the receiver has demonstrated adequacy by showing similar variation trends with the CODE GIM. Bernese also has showed the capability to produce high resolution IONosphere map EXchange (IONEX) data compared to the CODE GIM. At each station IPP, VTEC difference in two stations showed absolute maximum 3.3 and 2.3 Total Electron Content Unit (TECU) in Bernese and GIM, respectively. The ionospheric delay of the TWCP has showed maximum 5.69 and 2.54 ps from Bernese and CODE GIM, respectively. Bernese could correct up to 6.29 ps in ionospheric delay rather than using CODE GIM. The peak-to-peak value of the ionospheric delay for TWCP in Bernese is about 10 ps, and this has to be eliminated to get high precision TWCP results. The $10^{-16}$ level uncertainty of atomic clock corresponds to 10 ps for 1 day averaging time, so time synchronization performance needs less than 10 ps. Current time synchronization of a satellite and ground station is about 2 ns level, but the smaller required performance, like less than 1 ns, the better. In this perspective, since the ionospheric delay could exceed over 100 ps in a long baseline different from this short baseline case, the elimination of the ionospheric delay is thought to be important for more high precision time synchronization of a satellite and ground station. This paper showed detailed method how to eliminate ionospheric delay for TWCP, and a specific case is applied by using this technique. Anyone could apply this method to establish high precision TWCP capability, and it is possible to use other software such as GIPSYOASIS and GPSTk. This TWCP could be applied in the high precision atomic clocks and used in the ground stations of the future domestic satellite navigation system.

한 시멘트공장의 분진발생과 대기확산에 관한 조사연구 (A Study on the Emission and Dispersion of Particulate Matter from a Cement Plant)

  • 장남익;정용;권숙표
    • Journal of Preventive Medicine and Public Health
    • /
    • 제16권1호
    • /
    • pp.67-77
    • /
    • 1983
  • To investigate the an air pollution by particulate matter and its dispersion, a cement plant produceing portland cement 600,000 ton/year and its vicinity were surveyed from Obtober, 1980 to April, 1983. The survey was mainly focused on main stack emmission rate of the cement plant and particle size distribution in the dust, dustfall and total suspended particulate concentration in the area by month and distance from the stack. The results of the study were as follows; 1. The main stack emission rate was surveyed before and after the spray tower was additionally installed to the original E.P bag filter. Before the spray tower installed, the main stack emission rate was higher ($0.64g/Nm^3$) than the emission standard of Korean Environmental Preservation Law's ($0.59g/Nm^3$, amended to $0.4g/Nm^3$ on April 1983), but after the spray tower was installed, its main stack emission rate was markedly decreased to the standard ($0.43g/Nm^3$). 2. $2{\sim}3{\mu}m$ of the particle size was the largest portion (20.8%) of the dust particulate from the main stack and 50% of the frequency distribution was $1.5{\mu}m$ of the size. Most particle size was below $10{\mu}m$. 3. The spray tower reduced the dustfall to $37.81{\sim}9.76\;ton/km^2/month$ while dustfall appeared at $45.29-15.45ton/km^2/month$, in the vicinity of plant before spray tower installed 4. Mean concentrations of total suspended particulate for 24 hours of the various stations were determined in $20.6-200.0{\mu}g/m^3$, 3 stations of tham were higher than the value of Harry and William's arthmetic average standard $130{\mu}g/m^3$. 5. Linear regression between dustfall [X] and total suspended particulate[Y] concentration was an equation, Y=4.024X+11.479.[r=0.91] 6. During the whole seasons in the opposite area 100m apart from the omission source the prevailing wind direction was with estimated more than $30ton/km^2/month$, and the concentration of total suspended particulate for 24 hours averaging time was more than $140{\mu}g/m^3$ in the same area and direction. 7. Assuming the wind direction were constant through the day dustfalls for a day were estimated at $13.40ton/km^2/day,\;10.79ton/km^2/day$ and $4.55ton/km^2/day$ at various distances of 100m, 500m and 1,500m from the emission source respectively. 8. In the simutalion of dustfall and suspended dust by area, Gaussian dispersion model modified by size distribution of particulate matter was not applicated since the emission of dust were from multi sources other them stack. From the above results, it could be applied that the dispersion of dust from the cement plant is estimated and regulated for the purpose of environmental protection.

  • PDF

지역의 사회·경제적 인자와 용수공급체계를 고려한 가뭄 위험도 평가 (Drought risk assessment considering regional socio-economic factors and water supply system)

  • 김지은;김민지;최시중;이주헌;김태웅
    • 한국수자원학회논문집
    • /
    • 제55권8호
    • /
    • pp.589-601
    • /
    • 2022
  • 가뭄은 자연적 현상이지만, 지역의 물리적 및 사회적 요소와 결합되어 피해가 발생한다. 특히, 각종 용수 공급 및 수요과 연관되어 사회 경제적으로 큰 피해를 야기시킨다. 비슷한 심도의 기상학적 가뭄에도 지역의 특성과 용수공급체계에 따라 실제로 발생하는 가뭄 피해는 다르다. 본 연구에서는 지역의 사회·경제적 인자와 용수공급체계를 고려하여 가뭄 위험도를 평가하였다. 노출성은 용수공급 과부족량을 나타내는 결합가뭄관리지수(JDMI)를 등급화하여 평가하였다. 취약성은 가뭄에 영향을 받는 10개의 사회·경제적 인자에 엔트로피, PCA 및 GMM를 적용하여 가중평균하여 평가하였다. 대응능력은 지역의 용수능력을 나타내는 인자들을 베이지안 네트워크에 적용하여 평가하였다. 위험도는 노출성, 취약성 및 대응능력을 통합하여 결정하였다. 용수공급 실패 사상의 발생 가능성을 의미하는 가뭄 노출성을 평가한 결과, 괴산군이 0.81로 가장 높게 나타났다. 가뭄 취약성의 경우, 대전광역시가 0.61로 매우 취약한 것으로 나타났다. 지역의 용수공급체계가 고려된 가뭄 대응능력을 평가한 결과, 세종시가 가뭄 대응능력이 가장 낮은 것으로 나타났다. 마지막으로 위험도를 평가한 결과, 청주시가 가장 높게 나타났다. 이러한 결과를 통해 가뭄에 대한 위험 및 취약 원인을 파악하였으며, 향후 지역의 특성을 고려한 가뭄 피해 저감 정책 마련이 가능하다.

지지력 향상을 위한 확장형 강관말뚝에 관한 연구 (Development of Expandable Steel Pipe Piles to Improve Bearing Capacity)

  • 김의석;김정훈;김지윤;민병찬;최항석
    • 한국지반환경공학회 논문집
    • /
    • 제22권12호
    • /
    • pp.5-13
    • /
    • 2021
  • 기존의 뜬구조공법의 문제점을 개선하여 지하층 리모델링, 증축 시 안정성 확보와 공사비 절감이 가능한 공법을 개발하기 위하여 일반적인 마이크로파일에 비해 경제적이고 안정성 높은 확장형 강관말뚝을 개발하였다. 확장형 강관말뚝은 지중에서 강관을 확장하여 강관의 성능을 향상하는 공법으로 본 논문에서는 확장형 강관말뚝의 강관 형상에 따른 좌굴강도의 변화를 파악하고, 수치해석 모델을 개발하여 강관확장으로 인한 요철부의 턱효과를 규명하고, 재료시험을 통하여 최적 강관 확관량을 산정하였다. 강관의 확장 직경이 클수록 확관 턱 개수가 많을수록 좌굴강도가 커짐을 알 수 있었으며, 수치해석 결과에 따르면 확관률보다 확관 턱 개수가 좌굴강도에 큰 영향을 미침을 알 수 있었으며, 확관률는 1.2배 이상일 때, 확관 턱 개수는 증가할수록 좌굴강도 증가 효과가 크게 발생함을 알 수 있었다. 또한, 확장 각도가 45° 이하이고, 확관률이 1.3배 이상일 때 요철부의 턱효과가 크게 발생하는 것을 알 수 있었다. 강관이 파단되는 항복 시, 신율은 20~32%로 평균 25.4% 수준으로 확인되어 그 이상 변형은 강관의 성능을 발휘할 수 없었다. 재료시험을 통해 강관의 성능을 발휘하기 위한 최대 확관량은 항복 시 최저값으로 확인된 신율 20%에 안전측 80%를 고려하여 16%로 제한하는 것이 바람직한 것으로 분석되었다.

슈퍼커패시터를 적용한 에너지 저장시스템 설계에 관한 연구 (The Research On the Energy Storage System Using SuperCapacitor)

  • 김일송
    • 예술인문사회 융합 멀티미디어 논문지
    • /
    • 제8권11호
    • /
    • pp.215-222
    • /
    • 2018
  • 본 논문에서는 슈퍼커패시터를 적용한 에너지 저장장치에 대해서 연구하였다. 납축전지등을 사용하였던 기존방식과는 다르게 슈퍼커패시터를 사용하여 고출력을 달성할 수 있으며, 이러한 특징으로 인하여 연구된 저장장치는 비상전원용 전력을 공급할 수 있을 뿐만 아니라, 전력품질 개선이나 짧은 순간에 큰 출력을 요구하는 고속 기동용 모터등의 전원시스템으로도 사용이 가능한 장점이 있다. 에너지변환 시스템은 양방향 인버터와 컨버터로 구성되어 고속, 고출력의 충전, 방전을 수행할 수 있는 장점도 있다. 양방향 인버터를 이용한 two loop 제어기 설계 방법이 제시되었고, 실험과 제작을 통해 설계 방법의 정합성을 입증하였다. 설계 순서는 먼저 상태공간 평균화 기법과 Decoupling 기법을 이용하여 선형화된 전달함수를 계산하였다. 다음 20% 오버슈트와 안정시간을 만족시키는 전류 제어기 설계를 한 후에 오버슈트 없는 전압과 10배 이상 느린 안정시간을 가지는 전압제어기를 설계한 후 단위계단응답 커브를 통해 원하는 설계특성이 나오는 지를 확인하는 과정으로 이루어진다. 설계된 제어기들은 계통과 단위 역률을 이루며 제어되기 때문에 전력품질 향상을 이룰 수 있으며, 빠른 응답특성과 정상상태 오차를 0으로 만들 수 있는 우수한 제어특성을 가지고 있음을 확인할 수 있었다.

식생된 개수로에서 항력가중계수가 흐름에 미치는 영향 분석 (Impact of Drag-Related Weighting Coefficients in Vegetated Open-Channel Flows)

  • 강형식;최성욱
    • 대한토목학회논문집
    • /
    • 제26권5B호
    • /
    • pp.529-537
    • /
    • 2006
  • 본 논문에서는 식생된 개수로 흐름의 수치모의에 필요한 항력가중계수의 영향을 분석하였다. 이를 위해 시간 및 공간 평균기법을 이용하여 식생된 개수로 흐름에서 레이놀즈응력의 수송방정식을 유도하였다. 그 결과 총 레이놀즈응력은 시간의 변동 성분에 의한 레이놀즈응력과 공간상의 변동 성분에 의한 레이놀즈응력의 합이며, 총 레이놀즈응력의 수송방정식을 수치모의하기 위한 항력가중계수의 값은 $C_{fk}$ = 1.0인 것으로 나타났다. 그러나 시간의 변동 성분에 의한 레이놀즈응력을 수치모의하기 위해서는 거의 영에 가까운 항력가중계수를 갖는 것으로 나타났다. 이는 과거의 수치모의 연구에서 항력가중계수의 값이 거의 영에 가까울 때 실험결과와 잘 일치했는지에 대한 중요한 이유이다. 즉, 공간상의 변동성분에 의한 레이놀즈응력의 값은 실험을 통해 측정하기 매우 어렵기 때문에 식생된 개수로 흐름에서 측정된 레이놀즈응력은 대부분 시간상의 변동성분에 의한 레이놀즈응력이기 때문이다. 또한 레이놀즈응력모형을 이용하여 항력가중계수에 따른 식생된 개수로 흐름을 수치모의하고 기존의 실험 결과와 비교하였다. 그 결과 평균유속과 레이놀즈응력의 경우 항력가중계수의 영향은 작은 것으로 나타났으나, 난류강도 분포에서는 항력가중계수의 영향이 매우 크게 발생하였다. 또한 총 레이놀즈응력과 시간의 변동성분에 의한 레이놀즈응력의 수송방정식에서 각 항의 수지분석을 통하여 항력가중계수가 난류강도에 미치는 영향을 분석하였다.

가우시안 군집분석을 이용한 천리안 위성의 대기운동벡터 표적추적 알고리듬 개발 및 분석 (Development and Analysis of COMS AMV Target Tracking Algorithm using Gaussian Cluster Analysis)

  • 오유림;김재환;박형민;백강현
    • 대한원격탐사학회지
    • /
    • 제31권6호
    • /
    • pp.531-548
    • /
    • 2015
  • 위성영상을 이용하여 산출된 대기운동벡터(AMV)와 라디오존데 바람 관측 자료를 이용한 검증결과는 산출된 AMV가 지속적으로 관측 자료에 비해서 풍속이 약하게 나타나는 Slow Speed Bias(SSB)를 보여 주었다. 이러한 SSB는 표적추적, 표적선정, 그리고 고도할당 단계의 오차에 의해 야기될 수 있으며, 이 중 고도할당 단계의 오차는 SSB를 발생시키는 주된 요인으로 여겨진다. 그러나 최근 연구에서는 고도할당 단계의 개선만으로는 SSB 문제를 해결하는데 한계가 있음을 밝혔다. 그러므로 본 연구에서는 새로운 표적추적 알고리듬을 개발하여 SSB를 감소시킴으로서 기상청 현업 AMV 알고리듬의 성능을 개선하고자 하였다. 표적추적 단계의 오차는 표적 내에 다양한 시 공간 규모의 바람이 포함되어 벡터가 과도하게 평균된 움직임으로 계산되거나, 구름이 추적 시간동안 형태를 유지하지 못하고 변형되는 경우에 발생한다. 이러한 문제를 해결하기 위해 개발된 표적추적 알고리듬에서는 가우시안 군집분석(GMM)을 이용하여 변형이 적고 추적에 용이한 저온 군집을 표적으로 재선정하고, 이미지를 변형시켜 군집의 움직임을 보다 쉽게 추적할 수 있게 하였다. 또한 표적을 추적하기 위한 방법으로 거리제곱합 방법을 사용하였다. 개발된 알고리듬과 기존 COMS 알고리듬을 천리안 위성의 적외채널 영상에 적용하여 AMV를 산출하였으며, 이를 라디오존데 관측 자료와 비교 검증해 보았다. 제안된 알고리듬으로 산출된 AMV는 기존 알고리듬으로 산출된 AMV보다 평균 풍속이 $2.7ms^{-1}$증가함에 따라 SSB가 평균 29%까지 감소하는 개선된 결과를 보여주었다. 그러나 개발된 알고리듬으로 산출된 AMV는 중 하층의 정확도가 감소하였고, 기존 알고리듬에 비해 산출되는 AMV 벡터수가 약 40%까지 감소함을 보였다. 이에 따라 중 하층의 정확도 개선과 기존의 알고리듬과 비교하여 산출되는 벡터 개수가 감소하는 문제를 보완하기 위한 연구가 필요할 것으로 판단된다.

천부 굴절법 탄성파 탐사 자료의 정보정 (Statics corrections for shallow seismic refraction data)

  • Palmer Derecke;Nikrouz Ramin;Spyrou Andreur
    • 지구물리와물리탐사
    • /
    • 제8권1호
    • /
    • pp.7-17
    • /
    • 2005
  • 천부 탄성파 굴절법 탐사를 이용하여 굴절이 발생하는 지층의 속도를 산출하는 것은 ill-posed 문제이다. 계산된 시간 변수들에서의 작은 변화들이 이로부터 산출된 속도들에 커다란 수평적 변화를 가져올 수 있으며 이는 종종 역산 알고리듬의 인위적인 오차를 유발한다. 이러한 인위적인 오차들은 모델링을 통해 인지되거나 보정되지 않는다. 그러므로 만약 모델에 근거한 역산을 통해 정밀한 지하 굴절 모델을 얻고자 한다면 정확한 초기 모델이 필요하다. 탄성파속도에서 인위적인 오차의 원인은 일반적으로 불규칙한 굴절면에 있다. 대부분의 경우에 GRM 방법을 이용하면 불규칙한 굴절면을 다룰 수 있고 굴절면의 정밀한 초기 모델을 만들 수 있다. 하지만 지표에 매우 가까운 극천부 지역 또한 불규칙하다면 GRM 방법의 효능은 감소하고 풍화대 보정이 필요하다. 천부 불균질대에 대한 일반적인 보정방법들은 수평적 확장이 제한된 극천부지역의 불균질대의 경우 효과적이지 못하다. 이럴 경우 GRM 평활화 통계적 방법(Smoothing Statics Method; SSM)이 지층의 속도를 좀 더 정확하게 평가할 수 있는 간단하고 실용적인 방법이다. GRM SSM 방법은 제로 XY 값을 가지고 계산된 시간-심도값들로부터 실제 XY 값을 가지고 얻어진 시간-심도값들의 평균값을 빼줌으로써 평활화 정보정을 수행한다. 심도가 깊어질수록 시간-심도값들이 XY 값에 따라 크게 변하지 않으므로 이들의 평균값은 최적값과 훨씬 더 같아진다. 그러나 극천부의 불균질대에 대해 시간-심도값들은 XY 값들이 증가함에 따라 수평적으로 이동하고 평균화 과정을 통해 대폭 감소한다. 결과적으로, XY값들에 대해 평균화된 시간-심도단면도는 천부의 불균질대에 대한 보정에 효과적이다. 또한 제로 XY 값을 가지고 계산된 시간-심도값들은 천부 불균질대의 영향과 대상 굴절면에 대한 시간-심도값들의 합으로 주어지므로 그들의 차는 정보정을 위해 주시로부터 빼주어야 할 대략적인 값들을 제공한다. GRM SSM 방법은 결정론적인 풍화대에 대한 보정법이라기 보다는 평활화 과정이다. 이 방법은 수평적으로 확장이 매우 제한된 천부 불균질대에 대해 가장 효과적이다. 모델과 현장 적용 결과들을 통해 GRM SSM 방법을 이용하여 불규칙한 굴절면을 가진 지층들에 대해 좀 더 신뢰할 수 있는 정밀한 탄성파 속도를 산출할 수 있음을 보여주고 있다.

진해만의 빈산소 수괴 형성에 관한 수치실험 (Numerical Simulation of the Formation of Oxygen Deficient Water-masses in Jinhae Bay)

  • 최우정;박청길;이석모
    • 한국수산과학회지
    • /
    • 제27권4호
    • /
    • pp.413-433
    • /
    • 1994
  • 진해만은 빈산소로 인한 수산생물의 생산성 저하가 현저하다. 본 연구에서는 진해만의 수질환경개선 및 회복을 위해서 빈산소 수괴의 형성방지가 무엇보다도 중요하다고 생각되어 해수유동 및 물질순환 모델을 이용하여 빈산소 수괴 형성상태를 재현하였고, 주요 오염부하가 용존산소에 미치는 영향의 정도를 예측하여 어장환경 관리방안을 도출한 결과는 다음과 같다. 저층의 용존산소 농도분포는 수온 및 밀도성층과 밀접한 관계를 가지고 변동하였고 수온성층이 가장 강했던 $6{\sim}7$월에 양식어장이 밀집된 진해만 서부해역과 마산만 해역에서 2.0mg/l이하의 빈산소 수괴가 형성되었다. 해수유동 모델에 의한 $M_2$분조 계산결과 창조시 유향의 주류는 가덕수로를 통하여 서쪽으로 이동하였고 일부는 마산만으로 유입되었으며, 낙조시는 반대방향이었다. 최강유속은 가덕수로 부근에서 발생하였고 진동만, 고현성만, 원문만등의 진해만 서부해역과 마산만에서는 5cm/sec이하로 미약하였다. 항유는 만 중앙부와 칠천도 사이에서 반시계방향의 환유가 형성되는 것이 특징적이었으며, 마산만이나 당황만 입구 부근에서는 표층의 경우 남향의 흐름이 나타났고, 저층에서는 북향의 흐름이 형성 되었다. 물질순환 모델의 보정결과 상관성은 0.85이상을 보였고 상대오차는 $28\%$이하의 범위내에서 여름철의 빈산소 수괴를 재현하였다. 각종 오염부하가 용존산소에 미치는 영향의 정도와 범위를 시뮬레이션을 통하여 예측한 결과 SOD가 전지역에 걸쳐 가장 큰 영향을 나타내었고 ,마산만의 유입부하도 용존산소 분포에 큰 영향을 미쳤으나 마산만내 국한되어 나타났으며, 양식생물에 의한 부하는 영향도 적었고 진해만 서부해역에 국한되었다. 빈산소 수괴가 강하게 형성되는 마산만과 진해만 서부해역의 빈산소 수괴 형성방지를 위해 효율적인 오염저감 대책으로서 마산만의 경우 유입 COD부하와 SOD를 저감해야하며, 진해만 서부해역의 경우 SOD가 주가 되어야 함을 알 수 있었다. 회복시켜야 할 용존산소 농도를 해역II등급인 5.0mg/l의 농도를 유지시키기 위해서는 마산만의 경우 유입되는 COD 부하를 $50\%$, SOD를 $70\%$로 저감해야 하며, 진해만 서부해역의 경우 SOD를 $95\%$, 양식생물 오염부하를 $90\%$까지 저감해야 되는 것으로 나타났다.

  • PDF