• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.04a
• /
• pp.71-78
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity, the associated flow rule, and von-Mises yield criteria. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the center cracked plate under tensile loading. Those results are compared with the there cal solutions and the numerical solutions of ADINA software.

• Journal of the Korean Society for Precision Engineering
• /
• v.18 no.12
• /
• pp.88-94
• /
• 2001

The V-notched crack problem in dissimilar materials can be formulated as an eigenvalue problem. The RWCIM(Reciprocal Work Contour Integral Method) is applied to the determination of the eigenvector coefficients associated with eigenvalues for V-notched cracks in pseudo-isotropic and anisotropic dissimilar materials. The RWCIM algorithm is programed by the commercial numerical program, MATHEMATICA. The numerical results obtained are shown that the RWCIM is a useful method for determining the eigenvector coefficients of V-notched cracks in pseudo-isotropic and anisotropic dissimilar materials.

• Magazine of the Korea Concrete Institute
• /
• v.7 no.3
• /
• pp.175-186
• /
• 1995

The determination of the stress intensity factors is investigated by sur-face integral method around the crack tip of the nlass~vc: concrete structure. The surface integral met hod is naturally derived from the standal-ci path integral J. Howevcr. In the J integral method, pressure in the crack-face and body forces can not be considered, while this theory has advantage of ccmsidering many kind of forces, so t.his theory will be useful in investigating more accurate strt:ss states around crack tip. Furthermore. t h~s rrlethod can elerninate unntussary process of using singular elements and fine mesh around crack tip which is used 11; modelling the singularity around crack tip. A computer program for determming $K_I$, $K_{II}$ is tfcvulopcd by applying this theory. $K_I$, $K_{II}$ values usmg X noded isoparametric elements which was proved and variation of the stress intensity factor was investigated by application of darn structures.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.161-171
• /
• 1993

The aim of this study is to provide the stress-strain behavior of hardening geological materials such as rock or concrete on three dimensional spaces by using Desai model based on plastic theory. To validate proposed model, truly triaxial tests with high pressure under variety of stress paths in which three principal stresses were controlled independently using concrete materials were performed. The main results are summerized as follows: 1. Various stress paths for hardening materials used are satisfactorily explained by performing the truly triaxial test with high pressure. This is very important to investigate constitutive equations for materials like rock or concrete. 2. Since the proposed yield function is continuous, it avoids the singularity point at the intersection of two function in the previous models, thus, reducing the difficulties for computer implementation. 3. Analytic predictions for yielding behavior on $J_1-{\sqrt{J_{2D}}}$ octahedral and triaxial plane, as well as volumetric strain and stress-strain behavior agree well with experimental results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.2
• /
• pp.145-152
• /
• 2009

SGBEM(Symmetric Galerkin Boundary Element Method)-FEM alternating method has been proposed by Nikishkov, Park and Atluri. In the proposed method, arbitrarily shaped three-dimensional crack problems can be solved by alternating between the crack solution in an infinite body and the finite element solution without a crack. In the previous study, the SGBEM-FEM alternating method was extended further in order to solve elastic-plastic crack problems and to obtain elastic-plastic stress fields. For the elastic-plastic analysis the algorithm developed by Nikishkov et al. is used after modification. In the algorithm, the initial stress method is used to obtain elastic-plastic stress and strain fields. In this paper, elastic-plastic J integrals for three-dimensional cracks are obtained using the method. For that purpose, accurate values of displacement gradients and stresses are necessary on an integration path. In order to improve the accuracy of stress near crack surfaces, coordinate transformation and partitioning of integration domain are used. The coordinate transformation produces a transformation Jacobian, which cancels the singularity of the integrand. Using the developed program, simple three-dimensional crack problems are solved and elastic and elastic-plastic J integrals are obtained. The obtained J integrals are compared with the values obtained using a handbook solution. It is noted that J integrals obtained from the alternating method are close to the values from the handbook.

• Composites Research
• /
• v.19 no.4
• /
• pp.15-22
• /
• 2006

One of the primary factors limiting the application of composite-metal adhesively bonded joints in structural design is the lack of a good evaluation tool for the interfacial strength to predict the load bearing capacity of boned joints. In this paper composite-steel adhesion strength is evaluated in terms of stress intensity factor and fracture toughness of the interface corner. The load bearing capacity of double lap joints, fabricated by co-cured bonding of composite-steel adherends has been determined using fracture mechanical analysis. Bi-material interface comer stress singularity and its order are presented. Finally stress intensities and fracture toughness of the wedge shape bi-material interface corner are determined. Double lap joint failure locus and its mixed mode crack propagation criterion on $K_1-K_{11}$ plane have been developed by tension tests with different bond lengths.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.311-319
• /
• 2007

A Meshfree is a method used to establish algebraic equations of system for the whole problem domain without the use of a predefined mesh for the domain discretization. A point interpolation method is based on combining radial and polynomial basis functions. Involvement of radial basis functions overcomes possible singularity Furthermore, the interpolation function passes through all scattered points in an influence domain and thus shape functions are of delta function property. This makes the implementation of essential boundary conditions much easier than the meshfree methods based on the moving least-squares approximation. This study aims to investigate a stress analysis of structural element between a meshfree method and the finite element method. Examples on cantilever type plate, hollow cylinder and stress concentration problems show that the accuracy and convergence rate of the meshfree methods are high.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.3
• /
• pp.183-190
• /
• 2019

This paper introduces a near-tip grid refinement and explores its usefulness in the crack analysis by the natural element method(NEM). As a sort of local h-refinement in finite element method(FEM), a NEM grid is locally refined around the crack tip showing high stress singularity. This local grid refinement is completed in two steps in which grid points are added and Delaunay triangles sharing the crack tip node are divided. A plane strain rectangular plate with symmetric edge cracks is simulated to validate the proposed local grid refinement and to examine its usefulness in the crack analysis. The crack analysis is also simulated using a uniform NEM grid for comparison. Unlike the uniform grid, the refined grid provides near-tip stress distributions similar to the analytic solutions and the fine grid. In addition, the refined grid shows higher convergence than the uniform grid, the global relative error to the total number of grid points.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.3
• /
• pp.381-390
• /
• 2001

In this paper, a new improved crack analysis technique by Element-Free Galerkin(EFG) method is proposed, in which the singularity and the discontinuity of the crack successfully described by adding enrichment terms containing a singular basis function to the standard EFG approximation and a discontinuity function implemented in constructing the shape function across the crack surface. The standard EFG method requires considerable addition of nodes or modification of the model. In addition, the proposed method significantly decreases the size of system of equation compared to the previous enriched EFG method by using localized enrichment region near the crack tip. Numerical example show the improvement and th effectiveness of the previous method.

• 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.