• Structural Engineering and Mechanics
• /
• v.43 no.3
• /
• pp.321-336
• /
• 2012

As an improvement on the isoparametric element method, the derivation presented in this paper is close to that done by Wang (1990) for the 2-D finite element. We extend this idea to solve 3-D crack problems in this paper. A new displacement modelling is constructed with local solutions of three-dimensional cracks and a quasi-compatible isoparametric element for three-dimensional fracture mechanics analysis is presented. The stress intensity factors can be solved directly by means of the present method without any post-processing. A new method for calculating the stress intensity factors of three-dimensional cracks with complex geometries and loads is obtained. Numerical examples are given to demonstrate the validity of the present method. The accuracy of the results obtained by the proposed element is demonstrated by solving several crack problems. The results illustrate that this method not only saves much calculating time but also increases the accuracy of solutions. Because this quasi-compatible finite element of 3-D cracks contains any singularities and easily meets the requirement of compatibility, it can be easily implemented and incorporated into existing finite element codes.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.12
• /
• pp.2019-2027
• /
• 2003

Pitting wear is a dominant from of polyethylene surface damage in total knee replacements, and may originate from surface cracks that propagate under repeated tribological contact. In this study, stress intensity factors, K$\_$I/and $_{4}$, were calculated for a surface crack in a polyethylene-CoCr-bone system under the rolling and/or sliding contact pressures. Crack length and load location were considered in determination of probable crack propagation mechanisms and fracture modes. Positive K$\_$I/ values were obtained for shorter cracks in rolling contact and for all crack lengths when the sliding load was apart from the crack. $_{4}$ was the greatest when the load was directly adjacent to the crack (g/a=${\pm}$1). Sliding friction caused a substantial increase of both K$\_$I/$\^$max/ and $_{4}$$\^$max/. The effective Mode I stress intensity factors, K$\_$eff/, were the greatest at g/a=${\pm}$1, showing the significance of high shear stresses generated by loads adjacent to surface cracks. Such behavior of K$\_$eff/ suggests mechanisms for surface pitting by which surface cracks may propagate along their original plane under repeated rolling or sliding contact.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.10
• /
• pp.2398-2406
• /
• 1993

Abdi's numerical method(ref.13) for representing a stress singularity by shifting the mid-side nodes of isoparametric elements is reviewed. A simple technique to obtain the optimal position of the mid-side nodes in quadratic isoparametric finite element is presented. From this technique we can directly obtain the position of the side-nodes adjacent to the crack tip. It is also observed that the present technique provides good accuracy for the expression of the opening displacement and the determination of the mid-side nodes for more wide range of material properties than that obtained by Abdicant the finite element method is applied to determine stress intensity factors for pressurized crack perpendicular to and terminating at the interface of two bonded dissimilar materials. A proper definition for stress intensity factors of a crack perpendicular to bimaterial interface is provided. It is based upon a near-tip displacement solutions on the crack surface for interface crack between two dissimilar materials. Numerical testing is carried out with the eight-node and six-node elements. The results obtained are compared with the previous solutions.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.8
• /
• pp.147-155
• /
• 1995

Recent experimental studies have been shown that strain gages can be employed to determine either static or dynamic stress intensity factors $K_{I}$ wiht relatively simple experiments. However, it does not usually provide a reliable value of stress intensity factor because of local yielding and limited regions for strain gage placement at the vicinity of the crack tip. This paper attempted to define a valid region and to indicate procedures for locating and orienting the strain gage to determine static toughness $K_{Is}$ accurately form one strain gage readings with respect to varying loadings. The strain gage methods was used for compact tension specimens made of Polycarbonate and PMMA(polymethyl methacrylate). Series expansions of the static and dynamic strain fields are applied. Strain gage orientation and location are then studied to optimize the strain response. Especially, in the dynamic experiment, the specimen employed is an oversized Charpy V-notch specimen which has been modified to provide significant constraint with a large elevation of the flow stress. The impact behavior of the specimen is monitored by placing strain gage near the crack tip. The dynamic toughness $K_{Id}$ is determined from the strain time traces of this gage.e.

• Journal of Welding and Joining
• /
• v.30 no.6
• /
• pp.80-85
• /
• 2012

Fatigue cracks in structural components are the most common cause of structural failure when exposed to fatigue loading. In this respect, fatigue crack detection and structural health assessment are very important. Currently, various smart materials are used for detecting fatigue crack and measurement of SIFs(Stress Intensity Factors). So, this paper presented a measurement of SIFs using MFC(Micro Fiber Composite) sensor which is the one of the smart material. MFC sensor is more flexible, durable and reliable than other smart materials. The SIFs of Mode I(K I) as well as Mode II(K II) based on the piezoelectric constitutive law and fracture mechanics are calculated. In this study, the SIF values measured by MFC sensors are compared with the theoretical results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.1
• /
• pp.36-43
• /
• 1987

Stress intensity factors are computed by the boundary element method employing the multiregion technique along with the double-point concept. To demonstrate the validity of the current method, the stress intensity factors of the well-known simple models such as a slanted edge crack and an arcular crack are determined, in advanced, which are proved to be in good agreement within 5% with the pre-existing solutions. Z-shaped cracks are analyzed with various branch crack lengths and branching angles.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.8
• /
• pp.1875-1882
• /
• 1993

Variable thickness plates are commonly used as structural members in the majority of industrial sectors. Previous fracture mechanics researches on variable thickness plates were limited to mode I loading cases. In practice, however, cracks are usually located inclined to the loading direction. In this respect, combined mode(mode I/II) stress intensity factors $K_{I}$ and $K_{II}$ at the crack tip for a variable thickness plate were obtained by 3-dimensional finite element analysis. Variable thickness plates containing a slant edge crack were chosen. The parameters used in this study were dimensionless crack $length{\lambda}$, slant $angle{\alpha}$, thickness $ratio{\beta}$ and width ratio{\omega}$. Stress intensity factors were calculated by crack opening displacement(COD) and crack sliding displacement(CSD)method proposed by Ingraffea and Manu.

• Journal of the Korean Society of Safety
• /
• v.14 no.4
• /
• pp.28-36
• /
• 1999

In this paper, A circular arc crackered plate in biaxially stretched sheets was investigated in the boundary element method. The applications of fracture mechanics have traditionally concentrated on crack problems under an mode I, straight crack. However, many service failures occur from growth of cracks subjected to mixed mode loadings. A rectangular plate with arc crack or slanted central crack, under biaxial tensile loading, was treated analytically and also solved numerically. The Results from BEM applying different loading conditions, crack length (a/W), arc angle($\alpha$) are presented and discussed. The stress intensity factors are evaluated by the techniques of the J-integral. The decomposition method, used to decouple the stress intensity factors in mixed mode problems, is implemented by a considering a small circular contour path around each crack tip. The BIE method was successfully applied to a circular arc crackerd plate problem, also slanted centre cracked plate under mixed mode.

• Structural Engineering and Mechanics
• /
• v.58 no.2
• /
• pp.217-230
• /
• 2016

Effects of temperature dependent material properties on mixed mode fracture parameters of functionally graded materials subjected to thermal loading are investigated. A domain form of the $J_k$-integral method including temperature-dependent material properties and its numerical implementation using finite element analysis is presented. Temperature and displacement fields are calculated using finite element analysis and are used to compute mixed mode stress intensity factors using the $J_k$-integral. Numerical results indicate that temperature-dependency of material properties has considerable effect on the mixed-mode stress intensity factors of cracked functionally graded structures.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.276-286
• /
• 1999

A recently developed numerical method based on a volume integral formulation is developed for the effective accurate calculation of the stress intensity factors at the crack tips in unbounded isotropic solids in the presence of multiple anisotropic inclusions and cracks and subjected to external loads. In this paper, a detailed analysis of the stress intensity factors are carried out for an unbounded isotropic matrix containing an orthotropic cylindrical inclusion and a crack. The accuracy and effectiveness of the new method are examined through comparison with results obtained from analytical method and finite element method using ANSYS. It is demonstrated that this new method is very accurate and effective for solving plane elastostatic problems in unbounded solids containing anisotropic inclusions and cracks.