• Computers and Concrete
• /
• v.1 no.2
• /
• pp.131-150
• /
• 2004

We present in this work a coupled phenomenological chemo-mechanical model that represents the degradation of concrete-like materials. The chemical behaviour is described by the nowadays well known simplified calcium leaching approach. And the mechanical damage behaviour is described by a continuum damage model which involves the gradient of the damage quantity. The coupled nonlinear problem at hand is addressed within the context of the finite element method. For the equation governing the calcium dissolution-diffusion part of the problem, special care is taken to treat the highly nonlinear calcium conductivity and solid calcium functions. The algorithmic design is based on a Newton-type iterative scheme where use is made of a recently proposed relaxed linearization procedure. And for the equation governing the damage part of the problem, an augmented Lagrangian formulation is used to take into account the damage irreversibility constraint. Finally, numerical simulations are compared with experimental results on cement paste.

• Journal of Power System Engineering
• /
• v.12 no.1
• /
• pp.72-79
• /
• 2008

This paper designs a looper-tension controller for mass-flow stabilization in hot strip finishing mills. By Newton's 2nd law and Hooke's law, nonlinear dynamic equations on the looper-tension system are firstly derived, and linearized by a linearization algorithm using a Taylor's series expansion. Moreover, a tension calculation model is obtained from the nonlinear dynamic equations which is called as a soft sensor of strip tension between two neighboring stands. Next, a looper-tension servo controller is designed by an ILQ(Inverse Linear Quadratic optimal control) algorithm, and it is combined with a minimal disturbance observer which to attenuate speed disturbances by AGC and operator interventions, etc.. Finally, it is shown from by a computer simulation that the proposed ILQ controller with a disturbance observer is very effective in stabilizing the strip mass-flow under some disturbances, moreover it has a good command following performance.

• Structural Engineering and Mechanics
• /
• v.35 no.6
• /
• pp.677-697
• /
• 2010

This paper focuses on geometrically non-linear static analysis of a simply supported beam made of hyperelastic material subjected to a non-follower transversal uniformly distributed load. As it is known, the line of action of follower forces is affected by the deformation of the elastic system on which they act and therefore such forces are non-conservative. The material of the beam is assumed as isotropic and hyperelastic. Two types of simply supported beams are considered which have the following boundary conditions: 1) There is a pin at left end and a roller at right end of the beam (pinned-rolled beam). 2) Both ends of the beam are supported by pins (pinned-pinned beam). In this study, finite element model of the beam is constructed by using total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In order to use the solution procedures of Newton-Raphson type, there is need to linearized equilibrium equations, which can be achieved through the linearization of the principle of virtual work in its continuum form. In the study, the effect of the large deflections and rotations on the displacements and the normal stress and the shear stress distributions through the thickness of the beam is investigated in detail. It is known that in the failure analysis, the most important quantities are the principal normal stresses and the maximum shear stress. Therefore these stresses are investigated in detail. The convergence studies are performed for various numbers of finite elements. The effects of the geometric non-linearity and pinned-pinned and pinned-rolled support conditions on the displacements and on the stresses are investigated. By using a twelve-node quadratic element, the free boundary conditions are satisfied and very good stress diagrams are obtained. Also, some of the results of the total Lagrangian finite element model of two dimensional continuum for a twelve-node quadratic element are compared with the results of SAP2000 packet program. Numerical results show that geometrical nonlinearity plays very important role in the static responses of the beam.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Geotechnical Engineering
• /
• v.8 no.2
• /
• pp.5-20
• /
• 1992

The physical-based and lumped-parameter hydrologic groundwater flow model for predicting the rainfall-triggered rise of groundwater levels in hillside slopes is developed in this paper to assess the risk of landslides. The developed model consists of a vertical infiltration model for unsaturated zone linked to a linear storage reservoir model(LSRM) for saturated zone. The groundwater flow model has uncertain constants like soil depttL slope angle, saturated permeability, and potential evapotranspiration and four free model parameters like a, b, c, and K. The free model parameters could be estimated from known input-output records. The BARD algorithm is uses as the parameter estimation technique which is based on a linearization of the proposed model by Gauss -Newton method and Taylor series expansion. The application to examine the capacity of prediction shows that the developed model has a potential of use in forecast systems of predicting landslides and that the optimal estimate of potential 'a' in infiltration model is the most important in the global optimum analysis because small variation of it results in the large change of the objective function, the sum of squares of deviations of the observed and computed groundwater levels. 본 논문에서는 가파른 산사면에서 산사태의 발생을 예측하기 위한 수문학적 인 지하수 흐름 모델을 개발하였다. 이 모델은 물리적인 개념에 기본하였으며, Lumped-parameter를 이용하였다. 개발된 지하수 흐름 모델은 두 모델을 조합하여 구성되어 있으며, 비포화대 흐름을 위해서는 수정된 abcd 모델을, 포화대 흐름에 대해서는 시간 지체 효과를 고려할 수 있는 선형 저수지 모델을 이용하였다. 지하수 흐름 모델은 토층의 두께, 산사면의 경사각, 포화투수계수, 잠재 증발산 량과 같은 불확실한 상수들과 a, b, c, 그리고 K와 같은 자유모델변수들을 가진다. 자유모델변수들은 유입-유출 자료들로부터 평가할 수 있으며, 이를 위해서 본 논문에서는 Gauss-Newton 방법을 이용한 Bard 알고리즘을 사용하였다. 서울 구로구 시흥동 산사태 발생 지역의 산사면에 대하여 개발된 모델을 적용하여 예제 해석을 수행함으로써, 지하수 흐름 모델이 산사태 발생 예측을 위하여 이용할 수 있음을 입증하였다. 또한, 매개변수분석 연구를 통하여, 변수 a값은 작은 변화에 대하여 목적함수값에 큰 변화를 일으키므로 a의 값에 대한 최적값을 구하는 것이 가장 중요한 요소라는 결론을 얻었다.

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