• Computers and Concrete
• /
• v.25 no.6
• /
• pp.537-549
• /
• 2020

Geometrically nonlinear axisymmetric bending analysis of shear deformable circular plates on a nonlinear three-parameter elastic foundation was made. Plates ranging from "thin" to "moderately thick" were investigated for three types of material: isotropic, transversely isotropic, and orthotropic. The differential equations were discretized by means of the finite difference method (FDM) and the differential quadrature method (DQM). The Newton-Raphson method was applied to find the solution. A parametric investigation using seven unknowns per node was presented. The novelty of the paper is that detailed numerical simulations were made to highlight the combined effects of the material properties and the boundary conditions on (i) the deflection, (ii) the stress resultants, and (iii) the external load. The formulation was verified through comparison studies. It was observed that the results are highly influenced from the boundary conditions, and from the material properties.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.29-34
• /
• 2004

Numerical and experimental studies on pipeline laying for intake Deep Ocean Water are carried out. In the numerical study, an implicit finite difference algorithm is employed for three-dimensional pipe equations. Fluid non-linearity and bending stiffness are considered and solved by Newton-Raphson iteration. Seabed is modeled as elastic foundation with linear spring and damper. Top tension and general configuration of pipeline at a depth are predicted. It is found that control for tension to prevent being large curvature of pipeline is needed on th steep seabed and, it should be considered 23.5 ton of tension at a top of pipe on the process of pipeline laying at 400m of water depth The largest top tension of pipe on condition of the beam sea during pipe laying is shown from the experiment. The results of this study can be contributed to the design of pipeline laying for upwelling deep ocean water.

• Transactions on Electrical and Electronic Materials
• /
• v.18 no.2
• /
• pp.97-102
• /
• 2017

A current transformer (CT) is a type of sensor that consists of a combination of electric and magnetic circuits, and it measures large ac currents. When a large amount of current flows into the primary winding, the alternating magnetic flux in the iron core induces an electromotive force in the secondary winding. The characteristics of a CT are determined by the iron core design because the iron core is saturated above a certain magnetic flux density. In particular, when a large current, such as a current surge, is input into a CT, the iron core becomes saturated and the induced electromotive force in the secondary winding fluctuates severely. Under these conditions, the CT no longer functions as a sensor. In this study, the characteristics of the secondary winding were investigated using the time-difference finite element method when a current surge was provided as an input. The CT was modeled as a two-dimensional analysis object using constraints, and the saturation characteristics of the iron core were evaluated using the Newton-Rhapson method. The results of the calculation were compared with the experimental data. The results of this study will prove useful in the designs of the iron core and the windings of CTs.

• 제어로봇시스템학회:학술대회논문집
• /
• 2005.06a
• /
• pp.906-909
• /
• 2005

In this study, the dynamic flying characteristics of the worn head sliders are investigated theoretically due to the change in head geometry caused by head and disk contact. The film shapes can be approximated as taper- truncated cycloidal-flat film. Two-dimensional time dependent modified Reynolds equation included molecular slip effect are formulated with neglected the roughness effect. The motion of head slider was assumed to have two degree of freedom in this work. Finite difference approximation with Newton Raphson iterative technique and the fourth order Runge-Kutta method were implemented to obtain the transient response of the slider head with various change in head geometry numerically and compared with the transient response of the IBM3380 type head slider. The simulation results show the film shape has affects significantly on the static and dynamic characteristic of slider head in magnetic storage systems.

• KSTLE International Journal
• /
• v.10 no.1_2
• /
• pp.17-22
• /
• 2009

The numerical analysis of elastohydrodynamic lubrication for the ball and roller bearings is performed in order to study the effect of the number of rolling elements and the shaft load on the minimum film thickness. A finite difference method and the Newton-Raphson method are used in the analysis. For a given shaft load, the maximum load of rolling element is determined along with the number of rolling elements. And then the minimum film thickness is calculated for several rolling bearings. The shape of film thickness and the pressure distribution are also studied.

• Tribology and Lubricants
• /
• v.17 no.6
• /
• pp.435-440
• /
• 2001

A numerical analysis between the piston and cylinder in hydraulic piston pumps under reciprocating motion is presented. A finite difference method and the Newton-Raphson method are used simultaneously to solve the Reynolds equation in the clearance and the equation of motion for the piston. The tapered piston showed stable behaviors regardless of their initial eccentric positions in the clearance, and the reciprocating speed affect highly on the piston end trajectories. Therefore, the numerical methods and results of present study can be used in the lubrication study of other piston-cylinder type fluid machineries.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1999.11a
• /
• pp.129-135
• /
• 1999

A numerical analysis between the piston and cylinder in swash plate type hydraulic piston pumps under reciprocating motion is presented. A finite difference method and the Newton-Raphson method are used simultaneously to solve the Reynolds equation In the clearance and the equation of motion for the piston. The tapered piston showed stable behaviors regardless of their initial eccentric conditions in the clearance and the reciprocating speed affect highly on the piston end trajectories. Therefore, the results of present study can be used other types fluid machineries.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.433-451
• /
• 2015

This paper focuses on large deflection static behavior of edge cracked simple supported beams subjected to a non-follower transversal point load at the midpoint of the beam by using the total Lagrangian Timoshenko beam element approximation. The cross section of the beam is circular. The cracked beam is modeled as an assembly of two sub-beams connected through a massless elastic rotational spring. It is known that large deflection problems are geometrically nonlinear problems. The considered highly nonlinear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. There is no restriction on the magnitudes of deflections and rotations in contradistinction to von-Karman strain displacement relations of the beam. The beams considered in numerical examples are made of Aluminum. In the study, the effects of the location of crack and the depth of the crack on the non-linear static response of the beam are investigated in detail. The relationships between deflections, end rotational angles, end constraint forces, deflection configuration, Cauchy stresses of the edge-cracked beams and load rising are illustrated in detail in nonlinear case. Also, the difference between the geometrically linear and nonlinear analysis of edge-cracked beam is investigated in detail.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2002.10a
• /
• pp.191-196
• /
• 2002

The paper presents a numerical and experimental investigation on dynamic behaviors of a towed low tension cable. In the numerical study, an implicit finite difference algorithm is employed for three-dimensional cable equations. Fluid and geometric non-linearity and bending stiffness are considered and solved by Newton-Raphson iteration. Block tri-diagonal matrix method is applied for the fast calculation of the huge size of matrices. In order to verify the numerical results and to see real physical phenomena, an experiment is carried out for a 6m cable in a deep and long towing tank. The cable is towed in two different ways; one is towed at a constant speed and the other is towed at a constant speed with top end horizontal oscillations. Cable tension and shear forces are measured at the top end. Numerical and experimental results are compared with good agreements in most cases but with some differences in a few cases. The differences are due to drag coefficients caused by vortex shedding. In the numerical modeling, non-uniform element length needs to be employed to cope with the sharp variation of tension and shear forces at near top end.