• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.1
• /
• pp.13-23
• /
• 2008

In [7, 8], they proposed a new singular function(NSF) method to compute singular solutions of Poisson equations on a polygonal domain with re-entrant angles. Singularities are eliminated and only the regular part of the solution that is in $H^2$ is computed. The stress intensity factor and the solution can be computed as a post processing step. This method was extended to the interface problem and Poisson equations with the mixed boundary condition. In this paper, we give NSF method for the Helmholtz equations ${\Delta}u+Ku=f$ with homogeneous Dirichlet boundary condition. Examples with a singular point are given with numerical results.

• Structural Engineering and Mechanics
• /
• v.15 no.5
• /
• pp.563-578
• /
• 2003

Delaunay triangulation is combined with an adaptive finite element method for analysis of two-dimensional crack propagation problems. The content includes detailed descriptions of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around the crack tips and large elements in the other regions. Three examples for predicting the stress intensity factors of a center cracked plate, a compact tension specimen, a single edge cracked plate under mixed-mode loading, and an example for simulating crack growth behavior in a single edge cracked plate with holes, are used to evaluate the effectiveness of the procedure. These examples demonstrate that the proposed procedure can improve solution accuracy as well as reduce total number of unknowns and computational time.

• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.165-170
• /
• 2000

The weight function method is an efficient technique to calculate the stress intensity factors for various loading conditions in that only the stress analysis of an uncracked model is required. This paper analyzes the mixed-mode stress intensity factors of surface and deepest points for quarter elliptical surface cracks in mechanical joints by weight function method and the coefficients included in weight function are determined by finite element analyses for reference loadings. Results for the different number of terms in weight function are presented.

• Smart Structures and Systems
• /
• v.16 no.3
• /
• pp.401-414
• /
• 2015

A developed hybrid method for crack identification of beams is presented. Based on the Euler-Bernouli beam theory and concepts of fracture mechanics, governing equation of the cracked beams is reformulated. Finite element (FE) method as a powerful numerical tool is used to discritize the equation in space domain. After transferring the equations from time domain to frequency domain, frequencies and mode shapes of the beam are obtained. Efficiency of the governed equation for free vibration analysis of the beams is shown by comparing the results with those available in literature and via ANSYS software. The used equation yields to move the influence of cracks from the stiffness matrix to the mass matrix. For crack identification measured data are produced by applying random error to the calculated frequencies and mode shapes. An objective function is prepared as root mean square error between measured and calculated data. To minimize the function, hybrid genetic algorithms (GAs) and particle swarm optimization (PSO) technique is introduced. Efficiency, Robustness, applicability and usefulness of the mixed optimization numerical tool in conjunction with the finite element method for identification of cracks locations and depths are shown via solving different examples.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.9 s.180
• /
• pp.2344-2352
• /
• 2000

Mechanical joints such as bolted or riveted joints are widely used in mechanical components. The reliable determination of the stress intensity factors for cracks in bolted joints is needed to evaluate the safety and fatigue life of them. The weight function method is an efficient technique to calculate the stress intensity factors for various loading conditions because only the stress analysis of an uncracked model is required. In this paper the mixed-mode stress intensity factors for cracks in bolted joints are obtained by weight function method, in which the coefficients of weight function are determined by finite element analyses for reference loadings. Critical inclined angle that mode I stress intensity factor becomes maximum is determined and the effects of crack length and the magnitude of clearance on critical inclined angle are investigated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.322-329
• /
• 2001

A study of the shock wave turbulent boundary layer interaction is presented. The focus of the study is the interactions of the shock waves with the turbulent boundary layer on the falt plate. Three examples are investigated. The computations are performed, using mixed explicit-implicit generalized Galerkin finite element method. The linear equations at each time step are solved by a preconditioned GMRES algorithm. Numerical results indicate that the implicit scheme converges to the asymptotic steady state much faster than the explicit counterpart. The computed surface pressures and skin friction coefficients display good agreement with experimental data. The flowfield manifests a complex shock wave system and a pair of counter-rotating vortices.

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.621-639
• /
• 2011

The purpose of this paper is to introduce a general iterative process by viscosity approximation method with parallel method to ap-proximate a common element of the set of solutions of a mixed equilibrium problem and of the set of common fixed points of a finite family of $k_i$-strict pseudo-contractions in a Hilbert space. We obtain a strong convergence theorem of the proposed iterative method for a finite family of $k_i$-strict pseudo-contractions to the unique solution of variational inequality which is the optimality condition for a minimization problem under some mild conditions imposed on parameters. The results obtained in this paper improve and extend the corresponding results announced by Liu (2009), Plubtieng-Panpaeng (2007), Takahashi-Takahashi (2007), Peng et al. (2009) and some well-known results in the literature.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.395-405
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among four filling patterns at each triangular control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. Using the proposed numerical technique, the collapse of a water dam has been analyzed to predict flow phenomenon of fluid and the predicted front positions with respect to time have been compared with the reported experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.11
• /
• pp.1348-1358
• /
• 2004

The filling pattern and an adaptive grid refinement based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible transient viscous flow with free surfaces. The governing equation fur flow analysis is Navier-Stokes equation including inertia and gravity effects. The mixed FE formulation and predictor-corrector method are used effectively for unsteady numerical simulation. The flow front surface and the volume inflow rate are calculated using the filling pattern technique to select an adequate pattern among seven filling patterns at each tetrahedral control volume. By adaptive grid refinement, the new flow field that renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. In this domain the elements in the surface region are made finer than those in the remaining regions for more efficient computation. The collapse of a water dam and the filling of a fluidity spiral have been analyzed. The numerical results have been in good agreement with the experimental results and the efficiency of the adaptive grid refinement and filling pattern techniques have been verified.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.703-714
• /
• 2010

This paper proposes an approach for staged finite element modeling with coupled seepage and stress analysis. The stage modeling is based on the predefined inter-relationship between the base model and the unit stage models. A unit stage constitutes a complete finite element model, of which the geometries and attributes are subject to changes from stage to stage. The seepage analysis precedes the mechanical stress analysis at every stage. Division of the wet and dry zone and the pore pressures are evaluated from the seepage analysis and used in determining input data for the stress analysis. The results of the stress analysis may also be associated with the pore water pressures. For consolidation analysis, the pore pressure and the displacement variables are mixed in a coupled matrix equation. The time marching solution produces the dissipation of excess pore pressure and variation of stresses with passage of time. For undrained analysis, the excess pore pressures are computed from the stress increment due to loading applied in the unit stage and are used in revising the hydraulic head. The solution results of a unit stage are inherited and accumulated to the subsequent stages through the relationship of the base model and the individual unit stages. Implementation of the proposed approach is outlined on the basis of the core procedures, and numerical examples are presented for demonstration of its application.