• Journal of the Korean Wood Science and Technology
• /
• v.42 no.4
• /
• pp.367-375
• /
• 2014

A simplified material model, which was efficiently implemented in a three-dimensional finite solid element (3D FE) analysis for wood was developed. The bi-linear elasto-plastic anisotropic material theory was adopted to describe constitutive relations of wood in three major directions including longitudinal, radial and tangential direction. The assumption of transverse isotropy was made to reduce the requisite 27 material constants to 6 independent constants including elastic moduli, yield stresses and Poisson's ratios in the parallel, and perpendicular to grain directions. The results of Douglas fir compression tests in the three directions were compared to the 3D FE simulation incorporated with the wood material model developed in this study. Successful agreements of the results were found in the load-deformation curves and the permanent deformations. Future works and difficulties expected in the advanced application of the model were discussed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.7
• /
• pp.44-50
• /
• 2002

Low-velocity impact on composite sandwich panel has been investigated. Contact force is computed from a proposed modified Hertzian contact law. The Hertzian contact law is constructed by adjusting numerical value of the exponent and reducing the through-the- thickness elastic constant of honeycomb core. The equivalent transverse elastic constant is calculated from the rule of mixture. Nonlinear equation to calculate the contact force is solved by the Newton-Raphson method and time integration is done by the Newmark-beta method. A finite element program for the low-velocity impact analysis is coded by implementing these techniques and an 18-node assumed strain solid element. Behaviors of composite sandwich panels subjected to low-velocity impact are analyzed for various cases with different geometry and lay-ups. It has been found that the present code with the proposed contact law can predict measured contact forces and contact times for most cases within reasonable error bounds.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.10 no.9
• /
• pp.2183-2188
• /
• 2009

This paper describes an automatic mesh generation system for finite element analysis of three-dimensional cracks. It is consisting of fuzzy knowledge processing, bubble meshing and solid geometry modeler. This novel mesh generation process consists of three sub-processes: (a) definition of geometric model, i.e. analysis model, (b) generation of bubbles, and (c) generation of elements. One of commercial solid modelers is employed for three-dimensional crack structures. Bubble is generated if its distance from existing bubble points is similar to the bubble spacing function at the point. The bubble spacing function is well controlled by the fuzzy knowledge processing. The Delaunay method is introduced as a basic tool for element generation. Practical performances of the present system are demonstrated through several mesh generations for 3D cracks.

• Nuclear Engineering and Technology
• /
• v.53 no.8
• /
• pp.2696-2707
• /
• 2021

The purpose of this study is to investigate the efficiency of various structural modeling schemes for evaluating seismic performances and fragility of the reactor containment building (RCB) structure in the advanced power reactor 1400 (APR1400) nuclear power plant (NPP). Four structural modeling schemes, i.e. lumped-mass stick model (LMSM), solid-based finite element model (Solid FEM), multi-layer shell model (MLSM), and beam-truss model (BTM), are developed to simulate the seismic behaviors of the containment structure. A full three-dimensional finite element model (full 3D FEM) is additionally constructed to verify the previous numerical models. A set of input ground motions with response spectra matching to the US NRC 1.60 design spectrum is generated to perform linear and nonlinear time-history analyses. Floor response spectra (FRS) and floor displacements are obtained at the different elevations of the structure since they are critical outputs for evaluating the seismic vulnerability of RCB and secondary components. The results show that the difference in seismic responses between linear and nonlinear analyses gets larger as an earthquake intensity increases. It is observed that the linear analysis underestimates floor displacements while it overestimates floor accelerations. Moreover, a systematic assessment of the capability and efficiency of each structural model is presented thoroughly. MLSM can be an alternative approach to a full 3D FEM, which is complicated in modeling and extremely time-consuming in dynamic analyses. Specifically, BTM is recommended as the optimal model for evaluating the nonlinear seismic performance of NPP structures. Thereafter, linear and nonlinear BTM are employed in a series of time-history analyses to develop fragility curves of RCB for different damage states. It is shown that the linear analysis underestimates the probability of damage of RCB at a given earthquake intensity when compared to the nonlinear analysis. The nonlinear analysis approach is highly suggested for assessing the vulnerability of NPP structures.

• Coupled systems mechanics
• /
• v.8 no.1
• /
• pp.71-93
• /
• 2019

Finite element simulations of solid mechanics problems often involve the use of Non-Confirming Meshes (NCM) to increase accuracy in capturing nonlinear behavior, including damage and plasticity, in part of a solid domain without an undue increase in computational costs. In the presence of material nonlinearity and plasticity, higher-order variables are often needed to capture nonlinear behavior and material history on non-conforming interfaces. The most popular formulations for coupling non-conforming meshes are dual methods that involve the interpolation of a traction field on the interface. These methods are subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) stability condition, and are therefore limited in their implementation with the higher-order elements needed to capture nonlinear material behavior. Alternatively, the enriched discontinuous Galerkin approach (EDGA) (Haikal and Hjelmstad 2010) is a primal method that provides higher order kinematic fields on the interface, and in which interface tractions are computed from local finite element estimates, therefore facilitating its implementation with nonlinear material models. The inclusion of higher-order interface variables, however, presents the issue of preserving material history at integration points when a increase in integration order is needed. In this study, the enriched discontinuous Galerkin approach (EDGA) is extended to the case of small-deformation plasticity. An interface-driven Gauss-Kronrod integration rule is proposed to enable adaptive enrichment on the interface while preserving history-dependent material data at existing integration points. The method is implemented using classical J2 plasticity theory as well as the pressure-dependent Drucker-Prager material model. We show that an efficient treatment of interface variables can improve algorithmic performance and provide a consistent approach for coupling non-conforming meshes in inelasticity.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.26 no.2
• /
• pp.71-82
• /
• 2022

This paper presents the effect on the inelastic behavior and structural performance of concrete and filled steel pipe through a numerical method for reliable judgment under various load conditions of the CJS composite structural system. Variable values optimized for the CJS synthetic structural system and the effects of multiple variables used for finite element analysis to present analytical modeling were compared and analyzed with experimental results. The Winfrith concrete model was used as a concrete material model that describes the confinement effect well, and the concrete structure was modeled with solid elements. Through geometric analysis of shell and solid elements, rectangular steel pipe columns and steel elements were modeled as shell elements. In addition, the slip behavior of the joint between the concrete column and the rectangular steel pipe was described using the Surface-to-Surface function. After finite element analysis modeling, simulation was performed for cyclic loading after assuming that the lower part of the foundation was a pin in the same way as in the experiment. The analysis model was verified by comparing the calculated analysis results with the experimental results, focusing on initial stiffness, maximum strength, and energy dissipation capability.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.10a
• /
• pp.133-138
• /
• 1991

The aim of the present paper is to develop a computer program predicting ultimate fracture strength of initially cracked structure under monotonically increasing external loads. For this purpose, two kinds of 3-D isoparametric solid elements, one 6-node wedge element and another 8-node brick element are formulated along the small deformation theory. Plasticity in the element is checked using von Mises' yield criterion. Elasto-plastic stiffness matrix of the element is calculated taking account of strain hardening effect. If the principal strain at crack tip which is one nodal point exceeds the critical strain dependin on the material property, crack tip is supposed to be opened and the crack tip node which was previously constrained in the direction perpendicular to the crack line is released. After that, the crack lay be propagated to the adjacent node. Once a crack tip node is fractured, the energy of the newly fractured node should be released which is to be absorbed by the remaining part. The accumulated reaction force which was carried by the newly fractured node so far is then applied in the opposite direction. During the action of crack tip relief force, since unloading may be occured in the plastic element, unloading check should be made. If a plastic element unloads, elastic stress-strain equation is used in the calculation of the stiffness matrix of the element, while for a loading element, elasto-plastic stress-strain equation is continuously used. Verification of the computer program is made comparing with the experimental results for center cracked panel subjected to uniform tensile load. Also some factors affecting ultimate fracture strength of initially cracked plate are investigated. It is concluded that the computer program developed here gives an accurate solution and becomes useful tool for predicting ultimate fracture load of initially cracked structural system under monotonically increasing external loads.

• Geotechnical Engineering
• /
• v.12 no.3
• /
• pp.35-48
• /
• 1996

Reinforced soil structures may experience large local movements between soil and reinforcement. The failure modes of a reinforced structure depend on several factors which are governed by deformation and slipping of the reinforcement. In some cases, pulling out of the reinforcement may occur instead of rupturing, The growing use of geosynthetic liner system for storage of solid and liquid wastes has led to a number of slope instability problems where the synthetic liner may undergo a large amount of stretching and slipping as a result of the loading. The conventional finite element model for the soil-reinforcement interface uses a zero thickness joint element with normal and shear stiffnesses and can only accommodate a small amount of deformation. When a large slippage occurs, the model provides an i ncorrect mechanism for deformation. This paper presents a new interface finite element model which is able to simulate a large amount of slippage between soil and reinforcement. The formulation of the model is presented and the capability of the model is demonstrated using illustrative examples.

• Structural Engineering and Mechanics
• /
• v.39 no.6
• /
• pp.767-793
• /
• 2011

An efficient edge-based smoothed finite element method (ES-FEM) has been recently developed for solving solid mechanics problems. The ES-FEM uses triangular elements that can be generated easily for complicated domains. In this paper, the complexity study of the ES-FEM based on triangular elements is conducted in detail, which confirms the ES-FEM produces higher computational efficiency compared to the FEM. Therefore, the ES-FEM offers an excellent platform for adaptive analysis, and this paper presents an efficient adaptive procedure based on the ES-FEM. A smoothing domain based energy (SDE) error estimate is first devised making use of the features of the ES-FEM. The present error estimate differs from the conventional approaches and evaluates error based on smoothing domains used in the ES-FEM. A local refinement technique based on the Delaunay algorithm is then implemented to achieve high efficiency in the mesh refinement. In this refinement technique, each node is assigned a scaling factor to control the local nodal density, and refinement of the neighborhood of a node is accomplished simply by adjusting its scaling factor. Intensive numerical studies, including an actual engineering problem of an automobile part, show that the proposed adaptive procedure is effective and efficient in producing solutions of desired accuracy.

• Structural Engineering and Mechanics
• /
• v.12 no.2
• /
• pp.201-214
• /
• 2001

Four finite element (FE) models are examined to find the one that best estimates moment-rotation characteristics of top- and seat-angle with double web-angle connections. To efficiently simulate the real behavior of connections, finite element analyses are performed with following considerations: 1) all components of connection (beam, column, angles and bolts) are discretized by eight-node solid elements; 2) shapes of bolt shank, head, and nut are precisely taken into account in modeling; and 3) contact surface algorithm is applied as boundary condition. To improve accuracy in predicting moment-rotation behavior of a connection, bolt pretension is introduced before the corresponding connection moment being surcharged. The experimental results are used to investigate the applicability of FE method and to check the performance of three-parameter power model by making comparison among their moment-rotation behaviors and by assessment of deformation and stress distribution patterns at the final stage of loading. This research exposes two important features: (1) the FE method has tremendous potential for connection modeling for both monotonic and cyclic loading; and (2) the power model is able to predict moment-rotation characteristics of semi-rigid connections with acceptable accuracy.