• Geophysics and Geophysical Exploration
• /
• v.7 no.4
• /
• pp.234-244
• /
• 2004

In order to calculate 3-dimensional wavefield using finite-element method in frequency domain, we must factor so huge sparse impedance matrix. Because of difficulties of handling of this huge impedance matrix, 3-dimensional wave equation modeling is conducted mainly in time domain. In this study, we simulate the 3-D wavefield using finite-element method in Laplace domain by combining high-performance parallel finite-element solver and SWEET (Suppressed Wave Equation Estimation of Traveltime) algorithm which can calculate the traveltime and the amplitude. To verify this combination, we applied it to the SEG/EAGE 3D salt model in serial and parallel computing environments.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.04a
• /
• pp.172-175
• /
• 2005

The numerical simulation is performed for the acoustic emission and the wave propagation due to fiber breakage in single fiber composite plates by the finite element transient analysis. The acoustic emission and the following wave motions from a fiber breakage under a static loading is simulated to investigate the applicability of the explicit finite element method and the equivalent volume force model as a simulation tool of wave propagation and a modeling technique of an acoustic emission. For such a simple case of the damage event under static loading, various parameters affecting the wave motion are investigated for reliable simulations of the impact damage event. The high velocity and the small wave length of the acoustic emission require a refined analysis with dense distribution of the finite element and a small time step. In order to fulfill the requirement for capturing the exact wave propagation and to cover the 3-D simulation, we utilize the parallel FE transient analysis code and the parallel computing technology.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.12
• /
• pp.889-898
• /
• 2022

In this paper, parallel computing method on the three-dimensional electromagnetic field is proposed. The present electromagnetic scattering analysis is conducted based on the time-harmonic vector wave equation and the finite element method. The edge-based element and 2^nd -order absorbing boundary condition are used. Parallelization of the elemental numerical integration and the matrix assemblage is accomplished by allocating the partitioned finite element subdomain for each processor. The graph partitioning library, METIS, is employed for the subdomain generation. The large sparse matrix computation is conducted by MUMPS, which is the parallel computing library based on the multi-frontal method. The accuracy of the present program is validated by the comparison against the Mie-series analytical solution and the results by ANSYS HFSS. In addition, the scalability is verified by measuring the speed-up in terms of the number of processors used. The present electromagnetic scattering analysis is performed for a perfect electric conductor sphere, isotropic/anisotropic dielectric sphere, and the missile configuration. The algorithm of the present program will be applied to the finite element and tearing method, aiming for the further extended parallel computing performance.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.389-395
• /
• 2004

We consider the interface between the parallel distributed memory multifrontal solver and the finite element method. We give in detail the requirement and the data structure of parallel FEM interface which includes the element data and the node array. The full procedures of solving a large scale structural problem are assumed to have pre-post processors, of which algorithm is not considered in this paper. The main advantage of implementing the parallel FEM interface is shown up in the case that we use a distributed memory system with a large number of processors to solve a very large scale problem. The memory efficiency and the performance effect are examined by analyzing some examples on the Pegasus cluster system.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.8 no.4
• /
• pp.1-7
• /
• 1998

This paper describes a meshless of element-free method based on fuzzy knowledge processing. To efficiently simulate complicated physical phenomena with dynmics and non-linear ploblem using computational mechanics, special method is required such as parallel processing or adaptive analysis techniques. However, the conventional finite element method is too complicated to be employed in the above cases. In order to reduce the above complexity of the conventional finite element analysis systms, the so called meshles finite elements as an input information have been stuided. Node is generated if its distance form existing node points is similar to the node spacing fuction at the point. The node spacing function is well controlled by the fuzzy knowledge processing Practical performances of the present system are demonstrated through several three-dimensional(3D) problems.

• Communications of the Korean Mathematical Society
• /
• v.11 no.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.

• Computers and Concrete
• /
• v.11 no.6
• /
• pp.541-569
• /
• 2013

There is a wide variety of existing Genetic Algorithms (GA) operators and parameters in the literature. However, there is no unique technique that shows the best performance for different classes of optimization problems. Hence, the evaluation of these operators and parameters, which influence the effectiveness of the search process, must be carried out on a problem basis. This paper presents a comparison for the influence of GA operators and parameters on the performance of the damage identification problem using the finite element model updating method (FEMU). The damage is defined as reduction in bending rigidity of the finite elements of a reinforced concrete beam. A certain damage scenario is adopted and identified using different GA operators by minimizing the differences between experimental and analytical modal parameters. In this study, different selection, crossover and mutation operators are compared with each other based on the reliability, accuracy and efficiency criteria. The exploration and exploitation capabilities of different operators are evaluated. Also a comparison is carried out for the parallel and sequential GAs with different population sizes and the effect of the multiple use of some crossover operators is investigated. The results show that the roulettewheel selection technique together with real valued encoding gives the best results. It is also apparent that the Non-uniform Mutation as well as Parent Centric Normal Crossover can be confidently used in the damage identification problem. Nevertheless the parallel GAs increases both computation speed and the efficiency of the method.

• Structural Engineering and Mechanics
• /
• v.77 no.3
• /
• pp.369-381
• /
• 2021

The seismic responses of elevated tanks considering liquid-structure interaction are presented under horizontal earthquake. The scaled model tank is fabricated to study the dynamic responses of anchored tank and newly designed uplift tank with replaced dampers. The natural frequencies for structural mode are obtained by modal analysis. The dynamic responses of tanks are completed by finite element method, which are compared with the results from experiment. The displacement parallel and perpendicular to the excitation direction are both gained as well as structural acceleration. The strain of tank walls and the axial strain of columns are also obtained afterwards. The seismic responses of liquid storage tank can be calculated by the finite element model effectively and the results match well with the one measured by experiment. The aim is to provide a new type of tank system with vertical constraint relaxed which leads to lower stress level. With the liquid volume increasing, the structural fundamental frequency has a great reduction and the one of uplift tank are even smaller. Compared with anchored tank, the displacement of uplift tank is magnified, the strain for tank walls and columns parallel to excitation direction reduces obviously, while the one perpendicular to earthquake direction increases a lot, but the values are still small. The stress level of new tank seems to be more even due to uplift effect. The new type of tank can realize recoverable function by replacing dampers after earthquake.

• Structural Engineering and Mechanics
• /
• v.14 no.6
• /
• pp.625-647
• /
• 2002

Parallel execution of computational mechanics codes requires efficient mesh-partitioning techniques. These mesh-partitioning techniques divide the mesh into specified number of submeshes of approximately the same size and at the same time, minimise the interface nodes of the submeshes. This paper describes a new mesh partitioning technique, employing Genetic Algorithms. The proposed algorithm operates on the deduced graph (dual or nodal graph) of the given finite element mesh rather than directly on the mesh itself. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialise an optimisation of the graph partition problem. In practice, hierarchy of (usually more than two) graphs are used to obtain the final graph partition. The proposed partitioning algorithm is applied to graphs derived from unstructured finite element meshes describing practical engineering problems and also several example graphs related to finite element meshes given in the literature. The test results indicate that the proposed GA based graph partitioning algorithm generates high quality partitions and are superior to spectral and multilevel graph partitioning algorithms.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 2000.10a
• /
• pp.410-415
• /
• 2000

In this paper feasible direction method which is one of the optimization method is adopted to natural frequency optimization. In order to find the optimum design of structures that have characteristic natural frequency range, a numerical optimization method to solving eigenvalue problems is a widely used approach. However most cases, it is difficult to decide the accurate thickness and shape of structures that have allowable natural frequency in design constraints. Parallel analysis algorithm involving the feasible direction optimization method and Rayleight-Ritz eigenvalue solving method is developed. The method is implemented by using finite element method. It calculated the optimal thickness and the thickness ratio of each element of 2-D plane element through the parallel algorithm method which satisfy the design constraint of natural frequency.