• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.131-138
• /
• 1997

This paper concerns the stress singularity at the interface corner of the viscoelastic adhesive layer bonded between rigid adherends, subjected to a uniform transverse tensile strain. The characteristic equation is derived in the Laplace transformed space, following Williams, and the transformed characteristic equation is inverted analytically into real time space for the viscoelastic model considered here. The order of the singularity is obtained numerically. The time-domain boundary element method is employed to investigate the nature of stresses along the interface. Numerical results show that the order of the singularity increases with time while the free-edge stress intensity factors are relaxed with time.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.1-8
• /
• 2004

For the linearized differential algebraic equation of the nonlinear constrained system, exact initial values of the acceleration are needed to solve itself. It may be very troublesome to perform the inverse operation for obtaining the incremental quantities since the mass matrix contains the zero element in the diagonal. This fact makes the mass matrix impossible to be positive definite. To overcome this singularity phenomenon the mass matrix needs to be modified to allow the feasible application of predictor and corrector in the iterative computation. In this paper the proposed numerical algorithm based on the modified mass matrix combines the conventional implicit algorithm, Newton-Raphson method and Newmark method. The numerical example presents reliabilities for the proposed algorithm via comparisons of the 4th order Runge-kutta method. The proposed algorithm seems to be satisfactory even though the acceleration, Lagrange multiplier, and energy show unstable behaviour. Correspondingly, it provides one important clue to another algorithm for the enhancement of the numerical results.

• Communications for Statistical Applications and Methods
• /
• v.24 no.3
• /
• pp.303-315
• /
• 2017

Yoo (2015, Statistics and Probability Letters, 99, 109-113) derives theoretical results in an optimal sufficient dimension reduction with singular inner-product matrix. The results are promising, but Yoo (2015) only presents one simulation study. So, an evaluation of its practical usefulness is necessary based on numerical studies. This paper studies the asymptotic behaviors of Yoo (2015) through various simulation models and presents a real data example that focuses on ordinary least squares. Intensive numerical studies show that the $x^2$ test by Yoo (2015) outperforms the existing optimal sufficient dimension reduction method. The basis estimation by the former can be theoretically sub-optimal; however, there are no notable differences from that by the latter. This investigation confirms the practical usefulness of Yoo (2015).

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.2
• /
• pp.101-113
• /
• 2018

The dual singular function method(DSFM) is a numerical algorithm to get optimal solution including corner singularities for Poisson and Helmholtz equations. In this paper, we apply DSFM to solve heat equation which is a time dependent problem. Since the DSFM for heat equation is based on DSFM for Helmholtz equation, it also need to use Sherman-Morrison formula. This formula requires linear solver n + 1 times for elliptic problems on a domain including n reentrant corners. However, the DSFM for heat equation needs to pay only linear solver once per each time iteration to standard numerical method and perform optimal numerical accuracy for corner singularity problems. Because the Sherman-Morrison formula is rather complicated to apply computation, we introduce a simplified formula by reanalyzing the Sherman-Morrison method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.4
• /
• pp.281-292
• /
• 2009

The solution of the interface problem or Poisson problem with concave corner has singular perturbation at the interface corners or singular corners. The decoupled dual singular function method (DDSFM) which exploits the singular representations of the solutions was suggested in [3, 9] and estimated optimal accuracy in [10]. The convergence rates consist with theoretical results even for the problems with very strong singularity, with the efficiency depending on parameters used in the methods. Furthermore the errors in $L^2$ and $L^\infty$-spaces display some oscillation, in the cases with meshsize not small enough. In this paper, we present an answer to remove the oscillation via numerical experiments. We observe the effects of parameters in DDSFM, and show the consisting efficiency of the method over the strong singularity.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.2
• /
• pp.73-78
• /
• 1992

Numerical method by Abdi et al. for obtaining a stress singularity near a crack perpendicular to the interface between two elastic materials is reviewed. More efficient and simple method to obtain crack singularity by shifting a mid-side node of 8-node isoparametric element is presented. It is observed that the present analysis provides increased accuracy for the expression of the opening displacement and the determination of the optimal position of the mid-side node for a wide range of material properties.

• Journal of Institute of Control, Robotics and Systems
• /
• v.18 no.9
• /
• pp.823-829
• /
• 2012

An excitation table is commonly used for vibration and ride tests for parts or assemblies of automobiles, aircrafts, or other heavy systems. The authors have analyzed several kinematic properties of an excitation table that is under development for heavy transport vehicles. It consists of one table and 7 linear hydraulic actuators. The authors have performed mobility analysis, inverse kinematics, forward kinematics, and singularity analysis. Especially, we have proposed a fast forward kinematic solution considering the limited motion of the excitation table. On the assumption that the motion variables such as rotation angles and displacements are small, the forward kinematic problem is converted to the observer problem of a linear system. This provides a fast solution. Also we have verified that there are no singularity points in the working range by numerical analysis.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.1
• /
• pp.14-28
• /
• 1992

For a planar curve-sided paned with constant or linear density distributions of source or doublet in the singularity methods, Cantaloube and Rehbach(1986) show that the surface integral can be transformed into contour integral by using Stokes' formulas. As an extension of their formulations, this paper deals with a planar polygonal panel for which we derive the closed-forms of the potentials and the velocities induced by the singularity distributions. Test calculations show that the analytical evaluation of the closed-forms is superior to numerical integration(suggested by Cantaloube and Rehbach) of the contour integral. The compact and explicit expressions may produce accurate values of matrix elements of simultaneous linear equations in the singularity methods with much reduced computer tiome.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.4
• /
• pp.351-365
• /
• 2019

This article presents a new approach for conformal flattening with optimal cone singularity. The algorithm here takes an optimal control for selecting optimal cones and uses the Ricci flow to force the flattening. This work is considered as a modification to the work of Soliman et al. [1] in the sense that they make use of the Yamabe equation for the flattening, which is an approximation of the Ricci flow. We present a numerical algorithm based on the optimal control with the mathematical background. Several numerical results validate that our method is optimal in total cone angle and usage of the Ricci flow ensures the conformal flattening while selecting optimal cones.

• Honam Mathematical Journal
• /
• v.40 no.4
• /
• pp.785-794
• /
• 2018

In [8] they introduced a new finite element method for accurate numerical solutions of Poisson equations with corner singularities. They consider the Poisson equations with homogeneous Dirichlet boundary condition with one corner singularity at the origin, and compute the finite element solution using standard FEM and use the extraction formula to compute the stress intensity factor, then pose a PDE with a regular solution by imposing the nonhomogeneous boundary condition using the computed stress intensity factor, which converges with optimal speed. From the solution they could get an accurate solution just by adding the singular part. This approach uses the polar coordinate and the cut-off function to control the singularity and the boundary condition. In this paper we consider Poisson equations with multiple singular points, which involves different cut-off functions which might overlaps together and shows the way of cording in FreeFEM++ to control the singular functions and cut-off functions with numerical experiments.