• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• International Journal of Highway Engineering
• /
• v.11 no.1
• /
• pp.165-175
• /
• 2009

With the current move towards adopting mechanistic-empirical concepts in the design of pavement structures, state-of-the-art mechanistic analysis methodologies are needed to determine accurate pavement responses, such as stress, strain, and deformation. Previous laboratory studies of pavement foundation geomaterials, i.e., unbound granular materials used in base/subbase layers and fine-grained soils of a prepared subgrade, have shown that the resilient responses followed by nonlinear, stress-dependent behavior under repeated wheel loading. This nonlinear behavior is commonly characterized by stress-dependent resilient modulus material models that need to be incorporated into finite element (FE) based mechanistic pavement analysis methods to predict more realistically predict pavement responses for a mechanistic pavement analysis. Developed user material subroutine using aforementioned resilient model with nonlinear solution technique and convergence scheme with proven performance were successfully employed in general-purpose FE program, ABAQUS. This numerical analysis was investigated in predicted critical responses and domain selection with specific mesh generation was implemented to evaluate better prediction of pavement responses. Results obtained from both axisymmetric and three-dimensional (3D) nonlinear FE analyses were compared and remarkable findings were described for nonlinear FE analysis. The UMAT subroutine performance was also validated with the instrumented full scale pavement test section study results from the Federal Aviation Administration's National Airport Pavement Test Facility (FAA's NAPTF).

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.1
• /
• pp.36-43
• /
• 2002

Thermopiezoelastic snap-through phenomena of piezolaminated plates are investigated by applying an are-length scheme to Newton-Raphson method. Based on the layerwise displacement theory and von Karman strain-displacement relationships, nonlinear finite element formulations are derived for the thermopiezoelastic composite plates. From the static and dynamic viewpoint, nonlinear thermopierzoelastic behavior and vibration characteristicx are stuied for symmetric and eccentric structural models with various piezoelestric actuation modes. Present results show the possibility to enhance the performance, namely thermopiezoelastic snapping, induced by the excessive piezoelectric actuation in the active suppression of thermally buckled large deflection piezolaminated paltes.

• Journal of Mechanical Science and Technology
• /
• v.17 no.12
• /
• pp.1835-1843
• /
• 2003

This paper introduces the development of the caliper system for a geometry PIG (Pipeline Inspection Gauge). The objective of the caliper system is to detect and measure dents, wrinkles, and ovalities affect the pipe structural integrity. The developed caliper system consists of a finger arm, an anisotropic permanent magnet, a back yoke, pins, pinholes and a linear hall effect sensor. The angle displacement of the finger arm is measured by the change of the magnetic field in sensing module. Therefore the sensitivity of the caliper system mainly depends on the magnitude of the magnetic field inside the sensing module. In this research, the ring shaped anisotropic permanent magnet and linear hall effect sensors were used to produce and measure the magnetic field. The structure of the permanent magnet, the back yoke and pinhole positions were optimized that the magnitude of the magnetic field range between a high of 0.1020 Tesla and a low of zero by using three dimensional nonlinear finite element methods. A simulator was fabricated to prove the effectiveness of the developed caliper system and the computational scheme using the finite element method. The experimental results show that the developed caliper system is quite efficient for the geometry PIG with good performance.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.3
• /
• pp.403-414
• /
• 1994

Nonlinear analysis of RC containment structure under thermal load and pressure is presented to trace the behaviour after an assumed LOCA. The temperature distribution varying with time through the wall thickness is determined by transient finite element analysis with the two time level scheme in time domain. The layered shell finite elements are used to represent the containment structures in nuclear power plants. Both geometric and material nonlinearities are taken into account in the finite element formulation. The constitutive relation of concrete is modeled according to Drucker-Prager yield criteria in compression. Tension stiffening model is used to represent the tensile behaviour of concrete including bond effect. The reinforcing bars are modeled by smeared layer at the location of reinforcements accounting elasto-plastic axial behaviors. The steel liner model under Von Mises yield criteria is adopted to represent elastic-perfect plastic behaviour. Geometric nonlinearity is formulated to consider the large displacement effect. Thermal stress components are determined by the initial strain concept during each time step. The temperature differential between any two consecutive time steps is considered as a load incremental. The numerical results from this study reveal that nonlinear temperature gradient based on transient thermal analysis will produces excessive large displacement. Nonlinear behavior of containment structures up to ultimate stage can be traced reallistically. The present study allows more realistic analysis of concrete containment structures in nuclear power plants.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.247-250
• /
• 2007

Flexible media such as the paper, the film, etc. are thin, light and very flexible. They behave in geometrically nonlinear. Any of small force makes large deformation. So we must including aerodynamic effect when its behavior is predicted. Thus, it becomes fully coupled fluid-structure interaction(FSI) problem. In FSI problems, where the fluid mesh near the structure undergoes large deformations and becomes unacceptably distorted, which drive the time step to a very small value for explicit calculations, the arbitrary Lagrangian-Eulerian(ALE) methods or rezoning are used to create a new undistorted mesh for the fluid domain, which allows the calculations to continue. In this paper, FE sheet model considering geometric nonlinearity is formulated to simulate the behavior of the flexible media. Aerodynamic force to the media by surrounding air is calculated by solving the incompressible Navier-Stokes equations. Q2Q1(Taylor-Hood) element which means biquadratic for velocity and bilinear for pressure is used for fluid domain. Q2Q1 element satisfies LBB condition and any stabilization technique is not needed. In this paper, cantilevered sheet in the viscous incompressible Navier-Stokes flow is simulated to check the mesh motion and numerical integration scheme, and then falling paper in the air is simulated and the effects of some representative parameters are investigated.

• International Journal of Control, Automation, and Systems
• /
• v.3 no.2
• /
• pp.211-216
• /
• 2005

In electrical impedance tomography (EIT), various image reconstruction algorithms have been used in order to compute the internal resistivity distribution of the unknown object with its electric potential data at the boundary. Mathematically, the EIT image reconstruction algorithm is a nonlinear ill-posed inverse problem. This paper presents a simulated annealing technique as a statistical reconstruction algorithm for the solution of the static EIT inverse problem. Computer simulations with 32 channels synthetic data show that the spatial resolution of reconstructed images by the proposed scheme is improved as compared to that of the mNR algorithm at the expense of increased computational burden.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.159-164
• /
• 2004

In electrical impedance tomography (EIT), various image reconstruction algorithms have been used in order to compute the internal resistivity distribution of the unknown object with its electric potential data at the boundary. Mathematically the EIT image reconstruction algorithm is a nonlinear ill-posed inverse problem. This paper presents a simultaneous perturbation method as an image reconstruction algorithm for the solution of the static EIT inverse problem. Computer simulations with the 32 channels synthetic data show that the spatial resolution of reconstructed images by the proposed scheme is improved as compared to that of the mNR algorithm at the expense of increased computational burden.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.10 no.3
• /
• pp.176-184
• /
• 2002

The door stiffness is one of the important factors for the side impact. Generally, the researches have been conducted on the assembled door. A side impact door beam is installed in a door to protect occupants from the side impact. This research is only concentrated on the side impact beam and a side impact beam is designed. The cross section is defined to have an elliptic shape. An optimization problem is defined to find the design maximizing the intrusion stiffness within the specified weight. Design variables are the radii and the thickness of the ellipsoid. The analysis of the side impact is carried out by the nonlinear finite element method. The optimization problem is solved by two methods. One is the experimental design scheme using an orthogonal array. The other is the gradient-based optimization using the response surface method(RSM). Both methods have obtained the better designs than the current one.