• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.486-491
• /
• 2000

In this study, finite elements addition and removal method by stress range is applied to optimize shapes in structures, without using classical and numerical optimization methods and search methods. The program based on this algorithm is developed and compared to theoritial results with considerable accuracy. Classical methods need mesh generation for finite element analysis for every iteration, the developed method needs updated mesh data such as coordinates of nodes, elements connectivity, and loads on nodes. And other tools of finite element analysis can be in use as a black box to interface with this program.

• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.177-187
• /
• 2009

Machine or structural members subjected to fatigue loading will have a crack initiated during early part of their life. Therefore analysis of members with cracks and other discontinuities is very important. Finite element method has enjoyed widespread use in engineering, but it is not convenient for crack problems as the region very close to crack tip is to be discretized with very fine mesh. However, as the body force method (BFM), requires only the boundary of the discontinuity (crack or hole) to be discretized it is easy versatile technique to analyze such problems. In the present work fundamental solution for concentrated load x + iy acting in the semi-infinite plate at an arbitrary point $z_0=x_0+iy_0$ is considered. These fundamental solutions are in complex form ${\phi}(z)$ and ${\psi}(z)$ (England 1971). These potentials are known as Melan potentials (Ramakrishna 1994). A crack in the semi-infinite plate as shown in Fig. 1 is considered. This crack is divided into number of divisions. By applying pair of body forces on a division, the resultant forces on the remaining 'N'divisions are to be found for which ${\phi}_1(z)$ and ${\psi}_1(z)$ are derived. Body force method is applied to calculate stress intensity factor for crack in semi-infinite plate. Also for the case of crack emanating from circular hole in semi-infinite plate radial stress, hoop stress and shear stress are calculated around the hole and crack. Convergent results are obtained by body force method. These results are compared with FEM results.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.216-232
• /
• 1991

This Finite Element Program(TAFEM) has been developed to be able to carry out the structural analsis of tunnel section and simulate the surrounding ground behaviour due to New Austrian Tunnelling Method, of which main support is the surrounding ground, itself. The Elasto-plastic theory has been applied. The used finite elements are 8-noded isoparametric element(rock & shotcrete), 2 or 3-noded rod element(rock bolt) and infinite boundary element. The load incremental method and tangential stiffness method has been used. Associated flow rule was applied to plastic flow and yield criteria inclued not only Mohr-Coulomb but also Drucker-Prager. In this paper, Drucker-Prager yield criterion has been used. The relationship between plastic strain and stress is based on the incremental strain concept and stress-strain equation on the basis of the stress path of each gauss point has been adopted. It may be rational that rock is considered to be no-tension material, so that no-tension analysis has been adopted in accordance with the brittle fracture constitutive equation.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.8 no.3
• /
• pp.66-75
• /
• 1999

This paper is concerned with a analysis of noise by inverse problem available for analyzing the two and three-dimensional acoustic field problems. The noise of analysis considered in this study can be reduced to an optimum problem to find the optimal set of parameters defining the vibrating state of noise source surfaces. The optimal set of parameters are searched by the standard optimization procedure minimizing the square sum of the residuals between the measured and computed quantities of sound pressure at some points in the acoustic field. Computation is carried out for typical examples in which the noise sources are located on the infinite plane. It is demonstrated that the noise of analysis can be effectively made by using the sensitive reference data.

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.225-239
• /
• 2005

Applications of the Trefftz boundary element method (BEM) to anti-plane electroelastic problems are presented in this paper. Both direct and indirect methods with domain decomposition are discussed in details. Each crack is treated as semi-infinite thin slit defined in a subregion, for which a particular solution of the anti-plane problem, satisfying exactly the crack-face condition, is derived. The stress intensity factors defined at each crack tip can be directly computed from the coefficients of the particular solution. The performance of the proposed formulation is assessed by two examples and comparison is made with results obtained by other approaches. The Trefftz boundary element approach is demonstrated to be suitable for the analysis of the anti-plane problem of piezoelectric materials.

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

This paper presents a time domain method for soil-structure interaction analysis for seismic loadings. It is based on the finite element formulation incorporating analytical frequency-dependent infinite elements for the far field soil. The dynamic stiffness matrices of the far field region formulated using the present method in frequency domain can be easily transformed into the corresponding matrices in time domain. At first, the equivalent earthquake forces are evaluated along the interface between the near and the far fields from the free-field response analysis carried out in frequency domain, and the results are transformed into the time domain. An efficient procedure is developed for the convolution integrals to evaluate the interaction force along the interface, which depends on the response on the interface at the past time instances as well as the concurrent instance. Then, the dynamic responses are obtained for the equivalent earthquake force and the interaction force using Newmark direct integration technique. Since the response analysis is carried out in time domain, it can be easily extended to the nonlinear analysis. Example analysis has been carried out to verify the present method in a multi-layered half-space.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 2000.10a
• /
• pp.434-439
• /
• 2000

In this paper, in order to study the geometric factor effect of a circular hole near a crack tip in a semi-infinite plate, the Dimensionless Stress Intensity Factor, $F(=\frac K {\sigma {\sqrt{\pi a}}})$ is analyzed at the crack tip using a two Dimensional Boundary Element Method (BEM) program which is known as superior in Fracture Mechanics. Kelvin's solution was used as a fundamental solution in BEM analysis and displacement extrapolation method was used to determine Stress Intensity Factor.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.234-239
• /
• 2001

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isometric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.89-98
• /
• 1993

As the construction of nuclear power plants are increased, nuclear wastes disposal has been faced as a serious problem. If nuclear wastes are to be buried in the underground stratum, thermo-mechanical behavior of stratum must be analyzed, because high temperature distribution has a significant effect on tunnel and surrounding stratum. In this study, in order to analyze the structural behavior of the underground which is subject to concentrated heat sources, a coupling method of nonlinear finite elements and linear boundary elements is proposed. The nonlinear finite elements (NFE) are applied in the vicinity of nuclear depository where thermo-mechanical stress is concentrated. The boundary elements are also used in infinite domain where linear behavior is expected. Using the similar method as for the problem in mechanical field, the coupled nonlinear finite element-boundary element (NFEBE) is developed. It is found that NFEBE method is more efficient than NFE which considers nonlinearity in the whole domain for the nuclear wastes depository that is expected to exhibit local nonlinearity behavior. The effect of coefficients of the rock mass such as Poisson's ratio, elastic modulus, thermal diffusivity and thermal expansion coefficient is investigated through the developed method. As a result, it is revealed that the displacements around tunnel are largely dependent on the thermal expansion coefficients.

• Proceedings of the KSR Conference
• /
• 2010.06a
• /
• pp.1243-1248
• /
• 2010

Fatigue design and evaluation of welded joints are typically carried out by weld classification approach in which a family (theoretically infinite) of parallel nominal stress based S-N curves are used according to joint types and loading modes as well as extrapolation-based hot spot stress. Traditional finite element methods are not capable of consistently capturing the stress concentration effects on fatigue behavior due to their mesh-sensitivity in stress determination at welds resulted from notch stress singularity. The extrapolated hot spot stresses tend vary, depending on the element sizes, types, joint types, and loading mode. however, the equilibrium-equivalent structural stress method(E2S2) has been recently developed through several joint industry projects as a robust method to analyze welded components using finite element analysis. This method has been proven effective in correlating a large amount of published fatigue test results in the literature such as master S-N curve and has used for evaluating the fatigue life of welding components. In this study, fatigue analysis of the welding bogie frame is examined using E2S2 method with master S-N curve.