• Korean Journal of Computational Design and Engineering
• /
• v.21 no.2
• /
• pp.143-150
• /
• 2016

In shipbuilding, hull assemblies are manufactured by welding. The thermal deformation caused by the welding produces shape deformation. Counter-deformed design methods have been used in shipyards to cope with the weld-induced deformation of ship assembles. Finite element methods (FEMs) are frequently used to estimate welding distortion in the counter-deformed design. For the estimation of welding distortion, producing uniform rectangular elements is required to enter thermal loads on the welding line and obtain accurate analysis results. In this paper, a new automatic mesh generation method is proposed for prediction of welding deformation in FEM. Meshes are constructed for test cases to demonstrate the feasibility of the proposed mesh generation method.

• Structural Engineering and Mechanics
• /
• v.41 no.4
• /
• pp.441-458
• /
• 2012

The shear strength is an important factor in the design of prestressed concrete beams. Therefore, researchers have utilized various methods to determine the shear strength of these elements for the design purposes. To evaluate some of the proposed theoretical methods, numerous models of post-tensioned beams with or without vertical prestressing are selected and analyzed using the finite element method and assuming nonlinear behavior for the materials. In this regard the validity of modeling is evaluated based on some tests results. In the second part of the study two beam specimens are built and tested and their load-deformation curve and cracking pattern are studied. The analytical results consist of compressive strut slope and mid span load deflection are compared with some experimental results, and the results of some codes' formulas. Finally comparing the results of nonlinear analysis with the experimental values, a new formula is proposed for determining strut slopes in prestressed concrete beams.

• Journal of Technologic Dentistry
• /
• v.42 no.4
• /
• pp.348-354
• /
• 2020

Purpose: Zirconia is differentiated from other ceramics because of its high resistance to corrosion and wear, excellent flexural strength (900~1400 MPa), and high hardness. Dental zirconia with proven mechanical/biological stability is suitable for the manufacture of implants. However, there are limited in vivo studies evaluating stress distribution in zirconia compared with that in titanium implants and studies analyzing finite elements. This study was conducted to evaluate the stress distribution of the supporting bone surrounding zirconia and titanium implants using the finite element analysis method. Methods: For finite element analysis, a single implant-supported restoration was designed. Using a universal analysis program, eight occlusal points were set in the direction of the occlusal long axis. The occlusal load was simulated at 700 N. Results: The zirconia implant (47.7 MPa) von Mises stress decreased by 5.3% in the upper cortical bone compared with the titanium implant (50.2 MPa) von Mises stress. Similarly, the zirconia implant (20.8 MPa) von Mises stress decreased by almost 4% in the cancellous bone compared with the titanium implant (21.7 MPa) von Mises stress. The principal stress in the cortical and cancellous bone exhibited a similar propensity to von Mises stress. Conclusion: In the supporting bone, the zirconia implant is able to reduce bone resorption caused by mechanically transferred stress. It is believed that the zirconia implant can be a potential substitute for the titanium implant by reinforcing aesthetic characteristics and improving stress distribution.

• Nuclear Engineering and Technology
• /
• v.25 no.3
• /
• pp.353-360
• /
• 1993

Many types of pipe whip restraints are installed to protect the structural components from the anticipated pipe whip phenomena of high energy lines in nuclear power plants. It is necessary to investigate these phenomena accurately in order to design the pipe whip restraints properly and/or to evaluate the acceptability of the pipe whip restraint design. Various research programs have been conducted in many countries to develop analytical methods and to verify the validity of the methods. In this study, various types of finite elements in ANSYS[1], the general purpose finite element computer program, was used to simulate the postulated pipe whips to obtain impact loads and the calculated results were compared with the specific experimental results from the sample pipe whip test for the U-shaped pipe whip restraints. Some calculational models, having the gap element or the spring element between the pipe whip restraint and the pipe line, give reasonably good transient responses of the restraint forces compared with the experimental results, and could be useful in evaluating the acceptability of the pipe whip restraint design.

• Coupled systems mechanics
• /
• v.8 no.1
• /
• pp.71-93
• /
• 2019

Finite element simulations of solid mechanics problems often involve the use of Non-Confirming Meshes (NCM) to increase accuracy in capturing nonlinear behavior, including damage and plasticity, in part of a solid domain without an undue increase in computational costs. In the presence of material nonlinearity and plasticity, higher-order variables are often needed to capture nonlinear behavior and material history on non-conforming interfaces. The most popular formulations for coupling non-conforming meshes are dual methods that involve the interpolation of a traction field on the interface. These methods are subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) stability condition, and are therefore limited in their implementation with the higher-order elements needed to capture nonlinear material behavior. Alternatively, the enriched discontinuous Galerkin approach (EDGA) (Haikal and Hjelmstad 2010) is a primal method that provides higher order kinematic fields on the interface, and in which interface tractions are computed from local finite element estimates, therefore facilitating its implementation with nonlinear material models. The inclusion of higher-order interface variables, however, presents the issue of preserving material history at integration points when a increase in integration order is needed. In this study, the enriched discontinuous Galerkin approach (EDGA) is extended to the case of small-deformation plasticity. An interface-driven Gauss-Kronrod integration rule is proposed to enable adaptive enrichment on the interface while preserving history-dependent material data at existing integration points. The method is implemented using classical J2 plasticity theory as well as the pressure-dependent Drucker-Prager material model. We show that an efficient treatment of interface variables can improve algorithmic performance and provide a consistent approach for coupling non-conforming meshes in inelasticity.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.4
• /
• pp.75-86
• /
• 1985

A beam with web opening may reduce the cost of steel and the height of multistory steel buildings. Bower's analysis based on the theory of elasticity and Vierendeel analysis had evaluated the normal stresses around the holes, but these analyses have difficulties for practical uses because of complexity and the limitation for their application. In this study, it is shown that the finite element method, using smaller number of isoparametric elements by taking only a part of the beam which includes the hole, can diminish defects of the above two methods and it may represent more satisfactorily the distribution of the local stress concentration around the hole than the other methods which employed linear elements such as in the analysis by Samuel or Redwood. This study presents the effects of moments, shears, and eccentricities of a hole on the distribution of the normal stresses calculated by using the proposed finite element method. Consequently, it is found that the variations of shear force and hole depth give significant effects on the normal stresses around a hole, while the variations of eccentricities of the hole provide a little effect on them. The regression coefficients resulted from the multiple linear regression may be used for estimating the normal stresses around any arbitrary hole in the web of a beam, since the normal stresses guessed by this regression coefficient equation match well the results by the finite element method except the case of large eccentricity.

• Journal of Electrical Engineering and Technology
• /
• v.12 no.5
• /
• pp.1883-1890
• /
• 2017

Transformers are one of the most precious elements of the electric power system. Stability and reliability of the electric power network mainly depend on the working of the transformer. Leakage reactance of the transformer is one of the important factors and accurate calculation of the leakage reactance is necessary for the transformer designers and electric distributors. Leakage reactance of the transformer depends on the geometry of the transformer. There are many different methods for the calculations of the leakage reactance however mostly are usable when the axial heights of the high voltage and low voltage windings are equal. When the axial heights of high voltage and low voltage windings are asymmetric most of the analytical methods are not reliable. In this study, a new analytical method is introduced for the calculation of the leakage reactance. Fourteen different transformers are investigated in this study and four of them are presented in this paper. The results of the new analytical method are compared with the experimental results. Other analytical and numerical methods are also compared with this new method. Results show that this method is more reliable and accurate as compared to the other analytical methods. The maximum relative error between short-circuit test and proposed method for these fourteen transformers was less than 2.8%.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.12
• /
• pp.2121-2133
• /
• 1997

A finite element code is developed for 3-D self-contact problems, using continuum elements with a SRI(Selective Reduced Integration) scheme to prevent locking phenomenon by the incompressibility of rubber. Contact treatment is carried out in two ways : using the displacement constraints in case of rigid contact ; and imposing the same contact forces on two contact boundaries in case of self-contact. The finite element code developed is applied to the deformation analysis of C.V.joint boots which maintain lubrication conditions and protect the C.V.joint assembly from impact and dust. The boot accompanies large rotation depending on the rotation of the wheel axis and leading to the self-contact phenomena of the boot bellows. Since this contact phenomenon causes wear of the product and has great influence on the endurance life of the product, it is indispensable to carry out stress analysis of the rubber boots. In case of self-contact, various methods for determining contact forces have been suggested with an appropriate contact formulation. Especially, the types of penetration in self-contact are modularized to accelerate conputation with a contact algorithm.

• Nuclear Engineering and Technology
• /
• v.46 no.6
• /
• pp.783-796
• /
• 2014

The coarse mesh finite difference (CMFD) method is applied to the discontinuous finite element method based discrete ordinate calculation for source convergence acceleration. The three-dimensional (3-D) DFEM-Sn code FEDONA is developed for general geometry applications as a framework for the CMFD implementation. Detailed methods for applying the CMFD acceleration are established, such as the method to acquire the coarse mesh flux and current by combining unstructured tetrahedron elements to rectangular coarse mesh geometry, and the alternating calculation method to exchange the updated flux information between the CMFD and DFEM-Sn. The partial current based CMFD (p-CMFD) is also implemented for comparison of the acceleration performance. The modified p-CMFD method is proposed to correct the weakness of the original p-CMFD formulation. The performance of CMFD acceleration is examined first for simple two-dimensional multigroup problems to investigate the effect of the problem and coarse mesh sizes. It is shown that smaller coarse meshes are more effective in the CMFD acceleration and the modified p-CMFD has similar effectiveness as the standard CMFD. The effectiveness of CMFD acceleration is then assessed for three-dimensional benchmark problems such as the IAEA (International Atomic Energy Agency) and C5G7MOX problems. It is demonstrated that a sufficiently converged solution is obtained within 7 outer iterations which would require 175 iterations with the normal DFEM-Sn calculations for the IAEA problem. It is claimed that the CMFD accelerated DFEM-Sn method can be effectively used in the practical eigenvalue calculations involving general geometries.

• Journal of Korean Society of Steel Construction
• /
• v.12 no.3 s.46
• /
• pp.291-302
• /
• 2000

In recent years the development of high modulus, high strength and low density boron and graphite fibers bonded together has brought renewed interestes in structural elements. When a plate with arbitrarily oriented layers and clamped boundary conditions is subjected to uniform loading, it is difficult to analyze and apply, compared with isotropic and orthotropic cases. Therefore the numerical methods, such as finite difference method or finite element method, should be emloyed to analyse such problems. In this study the finite difference technique is used to formulate the bending analysis of symmetric composite laminated skew plates. When this technique is used to solve the problem, it is desirable to reduce the order of the derivatives in order to minimize the number of the pivotal points involved in each equation. The 4th order partial differential equations of laminated skew plates are converted to an equivalent three of 2nd order partial differential equations with three dependant variables.