• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.5
• /
• pp.141-153
• /
• 1999

A finite element procedure for the analysis of rubber-like hyperelastic material is developed. The volumetric incompressiblity conditions of the rubber deformation is included in the formulation by using penalty method. In this paper, the behavior of the rubber deformation is represented by hyperelastic constitutive relations based on a generalized Mooney-Rivlin model. The principle of virtual work is used to derive nonlinear finite element equation for the large displacement problem and presented in total-Lagrangian description. The finite element procedure using analytic differentiation resulted in very close solution to the result of the well known commercial packages NISAII AND ABAQUS. Numerical tests show that the results from the numerical differentiation method coincide very well with those from the analytic method and the well known commercial packages in static analysis. The convergency of rubber usingν iteration method is also discussed.

• Nuclear Engineering and Technology
• /
• v.27 no.3
• /
• pp.374-384
• /
• 1995

The mathematical adjoint solution of the Analytic Function Expansion (AFEN) method is found by solving the transposed matrix equation of AFEN nodal equation with only minor modification to the forward solution code AFEN. The perturbation calculations are then performed to estimate the change of reactivity by using the mathematical adjoint The adjoint calculational scheme in this study does not require the knowledge of the physical adjoint or the eigenvalue of the forward equation. Using the adjoint solutions, the exact and first-order perturbation calculations are peformed for the well-known benchmark problems (i.e., IAEA-2D benchmark problem and EPRI-9R benchmark problem). The results show that the mathematical adjoint flux calculated in the code is the correct adjoint solution of the AFEN method.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2230-2241
• /
• 2006

This paper presents a sole application of boundary element method to the conjugate heat transfer problem of thermally developing laminar flow in a thick walled pipe when the fluid velocities are fully developed. Due to the coupled mechanism of heat conduction in the solid region and heat convection in the fluid region, two separate solutions in the solid and fluid regions are sought to match the solid-fluid interface continuity condition. In this method, the dual reciprocity boundary element method (DRBEM) with the axial direction marching scheme is used to solve the heat convection problem and the conventional boundary element method (BEM) of axisymmetric model is applied to solve the heat conduction problem. An iterative and numerically stable BEM solution algorithm is presented, which uses the coupled interface conditions explicitly instead of uncoupled conditions. Both the local convective heat transfer coefficient at solid-fluid interface and the local mean fluid temperature are initially guessed and updated as the unknown interface thermal conditions in the iterative solution procedure. Two examples imposing uniform temperature and heat flux boundary conditions are tested in thermally developing region and compared with analytic solutions where available. The benchmark test results are shown to be in good agreement with the analytic solutions for both examples with different boundary conditions.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.11
• /
• pp.2773-2781
• /
• 1993

Laplace's equation is, perhaps, the most important equation, which governs various kinds of physical phenomena. Due to its importance, there have been several numerical techniques such as the finite element method, the finite difference method, and the boundary element method. However, these techniques do not appear very effective as they require a substantial amount of numerical calculation. In this paper, we develop a new most efficient technique based on analytic solutions for Laplace's equation in some convex polygons. Although a similar approach was used for the same problem, the present technique is unique as it solves directly Laplace's equation with the utilization of analytical solutions.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.4
• /
• pp.289-294
• /
• 2018

In this work, a recently developed semi-analytic technique, so called the residual power series method, is generalized to process higher-dimensional linear and nonlinear partial differential equations. The solutions obtained takes a form of an infinite power series which can, in turn, be expressed in a closed exact form. The results reveal that the proposed generalization is very effective, convenient and simple. This is achieved by handling the (m+1)-dimensional Burgers equation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.941-949
• /
• 1995

An efficient method based on analytic solutions is applied to solve anti-plane Mode III crack problems. The analytic technique developed earlier by the present authors for Laplace's equation in a simply-connected region is now extended to general Mode III crack problems. Unlike typical numerical methods which require fine meshing near crack tips, the present method divides the cracked bodies, typically non-convex or multiply-connected, into only a few super elements. In each super element, an element stiffness matrix, relating the series coefficients of the traction and displacement, is first formed. Then an assembly algorithm similar to that used in the finite elements, is first formed. Then an assembly algorithm similar to that used in the finite elements, is developed. A big advantage of the present method is that only the boundary conditions are to be satisfied in the solution procedure due to the use of analytic solutions. Several numerical results demonstrate the efficiency and accuracy of the present method.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.6
• /
• pp.461-469
• /
• 2017

A segment of the wave board has been expressed as a submerged line segment in the two dimensional wave flume. The lower end of the line segment could be extended to the bottom of the wave flume and the other opposite upper end of the board could be extended to the free surface. It is assumed that the motion of the wave board could be defined by the sinusoidal motion in horizontal direction on either end of the wave board. When the amplitude of sinusoidal motion of the wave board on lower and upper end are equal, the wave board motion could express the horizontally oscillating submerged segment of piston type wave generator. The submerged segment of flap type wave generator also could be expressed by taking the motion amplitude differently for the either end of the board. The pivot point of the segment motion could play a role of hinge point of the flap type wave generator. Simplified analytic solution of oscillating submerged wave board segment in water of finite depth has been derived through the first order perturbation method at two dimensional domain. The case study of the analytic solution has been carried out and it is found out that the solution could be utilized for the design of wave generator with arbitrary shape by linear superposition.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.10
• /
• pp.4706-4712
• /
• 2013

One of the efficient procedures for the solution of partial differential equations is the method of differential quadrature. This method has been applied to a large number of cases to circumvent the difficulties of programming complex algorithms for the computer, as well as excessive use of storage due to conditions of complex geometry and loading. Under in-plane uniform distributed load, the buckling of asymmetric curved beam with varying cross section is analyzed by using differential quadrature method (DQM). Critical load due to diverse cross section variation and opening angle is calculated. Analysis result of DQM is compared with the result of exact analytic solution. As DQM is used with small grid points, exact analysis result is shown. New result according to diverse cross section variation is also suggested.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.17-28
• /
• 1999

A new procedure of the boundary element method(BEM),say, singular BEM for the potential problems with singularities is presented. To obtain the numerical solution of which asymptotic behavior near the singularities is close to that of the analytic solution, we use particular elements on the boundary segments containing singularities. The Motz problem and the crack problem are taken as the typical examples, and numerical results of these cases show the efficiency of the present method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.6
• /
• pp.36-44
• /
• 2003

In this paper, an analytic solution method is proposed to overcome the numerical problems when the slewing dynamics of hybrid coordinate systems is investigated via time finite element analysis. It is shown that the dynamics of the hybrid coordinate systems is governed by the coupled dual differential equations for both slewing and structural modes. Structural modes are transformed into the time-based modal coordinates and analytic spatial propagation equations are derived for each space-dependent time mode. Slew angle history is obtained analytically by appropriate applications of the boundary conditions and structural propagation is re-calculated using the slew angle. Numerical examples are demonstrated to validate the proposed analytic method in comparison to the existing state transition matrix method.