• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.145-155
• /
• 1996

This paper suggests a new algorithm of automatic mesh generation over planar domains with arbitrarily shaped boundaries and control curves. The algorithm is based on the method of recursively subdividing the domain by the path connecting, with minimum penalty value, two points on the super-loop, which consists of the boundaries and the control curves, The algorithm is not subject to any limitation on the shape of the domain, and its process can be fully automated. Therefore, this algorithm can be implemented into computer programs which require minimal user intervention while generating finite element meshes over complicated domains. This algorithm can also be easily extended for application to the generation of meshes over curved surfaces, or to the adaptive mesh generation.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.9
• /
• pp.117-126
• /
• 1998

For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.

• Nuclear Engineering and Technology
• /
• v.52 no.12
• /
• pp.2918-2927
• /
• 2020

This paper proposes a simplified elastic-plastic analysis procedure using the penalty factors presented in the Code Case N-779 for strain-based fatigue assessment of nuclear safety class 1 components under severe seismic loads such as safety shutdown earthquake and beyond design-basis earthquake. First, a simplified elastic-plastic analysis procedure for strain-based fatigue assessment of nuclear safety class 1 components under the severe seismic loads was proposed based on the analysis result for the simplified elastic-plastic analysis procedure in the Code Case N-779 and the stress categories corresponding to normal operation and seismic loads. Second, total strain amplitude was calculated directly by performing finite element cyclic elastic-plastic seismic analysis for a hot leg nozzle in pressurizer surge line subject to combined loading including deadweight, pressure, seismic inertia load, and seismic anchor motion, as well as was derived indirectly by applying the proposed analysis procedure to the finite element elastic stress analysis result for each load. Third, strain-based fatigue assessment was implemented by applying the strain-based fatigue acceptance criteria in the ASME B&PV Code, Sec. III, Subsec. NB, Article NB-3200 and by using the total strain amplitude values calculated. Last, the total strain amplitude and the fatigue assessment result corresponding to the simplified elastic-plastic analysis were compared with those using the finite element elastic-plastic seismic analysis results. As a result of the comparison, it was identified that the proposed analysis procedure can derive reasonable and conservative results.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.168-175
• /
• 1996

This paper deals with the variable-node element for fluid flow and the adaptive h-version mesh refinement algorithm. The transient element has been formulated by the Galerkin approach in which the pressure term is replaced with the penalty function. The present element having variable mid-side node and is suitable for constructing a locally refined mesh avoiding the use of the highly distorted elements. A modified Gauss quadrature is needed to integrate the element matrices to solve the trouble associated with the discontinuity of derivatives of shape functions. Several numerical examples show that the proposed element can be effectively used in the h-version adapt ive mesh refinement

• Journal of applied mathematics & informatics
• /
• v.4 no.1
• /
• pp.167-178
• /
• 1997

In this paper we present an application to airfoil design of an optimum design method based on optimal control theory. The method used here transforms the design problem by way of a change of variable into an optimal control problem for a distributed system with Neumann boundary control. This results in a set of variational inequalities which is solved by adding a penalty term to the differential equation. This si inturn solved by a finite element method.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.245-248
• /
• 1999

Explicit dynamic finite element analysis has been used widely in the field of sheet metal forming. However in using the analysis technique there are some parameters which are not clearly defined so that engineers may obtain inaccurate solutions In the present study parameters such as time step damping ratio penalty constant and punch speed were investigated on their influence to the solution behavior. Considered forming processes are plane stain bending by a punch and axisymmetric deep drawing.

• Transactions of Materials Processing
• /
• v.5 no.4
• /
• pp.297-304
• /
• 1996

A process model is formulated considering the effect of crystallographic texture developed in forming process. The deformation induced plastic anisotropy can be predicted by capturing the evolution of texture during large deformation in the polycrystalline aggregate. The anisotropic stiffness matrix for the aggregate is derived and implemented in Eulerian finite element code using a Consistent Penalty method. As an application the evolution of texture in rolling drawing and extrusion processes are simulated. The numerical results show good agreements with report-ed experimental textures.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.3
• /
• pp.453-458
• /
• 1992

본 연구에서는 유한요소법에 바탕을 둔 최적공정설계법을 이용하여 압출 금형 형상의 최적설계를 실시하였다.설계의 결과를 참고문헌에서 발췌한 이론해와 비교 하였으며, 축대칭 압출공정에서 압하율, 마찰, 재료특성 등이 금형의 최적형상에 미치 는 영향에 대하여 조사하였다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.60-67
• /
• 1998

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Nuclear Engineering and Technology
• /
• v.49 no.4
• /
• pp.838-849
• /
• 2017

This paper proposes penalty factor equations that take into consideration the weld strength over-match given in the classified form similar to the revised equations presented in the Code Case N-779 via cyclic elastic-plastic finite element analysis. It was found that the $K_e$ analysis data reflecting elastic follow-up can be consolidated by normalizing the primary-plus-secondary stress intensity ranges excluding the nonlinear thermal stress intensity component, $S_n$ to over-match degree of yield strength, $M_F$. For the effect of over-match on $K_n{\times}K_{\nu}$, dispersion of the $K_n{\times}K_{\nu}$ analysis data can be sharply reduced by dividing total stress intensity range, excluding local thermal stresses, $S_{p-lt}$ by $M_F$. Finally, the proposed equations were applied to the weld between the safe end and the piping of a pressurizer surge nozzle in pressurized water reactors in order to calculate a cumulative usage factor. The cumulative usage factor was then compared with those derived by the previous $K_e$ factor equations. The result shows that application of the proposed equations can significantly reduce conservatism of fatigue assessment using the previous $K_e$ factor equations.