• 대한기계학회논문집A
• /
• 제26권5호
• /
• pp.904-913
• /
• 2002

In this paper, the dynamic stress intensity factors of fracture mechanics are numerically computed in time domain using the FEM. For which the finite element formulations are derived applying the Galerkin method to the equations of motion in convolution integral as has been presented in the previous paper. To assure the strain fields of r$^{-1}$ 2/ singularity near the crack tip, the triangular quarter-point singular elements are imbedded in the finite element mesh discretized by the isoparametric quadratic quadrilateral elements. Two-dimensional problems of the elastodynamic fracture mechanics under the impact load are solved and compared with the existing numerical and analytical solutions, being shown that numerical results of good accuracy are obtained by the presented method.

• 대한기계학회논문집A
• /
• 제28권10호
• /
• pp.1583-1589
• /
• 2004

In this paper, A new finite element method for the time domain analysis of the dynamic stress intensity factor of two-dimensional viscoelastic body with a stationary central crack under the transient dynamic load is presented, which is based on the intergrodifferential equations of motion in the isotropic linear viscoelasticity and the Galerkin's method. The vlscoelastic material is assumed to be elastic in dilatation and behaves like a standard linear solid in shear. As a numerical example, the Chen's problem in viscoelastodynamic version is solved for the parametric study about the effect of viscosity and relaxation time on the dynamic stress intensity factor.

• Journal of applied mathematics & informatics
• /
• 제32권5_6호
• /
• pp.585-598
• /
• 2014

In this paper we consider the nonlinear parabolic problems with mixed boundary condition. Under comparatively mild conditions of the coefficients related to the problem, we construct the discontinuous Galerkin approximation of the solution to the nonlinear parabolic problem. We discretize spatial variables and construct the finite element spaces consisting of discontinuous piecewise polynomials of which the semidiscrete approximations are composed. We present the proof of the convergence of the semidiscrete approximations in $L^{\infty}(H^1)$ and $L^{\infty}(L^2)$ normed spaces.

• Structural Engineering and Mechanics
• /
• 제3권1호
• /
• pp.61-74
• /
• 1995

A linearly tapered, doubly symmetric thin-walled closed member, such as power-transmission towers and tourist towers, are often characterized by local variation in mass and/or rigidity, due to additional mass and rigidity. On the preliminary stage of design the closed-form solution is more effective than the finite element method. In order to propose approximate solutions, the discontinuous and local variation in mass and/or rigidity is treated continuously by means of a usable function proposed by Takabatake(1988, 1991, 1993). Thus, a simplified analytical method and approximate solutions for the free and forced transverse vibrations in linear elasticity are demonstrated in general by means of the Galerkin method. The solutions proposed here are examined from the results obtained using the Galerkin method and Wilson-${\theta}$ method and from the results obtained using NASTRAN.

• 한국전산유체공학회지
• /
• 제13권4호
• /
• pp.13-23
• /
• 2008

This paper is an extension of previous study[1] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite (serendipity) cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires four degrees-of-freedom at each element corners. Those degrees-of-freedom are the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational basis functions from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[2].

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.33-41
• /
• 2008

This paper is an extension of previous study[9] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite serendipity cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires at each element corners four degrees-of-freedom such as the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational bases from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[7].

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.33-41
• /
• 2008

This paper is an extension of previous study[9] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite serendipity cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires at each element corners four degrees-of-freedom such as the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational bases from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[7].

• 한국전산유체공학회지
• /
• 제2권1호
• /
• pp.66-72
• /
• 1997

Optimum nozzle design exploiting the method of characteristic(M.O.C) has been in application as an efficient design methodology targeting a less weighted and short expansion nozzle. This paper treats the optimum nozzle design and the analysis of the inviscid compressible flow inside. Based on traditional Rao's method, the optimum nozzle design is coded with minor modifications for the identification of the control surface across which the mass flux should be conserved. Internal flow field is simulated numerically by M.O.C and implicit/explicit Taylor-Galerkin finite element method(F.E.M) with the aid of adaptive remeshing to capture the shock wave, hence improve the accuracy. Designed and calculated flow fields due to the separate analyses show that the mass flux predicted by optimum nozzle design with M.O.C is not conserved across the control surface and the sonic line should be located upstream of the nozzle throat. Rao's optimum nozzle design methodology exaggerates the momentum thrust and tends to overemphasize the engine performance loss.

• 한국전산구조공학회논문집
• /
• 제16권2호
• /
• pp.207-216
• /
• 2003

본 논문에서는 균열을 지닌 축대칭 문제를 해석하기 위하여 시간적분형 운동방정식을 바탕으로 한 유한요소 해법을 제시한다. 유한요소메쉬는 8절점 등매개변수 사변형 요소와 균열선단에서의 1/4절점 삼각형 특이요소로 구성되며, 동적 응력확대계수는 균열면상의 1/4절점의 y방향 변위로부터 구한다. 제시된 해법의 정확성과 타당성을 검증하기 위하여 내부에 원환균열을 지닌 무한 탄성체가 균열면상에서 충격하중을 받을 때의 동적 응력확대계수를 계산하고 타 수치결과와 비교 검토하였다. 응용 예제로서 원환균열과 원주균열을 지닌 중실축과 중공축의 동적 응력확대계수를 균열의 길이와 축의 길이에 따른 영향을 자세히 조사하였다. 균열길이가 커지면 동적 응력확대계수가 커지고, 축의 길이가 길어지면 동적 응력확대계수 곡선의 폭도 함께 증가됨을 확인하였다. 그리고 균열의 위치가 안쪽에 포함될 경우보다는 바깥쪽에 포함될 때 더 큰 동적 응력확대계수가 발생됨을 밝힌다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 춘계 학술대회논문집
• /
• pp.51-59
• /
• 2000

The Interactions of shock wave with turbulent boundary layers in high-speed flows cause complex flowfields which result in increased adverse pressure gradients, skin friction and temperature. Accurate and reliable prediction of such phenomena is needed in designing high-speed propulsion systems. Such analyses of the complex flowfields require sophisticated numerical scheme that can resolve interactions between shock wave and boundary layers accurately. Therefore the purpose of the present. article is to introduce an accurate and efficient mixed explicit-implicit generalized Galerkin finite element method. To demonstrate the validity of the theory and numerical procedure, several benchmark cases are investigated.