• 소음진동
• /
• 제9권1호
• /
• pp.103-112
• /
• 1999

A new modeling and analysis technique for one-dimensional distributed parameter systems is presented. First. discretized equations of motion in Laplace domain are derived by applying discretization methods for partial differential equations of a one-dimensional structure with respect to spatial coordinate. Secondly. the z and inverse z transformations are applied to the discretized equations of motion for obtaining a dynamic matrix for a uniform element. Four different discretization methods are tested with an example. Finally, taking infinite on the number of step for a uniform element leads to an exact dynamic matrix for the uniform element. A generalized modal analysis procedure for eigenvalue analysis and modal expansion is also presented. The resulting element dynamic matrix is tested with a numerical example. Another application example is provided to demonstrate the applicability of the proposed method.

• International Journal of Aeronautical and Space Sciences
• /
• 제13권2호
• /
• pp.127-153
• /
• 2012

This paper presents an advanced computational method for the prediction of the responses in the frequency domain of general linear dissipative structural-acoustic and fluid-structure systems, in the low-and medium-frequency domains and this includes uncertainty quantification. The system under consideration is constituted of a deformable dissipative structure that is coupled with an internal dissipative acoustic fluid. This includes wall acoustic impedances and it is surrounded by an infinite acoustic fluid. The system is submitted to given internal and external acoustic sources and to the prescribed mechanical forces. An efficient reduced-order computational model is constructed by using a finite element discretization for the structure and an internal acoustic fluid. The external acoustic fluid is treated by using an appropriate boundary element method in the frequency domain. All the required modeling aspects for the analysis of the medium-frequency domain have been introduced namely, a viscoelastic behavior for the structure, an appropriate dissipative model for the internal acoustic fluid that includes wall acoustic impedance and a model of uncertainty in particular for the modeling errors. This advanced computational formulation, corresponding to new extensions and complements with respect to the state-of-the-art are well adapted for the development of a new generation of software, in particular for parallel computers.

• 대한기계학회논문집B
• /
• 제23권2호
• /
• pp.243-253
• /
• 1999

The dual reciprocity boundary element method (DRBEM) is used to solve the Graetz problem of laminar flow inside circular duct. In this method the domain integral tenn of boundary integral equation resulting from source term of governing equation is transformed into equivalent boundary-only integrals by using the radial basis interpolation function, and therefore complicate domain discretization procedure Is completely removed. Velocity profile is obtained by solving the momentum equation first and then, using this velocities as Input data, energy equation Is solved to get the temperature profile by advancing from duct entrance through the axial direction marching scheme. DRBEM solution is tested for the uniform temperature and heat flux boundary condition cases. Local Nusselt number, mixed mean temperature and temperature profile inside duct at each dimensionless axial location are obtained and compared with exact solutions for the accuracy test Solutions arc in good agreement at the entry region as well as fully developed region of circular duct, and their accuracy are verified from error analysis.

• Korea-Australia Rheology Journal
• /
• 제19권4호
• /
• pp.233-242
• /
• 2007

We performed a numerical investigation to find out the optimal choice of the spatial discretization in the distributed-Lagrangian-multiplier/fictitious-domain (DLM/FD) method for the solid/fluid interaction problem. The elastic solid bar attached on the bottom in a pressure-driven channel flow of a Newtonian fluid was selected as a model problem. Our formulation is based on the scheme of Yu (2005) for the interaction between flexible bodies and fluid. A fixed regular rectangular discretization was applied for the description of solid and fluid domain by using the fictitious domain concept. The hydrodynamic interaction between solid and fluid was treated implicitly by the distributed Lagrangian multiplier method. Considering a simplified problem of the Stokes flow and the linearized elasticity, two numerical factors were investigated to clarify their effects and to find the optimum condition: the distribution of Lagrangian multipliers and the solid/fluid interfacial condition. The robustness of this method was verified through the mesh convergence and a pseudo-time step test. We found that the fluid stress in a fictitious solid domain can be neglected and that the Lagrangian multipliers are better to be applied on the entire solid domain. These results will be used to extend our study to systems of elastic particle in the Stokes flow, and of particles in the viscoelastic fluid.

• 대한수학회논문집
• /
• 제11권2호
• /
• pp.495-501
• /
• 1996

There is currently a regain of interest in ADI (Alternating Direction Implicit) method as a preconditioner for iterative Method for solving large sparse linear systems, because of its suitability for parallel computation. However the classical ADI is not applicable to FE(Finite Element) matrices. In this paper wer propose a Block-ADI method, which is applicable to Finite Element metrices. The new approach is a combination of classical ADI method and domain decompositi on. Also, we provide a partial proof of the convergence based on the results from the regular splittings, in case the discretization metrix is symmetric positive definite.

• 한국통신학회논문지
• /
• 제21권10호
• /
• pp.2724-2730
• /
• 1996

In this paper, the dispersive characteristics of the Finite-Difference Frqequency-Domain method based on the Multi-Resolution Technique(MR-FDFD) are numerically analyzed. A dispersion analysis of the MR-FDFD ority of the MR-FDFD method to the spatial discretization is shown. We expect that the multi-resoluation technique will improve the disavantage of the finite difference techqnique which needs the large comutational memory for accurate electromagnetic analysis.

• Structural Engineering and Mechanics
• /
• 제27권6호
• /
• pp.683-696
• /
• 2007

This article is concerned with the presentation of a time-domain BEM approach applied to the solution of the scalar wave equation for 2D problems. The basic idea is quite simple: the basic variables of the problem at time $t_n$ (potential and flux) are computed with the results related to the potential and to its time derivative at time $t_{n-1}$ playing the role of "initial conditions". This time-marching scheme needs the computation of the potential and its time derivative at all boundary nodes and internal points, as well as the entire discretization of the domain. The convolution integrals of the standard time-domain BEM formulation, however, are not computed; the matrices assembled, only at the initial time interval, are those related to the potential, flux and to the potential time derivative. Two examples are presented and discussed at the end of the article, in order to verify the accuracy and potentialities of the proposed formulation.

• 한국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.

• 한국전산구조공학회논문집
• /
• 제19권1호
• /
• pp.39-46
• /
• 2006

본 연구는 파 해석에 있어서 공간-시간 분할 개념을 도입하여 켈러킨 방법으로 해석하였다. 공간-시간 유한요소법은 오직 공간에 대해서만 분할하는 일반적인 유한요소법보다 간편하다. 비교적 큰 시간간격에 대해서 공간과 시간을 동시에 분할하는 방법을 제시하며 가중잔차법이 공간-시간 영역에서 유한요소 정식화에 이용되었다. 큰 시간 간격으로 인하여 문제의 해가 발산하는 경우가 동적인 문제에서 흔히 발생한다. 이러한 결점을 보완한 사각형 공간-시간 요소를 취하여 문제를 해석하고 해의 안정에 대해 기술하였다. 다수의 수치해석을 통하여 이 방법이 효과적 임을 알 수 있었다.

• Structural Engineering and Mechanics
• /
• 제75권6호
• /
• pp.685-699
• /
• 2020

Polygonal finite element provides a great flexibility in mesh generation of crack propagation problems where the topology of the domain changes significantly. However, the control of the discretization error in such problems is a main concern. In this paper, a polygonal-FEM is presented in modeling of crack propagation problems via an automatic adaptive mesh refinement procedure. The adaptive mesh refinement is accomplished based on the Zienkiewicz-Zhu error estimator in conjunction with a weighted SPR technique. Adaptive mesh refinement is employed in some steps for reduction of the discretization error and not for tracking the crack. In the steps that no adaptive mesh refinement is required, local modifications are applied on the mesh to prevent poor polygonal element shapes. Finally, several numerical examples are analyzed to demonstrate the efficiency, accuracy and robustness of the proposed computational algorithm in crack propagation problems.