• Coupled systems mechanics
• /
• 제7권6호
• /
• pp.669-690
• /
• 2018

This paper deals with numerical modeling of dynamic failure phenomena in rate-sensitive brittle and/or ductile materials. To this end, a two-dimensional continuum viscodamage-embedded discontinuity model, which is based on our previous work (see Do et al. 2017), is developed. More specifically, the pre-peak nonlinear and rate-sensitive hardening response of the material behavior, representing the fracture-process zone creation, is described by a rate-dependent continuum damage model. Meanwhile, an embedded displacement discontinuity model is used to formulate the post-peak response, involving the macro-crack creation accompanied by exponential softening. The numerical implementation in the context of the finite element method exploiting the second-order mid-point scheme is discussed in detail. In order to show the performance of the model several numerical examples are included.

• Journal of applied mathematics & informatics
• /
• 제36권3_4호
• /
• pp.181-203
• /
• 2018

In this paper, we have constructed an overlapping Schwarz method for singularly perturbed second order convection-diffusion equations. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region we use the central finite difference scheme on a uniform mesh while on the non-layer region we use the mid-point difference scheme on a uniform mesh. It is shown that the numerical approximations which converge in the maximum norm to the exact solution. When appropriate subdomains are used, the numerical approximations generated from the method are shown to be first order convergent. Furthermore it is shown that, two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantages of this method used with the proposed scheme is it reduces iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.

• 대한기계학회논문집A
• /
• 제29권2호
• /
• pp.253-261
• /
• 2005

In decades, a substantial body of work on a unified viscoplastic model which considers the mechanism of plastic deformation and creep deformation has developed. The systematic scheme for numerical analysis of unified model is necessary because the dominant failure mechanism is the defect growth and coalescence in materials. In the present study, the unified viscoplastic model for materials with defects suggested by Suquet and Michel was employed for numerical analysis. The constitutive equations are integrated based on the generalized mid-point rule and implemented into a finite element program (ABAQUS) by means of user-defined subroutine (UMAT). To evaluate the validity of the developed UMAT code and the assessment of the adopted viscoplastic model, the results obtained from the UMAT code was compared with the numerical reference solution and experimental data. The unit cell analysis also has been investigated to study the effect of strain rate, temperature, stress triaxiality and initial defect volume fraction on the growth and coalescence of the defect.

• Journal of applied mathematics & informatics
• /
• 제27권5_6호
• /
• pp.1001-1015
• /
• 2009

In this paper, two hybrid difference schemes on the Shishkin mesh are constructed for solving a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative. We prove that the schemes are almost second order convergence in the supremum norm independent of the diffusion parameter. Error bounds for the numerical solution and its derivative are established. Numerical results are provided to illustrate the theoretical results.

• 한국물환경학회지
• /
• 제29권4호
• /
• pp.514-522
• /
• 2013

The major 24 tributaries in Nam-River mid-watershed were monitored for discharge and water quality in order to understand the characteristics of the watershed and to select the tributary catchment for improving water quality. According to the analytical results of discharge and water quality monitoring data of 24 tributaries, the mean value of discharge below $0.1m^3/s$ was 62.5% among the monitored tributaries and it mostly exceeded the water quality standards of Nam-river mid-watershed ($BOD_5$ = 3 mg/L, T-P = 0.1 mg/L over). According to the stream grouping method and the water quality delivery load density ($kg/day/km^2$) based on the results of tributary discharge and water quality monitoring, the tributary watersheds for improving the water quality were selected. In the Nam-River mid-watershed, tributaries in the GaJwaCheon, HaChonCheon catchment (Group D, $BOD_5$ = 3 mg/L over) and in the UirYeongCheon, SeokGyoCheon catchment (Group A, T-P = 0.1 mg/L over), which have a small flow (and/or large flow) and a high concentrations of water pollutants. The various water quality improving scheme for tributaries, in accordance with the reduction of potential point source pollution by living sewage and livestock wastewater, should be established and implemented.

• 지구물리와물리탐사
• /
• 제3권2호
• /
• pp.61-69
• /
• 2000

지표 경계 조건은 유한 차분법을 이용한 탄성파 파동 방정식 모델링에서 수치해의 정확성을 떨어뜨리는 한편 포아송 비에 따른 해의 안정성을 제한하는 주 요인이 된다. 본 연구에서는 지표 경계 조건과 같은 Neumann 경계 조건의 처리에 효과적인 적분법(integration method)에 기반 하여 차분식을 유도하고, 이로부터 중간점 차분 기법을 제안하였다. 또한, 개발된 알고리즘을 Lamb의 문제에 적용하여 이론해와 비교함으로써 타당성을 검증하였다.

• Wind and Structures
• /
• 제9권4호
• /
• pp.315-330
• /
• 2006

Flow and scalar dispersion around Cheomseongdae are numerically investigated using a three-dimensional computational fluid dynamics (CFD) model with the renormalization group (RNG) $k-{\varepsilon}$ turbulence closure scheme. Cheomseongdae is an ancient astronomical observatory in Gyeongju, Korea, and is chosen as a model obstacle because of its unique shape, that is, a cylinder-shaped architectural structure with its radius varying with height. An interesting feature found is a mid-height saddle point behind Cheomseongdae. Different obstacle shapes and corresponding flow convergences help to explain the presence of the saddle point. The predicted size of recirculation zone formed behind Cheomseongdae increases with increasing ambient wind speed and decreases with increasing ambient turbulence intensity. The relative roles of inertial and eddy forces in producing cavity flow zones around an obstacle are conceptually presented. An increase in inertial force promotes flow separation. Consequently, cavity flow zones around the obstacle expand and flow reattachment occurs farther downwind. An increase in eddy force weakens flow separation by mixing momentum there. This results in the contraction of cavity flow zones and flow reattachment occurs less far downwind. An increase in ambient wind speed lowers predicted scalar concentration. An increase in ambient turbulence intensity lowers predicted maximum scalar concentration and acts to distribute scalars evenly.

• Steel and Composite Structures
• /
• 제9권2호
• /
• pp.131-158
• /
• 2009

This paper presents a modelling technique for the nonlinear analysis of composite steel-concrete beams with partial shear interaction. It extends the applicability of two stiffness elements previously derived by the authors using the direct stiffness method, i.e. the 6DOF and the 8DOF elements, to account for material nonlinearities. The freedoms are the vertical displacement, the rotation and the slip at both ends for the 6DOF stiffness element, as well as the axial displacement at the level of the reference axis for the 8DOF stiffness element. The solution iterative scheme is based on the secant method, with the convergence criteria relying on the ratios of the Euclidean norms of both forces and displacements. The advantage of the approach is that the displacement and force fields of the stiffness elements are extremely rich as they correspond to those required by the analytical solution of the elastic partial interaction problem, thereby producing a robust numerical technique. Experimental results available in the literature are used to validate the finite element proposed in the paper. For this purpose, those reported by Chapman and Balakrishnan (1964), Fabbrocino et al. (1998, 1999) and Ansourian (1981) are utilised; these consist of six simply supported beams with a point load applied at mid-span inducing positive bending moment in the beams, three simply supported beams with a point load applied at mid-span inducing negative bending moment in the beams, and six two-span continuous composite beams respectively. Based on these comparisons, a preferred degree of discretisation suitable for the proposed modelling technique expressed as a function of the ratio between the element length and depth is proposed, as is the number of Gauss stations needed. This allows for accurate prediction of the nonlinear response of composite beams.

• 한국항공우주학회지
• /
• 제49권1호
• /
• pp.21-30
• /
• 2021

탐색기를 탑재한 유도탄의 중기궤적은 탐색기 시야(FOV : Field-Of-View) 내에서 표적을 탐지하며, 전환 시점에서의 기동성을 최대화하도록 설계하는 것이 요구된다. 유도탄의 비행궤적 최적화 문제는 여러 구속조건이 적용된 비선형 문제로 일반적인 해석해를 도출하기 어렵기 때문에 그 동안 다양한 계산적인 방법들이 제시되어 왔다. 본 논문에서는 추진-활공 유도탄의 중기궤적 최적화 문제를 컨벡스 최적화 기법인 2차 원뿔 프로그래밍을 이용하여 산출하였다. 먼저, 운동방정식의 상태변수를 최소화하기 위해서 제어변수 구속조건이 추가된 제어변수 추가 형태의 운동방정식을 구성하였다. 또한, 자유 종말시간 문제와 추진시간 문제를 대처하기 위하여 정규화된 시간 변수를 독립 변수로 설정하였다. 그리고, 운동방정식과 제어변수 구속조건을 컨벡스 형태로 변환하기 위하여 각각 부분 선형화와 무손실 컨벡스 변환을 적용하였다. 마지막으로, 본 논문에서 제시된 방안의 타당성을 확인하기 위하여 비선형 최적화 프로그래밍 결과와 비교하였다.

• 한국수산과학회지
• /
• 제43권2호
• /
• pp.159-168
• /
• 2010

This paper presents a mathematical model and simulation method for investigating the performance of set net systems and fish cage systems influenced by currents and waves. Both systems consist of netting, mooring ropes, a floating collar and sinkers. The netting and ropes were considered flexible structures and the floating collar was considered an elastic structure. Both were modeled on a mass-spring model. The structures were divided into finite elements and mass points were placed at the mid-point of each element, and the mass points were connected by mass-less springs. Each mass point was subjected to external and internal forces and the total force was calculated at every integration step. An implicit integration scheme was used to solve the nonlinear dynamic system. The computation method was applied to dynamic simulation of actual systems simultaneously influenced by currents and waves in order to evaluate their practicality. The simulation results improved our understanding of the behavior of the structure and provided valuable information concerning the optimized design of set net and fish cage systems exposed to an open ocean environment.