• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2005년도 창립60주년 기념 추계 학술대회
• /
• pp.147.1-147.1
• /
• 2005

• 대한기계학회논문집A
• /
• 제30권6호
• /
• pp.705-712
• /
• 2006

A finite element method based on the penalized subdomain-interface framework is proposed for fully-coupled, nonlinear thermomechanical analyses with thermal contact anuor radiation boundaries. In the variational formulation, a well-known penalty functional scheme is adopted for connecting subdomains and interfaces that satisfy various continuity requirements. As a logical consequence, the whole domain can be arbitrarily divided into independently-modeled subdomains without considering the conformity of meshes along their interfaces. Since the nonlinearities due to the contact and radiation boundaries can be localized within a few subdomains, the computational efficiency of the present method is greatly increased with appropriate solution algorithms. By solving some numerical problems, these advantageous features are confirmed carefully.

• Journal of Electrical Engineering and Technology
• /
• 제12권2호
• /
• pp.778-789
• /
• 2017

This paper presents an improved subdomain method to predict the magnet field distributions and electromagnetic performance of the surface-mounted permanent magnet (SPM) machines with static or dynamic eccentricity. Conventional subdomain models are either based on the scalar magnet potential to predict the rotor eccentricity effect or dependent on the magnetic vector potential without considering the eccentric rotor. In this paper, both the magnetic vector potential and the perturbation theory are introduced in order to accurately calculate the effect of rotor eccentricity on the open-circuit and armature reaction performance. The calculation results are presented and validated by the corresponding finite-element method (FEM) results.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1995년도 추계학술대회 논문집
• /
• pp.625-628
• /
• 1995

A static/implicit finite element code for sheet forming (ITAS3D) is parallelized on IBM SP 6000 multi-processor computer. Computing-load-balanced domain decomposition method and the direct solution method at each subdomain (and interface) equation are developed. The system of equations for each subdomain are constructed by condensation and calculated on each processor. Approximated operation counts are calculated to set up the nonlinear equation system for balancing the compute load on each subdomain. Th esquare cup tests with several numbers of elements are used in demonstrating the performance of this parallel implementation. This procedure are proved to be efficient for moderate number of processors, especially for large number of elements.

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

Substructuring methods are usually used in finite element structural analyses. In this study a multi-level substructuring algorithm is developed and proposed as a possible candidate for incompressible fluid solves. Finite element formulation for incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et.al.[5]. The present algorithm consists of four stages such as a gathering stage, a condensing stage, a solving stage and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At highest level, each subdomain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each subdomain has been replaced by a sequential static condensation. The global algebraic system arising feom the assembly of each subdomains is solved using Conjugate Gradient Squared(CGS) method. In this case, pre-conditioning techniques usually accompanied by iterative solvers are not needed.

• 한국항공우주학회지
• /
• 제33권11호
• /
• pp.7-14
• /
• 2005

접촉 현상을 포함하는 완전 연계된 열기계학적 문제에 대한 정확하고 효율적인 해석을 위하여 부영역과 공유면에 근거한 유한요소 정식화 기법을 제안하였다. 부영역과 공유면을 결합하기 위한 등식 적합 조건을 벌칙 함수로 처리함으로써 모든 유효 강성 행렬이 양정치화되며, 역행렬과 같은 각종 수치 연산이 매우 간편하다. 또한 전체 구조 형상이 복잡하더라도 대상 영역을 임의의 부영역으로 분할한 후, 공유면에서의 절점 연속성 등을 고려하지 않고 각 부영역을 독립적으로 유한요소 모델링할 수 있다. 컴퓨터 코드의 개발 및 수치 예제의 해석을 통하여 본 기법에 대한 기본적인 특성을 확인하였다.

• Structural Engineering and Mechanics
• /
• 제65권4호
• /
• pp.359-368
• /
• 2018

Dynamic analysis of complex and large structures may be costly from a numerical point of view. For coupled vibroacoustic finite element models, the importance of reducing the size becomes obvious because the fluid degrees of freedom must be added to the structural ones. In this paper, a component mode synthesis method is proposed for large vibroacoustic interaction problems. This method couples fluid subdomains and dynamical substructuring of Craig and Bampton type. The acoustic formulation is written in terms of the velocity potential, which implies several advantages: coupled algebraic systems remain symmetric, and a potential formulation allows a direct extension of Craig and Bampton's method to acoustics. Those properties make the proposed method easy to implement in an existing finite element code because the local numerical treatment of substructures and fluid subdomains is undifferentiated. Test cases are then presented for axisymmetric geometries. Numerical results tend to prove the validity and the efficiency of the proposed method.

• 한국대기환경학회지
• /
• 제9권3호
• /
• pp.242-246
• /
• 1993

Theoretically, spectral method has the highest accuracy among present numerical methods, but it is generally difficult to apply to complex terrains because of complex boundary conditions. Recently, spectral-element method, basically divide the domain into a set of rectangular subdomain and solve the equation at each subdomain, has been introduced. However, boundary conditions become more complex and requires more computing time, thus spectral-element method is not powerful for all complex terrain problems. In this paper, potential flow theory was intorduced to solve the air flows and diffusion phenomenon in the presence of terrain obstacles. Using the velocity potential-stream line orthogonal coordinate space, the diffusion problems of hilly terrain by pseudospectral method were solved and compared those with no terrain real scale solutions.

• 대한조선학회지
• /
• 제24권3호
• /
• pp.1-8
• /
• 1987

In this paper, the localized finite element method(LFEM) is applied to 3-dimensional ship motion problems in water of infinite depth. The LFEM used here is based on the functional constructed by Bai & Yeung(1974). To test the present numerical scheme, a few vertical axisymmetric bodies are treated by general 3-dimensional formulation. The computed results of hydrodynamic coefficients for a few vertical spheroids and vertical circular cylinders show good agreement with results obtained by others. The advantages of the present numerical method compared with the method of integral equation are as follows; (i) The cumbersome existence of irregular frequencies in the method of conventional integral equation is removed. (ii) The final matrix is banded and symmetric and the computation of the matrix elements is comparatively easier, whereas the size of the matrix in the present scheme is much larger. (iii) In the future research, it is possible to accommodate with the nonlinear exact free surface boundary condition in the localized finite element subdomain, whereas the linear solution is assumed in the truncated(far field) subdomain.

• Journal of Mechanical Science and Technology
• /
• 제20권1호
• /
• pp.110-124
• /
• 2006

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.