• 대한기계학회논문집A
• /
• 제24권5호
• /
• 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
• /
• 제34권4호
• /
• pp.449-474
• /
• 2010

A smart hollow cylinder consisting of a host functionally graded elastic core layer and two surface homogeneous piezoelectric layers is presented in this paper. The bonding between the layers can be perfect or imperfect, depending on the parameters taken in the general linear spring-layer interface model. The effect of such weak interfaces on free vibration and steady-state response is then investigated. Piezoelectric layers at inner and outer surfaces are polarized axially or radially and act as a sensor and an actuator respectively. For a simply supported condition, the state equations with non-constant coefficients are obtained directly from the formulations of elasticity/piezoelasticity. An approximate laminated model is then introduced for the sake of solving the state equations conveniently. It is further assumed that the hollow cylinder is embedded in an elastic medium and is simultaneously filled with compressible fluid. The interaction between the structure and its surrounding media is taken into account. Numerical examples are finally given with discussions on the effect of some related parameters.

• 한국소음진동공학회논문집
• /
• 제14권5호
• /
• pp.424-431
• /
• 2004

An implicit direct-time integration method for obtaining transient responses of general dynamic systems is described. The conventional Newmark method cannot be directly applied to state-space first-order differential equations, which contain no explicit acceleration terms. The method proposed here is the state-space Newmark method that incorporates the average velocity concept, and can be applied to an analysis of general dynamic systems that are expressed by state-space first-order differential equations. It is also readily coded into a program. Stability and accuracy analyses indicate that the method is numerically unconditionally stable like the conventional Newmark method, and has a period error of 2nd-order accuracy for small damping and 4th-order for large damping and an amplitude error of 2nd-order, regardless of damping. In addition, its utility and validity are confirmed by two application examples. The results suggest that the proposed state-space Newmark method based on average velocity be generally applied to the analysis of transient responses of general dynamic systems with a high degree of reliability with respect to stability and accuracy.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2003년도 추계학술대회논문집
• /
• pp.465-470
• /
• 2003

In this study, the state-space Newmark method based on average velocity is presented to analyse the transient dynamic response for general dynamic system. The conventional Newmark method based on average acceleration cannot he directly to the first-order state-space differential equations introducing the state-space vector. To overcome this problem, the time-step integration algorithm, based on average velocity concept, suitable for the first-order state-space differential equations is proposed In results, the proposed method has %he numerical stability and order of accuracy, which is proved analytically, equal to those of the conventional Newmark method based on average acceleration. Also, the formulation for numerical solution is very simple and the calculation time Is nearly equal to that of the conventional Newmark method based on average acceleration in spite of an increase of two times over matrix size. This method will be look forward to applying the general dynamic system to calculate the transient dynamic response.

• Journal of Advanced Research in Ocean Engineering
• /
• 제1권1호
• /
• pp.1-13
• /
• 2015

Post-Katrina investigations revealed that most earthen levee damage occurred on the levee crest and landward-side slope as a result of either wave overtopping, storm surge overflow, or a combination of both. In this paper, combined wave overtopping and storm surge overflow of a levee embankment strengthened with high performance turf reinforcement mat (HPTRM) system was studied in a purely Lagrangian and meshless approach, two-dimensional smoothed particle hydrodynamics (SPH) model. After the SPH model is calibrated with full-scale overtopping test results, the overtopping discharge, flow thickness, flow velocity, average overtopping velocity, shear stress, and soil erosion rate are calculated. New equations are developed for average overtopping discharge. The shear stresses on landward-side slope are calculated and the characteristics of soil loss are given. Equations are also provided to estimate soil loss rate. The range of the application of these equations is discussed.

• 한국자동차공학회논문집
• /
• 제15권2호
• /
• pp.165-174
• /
• 2007

Variable displacement type piston pump uses various controllers for controlling more than one state quantity like pressure, flow, power, and so on. These controllers need the mathematical model closely expressing dynamic behavior of pump for analyzing the stability of control systems which usually use various kinds of state variables. This paper derives the nonlinear mathematical model for variable displacement type piston pump. This model consists of two 1st oder differential equations by the continuity equations and one 2nd oder differential equation by the motion equation. To simplify the model we obtain the linear state variable model by differentiating the three nonlinear equations. And we verify this linearized model by comparison of simulation with experimentation and analyze the stability for the flow/pressure control. Finally this paper suggests the design compensation to ensure the stability of the systems.

• Advances in Computational Design
• /
• 제5권3호
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• Geomechanics and Engineering
• /
• 제20권5호
• /
• pp.411-420
• /
• 2020

Numbers fitting-curve equations have been proposed to predict soil-water retention curve (SWRC) whose parameters have no definitude physical meaning. And these methods with precondition of measuring SWRC data is time-consuming. A simplified directly method to estimate SWRC without parameters obtained by fitting-curve is proposed. Firstly, the total SWRC can be discretized into linear segments respectively. Every segment can be represented by linear formulation and every turning point can be determined by the pore-size distribution (PSD) of Mercury Intrusion Porosimetry (MIP) tests. The pore diameters governing the air-entry condition (AEC) and residual condition (RC) can be determined by the PSDs of MIP test. The PSD changes significantly during drying in SWR test, so the determination of AEC and RC should use the PSD under corresponding suction conditions. Every parameter in proposed equations can be determined directly by PSD without curve-fitting procedure and has definitude physical meaning. The proposed equations give a good estimation of both unimodal and bimodal SWRCs.

• 한국소음진동공학회논문집
• /
• 제16권1호
• /
• pp.57-65
• /
• 2006

In this paper, we present the principle of virtual work, from which the exact non-linear equations of motion of a rotating ring can be derived, by using the theory of finite deformation. For a thin ring of which the effect of variation in curvature across the cross-section is neglected, the radial displacement and the extensional stress are determined from the principle of virtual work at the steady state where the ring is rotating with a constant angular velocity. And also we formulate systematically the governing equations concerned to the in-plane vibrations and the out-of-plane vibrations at the disturbed state by using the principle of virtual work which is expressed with the disturbed displacements about the steady state. The formulated governing equations are classified by four models along the cases of considering or neglecting all or partly the secondary effects of flexural shear, rotary inertia, circumferential extension, and twist inertia. The natural vibrations of thin rings are analyzed, and its results are compared and discussed.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2005년도 추계학술대회논문집
• /
• pp.585-592
• /
• 2005

In this paper, we present the principle of virtual work, from which the exact non-linear equations of motion of a rotating ring can be derived, by using the theory of finite deformation For a thin ring of which the effect of variation in curvature across the cross-section is neglected, the radial displacement and the extensional stress are determined from the principle of virtual work at the steady state where the ring is rotating with a constant angular velocity. And also we formulate systematically the governing equations concerned to the in-plane vibrations and the out-of-plane vibrations at the disturbed state by using the principle of virtual work which is expressed with the disturbed displacements about the steady state. The formulated governing equations are classified by four models along the cases of considering or neglecting all or partly the secondary effects of flexural shear, rotary inertia, circumferential extension, and twist inertia. The natural vibrations of thin rings are analyzed, and its results are compared and discussed.