• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.435-444
• /
• 2001

Even todays, accurate and efficient algorithms for the large deformation analysis of elastoplastic frame structures lack due to the complexities of kinematics, material nonlinearities and numerical methods to cater for. The author suggests appropriate beam element based upon the incremental formulation from the 3D rod theory where Cauchy stress and engineering strain are variables to incorporate plasticity equations so that objectivity may be satisfied. A rectum mapping methods which can integrate and satisfy yield criteria efficiently is suggested and a continuation method which has global convergency and quadratic speed is developed as well. leading-unloading example problems are tested and the ideas are proved to be valuable.

• Geomechanics and Engineering
• /
• v.13 no.4
• /
• pp.601-619
• /
• 2017

Cantilever retaining wall movements generally depend on the intensity and duration of ground motion, the response of the soil underlying the wall, the response of the backfill, the structural rigidity, and soil-structure interaction (SSI). This paper investigates the effect of material properties on seismic response of backfill-cantilever retaining wall-soil/foundation interaction system considering SSI. The material properties varied include the modulus of elasticity, Poisson's ratio, and mass density of the wall material. A series of nonlinear time history analyses with variation of material properties of the cantilever retaining wall are carried out by using the suggested finite element model (FEM). The backfill and foundation soil are modelled as an elastoplastic medium obeying the Drucker-Prager yield criterion, and the backfill-wall interface behavior is taken into consideration by using interface elements between the wall and soil to allow for de-bonding. The viscous boundary model is used in three dimensions to consider radiational effect of the seismic waves through the soil medium. In the seismic analyses, North-South component of the ground motion recorded during August 17, 1999 Kocaeli Earthquake in Yarimca station is used. Dynamic equations of motions are solved by using Newmark's direct step-by-step integration method. The response quantities incorporate the lateral displacements of the wall relative to the moving base and the stresses in the wall in all directions. The results show that while the modulus of elasticity has a considerable effect on seismic behavior of cantilever retaining wall, the Poisson's ratio and mass density of the wall material have negligible effects on seismic response.

• Geomechanics and Engineering
• /
• v.13 no.4
• /
• pp.585-600
• /
• 2017

This study deals with simple solutions for a spherical or circular opening excavated in elastic-brittle plastic rock mass compatible with a linear Mohr-Coulomb (M-C) or a nonlinear Hoek-Brown (H-B) yield criterion. Based on total strain approach, the closed-form solutions of stresses and displacement are derived simultaneously for circular and spherical openings using original H-B and M-C yield criteria. Two simple numerical procedures are proposed for the solution of generalized H-B and M-C yield criteria. Based on incremental approach, the similarity solution is derived for circular and spherical openings using generalized H-B and M-C yield criteria. The classical Runge-Kutta method is used to integrate the first-order ordinary differential equations. Using three data sets for M-C and H-B models, the results of the radial displacements, the spreading of the plastic radius with decreasing pressure, and the radial and circumferential stresses in the plastic region are compared. Excellent agreement among the solutions is obtained for all cases of spherical and circular openings. The importance of the use of proper initial values in the similarity solution is discussed.

• Smart Structures and Systems
• /
• v.14 no.6
• /
• pp.1131-1150
• /
• 2014

One of the issues in extending the range of applicable problems of real-time hybrid simulation is the computation speed of the simulator when large-scale computational models with a large number of DOF are used. In this study, functionality of real-time dynamic simulation of MDOF systems is achieved by creating a logic circuit that performs the step-by-step numerical time integration of the equations of motion of the system. The designed logic circuit can be implemented to an FPGA-based system; FPGA (Field Programmable Gate Array) allows large-scale parallel computing by implementing a number of arithmetic operators within the device. The operator splitting method is used as the numerical time integration scheme. The logic circuit consists of blocks of circuits that perform numerical arithmetic operations that appear in the integration scheme, including addition and multiplication of floating-point numbers, registers to store the intermediate data, and data busses connecting these elements to transmit various information including the floating-point numerical data among them. Case study on several types of linear and nonlinear MDOF system models shows that use of resource sharing in logic synthesis is crucial for effective application of FPGA to real-time dynamic simulation of structural response with time step interval of 1 ms.

• Journal of the Korea Concrete Institute
• /
• v.11 no.6
• /
• pp.13-24
• /
• 1999

The external, unbonded prestressed concrete(PSC) members exhibit very different structural behavior from that of internal bonded PSC members because of eccentricity change and slip occurrence during loading process. The purpose of the present study is to propose the ultimate failure stresses of prestressing (PS) steels for those external unbonded PSC members. To this end, a comprehensive analysis has been made using the nonlinear finite element analysis program developed recently for external unbonded PSC members by authors. A series of major influencing variables have been included in the analysis. It was found that the span-depth ratio, neutral axis depth-effective depth ratio, load geometry, amount of ordinary steel, and prestressing steel ration have great influence for the ultimate failue stress of PS steel is preposed and is compared with experimental dat as well as existing formulas for internal unbonded members. The Comparison indicates that the proposed equation agrees relatively well with experimental data and that existing formulas including ACI and AASHTO equations show some discrepancies from experimental ones. The present study allows more realistic analysis and design of prestressed concrete structures with external unbonded tendons.

• Proceedings of the KSME Conference
• /
• 2001.06d
• /
• pp.319-324
• /
• 2001

A finite element model for the process of squeeze casting for metal matrix composites (MMCs) in cylindrical mold is developed. The fluid flow and the heat transfer are the fundamental phenomena in the squeeze casing process. To describe heat transfer with solidification of molten aluminum, the energy equation in terms of temperature and enthalpy are applied to two dimensional axisymmetric model which is similar to the experimental system. And one dimensional flow model is employed to simulate the transient metal flow. The direct iteration technique was used to solve the resulting nonlinear algebraic equations. A computer program is developed to calculate the enthalpy, temperature and fluid velocity. Cooling curves and temperature distribution during infiltration and solidification are calculated for pure aluminum. The temperature is measured and recorded experimentally. At two points of the perform inside and one point of the mold outside, thermocouple wire are installed. The time-temperature data are compared with the calculated cooling curves. The experimental results show that the finite element model can estimate the solidification time and predict the cooling process.

• Journal of the Korean Mathematical Society
• /
• v.35 no.2
• /
• pp.465-490
• /
• 1998

In this paper, we consider the existence and asymptotic behavior of solutions of the following problem: $u_{tt}$ -(t, x) - (∥∇u(t, x)∥(equation omitted) ＋ ∥∇v(t, x) (equation omitted)$^{\gamma}$ $\Delta$u(t, x)＋$\delta$│ $u_{t}$ (t, x)│sup p-1/ $u_{t}$ (t, x) = $\mu$│u(t, x) $^{q-1}$u(t, x), x$\in$$\Omega$, t$\in$[0, T], $v_{tt}$ (t, x) - (∥∇uu(t, x) (equation omitted) ＋ ∥∇v(t, x) (equation omitted)sup ${\gamma}$/ $\Delta$v(t, x)＋$\delta$ │ $v_{t}$ (t, x)│sup p-1/ $u_{t}$ (t, x) = $\mu$ u(t, x) $^{q-1}$u(t, x), x$\in$$\Omega$, t$\in$[0, T], u(0, x) = $u_{0}$ (x), $u_{t}$ (0, x) = $u_1$(x), x$\in$$\Omega$, u(0, x) = $v_{0}$ (x), $v_{t}$ (0, x) = $v_1$(x), x$\in$$\Omega$, u│∂$\Omega$=v│∂$\Omega$=0 T > 0, q > 1, p $\geq$1, $\delta$ > 0, $\mu$ $\in$ R, ${\gamma}$ $\geq$ 1 and $\Delta$ is the Laplacian in $R^{N}$.X> N/.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.20 no.4
• /
• pp.1-8
• /
• 1983

This paper describes the application of the finite element method to the large deflection elastic plastic analysis and ultimate load calculation of axisymmetric shell of revolution with initial imperfection subjected to external pressure. The nonlinear equilibrium equations are linearized by the successive incremental method and are solved by the combination of load increment and iteration scheme with considering plastic deformation theory. To get the more realistic effect of large deflection, corrected coordinats and directions of applied load ar every load increment steps are used. The effects of the plasticity, initial imperfection and the shape of shells on the ultimate load of clamped circular cap under external pressure are investigated. Consequently, the following conclusions are obtained; (1) At same geometric parameter $\lambda$, each shape of clamped circular caps yield same elastic ultimate loads in both cases, i.e. with and without initial imperfections, whereas, in the case of elastic-plastic state the shell becomes thicker, the ultimate loads are getting smaller. (2) The effects of initial imperfection to ultimate load are most significant in the elastic case and are more senstive in the elastic-plastic state with the thinner shells.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.2
• /
• pp.14-26
• /
• 2002

A new model for degenerate semiconductor quantum well devices was developed. In this model, the multi-subband Boltzmann transport equation was formulated by applying the Pauli exclusion principle and coupled to the Schrodinger and Poisson equations. For the solution of the resulted nonlinear system, the finite difference method and the Newton-Raphson method was used and carrier energy distribution function was obtained for each subband. The model was applied to a Si MOSFET inversion layer. The results of the simulation showed the changes of the distribution function from Boltzmann like to Fermi-Dirac like depending on the electron density in the quantum well, which presents the appropriateness of this modeling, the effectiveness of the solution method, and the importance of the Pauli -exclusion principle according to the reduced size of semiconductor devices.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10b
• /
• pp.1904-1909
• /
• 1991

Software linear and exponential acceleration/deceleration algorithms for control of machine axes of motion in industrial robots and CNC machine tools are proposed. Typical hardware systems used to accelerate and decelerate axes of motion are mathematically analyzed. Discrete-time state equations are derived from the mathematical analyses for the development of software acceleration/deceleration algorithms. Synchronous control method of multiple axes of motion in industrial robots and CNC machine tools is shown to be easily obtained on the basis of the proposed acceleration/deceleration algorithms. The path error analyses are carried out for the case where the software linear and exponential acceleration/deceleration algorithms are applied to a circular interpolator. A motion control system based on a floating point digital signal processor (DSP) TMS 320C30 is developed in order to implement the proposed algorithms. Experimental results demonstrate that the developed algorithms and the motion control system are available for control of multiple axes and nonlinear motion composed of a combination of lines and circles which industrial robots and CNC machine tools require.