• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.4
• /
• pp.361-368
• /
• 2010

In this paper, the parallel algorithm using MPI(Message-Passing Interface) library is introduced in order to improve numerical efficiency for the staggered method for nearly incompressible and impermeable porous media which was introduced by Park and Tak(2010). The porous media theory and the staggered method are also briefly introduced in this paper. Moreover, we account for MPI library for blocking, non-blocking, and collective communication, and propose combined the staggered method with the blocking and nonblocking MPI library. And then, we present how to allocate CPUs on the staggered method and the MPI library, which is related with the numerical efficiency in order to solve unknown variables on nearly incompressible and impermeable porous media. Finally, the results comparing serial solution with parallel solution are verified by 2 dimensional saturated porous model according to the number of FEM meshes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.5
• /
• pp.475-482
• /
• 2010

In this paper, an efficient two-level parallel domain decomposition algorithm is suggested to solve large-DOF structural problems. Each subdomain is composed of the coarse problem and local problem. In the coarse problem, displacements at coarse nodes are computed by the iterative method that does not need to assemble a stiffness matrix for the whole coarse problem. Then displacements at local nodes are computed by Multi-Frontal Sparse Solver. A parallel version of PCG(Preconditioned Conjugate Gradient Method) is developed to solve the coarse problem iteratively, which minimizes the data communication amount between processors to increase the possible problem DOF size while maintaining the computational efficiency. The test results show that the suggested algorithm provides scalability on computing performance and an efficient approach to solve large-DOF structural problems.

• Structural Engineering and Mechanics
• /
• v.26 no.3
• /
• pp.343-361
• /
• 2007

This article presents a numerical tool for dimensioning two-threaded fasteners connecting prismatic parts subjected to fatigue tension loads that are coplanar with the screw axis. A simplified numerical model is developed from unidirectional finite elements, modeling the connected parts and screws with bent elements and the elastic contact layer between the parts with springs. An algorithm updating the contact stiffness matrix, calculating forces and displacements at each node of the structure and thus normal stresses in the screws in both static and fatigue is further developed using C language. An experimental study is also conducted in parallel with the numerical approach to validate the developed model assumptions, the numerical model and the 3D finite element results. Since stiffness values for the compressive zones in the parts are analytically difficult to determine, a statistical software method is used, from which a tuning factor is derived for identifying these stiffness values. The method is also applied to set out the influence of each parameter on the fatigue behaviour of each screw. Finally, the developed model will be used to establish a new, sophisticated, fast and accurate tool for dimensioning bolted mechanical structures.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.5
• /
• pp.482-490
• /
• 2009

In this paper, the author proposes a new method for acoustic radiation optimum design to minimize noise from a vibrating panel-like structure using a collaborative population-based search method called the particle swarm optimization algorithm(PSOA). The PSOA is a parallel evolutionary computation technique initially developed by Kennedy and Eberhart. The acoustic radiation optimization method based on the PSOA consists of two processes. In the first process, the acoustic radiation analysis by an integrated p-version FEM/BEM, which was developed by using MATLAB, is performed to evaluate the exterior acoustic radiation field of the panel. The second process is to search the optimum design variables: 1) Shape of Bezier curves and 2) Shape and position of ribs, to minimize noise from the panel using the PSOA. The optimization method based on the PSOA is compared to that based on the steady state genetic algorithm(SSGA) in order to verify the effectiveness and validity of the optimal solution by PSOA. Finally, it is shown that the optimal designs of the panel obtained by using the PSOA can achieve effective reductions in radiated sound power.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.1
• /
• pp.59-71
• /
• 2012

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in unbounded isotropic elastic solids containing multiple interacting isotropic or anisotropic diamond-shaped inclusions subject to remote uniaxial tension. The method is applied to two-dimensional problems involving long parallel diamond-shaped cylindrical inclusions. A detailed analysis of the stress field at the interface between the matrix and the central inclusion is carried out for square and hexagonal packing of the inclusions. The effects of the number of isotropic or anisotropic diamond-shaped inclusions and of the various fiber volume fractions for the circular inclusions circumscribing its respective diamond-shaped inclusion on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy and efficiency of the method are examined through comparison with results obtained using the finite element method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.555-564
• /
• 1994

This study has been made to investigate the stress distribution around defects and inclusions that behave as stress concentrators. The stress distribution and interation effects around defects and inclusions was analyzed using Finite Element Method. The results are as follows;(1) Maximum stress point in case of $E_I/E_M>1$($E_I$:elasticity modulus forthe inclusion, $E_M$/:elasticity modulus for the base material)is the vertical point with respect to force direction and in case of $E_I/E_M<1$ it is the parallel point along the hole edge. (2) Interaction effects of ${\sigma}_y$ for the inclusion side is larger than the defect side when the interval between inclusion and defect is near. (3) stress interation effects is large if the difference of ${\sigma}_y$ is small and it is small if the difference of ${\sigma}_y$ is large for the case that the interval between inclusion and defect whose size and property are different is near.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.5
• /
• pp.211-221
• /
• 1998

Thermal behavior on the cylinder block of a 4-cylinder, 4-stroke 2.0L SOHC gasoline engine was numerically and experimentally analyzed. The numerical calculation was performed using the finite element method. The cylinder block was modelled as a three dimensional finite element by considering its geometry. The physical domain was devided into hexahedron elements. 16 thermocouples were installed at points of 2mm inside from cylinder wall near top ring of piston in cylinder block, which points have suffered major thermal loads and suggested as proper measurement points for engine design by industrial engineers. Under full load and 9$0^{\circ}C$ coolant temperature condition, temperature behavior of cylinder block according to engine speed were analyzed. The results showed that temperature rose gradually to conform to a function of 2nd~4th order of engine speed at intake side, exhaust and siamese side, respectively. As engine load was changed from 100 to 50% by 25% step, temperature curve also conformed to 2nd~7th order function of engine speed. Temperature differences by load condition were similar among 100, 75% and 50%. Under full load and coolant temperature of 11$0^{\circ}C$, temperature behavior were also analyzed and the result also showed conformance to 2n d~7th order function of engine speed. Temperature curve was transferred in parallel upwards corresponding coolant temperature rise.

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.62 no.2
• /
• pp.72-78
• /
• 2013

In this study, samples from the site where there occurred unbalanced current when cable routing were analyzed, and the simulation program for electric power system analysis, PSCAD/EMTDC, was used to calculate the current unbalance on cable routing. Based on electromagnetic finite element analysis(FEA), electromagnetic parameters enabled the interlocking with COMSOL for the calculation of allowable current ampacity and magnetic filed distribution. This then led to modeling unbalanced current between common modes using the unbalanced current analysis program, thereby comparing and discussing the results from both. The analyzed model is a common mode 2 parallel circuit, which is a basic model for cable routing, and by arranging cables in various ways, the arrangement with the least current unbalance was suggested, which would, in the future, prevent earth faults and extend life for the whole cable.

• The korean journal of orthodontics
• /
• v.15 no.1
• /
• pp.7-22
• /
• 1985

The purpose of this study was to analyze the stress distribution and the displacement in the maxillary complex after the application of the reverse headgear. The direction of force was parallel to the occlusal plane. Orthopedic force,300gm, was applied to the maxilla of the dry human skull in a forward direction. The stress distribution and the displacement within the maxillary Complex was analyzed by a 3-dimensional finite element method. The results were as follows: 1. The stress distribution at the molar region was greater than that at the anterior. 2. The stress distribution at the lateral side of the premaxilla was greater than that at the middle aide, especially high stress was noted at the canine eminence. 9. Compressive stress was noted only at the frontozygomatic suture of the zygomatic arch. 4. A forward, upward, and sideward displacement was noted at the entire nodal points of the zygomaticomaxillary suture portion. A displacement with a slight rotation was observed on the transverse palatine suture. 5. The maximum stress was observed at the lateral side of the maxillary tuberosity area, and generally the forward and downward displacement was noted at all this area.

• Steel and Composite Structures
• /
• v.11 no.1
• /
• pp.39-58
• /
• 2011

The subject of the ongoing research work is to analyze the composite action of the structural elements of composite slabs with profiled steel decking by experimental and numerical studies. The mechanical and frictional interlocks result in a complex behaviour and failure under horizontal shear action. This is why the design characteristics can be determined only by standardized experiments. The aim of the current research is to develop a computational method which can predict the behaviour of embossed mechanical bond under shear actions, in order to derive the design characteristics of composite slabs with profiled steel decking. In the first phase of the research a novel experimental analysis is completed on an individual concrete encased embossment of steel strip under shear action. The experimental behaviour modes and failure mechanisms are determined. In parallel with the tests a finite element model is developed to follow the ultimate behaviour of this type of embossment, assuming that the phenomenon is governed by the failure of the steel part. The model is verified and applied to analyse the effect of embossment's parameters on the behaviour. In the extended investigation different friction coefficients, plate thicknesses, heights and the size effects are studied. On the basis of the results the tendencies of the ultimate behaviour and resistance by the studied embossment's characteristics are concluded.