• Steel and Composite Structures
• /
• v.47 no.5
• /
• pp.569-585
• /
• 2023

This paper proposes an efficient approach for the structural topology optimization of bi-directional functionally graded structures by incorporating popular radial basis functions (RBFs) into an implicit level set (ILS) method. Compared to traditional element density-based methods, a level set (LS) description of material boundaries produces a smoother boundary description of the design. The paper develops RBF implicit modeling with multiquadric (MQ) splines, thin-plate spline (TPS), exponential spline (ES), and Gaussians (GS) to define the ILS function with high accuracy and smoothness. The optimization problem is formulated by considering RBF-based nodal densities as design variables and minimizing the compliance objective function. A LS-RBF optimization method is proposed to transform a Hamilton-Jacobi partial differential equation (PDE) into a system of coupled non-linear ordinary differential equations (ODEs) over the entire design domain using a collocation formulation of the method of lines design variables. The paper presents detailed mathematical expressions for BiDFG beams topology optimization with two different material models: continuum functionally graded (CFG) and mechanical functionally graded (MFG). Several numerical examples are presented to verify the method's efficiency, reliability, and success in accuracy, convergence speed, and insensitivity to initial designs in the topology optimization of two-dimensional (2D) structures. Overall, the paper presents a novel and efficient approach to topology optimization that can handle bi-directional functionally graded structures with complex geometries.

• Nuclear Engineering and Technology
• /
• v.51 no.6
• /
• pp.1504-1513
• /
• 2019

Recently, ocean nuclear reactors have received attention due to enhanced safety features. The movable and transportable characteristics distinguish ocean nuclear reactors from land-based nuclear reactors. Therefore, for safety/design analysis of the ocean reactor, the thermos-hydraulics must be investigated in the moving system. However, there are no studies reporting the general two-fluid equations that can be used for multi-dimensional simulations of two-phase flows in moving systems. This study is to systematically formulate the multi-dimensional two-fluid equations in the non-inertial frame of reference. To demonstrate the applicability of the formulated equations, we perform a total of six different simulations in 2D tanks with translational and/or rotational motions.

• Nuclear Engineering and Technology
• /
• v.45 no.4
• /
• pp.439-458
• /
• 2013

REX-10 is a fully-passive small modular reactor in which the coolant flow is driven by natural circulation, the RCS is pressurized by a steam-gas pressurizer, and the decay heat is removed by the PRHRS. To confirm design decisions and analyze the transient responses of an integral PWR such as REX-10, a thermal-hydraulic system code named TAPINS (Thermal-hydraulic Analysis Program for INtegral reactor System) is developed in this study. Based on a one-dimensional four-equation drift-flux model, TAPINS incorporates mathematical models for the core, the helical-coil steam generator, and the steam-gas pressurizer. The system of difference equations derived from the semi-implicit finite-difference scheme is numerically solved by the Newton Block Gauss Seidel (NBGS) method. TAPINS is characterized by applicability to transients with non-equilibrium effects, better prediction of the transient behavior of a pressurizer containing non-condensable gas, and code assessment by using the experimental data from the autonomous integral effect tests in the RTF (REX-10 Test Facility). Details on the hydrodynamic models as well as a part of validation results that reveal the features of TAPINS are presented in this paper.

• Structural Engineering and Mechanics
• /
• v.57 no.6
• /
• pp.989-1014
• /
• 2016

This paper deals with the analysis of the natural frequencies, mode shapes of an axially loaded beam system carrying ends consisting of non-concentrated tip masses and three spring-two mass sub-systems. The influence of system design and sub-system parameters on the combined system characteristics is the major part of this investigation. The effect of material properties, rotary inertia and shear deformation of the beam system is included. The end masses are elastically supported against rotation and translation at an offset point from the point of attachment. Sub-systems are attached to center of gravity eccentric points out of the beam span. The boundary conditions of the ordinary differential equation governing the lateral deflections and slope due to bending of the beam system including developed shear force frequency dependent terms, due to the sub.system suspension, have been formulated. Exact formulae for the modal frequencies and the modal shapes have been derived. Based on these formulae, detailed parametric studies are carried out. The geometrical and mechanical parameters of the system under study have been presented in non-dimensional analysis. The applied mathematical model is presented to cover wide range of mechanical, naval and structural engineering applications.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.4 no.2
• /
• pp.50-56
• /
• 2005

To demonstrate the feasibility of utilizing FCAD chart proposed in our previous work, series of crack growth/arrest behavior in the integrally stiffened panels were simulated by numerical method using upper mentioned FCAD charts and a new crack growth rate equation. It is concluded that proposed family of FCAD curves, in the form of non-dimensional arrest load ranges, are reliable indicators of fatigue crack growth/arrest behavior of integrally stiffened panels considered here.

• Journal of Engineering Education Research
• /
• v.15 no.5
• /
• pp.93-97
• /
• 2012

Most of ubiquitous computing devices such as stereo camera, ultrasonic sensor based MIT cricket system and other wireless sensor network devices are widely applied to the 2 Dimensional(2D) localization system in today. Because stereo camera cannot estimate the optimal location between moving node and beacon node in Wireless Personal Area Network(WPAN) under Non Line Of Sight(NLOS) environment, it is a great weakness point to the design of the 2D localization system in indoor environment. But the conventional 2D triangulation scheme that is adapted to the MIT cricket system cannot estimate the 3 Dimensional(3D) coordinate values for estimation of the optimal location of the moving node generally. Therefore, the 3D triangulation scheme based on the space segmentation in WPAN is suggested in this paper. The measuring data in the suggested scheme by computer simulation is compared with that of the geographic measuring data in the AutoCAD software system. The average error of coordinates values(x,y,z) of the moving node is calculated to 0.008m by the suggested scheme. From the results, it can be seen that the location correctness of the suggested scheme is very excellent for using the localization system in WPAN.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.8 no.1
• /
• pp.59-65
• /
• 2004

Recently, as steel structures become higher and more long-spanned, application of high strength steel. SM570, is considered, For accurate seismic design, seismic analysis of steel structures needs a constitutive equation describing the characteristic of steel materials under non-proportional cyclic loading, While the use of SM570 material is much increased these days, research for description and generalization of cyclic plasticity behavior are insufficient, In this study, a cyclic plasticity model is proposed by results of material tests, i.e, monotonic and low cycle tests, Proposed cyclic plasticity model is applied to 3-Dimensional FE program and we carried out seismic analysis of pipe-section steel pier using SM570, Comparison between experiment and analysis results shows that the proposed constitutive equation is able to describe exactly the complicated plastic behavior of steel structure using SM570.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1457-1463
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a multi-blade centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k - c turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in this centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2157-2161
• /
• 2003

In this paper, the response surface method using three-dimensional Navier-Stokes analysis to optimize the shape of a forward-curved blades centrifugal fan, is described. For numerical analysis, Reynolds-averaged Navier-Stokes equations with standard k-e turbulence model are transformed into non-orthogonal curvilinear coordinate system, and are discretized with finite volume approximations. Due to the large number of blades in forward-curved blades centrifugal fan, the flow inside of the fan is regarded as steady flow by introducing the impeller force models for economic calculations. Linear Upwind Differencing Scheme(LUDS) is used to approximate the convection terms in the governing equations. SIMPLEC algorithm is used as a velocity-pressure correction procedure. Design variables, location of cur off, radius of cut off, expansion angle of scroll and width of impeller were selected to optimize the shapes of scroll and blades. Data points for response evaluations were selected by D-optimal design, and linear programming method was used for the optimization on the response surface. As a main result of the optimization, the efficiency was successfully improved. It was found that the optimization process provides reliable design of this kind of fans with reasonable computing time

• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.