• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.1
• /
• pp.10-17
• /
• 1985

A new method is suggested for the solution of plane elasticity problems. With use of the dislocation model in the crack problems, the basic scheme of this method is to find equilibrium Burgers vectors of dislocations which are distributed along the boundary of the first fundamental boundary value problems. The stress distribution in the region can be found by superposition of the contributions of each dislocation. The method is applied to three cases with known analytical solutions, and to a V-notched specimen under uniaxial tension. The numerical results are compared with other available solutions. This method is effective and simple in its use, compared with other numerical methods. The method also provides very accurate solutions in the region except near the boundary where the discretization error is significant. The extrapolation method is suggested for the stresses in the boundary region. Extensive application are also suggested for a general estimate of the computational efficiency of the method.

• Communications of the Korean Mathematical Society
• /
• v.20 no.3
• /
• pp.563-578
• /
• 2005

In [Z. Cai and B. C. Shin, SIAM J. Numer. Anal. 40 (2002), 307-318], we developed the discrete first-order system least squares method for the second-order elliptic boundary value problem by directly approximating $H(div){\cap}H(curl)-type$ space based on the Helmholtz decomposition. Under general assumptions, error estimates were established in the $L^2\;and\;H^1$ norms for the vector and scalar variables, respectively. Such error estimates are optimal with respect to the required regularity of the solution. In this paper, we study solution methods for solving the system of linear equations arising from the discretization of variational formulation which possesses discrete biharmonic term and focus on numerical results including the performances of multigrid preconditioners and the finite element accuracy.

• Journal of KSNVE
• /
• v.8 no.4
• /
• pp.632-642
• /
• 1998

A method based on finite element discretization is developed for optimizing the polarization profile of PVDF film to create the modal transducer for specific modes. Using this concept, one can design the modal transducer in two-dimensional structure having arbitrary geometry and boundary conditions. As a practical means for implementing this polarization profile without repoling the PVDF film the polarization profile is approximated by optimizing electrode patterns, lamination angles, and poling directions of the multi-layered PVDF transducer. This corresponds to the approximation of a continuous function using discrete values. The electrode pattern of each PVDF layer is optimized by deciding the electrode of each finite element to be used or not. Genetic algorithm, suitable for discrete problems, is used as an optimization scheme. For the optimization of each layers lamination angle, the continuous lamination angle is encoded into discrete value using binary 5 bit string. For the experimental demonstration, a modal sensor for first and second modes of cantilevered composite plate is designed using two layers of PVDF films. The actuator is designed based on the criterion of minimizing the system energy in the control modes under a given initial condition. Experimental results show that the signals from residual modes are successfully reduced using the optimized multi-layered PVDF sensor. Using discrete LQG control law, the modal peaks of first and second modes are reduced in the amount of 12 dB and 4 dB, resepctively.

• Journal of the Korean Mathematical Society
• /
• v.50 no.6
• /
• pp.1291-1310
• /
• 2013

This paper develops least-squares pseudo-spectral collocation methods for elliptic boundary value problems having interface conditions given by discontinuous coefficients and singular source term. From the discontinuities of coefficients and singular source term, we derive the interface conditions and then we impose such interface conditions to solution spaces. We define two types of discrete least-squares functionals summing discontinuous spectral norms of the residual equations over two sub-domains. In this paper, we show that the homogeneous least-squares functionals are equivalent to appropriate product norms and the proposed methods have the spectral convergence. Finally, we present some numerical results to provide evidences for analysis and spectral convergence of the proposed methods.

• Smart Structures and Systems
• /
• v.2 no.3
• /
• pp.209-224
• /
• 2006

The use of smart dampers to optimally control the response of structures is on the increase. To maximize the potential use of such damper systems, their accurate modeling and assessment of their performance is of vital interest. In this study, the performance of a controllable fluid dashpot damper, in terms of damper forces, damper dynamic range and damping force hysteretic loops, respectively, is studied mathematically. The study employs a damper Bingham-Maxwell (BingMax) model whose mathematical formulation is developed using a Fourier series technique. The technique treats this one-dimensional Navier-Stokes's momentum equation as a linear superposition of initial-boundary value problems (IBVPs): boundary conditions, viscous term, constant Direct Current (DC) induced fluid plug and fluid inertial term. To hold the formulation applicable, the DC current level to the damper is supplied as discrete constants. The formulation and subsequent simulation are validated with experimental results of a commercially available magneto rheological (MR) dashpot damper (Lord model No's RD-1005-3) subjected to a sinusoidal stroke motion using a 'SCHENK' material testing machine in the Materials Laboratory at the University of Technology, Sydney.

• Korea-Australia Rheology Journal
• /
• v.15 no.3
• /
• pp.131-150
• /
• 2003

A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.