• Structural Engineering and Mechanics
• /
• v.30 no.5
• /
• pp.619-645
• /
• 2008

This study presents the efficiency of simulating structural systems using a method that combines a simplified component model (SCM) and rigorous component model (RCM). To achieve a realistic simulation of structural systems, a numerical model must be adequately capturing the detailed behaviors of real systems at various scales. However, capturing all details represented within an entire structural system by very fine meshes is practically impossible due to technological limitations on computational engineering. Therefore, this research develops an approach to simulate large-scale structural systems that combines a simplified global model with multiple detailed component models adjusted to various scales. Each correlated multi-scale simulation model is linked to others using a multi-level hierarchical modeling simulation method. Simulations are performed using nonlinear finite element analysis. The proposed method is applied in an analysis of a simple reinforced concrete structure and the Reuipu Elementary School (an existing structure), with analysis results then compared to actual onsite observations. The proposed method obtained results very close to onsite observations, indicating the efficiency of the proposed model in simulating structural system behavior.

• Journal of Korea Water Resources Association
• /
• v.30 no.1
• /
• pp.63-73
• /
• 1997

Finite element model is established for the simulation of groundwater flow due to water curtain around underground oil storage Choleski decomposition method. The symmetric global conductance matrix is solved by vector storage Choleski decomposition method. The model is verified through comparison with the results of electric analogy. For the application of this model to real site, the finite element meshes are constructed according to representative vertical cross and longitudinal sections. In cross-sectional analysis, potential and flow distributions are compared based on the cavern pressure and horizontal water curtain. For longitudinal section, effects between nearly located caverns with or without vertical water curtain are analyzed. These results prove that the established model can be used as a tool for flow analysis around underground caverns.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.99-110
• /
• 2005

This work presents results of finite element analysis of isothermal incompressible creeping viscoelastic flows with the tensor-logarithmic formulation of the Leonov model especially for the planar geometry with singular comers in the domain. In the case of 4:1 contraction flow, for all 5 meshes we have obtained solutions over the Deborah number of 100, even though there exists slight decrease of convergence limit as the mesh becomes finer. From this analysis, singular behavior of the comer vortex has been clearly seen and proper interpolation of variables in terms of the logarithmic transformation is demonstrated. Solutions of 4:1:4 contraction/expansion flow are also presented, where there exists 2 singular comers. 5 different types spatial resolutions are also employed, in which convergent solutions are obtained over the Deborah number of 10. Although the convergence limit is rather low in comparison with the result of the contraction flow, the results presented herein seem to be the only numerical outcome available for this flow type. As the flow rate increases, the upstream vortex increases, but the downstream vortex decreases in their size. In addition, peculiar deflection of the streamlines near the exit comer has been found. When the spatial resolution is fine enough and the Deborah number is high, small lip vortex just before the exit comer has been observed. It seems to occur due to abrupt expansion of the elastic liquid through the constriction exit that accompanies sudden relaxation of elastic deformation.

• Korea-Australia Rheology Journal
• /
• v.13 no.1
• /
• pp.37-45
• /
• 2001

Numerical simulation of coextrusion process of viscoelastic fluids within a die has been carried out. In the coextrusion process velocity profile at the outflow boundary is not known a priori, which makes it difficult to impose the proper boundary condition at the outflow boundary. This difficulty has been avoided by using the open boundary condition (OBC) method. In this study, elastic viscous stress splitting (EVSS) formulation with streamline upwind (SU) method has been used in the finite element method. In order to test the validity of the OBC method, comparison between the results of fully developed condition at the outlet and those of OBC has been made for a Newtonian fluid. In the case of upper convected Maxwell (UCM) fluid, the effect of outflow boundary condition on the interface position has been investigated by using two meshes having different downstream lengths. In both cases, the results with the OBC method showed reasonable interface shape. In particular, for the UCM fluid the interface shape calculated with OBC was independent of the downstream length, while the results with the zero traction condition showed oscillation of interface position close to the outlet. Viscosity difference was found to be more important than elasticity difference in determining the final interface position. However, the overshoot of interface position near the con-fluent point increased with elasticity.

• Journal of the Korean Society of Safety
• /
• v.21 no.2 s.74
• /
• pp.80-88
• /
• 2006

In the present study damage behavior of resilient sleepers on concrete ballasts is analyzed. Cracks of resilient sleepers in a railway track system are concentrated on inside of blocks to which the tie bars are connected. Finite element analysis is performed by dividing a block into the straight section and the curved section according to the load condition of the resilient sleeper, and limited the interpretation within the range of resilience. In addition, the value of stress obtained from the interpretation was compared with the allowable stress of concrete to determine the safety. According to the result of numerical analysis, compared with the stress before unequal settlement, the tensile stress of the inside of the block increased significantly after the settlement considering the entire block, and the tensile stress of this part exceeded the allowable stress of concrete, so was undesirable in terms of safety. In reality, the arrangement of tensile stiffeners inside blocks connected to tie bars is improper in the design of resilient sleepers, and when unequal settlement occurs, tensile stress increases on this part and consequently causes cracking damage. It is necessary to arrange wire meshes or tensile reinforcing bars in a structurally safe way to reinforce the inside of blocks on which cracks are concentrated.

• Proceedings of the KIEE Conference
• /
• 2006.07c
• /
• pp.1494-1495
• /
• 2006

A complete finite element analysis method for discharge onset process, which is governed and coupled by charge transport equation and electric field equation, was presented. The charge transport equation of first order was transformed into a second-order one by utilizing the artificial diffusion scheme. The two second-order equations were analyzed by the finite element formulation which is well-developed for second-order ones. The Fowler-Nordheim injection boundary condition was adopted for charge transport equation. After verifying the numerical results by comparing to the analytic solutions using parallel plane electrodes with one carrier system, we extended the result to blade-plane electrodes in 2D xy geometry with three carriers system. Radius of the sharp tip was taken to be 50 ${\mu}m$. When this sharp geometry was solved by utilizing the space discretizing methods, the very sharp tip was found to cause a singularity in electric field and space charge distribution around the tip. To avoid these numerical difficulties in the FEM, finer meshes, a higher order shape function, and artificial diffusion scheme were employed.

• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.5 s.170
• /
• pp.205-213
• /
• 2005

This paper addresses a method for numerical simulation in the pressing process of hot glass. Updated Lagrangian finite element formulations are employed for the flow and energy equations to accommodate moving meshes. The model is assumed axi-symmetric and creep flow is assumed due to the high viscosity. Commercial software ANSYS is used to solve the coupled flow and energy equations. Moving contact points as well as free surface during the pressing are effectively calculated and updated by utilizing API functions of CAD software Unigraphics. The mesh distortion problem near the wall is overcome by automatic remeshing, and the temperatures of the new mesh are conveniently interpolated by using a unique function of ANSYS. The developed model is applied to the pressing process of TV glasses. In conclusion, the presented method shows that the pressing process accompanying moving boundary can be simulated by effectively combining general purpose software without resorting to special dedicated codes.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.6
• /
• pp.1004-1014
• /
• 2001

The fully automatic algorithm from initial finite element mesh generation to remeshing in two dimensional geometry is introduced using bubble packing method (BPM) for finite element analysis. BPM determines the node placement by force-balancing configuration of bubbles and the triangular meshes are made by Delaunay triangulation with advancing front concept. In BPM, we suggest two node-search algorithms and the adaptive/recursive bubble controls to search the optimal nodal position. To use the automatically generated mesh information in FEA, the new enhanced bandwidth minimization scheme with high efficiency in CPU time is developed. In the remeshing stage, the mesh refinement is incorporated by the control of bubble size using two parameters. And Superconvergent Patch Recovery (SPR) technique is used for error estimation. To verify the capability of this algorithm, we consider two elasticity problems, one is the bending problem of short cantilever beam and the tension problem of infinite plate with hole. The numerical results indicate that the algorithm by BPM is able to refine the mesh based on a posteriori error and control the mesh size easily by two parameters.

• Nuclear Engineering and Technology
• /
• v.41 no.9
• /
• pp.1171-1180
• /
• 2009

The internal core baffle structure of a French Pressurized Water Reactor (PWR) consists of a collection of baffles and formers that are attached to the barrel. The connections are done thanks to a large number of bolts (about 1500). After inspection, some of the bolts have been found cracked. This has been attributed to the Irradiation Assisted Stress Corrosion Cracking (IASCC). The $Electricit\acute{e}$ De France (EDF) has set up a research program to gain better knowledge of the temperature distribution, which may affect the bolts and the whole structure. The temperature distribution in the structure was calculated thanks to the thermal code SYRTHES that used a finite element approach. The heat transfer between the by-pass flow inside the cavities of the core baffle and the structure was accounted for thanks to a strong thermal coupling between the thermal code SYRTHES and the CFD code named Code_Saturne. The results for the CP0 plant design show that both the high temperature and strong temperature gradients could potentially induce mechanical stresses. The CPY design, where each bolt is individually cooled, had led to a reduction of temperatures inside the structures. A new parallel version of SYRTHES, for calculations on very large meshes and based on MPI, has been developed. A demonstration test on the complete structure that has led to about 1.1 billion linear tetraedra has been calculated on 2048 processors of the EDF Blue Gene computer.

• Steel and Composite Structures
• /
• v.51 no.4
• /
• pp.363-375
• /
• 2024

Within the optimization field, addressing the intricate posed by fluidic pressure loads on functionally graded structures with frequency-related designs is a kind of complex design challenges. This paper thus introduces an innovative density-based topology optimization strategy for frequency-constraint functionally graded structures incorporating Darcy's law and a drainage term. It ensures consistent treatment of design-dependent fluidic pressure loads to frequency-related structures that dynamically adjust their direction and location throughout the design evolution. The porosity of each finite element, coupled with its drainage term, is intricately linked to its density variable through a Heaviside function, ensuring a seamless transition between solid and void phases. A design-specific pressure field is established by employing Darcy's law, and the associated partial differential equation is solved using finite element analysis. Subsequently, this pressure field is utilized to ascertain consistent nodal loads, enabling an efficient evaluation of load sensitivities through the adjoint-variable method. Moreover, this novel approach incorporates load-dependent structures, frequency constraints, functionally graded material models, and polygonal meshes, expanding its applicability and flexibility to a broader range of engineering scenarios. The proposed methodology's effectiveness and robustness are demonstrated through numerical examples, including fluidic pressure-loaded frequency-constraint structures undergoing small deformations, where compliance is minimized for structures optimized within specified resource constraints.