• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2002년도 가을 학술발표회 논문집
• /
• pp.575-582
• /
• 2002

A new coupling method of Element-Free Galerkin Method(EFGM) and Boundary Element Method(BEM) using the domain decomposition method is presented in this paper. This proposed methodology is that the problem domain is decomposed into sub-domains being modeled by the EFGM and BEM respectively and the respective EFGM and BEM domains share a partially overlapped region over an entire domain. Then, the each sub-domain is separately computed and the variables on common region are iteratively updated until converging. It is an important note that in the developed coupling method, there is no need to combine the coefficient matrices of EFGM and BEM sub-domains, in contrast with the other conventional coupling methods. In the first part of this paper, a theory of EFGM and BEM is summarized, and then a brief introduction of domain decomposition method is described. Then, a new coupling method is presented. Also, patch test and Some numerical examples are studied to verify stability, accuracy and efficiency of the proposed method, in which numerical performance of the method is compared with that of conventional method such as EFGM-BEM variational coupling method, EFGM and BEM.

• 한국CDE학회논문집
• /
• 제12권4호
• /
• pp.255-262
• /
• 2007

In the finite element analysis of forming process, objects are described with a finite number of elements and nodes and the approximated solutions can be obtained by the variational principle. One of the shortcomings of a finite element analysis is that the structure of mesh has become inefficient and unusable because discretization error increases as deformation proceeds due to severe distortion of elements. If the state of current mesh satisfies a certain remeshing criterion, analysis is stopped instantly and resumed with a reconstructed mesh. In the study, a new remeshing algorithm using tetrahedral elements has been developed, which is adapted to the desired mesh density. In order to reduce the discretization error, desired mesh sizes in each lesion of the workpiece are calculated using the Zinkiewicz and Zhu's a-posteriori error estimation scheme. The pre-constructed mesh is constructed based on the modified point insertion technique which is adapted to the density function. The object domain is divided into uniformly-sized sub-domains and the numbers of nodes in each sub-domain are redistributed, respectively. After finishing the redistribution process of nodes, a tetrahedral mesh is reconstructed with the redistributed nodes, which is adapted to the density map and resulting in good mesh quality. A goodness and adaptability of the constructed mesh is verified with a testing measure. The proposed remeshing technique is applied to the finite element analyses of forging processes.

• Structural Engineering and Mechanics
• /
• 제38권3호
• /
• pp.261-281
• /
• 2011

A new method of formulation of a class of elements that are immune to mesh distortion effects is proposed here. The simple three-noded bar element with an offset of the internal node from the element center is employed here to demonstrate the method and the principles on which it is founded upon. Using the function space approach, the modified formulation is shown here to be superior to the conventional isoparametric version of the element since it satisfies the completeness requirement as the metric formulation, and yet it is in agreement with the best-fit paradigm in both the metric and the parametric domains. Furthermore, the element error is limited to only those that are permissible by the classical projection theorem of strains and stresses. Unlike its conventional counterpart, the modified element is thus not prone to any errors from mesh distortion. The element formulation is symmetric and thus satisfies the requirement of the conservative nature of problems associated with all self-adjoint differential operators. The present paper indicates that a proper mapping set for distortion immune elements constitutes geometric and displacement interpolations through parametric and metric shape functions respectively, with the metric components in the displacement/strain replaced by the equivalent geometric interpolation in parametric co-ordinates.

• 대한기계학회논문집
• /
• 제16권8호
• /
• pp.1429-1437
• /
• 1992

본 연구에서는 이완형 물성방정식을 바탕으로 하며 프와송 비가 일정하다는 가정을 하지 않는다. 또한 점탄성 지배방정식에 변분원리를 적용하고 유도되어진 식 에 유한요소해법을 사용하여 시스템 기본해석을 위한 연립방정식을 유도한다. 이와 함께 점탄성 물성함수의 유도 및 응력계산을 위한 공식화 과정도 설명한다. 제시된 방법론의 타당성 및 정확성을 보이기 위해서 평면응력 및 평면변형 문제의 변위 및 응력을 수치해석하여 이론해와 비교 검토하며, 아울러 시간증분의 변화와 Gauss poi- nts수가 수치정확도에 끼치는 영향을 조사한다.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• 제16권9호
• /
• pp.2927-2941
• /
• 2022

The primary task of machine fault diagnosis is to judge whether the current state is normal or damaged, so it is a typical binary classification problem with mutual exclusion. Mutually exclusive events and out-of-domain detection have one thing in common: there are two types of data and no intersection. We proposed a fusion model method to improve the accuracy of machine fault diagnosis, which is based on the mutual exclusivity of events and the commonality of out-of-distribution detection, and finally generalized to all binary classification problems. It is reported that the performance of a convolutional neural network (CNN) will decrease as the recognition type increases, so the variational auto-encoder (VAE) is used as the primary model. Two VAE models are used to train the machine's normal and fault sound data. Two reconstruction probabilities will be obtained during the test. The smaller value is transformed into a correction value of another value according to the mutually exclusive characteristics. Finally, the classification result is obtained according to the fusion algorithm. Filtering normal data features from fault data features is proposed, which shields the interference and makes the fault features more prominent. We confirm that good performance improvements have been achieved in the machine fault detection data set, and the results are better than most mainstream models.

• 한국전산구조공학회논문집
• /
• 제20권3호
• /
• pp.339-345
• /
• 2007

본 논문에서는 등기하 해석법을 이용하여 평면 탄성문제의 변분식을 유도하였다. 등기하 해석법은 새로이 부각되고 있는 해석법으로서 기저 함수가 NURBS(Non-Uniform Rational B-Splines) 로부터 직접 생성되므로 해 공간은 CAD 모델을 구성하는 함수로써 표현된다. 또한 CAD 모델의 B-Spline 기저 함수를 직접 사용하므로 기하학적으로 엄밀한 형상을 표현할 수 있고 요소망의 재구성 없이 해석모델을 정밀화(Refinement)할 수 있는 강점이 있다. 본 논문에서는 이를 확장하여 연속체 기반의 애드조인트 설계 민감도 해석법을 사용하는 등기하 설계민감도 해석법을 유도하였다. 기존의 유한요소 기반형상 최적설계는 형상의 매개화에 어려움을 겪었으나 등기하 기반 최적설계에서는 기하학적 정보가 이미 B-spline 기저함수와 조정점에 포함되어 있으므로 이러한 어려움을 피할 수 있는 잠재력을 가지고 있다. 몇몇 수치 예제를 통해서 등기하 해석법을 사용한 설계 민감도 해석을 수행하였으며 유한차분 민감도와 비교하여 정확성을 확인하였다.

• Journal of Advanced Research in Ocean Engineering
• /
• 제1권3호
• /
• pp.157-167
• /
• 2015

In order to provide simple and accurate wave theory in design of offshore structure, an analytical approximation is introduced in this paper. The solution is limited to flat bottom having a constant water depth. Water is considered as inviscid, incompressible and irrotational. The solution satisfies the continuity equation, bottom boundary condition and non-linear kinematic free surface boundary condition exactly. Error for dynamic condition is quite small. The solution is suitable in description of breaking waves. The solution is presented with closed form and dispersion relation is also presented with closed form. In the last century, there have been two main approaches to the nonlinear problems. One of these is perturbation method. Stokes wave and Cnoidal wave are based on the method. The other is numerical method. Dean's stream function theory is based on the method. In this paper, power series method was considered. The power series method can be applied to certain nonlinear differential equations (initial value problems). The series coefficients are specified by a nonlinear recurrence inherited from the differential equation. Because the non-linear wave problem is a boundary value problem, the power series method cannot be applied to the problem in general. But finite number of coefficients is necessary to describe the wave profile, truncated power series is enough. Therefore the power series method can be applied to the problem. In this case, the series coefficients are specified by a set of equations instead of recurrence. By using the set of equations, the nonlinear wave problem has been solved in this paper.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• 대한기계학회논문집B
• /
• 제26권12호
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• 한국표면공학회지
• /
• 제53권3호
• /
• pp.95-103
• /
• 2020

In this study, the material used in the hot dip galvanizing equipment was poorly corrosion-resistant, so it was performed to solve the cost and time problems caused by equipment replacement. The theoretical calculation was performed using the DV-Xα method(Discrete Variational Local-density approximation method). The alloy (STS4XX series) of the equipment currently used has a martensite phase. Therefore, the theoretical calculation was performed by applying P4 / mmm, which is a tetragonal structure. The new alloy was chosen by designing theoretical values close to existing materials. Considering elements that contribute to corrosion, most have high prices. Therefore, the design was completed by adjusting the content using only the components of the reference material in the theoretical design. The final design alloys were chosen as D6 and D9. Designed D6 and D9 were dissolved and prepared using an induction furnace. After the heat treatment process was completed, the corrosion rate of the alloys was confirmed by using the potentiodynamic polarization test. The surface of the prepared alloys were processed horizontally and then polished to # 1200 using sand paper to perform potentiodynamic polarization test. Domestic products: 4.735 mpy (mils / year), D6: 0.9166 mpy, D9: 0.3372 mpy, alloys designed than domestic products had a lower corrosion rate. Therefore, the designed alloy was expected to have better erosion resistance.