• Steel and Composite Structures
• /
• 제47권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
• /
• 제51권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
• /
• 제45권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
• /
• 제57권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.

• 한국기계가공학회지
• /
• 제4권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.

• 공학교육연구
• /
• 제15권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.

• 한국지진공학회논문집
• /
• 제8권1호
• /
• pp.59-65
• /
• 2004

최근 강구조물의 고층화 및 장경간화로 인하여 SM570강재와 같은 고강도 강재의 적용을 필요로 하고 있다. 강구조물의 정확한 내진설계를 위한 내진구조해석시 비선형 반복하중을 받는 강재의 특성을 명확히 포현할 수 있는 구성식이 필요하다. SM570는 최근 그 사용이 증가하고 있으나 아직 반복소성거동의 구현 및 정식화에 관한 연구는 아직 미진하다. 본 연구에서는 인장 및 저싸이클 피로 실험을 통하여 SM570 강재의 반복소성모델을 제안하였다. 제안된 반복소성모델을 3차원 유한요소에 적용하여 SM570이 사용된 원형 강교각의 내진해석을 수행하였다. 실험결과와 내진해석을 통하여 본 연구에서 제안한 구성식은 SM570이 사용된 강구조물의 복잡한 소성거동을 정도 높게 구현함을 알 수 있었다.

• 대한기계학회논문집B
• /
• 제27권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• 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
• /
• 제5권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.