• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1075-1083
• /
• 2000

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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. 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. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented 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 line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Structural Engineering and Mechanics
• /
• v.82 no.1
• /
• pp.121-131
• /
• 2022

The main goal of this study is to prepare a program for analyzing High Strength Steel Fibrous Reinforced Concrete (HSSFRC) slabs and predict the response and strength of the slab instead of preparing a prototype and testing it in the laboratory. For this purpose, new equations are proposed to represent the material properties of High Strength Steel Fibrous Reinforced Concrete. The proposed equations obtained from performing regression analysis on many experimental results using statistical programs. The finite element method is adopted for non-linear analysis of the slabs. The eight-node "Serendipity element" (3 DoF) is chosen to represent the concrete. The layered approach is adopted for concrete elements and the steel reinforcement is represented by a smeared layer. The compression properties of the concrete are modeled by a work hardening plasticity approach and the yield condition is determined depending on the first two stress invariants. A tensile strength criterion is adopted in order to estimate the cracks propagation. many experimental results for testing slabs are compared with the numerical results of the present study and a good agreement is achieved regarding load-deflection curves and crack pattern. The response of the load deflection curve is slightly stiff at the beginning because the creep effect is not considered in this study and for assuming perfect bond between the steel reinforcement and the concrete, however, a great agreement is achieved between the ultimate load from the present study and experimental results. For the models of the tension stiffening and cracked shear modulus, the value of Bg and Bt (Where Bg and Bt are the curvature factor for the cracked shear modulus and tension stiffening models respectively) equal to 0.005 give good results compared with experimental result.

• Tunnel and Underground Space
• /
• v.20 no.2
• /
• pp.97-104
• /
• 2010

The materials failure relation $\ddot{\Omega}=A{(\dot{\Omega})}^\alpha$ where $\Omega$ is a measurable quantity such as displacement and the dot superscript is the time derivative, may be used to analyze the accelerating creep of materials. Coefficients, A and $\alpha$, are determined by fitting given data sets. In this study, it is tried to predict the failure time of tunnel using the materials failure relation. Four fitting techniques of applying the materials failure relation are attempted to forecast a failure time. Log velocity versus log acceleration technique, log time versus log velocity technique, inverse velocity technique are based on the linear least squares fits and non-linear least squares technique utilizes the Levenberg-Marquardt algorithm. Since the log velocity versus log acceleration technique utilizes a logarithmic representation of the materials failure relation, it indicates the suitability of the materials failure relation applied to predict a failure time of tunnel. A linear correlation between log velocity and log acceleration appears satisfactory(R=0.84) and this represents that the materials failure relation is a suitable model for predicting a failure time of tunnel. Through comparing the real failure time of tunnel with the predicted failure times from four curve fittings, it is shown that the log time versus log velocity technique results in the best prediction.

• JSTS:Journal of Semiconductor Technology and Science
• /
• v.2 no.4
• /
• pp.246-252
• /
• 2002

Two kinds of PZT cantilevers integrated with a piezoresistor have been newly designed, fabricated, and characterized for high speed AFM. In first cantilever, a piezoresistor is used to sense atomic force acting on tip, while in second cantilever, a piezoresistor is integrated to calibrate hysteresis and creep phenomena of the PZT cantilever. The fabricated PZT cantilevers provide high tip displacement of $0.55\mu\textrm{m}/V$ and high resonant frequency of 73 KHz. A new cantilever structure has been designed to prevent electrical coupling between sensor and PZT actuator and the proposed cantilever shows 5 times lower coupling voltage than that of the previous cantilever. The fabricated PZT cantilever shows a crisp scanned image at 1mm/sec, while the conventional piezo-tube scanner shows blurred image even at $180\mu\textrm{m}/sec$. The non-linear properties of the PZT actuator are also well calibrated using the piezoresistive sensor for calibration.

• Journal of the Korea Concrete Institute
• /
• v.25 no.5
• /
• pp.497-508
• /
• 2013

This paper introduces a numerical model which can evaluate the fire-resistant capacity of reinforced concrete members. On the basis of the transient heat transfer considering the heat conduction, convection and radiation, time-dependent temperature distribution across a section is determined. A layered fiber section method is adopted to consider non-linear material properties depending on the temperature and varying with the position of a fiber. Furthermore, effects of non-mechanical strains of each fiber like thermal expansion, transient strain and creep strain are reflected on the non-linear structural analysis to take into account the extreme temperature variation induced by the fire. Analysis results by the numerical model are compared with experimental data from the standard fire tests to validate an exactness of the introduced numerical model. Also, time-dependent changes in the resisting capacities of reinforced concrete members exposed to fire are investigated through the analyses and, the resisting capacities evaluated are compared with those determined by the design code.

• 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.

• Nuclear Engineering and Technology
• /
• v.20 no.2
• /
• pp.88-104
• /
• 1988

The FUel Rod Analysis(FURA) code is developed using two-dimensional finite element methods for axisymmetric and plane stress analysis of fuel rod. It predicts the thermal and mechanical behavior of fuel rod during normal and load follow operations. To evaluate the exact temperature distribution and the inner gas pressure, the radial deformation of pellet and clad, the fission gas release are considered over the full-length of fuel rod. The thermal element equation is derived using Galerkin's techniques. The displacement element equation is derived using the principle of virtual works. The mechanical analysis can accommodate various components of strain: elastic, plastic, creep and thermal strain as well as strain due to swelling, relocation and densification. The 4-node quadratic isoparametric elements are adopted, and the geometric model is confined to a half-pellet-height region with the assumption that pellet-pellet interaction is symmetrical. The pellet cracking and crack healing, pellet-cladding interaction are modelled. The Newton-Raphson iteration with an implicit algorithm is applied to perform the analysis of non-linear material behavior accurately and stably. The pellet and cladding model has been compared with both analytical solutions and experimental results. The observed and predicted results are in good agreement. The general behavior of fuel rod is calculated by axisymmetric system and the cladding behavior against radial crack is used by plane stress system. The sensitivity of strain aging of PWR fuel cladding tube due to load following is evaluated in terms of linear power, load cycle frequency and amplitude.