• Korea-Australia Rheology Journal
• /
• 제13권1호
• /
• pp.37-45
• /
• 2001

Numerical simulation of coextrusion process of viscoelastic fluids within a die has been carried out. In the coextrusion process velocity profile at the outflow boundary is not known a priori, which makes it difficult to impose the proper boundary condition at the outflow boundary. This difficulty has been avoided by using the open boundary condition (OBC) method. In this study, elastic viscous stress splitting (EVSS) formulation with streamline upwind (SU) method has been used in the finite element method. In order to test the validity of the OBC method, comparison between the results of fully developed condition at the outlet and those of OBC has been made for a Newtonian fluid. In the case of upper convected Maxwell (UCM) fluid, the effect of outflow boundary condition on the interface position has been investigated by using two meshes having different downstream lengths. In both cases, the results with the OBC method showed reasonable interface shape. In particular, for the UCM fluid the interface shape calculated with OBC was independent of the downstream length, while the results with the zero traction condition showed oscillation of interface position close to the outlet. Viscosity difference was found to be more important than elasticity difference in determining the final interface position. However, the overshoot of interface position near the con-fluent point increased with elasticity.

• 한국기계가공학회지
• /
• 제2권4호
• /
• pp.61-66
• /
• 2003

An excellent mesh construction is of Importance in yielding good results of finite element analysis. The new mesh generation algorithm, which offsets boundaries inward, was developed on the basis of a looping method. An user interface technique and automatic splitting lines which both divide a given domain into subdomains manually or automatically, were used. In addition, the separation method has advantages to prevent the large scale of element size and to control numbers of nodes and elements. This new mesh generation algorithm was proved in practice.

• 콘크리트학회논문집
• /
• 제12권2호
• /
• pp.63-69
• /
• 2000

This study investigates a new method of using a concrete disk to calculate stress intensity factor (SIF) for mixed mode cases. The results indicate that the disk method is more accurate than three point bending test (TPB) in obtaining correct SIF values for mixed mode fracture propagation. Stress intensity factors $K_{I}$ and $K_{II}$ are calculated using a center notched disk subjected to splitting load. The notch angle is calculated by finite element (FEM). Fracture toughness $K_\textsc{k}$ of the concrete is obtained from the load intensities at the initiation of crack propagation. According to the finite element analysis(FEA) and disk test, the results show that mode I and mixed mode cracks propagate toward the directions of crack face and loading point, respectively. The results from FEA with maximum stress theory compare well with the experimental date. Unlike TPB method where an accurate fracture toughness value is difficult to obtain due to the irregular shape of load deflection curve and delayed final crack propagation (following slow stable cracking). fracture toughness value is easily measured in the disk test from the crack initial load. Therefore, it is safe to conclude that disk method is more advantageous than TPB method in analyzing combined mode fracture problems.

• Wind and Structures
• /
• 제20권3호
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• Wind and Structures
• /
• 제17권3호
• /
• pp.263-274
• /
• 2013

A numerical technique for fluid-structure interaction, which is based on the finite element method (FEM) and computational fluid dynamics (CFD), was developed for application to an industrial chimney equipped with a pendulum tuned mass damper (TMD). In order to solve the structural problem, a one-dimensional beam model (Navier-Bernoulli) was considered and, for the dynamical problem, the standard second-order Newmark method was used. Navier-Stokes equations for incompressible flow are solved in several horizontal planes to determine the pressure in the boundary of the corresponding cross-section of the chimney. Forces per unit length were obtained by integrating the pressure and are introduced in the structure using standard FEM interpolation techniques. For the fluid problem, a fractional step scheme based on a second order pressure splitting has been used. In each fluid plane, the displacements have been taken into account considering an Arbitrary Lagrangian Eulerian approach. The stabilization of convection and diffusion terms is achieved by means of quasi-static orthogonal subscales. For each period of time, the fluid problem was solved and the geometry of the mesh of each fluid plane is updated according to the structure displacements. Using this technique, along-wind and across-wind effects have been properly explained. The method was applied to an industrial chimney in three scenarios (with or without TMD and for different damping values) and for two wind speeds, showing different responses.

• 대한기계학회논문집B
• /
• 제27권8호
• /
• pp.1182-1190
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for large-scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady flow computation of the incompressible Navier-Stokes equation, and Smagorinsky LES model is chosen for turbulent flow computation. Both METIS and MPI Libraries are used for domain partitioning and data communication between processors, respectively. Tiburon model of Hyundai Motor Company is chosen as the computational model at Re=7.5 $\times$ 10$^{5}$ , which is based on the car height. The calculation is carried out under both the wind tunnel condition and the road condition using IBM SP parallel architecture at KISTI Super Computing Center. Compared with the existing experimental data, both the velocity and pressure fields are predicted reasonably well and the drag coefficient is in good agreement. Furthermore, it is confirmed that the drag under the road condition is smaller than that under the wind-tunnel condition.

• Nuclear Engineering and Technology
• /
• 제53권6호
• /
• pp.1846-1857
• /
• 2021

Ultrasonic nondestructive testing is important for monitoring the structural integrity of dissimilar metal welds (DMWs) in pressure vessels and piping in nuclear power plants. However, there is a low probability of crack detection via inspection of DMWs using ultrasonic waves because the grain structures (grain orientations) of the weld area cause distortion and splitting of ultrasonic beams propagating in anisotropic media. To overcome this issue, the grain orientation should be known, and a precise ultrasonic wave simulation technique in anisotropic media is required to model the distortion and splitting of the waves accurately. In this study, a method for nondestructive prediction of the DMW grain orientations is presented for accurate simulation of ultrasonic wave propagation behavior in the weld area. The ultrasonic wave propagation behavior in anisotropic media is simulated via finite-element analysis when ultrasonic waves propagate in a transversely isotropic material. In addition, a methodology to predict the DMW grain orientation is proposed that employs a simulation technique for ultrasonic wave propagation behavior calculation and an optimization technique. The simulated ultrasonic wave behaviors with the grain orientations predicted via the proposed method demonstrate its usefulness. Moreover, the method can be used to determine the focal law in DMWs.

• 터널과지하공간
• /
• 제29권6호
• /
• pp.523-531
• /
• 2019

암석역학이나 암반공학의 여러 활용 분야 중 발파, 천공 및 기계굴착의 경우에는 파괴인성이 중요한 요소이며, 암석의 파괴인성을 측정하기 위한 여러 가지 방법들이 제안되었다. 본 연구에서는 간단히 시험편을 제작할 수 있고, 압축 하중으로 암석의 모드 I 파괴인성을 구할 수 있도록 쐐기 분열 시험편(Wedge Splitting Test Specimen)을 이용하였다. 쐐기의 작용하는 연직하중에 의한 수직 및 수평하중 모두를 고려한 균열선단의 응력상태로부터 수치해석을 통한 응력확대계수의 식을 제안하였다. 쐐기 분열 시험편에 의해 구해진 모드 I 파괴인성값을 GD와 SENB 시험편의 시험 결과와 비교하여 시험법의 타당성을 확인하였다.

• Geomechanics and Engineering
• /
• 제10권6호
• /
• pp.775-791
• /
• 2016

The indirect tensile strengths (ITSs) of different cemented paste backfill mixes with different curing times were determined by considering crack initiation and fracture toughness concepts under different loading conditions of steel loading arcs with various contact angles, flat platens and the standard Brazilian test jaw. Because contact area of the ITS test discs developes rapidly and varies in accordance with the deformability, ITSs of curing materials were not found convenient to determine under the loading apparatus with indefinite contact angle. ITS values increasing with an increase in contact angle can be measured to be excessively high because of the high contact angles resulted from the deformable characteristics of the soft paste backfill materials. As a result of the change of deformation characteristics with the change of curing time, discs have different contact conditions causing an important disadvantage to reflect the strength change due to the curing reactions. In addition to the experimental study, finite element analyses were performed on several types of disc models under various loading conditions. As a result, a comparison between all loading conditions was made to determine the best ITSs of the cemented paste backfill materials. Both experimental and numerical analyses concluded that loading arcs with definite contact angles gives better results than those obtained with the other loading apparatus without a definite contact angle. Loading arcs with the contact angle of $15^{\circ}$ was found the most convenient loading apparatus for the typical cemented paste backfill materials, although it should be used carefully considering the failure cracks for a valid test.

• 콘크리트학회논문집
• /
• 제14권1호
• /
• pp.58-66
• /
• 2002

본 연구의 목적은 유효탄성균열모델과 점성균열모델의 개념에 기초한 임계응력확대계수, 임계균열단개구변위와 파괴에너지, 이선형 연화 곡선같은 콘크리트의 파괴특성들을 초기재령 콘크리트에 관해 구명하는 것이다. 이를 위해 모드 I의 파괴를 일으킬 수 있는 쐐기쪼갬시험이 노치가 있는 육각형의 쐐기 시험체에 대하여 수행되었다. 강도와 재령의 변화에 따라 하중-균열입구변위 곡선이 얻어졌으며, 이것은 선형탄성 파괴역학과 유한요소법에 의해 분석되었다. 실험 결과를 분석한 결과, 재령 1일부터 재령 28일까지의 임계응력확대계수와 파괴에너지는 증가하였으며, 임계균열단개구변위는 감소하였다. 또한 수치해석을 통하여 재령 1일부터 재령 28일까지의 이선형 연화 곡선의 네 파라미터를 구할 수 있었다. 이렇게 얻어진 초개재령 콘크리트의 파괴특성치와 이선형 연화 곡선은 초기재령 콘크리트의 파괴 기준과 유한요소해석시의 입력 상수로서 사용될 수 있을 것이다.