• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.1-7
• /
• 2011

In this paper, based on an isogeometric approach, we have developed a shape design optimization method for plane elasticity problems subjected to design-dependent loads. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry. In an isogeometric analysis, however, the geometric properties are already embedded in the B-spline basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. The solution space for the response analysis can be represented in terms of the same NURBS basis functions to represent the geometry, which yields a precise analysis model that exactly represents the normal and curvature depending on the applied loads. A continuum-based isogeometric adjoint sensitivity is extensively derived for the plane elasticity problems under the design-dependent loads. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Journal of the Korean GEO-environmental Society
• /
• v.16 no.5
• /
• pp.35-42
• /
• 2015

This is conducted to observe underground temperature and to analyze its change affected by climate condition and soil infiltration in the mountainous area, Yesan region, Chungcheong-namdo province. Additionally, underground temperature change is also simulated using air temperature and soil thermal properties with a numerical model. Soil temperature monitoring data acquired from each depth, 20 cm, 50 cm, and 100 cm, indicates that the data within 50 cm in depth shows peak-shaped big fluctuation directly affected by air temperature and it at 100 cm has open-shaped small fluctuation. Underground temperature variation, a difference between high and low values, during monitoring period is weakly proportional to hydraulic conductivity of the sediment and it is assumed that water plays a part in delivering air temperature in soil. The underground temperature estimated by a numerical model is very similar to the observed data with an average value of 0.99 cross-correlation coefficient. From the result of this study, the aquifer unsaturated hydraulic conductivity of the soil and the groundwater recharge is likely to be able to estimate with underground temperature profile calculated using a numerical model.

• Journal of the Korean earth science society
• /
• v.36 no.3
• /
• pp.210-221
• /
• 2015

The terrain features and surface characteristics are the most important elements not only in meteorological modeling but also in air quality modeling. The diurnal evolution of local climate over complex terrain may be significantly controlled by the ground irregularities. Such topographic features can affect a thermally driven flow, either directly by causing changes in the wind direction or indirectly, by inducing significant variations in the ground temperature. Over a complex terrain, these variations are due to the nonuniform distribution of solar radiation, which is highly determined by the ground geometrical characteristics, i.e. slope and orientation. Therefore, the accuracy of prediction of regional scale circulation is strong associated with the accuracy of land-use and topographic information in meso-scale circulation assessment. The objective of this work is a numerical simulation using MM5-A2C model with the detailed topography and land-use information as the surface boundary conditions of the air flow field in mountain regions. Meteorological conditions estimated by MM5-A2C command a great influence on the dispersion of mountain areas with the reasonable feature of topography where there is an important difference in orographic forcing.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.337-343
• /
• 2014

This paper presents a continuum-based shape design sensitivity analysis(DSA) method for crack propagation problems using a reproducing kernel method(RKM), which facilitates the remeshing problem required for finite element analysis(FEA) and provides the higher order shape functions by increasing the continuity of the kernel functions. A linear elasticity is considered to obtain the required stress field around the crack tip for the evaluation of J-integral. The sensitivity of displacement field and stress intensity factor(SIF) with respect to shape design variables are derived using a material derivative approach. For efficient computation of design sensitivity, an adjoint variable method is employed tather than the direct differentiation method. Through numerical examples, The mesh-free and the DSA methods show excellent agreement with finite difference results. The DSA results are further extended to a shape optimization of crack propagation problems to control the propagation path.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.181-190
• /
• 2007

The linkage framework of surface geometric modeling based on the NURBS and shell finite element analysis is developed in this study. In the geometrically exact shell finite element analysis, the accuracy of the analysis strongly depends upon the accurate computation of the surface geometric quantities. Therefore if we obtain the necessary geometric quantities from the NVRBS surface equation, it's possible to construct the effective linkage framework of surface modeling in the CAD systems and shell finite element analysis using geometrically exact shell finite element. Besides, the linkage framework can be applied to the analysis of general and complex surfaces as well as simple surfaces. In this study, the shell surfaces are generated by interpolating given set of data points based on the NURBS surfaces. These data points usually can be obtained from surface scanning. But the representations of the generated NURBS surface are not same to one another. The accuracy depends on the chosen parameterization methods used in NURBS. Therefore, it is needed to select the suitable parameterization method according to the geometry of the surfaces. To verify the performance and accuracy of our developed linkage framework, we solve several well-known benchmark problems and assess the performance of the developed method.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1342-1347
• /
• 2007

최근 수자원 분야에서 GIS의 활용은 유역의 수리, 수문분석을 위한 모형의 입력자료 생성 및 모의 결과를 가시화, 이에 따른 유역관리 시스템의 구축 등 폭넓게 활용되고 있다. 또한 국가지리정보시스템에 조사를 통해서 수치지형도 및 주제도를 구축, 구축된 자료의 표준화를 실시하고 있는 실정이다. 그러나 수자원 분야에서 GIS를 활용하기 위한 기술은 주로 선진국을 중심으로 발전하여 왔기 때문에 우리의 실정에 맞지 않아 활용측면에 있어 신뢰할 만한 결과를 얻지 못하고 있는 실정이다. 따라서 국내 상황을 고려하면서 수자원이라는 전문분야 적용에 적합한 GIS(HyGIS)를 개발하고 여기에 수리, 수문분석모형을 연계하여 국내 실무분야에 적용함에 있어 편의성과 실용성을 구비한 모형개발이 이루어져야 한다. 따라서 본 연구에서는 국내 소프트웨어 GeoMania(HyGIS)에 의한 GIS 정보처리의 자동화를 기반으로 하여 실무에서 활용도가 높은 수문모형인 HEC-HMS 및 HEC-1과 수리모형인 HEC-RAS를 연계 및 통합하기 위해서 HyGIS에서 DLL형태로 제공되도록 하였다. HyGIS에서는 수문학적 DEM 분석 및 공간정보 생성, 선형참조가 가능한 하천 네트워크 생성, 유역 시설물 관리 등의 기능을 제공하며 COM(Component Object Model)을 기반으로 개발된 시스템으로 다른 소스로부터 개발된 컴포넌트를 연계하여 시스템의 기능 확장을 손쉽게 수행할 수 있도록 하였을 뿐만 아니라 공간 DB는 GeoMania의 고유 DB인 GSS를 이용한다. HyGIS-HMS는 HyGIS를 통한 국내 유역의 지리정보를 활용하여 HEC-HMS 뿐만 아니라 HEC-1을 추가하여 사용자의 기호와 편의에 따라 모형을 선택할 수 있도록 하였으며 HEC-1의 결과를 가시화하기 위해서 챠트 기능을 추가하였으며 매개변수를 자동으로 산정할 수 있도록 시스템을 구축하였다. HyGIS-RAS는 국내 하천유역에 대해서 기구축 되어있는 하천관리지리정보시스템(RIMGIS)자료를 직접 활용하도록 구성되어있고 자료를 활용하여 제내지와 제외지를 통합하여 TIN분석을 실시하여 범람 홍수해석에 활용할 수 있도록 하였다. 하천수리해석의 기능을 보강하기 위해 역산조도계수 산정모형, 상류-사류 천이류 구간에 대한 부등류 해석모형, 범람 홍수류에 대한 홍수위 산정모형, 하천수리계산시의 불확실도 해석모형 등의 새로운 기능을 추가하여 제시하였다. 모든 입출력자료는 프로젝트 단위별로 운영되어 data의 관리가 손쉽도록 하였으며 결과를 DB에 저장하여 다른 모형에서도 적용할 수 있도록 하였다. 그리고 HyGIS-HMS 및 HyGIS-RAS 모형에서 강우-유출-하도 수리해석-범람해석 등이 일괄되게 하나의 시스템 내에서 구현될 수 있도록 하였다. 따라서 HyGIS와 통합된 수리, 수문모형은 국내 하천 및 유역에 적합한 시스템으로서 향후 HydroInformatics 구현을 염두에 둔 특화된 국내 수자원 분야 소프트웨어의 개발에 기본 토대를 제공할 것으로 판단된다.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.121-128
• /
• 2009

During a hypothetical high-pressure accident in a nuclear power plant (NPP), molten corium can be ejected through a breach of a reactor pressure vessel (RPV) and dispersed by a following jet of a high-pressure steam in the RPV. The dispersed corium is fragmented into smaller droplets in a reactor cavity of the NPP by the steam jet and released into other compartments of the NPP by a overpressure in the cavity. The fragments of the corium transfer thermal energy to the ambient air in the containment or interact chemically with steam and generate hydrogen which may be burnt in the containment. The thermal loads from the ejected molten corium on the containment which is called direct containment heating (DCH) can threaten the integrity of the containment. DCH in a NPP containment is related to many physical phenomena such as multi-phase hydrodynamics, thermodynamics and chemical process. In the evaluation of the DCH load, the melt dispersion rates depending on the RPV pressure are the most important parameter. Mostly, DCH was evaluated by using lumped-analysis codes with some correlations obtained from experiments for the dispersion rates. In this study, MC3D code was used to evaluate the dispersion rates in the APR1400 NPP during the high-pressure accidents. MC3D is a two-phase analysis code based on Eulerian four-fields for melt jet, melt droplets, gas and water. The dispersion rates of the corium melt depending on the RPV pressure were obtained from the MC3D analyses and the values specific to the APR1400 cavity geometry were compared to a currently available correlation.

• Journal of the Korean Geosynthetics Society
• /
• v.14 no.3
• /
• pp.43-51
• /
• 2015

In seasonal frozen areas which have a temperature difference in the winter and spring season like south korea, if stiffness reduction by repeated freezing and thawing occurs to slopes adjacent to private facilities or mountain slopes, safety factor is insufficient to design criteria and landslide could be occurs due to rainfall or snowfall. It can lead to large damage of human life and property. In this study, in order to examine the safety changes of mountain slopes by repeated freezing and thawing, soil samples series of SP and SM by USCS distributed in surface soil of mountain slopes were collected for specimens. Through the direct shear test, the characterestics of frozen soil shear strength were analyzed and by utilizing numerical methods, chracteristics of strength reduction of weathered granite soil according to repeated action of freezing and thawing, changes in the stability of the slopes when applying freezing and thawing of the soil samples were examined. As a result, the maximum shear stress decreased approximately 10%, and slope stability analysis confirmed that required safety factor is less than compare with the non-frozen samples.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.75-82
• /
• 1987

It is the objective of this study to estimate J-integral by numerical analysis, in which J-integral as aparameters in fracture mechanics can be used to evaluate unstable ductile fracture which is a important problem with respect to structural stability when the scope is beyond small scale yielding criteria. For this, 8-node isoparametric singular element as crack tip element of a cracked material was used to solve plastic blunting phenomenon at crack tip, and crack opening was forced to start when J-value was exceeding fracture toughness $J_{IC}$ of the material. And crack propagation behaviour was treated by using crack opening angle. From this study, it was shown that crack opening, stable propagation and unstable opening point of the cracked material found by using J-value obtained from this study were accord with the other study, so think, J-value obtained from this study can be directly used as a parameter in fracture mechanics to deal with the problem of stable propagation of crack and unstable ductile fracture.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.6 no.6
• /
• pp.542-546
• /
• 2005

The J-integral used as a ductile crack initiation criterion has been discussed for the impact loaded elastic-plastic 3PB specimens. The experimental method to measure or estimate the J-integral history has been investigated and its result has been compared to the obtained elastic-plastic values by the finite element model of this study. These numerical results and the experimental curves are found to agree closely. J-integral can be calculated by only numerical analysis with the finite element model. It is proved that simple calculation can be made in order to find the possible value of J-integral by crack mouth opening displacement(CMOD) in the dynamic nonlinear fracture experiment of 3-point bend(3PB) specimen. The property of elastic-plastic material is considered at different impact velocities. The J-integral may be estimated from the crack mouth opening displacement which can be measured directly kom photographs taken during impact experiments.