• Journal of Korean Society of Steel Construction
• /
• v.12 no.5 s.48
• /
• pp.569-579
• /
• 2000

An improved multi-level (IML) optimization algorithm using automatic differentiation (AD) of framed structures is proposed in this paper. For the efficiency of the proposed algorithm, multi-level optimization techniques using a decomposition method that separates both system-level and element-level optimizations, that utilizes and an artificial constraint deletion technique, are incorporated in the algorithm. And also to save the numerical efforts, an efficient reanalysis technique through approximated structural responses such as moments and frequencies with respect to intermediate variables is proposed in the paper. Sensitivity analysis of dynamic structural response is executed by AD that is a powerful technique for computing complex or implicit derivatives accurately and efficiently with minimal human effort. The efficiency and robustness of the IML algorithm, compared with a plain multi-level (PML) algorithm, is successfully demonstrated in the numerical examples.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.3
• /
• pp.211-219
• /
• 2008

A three-dimensional inviscid flow solver has been developed based on unstructured meshes for the investigation of the effect of the external stores on the tail surfaces of a generic fighter aircraft. The numerical method is based on a vertex-centered finite-volume discretization and an implicit point Gauss-Seidel time integration. The calculations were made for a steady flow and the computed results were compared with experimental data to validate the flow solver. An unsteady time-accurate computation of the generic fighter aircraft with external stores at transonic flight conditions showed that the external stores cause unsteady loading on the horizontal tail surface due to the mutual interference between their wake and the horizontal tail surface. It was shown that downward deflection of the trailing edge flap significantly reduces the undesirable interference effect.

• Tunnel and Underground Space
• /
• v.12 no.4
• /
• pp.304-311
• /
• 2002

This study reports the influence of vocid geometry on fracture closure and permeability from numerical experiments. As the aperture distributions of rock fractures are characterized by statistical methods, synthetic fractures have successfully been simulated in this way. Based on the generated fracture models, models for fracture closure and flow calculation have been developed. A fracture closure model has been developed by considering the asperity compression and half-space deformation, and flow calculations have been performed using a finite difference method adopting a local cubic law. The results of numerical experiments have shown that the increase in the aperture spatial correlation leads the fracture closure and the decrease in fracture permeability to increase. Also, it has been indicated that there is an implicit relation between fracture normal stiffness and permeability. The importance of this study is to enhance the understanding the hydro-mechanical behavior of fractured rock massed due to engineering projects.

• Computers and Concrete
• /
• v.25 no.1
• /
• pp.67-73
• /
• 2020

Traditionally used analytical approach to predict the fatigue failure of reinforced concrete (RC) structure is generally conservative and has certain limitations. The nonlinear finite element method (FEM) offers less expensive solution for fatigue analysis with sufficient accuracy. However, the conventional implicit dynamic analysis is very expensive for high level computation. Whereas, an explicit dynamic analysis approach offers a computationally operative modelling to predict true responses of a structural element under periodic loading and might be perfectly matched to accomplish long life fatigue computations. Hence, this study simulates the fatigue behaviour of RC beams with finite element (FE) assemblage presenting a simplified explicit dynamic numerical solution to show computer aided fatigue behaviour of RC beam. A commercial FEM package, ABAQUS has been chosen for this complex modelling. The concrete has been modelled as a 8-node solid element providing competent compression hardening and tension stiffening. The steel reinforcements are simulated as two-node truss elements comprising elasto-plastic stress-strain behaviour. All the possible nonlinearities are duly incorporated. Time domain analysis has been adopted through an automatic Newmark-β time incremental technique. The program consists of twelve RC beams to visualize the real behaviour during fatigue process and to obtain the reliability of the study. Both the numerical and experimental results indicate a redistribution of stresses along the time and damage accumulation of beam which severely affect the serviceability and ultimate capacity of RC beam. The output of the FEM analysis demonstrates good match with the experimental consequences which affirm the efficacy of the computer aided model. The controlled fatigue damage evolution at service fatigue load limits makes the FE model an efficient tool in predicting high cycle fatigue behaviour of RC structures.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.1
• /
• pp.35-45
• /
• 2005

Numerical study is carried out to investigate the effects of chemistry and thermal radiation on the rocket plume flow field at various altitudes. Navier-Stokes equations for compressible flows were solved by a fully-implicit TVD code based on the finite volume method. An infinitely fast chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thick media were incorporated with the fluid dynamics code. The plume flow fields of a kerosene-fueled rocket flying at Mach number zero at sea-level, 1.16 at altitude of 5.06 km and 2.90 at 17.34 km were numerically analyzed. Results showed the plume structures at different altitude conditions with the effects of chemistry and radiation. It is understood that the excess temperature by the chemical reactions in the exhaust gas may not be ignored in the view point of propulsion performance and thermal protection of the rocket base, especially at higher altitude conditions.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.19 no.1
• /
• pp.29-36
• /
• 2015

Base isolation is considered as a seismic protective system in the design of next generation Nuclear Power Plants (NPPs). If seismic isolation devices are installed in nuclear power plants then the safety under a seismic load of the power plant may be improved. However, with respect to some equipment, seismic risk may increase because displacement may become greater than before the installation of a seismic isolation device. Therefore, it is estimated to be necessary to select equipment in which the seismic risk increases due to an increase in the displacement by the installation of a seismic isolation device, and to perform research on the seismic performance of each piece of equipment. In this study, modified NRC-BNL benchmark models were used for seismic analysis. The numerical models include representations of isolation devices. In order to validate the numerical piping system model and to define the failure mode, a quasi-static loading test was conducted on the piping components before the analysis procedures. The fragility analysis was performed by using the results of the inelastic seismic response analysis. Inelastic seismic response analysis was carried out by using the shell finite element model of a piping system considering internal pressure. The implicit method was used for the direct integration time history analysis. In addition, the collapse load point was used for the failure mode for the fragility analysis.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.88-99
• /
• 1996

The explicit scheme for finite element analysis of sheet metal forming problems has been widely used for providing practical solutions since it improves the convergency problem, memory size and computational time especially for the case of complicated geometry and large element number. The explicit schemes in general use are based on the elastic-plastic modeling of material requiring large computataion time. In the present work, a basic formulation for rigid-plastic explicit finite element analysis of plain strain sheet metal forming problems has been proposed. The effect of some basic parameters involved in the dynamic analysis has been studied in detail. Thus, the effective ranges of parameters have been proposed for numerical simultion by the rigid-plastic explicit finite element method. A direct trial-and-error method is introduced to treat contact and friction. In computation, sheet material is assumed to possess normal anisotropy and rigid-plastic workhardening characteristics. In order to show the validity and effectiveness of the proposed explicit scheme, computations are carried out for cylindrical punch stretching and the computational results are compared with those by the implicit scheme as well as with a commercial code. The proposed rigid-plastic exlicit finite element method can be used as a robust and efficient computational method for analysis of sheet metal forming.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.162-176
• /
• 2022

In Lead-based reactor (LBR) severe accident, the meltdown and migration inside the reactor core will lead to fuel fragment concentration, which may further cause re-criticality and even core disintegration. Accurately predicting the migration and solidification behavior of melt in LBR severe accidents is of prime importance for safety analysis of LBR. In this study, the Moving Particle Semi-implicit (MPS) method is validated and used to simulate the migration and solidification behavior. Two main surface tension models are validated and compared. Meanwhile, the MPS method is validated by the L-plate solidification test. Based on the improved MPS method, the migration and solidification behavior of melt in LBR severe accident was studied furthermore. In the Pb-Bi coolant, the melt flows upward due to density difference. The migration and solidification behavior are greatly affected by the surface tension and viscous resistance varying with enthalpy. The whole movement process can be divided into three stages depending on the change in velocity. The heat transfer of core melt is determined jointly by two heat transfer modes: flow heat transfer and solid conductivity. Generally, the research results indicate that the MPS method has unique advantage in studying the migration and solidification behavior in LBR severe accident.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.567-578
• /
• 2002

The main purpose of this study is to investigate the dynamic behaviour of containment building in nuclear power plant excited by aircraft impact loading using a lower order 8-node solid element. The yield and failure surfaces for concrete material model is formulated on the basis of Drucker-Prager yield criteria and are assumed to be varied by taking account of the visco-plastic energy dissipation. The standard 8-node solid element has prone to exhibit the element deficiencies and the so-called B bar method proposed by Hughes is therefore adopted in this study. The implicit Newmark method is adopted to ensure the numerical stability during the analysis. Finally, the effect of different levels of cracking strain and several types of aircraft loading are examined on the dynamic behaviour of containment building and the results are quantitatively summarized as a future benchmark.