• Structural Engineering and Mechanics
• /
• v.4 no.1
• /
• pp.25-35
• /
• 1996

This paper describes a continuum variational formulation for design optimization of nonlinear structures in the elastic-plastic domain, where unloading and reloading of the structures are allowed to occur. The Total Lagrangian procedure is used for the description of the structural deformation. The direct differentiation approach is used to derive the sensitivities of the various structural response measures with respect to the design parameters. Since the material goes into the inelastic range and unloading and reloading of the structure are allowed to occur, the structural response is path dependent and an additional step is needed to integrate the constitutive equations. It can be shown, consequently, that design sensitivity analysis is also path-dependent. The theory has been discretized by the finite element technique and implemented in a structural analysis code. Mathematical programming approach is used for the optimization process. Numerical applications on trusses are performed, where cross-sectional areas and nodal point coordinates are treated as design variables. Optimal designs have been obtained and compared by using two different strategies: a two level strategy where the levels are defined accordingly the type of design variables, cross sectional areas or node coordinates, and optimizing simultaneously with respect to both types of design variables.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.22 no.6
• /
• pp.763-772
• /
• 2019

There is a need to use sheath flow nozzle to detect bioaerosol such as virus and bacteria due to their characteristics. In order to enhance the detection performance depending on nozzle parameters, numerical analysis was carried out using a commercial code, ANSYS CFX. Eulerian-lagrangian approach method is used in this simulation. Multiphase flow characteristics between primary fluid and solid were considered. The detection performance was evaluated based on the results of flow field in nozzle chamber such as focusing efficiency and swirl strength. In addition, Latin hypercube sampling(LHS) of design of experiment(DOE) was used for generating a near-random sampling. Then, the acceleration section is optimized using response surface method(RSM). Results show that the optimized model achieved a 6.13 % in a focusing efficiency and 11.47 % increase in swirl strength over the reference model.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2008.05a
• /
• pp.235-237
• /
• 2008

A unified hydrocode, ExLO, in which Largrangian, ALE and Eulerian solvers are incorporated into a single framework, has recently been developed in Korea. It is based on the three dimensional explicit finite element method and written in C++. ExLO is mainly designed for the calculation of structural responses to highly transient loading conditions, such as high-speed impacts, high-speed machining, high speed forming and explosions. In this paper the numerical schemes are described. Some improvements of the material interface and advection scheme are included. Details and issues of the momentum advection scheme are provided. In this paper the modeling capability of ExLO has been described for two extreme loading events; high-speed impacts and explosions. Numerical predictions are in good agreement with the existing experimental data. Specific applications of the code are discussed in a separate paper in this journal. Eventually ExLO will be providing an optimum simulation environment to engineering problems including the fluid-structure interaction problems, since it allows regions of a problem to be modeled with Lagrangian, ALE or Eulerian schemes in a single framework.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.9
• /
• pp.117-126
• /
• 1998

For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.

• Structural Engineering and Mechanics
• /
• v.68 no.6
• /
• pp.721-733
• /
• 2018

The modal frequency responses of functionally graded (FG) sandwich doubly curved shell panels are investigated using a higher-order finite element formulation. The system of equations of the panel structure derived using Hamilton's principle for the evaluation of natural frequencies. The present shell panel model is discretised using the isoparametric Lagrangian element (nine nodes and nine degrees of freedom per node). An in-house MATLAB code is prepared using higher-order kinematics in association with the finite element scheme for the calculation of modal values. The stability of the opted numerical vibration frequency solutions for the various shell geometries i.e., single and doubly curved FG sandwich structure are proven via the convergence test. Further, close conformance of the finite element frequency solutions for the FG sandwich structures is found when compared with the published theoretical predictions (numerical, analytical and 3D elasticity solutions). Subsequently, appropriate numerical examples are solved pertaining to various design factors (curvature ratio, core-face thickness ratio, aspect ratio, support conditions, power-law index and sandwich symmetry type) those have the significant influence on the free vibration modal data of the FG sandwich curved structure.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.749-758
• /
• 2018

This article derived a hybrid coupling technique using the higher-order displacement polynomial and three soft computing techniques (teaching learning-based optimization, particle swarm optimization, and artificial bee colony) to predict the optimal stacking sequence of the layered structure and the corresponding frequency values. The higher-order displacement kinematics is adopted for the mathematical model derivation considering the necessary stress and stain continuity and the elimination of shear correction factor. A nine noded isoparametric Lagrangian element (eighty-one degrees of freedom at each node) is engaged for the discretisation and the desired model equation derived via the classical Hamilton's principle. Subsequently, three soft computing techniques are employed to predict the maximum natural frequency values corresponding to their optimum layer sequences via a suitable home-made computer code. The finite element convergence rate including the optimal solution stability is established through the iterative solutions. Further, the predicted optimal stacking sequence including the accuracy of the frequency values are verified with adequate comparison studies. Lastly, the derived hybrid models are explored further to by solving different numerical examples for the combined structural parameters (length to width ratio, length to thickness ratio and orthotropicity on frequency and layer-sequence) and the implicit behavior discuss in details.

• Structural Engineering and Mechanics
• /
• v.77 no.1
• /
• pp.57-74
• /
• 2021

The thermal eigenvalue responses of the graded sandwich shell structure are evaluated numerically under the variable thermal loadings considering the temperature-dependent properties. The polynomial type rule-based sandwich panel model is derived using higher-order type kinematics considering the shear deformation in the framework of the equivalent single-layer theory. The frequency values are computed through an own home-made computer code (MATLAB environment) prepared using the finite element type higher-order formulation. The sandwich face-sheets and the metal core are discretized via isoparametric quadrilateral Lagrangian element. The model convergence is checked by solving the similar type published numerical examples in the open domain and extended for the comparison of natural frequencies to have the final confirmation of the model accuracy. Also, the influence of each variable structural parameter, i.e. the curvature ratios, core-face thickness ratios, end-support conditions, the power-law indices and sandwich types (symmetrical and unsymmetrical) on the thermal frequencies of FG sandwich curved shell panel model. The solutions are helping to bring out the necessary influence of one or more parameters on the frequencies. The effects of individual and the combined parameters as well as the temperature profiles (uniform, linear and nonlinear) are examined through several numerical examples, which affect the structural strength/stiffness values. The present study may help in designing the future graded structures which are under the influence of the variable temperature loading.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.72.3-73
• /
• 2021

Baryon Acoustic Oscillations (BAO) are caused by acoustic density waves in the early universe and act as a standard ruler in the clustering pattern of galaxies in the late Universe. Measuring the BAO feature in the 2-point correlation function of a sample of galaxies allows us to estimate cosmological distances to the galaxies mean redshift, , which is important for testing and constraining the cosmology model. The BAO feature is also expected to appear in the higher order statistics. In this work we measure the generalized spatial N-point point correlation functions up to 4th order. We made measurements of the 2, 3, and 4-point correlation functions in the SDSS-III DR12 CMASS data, comprising of 777,202 galaxies. The errors and covariances matrices were estimated from 500 mock catalogues. We created a theoretical model for these statistics by measuring the N-point functions in halo catalogues produced by the approximate Lagrangian perturbation theory based simulation code, PINOCCHIO. We created simulations using initial conditions with and without the BAO feature. We find that the BAO is detected to high significance up to the 4-point correlation function.

• Journal of Korean Society on Water Environment
• /
• v.32 no.3
• /
• pp.243-252
• /
• 2016

In cases of water pollution accidents, accurate prediction for arrival time and concentration of contaminants in a river is essential to take proper measures and minimize their impact on downstream water intake facilities. It is critical to fully understand the behavior characteristics of contaminants on river surface, especially in case of oil spill accidents. Therefore, in this study, the effects of main parameters of advection and diffusion of contaminants were analyzed and validated by comparing the results of Lagrangian particle tracking (LPT) simulation of Environmental Fluid Dynamic Code (EFDC) model with those of Global Position System (GPS)-equipped drifter experiment. Prevention scenario modeling was accomplished by taking cases of movable weir operation into account. The simulated water level and flow velocity fluctuations agreed well with observations. There was no significant difference in the speed of surface particle movement between 5 and 10 layer modeling. Therefore, 5 layer modeling could be chosen to reduce computational time. It was found that full three dimensional modeling simulated wind effects on surface particle movements more sensitively than depth-averaged two dimensional modeling. The diffusion range of particles was linearly proportional to horizontal diffusivity by sensitivity analysis. Horizontal diffusivity estimated from the results of GPS-equipped drifter experiment was 0.096 m²/sec, which was considered to be valid for applying the LPT module in this area. Finally, the scenario analysis results showed that particle movements could be stagnant when discharge from the upstream weir was reduced, implying the possibility of securing time for mitigation actions such as oil boom installation and wiping oil contaminants. The outcomes of this study can help improve the prediction accuracy of particle tracking simulation to establish the most suitable mitigation plan considering the combination of movable weir operation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.10
• /
• pp.1793-1799
• /
• 2003

An optimal design of a multi-layered plate structure to endure high-velocity impact has been suggested by using size optimization after numerical simulations. The NET2D, a Lagrangian explicit time-integration finite element code for analyzing high-velocity impact, was used to find the parameters for the optimization. Three different materials such as mild steel, aluminum for a multi-layered plate structure and die steel for the pellet, were assumed. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, Johnson-Cook model and Phenomenological Material Model were used as constitutive models for the simulation. It was carried out with several different gaps and thickness of layers to figure out the trend in terms of those parameters' changes under the constraint, which is against complete penetration. Also, the measuring domain has been shrunk with several elements to reduce the analyzing time. The response surface method based on the design of experiments was used as optimization algorithms. The optimized thickness of each layer in which perforation does not occur has been obtained at a constant velocity and a designated total thickness. The result is quite acceptable satisfying both the minimized deformation energy and the weight criteria. Furthermore, a conceptual idea for topology optimization was suggested for the future work.