• Structural Engineering and Mechanics
• /
• v.28 no.4
• /
• pp.461-477
• /
• 2008

An adaptive finite element method for analyzing two-dimensional and axisymmetric nonlinear elastic fracture mechanics problems with cracks is presented. The J-integral is used as a parameter to characterize the severity of stresses and deformation near crack tips. The domain integral technique, for which all relevant quantities are integrated over any arbitrary element areas around the crack tips, is utilized as the J-integral solution scheme with 9-node degenerated crack tip elements. The solution accuracy is further improved by incorporating an error estimation procedure onto a remeshing algorithm with a solution mapping scheme to resume the analysis at a particular load level after the adaptive remeshing technique has been applied. Several benchmark problems are analyzed to evaluate the efficiency of the combined domain integral technique and the adaptive finite element method.

• Kyungpook Mathematical Journal
• /
• v.60 no.4
• /
• pp.869-881
• /
• 2020

Most first kind integral equations are ill-posed, and obtaining their numerical solution often requires solving a linear system of algebraic equations of large condition number, which may be difficult or impossible. This article proposes a regularization-direct method to numerically solve first kind Fredholm integral equations. The vector forms of block-pulse functions and related properties are applied to formulate the direct method and reduce the integral equation to a linear system of algebraic equations. We include a regularization scheme to overcome the ill-posedness of integral equation and obtain a stable numerical solution. Some test problems are solved using the proposed regularization-direct method to illustrate its efficiency for solving first kind Fredholm integral equations.

• Structural Engineering and Mechanics
• /
• v.34 no.1
• /
• pp.85-95
• /
• 2010

In this paper, numerical solution of the singular integral equation for the multiple curved branch-cracks is investigated. If some quadrature rule is used, one difficult point in the problem is to balance the number of unknowns and equations in the solution. This difficult point was overcome by taking the following steps: (a) to place a point dislocation at the intersecting point of branches, (b) to use the curve length method to covert the integral on the curve to an integral on the real axis, (c) to use the semi-open quadrature rule in the integration. After taking these steps, the number of the unknowns is equal to the number of the resulting algebraic equations. This is a particular advantage of the suggested method. In addition, accurate results for the stress intensity factors (SIFs) at crack tips have been found in a numerical example. Finally, several numerical examples are given to illustrate the efficiency of the method presented.

• Journal of Ship and Ocean Technology
• /
• v.3 no.4
• /
• pp.1-14
• /
• 1999

The improved Green integral equation using the Kelvin-type Green function in known free of irregular frequencies where the integral over the inner free surface integral is removed from the integral equation, resulting in an overdetermined integral equation. The solution of the overdetermined Green integral equation is shown identical with the solution of the improved Green integral equation Using the B-spline higher order panel method, the overdetermined equation is discretized in two different ways; one of the resulting linear system is square and the other is redundant. Numerical experiments show that the solutions of both are identical. Using the present methods, the exact values and higher derivatives of the potential at any place over the wetted surface of the body can be found with much fewer panels than low order panel method.

• Structural Engineering and Mechanics
• /
• v.72 no.4
• /
• pp.491-502
• /
• 2019

This paper explores the analytical buckling solution of rectangular thin plates by the finite integral transform method. Although several analytical and numerical developments have been made, a benchmark analytical solution is still very few due to the mathematical complexity of solving high order partial differential equations. In solution procedure, the governing high order partial differential equation with specified boundary conditions is converted into a system of linear algebraic equations and the analytical solution is obtained classically. The primary advantage of the present method is its simplicity and generality and does not need to pre-determine the deflection function which makes the solving procedure much reasonable. Another advantage of the method is that the analytical solutions obtained converge rapidly due to utilization of the sum functions. The application of the method is extensive and can also handle moderately thick and thick elastic plates as well as bending and vibration problems. The present results are validated by extensive numerical comparison with the FEA using (ABAQUS) software and the existing analytical solutions which show satisfactory agreement.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.52 no.4
• /
• pp.185-191
• /
• 2003

In this Paper, we propose an accurate and stable solution of the transient electromagnetic response from three-dimensional arbitrarily shaped conducting objects by using a time domain magnetic field integral equation. This method does not utilize the conventional marching-on in time (MOT) solution. Instead we solve the time domain integral equation by expressing the transient behavior of the induced current in terms of temporal expansion functions with decaying exponential functions and Laguerre·polynomials. Since these temporal expansion functions converge to zero as time progresses, the transient response of the induced current does not have a late time oscillation and converges to zero unconditionally. To show the validity of the proposed method, we solve a time domain magnetic field integral equation for three closed conducting objects and compare the results of Mie solution and the inverse discrete Fourier transform (IDFT) of the solution obtained in the frequency domain.

• The Pure and Applied Mathematics
• /
• v.14 no.4
• /
• pp.279-282
• /
• 2007

Using the center instead of the Lipschitz condition we show how to provide weaker sufficient convergence conditions of the modified Newton Kantorovich method for the solution of nonlinear singular integral equations with Curleman shift (NLSIES). Finer error bounds on the distances involved and a more precise information on the location of the solution are also obtained and under the same computational cost than in [1].

• East Asian mathematical journal
• /
• v.23 no.2
• /
• pp.257-260
• /
• 2007

Using the center instead of the Lipschitz condition we show how to provide weaker sufficient convergence conditions of the modified Newton Kantorovich method for the solution of nonlinear singular integral equations with Curleman shift (NLSIES). Finer error bounds on the distances involved and a more precise information on the location of the solution are also obtained and under the same computational cost than in [1].

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.8
• /
• pp.2524-2539
• /
• 1996

A design process and an axisymmetric boundary integral equation method for precise structural analysis of the aircraft gas turbine rotor disk were developed. This axisymmetric boundary integral equation method for stress and steady-state thermal analysis was improved in solution accuracy by appling an implicit technique for Cauchy principal value evaluation, a subelement technique for weak singular integral evaluation and a double exponentical integral technoque for internal point solution near boundary surfaces. Stresses, temperatures, low cycle fatigue lifes and critical speeds for the turbine rotor disk of the thrust 1421 N class turbojet engine were analysed in a pratical calculation model problem.

• Journal of Ship and Ocean Technology
• /
• v.6 no.1
• /
• pp.34-53
• /
• 2002

Numerical solution of the forward-speed radiation problem for a half-immersed cylinder advancing in regular waves is presented by making use of the improved Green integral equation in the frequency domain. The B-spline higher order panel method is employed stance the potential and its derivative are unknown at the same time. The present numerical solution of the improved Green integral equation by the B-spline higher order Kelvin panel method is shown to be free of irregular frequencies which are present in the Green integral equation using the forward-speed Kelvin-type Green function.