• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.2
• /
• pp.135-143
• /
• 2017

A customized CFD simulator to perform thermo-fluid dynamic simulations of an HVAC for an indoor space is presented. This simulation system has been developed for engineers studying architectural engineering, as the HVAC mechanical systems used in housings and buildings. Hence, all functions and options are so designed to be suitable that they are suitable for non-CFD experts as well as CFD engineers. A Computational mesh is generated by open-source libraries, FEMM (Finite Element Method Magnetics), and OpenFOAM. Once the boundary conditions are set, the fluid dynamic calculations are performed using the OpenFOAM solver. Numerical results are validated by comparing them with the experimental data for a simple indoor air flow case. In this paper, an entirely new calculation process is introduced, and the flow simulation results for a sample office room are also discussed.

• Steel and Composite Structures
• /
• v.23 no.4
• /
• pp.409-420
• /
• 2017

This paper theoretically studies the cyclic behavior of hybrid connections between steel coupling beams and concrete shear walls with embedded steel columns. Finite element models of connections with long and short embedded steel columns are built in ABAQUS and validated against the test results in the companion paper. Parametric studies are carried out using the validated FE model to determine the key influencing factors on the load-bearing capacity of connections. A close-form solution of the load-bearing capacity of connections is proposed by considering the contributions from the compressive strength of concrete at the interface between the embedded beam and concrete, shear yielding of column web in the tensile region, and shear capacity of column web and concrete in joint zone. The results show that the bond slip between embedded steel members and concrete should be considered which can be simulated by defining contact boundary conditions. It is found that the loadbearing capacity of connections strongly depends on the section height, flange width and web thickness of the embedded column. The accuracy of the proposed calculation method is validated against test results and also verified against FE results (with differences within 10%). It is recommended that embedded steel columns should be placed along the entire height of shear walls to facilitate construction and enhance the ductility. The thickness and section height of embedded columns should be increased to enhance the load-bearing capacity of connections. The stirrups in the joint zone should be strengthened and embedded columns with very small section height should be avoided.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.223-231
• /
• 2019

A continuum-based design sensitivity analysis(DSA) method is developed for electrostatic problems. To consider high order objective functions, we use 9-node finite element basis functions for analysis and DSA methods. As the design variables are parameterized with B-spline functions, smooth boundary variations are naturally obtained. To solve mesh entanglement problems during the optimization process, a mesh regularization scheme is employed. By minimizing the Dirichlet energy functional, mesh uniformity can be automatically achieved. In numerical examples for maximizing dielectrophoresis forces, the numerical results are compared with well-known electrode geometries and the obtained characteristics are discussed.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.7
• /
• pp.265-274
• /
• 2018

The increasing use of curved beams in buildings, vehicles, ships, and aircraft has prompted studies directed toward the development of an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of a large number of investigations. Solutions of the relevant differential equations have been obtained traditionally using standard finite difference or finite element methods. These techniques require a great deal of computer time as the number of discrete nodes becomes relatively large under the conditions of complex geometry and loading. One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. The differential quadrature method (DQM) has been applied to a large number of cases to overcome the difficulties of the complex algorithms of programming for the computer, as well as the excessive use of storage due to the conditions of complex geometry and loading. The in-plane buckling of curved beams considering the extensibility of the arch axis was analyzed under uniformly distributed radial loads using the DQM. The critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with the precise results by other methods for cases, in which they were available. The DQM, using only a limited number of grid points, provided results that agreed very well (less than 0.3%) with the exact ones. New results according to diverse variations were obtained, showing the important roles in the buckling behavior of curved beams, and can be used in comparisons with other numerical solutions or with experimental test data.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.10
• /
• pp.300-309
• /
• 2016

Curved beams are increasingly used in buildings, vehicles, ships, and aircraft, which has resulted in considerable effort towards developing an accurate method for analyzing the dynamic behavior of such structures. The stability behavior of elastic curved beams has been the subject of many investigations. Solutions to the relevant differential equations have traditionally been obtained by the standard finite difference or finite element methods. However, these techniques require a great deal of computer time for a large number of discrete nodes with conditions of complex geometry and loading. One efficient procedure for the solution of partial differential equations is the differential quadrature method (DQM). This method has been applied to many cases to overcome the difficulties of complex algorithms and high storage requirements for complex geometry and loading conditions. Out-of-plane buckling of curved beams with rotatory inertia were analyzed using DQM under uniformly distributed radial loads. Critical loads were calculated for the member with various parameter ratios, boundary conditions, and opening angles. The results were compared with exact results from other methods for available cases. The DQM used only a limited number of grid points and shows very good agreement with the exact results (less than 0.3% error). New results according to diverse variation are also suggested, which show important roles in the buckling behavior of curved beams and can be used for comparisons with other numerical solutions or experimental test data.

• Wind and Structures
• /
• v.10 no.4
• /
• pp.367-382
• /
• 2007

A nonlinear numerical method was developed to assess the stability of suspension bridge catwalks under a wind load. A section model wind tunnel test was used to obtain a catwalk's aerostatic coefficients, from which the displacement-dependent wind loads were subsequently derived. The stability of a suspension bridge catwalk was analyzed on the basis of the geometric nonlinear behavior of the structure. In addition, a full model test was conducted on the catwalk, which spanned 960 m. A comparison of the displacement values between the test and the numerical simulation shows that a numerical method based on a section model test can be used to effectively and accurately evaluate the stability of a catwalk. A case study features the stability of the catwalk of the Runyang Yangtze suspension bridge, the main span of which is 1490 m. Wind can generally attack the structure from any direction. Whenever the wind comes at a yaw angle, there are six wind load components that act on the catwalk. If the yaw angle is equal to zero, the wind is normal to the catwalk (called normal wind) and the six load components are reduced to three components. Three aerostatic coefficients of the catwalk can be obtained through a section model test with traditional test equipment. However, six aerostatic coefficients of the catwalk must be acquired with the aid of special section model test equipment. A nonlinear numerical method was used study the stability of a catwalk under a yaw wind, while taking into account the six components of the displacement-dependent wind load and the geometric nonlinearity of the catwalk. The results show that when wind attacks with a slight yaw angle, the critical velocity that induces static instability of the catwalk may be lower than the critical velocity of normal wind. However, as the yaw angle of the wind becomes larger, the critical velocity increases. In the atmospheric boundary layer, the wind is turbulent and the velocity history is a random time history. The effects of turbulent wind on the stability of a catwalk are also assessed. The wind velocity fields are regarded as stationary Gaussian stochastic processes, which can be simulated by a spectral representation method. A nonlinear finite-element model set forepart and the Newmark integration method was used to calculate the wind-induced buffeting responses. The results confirm that the turbulent character of wind has little influence on the stability of the catwalk.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.33 no.5
• /
• pp.339-349
• /
• 2020

Double-leaf blast-resistant doors consisting of steel box and slab are application-specific structures installed at the entrances of protective facilities. In these structural systems, certain spacing is provided between the door and wall. However, variation in the boundary condition and structural behavior due to this spacing are not properly considered in the explosion analysis and design. In this study, the structural response and failure behavior based on two variables such as the spacing and blast pressure were analyzed using the finite element method. The results revealed that the two variables affected the overall structural behavior such as the maximum and permanent deflections. The degree of contact due to collision between the door and wall and the impact force applied to the door varied according to the spacing. Hence, the shear-failure behavior of the concrete slab was affected by this impact force. Doors with spacing of less than 10 mm were vulnerable to shear failure, and the case of approximately 15-mm spacing was more reasonable for increasing the flexural performance. For further study, tests and numerical research on the structural behavior are needed by considering other variables such as specifications of the structural members and details of the slab shear design.

• The Journal of Engineering Geology
• /
• v.7 no.2
• /
• pp.139-149
• /
• 1997

Electrical resistivity monitoring methods were adopted to detect groundwater table change in alluvium. Numerical modelling test using finite element method(FEM) and field resisfivity monitoring were conducted in the study. The field monitoring data were acquired in the alluvium deposit site in Jeong-Dong Ri, Geum River where pumping test had been conducted continuously for 20 days to make artificial changes of groundwater table. The unit distance of the electrode array was 4m and 21 fixed electrodes were applied in numerical calculation and field data acquisition. "Modified Wenner" and dipole-dipole array configurations were used in the study. The models used in two-dimensional numerical test were designed on the basis of the simplifving geological model of the alluvium in Jeong Dong Ri, Geum River. Numerical test results show that the apparent resistivity pseudosections were changed in the vicinity of the pootion where groundwater table was changed. Furthermore, there are some apparent resistivity changes in the boundary between aquifer and crystalline basement rock which overlays the aquifer. The field monitoring data also give similar results which were observed in numerical tests. From the numerical test using FEM and field resistivity monitoring observations in alluvium site of Geum River, the electrical monitoring method is proved to be a useful tool for detecting groundwater behavior including groundwater table change. There are some limitations, however, in the application of the resistivity method only because the change of groundwater table does not give enough variations in the apparent resistivity pseudosections to estimate the amount of groundwater table change. For the improved detection of groundwater table changes, it is desirable to combine the resistivity method with other geophysical methods that reveal the underground image such as high-resolution seismic and/or ground penetrating radar surveys.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.19 no.3
• /
• pp.194-202
• /
• 2016

This study aims to review a topology optimization based on finite element analysis (FEA) for conceptual design of platform in the 10MW class floating type wave-wind hybrid power generation system (WHPGS). Two topology optimization theories, density method (DM) and homogenization design method (HDM) were used to check which one is more effective for a simplified structural design problem prior to the topology optimization of platform of WHPGS. From the results of the simplified design problem, the HDM was applied to the topology optimization of platform of WHPGS. For the conceptual platform design of WHPGS, FEA model was created and then the structural analysis was performed considering offshore environmental loads at installation site. Hydrodynamics analysis was carried out to calculate pressure on platform and tension forces in mooring lines induced from the offshore environmental loads such as design wave and current. Loading conditions for the structural analysis included the analysis results from the hydrodynamic analysis and the weights of WHPGS. Boundary condition was realized using inertia relief method. The topology optimization of WHPGS platform was performed using the HDM, and then the conceptual arrangement of main structural members was suggested. From the results, it was confirmed that the topology optimization might be a useful tool to design the conceptual arrangement of main structural members for a newly developed offshore structure such as the floating type WHPGS.

• The Journal of Korean Academy of Prosthodontics
• /
• v.49 no.2
• /
• pp.168-176
• /
• 2011

Purpose: This study aims at investigating the influence of various insertion torques on thermal changes of bone. A proper insertion torque is derived based on the thermal analysis with two different implant designs. Materials and methods: For implant materials, bovine scapula bone of 15 - 20 mm thickness was cut into 35 mm by 40 - 50 mm pieces. Of these, the pieces having 2 - 3 mm thickness cortical bone were used as samples. Then, the half of the sample was immersed in a bath of $36.5^{\circ}C$ and the other half was exposed to ambient temperature of $25^{\circ}C$, so that the inner and surface temperatures reached $36.5^{\circ}C$ and $28^{\circ}C$, respectively. Two types of implants ($4.5{\times}10\;mm$ Br${\aa}$nemark type, $4.8{\times}10\;mm$ Microthread type) were inserted into bovine scapula bone and the temperature was measured by a thermocouple at 0.2 mm from the measuring point. Finite element method (FEM) was used to analyze the thermal changes at contacting surface assuming that the sample is a cube of $4\;cm{\times}4\;cm{\times}2\;cm$ and a layer up to 2 mm from the top is cortical bone and below is a cancellous bone. Boundary conditions were set on the basis of the shape of cavity after implants. SolidWorks was used as a CAD program with the help of Abaqus 6.9-1. Results: In the in-vitro experiment, the Microhead type implant gives a higher maximum temperature than that of the Br${\aa}$nemark type, which is attributed to high frictional heat that is associated with the implant shape. In both types, an Eriksson threshold was observed at torques of 50 Ncm (Br${\aa}$nemark) and 35 Ncm (Microthread type), respectively. Based on these findings, the Microthread type implant is more affected by insertion torques. Conclusion: This study demonstrate that a proper choice of insertion torque is important when using a specific type of implant. In particular, for the Microthread type implant, possible bone damage may be expected as a result of frictional heat, which compensates for initial high success rate of fixation. Therefore, the insertion torque should be adjusted for each implant design. Furthermore, the operation skills should be carefully chosen for each implant type and insertion torque.