• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.2113-2122
• /
• 1994

The 3-D FEM analysis for simulating the stamping operation of planar anisotropic sheet metals with arbitrarily-shaped tools is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The consistent full set of governing relations, comprising equilibrium equation and mesh-normal geometric constraints, is appropriately linearized. The linear triangular elements are used for depicting the formed sheet, based on membrane approximation. Barlat's non-quadratic anisotropic yield criterion(strain-rate potential) is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and non-quadratic function parameter. The planar anisotropic finite element formulation is tested with the numerical simulations of the stamping of an automotive hood inner panel and the drawing of a hemispherical punch. The in-plane anisotropic effects on the formability of both mild steel and aluminum alloy sheet metals are examined.

• Journal of Ocean Engineering and Technology
• /
• v.8 no.1
• /
• pp.33-40
• /
• 1994

Welding of structure produces welding residual stresses which influence buckling strength, brittle fracture strength and cold crack on the weld parts. Therefore, it is very important to accurately analyze the residual stress before welding in order to guarantee the safety of weldment. If the weld length is long enough compared to the thickness and the breadth of plate, thermal and mechanical behaviors in the middle portion of the plate are assumed to be uniform along the thickness direction(z-axis). Thus, the following conditions(so-called plane deformation) can be assumed for the plate except near its end;1) distributions of stress and strain are independent on the z-axis;2) plane normal to z-axis before deformation remains plane during and after deformation. In this paper, plane-deformation thermal elasto-plastic problem is formulated by being based on the finite element method. Moreover special regards and paid to the fact that material properties in elastic and plastic region are temperature-dependence. And the method to solve the plane-deformation thermal elasto-plastic problem is shown by using the incremental technique. From the results of analysis, the characterisics of distribution of welding residual stress and plastic strain with the production mechanism are clarified.

• Structural Engineering and Mechanics
• /
• v.10 no.1
• /
• pp.67-79
• /
• 2000

Critical loads and load-carrying capacities for steel scaffolds used as shoring systems were compared using computational and experimental methods in Part I of this paper. In that paper, a simple 2-D model was established for use in evaluating the structural behavior of scaffold-shoring systems. This 2-D model was derived using an incremental finite element analysis (FEA) of a typical complete scaffold-shoring system. Although the simplified model is only two-dimensional, it predicts the critical loads and failure modes of the complete system. The objective of this paper is to present a closed-form solution to the 2-D model. To simplify the analysis, a simpler model was first established to replace the 2-D model. Then, a closed-form solution for the critical loads and failure modes based on this simplified model were derived using a bifurcation (eigenvalue) approach to the elastic-buckling problem. In this closed-form equation, the critical loads are shown to be function of the number of stories, material properties, and section properties of the scaffolds. The critical loads and failure modes obtained from the analytical (closed-form) solution were compared with the results from the 2-D model. The comparisons show that the critical loads from the analytical solution (simplified model) closely match the results from the more complex model, and that the predicted failure modes are nearly identical.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.170-177
• /
• 2003

This paper proposes a robust method for the Ramberg-Osgood (R-O) fit to accurately estimate elastic-plastic J from engineering fracture mechanics analysis based on deformation plasticity. The proposal is based on engineering stress-strain data to determine the R-O parameters, instead of true stress-strain data. Moreover, for practical applications, the method is given not only for the case when full stress-strain data are available but also for the case when only yield and tensile strengths are available. Reliability of the proposed method for the R-O fit is validated against detailed 3-D Finite Element (FE) analyses for circumferential through-wall cracked pipes under global bending using five different materials, three stainless steels and two ferritic steels. Taking the FE J results based on incremental plasticity using actual stress-strain data as reference, the FE J results based on deformation plasticity using various R-O fits are compared with reference J values. Comparisons show that the proposed R-O fit provides more accurate J values for all cases, compared to existing methods for the R-O fit. Advantages of the proposed R-O fit in practical applications are discussed, together with its accuracy.

• Steel and Composite Structures
• /
• v.25 no.5
• /
• pp.557-569
• /
• 2017

Moment-thrust-curvatures ($M-P-{\Phi}$ curves) are fundamental quantities for detailed descriptions of basic properties such as stiffness and strength of a section under axial loads required for accurate computation of the deformations of reinforced concrete or composite columns. Currently, the finite-element-based methods adopting small fibers for analyzing a section are commonly used for generating the $M-P-{\Phi}$ curves and they require large amounts of computational time and effort. Further, the conventional numerical procedure using the force-control method might encounter divergence problems under high compression or tension. Therefore, this paper proposes a divergence-free approach, combining the use of the displacement-control and the Quasi-Newton scheme in the incremental-iterative procedure, for generating the $M-P-{\Phi}$ curves of arbitrary sections. An efficient method for computing the strength from concrete components is employed, where the stress integration is executed by layer-based algorithms. For easy modeling of residual stress, cross sections of structural steel components are meshed into fibers for strength resultants. The numerical procedure is elaborated in detail with flowcharts. Finally, extensive validating examples from previously published research are given for verifying the accuracy of the proposed method.

• The Journal of the Korea Contents Association
• /
• v.15 no.12
• /
• pp.604-611
• /
• 2015

When the moisture in concrete member evaporates by high temperature, the evaporation heat which absorbs surrounding temperature occurs. The incremental rate of the internal temperature in concrete is reduced due to the evaporation heat in spite of continuously increasing external temperature. Therefore, this paper has proposed the evaluation algorithm for predicting the internal temperature of concrete members considering the evaporation heat under the high temperature. Finite element method is employed to facilitate thermal analysis for any position of member. To demonstrate the validity of this numerical procedure, the prediction by the proposed algorithm is compared with the test results of other researchers. The proposed algorithm shows a good agreement with the experimental results including the phenomenon that temperature is lost by the evaporation heat.

• Computational Structural Engineering
• /
• v.11 no.4
• /
• pp.155-168
• /
• 1998

보강된 판 및 쉘구조의 안정성 및 후좌굴을 포함하는 기하학적 비선형 해석을 수행하기 위하여, total Lagrangian formulation에 근거한 연속체의 증분평형방정식으로부터 변형된 쉘요소인 유한요소이론을 제시하였다. 쉘구조의 곡률이 불연속적으로 변하거나 쉘부재들이 유한한 각도로 만나는 보강된 판 및 쉘구조의 비선형 해석이 가능하도록 주부재와 보강재 간의 연결점에 대한 일반적인 변환관계를 제시하였으며 좌굴해석 및 기하학적 비선형해석의 경우에 해의 정확성 및 수렴성을 개선시키기 위하여 접선강도행렬 산정시 회전각의 2차항을 포함시켰다. 또한, shear locking 현상을 극복하기 위하여 감차적분을 적용하였고 쉘구조의 좌굴해석에서는 power method를 적용하여 해석의 효율을 높였으며, 후좌굴해석에서는 변위 및 하중증분법을 적절히 결합시켜 보강된 쉘구조의 후좌굴 거동추적이 용이하였다. 또한, 입력자료를 손쉽게 준비하고 좌굴모드 및 후좌굴거동을 효율적으로 분석하기 위하여 전, 후 처리 프로그램을 개발하였고 다양한 해석예제를 통하여 다른 문헌의 해석결과를 비교함으로써 본 연구에서 개발된 유한요소 해석프로그램의 타당성 및 정확성을 입증하였다.

• Computational Structural Engineering
• /
• v.2 no.3
• /
• pp.97-103
• /
• 1989

In this study, a usefulness of the iterative solution procedures is reviewed for the elasto-plastic large deflection analysis of imperfect plates by finite element method. Three typical solution techniques such as simple incremental(SI) method, Newton-Raphson(NR) method and modified Newton-Raphson (mNR) method are compared. It is concluded that for thin plates which are given rise to the large deflection, iteration for the convergence of the unbalance force should be performed and in this case mNR method is more useful than NR method since the computing time of the former becomes to be a half of the latter, in which the accuracy of the result remains same. For thick plates or thin plates with large initial deflection, however, the use of SI method is quite better since the unbalance force may be negligible.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1571-1578
• /
• 2003

This paper proposes a robust method for the Ramberg-Osgood(R-O)fit to accurately estimate elastic-plastic J from engineering fracture mechanics analysis based on deformation plasticity. The proposal is based on engineering stress-strain data to determine the R-O parameters, instead of true stress-strain data. Moreover, for practical applications, the method is given not only for the case when full stress-strain data are available but also for the case when only yield and tensile strengths are available. Reliability of the proposed method for the R-O fit is validated against detailed 3-D Finite Element (FE) analyses for circumferential through-wall cracked pipes under global bending using five different materials, three stainless steels and two ferritic steels. Taking the FE J results based on incremental plasticity using actual stress-strain data as reference, the FE J results based on deformation plasticity using various R-O fits are compared with reference J values. Comparisons show that the proposed R-O fit provides more accurate J values for all cases, compared to existing methods for the R-O fit. Advantages of the proposed R-O fit in practical applications are discussed, together with its accuracy.

• Proceedings of the KSR Conference
• /
• 2003.05a
• /
• pp.376-381
• /
• 2003

This study is focused on the dynamic response of curved bridge when the rubber tired AGT vehicles is running with alternative articulations. For the analytic approach, there is necessary for the three dimensional vehicle model with 11 degree of freedom and the three dimensional curved bridge model by means of finite element method. It can be described by conventional Lagrangian formula with respect to the dynamic interactions between vehicles and its met bridge. The formula is implemented by Fortran language on the simulation program designated BADIA II(Bridge-AGT Dynamic Interaction Analysis II). The solutions of the formula are derived by Newmark- ${\beta}$ method. The BADIA II is for the dynamic interactions between vehicle and curved bridge in terms of the roughness of running surface and guide rail. The applicability of the BADIA II is verified in terms of displacement and modal frequency. This study is described that the dynamic interactive behaviors between the rubber tired AGT vehicle and curved bridge in terms of the radius of curvatures of curved bridge, vehicle articulations, vehicle speeds, vehicle weights, flatness of running surface and roughness of guide rail using BADIA II.