• Transactions of Materials Processing
• /
• v.9 no.6
• /
• pp.590-597
• /
• 2000

A finite element formulation lot the simulation of axisymmetric sheet hydroforming is proposed, and an implicit program is coded. In order to describe normal anisotropy of steel sheet, Hill's non-quadratic yield function (Hill, 1979) is employed. Frictional contacts among sheet surface, rigid tool surface, and flexible hydrostatic pressure are considered using mesh normal vectors based on finite element of the sheet. Applied hydraulic pressure is also considered as a function of forming rate and time and treated as an external loading. The complete set of the governing relations comprising equilibrium and interfacial equations is approximately linearized for Newton-Raphson algorithm. In order to verify the validity of the developed finite element formulation, the axisymmetric bulge test is simulated. Simulation results are compared with other FEM results and experimental measurements and showed good agreements. In axisymmetric hydroforming processes of a disk cover, formability changes are observed according to the hydraulic pressure curve changes.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2009.05a
• /
• pp.172-175
• /
• 2009

The purpose of this study is to predict the forming limit diagram (FLD) of strain-rate sensitive materials on the basis of the Marciniak and Kuczynski (M-K) theory. The strain-rate effect is taken into consideration in such a way that the stress-strain curves for various strain-rates are inputted into the formulation as point data, not as curve-fitted models such as power function. To solve the nonlinear system of equations derived from the equilibrium and constraints in the groove region and the safe zone, the Newton-Raphson method is used. The theoretical FLDs using four different yield criteria, that are von Mises, Hill (1948), Hill (1979), Logan and Hosford, are compared with the experimental, numerical (FEA) and other theoretical results. A new trial is made where a modified M-K model having n-step grooves is introduced to describe a real localized neck.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Journal of the Korean Solar Energy Society
• /
• v.21 no.3
• /
• pp.51-60
• /
• 2001

A greenhouse equipped with latent heat storage system was built to obtain various thermal properties, such as greenhouse air temperature, soil surface temperature, energy flow in latent heat storage, etc., which could be used in validation of greenhouse numerical model to be developed in this study. This numerical model expressed with Newton-Raphson method was programed by C-language and utilized to simulate greenhouse thermal behavior. Greenhouse air temperature and soil surface temperature predicted by the greenhouse model developed in this study were very close to the measured data obtained through almost 3 years of experiment. Therefore, it is concluded that the greenhouse model developed and verified by measured data could be utilized for simulating various thermal behaviors of greenhouses equipped with latent heat storage system to be used for energy saving purposes.

• Proceedings of the KIEE Conference
• /
• 1997.07a
• /
• pp.306-308
• /
• 1997

In this paper, We present the transient analysis method of induction motor by TDFE(Time Domain Finite Element) method. For simulation of transient performance, Maxwell's equations are solved using 2-Dimensional TDFE method, and the circuit equations from the stator and rotor are solved simultaneously. The time derivatives are discretized with Euler scheme and the Newton-Raphson iteration method is applied to a large system of equations which are representing the whole magnetic and feeding circuit equations because of the magnetic nonlinearity of the stator and rotor core. The presented method is applied to three phase induction motor. And we obtained the phase currents, torque and rotor position until the steady state.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.15 no.5
• /
• pp.389-396
• /
• 2003

Systematic methods to calculate thermodynamic properties of carbon dioxide, HFC-l34a and HCFC-22 are presented. First, application of a basic method to identify the saturation state with given temperature or pressure is attempted and the feasibility of auxiliary equations is tested. Next, detailed procedures are suggested to tell a phase when temperature/pressure and another property are specified. Finally Newton-Raphson method is applied to calculate unknown thermodynamic properties fixing the state with the two independent properties specified. The procedures described here are utilized to develop a computer program, which is used to find the relation between temperature and pressure with maximum isobaric heat capacity for super-critical carbon dioxide.

• Journal of Ocean Engineering and Technology
• /
• v.2 no.2
• /
• pp.71-77
• /
• 1988

For the collapse of structures due to the variable repeated load, two types of collapse mechanisms, i.e., incremental collapse and alternating plasticity, exist. Under the similar variable repeated loading conditions there exists shakedown state in the structures. In shakedown state, the number of plastic hinges are not increased and all further loading will be resulted in the elastic moment changes. Namely, under the shakedown state, structures do not collapse. In this investigation, shakedown analysis are performed by composing new computer programs. Basic theories employed to compose the programs are as follows. 1. Newton-Raphson methods are added to the existing matrix method for the plastic analysis. 2. An effort to construct the stiffness of axial and bending springs attached at both ends of the member has been made. By using the programs developed, it is possible to anticipate the collapse mechanisms (Incremental collapse, alternating plasticity). Lastly for the verification of performance of the program, demonstration examples have been solved and the results are compared with other sources.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.54 no.6
• /
• pp.267-273
• /
• 2005

This paper proposes a line-to-ground fault location algorithm for underground cable system. A feature of the proposed method is a new algorithm based on the analytic research which has not been tried until now. The proposed method firstly makes voltage and current equations using analysis of distributed parameter circuit for each of cores and sheathes respectively, and then establishes an equation of the fault distance according to the analysis of the fault conditions. Finally the solution of this equation is calculated by Newton-Raphson iteration method. The effectiveness of this proposed algorithm has been proven through PSCAD/EMTDC(Ver. 4.1) simulations.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.49 no.11
• /
• pp.542-548
• /
• 2000

In this paper, we proposed a new load flow analysis algorithm. In order to develop it, we converted typical power flow problem into optimal problem. This problem has the objective function that minimize the difference between calculated values and specified values of bus powers and subject to bus power equations of P and Q. Using it, we solved the divergence by singularity of Jacobian matrix, the divergence by initial value in the typical power flow study. In the study of a sample system, we verified the superiority of proposed algorithm.

• Proceedings of the KIEE Conference
• /
• 1999.07a
• /
• pp.70-72
• /
• 1999

일반적으로 비선형 정자계 문제를 해석하기 위해서 뉴튼-�N슨(Newton -Raphson : NR)법이 이용된다. 하지만 뉴튼-�N슨법의 경우 각 반복계산 때마다 새로운 선형 시스템의 해를 구하기 위해서 LU-decomposition과 같은 과정을 매 반복계산 때마다 시행해야 하므로 절점(node)의 수가 증가할 경우 계산시간이 증가한다는 단점이 있다. 이러한 단점을 보완하기 위해서 최근 TLM (Transmission Line Modeling)법이 새로운 반복계산법으로 비선형 유한 요소 해석에 적용되었으며 뉴튼-�N슨법에 비해 훨씬 우수한 특성을 보여주었다. 하지만 지금까지의 TLM법은 2차원의 정식화만 이루어졌고 3차원에는 적용되지 못한 것이 사실이다. 본 논문에서는 3차원의 비선형 정자계 문제에 TLM법을 적용할 수 있는 수식을 최초로 제안하며 3차원 코어(core)모델에 대해 TLM법을 적용하여 그 타당성을 검증하기로 한다. 또한 3차원 비선형 TLM법을 이용한 해석 결과가 뉴튼-�N슨법에 의한 결과와 완전히 일치하며 수렴 속도에 있어서도 훨씬 향상된 결과를 나타냄을 보이도록 하겠다.