• Journal of electromagnetic engineering and science
• /
• v.19 no.1
• /
• pp.6-12
• /
• 2019

An IE-FFT algorithm is implemented and applied to the electromagnetic (EM) solution of perfect electric conducting (PEC) scattering problems. The solution of the method of moments (MoM), based on the magnetic field integral equation (MFIE), is obtained for PEC objects with closed surfaces. The IE-FFT algorithm uses a uniform Cartesian grid to apply a global fast Fourier transform (FFT), which leads to significantly reduce memory requirement and speed up CPU with an iterative solver. The IE-FFT algorithm utilizes two discretizations, one for the unknown induced surface current on the planar triangular patches of 3D arbitrary geometries and the other on a uniform Cartesian grid for interpolating the free-space Green's function. The uniform interpolation of the Green's functions allows for a global FFT for far-field interaction terms, and the near-field interaction terms should be adequately corrected. A 3D block-Toeplitz structure for the Lagrangian interpolation of the Green's function is proposed. The MFIE formulation with the IE-FFT algorithm, without the help of a preconditioner, is converged in certain iterations with a generalized minimal residual (GMRES) method. The complexity of the IE-FFT is found to be approximately $O(N^{1.5})$and $O(N^{1.5}logN)$ for memory requirements and CPU time, respectively.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.143-150
• /
• 2000

A physical lumped parameter model is proposed for the time domain analysis of dam-reservoir system. The exact solution of transmitting boundary is derived for a semi-infinite 2-D reservoir of constant depth. The characteristics of the solution are examined in both frequency and the domains. Mass and damping coefficient are obtained from asymptotic behavior of the frequency domain solution. Further refinement to the lumped model is made by approximating the kernel function of the convolution integral in the exact solution. Finally a new physical lumped parameter model is proposed that consists of two masses, a spring and two dampers for each mode. It is demonstrated that new lumped parameter model of transmitting boundary can give excellent results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.8 s.251
• /
• pp.949-956
• /
• 2006

The method of Greenwood and Williamson is extended to obtain a solution to the coupled non-linear problem of steady-state electrical and thermal conduction across a crack in a conductive layer, for which the electrical resistivity and thermal conductivity are functions of temperature. The problem can be decomposed into the solution of a pair of non-linear algebraic equations involving boundary values and material properties. The new mixed-boundary value problem given from the thermal and electrical boundary conditions for the crack in the conductive layer is reduced in order to solve a singular integral equation of the first kind, the solution of which can be expressed in terms of the product of a series of the Chebyshev polynomials and their weight function. The non-existence of the solution for an infinite conductor in electrical and thermal conduction is shown. Numerical results are given showing the temperature field around the crack.

• 제어로봇시스템학회:학술대회논문집
• /
• 1990.10b
• /
• pp.1145-1149
• /
• 1990

This paper proposes a numerical method for redundant manipulators using predetermined optimal resolution. In order to obtain optimal joint trajectories, it is desirable to formulate redundancy resolution as an optimization problem having an integral cost criterion. We predetermine the trajectories of redundant joints in terms of the Nth partial sum of the Fourier series, which lead to the solution in the desirable homotopy class. Then optimal coefficients of the Fourier series, which yield the optimal solution within the predetermined class, are searched by the Powell's method. The proposed method is applied to a 3-link planar manipulator for cyclic tasks in Cartesian space. As the results, we can obtain the optimal solution in the desirable homotopy class without topological liftings of the solution. To show the validity of the proposed method, we analyze both optimal and extremal solutions by the Fast Fourier Transform (FFT) and discuss joint trajectories on the phase plane.

• Structural Engineering and Mechanics
• /
• v.38 no.5
• /
• pp.593-618
• /
• 2011

This paper deals with the problem of the determination of the response of a viscoelastic Bernoulli-Navier beam, which is resting on an elastic medium. Assuming uniaxial bending, the displacement of the beam axis is governed by an integro-differential equation. The compatibility of the displacements between the beam and the elastic medium is imposed through an integral equation. In general and in particular in the case of a Boussinesq medium, the solution has to be pursued numerically. On the contrary, in the case of a Winkler's medium the compatibility equation becomes a linear finite relationship, which allows finding an original analytical solution of the problem for both hereditary and aging behavior of the beam. Some numerical examples complete the paper, in which a comparison is made between the hereditary and the aging model for the creep of the beam.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.12
• /
• pp.945-953
• /
• 2008

Laminar film condensation of saturated vapor in forced flow over a flat plate is analysed. The problem is formulated as exact boundary-layer solution and integral approximate solution. From numerical solutions of the governing equations, it is found that the energy transfer by convection and the effect of inertia term in the momentum equation in negligibly small for low pressure but quite important for high pressure. The condensate rate, liquid-vapor interfacial shear stress and local heat transfer are strongly dependent on the reduced pressure $P_r$ and the modified Jacob number Ja/Pr.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1326-1338
• /
• 2017

The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of Consortium for Advanced Simulation of Light Water Reactors (CASL) tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated good agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle hot full power equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics/subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of $0.99994{\pm}6.8E-6$ (95% confidence interval). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.7
• /
• pp.759-768
• /
• 2012

This paper proposes a method to determine the elastic follow-up factors for the $C(t)$-integral under secondary stress. The rate of creep crack growth for transient creep is correlated with the $C(t)$-integral. Elastic follow-up behavior, which occurs in structures under secondary loading, prevents a relaxation of stress during transient creep. Thus, both the values of $C(t)$ and creep crack growth increase as increasing elastic follow-up. An estimation solution for $C(t)$ was proposed by Ainsworth and Dean based on the reference stress method. To predict the value of $C(t)$ using this solution, an independent method to determine the elastic follow-up factors for cracked bodies is needed. This paper proposed that the elastic follow-up factors for $C(t)$ can be determined by elastic-plastic analyses using the plastic-creep analogy. Finite element analyses were performed to verify this method.

• Steel and Composite Structures
• /
• v.34 no.4
• /
• pp.511-524
• /
• 2020

In this research, a simple four-variable trigonometric integral shear deformation model is proposed for the static behavior of advanced functionally graded (AFG) ceramic-metal plates supported by a two-parameter elastic foundation and subjected to a nonlinear hygro-thermo-mechanical load. The elastic properties, including both the thermal expansion and moisture coefficients of the plate, are also supposed to be varied within thickness direction by following a power law distribution in terms of volume fractions of the components of the material. The interest of the current theory is seen in its kinematics that use only four independent unknowns, while first-order plate theory and other higher-order plate theories require at least five unknowns. The "in-plane displacement field" of the proposed theory utilizes cosine functions in terms of thickness coordinates to calculate out-of-plane shear deformations. The vertical displacement includes flexural and shear components. The elastic foundation is introduced in mathematical modeling as a two-parameter Winkler-Pasternak foundation. The virtual displacement principle is applied to obtain the basic equations and a Navier solution technique is used to determine an analytical solution. The numerical results predicted by the proposed formulation are compared with results already published in the literature to demonstrate the accuracy and efficiency of the proposed theory. The influences of "moisture concentration", temperature, stiffness of foundation, shear deformation, geometric ratios and volume fraction variation on the mechanical behavior of AFG plates are examined and discussed in detail.

• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.