• Korea-Australia Rheology Journal
• /
• v.18 no.4
• /
• pp.171-181
• /
• 2006

We present a short review for authors' previous work on direct numerical simulations for inertialess hard particle suspensions formulated either with a Newtonian fluid or with viscoelastic polymeric fluids to understand the microstructural evolution and the bulk material behavior. We employ two well-defined bi-periodic domain concepts such that a single cell problem with a small number of particles may represent a large number of repeated structures: one is the sliding bi-periodic frame for simple shear flow and the other is the extensional bi-periodic frame for planar elongational flow. For implicit treatment of hydrodynamic interaction between particle and fluid, we use the finite-element/fictitious-domain method similar to the distributed Lagrangian multiplier (DLM) method together with the rigid ring description. The bi-periodic boundary conditions can be effectively incorportated as constraint equations and implemented by Lagrangian multipliers. The bulk stress can be evaluated by simple boundary integrals of stresslets on the particle boundary in such formulations. Some 2-D example results are presented to show effects of the solid fraction and the particle configuration on the shear and elongational viscosity along with the micro-structural evolution for both particles and fluid. Effects of the fluid elasticity has been also presented.

• Journal of the korean Society of Automotive Engineers
• /
• v.12 no.1
• /
• pp.26-32
• /
• 1990

This study investigated a method for the determination of the J-integral for a tough glassy polymer such as polycarbonate plates by using the method of caustics. Comparing the values of J-integral determined by a numerical analysis and by the method of caustics, the method of caustics was found to be an effective experimental technique for the determination of the J-integral. The ratio between two J-integrals determined by the method of caustics and by finite element method converged into 1 within the limit of low load. However, it was noticed that the greater the plastic zone at the crack tip was, the lower the J-integral obtained by the reflect method of caustics. This difference may be deduced from the damage at the crack tip such as craze appeared in the polycarbonate plate. It was confirmed that the ratio of longitudinal diameter( $D_{l}$ ) to transverse diameter ( $D_{t}$) of caustics generally converged into 1 at the low load. The transition of the state of stress at the vicinity of a crack tip from plane strain to plane stress was deduced by noticing that the longitudinal diameter( $D_{l}$ ) grew faster than the transverse diameter( $D_{t}$) of caustics within the higher load range.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.1
• /
• pp.13-26
• /
• 2006

When the stagnation temperature of a perfect gas increases, the specific heat for constant pressure and ratio of the specefic heats do not remain constant any more and start to vary with this temperature. The gas remains perfect: its state equation remains always valid, with exception that it will be named by calorically imperfect gas. The aim of this research is to develop the relations of the necessary thermodynamics and geometrical ratios. and to study the supersonic flow at high temperature. lower than the threshold of dissociation. The results are found by the resolution of nonlinear algebraic equations and integration of complex analytical functions where the exact calculation is impossible. The dichotomy method is used to solve the nonlinear equation. and the Simpson algorithm for the numerical integration of the found integrals. A condensation of the nodes is used. Since. the functions to be integrated have a high gradient at the extremity of the interval of integration. The comparison is made with the calorifcally perfect gas to determine the error made by this last. The application is made for the air in a supersonic nozzle.

• Wind and Structures
• /
• v.7 no.6
• /
• pp.421-447
• /
• 2004

This paper presents a time domain approach for predicting buffeting response of long suspension bridges under skew winds. The buffeting forces on an oblique strip of the bridge deck in the mean wind direction are derived in terms of aerodynamic coefficients measured under skew winds and equivalent fluctuating wind velocities with aerodynamic impulse functions included. The time histories of equivalent fluctuating wind velocities and then buffeting forces along the bridge deck are simulated using the spectral representation method based on the Gaussian distribution assumption. The self-excited forces on an oblique strip of the bridge deck are represented by the convolution integrals involving aerodynamic impulse functions and structural motions. The aerodynamic impulse functions of self-excited forces are derived from experimentally measured flutter derivatives under skew winds using rational function approximations. The governing equation of motion of a long suspension bridge under skew winds is established using the finite element method and solved using the Newmark numerical method. The proposed time domain approach is finally applied to the Tsing Ma suspension bridge in Hong Kong. The computed buffeting responses of the bridge under skew winds during Typhoon Sam are compared with those obtained from the frequency domain approach and the field measurement. The comparisons are found satisfactory for the bridge response in the main span.

• Journal of the Korean Society for Railway
• /
• v.18 no.3
• /
• pp.194-202
• /
• 2015

In this paper, dynamic modeling with viscoelastic properties of a human body resting on a seat is presented to quantitatively analyze ride quality of passengers exposed to vertical vibrations. In describing the motions of a seated body, a 5 degree-of-freedom multibody model from the literature is investigated. The viscoelastic characteristics of seats used in railway vehicles are mathematically formulated with nonlinear stiffness characteristics and convolution integrals representing time delay terms. Transfer functions for the floor input are investigated and it is found that these are different in accordance with the input magnitude due to nonlinear characteristics of the seat. Measured floor input at the railway vehicle is used to analyze realistic human vibration characteristics. Frequency weighted RMS acceleration values are calculated and the effects of the seat design parameters on the frequency weighted RMS acceleration values are presented.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.487-492
• /
• 2007

The enhancement of the service life of damaged or cracked structures is a major issue for researchers and engineers. The hierarchical void element with the integrals of Legend polynomials is used to characterize the fracture behavior of unpatched crack as well as repaired crack with bonded composite patches by computing the stress intensity factors and stress contours at the crack tip. The numerical approach is based on the v-version degenerate shell element including the theory of anisotropic laminated composites. Since the equivalent single layer approach is adopted in this study, the proposed element is necessary to represent a discontinuous crack part as a continuum body with zero stiffness of materials. Thus the aspect ratio of this element to represent the crack should be extremely slender. The sensitivity of numerical solution with respect to energy release rate, displacement and stress has been tested to show the robustness of hierarchical void element as the aspect ratio is increased up to 2000. The stiffness derivative method and displacement extrapolation method have been applied to calculate the stress intensity factors of Mode I problem.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.481-486
• /
• 2007

The Zienkiewicz-Zhu(Z/Z) error estimate is slightly modified for the hierarchical p-refinement, and is then applied to L-shaped plates subjected to bending to demonstrate its effectiveness. An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the superconvergent patch recovery(SPR) technique. The modified Z/Z error estimate p-refinement is different from the conventional approach because the high order shape functions based on integrals of Legendre polynomials are used to interpolate displacements within an element, on the other hand, the same order of basis function based on Pascal's triangle tree is also used to interpolate recovered stresses. The least-square method is used to fit a polynomial to the stresses computed at the sampling points. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly or selectively. It is noted that the error decreases rapidly with an increase in the number of degrees of freedom and the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.375-382
• /
• 2008

The particular integral formulation for two(2D) and three(3D) dimensional inelastic transient dynamic stress analysis is presented. The elastostatic equation is used for the complementary solution. Using the concept of global shape function, the particular integrals for displacement and traction rates are obtained to approximate acceleration of the inhomogeneous equation. The Houbolt time integration scheme is used for the time-marching process. The Newton-Raphson algorithm for plastic multiplier is used to solve the system equation. Numerical results of four example problems are given to demonstrate the validity and accuracy of the present formulation.

• Journal of the Institute of Convergence Signal Processing
• /
• v.5 no.4
• /
• pp.317-324
• /
• 2004

Electromagnetically coupled coaxial dipole (ECCD) array antenna with and without short-ended termination is investigated theoretically. The integral equations are derived for the structure by use of the Fourier transform and mode expansion of radial waveguide. The integrals appearing in the integral equations are evaluated along the branch cut instead of real axis for a faster convergent integral. The effects of slots and dipoles, short-ended termination length, and dielectric coating on the radiation characteristics are presented. Radiation pattern of the structure is also investigated. The results of the present method are compared with those of the commercial EM simulator and good agreement is found.

• Geophysics and Geophysical Exploration
• /
• v.3 no.1
• /
• pp.13-18
• /
• 2000

We describe the concept of the singularity in the integral equation of electromagnetic (EM) modeling in comparison with that in the integral representation of electric fields in EM theory, which would clarify the singular integral problems of the Green's tensor. We have also derived and classified the singular integrals of the Green's tensors in 3-D, 2.5-D and 2-D as well as in the thin sheet integral equations of the EM scattering problem, which have the most important effect on the accuracy of the numerical solution of the problems.